Differential regulation and interaction of homoeologous WRKY18 and WRKY40 in Arabidopsis allotetraploids and biotic stress responses.

WRKY transcription factors (TFs) belong to a large family of regulatory proteins in plants that modulate many plant processes. Extensive studies have been conducted on WRKY-mediated defense response in Arabidopsis thaliana and several crop species. Here, we aimed to investigate the potential roles and contributions of WRKY TFs in improving the defense response in the resynthesized Arabidopsis allotetraploids (Arabidopsis suecica) derived from two related autotetraploid progenitors, Arabidopsis thaliana (At4) and Arabidopsis arenosa (Aa). Rapid and differential induction of WRKY18 and WRKY40 expression was evident in response to Pseudomonas syringae and salicylic acid (SA) treatments in the allotetraploids. Selected direct targets of the WRKYs and PR1 also showed altered induction kinetics in the allotetraploids. Cleaved amplified polymorphic sequence analysis further revealed the accumulation of preferential homoeologous alleles (AtWRKY18, AaWRKY40, and AtWRKY60) in the allotetraploids, suggesting the potential for altered protein-protein interaction networks in the hybrids. Indeed, results showed that the cis-interacting AtWRKY18/AtWRKY18 homodimer or trans-interacting AtWRKY18/AaWRKY40 heterodimer exists as the preferred dimer interaction. Moreover, differential affinities of WRKY18 and WRKY40 homo- and heterodimers toward the W-boxes in the WRKY60 promoter were observed. Transient and stable expression of the selected WRKYs in transgenic Arabidopsis further supported the idea that differential interactions lead to changes in PR1 induction and direct target expression under stress, respectively. Our data suggest that differential expression as well as differences in the strength of protein-protein and/or protein-DNA interactions among the WRKY homoeologs could lead to altered regulatory networks of defense genes, contributing to improved defense in allotetraploids.

Regulating the Regulators: The Control of Transcription Factors in Plant Defense Signaling.

Being sessile, plants rely on intricate signaling pathways to mount an efficient defense against external threats while maintaining the cost balance for growth. Transcription factors (TFs) form a repertoire of master regulators in controlling various processes of plant development and responses against external stimuli. There are about 58 families of TFs in plants and among them, six major TF families (AP2/ERF (APETALA2/ethylene responsive factor), bHLH (basic helix-loop-helix), MYB (myeloblastosis related), NAC (no apical meristem (NAM), Arabidopsis transcription activation factor (ATAF1/2), and cup-shaped cotyledon (CUC2)), WRKY, and bZIP (basic leucine zipper)) are found to be involved in biotic and abiotic stress responses. As master regulators of plant defense, the expression and activities of these TFs are subjected to various transcriptional and post-transcriptional controls, as well as post-translational modifications. Many excellent reviews have discussed the importance of these TFs families in mediating their downstream target signaling pathways in plant defense. In this review, we summarize the molecular regulatory mechanisms determining the expression and activities of these master regulators themselves, providing insights for studying their variation and regulation in crop wild relatives (CWR). With the advance of genome sequencing and the growing collection of re-sequencing data of CWR, now is the time to re-examine and discover CWR for the lost or alternative alleles of TFs. Such approach will facilitate molecular breeding and genetic improvement of domesticated crops, especially in stress tolerance and defense responses, with the aim to address the growing concern of climate change and its impact on agriculture crop production.

A Role for CHH Methylation in the Parent-of-Origin Effect on Altered Circadian Rhythms and Biomass Heterosis in Arabidopsis Intraspecific Hybrids.

Hybrid plants and animals often show increased levels of growth and fitness, a phenomenon known as hybrid vigor or heterosis. Circadian rhythms optimize physiology and metabolism in plants and animals. In plant hybrids and polyploids, expression changes of the genes within the circadian regulatory network, such as CIRCADIAN CLOCK ASSOCIATED1 (CCA1), lead to heterosis. However, the relationship between allelic CCA1 expression and heterosis has remained elusive. Here, we show a parent-of-origin effect on altered circadian rhythms and heterosis in Arabidopsis thaliana F1 hybrids. This parent-of-origin effect on biomass heterosis correlates with altered CCA1 expression amplitudes, which are associated with methylation levels of CHH (where H = A, T, or C) sites in the promoter region. The direction of rhythmic expression and hybrid vigor is reversed in reciprocal F1 crosses involving mutants that are defective in the RNA-directed DNA methylation pathway (argonaute4 and nuclear RNA polymerase D1a) but not in the maintenance methylation pathway (methyltransferase1 and decrease in DNA methylation1). This parent-of-origin effect on circadian regulation and heterosis is established during early embryogenesis and maintained throughout growth and development.

Coordinated histone modifications are associated with gene expression variation within and between species.

Histone modifications regulate gene expression in eukaryotes, but their effects on transcriptomes of a multicellular organism and on transcriptomic divergence between species are poorly understood. Here we present the first nucleotide-resolution maps of histone acetylation, methylation, and core histone in Arabidopsis thaliana and a comprehensive analysis of these and all other available maps with gene expression data in A. thaliana, Arabidopsis arenosa, and allotetraploids. H3K9 acetylation (H3K9ac) and H3K4 trimethylation (H3K4me3) are correlated, and their distribution patterns are associated with Gene Ontology (GO) functional classifications. Highly dense and narrow distributions of these modifications near transcriptional start sites are associated with constitutive expression of genes involved in translation, whereas broad distributions toward coding regions correlate with expression variation of the genes involved in photosynthesis, carbohydrate metabolism, and defense responses. Compared to animal stem cells, dispersed distributions of H3K27me3 without bivalent H3K4me3 and H3K9ac marks correlate with developmentally repressed genes in Arabidopsis. Finally, genes affected by A. thaliana histone deacetylase 1 mutation tend to show high levels of expression variation within and between species. The data suggest that genome-wide coordinated modifications of histone acetylation and methylation provide a general mechanism for gene expression changes within and between species and in allopolyploids.

cis- and trans-Regulation of miR163 and target genes confers natural variation of secondary metabolites in two Arabidopsis species and their allopolyploids.

MicroRNAs (miRNAs) play essential roles in plant and animal development, but the cause and effect of miRNA expression divergence between closely related species and in interspecific hybrids or allopolyploids are unknown. Here, we show differential regulation of a miR163-mediated pathway in allotetraploids and their progenitors, Arabidopsis thaliana and Arabidopsis arenosa. miR163 is a recently evolved miRNA in A. thaliana and highly expressed in A. thaliana, but its expression was undetectable in A. arenosa and repressed in resynthesized allotetraploids. Repression of A. arenosa MIR163 (Aa MIR163) is caused by a weak cis-acting promoter and putative trans-acting repressor(s) present in A. arenosa and allotetraploids. Moreover, ectopic Aa MIR163 precursors were processed more efficiently in A. thaliana than in resynthesized allotetraploids, suggesting a role of posttranscriptional regulation in mature miR163 abundance. Target genes of miR163 encode a family of small molecule methyltransferases involved in secondary metabolite biosynthetic pathways that are inducible by a fungal elicitor, alamethicin. Loss of miR163 or overexpression of miR163 in mir163 mutant plants alters target transcript and secondary metabolite profiles. We suggest that cis- and trans-regulation of miRNA and other genes provides a molecular basis for natural variation of biochemical and metabolic pathways that are important to growth vigor and stress responses in Arabidopsis-related species and allopolyploids.

Roles of target site location and sequence complementarity in trans-acting siRNA formation in Arabidopsis.

In plants, many mRNAs and non-coding RNAs are cleaved by RNA-induced silencing complexes. After cleavage, only a limited number of RNAs are processed into trans-acting siRNAs (tasiRNAs). One reason is that 22 nt small RNAs, but not the more common 21 nt small RNAs, can efficiently trigger tasiRNA formation. The characteristics of the target transcripts may also affect tasiRNA production. Here we report the effects of target site location and sequence complementarity on tasiRNA formation. A synthetic sequence that included a miR173 target site and two siRNAs targeting an endogenous mRNA encoding PHYTOENE DESATURASE3 was introduced into a protein-coding (GFP) gene in the coding region or 3' UTR. tasiRNAs were generated in the transgenic seedlings, and the PDS3 mRNA level was reduced, leading to a photobleaching phenotype. It was found that tasiRNAs were most efficiently produced when the miR173 target site was placed immediately after the stop codon. Introducing premature stop codons caused a dramatic reduction of tasiRNAs and over-accumulation of 3' cleavage products, suggesting positive effects of translation on processing the 3' cleavage products into tasiRNAs. By systematically mutating the miR173 target site, we found that perfect complementarity between the 3' end of miR173 and the 5' end of the target sequence was crucial. Mismatches at that position abolished tasiRNA formation, but mismatches at the 5' end of miR173 had less effect. These data suggest important roles for translation and specific sequence complementarity in tasiRNA formation, providing new insights into tasiRNA biogenesis as well as a strategy for improving the efficiency of RNA interference (RNAi) using tasiRNAs.

Production of the 42-kDa fragment of Plasmodium falciparum merozoite surface protein 1, a leading malaria vaccine antigen, in Arabidopsis thaliana seeds.

Malaria is widely associated with poverty, and a low-cost vaccine against malaria is highly desirable for implementing comprehensive vaccination programmes in developing countries. Production of malaria antigens in plants is a promising approach, but its development has been hindered by poor expression of the antigens in plant cells. In the present study, we targeted plant seeds as a low-cost vaccine production platform and successfully expressed the Plasmodium falciparum 42-kDa fragment of merozoite surface protein 1 (MSP142), a leading malaria vaccine candidate, at a high level in transgenic Arabidopsis seeds. We overcame hurdles of transcript and protein instabilities of MSP142 in plants by synthesizing a plant-optimized MSP142 cDNA and either targeting the recombinant protein to protein storage vacuoles or fusing it with a stable plant storage protein. An exceptional improvement in MSP142 expression, from an undetectable level to 5% of total extractable protein, was achieved with these combined strategies. Importantly, the plant-derived MSP142 maintains its natural antigenicity and can be recognized by immune sera from malaria-infected patients. Our results provide a strong basis for the development of a plant-based, low-cost malaria vaccine.

RNAi-mediated down-regulation of DCL1 and AGO1 induces developmental changes in resynthesized Arabidopsis allotetraploids.

Both natural and newly synthesized allopolyploids display nonadditive gene expression changes through genetic and epigenetic mechanisms. The nonadditively expressed genes include many microRNA (miRNA) targets, suggesting a role for miRNAs and their targets in morphological variation in the allopolyploids and their progenitors. We produced dominant-negative transgenic allotetraploid plants in Arabidopsis using RNA interference (RNAi) that downregulates the expression of miRNA biogenesis genes, including DCL1 and AGO1. RNAi of DCL1 and AGO1 led to dominant negative phenotypes and decreased accumulation of several miRNAs and a tasiRNA tested in the transgenic resynthesized allotetraploids. The results demonstrated that miRNA biogenesis genes are effectively downregulated in the resynthesized allotetraploids containing redundant homoeologous genes that are difficult to be manipulated by conventional mutation screens. These lines will be useful for studying the effects of miRNA biogenesis genes on growth and developmental variation in the allopolyploids.

Plant SET domain-containing proteins: structure, function and regulation.

Modification of the histone proteins that form the core around which chromosomal DNA is looped has profound epigenetic effects on the accessibility of the associated DNA for transcription, replication and repair. The SET domain is now recognized as generally having methyltransferase activity targeted to specific lysine residues of histone H3 or H4. There is considerable sequence conservation within the SET domain and within its flanking regions. Previous reviews have shown that SET proteins from Arabidopsis and maize fall into five classes according to their sequence and domain architectures. These classes generally reflect specificity for a particular substrate. SET proteins from rice were found to fall into similar groupings, strengthening the merit of the approach taken. Two additional classes, VI and VII, were established that include proteins with truncated/interrupted SET domains. Diverse mechanisms are involved in shaping the function and regulation of SET proteins. These include protein-protein interactions through both intra- and inter-molecular associations that are important in plant developmental processes, such as flowering time control and embryogenesis. Alternative splicing that can result in the generation of two to several different transcript isoforms is now known to be widespread. An exciting and tantalizing question is whether, or how, this alternative splicing affects gene function. For example, it is conceivable that one isoform may debilitate methyltransferase function whereas the other may enhance it, providing an opportunity for differential regulation. The review concludes with the speculation that modulation of SET protein function is mediated by antisense or sense-antisense RNA.

Regulation of miR163 and its targets in defense against Pseudomonas syringae in Arabidopsis thaliana.

Small RNAs are important regulators for a variety of biological processes, including leaf development, flowering-time, embryogenesis and defense responses. miR163 is a non-conserved miRNA and its locus has evolved recently through inverted duplication of its target genes to which they belong to the SABATH family of related small-molecule methyltransferases (MTs). In Arabidopsis thaliana, previous study demonstrated that miR163 accumulation was induced by alamethicin treatment, suggesting its roles in defense response pathways. Enhanced resistance against Pseudomonas syringae pv. tomato (Pst) was observed in the mir163 mutant, whereas transgenic lines overexpressing miR163 showed increase sensitivity to Pst, suggesting that miR163 is a negative regulator of defense response. Elevated level of miR163 and its targets in A. thaliana were observed upon Pst treatment, suggesting a modulating relationship between miR163 and its targets. In addition, miR163 and histone deacetylase were found to act cooperatively in mediating defense against Pst. Transgenic plants overexpressing miR163-resistant targets suggested their different contributions in defense. Results from this study revealed that the stress-inducible miR163 and its targets act in concert to modulate defense responses against bacterial pathogen in A. thaliana.

Heterologous protein-DNA interactions lead to biased allelic expression of circadian clock genes in interspecific hybrids.

Genomic interactions in allopolyploids create expression variation of homoeologous alleles through protein-protein and protein-DNA interactions. However, the molecular basis for this is largely unknown. Here we investigated the protein-protein and protein-DNA interactions among homoeologous transcription factors in the circadian-clock feedback loop, consisting of CCA1 HIKING EXPEDITION (CHE), CIRCADIAN CLOCK ASSOCIATED1 (CCA1), and TIMING OF CAB EXPRESSION1 (TOC1), plus the interaction with a chromatin factor, HISTONE DEACETYLASE1 (HD1). In the allotetraploids formed between A. thaliana (At) and Arabidopsis arenosa (Aa), AtCCA1 is expressed at lower levels than AaCCA1, which could alter clock output traits. The reduced AtCCA1 expressions in the allotetraploids are consistent with the biochemical data that AaCHE showed preferential binding to the AtCCA1 promoter, in which AaCHE interacts with a higher affinity to AtHD1 than AtCHE. AaCHE also showed a higher affinity to TOC1 than AtCHE, consistent with the effect of TOC1 on repressing CCA1. Thus, stronger AaCHE-TOC1 and AaCHE-AtHD1 interactions reduce AtCC1 allelic expression. Our current data suggest a biochemical basis for protein interactions in trans with a preference to the cis-acting elements in heterologous combinations to reduce AtCCA1 expression, while altered CCA1 expression has been shown to affect metabolic and biomass heterosis in interspecific hybrids or allotetraploids.

High-throughput RNA-seq for allelic or locus-specific expression analysis in Arabidopsis-related species, hybrids, and allotetraploids.

With the next generation sequencing technology, RNA-Seq (RNA sequencing) becomes one of the most powerful tools in quantification of global transcriptomes, discovery of new transcripts and alternative isoforms, as well as detection of single nucleotide polymorphisms (SNPs). RNA-Seq is advantageous over hybridization-based gene quantification methods: (1) it does not require prior information about genomic sequences, (2) it avoids high background problem caused by cross-hybridization, and (3) it is highly sensitive and avoids background and saturation of signals; and finally it is capable of detecting allelic expression differences in hybrids and allopolyploids. We used the RNA-Seq method to determine the genome-wide transcriptome changes in Arabidopsis allotetraploids and their parents, A. thaliana and A. arenosa. The use of this approach allows us to quantify transcriptome from these species and more importantly, to identify allelic or homoeologous-specific gene expression that plays a role in morphological evolution of allopolyploids. The computational pipelines developed are also applicable to the analysis of chromatin immunoprecipitation sequencing (ChIP-seq) data in Arabidopsis-related species, hybrids, and allopolyploids. Comparative analysis of RNA-Seq and ChIP-Seq data will allow us to determine the effects of chromatin modifications on nonadditive gene expression in hybrids and allopolyploids.

Big roles for small RNAs in polyploidy, hybrid vigor, and hybrid incompatibility.

Small RNAs, including microRNAs (miRNAs), small interfering RNAs (siRNAs), and trans-acting siRNAs (ta-siRNAs), mediate gene expression and epigenetic regulation. While siRNAs are highly diverged, miRNAs and ta-siRNAs are generally conserved but many are differentially expressed between related species and in interspecific hybrids and allopolyploids. On one hand, combination of diverged maternal and paternal siRNAs in the same nucleus may exert cis-acting and trans-acting effects on transposable elements (TEs) and TE-associated genes, leading to genomic instability and endosperm and embryo failures, constituting a bottleneck for the evolution of hybrids and polyploids. On the other hand, cis and trans-acting small RNAs induce quantitative and qualitative changes in epigenetic regulation, leading to morphological variation and hybrid vigor in F1 hybrids and stable allopolyploids as well as transgressive phenotypes in the progeny, increasing a potential for adaptive evolution.

Cis- and trans-regulatory divergence between progenitor species determines gene-expression novelty in Arabidopsis allopolyploids.

Gene-expression divergence between species shapes morphological evolution, but the molecular basis is largely unknown. Here we show cis- and trans-regulatory elements and chromatin modifications on gene-expression diversity in genetically tractable Arabidopsis allotetraploids. In Arabidopsis thaliana and Arabidopsis arenosa, both cis and trans with predominant cis-regulatory effects mediate gene-expression divergence. The majority of genes with both cis- and trans-effects are subjected to compensating interactions and stabilizing selection. Interestingly, cis- and trans-regulation is associated with chromatin modifications. In F1 allotetraploids, Arabidopsis arenosa trans factors predominately affect allelic expression divergence. Arabidopsis arenosa trans factors tend to upregulate Arabidopsis thaliana alleles, whereas Arabidopsis thaliana trans factors up- or down-regulate Arabidopsis arenosa alleles. In resynthesized and natural allotetraploids, trans effects drive expression of both homoeologous loci into the same direction. We provide evidence for natural selection and chromatin regulation in shaping gene-expression diversity during plant evolution and speciation.

Ordered histone modifications are associated with transcriptional poising and activation of the phaseolin promoter.

The phaseolin (phas) promoter drives copious production of transcripts encoding the protein phaseolin during seed embryogenesis but is silent in vegetative tissues, in which a nucleosome is positioned over its three-phased TATA boxes. Transition from the inactive state in transgenic Arabidopsis thaliana leaves was accomplished by ectopic expression of the transcription factor Phaseolus vulgaris ABI3-like factor (ALF) and application of abscisic acid (ABA). Placement of hemagglutinin-tagged ALF expression under the control of an estradiol-inducible promoter permitted chromatin immunoprecipitation analysis of chronological changes in histone modifications, notably increased acetylation of H3-K9 and H4-K12, as phas chromatin was remodeled (potentiated). A different array of changes, including acetylation of H3-K14 and methylation of H3-K4, was found to be associated with ABA-mediated activation. Thus, temporal separation of phas potentiation from activation revealed that histone H3 and H4 Lys residues are not globally hyperacetylated during phas expression. Whereas decreases in histone H3 and H4 levels were detected during ALF-mediated remodeling, slight increases occurred after ABA-mediated activation, suggesting the restoration of histone-phas interactions or the replacement of histones in the phas chromatin. The observed histone modifications provide insight into factors involved in the euchromatinization and activation of a plant gene and expand the evidence for histone code conservation among eukaryotes.

The 5' UTR negatively regulates quantitative and spatial expression from the ABI3 promoter.

The involvement of transcription factors Arabidopsis abscisic acid-insensitive3 (ABI3), maize viviparous1 (VP1) and Phaseolus vulgaris ABI3-like factor (PvALF) in the spatial control of storage protein gene expression is well established. However, little insight exists as to how they are themselves regulated. To address this, a 5.15 kb ABI3 upstream sequence including a 4.6 kb full-length promoter and 519 bp of 5'-untranslated region (UTR) was used to drive either beta-glucuronidase (GUS) or green fluorescent protein (GFP) expression in Arabidopsis. Expression from the full-length (- 4630/ + 519ABI3 ) and various 5'-truncated promoters was detected during embryogenesis in all lines, except those transgenic for promoter elements shorter than 364 bp. Two upstream activating regions, -3600 to -2033 and -2033 to -882, enhanced GUS expression in seeds. The -882 to -364 region was sufficient to confer seed-specific expression of GUS when fused to a - 64/ + 6CaMV 35S minimal promoter. Expression from the ABI3 promoter constructs was seed-specific, except in the presence of exogenous abscisic acid (ABA) (>0.3 microM), when GUS expression was detected in seedling roots. Excision of a 405 bp region containing three upstream open reading frames (uORFs) from the 5'-UTR dramatically increased GUS expression and debilitated constraint of reporter expression in roots. Negative regulation of ABI3 expression by the 5'-UTR may involve a post-transcriptional mechanism analogous to that of tumor suppressor genes which also bear long, uORF-containing, 5'-UTRs, or through interactions with RNA-binding proteins.