3D Chromatin Architecture of Large Plant Genomes Determined by Local A/B Compartments

A method for chromatin domain partitioning based on hypergraph clustering

For many years, multi-scale models of chromatin domains, such as A/B compartments, sub-compartments, topologically associated domains (TADs), sub-TADs, and loops have been popular. However, existing methods can only identify structures at a single scale and cannot partition multi-scale structures. In this paper, we proposed a method (TORNADOES) for chromatin domain partitioning based on hypergraph clustering. First, we use a density clustering algorithm to identify TADs at different scales based on Hi-C data with different resolutions. Then, by combining ChIP-seq data features and TAD results at different scales, we generate a hypergraph based on these TADs. Finally, we partition the chromatin domain structure at different scales, including A/B, A1, A2, B1, B2, and B3 based on the Laplacian matrix feature of the hypergraph. Similarity comparison experiments and ChIP-seq signal enrichment analysis are performed on the A/B region and sub-TAD levels, respectively, demonstrating that our method can identify chromatin domains with distinct features and provide a deeper understanding of the organizational patterns and functional differences in TADs at the genomic hierarchical structure. Comparative analysis of multiple cell line data shows that TORNADOES can better classify different numbers and types of compartments by changing the factors ChIP-seq data and clustering number used to characterize TAD compared to other methods. Source code for the TORNADOES method can be found at https://github.com/ghaiyan/TORNADOES.

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.

The jet-like chromatin structure defines active secondary metabolism in fungi

Eukaryotic genomes are spatially organized within the nucleus in a nonrandom manner. However, fungal genome arrangement and its function in development and adaptation remain largely unexplored. Here, we show that the high-order chromosome structure of Fusarium graminearum is sculpted by both H3K27me3 modification and ancient genome rearrangements. Active secondary metabolic gene clusters form a structure resembling chromatin jets. We demonstrate that these jet-like domains, which can propagate symmetrically for 54 kb, are prevalent in the genome and correlate with active gene transcription and histone acetylation. Deletion of GCN5, which encodes a core and functionally conserved histone acetyltransferase, blocks the formation of the domains. Insertion of an exogenous gene within the jet-like domain significantly augments its transcription. These findings uncover an interesting link between alterations in chromatin structure and the activation of fungal secondary metabolism, which could be a general mechanism for fungi to rapidly respond to environmental cues, and highlight the utility of leveraging three-dimensional genome organization in improving gene transcription in eukaryotes.

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.

Genome organization and botanical diversity

The rich diversity of angiosperms, both the planet's dominant flora and the cornerstone of agriculture, is integrally intertwined with a distinctive evolutionary history. Here, we explore the interplay between angiosperm genome organization and botanical diversity, empowered by genomic approaches ranging from genetic linkage mapping to analysis of gene regulation. Commonality in the genetic hardware of plants has enabled robust comparative genomics that has provided a broad picture of angiosperm evolution and implicated both general processes and specific elements in contributing to botanical diversity. We argue that the hardware of plant genomes-both in content and in dynamics-has been shaped by selection for rather substantial differences in gene regulation between plants and animals such as maize and human, organisms of comparable genome size and gene number. Their distinctive genome content and dynamics may reflect in part the indeterminate development of plants that puts strikingly different demands on gene regulation than in animals. Repeated polyploidization of plant genomes and multiplication of individual genes together with extensive rearrangement and differential retention provide rich raw material for selection of morphological and/or physiological variations conferring fitness in specific niches, whether natural or artificial. These findings exemplify the burgeoning information available to employ in increasing knowledge of plant biology and in modifying selected plants to better meet human needs.

Chromosome-scale genome assembly of bread wheat's wild relative Triticum timopheevii

Wheat (Triticum aestivum) is one of the most important food crops with an urgent need for increase in its production to feed the growing world. Triticum timopheevii (2n = 4x = 28) is an allotetraploid wheat wild relative species containing the At and G genomes that has been exploited in many pre-breeding programmes for wheat improvement. In this study, we report the generation of a chromosome-scale reference genome assembly of T. timopheevii accession PI 94760 based on PacBio HiFi reads and chromosome conformation capture (Hi-C). The assembly comprised a total size of 9.35 Gb, featuring a contig N50 of 42.4 Mb and included the mitochondrial and plastid genome sequences. Genome annotation predicted 166,325 gene models including 70,365 genes with high confidence. DNA methylation analysis showed that the G genome had on average more methylated bases than the At genome. In summary, the T. timopheevii genome assembly provides a valuable resource for genome-informed discovery of agronomically important genes for food security.

Integrative multiomics profiling of passion fruit reveals the genetic basis for fruit color and aroma

Passion fruit (Passiflora edulis) possesses a complex aroma and is widely grown in tropical and subtropical areas. Here, we conducted the de novo assembly, annotation, and comparison of PPF (P. edulis Sims) and YPF (P. edulis f. flavicarpa) reference genomes using PacBio, Illumina, and Hi-C technologies. Notably, we discovered evidence of recent whole-genome duplication events in P. edulis genomes. Comparative analysis revealed 7.6∼8.1 million single nucleotide polymorphisms, 1 million insertions/deletions, and over 142 Mb presence/absence variations among different P. edulis genomes. During the ripening of yellow passion fruit, metabolites related to flavor, aroma, and color were substantially accumulated or changed. Through joint analysis of genomic variations, differentially expressed genes, and accumulated metabolites, we explored candidate genes associated with flavor, aroma, and color distinctions. Flavonoid biosynthesis pathways, anthocyanin biosynthesis pathways, and related metabolites are pivotal factors affecting the coloration of passion fruit, and terpenoid metabolites accumulated more in PPF. Finally, by heterologous expression in yeast (Saccharomyces cerevisiae), we functionally characterized 12 terpene synthases. Our findings revealed that certain TPS homologs in both YPF and PPF varieties produce identical terpene products, while others yield distinct compounds or even lose their functionality. These discoveries revealed the genetic and metabolic basis of unique characteristics in aroma and flavor between the 2 passion fruit varieties. This study provides resources for better understanding the genome architecture and accelerating genetic improvement of passion fruits.

Cohesin mediated loop extrusion from active enhancers form chromatin jets in C. elegans

In mammals, cohesin and CTCF organize the 3D genome into topologically associated domains (TADs) to regulate communication between cis-regulatory elements. However, many organisms, including S. cerevisiae , C. elegans , and A. thaliana lack CTCF. Here, we use C. elegans as a model to investigate the function of cohesin in 3D genome organization in an animal without CTCF. We use auxin-inducible degradation to acutely deplete SMC-3 or its negative regulator WAPL-1 from somatic cells. Using Hi-C data, we identify a cohesin-dependent 3D genome organization feature called chromatin jets (aka fountains). These are population average reflections of DNA loops that are ∼20-40 kb in scale and often cover a few transcribed genes. The jets emerge from NIPBL occupied segments, and the trajectory of the jets coincides with cohesin binding. Cohesin translocation from jet origins depends on a fully intact complex and is extended upon WAPL-1 depletion. Hi-C results support the idea that cohesin is preferentially loaded at NIPBL occupied sites and loop extrudes in an effectively two-sided manner. The location of putative loading sites coincide with active enhancers and the strength of chromatin jet pattern correlates with transcription. Hi-C analyses upon WAPL-1 depletion reveal unequal loop extrusion processivity on each side and stalling due to cohesin molecules colliding. Compared to mammalian systems, average processivity of C. elegans cohesin is ∼10-fold shorter and NIPBL binding does not depend on cohesin. We conclude that the processivity of cohesin scales with genome size; and regardless of CTCF presence, preferential loading of cohesin at enhancers is a conserved mechanism of genome organization that regulates the interaction of gene regulatory elements in 3D.

Chromosome compartmentalization: causes, changes, consequences, and conundrums

The spatial segregation of the genome into compartments is a major feature of 3D genome organization. New data on mammalian chromosome organization across different conditions reveal important information about how and why these compartments form and change. A combination of epigenetic state, nuclear body tethering, physical forces, gene expression, and replication timing (RT) can all influence the establishment and alteration of chromosome compartments. We review the causes and implications of genomic regions undergoing a 'compartment switch' that changes their physical associations and spatial location in the nucleus. About 20-30% of genomic regions change compartment during cell differentiation or cancer progression, whereas alterations in response to a stimulus within a cell type are usually much more limited. However, even a change in 1-2% of genomic bins may have biologically relevant implications. Finally, we review the effects of compartment changes on gene regulation, DNA damage repair, replication, and the physical state of the cell.

HiCrayon reveals distinct layers of multi-state 3D chromatin organization

The co-visualization of chromatin conformation with 1D 'omics data is key to the multi-omics driven data analysis of 3D genome organization. Chromatin contact maps are often shown as 2D heatmaps and visually compared to 1D genomic data by simple juxtaposition. While common, this strategy is imprecise, placing the onus on the reader to align features with each other. To remedy this, we developed HiCrayon, an interactive tool that facilitates the integration of 3D chromatin organization maps and 1D datasets. This visualization method integrates data from genomic assays directly into the chromatin contact map by coloring interactions according to 1D signal. HiCrayon is implemented using R shiny and python to create a graphical user interface (GUI) application, available in both web or containerized format to promote accessibility. HiCrayon is implemented in R, and includes a graphical user interface (GUI), as well as a slimmed-down web-based version that lets users quickly produce publication-ready images. We demonstrate the utility of HiCrayon in visualizing the effectiveness of compartment calling and the relationship between ChIP-seq and various features of chromatin organization. We also demonstrate the improved visualization of other 3D genomic phenomena, such as differences between loops associated with CTCF/cohesin vs. those associated with H3K27ac. We then demonstrate HiCrayon's visualization of organizational changes that occur during differentiation and use HiCrayon to detect compartment patterns that cannot be assigned to either A or B compartments, revealing a distinct 3rd chromatin compartment. Overall, we demonstrate the utility of co-visualizing 2D chromatin conformation with 1D genomic signals within the same matrix to reveal fundamental aspects of genome organization. Local version: https://github.com/JRowleyLab/HiCrayon Web version: https://jrowleylab.com/HiCrayon.

Increased long-distance and homo-trans interactions related to H3K27me3 in Arabidopsis hybrids

In plants, the genome structure of hybrids changes compared with their parents, but the effects of these changes in hybrids remain elusive. Comparing reciprocal crosses between Col × C24 and C24 × Col in Arabidopsis using high-throughput chromosome conformation capture assay (Hi-C) analysis, we found that hybrid three-dimensional (3D) chromatin organization had more long-distance interactions relative to parents, and this was mainly located in promoter regions and enriched in genes with heterosis-related pathways. The interactions between euchromatin and heterochromatin were increased, and the compartment strength decreased in hybrids. In compartment domain (CD) boundaries, the distal interactions were more in hybrids than their parents. In the hybrids of CURLY LEAF (clf) mutants clfCol  × clfC24 and clfC24  × clfCol , the heterosis phenotype was damaged, and the long-distance interactions in hybrids were fewer than in their parents with lower H3K27me3. ChIP-seq data revealed higher levels of H3K27me3 in the region adjacent to the CD boundary and the same interactional homo-trans sites in the wild-type (WT) hybrids, which may have led to more long-distance interactions. In addition, the differentially expressed genes (DEGs) located in the boundaries of CDs and loop regions changed obviously in WT, and the functional enrichment for DEGs was different between WT and clf in the long-distance interactions and loop regions. Our findings may therefore propose a new epigenetic explanation of heterosis in the Arabidopsis hybrids and provide new insights into crop breeding and yield increase.

Conserved chromatin and repetitive patterns reveal slow genome evolution in frogs

Frogs are an ecologically diverse and phylogenetically ancient group of anuran amphibians that include important vertebrate cell and developmental model systems, notably the genus Xenopus. Here we report a high-quality reference genome sequence for the western clawed frog, Xenopus tropicalis, along with draft chromosome-scale sequences of three distantly related emerging model frog species, Eleutherodactylus coqui, Engystomops pustulosus, and Hymenochirus boettgeri. Frog chromosomes have remained remarkably stable since the Mesozoic Era, with limited Robertsonian (i.e., arm-preserving) translocations and end-to-end fusions found among the smaller chromosomes. Conservation of synteny includes conservation of centromere locations, marked by centromeric tandem repeats associated with Cenp-a binding surrounded by pericentromeric LINE/L1 elements. This work explores the structure of chromosomes across frogs, using a dense meiotic linkage map for X. tropicalis and chromatin conformation capture (Hi-C) data for all species. Abundant satellite repeats occupy the unusually long (~20 megabase) terminal regions of each chromosome that coincide with high rates of recombination. Both embryonic and differentiated cells show reproducible associations of centromeric chromatin and of telomeres, reflecting a Rabl-like configuration. Our comparative analyses reveal 13 conserved ancestral anuran chromosomes from which contemporary frog genomes were constructed.

3D genome structural variations play important roles in regulating seed oil content of Brassica napus

Dissecting the complex regulatory mechanism of seed oil content (SOC) is one of the main research goals in Brassica napus. Increasing evidence suggests that genome architecture is linked to multiple biological functions. However, the effect of genome architecture on SOC regulation remains unclear. Here, we used high-throughput chromatin conformation capture to characterize differences in the three-dimensional (3D) landscape of genome architecture of seeds from two B. napus lines, N53-2 (with high SOC) and Ken-C8 (with low SOC). Bioinformatics analysis demonstrated that differentially accessible regions and differentially expressed genes between N53-2 and Ken-C8 were preferentially enriched in regions with quantitative trait loci (QTLs)/associated genomic regions (AGRs) for SOC. A multi-omics analysis demonstrated that expression of SOC-related genes was tightly correlated with genome structural variations in QTLs/AGRs of B. napus. The candidate gene BnaA09g48250D, which showed structural variation in a QTL/AGR on chrA09, was identified by fine-mapping of a KN double-haploid population derived from hybridization of N53-2 and Ken-C8. Overexpression and knockout of BnaA09g48250D led to significant increases and decreases in SOC, respectively, in the transgenic lines. Taken together, our results reveal the 3D genome architecture of B. napus seeds and the roles of genome structural variations in SOC regulation, enriching our understanding of the molecular mechanisms of SOC regulation from the perspective of spatial chromatin structure.

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.

MINI-AC: inference of plant gene regulatory networks using bulk or single-cell accessible chromatin profiles

Gene regulatory networks (GRNs) represent the interactions between transcription factors (TF) and their target genes. Plant GRNs control transcriptional programs involved in growth, development, and stress responses, ultimately affecting diverse agricultural traits. While recent developments in accessible chromatin (AC) profiling technologies make it possible to identify context-specific regulatory DNA, learning the underlying GRNs remains a major challenge. We developed MINI-AC (Motif-Informed Network Inference based on Accessible Chromatin), a method that combines AC data from bulk or single-cell experiments with TF binding site (TFBS) information to learn GRNs in plants. We benchmarked MINI-AC using bulk AC datasets from different Arabidopsis thaliana tissues and showed that it outperforms other methods to identify correct TFBS. In maize, a crop with a complex genome and abundant distal AC regions, MINI-AC successfully inferred leaf GRNs with experimentally confirmed, both proximal and distal, TF-target gene interactions. Furthermore, we showed that both AC regions and footprints are valid alternatives to infer AC-based GRNs with MINI-AC. Finally, we combined MINI-AC predictions from bulk and single-cell AC datasets to identify general and cell-type specific maize leaf regulators. Focusing on C4 metabolism, we identified diverse regulatory interactions in specialized cell types for this photosynthetic pathway. MINI-AC represents a powerful tool for inferring accurate AC-derived GRNs in plants and identifying known and novel candidate regulators, improving our understanding of gene regulation in plants.

Evolutionary insights into 3D genome organization and epigenetic landscape of Vigna mungo

Eukaryotic genomes show an intricate three-dimensional (3D) organization within the nucleus that regulates multiple biological processes including gene expression. Contrary to animals, understanding of 3D genome organization in plants remains at a nascent stage. Here, we investigate the evolution of 3D chromatin architecture in legumes. By using cutting-edge PacBio, Illumina, and Hi-C contact reads, we report a gap-free, chromosome-scale reference genome assembly of Vigna mungo, an important minor legume cultivated in Southeast Asia. We spatially resolved V. mungo chromosomes into euchromatic, transcriptionally active A compartment and heterochromatic, transcriptionally-dormant B compartment. We report the presence of TAD-like-regions throughout the diagonal of the HiC matrix that resembled transcriptional quiescent centers based on their genomic and epigenomic features. We observed high syntenic breakpoints but also high coverage of syntenic sequences and conserved blocks in boundary regions than in the TAD-like region domains. Our findings present unprecedented evolutionary insights into spatial 3D genome organization and epigenetic patterns and their interaction within the V. mungo genome. This will aid future genomics and epigenomics research and breeding programs of V. mungo.

Genome-wide DNase I-hypersensitive site assay reveals distinct genomic distributions and functional features of open chromatin in autopolyploid sugarcane

The characterization of cis-regulatory DNA elements (CREs) is essential for deciphering the regulation of gene expression in eukaryotes. Although there have been endeavors to identify CREs in plants, the properties of CREs in polyploid genomes are still largely unknown. Here, we conducted the genome-wide identification of DNase I-hypersensitive sites (DHSs) in leaf and stem tissues of the auto-octoploid species Saccharum officinarum. We revealed that DHSs showed highly similar distributions in the genomes of these two S. officinarum tissues. Notably, we observed that approximately 74% of DHSs were located in distal intergenic regions, suggesting considerable differences in the abundance of distal CREs between S. officinarum and other plants. Leaf- and stem-dependent transcriptional regulatory networks were also developed by mining the binding motifs of transcription factors (TFs) from tissue-specific DHSs. Four TEOSINTE BRANCHED 1, CYCLOIDEA, and PCF1 (TCP) TFs (TCP2, TCP4, TCP7, and TCP14) and two ethylene-responsive factors (ERFs) (ERF109 and ERF03) showed strong causal connections with short binding distances from each other, pointing to their possible roles in the regulatory networks of leaf and stem development. Through functional validation in transiently transgenic protoplasts, we isolate a set of tissue-specific promoters. Overall, the DHS maps presented here offer a global view of the potential transcriptional regulatory elements in polyploid sugarcane and can be expected to serve as a valuable resource for both transcriptional network elucidation and genome editing in sugarcane breeding.

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.

HSFA1a: the quarterback of heat stress response and 3D-chromatin organization

In eukaryotes, the highly complex chromatin exhibits structural dynamism to modulate cellular responses. Recently, Huang et al. have shown a novel role of the heat stress master regulator, HSFA1a, in transiently reorganizing the 3D chromatin to promote distal and proximal enhancer-promoter contacts to regulate stress-responsive genes in tomato.

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.

Chromatin architecture of two different strains of Artemisia annua reveals the alterations in interaction and gene expression

MAIN CONCLUSION

The hierarchical architecture of chromatins affects the gene expression level of glandular secreting trichomes and the artemisinin biosynthetic pathway-related genes, consequently bringing on huge differences in the content of artemisinin and its derivatives of A. annua. The plant of traditional Chinese medicine "Qinghao" is called Artemisia annua L. in Chinese Pharmacopoeia. High content and the total amount of artemisinin is the main goal of A. annua breeding, nevertheless, the change of chromatin organization during the artemisinin synthesis process has not been discovered yet. This study intended to find the roles of chromatin structure in the production of artemisinin through bioinformatics and experimental validation. Chromosome conformation capture analysis was used to scrutinize the interactions among chromosomes and categorize various scales of chromatin during artemisinin synthesis in A. annua. To confirm the effect of the changes in chromatin structure, Hi-C and RNA-sequencing were performed on two different strains to find the correlation between chromatin structure and gene expression levels on artemisinin synthesis progress and regulation. Our results revealed that the frequency of intra-chromosomal interactions was higher in the inter-chromosomal interactions between the root and leaves on a high artemisinin production strain (HAP) compared to a low artemisinin production strain (LAP). We found that compartmental transition was connected with interactions among different chromatins. Interestingly, glandular secreting trichomes (GSTs) and the artemisinin biosynthetic pathway (ABP) related genes were enriched in the areas which have the compartmental transition, reflecting the regulation of artemisinin synthesis. Topologically associated domain boundaries were associated with various distributions of genes and expression levels. Genes associated with ABP and GST in the adjacent loop were highly expressed, suggesting that epigenetic regulation plays an important role during artemisinin synthesis and glandular secreting trichomes production process. Chromatin structure could show an important status in the mechanisms of artemisinin synthesis process in A. annua.

A multi-omic Nicotiana benthamiana resource for fundamental research and biotechnology

Nicotiana benthamiana is an invaluable model plant and biotechnology platform with a ~3 Gb allotetraploid genome. To further improve its usefulness and versatility, we have produced high-quality chromosome-level genome assemblies, coupled with transcriptome, epigenome, microRNA and transposable element datasets, for the ubiquitously used LAB strain and a related wild accession, QLD. In addition, single nucleotide polymorphism maps have been produced for a further two laboratory strains and four wild accessions. Despite the loss of five chromosomes from the ancestral tetraploid, expansion of intergenic regions, widespread segmental allopolyploidy, advanced diploidization and evidence of recent bursts of Copia pseudovirus (Copia) mobility not seen in other Nicotiana genomes, the two subgenomes of N. benthamiana show large regions of synteny across the Solanaceae. LAB and QLD have many genetic, metabolic and phenotypic differences, including disparate RNA interference responses, but are highly interfertile and amenable to genome editing and both transient and stable transformation. The LAB/QLD combination has the potential to be as useful as the Columbia-0/Landsberg errecta partnership, utilized from the early pioneering days of Arabidopsis genomics to today.

Low temperature-induced regulatory network rewiring via WRKY regulators during banana peel browning

Banana (Musa spp.) fruits, as typical tropical fruits, are cold sensitive, and lower temperatures can disrupt cellular compartmentalization and lead to severe browning. How tropical fruits respond to low temperature compared to the cold response mechanisms of model plants remains unknown. Here, we systematically characterized the changes in chromatin accessibility, histone modifications, distal cis-regulatory elements, transcription factor binding, and gene expression levels in banana peels in response to low temperature. Dynamic patterns of cold-induced transcripts were generally accompanied by concordant chromatin accessibility and histone modification changes. These upregulated genes were enriched for WRKY binding sites in their promoters and/or active enhancers. Compared to banana peel at room temperature, large amounts of banana WRKYs were specifically induced by cold and mediated enhancer-promoter interactions regulating critical browning pathways, including phospholipid degradation, oxidation, and cold tolerance. This hypothesis was supported by DNA affinity purification sequencing, luciferase reporter assays, and transient expression assay. Together, our findings highlight widespread transcriptional reprogramming via WRKYs during banana peel browning at low temperature and provide an extensive resource for studying gene regulation in tropical plants in response to cold stress, as well as potential targets for improving cold tolerance and shelf life of tropical fruits.

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.

Remodeling of the 3D chromatin architecture in the marine microalga Nannochloropsis oceanica during lipid accumulation

BACKGROUND

Genomic three-dimensional (3D) spatial organization plays a key role in shaping gene expression and associated chromatin modification, and it is highly sensitive to environmental stress conditions. In microalgae, exposure to nitrogen stress can drive lipid accumulation, yet the associated functional alterations in the spatial organization of the microalgal genome have yet to be effectively characterized.

RESULTS

Accordingly, the present study employed RNA-seq, Hi-C, and ChIP-seq approaches to explore the relationship between 3D chromosomal architecture and gene expression during lipid accumulation in the marine microalga Nannochloropsis oceanica in response to nitrogen deprivation (ND). These analyses revealed that ND resulted in various changes in chromosomal organization, including A/B compartment transitions, topologically associating domain (TAD) shifts, and the disruption of short-range interactions. Significantly higher levels of gene expression were evident in A compartments and TAD boundary regions relative to B compartments and TAD interior regions, consistent with observed histone modification enrichment in these areas. ND-induced differentially expressed genes (DEGs) were notably enriched in altered TAD-associated regions and regions exhibiting differential genomic contact. These DEGs were subjected to Gene Ontology (GO) term analyses that indicated they were enriched in the 'fatty acid metabolism', 'response to stress', 'carbon fixation' and 'photosynthesis' functional categories, in line with the ND treatment conditions used to conduct this study. These data indicate that Nannochloropsis cells exhibit a clear association between chromatin organization and transcriptional activity under nitrogen stress conditions. Pronounced and extensive histone modifications were evident in response to ND. Observed changes in chromatin architecture were linked to shifts in histone modifications and gene expression.

CONCLUSIONS

Overall, the reprogramming of many lipid metabolism-associated genes was evident under nitrogen stress conditions with respect to both histone modifications and chromosomal organization. Together these results revealed that higher-order chromatin architecture represents a new layer that can guide efforts to understand the transcriptional regulation of lipid metabolism in nitrogen-deprived microalgae.

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.

Epigenetic regulation of nuclear processes in fungal plant pathogens

Through the association of protein complexes to DNA, the eukaryotic nuclear genome is broadly organized into open euchromatin that is accessible for enzymes acting on DNA and condensed heterochromatin that is inaccessible. Chemical and physical alterations to chromatin may impact its organization and functionality and are therefore important regulators of nuclear processes. Studies in various fungal plant pathogens have uncovered an association between chromatin organization and expression of in planta-induced genes that are important for pathogenicity. This review discusses chromatin-based regulation mechanisms as determined in the fungal plant pathogen Verticillium dahliae and relates the importance of epigenetic transcriptional regulation and other nuclear processes more broadly in fungal plant pathogens.

Shifting landscapes: the role of 3D genomic organizations in gene regulatory strategies

3D genome folding enables the physical storage of chromosomes into the compact volume of a cell's nucleus, allows for the accurate segregation of chromatin to daughter cells, and has been shown to be tightly coupled to the way in which genetic information is converted into transcriptional programs [1-3]. Importantly, this link between chromatin architecture and gene regulation is a selectable feature in which modifications to chromatin organization accompany, or perhaps even drive the establishment of new regulatory strategies with enduring impacts on animal body plan complexity. Here, we discuss the nature of different 3D genome folding systems found across the tree of life, with particular emphasis on metazoans, and the relative influence of these systems on gene regulation. We suggest how the properties of these folding systems have influenced regulatory strategies employed by different lineages and may have catalyzed the partitioning and specialization of genetic programs that enabled multicellularity and organ-grade body plan complexity.

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.

Structural diversity and stress regulation of the plant immunity-associated CALMODULIN-BINDING PROTEIN 60 (CBP60) family of transcription factors in Solanum lycopersicum (tomato)

Cellular signaling generates calcium (Ca2+) ions, which are ubiquitous secondary messengers decoded by calcium-dependent protein kinases, calcineurins, calreticulin, calmodulins (CAMs), and CAM-binding proteins. Previous studies in the model plant Arabidopsis thaliana have shown the critical roles of the CAM-BINDING PROTEIN 60 (CBP60) protein family in plant growth, stress responses, and immunity. Certain CBP60 factors can regulate plant immune responses, like pattern-triggered immunity, effector-triggered immunity, and synthesis of major plant immune-activating metabolites salicylic acid (SA) and N-hydroxypipecolic acid (NHP). Although homologous CBP60 sequences have been identified in the plant kingdom, their function and regulation in most species remain unclear. In this paper, we specifically characterized 11 members of the CBP60 family in the agriculturally important crop tomato (Solanum lycopersicum). Protein sequence analyses revealed that three CBP60 homologs have the closest amino acid identity to Arabidopsis CBP60g and SARD1, master transcription factors involved in plant immunity. Strikingly, AlphaFold deep learning-assisted prediction of protein structures highlighted close structural similarity between these tomato and Arabidopsis CBP60 homologs. Conserved domain analyses revealed that they possess CAM-binding domains and DNA-binding domains, reflecting their potential involvement in linking Ca2+ signaling and transcriptional regulation in tomato plants. In terms of their gene expression profiles under biotic (Pseudomonas syringae pv. tomato DC3000 pathogen infection) and/or abiotic stress (warming temperatures), five tomato CBP60 genes were pathogen-responsive and temperature-sensitive, reminiscent of Arabidopsis CBP60g and SARD1. Overall, we present a genome-wide identification of the CBP60 gene/protein family in tomato plants, and we provide evidence on their regulation and potential function as Ca2+-sensing transcriptional regulators.

An upgraded method of high-throughput chromosome conformation capture (Hi-C 3.0) in cotton (Gossypium spp.)

High-throughput chromosome conformation capture (Hi-C) technology has been applied to explore the chromatin interactions and shed light on the biological functions of three-dimensional genomic features. However, it remains challenging to guarantee the high quality of Hi-C library in plants and hence the reliable capture of chromatin structures, especially loops, due to insufficient fragmentation and low efficiency of proximity ligations. To overcome these deficiencies, we optimized the parameters of the Hi-C protocol, principally the cross-linking agents and endonuclease fragmentation strategy. The double cross-linkers (FA+DSG) and double restriction enzymes (DpnII+DdeI) were utilized. Thus, a systematic in situ Hi-C protocol was designed using plant tissues embedded with comprehensive quality controls to monitor the library construction. This upgraded method, termed Hi-C 3.0, was applied to cotton leaves for trial. In comparison with the conventional Hi-C 2.0, Hi-C 3.0 can obtain more than 50% valid contacts at a given sequencing depth to improve the signal-to-noise ratio. Hi-C 3.0 can furthermore enhance the capturing of loops almost as twice as that of Hi-C 2.0. In addition, Hi-C 3.0 showed higher efficiency of compartment detection and identified compartmentalization more accurately. In general, Hi-C 3.0 contributes to the advancement of the Hi-C method in plants by promoting its capability on decoding the chromatin organization.

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.

Assignment of the somatic A/B compartments to chromatin domains in giant transcriptionally active lampbrush chromosomes

BACKGROUND

The three-dimensional configuration of the eukaryotic genome is an emerging area of research. Chromosome conformation capture outlined genome segregation into large scale A and B compartments corresponding mainly to transcriptionally active and repressive chromatin. It remains unknown how the compartmentalization of the genome changes in growing oocytes of animals with hypertranscriptional type of oogenesis. Such oocytes are characterized by highly elongated chromosomes, called lampbrush chromosomes, which acquire a typical chromomere-loop appearance, representing one of the classical model systems for exploring the structural and functional organization of chromatin domains.

RESULTS

Here, we compared the distribution of A/B compartments in chicken somatic cells with chromatin domains in lampbrush chromosomes. We found that in lampbrush chromosomes, the extended chromatin domains, restricted by compartment boundaries in somatic cells, disintegrate into individual chromomeres. Next, we performed FISH-mapping of the genomic loci, which belong to A or B chromatin compartments as well as to A/B compartment transition regions in embryonic fibroblasts on isolated lampbrush chromosomes. We found, that in chicken lampbrush chromosomes, clusters of dense compact chromomeres bearing short lateral loops and enriched with repressive epigenetic modifications generally correspond to constitutive B compartments in somatic cells. A compartments align with lampbrush chromosome segments with smaller, less compact chromomeres, longer lateral loops, and a higher transcriptional status. Clusters of small loose chromomeres with relatively long lateral loops show no obvious correspondence with either A or B compartment identity. Some genes belonging to facultative B (sub-) compartments can be tissue-specifically transcribed during oogenesis, forming distinct lateral loops.

CONCLUSIONS

Here, we established a correspondence between the A/B compartments in somatic interphase nucleus and chromatin segments in giant lampbrush chromosomes from diplotene stage oocytes. The chromomere-loop structure of the genomic regions corresponding to interphase A and B compartments reveals the difference in how they are organized at the level of chromatin domains. The results obtained also suggest that gene-poor regions tend to be packed into chromomeres.

Chromatin alternates between A and B compartments at kilobase scale for subgenic organization

Nuclear compartments are prominent features of 3D chromatin organization, but sequencing depth limitations have impeded investigation at ultra fine-scale. CTCF loops are generally studied at a finer scale, but the impact of looping on proximal interactions remains enigmatic. Here, we critically examine nuclear compartments and CTCF loop-proximal interactions using a combination of in situ Hi-C at unparalleled depth, algorithm development, and biophysical modeling. Producing a large Hi-C map with 33 billion contacts in conjunction with an algorithm for performing principal component analysis on sparse, super massive matrices (POSSUMM), we resolve compartments to 500 bp. Our results demonstrate that essentially all active promoters and distal enhancers localize in the A compartment, even when flanking sequences do not. Furthermore, we find that the TSS and TTS of paused genes are often segregated into separate compartments. We then identify diffuse interactions that radiate from CTCF loop anchors, which correlate with strong enhancer-promoter interactions and proximal transcription. We also find that these diffuse interactions depend on CTCF's RNA binding domains. In this work, we demonstrate features of fine-scale chromatin organization consistent with a revised model in which compartments are more precise than commonly thought while CTCF loops are more protracted.

Plant Tolerance to Drought Stress with Emphasis on Wheat

Environmental stresses, such as drought, have negative effects on crop yield. Drought is a stress whose impact tends to increase in some critical regions. However, the worldwide population is continuously increasing and climate change may affect its food supply in the upcoming years. Therefore, there is an ongoing effort to understand the molecular processes that may contribute to improving drought tolerance of strategic crops. These investigations should contribute to delivering drought-tolerant cultivars by selective breeding. For this reason, it is worthwhile to review regularly the literature concerning the molecular mechanisms and technologies that could facilitate gene pyramiding for drought tolerance. This review summarizes achievements obtained using QTL mapping, genomics, synteny, epigenetics, and transgenics for the selective breeding of drought-tolerant wheat cultivars. Synthetic apomixis combined with the msh1 mutation opens the way to induce and stabilize epigenomes in crops, which offers the potential of accelerating selective breeding for drought tolerance in arid and semi-arid regions.

Insights into the convergent evolution of fructan biosynthesis in angiosperms from the highly characteristic chicory genome

Fructans in angiosperms play essential roles in physiological functions and environmental adaptations. As a major source of industrial fructans (especially inulin-type), chicory (Cichorium intybus L.) is a model species for studying fructan biosynthesis. However, the genes underlying this process and their evolutionary history in angiosperms remain elusive. We combined multiple sequencing technologies to assemble and annotate the chicory genome and scan its (epi)genomic features, such as genomic components, DNA methylation, and three-dimensional (3D) structure. We also performed a comparative genomics analysis to uncover the associations between key traits and gene families. We achieved a nearly complete chicory genome assembly and found that continuous bursts of a few highly active retrotransposon families largely shaped the (epi)genomic characteristics. The highly methylated genome with its unique 3D structure potentially influences critical biological processes. Our comprehensive comparative genomics analysis deciphered the genetic basis for the rich sesquiterpene content in chicory and indicated that the fructan-accumulating trait resulted from convergent evolution in angiosperms due to shifts in critical sites of fructan-active enzymes. The highly characterized chicory genome provides insight into Asteraceae evolution and fructan biosynthesis in angiosperms.

DNA methylation signatures of duplicate gene evolution in angiosperms

Gene duplication is a source of evolutionary novelty. DNA methylation may play a role in the evolution of duplicate genes (paralogs) through its association with gene expression. While this relationship has been examined to varying extents in a few individual species, the generalizability of these results at either a broad phylogenetic scale with species of differing duplication histories or across a population remains unknown. We applied a comparative epigenomics approach to 43 angiosperm species across the phylogeny and a population of 928 Arabidopsis (Arabidopsis thaliana) accessions, examining the association of DNA methylation with paralog evolution. Genic DNA methylation was differentially associated with duplication type, the age of duplication, sequence evolution, and gene expression. Whole genome duplicates were typically enriched for CG-only gene-body methylated or unmethylated genes, while single-gene duplications were typically enriched for non-CG methylated or unmethylated genes. Non-CG methylation, in particular, was characteristic of more recent single-gene duplicates. Core angiosperm gene families differentiated into those which preferentially retain paralogs and 'duplication-resistant' families, which convergently reverted to singletons following duplication. Duplication-resistant families that still have paralogous copies were, uncharacteristically for core angiosperm genes, enriched for non-CG methylation. Non-CG methylated paralogs had higher rates of sequence evolution, higher frequency of presence-absence variation, and more limited expression. This suggests that silencing by non-CG methylation may be important to maintaining dosage following duplication and be a precursor to fractionation. Our results indicate that genic methylation marks differing evolutionary trajectories and fates between paralogous genes and have a role in maintaining dosage following duplication.

Considerations and caveats for analyzing chromatin compartments

Genomes are organized into nuclear compartments, separating active from inactive chromatin. Chromatin compartments are readily visible in a large number of species by experiments that map chromatin conformation genome-wide. When analyzing these maps, a common step is the identification of genomic intervals that interact within A (active) and B (inactive) compartments. It has also become increasingly common to identify and analyze subcompartments. We review different strategies to identify A/B and subcompartment intervals, including a discussion of various machine-learning approaches to predict these features. We then discuss the strengths and limitations of current strategies and examine how these aspects of analysis may have impacted our understanding of chromatin compartments.

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.

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.

Moderation of Structural DNA Properties by Coupled Dinucleotide Contents in Eukaryotes

Dinucleotides are known as determinants for various structural and physiochemical properties of DNA and for binding affinities of proteins to DNA. These properties (e.g., stiffness) and bound proteins (e.g., transcription factors) are known to influence important biological functions, such as transcription regulation and 3D chromatin organization. Accordingly, the question arises of how the considerable variations in dinucleotide contents of eukaryotic chromosomes could still provide consistent DNA properties resulting in similar functions and 3D conformations. In this work, we investigate the hypothesis that coupled dinucleotide contents influence DNA properties in opposite directions to moderate each other's influences. Analyzing all 2478 chromosomes of 155 eukaryotic species, considering bias from coding sequences and enhancers, we found sets of correlated and anti-correlated dinucleotide contents. Using computational models, we estimated changes of DNA properties resulting from this coupling. We found that especially pure A/T dinucleotides (AA, TT, AT, TA), known to influence histone positioning and AC/GT contents, are relevant moderators and that, e.g., the Roll property, which is known to influence histone affinity of DNA, is preferably moderated. We conclude that dinucleotide contents might indirectly influence transcription and chromatin 3D conformation, via regulation of histone occupancy and/or other mechanisms.

High-Resolution 3D Genome Map of Brucella Chromosomes in Exponential and Stationary Phases

The three-dimensional (3D) genome structure of an organism or cell is highly relevant to its biological activities, but the availability of 3D genome information for bacteria, especially intracellular pathogens, is still limited. Here, we used Hi-C (high-throughput chromosome conformation capture) technology to determine the 3D chromosome structures of exponential- and stationary-phase Brucella melitensis at a 1-kb resolution. We observed that the contact heat maps of the two B. melitensis chromosomes contain a prominent diagonal and a secondary diagonal. Then, 79 chromatin interaction domains (CIDs) were detected at an optical density at 600 nm (OD600) of 0.4 (exponential phase), with the longest CID being 106 kb and the shortest being 12 kb. Moreover, we obtained 49,363 significant cis-interaction loci and 59,953 significant trans-interaction loci. Meanwhile, 82 CIDs of B. melitensis at an OD600 of 1.5 (stationary phase) were detected, with the longest CID being 94 kb and the shortest being 16 kb. In addition, 25,965 significant cis-interaction loci and 35,938 significant trans-interaction loci were obtained in this phase. Furthermore, we found that as the B. melitensis cells grew from the logarithmic to the plateau phase, the frequency of short-range interactions increased, while that of long-range interactions decreased. Finally, combined analysis of 3D genome and whole-genome transcriptome (RNA-seq) data revealed that the strength of short-range interactions in Chr1 is specifically and strongly correlated with gene expression. Overall, our study provides a global view of the chromatin interactions in the B. melitensis chromosomes, which will serve as a resource for further study of the spatial regulation of gene expression in Brucella. IMPORTANCE The spatial structure of chromatin plays important roles in normal cell functions and in the regulation of gene expression. Three-dimensional genome sequencing has been performed in many mammals and plants, but the availability of such data for bacteria, especially intracellular pathogens, is still limited. Approximately 10% of sequenced bacterial genomes contain more than one replicon. However, how multiple replicons are organized within bacterial cells, how they interact, and whether these interactions help to maintain or segregate these multipartite genomes are unresolved issues. Brucella is a Gram-negative, facultative intracellular, and zoonotic bacterium. Except for Brucella suis biovar 3, Brucella species have two chromosomes. Here, we applied Hi-C technology to determine the 3D genome structures of exponential- and stationary-phase Brucella melitensis chromosomes at a 1-kb resolution. Combined analysis of the 3D genome and RNA-seq data indicated that the strength of short-range interactions in B. melitensis Chr1 is specifically and strongly correlated with gene expression. Our study provides a resource to achieve a deeper understanding of the spatial regulation of gene expression in Brucella.

An improved assembly of the "Cascade" hop (Humulus lupulus) genome uncovers signatures of molecular evolution and refines time of divergence estimates for the Cannabaceae family

We present a chromosome-level assembly of the Cascade hop (Humulus lupulus L. var. lupulus) genome. The hop genome is large (2.8 Gb) and complex, and early attempts at assembly were fragmented. Recent advances have made assembly of the hop genome more tractable, transforming the extent of investigation that can occur. The chromosome-level assembly of Cascade was developed by scaffolding the previously reported Cascade assembly generated with PacBio long-read sequencing and polishing with Illumina short-read DNA sequencing. We developed gene models and repeat annotations and used a controlled bi-parental mapping population to identify significant sex-associated markers. We assessed molecular evolution in gene sequences, gene family expansion and contraction, and time of divergence from Cannabis sativa and other closely related plant species using Bayesian inference. We identified the putative sex chromosome in the female genome based on significant sex-associated markers from the bi-parental mapping population. While the estimate of repeat content (~64%) is similar to the estimate for the hemp genome, syntenic blocks in hop contain a greater percentage of LTRs. Hop is enriched for disease resistance-associated genes in syntenic gene blocks and expanded gene families. The Cascade chromosome-level assembly will inform cultivation strategies and serve to deepen our understanding of the hop genomic landscape, benefiting hop researchers and the Cannabaceae genomics community.

Evolution of 3D chromatin organization at different scales

Most animal genomes fold in 3D chromatin domains called topologically associated domains (TADs) that facilitate interactions between cis-regulatory elements (CREs) and promoters. Owing to their critical role in the control of developmental gene expression, we explore how TADs have shaped animal evolution. In the light of recent studies that profile TADs in disparate animal lineages, we discuss their phylogenetic distribution and the mechanisms that underlie their formation. We present evidence indicating that TADs are plastic entities composed of genomic strata of different ages: ancient cores are combined with newer regions and brought into extant TADs through genomic rearrangements. We highlight that newly incorporated TAD strata enable the establishment of new CRE-promoter interactions and in turn new expression patterns that can drive phenotypical innovation. We further highlight how subtle changes in chromatin folding may fine-tune the expression levels of developmental genes and hold a potential for evolutionary significance.

HSFA1a modulates plant heat stress responses and alters the 3D chromatin organization of enhancer-promoter interactions

The complex and dynamic three-dimensional organization of chromatin within the nucleus makes understanding the control of gene expression challenging, but also opens up possible ways to epigenetically modulate gene expression. Because plants are sessile, they evolved sophisticated ways to rapidly modulate gene expression in response to environmental stress, that are thought to be coordinated by changes in chromatin conformation to mediate specific cellular and physiological responses. However, to what extent and how stress induces dynamic changes in chromatin reorganization remains poorly understood. Here, we comprehensively investigated genome-wide chromatin changes associated with transcriptional reprogramming response to heat stress in tomato. Our data show that heat stress induces rapid changes in chromatin architecture, leading to the transient formation of promoter-enhancer contacts, likely driving the expression of heat-stress responsive genes. Furthermore, we demonstrate that chromatin spatial reorganization requires HSFA1a, a transcription factor (TF) essential for heat stress tolerance in tomato. In light of our findings, we propose that TFs play a key role in controlling dynamic transcriptional responses through 3D reconfiguration of promoter-enhancer contacts.

Pan-3D genome analysis reveals structural and functional differentiation of soybean genomes

BACKGROUND

High-order chromatin structure plays important roles in gene regulation. However, the diversity of the three-dimensional (3D) genome across plant accessions are seldom reported.

RESULTS

Here, we perform the pan-3D genome analysis using Hi-C sequencing data from 27 soybean accessions and comprehensively investigate the relationships between 3D genomic variations and structural variations (SVs) as well as gene expression. We find that intersection regions between A/B compartments largely contribute to compartment divergence. Topologically associating domain (TAD) boundaries in A compartments exhibit significantly higher density compared to those in B compartments. Pan-3D genome analysis shows that core TAD boundaries have the highest transcription start site (TSS) density and lowest GC content and repeat percentage. Further investigation shows that non-long terminal repeat (non-LTR) retrotransposons play important roles in maintaining TAD boundaries, while Gypsy elements and satellite repeats are associated with private TAD boundaries. Moreover, presence and absence variation (PAV) is found to be the major contributor to 3D genome variations. Nevertheless, approximately 55% of 3D genome variations are not associated with obvious genetic variations, and half of them affect the flanking gene expression. In addition, we find that the 3D genome may also undergo selection during soybean domestication.

CONCLUSION

Our study sheds light on the role of 3D genomes in plant genetic diversity and provides a valuable resource for studying gene regulation and genome evolution.

Emerging mechanistic insights into the regulation of specialized metabolism in plants

Plants biosynthesize a broad range of natural products through specialized and species-specific metabolic pathways that are fuelled by core metabolism, together forming a metabolic network. Specialized metabolites have important roles in development and adaptation to external cues, and they also have invaluable pharmacological properties. A growing body of evidence has highlighted the impact of translational, transcriptional, epigenetic and chromatin-based regulation and evolution of specialized metabolism genes and metabolic networks. Here we review the forefront of this research field and extrapolate to medicinal plants that synthetize rare molecules. We also discuss how this new knowledge could help in improving strategies to produce useful plant-derived pharmaceuticals.