SANT proteins modulate gene expression by coordinating histone H3KAc and Khib levels and regulate plant heat tolerance.

Histone post-translational modifications (PTMs), such as acetylation and recently identified lysine 2-hydroxyisobutyrylation (Khib), act as active epigenomic marks in plants. SANT domain-containing proteins SANT1, SANT2, SANT3 and SANT4 (SANT1/2/3/4), derived from PIF/Harbinger transposases, form a complex with HISTONE DEACETYLASE 6 (HDA6) to regulate gene expression via histone deacetylation. However, whether SANT1/2/3/4 coordinate different types of PTMs to regulate transcription and mediate responses to specific stresses in plants remains unclear. Here, in addition to modulating histone deacetylation, we found that SANT1/2/3/4 proteins acted like HDA6 or HDA9 in regulating the removal of histone Khib in Arabidopsis (Arabidopsis thaliana). Histone H3 lysine acetylation (H3KAc) and histone Khib were coordinated by SANT1/2/3/4 to regulate gene expression, with H3KAc playing a predominant role and Khib acting complementarily to H3KAc. SANT1/2/3/4 mutation significantly increased the expression of heat-inducible genes with concurrent change of H3KAc levels under normal and heat stress conditions, resulting in enhanced thermotolerance. This study revealed the critical roles of Harbinger transposon-derived SANT domain-containing proteins in transcriptional regulation by coordinating different types of histone PTMs and in the regulation of plant thermotolerance by mediating histone acetylation modification.

Subtilisin-like proteases from Fusarium graminearum induce plant cell death and contribute to virulence.

Fusarium head blight (FHB), caused by Fusarium graminearum, causes huge annual economic losses in cereal production. To successfully colonize host plants, pathogens secrete hundreds of effectors that interfere with plant immunity and facilitate infection. However, the roles of most secreted effectors of F. graminearum in pathogenesis remain unclear. We analyzed the secreted proteins of F. graminearum and identified 255 candidate effector proteins by liquid chromatography-mass spectrometry (LC-MS). Five subtilisin-like family proteases (FgSLPs) were identified that can induce cell death in Nicotiana benthamiana leaves. Further experiments showed that these FgSLPs induced cell death in cotton (Gossypium barbadense) and Arabidopsis (Arabidopsis thaliana). A signal peptide and light were not essential for the cell death-inducing activity of FgSLPs. The I9 inhibitor domain and the entire C-terminus of FgSLPs were indispensable for their self-processing and cell death-inducing activity. FgSLP-induced cell death occurred independent of the plant signal transduction components BRI-ASSOCIATED KINASE 1 (BAK1), SUPPRESSOR OF BIR1 1 (SOBIR1), ENHANCED DISEASE SUSCEPTIBILITY 1 (EDS1), and PHYTOALEXIN DEFICIENT 4 (PAD4). Reduced virulence was observed when FgSLP1 and FgSLP2 were simultaneously knocked out. This study reveals a class of secreted toxic proteins essential for F. graminearum virulence.

Transcriptome Profiling Reveals the Gene Network Responding to Low Nitrogen Stress in Wheat.

As one of the essential nutrients for plants, nitrogen (N) has a major impact on the yield and quality of wheat worldwide. Due to chemical fertilizer pollution, it has become increasingly important to improve crop yield by increasing N use efficiency (NUE). Therefore, understanding the response mechanisms to low N (LN) stress is essential for the regulation of NUE in wheat. In this study, LN stress significantly accelerated wheat root growth, but inhibited shoot growth. Further transcriptome analysis showed that 8468 differentially expressed genes (DEGs) responded to LN stress. The roots and shoots displayed opposite response patterns, of which the majority of DEGs in roots were up-regulated (66.15%; 2955/4467), but the majority of DEGs in shoots were down-regulated (71.62%; 3274/4565). GO and KEGG analyses showed that nitrate reductase activity, nitrate assimilation, and N metabolism were significantly enriched in both the roots and shoots. Transcription factor (TF) and protein kinase analysis showed that genes such as MYB-related (38/38 genes) may function in a tissue-specific manner to respond to LN stress. Moreover, 20 out of 107 N signaling homologous genes were differentially expressed in wheat. A total of 47 transcriptome datasets were used for weighted gene co-expression network analysis (17,840 genes), and five TFs were identified as the potential hub regulatory genes involved in the response to LN stress in wheat. Our findings provide insight into the functional mechanisms in response to LN stress and five candidate regulatory genes in wheat. These results will provide a basis for further research on promoting NUE in wheat.

Identification of natural allelic variation in TTL1 controlling thermotolerance and grain size by a rice super pan-genome.

Continuously increasing global temperatures present great challenges to food security. Grain size, one of the critical components determining grain yield in rice (Oryza sativa L.), is a prime target for genetic breeding. Thus, there is an immediate need for genetic improvement in rice to maintain grain yield under heat stress. However, quantitative trait loci (QTLs) endowing heat stress tolerance and grain size in rice are extremely rare. Here, we identified a novel negative regulator with pleiotropic effects, Thermo-Tolerance and grain Length 1 (TTL1), from the super pan-genomic and transcriptomic data. Loss-of-function mutations in TTL1 enhanced heat tolerance, and caused an increase in grain size by coordinating cell expansion and proliferation. TTL1 was shown to function as a transcriptional regulator and localized to the nucleus and cell membrane. Furthermore, haplotype analysis showed that hapL and hapS of TTL1 were obviously correlated with variations of thermotolerance and grain size in a core collection of cultivars. Genome evolution analysis of available rice germplasms suggested that TTL1 was selected during domestication of the indica and japonica rice subspecies, but still had much breeding potential for increasing grain length and thermotolerance. These findings provide insights into TTL1 as a novel potential target for the development of high-yield and thermotolerant rice varieties.

Transcriptome Profile of Fusarium graminearum Treated by Putrescine.

Fusarium graminearum (F. graminearum) is the main pathogen of Fusarium head blight (FHB) in wheat, barley, and corn. Deoxynivalenol (DON), produced by F. graminearum, is the most prevalent toxin associated with FHB. The wheat defense compound putrescine can promote DON production during F. graminearum infection. However, the underlying mechanisms of putrescine-induced DON synthesis are not well-studied. To investigate the effect of putrescine on the global transcriptional regulation of F. graminearum, we treated F. graminearum with putrescine and performed RNA deep sequencing. We found that putrescine can largely affect the transcriptome of F. graminearum. Gene ontology (GO) and KEGG enrichment analysis revealed that having a large amount of DEGs was associated with ribosome biogenesis, carboxylic acid metabolism, glycolysis/gluconeogenesis, and amino acid metabolism pathways. Co-expression analysis showed that 327 genes had similar expression patterns to FgTRI genes and were assigned to the same module. In addition, three transcription factor genes were identified as hub genes in this module, indicating that they may play important roles in DON synthesis. These results provide important clues for further analysis of the molecular mechanisms of putrescine-induced DON synthesis and will facilitate the study of the pathogenic mechanisms of FHB.

A domesticated Harbinger transposase forms a complex with HDA6 and promotes histone H3 deacetylation at genes but not TEs in Arabidopsis.

In eukaryotes, histone acetylation is a major modification on histone N-terminal tails that is tightly connected to transcriptional activation. HDA6 is a histone deacetylase involved in the transcriptional regulation of genes and transposable elements (TEs) in Arabidopsis thaliana. HDA6 has been shown to participate in several complexes in plants, including a conserved SIN3 complex. Here, we uncover a novel protein complex containing HDA6, several Harbinger transposon-derived proteins (HHP1, SANT1, SANT2, SANT3, and SANT4), and MBD domain-containing proteins (MBD1, MBD2, and MBD4). We show that mutations of all four SANT genes in the sant-null mutant cause increased expression of the flowering repressors FLC, MAF4, and MAF5, resulting in a late flowering phenotype. Transcriptome deep sequencing reveals that while the SANT proteins and HDA6 regulate the expression of largely overlapping sets of genes, TE silencing is unaffected in sant-null mutants. Our global histone H3 acetylation profiling shows that SANT proteins and HDA6 modulate gene expression through deacetylation. Collectively, our findings suggest that Harbinger transposon-derived SANT domain-containing proteins are required for histone deacetylation and flowering time control in plants.

Chloroplast genomes in Populus (Salicaceae): comparisons from an intensively sampled genus reveal dynamic patterns of evolution.

The chloroplast is one of two organelles containing a separate genome that codes for essential and distinct cellular functions such as photosynthesis. Given the importance of chloroplasts in plant metabolism, the genomic architecture and gene content have been strongly conserved through long periods of time and as such are useful molecular tools for evolutionary inferences. At present, complete chloroplast genomes from over 4000 species have been deposited into publicly accessible databases. Despite the large number of complete chloroplast genomes, comprehensive analyses regarding genome architecture and gene content have not been conducted for many lineages with complete species sampling. In this study, we employed the genus Populus to assess how more comprehensively sampled chloroplast genome analyses can be used in understanding chloroplast evolution in a broadly studied lineage of angiosperms. We conducted comparative analyses across Populus in order to elucidate variation in key genome features such as genome size, gene number, gene content, repeat type and number, SSR (Simple Sequence Repeat) abundance, and boundary positioning between the four main units of the genome. We found that some genome annotations were variable across the genus owing in part from errors in assembly or data checking and from this provided corrected annotations. We also employed complete chloroplast genomes for phylogenetic analyses including the dating of divergence times throughout the genus. Lastly, we utilized re-sequencing data to describe the variations of pan-chloroplast genomes at the population level for P. euphratica. The analyses used in this paper provide a blueprint for the types of analyses that can be conducted with publicly available chloroplast genomes as well as methods for building upon existing datasets to improve evolutionary inference.

Dual Recognition of H3K4me3 and DNA by the ISWI Component ARID5 Regulates the Floral Transition in Arabidopsis.

Chromatin remodeling and histone modifications are important for development and floral transition in plants. However, it is largely unknown whether and how these two epigenetic regulators coordinately regulate the important biological processes. Here, we identified three types of Imitation Switch (ISWI) chromatin-remodeling complexes in Arabidopsis (Arabidopsis thaliana). We found that AT-RICH INTERACTING DOMAIN5 (ARID5), a subunit of a plant-specific ISWI complex, can regulate development and floral transition. The ARID-PHD dual domain cassette of ARID5 recognizes both the H3K4me3 histone mark and AT-rich DNA. We determined the ternary complex structure of the ARID5 ARID-PHD cassette with an H3K4me3 peptide and an AT-containing DNA. The H3K4me3 peptide is combinatorially recognized by the PHD and ARID domains, while the DNA is specifically recognized by the ARID domain. Both PHD and ARID domains are necessary for the association of ARID5 with chromatin. The results suggest that the dual recognition of AT-rich DNA and H3K4me3 by the ARID5 ARID-PHD cassette may facilitate the association of the ISWI complex with specific chromatin regions to regulate development and floral transition.

FHA2 is a plant-specific ISWI subunit responsible for stamen development and plant fertility.

Imitation Switch (ISWI) chromatin remodelers are known to function in diverse multi-subunit complexes in yeast and animals. However, the constitution and function of ISWI complexes in Arabidopsis thaliana remain unclear. In this study, we identified forkhead-associated domain 2 (FHA2) as a plant-specific subunit of an ISWI chromatin-remodeling complex in Arabidopsis. By in vivo and in vitro analyses, we demonstrated that FHA2 directly binds to RLT1 and RLT2, two redundant subunits of the ISWI complex in Arabidopsis. The stamen filament is shorter in the fha2 and rlt1/2 mutants than in the wild type, whereas their pistil lengths are comparable. The shorter filament, which is due to reduced cell size, results in insufficient pollination and reduced fertility. The rlt1/2 mutant shows an early-flowering phenotype, whereas the phenotype is not shared by the fha2 mutant. Consistent with the functional specificity of FHA2, our RNA-seq analysis indicated that the fha2 mutant affects a subset of RLT1/2-regulated genes that does not include genes involved in the regulation of flowering time. This study demonstrates that FHA2 functions as a previously uncharacterized subunit of the Arabidopsis ISWI complex and is exclusively involved in regulating stamen development and plant fertility.

Large-scale identification of expression quantitative trait loci in Arabidopsis reveals novel candidate regulators of immune responses and other processes.

The extensive phenotypic diversity within natural populations of Arabidopsis is associated with differences in gene expression. Transcript levels can be considered as inheritable quantitative traits, and used to map expression quantitative trait loci (eQTL) in genome-wide association studies (GWASs). In order to identify putative genetic determinants for variations in gene expression, we used publicly available genomic and transcript variation data from 665 Arabidopsis accessions and applied the single nucleotide polymorphism-set (Sequence) Kernel Association Test (SKAT) method for the identification of eQTL. Moreover, we used the penalized orthogonal-components regression (POCRE) method to increase the power of statistical tests. Then, gene annotations were used as test units to identify genes that are associated with natural variations in transcript accumulation, which correspond to candidate regulators, some of which may have a broad impact on gene expression. Besides increasing the chances to identify real associations, the analysis using POCRE and SKAT significantly reduced the computational cost required to analyze large datasets. As a proof of concept, we used this approach to identify eQTL that represent novel candidate regulators of immune responses. The versatility of this approach allows its application to any process that is subjected to natural variation among Arabidopsis accessions.

A plant-specific SWR1 chromatin-remodeling complex couples histone H2A.Z deposition with nucleosome sliding.

Deposition of H2A.Z in chromatin is known to be mediated by a conserved SWR1 chromatin-remodeling complex in eukaryotes. However, little is known about whether and how the SWR1 complex cooperates with other chromatin regulators. Using immunoprecipitation followed by mass spectrometry, we found all known components of the Arabidopsis thaliana SWR1 complex and additionally identified the following three classes of previously uncharacterized plant-specific SWR1 components: MBD9, a methyl-CpG-binding domain-containing protein; CHR11 and CHR17 (CHR11/17), ISWI chromatin remodelers responsible for nucleosome sliding; and TRA1a and TRA1b, accessory subunits of the conserved NuA4 histone acetyltransferase complex. MBD9 directly interacts with CHR11/17 and the SWR1 catalytic subunit PIE1, and is responsible for the association of CHR11/17 with the SWR1 complex. MBD9, TRA1a, and TRA1b function as canonical components of the SWR1 complex to mediate H2A.Z deposition. CHR11/17 are not only responsible for nucleosome sliding but also involved in H2A.Z deposition. These results indicate that the association of the SWR1 complex with CHR11/17 may facilitate the coupling of H2A.Z deposition with nucleosome sliding, thereby co-regulating gene expression, development, and flowering time.

Histone acetylation recruits the SWR1 complex to regulate active DNA demethylation in Arabidopsis.

Active DNA demethylation is critical for controlling the DNA methylomes in plants and mammals. However, little is known about how DNA demethylases are recruited to target loci, and the involvement of chromatin marks in this process. Here, we identify 2 components of the SWR1 chromatin-remodeling complex, PIE1 and ARP6, as required for ROS1-mediated DNA demethylation, and discover 2 SWR1-associated bromodomain-containing proteins, AtMBD9 and nuclear protein X1 (NPX1). AtMBD9 and NPX1 recognize histone acetylation marks established by increased DNA methylation 1 (IDM1), a known regulator of DNA demethylation, redundantly facilitating H2A.Z deposition at IDM1 target loci. We show that at some genomic regions, H2A.Z and DNA methylation marks coexist, and H2A.Z physically interacts with ROS1 to regulate DNA demethylation and antisilencing. Our results unveil a mechanism through which DNA demethylases can be recruited to specific target loci exhibiting particular histone marks, providing a conceptual framework to understand how chromatin marks regulate DNA demethylation.

Arabidopsis AGDP1 links H3K9me2 to DNA methylation in heterochromatin.

Heterochromatin is a tightly packed form of chromatin that is associated with DNA methylation and histone 3 lysine 9 methylation (H3K9me). Here, we identify an H3K9me2-binding protein, Agenet domain (AGD)-containing p1 (AGDP1), in Arabidopsis thaliana. Here we find that AGDP1 can specifically recognize the H3K9me2 mark by its three pairs of tandem AGDs. We determine the crystal structure of the Agenet domain 1 and 2 cassette (AGD12) of Raphanus sativus AGDP1 in complex with an H3K9me2 peptide. In the complex, the histone peptide adopts a unique helical conformation. AGD12 specifically recognizes the H3K4me0 and H3K9me2 marks by hydrogen bonding and hydrophobic interactions. In addition, we find that AGDP1 is required for transcriptional silencing, non-CG DNA methylation, and H3K9 dimethylation at some loci. ChIP-seq data show that AGDP1 preferentially occupies long transposons and is associated with heterochromatin marks. Our findings suggest that, as a heterochromatin-binding protein, AGDP1 links H3K9me2 to DNA methylation in heterochromatin regions.

The PEAT protein complexes are required for histone deacetylation and heterochromatin silencing.

In eukaryotes, heterochromatin regions are typically subjected to transcriptional silencing. DNA methylation has an important role in such silencing and has been studied extensively. However, little is known about how methylated heterochromatin regions are subjected to silencing. We conducted a genetic screen and identified an epcr (enhancer of polycomb-related) mutant that releases heterochromatin silencing in Arabidopsis thaliana We demonstrated that EPCR1 functions redundantly with its paralog EPCR2 and interacts with PWWP domain-containing proteins (PWWPs), AT-rich interaction domain-containing proteins (ARIDs), and telomere repeat binding proteins (TRBs), thus forming multiple functionally redundant protein complexes named PEAT (PWWPs-EPCRs-ARIDs-TRBs). The PEAT complexes mediate histone deacetylation and heterochromatin condensation and thereby facilitate heterochromatin silencing. In heterochromatin regions, the production of small interfering RNAs (siRNAs) and DNA methylation is repressed by the PEAT complexes. The study reveals how histone deacetylation, heterochromatin condensation, siRNA production, and DNA methylation interplay with each other and thereby maintain heterochromatin silencing.

The Arabidopsis acetylated histone-binding protein BRAT1 forms a complex with BRP1 and prevents transcriptional silencing.

Transposable elements and other repetitive DNA sequences are usually subject to DNA methylation and transcriptional silencing. However, anti-silencing mechanisms that promote transcription in these regions are not well understood. Here, we describe an anti-silencing factor, Bromodomain and ATPase domain-containing protein 1 (BRAT1), which we identified by a genetic screen in Arabidopsis thaliana. BRAT1 interacts with an ATPase domain-containing protein, BRP1 (BRAT1 Partner 1), and both prevent transcriptional silencing at methylated genomic regions. Although BRAT1 mediates DNA demethylation at a small set of loci targeted by the 5-methylcytosine DNA glycosylase ROS1, the involvement of BRAT1 in anti-silencing is largely independent of DNA demethylation. We also demonstrate that the bromodomain of BRAT1 binds to acetylated histone, which may facilitate the prevention of transcriptional silencing. Thus, BRAT1 represents a potential link between histone acetylation and transcriptional anti-silencing at methylated genomic regions, which may be conserved in eukaryotes.

The SET domain proteins SUVH2 and SUVH9 are required for Pol V occupancy at RNA-directed DNA methylation loci.

RNA-directed DNA methylation (RdDM) is required for transcriptional silencing of transposons and other DNA repeats in Arabidopsis thaliana. Although previous research has demonstrated that the SET domain-containing SU(VAR)3-9 homologs SUVH2 and SUVH9 are involved in the RdDM pathway, the underlying mechanism remains unknown. Our results indicated that SUVH2 and/or SUVH9 not only interact with the chromatin-remodeling complex termed DDR (DMS3, DRD1, and RDM1) but also with the newly characterized complex composed of two conserved Microrchidia (MORC) family proteins, MORC1 and MORC6. The effect of suvh2suvh9 on Pol IV-dependent siRNA accumulation and DNA methylation is comparable to that of the Pol V mutant nrpe1 and the DDR complex mutant dms3, suggesting that SUVH2 and SUVH9 are functionally associated with RdDM. Our CHIP assay demonstrated that SUVH2 and SUVH9 are required for the occupancy of Pol V at RdDM loci and facilitate the production of Pol V-dependent noncoding RNAs. Moreover, SUVH2 and SUVH9 are also involved in the occupancy of DMS3 at RdDM loci. The putative catalytic active site in the SET domain of SUVH2 is dispensable for the function of SUVH2 in RdDM and H3K9 dimethylation. We propose that SUVH2 and SUVH9 bind to methylated DNA and facilitate the recruitment of Pol V to RdDM loci by associating with the DDR complex and the MORC complex.

The PRP6-like splicing factor STA1 is involved in RNA-directed DNA methylation by facilitating the production of Pol V-dependent scaffold RNAs.

DNA methylation is a conserved epigenetic marker in plants and animals. In Arabidopsis, DNA methylation can be established through an RNA-directed DNA methylation (RdDM) pathway. By screening for suppressors of ros1, we identified STA1, a PRP6-like splicing factor, as a new RdDM regulator. Whole-genome bisulfite sequencing suggested that STA1 and the RdDM pathway share a large number of common targets in the Arabidopsis genome. Small RNA deep sequencing demonstrated that STA1 is predominantly involved in the accumulation of the siRNAs that depend on both Pol IV and Pol V. Moreover, the sta1 mutation partially reduces the levels of Pol V-dependent RNA transcripts. Immunolocalization assay indicated that STA1 signals are exclusively present in the Cajal body and overlap with AGO4 in most nuclei. STA1 signals are also partially overlap with NRPE1. Localization of STA1 to AGO4 and NRPE1 signals is probably related to the function of STA1 in the RdDM pathway. Based on these results, we propose that STA1 acts downstream of siRNA biogenesis and facilitates the production of Pol V-dependent RNA transcripts in the RdDM pathway.

DTF1 is a core component of RNA-directed DNA methylation and may assist in the recruitment of Pol IV.

DNA methylation is an important epigenetic mark in many eukaryotic organisms. De novo DNA methylation in plants can be achieved by the RNA-directed DNA methylation (RdDM) pathway, where the plant-specific DNA-dependent RNA polymerase IV (Pol IV) transcribes target sequences to initiate 24-nt siRNA production and action. The putative DNA binding protein DTF1/SHH1 of Arabidopsis has been shown to associate with Pol IV and is required for 24-nt siRNA accumulation and transcriptional silencing at several RdDM target loci. However, the extent and mechanism of DTF1 function in RdDM is unclear. We show here that DTF1 is necessary for the accumulation of the majority of Pol IV-dependent 24-nt siRNAs. It is also required for a large proportion of Pol IV-dependent de novo DNA methylation. Interestingly, there is a group of RdDM target loci where 24-nt siRNA accumulation but not DNA methylation is dependent on DTF1. DTF1 interacts directly with the chromatin remodeling protein CLASSY 1 (CLSY1), and both DTF1 and CLSY1 are associated in vivo with Pol IV but not Pol V, which functions downstream in the RdDM effector complex. DTF1 and DTF2 (a DTF1-like protein) contain a SAWADEE domain, which was found to bind specifically to histone H3 containing H3K9 methylation. Taken together, our results show that DTF1 is a core component of the RdDM pathway, and suggest that DTF1 interacts with CLSY1 to assist in the recruitment of Pol IV to RdDM target loci where H3K9 methylation may be an important feature. Our results also suggest the involvement of DTF1 in an important negative feedback mechanism for DNA methylation at some RdDM target loci.

The splicing machinery promotes RNA-directed DNA methylation and transcriptional silencing in Arabidopsis.

DNA methylation in transposons and other DNA repeats is conserved in plants as well as in animals. In Arabidopsis thaliana, an RNA-directed DNA methylation (RdDM) pathway directs de novo DNA methylation. We performed a forward genetic screen for suppressors of the DNA demethylase mutant ros1 and identified a novel Zinc-finger and OCRE domain-containing Protein 1 (ZOP1) that promotes Pol IV-dependent siRNA accumulation, DNA methylation, and transcriptional silencing. Whole-genome methods disclosed the genome-wide effects of zop1 on Pol IV-dependent siRNA accumulation and DNA methylation, suggesting that ZOP1 has both RdDM-dependent and -independent roles in transcriptional silencing. We demonstrated that ZOP1 is a pre-mRNA splicing factor that associates with several typical components of the splicing machinery as well as with Pol II. Immunofluorescence assay revealed that ZOP1 overlaps with Cajal body and is partially colocalized with NRPE1 and DRM2. Moreover, we found that the other development-defective splicing mutants tested including mac3a3b, mos4, mos12 and mos14 show defects in RdDM and transcriptional silencing. We propose that the splicing machinery rather than specific splicing factors is involved in promoting RdDM and transcriptional silencing.