Nuclear lamina component KAKU4 regulates chromatin states and transcriptional regulation in the Arabidopsis genome.

BACKGROUND: The nuclear lamina links the nuclear membrane to chromosomes and plays a crucial role in regulating chromatin states and gene expression. However, current knowledge of nuclear lamina in plants is limited compared to animals and humans. RESULTS: This study mainly focused on elucidating the mechanism through which the putative nuclear lamina component protein KAKU4 regulates chromatin states and gene expression in Arabidopsis leaves. Thus, we constructed a network using the association proteins of lamin-like proteins, revealing that KAKU4 is strongly associated with chromatin or epigenetic modifiers. Then, we conducted ChIP-seq technology to generate global epigenomic profiles of H3K4me3, H3K27me3, and H3K9me2 in Arabidopsis leaves for mutant (kaku4-2) and wild-type (WT) plants alongside RNA-seq method to generate gene expression profiles. The comprehensive chromatin state-based analyses indicate that the knockdown of KAKU4 has the strongest effect on H3K27me3, followed by H3K9me2, and the least impact on H3K4me3, leading to significant changes in chromatin states in the Arabidopsis genome. We discovered that the knockdown of the KAKU4 gene caused a transition between two types of repressive epigenetics marks, H3K9me2 and H3K27me3, in some specific PLAD regions. The combination analyses of epigenomic and transcriptomic data between the kaku4-2 mutant and WT suggested that KAKU4 may regulate key biological processes, such as programmed cell death and hormone signaling pathways, by affecting H3K27me3 modification in Arabidopsis leaves. CONCLUSIONS: In summary, our results indicated that KAKU4 is directly and/or indirectly associated with chromatin/epigenetic modifiers and demonstrated the essential roles of KAKU4 in regulating chromatin states, transcriptional regulation, and diverse biological processes in Arabidopsis.

The chromosome-level genome and functional database accelerate research about biosynthesis of secondary metabolites in Rosa roxburghii.

Rosa roxburghii Tratt, a valuable plant in China with long history, is famous for its fruit. It possesses various secondary metabolites, such as L-ascorbic acid (vitamin C), alkaloids and poly saccharides, which make it a high nutritional and medicinal value. Here we characterized the chromosome-level genome sequence of R. roxburghii, comprising seven pseudo-chromosomes with a total size of 531 Mb and a heterozygosity of 0.25%. We also annotated 45,226 coding gene loci after masking repeat elements. Orthologs for 90.1% of the Complete Single-Copy BUSCOs were found in the R. roxburghii annotation. By aligning with protein sequences from public platform, we annotated 85.89% genes from R. roxburghii. Comparative genomic analysis revealed that R. roxburghii diverged from Rosa chinensis approximately 5.58 to 13.17 million years ago, and no whole-genome duplication event occurred after the divergence from eudicots. To fully utilize this genomic resource, we constructed a genomic database RroFGD with various analysis tools. Otherwise, 69 enzyme genes involved in L-ascorbate biosynthesis were identified and a key enzyme in the biosynthesis of vitamin C, GDH (L-Gal-1-dehydrogenase), is used as an example to introduce the functions of the database. This genome and database will facilitate the future investigations into gene function and molecular breeding in R. roxburghii.

KAKU4 regulates leaf senescence through modulation of H3K27me3 deposition in the Arabidopsis genome.

Lamins are the major components of the nuclear lamina, which regulate chromatin structure and gene expression. KAKU4 is a unique nuclear lamina component in the nuclear periphery, modulates nuclear shape and size in Arabidopsis. The knowledge about the regulatory role of KAKU4 in leaf development remains limited. Here we found that knockdown of KAKU4 resulted in an accelerated leaf senescence phenotype, with elevated levels of H2O2 and hormones, particularly SA, JA, and ABA. Our results demonstrated the importance of KAKU4 as a potential negative regulator in age-triggered leaf senescence in Arabidopsis. Furthermore, we conducted combination analyses of transcriptomic and epigenomic data for the kaku4 mutant and WT leaves. The knockdown of KAKU4 lowered H3K27me3 deposition in the up-regulated genes associated with hormone pathways, programmed cell death, and leaf senescence, including SARD1, SAG113/HAI1, PR2, and so forth. In addition, we found the functional crosstalks between KAKU4 and its associated proteins (CRWN1/4, PNET2, GBPL3, etc.) through comparing multiple transcriptome datasets. Overall, our results indicated that KAKU4 may inhibit the expression of a series of genes related to hormone signals and H2O2 metabolism by affecting the deposition of H3K27me3, thereby suppressing leaf senescence.

PlantGSAD: a comprehensive gene set annotation database for plant species.

With the accumulation of massive data sets from high-throughput experiments and the rapid emergence of new types of omics data, gene sets have become more diverse and essential for the refinement of gene annotation at multidimensional levels. Accordingly, we collected and defined 236 007 gene sets across different categories for 44 plant species in the Plant Gene Set Annotation Database (PlantGSAD). These gene sets were divided into nine main categories covering many functional subcategories, such as trait ontology, co-expression modules, chromatin states, and liquid-liquid phase separation. The annotations from the collected gene sets covered all of the genes in the Brassicaceae species Arabidopsis and Poaceae species Oryza sativa. Several GSEA tools are implemented in PlantGSAD to improve the efficiency of the analysis, including custom SEA for a flexible strategy based on customized annotations, SEACOMPARE for the cross-comparison of SEA results, and integrated visualization features for ontological analysis that intuitively reflects their parent-child relationships. In summary, PlantGSAD provides numerous gene sets for multiple plant species and highly efficient analysis tools. We believe that PlantGSAD will become a multifunctional analysis platform that can be used to predict and elucidate the functions and mechanisms of genes of interest. PlantGSAD is publicly available at http://systemsbiology.cau.edu.cn/PlantGSEAv2/.

Cardiac involvement in muscular dystrophies: Role of myocardial perfusion imaging.

Technetium-99m-methoxy isobutyl isonitrile (99mTc-MIBI) myocardial perfusion imaging (MPI) is a functional imaging method with relatively poor specificity but high sensitivity. We present 48-year-old man with cardiac involvement due to muscular dystrophies (MD). Myocardial perfusion imaging rest images revealed regional myocardial perfusion decrease in multiple walls, enlarged heart and decreased left ventricular systolic function. The lesion location of MPI was consistent with that seen on CMR. Our case showed MPI was useful for detection and evaluation of the MD patient with cardiac involvement. In addition, imaging findings in combination with clinical history and other data are important. The case highlight is thevalue of MPI in myocardiopathy.

Green leaf volatile (Z)-3-hexeny-1-yl acetate reduces salt stress in peanut by affecting photosynthesis and cellular redox homeostasis.

Green leaf volatiles (GLVs) are released by plants when they encounter biotic stress, but their functions in the response to abiotic stress have not been determined. We have previously shown that exogenous application of (Z)-3-hexeny-1-yl acetate (Z-3-HAC), a kind of GLV, could alleviate salt stress in peanut (Arachis hypogaea L.) seedlings; however, notably little is known concerning the transcription regulation mechanisms of Z-3-HAC. In this study, we comprehensively characterized the transcriptomes and physiological indices of peanut seedlings exposed to Z-3-HAC and/or salt stress. Analysis of transcriptome data showed that 1420 genes were upregulated in the seedlings primed with Z-3-HAC under salt stress compared with the non-primed treatment. Interestingly, these genes were significantly enriched in the photosynthetic and ascorbate metabolism-related categories, as well as several plant hormone metabolism pathways. The physiological data revealed that Z-3-HAC significantly increased the net photosynthetic rate, SPAD value, plant height and shoot biomass compared with the non-primed peanut seedlings under salt stress. A significantly higher ratio of K+ :Na+ , reduced-to-oxidized glutathione (GSH:GSSG), and ascorbate-to-dehydroascorbate (AsA:DHA) were also observed for the plants primed with Z-3-HAC compared with the salt stress control. Meanwhile, Z-3-HAC significantly increased the activity of enzymes in the AsA-GSH cycle. Taken together, these results highlight the importance of Z-3-HAC in protecting peanut seedlings against salt stress by affecting photosynthesis, cellular redox homeostasis, K+ :Na+ homeostasis, and phytohormones.

PCSD: a plant chromatin state database.

Genome-wide maps of chromatin states have become a powerful representation of genome annotation and regulatory activity. We collected public and in-house plant epigenomic data sets and applied a Hidden Markov Model to define chromatin states, which included 290 553 (36 chromatin states), 831 235 (38 chromatin states) and 3 936 844 (26 chromatin states) segments across the whole genome of Arabidopsis thaliana, Oryza sativa and Zea mays, respectively. We constructed a Plant Chromatin State Database (PCSD, http://systemsbiology.cau.edu.cn/chromstates) to integrate detailed information about chromatin states, including the features and distribution of states, segments in states and related genes with segments. The self-organization mapping (SOM) results for these different chromatin signatures and UCSC Genome Browser for visualization were also integrated into the PCSD database. We further provided differential SOM maps between two epigenetic marks for chromatin state comparison and custom tools for new data analysis. The segments and related genes in SOM maps can be searched and used for motif and GO analysis, respectively. In addition, multi-species integration can be used to discover conserved features at the epigenomic level. In summary, our PCSD database integrated the identified chromatin states with epigenetic features and may be beneficial for communities to discover causal functions hidden in plant chromatin.

agriGO v2.0: a GO analysis toolkit for the agricultural community, 2017 update.

The agriGO platform, which has been serving the scientific community for >10 years, specifically focuses on gene ontology (GO) enrichment analyses of plant and agricultural species. We continuously maintain and update the databases and accommodate the various requests of our global users. Here, we present our updated agriGO that has a largely expanded number of supporting species (394) and datatypes (865). In addition, a larger number of species have been classified into groups covering crops, vegetables, fish, birds and insects closely related to the agricultural community. We further improved the computational efficiency, including the batch analysis and P-value distribution (PVD), and the user-friendliness of the web pages. More visualization features were added to the platform, including SEACOMPARE (cross comparison of singular enrichment analysis), direct acyclic graph (DAG) and Scatter Plots, which can be merged by choosing any significant GO term. The updated platform agriGO v2.0 is now publicly accessible at http://systemsbiology.cau.edu.cn/agriGOv2/.

A systemic identification approach for primary transcription start site of Arabidopsis miRNAs from multidimensional omics data.

The 22-nucleotide non-coding microRNAs (miRNAs) are mostly transcribed by RNA polymerase II and are similar to protein-coding genes. Unlike the clear process from stem-loop precursors to mature miRNAs, the primary transcriptional regulation of miRNA, especially in plants, still needs to be further clarified, including the original transcription start site, functional cis-elements and primary transcript structures. Due to several well-characterized transcription signals in the promoter region, we proposed a systemic approach integrating multidimensional "omics" (including genomics, transcriptomics, and epigenomics) data to improve the genome-wide identification of primary miRNA transcripts. Here, we used the model plant Arabidopsis thaliana to improve the ability to identify candidate promoter locations in intergenic miRNAs and to determine rules for identifying primary transcription start sites of miRNAs by integrating high-throughput omics data, such as the DNase I hypersensitive sites, chromatin immunoprecipitation-sequencing of polymerase II and H3K4me3, as well as high throughput transcriptomic data. As a result, 93% of refined primary transcripts could be confirmed by the primer pairs from a previous study. Cis-element and secondary structure analyses also supported the feasibility of our results. This work will contribute to the primary transcriptional regulatory analysis of miRNAs, and the conserved regulatory pattern may be a suitable miRNA characteristic in other plant species.

PlantCHRs: A comprehensive database of plant chromatin remodeling factors.

The Snf2 protein family is a group of ATP-dependent chromatin remodeling factors (CHRs) that play an essential role in gene expression regulation. In plants, Snf2 is involved in growth, development, as well as stress resistance. However, only a very limited number of experimentally validated Snf2 have been identified and reported, while the majority remaining undiscovered in most species . In this study, we predicted 3135 Snf2 proteins and 8398 chromatin remodeling complex (CRC) subunits in diverse plant species, and constructed the Plant Chromatin Remodeling Factors Database (PlantCHRs, http://www.functionalgenomics.cn/PlantCHRs/), which provide a comprehensive resource for researchers to access information about plant CHRs. We also developed an online tool capable of predicting CHRs and CRC subunits. Moreover, we investigated the distribution of Snf2 proteins in different species and observed a significant increase in the number of Snf2 proteins and the diversity of the Snf2 subfamily during the evolution, highlighting their evolutionary importance. By analyzing the expression patterns of the Snf2 genes in different tissues of maize and Arabidopsis, we found that the Snf2 proteins may show some conservation across different species in regulating plant growth and development. Over the all, we established a comprehensive database for plant CHRs, which will facilitate the researches on plant chromatin remodeling.

GURFAP: A Platform for Gene Function Analysis in Glycyrrhiza Uralensis.

Glycyrrhiza uralensis (Licorice), which belongs to Leguminosae, is famous for the function of pharmacologic action and natural sweetener with its dried roots and rhizomes. In recent years, the whole-genome sequence of G. uralensis has been completed, which will help to lay the foundation for the study of gene function. Here, we integrated the available genomic and transcriptomic data of G. uralensis and constructed the G. uralensis gene co-expression network. We then annotated gene functions of G. uralensis via aligning with public databases. Furthermore, gene families of G. uralensis were predicted by tools including iTAK (Plant Transcription factor and Protein kinase Identifier and Classifier), HMMER (hidden Markov models), InParanoid, and PfamScan. Finally, we constructed a platform for gene function analysis in G. uralensis (GURFAP, www.gzybioinfoormatics.cn/GURFAP). For analyzed and predicted gene function, we introduced various tools including BLAST (Basic local alignment search tool), GSEA (Gene set enrichment analysis), Motif, Heatmap, and JBrowse. Our analysis based on this platform indicated that the biosynthesis of glycyrrhizin might be regulated by MYB and bHLH. We also took CYP88D6, CYP72A154, and bAS gene in the synthesis pathway of glycyrrhizin as examples to demonstrate the reliability and availability of our platform. Our platform GURFAP will provide convenience for researchers to mine the gene function of G. uralensis and thus discover more key genes involved in the biosynthetic pathway of active ingredients.

Refinement of bamboo genome annotations through integrative analyses of transcriptomic and epigenomic data.

Bamboo, one of the most crucial nontimber forest resources worldwide, has the capacity for rapid growth. In recent years, the genome of moso bamboo (Phyllostachys edulis) has been decoded, and a large amount of transcriptome data has been published. In this study, we generated the genome-wide profiles of the histone modification H3K4me3 in leaf, stem, and root tissues of bamboo. The trends in the distribution patterns were similar to those in rice. We developed a processing pipeline for predicting novel transcripts to refine the structural annotation of the genome using H3K4me3 ChIP-seq data and 29 RNA-seq datasets. As a result, 12,460 novel transcripts were predicted in the bamboo genome. Compared with the transcripts in the newly released version 2.0 of the bamboo genome, these novel transcripts are tissue-specific and shorter, and most have a single exon. Some representative novel transcripts were validated by semiquantitative RT-PCR and qRT-PCR analyses. Furthermore, we put these novel transcripts back into the ChIP-seq analysis pipeline and discovered that the percentages of H3K4me3 in genic elements were increased. Overall, this work integrated transcriptomic data and epigenomic data to refine the annotation of the genome in order to discover more functional genes and study bamboo growth and development, and the application of this predicted pipeline may help refine the structural annotation of the genome in other species.

Systematic Analysis of Differential H3K27me3 and H3K4me3 Deposition in Callus and Seedling Reveals the Epigenetic Regulatory Mechanisms Involved in Callus Formation in Rice.

Plant growth and development occurs through meristematic cell activity, and cell fate transition is accompanied by epigenetic modifications. Callus with cell pluripotency exhibits the ability to undergo continued cell division, and is ideal for studying plant meristematic differentiation. By comparing the differential epigenetic modifications between callus and seedling, the changes in chromatin state and effects of various epigenetic modifications on the growth and development of plants can be revealed, and the key genes related to plant growth and development can be identified, providing novel insights into the regulation of plant growth and development. In this study, we performed ChIP assays using various antibodies in rice seed-induced callus and seedlings grown for about 15 days to examine the differential deposition of H3K27me3 and H3K4me3. Furthermore, data for DNase I-hypersensitive sites in the corresponding tissues were downloaded from National Center for Biotechnology Information. We analyzed 4,562 callus H3K27me3-decreased genes especially those encoding transcription factors in callus, and found that most of the transcription factors, including AP2-ERREBP, NAC, and HB gene families, were related to growth and development. Genes related to meristemization, such as OsWOX9, OsWOX11, OsPLT4, OsPLT5, and OsSHR, were also included. In contrast, H3K4me3 positively regulated callus characteristics through its higher deposition in the callus than in the seedling. We further performed transcriptomic analysis on 45 sets of Affymetrix GeneChip arrays and identified 1,565 genes preferentially expressed in the callus. Callus development and root development in rice were found to share a common regulatory mechanism. We found that these genes, which are associated with meristems, require the removal of H3K27me3 and the deposition of H3K4me3, and DNase I-hypersensitive sites to maintain a relatively active state in the callus than in the seedling. The present study provides novel data about the epigenetic mechanisms involved in callus formation and additional resources for the study of cell division and differentiation in plants.

croFGD: Catharanthus roseus Functional Genomics Database.

Catharanthus roseus is a medicinal plant, which can produce monoterpene indole alkaloid (MIA) metabolites with biological activity and is rich in vinblastine and vincristine. With release of the scaffolded genome sequence of C. roseus, it is necessary to annotate gene functions on the whole-genome level. Recently, 53 RNA-seq datasets are available in public with different tissues (flower, root, leaf, seedling, and shoot) and different treatments (MeJA, PnWB infection and yeast elicitor). We used in-house data process pipeline with the combination of PCC and MR algorithms to construct a co-expression network exploring multi-dimensional gene expression (global, tissue preferential, and treat response) through multi-layered approaches. In the meanwhile, we added miRNA-target pairs, predicted PPI pairs into the network and provided several tools such as gene set enrichment analysis, functional module enrichment analysis, and motif analysis for functional prediction of the co-expression genes. Finally, we have constructed an online croFGD database (http://bioinformatics.cau.edu.cn/croFGD/). We hope croFGD can help the communities to study the C. roseus functional genomics and make novel discoveries about key genes involved in some important biological processes.

Chromatin State-Based Analysis of Epigenetic H3K4me3 Marks of Arabidopsis in Response to Dark Stress.

Light is essential to plant growth and development. Extended darkness causes dramatic gene expression changes, leading to leaf senescence, hypocotyl growth, petiole elongation, reduced leaf area, and early flowering, etc. However, the underlying mechanism of response to darkness at epigenetic levels remains largely unknown. In this study, we conducted ChIP-seq to generate global epigenomic profiles of H3K4me3 under 3-day extended darkness and normal light conditions in Arabidopsis. We applied chromatin state analysis together with self-organization mapping (SOM) to study the combination of epigenetic regulation under dark stress. The SOM map clusters the segments on the genome according to multiple diverse epigenomic datasets, which breaks the limit of dispersed distribution of epigenetic marks on the genome. Through SOM analysis, we also found that the signals of H3K4me3 were mainly increased after darkness. Analysis of H3K4me3-changed genes together with differentially expressed genes indicated that the genes showing dark-increased H3K4me3 were most involved in senescence and autophagy, and cross-talk existed between dark-induced and natural senescence. In summary, we studied the regulation of the epigenetic H3K4me3 marks of Arabidopsis in response to dark stress using chromatin state and SOM analyses. Our study revealed the regulatory mechanisms of the epigenome in response to dark stress, and SOM analysis based on chromatin states used in our study will also be helpful for other studies on dynamic changes of multiple epigenetic marks.

AtSPX1-mediated transcriptional regulation during leaf senescence in Arabidopsis thaliana.

Leaf senescence is the final stage of leaf growth, a highly coordinated and complicated process. Phosphorus as an essential macronutrient for plant growth is remobilized from senescing leaves to other vigorous parts of the plant. In this study, through data mining, we found some phosphate starvation induced genes such as AtSPX1, were significantly induced in aging leaves in Arabidopsis. We applied a reverse genetics approach to investigate the phenotypes of transgenic plants and mutant plants, and the results showed that the overexpression of AtSPX1 accelerated leaf senescence, suppressed Pi accumulation, promoted SA production and H2O2 levels in leaves, while the mutant lines of AtSPX1 showed slightly delayed leaf senescence. We conducted RNA-seq-based transcriptome analysis together with GO and GSEA enrichment analyses for transgenic vs. wild-type plants to elucidate the possible underlying regulatory mechanism. The 558 genes that were up-regulated in the overexpression plants 35S::AtSPX1/WT, were significantly enriched in the process of leaf senescence, Pi starvation responses and SA signaling pathways, as were the target genes of some transcription factors such as WRKYs and NACs. In a word, we characterized AtSPX1 as a key regulator, which mediated the crosstalks among leaf senescence, Pi starvation and SA signaling pathways in Arabidopsis thaliana.

MCENet: A database for maize conditional co-expression network and network characterization collaborated with multi-dimensional omics levels.

Maize (Zea mays) is the most widely grown grain crop in the world, playing important roles in agriculture and industry. However, the functions of maize genes remain largely unknown. High-quality genome-wide transcriptome datasets provide important biological knowledge which has been widely and successfully used in plants not only by measuring gene expression levels but also by enabling co-expression analysis for predicting gene functions and modules related to agronomic traits. Recently, thousands of maize transcriptomic data are available across different inbred lines, development stages, tissues, and treatments, or even across different tissue sections and cell lines. Here, we integrated 701 transcriptomic and 108 epigenomic data and studied the different conditional networks with multi-dimensional omics levels. We constructed a searchable, integrative, one-stop online platform, the maize conditional co-expression network (MCENet) platform. MCENet provides 10 global/conditional co-expression networks, 5 network accessional analysis toolkits (i.e., Network Search, Network Remodel, Module Finder, Network Comparison, and Dynamic Expression View) and multiple network functional support toolkits (e.g., motif and module enrichment analysis). We hope that our database might help plant research communities to identify maize functional genes or modules that regulate important agronomic traits. MCENet is publicly accessible at http://bioinformatics.cau.edu.cn/MCENet/.

PlantEAR: Functional Analysis Platform for Plant EAR Motif-Containing Proteins.

The Ethylene-responsive element binding factor-associated Amphiphilic Repression (EAR) motifs, which were initially identified in members of the Arabidopsis ethylene response factor (ERF) family, are transcriptional repression motifs in plants and are defined by the consensus sequence patterns of either LxLxL or DLNxxP. EAR motif-containing proteins can function as transcription repressors, thus interacting with co-repressors, such as TOPLESS and AtSAP18, affecting the structure of chromatin by histone modifications and thereby repressing gene transcription. EAR motif-containing proteins are highly conserved across diverse plant species and play important roles in hormone signal transduction, stress responses and development, but they have not been identified in most plants. In this study, we identified 20,542 EAR motif-containing proteins from 71 plant species based on a Hidden Markov Model and orthologous gene search, and then we constructed a functional analysis platform for plant EAR motif-containing proteins (PlantEAR, http://structuralbiology.cau.edu.cn/plantEAR) by integrating a variety of functional annotations and processed data. Several tools were provided as functional support for EAR motif-containing proteins, such as browse, search, co-expression and protein-protein interaction (PPI) network analysis as well as cis-element analysis and gene set enrichment analysis (GSEA). In addition, basing on the identified EAR motif-containing proteins, we also explored their distribution in various species and found that the numbers of EAR motif-containing proteins showed an increasing trend in evolution from algae to angiosperms.

Co-expression Gene Network Analysis and Functional Module Identification in Bamboo Growth and Development.

Bamboo is one of the fastest-growing non-timber forest plants. Moso bamboo (Phyllostachys edulis) is the most economically valuable bamboo in Asia, especially in China. With the release of the whole-genome sequence of moso bamboo, there are increasing demands for refined annotation of bamboo genes. Recently, large amounts of bamboo transcriptome data have become available, including data on the multiple growth stages of tissues. It is now feasible for us to construct co-expression networks to improve bamboo gene annotation and reveal the relationships between gene expression and growth traits. We integrated the genome sequence of moso bamboo and 78 transcriptome data sets to build genome-wide global and conditional co-expression networks. We overlaid the gene expression results onto the network with multiple dimensions (different development stages). Through combining the co-expression network, module classification and function enrichment tools, we identified 1,896 functional modules related to bamboo development, which covered functions such as photosynthesis, hormone biosynthesis, signal transduction, and secondary cell wall biosynthesis. Furthermore, an online database (http://bioinformatics.cau.edu.cn/bamboo) was built for searching the moso bamboo co-expression network and module enrichment analysis. Our database also includes cis-element analysis, gene set enrichment analysis, and other tools. In summary, we integrated public and in-house bamboo transcriptome data sets and carried out co-expression network analysis and functional module identification. Through data mining, we have yielded some novel insights into the regulation of growth and development. Our established online database might be convenient for the bamboo research community to identify functional genes or modules with important traits.

Genome-wide mapping of DNase I hypersensitive sites reveals chromatin accessibility changes in Arabidopsis euchromatin and heterochromatin regions under extended darkness.

Light, as the energy source in photosynthesis, is essential for plant growth and development. Extended darkness causes dramatic gene expression changes. In this study, we applied DNase-seq (DNase I hypersensitive site sequencing) to study changes of chromatin accessibility in euchromatic and heterochromatic regions under extended darkness in Arabidopsis. We generated 27 Gb DNase-seq and 67.6 Gb RNA-seq data to investigate chromatin accessibility changes and global gene expression under extended darkness and control condition in Arabidopsis. We found that ~40% DHSs (DNaseI hypersensitive sites) were diminished under darkness. In non-TE regions, the majority of DHS-changed genes were DHS-diminished under darkness. A total of 519 down-regulated genes were associated with diminished DHSs under darkness, mainly involved in photosynthesis process and retrograde signaling, and were regulated by chloroplast maintenance master regulators such as GLK1. In TE regions, approximately half of the DHS-changed TEs were DHS-increased under darkness and were primarily associated with the LTR/Gypsy retrotransposons in the heterochromatin flanking the centromeres. In contrast, DHS-diminished TEs under darkness were enriched in Copia, LINE, and MuDR dispersed across chromosomes. Together, our results indicated that extended darkness resulted in more increased chromatin compaction in euchromatin and decompaction in heterochromatin, thus further leading to gene expression changes in Arabidopsis.