Characterization of functional relationships of R-loops with gene transcription and epigenetic modifications in rice

Non-B DNA in plant genomes: prediction, mapping, and emerging roles

Regulating gene expression in plant development and environmental responses is vital for mitigating the effects of climate change on crop growth and productivity. The eukaryotic genome largely shows the canonical B-DNA structure that is organized into nucleosomes with histone modifications shaping the epigenome. Nuclear proteins and RNA interactions influence chromatin conformations and dynamically modulate gene activity. Non-B DNA conformations and their transitions introduce novel aspects to gene expression modulation, particularly in response to environmental shifts. We explore the current understanding of non-B DNA structures in plant genomes, their interplay with epigenomics and gene expression, and advances in methods for their mapping and characterization. The exploration of so far uncharacterized non-B DNA structures remains an intriguing area in plant chromatin research and offers insights into their potential role in gene regulation.

Genome-wide characterization of single-stranded DNA in rice

Single-stranded DNA (ssDNA) is essential for various DNA-templated processes in both eukaryotes and prokaryotes. However, comprehensive characterizations of ssDNA still lag in plants compared to nonplant systems. Here, we conducted in situ S1-sequencing, with starting gDNA ranging from 5 µg to 250 ng, followed by comprehensive characterizations of ssDNA in rice (Oryza sativa L.). We found that ssDNA loci were substantially associated with a subset of non-B DNA structures and functional genomic loci. Subtypes of ssDNA loci had distinct epigenetic features. Importantly, ssDNA may act alone or partly coordinate with non-B DNA structures, functional genomic loci, or epigenetic marks to actively or repressively modulate gene transcription, which is genomic region dependent and associated with the distinct accumulation of RNA Pol II. Moreover, distinct types of ssDNA had differential impacts on the activities and evolution of transposable elements (TEs) (especially common or conserved TEs) in the rice genome. Our study showcases an antibody-independent technique for characterizing non-B DNA structures or functional genomic loci in plants. It lays the groundwork and fills a crucial gap for further exploration of ssDNA, non-B DNA structures, or functional genomic loci, thereby advancing our understanding of their biology in plants.

Unveiling the distinctive traits of functional rye centromeres: minisatellites, retrotransposons, and R-loop formation

Centromeres play a vital role in cellular division by facilitating kinetochore assembly and spindle attachments. Despite their conserved functionality, centromeric DNA sequences exhibit rapid evolution, presenting diverse sizes and compositions across species. The functional significance of rye centromeric DNA sequences, particularly in centromere identity, remains unclear. In this study, we comprehensively characterized the sequence composition and organization of rye centromeres. Our findings revealed that these centromeres are primarily composed of long terminal repeat retrotransposons (LTR-RTs) and interspersed minisatellites. We systematically classified LTR-RTs into five categories, highlighting the prevalence of younger CRS1, CRS2, and CRS3 of CRSs (centromeric retrotransposons of Secale cereale) were primarily located in the core centromeres and exhibited a higher association with CENH3 nucleosomes. The minisatellites, mainly derived from retrotransposons, along with CRSs, played a pivotal role in establishing functional centromeres in rye. Additionally, we observed the formation of R-loops at specific regions of CRS1, CRS2, and CRS3, with both rye pericentromeres and centromeres exhibiting enrichment in R-loops. Notably, these R-loops selectively formed at binding regions of the CENH3 nucleosome in rye centromeres, suggesting a potential role in mediating the precise loading of CENH3 to centromeres and contributing to centromere specification. Our work provides insights into the DNA sequence composition, distribution, and potential function of R-loops in rye centromeres. This knowledge contributes valuable information to understanding the genetics and epigenetics of rye centromeres, offering implications for the development of synthetic centromeres in future plant modifications and beyond.

LncRNAs: the art of being influential without protein

The plant long noncoding (lnc)RNA field is on the brink of transitioning from large-scale identification of lncRNAs to their functional characterization. Due to the cross-kingdom conservation of interaction types and molecular functions, there is much to be learned from mammalian lncRNA research. Here, we discuss the different molecular processes involving lncRNAs from the regulation of chromatin to splicing. Furthermore, we discuss the lncRNA interactome, which includes proteins, other RNAs, and DNA. We explore and discuss how mammalian lncRNA functionalities could be reflected in similar pathways in plants and hypothesize that several breakthroughs in mammalian research could lead to the discovery of novel plant lncRNA molecular functions. Expanding our knowledge of the biological role of lncRNAs and their multiple applications paves the way for future agricultural applications.

Non-B-form DNA is associated with centromere stability in newly-formed polyploid wheat

Non-B-form DNA differs from the classic B-DNA double helix structure and plays a crucial regulatory role in replication and transcription. However, the role of non-B-form DNA in centromeres, especially in polyploid wheat, remains elusive. Here, we systematically analyzed seven non-B-form DNA motif profiles (A-phased DNA repeat, direct repeat, G-quadruplex, inverted repeat, mirror repeat, short tandem repeat, and Z-DNA) in hexaploid wheat. We found that three of these non-B-form DNA motifs were enriched at centromeric regions, especially at the CENH3-binding sites, suggesting that non-B-form DNA may create a favorable loading environment for the CENH3 nucleosome. To investigate the dynamics of centromeric non-B form DNA during the alloploidization process, we analyzed DNA secondary structure using CENH3 ChIP-seq data from newly formed allotetraploid wheat and its two diploid ancestors. We found that newly formed allotetraploid wheat formed more non-B-form DNA in centromeric regions compared with their parents, suggesting that non-B-form DNA is related to the localization of the centromeric regions in newly formed wheat. Furthermore, non-B-form DNA enriched in the centromeric regions was found to preferentially form on young LTR retrotransposons, explaining CENH3's tendency to bind to younger LTR. Collectively, our study describes the landscape of non-B-form DNA in the wheat genome, and sheds light on its potential role in the evolution of polyploid centromeres.

Approaches for Mapping and Analysis of R-loops

R-loops are nucleic acid structures composed of a DNA:RNA hybrid with a displaced non-template single-stranded DNA. Current approaches to identify and map R-loop formation across the genome employ either an antibody targeted against R-loops (S9.6) or a catalytically inactivated form of RNase H1 (dRNH1), a nuclease that can bind and resolve DNA:RNA hybrids via RNA exonuclease activity. This overview article outlines several ways to map R-loops using either methodology, explaining the differences and similarities among the approaches. Bioinformatic analysis of R-loops involves several layers of quality control and processing before visualizing the data. This article provides resources and tools that can be used to accurately process R-loop mapping data and explains the advantages and disadvantages of the resources as compared to one another. © 2024 Wiley Periodicals LLC.

The roles of DNA methylation on pH dependent i-motif (iM) formation in rice

I-motifs (iMs) are four-stranded non-B DNA structures containing C-rich DNA sequences. The formation of iMs is sensitive to pH conditions and DNA methylation, although the extent of which is still unknown in both humans and plants. To investigate this, we here conducted iMab antibody-based immunoprecipitation and sequencing (iM-IP-seq) along with bisulfite sequencing using CK (original genomic DNA without methylation-related treatments) and hypermethylated or demethylated DNA at both pH 5.5 and 7.0 in rice, establishing a link between pH, DNA methylation and iM formation on a genome-wide scale. We found that iMs folded at pH 7.0 displayed higher methylation levels than those formed at pH 5.5. DNA demethylation and hypermethylation differently influenced iM formation at pH 7.0 and 5.5. Importantly, CG hypo-DMRs (differentially methylated regions) and CHH (H = A, C and T) hyper-DMRs alone or coordinated with CG/CHG hyper-DMRs may play determinant roles in the regulation of pH dependent iM formation. Thus, our study shows that the nature of DNA sequences alone or combined with their methylation status plays critical roles in determining pH-dependent formation of iMs. It therefore deepens the understanding of the pH and methylation dependent modulation of iM formation, which has important biological implications and practical applications.

The Epigenomic Features and Potential Functions of PEG- and PDS-Favorable DNA G-Quadruplexes in Rice

A G-quadruplex (G4) is a typical non-B DNA structure and involved in various DNA-templated events in eukaryotic genomes. PEG and PDS chemicals have been widely applied for promoting the folding of in vivo or in vitro G4s. However, how PEG and PDS preferentially affect a subset of G4 formation genome-wide is still largely unknown. We here conducted a BG4-based IP-seq in vitro under K++PEG or K++PDS conditions in the rice genome. We found that PEG-favored IP-G4s+ have distinct sequence features, distinct genomic distributions and distinct associations with TEGs, non-TEGs and subtypes of TEs compared to PDS-favored ones. Strikingly, PEG-specific IP-G4s+ are associated with euchromatin with less enrichment levels of DNA methylation but with more enriched active histone marks, while PDS-specific IP-G4s+ are associated with heterochromatin with higher enrichment levels of DNA methylation and repressive marks. Moreover, we found that genes with PEG-specific IP-G4s+ are more expressed than those with PDS-specific IP-G4s+, suggesting that PEG/PDS-specific IP-G4s+ alone or coordinating with epigenetic marks are involved in the regulation of the differential expression of related genes, therefore functioning in distinct biological processes. Thus, our study provides new insights into differential impacts of PEG and PDS on G4 formation, thereby advancing our understanding of G4 biology.

R-loops act as regulatory switches modulating transcription of COLD-responsive genes in rice

COLD is a major naturally occurring stress that usually causes complex symptoms and severe yield loss in crops. R-loops function in various cellular processes, including development and stress responses, in plants. However, how R-loops function in COLD responses is largely unknown in COLD susceptible crops like rice (Oryza sativa L.). We conducted DRIP-Seq along with other omics data (RNA-Seq, DNase-Seq and ChIP-Seq) in rice with or without COLD treatment. COLD treatment caused R-loop reprogramming across the genome. COLD-biased R-loops had higher GC content and novel motifs for the binding of distinct transcription factors (TFs). Moreover, R-loops can directly/indirectly modulate the transcription of a subset of COLD-responsive genes, which can be mediated by R-loop overlapping TF-centered or cis-regulatory element-related regulatory networks and lncRNAs, accounting for c. 60% of COLD-induced expression of differential genes in rice, which is different from the findings in Arabidopsis. We validated two R-loop loci with contrasting (negative/positive) roles in the regulation of two individual COLD-responsive gene expression, as potential targets for enhanced COLD resistance. Our study provides detailed evidence showing functions of R-loop reprogramming during COLD responses and provides some potential R-loop loci for genetic and epigenetic manipulation toward breeding of rice varieties with enhanced COLD tolerance.

Ribonucleotide and R-Loop Damage in Plastid DNA and Mitochondrial DNA during Maize Development

Although the temporary presence of ribonucleotides in DNA is normal, their persistence represents a form of DNA damage. Here, we assess such damage and damage defense to DNA in plastids and mitochondria of maize. Shoot development proceeds from meristematic, non-pigmented cells containing proplastids and promitochondria at the leaf base to non-dividing green cells in the leaf blade containing mature organelles. The organellar DNAs (orgDNAs) become fragmented during this transition. Previously, orgDNA damage and damage defense of two types, oxidative and glycation, was described in maize, and now a third type, ribonucleotide damage, is reported. We hypothesized that ribonucleotide damage changes during leaf development and could contribute to the demise of orgDNAs. The levels of ribonucleotides and R-loops in orgDNAs and of RNase H proteins in organelles were measured throughout leaf development and in leaves grown in light and dark conditions. The data reveal that ribonucleotide damage to orgDNAs increased by about 2- to 5-fold during normal maize development from basal meristem to green leaf and when leaves were grown in normal light conditions compared to in the dark. During this developmental transition, the levels of the major agent of defense, RNase H, declined. The decline in organellar genome integrity during maize development may be attributed to oxidative, glycation, and ribonucleotide damages that are not repaired.

A DNA adenine demethylase impairs PRC2-mediated repression of genes marked by a specific chromatin signature

BACKGROUND

The Fe (II)- and α-ketoglutarate-dependent AlkB family dioxygenases are implicated in nucleotide demethylation. AlkB homolog1 (ALKBH1) is shown to demethylate DNA adenine methylation (6mA) preferentially from single-stranded or unpaired DNA, while its demethylase activity and function in the chromatin context are unclear.

RESULTS

Here, we find that loss-of-function of the rice ALKBH1 gene leads to increased 6mA in the R-loop regions of the genome but has a limited effect on the overall 6mA level. However, in the context of mixed tissues, rather than on individual loci, the ALKBH1 mutation or overexpression mainly affects the expression of genes with a specific combination of chromatin modifications in the body region marked with H3K4me3 and H3K27me3 but depleted of DNA CG methylation. In the similar context of mixed tissues, further analysis reveals that the ALKBH1 protein preferentially binds to genes marked by the chromatin signature and has a function to maintain a high H3K4me3/H3K27me3 ratio by impairing the binding of Polycomb repressive complex 2 (PRC2) to the targets, which is required for both the basal and stress-induced expression of the genes.

CONCLUSION

Our findings unravel a function of ALKBH1 to control the balance between the antagonistic histone methylations for gene activity and provide insight into the regulatory mechanism of PRC2-mediated H3K27me3 deposition within the gene body region.

Side-by-side comparison of G-quadruplex (G4) capture efficiency of the antibody BG4 versus the small-molecule ligands TASQs

The search for G-quadruplex (G4)-forming sequences across the genome is motivated by their involvement in key cellular processes and their putative roles in dysregulations underlying human genetic diseases. Sequencing-based methods have been developed to assess the prevalence of DNA G4s genome wide, including G4-seq to detect G4s in purified DNA (in vitro) using the G4 stabilizer PDS, and G4 chromatin immunoprecipitation sequencing (G4 ChIP-seq) to detect G4s in in situ fixed chromatin (in vivo) using the G4-specific antibody BG4. We recently reported on G4-RNA precipitation and sequencing (G4RP-seq) to assess the in vivo prevalence of RNA G4 landscapes transcriptome wide using the small molecule BioTASQ. Here, we apply this technique for mapping DNA G4s in plants (rice) and compare the efficiency of this new technique, G4-DNA precipitation and sequencing, G4DP-seq, to that of BG4-DNA-IP-seq that we developed for mapping of DNA G4s in rice using BG4. By doing so, we compare the G4 capture ability of small-sized ligands (BioTASQ and BioCyTASQ) versus the antibody BG4.

Non-B-form DNA tends to form in centromeric regions and has undergone changes in polyploid oat subgenomes

Centromeres are the specialized regions of the chromosomes that direct faithful chromosome segregation during cell division. Despite their functional conservation, centromeres display features of rapidly evolving DNA and wide evolutionary diversity in size and organization. Previous work found that the noncanonical B-form DNA structures are abundant in the centromeres of several eukaryotic species with a possible implication for centromere specification. Thus far, systematic studies into the organization and function of non-B-form DNA in plants remain scarce. Here, we applied the oat system to investigate the role of non-B-form DNA in centromeres. We conducted chromatin immunoprecipitation sequencing using an antibody to the centromere-specific histone H3 variant (CENH3); this accurately positioned oat centromeres with different ploidy levels and identified a series of centromere-specific sequences including minisatellites and retrotransposons. To define genetic characteristics of oat centromeres, we surveyed the repeat sequences and found that dyad symmetries were abundant in oat centromeres and were predicted to form non-B-DNA structures in vivo. These structures including bent DNA, slipped DNA, Z-DNA, G-quadruplexes, and R-loops were prone to form within CENH3-binding regions. Dynamic conformational changes of predicted non-B-DNA occurred during the evolution from diploid to tetraploid to hexaploid oat. Furthermore, we applied the single-molecule technique of AFM and DNA:RNA immunoprecipitation with deep sequencing to validate R-loop enrichment in oat centromeres. Centromeric retrotransposons exhibited strong associations with R-loop formation. Taken together, our study elucidates the fundamental character of non-B-form DNA in the oat genome and reveals its potential role in centromeres.

R-loop: The new genome regulatory element in plants

An R-loop is a three-stranded chromatin structure that consists of a displaced single strand of DNA and an RNA:DNA hybrid duplex, which was thought to be a rare by-product of transcription. However, recent genome-wide data have shown that R-loops are widespread and pervasive in a variety of genomes, and a growing body of experimental evidence indicates that R-loops have both beneficial and harmful effects on an organism. To maximize benefit and avoid harm, organisms have evolved several means by which they tightly regulate R-loop levels. Here, we summarize our current understanding of the biogenesis and effects of R-loops, the mechanisms that regulate them, and methods of R-loop profiling, reviewing recent research advances on R-loops in plants. Furthermore, we provide perspectives on future research directions for R-loop biology in plants, which might lead to a more comprehensive understanding of R-loop functions in plant genome regulation and contribute to future agricultural improvements.

Light, chromatin, action: nuclear events regulating light signaling in Arabidopsis

The plant nucleus provides a major hub for environmental signal integration at the chromatin level. Multiple light signaling pathways operate and exchange information by regulating a large repertoire of gene targets that shape plant responses to a changing environment. In addition to the established role of transcription factors in triggering photoregulated changes in gene expression, there are eminent reports on the significance of chromatin regulators and nuclear scaffold dynamics in promoting light-induced plant responses. Here, we report and discuss recent advances in chromatin-regulatory mechanisms modulating plant architecture and development in response to light, including the molecular and physiological roles of key modifications such as DNA, RNA and histone methylation, and/or acetylation. The significance of the formation of biomolecular condensates of key light signaling components is discussed and potential applications to agricultural practices overviewed.

Epigenomic Features and Potential Functions of K+ and Na+ Favorable DNA G-Quadruplexes in Rice

DNA G-quadruplexes (G4s) are non-canonical four-stranded DNA structures involved in various biological processes in eukaryotes. Molecularly crowded solutions and monovalent cations have been reported to stabilize in vitro and in vivo G4 formation. However, how K⁺ and Na⁺ affect G4 formation genome-wide is still unclear in plants. Here, we conducted BG4-DNA-IP-seq, DNA immunoprecipitation with anti-BG4 antibody coupled with sequencing, under K⁺ and Na⁺ + PEG conditions in vitro. We found that K⁺-specific IP-G4s had a longer peak size, more GC and PQS content, and distinct AT and GC skews compared to Na⁺-specific IP-G4s. Moreover, K⁺- and Na⁺-specific IP-G4s exhibited differential subgenomic enrichment and distinct putative functional motifs for the binding of certain trans-factors. More importantly, we found that K⁺-specific IP-G4s were more associated with active marks, such as active histone marks, and low DNA methylation levels, as compared to Na⁺-specific IP-G4s; thus, K⁺-specific IP-G4s in combination with active chromatin features facilitate the expression of overlapping genes. In addition, K⁺- and Na⁺-specific IP-G4 overlapping genes exhibited differential GO (gene ontology) terms, suggesting they may have distinct biological relevance in rice. Thus, our study, for the first time, explores the effects of K⁺ and Na⁺ on global G4 formation in vitro, thereby providing valuable resources for functional G4 studies in rice. It will provide certain G4 loci for the biotechnological engineering of rice in the future.

Splicing-mediated activation of SHAGGY-like kinases underpinning carbon partitioning in Arabidopsis seeds

Glycogen synthase kinase 3 (GSK3) family members serve as signaling hubs for plant development and stress responses, yet the underlying mechanism of their transcriptional regulation remains a long-standing mystery. Here we show that the transcription of SHAGGY-like kinase 11/12 (SK11/12), two members of the GSK3 gene family, is promoted by the splicing factor SmD1b, which is essential for distributing carbon sources into storage and protective components in Arabidopsis seeds. The chromatin recruitment of SmD1b at the SK11/12 loci promotes their transcription associated with co-transcriptional splicing of the first introns in the 5'-untranslated region of SK11/12. The loss of SmD1b function generates transcripts with unspliced introns that create disruptive R-loops to hamper the transcriptional elongation of SK11/12, in addition to compromising the recruitment of RNA polymerase II to the SK11/12 genomic regions. These effects imposed by SmD1b determine the transcription of SK11/12 to confer a key switch of carbon flow among metabolic pathways in zygotic and maternal tissues in seeds.

Non-B DNA structures emerging from plant genomes

Noncanonical (non-B) DNA structures, including G-quadruplexes, R-loops, and various DNA alternate conformations, affect genome organization, stability, and activity. Very few non-B structures have been mapped across plant genomes. Recently, Ma et al. identified intercalated-motifs (cytosine-rich, four-stranded quadruplex DNA) in rice, expanding how non-B DNA may contribute to plant gene regulation.

R-loop mapping and characterization during Drosophila embryogenesis reveals developmental plasticity in R-loop signatures

R-loops are involved in transcriptional regulation, DNA and histone post-translational modifications, genome replication and genome stability. To what extent R-loop abundance and genome-wide localization is actively regulated during metazoan embryogenesis is unknown. Drosophila embryogenesis provides a powerful system to address these questions due to its well-characterized developmental program, the sudden onset of zygotic transcription and available genome-wide data sets. Here, we measure the overall abundance and genome localization of R-loops in early and late-stage embryos relative to Drosophila cultured cells. We demonstrate that absolute R-loop levels change during embryogenesis and that RNaseH1 catalytic activity is critical for embryonic development. R-loop mapping by strand-specific DRIP-seq reveals that R-loop localization is plastic across development, both in the genes which form R-loops and where they localize relative to gene bodies. Importantly, these changes are not driven by changes in the transcriptional program. Negative GC skew and absolute changes in AT skew are associated with R-loop formation in Drosophila. Furthermore, we demonstrate that while some chromatin binding proteins and histone modifications such as H3K27me3 are associated with R-loops throughout development, other chromatin factors associated with R-loops in a developmental specific manner. Our findings highlight the importance and developmental plasticity of R-loops during Drosophila embryogenesis.

R-loops at microRNA encoding loci promote co-transcriptional processing of pri-miRNAs in plants

In most organisms, the maturation of nascent RNAs is coupled to transcription. Unlike in animals, the RNA polymerase II (RNAPII) transcribes microRNA genes (MIRNAs) as long and structurally variable pri-miRNAs in plants. Current evidence suggests that the miRNA biogenesis complex assembly initiates early during the transcription of pri-miRNAs in plants. However, it is unknown whether miRNA processing occurs co-transcriptionally. Here, we used native elongating transcript sequencing data and imaging techniques to demonstrate that plant miRNA biogenesis occurs coupled to transcription. We found that the entire biogenesis occurs co-transcriptionally for pri-miRNAs processed from the loop of the hairpin but requires a second nucleoplasmic step for those processed from the base. Furthermore, we found that co- and post-transcriptional miRNA processing mechanisms co-exist for most miRNAs in a dynamic balance. Notably, we discovered that R-loops, formed near the transcription start site region of MIRNAs, promote co-transcriptional pri-miRNA processing. Furthermore, our results suggest the neofunctionalization of co-transcriptionally processed miRNAs, boosting countless regulatory scenarios.

Role of Nuclear Non-Canonical Nucleic Acid Structures in Organismal Development and Adaptation to Stress Conditions

Over the last decades, numerous examples have involved nuclear non-coding RNAs (ncRNAs) in the regulation of gene expression. ncRNAs can interact with the genome by forming non-canonical nucleic acid structures such as R-loops or DNA:RNA triplexes. They bind chromatin and DNA modifiers and transcription factors and favor or prevent their targeting to specific DNA sequences and regulate gene expression of diverse genes. We review the function of these non-canonical nucleic acid structures in regulating gene expression of multicellular organisms during development and in response to different stress conditions and DNA damage using examples described in several organisms, from plants to humans. We also overview recent techniques developed to study where R-loops or DNA:RNA triplexes are formed in the genome and their interaction with proteins.

Epigenomic features of DNA G-quadruplexes and their roles in regulating rice gene transcription

A DNA G-quadruplex (G4) is a non-canonical four-stranded nucleic acid structure involved in many biological processes in mammals. The current knowledge on plant DNA G4s, however, is limited; whether and how DNA G4s impact gene expression in plants is still largely unknown. Here, we applied a protocol referred to as BG4-DNA-IP-seq followed by a comprehensive characterization of DNA G4s in rice (Oryza sativa L.); we next integrated dG4s (experimentally detectable G4s) with existing omics data and found that dG4s exhibited differential DNA methylation between transposable element (TE) and non-TE genes. dG4 regions displayed genic-dependent enrichment of epigenomic signatures; finally, we showed that these sites displayed a positive association with expression of DNA G4-containing genes when located at promoters, and a negative association when located in the gene body, suggesting localization-dependent promotional/repressive roles of DNA G4s in regulating gene transcription. This study reveals interrelations between DNA G4s and epigenomic signatures, as well as implicates DNA G4s in modulating gene transcription in rice. Our study provides valuable resources for the functional characterization or bioengineering of some of key DNA G4s in rice.

Genome-wide characterization of i-motifs and their potential roles in the stability and evolution of transposable elements in rice

I-motifs (iMs) are non-canonical DNA secondary structures that fold from cytosine (C)-rich genomic DNA regions termed putative i-motif forming sequences (PiMFSs). The structure of iMs is stabilized by hemiprotonated C-C base pairs, and their functions are now suspected in key cellular processes in human cells such as genome stability and regulation of gene transcription. In plants, their biological relevance is still largely unknown. Here, we characterized PiMFSs with high potential for i-motif formation in the rice genome by developing and applying a protocol hinging on an iMab antibody-based immunoprecipitation (IP) coupled with high-throughput sequencing (seq), consequently termed iM-IP-seq. We found that PiMFSs had intrinsic subgenomic distributions, cis-regulatory functions and an intricate relationship with DNA methylation. We indeed found that the coordination of PiMFSs with DNA methylation may affect dynamics of transposable elements (TEs) among different cultivated Oryza subpopulations or during evolution of wild rice species. Collectively, our study provides first and unique insights into the biology of iMs in plants, with potential applications in plant biotechnology for improving important agronomic rice traits.

The landscape of promoter-centred RNA-DNA interactions in rice

Chromatin-associated RNAs play key roles in various biological processes. However, both their repository and conjugation genomic loci and potential functions remain largely unclear. Here, we develop an effective method for mapping of chromatin-associated RNA-DNA interactions, followed by paired-end-tag sequencing (ChRD-PET) in rice. We present a comprehensive interaction map between RNAs and H3K4me3-marked regions based on H3K4me3 ChRD-PET data, showing three types of RNA-DNA interactions-local, proximal and distal. We further characterize the origin and composition of the RNA strand in R-loop RNA-DNA hybrids and identify that extensive cis and trans RNAs, including trans-non-coding RNAs, are prevalently involved in the R-loop. Integrative analysis of rice epigenome and three-dimensional genome data suggests that both coding and non-coding RNAs engage extensively in the formation of chromatin loops and chromatin-interacting domains. In summary, ChRD-PET is an efficient method for studying the features of RNA-chromatin interactions, and the resulting datasets constitute a valuable resource for the study of RNAs and their biological functions.

Prospects and challenges of epigenomics in crop improvement

BACKGROUND

The advent of high-throughput epigenome mapping techniques has ushered in a new era of multiomics with powerful tools now available to map and record genomic output at different levels. Integrating the different components of the epigenome from these multiomics measures allows investigations of cis-regulatory elements on a genome-scale. Mapping of chromatin state, chromatin accessibility dynamics, and higher-order chromatin structure enables a new level of understanding of cell fate determination, identity and function in normal growth and development, disease resistance, and yield.

OBJECTIVE

In this paper, the recent advances in epigenomics research of rice, maize, and wheat are reviewed, and the development trends of epigenomics of major crops in the coming years are projected.

METHODS

We highlight the role of epigenomics in regulating growth and development and identifying potential distal cis-regulatory elements in three major crops, and discuss the prospects and challenges for new epigenetics-mediated breeding technologies in crop improvement.

CONCLUSION

In this review, we summarize and analyze recent epigenomic advances in three major crops epigenomics and discuss possibilities and challenges for future research in the field.

R-loop and its functions at the regulatory interfaces between transcription and (epi)genome

R-loop represents a prevalent and specialized chromatin structure critically involved in a wide range of biological processes. In particular, co-transcriptional R-loops, produced often due to RNA polymerase pausing or RNA biogenesis malfunction, can initiate molecular events to context-dependently regulate local gene transcription and crosstalk with chromatin modifications. Cellular "readers" of R-loops are identified, exerting crucial impacts on R-loop homeostasis and gene regulation. Mounting evidence also supports R-loop deregulation as a frequent, sometimes initiating, event during the development of human pathologies, notably cancer and neurological disorder. The purpose of this review is to cover recent advances in understanding the fundamentals of R-loop biology, which have started to unveil complex interplays of R-loops with factors involved in various biological processes such as transcription, RNA processing and epitranscriptomic modification (such as N6-methyladenosine), DNA damage sensing and repair, and epigenetic regulation.

Toward an understanding of the detection and function of R-loops in plants

Although lagging behind studies in humans and other mammals, studies of R-loops in plants have recently entered an exciting stage in which the roles of R-loops in gene expression, genome stability, epigenomic signatures, and plant development and stress responses are being elucidated. Here, we review the strengths and weaknesses of existing methodologies, which were largely developed for R-loop studies in mammals, and then discuss the potential challenges of applying these methodologies to R-loop studies in plants. We then focus on recent advances in the functional characterization of R-loops in Arabidopsis thaliana and rice. Recent studies in plants indicate that there are coordinated relationships between R-loops and gene expression, and between R-loops and epigenomic signatures that depend, in part, on the types of R-loops involved. Finally, we discuss the emerging roles of R-loops in plants and directions for future research.

Genome-wide mapping reveals R-loops associated with centromeric repeats in maize

R-loops are stable chromatin structures comprising a DNA:RNA hybrid and a displaced single-stranded DNA. R-loops have been implicated in gene expression and chromatin structure, as well as in replication blocks and genome instability. Here, we conducted a genome-wide identification of R-loops and identified more than 700,000 R-loop peaks in the maize (Zea mays) genome. We found that sense R-loops were mainly enriched in promoters and transcription termination sites and relatively less enriched in gene bodies, which is different from the main gene-body localization of sense R-loops in Arabidopsis and Oryza sativa. At the chromosome scale, maize R-loops were enriched in pericentromeric heterochromatin regions, and a significant portion of R-loops were derived from transposable elements. In centromeres, R-loops preferentially formed within the binding regions of the centromere-specific histone CENH3, and centromeric retrotransposons were strongly associated with R-loop formation. Furthermore, centromeric retrotransposon R-loops were observed by applying the single-molecule imaging technique of atomic force microscopy. These findings elucidate the fundamental character of R-loops in the maize genome and reveal the potential role of R-loops in centromeres.

Interplay of DNA and RNA N 6 -methyladenine with R-loops in regulating gene transcription in Arabidopsis

R-loops and covalent modifications of N ⁶ -methyladenine on DNA (D-6 mA) or RNA (R-m⁶A) have been documented to function in various cellular processes in eukaryotes. However, the relationships between R-loops and both covalent modifications are still elusive in plants. Here, we integrated existing ssDRIP-seq with D-6 mA and R-m⁶A data from Arabidopsis thaliana. We found that the presence of either of both modifications facilitates R-loop formation and transcription of overlapping genes. Interestingly, our study suggests that the presence of R-m⁶A is key to affect R-loop intensity and positively regulate gene transcription. Moreover, the presence of D-6 mA plays an additive role to facilitate the effect of R-m⁶A on R-loop intensity, however, D-6 mA may negatively regulate gene transcription when coexisted with R-m⁶A. Our analyses indicate that D-6 mA, R-m⁶A, or histone marks may act individually and cooperatively with R-loops in regulating gene transcription. Our study is the first to link R-loops with D-6 mA and R-m⁶A in plants, thereby providing new insights into interactions between R-loops with D-6 mA, R-m⁶A, and histone marks for regulating gene transcription.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01010-5.

Identification of R-loop-forming Sequences in Drosophila melanogaster Embryos and Tissue Culture Cells Using DRIP-seq

R-loops are non-canonical nucleic structures composed of an RNA-DNA hybrid and a displaced ssDNA. Originally identified as a source of genomic instability, R-loops have been shown over the last decade to be involved in the targeting of proteins and to be associated with different histone modifications, suggesting a regulatory function. In addition, R-loops have been demonstrated to form differentially during the development of different tissues in plants and to be associated with diseases in mammals. Here, we provide a single-strand DRIP-seq protocol to identify R-loop-forming sequences in Drosophila melanogaster embryos and tissue culture cells. This protocol differs from earlier DRIP protocols in the fragmentation step. Sonication, unlike restriction enzymes, generates a homogeneous and highly reproducible nucleic acid fragment pool. In addition, it allows the use of this protocol in any organism with minimal optimization. This protocol integrates several steps from published protocols to identify R-loop-forming sequences with high stringency, suitable for de novo characterization. Graphic abstract: Figure 1.Overview of the strand-specific DRIP-seq protocol.

Emerging roles of centromeric RNAs in centromere formation and function

BACKGROUND

Centromeres are specialized chromosomal domains involved in kinetochore formation and faithful chromosome segregation. Despite a high level of functional conservation, centromeres are not identified by DNA sequences, but by epigenetic means. Universally, centromeres are typically formed on highly repetitive DNA, which were previously considered to be silent. However, recent studies have shown that transcription occurs in this region, known as centromeric-derived RNAs (cenRNAs). CenRNAs that contribute to fundamental aspects of centromere function have been recently investigated in detail. However, the distribution, behavior and contributions of centromeric transcripts are still poorly understood.

OBJECTIVE

The aim of this article is to provide an overview of the roles of cenRNAs in centromere formation and function.

METHODS

We describe the structure and DNA sequence of centromere from yeast to human. In addition, we briefly introduce the roles of cenRNAs in centromere formation and function, kinetochore structure, accurate chromosome segregation, and pericentromeric heterochromatin assembly. Centromeric circular RNAs (circRNAs) and R-loops are rising stars in centromere function. CircRNAs have been successfully identified in various species with the assistance of high-throughput sequencing and novel computational approaches for non-polyadenylated RNA transcripts. Centromeric R-loops can be identified by the single-strand DNA ligation-based library preparation technique. But the molecular features and function of these centromeric R-loops and circRNAs are still being investigated.

CONCLUSION

In this review, we summarize recent findings on the epigenetic regulation of cenRNAs across species, which would provide useful information about cenRNAs and interesting hints for further studies.

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

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

Decoding the plant genome: From epigenome to 3D organization

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

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

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