Genome-wide analysis of local chromatin packing in Arabidopsis thaliana

Chromosome-level reference genome and resequencing of 322 accessions reveal evolution, genomic imprint and key agronomic traits in adzuki bean

Adzuki bean (Vigna angularis) is an important legume crop cultivated in over 30 countries worldwide. We developed a high-quality chromosome-level reference genome of adzuki bean cultivar Jingnong6 by combining PacBio Sequel long-read sequencing with short-read and Hi-C technologies. The assembled genome covers 97.8% of the adzuki bean genome with a contig N50 of approximately 16 Mb and a total of 32 738 protein-coding genes. We also generated a comprehensive genome variation map of adzuki bean by whole-genome resequencing (WGRS) of 322 diverse adzuki beans accessions including both wild and cultivated. Furthermore, we have conducted comparative genomics and a genome-wide association study (GWAS) on key agricultural traits to investigate the evolution and domestication. GWAS identified several candidate genes, including VaCycA3;1, VaHB15, VaANR1 and VaBm, that exhibited significant associations with domestication traits. Furthermore, we conducted functional analyses on the roles of VaANR1 and VaBm in regulating seed coat colour. We provided evidence for the highest genetic diversity of wild adzuki (Vigna angularis var. nipponensis) in China with the presence of the most original wild adzuki bean, and the occurrence of domestication process facilitating transition from wild to cultigen. The present study elucidates the genetic basis of adzuki bean domestication traits and provides crucial genomic resources to support future breeding efforts in adzuki bean.

Chromosome-Scale Genome Assembly of the Sheep-Biting Louse Bovicola ovis Using Nanopore Sequencing Data and Pore-C Analysis

Bovicola ovis, commonly known as the sheep-biting louse, is an ectoparasite that adversely affects the sheep industry. Sheep louse infestation lowers the quality of products, including wool and leather, causing a loss of approximately AUD 123M per annum in Australia alone. The lack of a high-quality genome assembly for the sheep-biting louse, as well as any closely related livestock lice, has hindered the development of louse research and management control tools. In this study, we present the assembly of B. ovis with a genome size of ~123 Mbp based on a nanopore long-read sequencing library and Illumina RNA sequencing, complemented with a chromosome-level scaffolding using the Pore-C multiway chromatin contact dataset. Combining multiple alignment and gene prediction tools, a comprehensive annotation on the assembled B. ovis genome was conducted and recalled 11,810 genes as well as other genomic features including orf, ssr, rRNA and tRNA. A manual curation using alignment with the available closely related louse species, Pediculus humanus, increased the number of annotated genes to 16,024. Overall, this study reported critical genetic resources and biological insights for the advancement of sheep louse research and the development of sustainable control strategies in the sheep industry.

The MVPs (masterful versatile players): Chromatin factors as pivotal mediators between 3D genome organization and the response to environment

In recent years, the study of genome dynamics has become a prominent research field due to its influence on understanding the control of gene expression. The study of 3D genome organization has unveiled multiple mechanisms in orchestrating chromosome folding. Growing evidence reveals that these mechanisms are not only important for genome organization, but play a pivotal role in enabling plants to adapt to environmental stimuli. In this review, we provide an overview of the current knowledge concerning epigenetic factors and regulatory elements driving 3D genome dynamics and their responses to external stimuli. We discuss the most recent findings, previous evidence, and explore their implications for future research.

An atlas of the tomato epigenome reveals that KRYPTONITE shapes TAD-like boundaries through the control of H3K9ac distribution

In recent years, the exploration of genome three-dimensional (3D) conformation has yielded profound insights into the regulation of gene expression and cellular functions in both animals and plants. While animals exhibit a characteristic genome topology defined by topologically associating domains (TADs), plants display similar features with a more diverse conformation across species. Employing advanced high-throughput sequencing and microscopy techniques, we investigated the landscape of 26 histone modifications and RNA polymerase II distribution in tomato (Solanum lycopersicum). Our study unveiled a rich and nuanced epigenetic landscape, shedding light on distinct chromatin states associated with heterochromatin formation and gene silencing. Moreover, we elucidated the intricate interplay between these chromatin states and the overall topology of the genome. Employing a genetic approach, we delved into the role of the histone modification H3K9ac in genome topology. Notably, our investigation revealed that the ectopic deposition of this chromatin mark triggered a reorganization of the 3D chromatin structure, defining different TAD-like borders. Our work emphasizes the critical role of H3K9ac in shaping the topology of the tomato genome, providing valuable insights into the epigenetic landscape of this agriculturally significant crop species.

An improved genome assembly of Chrysanthemum nankingense reveals expansion and functional diversification of terpene synthase gene family

BACKGROUND

Terpenes are important components of plant aromas, and terpene synthases (TPSs) are the key enzymes driving terpene diversification. In this study, we characterized the volatile terpenes in five different Chrysanthemum nankingense tissues. In addition, genome-wide identification and expression analysis of TPS genes was conducted utilizing an improved chromosome-scale genome assembly and tissue-specific transcriptomes. The biochemical functions of three representative TPSs were also investigated.

RESULTS

We identified tissue-specific volatile organic compound (VOC) and volatile terpene profiles. The improved Chrysanthemum nankingense genome assembly was high-quality, including a larger assembled size (3.26 Gb) and a better contig N50 length (3.18 Mb) compared to the old version. A total of 140 CnTPS genes were identified, with the majority representing the TPS-a and TPS-b subfamilies. The chromosomal distribution of these TPS genes was uneven, and 26 genes were included in biosynthetic gene clusters. Closely-related Chrysanthemum taxa were also found to contain diverse TPS genes, and the expression profiles of most CnTPSs were tissue-specific. The three investigated CnTPS enzymes exhibited versatile activities, suggesting multifunctionality.

CONCLUSIONS

We systematically characterized the structure and diversity of TPS genes across the Chrysanthemum nankingense genome, as well as the potential biochemical functions of representative genes. Our results provide a basis for future studies of terpene biosynthesis in chrysanthemums, as well as for the breeding of improved chrysanthemum varieties.

Haplotype-resolved chromosome-level genome assembly of Ehretia macrophylla

Ehretia macrophylla Wall, known as wild loquat, is an ecologically, economically, and medicinally significant tree species widely grown in China, Japan, Vietnam, and Nepal. In this study, we have successfully generated a haplotype-resolved chromosome-scale genome assembly of E. macrophylla by integrating PacBio HiFi long-reads, Illumina short-reads, and Hi-C data. The genome assembly consists of two haplotypes, with sizes of 1.82 Gb and 1.58 Gb respectively, and contig N50 lengths of 28.11 Mb and 21.57 Mb correspondingly. Additionally, 99.41% of the assembly was successfully anchored into 40 pseudo-chromosomes. We predicted 58,886 protein-coding genes, of which 99.60% were functionally annotated from databases. We furthermore detected 2.65 Gb repeat sequences, 659,290 rRNAs, 4,931 tRNAs and 4,688 other ncRNAs. The high-quality assembly of the genome offers a solid basis for furthering the fields of molecular breeding and functional genomics of E. macrophylla.

Evolutionary insights into the organization of chromatin structure and landscape of transcriptional regulation in plants

Development of complex traits necessitates the functioning and coordination of intricate regulatory networks involving multiple genes. Understanding 3D chromatin structure can facilitate insight into the regulation of gene expression by regulatory elements. This potential, of visualizing the role of chromatin organization in the evolution and function of regulatory elements, remains largely unexplored. Here, we describe new perspectives that arise from the dual considerations of sequence variation of regulatory elements and chromatin structure, with a special focus on whole-genome doubling or polyploidy. We underscore the significance of hierarchical chromatin organization in gene regulation during evolution. In addition, we describe strategies for exploring chromatin organization in future investigations of regulatory evolution in plants, enabling insights into the evolutionary influence of regulatory elements on gene expression and, hence, phenotypes.

Multiple strategies, including 6mA methylation, affecting plant alternative splicing in allopolyploid peanut

Alternative splicing (AS), an important post-transcriptional regulation mechanism in eukaryotes, can significantly increase transcript diversity and contribute to gene expression regulation and many other complicated developmental processes. While plant gene AS events are well described, few studies have investigated the comprehensive regulation machinery of plant AS. Here, we use multi-omics to analyse peanut AS events. Using long-read isoform sequencing, 146 464 full-length non-chimeric transcripts were obtained, resulting in annotation corrections for 1782 genes and the identification of 4653 new loci. Using Iso-Seq RNA sequences, 271 776 unique splice junctions were identified, 82.49% of which were supported by transcriptome data. We characterized 50 977 polyadenylation sites for 23 262 genes, 12 369 of which had alternative polyadenylation sites. AS allows differential regulation of the same gene by miRNAs at the isoform level coupled with polyadenylation. In addition, we identified many long non-coding RNAs and fusion transcripts. There is a suppressed effect of 6mA on AS and gene expression. By analysis of chromatin structures, the genes located in the boundaries of topologically associated domains, proximal chromosomal telomere regions, inter- or intra-chromosomal loops were found to have more unique splice isoforms, higher expression, lower 6mA and more transposable elements (TEs) in their gene bodies than the other genes, indicating that chromatin interaction, 6mA and TEs play important roles in AS and gene expression. These results greatly refine the peanut genome annotation and contribute to the study of gene expression and regulation in peanuts. This work also showed AS is associated with multiple strategies for gene regulation.

Insight into chromatin compaction and spatial organization in rice interphase nuclei

Chromatin organization and its interactions are essential for biological processes, such as DNA repair, transcription, and DNA replication. Detailed cytogenetics data on chromatin conformation, and the arrangement and mutual positioning of chromosome territories in interphase nuclei are still widely missing in plants. In this study, level of chromatin condensation in interphase nuclei of rice (Oryza sativa) and the distribution of chromosome territories (CTs) were analyzed. Super-resolution, stimulated emission depletion (STED) microscopy showed different levels of chromatin condensation in leaf and root interphase nuclei. 3D immuno-FISH experiments with painting probes specific to chromosomes 9 and 2 were conducted to investigate their spatial distribution in root and leaf nuclei. Six different configurations of chromosome territories, including their complete association, weak association, and complete separation, were observed in root meristematic nuclei, and four configurations were observed in leaf nuclei. The volume of CTs and frequency of their association varied between the tissue types. The frequency of association of CTs specific to chromosome 9, containing NOR region, is also affected by the activity of the 45S rDNA locus. Our data suggested that the arrangement of chromosomes in the nucleus is connected with the position and the size of the nucleolus.

Light control of three-dimensional chromatin organization in soybean

Higher-order chromatin structure is critical for regulation of gene expression. In plants, light profoundly affects the morphogenesis of emerging seedlings as well as global gene expression to ensure optimal adaptation to environmental conditions. However, the changes and functional significance of chromatin organization in response to light during seedling development are not well documented. We constructed Hi-C contact maps for the cotyledon, apical hook and hypocotyl of soybean subjected to dark and light conditions. The resulting high-resolution Hi-C contact maps identified chromosome territories, A/B compartments, A/B sub-compartments, TADs (Topologically Associated Domains) and chromatin loops in each organ. We observed increased chromatin compaction under light and we found that domains that switched from B sub-compartments in darkness to A sub-compartments under light contained genes that were activated during photomorphogenesis. At the local scale, we identified a group of TADs constructed by gene clusters consisting of different numbers of Small Auxin-Upregulated RNAs (SAURs), which exhibited strict co-expression in the hook and hypocotyl in response to light stimulation. In the hypocotyl, RNA polymerase II (RNAPII) regulated the transcription of a SAURs cluster under light via TAD condensation. Our results suggest that the 3D genome is involved in the regulation of light-related gene expression in a tissue-specific manner.

Maintenance and dynamic reprogramming of chromatin organization during development

Controlled transcription of genes is critical for cell differentiation and development. Gene expression regulation therefore involves a multilayered control from nucleosome composition in histone variants and their post-translational modifications to higher-order folding of chromatin fibers and chromatin interactions in nuclear space. Recent technological advances have allowed gaining insight into these mechanisms, the interplay between local and higher-order chromatin organization, and the dynamic changes that occur during stress response and developmental transitions. In this review, we will discuss chromatin organization from the nucleosome to its three-dimensional structure in the nucleus, and consider how these different layers of organization are maintained during the cell cycle or rapidly reprogrammed during development.

Linker histone H1 drives heterochromatin condensation via phase separation in Arabidopsis

In the eukaryotic nucleus, heterochromatin forms highly condensed, visible foci known as heterochromatin foci (HF). These HF are enriched with linker histone H1, a key player in heterochromatin condensation and silencing. However, it is unknown how H1 aggregates HF and condenses heterochromatin. In this study, we established that H1 facilitates heterochromatin condensation by enhancing inter- and intrachromosomal interactions between and within heterochromatic regions of the Arabidopsis (Arabidopsis thaliana) genome. We demonstrated that H1 drives HF formation via phase separation, which requires its C-terminal intrinsically disordered region (C-IDR). A truncated H1 lacking the C-IDR fails to form foci or recover HF in the h1 mutant background, whereas C-IDR with a short stretch of the globular domain (18 out of 71 amino acids) is sufficient to rescue both defects. In addition, C-IDR is essential for H1's roles in regulating nucleosome repeat length and DNA methylation in Arabidopsis, indicating that phase separation capability is required for chromatin functions of H1. Our data suggest that bacterial H1-like proteins, which have been shown to condense DNA, are intrinsically disordered and capable of mediating phase separation. Therefore, we propose that phase separation mediated by H1 or H1-like proteins may represent an ancient mechanism for condensing chromatin and DNA.

Mapping nucleosome-resolution chromatin organization and enhancer-promoter loops in plants using Micro-C-XL

Although chromatin organizations in plants have been dissected at the scales of compartments and topologically associating domain (TAD)-like domains, there remains a gap in resolving fine-scale structures. Here, we use Micro-C-XL, a high-throughput chromosome conformation capture (Hi-C)-based technology that involves micrococcal nuclease (instead of restriction enzymes) and long cross-linkers, to dissect single nucleosome-resolution chromatin organization in Arabidopsis. Insulation analysis reveals more than 14,000 boundaries, which mostly include chromatin accessibility, epigenetic modifications, and transcription factors. Micro-C-XL reveals associations between RNA Pols and local chromatin organizations, suggesting that gene transcription substantially contributes to the establishment of local chromatin domains. By perturbing Pol II both genetically and chemically at the gene level, we confirm its function in regulating chromatin organization. Visible loops and stripes are assigned to super-enhancers and their targeted genes, thus providing direct insights for the identification and mechanistic analysis of distal CREs and their working modes in plants. We further investigate possible factors regulating these chromatin loops. Subsequently, we expand Micro-C-XL to soybean and rice. In summary, we use Micro-C-XL for analyses of plants, which reveal fine-scale chromatin organization and enhancer-promoter loops and provide insights regarding three-dimensional genomes in plants.

Chromatin remodeling analysis reveals the RdDM pathway responds to low-phosphorus stress in maize

Chromatin in eukaryotes folds into a complex three-dimensional (3D) structure that is essential for controlling gene expression and cellular function and is dynamically regulated in biological processes. Studies on plant phosphorus signaling have concentrated on single genes and gene interactions. It is critical to expand the existing signaling pathway in terms of its 3D structure. In this study, low-Pi treatment led to greater chromatin volume. Furthermore, low-Pi stress increased the insulation score and the number of TAD-like domains, but the effects on the A/B compartment were not obvious. The methylation levels of target sites (hereafter as RdDM levels) peaked at specific TAD-like boundaries, whereas RdDM peak levels at conserved TAD-like boundaries shifted and decreased sharply. The distribution pattern of RdDM sites originating from the Helitron transposons matched that of genome-wide RdDM sites near TAD-like boundaries. RdDM pathway genes were upregulated in the middle or early stages and downregulated in the later stages under low-Pi conditions. The RdDM pathway mutant ddm1a showed increased tolerance to low-Pi stress, with shortened and thickened roots contributing to higher Pi uptake from the shallow soil layer. ChIP-seq results revealed that ZmDDM1A could bind to Pi- and root development-related genes. Strong associations were found between interacting genes in significantly different chromatin-interaction regions and root traits. These findings not only expand the mechanisms by which plants respond to low-Pi stress through the RdDM pathway but also offer a crucial framework for the analysis of biological issues using 3D genomics.

Three-dimensional chromatin architecture in plants - General features and novelties

Research on the three-dimensional (3D) structure of the genome and its distribution within the nuclear space has made a big leap in the last two decades. Work in the animal field has led to significant advances in our general understanding on eukaryotic genome organization. This did not only bring along insights into how the 3D genome interacts with the epigenetic landscape and the transcriptional machinery but also how 3D genome architecture is relevant for fundamental developmental processes, such as cell differentiation. In parallel, the 3D organization of plant genomes have been extensively studied, which resulted in both congruent and novel findings, contributing to a more complete view on how eukaryotic genomes are organized in multiple dimensions. Plant genomes are remarkably diverse in size, composition, and ploidy. Furthermore, as intrinsically sessile organisms without the possibility to relocate to more favorable environments, plants have evolved an elaborate epigenetic repertoire to rapidly respond to environmental challenges. The diversity in genome organization and the complex epigenetic programs make plants ideal study subjects to acquire a better understanding on universal features and inherent constraints of genome organization. Furthermore, considering a wide range of species allows us to study the evolutionary crosstalk between the various levels of genome architecture. In this article, we aim at summarizing important findings on 3D genome architecture obtained in various plant species. These findings cover many aspects of 3D genome organization on a wide range of levels, from gene loops to topologically associated domains and to global 3D chromosome configurations. We present an overview on plant 3D genome organizational features that resemble those in animals and highlight facets that have only been observed in plants to date.

Accessible gene borders establish a core structural unit for chromatin architecture in Arabidopsis

Three-dimensional (3D) chromatin structure is linked to transcriptional regulation in multicellular eukaryotes including plants. Taking advantage of high-resolution Hi-C (high-throughput chromatin conformation capture), we detected a small structural unit with 3D chromatin architecture in the Arabidopsis genome, which lacks topologically associating domains, and also in the genomes of tomato, maize, and Marchantia polymorpha. The 3D folding domain unit was usually established around an individual gene and was dependent on chromatin accessibility at the transcription start site (TSS) and transcription end site (TES). We also observed larger contact domains containing two or more neighboring genes, which were dependent on accessible border regions. Binding of transcription factors to accessible TSS/TES regions formed these gene domains. We successfully simulated these Hi-C contact maps via computational modeling using chromatin accessibility as input. Our results demonstrate that gene domains establish basic 3D chromatin architecture units that likely contribute to higher-order 3D genome folding in plants.

Chromatin dynamics and subnuclear gene positioning for transcriptional regulation

Plants have been found to exhibit diverse characteristics and functions of chromatin organization, showing both similarities and differences to animals. It is becoming clear how chromatin organization is linked to transcriptional regulation in response to environmental stresses. Regulation of specific chromatin positions in the nuclear space is important for transcription, and the mechanisms that enable such chromatin dynamics are gradually being unveiled. Genes move between subdomains responsible for transcriptional activation or suppression in the subnuclear space in a gene repositioning cycle. We propose a model of localized chromatin interaction in nuclear subdomains, in which the dynamics of local chromatin interactions have a more important impact on the regulation of gene expression than large-scale chromatin organization. In this mini-review, we highlight recent findings on chromatin dynamics, particularly involving transcriptional regulation, and discuss future directions in the study of chromatin organization in plants.

Genome-wide high-throughput chromosome conformation capture analysis reveals hierarchical chromatin interactions during early somatic embryogenesis

Somatic embryogenesis (SE), like zygotic embryo development, is a progressive process. Early SE is the beginning of a switch from a somatic to an embryogenic state and is an important stage for initiating chromatin reprogramming of SE. Previous studies suggest that changes in chromatin accessibility occur during early SE, although information on the 3D structure of chromatin is not yet available. Here, we present a chromosome-level genome assembly of longan (Dimocarpus longan) using PacBio combined with high-through chromosome conformation capture scaffolding, which resulted in a 446 Mb genome assembly anchored onto 15 scaffolds. During early SE, chromatin was concentrated and then decondensed, and a large number of long terminal repeat retrotransposons (LTR-RTs) were enriched in the local chromatin interaction region, suggesting LTR-RTs were involved in chromatin reorganization. Early SE was accompanied by the transformation from A to B compartments, and the interactions between B compartments were enhanced. Results from chromatin accessibility, monomethylation of histone H3 at lysine 4 (H3K4me1) modification, and transcription analyses further revealed a gene regulatory network for cell wall thickening during SE. Particularly, we found that the H3K4me1 differential peak binding motif showed abnormal activation of ethylene response factor transcription factors and participation in SE. The chromosome-level genomic and multiomics analyses revealed the 3D conformation of chromatin during early SE, providing insight into the molecular mechanisms underlying cell wall thickening and the potential regulatory networks of TFs during early SE in D. longan. These results provide additional clues for revealing the molecular mechanisms of plant SE.

Histone H1 protects telomeric repeats from H3K27me3 invasion in Arabidopsis

While the pivotal role of linker histone H1 in shaping nucleosome organization is well established, its functional interplays with chromatin factors along the epigenome are just starting to emerge. Here we show that, in Arabidopsis, as in mammals, H1 occupies Polycomb Repressive Complex 2 (PRC2) target genes where it favors chromatin condensation and H3K27me3 deposition. We further show that, contrasting with its conserved function in PRC2 activation at genes, H1 selectively prevents H3K27me3 accumulation at telomeres and large pericentromeric interstitial telomeric repeat (ITR) domains by restricting DNA accessibility to Telomere Repeat Binding (TRB) proteins, a group of H1-related Myb factors mediating PRC2 cis recruitment. This study provides a mechanistic framework by which H1 avoids the formation of gigantic H3K27me3-rich domains at telomeric sequences and contributes to safeguard nucleus architecture.

3D organization of regulatory elements for transcriptional regulation in Arabidopsis

BACKGROUND

Although spatial organization of compartments and topologically associating domains at large scale is relatively well studied, the spatial organization of regulatory elements at fine scale is poorly understood in plants.

RESULTS

Here we perform high-resolution chromatin interaction analysis using paired-end tag sequencing approach. We map chromatin interactions tethered with RNA polymerase II and associated with heterochromatic, transcriptionally active, and Polycomb-repressive histone modifications in Arabidopsis. Analysis of the regulatory repertoire shows that distal active cis-regulatory elements are linked to their target genes through long-range chromatin interactions with increased expression of the target genes, while poised cis-regulatory elements are linked to their target genes through long-range chromatin interactions with depressed expression of the target genes. Furthermore, we demonstrate that transcription factor MYC2 is critical for chromatin spatial organization, and propose that MYC2 occupancy and MYC2-mediated chromatin interactions coordinately facilitate transcription within the framework of 3D chromatin architecture. Analysis of functionally related gene-defined chromatin connectivity networks reveals that genes implicated in flowering-time control are functionally compartmentalized into separate subdomains via their spatial activity in the leaf or shoot apical meristem, linking active mark- or Polycomb-repressive mark-associated chromatin conformation to coordinated gene expression.

CONCLUSION

The results reveal that the regulation of gene transcription in Arabidopsis is not only by linear juxtaposition, but also by long-range chromatin interactions. Our study uncovers the fine scale genome organization of Arabidopsis and the potential roles of such organization in orchestrating transcription and development.

Conserved H3K27me3-associated chromatin looping mediates physical interactions of gene clusters in plants

Higher-order chromatin organization is essential for transcriptional regulation, genome stability maintenance, and other genome functions. Increasing evidence has revealed significant differences in 3D chromatin organization between plants and animals. However, the extent, pattern, and rules of chromatin organization in plants are still unclear. In this study, we systematically identified and characterized long-range chromatin loops in the Arabidopsis 3D genome. We identified hundreds of long-range cis chromatin loops and found their anchor regions are closely associated with H3K27me3 epigenetic modifications. Furthermore, we demonstrated that these chromatin loops are dependent on Polycomb group (PcG) proteins, suggesting that the Polycomb repressive complex 2 (PRC2) complex is essential for establishing and maintaining these novel loops. Although most of these PcG-medicated chromatin loops are stable, many of these loops are tissue-specific or dynamically regulated by different treatments. Interestingly, tandemly arrayed gene clusters and metabolic gene clusters are enriched in anchor regions. Long-range H3K27me3-marked chromatin interactions are associated with the coregulation of specific gene clusters. Finally, we also identified H3K27me3-associated chromatin loops associated with gene clusters in Oryza sativa and Glycine max, indicating that these long-range chromatin loops are conserved in plants. Our results provide novel insights into genome evolution and transcriptional coregulation in plants.

Binding by the Polycomb complex component BMI1 and H2A monoubiquitination shape local and long-range interactions in the Arabidopsis genome

Three-dimensional (3D) chromatin organization is highly dynamic during development and seems to play a crucial role in regulating gene expression. Self-interacting domains, commonly called topologically associating domains (TADs) or compartment domains (CDs), have been proposed as the basic structural units of chromatin organization. Surprisingly, although these units have been found in several plant species, they escaped detection in Arabidopsis (Arabidopsis thaliana). Here, we show that the Arabidopsis genome is partitioned into contiguous CDs with different epigenetic features, which are required to maintain appropriate intra-CD and long-range interactions. Consistent with this notion, the histone-modifying Polycomb group machinery is involved in 3D chromatin organization. Yet, while it is clear that Polycomb repressive complex 2 (PRC2)-mediated trimethylation of histone H3 on lysine 27 (H3K27me3) helps establish local and long-range chromatin interactions in plants, the implications of PRC1-mediated histone H2A monoubiquitination on lysine 121 (H2AK121ub) are unclear. We found that PRC1, together with PRC2, maintains intra-CD interactions, but it also hinders the formation of H3K4me3-enriched local chromatin loops when acting independently of PRC2. Moreover, the loss of PRC1 or PRC2 activity differentially affects long-range chromatin interactions, and these 3D changes differentially affect gene expression. Our results suggest that H2AK121ub helps prevent the formation of transposable element/H3K27me1-rich long loops and serves as a docking point for H3K27me3 incorporation.

Toward cis-regulation in soybean: a 3D genome scope

In eukaryotic cells, 3D genome plays an important role in the regulation of gene spatiotemporal expression, which is essential for the biological and developmental processes in a life cycle. In the past decade, the development of high-throughput technologies greatly enhances our ability to map the 3D genome organization, identifies multiple 3D genome structures, and investigates the functional role of 3D genome organization in gene regulation, which facilitates our understandings of cis-regulatory landscape and biological development. Comparing with the comprehensive analyses of 3D genome in mammals and model plants, the progress in soybean is much less. Future development and application of tools to precisely manipulate 3D genome structure at different levels will significantly strengthen the functional genome study and molecular breeding in soybean. Here, we review the recent progresses in 3D genome study and discuss future directions, which may help to improve soybean 3D functional genome study and molecular breeding.

Reorganization of three-dimensional chromatin architecture in Medicago truncatula under phosphorus deficiency

Emerging evidence reveals that the three-dimensional (3D) chromatin architecture plays a key regulatory role in various biological processes of plants. However, information on the 3D chromatin architecture of the legume model plant Medicago truncatula and its potential roles in the regulation of response to mineral nutrient deficiency are very limited. Using high-resolution chromosome conformation capture sequencing, we identified the 3D genome structure of M. truncatula in terms of A/B compartments, topologically associated domains (TADs) and chromatin loops. The gene density, expressional level, and active histone modification were higher in A compartments than in B compartments. Moreover, we analysed the 3D chromatin architecture reorganization in response to phosphorus (P) deficiency. The intra-chromosomal cis-interaction proportion was increased by P deficiency, and a total of 748 A/B compartment switch regions were detected. In these regions, density changes in H3K4me3 and H3K27ac modifications were associated with expression of P deficiency-responsive genes involved in root system architecture and hormonal responses. Furthermore, these genes enhanced P uptake and mobilization by increasing root surface area and strengthening signal transduction under P deficiency. These findings advance our understanding of the potential roles of 3D chromatin architecture in responses of plants in general, and in particular in M. truncatula, to P deficiency.

Genomic rearrangements and evolutionary changes in 3D chromatin topologies in the cotton tribe (Gossypieae)

BACKGROUND

Analysis of the relationship between chromosomal structural variation (synteny breaks) and 3D-chromatin architectural changes among closely related species has the potential to reveal causes and correlates between chromosomal change and chromatin remodeling. Of note, contrary to extensive studies in animal species, the pace and pattern of chromatin architectural changes following the speciation of plants remain unexplored; moreover, there is little exploration of the occurrence of synteny breaks in the context of multiple genome topological hierarchies within the same model species.

RESULTS

Here we used Hi-C and epigenomic analyses to characterize and compare the profiles of hierarchical chromatin architectural features in representative species of the cotton tribe (Gossypieae), including Gossypium arboreum, Gossypium raimondii, and Gossypioides kirkii, which differ with respect to chromosome rearrangements. We found that (i) overall chromatin architectural territories were preserved in Gossypioides and Gossypium, which was reflected in their similar intra-chromosomal contact patterns and spatial chromosomal distributions; (ii) the non-random preferential occurrence of synteny breaks in A compartment significantly associate with the B-to-A compartment switch in syntenic blocks flanking synteny breaks; (iii) synteny changes co-localize with open-chromatin boundaries of topologically associating domains, while TAD stabilization has a greater influence on regulating orthologous expression divergence than do rearrangements; and (iv) rearranged chromosome segments largely maintain ancestral in-cis interactions.

CONCLUSIONS

Our findings provide insights into the non-random occurrence of epigenomic remodeling relative to the genomic landscape and its evolutionary and functional connections to alterations of hierarchical chromatin architecture, on a known evolutionary timescale.

Conserved noncoding sequences correlate with distant gene contacts in Arabidopsis and Brassica

Physical contact between genes distant on chromosomes is a potentially important way for genes to coordinate their expressions. To investigate the potential importance of distant contacts, we performed high-throughput chromatin conformation capture (Hi-C) experiments on leaf nuclei isolated from Brassica rapa and Brassica oleracea. We then combined our results with published Hi-C data from Arabidopsis thaliana. We found that distant genes come into physical contact and do so preferentially between the proximal promoter of one gene and the downstream region of another gene. Genes with higher numbers of conserved noncoding sequences (CNSs) nearby were more likely to have contact with distant genes. With more CNSs came higher numbers of transcription factor binding sites and more histone modifications associated with the activity. In addition, for the genes we studied, distant contacting genes with CNSs were more likely to be transcriptionally coordinated. These observations suggest that CNSs may enrich active histone modifications and recruit transcription factors, correlating with distant contacts to ensure coordinated expression. This study advances our knowledge of gene contacts and provides insights into the relationship between CNSs and distant gene contacts in plants.

A Comparison of Topologically Associating Domain Callers Based on Hi-C Data

Topologically associating domains (TADs) are local chromatin interaction domains, which have been shown to play an important role in gene expression regulation. TADs were originally discovered in the investigation of 3D genome organization based on High-throughput Chromosome Conformation Capture (Hi-C) data. Continuous considerable efforts have been dedicated to developing methods for detecting TADs from Hi-C data. Different computational methods for TADs identification vary in their assumptions and criteria in calling TADs. As a consequence, the TADs called by these methods differ in their similarities and biological features they are enriched in. In this work, we performed a systematic comparison of twenty-six TAD callers. We first compared the TADs and gaps between adjacent TADs across different methods, resolutions, and sequencing depths. We then assessed the quality of TADs and TAD boundaries according to three criteria: the decay of contact frequencies over the genomic distance, enrichment and depletion of regulatory elements around TAD boundaries, and reproducibility of TADs and TAD boundaries in replicate samples. Last, due to the lack of a gold standard of TADs, we also evaluated the performance of the methods on synthetic datasets. We discussed the key principles of TAD callers, and pinpointed current situation in the detection of TADs. We provide a concise, comprehensive, and systematic framework for evaluating the performance of TAD callers, and expect our work will provide useful guidance in choosing suitable approaches for the detection and evaluation of TADs.

Chromosome Conformation Capture of Mitotic Chromosomes

Despite more than a century of intensive study of mitotic chromosomes, their three-dimensional organization remains enigmatic. The last decade established Hi-C as a method of choice for study of spatial genome-wide interactions. Although its utilization has been focused mainly on studying genomic interactions in interphase nuclei, the method can be also successfully applied to study 3D architecture and genome folding in mitotic chromosomes. However, obtaining sufficient number of mitotic chromosomes as an input material and effective coupling with Hi-C method is challenging in plant species. An elegant way to overcome hindrances with obtaining a pure mitotic chromosome fraction is their isolation via flow cytometric sorting. This chapter presents a protocol describing plant sample preparation for chromosome conformation studies, for flow-sorting of plant mitotic metaphase chromosomes and for the Hi-C procedure.

CW198 acts as a genetic insulator to block enhancer-promoter interaction in plants

Insulators in vertebrates play a role in genome architecture and orchestrate temporo-spatial enhancer-promoter interactions. In plants, insulators and their associated binding factors have not been documented as of yet, largely as a result of a lack of characterized insulators. In this study, we took a comprehensive strategy to identify and validate the enhancer-blocking insulator CW198. We show that a 1.08-kb CW198 fragment from Arabidopsis can, when interposed between an enhancer and a promoter, efficiently abrogate the activation function of both constitutive and floral organ-specific enhancers in transgenic Arabidopsis and tobacco plants. In plants, both transcriptional crosstalk and spreading of histone modifications were rarely detectable across CW198, which resembles the insulation property observed across the CTCF insulator in the mammalian genome. Taken together, our findings support that CW198 acts as an enhancer-blocking insulator in both Arabidopsis and tobacco. The significance of the present findings and their relevance to the mitigation of mutual interference between enhancers and promoters, as well as multiple promoters in transgenes, is discussed.

Chromatin architectural alterations due to null mutation of a major CG methylase in rice

Associations between 3D chromatin architectures and epigenetic modifications have been characterized in animals. However, any impact of DNA methylation on chromatin architecture in plants is understudied, which is confined to Arabidopsis thaliana. Because plant species differ in genome size, composition, and overall chromatin packing, it is unclear to what extent findings from A. thaliana hold in other species. Moreover, the incomplete chromatin architectural profiles and the low-resolution high-throughput chromosome conformation capture (Hi-C) data from A. thaliana have hampered characterizing its subtle chromatin structures and their associations with DNA methylation. We constructed a high-resolution Hi-C interaction map for the null OsMET1-2 (the major CG methyltransferase in rice) mutant (osmet1-2) and isogenic wild-type rice (WT). Chromatin structural changes occurred in osmet1-2, including intra-/inter-chromosomal interactions, compartment transition, and topologically associated domains (TAD) variations. Our findings provide novel insights into the potential function of DNA methylation in TAD formation in rice and confirmed DNA methylation plays similar essential roles in chromatin packing in A. thaliana and rice.

Identification and functional validation of super-enhancers in Arabidopsis thaliana

Super-enhancers (SEs) are exceptionally large enhancers and are recognized to play prominent roles in cell identity in mammalian species. We surveyed the genomic regions containing large clusters of accessible chromatin regions (ACRs) marked by deoxyribonuclease (DNase) I hypersensitivity in Arabidopsis thaliana. We identified a set of 749 putative SEs, which have a minimum length of 1.5 kilobases and represent the top 2.5% of the largest ACR clusters. We demonstrate that the genomic regions associating with these SEs were more sensitive to DNase I than other nonpromoter ACRs. The SEs were preferentially associated with topologically associating domains. Furthermore, the SEs and their predicted cognate genes were frequently associated with organ development and tissue identity in A. thaliana. Therefore, the A. thaliana SEs and their cognate genes mirror the functional characteristics of those reported in mammalian species. We developed CRISPR/Cas-mediated deletion lines of a 3,578-bp SE associated with the thalianol biosynthetic gene cluster (BGC). Small deletions (131-157 bp) within the SE resulted in distinct phenotypic changes and transcriptional repression of all five thalianol genes. In addition, T-DNA insertions in the SE region resulted in transcriptional alteration of all five thalianol genes. Thus, this SE appears to play a central role in coordinating the operon-like expression pattern of the thalianol BGC.

Histone H2B.8 compacts flowering plant sperm through chromatin phase separation

Sperm chromatin is typically transformed by protamines into a compact and transcriptionally inactive state^1,2. Sperm cells of flowering plants lack protamines, yet they have small, transcriptionally active nuclei with chromatin condensed through an unknown mechanism^3,4. Here we show that a histone variant, H2B.8, mediates sperm chromatin and nuclear condensation in Arabidopsis thaliana. Loss of H2B.8 causes enlarged sperm nuclei with dispersed chromatin, whereas ectopic expression in somatic cells produces smaller nuclei with aggregated chromatin. This result demonstrates that H2B.8 is sufficient for chromatin condensation. H2B.8 aggregates transcriptionally inactive AT-rich chromatin into phase-separated condensates, which facilitates nuclear compaction without reducing transcription. Reciprocal crosses show that mutation of h2b.8 reduces male transmission, which suggests that H2B.8-mediated sperm compaction is important for fertility. Altogether, our results reveal a new mechanism of nuclear compaction through global aggregation of unexpressed chromatin. We propose that H2B.8 is an evolutionary innovation of flowering plants that achieves nuclear condensation compatible with active transcription.

Exploiting Genomic Features to Improve the Prediction of Transcription Factor-Binding Sites in Plants

The identification of transcription factor (TF) target genes is central in biology. A popular approach is based on the location by pattern matching of potential cis-regulatory elements (CREs). During the last few years, tools integrating next-generation sequencing data have been developed to improve the performance of pattern matching. However, such tools have not yet been comprehensively evaluated in plants. Hence, we developed a new streamlined method aiming at predicting CREs and target genes of plant TFs in specific organs or conditions. Our approach implements a supervised machine learning strategy, which allows decision rule models to be learnt using TF ChIP-chip/seq experimental data. Different layers of genomic features were integrated in predictive models: the position on the gene, the DNA sequence conservation, the chromatin state and various CRE footprints. Among the tested features, the chromatin features were crucial for improving the accuracy of the method. Furthermore, we evaluated the transferability of predictive models across TFs, organs and species. Finally, we validated our method by correctly inferring the target genes of key TFs controlling metabolite biosynthesis at the organ level in Arabidopsis. We developed a tool-Wimtrap-to reproduce our approach in plant species and conditions/organs for which ChIP-chip/seq data are available. Wimtrap is a user-friendly R package that supports an R Shiny web interface and is provided with pre-built models that can be used to quickly get predictions of CREs and TF gene targets in different organs or conditions in Arabidopsis thaliana, Solanum lycopersicum, Oryza sativa and Zea mays.

Mysteries of gene regulation: Promoters are not the sole triggers of gene expression

•

TF binding peaks were widely distributed in nonpromoters, especially downstream regions of transcription termination sites.

•

Exons of non-coding regions were the prominent regions of TF binding.

•

TAD boundaries were colocalized with activating histone marks and TF binding.

•

Genes with distinct functions demonstrated substantially different behaviors in cis-regulation and epigenetic signatures.

Cis-regulatory elements of promoters are essential for gene regulation by transcription factors (TFs). However, the regulatory roles of nonpromoter regions, TFs, and epigenetic marks remain poorly understood in plants. In this study, we characterized the cis-regulatory regions of 53 TFs and 19 histone marks in 328 chromatin immunoprecipitation (ChIP-seq) datasets from Arabidopsis. The genome-wide maps indicated that both promoters and regions around the transcription termination sites of protein-coding genes recruit the most TFs. The maps also revealed a diverse of histone combinations. The analysis suggested that exons play critical roles in the regulation of non-coding genes. Additionally, comparative analysis between heat-stress-responsive and nonresponsive genes indicated that the genes with distinct functions also exhibited substantial differences in cis-regulatory regions, histone regulation, and topologically associating domain (TAD) boundary organization. By integrating multiple high-throughput sequencing datasets, this study generated regulatory models for protein-coding genes, non-coding genes, and TAD boundaries to explain the complexity of transcriptional regulation.

DNA methylation in transposable elements buffers the connection between three-dimensional chromatin organization and gene transcription upon rice genome duplication

INTRODUCTION

Polyploidy is a major force in plant evolution and the domestication of cultivated crops.

OBJECTIVES

The study aimed to explore the relationship and underlying mechanism between three-dimensional (3D) chromatin organization and gene transcription upon rice genome duplication.

METHODS

The 3D chromatin structures between diploid (2C) and autotetraploid (4C) rice were compared using high-throughput chromosome conformation capture (Hi-C) analysis. The study combined genetics, transcriptomics, whole-genome bisulfite sequencing (WGBS-seq) and 3D genomics approaches to uncover the mechanism for DNA methylation in modulating gene transcription through 3D chromatin architectures upon rice genome duplication.

RESULTS

We found that 4C rice presents weakened intra-chromosomal interactions compared to its 2C progenitor in some chromosomes. In addition, we found that changes of 3D chromatin organizations including chromatin compartments, topologically associating domains (TADs), and loops, are uncorrelated with gene transcription. Moreover, DNA methylations in the regulatory sequences of genes in compartment A/B switched regions and TAD boundaries are unrelated to their expression. Importantly, although there was no significant difference in the methylation levels in transposable elements (TEs) in differentially expressed gene (DEG) and non-DEG promoters between 2C and 4C rice, we found that the hypermethylated TEs across genes in compartment A/B switched regions and TAD boundaries may suppress the expression of these genes.

CONCLUSION

The study proposed that the rice genome doubling might modulate TE methylation to buffer the effects of chromatin architecture on gene transcription in compartment A/B switched regions and TAD boundaries, resulting in the disconnection between 3D chromatin structure alteration and gene transcription upon rice genome duplication.

Plant biosynthetic gene clusters in the context of metabolic evolution

Covering: up to 2022Plants produce a wide range of structurally and biosynthetically diverse natural products to interact with their environment. These specialised metabolites typically evolve in limited taxonomic groups presumably in response to specific selective pressures. With the increasing availability of sequencing data, it has become apparent that in many cases the genes encoding biosynthetic enzymes for specialised metabolic pathways are not randomly distributed on the genome. Instead they are physically linked in structures such as arrays, pairs and clusters. The exact function of these clusters is debated. In this review we take a broad view of gene arrangement in plant specialised metabolism, examining types of structures and variation. We discuss the evolution of biosynthetic gene clusters in the wider context of metabolism, populations and epigenetics. Finally, we synthesise our observations to propose a new hypothesis for biosynthetic gene cluster formation in plants.

Topoisomerase VI participates in an insulator-like function that prevents H3K9me2 spreading

The organization of the genome into transcriptionally active and inactive chromatin domains requires well-delineated chromatin boundaries and insulator functions in order to maintain the identity of adjacent genomic loci with antagonistic chromatin marks and functionality. In plants that lack known chromatin insulators, the mechanisms that prevent heterochromatin spreading into euchromatin remain to be identified. Here, we show that DNA Topoisomerase VI participates in a chromatin boundary function that safeguards the expression of genes in euchromatin islands within silenced heterochromatin regions. While some transposable elements are reactivated in mutants of the Topoisomerase VI complex, genes insulated in euchromatin islands within heterochromatic regions of the Arabidopsis thaliana genome are specifically down-regulated. H3K9me2 levels consistently increase at euchromatin island loci and decrease at some transposable element loci. We further show that Topoisomerase VI physically interacts with S-adenosylmethionine synthase methionine adenosyl transferase 3 (MAT3), which is required for H3K9me2. A Topoisomerase VI defect affects MAT3 occupancy on heterochromatic elements and its exclusion from euchromatic islands, thereby providing a possible mechanistic explanation to the essential role of Topoisomerase VI in the delimitation of chromatin domains.

Chromatin remodeling complexes regulate genome architecture in Arabidopsis

In eukaryotes, three-dimensional (3D) chromatin architecture maintains genome stability and is important in regulating gene transcription. However, little is known about the mechanisms by which diverse ATP-dependent chromatin remodeling complexes regulate the 3D chromatin structure in plants. We examined the 3D chromatin structure within the ATPase subunit of the SWI/SNF, ISWI, INO80, and CHD remodeling complexes in wild-type (WT) and mutant Arabidopsis thaliana plants by combining high-throughput sequencing with in situ Hi-C, the enrichment of histone marks, nucleosome density, and gene expression. We found that compartment regions switched and compartmental strength was significantly weakened in all four enzyme mutants. Chromatin remodeling complexes differentially regulated the nucleosome distribution pattern and density within the switching compartments. Alterations of nucleosome distribution pattern and density were associated with a reduction in H3K27me3 levels in the chromatin remodeling enzyme mutants and led to compartment switching. Our data show that chromatin remodeling complexes regulate the linear nucleosome distribution pattern and density to promote H3K27me3 deposition, which in turn regulates 3D chromatin structure.

ChIP-Hub provides an integrative platform for exploring plant regulome

Plant genomes encode a complex and evolutionary diverse regulatory grammar that forms the basis for most life on earth. A wealth of regulome and epigenome data have been generated in various plant species, but no common, standardized resource is available so far for biologists. Here, we present ChIP-Hub, an integrative web-based platform in the ENCODE standards that bundles >10,000 publicly available datasets reanalyzed from >40 plant species, allowing visualization and meta-analysis. We manually curate the datasets through assessing ~540 original publications and comprehensively evaluate their data quality. As a proof of concept, we extensively survey the co-association of different regulators and construct a hierarchical regulatory network under a broad developmental context. Furthermore, we show how our annotation allows to investigate the dynamic activity of tissue-specific regulatory elements (promoters and enhancers) and their underlying sequence grammar. Finally, we analyze the function and conservation of tissue-specific promoters, enhancers and chromatin states using comparative genomics approaches. Taken together, the ChIP-Hub platform and the analysis results provide rich resources for deep exploration of plant ENCODE. ChIP-Hub is available at https://biobigdata.nju.edu.cn/ChIPHub/.

A comprehensive data portal to explore plant regulomes is still unavailable. Here, the authors develop a web-based platform ChIP-Hub in the ENCODE standards and demonstrate its applications in the identification of hierarchical regulatory network, tissue-specific chromatin dynamics, putative enhancers and chromatin states.

Spatial Features and Functional Implications of Plant 3D Genome Organization

The advent of high-throughput sequencing-based methods for chromatin conformation, accessibility, and immunoprecipitation assays has been a turning point in 3D genomics. Altogether, these new tools have been pushing upward the interpretation of pioneer cytogenetic evidence for a higher order in chromatin packing. Here, we review the latest development in our understanding of plant spatial genome structures and different levels of organization and discuss their functional implications. Then, we spotlight the complexity of organellar (i.e., mitochondria and plastids) genomes and discuss their 3D packing into nucleoids. Finally, we propose unaddressed research axes to investigate functional links between chromatin-like dynamics and transcriptional regulation within organellar nucleoids.

A long noncoding RNA functions in high-light-induced anthocyanin accumulation in apple by activating ethylene synthesis

Anthocyanin production in apple (Malus domestica) fruit and their consequent coloration can be induced by high-light treatment. The hormone ethylene is also essential for this coloration, but the regulatory relationships that link ethylene and light with anthocyanin-associated coloration are not well defined. In this study, we observed that high-light treatment of apple fruit increased anthocyanin accumulation more than moderate-light treatment did and was the main contributor of induced ethylene production and activation of anthocyanin biosynthesis. A transcriptome study of light-treated apple fruit suggested that a long noncoding RNA (lncRNA), MdLNC610, the corresponding gene of which is physically located downstream from the 1-aminocyclopropane-1-carboxylate oxygenase (ACO) ethylene biosynthesis gene MdACO1, likely affects anthocyanin biosynthesis under high-light treatment. Expression and promoter β-glucuronidase reporter analyses further showed that MdLNC610 upregulates expression of MdACO1 and so likely participates in high-light-induced ethylene biosynthesis. Overexpression of MdACO1 and MdLNC610 in apple fruit and calli indicated that a major increase in MdLNC610 expression activates MdACO1 expression, thereby causing an increase in ethylene production and anthocyanin levels. These results suggest that MdLNC610 participates in the regulation of high-light-induced anthocyanin production by functioning as a positive regulator to promote MdACO1 expression and ethylene biosynthesis. Our study provides insights into the relationship between mRNA and lncRNA networks in the ethylene biosynthetic pathway and anthocyanin accumulation in apple fruit.

Redox status of the plant cell determines epigenetic modifications under abiotic stress conditions and during developmental processes

BACKGROUND

The oxidation-reduction (redox) status of the cell influences or regulates transcription factors and enzymes involved in epigenetic changes, such as DNA methylation, histone protein modifications, and chromatin structure and remodeling. These changes are crucial regulators of chromatin architecture, leading to differential gene expression in eukaryotes. But the cell's redox homeostasis is difficult to sustain since the production of reactive oxygen species (ROS) and reactive nitrogen species (RNS) is not equal in plants at different developmental stages and under abiotic stress conditions. Exceeding optimum ROS and RNS levels leads to oxidative stress and thus alters the redox status of the cell. Consequently, this alteration modulates intracellular epigenetic modifications that either mitigate or mediate the plant growth and stress response.

AIM OF REVIEW

Recent studies suggest that the altered redox status of the cell reform the cellular functions and epigenetic changes. Recent high-throughput techniques have also greatly advanced redox-mediated gene expression discovery, but the integrated view of the redox status, and its associations with epigenetic changes and subsequent gene expression in plants are still scarce. In this review, we accordingly focus on how the redox status of the cell affects epigenetic modifications in plants under abiotic stress conditions and during developmental processes. This is a first comprehensive review on the redox status of the cell covering the redox components and signaling, redox status alters the post-translational modification of proteins, intracellular epigenetic modifications, redox interplay during DNA methylation, redox regulation of histone acetylation and methylation, redox regulation of miRNA biogenesis, redox regulation of chromatin structure and remodeling and conclusion, future perspectives and biotechnological opportunities for the future development of the plants.

KEY SCIENTIFIC CONCEPTS OF REVIEW

The interaction of redox mediators such as ROS, RNS and antioxidants regulates redox homeostasis and redox-mediated epigenetic changes. We discuss how redox mediators modulate epigenetic changes and show the opportunities for smart use of the redox status of the cell in plant development and abiotic stress adaptation. However, how a redox mediator triggers epigenetic modification without activating other redox mediators remains yet unknown.

A comparison of topologically associating domain callers over mammals at high resolution

BACKGROUND

Topologically associating domains (TADs) are locally highly-interacting genome regions, which also play a critical role in regulating gene expression in the cell. TADs have been first identified while investigating the 3D genome structure over High-throughput Chromosome Conformation Capture (Hi-C) interaction dataset. Substantial degree of efforts have been devoted to develop techniques for inferring TADs from Hi-C interaction dataset. Many TAD-calling methods have been developed which differ in their criteria and assumptions in TAD inference. Correspondingly, TADs inferred via these callers vary in terms of both similarities and biological features they are enriched in.

RESULT

We have carried out a systematic comparison of 27 TAD-calling methods over mammals. We use Micro-C, a recent high-resolution variant of Hi-C, to compare TADs at a very high resolution, and classify the methods into 3 categories: feature-based methods, Clustering methods, Graph-partitioning methods. We have evaluated TAD boundaries, gaps between adjacent TADs, and quality of TADs across various criteria. We also found particularly CTCF and Cohesin proteins to be effective in formation of TADs with corner dots. We have also assessed the callers performance on simulated datasets since a gold standard for TADs is missing. TAD sizes and numbers change remarkably between TAD callers and dataset resolutions, indicating that TADs are hierarchically-organized domains, instead of disjoint regions. A core subset of feature-based TAD callers regularly perform the best while inferring reproducible domains, which are also enriched for TAD related biological properties.

CONCLUSION

We have analyzed the fundamental principles of TAD-calling methods, and identified the existing situation in TAD inference across high resolution Micro-C interaction datasets over mammals. We come up with a systematic, comprehensive, and concise framework to evaluate the TAD-calling methods performance across Micro-C datasets. Our research will be useful in selecting appropriate methods for TAD inference and evaluation based on available data, experimental design, and biological question of interest. We also introduce our analysis as a benchmarking tool with publicly available source code.

Reshuffling of the ancestral core-eudicot genome shaped chromatin topology and epigenetic modification in Panax

All extant core-eudicot plants share a common ancestral genome that has experienced cyclic polyploidizations and (re)diploidizations. Reshuffling of the ancestral core-eudicot genome generates abundant genomic diversity, but the role of this diversity in shaping the hierarchical genome architecture, such as chromatin topology and gene expression, remains poorly understood. Here, we assemble chromosome-level genomes of one diploid and three tetraploid Panax species and conduct in-depth comparative genomic and epigenomic analyses. We show that chromosomal interactions within each duplicated ancestral chromosome largely maintain in extant Panax species, albeit experiencing ca. 100–150 million years of evolution from a shared ancestor. Biased genetic fractionation and epigenetic regulation divergence during polyploidization/(re)diploidization processes generate remarkable biochemical diversity of secondary metabolites in the Panax genus. Our study provides a paleo-polyploidization perspective of how reshuffling of the ancestral core-eudicot genome leads to a highly dynamic genome and to the metabolic diversification of extant eudicot plants.

The role of polyploidization generated genomic diversity in shaping the hierarchical genome architecture remains unclear. Here, the authors show that repatterning of the ancestral eudicot genome has resulted in multi-dimensional genome plasticity and secondary metabolite diversification via comparisons of Panax genomes.

Chromosome-Scale Assembly of the Dendrobium nobile Genome Provides Insights Into the Molecular Mechanism of the Biosynthesis of the Medicinal Active Ingredient of Dendrobium

Orchids constitute approximately 10% of flowering plant species. However, only about 10 orchid genomes have been published. Metabolites are the main way through which orchids respond to their environment. Dendrobium nobile, belonging to Dendrobium, the second largest genus in Orchidaceae, has high ornamental, medicinal, and ecological value. D. nobile is the source of many popular horticultural varieties. Among the Dendrobium species, D. nobile has the highest amount of dendrobine, which is regarded as one of the criteria for evaluating medicinal quality. Due to lack of data and analysis at the genomic level, the biosynthesis pathways of dendrobine and other related medicinal ingredients in D. nobile are unknown. In this paper, we report a chromosome-scale reference genome of D. nobile to facilitate the investigation of its genomic characteristics for comparison with other Dendrobium species. The assembled genome size of D. nobile was 1.19 Gb. Of the sequences, 99.45% were anchored to 19 chromosomes. Furthermore, we identified differences in gene number and gene expression patterns compared with two other Dendrobium species by integrating whole-genome sequencing and transcriptomic analysis [e.g., genes in the polysaccharide biosynthesis pathway and upstream of the alkaloid (dendrobine) biosynthesis pathway]. Differences in the TPS and CYP450 gene families were also found among orchid species. All the above differences might contribute to the species-specific medicinal ingredient biosynthesis pathways. The metabolic pathway-related analysis will provide further insight into orchid responses to the environment. Additionally, the reference genome will provide important insights for further molecular elucidation of the medicinal active ingredients of Dendrobium and enhance the understanding of orchid evolution.

Rice RS2-9, which is bound by transcription factor OSH1, blocks enhancer-promoter interactions in plants

Insulators characterized in Drosophila and mammals have been shown to play a key role in the restriction of promiscuous enhancer-promoter interactions, as well as reshaping the topological landscape of chromosomes. Yet the role of insulators in plants remains poorly understood, in large part because of a lack of well-characterized insulators and binding factor(s). In this study, we isolated a 1.2-kb RS2-9 insulator from the Oryza sativa (rice) genome that can, when interposed between an enhancer and promoter, efficiently block the activation function of both constitutive and floral organ-specific enhancers in transgenic Arabidopsis and Nicotiana tabacum (tobacco). In the rice genome, the genes flanking RS2-9 exhibit an absence of mutual transcriptional interactions, as well as a lack of histone modification spread. We further determined that O. sativa Homeobox 1 (OSH1) bound two regions of RS2-9, as well as over 50 000 additional sites in the rice genome, the majority of which resided in intergenic regions. Mutation of one of the two OSH1-binding sites in RS2-9 impaired insulation activity by up to 60%, whereas the mutation of both binding sites virtually abolished insulator function. We also demonstrated that OSH1 binding sites were associated with 72% of the boundaries of topologically associated domains (TADs) identified in the rice genome, which is comparable to the 77% of TAD boundaries bound by the insulator CCCTC-binding factor (CTCF) in mammals. Taken together, our findings indicate that OSH1-RS2-9 acts as a true insulator in plants, and highlight a potential role for OSH1 in gene insulation and topological organization in plant genomes.

Diversification of chromatin organization in eukaryotes

Our knowledge of the chromatin landscape and its regulation was originally discovered using yeast and a limited number of animals models. A wealth of studies in model plants now strongly demonstrates the conservation of certain features while illuminating the diversification of others. Here we summarize recent advances that describe specific features of chromatin organization of transcriptional units and specific regulation of heterochromatin in flowering plants. We highlight the importance of transcriptional regulation in plant chromatin organization and the need to investigate a more diverse range of species to understand the chromatin landscape in eukaryotes.

Conformational Changes in Three-Dimensional Chromatin Structure in Paulownia fortunei After Phytoplasma Infection

Higher-order chromatin structures play important roles in regulating multiple biological processes such as growth and development as well as biotic and abiotic stress response. However, little is known about three-dimensional chromatin structures in Paulownia or about whole-genome chromatin conformational changes that occur in response to Paulownia witches' broom (PaWB) disease. We used high-throughput chromosome conformation capture (Hi-C) to obtain genome-wide profiles of chromatin conformation in both healthy and phytoplasma-infected Paulownia fortunei genome. The heat map results indicated that the strongest interactions between chromosomes were in the telomeres. We confirmed that the main structural characteristics of A/B compartments, topologically associated domains, and chromatin loops were prominent in the Paulownia genome and were clearly altered in phytoplasma-infected plants. By combining chromatin immunoprecipitation sequencing, Hi-C signals, and RNA sequencing data, we inferred that the chromatin structure changed and the modification levels of three histones (H3K4me3/K9ac/K36me3) increased in phytoplasma-infected P. fortunei, which was associated with changes of transcriptional activity. We concluded that for epigenetic modifications, transcriptional activity might function in combination to shape chromatin packing in healthy and phytoplasm-infected Paulownia. Finally, 11 genes (e.g., RPN6, Sec61 subunit-α) that were commonly located at specific topologically associated domain boundaries, A/B compartment switching and specific loops, and had been associated with histone marks were identified and considered as closely related to PaWB stress. Our results provide new insights into the nexus between gene regulation and chromatin conformational alterations in nonmodel plants upon phytopathogen infection and plant disease resistance.

Molecular properties of epimutation hotspots

Mistakes in the maintenance of CG methylation are a source of heritable epimutations in plants. Multigenerational surveys indicate that the rate of these stochastic events varies substantially across the genome, with some regions harbouring localized 'epimutation hotspots'. Using Arabidopsis as a model, we show that epimutation hotspots are indexed by a specific set of chromatin states that map to subregions of gene body methylation genes. Although these regions comprise only ~12% of all CGs in the genome, they account for ~63% of all epimutation events per unit time. Molecular profiling revealed that these regions contain unique sequence features, harbour steady-state intermediate methylation levels and act as putative targets of antagonistic DNA methylation pathways. We further demonstrate that experimentally induced shifts in steady-state methylation in these hotspot regions are sufficient to significantly alter local epimutation intensities. Our work lays the foundation for dissecting the molecular mechanisms and evolutionary consequences of epimutation hotspots in plants.