Erratum: Acetate-mediated novel survival strategy against drought in plants

This corrects the article DOI: 10.1038/nplants.2017.97.

Acetate-mediated novel survival strategy against drought in plants

Water deficit caused by global climate changes seriously endangers the survival of organisms and crop productivity, and increases environmental deterioration^1,2. Plants' resistance to drought involves global reprogramming of transcription, cellular metabolism, hormone signalling and chromatin modification^3-8. However, how these regulatory responses are coordinated via the various pathways, and the underlying mechanisms, are largely unknown. Herein, we report an essential drought-responsive network in which plants trigger a dynamic metabolic flux conversion from glycolysis into acetate synthesis to stimulate the jasmonate (JA) signalling pathway to confer drought tolerance. In Arabidopsis, the ON/OFF switching of this whole network is directly dependent on histone deacetylase HDA6. In addition, exogenous acetic acid promotes de novo JA synthesis and enrichment of histone H4 acetylation, which influences the priming of the JA signalling pathway for plant drought tolerance. This novel acetate function is evolutionarily conserved as a survival strategy against environmental changes in plants. Furthermore, the external application of acetic acid successfully enhanced the drought tolerance in Arabidopsis, rapeseed, maize, rice and wheat plants. Our findings highlight a radically new survival strategy that exploits an epigenetic switch of metabolic flux conversion and hormone signalling by which plants adapt to drought.

The expression of miR-125b-5p is increased in the serum of patients with chronic hepatitis B infection and inhibits the detection of hepatitis B virus surface antigen

MicroRNAs were first discovered as small endogenous RNA molecules and some viruses have been reported to interact with host miRNAs. By investigating miRNA expression in serum derived from HBV-infected patients, we have clarified the relationship between miRNA expression and chronic HBV infection. Additionally, we demonstrate the use of miRNAs as both novel biomarkers and new therapies against HBV. We included the sera of 20 patients with chronic HBV infection, sera of 20 patients with HCV infection and sera of 10 healthy controls in this study. The miRNA libraries were sequenced using a 32-mer single end sequence. The validation study of circulating miRNA in serum was conducted by qRT-PCR. The HBV genomic regions of genotype B and genotype C that were speculated to be targeted by miRNA were constructed using complementary oligonucleotides in the vectors. Reporter assays were performed 48 h after transfection. The expression levels of 21 miRNAs were found to be differentially expressed in the three groups. 10 miRNAs (hsa-miR-100-5p, miR-125b-5p, miR-193b-3p, miR-194-3p, miR-30a-3p, miR-30c-2-3p, miR-3591-5p, miR-4709-3p, miR-574-3p and miR-99a-5p) were found to be upregulated in CH-B by deep sequence analysis. The computer analysis showed that two regions of HBsAg are potential targets of miR-125b-5p and miR-30c-2-3p and that these miRNAs may downregulate the expression of HBV-S. The HBV genotype C segment speculated to be targeted by hsa-miR-125b-5p significantly decreased the expression of the reporter. This study indicated that expression of miR-125b-5p was related to the etiology of chronic hepatitis B infection and regulated the expression of HBsAg.

Sm-Like Protein-Mediated RNA Metabolism Is Required for Heat Stress Tolerance in Arabidopsis

Sm-like proteins play multiple functions in RNA metabolism, which is essential for biological processes such as stress responses in eukaryotes. The Arabidopsis thaliana sad1 mutant has a mutation of sm-like protein 5 (LSM5) and shows impaired drought and salt stress tolerances. The lsm5/sad1 mutant also showed hypersensitivity to heat stress. GFP-fused LSM5/SAD1 was localized in the nucleus under optimal growth conditions. After heat stress treatment, GFP-fused LSM5/SAD1 fluorescence was also observed as small cytoplasmic dots, in addition to nuclear localization. Whole genome transcriptome analysis revealed that many genes in Arabidopsis were drastically changed in response to heat stress. More heat-responsive genes were highly expressed in lsm5/sad1 mutant at both 2 and 6 h after heat stress treatment. Additionally, intron-retained and capped transcripts accumulated in the lsm5/sad1 mutant after heat stress treatment. In this study, we also identified non-Arabidopsis Genome Initiative transcripts that were expressed from unannotated regions. Most of these transcripts were antisense transcripts, and many capped non-AGI transcripts accumulated in the lsm5/sad1 mutant during heat stress treatment. These results indicated that LSM5/SAD1 functions to degrade aberrant transcripts through appropriate mRNA splicing and decapping, and precise RNA metabolic machinery is required for heat stress tolerance.

Analysis of differential expression patterns of mRNA and protein during cold-acclimation and de-acclimation in Arabidopsis

Overwintering plants are capable of exhibiting high levels of cold tolerance, which is acquired through the process of cold acclimation (CA). In contrast to CA, the acquired freezing tolerance is rapidly reduced during cold de-acclimation (DA) and plants resume growth after sensing warm temperatures. In order to better understand plant growth and development, and to aid in the breeding of cold-tolerant plants, it is important to decipher the functional mechanisms of the DA process. In this study, we performed comparative transcriptomic and proteomic analyses during CA and DA. As revealed by shotgun proteomics, we identified 3987 peptides originating from 1569 unique proteins and the corresponding mRNAs were analyzed. Among the 1569 genes, 658 genes were specifically induced at the transcriptional level during the process of cold acclimation. In order to investigate the relationship between mRNA and the corresponding protein expression pattern, a Pearson correlation was analyzed. Interestingly, 199 genes showed a positive correlation of mRNA and protein expression pattern, indicating that both their transcription and translation occurred during CA. However, 226 genes showed a negative correlation of mRNA and protein expression pattern, indicating that their mRNAs were transcribed during CA and were stored for the subsequent DA step. Under this scenario, those proteins were specifically increased during DA without additional transcription of mRNA. In order to confirm the negative correlation of mRNA and protein expression patterns, qRT-PCR and western blot analyses were performed. Mitochondrial malate dehydrogenase 1 (mMDH1) exhibited a negative correlation of mRNA and protein levels, which was characterized by CA-specific mRNA induction and protein accumulation specifically during DA. These data indicate that the expression of specific mRNAs and subsequent accumulation of corresponding proteins are not always in accordance under low temperature stress conditions in plants.

Transcriptome analysis using a high-density oligomicroarray under drought stress in various genotypes of cassava: an important tropical crop

Cassava is an important crop that provides food security and income generation in many tropical countries and is known for its adaptability to various environmental conditions. Despite its global importance, the development of cassava microarray tools has not been well established. Here, we describe the development of a 60-mer oligonucleotide Agilent microarray representing ∼20,000 cassava genes and how it can be applied to expression profiling under drought stress using three cassava genotypes (MTAI16, MECU72 and MPER417-003). Our results identified about 1300 drought stress up-regulated genes in cassava and indicated that cassava has similar mechanisms for drought stress response and tolerance as other plant species. These results demonstrate that our microarray is a useful tool for analysing the cassava transcriptome and that it is applicable for various cassava genotypes.

Arabidopsis HDA6 regulates locus-directed heterochromatin silencing in cooperation with MET1

Heterochromatin silencing is pivotal for genome stability in eukaryotes. In Arabidopsis, a plant-specific mechanism called RNA–directed DNA methylation (RdDM) is involved in heterochromatin silencing. Histone deacetylase HDA6 has been identified as a component of such machineries; however, its endogenous targets and the silencing mechanisms have not been analyzed globally. In this study, we investigated the silencing mechanism mediated by HDA6. Genome-wide transcript profiling revealed that the loci silenced by HDA6 carried sequences corresponding to the RDR2-dependent 24-nt siRNAs, however their transcript levels were mostly unaffected in the rdr2 mutant. Strikingly, we observed significant overlap of genes silenced by HDA6 to those by the CG DNA methyltransferase MET1. Furthermore, regardless of dependence on RdDM pathway, HDA6 deficiency resulted in loss of heterochromatic epigenetic marks and aberrant enrichment for euchromatic marks at HDA6 direct targets, along with ectopic expression of these loci. Acetylation levels increased significantly in the hda6 mutant at all of the lysine residues in the H3 and H4 N-tails, except H4K16. Interestingly, we observed two different CG methylation statuses in the hda6 mutant. CG methylation was sustained in the hda6 mutant at some HDA6 target loci that were surrounded by flanking DNA–methylated regions. In contrast, complete loss of CG methylation occurred in the hda6 mutant at the HDA6 target loci that were isolated from flanking DNA methylation. Regardless of CG methylation status, CHG and CHH methylation were lost and transcriptional derepression occurred in the hda6 mutant. Furthermore, we show that HDA6 binds only to its target loci, not the flanking methylated DNA, indicating the profound target specificity of HDA6. We propose that HDA6 regulates locus-directed heterochromatin silencing in cooperation with MET1, possibly recruiting MET1 to specific loci, thus forming the foundation of silent chromatin structure for subsequent non-CG methylation.

Eukaryotes are defended from potentially harmful DNA elements, such as transposons, by forming inactive genomic structure. Chromatin, which consists of DNA and histone proteins, is densely packed in the silent structure, and chromatin chemical modifications such as DNA methylation and histone modifications are known to be essential for this packing. In plants, small RNA molecules have been thought to trigger DNA methylation and resulting silent chromatin formation. We revealed that elimination of specific histone modifications concomitant with DNA methylation is pivotal for the silent chromatin. Furthermore, the histone deacetylase was shown to have more profound target specificity than the DNA methyltransferase and is required for locus-directed DNA methylation, implying the involvement of the histone deacetylase for targeting the DNA methyltransferase to specific places on the genome. These proteins and their functions for gene silencing are evolutionarily conserved in higher eukaryotes, and several proteins involved in small RNA production are plant-specific. Thus, we present a hypothesis that the plant genome may build the protecting foundation by the conserved genome surveillance in eukaryotes, and the reinforcing machinery involving small RNAs could be evolutionarily added to the plant heterochromatin silencing system.

ARTADE2DB: improved statistical inferences for Arabidopsis gene functions and structure predictions by dynamic structure-based dynamic expression (DSDE) analyses

Recent advances in technologies for observing high-resolution genomic activities, such as whole-genome tiling arrays and high-throughput sequencers, provide detailed information for understanding genome functions. However, the functions of 50% of known Arabidopsis thaliana genes remain unknown or are annotated only on the basis of static analyses such as protein motifs or similarities. In this paper, we describe dynamic structure-based dynamic expression (DSDE) analysis, which sequentially predicts both structural and functional features of transcripts. We show that DSDE analysis inferred gene functions 12% more precisely than static structure-based dynamic expression (SSDE) analysis or conventional co-expression analysis based on previously determined gene structures of A. thaliana. This result suggests that more precise structural information than the fixed conventional annotated structures is crucial for co-expression analysis in systems biology of transcriptional regulation and dynamics. Our DSDE method, ARabidopsis Tiling-Array-based Detection of Exons version 2 and over-representation analysis (ARTADE2-ORA), precisely predicts each gene structure by combining two statistical analyses: a probe-wise co-expression analysis of multiple transcriptome measurements and a Markov model analysis of genome sequences. ARTADE2-ORA successfully identified the true functions of about 90% of functionally annotated genes, inferred the functions of 98% of functionally unknown genes and predicted 1,489 new gene structures and functions. We developed a database ARTADE2DB that integrates not only the information predicted by ARTADE2-ORA but also annotations and other functional information, such as phenotypes and literature citations, and is expected to contribute to the study of the functional genomics of A. thaliana. URL: http://artade.org.

Genome-wide analysis of endogenous abscisic acid-mediated transcription in dry and imbibed seeds of Arabidopsis using tiling arrays

The phytohormone abscisic acid (ABA) plays important roles in the induction and maintenance of seed dormancy. Although application of exogenous ABA inhibits germination, the effects of exogenous ABA on ABA-mediated gene transcription differ from those of endogenous ABA. To understand how endogenous ABA regulates the transcriptomes in seeds, we performed comprehensive expression analyses using whole-genome Affymetrix tiling arrays in two ABA metabolism mutants - an ABA-deficient mutant (aba2) and an ABA over-accumulation mutant (cyp707a1a2a3 triple mutant). Hierarchical clustering and principal components analyses showed that differences in endogenous ABA levels do not influence global expression of stored mRNA in dry seeds. However, the transcriptome after seed imbibition was related to endogenous ABA levels in both types of mutant. Endogenous ABA-regulated genes expressed in imbibed seeds included those encoding key ABA signaling factors and gibberellin-related components. In addition, cohorts of ABA-upregulated genes partially resembled those of dormant genes, whereas ABA-downregulated genes were partially overlapped with after-ripening-regulated genes. Bioinformatic analyses revealed that 6105 novel genes [non-Arabidopsis Genome Initiative (AGI) transcriptional units (TUs)] were expressed from unannotated regions. Interestingly, approximately 97% of non-AGI TUs possibly encoded hypothetical non-protein-coding RNAs, including a large number of antisense RNAs. In dry and imbibed seeds, global expression profiles of non-AGI TUs were similar to those of AGI genes. For both non-AGI TUs and AGI code genes, we identified those that were regulated differently in embryo and endosperm tissues. Our results suggest that transcription in Arabidopsis seeds is more complex and dynamic than previously thought.

Transcriptome analyses revealed diverse expression changes in ago1 and hyl1 Arabidopsis mutants

MicroRNAs (miRNAs) are 20-24 nucleotide endogenous regulatory molecules conserved in higher eukaryotes. In Arabidopsis, miRNAs are produced through step-wise cleavages of primary miRNA precursors (pri-miRNAs) by DICER-LIKE1 (DCL1). This cleavage step is also supported by a double-stranded RNA-binding protein, HYPONASTIC LEAVES1 (HYL1). In many cases, mature miRNA is predominantly incorporated into an endonuclease, ARGONAUTE1 (AGO1), which degrades miRNA-targeted mRNAs. Here, we examined and revealed whole genome transcriptomes in ago1-25 and hyl1-2 mutants using tiling arrays. The data in this paper are valuable for understanding the relationship between the miRNA pathway and its effect on transcriptomes.

Arabidopsis transcriptome analysis under drought, cold, high-salinity and ABA treatment conditions using a tiling array

Plants respond and adapt to drought, cold and high-salinity stresses in order to survive. In this study, we applied Arabidopsis Affymetrix tiling arrays to study the whole genome transcriptome under drought, cold, high-salinity and ABA treatment conditions. The bioinformatic analysis using the tiling array data showed that 7,719 non-AGI transcriptional units (TUs) exist in the unannotated "intergenic" regions of Arabidopsis genome. These include 1,275 and 181 TUs that are induced and downregulated, respectively, by the stress or ABA treatments. Most of the non-AGI TUs are hypothetical non-protein-coding RNAs. About 80% of the non-AGI TUs belong to pairs of the fully overlapping sense-antisense transcripts (fSATs). Significant linear correlation between the expression ratios (treated/untreated) of the sense TUs and the ratios of the antisense TUs was observed in the SATs of AGI code/non-AGI TU. We studied the biogenesis mechanisms of the stress- or ABA-inducible antisense RNAs and found that the expression of sense TUs is necessary for the stress- or ABA-inducible expression of the antisense TUs in the fSATs (AGI code/non-AGI TU).