Evaluation of experimental design and computational parameter choices affecting analyses of ChIP-seq and RNA-seq data in undomesticated poplar trees

Genome-wide profiling of histone (H3) lysine 4 (K4) tri-methylation (me3) under drought, heat, and combined stresses in switchgrass

BACKGROUND

Switchgrass (Panicum virgatum L.) is a warm-season perennial (C4) grass identified as an important biofuel crop in the United States. It is well adapted to the marginal environment where heat and moisture stresses predominantly affect crop growth. However, the underlying molecular mechanisms associated with heat and drought stress tolerance still need to be fully understood in switchgrass. The methylation of H3K4 is often associated with transcriptional activation of genes, including stress-responsive. Therefore, this study aimed to analyze genome-wide histone H3K4-tri-methylation in switchgrass under heat, drought, and combined stress.

RESULTS

In total, ~ 1.3 million H3K4me3 peaks were identified in this study using SICER. Among them, 7,342; 6,510; and 8,536 peaks responded under drought (DT), drought and heat (DTHT), and heat (HT) stresses, respectively. Most DT and DTHT peaks spanned 0 to + 2000 bases from the transcription start site [TSS]. By comparing differentially marked peaks with RNA-Seq data, we identified peaks associated with genes: 155 DT-responsive peaks with 118 DT-responsive genes, 121 DTHT-responsive peaks with 110 DTHT-responsive genes, and 175 HT-responsive peaks with 136 HT-responsive genes. We have identified various transcription factors involved in DT, DTHT, and HT stresses. Gene Ontology analysis using the AgriGO revealed that most genes belonged to biological processes. Most annotated peaks belonged to metabolite interconversion, RNA metabolism, transporter, protein modifying, defense/immunity, membrane traffic protein, transmembrane signal receptor, and transcriptional regulator protein families. Further, we identified significant peaks associated with TFs, hormones, signaling, fatty acid and carbohydrate metabolism, and secondary metabolites. qRT-PCR analysis revealed the relative expressions of six abiotic stress-responsive genes (transketolase, chromatin remodeling factor-CDH3, fatty-acid desaturase A, transmembrane protein 14C, beta-amylase 1, and integrase-type DNA binding protein genes) that were significantly (P < 0.05) marked during drought, heat, and combined stresses by comparing stress-induced against un-stressed and input controls.

CONCLUSION

Our study provides a comprehensive and reproducible epigenomic analysis of drought, heat, and combined stress responses in switchgrass. Significant enrichment of H3K4me3 peaks downstream of the TSS of protein-coding genes was observed. In addition, the cost-effective experimental design, modified ChIP-Seq approach, and analyses presented here can serve as a prototype for other non-model plant species for conducting stress studies.

PSDX: A Comprehensive Multi-Omics Association Database of Populus trichocarpa With a Focus on the Secondary Growth in Response to Stresses

Populus trichocarpa (P. trichocarpa) is a model tree for the investigation of wood formation. In recent years, researchers have generated a large number of high-throughput sequencing data in P. trichocarpa. However, no comprehensive database that provides multi-omics associations for the investigation of secondary growth in response to diverse stresses has been reported. Therefore, we developed a public repository that presents comprehensive measurements of gene expression and post-transcriptional regulation by integrating 144 RNA-Seq, 33 ChIP-seq, and six single-molecule real-time (SMRT) isoform sequencing (Iso-seq) libraries prepared from tissues subjected to different stresses. All the samples from different studies were analyzed to obtain gene expression, co-expression network, and differentially expressed genes (DEG) using unified parameters, which allowed comparison of results from different studies and treatments. In addition to gene expression, we also identified and deposited pre-processed data about alternative splicing (AS), alternative polyadenylation (APA) and alternative transcription initiation (ATI). The post-transcriptional regulation, differential expression, and co-expression network datasets were integrated into a new P. trichocarpa Stem Differentiating Xylem (PSDX) database (http://forestry.fafu.edu.cn/db/SDX), which further highlights gene families of RNA-binding proteins and stress-related genes. The PSDX also provides tools for data query, visualization, a genome browser, and the BLAST option for sequence-based query. Much of the data is also available for bulk download. The availability of PSDX contributes to the research related to the secondary growth in response to stresses in P. trichocarpa, which will provide new insights that can be useful for the improvement of stress tolerance in woody plants.

Comparative transcriptome analyses define genes and gene modules differing between two Populus genotypes with contrasting stem growth rates

BACKGROUND

Wood provides an important biomass resource for biofuel production around the world. The radial growth of tree stems is central to biomass production for forestry and biofuels, but it is challenging to dissect genetically because it is a complex trait influenced by many genes. In this study, we adopted methods of physiology, transcriptomics and genetics to investigate the regulatory mechanisms of tree radial growth and wood development.

RESULTS

Physiological comparison showed that two Populus genotypes presented different rates of radial growth of stems and accumulation of woody biomass. A comparative transcriptional network approach was used to define and characterize functional differences between two Populus genotypes. Analyses of transcript profiles from wood-forming tissue of the two genotypes showed that 1542, 2295 and 2110 genes were differentially expressed in the pre-growth, fast-growth and post-growth stages, respectively. The co-expression analyses identified modules of co-expressed genes that displayed distinct expression profiles. Modules were further characterized by correlating transcript levels with genotypes and physiological traits. The results showed enrichment of genes that participated in cell cycle and division, whose expression change was consistent with the variation of radial growth rates. Genes related to secondary vascular development were up-regulated in the faster-growing genotype in the pre-growth stage. We characterized a BEL1-like (BELL) transcription factor, PeuBELL15, which was up-regulated in the faster-growing genotype. Analyses of transgenic Populus overexpressing as well as CRISPR/Cas9-induced mutants for BELL15 showed that PeuBELL15 improved accumulation of glucan and lignin, and it promoted secondary vascular growth by regulating the expression of genes relevant for cellulose synthases and lignin biosynthesis.

CONCLUSIONS

This study illustrated that active division and expansion of vascular cambium cells and secondary cell wall deposition of xylem cells contribute to stem radial increment and biomass accumulation, and it identified relevant genes for these complex growth traits, including a BELL transcription factor gene PeuBELL15. This provides genetic resources for improving and breeding elite genotypes with fast growth and high wood biomass.

The developmental dynamics of the Populus stem transcriptome

The Populus shoot undergoes primary growth (longitudinal growth) followed by secondary growth (radial growth), which produces biomass that is an important source of energy worldwide. We adopted joint PacBio Iso-Seq and RNA-seq analysis to identify differentially expressed transcripts along a developmental gradient from the shoot apex to the fifth internode of Populus Nanlin895. We obtained 87 150 full-length transcripts, including 2081 new isoforms and 62 058 new alternatively spliced isoforms, most of which were produced by intron retention, that were used to update the Populus annotation. Among these novel isoforms, there are 1187 long non-coding RNAs and 356 fusion genes. Using this annotation, we found 15 838 differentially expressed transcripts along the shoot developmental gradient, of which 1216 were transcription factors (TFs). Only a few of these genes were reported previously. The differential expression of these TFs suggests that they may play important roles in primary and secondary growth. AP2, ARF, YABBY and GRF TFs are highly expressed in the apex, whereas NAC, bZIP, PLATZ and HSF TFs are likely to be important for secondary growth. Overall, our findings provide evidence that long-read sequencing can complement short-read sequencing for cataloguing and quantifying eukaryotic transcripts and increase our understanding of the vital and dynamic process of shoot development.

Transcriptional and temporal response of Populus stems to gravi-stimulation

Plants modify development in response to external stimuli, to produce new growth that is appropriate for environmental conditions. For example, gravi-stimulation of leaning branches in angiosperm trees results in modifications of wood development, to produce tension wood that pulls leaning stems upright. Here, we use gravi-stimulation and tension wood response to dissect the temporal changes in gene expression underlying wood formation in Populus stems. Using time-series analysis of seven time points over a 14-d experiment, we identified 8,919 genes that were differentially expressed between tension wood (upper) and opposite wood (lower) sides of leaning stems. Clustering of differentially expressed genes showed four major transcriptional responses, including gene clusters whose transcript levels were associated with two types of tissue-specific impulse responses that peaked at about 24-48 h, and gene clusters with sustained changes in transcript levels that persisted until the end of the 14-d experiment. Functional enrichment analysis of those clusters suggests they reflect temporal changes in pathways associated with hormone regulation, protein localization, cell wall biosynthesis and epigenetic processes. Time-series analysis of gene expression is an underutilized approach for dissecting complex developmental responses in plants, and can reveal gene clusters and mechanisms influencing development.

Identifying gene coexpression networks underlying the dynamic regulation of wood-forming tissues in Populus under diverse environmental conditions

Trees modify wood formation through integration of environmental and developmental signals in complex but poorly defined transcriptional networks, allowing trees to produce woody tissues appropriate to diverse environmental conditions. In order to identify relationships among genes expressed during wood formation, we integrated data from new and publically available datasets in Populus. These datasets were generated from woody tissue and include transcriptome profiling, transcription factor binding, DNA accessibility and genome-wide association mapping experiments. Coexpression modules were calculated, each of which contains genes showing similar expression patterns across experimental conditions, genotypes and treatments. Conserved gene coexpression modules (four modules totaling 8398 genes) were identified that were highly preserved across diverse environmental conditions and genetic backgrounds. Functional annotations as well as correlations with specific experimental treatments associated individual conserved modules with distinct biological processes underlying wood formation, such as cell-wall biosynthesis, meristem development and epigenetic pathways. Module genes were also enriched for DNase I hypersensitivity footprints and binding from four transcription factors associated with wood formation. The conserved modules are excellent candidates for modeling core developmental pathways common to wood formation in diverse environments and genotypes, and serve as testbeds for hypothesis generation and testing for future studies.

DNA-Binding Factor Target Identification by Chromatin Immunoprecipitation (ChIP) in Plants

Chromatin immunoprecipitation (ChIP) allows the precise identification of genomic loci that physically interact with a protein of interest, whether that protein is a transcription factor, a core polymerase, a histone, or other chromatin-associated protein. In short, tissue is first cross-linked to freeze a population of DNA-protein interactions at a stage of interest. Chromatin is then extracted, fragmented, and incubated with a specific antibody against the protein of interest. Next, the resultant DNA-protein complexes are immunoprecipitated and captured using beads that bind to the antibody constant region. Samples are finally reverse cross-linked to separate the bound fragments and the DNA is purified. This DNA is analyzed by quantitative PCR for enrichment of genomic regions expected to be bound by the protein under study. The protocol detailed in this chapter has been successfully applied in the identification of target genes for seven transcriptional regulators of diverse classes involved in Arabidopsis thaliana floral transition.

A resource for characterizing genome-wide binding and putative target genes of transcription factors expressed during secondary growth and wood formation in Populus

Identifying transcription factor target genes is essential for modeling the transcriptional networks underlying developmental processes. Here we report a chromatin immunoprecipitation sequencing (ChIP-seq) resource consisting of genome-wide binding regions and associated putative target genes for four Populus homeodomain transcription factors expressed during secondary growth and wood formation. Software code (programs and scripts) for processing the Populus ChIP-seq data are provided within a publically available iPlant image, including tools for ChIP-seq data quality control and evaluation adapted from the human Encyclopedia of DNA Elements (ENCODE) project. Basic information for each transcription factor (including members of Class I KNOX, Class III HD ZIP, BEL1-like families) binding are summarized, including the number and location of binding regions, distribution of binding regions relative to gene features, associated putative target genes, and enriched functional categories of putative target genes. These ChIP-seq data have been integrated within the Populus Genome Integrative Explorer (PopGenIE) where they can be analyzed using a variety of web-based tools. We present an example analysis that shows preferential binding of transcription factor ARBORKNOX1 to the nearest neighbor genes in a pre-calculated co-expression network module, and enrichment for meristem-related genes within this module including multiple orthologs of Arabidopsis KNOTTED-like Arabidopsis 2/6.

Differential gene expression and alternative splicing between diploid and tetraploid watermelon

The exploitation of synthetic polyploids for producing seedless fruits is well known in watermelon. Tetraploid progenitors of triploid watermelon plants, compared with their diploid counterparts, exhibit wide phenotypic differences. Although many factors modulate alternative splicing (AS) in plants, the effects of autopolyploidization on AS are still unknown. In this study, we used tissues of leaf, stem, and fruit of diploid and tetraploid sweet watermelon to understand changes in gene expression and the occurrence of AS. RNA-sequencing analysis was performed along with reverse transcription quantitative PCR and rapid amplification of cDNA ends (RACE)-PCR to demonstrate changes in expression and splicing. All vegetative tissues except fruit showed an increased level of AS in the tetraploid watermelon throughout the growth period. The ploidy levels of diploids and the tetraploid were confirmed using a ploidy analyser. We identified 5362 and 1288 genes that were up- and downregulated, respectively, in tetraploid as compared with diploid plants. We further confirmed that 22 genes underwent AS events across tissues, indicating possibilities of generating different protein isoforms with altered functions of important transcription factors and transporters. Arginine biosynthesis, chlorophyllide synthesis, GDP mannose biosynthesis, trehalose biosynthesis, and starch and sucrose degradation pathways were upregulated in autotetraploids. Phloem protein 2, chloroplastic PGR5-like protein, zinc-finger protein, fructokinase-like 2, MYB transcription factor, and nodulin MtN21 showed AS in fruit tissues. These results should help in developing high-quality seedless watermelon and provide additional transcriptomic information related to other cucurbits.

The Populus ARBORKNOX1 homeodomain transcription factor regulates woody growth through binding to evolutionarily conserved target genes of diverse function

The class I KNOX homeodomain transcription factor ARBORKNOX1 (ARK1) is a key regulator of vascular cambium maintenance and cell differentiation in Populus. Currently, basic information is lacking concerning the distribution, functional characteristics, and evolution of ARK1 binding in the Populus genome. Here, we used chromatin immunoprecipitation sequencing (ChIP-seq) technology to identify ARK1 binding loci genome-wide in Populus. Computational analyses evaluated the distribution of ARK1 binding loci, the function of genes associated with bound loci, the effect of ARK1 binding on transcript levels, and evolutionary conservation of ARK1 binding loci. ARK1 binds to thousands of loci which are highly enriched proximal to the transcriptional start sites of genes of diverse functions. ARK1 target genes are significantly enriched in paralogs derived from the whole-genome salicoid duplication event. Both ARK1 and a maize (Zea mays) homolog, KNOTTED1, preferentially target evolutionarily conserved genes. However, only a small portion of ARK1 target genes are significantly differentially expressed in an ARK1 over-expression mutant. This study describes the functional characteristics and evolution of DNA binding by a transcription factor in an undomesticated tree, revealing complexities similar to those shown for transcription factors in model animal species.