Single-molecule regulatory architectures captured by chromatin fiber sequencing

Histone modification inhibitors: An emerging frontier in thyroid Cancer therapy

Thyroid cancer (TC) is the most common endocrine cancer and is a serious health concern due to its aggressiveness and high incidence. Histone modifications affect DNA accessibility and gene transcriptional activity by altering the structure of chromatin. Abnormal histone modifications may affect genome stability and disrupt gene expression patterns, leading to many diseases, including cancer. A growing body of research suggests that histone modifications and TC progression are inextricably linked. This article discusses the impact of aberrant histone modification patterns on TC. By targeting specific histone-modifying enzymes, it may be possible to regulate gene expression and inhibit the growth of TC. Finally, we summarize the relevant histone modification inhibitors to better understand the development stage of the use of these drugs to inhibit histone-modifying enzymes in cancer treatment.

Quantification and potential functional relevance of binding cooperativity of adjacent transcription factors on DNA

In eukaryotes, the expression of specific genes is regulated by a combination of transcription factors (TFs) bound on regulatory regions of the genomic DNA (promoters and enhancers). Recent advances in genomic sequencing technology have enabled the measurements of TFs' footprints and binding affinities on DNA at the single-molecule level, facilitating the probing of binding cooperativity among adjacent TFs. This necessitates quantitative descriptions of TFs' binding cooperativity and understanding of its potential functional relevance. In this study, we show that the binding cooperativities between two adjacent TFs can be quantified by the [Formula: see text] coefficient, which can be experimentally determined. Under thermodynamic equilibrium, the binding affinities of two TFs can either increase together (positive cooperativity) or decrease together (negative cooperativity), but not in opposing directions (one increases while the other decreases). Within the framework of thermodynamics, we investigate the functional relevance of cooperativity. The functional relevance of positive cooperativity, which has been extensively discussed in the literature, is the sigmoidal binding curve around a TF concentration threshold (analogous to oxygen binding to hemoglobin), whereas the functional relevance of negative cooperativity is twofold. First, mutual exclusion of the two TFs enables bidirectional gene switching, akin to the CI-Cro system in phage [Formula: see text]. Second, while TFs often exhibit intranuclear concentration fluctuations, negative binding cooperativity assures fast TF dissociation from DNA and hence rapid response for gene expression regulation. Furthermore, the nonequilibrium steady states of living cells can lead to either positive or negative cooperativity, which can also be quantified by the [Formula: see text] coefficient.

Mouse promoters are characterised by low occupancy and high turnover of RNA polymerase II

The general transcription machinery and its occupancy at promoters are highly conserved across metazoans. This contrasts with the kinetics of mRNA production that considerably differ between model species such as Drosophila and mouse. The molecular basis for these kinetic differences is currently unknown. Here, we used Single-Molecule Footprinting to measure RNA Polymerase II (Pol II) occupancy, the fraction of DNA molecules bound, at promoters in mouse and Drosophila cell lines. Single-molecule data reveals that Pol II occupancy is on average 3-5 times more frequent at transcriptionally active Drosophila promoters than active mouse promoters. Kinetic modelling of the occupancy states suggests that these differences in Pol II occupancy are determined by the ratio between the transcription initiation and Pol II turnover rates. We used chemical perturbation of transcription initiation to determine Pol II turnover rate in both species. Integration of these data into the model shows that infrequent Pol II occupancy in mouse is explained by the combination of high Pol II turnover and low transcription initiation rates.

Always on the Move: Overview on Chromatin Dynamics within Nuclear Processes

Our genome is organized into chromatin, a dynamic and modular structure made of nucleosomes. Chromatin organization controls access to the DNA sequence, playing a fundamental role in cell identity and function. How nucleosomes enable these processes is an active area of study. In this review, we provide an overview of chromatin dynamics, its properties, mechanisms, and functions. We highlight the diverse ways by which chromatin dynamics is controlled during transcription, DNA replication, and repair. Recent technological developments have promoted discoveries in this area, to which we provide an outlook on future research directions.

Centromeres in the thermotolerant yeast K. marxianus mediate attachment to a single microtubule

Eukaryotic chromosome segregation requires spindle microtubules to attach to chromosomes through kinetochores. The chromosomal locus that mediates kinetochore assembly is the centromere and is epigenetically specified in most organisms by a centromeric histone H3 variant called CENP-A. An exception to this is budding yeast which have short, sequenced-defined point centromeres. In S. cerevisiae , a single CENP-A nucleosome is formed at the centromere and is sufficient for kinetochore assembly. The thermophilic budding yeast Kluyveromyces marxianus also has a point centromere but its length is nearly double the S. cerevisiae centromere and the number of centromeric nucleosomes and kinetochore attachment sites is unknown. Purification of native kinetochores from K. marxianus yielded a mixed population, with one subpopulation that appeared to consist of doublets, making it unclear whether K. marxianus shares the same attachment architecture as S. cerevisiae. Here, we demonstrate that though the doublet kinetochores have a functional impact on kinetochore strength, kinetochore localization throughout the cell cycle appears conserved between these two yeasts. In addition, whole spindle electron tomography demonstrates that a single microtubule binds to each chromosome. Single-molecule nucleosome mapping analysis suggests the presence of a single centromeric nucleosome. Taken together, we propose that the K. marxianus point centromere assembles a single centromeric nucleosome that mediates attachment to one microtubule.

Zygotic activation of transposable elements during zebrafish early embryogenesis

Although previous studies have shown that transposable elements (TEs) are conservatively activated to play key roles during early embryonic development, the details of zygotic TE activation (ZTA) remain poorly understood. Here, we employ long-read sequencing to precisely identify that only a small subset of TE loci are activated among numerous copies, allowing us to map their hierarchical transcriptional cascades at the single-locus and single-transcript level. Despite the heterogeneity of ZTA across family, subfamily, locus, and transcript levels, our findings reveal that ZTA follows a markedly different pattern from conventional zygotic gene activation (ZGA): ZTA occurs significantly later than ZGA and shows a pronounced bias for nuclear localization of TE transcripts. This study advances our understanding of TE activation by providing a high-resolution view of TE copies and creating a comprehensive catalog of thousands of previously unannotated transcripts and genes that are activated during early zebrafish embryogenesis. Among these genes, we highlight two that are essential for zebrafish development.

Unraveling the hidden complexity of cancer through long-read sequencing

Cancer is fundamentally a disease of the genome, characterized by extensive genomic, transcriptomic, and epigenomic alterations. Most current studies predominantly use short-read sequencing, gene panels, or microarrays to explore these alterations; however, these technologies can systematically miss or misrepresent certain types of alterations, especially structural variants, complex rearrangements, and alterations within repetitive regions. Long-read sequencing is rapidly emerging as a transformative technology for cancer research by providing a comprehensive view across the genome, transcriptome, and epigenome, including the ability to detect alterations that previous technologies have overlooked. In this Perspective, we explore the current applications of long-read sequencing for both germline and somatic cancer analysis. We provide an overview of the computational methodologies tailored to long-read data and highlight key discoveries and resources within cancer genomics that were previously inaccessible with prior technologies. We also address future opportunities and persistent challenges, including the experimental and computational requirements needed to scale to larger sample sizes, the hurdles in sequencing and analyzing complex cancer genomes, and opportunities for leveraging machine learning and artificial intelligence technologies for cancer informatics. We further discuss how the telomere-to-telomere genome and the emerging human pangenome could enhance the resolution of cancer genome analysis, potentially revolutionizing early detection and disease monitoring in patients. Finally, we outline strategies for transitioning long-read sequencing from research applications to routine clinical practice.

Evolution of genome-wide methylation profiling technologies

In this mini-review, we explore the advancements in genome-wide DNA methylation profiling, tracing the evolution from traditional methods such as methylation arrays and whole-genome bisulfite sequencing to the cutting-edge single-molecule profiling enabled by long-read sequencing (LRS) technologies. We highlight how LRS is transforming clinical and translational research, particularly by its ability to simultaneously measure genetic and epigenetic information, providing a more comprehensive understanding of complex disease mechanisms. We discuss current challenges and future directions in the field, emphasizing the need for innovative computational tools and robust, reproducible approaches to fully harness the capabilities of LRS in molecular diagnostics.

A Hitchhiker's Guide to long-read genomic analysis

Over the past decade, long-read sequencing has evolved into a pivotal technology for uncovering the hidden and complex regions of the genome. Significant cost efficiency, scalability, and accuracy advancements have driven this evolution. Concurrently, novel analytical methods have emerged to harness the full potential of long reads. These advancements have enabled milestones such as the first fully completed human genome, enhanced identification and understanding of complex genomic variants, and deeper insights into the interplay between epigenetics and genomic variation. This mini-review provides a comprehensive overview of the latest developments in long-read DNA sequencing analysis, encompassing reference-based and de novo assembly approaches. We explore the entire workflow, from initial data processing to variant calling and annotation, focusing on how these methods improve our ability to interpret a wide array of genomic variants. Additionally, we discuss the current challenges, limitations, and future directions in the field, offering a detailed examination of the state-of-the-art bioinformatics methods for long-read sequencing.

Transformer-based deep learning for accurate detection of multiple base modifications using single molecule real-time sequencing

We had previously reported a convolutional neural network (CNN) based approach, called the holistic kinetic model (HK model 1), for detecting 5-methylcytosine (5mC) by single molecule real-time sequencing (Pacific Biosciences). In this study, we constructed a hybrid model with CNN and transformer layers, named HK model 2. We improve the area under the receiver operating characteristic curve (AUC) for 5mC detection from 0.91 for HK model 1 to 0.99 for HK model 2. We further demonstrate that HK model 2 can detect other types of base modifications, such as 5-hydroxymethylcytosine (5hmC) and N6-methyladenine (6mA). Using HK model 2 to analyze 5mC patterns of cell-free DNA (cfDNA) molecules, we demonstrate the enhanced detection of patients with hepatocellular carcinoma, with an AUC of 0.97. Moreover, HK model 2-based detection of 6mA enables the detection of jagged ends of cfDNA and the delineation of cellular chromatin structures. HK model 2 is thus a versatile tool expanding the applications of single molecule real-time sequencing in liquid biopsies.

Conservation of dichromatin organization along regional centromeres

The attachment of the kinetochore to the centromere is essential for genome maintenance, yet the highly repetitive nature of satellite regional centromeres limits our understanding of their chromatin organization. We demonstrate that single-molecule chromatin fiber sequencing (Fiber-seq) can uniquely co-resolve kinetochore and surrounding chromatin architectures along point centromeres, revealing largely homogeneous single-molecule kinetochore occupancy. In contrast, the application of Fiber-seq to regional centromeres exposed marked per-molecule heterogeneity in their chromatin organization. Regional centromere cores uniquely contain a dichotomous chromatin organization (dichromatin) composed of compacted nucleosome arrays punctuated with highly accessible chromatin patches. CENP-B occupancy phases dichromatin to the underlying alpha-satellite repeat within centromere cores but is not necessary for dichromatin formation. Centromere core dichromatin is conserved between humans and primates, including along regional centromeres lacking satellite repeats. Overall, the chromatin organization of regional centromeres is defined by marked per-molecule heterogeneity, buffering kinetochore attachment against sequence and structural variability within regional centromeres.

Computational analysis of DNA methylation from long-read sequencing

DNA methylation is a critical epigenetic mechanism in numerous biological processes, including gene regulation, development, ageing and the onset of various diseases such as cancer. Studies of methylation are increasingly using single-molecule long-read sequencing technologies to simultaneously measure epigenetic states such as DNA methylation with genomic variation. These long-read data sets have spurred the continuous development of advanced computational methods to gain insights into the roles of methylation in regulating chromatin structure and gene regulation. In this Review, we discuss the computational methods for calling methylation signals, contrasting methylation between samples, analysing cell-type diversity and gaining additional genomic insights, and then further discuss the challenges and future perspectives of tool development for DNA methylation research.

Genetic diversity and regulatory features of human-specific NOTCH2NL duplications

NOTCH2NL (NOTCH2-N-terminus-like) genes arose from incomplete, recent chromosome 1 segmental duplications implicated in human brain cortical expansion. Genetic characterization of these loci and their regulation is complicated by the fact they are embedded in large, nearly identical duplications that predispose to recurrent microdeletion syndromes. Using nearly complete long-read assemblies generated from 67 human and 12 ape haploid genomes, we show independent recurrent duplication among apes with functional copies emerging in humans ~2.1 million years ago. We distinguish NOTCH2NL paralogs present in every human haplotype (NOTCH2NLA) from copy number variable ones. We also characterize large-scale structural variation, including gene conversion, for 28% of haplotypes leading to a previously undescribed paralog, NOTCH2tv. Finally, we apply Fiber-seq and long-read transcript sequencing to human cortical neurospheres to characterize the regulatory landscape and find that the most fixed paralogs, NOTCH2 and NOTCH2NLA, harbor the greatest number of paralog-specific elements potentially driving their regulation.

DiMeLo-cito: a one-tube protocol for mapping protein-DNA interactions reveals CTCF bookmarking in mitosis

Genome regulation relies on complex and dynamic interactions between DNA and proteins. Recently, powerful methods have emerged that leverage third-generation sequencing to map protein-DNA interactions genome-wide. For example, Directed Methylation with Long-read sequencing (DiMeLo-seq) enables mapping of protein-DNA interactions along long, single chromatin fibers, including in highly repetitive genomic regions. However, DiMeLo-seq involves lossy centrifugation-based wash steps that limit its applicability to many sample types. To address this, we developed DiMeLo-cito, a single-tube, wash-free protocol that maximizes the yield and quality of genomic DNA obtained for long-read sequencing. This protocol enables the interrogation of genome-wide protein binding with as few as 100,000 cells and without the requirement of a nuclear envelope, enabling confident measurement of protein-DNA interactions during mitosis. Using this protocol, we detected strong binding of CTCF to mitotic chromosomes in diploid human cells, in contrast with earlier studies in karyotypically unstable cancer cell lines, suggesting that CTCF "bookmarks" specific sites critical for maintaining genome architecture across cell divisions. By expanding the capabilities of DiMeLo-seq to a broader range of sample types, DiMeLo-cito can provide new insights into genome regulation and organization.

SMAC: identifying DNA N6-methyladenine (6mA) at the single-molecule level using SMRT CCS data

DNA modifications, such as N6-methyladenine (6mA), play important roles in various processes in eukaryotes. Single-molecule, real-time (SMRT) sequencing enables the direct detection of DNA modifications without requiring special sample preparation. However, most SMRT-based studies of 6mA rely on ensemble-level consensus by combining multiple reads covering the same genomic position, which misses the single-molecule heterogeneity. While recent methods have aimed at single-molecule level detection of 6mA, limitations in sequencing platforms, resolution, accuracy, and usability restrict their application in comprehensive epigenetic studies. Here, we present SMAC (single-molecule 6mA analysis of CCS reads), a novel framework for accurately detecting 6mA at the single-molecule level using SMRT circular consensus sequencing (CCS) data from the Sequel II system. It is an automated method that streamlines the entire workflow by packaging both existing softwares and built-in scripts, with user-defined parameters to allow easy adaptation for various studies. By utilizing the statistical distribution characteristics of enzyme kinetic indicators on single DNA molecules rather than a fixed cutoff, SMAC significantly improves 6mA detection accuracy at the single-nucleotide and single-molecule levels. It simplifies analysis by providing comprehensive information, including quality control, statistical analysis, and site visualization, directly from raw sequencing data. SMAC is a powerful new tool that enables de novo detection of 6mA and empowers investigation of its functions in modulating physiological processes.

Stable centromere association of the yeast histone variant Cse4 requires its essential N-terminal domain

Chromosome segregation relies on kinetochores that assemble on specialized centromeric chromatin containing a histone H3 variant. In budding yeast, a single centromeric nucleosome containing Cse4 assembles at a sequence-defined 125 bp centromere. Yeast centromeric sequences are poor templates for nucleosome formation in vitro, suggesting the existence of mechanisms that specifically stabilize Cse4 nucleosomes in vivo. The extended Cse4 N-terminal tail binds to the chaperone Scm3, and a short essential region called END within the N-terminal tail binds the inner kinetochore complex Okp1/Ame1. To address the roles of these interactions, we utilized single-molecule fluorescence assays to monitor Cse4 during kinetochore assembly. We found that Okp1/Ame1 and Scm3 independently stabilize Cse4 at centromeres via their END interaction. Scm3 and Cse4 stability at the centromere are enhanced by Ipl1/Aurora B phosphorylation of the Cse4 END, identifying a previously unknown role for Ipl1 in ensuring Cse4 stability. Strikingly, a phosphomimetic mutation in the Cse4 END restores Cse4 recruitment in mutants defective in Okp1/Ame1 binding. Together, these data suggest that a key function of the essential Cse4 N-terminus is to ensure Cse4 localization at centromeres.

Synchronized long-read genome, methylome, epigenome and transcriptome profiling resolve a Mendelian condition

Resolving the molecular basis of a Mendelian condition remains challenging owing to the diverse mechanisms by which genetic variants cause disease. To address this, we developed a synchronized long-read genome, methylome, epigenome and transcriptome sequencing approach, which enables accurate single-nucleotide, insertion-deletion and structural variant calling and diploid de novo genome assembly. This permits the simultaneous elucidation of haplotype-resolved CpG methylation, chromatin accessibility and full-length transcript information in a single long-read sequencing run. Application of this approach to an Undiagnosed Diseases Network participant with a chromosome X;13-balanced translocation of uncertain significance revealed that this translocation disrupted the functioning of four separate genes (NBEA, PDK3, MAB21L1 and RB1) previously associated with single-gene Mendelian conditions. Notably, the function of each gene was disrupted via a distinct mechanism that required integration of the four 'omes' to resolve. These included fusion transcript formation, enhancer adoption, transcriptional readthrough silencing and inappropriate X-chromosome inactivation of autosomal genes. Overall, this highlights the utility of synchronized long-read multi-omic profiling for mechanistically resolving complex phenotypes.

Multiscale footprints reveal the organization of cis-regulatory elements

Cis-regulatory elements (CREs) control gene expression and are dynamic in their structure and function, reflecting changes in the composition of diverse effector proteins over time1. However, methods for measuring the organization of effector proteins at CREs across the genome are limited, hampering efforts to connect CRE structure to their function in cell fate and disease. Here we developed PRINT, a computational method that identifies footprints of DNA-protein interactions from bulk and single-cell chromatin accessibility data across multiple scales of protein size. Using these multiscale footprints, we created the seq2PRINT framework, which uses deep learning to allow precise inference of transcription factor and nucleosome binding and interprets regulatory logic at CREs. Applying seq2PRINT to single-cell chromatin accessibility data from human bone marrow, we observe sequential establishment and widening of CREs centred on pioneer factors across haematopoiesis. We further discover age-associated alterations in the structure of CREs in murine haematopoietic stem cells, including widespread reduction of nucleosome footprints and gain of de novo identified Ets composite motifs. Collectively, we establish a method for obtaining rich insights into DNA-binding protein dynamics from chromatin accessibility data, and reveal the architecture of regulatory elements across differentiation and ageing.

Interpreting the CTCF-mediated sequence grammar of genome folding with AkitaV2

Interphase mammalian genomes are folded in 3D with complex locus-specific patterns that impact gene regulation. CTCF (CCCTC-binding factor) is a key architectural protein that binds specific DNA sites, halts cohesin-mediated loop extrusion, and enables long-range chromatin interactions. There are hundreds of thousands of annotated CTCF-binding sites in mammalian genomes; disruptions of some result in distinct phenotypes, while others have no visible effect. Despite their importance, the determinants of which CTCF sites are necessary for genome folding and gene regulation remain unclear. Here, we update and utilize Akita, a convolutional neural network model, to extract the sequence preferences and grammar of CTCF contributing to genome folding. Our analyses of individual CTCF sites reveal four predictions: (i) only a small fraction of genomic sites are impactful; (ii) impact is highly dependent on sequences flanking the core CTCF binding motif; (iii) core and flanking nucleotides contribute largely additively to the overall impact of a site; (iv) sites created as combinations of different core and flanking sequences have impacts proportional to the product of their average impacts, i.e. they are broadly compatible. Our analysis of collections of CTCF sites make two predictions for multi-motif grammar: (i) insulation strength depends on the number of CTCF sites within a cluster, and (ii) pattern formation is governed by the orientation and spacing of these sites, rather than any inherent specialization of the CTCF motifs themselves. In sum, we present a framework for using neural network models to probe the sequences instructing genome folding and provide a number of predictions to guide future experimental inquiries.

The regulatory potential of transposable elements in maize

The genomes of flowering plants consist largely of transposable elements (TEs), some of which modulate gene regulation and function. However, the repetitive nature of TEs and difficulty of mapping individual TEs by short-read-sequencing have hindered our understanding of their regulatory potential. We demonstrate that long-read chromatin fiber sequencing (Fiber-seq) comprehensively identifies accessible chromatin regions (ACRs) and CpG methylation across the maize genome. We uncover stereotypical ACR patterns at young TEs that degenerate with evolutionary age, resulting in TE-enhancers preferentially marked by a novel plant-specific epigenetic feature: simultaneous hyper-CpG methylation and chromatin accessibility. We show that TE ACRs are co-opted as gene promoters and that ACR-containing TEs can facilitate gene amplification. Lastly, we uncover a pervasive epigenetic signature - hypo-5mCpG methylation and diffuse chromatin accessibility - directing TEs to specific loci, including the loci that sparked McClintock's discovery of TEs.

Centromeres in the thermotolerant yeast K. marxianus mediate attachment to a single microtubule

Eukaryotic chromosome segregation requires spindle microtubules to attach to chromosomes through kinetochores. The chromosomal locus that mediates kinetochore assembly is the centromere and is epigenetically specified in most organisms by a centromeric histone H3 variant called CENP-A. An exception to this is budding yeast which have short, sequenced-defined point centromeres. In S. cerevisiae, a single CENP-A nucleosome is formed at the centromere and is sufficient for kinetochore assembly. The thermophilic budding yeast Kluyveromyces marxianus also has a point centromere but its length is nearly double the S. cerevisiae centromere and the number of centromeric nucleosomes and kinetochore attachment sites is unknown. Purification of native kinetochores from K. marxianus yielded a mixed population, with one subpopulation that appeared to consist of doublets, making it unclear whether K. marxianus shares the same attachment architecture as S. cerevisiae. Here, we demonstrate that though the doublet kinetochores have a functional impact on kinetochore strength, kinetochore localization throughout the cell cycle appears conserved between these two yeasts. In addition, whole spindle electron tomography demonstrates that a single microtubule binds to each chromosome. Single-molecule nucleosome mapping analysis suggests the presence of a single centromeric nucleosome. Taken together, we propose that the K. marxianus point centromere assembles a single centromeric nucleosome that mediates attachment to one microtubule.

Pervasive and programmed nucleosome distortion patterns on single mammalian chromatin fibers

We present a genome-scale method to map the single-molecule co-occupancy of structurally distinct nucleosomes, subnucleosomes, and other protein-DNA interactions via long-read high-resolution adenine methyltransferase footprinting. Iteratively Defined Lengths of Inaccessibility (IDLI) classifies nucleosomes on the basis of shared patterns of intranucleosomal accessibility, into: i.) minimally-accessible chromatosomes; ii.) octasomes with stereotyped DNA accessibility from superhelical locations (SHLs) ±1 through ±7; iii.) highly-accessible unwrapped nucleosomes; and iv.) subnucleosomal species, such as hexasomes, tetrasomes, and other short DNA protections. Applying IDLI to mouse embryonic stem cell (mESC) chromatin, we discover widespread nucleosomal distortion on individual mammalian chromatin fibers, with >85% of nucleosomes surveyed displaying degrees of intranucleosomally accessible DNA. We observe epigenomic-domain-specific patterns of distorted nucleosome co-occupancy and positioning, including at enhancers, promoters, and mouse satellite repeat sequences. Nucleosome distortion is programmed by the presence of bound transcription factors (TFs) at cognate motifs; occupied TF binding sites are differentially decorated by distorted nucleosomes compared to unbound sites, and degradation experiments establish direct roles for TFs in structuring binding-site proximal nucleosomes. Finally, we apply IDLI in the context of primary mouse hepatocytes, observing evidence for pervasive nucleosomal distortion in vivo. Further genetic experiments reveal a role for the hepatocyte master regulator FOXA2 in directly impacting nucleosome distortion at hepatocyte-specific regulatory elements in vivo. Our work suggests extreme-but regulated-plasticity in nucleosomal DNA accessibility at the single-molecule level. Further, our study offers an essential new framework to model transcription factor binding, nucleosome remodeling, and cell-type specific gene regulation across biological contexts.

The single-molecule accessibility landscape of newly replicated mammalian chromatin

We present replication-aware single-molecule accessibility mapping (RASAM), a method to nondestructively measure replication status and protein-DNA interactions on chromatin genome-wide. Using RASAM, we uncover a genome-wide state of single-molecule "hyperaccessibility" post-replication that resolves over several hours. Combining RASAM with cellular models for rapid protein degradation, we demonstrate that histone chaperone CAF-1 reduces nascent chromatin accessibility by filling single-molecular "gaps" and generating closely spaced dinucleosomes on replicated DNA. At cis-regulatory elements, we observe unique modes by which nascent chromatin hyperaccessibility resolves: at CCCTC-binding factor (CTCF)-binding sites, CTCF and nucleosomes compete, reducing CTCF occupancy and motif accessibility post-replication; at active transcription start sites, high chromatin accessibility is maintained, implying rapid re-establishment of nucleosome-free regions. Our study introduces a new paradigm for studying replicated chromatin fiber organization. More broadly, we uncover a unique organization of newly replicated chromatin that must be reset by active processes, providing a substrate for epigenetic reprogramming.

Chromatin architecture mapping by multiplex proximity tagging

Chromatin plays a pivotal role in genome expression, maintenance, and replication. To better understand chromatin organization, we developed a novel proximity-tagging method which assigns unique DNA barcodes to molecules that associate in 3D space. Using this method - Proximity Copy Paste (PCP) - we mapped the connectivity of individual nucleosomes in Saccharomyces cerevisiae. By analyzing nucleosome positions and spacing on single molecule fibers, we show that chromatin is predominantly organized into regularly spaced nucleosome arrays that can be positioned or delocalized. Basic features of nucleosome arrays are generally explained by gene size and transcription. PCP can also map long-range, multi-way interactions and we provide the first direct evidence supporting a model that metaphase chromosomes are compacted by cohesin loop clustering. Analyzing single-molecule nuclease footprinting data we define distinct chromatin states within a mixed population to show that non-canonical nucleosomes, notably Overlapping-Di-Nucleosomes (OLDN) are a stable feature of chromatin. PCP is a versatile method allowing the detection of the connectivity of individual molecules locally and over large distance to be mapped at high-resolution in a single experiment.

Profiling the epigenome using long-read sequencing

The advent of single-molecule, long-read sequencing (LRS) technologies by Oxford Nanopore Technologies and Pacific Biosciences has revolutionized genomics, transcriptomics and, more recently, epigenomics research. These technologies offer distinct advantages, including the direct detection of methylated DNA and simultaneous assessment of DNA sequences spanning multiple kilobases along with their modifications at the single-molecule level. This has enabled the development of new assays for analyzing chromatin states and made it possible to integrate data for DNA methylation, chromatin accessibility, transcription factor binding and histone modifications, thereby facilitating comprehensive epigenomic profiling. Owing to recent advancements, alternative, nascent and translating transcripts can be detected using LRS approaches. This Review discusses LRS-based experimental and computational strategies for characterizing chromatin states and highlights their advantages over short-read sequencing methods. Furthermore, we demonstrate how various long-read methods can be integrated to design multi-omics studies to investigate the relationship between chromatin states and transcriptional dynamics.

Genome-wide single-cell and single-molecule footprinting of transcription factors with deaminase

Decades of research have established that mammalian transcription factors (TFs) bind to each gene's regulatory regions and cooperatively control tissue specificity, timing, and intensity of gene transcription. Mapping the combination of TF binding sites genome wide is critically important for understanding functional genomics. Here, we report a technique to measure TFs' binding sites on the human genome with a near single-base resolution by footprinting with deaminase (FOODIE) on a single-molecule and single-cell basis. Single-molecule sequencing reads after enzymatic deamination allow detection of the TF binding fraction on a particular footprint and the binding cooperativity of any two adjacent TFs, which can be either positive or negative. As a newcomer of single-cell genomics, single-cell FOODIE enables the detection of cell-type-specific TF footprints in a pure cell population in a heterogeneous tissue, such as the brain. We found that genes carrying out a certain biological function together in a housing-keeping correlated gene module (CGM) or a tissues-specific CGM are coordinated by shared TFs in the gene's promoters and enhancers, respectively. Scalable and cost-effective, FOODIE allows us to create an open FOODIE database for cell lines, with applicability to human tissues and clinical samples.

Resolving the chromatin impact of mosaic variants with targeted Fiber-seq

Accurately quantifying the functional consequences of noncoding mosaic variants requires the pairing of DNA sequences with both accessible and closed chromatin architectures along individual DNA molecules-a pairing that cannot be achieved using traditional fragmentation-based chromatin assays. We demonstrate that targeted single-molecule chromatin fiber sequencing (Fiber-seq) achieves this, permitting single-molecule, long-read genomic, and epigenomic profiling across targeted >100 kb loci with ∼10-fold enrichment over untargeted sequencing. Targeted Fiber-seq reveals that pathogenic expansions of the DMPK CTG repeat that underlie Myotonic Dystrophy 1 are characterized by somatic instability and disruption of multiple nearby regulatory elements, both of which are repeat length-dependent. Furthermore, we reveal that therapeutic adenine base editing of the segmentally duplicated γ-globin (HBG1/HBG2) promoters in primary human hematopoietic cells induced toward an erythroblast lineage increases the accessibility of the HBG1 promoter as well as neighboring regulatory elements. Overall, we find that these non-protein coding mosaic variants can have complex impacts on chromatin architectures, including extending beyond the regulatory element harboring the variant.

Single chromatin fiber profiling and nucleosome position mapping in the human brain

We apply a single-molecule chromatin fiber sequencing (Fiber-seq) protocol designed for amplification-free cell-type-specific mapping of the regulatory architecture at nucleosome resolution along extended ∼10-kb chromatin fibers to neuronal and non-neuronal nuclei sorted from human brain tissue. Specifically, application of this method enables the resolution of cell-selective promoter and enhancer architectures on single fibers, including transcription factor footprinting and position mapping, with sequence-specific fixation of nucleosome arrays flanking transcription start sites and regulatory motifs. We uncover haplotype-specific chromatin patterns, multiple regulatory elements cis-aligned on individual fibers, and accessible chromatin at 20,000 unique sites encompassing retrotransposons and other repeat sequences hitherto "unmappable" by short-read epigenomic sequencing. Overall, we show that Fiber-seq is applicable to human brain tissue, offering sharp demarcation of nucleosome-depleted regions at sites of open chromatin in conjunction with multi-kilobase nucleosomal positioning at single-fiber resolution on a genome-wide scale.

Opportunities and challenges of single-cell and spatially resolved genomics methods for neuroscience discovery

Over the past decade, single-cell genomics technologies have allowed scalable profiling of cell-type-specific features, which has substantially increased our ability to study cellular diversity and transcriptional programs in heterogeneous tissues. Yet our understanding of mechanisms of gene regulation or the rules that govern interactions between cell types is still limited. The advent of new computational pipelines and technologies, such as single-cell epigenomics and spatially resolved transcriptomics, has created opportunities to explore two new axes of biological variation: cell-intrinsic regulation of cell states and expression programs and interactions between cells. Here, we summarize the most promising and robust technologies in these areas, discuss their strengths and limitations and discuss key computational approaches for analysis of these complex datasets. We highlight how data sharing and integration, documentation, visualization and benchmarking of results contribute to transparency, reproducibility, collaboration and democratization in neuroscience, and discuss needs and opportunities for future technology development and analysis.

Mapping protein-DNA interactions with DiMeLo-seq

We recently developed directed methylation with long-read sequencing (DiMeLo-seq) to map protein-DNA interactions genome wide. DiMeLo-seq is capable of mapping multiple interaction sites on single DNA molecules, profiling protein binding in the context of endogenous DNA methylation, identifying haplotype-specific protein-DNA interactions and mapping protein-DNA interactions in repetitive regions of the genome that are difficult to study with short-read methods. With DiMeLo-seq, adenines in the vicinity of a protein of interest are methylated in situ by tethering the Hia5 methyltransferase to an antibody using protein A. Protein-DNA interactions are then detected by direct readout of adenine methylation with long-read, single-molecule DNA sequencing platforms such as Nanopore sequencing. Here we present a detailed protocol and practical guidance for performing DiMeLo-seq. This protocol can be run on nuclei from fresh, lightly fixed or frozen cells. The protocol requires 1-2 d for performing in situ targeted methylation, 1-5 d for library preparation depending on desired fragment length and 1-3 d for Nanopore sequencing depending on desired sequencing depth. The protocol requires basic molecular biology skills and equipment, as well as access to a Nanopore sequencer. We also provide a Python package, dimelo, for analysis of DiMeLo-seq data.

Better together: how cooperativity influences transcriptional bursting

Transcriptional bursting refers to the stochastic transition of a promoter between transcriptionally active and inactive states. This dynamic process is highly regulated by the dynamics of transcription factor binding to DNA, their interactions with coactivators, and the 3D interactions between promoters, condensates, and enhancers. In this mini-review, we discuss recent insights into the kinetics of transcription factors and cofactors in both simple and complex regulatory environments to understand their impact on transcriptional bursting. We examine the novel concept of transcription factor exchange and relate it to different cooperativity models. Finally, we discuss recent live-cell imaging studies on the regulation of transcriptional bursting by enhancers and transcriptional condensates.

Single-molecule states link transcription factor binding to gene expression

The binding of multiple transcription factors (TFs) to genomic enhancers drives gene expression in mammalian cells1. However, the molecular details that link enhancer sequence to TF binding, promoter state and transcription levels remain unclear. Here we applied single-molecule footprinting2,3 to measure the simultaneous occupancy of TFs, nucleosomes and other regulatory proteins on engineered enhancer-promoter constructs with variable numbers of TF binding sites for both a synthetic TF and an endogenous TF involved in the type I interferon response. Although TF binding events on nucleosome-free DNA are independent, activation domains recruit cofactors that destabilize nucleosomes, driving observed TF binding cooperativity. Average TF occupancy linearly determines promoter activity, and we decompose TF strength into separable binding and activation terms. Finally, we develop thermodynamic and kinetic models that quantitatively predict both the enhancer binding microstates and gene expression dynamics. This work provides a template for the quantitative dissection of distinct contributors to gene expression, including TF activation domains, concentration, binding affinity, binding site configuration and recruitment of chromatin regulators.

Independent expansion, selection, and hypervariability of the TBC1D3 gene family in humans

TBC1D3 is a primate-specific gene family that has expanded in the human lineage and has been implicated in neuronal progenitor proliferation and expansion of the frontal cortex. The gene family and its expression have been challenging to investigate because it is embedded in high-identity and highly variable segmental duplications. We sequenced and assembled the gene family using long-read sequencing data from 34 humans and 11 nonhuman primate species. Our analysis shows that this particular gene family has independently duplicated in at least five primate lineages, and the duplicated loci are enriched at sites of large-scale chromosomal rearrangements on Chromosome 17. We find that all human copy-number variation maps to two distinct clusters located at Chromosome 17q12 and that humans are highly structurally variable at this locus, differing by as many as 20 copies and ∼1 Mbp in length depending on haplotypes. We also show evidence of positive selection, as well as a significant change in the predicted human TBC1D3 protein sequence. Last, we find that, despite multiple duplications, human TBC1D3 expression is limited to a subset of copies and, most notably, from a single paralog group: TBC1D3-CDKL These observations may help explain why a gene potentially important in cortical development can be so variable in the human population.

DNA-m6A calling and integrated long-read epigenetic and genetic analysis with fibertools

Long-read DNA sequencing has recently emerged as a powerful tool for studying both genetic and epigenetic architectures at single-molecule and single-nucleotide resolution. Long-read epigenetic studies encompass both the direct identification of native cytosine methylation and the identification of exogenously placed DNA N 6 -methyladenine (DNA-m6A). However, detecting DNA-m6A modifications using single-molecule sequencing, as well as coprocessing single-molecule genetic and epigenetic architectures, is limited by computational demands and a lack of supporting tools. Here, we introduce fibertools, a state-of-the-art toolkit that features a semisupervised convolutional neural network for fast and accurate identification of m6A-marked bases using Pacific Biosciences (PacBio) single-molecule long-read sequencing, as well as the coprocessing of long-read genetic and epigenetic data produced using either the PacBio or Oxford Nanopore Technologies (ONT) sequencing platforms. We demonstrate accurate DNA-m6A identification (>90% precision and recall) along >20 kb long DNA molecules with an ∼1000-fold improvement in speed. In addition, we demonstrate that fibertools can readily integrate genetic and epigenetic data at single-molecule resolution, including the seamless conversion between molecular and reference coordinate systems, allowing for accurate genetic and epigenetic analyses of long-read data within structurally and somatically variable genomic regions.

Epigenetic control and inheritance of rDNA arrays

Ribosomal RNA (rRNA) genes exist in multiple copies arranged in tandem arrays known as ribosomal DNA (rDNA). The total number of gene copies is variable, and the mechanisms buffering this copy number variation remain unresolved. We surveyed the number, distribution, and activity of rDNA arrays at the level of individual chromosomes across multiple human and primate genomes. Each individual possessed a unique fingerprint of copy number distribution and activity of rDNA arrays. In some cases, entire rDNA arrays were transcriptionally silent. Silent rDNA arrays showed reduced association with the nucleolus and decreased interchromosomal interactions, indicating that the nucleolar organizer function of rDNA depends on transcriptional activity. Methyl-sequencing of flow-sorted chromosomes, combined with long read sequencing, showed epigenetic modification of rDNA promoter and coding region by DNA methylation. Silent arrays were in a closed chromatin state, as indicated by the accessibility profiles derived from Fiber-seq. Removing DNA methylation restored the transcriptional activity of silent arrays. Array activity status remained stable through the iPS cell re-programming. Family trio analysis demonstrated that the inactive rDNA haplotype can be traced to one of the parental genomes, suggesting that the epigenetic state of rDNA arrays may be heritable. We propose that the dosage of rRNA genes is epigenetically regulated by DNA methylation, and these methylation patterns specify nucleolar organizer function and can propagate transgenerationally.

RNA polymerases reshape chromatin architecture and couple transcription on individual fibers

RNA polymerases must initiate and pause within a complex chromatin environment, surrounded by nucleosomes and other transcriptional machinery. This environment creates a spatial arrangement along individual chromatin fibers ripe for both competition and coordination, yet these relationships remain largely unknown owing to the inherent limitations of traditional structural and sequencing methodologies. To address this, we employed long-read chromatin fiber sequencing (Fiber-seq) in Drosophila to visualize RNA polymerase (Pol) within its native chromatin context with single-molecule precision along up to 30 kb fibers. We demonstrate that Fiber-seq enables the identification of individual Pol II, nucleosome, and transcription factor footprints, revealing Pol II pausing-driven destabilization of downstream nucleosomes. Furthermore, we demonstrate pervasive direct distance-dependent transcriptional coupling between nearby Pol II genes, Pol III genes, and transcribed enhancers, modulated by local chromatin architecture. Overall, transcription initiation reshapes surrounding nucleosome architecture and couples nearby transcriptional machinery along individual chromatin fibers.

Determining mesoscale chromatin structure parameters from spatially correlated cleavage data using a coarse-grained oligonucleosome model

The three-dimensional structure of chromatin has emerged as an important feature of eukaryotic gene regulation. Recent technological advances in DNA sequencing-based assays have revealed locus- and chromatin state-specific structural patterns at the length scale of a few nucleosomes (~1 kb). However, interpreting these data sets remains challenging. Radiation-induced correlated cleavage of chromatin (RICC-seq) is one such chromatin structure assay that maps DNA-DNA-contacts at base pair resolution by sequencing single-stranded DNA fragments released from irradiated cells. Here, we develop a flexible modeling and simulation framework to enable the interpretation of RICC-seq data in terms of oligonucleosome structure ensembles. Nucleosomes are modeled as rigid bodies with excluded volume and adjustable DNA wrapping, connected by linker DNA modeled as a worm-like chain. We validate the model's parameters against cryo-electron microscopy and sedimentation data. Our results show that RICC-seq is sensitive to nucleosome spacing, nucleosomal DNA wrapping, and the strength of inter-nucleosome interactions. We show that nucleosome repeat lengths consistent with orthogonal assays can be extracted from experimental RICC-seq data using a 1D convolutional neural net trained on RICC-seq signal predicted from simulated ensembles. We thus provide a suite of analysis tools that add quantitative structural interpretability to RICC-seq experiments.

Exploring the reciprocity between pioneer factors and development

Development is regulated by coordinated changes in gene expression. Control of these changes in expression is largely governed by the binding of transcription factors to specific regulatory elements. However, the packaging of DNA into chromatin prevents the binding of many transcription factors. Pioneer factors overcome this barrier owing to unique properties that enable them to bind closed chromatin, promote accessibility and, in so doing, mediate binding of additional factors that activate gene expression. Because of these properties, pioneer factors act at the top of gene-regulatory networks and drive developmental transitions. Despite the ability to bind target motifs in closed chromatin, pioneer factors have cell type-specific chromatin occupancy and activity. Thus, developmental context clearly shapes pioneer-factor function. Here, we discuss this reciprocal interplay between pioneer factors and development: how pioneer factors control changes in cell fate and how cellular environment influences pioneer-factor binding and activity.

A common flanking variant is associated with enhanced stability of the FGF14-SCA27B repeat locus

The factors driving or preventing pathological expansion of tandem repeats remain largely unknown. Here, we assessed the FGF14 (GAA)·(TTC) repeat locus in 2,530 individuals by long-read and Sanger sequencing and identified a common 5'-flanking variant in 70.34% of alleles analyzed (3,463/4,923) that represents the phylogenetically ancestral allele and is present on all major haplotypes. This common sequence variation is present nearly exclusively on nonpathogenic alleles with fewer than 30 GAA-pure triplets and is associated with enhanced stability of the repeat locus upon intergenerational transmission and increased Fiber-seq chromatin accessibility.

Semiconservative transmission of DNA N 6-adenine methylation in a unicellular eukaryote

Although DNA N 6-adenine methylation (6mA) is best known in prokaryotes, its presence in eukaryotes has recently generated great interest. Biochemical and genetic evidence supports that AMT1, an MT-A70 family methyltransferase (MTase), is crucial for 6mA deposition in unicellular eukaryotes. Nonetheless, the 6mA transmission mechanism remains to be elucidated. Taking advantage of single-molecule real-time circular consensus sequencing (SMRT CCS), here we provide definitive evidence for semiconservative transmission of 6mA in Tetrahymena thermophila In wild-type (WT) cells, 6mA occurs at the self-complementary ApT dinucleotide, mostly in full methylation (full-6mApT); after DNA replication, hemi-methylation (hemi-6mApT) is transiently present on the parental strand, opposite to the daughter strand readily labeled by 5-bromo-2'-deoxyuridine (BrdU). In ΔAMT1 cells, 6mA predominantly occurs as hemi-6mApT. Hemi-to-full conversion in WT cells is fast, robust, and processive, whereas de novo methylation in ΔAMT1 cells is slow and sporadic. In Tetrahymena, regularly spaced 6mA clusters coincide with the linker DNA of nucleosomes arrayed in the gene body. Importantly, in vitro methylation of human chromatin by the reconstituted AMT1 complex recapitulates preferential targeting of hemi-6mApT sites in linker DNA, supporting AMT1's intrinsic and autonomous role in maintenance methylation. We conclude that 6mA is transmitted by a semiconservative mechanism: full-6mApT is split by DNA replication into hemi-6mApT, which is restored to full-6mApT by AMT1-dependent maintenance methylation. Our study dissects AMT1-dependent maintenance methylation and AMT1-independent de novo methylation, reveals a 6mA transmission pathway with a striking similarity to 5-methylcytosine (5mC) transmission at the CpG dinucleotide, and establishes 6mA as a bona fide eukaryotic epigenetic mark.

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.

Direct transposition of native DNA for sensitive multimodal single-molecule sequencing

Concurrent readout of sequence and base modifications from long unamplified DNA templates by Pacific Biosciences of California (PacBio) single-molecule sequencing requires large amounts of input material. Here we adapt Tn5 transposition to introduce hairpin oligonucleotides and fragment (tagment) limiting quantities of DNA for generating PacBio-compatible circular molecules. We developed two methods that implement tagmentation and use 90-99% less input than current protocols: (1) single-molecule real-time sequencing by tagmentation (SMRT-Tag), which allows detection of genetic variation and CpG methylation; and (2) single-molecule adenine-methylated oligonucleosome sequencing assay by tagmentation (SAMOSA-Tag), which uses exogenous adenine methylation to add a third channel for probing chromatin accessibility. SMRT-Tag of 40 ng or more human DNA (approximately 7,000 cell equivalents) yielded data comparable to gold standard whole-genome and bisulfite sequencing. SAMOSA-Tag of 30,000-50,000 nuclei resolved single-fiber chromatin structure, CTCF binding and DNA methylation in patient-derived prostate cancer xenografts and uncovered metastasis-associated global epigenome disorganization. Tagmentation thus promises to enable sensitive, scalable and multimodal single-molecule genomics for diverse basic and clinical applications.

Single-Macromolecule Studies of Eukaryotic Genomic Maintenance

Genomes are self-organized and self-maintained as long, complex macromolecules of chromatin. The inherent heterogeneity, stochasticity, phase separation, and chromatin dynamics of genome operation make it challenging to study genomes using ensemble methods. Various single-molecule force-, fluorescent-, and sequencing-based techniques rooted in different disciplines have been developed to fill critical gaps in the capabilities of bulk measurements, each providing unique, otherwise inaccessible, insights into the structure and maintenance of the genome. Capable of capturing molecular-level details about the organization, conformational changes, and packaging of genetic material, as well as processive and stochastic movements of maintenance factors, a single-molecule toolbox provides an excellent opportunity for collaborative research to understand how genetic material functions in health and malfunctions in disease. In this review, we discuss novel insights brought to genomic sciences by single-molecule techniques and their potential to continue to revolutionize the field-one molecule at a time.

Examining chromatin heterogeneity through PacBio long-read sequencing of M.EcoGII methylated genomes: an m6A detection efficiency and calling bias correcting pipeline

Recent studies have combined DNA methyltransferase footprinting of genomic DNA in nuclei with long-read sequencing, resulting in detailed chromatin maps for multi-kilobase stretches of genomic DNA from one cell. Theoretically, nucleosome footprints and nucleosome-depleted regions can be identified using M.EcoGII, which methylates adenines in any sequence context, providing a high-resolution map of accessible regions in each DNA molecule. Here, we report PacBio long-read sequence data for budding yeast nuclei treated with M.EcoGII and a bioinformatic pipeline which corrects for three key challenges undermining this promising method. First, detection of m6A in individual DNA molecules by the PacBio software is inefficient, resulting in false footprints predicted by random gaps of seemingly unmethylated adenines. Second, there is a strong bias against m6A base calling as AT content increases. Third, occasional methylation occurs within nucleosomes, breaking up their footprints. After correcting for these issues, our pipeline calculates a correlation coefficient-based score indicating the extent of chromatin heterogeneity within the cell population for every gene. Although the population average is consistent with that derived using other techniques, we observe a wide range of heterogeneity in nucleosome positions at the single-molecule level, probably reflecting cellular chromatin dynamics.

A data-adaptive methods in detecting exogenous methyltransferase accessible chromatin in human genome using nanopore sequencing

MOTIVATION

Identifying chromatin accessibility is one of the key steps in studying the regulation of eukaryotic genomes. The combination of exogenous methyltransferase and nanopore sequencing provides an strategy to identify open chromatin over long genomic ranges at the single-molecule scale. However, endogenous methylation, non-open-chromatin-specific exogenous methylation and base-calling errors limit the accuracy and hinders its application to complex genomes.

RESULTS

We systematically evaluated the impact of these three influence factors, and developed a model-based computational method, methyltransferase accessible genome region finder (MAGNIFIER), to address the issues. By incorporating control data, MAGNIFIER attenuates the three influence factors with data-adaptive comparison strategy. We demonstrate that MAGNIFIER is not only sensitive to identify the open chromatin with much improved accuracy, but also able to detect the chromatin accessibility of repetitive regions that are missed by NGS-based methods. By incorporating long-read RNA-seq data, we revealed the association between the accessible Alu elements and non-classic gene isoforms.

AVAILABILITY AND IMPLEMENTATION

Freely available on web at https://github.com/Goatofmountain/MAGNIFIER.

STR mutations on chromosome 15q cause thyrotropin resistance by activating a primate-specific enhancer of MIR7-2/MIR1179

Thyrotropin (TSH) is the master regulator of thyroid gland growth and function. Resistance to TSH (RTSH) describes conditions with reduced sensitivity to TSH. Dominantly inherited RTSH has been linked to a locus on chromosome 15q, but its genetic basis has remained elusive. Here we show that non-coding mutations in a (TTTG)4 short tandem repeat (STR) underlie dominantly inherited RTSH in all 82 affected participants from 12 unrelated families. The STR is contained in a primate-specific Alu retrotransposon with thyroid-specific cis-regulatory chromatin features. Fiber-seq and RNA-seq studies revealed that the mutant STR activates a thyroid-specific enhancer cluster, leading to haplotype-specific upregulation of the bicistronic MIR7-2/MIR1179 locus 35 kb downstream and overexpression of its microRNA products in the participants' thyrocytes. An imbalance in signaling pathways targeted by these micro-RNAs provides a working model for this cause of RTSH. This finding broadens our current knowledge of genetic defects altering pituitary-thyroid feedback regulation.

Genome organization and botanical diversity

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

Insight on ecDNA-mediated tumorigenesis and drug resistance

Extrachromosomal DNAs (ecDNAs) are a pervasive feature found in cancer and contain oncogenes and their corresponding regulatory elements. Their unique structural properties allow a rapid amplification of oncogenes and alter chromatin accessibility, leading to tumorigenesis and malignant development. The uneven segregation of ecDNA during cell division enhances intercellular genetic heterogeneity, which contributes to tumor evolution that might trigger drug resistance and chemotherapy tolerance. In addition, ecDNA has the ability to integrate into or detach from chromosomal DNA, such progress results into structural alterations and genomic rearrangements within cancer cells. Recent advances in multi-omics analysis revealing the genomic and epigenetic characteristics of ecDNA are anticipated to make valuable contributions to the development of precision cancer therapy. Herein, we conclud the mechanisms of ecDNA generation and the homeostasis of its dynamic structure. In addition to the latest techniques in ecDNA research including multi-omics analysis and biochemical validation methods, we also discuss the role of ecDNA in tumor development and treatment, especially in drug resistance, and future challenges of ecDNA in cancer therapy.