Plants develop tolerance to drought by activating genes with altered levels of epigenetic modifications. Specific transcription factors are involved in this activation, but the molecular connections within the regulatory system are unclear. Here, we analyzed genome-wide acetylated lysine residue 9 of histone H3 (H3K9ac) enrichment and examined its association with transcriptomes in Populus trichocarpa under drought stress. We revealed that abscisic acid-Responsive Element (ABRE) motifs in promoters of the drought-responsive genes PtrNAC006, PtrNAC007, and PtrNAC120 are involved in H3K9ac enhancement and activation of these genes. Overexpressing these PtrNAC genes in P. trichocarpa resulted in strong drought-tolerance phenotypes. We showed that the ABRE binding protein PtrAREB1-2 binds to ABRE motifs associated with these PtrNAC genes and recruits the histone acetyltransferase unit ADA2b-GCN5, forming AREB1-ADA2b-GCN5 ternary protein complexes. Moreover, this recruitment enables GCN5-mediated histone acetylation to enhance H3K9ac and enrich RNA polymerase II specifically at these PtrNAC genes for the development of drought tolerance. CRISPR editing or RNA interference-mediated downregulation of any of the ternary members results in highly drought-sensitive P. trichocarpa. Thus, the combinatorial function of the ternary proteins establishes a coordinated histone acetylation and transcription factor-mediated gene activation for drought response and tolerance in Populus species.

Wood remains the world’s most abundant and renewable resource for timber and pulp and is an alternative to fossil fuels. Understanding the molecular regulation of wood formation can advance the engineering of wood for more efficient material and energy productions. We integrated a black cottonwood (Populus trichocarpa) wood-forming cell system with quantitative transcriptomics and chromatin binding assays to construct a transcriptional regulatory network (TRN) directed by a key transcription factor (TF), PtrSND1-B1 (secondary wall-associated NAC-domain protein). The network consists of four layers of TF–target gene interactions with quantitative regulatory effects, describing the specificity of how the regulation is transduced through these interactions to activate cell wall genes (effector genes) for wood formation. PtrSND1-B1 directs 57 TF–DNA interactions through 17 TFs transregulating 27 effector genes. Of the 57 interactions, 55 are novel. We tested 42 of these 57 interactions in 30 genotypes of transgenic P. trichocarpa and verified that ∼90% of the tested interactions function in vivo. The TRN reveals common transregulatory targets for distinct TFs, leading to the discovery of nine TF protein complexes (dimers and trimers) implicated in regulating the biosynthesis of specific types of lignin. Our work suggests that wood formation may involve regulatory homeostasis determined by combinations of TF–DNA and TF–TF (protein–protein) regulations.

Bacterial spot (BS) is one of the most devastating foliar bacterial diseases of tomato and is caused by multiple species of Xanthomonas. We performed the RNA sequencing (RNA-Seq) analysis of three tomato lines with different levels of resistance to Xanthomonas perforans race T4 to study the differentially expressed genes (DEGs) and transcript-based sequence variations. Analysis between inoculated and control samples revealed that resistant genotype Solanum pimpinellifolium accession PI 270443 had more DEGs (834), followed by susceptible genotype tomato (S. lycopersicum L) breeding line NC 714 (373), and intermediate genotype tomato breeding line NC 1CELBR (154). Gene ontology (GO) terms revealed that more GO terms (51) were enriched for upregulated DEGs in the resistant genotype PI 270443, and more downregulated DEGs (67) were enriched in the susceptible genotype NC 714. DEGs in the biotic stress pathway showed more upregulated biotic stress pathway DEGs (67) for PI 270443 compared to more downregulated DEGs (125) for the susceptible NC 714 genotype. Resistant genotype PI 270443 has three upregulated DEGs for pathogenesis-related (PR) proteins, and susceptible genotype NC 714 has one downregulated R gene. Sequence variations called from RNA-Seq reads against the reference genome of susceptible Heinz 1706 showed that chr11, which has multiple reported resistance quantitative trait loci (QTLs) to BS race T4, is identical between two resistant lines, PI 270443 and NC 1CELBR, suggesting that these two lines share the same resistance QTLs on this chromosome. Several loci for PR resistance proteins with sequence variation between the resistant and susceptible tomato lines were near the known Rx4 resistance gene on chr11, and additional biotic stress associated DEGs near to the known Rx4 resistance gene were also identified from the susceptible NC 714 line.

Summary: Crop plant partial resistance to plant pathogens controlled by quantitative trait loci (QTL) is desirable in cultivar development programmes because of its increased durability. Mechanisms underlying such resistance are difficult to study. We performed RNA‐seq analyses for tobacco (Nicotiana tabacum) nearly isogenic lines (NILs) with and without favourable allele(s) at Phn7.1, a major QTL influencing partial resistance to the soil‐borne pathogens Phytophthora nicotianae and Ralstonia solanacearum. Based upon combined analyses of transcriptome‐based sequence variation and gene expression profiles, we concluded that allelic variability at the Phn7.1 locus was likely generated from homoeologous exchange, which led to deletion of low‐expressing members of the SAR8.2 gene family and duplication of high‐expressing SAR8.2 genes from a different subgenome of allotetraploid tobacco. The high expression of endogenous Phn7.1‐associated SAR8.2 genes was correlated with observed resistance to P. nicotianae. Our findings suggest a role for genomic rearrangements in the generation of favourable genetic variability affecting resistance to pathogens in plants.