Cistrome and Epicistrome Features Shape the Regulatory DNA Landscape.

The cistrome is the complete set of transcription factor (TF) binding sites (cis-elements) in an organism, while an epicistrome incorporates tissue-specific DNA chemical modifications and TF-specific chemical sensitivities into these binding profiles. Robust methods to construct comprehensive cistrome and epicistrome maps are critical for elucidating complex transcriptional networks that underlie growth, behavior, and disease. Here, we describe DNA affinity purification sequencing (DAP-seq), a high-throughput TF binding site discovery method that interrogates genomic DNA with in-vitro-expressed TFs. Using DAP-seq, we defined the Arabidopsis cistrome by resolving motifs and peaks for 529 TFs. Because genomic DNA used in DAP-seq retains 5-methylcytosines, we determined that >75% (248/327) of Arabidopsis TFs surveyed were methylation sensitive, a property that strongly impacts the epicistrome landscape. DAP-seq datasets also yielded insight into the biology and binding site architecture of numerous TFs, demonstrating the value of DAP-seq for cost-effective cistromic and epicistromic annotation in any organism.

Dynamics and function of DNA methylation in plants.

DNA methylation is a conserved epigenetic modification that is important for gene regulation and genome stability. Aberrant patterns of DNA methylation can lead to plant developmental abnormalities. A specific DNA methylation state is an outcome of dynamic regulation by de novo methylation, maintenance of methylation and active demethylation, which are catalysed by various enzymes that are targeted by distinct regulatory pathways. In this Review, we discuss DNA methylation in plants, including methylating and demethylating enzymes and regulatory factors, and the coordination of methylation and demethylation activities by a so-called methylstat mechanism; the functions of DNA methylation in regulating transposon silencing, gene expression and chromosome interactions; the roles of DNA methylation in plant development; and the involvement of DNA methylation in plant responses to biotic and abiotic stress conditions.

The complex polyploid genome architecture of sugarcane.

Sugarcane, the world's most harvested crop by tonnage, has shaped global history, trade and geopolitics, and is currently responsible for 80% of sugar production worldwide1. While traditional sugarcane breeding methods have effectively generated cultivars adapted to new environments and pathogens, sugar yield improvements have recently plateaued2. The cessation of yield gains may be due to limited genetic diversity within breeding populations, long breeding cycles and the complexity of its genome, the latter preventing breeders from taking advantage of the recent explosion of whole-genome sequencing that has benefited many other crops. Thus, modern sugarcane hybrids are the last remaining major crop without a reference-quality genome. Here we take a major step towards advancing sugarcane biotechnology by generating a polyploid reference genome for R570, a typical modern cultivar derived from interspecific hybridization between the domesticated species (Saccharum officinarum) and the wild species (Saccharum spontaneum). In contrast to the existing single haplotype ('monoploid') representation of R570, our 8.7 billion base assembly contains a complete representation of unique DNA sequences across the approximately 12 chromosome copies in this polyploid genome. Using this highly contiguous genome assembly, we filled a previously unsized gap within an R570 physical genetic map to describe the likely causal genes underlying the single-copy Bru1 brown rust resistance locus. This polyploid genome assembly with fine-grain descriptions of genome architecture and molecular targets for biotechnology will help accelerate molecular and transgenic breeding and adaptation of sugarcane to future environmental conditions.

Molecular mechanisms of the RNA polymerases in plant RNA-directed DNA methylation.

In plants, two atypical DNA-dependent RNA polymerases, RNA polymerase IV (Pol IV) and Pol V, and an RNA-DEPENDENT RNA POLYMERASE 2 (RDR2) together produce noncoding RNAs (ncRNAs) to guide the plant-specific RNA-directed DNA methylation (RdDM). Although both Pol IV and Pol V have evolved from the canonical Pol II, they have adapted to different roles in RdDM. The mechanisms of their adaptation are key to understanding plant DNA methylation and the divergent evolution of polymerases. In this review, we summarize insights that have emerged from recent structural studies of Pol IV, Pol V, and RDR2 and discuss their structural features critical for efficient ncRNA production in RdDM.

Brassinosteroid-dependent phosphorylation of PHOSPHATE STARVATION RESPONSE2 reduces its DNA-binding ability in rice.

Brassinosteroids (BRs) are widely used as plant growth regulators in modern agriculture. Understanding how BRs regulate nutrient signaling is crucial for reducing fertilizer usage. Here we elucidate that the central BR signaling inhibitor GSK3/SHAGGY-LIKE KINASE2 (GSK2) interacts directly with and phosphorylates PHOSPHATE STARVATION RESPONSE2 (OsPHR2), the key regulator of phosphate (Pi) signaling, to suppress its transcription factor activity in rice (Oryza sativa). We identify a critical phosphorylation site at serine residue S269 of OsPHR2 and demonstrate that phosphorylation by GSK2 or phosphor-mimic mutation of S269 substantially impairs the DNA-binding activity of OsPHR2, and thus diminishes expression of OsPHR2-induced genes and reduces Pi levels. Like BRs, Pi starvation noticeably induces GSK2 instability. We further show that this site-specific phosphorylation event is conserved in Arabidopsis (Arabidopsis thaliana), but varies among the PHR-family members, being present only in most land plants. These results unveil a distinctive post-transcriptional regulatory mechanism in Pi signaling by which BRs promote Pi acquisition, with a potential contribution to the environmental adaptability of plants during their evolution.

Dual binding of chromomethylase domains to H3K9me2-containing nucleosomes directs DNA methylation in plants.

DNA methylation and histone modification exert epigenetic control over gene expression. CHG methylation by CHROMOMETHYLASE3 (CMT3) depends on histone H3K9 dimethylation (H3K9me2), but the mechanism underlying this relationship is poorly understood. Here, we report multiple lines of evidence that CMT3 interacts with H3K9me2-containing nucleosomes. CMT3 genome locations nearly perfectly correlated with H3K9me2, and CMT3 stably associated with H3K9me2-containing nucleosomes. Crystal structures of maize CMT3 homolog ZMET2, in complex with H3K9me2 peptides, showed that ZMET2 binds H3K9me2 via both bromo adjacent homology (BAH) and chromo domains. The structures reveal an aromatic cage within both BAH and chromo domains as interaction interfaces that capture H3K9me2. Mutations that abolish either interaction disrupt CMT3 binding to nucleosomes and show a complete loss of CMT3 activity in vivo. Our study establishes dual recognition of H3K9me2 marks by BAH and chromo domains and reveals a distinct mechanism of interplay between DNA methylation and histone modification.

Reaction Mechanisms of Pol IV, RDR2, and DCL3 Drive RNA Channeling in the siRNA-Directed DNA Methylation Pathway.

In eukaryotes with multiple small RNA pathways, the mechanisms that channel RNAs within specific pathways are unclear. Here, we reveal the reactions that account for channeling in the small interfering RNA (siRNA) biogenesis phase of the Arabidopsis RNA-directed DNA methylation pathway. The process begins with template DNA transcription by NUCLEAR RNA POLYMERASE IV (Pol IV), whose atypical termination mechanism, induced by nontemplate DNA base-pairing, channels transcripts to the associated RNA-dependent RNA polymerase RDR2. RDR2 converts Pol IV transcripts into double-stranded RNAs and then typically adds an extra untemplated 3' terminal nucleotide to the second strands. The dicer endonuclease DCL3 cuts resulting duplexes to generate 24- and 23-nt siRNAs. The 23-nt RNAs bear the untemplated terminal nucleotide of the RDR2 strand and are underrepresented among ARGONAUTE4-associated siRNAs. Collectively, our results provide mechanistic insights into Pol IV termination, Pol IV-RDR2 coupling, and RNA channeling, from template DNA transcription to siRNA strand discrimination.

Transgenerational increases in DNA methylation in Arabidopsis plants defective in active DNA demethylation.

Spontaneous gain or loss of DNA methylation occurs in plant and animal genomes, and DNA methylation changes can lead to meiotically stable epialleles that generate heritable phenotypic diversity. However, it is unclear whether transgenerational epigenetic stability may be regulated by any cellular factors. Here, we examined spontaneously occurring variations in DNA methylation in wild-type and ros1 mutant Arabidopsis plants that were propagated for ten generations from single-seed descent. We found that the ros1 mutant, which is defective in active DNA demethylation, showed an increased transgenerational epimutation rate. The ros1 mutation led to more spontaneously gained methylation than lost methylation at individual cytosines, compared to the wild type which had similar numbers of spontaneously gained and lost methylation cytosines. Consistently, transgenerational differentially methylated regions were also biased toward hypermethylation in the ros1 mutant. Our results reveal a genetic contribution of the ROS1 DNA demethylase to transgenerational epigenetic stability and suggest that ROS1 may have an unexpected surveillance function in preventing transgenerational DNA methylation increases.

Genome copy number predicts extreme evolutionary rate variation in plant mitochondrial DNA.

Nuclear and organellar genomes can evolve at vastly different rates despite occupying the same cell. In most bilaterian animals, mitochondrial DNA (mtDNA) evolves faster than nuclear DNA, whereas this trend is generally reversed in plants. However, in some exceptional angiosperm clades, mtDNA substitution rates have increased up to 5,000-fold compared with closely related lineages. The mechanisms responsible for this acceleration are generally unknown. Because plants rely on homologous recombination to repair mtDNA damage, we hypothesized that mtDNA copy numbers may predict evolutionary rates, as lower copy numbers may provide fewer templates for such repair mechanisms. In support of this hypothesis, we found that copy number explains 47% of the variation in synonymous substitution rates of mtDNA across 60 diverse seed plant species representing ~300 million years of evolution. Copy number was also negatively correlated with mitogenome size, which may be a cause or consequence of mutation rate variation. Both relationships were unique to mtDNA and not observed in plastid DNA. These results suggest that homologous recombinational repair plays a role in driving mtDNA substitution rates in plants and may explain variation in mtDNA evolution more broadly across eukaryotes. Our findings also contribute to broader questions about the relationships between mutation rates, genome size, selection efficiency, and the drift-barrier hypothesis.

Simultaneous profiling of chromatin accessibility and DNA methylation in complete plant genomes using long-read sequencing.

Epigenetic regulations, including chromatin accessibility, nucleosome positioning and DNA methylation intricately shape genome function. However, current chromatin profiling techniques relying on short-read sequencing technologies fail to characterise highly repetitive genomic regions and cannot detect multiple chromatin features simultaneously. Here, we performed Simultaneous Accessibility and DNA Methylation Sequencing (SAM-seq) of purified plant nuclei. Thanks to the use of long-read nanopore sequencing, SAM-seq enables high-resolution profiling of m6A-tagged chromatin accessibility together with endogenous cytosine methylation in plants. Analysis of naked genomic DNA revealed significant sequence preference biases of m6A-MTases, controllable through a normalisation step. By applying SAM-seq to Arabidopsis and maize nuclei we obtained fine-grained accessibility and DNA methylation landscapes genome-wide. We uncovered crosstalk between chromatin accessibility and DNA methylation within nucleosomes of genes, TEs, and centromeric repeats. SAM-seq also detects DNA footprints over cis-regulatory regions. Furthermore, using the single-molecule information provided by SAM-seq we identified extensive cellular heterogeneity at chromatin domains with antagonistic chromatin marks, suggesting that bivalency reflects cell-specific regulations. SAM-seq is a powerful approach to simultaneously study multiple epigenetic features over unique and repetitive sequences, opening new opportunities for the investigation of epigenetic mechanisms.

The clade-specific target recognition mechanisms of plant RISCs.

Eukaryotic Argonaut proteins (AGOs) assemble RNA-induced silencing complexes (RISCs) with guide RNAs that allow binding to complementary RNA sequences and subsequent silencing of target genes. The model plant Arabidopsis thaliana encodes 10 different AGOs, categorized into three distinct clades based on amino acid sequence similarity. While clade 1 and 2 RISCs are known for their roles in post-transcriptional gene silencing, and clade 3 RISCs are associated with transcriptional gene silencing in the nucleus, the specific mechanisms of how RISCs from each clade recognize their targets remain unclear. In this study, I conducted quantitative binding analyses between RISCs and target nucleic acids with mismatches at various positions, unveiling distinct target binding characteristics unique to each clade. Clade 1 and 2 RISCs require base pairing not only in the seed region but also in the 3' supplementary region for stable target RNA binding, with clade 1 exhibiting a higher stringency. Conversely, clade 3 RISCs tolerate dinucleotide mismatches beyond the seed region. Strikingly, they bind to DNA targets with an affinity equal to or surpassing that of RNA, like prokaryotic AGO complexes. These insights challenge existing views on plant RNA silencing and open avenues for exploring new functions of eukaryotic AGOs.

Diversity of plant DNA in stool is linked to dietary quality, age, and household income.

Eating a varied diet is a central tenet of good nutrition. Here, we develop a molecular tool to quantify human dietary plant diversity by applying DNA metabarcoding with the chloroplast trnL-P6 marker to 1,029 fecal samples from 324 participants across two interventional feeding studies and three observational cohorts. The number of plant taxa per sample (plant metabarcoding richness or pMR) correlated with recorded intakes in interventional diets and with indices calculated from a food frequency questionnaire in typical diets (ρ = 0.40 to 0.63). In adolescents unable to collect validated dietary survey data, trnL metabarcoding detected 111 plant taxa, with 86 consumed by more than one individual and four (wheat, chocolate, corn, and potato family) consumed by >70% of individuals. Adolescent pMR was associated with age and household income, replicating prior epidemiologic findings. Overall, trnL metabarcoding promises an objective and accurate measure of the number and types of plants consumed that is applicable to diverse human populations.

The interplay of homology-directed repair pathways in the repair of zebularine-induced DNA-protein crosslinks in Arabidopsis.

DNA-protein crosslinks (DPCs) are highly toxic DNA lesions represented by proteins covalently bound to the DNA. Persisting DPCs interfere with fundamental genetic processes such as DNA replication and transcription. Cytidine analog zebularine (ZEB) has been shown to crosslink DNA METHYLTRANSFERASE1 (MET1). Recently, we uncovered a critical role of the SMC5/6-mediated SUMOylation in the repair of DPCs. In an ongoing genetic screen, we identified two additional candidates, HYPERSENSITIVE TO ZEBULARINE 2 and 3, that were mapped to REGULATOR OF TELOMERE ELONGATION 1 (RTEL1) and polymerase TEBICHI (TEB), respectively. By monitoring the growth of hze2 and hze3 plants in response to zebularine, we show the importance of homologous recombination (HR) factor RTEL1 and microhomology-mediated end-joining (MMEJ) polymerase TEB in the repair of MET1-DPCs. Moreover, genetic interaction and sensitivity assays showed the interdependency of SMC5/6 complex, HR, and MMEJ in the homology-directed repair of MET1-DPCs in Arabidopsis. Altogether, we provide evidence that MET1-DPC repair in plants is more complex than originally expected.

How diverse a monocentric chromosome can be? Repeatome and centromeric organization of Juncus effusus (Juncaceae).

Juncus is the largest genus of Juncaceae and was considered holocentric for a long time. Recent findings, however, indicated that 11 species from different clades of the genus have monocentric chromosomes. Thus, the Juncus centromere organization and evolution need to be reassessed. We aimed to investigate the major repetitive DNA sequences of two accessions of Juncus effusus and its centromeric structure by employing whole-genome analyses, fluorescent in situ hybridization, CENH3 immunodetection, and chromatin immunoprecipitation sequencing. We showed that the repetitive fraction of the small J. effusus genome (~270 Mbp/1C) is mainly composed of Class I and Class II transposable elements (TEs) and satellite DNAs. Three identified satellite DNA families were mainly (peri)centromeric, with two being associated with the centromeric protein CENH3, but not strictly centromeric. Two types of centromere organization were discerned in J. effusus: type 1 was characterized by a single CENH3 domain enriched with JefSAT1-155 or JefSAT2-180, whereas type 2 showed multiple CENH3 domains interrupted by other satellites, TEs or genes. Furthermore, while type 1 centromeres showed a higher degree of satellite identity along the array, type 2 centromeres had less homogenized arrays along the multiple CENH3 domains per chromosome. Although the analyses confirmed the monocentric organization of J. effusus chromosomes, our data indicate a more dynamic arrangement of J. effusus centromeres than observed for other plant species, suggesting it may constitute a transient state between mono- and holocentricity.

Degenerate oligonucleotide primer MIG-seq: an effective PCR-based method for high-throughput genotyping.

Next-generation sequencing (NGS) library construction often involves using restriction enzymes to decrease genome complexity, enabling versatile polymorphism detection in plants. However, plant leaves frequently contain impurities, such as polyphenols, necessitating DNA purification before enzymatic reactions. To overcome this problem, we developed a PCR-based method for expeditious NGS library preparation, offering flexibility in number of detected polymorphisms. By substituting a segment of the simple sequence repeat sequence in the MIG-seq primer set (MIG-seq being a PCR method enabling library construction with low-quality DNA) with degenerate oligonucleotides, we introduced variability in detectable polymorphisms across various crops. This innovation, named degenerate oligonucleotide primer MIG-seq (dpMIG-seq), enabled a streamlined protocol for constructing dpMIG-seq libraries from unpurified DNA, which was implemented stably in several crop species, including fruit trees. Furthermore, dpMIG-seq facilitated efficient lineage selection in wheat and enabled linkage map construction and quantitative trait loci analysis in tomato, rice, and soybean without necessitating DNA concentration adjustments. These findings underscore the potential of the dpMIG-seq protocol for advancing genetic analyses across diverse plant species.

Population epigenetics: DNA methylation in the plant omics era.

DNA methylation plays an important role in many biological processes. The mechanisms underlying the establishment and maintenance of DNA methylation are well understood thanks to decades of research using DNA methylation mutants, primarily in Arabidopsis (Arabidopsis thaliana) accession Col-0. Recent genome-wide association studies (GWASs) using the methylomes of natural accessions have uncovered a complex and distinct genetic basis of variation in DNA methylation at the population level. Sequencing following bisulfite treatment has served as an excellent method for quantifying DNA methylation. Unlike studies focusing on specific accessions with reference genomes, population-scale methylome research often requires an additional round of sequencing beyond obtaining genome assemblies or genetic variations from whole-genome sequencing data, which can be cost prohibitive. Here, we provide an overview of recently developed bisulfite-free methods for quantifying methylation and cost-effective approaches for the simultaneous detection of genetic and epigenetic information. We also discuss the plasticity of DNA methylation in a specific Arabidopsis accession, the contribution of DNA methylation to plant adaptation, and the genetic determinants of variation in DNA methylation in natural populations. The recently developed technology and knowledge will greatly benefit future studies in population epigenomes.

The Arabidopsis ATR-SOG1 signaling module regulates pleiotropic developmental adjustments in response to 3'-blocked DNA repair intermediates.

Base excision repair and active DNA demethylation produce repair intermediates with DNA molecules blocked at the 3'-OH end by an aldehyde or phosphate group. However, both the physiological consequences of these accumulated single-strand DNAs break with 3'-blocked ends (DNA 3'-blocks) and the signaling pathways responding to unrepaired DNA 3'-blocks remain unclear in plants. Here, we investigated the effects of DNA 3'-blocks on plant development using the zinc finger DNA 3'-phosphoesterase (zdp) AP endonuclease2 (ape2) double mutant, in which 3'-blocking residues are poorly repaired. The accumulation of DNA 3'-blocked triggered diverse developmental defects that were dependent on the ATM and RAD3-related (ATR)-suppressor of gamma response 1 (SOG1) signaling module. SOG1 mutation rescued the developmental defects of zdp ape2 leaves by preventing cell endoreplication and promoting cell proliferation. However, SOG1 mutation caused intensive meristematic cell death in the radicle of zdp ape2 following germination, resulting in rapid termination of radicle growth. Notably, mutating FORMAMIDOPYRIMIDINE DNA GLYCOSYLASE (FPG) in zdp ape2 sog1 partially recovered its radicle growth, demonstrating that DNA 3'-blocks generated by FPG caused the meristematic defects. Surprisingly, despite lacking a functional radicle, zdp ape2 sog1 mutants compensated the lack of root growth by generating anchor roots having low levels of DNA damage response. Our results reveal dual roles of SOG1 in regulating root establishment when seeds germinate with excess DNA 3'-blocks.

A diRNA-protein scaffold module mediates SMC5/6 recruitment in plant DNA repair.

In eukaryotes, the STRUCTURAL MAINTENANCE OF CHROMOSOME 5/6 (SMC5/6) complex is critical to maintaining chromosomal structures around double-strand breaks (DSBs) in DNA damage repair. However, the recruitment mechanism of this conserved complex at DSBs remains unclear. In this study, using Arabidopsis thaliana as a model, we found that SMC5/6 localization at DSBs is dependent on the protein scaffold containing INVOLVED IN DE NOVO 2 (IDN2), CELL DIVISION CYCLE 5 (CDC5), and ALTERATION/DEFICIENCY IN ACTIVATION 2B (ADA2b), whose recruitment is further mediated by DNA-damage-induced RNAs (diRNAs) generated from DNA regions around DSBs. The physical interactions of protein components including SMC5-ADA2b, ADA2b-CDC5, and CDC5-IDN2 result in formation of the protein scaffold. Further analysis indicated that the DSB localization of IDN2 requires its RNA-binding activity and ARGONAUTE 2 (AGO2), indicating a role for the AGO2-diRNA complex in this process. Given that most of the components in the scaffold are conserved, the mechanism presented here, which connects SMC5/6 recruitment and small RNAs, will improve our understanding of DNA repair mechanisms in eukaryotes.

H2A.Z-containing nucleosomes mediate the thermosensory response in Arabidopsis.

Plants are highly sensitive to temperature and can perceive a difference of as little as 1 degrees C. How temperature is sensed and integrated in development is unknown. In a forward genetic screen in Arabidopsis, we have found that nucleosomes containing the alternative histone H2A.Z are essential to perceiving ambient temperature correctly. Genotypes deficient in incorporating H2A.Z into nucleosomes phenocopy warm grown plants, and show a striking constitutive warm temperature transcriptome. We show that nucleosomes containing H2A.Z display distinct responses to temperature in vivo, independently of transcription. Using purified nucleosomes, we are able to show that H2A.Z confers distinct DNA-unwrapping properties on nucleosomes, indicating a direct mechanism for the perception of temperature through DNA-nucleosome fluctuations. Our results show that H2A.Z-containing nucleosomes provide thermosensory information that is used to coordinate the ambient temperature transcriptome. We observe the same effect in budding yeast, indicating that this is an evolutionarily conserved mechanism.

Evolutionarily Conserved Principles Predict 3D Chromatin Organization.

Topologically associating domains (TADs), CTCF loop domains, and A/B compartments have been identified as important structural and functional components of 3D chromatin organization, yet the relationship between these features is not well understood. Using high-resolution Hi-C and HiChIP, we show that Drosophila chromatin is organized into domains we term compartmental domains that correspond precisely with A/B compartments at high resolution. We find that transcriptional state is a major predictor of Hi-C contact maps in several eukaryotes tested, including C. elegans and A. thaliana. Architectural proteins insulate compartmental domains by reducing interaction frequencies between neighboring regions in Drosophila, but CTCF loops do not play a distinct role in this organism. In mammals, compartmental domains exist alongside CTCF loop domains to form topological domains. The results suggest that compartmental domains are responsible for domain structure in all eukaryotes, with CTCF playing an important role in domain formation in mammals.