
==== Front
Plant Cell Rep
Plant Cell Rep
Plant Cell Reports
0721-7714
1432-203X
Springer Berlin Heidelberg Berlin/Heidelberg

3322
10.1007/s00299-024-03322-8
Original Article
Connecting high-resolution 3D chromatin maps with cell division and cell differentiation at the root apical meristem
Caballero Lara
http://orcid.org/0000-0001-5061-4076
Pasternak Taras
http://orcid.org/0000-0003-2580-2811
Riyazuddin Riyazuddin
http://orcid.org/0000-0003-2848-4919
Pérez-Pérez José Manuel jmperez@umh.es

https://ror.org/01azzms13 grid.26811.3c 0000 0001 0586 4893 Instituto de Bioingeniería, Universidad Miguel Hernández, 03202 Elche, Spain
Communicated by Attila Feher.

16 9 2024
16 9 2024
2024
43 10 23229 7 2024
9 9 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Key message

We used marker-free technologies to study chromatin at cellular resolution. Our results show asymmetric chromatin distribution, explore chromatin dynamics during mitosis, and reveal structural differences between trichoblast and atrichoblast cell.

Abstract

The shapes, sizes, and structural organizations of plant nuclei vary considerably among cell types, tissues, and species. This diversity is dependent on various factors, including cellular function, developmental stage, and environmental or physiological conditions. The differences in nuclear structure reflect the state of chromatin, which, in turn, controls gene expression and regulates cell fate. To examine the interrelationship between nuclear structure, cell morphology, and tissue-specific cell proliferation and differentiation processes, we conducted multiple visualizations of H3K4me1, H3K9me2, 4′,6-diamidino-2-phenylindole, 5-ethynyl 2′-deoxyuridine, and SCRI Renaissance 2200, followed by subsequent quantitative analysis of individual cells and nuclei. By assigning cylindrical coordinates to the nuclei in the iRoCS toolbox, we were able to construct in situ digital three-dimensional chromatin maps for all the tissue layers of individual roots. A detailed analysis of the nuclei features of H3K4me1 and H3K9me2 in the mitotic and the elongation zones in trichoblast and atrichoblast cells at the root apical meristem revealed cell type-specific chromatin dynamics with asymmetric distribution of euchromatin and heterochromatin marks that may be associated with cell cycle and cell differentiation characteristics of specific cells. Furthermore, the spatial distribution of nuclei stained with 5-ethynyl 2′-deoxyuridine in the epidermis and cortex tissues suggests short-range coordination of cell division and nuclear migration in a linear sequence through an unknown regulatory mechanism.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00299-024-03322-8.

Keywords

Root apical meristem
Cell division
DNA replication
Cell differentiation
Euchromatin
Heterochromatin
Nuclear structure
http://dx.doi.org/10.13039/501100004837 Ministerio de Ciencia e Innovación PID2021-126840OB-I00 Pérez-Pérez José Manuel http://dx.doi.org/10.13039/501100003359 Generalitat Valenciana INVEST/2022/247 Caballero Lara Universidad Miguel HernándezOpen Access funding provided thanks to the CRUE-CSIC agreement with Springer Nature.

issue-copyright-statement© Springer-Verlag GmbH Germany, part of Springer Nature 2024
==== Body
pmcIntroduction

The nucleus of plant cells is a complex structure that plays a critical role not only in DNA storage but also in DNA replication, DNA repair, and regulation of gene expression in response to physiological and environmental signals (Goto et al. 2021). The principles of nuclear organization are conserved throughout evolution (Misteli 2020). The distribution of histones and histone marks, the organization of chromosomal territories, and the formation of higher-order chromosomal domains have been extensively studied using advanced confocal laser scanning, and super-resolution microscopy (Dumur et al. 2019; Amiad-Pavlov et al. 2021; Zhang et al. 2024). Several studies in different species, including the model plant Arabidopsis thaliana, have established three-dimensional (3D) models of chromatin organization (Pontvianne and Grob 2020). However, most of these studies have been performed on isolated nuclei, without taking into account cell types and positions.

Plant cells exhibit a great diversity of nuclear shapes, ploidy levels, and nuclear structural organization that varies among plant species and cell types within an organism. This diversity is influenced by cellular function, developmental stage, and physiological conditions (Carotenuto et al. 2019; Cantwell & Dey 2022). Differences in nuclear structure reflect the state of chromatin, which is known to regulate gene expression and thus cell fate (Mattout et al. 2015). Although these differences are known and have been described in previous work (Pavlova et al. 2022), they have not yet been precisely quantified at the whole organ level. Therefore, there is a growing need to better understand chromatin structure in plants in situ. Comparing nuclear structure in different plant species and at different developmental stages has the potential to reveal evolutionarily significant functional features of nuclear organization.

The INDEPTH project, which stands for Impact of Nuclear Domains on Gene Expression and Plant Traits, investigates the function of chromatin domains in nuclear processes and gene expression control (Tatout et al. 2022). Several semi-automated tools have been provided by the participating groups to characterize plant chromatin architecture, including NucleusJ (Desset et al. 2018), NucleusJ2.0 (Dubos et al. 2020), and NODeJ (Dubos et al. 2022). Standardized image acquisition and analysis procedures were developed to examine the spatial distribution of nuclear and chromosomal signals from 3D image stacks (Randall et al. 2022). Two complementary segmentation methods have been developed, one semi-automated (iCRAQ) and the other based on deep learning (Nucl.Eye.D). These methods allow the software to be trained with user-specific datasets and purposes, thus expanding the available toolbox in the field of plant cytogenetics (Johann to Berens et al. 2022). In addition, the 3D Object Counter plugin (Bolte and Cordelières 2006) allows the extraction of at least 25 parameters related to nuclear organization, including the quantification of ploidy level, chromatin abundance and distribution, nuclear shape, and nuclear asymmetry. The same plugin can be used to quantify euchromatin, heterochromatin, and DNA replication levels. In addition, the introduction of a spherical coordinate system for each individual nucleus may allow the study of their internal asymmetry in relation to its cell fate.

In our study, we applied nuclear DNA labeling in combination with whole-mount in situ localization of euchromatin and heterochromatin marks, quantitative analysis of DNA replication with 5-ethynyl 2´-deoxyuridine (EdU), and cell wall staining with SCRI Renaissance 2200 (SR2200) in each cell within the root apical meristem (RAM) of primary roots (PR) of young A. thaliana and Solanum lycopersicum seedlings. Subsequently, deep image analysis was performed to examine cellular morphology and nuclear organization in these root apices with 3D spatial resolution. Our experimental design uses conventional confocal laser scanning microscopy to examine the nuclear organization within the tissue context of whole organs. The objective of this study was to examine the dynamic regulation of euchromatin and heterochromatin in specific cell types within the proliferation domain of the RAM region. This approach recognizes the complexity and diversity of nuclear organization in plant organs and aims to capture it comprehensively. We examine here the diversity of nuclear and chromatin organization and its potential impact on nuclear and chromatin morphology, as well as DNA replication. Our results are consistent with the hypothesis that there is no universal model for nuclear organization and chromatin modification, and that variation exists between different cell types within plant organs and between species.

Materials and methods

Plant material and growth conditions

The plant genotypes used in these experiments were wild-type A. thaliana (L.) Heyhn. Columbia-0, FAN-GFP (Liu et al. 2023), and S. lycopersicum L. var. M82. Surface-sterilized A. thaliana seeds were distributed in square Petri dishes with TK1 medium (Pasternak et al. 2023) supplemented with 150 mg/L Bacto Tryptone (Roth, Karlsruhe, Germany), 1% sucrose, 5 mM (N- morpholino) ethanesulfonic acid (MES), and 1.1% agar (Roth) adjusted to a pH of 5.6. These plates were maintained at room temperature for 4 h to allow seed imbibition. After sowing, the plates were maintained at 4 °C for 12 h. The seedlings were cultivated on vertically oriented plates within a plant growth chamber at 22 °C under continuous white light. S. lycopersicum seeds were surface-sterilized in 25% hypochlorite for 30 min and sown in TK4 medium supplemented with 2% sucrose (Pasternak et al. 2021).

5-ethynyl 2´-deoxyuridine (EdU) labeling and immunolocalization

For labeling with EdU, three to five four-day-old A. thaliana seedlings were pre-incubated in TK1 liquid medium for 120 min in several vials as technical repeats, and further processed as described previously (Pasternak et al. 2022). In brief, 10 μM EdU (final concentration) was added to each vial at different time intervals (120, 310, and 355 min), and the seedlings were fixed in microtubule-stabilizing buffer (MTSB) 420 min after transfer to TK1 liquid medium. Following fixation of seedlings was carried out in a 4% formaldehyde solution in MTSB with a short vacuum pulse (2–3 min) followed by shaking for 45 min. Subsequently, the seedlings were treated with warm methanol (35–40 °C) for 15–20 min, followed by rehydration. Next, the seedlings were subjected to a series of treatments including cell wall digestion (35 min), membrane permeabilization (35 min), blocking with 2% BSA (30 min), primary antibody incubation (16 h at 4 °C), incubation with the EdU detection solution (45 min) with Alexa 647 as a marker, secondary antibody incubation (45 min), washing, and labeling with 0.3 mg/L of 4′,6-diamidine-2′-phenylindole (DAPI) or with SCRI Renaissance 2200 (SR2200) (dilution 1:10,000), before being mounted on a microscopic slide (Pasternak et al. 2022).

Immunolocalization was performed in A. thaliana Col-0 and S. lycopersicum M82 seedlings using a previously published whole mount in situ protocol (Pasternak et al. 2015). The affinity-purified anti-histone H3 antibodies utilized are detailed in Table 1. The Alexa-488/Alexa-555-conjugated anti-mouse or anti-rabbit secondary antibodies were diluted 1:400.Table 1 Antibodies and dilutions used in this work

Antibodies	Host species	Clonality	Dilution used (WB; ICH)	Source	Product no	
Anti-H3K4me1 (monomethyl Lys4)	Rabbit	Polyclonal	1:1000; 1:200	Abcam	ab1791	
Anti-H3K4me2 (dimethyl Lys4)	Rabbit	Monoclonal	1:5000; 1:200	Abcam	ab176882	
Anti-H3K4me3 (trimethyl Lys4)	Rabbit	Polyclonal	1:1000; 1:200	Agrisera	AS16 3190	
Anti-H3K9me1 (monomethyl Lys9)	Rabbit	Polyclonal	1:1000; 1:200	Abcam	ab9045	
Anti-H3K9me2 (dimethyl Lys9)	Mouse	Monoclonal	1:1000; 1:300	Abcam	ab1220	
Anti-H3K9me3 (trimethyl Lys9)	Rabbit	Monoclonal	1:1000; 1:200	Abcam	ab176916	
Anti-H3K27me1 (monomethyl Lys27)	Rabbit	Polyclonal	1:1000; 1:200	Sigma Aldrich	07–448	
Anti-H3K27me3 (trimethyl Lys27)	Rabbit	Monoclonal	1:1000; 1:200	Sigma Aldrich	SAB5600009	
Anti-H3K36me2 (dimethyl Lys36)	Rabbit	Polyclonal	1:1000; 1:200	GeneTex	GTX630555	
Anti-H3K36me3 (trimethyl Lys36)	Rabbit	Polyclonal	1:1000; 1:300	Abcam	ab9050	
Anti-H3K9Ac (acetyl Lys9)	Rabbit	Polyclonal	1:5000; 1:200	Sigma Aldrich	07–352	
WB: Western blotting, ICH: immunohistochemistry

Histone purification, ELISA, and Western blotting

Nuclear proteins were extracted from two-week-old seedlings of A. thaliana Col-0 and S. lycopersicum M82 under standard growing conditions by using the EpiQuik Total Histone Extraction Kit (EpigenTek, Farmingdale, NY, USA). The colorimetric assay EpiQuik Histone H3 Modification Multiplex Assay Kit (EpigenTek) was employed to screen for several histone H3 modification sites in the seedlings of A. thaliana Col-0 and S. lycopersicum M82 samples. This experiment employed an enzyme-linked immunosorbent assay (ELISA), in which several histone H3 modifications were captured simultaneously with their respective antibodies in a 96-well plate. Absorbance at 450 nm was then measured after colour development. Quantification and normalization were conducted in accordance with the manufacturer’s protocol. The purified nuclear proteins were separated on a 12% SDS–polyacrylamide gel and subsequently transferred to PVDF blotting membranes (GE HealthCare Technologies, Chicago, IL, USA) as indicated elsewhere (Kenesi et al. 2023). The PVDF membrane was incubated in the blocking solution overnight. The blocking solution employed was TBS-T (150 mM NaCl, 50 mM Tris–HCl, 0.2% Tween; pH 7.5) with 5% non-fat milk (Kenesi et al. 2023). Subsequently, the membranes were washed with TBS-T at least three times. Subsequently, the membranes were incubated with histone-specific primary antibodies (Table 1) for 4 h at room temperature. The membranes were then washed with TBS-T and incubated for 1 h with secondary antibodies (Goat-Anti-Rabbit-HRP: 1:10,000 [Agrisera, AS09602], Goat-Anti-mouse-HRP: 1:5000 [GE HealthCare, NA 931VS]) at room temperature. Subsequently, the immunoreactive signal was developed using WesternSure PREMIUM chemiluminescent substrate, and detection was performed using the ChemiDoc XRS + imaging system (Bio-Rad, Hercules, CA, USA).

Image acquisition and processing

The TCS SP8 confocal microscope (Leica Microsystems, Mannheim, Germany) controlled by LAS X v3 software was used for image acquisition. The samples were imaged using a 20 × /0.75 air immersion objective for observation and an 86 × /1.20 water immersion objective (HC PL APO CS2) for image acquisition. Each acquisition consisted of a Z-stack capture of a single root with a 0.3−0.5 µm Z-slice distance, comprising 2–4 tiles. The excitation (Ex.) and emission (Em.) wavelengths used were as follows: DAPI and SR2200 (Ex./Em.): 405/424−503; Alexa 488 (Ex./Em.): 488/504−555; Alexa 555 (Ex./Em.): 555/558−643; Alexa 647 (Ex./Em.): 651/658−834. The images were converted to HDF5, stitched using the XUV tool, and saved as HDF5. The HDF5 files were utilized directly for root analysis or converted to TIFF for nuclei analysis.

Nuclei/chromatin/replication analysis

Whole DNA staining was carried out by employing either DAPI or propidium iodide (PI) for nuclei labelling. Besides, the SCRI Renaissance 2200 (SR2200) or calcofluor white stain was also applied to visualize the cell walls (Robert et al. 2013). Nuclei detection and cell segmentation were performed according to a previous protocol (Pasternak and Pérez-Pérez 2021). Nuclei/chromatin analysis was conducted using the 3D Object Detection plugin in Fiji (Schindelin et al. 2012). The parameters analyzed for the selected nuclei included their relative XYZ position, volume occupied by the DNA, surface area, total fluorescent signal from the nuclei, nuclei width, height, depth (W, H, D) values, and relative distance between the center of mass and centroid, which were converted to Euclidian vectors based on the identified parameters in the 3D Object Counter plug-in (Bolte and Cordelières 2006) in Fiji (Schindelin et al. 2012). The measured dataset was saved as CSV files and exported to XLS files. The data was manually curated to exclude false positives (fused nuclei) and false negatives (background noise). The Cartesian coordinates were converted to cylindrical coordinates for the entire root, while individual nuclei were converted to spherical coordinates. Subsequently, cell fate was assigned to each nucleus based on its relative position along the longitudinal and radial axes of the root. Alternatively, the image file was exported to HDF5 and loading it into the iRocS toolbox (Schmidt et al. 2014; Pasternak et al. 2024). For the individual nuclei analysis, multiple nuclei sequences of the same type were selected, with each channel analyzed independently. The data were extracted according to the 3D Object Counter plugin (Bolte and Cordelières 2006) in Fiji (Schindelin et al. 2012).

To ascertain the extent of overlap between the DAPI and H3K9me1/2 signals, nuclei from a range of cell types and regions within the RAM were analyzed. The ImageJ program was employed to analyze each nucleus, selecting the equatorial plane, and creating a projection on the z-axis of a stack size. For each channel (DAPI and H3K9me1), a reference line is defined at the periphery of the nucleus. This enables the recovery of the coordinates and intensity values of the drawn line. To quantify the co-localization of the H3K4me1 and H3K9me2 signals, the ImageJ threshold color tool was employed. This tool allows for the differentiation of pixels with different color ranges within an RGB image. The nuclei of interest were isolated from the global image using tools such as the rectangle and the mirror. For the colocalization analysis, the composition of channels whose signal is of interest was maintained, and a range of colors was established to determine each region of interest. The total signal area is defined as ranging from 0 to 255, the exclusive red signal area is defined as ranging from 0 to 15, the signal area corresponding to the overlay is defined as ranging from 16 to 58, and the exclusive green signal area is defined as ranging from 59 to 255.

DNA replication analysis

To analyze DNA replication, the labeled roots were sequentially scanned in four channels (nuclei, H3K4me3, H3K9me1/2, and EdU). To prevent saturation, the signal intensity in all channels was adjusted. Subsequently, the resulting images were exported to ImageJ, where the channels were split, and independent nuclei analysis was performed for each channel. The quantitative data pertaining to the level of DNA replication, as indicated by the integral EdU signal intensity for each nucleus, have been extracted. While in the case of heterochromatin replication, the intensity of each EdU spot was quantified for each nucleus.

Results

Establishing an experimental pipeline for a high-throughput study of chromatin structure in roots

An experimental procedure for DNA and cell wall staining was followed to characterize in detail the quantitative nuclear morphology and cell properties of root tips, including the identification of various euchromatin and heterochromatin marks. Serial stacks and tiles of the root tip from the columella and lateral root cap (LRC) cells to the elongation zone of the PRs, were obtained in both A. thaliana (Fig. 1) and tomato (Suppl. Fig. S1) seedlings. In the proliferation domain of the RAM region of both A. thaliana and tomato roots, nuclei exhibited intense DNA staining in discrete spots or chromocenters near the nuclear periphery (Fig. 1A, E and Suppl. Fig. S1). Moreover, no DNA staining was observed in the innermost region of the nuclei, which resembles the nucleoli. This was confirmed by the absence of overlap between the nucleolus-localized FAN-GFP fusion protein (Liu et al. 2023) and DAPI staining (Fig. 1A, B).Fig. 1 Methodology used to characterize nuclear heterogeneity in the RAM. A A representative A. thaliana primary root (PR) stained with DAPI for nuclei visualization, showing a longitudinal zonation pattern (Ivanov and Dubrovsky 2013). PD cell proliferation domain, TD transition domain, EZ elongation zone, RAM root apical meristem. B Nucleoli labeling with the FAN-GFP protein (magenta) and counterstained with DAPI (gray). C, D Immunolocalization was performed using the antibodies against the histone mark H3K4me1 (C) and H3K9me1 (D) in PRs of A. thaliana. E, F Detail of confocal images displaying the cells surrounding the QC that exhibit specific patterns of euchromatin (H3K4me1, E) and heterochromatin (H3K9me1, F). Yellow asterisks in A, E, and F indicate the quiescent center (QC) cells

A total of 21 distinct histone H3 modifications were screened and quantified in nuclear extracts of young A. thaliana and tomato seedlings using a commercially available ELISA kit (Suppl. Fig. S2A). Some of the H3 modifications that have previously been correlated with different chromatin domains (Wang et al. 2015; Bi et al. 2017) were selected for further investigation (Table 1). In particular, H3K9me2/me3 is known to be associated with constitutive heterochromatin, whereas H3K4me3 and H3K4me1 are linked to euchromatin, given that they are involved in gene activation and transcription elongation, respectively (Hu and Du 2022), Among the H3 modifications that were tested, negative results were obtained for antibodies against H3K9me3, H3K27me3, and H3K36me2. Furthermore, the polyclonal antibodies against H3K27me1, and H3K36me3 showed non-specific bands of larger molecular weight and cannot be used for in situ maps as they produced a strong background during immunolocalization in root tissues (Suppl. Fig. S2B, C).

In our immunolocalization studies, H3K4me1, H3K4me2, H3K9me1, and H3K9me2 exhibited a consistent and precise nuclear pattern in all the roots studied for both species (Suppl. Fig. S2C). The euchromatic H3K4me1 mark (Fuchs et al. 2006) was observed to be homogenously distributed in all detectable nuclei within the A. thaliana and S. lycopersicum root tips, with no presence in the nucleolus (Fig. 1C, E). In contrast, H3K9me1 and H3K9me2 marks were observed to accumulate in discrete speckles around the nuclear periphery (Fig. 1D, E), which is a characteristic of heterochromatin (Baroux et al. 2007). The extent of colocalization between the H3K9me1, H3K9me2 signal and the chromocenters observed in DAPI staining was determined by comparing the signal intensity at the periphery of the central section of the nucleus. A strong colocalization between H3K9me1 and H3K9me2 was observed in the outer cell layers of the proliferation domain of the RAM region, irrespective of tissue type (Suppl. Fig. S3). Moreover, the chromatin staining profiles of adjacent nuclei of the same cell type in the proliferation domain of the RAM region revealed significant heterogeneity (Suppl. Fig. S3). Therefore, this nuclear heterogeneity as regards H3K9me1 and H3K9me2 may be indicative of the differentiation status of individual nuclei.

To ascertain the proportion of nuclei exhibiting replicating DNA, EdU labeling was employed at varying incubation times from 45 to 300 min, as described elsewhere (Pasternak et al. 2022). We focused both on the proliferation domain and the transition domain of the RAM region, which is where active DNA replication occurred. As an illustration, 45 min incubation with EdU is shown in Fig. 2A. A distinctive pattern of EdU staining was observed, extending from the proliferation domain to the transition domain of the RAM region (Pasternak et al. 2022). The data revealed the presence of multiple patterns of EdU incorporation. These included a homogeneous EdU staining pattern (S1), a spotted, incomplete EdU staining pattern that was colocalized with heterochromatin (S2), and combinations of both (Fig. 2B, C).Fig. 2 Detection of DNA replication in the RAM. A A representative A. thaliana root tip stained with EdU for measuring DNA replication. B Three EdU staining patterns were distinguished: S0, nuclei exclusively stained with DAPI; S1, homogeneous EdU staining; S2, characterized by the presence of EdU spots. C The study analyzed the proportions of S0, S1, and S2 nuclei at 45, 90, and 300 min of EdU incubation. Letters indicate statistically significant differences in the proportion of cells of each pattern depending on the incubation time with EdU (P-value < 0.05; Chi-squared test)

We conducted double, triple and quadruple labeling with histone marks (H3K9me1, H3K9me2, H3K4me1) in conjunction with DAPI for DNA staining, SR2200 for cell wall visualization (Fig. 3 and Suppl. Movie 1), and EdU for visualization of DNA replication. The results of our analyses demonstrated that the triple labeling of H3K9me1, H3K9me2, and DAPI in tomato roots exhibited a tissue-specific localization pattern for heterochromatin marks in the LRC cells (Fig. 3A). Moreover, triple labeling with H3K4me1, H3K9me2, and EdU was successfully achieved for the entire A. thaliana RAM region (Fig. 3E). This approach enabled us to characterize the cellular morphology and nuclear heterogeneity in different cell types and throughout cell division with topological precision. Furthermore, subsequent studies were conducted on H3K9me2, given its precise localization to chromocenters in both species. We have elected to concentrate our efforts on A. thaliana seedlings, given their relatively low genome complexity and the extensive availability of genetic tools.Fig. 3 Multiple combinations for the characterization of nuclear heterogeneity. A Triple labeling of a tomato lateral root (LR) with DAPI (gray), H3K9me1 (cyan), and H3K9me2 (magenta). Lower panel is a magnification of the yellow inset from upper panel. B Triple staining with DAPI (gray), H3K4me1 (green), and H3K9me2 (magenta). C Double staining with H3K4me1 (green), and H3K9me2 (magenta). D Triple staining with SR2200 (gray), H3K4me1 (green), and H3K9me2 (magenta). E Quadruple staining with DAPI (gray), H3K4me1 (green), H3K9me2 (magenta), and EdU (cyan). A-D was performed in tomato LRs, and E was performed in A. thaliana PRs

Nuclear morphology and cellular parameters estimation along the axial dimension of the RAM

As illustrated in Fig. 1, there are notable differences in the morphology of the nuclei and the internal organization of chromatin among cell types within the RAM of A. thaliana, including the quiescent center (QC) region and the proliferation domain. To quantify these variations, a 3D-object analysis was conducted on all nuclei in single roots up to 600−700 µm (Fig. 4A, B). A regression analysis of the results from three representative cell types, including the LRC, epidermis, and pericycle, revealed significant differences in DNA volume and nuclear shape (Fig. 4C, D). The study discovered that the volume occupied by DNA of the LRC and pericycle cells remained stable, while the DNA volume of the epidermal cells increased progressively (Fig. 4C). Furthermore, the LRC nuclei elongated significantly along the root’s longitudinal axis, while the pericycle nuclei remained almost spherical (Fig. 4D). Moreover, based on the integral intensity of DNA labeling, we constructed maps of the relative nuclear ploidy distribution (Fig. 4E, F). A correlation was observed between the level of DNA fluorescence and the position of nuclei along the longitudinal axis of the root in the outer cell layers (Fig. 4E, F). This correlation defines the proliferation domain and the transition domain within the RAM region, thereby corroborating the findings of conventional ploidy studies but with otherwise spatial resolution.Fig. 4 Nuclear parameter studies in different tissues. A A representative PR of A. thaliana stained with DAPI used for ploidy estimation. 3D view from 12 central sections is shown. White arrowheads show the position of the last mitosis of different tissues that define the proliferation domain (PD). B nuclei segmentation with 3D Object Counter plugin. C DNA volume in the indicated cell layers was estimated from the DAPI signal, and D the nuclei shape factor was determined by the width-to-height ratio of the fluorescent DAPI area. E In situ DNA labeling intensity map. About 2,500 nuclei were represented according to the radial axis of the root and their distance from the QC. The fluorescence of the DAPI signal is represented by the color intensity bar in the right side. F Frequency histogram based on DAPI fluorescence analysis of these nuclei

The root epidermis is composed of two distinct cell types, trichoblasts (T) and atrichoblasts (AT), which exhibit varying cellular attributes along the longitudinal and radial growth axes (Fig. 5A, B). In longitudinal confocal sections along the root surface covering the proliferation domain and the transition domain of the RAM region (Fig. 5C and Suppl. Movie 2), it was observed that T cells in the proliferation domain were approximately 1.8 times thinner than neighboring AT cells (Fig. 5D, distance between nuclei). Although the nuclei size in T and AT cells was similar (Fig. 5D, nuclei area), the nuclei of T cells exhibited a flattened morphology along the longitudinal axis of the root compared to those of AT cells. This was confirmed by the pronounced disparities in the shape factor values between T and AT nuclei (Fig. 5D, ratio width/height). The nucleoli area of T cells was found to be significantly larger than that of neighboring AT cells (Fig. 5D, nucleoli area).Fig. 5 Spatial distribution of nuclear parameters in the root epidermis. A Triple labelling and staining with H3K4me1 (green), H3K9me2 (magenta), and DAPI (white) in a representative file of T and AT cells. B, C Spatial distribution of AT and T and their underlying cortical (Co) cells in a representative A. thaliana root; cross-section (B) and longitudinal view (C) of the proliferation domain (PD) and the elongation domain (EZ) of the RAM region. D Analysis of nuclear parameters related to nuclei shape and size in T and AT cells. E–H Morphometric analysis of subcellular parameters analyzed in cortex and overlying AT cells, including the relative position of the nucleus in the cell (F), measured by the distances from the edge of the nucleus to the outer cell wall (L1) and the inner cell wall (L2), as well as the number (G) and diameter (H) of vacuoles per cell (yellow asterisks). Letters in F–H indicate statistically significant differences between cell types (P-value < 0.01; LSD)

At the proliferation domain of the root, the epidermal and cortex nuclei remained in a central position within the cell (Fig. 5E). However, as the nuclei proceeded through the transition domain and into the elongation zone, the epidermal nuclei migrated within the cell towards the center on the radial root axis, while the cortex nuclei moved away from it in an opposite, centripetal direction (Fig. 5F). The observed nuclear migration within epidermal and cortical cells may be dependent on vacuole morphology, as has been demonstrated for polar nuclear migration within T cells during root hair elongation (Nakamura et al. 2018). To test this hypothesis, we examined the number and size of vacuoles in AT and cortical cells located in the proliferation domain and the transition domain of the RAM region. It was observed that the number of vacuoles in both the AT and cortex cells was greater in the proliferation domain than in the transition domain (Fig. 5G). As the nuclei proceeded through the transition domain and into the elongation zone of the root, the number of vacuoles diminished, and their diameter increased (Fig. 5H). This observation is consistent with the migration patterns observed in the aforementioned cell types.

Spatial distribution of euchromatin and heterochromatin marks during cell division and cell differentiation

The objective of this study was to investigate the dynamic regulation of euchromatin and heterochromatin marks during the mitotic stages in AT cells at the proliferation domain of the RAM region in A. thaliana seedlings. To achieve this, double staining with H3K4me1 and H3K9me2 was employed. A mere 6.1% of AT cells exhibited visible mitotic phases, spanning from metaphase to early telophase, within the proliferation domain of the RAM region. The analysis of 2D and 3D images of interphase nuclei (n = 170) demonstrated that the H3K9me2 mark is predominantly located in the peripheral region of the nucleus in the form of speckles of varying sizes, whereas the H3K4me1 mark exhibits a uniform distribution throughout the DNA and is excluded from the nucleolus (Fig. 6A and Supp. Movie 3). Although some degree of colocalization was observed between the two marks, the majority of their signal did not coincide (Fig. 6A and Suppl Fig. S4A).Fig. 6 Distribution of H3K9me2 and H3K4me1 in AT cells. 2D and 3D views of several nuclei in interphase (A), metaphase (B), and late anaphase (C) are shown. H3K9me2 is shown in magenta/red, and H3K4me1 is depicted in green. Dashed yellow lines in B and C show the equatorial plane of the cell division. D, E The signal intensity for different chromosomic regions during anaphase was measured for H3K9me2 (D) and H3K4me1 (E). Letters indicate statistically significant differences between regions (P-value < 0.01; LSD). F Possible signal distributions considered in the individual nucleus for each mark. Numbers indicate the orientation of the signal in the axial direction. In H, the first digit indicates H3K4me1 orientation and the second, the orientation of H3K9me2. G A representative example of contiguous AT cells with an indication of the signal distribution patterns of H3K4me1 and H3K9me2 assigned. White squares indicate a contiguous nucleus that recently completed mitosis (1_2 pattern). H Frequency histogram of all possible H3K4me1 and H3K9me2 signal distributions in individual nucleus. The statistically significant differences are represented by the asterisk (P-value < 0.01; LSD); n.f.: not found (n = 230 nuclei, 3 roots)

At the metaphase plate, the H3K9me2 heterochromatin mark was observed to be internally located, while the H3K4me1 euchromatin mark was situated in a peripheral position in relation to it (Fig. 6B and Supp. Movie 4). Images captured during late anaphase revealed an uneven distribution of heterochromatin in the two daughter nuclei (Fig. 6C, Supp. Movie 5, and Suppl Fig. S4A), with a greater proportion of the H3K9me2 signal observed in regions situated further away from the equatorial plane of the mitotic plate (Fig. 6D). In contrast, H3K4me1 was distributed uniformly within the two daughter nuclei, exhibiting a slight increase in the signal in regions closer to the equatorial plane (Fig. 6E). It has been proposed that contiguous nuclei with opposing H3K4me1 and H3K9me2 signals along the longitudinal axis of the root, may indicate cells that have recently undergone mitosis. Given the heterogeneous distribution of H3K4me1 and H3K9me2 signals across contiguous AT nuclei, we have devised a nomenclature to delineate the spatial relationship between these two marks in individual nuclei (Fig. 6F). It was postulated that the probability of occurrence for all potential configurations of uneven H3K4me1 and H3K9me2 in single nuclei would be equal if a random chromatin rearrangement occurred following mitosis. However, we observed a higher proportion than expected of interphase nuclei with the opposing distribution of H3K4me1 and H3K9me2, which we have designated the 1_2 configuration (Fig. 6G, H). Accordingly, this configuration indicates that H3K4me1 (euchromatin) and H3K9me2 (heterochromatin) were in opposite compartments of the nucleus. This suggests that the opposite euchromatin and heterochromatin distribution (1_2 configuration) may be more stable in interphase AT nuclei. It is noteworthy that chromatin configurations in which H3K4me1 and H3K9me2 signals colocalize completely in regions orthogonal to the equatorial plate of cell division (designated as 3_3 and 4_4 configurations), were not observed (Fig. 6H). In conclusion, the data indicated that there was a degree of chromatin asymmetry in the nuclei, with a clear separation of euchromatin and heterochromatin marks on opposite sides of the nucleus shortly after mitosis.

A significant difference was identified in the colocalization degree between H3K4me1 and H3K9me2 marks in T and AT cells within the proliferation domain of the RAM region (Suppl. Fig. S4B, C). Notably, the H3K9me2 signal exhibited a higher degree of variability in nuclei at the upper elongation zone and adjacent to the differentiation zone of the root when compared to the H3K4me1 signal in this region (Suppl. Fig. S4D–F and Suppl. Movie 6). The ratio of H3K9me2 to H3K4me1 is higher in the outer radial layers, such as the epidermis and cortex, compared to the internal cells of the pericycle and vascular region, which have a lower ratio of H3K9me1 to H3K4me1 (Suppl. Fig. S4G–I). This may be regarded as an epigenetic signature of the transcription of genes associated with the topology of individual cells within the root tip. Such significant differences in euchromatin topology led to epigenetic and geometric constraints and serve as a reason for root abnormalities under different environmental perturbations, such as those evidenced by carbon starvation of the PRs (Supp. Fig. S5). An intriguing tissue-specific pattern of H3K9me2-enriched heterochromatin was observed in the elongation zone of the RAM. In this region, the levels of H3K9me2 were significantly higher in the epidermis and cortex cells and markedly lower in the inner tissues of the pericycle and vascular cells. These distinctive features, including nuclear enlargement and large-scale chromatin rearrangements, are consistent with cells that have exited the cell cycle exit prior to terminal differentiation. However, this remodeling was not observed in the internal cell layers, including the pericycle, which exhibited a lower ratio of H3K9me2 to H3K4me1 in their nuclei. These findings are in accordance with the hypothesis of a dome-shaped structure of the RAM, whereby the outer layers undergo differentiation at a higher rate than the inner tissues, as previously proposed (Lavrekha et al. 2017).

Staining of DNA replication as a tool to measure cell cycle synchronization

To evaluate DNA replication patterns, EdU was employed on A. thaliana roots following distinct incubation periods from 45 to 300 min (Suppl. Fig. S6A). Nuclei were classified based on the distinct EdU staining patterns observed as homogeneous EdU staining pattern (S1) and spotted, incomplete EdU staining pattern (S2) (Fig. 7A and Supp. Movie 7). Following a brief 45-min incubation period with EdU allowed us to visualize only cells during DNA replication or shortly after (Pasternak et al 2022). It was observed that between 10% (LRC cells, n = 568) and 20% (epidermis [n = 297], and cortex [n = 122] cells) of the nuclei exhibited EdU staining (Suppl. Fig. S6B). Additionally, mitotic chromosomes were observed in only 1.6% of the analyzed nuclei (n = 987). Following an extended incubation period with EdU (300 min), the proportion of stained nuclei increased to 60% in epidermis (n = 483) and cortex (n = 203) tissues, while in LRC cells (n = 621), the proportion reached 30% after 300 min of EdU incubation (Fig. 7A and Suppl. Fig. S6A, B). Moreover, the proportion of nuclei undergoing mitosis exhibited variability across different cell types, with 0.63% in LRC (n = 1,595) and 1.58% in cortex (n = 569), indicating a differential division rate between tissues, as previously reported (Pasternak et al. 2022). In these conditions, a progressive increase in the proportion of nuclei in pattern S1 was observed in all cell types from 45 to 300 min of EdU incubation. However, the proportion of nuclei in pattern S2 remained constant even after prolonged incubation with EdU (Suppl. Fig. S6B). This suggests that pattern S1 identifies nuclei in the early S phase, whereas pattern S2 identifies nuclei in the late S phase. Furthermore, to ascertain which type of chromatin replicates DNA during the early or late S phase, we conducted an analysis of the colocalization of the H3K4me1 or H3K9me2 marks with the S1 and S2 patterns of EdU in the epidermis (Suppl. Fig. S6C, D). In both T and AT cell types, a greater degree of colocalization of the euchromatin mark H3K4me1 was observed in the nuclei of the pattern S1, whereas the pattern S2 exhibited a greater degree of colocalization with the H3K9me2 mark, which is characteristic of heterochromatin (Suppl. Fig. S6E, F). These findings align with the proposed hypothesis that DNA replication initiates concurrently in broad euchromatic regions within the nucleus during the early S phase (S1 pattern), while heterochromatin replication occurs later (S2 pattern) (Costas et al. 2011; Concia et al. 2018).Fig. 7 Quantification of DNA replication in the RAM. A A representative A. thaliana root tip used for the quantification of DNA replication. Seedlings were labelled with EdU for 90 min, scanned with a confocal microscope, stitched and subjected to 3D segmentation. B The analyzed nuclei were represented along the longitudinal and radial growth axes. Bubble size indicates the intensity of fluorescence per nuclei. C Triple labelling and staining with H3K4me1 (green), H3K9me2 (magenta), and EdU (white) in a representative file of T cells along the region between the PD and the EZ, dubbed as the transition domain (TD). The third panel represents the results from EdU segmentation. White asterisks indicate nuclei that were not correctly segmented or with low/no fluorescence signal. D Histograms with relative fluorescence of the EdU (grey, logarithmic scale) and H3K4me2 (green, linear scale) on the nuclei depicted in C. S0 refers to cells without EdU staining. Orange lines in C and D mark cells in the EZ. E A representative unrolled root epidermis showing AT nuclei in S0 (blue), S1 (orange), and S2 (grey) patterns. Nuclei with mitotic figures are depicted as green diamonds. The orange dashed line indicates the PD/TD boundary, which was defined by cell size increase. F The distribution of clusters of nuclei in the analyzed tissues of AT cells (top) and cortex cells (bottom) for the S0 and S1/S2 patterns of EdU staining are represented in histograms. Letters indicate significant differences between grouping in S0 and S1/S2 nuclei

We postulated that the fluorescence intensity of H3K4me1, H3K9me2, and EdU could be employed to estimate different parameters (volume, surface, etc.,) of euchromatin, heterochromatin, and replicating DNA, respectively, in individual nuclei with spatial resolution. Following a 90-min period of EdU incubation, a notable variation in EdU fluorescence was observed along the longitudinal and radial growth axes of the root tip (Fig. 7B). To further investigate the relationship between chromatin structure and DNA replication, we conducted triple staining after 45 min of incubation with EdU, in defined regions of the RAM. The quantification of the EdU signal was employed to identify nuclei that were not engaged in DNA replication (i.e., exhibiting the S0 pattern) and to determine the extent of DNA replication in neighboring nuclei. In a representative sample of T nuclei from the region between the proliferation domain and the elongation zone, we observed that H3K4me1 levels were similar in the proliferation domain region and increased 2.5-fold in the elongation zone (Fig. 7C, D), which is consistent with the observed increase in ploidy levels in this cell type (Fig. 3A). Regular patterns and levels of EdU staining were observed in cells in the proliferation domain, with contiguous cells displaying the same pattern and similar levels of fluorescence localization. In the elongation zone, nuclei mostly exhibited high levels of EdU signal (Fig. 7C, D), which colocalize with euchromatin, indicating coordinated entry into the endoreduplication phase. Following a brief pulse labeling with EdU, patches of nuclei with similar staining patterns were predominantly observed along longitudinal cell files, including the epidermis and cortex (Fig. 7E). The nuclei that were not stained with EdU (pattern S0) exhibited an inverse Gaussian distribution in both tissues analyzed, AT cells and cortex cells (Fig. 7F), with no significant differences between them (P-value = 0.532). In contrast, the nuclei stained with EdU were most frequently observed in groups of one (42.2%) or two (31.8%) nuclei, with no significant differences between AT and cortex nuclei (P-value = 0.916). No significant differences were observed in the grouping pattern of non-EdU stained nuclei between the proliferation domain and the elongation zone in any of the analyzed tissues (P-value = 0.118). The average number of EdU-stained nuclei was 2.3 in the proliferation domain and 1.7 in the elongation zone for both tissues (P-value = 0.012). A comprehensive examination of the alterations in nuclear morphology and the distribution of euchromatin and heterochromatin marks in diverse cell types within the RAM has enabled us to discern significant discrepancies in the growth and nuclear organization of cells in adjacent layers. Furthermore, the intracellular migration of their nuclei suggests the potential existence of regulatory mechanisms that govern these processes. Our long-term objective is to adapt this approach to other crops to elucidate the dynamics of chromatin maps and the subsequent impact on gene expression associated with development and stress tolerance.

Discussion

The present work outlines a methodology for the establishment of an in situ chromatin profile with cellular resolution. The procedure entailed staining the cell border, nucleus, euchromatin, and heterochromatin, as well as the detection of DNA replication events. The combination of a cylindrical coordinate system with the iRoCS toolbox (Schmidt et al. 2014) and graphical visualization tools like ggPlantmap (Jo and Kajala 2024) enables the creation of digital 3D maps of the RAM, encompassing cellular, nuclear, and chromatin features for each cell.

Histone modifications, including acetylation, methylation, and phosphorylation, exert a direct influence on the structure of chromatin and, consequently, on nuclear morphology. This influence is exerted either directly or by regulating the binding of chromatin-associated factors (Bannister and Kouzarides 2011; Tessarz and Kouzarides 2014). The present study has confirmed the specificity of several chromatin methylation marks in H3K4 and H3K9 in the roots of A. thaliana and tomato. Two of these histone modifications are associated with heterochromatin while the other two are associated with euchromatin. The H3K4me1 mark was asymmetrically distributed throughout the nucleus, exhibiting a dense and highly organized distribution in the proliferation domain of the RAM region and being absent from the nucleoli in all studied cells. This finding is consistent with its known location in active gene enhancers situated within the euchromatic regions of the genome (Zhang et al. 2009). The chromocenters observed in the DAPI staining of the nuclear DNA were found to correspond to the heterochromatin regions marked with H3K9me1 and H3K9me2 (Fig. 5). Chromocenters are typically situated in the outermost region of the nucleus, in close proximity to the nuclear envelope. A noteworthy diversity in chromocenter number, size, and spatial distribution was observed across different cell types and along the longitudinal axis of the root, which were found to be associated with specific cell fate and/or cell cycle stages. Specifically, the H3K9me2 signal displayed increased intensity in the outer cell layers of both A. thaliana and tomato roots.

H3K9me2 is a conserved mark of peripheral heterochromatin that is closely associated with the nuclear lamina in animal cells (Poleshko et al. 2019). In A. thaliana, the H3K9me2 mark is present in extensive blocks of pericentromeric heterochromatin, as well as in specific regions comprising euchromatin repeats and transposons (Roudier et al. 2011; Ramakrishnan Chandra et al. 2024). In interphase nuclei, the H3K9me2 heterochromatin mark is located in specific speckles along the peripheral region of the nucleus and in close proximity to the nuclear envelope, with minimal overlap with the euchromatin H3K4me1 mark (Pei et al. 2021). The distribution of the latter is homogeneous within the interphase nucleus. The observed dynamic rearrangement of H3K9me2 and H3K4me1 during metaphase and anaphase is consistent with the hypothesis that H3K9me2-enriched pericentromeric heterochromatin plays a role in anchoring the mitotic spindle during mitosis (Yamada and Goshima 2017; Ramakrishnan Chandra et al. 2024). During metaphase, the mitotic spindle assembles with centromeric heterochromatin. In late anaphase, kinetochore microtubules undergo a further shortening, which contributes to the separation of heterochromatin towards the poles of the equatorial plane of mitosis. Euchromatin, which is less condensed, remains in the region opposite to that of the condensed chromatin, situated in close proximity to the equatorial plane. Following the conclusion of mitosis and the reconstruction of the nuclear envelope, the two contiguous daughter cells exhibit concave and convex patterns, respectively, of H3K4me1 nuclear distribution (white square in Fig. 6G). Moreover, following the formation of the cell plate, H3K9me2 exhibited contrasting distribution patterns (1_2 configuration). Our study identified approximately 40% of AT nuclei exhibiting this pattern, indicating that these cells have completed mitosis and are in the early G1 phase. It is noteworthy that the spatial distribution patterns of the H3K4me1 and H3K9me2 marks in the remaining 60% of the analyzed AT nuclei exhibited an independent and random distribution, with the exception of the absence of patterns in which these two marks were present in the region of the nucleus closest to the T cells. These findings suggest the possibility of spatial constraints on chromatin distribution between cells across different rows, but not within the same row. Further experimentation is required to extend this observation to other tissues, including the cortex and endodermis.

The relative DNA content of each nucleus can be estimated from the provided dataset by integrating the H3K4me1 signal and/or DAPI fluorescence data. It was assumed that the nuclei of the initial cells surrounding the QC exhibited a DNA level of 2C (Bhosale et al. 2018), which can be used as a reference for estimating the ploidy levels of the other nuclei under study. Cylindrical coordinates were assigned to individual nuclei in the iRoCS toolbox (Schmidt et al. 2014) for the purpose of constructing in situ ploidy maps for all the tissue layers within a single root. The maps can be utilized to characterize known mutants affected in endoreduplication and cell cycle progression, such as kip-related protein 5 (Jégu et al. 2013) and cyclin D3 triple mutants (Dewitte et al. 2007). The DNA volume and shape factor of the nucleus exhibited variation depending on the cell type and differentiation status along the longitudinal root axis. The nuclei of differentiating LRC cells underwent a transformation from a spherical shape with homogeneous DNA distribution to an elongated shape with marked chromocenters, accompanied by minor fluctuations in the amount of DNA per nucleus. In contrast, the nuclei observed at the root epidermis exhibited a notable increase in DNA levels, extending beyond the RAM region of the root. This observation is consistent with the endoreduplication of DNA that occurs in this region for the epidermal cells and that precedes their rapid cell expansion at the elongation zone of the root (Hayashi et al. 2013).

Mitosis is a relatively brief process, resulting in a low percentage of cells undergoing mitosis within the RAM region at a given time (Pasternak et al. 2022). By incubating young seedlings with EdU for increasing periods of time, we were able to conduct a comprehensive examination of the dynamics of DNA replication during the S-phase of the mitotic cell cycle and the endocycle at a cellular resolution in the RAM. This process is dependent on the activation of the replication origins through ORC1b (Vergara et al. 2023). The staining patterns observed for EdU, which are homogeneous (S1) and spotted (S2), suggest that DNA replication initiates concurrently in broad euchromatic regions within the nucleus during the early S phase (S1 pattern), while heterochromatin replication occurs later (S2 pattern) (Costas et al. 2011; Concia et al. 2018). These findings align with the proposed hypothesis, as evidenced by a stronger overlap between H3K9me2 and EdU signals in the S2 pattern, and a higher degree of overlap of H3K4me1 with EdU in the S2 pattern. The relative abundance of EdU over H3K4me1 was quantified to assess the extent of DNA replication occurring during the mitotic cell cycle and the endocycle. The observed differences in the pattern of EdU-stained nuclei of cells at the transition domain suggest the existence of tissue-specific S-phase durations, which require further investigation.

In the roots of A. thaliana and tomato, the outermost layer of the epidermis is composed of alternating T and AT cells with distinct geometry and cell fate (Berger et al. 1998). T cells are radially connected to two inner cells of the exodermis layer (in the case of the tomato) or the cortex layer (in the case of A. thaliana), while AT cells are connected to a single underlying exodermal or cortical cell. In the EZ, the number of T cells in each longitudinal row is approximately twice that of AT cells. The nuclei of T cells are flattened along the longitudinal axis of the root and exhibit larger nucleoli than those of AT cells, which may result in elevated mechanical stress on the nuclei of the T cells. In both T and AT cells within the proliferation domain of the RAM region, the nuclei are positioned centrally within the cell. These observations are consistent with the hypothesis that T cell undergo a higher frequency of formative divisions than AT cells (Berger et al. 1998). In consideration of the three-dimensional structure of the epidermis and the subepidermal layer, it becomes evident that a supracellular mechanism is indispensable for the coordination of the expansion of AT cells and their underlying cortex/exodermis cells with the cell division of the T cells with which they are in contact. It has been proposed that in root cells, brassinosteroid-mediated microtubule arrangement regulates cell shape through differential regulation of directional growth between the inner and outer tissues (Fridman et al. 2021). Given that brassinosteroid insensitive 1 mutants exhibit a reduced number of divisions in the epidermis (Fridman et al. 2021), it is imperative to conduct a comprehensive examination of the cellular and nuclear characteristics of T and AT cells in these mutants to elucidate the mechanism of growth coordination in roots.

In plant cells, the nuclei undergo a change in position in response to alterations in the surrounding environment, the polar growth of specific cells, and during cell division. External factors, such as the direction of light or the presence of other cells, have been demonstrated to influence this positioning (Wada 2018). In the leaves of A. thaliana, the nuclei of mesophyll cells are observed to relocate to the side walls following continuous exposure to intense blue light. This relocation is proposed to serve as a protective mechanism against ultraviolet (UV)-induced DNA damage (Iwabuchi et al. 2016). The polar nuclear movement observed during root hair growth in T cells at the differentiation zone of the root is driven by an as yet unidentified actin-based mechanism that is dependent on auxin signaling (Nakamura et al. 2018; Brueggeman et al. 2022). Our study revealed a previously unidentified phenomenon: the migration of the nuclei of AT cells and cortex cells in the transition domain of the A. thaliana roots. The results indicate that AT nuclei rapidly migrate towards the inner side of the cell, while nuclei in the cortex cells exhibited a centripetal movement upon reaching the transition domain of the RAM. It is conceivable that the concomitant migration of AT and subjacent cortex nuclei may facilitate the coordination of growth responses among neighboring cells through the mediation of specific short-range non-cell autonomous signals via the plasmodesmata. The results of our study indicate that vacuole biogenesis may be a crucial factor in the migration of AT and cortex nuclei. The shoot gravitropism 2 mutant, which displays abnormalities in its vacuolar structure, exhibited impeded nuclear migration during zygote polarization (Kimata et al. 2019). It would be interesting to analyze whether this mutant and others affected in vacuole biogenesis (Cui et al. 2018) exhibit evident alterations in the migration of the AT nuclei and cortex in their roots.

Conclusions

A comprehensive investigation has been conducted into the distribution of euchromatin and constitutive heterochromatin in the roots of A. thaliana, in conjunction with the analysis of DNA replication. Our findings illustrate that chromatin reorganization is a highly dynamic process occurring during the distinct phases of mitosis and immediately following cytokinesis. Furthermore, we observed an intriguing phenomenon of nuclear migration in certain cells at the transition domain of the RAM region and towards the elongation zone of the root. This phenomenon is associated with vacuole biogenesis and appears to be coordinated between the epidermis and the cells of the underlying cortex. The results demonstrate the structural dissimilarities of the nucleus and its internal organization between disparate cell types (T and AT), and illustrate a correlation between these structural dissimilarities and the processes of cell division and cell differentiation that occur in the RAM region of A. thaliana roots. This constitutes an invaluable tool for the study of mutants affected in the coordination of growth and the investigation of the impact of stress factors on chromatin organization, which will be the focus of future studies.

Supplementary Information

Below is the link to the electronic supplementary material.Supplementary file1 (PDF 3701 KB)

Supplementary file2 (PPTX 47723 KB)

Acknowledgements

We thank Dr. Eduardo Larriba (Instituto de Bioingeniería, UMH, Spain) for his critical comments and valuable suggestions and Profs. José Luis Micol and María Rosa Ponce (Instituto de Bioingeniería, UMH, Spain) for sharing their equipment.

Author contributions

LC, TP, and RR designed and performed all the experiments. LC, TP, and RR analyzed the data. The first draft of the manuscript was written by LC and TP. JMP-P wrote the final manuscript. All authors read and approved the final manuscript.

Funding

Open Access funding provided thanks to the CRUE-CSIC agreement with Springer Nature. This research is supported by the grant PID2021-126840OB-I00 funded by MCIN/AEI/10.13039/501100011033, and by the “ERDF A way of making Europe”. LC holds a “Programa INVESTIGO” contract (INVEST/2022/247) from the Generalitat Valenciana. TP holds a ‘María Zambrano’ distinguished researcher contract at the UMH.

Data availability

All data generated or analyzed during this study are provided in this published article and its supplementary data files or it will be provided upon a reasonable request.

Declarations

Conflict of interest

The authors declare no conflicts of interest.

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Lara Caballero and Taras Pasternak contributed equally to this work.
==== Refs
References

Amiad-Pavlov D Lorber D Bajpai G Live imaging of chromatin distribution reveals novel principles of nuclear architecture and chromatin compartmentalization Sci Adv 2021 7 6251 6253 10.1126/SCIADV.ABF6251
Amiad-Pavlov D, Lorber D, Bajpai G et al (2021) Live imaging of chromatin distribution reveals novel principles of nuclear architecture and chromatin compartmentalization. Sci Adv 7:6251–6253. 10.1126/SCIADV.ABF625110.1126/SCIADV.ABF6251
Bannister AJ Kouzarides T Regulation of chromatin by histone modifications Cell Res 2011 21 381 395 10.1038/cr.2011.22 21321607
Bannister AJ, Kouzarides T (2011) Regulation of chromatin by histone modifications. Cell Res 21:381–395. 10.1038/cr.2011.2221321607 10.1038/cr.2011.22
Baroux C Pecinka A Fuchs J The Triploid endosperm genome of arabidopsis adopts a peculiar, parental-dosage-dependent chromatin organization Plant Cell 2007 19 1782 1794 10.1105/TPC.106.046235 17557811
Baroux C, Pecinka A, Fuchs J et al (2007) The Triploid endosperm genome of arabidopsis adopts a peculiar, parental-dosage-dependent chromatin organization. Plant Cell 19:1782–1794. 10.1105/TPC.106.04623517557811 10.1105/TPC.106.046235
Berger F Hung CY Dolan L Schiefelbein J Control of cell division in the root epidermis of arabidopsis thaliana Dev Biol 1998 194 235 245 10.1006/DBIO.1997.8813 9501025
Berger F, Hung CY, Dolan L, Schiefelbein J (1998) Control of cell division in the root epidermis of arabidopsis thaliana. Dev Biol 194:235–245. 10.1006/DBIO.1997.88139501025 10.1006/DBIO.1997.8813
Bhosale R Boudolf V Cuevas F A Spatiotemporal DNA endoploidy map of the Arabidopsis root reveals roles for the endocycle in root development and stress adaptation Plant Cell 2018 30 2330 2351 10.1105/TPC.17.00983 30115738
Bhosale R, Boudolf V, Cuevas F et al (2018) A Spatiotemporal DNA endoploidy map of the Arabidopsis root reveals roles for the endocycle in root development and stress adaptation. Plant Cell 30:2330–2351. 10.1105/TPC.17.0098330115738 10.1105/TPC.17.00983
Bi X Cheng YJ Hu B Nonrandom domain organization of the Arabidopsis genome at the nuclear periphery Genome Res 2017 27 1162 1173 10.1101/GR.215186.116 28385710
Bi X, Cheng YJ, Hu B et al (2017) Nonrandom domain organization of the Arabidopsis genome at the nuclear periphery. Genome Res 27:1162–1173. 10.1101/GR.215186.11628385710 10.1101/GR.215186.116
Bolte S Cordelières FP A guided tour into subcellular colocalization analysis in light microscopy J Microsc 2006 224 213 232 10.1111/J.1365-2818.2006.01706.X 17210054
Bolte S, Cordelières FP (2006) A guided tour into subcellular colocalization analysis in light microscopy. J Microsc 224:213–232. 10.1111/J.1365-2818.2006.01706.X17210054 10.1111/J.1365-2818.2006.01706.X
Brueggeman JM Windham IA Nebenführ A Nuclear movement in growing Arabidopsis root hairs involves both actin filaments and microtubules J Exp Bot 2022 73 5388 5399 10.1093/JXB/ERAC207 35554524
Brueggeman JM, Windham IA, Nebenführ A (2022) Nuclear movement in growing Arabidopsis root hairs involves both actin filaments and microtubules. J Exp Bot 73:5388–5399. 10.1093/JXB/ERAC20735554524 10.1093/JXB/ERAC207
Cantwell H Dey G Nuclear size and shape control Semin Cell Dev Biol 2022 130 90 97 10.1016/j.semcdb.2021.10.013 34776332
Cantwell H, Dey G (2022) Nuclear size and shape control. Semin Cell Dev Biol 130:90–97. 10.1016/j.semcdb.2021.10.01334776332 10.1016/j.semcdb.2021.10.013
Carotenuto G Sciascia I Oddi L Size matters: three methods for estimating nuclear size in mycorrhizal roots of Medicago truncatula by image analysis BMC Plant Biol 2019 19 1 14 10.1186/S12870-019-1791-1 30606102
Carotenuto G, Sciascia I, Oddi L et al (2019) Size matters: three methods for estimating nuclear size in mycorrhizal roots of Medicago truncatula by image analysis. BMC Plant Biol 19:1–14. 10.1186/S12870-019-1791-130606102 10.1186/S12870-019-1791-1
Concia L Brooks AM Wheeler E Genome-wide analysis of the Arabidopsis replication timing program Plant Physiol 2018 176 2166 2185 10.1104/PP.17.01537 29301956
Concia L, Brooks AM, Wheeler E et al (2018) Genome-wide analysis of the Arabidopsis replication timing program. Plant Physiol 176:2166–2185. 10.1104/PP.17.0153729301956 10.1104/PP.17.01537
Costas C De La Paz SM Stroud H Genome-wide mapping of Arabidopsis thaliana origins of DNA replication and their associated epigenetic marks Nat Struct Mol Biol 2011 18 395 400 10.1038/nsmb.1988 21297636
Costas C, De La Paz SM, Stroud H et al (2011) Genome-wide mapping of Arabidopsis thaliana origins of DNA replication and their associated epigenetic marks. Nat Struct Mol Biol 18:395–400. 10.1038/nsmb.198821297636 10.1038/nsmb.1988
Cui Y Cao W He Y A whole-cell electron tomography model of vacuole biogenesis in Arabidopsis root cells Nat Plants 2018 5 95 105 10.1038/s41477-018-0328-1 30559414
Cui Y, Cao W, He Y et al (2018) A whole-cell electron tomography model of vacuole biogenesis in Arabidopsis root cells. Nat Plants 5:95–105. 10.1038/s41477-018-0328-130559414 10.1038/s41477-018-0328-1
Desset S Poulet A Tatout C Quantitative 3D analysis of nuclear morphology and heterochromatin organization from whole-mount plant tissue using NucleusJ Methods Mol Biol 2018 1675 615 632 10.1007/978-1-4939-7318-7_33 29052214
Desset S, Poulet A, Tatout C (2018) Quantitative 3D analysis of nuclear morphology and heterochromatin organization from whole-mount plant tissue using NucleusJ. Methods Mol Biol 1675:615–632. 10.1007/978-1-4939-7318-7_3329052214 10.1007/978-1-4939-7318-7_33
Dewitte W Scofield S Alcasabas AA Arabidopsis CYCD3 D-type cyclins link cell proliferation and endocycles and are rate-limiting for cytokinin responses Proc Natl Acad Sci U S A 2007 104 14537 14542 10.1073/PNAS.0704166104 17726100
Dewitte W, Scofield S, Alcasabas AA et al (2007) Arabidopsis CYCD3 D-type cyclins link cell proliferation and endocycles and are rate-limiting for cytokinin responses. Proc Natl Acad Sci U S A 104:14537–14542. 10.1073/PNAS.070416610417726100 10.1073/PNAS.0704166104
Dubos T Poulet A Gonthier-Gueret C Automated 3D bio-imaging analysis of nuclear organization by NucleusJ 2.0 Nucleus 2020 11 315 329 10.1080/19491034.2020.1845012 33153359
Dubos T, Poulet A, Gonthier-Gueret C et al (2020) Automated 3D bio-imaging analysis of nuclear organization by NucleusJ 2.0. Nucleus 11:315–329. 10.1080/19491034.2020.184501233153359 10.1080/19491034.2020.1845012
Dubos T Poulet A Thomson G NODeJ: an ImageJ plugin for 3D segmentation of nuclear objects BMC Bioinformatics 2022 23 1 11 10.1186/S12859-022-04743-6 34983366
Dubos T, Poulet A, Thomson G et al (2022) NODeJ: an ImageJ plugin for 3D segmentation of nuclear objects. BMC Bioinformatics 23:1–11. 10.1186/S12859-022-04743-634983366 10.1186/S12859-022-04743-6
Dumur T Duncan S Graumann K Probing the 3D architecture of the plant nucleus with microscopy approaches: challenges and solutions Nucleus 2019 10 181 10.1080/19491034.2019.1644592 31362571
Dumur T, Duncan S, Graumann K et al (2019) Probing the 3D architecture of the plant nucleus with microscopy approaches: challenges and solutions. Nucleus 10:181. 10.1080/19491034.2019.164459231362571 10.1080/19491034.2019.1644592
Fridman Y Strauss S Horev G The root meristem is shaped by brassinosteroid control of cell geometry Nat Plants 2021 7 1475 1484 10.1038/s41477-021-01014-9 34782771
Fridman Y, Strauss S, Horev G et al (2021) The root meristem is shaped by brassinosteroid control of cell geometry. Nat Plants 7:1475–1484. 10.1038/s41477-021-01014-934782771 10.1038/s41477-021-01014-9
Fuchs J Demidov D Houben A Schubert I Chromosomal histone modification patterns—from conservation to diversity Trends Plant Sci 2006 11 199 208 10.1016/J.TPLANTS.2006.02.008 16546438
Fuchs J, Demidov D, Houben A, Schubert I (2006) Chromosomal histone modification patterns—from conservation to diversity. Trends Plant Sci 11:199–208. 10.1016/J.TPLANTS.2006.02.00816546438 10.1016/J.TPLANTS.2006.02.008
Goto C Hara-Nishimura I Tamura K Regulation and physiological significance of the nuclear shape in plants Front Plant Sci 2021 12 673905 10.3389/FPLS.2021.673905 34177991
Goto C, Hara-Nishimura I, Tamura K (2021) Regulation and physiological significance of the nuclear shape in plants. Front Plant Sci 12:673905. 10.3389/FPLS.2021.67390534177991 10.3389/FPLS.2021.673905
Hayashi K Hasegawa J Matsunaga S The boundary of the meristematic and elongation zones in roots: endoreduplication precedes rapid cell expansion Sci Rep 2013 3 2723 10.1038/SREP02723 24121463
Hayashi K, Hasegawa J, Matsunaga S (2013) The boundary of the meristematic and elongation zones in roots: endoreduplication precedes rapid cell expansion. Sci Rep 3:2723. 10.1038/SREP0272324121463 10.1038/SREP02723
Hu H Du J Structure and mechanism of histone methylation dynamics in Arabidopsis Curr Opin Plant Biol 2022 67 102211 10.1016/j.pbi.2022.102211 35452951
Hu H, Du J (2022) Structure and mechanism of histone methylation dynamics in Arabidopsis. Curr Opin Plant Biol 67:102211. 10.1016/j.pbi.2022.10221135452951 10.1016/j.pbi.2022.102211
Ivanov VB Dubrovsky JG Longitudinal zonation pattern in plant roots: conflicts and solutions Trends Plant Sci 2013 18 237 243 10.1016/j.tplants.2012.10.002 23123304
Ivanov VB, Dubrovsky JG (2013) Longitudinal zonation pattern in plant roots: conflicts and solutions. Trends Plant Sci 18:237–243. 10.1016/j.tplants.2012.10.00223123304 10.1016/j.tplants.2012.10.002
Iwabuchi K Hidema J Tamura K Plant nuclei move to escape ultraviolet-induced DNA damage and cell death Plant Physiol 2016 170 678 685 10.1104/PP.15.01400 26681797
Iwabuchi K, Hidema J, Tamura K et al (2016) Plant nuclei move to escape ultraviolet-induced DNA damage and cell death. Plant Physiol 170:678–685. 10.1104/PP.15.0140026681797 10.1104/PP.15.01400
Jégu T Latrasse D Delarue M Multiple functions of kip-related Protein5 connect endoreduplication and cell elongation Plant Physiol 2013 161 1694 1705 10.1104/PP.112.212357 23426196
Jégu T, Latrasse D, Delarue M et al (2013) Multiple functions of kip-related Protein5 connect endoreduplication and cell elongation. Plant Physiol 161:1694–1705. 10.1104/PP.112.21235723426196 10.1104/PP.112.212357
Jo L Kajala K ggPlantmap: an open-source R package for the creation of informative and quantitative ggplot maps derived from plant images J Exp Bot 2024 10.1093/JXB/ERAE043 38329371
Jo L, Kajala K (2024) ggPlantmap: an open-source R package for the creation of informative and quantitative ggplot maps derived from plant images. J Exp Bot. 10.1093/JXB/ERAE04338329371 10.1093/JXB/ERAE043
Johann to Berens P Schivre G Theune M Advanced image analysis methods for automated segmentation of subnuclear chromatin domains Epigenomes 2022 6 34 10.3390/epigenomes6040034 36278680
Johann to Berens P, Schivre G, Theune M et al (2022) Advanced image analysis methods for automated segmentation of subnuclear chromatin domains. Epigenomes 6:34. 10.3390/epigenomes604003436278680 10.3390/epigenomes6040034
Kenesi E Kolbert Z Kaszler N The ROP2 GTPase participates in nitric oxide (NO)-induced root shortening in Arabidopsis Plants 2023 12 750 10.3390/plants12040750 36840099
Kenesi E, Kolbert Z, Kaszler N et al (2023) The ROP2 GTPase participates in nitric oxide (NO)-induced root shortening in Arabidopsis. Plants 12:750. 10.3390/plants1204075036840099 10.3390/plants12040750
Kimata Y Kato T Higaki T Polar vacuolar distribution is essential for accurate asymmetric division of Arabidopsis zygotes Proc Natl Acad Sci U S A 2019 116 2338 2343 10.1073/pnas.1814160116 30651313
Kimata Y, Kato T, Higaki T et al (2019) Polar vacuolar distribution is essential for accurate asymmetric division of Arabidopsis zygotes. Proc Natl Acad Sci U S A 116:2338–2343. 10.1073/pnas.181416011630651313 10.1073/pnas.1814160116
Lavrekha VV Pasternak T Ivanov VB 3D analysis of mitosis distribution highlights the longitudinal zonation and diarch symmetry in proliferation activity of the Arabidopsis thaliana root meristem Plant J 2017 92 834 845 10.1111/TPJ.13720 28921702
Lavrekha VV, Pasternak T, Ivanov VB et al (2017) 3D analysis of mitosis distribution highlights the longitudinal zonation and diarch symmetry in proliferation activity of the Arabidopsis thaliana root meristem. Plant J 92:834–845. 10.1111/TPJ.1372028921702 10.1111/TPJ.13720
Liu F Wang B Wang X FAN, the homolog of mammalian apoptosis antagonizing transcription factor AATF/Che-1 protein, is involved in safeguarding genome stability through the ATR induced pathway in Arabidopsis bioRxiv 2023 10.1101/2023.03.10.531693 38405786
Liu F, Wang B, Wang X et al (2023) FAN, the homolog of mammalian apoptosis antagonizing transcription factor AATF/Che-1 protein, is involved in safeguarding genome stability through the ATR induced pathway in Arabidopsis. bioRxiv. 10.1101/2023.03.10.53169338405786 10.1101/2023.03.10.531693
Mattout A Cabianca DS Gasser SM Chromatin states and nuclear organization in development—a view from the nuclear lamina Genome Biol 2015 16 174 10.1186/s13059-015-0747-5 26303512
Mattout A, Cabianca DS, Gasser SM (2015) Chromatin states and nuclear organization in development—a view from the nuclear lamina. Genome Biol 16:174. 10.1186/s13059-015-0747-526303512 10.1186/s13059-015-0747-5
Misteli T The self-organizing genome: principles of genome architecture and function Cell 2020 183 28 45 10.1016/J.CELL.2020.09.014 32976797
Misteli T (2020) The self-organizing genome: principles of genome architecture and function. Cell 183:28–45. 10.1016/J.CELL.2020.09.01432976797 10.1016/J.CELL.2020.09.014
Nakamura M Claes AR Grebe T Auxin and ROP GTPase signaling of polar nuclear migration in root epidermal hair cells Plant Physiol 2018 176 378 391 10.1104/PP.17.00713 29084900
Nakamura M, Claes AR, Grebe T et al (2018) Auxin and ROP GTPase signaling of polar nuclear migration in root epidermal hair cells. Plant Physiol 176:378–391. 10.1104/PP.17.0071329084900 10.1104/PP.17.00713
Pasternak T Pérez-Pérez JM Methods of in situ quantitative root biology Plants 2021 10 2399 10.3390/PLANTS10112399 34834762
Pasternak T, Pérez-Pérez JM (2021) Methods of in situ quantitative root biology. Plants 10:2399. 10.3390/PLANTS1011239934834762 10.3390/PLANTS10112399
Pasternak T Tietz O Rapp K Protocol: an improved and universal procedure for whole-mount immunolocalization in plants Plant Methods 2015 11 1 10 10.1186/s13007-015-0094-2 25649124
Pasternak T, Tietz O, Rapp K et al (2015) Protocol: an improved and universal procedure for whole-mount immunolocalization in plants. Plant Methods 11:1–10. 10.1186/s13007-015-0094-225649124 10.1186/s13007-015-0094-2
Pasternak T Paponov IA Kondratenko S Optimizing protocols for Arabidopsis shoot and root protoplast cultivation Plants 2021 10 375 10.3390/PLANTS10020375 33672063
Pasternak T, Paponov IA, Kondratenko S (2021) Optimizing protocols for Arabidopsis shoot and root protoplast cultivation. Plants 10:375. 10.3390/PLANTS1002037533672063 10.3390/PLANTS10020375
Pasternak T Kircher S Pérez-Pérez JM Palme K A simple pipeline for cell cycle kinetic studies in the root apical meristem J Exp Bot 2022 73 4683 4695 10.1093/JXB/ERAC123 35312781
Pasternak T, Kircher S, Pérez-Pérez JM, Palme K (2022) A simple pipeline for cell cycle kinetic studies in the root apical meristem. J Exp Bot 73:4683–4695. 10.1093/JXB/ERAC12335312781 10.1093/JXB/ERAC123
Pasternak T Pérez-Pérez JM Ruperti B A new in vitro growth system for phenotypic characterization and seed propagation of Arabidopsis thaliana J Plant Growth Regul 2023 43 652 658 10.1007/s00344-023-11093-x
Pasternak T, Pérez-Pérez JM, Ruperti B et al (2023) A new in vitro growth system for phenotypic characterization and seed propagation of Arabidopsis thaliana. J Plant Growth Regul 43:652–658. 10.1007/s00344-023-11093-x10.1007/s00344-023-11093-x
Pasternak T, Falk T, Paponov I (2024) Deep-resolution plant phenotyping platform description. In: Protocols.Io. https://www.protocols.io/view/deep-resolution-plant-phenotyping-platform-descrip-8epv5z2y4v1b/v1. Accessed 20 Feb 2024
Pavlova P van Zanten M Snoek BL 2D morphometric analysis of Arabidopsis thaliana nuclei reveals characteristic profiles of different cell types and accessions Chromosom Res 2022 30 5 24 10.1007/S10577-021-09673-2
Pavlova P, van Zanten M, Snoek BL et al (2022) 2D morphometric analysis of Arabidopsis thaliana nuclei reveals characteristic profiles of different cell types and accessions. Chromosom Res 30:5–24. 10.1007/S10577-021-09673-210.1007/S10577-021-09673-2
Pei L Li G Lindsey K Plant 3D genomics: the exploration and application of chromatin organization New Phytol 2021 230 1772 1786 10.1111/NPH.17262 33560539
Pei L, Li G, Lindsey K et al (2021) Plant 3D genomics: the exploration and application of chromatin organization. New Phytol 230:1772–1786. 10.1111/NPH.1726233560539 10.1111/NPH.17262
Poleshko A Smith CL Nguyen SC H3k9me2 orchestrates inheritance of spatial positioning of peripheral heterochromatin through mitosis Elife 2019 8 e49278 10.7554/ELIFE.49278 31573510
Poleshko A, Smith CL, Nguyen SC et al (2019) H3k9me2 orchestrates inheritance of spatial positioning of peripheral heterochromatin through mitosis. Elife 8:e49278. 10.7554/ELIFE.4927831573510 10.7554/ELIFE.49278
Pontvianne F Grob S Three-dimensional nuclear organization in Arabidopsis thaliana J Plant Res 2020 133 479 488 10.1007/S10265-020-01185-0 32240449
Pontvianne F, Grob S (2020) Three-dimensional nuclear organization in Arabidopsis thaliana. J Plant Res 133:479–488. 10.1007/S10265-020-01185-032240449 10.1007/S10265-020-01185-0
Ramakrishnan Chandra J Kalidass M Demidov D The role of centromeric repeats and transcripts in kinetochore assembly and function Plant J 2024 118 982 996 10.1111/TPJ.16445 37665331
Ramakrishnan Chandra J, Kalidass M, Demidov D et al (2024) The role of centromeric repeats and transcripts in kinetochore assembly and function. Plant J 118:982–996. 10.1111/TPJ.1644537665331 10.1111/TPJ.16445
Randall RS Jourdain C Nowicka A Image analysis workflows to reveal the spatial organization of cell nuclei and chromosomes Nucleus 2022 13 277 299 10.1080/19491034.2022.2144013 36447428
Randall RS, Jourdain C, Nowicka A et al (2022) Image analysis workflows to reveal the spatial organization of cell nuclei and chromosomes. Nucleus 13:277–299. 10.1080/19491034.2022.214401336447428 10.1080/19491034.2022.2144013
Robert HS Grones P Stepanova AN Local auxin sources orient the apical-basal axis in arabidopsis embryos Curr Biol 2013 23 2506 2512 10.1016/j.cub.2013.09.039 24291089
Robert HS, Grones P, Stepanova AN et al (2013) Local auxin sources orient the apical-basal axis in arabidopsis embryos. Curr Biol 23:2506–2512. 10.1016/j.cub.2013.09.03924291089 10.1016/j.cub.2013.09.039
Roudier F Ahmed I Bérard C Integrative epigenomic mapping defines four main chromatin states in Arabidopsis EMBO J 2011 30 1928 1938 10.1038/EMBOJ.2011.103 21487388
Roudier F, Ahmed I, Bérard C et al (2011) Integrative epigenomic mapping defines four main chromatin states in Arabidopsis. EMBO J 30:1928–1938. 10.1038/EMBOJ.2011.10321487388 10.1038/EMBOJ.2011.103
Schindelin J Arganda-Carreras I Frise E Fiji: an open-source platform for biological-image analysis Nat Methods 2012 9 676 682 10.1038/nmeth.2019 22743772
Schindelin J, Arganda-Carreras I, Frise E et al (2012) Fiji: an open-source platform for biological-image analysis. Nat Methods 9:676–682. 10.1038/nmeth.201922743772 10.1038/nmeth.2019
Schmidt T Pasternak T Liu K The iRoCS toolbox—3D analysis of the plant root apical meristem at cellular resolution Plant J 2014 77 806 814 10.1111/TPJ.12429 24417645
Schmidt T, Pasternak T, Liu K et al (2014) The iRoCS toolbox—3D analysis of the plant root apical meristem at cellular resolution. Plant J 77:806–814. 10.1111/TPJ.1242924417645 10.1111/TPJ.12429
Tatout C Mougeot G Parry G The INDEPTH (impact of nuclear domains on gene expression and plant traits) academy: a community resource for plant science J Exp Bot 2022 73 1926 1933 10.1093/JXB/ERAC005 35090020
Tatout C, Mougeot G, Parry G et al (2022) The INDEPTH (impact of nuclear domains on gene expression and plant traits) academy: a community resource for plant science. J Exp Bot 73:1926–1933. 10.1093/JXB/ERAC00535090020 10.1093/JXB/ERAC005
Tessarz P Kouzarides T Histone core modifications regulating nucleosome structure and dynamics Nat Rev Mol Cell Biol 2014 15 703 708 10.1038/nrm3890 25315270
Tessarz P, Kouzarides T (2014) Histone core modifications regulating nucleosome structure and dynamics. Nat Rev Mol Cell Biol 15:703–708. 10.1038/nrm389025315270 10.1038/nrm3890
Vergara Z Gomez MS Desvoyes B Distinct roles of Arabidopsis ORC1 proteins in DNA replication and heterochromatic H3K27me1 deposition Nat Commun 2023 14 1 13 10.1038/s41467-023-37024-8 36596776
Vergara Z, Gomez MS, Desvoyes B et al (2023) Distinct roles of Arabidopsis ORC1 proteins in DNA replication and heterochromatic H3K27me1 deposition. Nat Commun 14:1–13. 10.1038/s41467-023-37024-836596776 10.1038/s41467-023-37024-8
Wada M Nuclear movement and positioning in plant cells Semin Cell Dev Biol 2018 82 17 24 10.1016/J.SEMCDB.2017.10.001 29054391
Wada M (2018) Nuclear movement and positioning in plant cells. Semin Cell Dev Biol 82:17–24. 10.1016/J.SEMCDB.2017.10.00129054391 10.1016/J.SEMCDB.2017.10.001
Wang C Liu C Roqueiro D Genome-wide analysis of local chromatin packing in Arabidopsis thaliana Genome Res 2015 25 246 256 10.1101/GR.170332.113 25367294
Wang C, Liu C, Roqueiro D et al (2015) Genome-wide analysis of local chromatin packing in Arabidopsis thaliana. Genome Res 25:246–256. 10.1101/GR.170332.11325367294 10.1101/GR.170332.113
Yamada M Goshima G Mitotic spindle assembly in land plants: molecules and mechanisms Biol 2017 6 6 10.3390/BIOLOGY6010006
Yamada M, Goshima G (2017) Mitotic spindle assembly in land plants: molecules and mechanisms. Biol 6:6. 10.3390/BIOLOGY601000610.3390/BIOLOGY6010006
Zhang X Bernatavichute YV Cokus S Genome-wide analysis of mono-, di- and trimethylation of histone H3 lysine 4 in Arabidopsis thaliana Genome Biol 2009 10 R62 10.1186/GB-2009-10-6-R62 19508735
Zhang X, Bernatavichute YV, Cokus S et al (2009) Genome-wide analysis of mono-, di- and trimethylation of histone H3 lysine 4 in Arabidopsis thaliana. Genome Biol 10:R62. 10.1186/GB-2009-10-6-R6219508735 10.1186/GB-2009-10-6-R62
Zhang Y Dong Q Wang Z Liu Q Yu H Sun W Zhang H A fine-scale Arabidopsis chromatin landscape reveals chromatin conformation-associated transcriptional dynamics Nat Commun 2024 15 3253 10.1038/s41467-024-47678-7 38627396
Zhang Y, Dong Q, Wang Z, Liu Q, Yu H, Sun W, Zhang H (2024) A fine-scale Arabidopsis chromatin landscape reveals chromatin conformation-associated transcriptional dynamics. Nat Commun 15:3253. 10.1038/s41467-024-47678-738627396 10.1038/s41467-024-47678-7
