Genome-wide analysis of local chromatin packing in Arabidopsis thaliana

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.

Stress Conditions Modulate the Chromatin Interactions Network in Arabidopsis

Stresses have been known to cause various responses like cellular physiology, gene regulation, and genome remodeling in the organism to cope and survive. Here, we assessed the impact of stress conditions on the chromatin-interactome network of Arabidopsis thaliana. We identified thousands of chromatin interactions in native as well as in salicylic acid treatment and high temperature conditions in a genome-wide fashion. Our analysis revealed the definite pattern of chromatin interactions and stress conditions could modulate the dynamics of chromatin interactions. We found the heterochromatic region of the genome actively involved in the chromatin interactions. We further observed that the establishment or loss of interactions in response to stress does not result in the global change in the expression profile of interacting genes; however, interacting regions (genes) containing motifs for known TFs showed either lower expression or no difference than non-interacting genes. The present study also revealed that interactions preferred among the same epigenetic state (ES) suggest interactions clustered the same ES together in the 3D space of the nucleus. Our analysis showed that stress conditions affect the dynamics of chromatin interactions among the chromatin loci and these interaction networks govern the folding principle of chromatin by bringing together similar epigenetic marks.

Genome assembly of alfalfa cultivar zhongmu-4 and identification of SNPs associated with agronomic traits

Alfalfa (Medicago sativa L.) is the most important legume forage crop worldwide with high nutritional value and yield. For a long time, the breeding of alfalfa was hampered by lacking reliable information on the autotetraploid genome and molecular markers linked to important agronomic traits. We herein reported the de novo assembly of the allele-aware chromosome-level genome of Zhongmu-4, a cultivar widely cultivated in China, and a comprehensive database of genomic variations based on resequencing of 220 germplasms. Approximate 2.74 Gb contigs (N50 of 2.06 Mb), accounting for 88.39% of the estimated genome, were assembled, and 2.56 Gb contigs were anchored to 32 pseudo-chromosomes. A total of 34,922 allelic genes were identified from the allele-aware genome. We observed the expansion of gene families, especially those related to the nitrogen metabolism, and the increase of repetitive elements including transposable elements, which probably resulted in the increase of Zhongmu-4 genome compared with Medicago truncatula. Population structure analysis revealed that the accessions from Asia and South America had relatively lower genetic diversity than those from Europe, suggesting that geography may influence alfalfa genetic divergence during local adaption. Genome-wide association studies identified 101 single nucleotide polymorphisms (SNPs) associated with 27 agronomic traits. Two candidate genes were predicted to be correlated with fall dormancy and salt response. We believe that the allele-aware chromosome-level genome sequence of Zhongmu-4 combined with the resequencing data of the diverse alfalfa germplasms will facilitate genetic research and genomics-assisted breeding in variety improvement of alfalfa.

Plant In Situ Hi-C Experimental Protocol and Bioinformatic Analysis

Hi-C enables the characterization of the 0conformation of the genome in the three-dimensional nuclear space. This technique has revolutionized our ability to detect interactions between linearly distant genomic sites on a genome-wide scale. Here, we detail a protocol to carry out in situ Hi-C in plants and describe a straightforward bioinformatics pipeline for the analysis of such data, in particular for comparing samples from different organs or conditions.

Interrogating Global Chromatin Interaction Network by High-Throughput Chromosome Conformation Capture (Hi-C) in Plants

High-throughput chromosome conformation capture (Hi-C) enables the global quantification of chromatin interaction frequency in eukaryotic nuclei. This method is based on in situ Hi-C, in which chromatin is cross-linked with formaldehyde, then digested with restriction enzyme. Biotin-labeled nucleotide is incorporated before the spatially adjacent DNA ends are ligated, making it possible to enrich specifically the chimeric ligation products for deep sequencing. In this chapter, we describe a modified in situ Hi-C protocol for the global chromatin interaction analysis in plants.

HiCORE: Hi-C Analysis for Identification of Core Chromatin Looping Regions with Higher Resolution

Genome-wide chromosome conformation capture (3C)- based high-throughput sequencing (Hi-C) has enabled identification of genome-wide chromatin loops. Because the Hi-C map with restriction fragment resolution is intrinsically associated with sparsity and stochastic noise, Hi-C data are usually binned at particular intervals; however, the binning method has limited reliability, especially at high resolution. Here, we describe a new method called HiCORE, which provides simple pipelines and algorithms to overcome the limitations of single-layered binning and predict core chromatin regions with three-dimensional physical interactions. In this approach, multiple layers of binning with slightly shifted genome coverage are generated, and interacting bins at each layer are integrated to infer narrower regions of chromatin interactions. HiCORE predicts chromatin looping regions with higher resolution, both in human and Arabidopsis genomes, and contributes to the identification of the precise positions of potential genomic elements in an unbiased manner.

Understanding 3D Genome Organization and Its Effect on Transcriptional Gene Regulation Under Environmental Stress in Plant: A Chromatin Perspective

The genome of a eukaryotic organism is comprised of a supra-molecular complex of chromatin fibers and intricately folded three-dimensional (3D) structures. Chromosomal interactions and topological changes in response to the developmental and/or environmental stimuli affect gene expression. Chromatin architecture plays important roles in DNA replication, gene expression, and genome integrity. Higher-order chromatin organizations like chromosome territories (CTs), A/B compartments, topologically associating domains (TADs), and chromatin loops vary among cells, tissues, and species depending on the developmental stage and/or environmental conditions (4D genomics). Every chromosome occupies a separate territory in the interphase nucleus and forms the top layer of hierarchical structure (CTs) in most of the eukaryotes. While the A and B compartments are associated with active (euchromatic) and inactive (heterochromatic) chromatin, respectively, having well-defined genomic/epigenomic features, TADs are the structural units of chromatin. Chromatin architecture like TADs as well as the local interactions between promoter and regulatory elements correlates with the chromatin activity, which alters during environmental stresses due to relocalization of the architectural proteins. Moreover, chromatin looping brings the gene and regulatory elements in close proximity for interactions. The intricate relationship between nucleotide sequence and chromatin architecture requires a more comprehensive understanding to unravel the genome organization and genetic plasticity. During the last decade, advances in chromatin conformation capture techniques for unravelling 3D genome organizations have improved our understanding of genome biology. However, the recent advances, such as Hi-C and ChIA-PET, have substantially increased the resolution, throughput as well our interest in analysing genome organizations. The present review provides an overview of the historical and contemporary perspectives of chromosome conformation capture technologies, their applications in functional genomics, and the constraints in predicting 3D genome organization. We also discuss the future perspectives of understanding high-order chromatin organizations in deciphering transcriptional regulation of gene expression under environmental stress (4D genomics). These might help design the climate-smart crop to meet the ever-growing demands of food, feed, and fodder.

The Chinese pine genome and methylome unveil key features of conifer evolution

Conifers dominate the world's forest ecosystems and are the most widely planted tree species. Their giant and complex genomes present great challenges for assembling a complete reference genome for evolutionary and genomic studies. We present a 25.4-Gb chromosome-level assembly of Chinese pine (Pinus tabuliformis) and revealed that its genome size is mostly attributable to huge intergenic regions and long introns with high transposable element (TE) content. Large genes with long introns exhibited higher expressions levels. Despite a lack of recent whole-genome duplication, 91.2% of genes were duplicated through dispersed duplication, and expanded gene families are mainly related to stress responses, which may underpin conifers' adaptation, particularly in cold and/or arid conditions. The reproductive regulation network is distinct compared with angiosperms. Slow removal of TEs with high-level methylation may have contributed to genomic expansion. This study provides insights into conifer evolution and resources for advancing research on conifer adaptation and development.

Plant 3D Chromatin Organization: Important Insights from Chromosome Conformation Capture Analyses of the Last 10 Years

Over the past few decades, eukaryotic linear genomes and epigenomes have been widely and extensively studied for understanding gene expression regulation. More recently, the three-dimensional (3D) chromatin organization was found to be important for determining genome functionality, finely tuning physiological processes for appropriate cellular responses. With the development of visualization techniques and chromatin conformation capture (3C)-based techniques, increasing evidence indicates that chromosomal architecture characteristics and chromatin domains with different epigenetic modifications in the nucleus are correlated with transcriptional activities. Subsequent studies have further explored the intricate interplay between 3D genome organization and the function of interacting regions. In this review, we summarize spatial distribution patterns of chromatin, including chromatin positioning, configurations and domains, with a particular focus on the effect of a unique form of interaction between varieties of factors that shape the 3D genome conformation in plants. We further discuss the methods, advantages and limitations of various 3C-based techniques, highlighting the applications of these technologies in plants to identify chromatin domains, and address their dynamic changes and functional implications in evolution, and adaptation to development and changing environmental conditions. Moreover, the future implications and emerging research directions of 3D genome organization are discussed.

Centromere-Specific Retrotransposons and Very-Long-Chain Fatty Acid Biosynthesis in the Genome of Yellowhorn (Xanthoceras sorbifolium, Sapindaceae), an Oil-Producing Tree With Significant Drought Resistance

In-depth genome characterization is still lacking for most of biofuel crops, especially for centromeres, which play a fundamental role during nuclear division and in the maintenance of genome stability. This study applied long-read sequencing technologies to assemble a highly contiguous genome for yellowhorn (Xanthoceras sorbifolium), an oil-producing tree, and conducted extensive comparative analyses to understand centromere structure and evolution, and fatty acid biosynthesis. We produced a reference-level genome of yellowhorn, ∼470 Mb in length with ∼95% of contigs anchored onto 15 chromosomes. Genome annotation identified 22,049 protein-coding genes and 65.7% of the genome sequence as repetitive elements. Long terminal repeat retrotransposons (LTR-RTs) account for ∼30% of the yellowhorn genome, which is maintained by a moderate birth rate and a low removal rate. We identified the centromeric regions on each chromosome and found enrichment of centromere-specific retrotransposons of LINE1 and Gypsy in these regions, which have evolved recently (∼0.7 MYA). We compared the genomes of three cultivars and found frequent inversions. We analyzed the transcriptomes from different tissues and identified the candidate genes involved in very-long-chain fatty acid biosynthesis and their expression profiles. Collinear block analysis showed that yellowhorn shared the gamma (γ) hexaploidy event with Vitis vinifera but did not undergo any further whole-genome duplication. This study provides excellent genomic resources for understanding centromere structure and evolution and for functional studies in this important oil-producing plant.

Chromatin spatial organization of wild type and mutant peanuts reveals high-resolution genomic architecture and interaction alterations

BACKGROUND

Three-dimensional (3D) chromatin organization provides a critical foundation to investigate gene expression regulation and cellular homeostasis.

RESULTS

Here, we present the first 3D genome architecture maps in wild type and mutant allotetraploid peanut lines, which illustrate A/B compartments, topologically associated domains (TADs), and widespread chromatin interactions. Most peanut chromosomal arms (52.3%) have active regions (A compartments) with relatively high gene density and high transcriptional levels. About 2.0% of chromosomal regions switch from inactive to active (B-to-A) in the mutant line, harboring 58 differentially expressed genes enriched in flavonoid biosynthesis and circadian rhythm functions. The mutant peanut line shows a higher number of genome-wide cis-interactions than its wild-type. The present study reveals a new TAD in the mutant line that generates different chromatin loops and harbors a specific upstream AP2EREBP-binding motif which might upregulate the expression of the GA2ox gene and decrease active gibberellin (GA) content, presumably making the mutant plant dwarf.

CONCLUSIONS

Our findings will shed new light on the relationship between 3D chromatin architecture and transcriptional regulation in plants.

The 3D Genome: From Structure to Function

The genome is the most functional part of a cell, and genomic contents are organized in a compact three-dimensional (3D) structure. The genome contains millions of nucleotide bases organized in its proper frame. Rapid development in genome sequencing and advanced microscopy techniques have enabled us to understand the 3D spatial organization of the genome. Chromosome capture methods using a ligation approach and the visualization tool of a 3D genome browser have facilitated detailed exploration of the genome. Topologically associated domains (TADs), lamin-associated domains, CCCTC-binding factor domains, cohesin, and chromatin structures are the prominent identified components that encode the 3D structure of the genome. Although TADs are the major contributors to 3D genome organization, they are absent in Arabidopsis. However, a few research groups have reported the presence of TAD-like structures in the plant kingdom.

The long noncoding RNA MdLNC499 bridges MdWRKY1 and MdERF109 function to regulate early-stage light-induced anthocyanin accumulation in apple fruit

Anthocyanin pigments contribute to plant coloration and are valuable sources of antioxidants in the human diet as components of fruits and vegetables. Their production is known to be induced by light in apple fruit (Malus domestica); however, the underlying molecular mechanism responsible for early-stage light-induced anthocyanin biosynthesis remains unclear. Here, we identified an ethylene response factor (ERF) protein, ERF109, involved in light-induced anthocyanin biosynthesis and found that it promotes coloration by directly binding to anthocyanin-related gene promoters. Promoter::β-glucuronidase reporter analysis and Hi-C sequencing showed that a long noncoding RNA, MdLNC499, located nearby MdERF109, induces the expression of MdERF109. A W-box cis-element in the MdLNC499 promoter was found to be regulated by a transcription factor, MdWRKY1. Transient expression in apple fruit and stable transformation of apple calli allowed us to reconstruct a MdWRKY1-MdLNC499-MdERF109 transcriptional cascade in which MdWRKY1 is activated by light to increase the transcription of MdLNC499, which in turn induces MdERF109. The MdERF109 protein induces the expression of anthocyanin-related genes and the accumulation of anthocyanins in the early stages of apple coloration. Our results provide a platform for better understanding the various regulatory mechanisms involved in light-induced apple fruit coloration.

Global profiling of RNA-chromatin interactions reveals co-regulatory gene expression networks in Arabidopsis

It is increasingly evident that various RNAs can bind chromatin to regulate gene expression and genome organization. Here we adapted a sequencing-based technique to profile RNA-chromatin interactions at a genome-wide scale in Arabidopsis seedlings. We identified more than 10,000 RNA-chromatin interactions mediated by protein-coding RNAs and non-coding RNAs. Cis and intra-chromosomal interactions are mainly mediated by protein-coding RNAs, whereas inter-chromosomal interactions are primarily mediated by non-coding RNAs. Many RNA-chromatin interactions tend to positively correlate with DNA-DNA interactions, suggesting their mutual influence and reinforcement. We further show that some RNA-chromatin interactions undergo alterations in response to biotic and abiotic stresses and that altered RNA-chromatin interactions form co-regulatory networks. Our study provides a global view on RNA-chromatin interactions in Arabidopsis and a rich resource for future investigations of regulatory roles of RNAs in gene expression and genome organization.

Insulators in Plants: Progress and Open Questions

The genomes of higher eukaryotes are partitioned into topologically associated domains or TADs, and insulators (also known as boundary elements) are the key elements responsible for their formation and maintenance. Insulators were first identified and extensively studied in Drosophila as well as mammalian genomes, and have also been described in yeast and plants. In addition, many insulator proteins are known in Drosophila, and some have been investigated in mammals. However, much less is known about this important class of non-coding DNA elements in plant genomes. In this review, we take a detailed look at known plant insulators across different species and provide an overview of potential determinants of plant insulator functions, including cis-elements and boundary proteins. We also discuss methods previously used in attempts to identify plant insulators, provide a perspective on their importance for research and biotechnology, and discuss areas of potential future research.

A novel bivalent chromatin associates with rapid induction of camalexin biosynthesis genes in response to a pathogen signal in Arabidopsis

Temporal dynamics of gene expression underpin responses to internal and environmental stimuli. In eukaryotes, regulation of gene induction includes changing chromatin states at target genes and recruiting the transcriptional machinery that includes transcription factors. As one of the most potent defense compounds in Arabidopsis thaliana, camalexin can be rapidly induced by bacterial and fungal infections. Though several transcription factors controlling camalexin biosynthesis genes have been characterized, how the rapid activation of genes in this pathway upon a pathogen signal is enabled remains unknown. By combining publicly available epigenomic data with in vivo chromatin modification mapping, we found that camalexin biosynthesis genes are marked with two epigenetic modifications with opposite effects on gene expression, trimethylation of lysine 27 of histone 3 (H3K27me3) (repression) and acetylation of lysine 18 of histone 3 (H3K18ac) (activation), to form a previously uncharacterized type of bivalent chromatin. Mutants with reduced H3K27me3 or H3K18ac suggested that both modifications were required to determine the timing of gene expression and metabolite accumulation at an early stage of the stress response. Our study indicates that the H3K27me3-H3K18ac bivalent chromatin, which we name as kairostat, plays an important role in controlling the timely induction of gene expression upon stress stimuli in plants.

Open chromatin in grapevine marks candidate CREs and with other chromatin features correlates with gene expression

Vitis vinifera is an economically important crop and a useful model in which to study chromatin dynamics. In contrast to the small and relatively simple genome of Arabidopsis thaliana, grapevine contains a complex genome of 487 Mb that exhibits extensive colonization by transposable elements. We used Hi-C, ChIP-seq and ATAC-seq to measure how chromatin features correlate to the expression of 31 845 grapevine genes. ATAC-seq revealed the presence of more than 16 000 open chromatin regions, of which we characterize nearly 5000 as possible distal enhancer candidates that occur in intergenic space > 2 kb from the nearest transcription start site (TSS). A motif search identified more than 480 transcription factor (TF) binding sites in these regions, with those for TCP family proteins in greatest abundance. These open chromatin regions are typically within 15 kb from their nearest promoter, and a gene ontology analysis indicated that their nearest genes are significantly enriched for TF activity. The presence of a candidate cis-regulatory element (cCRE) > 2 kb upstream of the TSS, location in the active nuclear compartment as determined by Hi-C, and the enrichment of H3K4me3, H3K4me1 and H3K27ac at the gene are correlated with gene expression. Taken together, these results suggest that regions of intergenic open chromatin identified by ATAC-seq can be considered potential candidates for cis-regulatory regions in V. vinifera. Our findings enhance the characterization of a valuable agricultural crop, and help to clarify the understanding of unique plant biology.

Comparative Genome Analyses Highlight Transposon-Mediated Genome Expansion and the Evolutionary Architecture of 3D Genomic Folding in Cotton

Transposable element (TE) amplification has been recognized as a driving force mediating genome size expansion and evolution, but the consequences for shaping 3D genomic architecture remains largely unknown in plants. Here, we report reference-grade genome assemblies for three species of cotton ranging 3-fold in genome size, namely Gossypium rotundifolium (K2), G. arboreum (A2), and G. raimondii (D5), using Oxford Nanopore Technologies. Comparative genome analyses document the details of lineage-specific TE amplification contributing to the large genome size differences (K2, 2.44 Gb; A2, 1.62 Gb; D5, 750.19 Mb) and indicate relatively conserved gene content and synteny relationships among genomes. We found that approximately 17% of syntenic genes exhibit chromatin status change between active ("A") and inactive ("B") compartments, and TE amplification was associated with the increase of the proportion of A compartment in gene regions (∼7,000 genes) in K2 and A2 relative to D5. Only 42% of topologically associating domain (TAD) boundaries were conserved among the three genomes. Our data implicate recent amplification of TEs following the formation of lineage-specific TAD boundaries. This study sheds light on the role of transposon-mediated genome expansion in the evolution of higher-order chromatin structure in plants.

Haplotype-resolved genome of diploid ginger (Zingiber officinale) and its unique gingerol biosynthetic pathway

Ginger (Zingiber officinale), the type species of Zingiberaceae, is one of the most widespread medicinal plants and spices. Here, we report a high-quality, chromosome-scale reference genome of ginger 'Zhugen', a traditionally cultivated ginger in Southwest China used as a fresh vegetable, assembled from PacBio long reads, Illumina short reads, and high-throughput chromosome conformation capture (Hi-C) reads. The ginger genome was phased into two haplotypes, haplotype 1 (1.53 Gb with a contig N50 of 4.68 M) and haplotype 0 (1.51 Gb with a contig N50 of 5.28 M). Homologous ginger chromosomes maintained excellent gene pair collinearity. In 17,226 pairs of allelic genes, 11.9% exhibited differential expression between alleles. Based on the results of ginger genome sequencing, transcriptome analysis, and metabolomic analysis, we proposed a backbone biosynthetic pathway of gingerol analogs, which consists of 12 enzymatic gene families, PAL, C4H, 4CL, CST, C3'H, C3OMT, CCOMT, CSE, PKS, AOR, DHN, and DHT. These analyses also identified the likely transcription factor networks that regulate the synthesis of gingerol analogs. Overall, this study serves as an excellent resource for further research on ginger biology and breeding, lays a foundation for a better understanding of ginger evolution, and presents an intact biosynthetic pathway for species-specific gingerol biosynthesis.

Disruption of topologically associating domains by structural variations in tetraploid cottons

Structural variations (SVs) are recognized to have an important role in transcriptional regulation, especially in the light of resolved 3D genome structure using high-throughput chromosome conformation capture (Hi-C) technology in mammals. However, the effect of SVs on 3D genome organization in plants remains rarely understood. In this study, we identified 295,496 SVs and 5251 topologically associating domains (TADs) in two diploid and two tetraploid cottons. We observed that approximately 16% of SVs occurred in TAD boundary regions that were called boundary affecting-structural variations (BA-SVs), and had a large effect on disrupting TAD organization. Nevertheless, SVs preferred occurring in TAD interior instead of TAD boundary, probably associated with the relaxed evolutionary selection pressure. We noticed the biased evolution of the At and Dt subgenomes of tetraploid cottons, in terms of SV-mediated disruption of 3D genome structure relative to diploids. In addition, we provide evidence showing that both SVs and TAD disruption could lead to expression difference of orthologous genes. This study advances our understanding of the effect of SVs on 3D genome organization and gene expression regulation in plants.

Three-dimensional chromosome organization in flowering plants

Research on plant three-dimensional (3D) genome architecture made rapid progress over the past 5 years. Numerous Hi-C interaction data sets were generated in a wide range of plant species, allowing for a comprehensive overview on 3D chromosome folding principles in the plant kingdom. Plants lack important genes reported to be vital for chromosome folding in animals. However, similar 3D structures such as topologically associating domains and chromatin loops were identified. Recent studies in Arabidopsis thaliana revealed how chromosomal regions are positioned within the nucleus by determining their association with both, the nuclear periphery and the nucleolus. Additionally, many plant species exhibit high-frequency interactions among KNOT entangled elements, which are associated with safeguarding the genome from invasive DNA elements. Many of the recently published Hi-C data sets were generated to aid de novo genome assembly and remain to date little explored. These data sets represent a valuable resource for future comparative studies, which may lead to a more profound understanding of the evolution of 3D chromosome organization in plants.

Altered chromatin architecture and gene expression during polyploidization and domestication of soybean

Polyploidy or whole-genome duplication (WGD) is widespread in plants and is a key driver of evolution and speciation, accompanied by rapid and dynamic changes in genomic structure and gene expression. The 3D structure of the genome is intricately linked to gene expression, but its role in transcription regulation following polyploidy and domestication remains unclear. Here, we generated high-resolution (∼2 kb) Hi-C maps for cultivated soybean (Glycine max), wild soybean (Glycine soja), and common bean (Phaseolus vulgaris). We found polyploidization in soybean may induce architecture changes of topologically associating domains and subsequent diploidization led to chromatin topology alteration around chromosome-rearrangement sites. Compared with single-copy and small-scale duplicated genes, WGD genes displayed more long-range chromosomal interactions and were coupled with higher levels of gene expression and chromatin accessibilities but void of DNA methylation. Interestingly, chromatin loop reorganization was involved in expression divergence of the genes during soybean domestication. Genes with chromatin loops were under stronger artificial selection than genes without loops. These findings provide insights into the roles of dynamic chromatin structures on gene expression during polyploidization, diploidization, and domestication of soybean.

Polycomb-dependent differential chromatin compartmentalization determines gene coregulation in Arabidopsis

In animals, distant H3K27me3-marked Polycomb targets can establish physical interactions forming repressive chromatin hubs. In plants, growing evidence suggests that H3K27me3 acts directly or indirectly to regulate chromatin interactions, although how this histone modification modulates 3D chromatin architecture remains elusive. To decipher the impact of the dynamic deposition of H3K27me3 on the Arabidopsis thaliana nuclear interactome, we combined genetics, transcriptomics, and several 3D epigenomic approaches. By analyzing mutants defective for histone H3K27 methylation or demethylation, we uncovered the crucial role of this chromatin mark in short- and previously unnoticed long-range chromatin loop formation. We found that a reduction in H3K27me3 levels led to a decrease in the interactions within Polycomb-associated repressive domains. Regions with lower H3K27me3 levels in the H3K27 methyltransferase clf mutant established new interactions with regions marked with H3K9ac, a histone modification associated with active transcription, indicating that a reduction in H3K27me3 levels induces a global reconfiguration of chromatin architecture. Altogether, our results reveal that the 3D genome organization is tightly linked to reversible histone modifications that govern chromatin interactions. Consequently, nuclear organization dynamics shapes the transcriptional reprogramming during plant development and places H3K27me3 as a key feature in the coregulation of distant genes.

The nearly complete genome of Ginkgo biloba illuminates gymnosperm evolution

Gymnosperms are a unique lineage of plants that currently lack a high-quality reference genome due to their large genome size and high repetitive sequence content. Here, we report a nearly complete genome assembly for Ginkgo biloba with a genome size of 9.87 Gb, an N50 contig size of 1.58 Mb and an N50 scaffold size of 775 Mb. We were able to accurately annotate 27,832 protein-coding genes in total, superseding the inaccurate annotation of 41,840 genes in a previous draft genome assembly. We found that expansion of the G. biloba genome, accompanied by the notable extension of introns, was mainly caused by the insertion of long terminal repeats rather than the recent occurrence of whole-genome duplication events, in contrast to the findings of a previous report. We also identified candidate genes in the central pair, intraflagellar transport and dynein protein families that are associated with the formation of the spermatophore flagellum, which has been lost in all seed plants except ginkgo and cycads. The newly obtained Ginkgo genome provides new insights into the evolution of the gymnosperm genome.

Plant 3D genomics: the exploration and application of chromatin organization

Eukaryotic genomes are highly folded for packing into higher-order chromatin structures in the nucleus. With the emergence of state-of-the-art chromosome conformation capture methods and microscopic imaging techniques, the spatial organization of chromatin and its functional implications have been interrogated. Our knowledge of 3D chromatin organization in plants has improved dramatically in the past few years, building on the early advances in animal systems. Here, we review recent advances in 3D genome mapping approaches, our understanding of the sophisticated organization of spatial structures, and the application of 3D genomic principles in plants. We also discuss directions for future developments in 3D genomics in plants.

Current progress and challenges in crop genetic transformation

Plant transformation remains the most sought-after technology for functional genomics and crop genetic improvement, especially for introducing specific new traits and to modify or recombine already existing traits. Along with many other agricultural technologies, the global production of genetically engineered crops has steadily grown since they were first introduced 25 years ago. Since the first transfer of DNA into plant cells using Agrobacterium tumefaciens, different transformation methods have enabled rapid advances in molecular breeding approaches to bring crop varieties with novel traits to the market that would be difficult or not possible to achieve with conventional breeding methods. Today, transformation to produce genetically engineered crops is the fastest and most widely adopted technology in agriculture. The rapidly increasing number of sequenced plant genomes and information from functional genomics data to understand gene function, together with novel gene cloning and tissue culture methods, is further accelerating crop improvement and trait development. These advances are welcome and needed to make crops more resilient to climate change and to secure their yield for feeding the increasing human population. Despite the success, transformation remains a bottleneck because many plant species and crop genotypes are recalcitrant to established tissue culture and regeneration conditions, or they show poor transformability. Improvements are possible using morphogenetic transcriptional regulators, but their broader applicability remains to be tested. Advances in genome editing techniques and direct, non-tissue culture-based transformation methods offer alternative approaches to enhance varietal development in other recalcitrant crops. Here, we review recent developments in plant transformation and regeneration, and discuss opportunities for new breeding technologies in agriculture.

Chromatin architectural proteins regulate flowering time by precluding gene looping

Chromatin structure is critical for gene expression and many other cellular processes. In Arabidopsis thaliana, the floral repressor FLC adopts a self-loop chromatin structure via bridging of its flanking regions. This local gene loop is necessary for active FLC expression. However, the molecular mechanism underlying the formation of this class of gene loops is unknown. Here, we report the characterization of a group of linker histone-like proteins, named the GH1-HMGA family in Arabidopsis, which act as chromatin architecture modulators. We demonstrate that these family members redundantly promote the floral transition through the repression of FLC A genome-wide study revealed that this family preferentially binds to the 5' and 3' ends of gene bodies. The loss of this binding increases FLC expression by stabilizing the FLC 5' to 3' gene looping. Our study provides mechanistic insights into how a family of evolutionarily conserved proteins regulates the formation of local gene loops.

INT-Hi-C reveals distinct chromatin architecture in endosperm and leaf tissues of Arabidopsis

Higher-order chromatin structure undergoes striking changes in response to various developmental and environmental signals, causing distinct cell types to adopt specific chromatin organization. High throughput chromatin conformation capture (Hi-C) allows studying higher-order chromatin structure; however, this technique requires substantial amounts of starting material, which has limited the establishment of cell type-specific higher-order chromatin structure in plants. To overcome this limitation, we established a protocol that is applicable to a limited amount of nuclei by combining the INTACT (isolation of nuclei tagged in specific cell types) method and Hi-C (INT-Hi-C). Using this INT-Hi-C protocol, we generated Hi-C data from INTACT purified endosperm and leaf nuclei. Our INT-Hi-C data from leaf accurately reiterated chromatin interaction patterns derived from conventional leaf Hi-C data. We found that the higher-order chromatin organization of mixed leaf tissues and endosperm differs and that DNA methylation and repressive histone marks positively correlate with the chromatin compaction level. We furthermore found that self-looped interacting genes have increased expression in leaves and endosperm and that interacting intergenic regions negatively impact on gene expression in the endosperm. Last, we identified several imprinted genes involved in long-range and trans interactions exclusively in endosperm. Our study provides evidence that the endosperm adopts a distinct higher-order chromatin structure that differs from other cell types in plants and that chromatin interactions influence transcriptional activity.

Comparative genome analyses highlight transposon-mediated genome expansion and the evolutionary architecture of 3D genomic folding in cotton

Transposable element (TE) amplification has been recognized as a driving force mediating genome size expansion and evolution, but the consequences for shaping 3D genomic architecture remains largely unknown in plants. Here, we report reference-grade genome assemblies for three species of cotton ranging 3-fold in genome size, namely Gossypium rotundifolium (K₂), G. arboreum (A₂), and G. raimondii (D₅), using Oxford Nanopore Technologies. Comparative genome analyses document the details of lineage-specific TE amplification contributing to the large genome size differences (K₂, 2.44 Gb; A₂, 1.62 Gb; D₅, 750.19 Mb) and indicate relatively conserved gene content and synteny relationships among genomes. We found that approximately 17% of syntenic genes exhibit chromatin status change between active (“A”) and inactive (“B”) compartments, and TE amplification was associated with the increase of the proportion of A compartment in gene regions (∼7,000 genes) in K₂ and A₂ relative to D₅. Only 42% of topologically associating domain (TAD) boundaries were conserved among the three genomes. Our data implicate recent amplification of TEs following the formation of lineage-specific TAD boundaries. This study sheds light on the role of transposon-mediated genome expansion in the evolution of higher-order chromatin structure in plants.

Subtelomeric assembly of a multi-gene pathway for antimicrobial defense compounds in cereals

Non-random gene organization in eukaryotes plays a significant role in genome evolution. Here, we investigate the origin of a biosynthetic gene cluster for production of defence compounds in oat-the avenacin cluster. We elucidate the structure and organisation of this 12-gene cluster, characterise the last two missing pathway steps, and reconstitute the entire pathway in tobacco by transient expression. We show that the cluster has formed de novo since the divergence of oats in a subtelomeric region of the genome that lacks homology with other grasses, and that gene order is approximately colinear with the biosynthetic pathway. We speculate that the positioning of the late pathway genes furthest away from the telomere may mitigate against a 'self-poisoning' scenario in which toxic intermediates accumulate as a result of telomeric gene deletions. Our investigations reveal a striking example of adaptive evolution underpinned by remarkable genome plasticity.

Large-scale reconstruction of chromatin structures of maize temperate and tropical inbred lines

Maize is a model plant species often used for genetics and genomics research because of its genetic diversity. There are prominent morphological, genetic, and epigenetic variations between tropical and temperate maize lines. However, the genome-wide chromatin conformations of these two maize types remain unexplored. We applied a Hi-C approach to compare the genome-wide chromatin interactions between temperate inbred line D132 and tropical line CML288. A reconstructed maize three-dimensional genome model revealed the spatial segregation of the global A and B compartments. The A compartments contain enriched genes and active epigenome marks, whereas the B compartments are gene-poor, transcriptionally silent chromatin regions. Whole-genome analyses indicated that the global A compartment content of CML288 was 3.12% lower than that of D132. Additionally, global and A/B sub-compartments were associated with differential gene expression and epigenetic changes between two inbred lines. About 25.3% of topologically associating domains (TADs) were determined to be associated with complex domain-level modifications that induced transcriptional changes, indicative of a large-scale reorganization of chromatin structures between the inbred maize lines. Furthermore, differences in chromatin interactions between the two lines correlated with epigenetic changes. These findings provide a solid foundation for the wider plant community to further investigate the genome-wide chromatin structures in other plant species.

SPADIS: An Algorithm for Selecting Predictive and Diverse SNPs in GWAS

Phenotypic heritability of complex traits and diseases is seldom explained by individual genetic variants identified in genome-wide association studies (GWAS). Many methods have been developed to select a subset of variant loci, which are associated with or predictive of the phenotype. Selecting connected SNPs on SNP-SNP networks have been proven successful in finding biologically interpretable and predictive SNPs. However, we argue that the connectedness constraint favors selecting redundant features that affect similar biological processes and therefore does not necessarily yield better predictive performance. In this paper, we propose a novel method called SPADIS that favors the selection of remotely located SNPs in order to account for their complementary effects in explaining a phenotype. SPADIS selects a diverse set of loci on a SNP-SNP network. This is achieved by maximizing a submodular set function with a greedy algorithm that ensures a constant factor approximation to the optimal solution. We compare SPADIS to the state-of-the-art method SConES, on a dataset of Arabidopsis Thaliana with continuous flowering time phenotypes. SPADIS has better average phenotype prediction performance in 15 out of 17 phenotypes when the same number of SNPs are selected and provides consistent improvements across multiple networks and settings on average. Moreover, it identifies more candidate genes and runs faster.

Correlated Evolution of Large DNA Fragments in the 3D Genome of Arabidopsis thaliana

In eukaryotes, the three-dimensional (3D) conformation of the genome is far from random, and this nonrandom chromatin organization is strongly correlated with gene expression and protein function, which are two critical determinants of the selective constraints and evolutionary rates of genes. However, whether genes and other elements that are located close to each other in the 3D genome evolve in a coordinated way has not been investigated in any organism. To address this question, we constructed chromatin interaction networks (CINs) in Arabidopsis thaliana based on high-throughput chromosome conformation capture data and demonstrated that adjacent large DNA fragments in the CIN indeed exhibit more similar levels of polymorphism and evolutionary rates than random fragment pairs. Using simulations that account for the linear distance between fragments, we proved that the 3D chromosomal organization plays a role in the observed correlated evolution. Spatially interacting fragments also exhibit more similar mutation rates and functional constraints in both coding and noncoding regions than the random expectations, indicating that the correlated evolution between 3D neighbors is a result of combined evolutionary forces. A collection of 39 genomic and epigenomic features can explain much of the variance in genetic diversity and evolutionary rates across the genome. Moreover, features that have a greater effect on the evolution of regional sequences tend to show higher similarity between neighboring fragments in the CIN, suggesting a pivotal role of epigenetic modifications and chromatin organization in determining the correlated evolution of large DNA fragments in the 3D genome.

KINESIN-12E regulates metaphase spindle flux and helps control spindle size in Arabidopsis

The bipolar mitotic spindle is a highly conserved structure among eukaryotes that mediates chromosome alignment and segregation. Spindle assembly and size control are facilitated by force-generating microtubule-dependent motor proteins known as kinesins. In animals, kinesin-12 cooperates with kinesin-5 to produce outward-directed forces necessary for spindle assembly. In plants, the relevant molecular mechanisms for spindle formation are poorly defined. While an Arabidopsis thaliana kinesin-5 ortholog has been identified, the kinesin-12 ortholog in plants remains elusive. In this study, we provide experimental evidence for the function of Arabidopsis KINESIN-12E in spindle assembly. In kinesin-12e mutants, a delay in spindle assembly is accompanied by the reduction of spindle size, demonstrating that KINESIN-12E contributes to mitotic spindle architecture. Kinesin-12E localization is mitosis-stage specific, beginning with its perinuclear accumulation during prophase. Upon nuclear envelope breakdown, KINESIN-12E decorates subpopulations of microtubules in the spindle and becomes progressively enriched in the spindle midzone. Furthermore, during cytokinesis, KINESIN-12E shares its localization at the phragmoplast midzone with several functionally diversified Arabidopsis KINESIN-12 members. Changes in the kinetochore and in prophase and metaphase spindle dynamics occur in the absence of KINESIN-12E, suggest it might play an evolutionarily conserved role during spindle formation similar to its spindle-localized animal kinesin-12 orthologs.

Polyploidy: an evolutionary and ecological force in stressful times

Polyploidy has been hypothesized to be both an evolutionary dead-end and a source for evolutionary innovation and species diversification. Although polyploid organisms, especially plants, abound, the apparent nonrandom long-term establishment of genome duplications suggests a link with environmental conditions. Whole-genome duplications seem to correlate with periods of extinction or global change, while polyploids often thrive in harsh or disturbed environments. Evidence is also accumulating that biotic interactions, for instance, with pathogens or mutualists, affect polyploids differently than nonpolyploids. Here, we review recent findings and insights on the effect of both abiotic and biotic stress on polyploids versus nonpolyploids and propose that stress response in general is an important and even determining factor in the establishment and success of polyploidy.

Bioinformatic analysis of chromatin organization and biased expression of duplicated genes between two poplars with a common whole-genome duplication

The nonrandom three-dimensional organization of chromatin plays an important role in the regulation of gene expression. However, it remains unclear whether this organization is conserved and whether it is involved in regulating gene expression during speciation after whole-genome duplication (WGD) in plants. In this study, high-resolution interaction maps were generated using high-throughput chromatin conformation capture (Hi-C) techniques for two poplar species, Populus euphratica and Populus alba var. pyramidalis, which diverged ~14 Mya after a common WGD. We examined the similarities and differences in the hierarchical chromatin organization between the two species, including A/B compartment regions and topologically associating domains (TADs), as well as in their DNA methylation and gene expression patterns. We found that chromatin status was strongly associated with epigenetic modifications and gene transcriptional activity, yet the conservation of hierarchical chromatin organization across the two species was low. The divergence of gene expression between WGD-derived paralogs was associated with the strength of chromatin interactions, and colocalized paralogs exhibited strong similarities in epigenetic modifications and expression levels. Thus, the spatial localization of duplicated genes is highly correlated with biased expression during the diploidization process. This study provides new insights into the evolution of chromatin organization and transcriptional regulation during the speciation process of poplars after WGD.

Polymer modelling unveils the roles of heterochromatin and nucleolar organizing regions in shaping 3D genome organization in Arabidopsis thaliana

The 3D genome is characterized by a complex organization made of genomic and epigenomic layers with profound implications on gene regulation and cell function. However, the understanding of the fundamental mechanisms driving the crosstalk between nuclear architecture and (epi)genomic information is still lacking. The plant Arabidopsis thaliana is a powerful model organism to address these questions owing to its compact genome for which we have a rich collection of microscopy, chromosome conformation capture (Hi-C) and ChIP-seq experiments. Using polymer modelling, we investigate the roles of nucleolus formation and epigenomics-driven interactions in shaping the 3D genome of A. thaliana. By validation of several predictions with published data, we demonstrate that self-attracting nucleolar organizing regions and repulsive constitutive heterochromatin are major mechanisms to regulate the organization of chromosomes. Simulations also suggest that interphase chromosomes maintain a partial structural memory of the V-shapes, typical of (sub)metacentric chromosomes in anaphase. Additionally, self-attraction between facultative heterochromatin regions facilitates the formation of Polycomb bodies hosting H3K27me3-enriched gene-clusters. Since nucleolus and heterochromatin are highly-conserved in eukaryotic cells, our findings pave the way for a comprehensive characterization of the generic principles that are likely to shape and regulate the 3D genome in many species.

Homology-mediated inter-chromosomal interactions in hexaploid wheat lead to specific subgenome territories following polyploidization and introgression

BACKGROUND

Polyploidization and introgression are major events driving plant genome evolution and influencing crop breeding. However, the mechanisms underlying the higher-order chromatin organization of subgenomes and alien chromosomes are largely unknown.

RESULTS

We probe the three-dimensional chromatin architecture of Aikang 58 (AK58), a widely cultivated allohexaploid wheat variety in China carrying the 1RS/1BL translocation chromosome. The regions involved in inter-chromosomal interactions, both within and between subgenomes, have highly similar sequences. Subgenome-specific territories tend to be connected by subgenome-dominant homologous transposable elements (TEs). The alien 1RS chromosomal arm, which was introgressed from rye and differs from its wheat counterpart, has relatively few inter-chromosome interactions with wheat chromosomes. An analysis of local chromatin structures reveals topologically associating domain (TAD)-like regions covering 52% of the AK58 genome, the boundaries of which are enriched with active genes, zinc-finger factor-binding motifs, CHH methylation, and 24-nt small RNAs. The chromatin loops are mostly localized around TAD boundaries, and the number of gene loops is positively associated with gene activity.

CONCLUSIONS

The present study reveals the impact of the genetic sequence context on the higher-order chromatin structure and subgenome stability in hexaploid wheat. Specifically, we characterized the sequence homology-mediated inter-chromosome interactions and the non-canonical role of subgenome-biased TEs. Our findings may have profound implications for future investigations of the interplay between genetic sequences and higher-order structures and their consequences on polyploid genome evolution and introgression-based breeding of crop plants.

Unraveling the 3D Genome Architecture in Plants: Present and Future

The eukaryotic genome has a hierarchical three-dimensional (3D) organization with functional implications for DNA replication, DNA repair, and transcriptional regulation. Over the past decade, scientists have endeavored to elucidate the spatial characteristics and functions of plant genome architecture using high-throughput chromatin conformation capturing technologies such as Hi-C, ChIA-PET, and HiChIP. Here, we systematically review current understanding of chromatin organization in plants at multiple scales. We also discuss the emerging opinions and concepts in 3D genome research, focusing on state-of-the-art 3D genome techniques, RNA-chromatin interactions, liquid-liquid phase separation, and dynamic chromatin alterations. We propose the application of single-cell/single-molecule multi-omics, multiway (DNA-DNA, DNA-RNA, and RNA-RNA interactions) chromatin conformation capturing methods, and proximity ligation-independent 3D genome-mapping technologies to explore chromatin organization structure and function in plants. Such methods could reveal the spatial interactions between trait-related SNPs and their target genes at various spatiotemporal resolutions, and elucidate the molecular mechanisms of the interactions among DNA elements, RNA molecules, and protein factors during the formation of key traits in plants.

High-quality genome assembly of Metaphire vulgaris

Earthworms enrich the soil and protect the health of their ecological environment. Previous studies on these invertebrates determined their protein content, hormone secretions, medicinal value, and ecological habits, but their whole genomic sequence remains incomplete. We performed whole genome sequencing of Metaphire vulgaris (Chen, 1930), which belongs to the genus Metaphire of the family Megascolecidae. The genome assembly was 729 Mb, with a N50 contig size of 4.2 Mb. In total, 559 contigs were anchored to 41 chromosomes according to the results of Hi-C (High-throughput Chromosome Conformation Capture) technology, which was confirmed by karyological analysis. A comparison of the genomic sequences and genes indicated that there was a whole-genome duplication in M. vulgaris followed by several chromosome fusion events. Hox genes and lumbrokinase genes were identified as partial clusters surrounding the genome. Our high-quality genome assembly of M. vulgaris will provide valuable information for gene function and evolutionary studies in earthworms.

Comparative Genome Analysis of Scutellaria baicalensis and Scutellaria barbata Reveals the Evolution of Active Flavonoid Biosynthesis

Scutellaria baicalensis (S. baicalensis) and Scutellaria barbata (S. barbata) are common medicinal plants of the Lamiaceae family. Both produce specific flavonoid compounds, including baicalein, scutellarein, norwogonin, and wogonin, as well as their glycosides, which exhibit antioxidant and antitumor activities. Here, we report chromosome-level genome assemblies of S. baicalensis and S. barbata with quantitative chromosomal variation (2n = 18 and 2n = 26, respectively). The divergence of S. baicalensis and S. barbata occurred far earlier than previously reported, and a whole-genome duplication (WGD) event was identified. The insertion of long terminal repeat elements after speciation might be responsible for the observed chromosomal expansion and rearrangement. Comparative genome analysis of the congeneric species revealed the species-specific evolution of chrysin and apigenin biosynthetic genes, such as the S. baicalensis-specific tandem duplication of genes encoding phenylalanine ammonia lyase and chalcone synthase, and the S. barbata-specific duplication of genes encoding 4-CoA ligase. In addition, the paralogous duplication, colinearity, and expression diversity of CYP82D subfamily members revealed the functional divergence of genes encoding flavone hydroxylase between S. baicalensis and S. barbata. Analyzing these Scutellaria genomes reveals the common and species-specific evolution of flavone biosynthetic genes. Thus, these findings would facilitate the development of molecular breeding and studies of biosynthesis and regulation of bioactive compounds.

Mapping Regulatory Determinants in Plants

The domestication and improvement of many plant species have frequently involved modulation of transcriptional outputs and continue to offer much promise for targeted trait engineering. The cis-regulatory elements (CREs) controlling these trait-associated transcriptional variants however reside within non-coding regions that are currently poorly annotated in most plant species. This is particularly true in large crop genomes where regulatory regions constitute only a small fraction of the total genomic space. Furthermore, relatively little is known about how CREs function to modulate transcription in plants. Therefore understanding where regulatory regions are located within a genome, what genes they control, and how they are structured are important factors that could be used to guide both traditional and synthetic plant breeding efforts. Here, we describe classic examples of regulatory instances as well as recent advances in plant regulatory genomics. We highlight valuable molecular tools that are enabling large-scale identification of CREs and offering unprecedented insight into how genes are regulated in diverse plant species. We focus on chromatin environment, transcription factor (TF) binding, the role of transposable elements, and the association between regulatory regions and target genes.

DNA replication and chromosome positioning throughout the interphase in three-dimensional space of plant nuclei

Despite much recent progress, our understanding of the principles of plant genome organization and its dynamics in three-dimensional space of interphase nuclei remains surprisingly limited. Notably, it is not clear how these processes could be affected by the size of a plant's nuclear genome. In this study, DNA replication timing and interphase chromosome positioning were analyzed in seven Poaceae species that differ in their genome size. To provide a comprehensive picture, a suite of advanced, complementary methods was used: labeling of newly replicated DNA by ethynyl-2'-deoxyuridine, isolation of nuclei at particular cell cycle phases by flow cytometric sorting, three-dimensional immunofluorescence in situ hybridization, and confocal microscopy. Our results revealed conserved dynamics of DNA replication in all species, and a similar replication timing order for telomeres and centromeres, as well as for euchromatin and heterochromatin regions, irrespective of genome size. Moreover, stable chromosome positioning was observed while transitioning through different stages of interphase. These findings expand upon earlier studies in suggesting that a more complex interplay exists between genome size, organization of repetitive DNA sequences along chromosomes, and higher order chromatin structure and its maintenance in interphase, albeit controlled by currently unknown factors.

Epigenetic silencing of clustered tRNA genes in Arabidopsis

Beyond their key role in translation, cytosolic transfer RNAs (tRNAs) are involved in a wide range of other biological processes. Nuclear tRNA genes (tDNAs) are transcribed by the RNA polymerase III (RNAP III) and cis-elements, trans-factors as well as genomic features are known to influence their expression. In Arabidopsis, besides a predominant population of dispersed tDNAs spread along the 5 chromosomes, some clustered tDNAs have been identified. Here, we demonstrate that these tDNA clusters are transcriptionally silent and that pathways involved in the maintenance of DNA methylation play a predominant role in their repression. Moreover, we show that clustered tDNAs exhibit repressive chromatin features whilst their dispersed counterparts contain permissive euchromatic marks. This work demonstrates that both genomic and epigenomic contexts are key players in the regulation of tDNAs transcription. The conservation of most of these regulatory processes suggests that this pioneering work in Arabidopsis can provide new insights into the regulation of RNA Pol III transcription in other organisms, including vertebrates.

Decoding the plant genome: From epigenome to 3D organization

The linear genome of eukaryotes is partitioned into diverse chromatin states and packaged into a three-dimensional (3D) structure, which has functional implications in DNA replication, DNA repair, and transcriptional regulation. Over the past decades, research on plant functional genomics and epigenomics has made great progress, with thousands of genes cloned and molecular mechanisms of diverse biological processes elucidated. Recently, 3D genome research has gradually attracted great attention of many plant researchers. Herein, we briefly review the progress in genomic and epigenomic research in plants, with a focus on Arabidopsis and rice, and summarize the currently used technologies and advances in plant 3D genome organization studies. We also discuss the relationships between one-dimensional linear genome sequences, epigenomic states, and the 3D chromatin architecture. This review provides basis for future research on plant 3D genomics.

Marchantia TCP transcription factor activity correlates with three-dimensional chromatin structure

Information in the genome is not only encoded within sequence or epigenetic modifications, but is also found in how it folds in three-dimensional space. The formation of self-interacting genomic regions, named topologically associated domains (TADs), is known as a key feature of genome organization beyond the nucleosomal level. However, our understanding of the formation and function of TADs in plants is extremely limited. Here we show that the genome of Marchantia polymorpha, a member of a basal land plant lineage, exhibits TADs with epigenetic features similar to those of higher plants. By analysing various epigenetic marks across Marchantia TADs, we find that these regions generally represent interstitial heterochromatin and their borders are enriched with Marchantia transcription factor TCP1. We also identify a type of TAD that we name 'TCP1-rich TAD', in which genomic regions are highly accessible and are densely bound by TCP1 proteins. Transcription of TCP1 target genes differs on the basis gene location, and those in TCP1-rich TADs clearly show a lower expression level. In tcp1 mutant lines, neither TCP1-bound TAD borders nor TCP1-rich TADs display drastically altered chromatin organization patterns, suggesting that, in Marchantia, TCP1 is dispensable for TAD formation. However, we find that in tcp1 mutants, genes residing in TCP1-rich TADs have a greater extent of expression fold change as opposed to genes that do not belong to these TADs. Our results suggest that, besides standing as spatial chromatin-packing modules, plant TADs function as nuclear microcompartments associated with transcription factor activities.

Weak interactions in higher-order chromatin organization

The detailed principles of the hierarchical folding of eukaryotic chromosomes have been revealed during the last two decades. Along with structures composing three-dimensional (3D) genome organization (chromatin compartments, topologically associating domains, chromatin loops, etc.), the molecular mechanisms that are involved in their establishment and maintenance have been characterized. Generally, protein-protein and protein-DNA interactions underlie the spatial genome organization in eukaryotes. However, it is becoming increasingly evident that weak interactions, which exist in biological systems, also contribute to the 3D genome. Here, we provide a snapshot of our current understanding of the role of the weak interactions in the establishment and maintenance of the 3D genome organization. We discuss how weak biological forces, such as entropic forces operating in crowded solutions, electrostatic interactions of the biomolecules, liquid-liquid phase separation, DNA supercoiling, and RNA environment participate in chromosome segregation into structural and functional units and drive intranuclear functional compartmentalization.

The matrix revolutions: towards the decoding of the plant chromatin three-dimensional reality

In recent years, we have witnessed a significant increase in studies addressing the three-dimensional (3D) chromatin organization of the plant nucleus. Important advances in chromatin conformation capture (3C)-derived and related techniques have allowed the exploration of the nuclear topology of plants with large and complex genomes, including various crops. In addition, the increase in their resolution has permitted the depiction of chromatin compartmentalization and interactions at the gene scale. These studies have revealed the highly complex mechanisms governing plant nuclear architecture and the remarkable knowledge gaps in this field. Here we discuss the state-of-the-art in plant chromosome architecture, including our knowledge of the hierarchical organization of the genome in 3D space and regarding other nuclear components. Furthermore, we highlight the existence in plants of topologically associated domain (TAD)-like structures that display striking differences from their mammalian counterparts, proposing the concept of ICONS-intergenic condensed spacers. Similarly, we explore recent advances in the study of chromatin loops and R-loops, and their implication in the regulation of gene activity. Finally, we address the impact that polyploidization has had on the chromatin topology of modern crops, and how this is related to phenomena such as subgenome dominance and biased gene retention in these organisms.

Chromatin regulation of somatic abiotic stress memory

In nature, plants are often subjected to periods of recurrent environmental stress that can strongly affect their development and productivity. To cope with these conditions, plants can remember a previous stress, which allows them to respond more efficiently to a subsequent stress, a phenomenon known as priming. This ability can be maintained at the somatic level for a few days or weeks after the stress is perceived, suggesting that plants can store information of a past stress during this recovery phase. While the immediate responses to a single stress event have been extensively studied, knowledge on priming effects and how stress memory is stored is still scarce. At the molecular level, memory of a past condition often involves changes in chromatin structure and organization, which may be maintained independently from transcription. In this review, we will summarize the most recent developments in the field and discuss how different levels of chromatin regulation contribute to priming and plant abiotic stress memory.