Study / S7LAJ71V02023-12-01

MAE-seq refines regulatory elements across the genome

Xiusheng Zhu, Qitong Huang, Lei Huang, Jing Luo, Qing Li et al.

About this study

Proper cell fate determination relies on precise spatial and temporal genome-wide cooperation between regulatory elements (REs) and their targeted genes. However, the lengths of REs defined using different methods vary, which indicates that there is sequence redundancy and that the context of the genome may be unintelligible. We developed a method called MAE-seq (Massive Active Enhancers by Sequencing) to experimentally identify functional REs at a 25-bp scale. In this study, MAE-seq was used to identify 626879, 541617 and 554826 25-bp enhancers in mouse embryonic stem cells (mESCs), C2C12 and HEK 293T, respectively. Using ∼1.6 trillion 25 bp DNA fragments and screening 12 billion cells, we identified 626879 as active enhancers in mESCs as an example. Comparative analysis revealed that most of the histone modification datasets were annotated by MAE-Seq loci. Furthermore, 33.85% (212195) of the identified enhancers were identified as de novo ones with no epigenetic modification. Intriguingly, distinct chromatin states dictate the requirement for dissimilar cofactors in governing novel and known enhancers. Validation results show that these 25-bp sequences could act as a functional unit, which shows identical or similar expression patterns as the previously defined larger elements, Enhanced resolution facilitated the identification of numerous cell-specific enhancers and their accurate annotation as super enhancers. Moreover, we characterized novel elements capable of augmenting gene activity. By integrating with high-resolution Hi-C data, over 55.64% of novel elements may have a distal association with different targeted genes. For example, we found that the Cdh1 gene interacts with one novel and two known REs in mESCs. The biological effects of these interactions were investigated using CRISPR-Cas9, revealing their role in coordinating Cdh1 gene expression and mESC proliferation. Our study presents an experimental approach to refine the REs at 25-bp resolution, advancing the precision of genome annotation and unveiling the underlying genome context. This novel approach not only advances our understanding of gene regulation but also opens avenues for comprehensive exploration of the genomic landscape.

Full author list & citation

Xiusheng Zhu, Qitong Huang, Lei Huang, Jing Luo, Qing Li, Dashuai Kong, Biao Deng, Yi Gu, Xueyan Wang, Chenying Li, Siyuan Kong, Yubo Zhang. MAE-seq refines regulatory elements across the genome. 2023-12-01. https://doi.org/10.1093/nar/gkad1129

Experiments 3

E3QMXUWED

mESC MAE-seq 25-bp enhancer screen

Transient pMX-mP-mCherry plasmid libraries carrying random 25-bp genomic fragments were transfected into E14 mouse embryonic stem cells; fluorescent cells were flow-sorted and sequenced against an input library. Four deposited output libraries (outputA-outputD) are represented as one mESC experiment.

Sort-Seq / Flow-Seq MPRAMousemm9
Explore data
E4TQ53FZ5

C2C12 MAE-seq 25-bp enhancer catalog

Supplementary Table S7 provides the final MAE-seq enhancer calls for C2C12 mouse myoblasts as 25-bp genomic coordinates classified as known or novel. The accessible GEO series does not contain a C2C12 per-locus count table, so this experiment packages the genuine MPRA-derived coordinate/classification endpoint available in the supplement.

Sort-Seq / Flow-Seq MPRAMousemm9
Explore data
E68UK392T

HEK293T MAE-seq 25-bp enhancer screen

Transient pMX-mP-mCherry plasmid libraries carrying random 25-bp genomic fragments were transfected into HEK293T cells; fluorescent cells were flow-sorted and sequenced against an input library. Three deposited output libraries (output1-output3) are represented as one HEK293T experiment.

Sort-Seq / Flow-Seq MPRAHumanhg19
Explore data

Raw source data 16 files

Original supplemental and deposited inputs retained for this study. Download files individually or together as a ZIP; nested folders are preserved. Source reuse terms apply, and sequencing reads may be omitted.

Download all 16 files (ZIP)gkad1129_supplemental_files.zipGSE193494_family.soft.gzGSE193494_filelist.txtGSM5808469_input_hg19_normalized.txt.gzGSM5808469_input_mm9_normalized.txt.gzGSM5808470_outputA_effective_cv1nc_q0.01.txt.gzGSM5808471_outputB_effective_cv1nc_q0.01.txt.gzGSM5808472_outputC_effective_cv1nc_q0.01.txt.gzGSM5808473_outputD_effective_cv1nc_q0.01.txt.gzGSM5808474_output1_effective_cv1nc_q0.01.txt.gzGSM5808475_output2_effective_cv1nc_q0.01.txt.gzGSM5808476_output3_effective_cv1nc_q0.01.txt.gzsource_manifest.txtTableS5_mESC_mm9.xlsxTableS6_HEK293T_hg19.csvTableS7_C2C12_mm9.csv

Cite OpenMPRA

Cite the OpenMPRA database. Include your access date because the collection changes over time.

Please also cite the source studies when using their data.