The study discloses the application of transformer-based deep learning models for the task of super-enhancers prediction in human tumor cell lines with a specific focus on sequence-only features within studied entities of super-enhancer and enhancer elements in the human genome. The proposed SE-prediction method included the GENA-LM application at handling long DNA sequences with the classification task, distinguishing super-enhancers from enhancers using H3K36me, H3K4me1, H3K4me3 and H3K27ac landscape datasets from HeLa, HEK293, H2171, Jurkat, K562, MM1S and U87 cell lines. The model was fine-tuned on relevant sequence data, allowing for the analysis of extended genomic sequences without the need for epigenetic markers as proposed in early approaches. The study achieved balanced accuracy metrics, surpassing previous models like SENet, particularly in HEK293 and K562 cell lines. Also, it was shown that super-enhancers frequently co-localize with epigenetic marks such as H3K4me3 and H3K27ac. Therefore, the attention mechanism of the model provided insights into the sequence features contributing to SE classification, indicating a correlation between sequence-only features and mentioned epigenetic landscapes. These findings support the potential transformer models use in further genomic sequence analysis for bioinformatics applications in enhancer/super-enhancer characterization and gene regulation studies.

Super-enhancers (SEs) orchestrate high-level transcription by integrating multiple regulatory elements and signals. Although chromatin accessibility and transcription factor binding within SEs are extensively studied, the role of non-canonical DNA structures, particularly Z-DNA, remains underexplored. In this study, genome-wide predictions of Z-DNA-forming sequences (generated by the Z-DNA-BERT model) were applied to systematically investigate their distribution within typical enhancers and SEs across multiple human cancer cell lines. Statistically significant enrichment of Z-DNA sequences within SE regions, compared to random genomic controls, was observed. Furthermore, genetic variants overlapping these Z-DNA regions, identified using data from the 1000 Genomes Project, were found to alter binding motifs of the SP/KLF transcription factor family. These mutations exhibited population-specific clustering and overlapped previously reported pathogenic copy-number variations (CNVs) associated with neurodevelopmental disorders, potentially affecting transcription factor binding motifs related to neuronal growth and differentiation pathways. Population-level phylogenetic analysis revealed distinct clustering patterns of these variants, suggesting frequency-specific genetic architecture. Overall, the computational findings indicate that Z-DNA structures within super-enhancers might play regulatory roles and potentially influence population-specific genetic variation, highlighting specific genomic targets and providing new avenues for future experimental research.

Aim. To study the relationship of gut microbiota (GM) composition with clinical and paraclinical characteristics of patients with coronary artery disease (CAD) and heart failure with reduced ejection fraction (HFrEF). Material and methods . The study included 64 patients (men, 71,9%) with CAD, 33 (men, 87,9%) with HFrEF and 38 (men, 63,2%) with various cardiovascular risk factors (control group). GM composition was assessed using real-time polymerase chain reaction and next- generation sequencing. To identify correlations of GM composition with clinical and paraclinical characteristics, the Spearman rank correlation coefficient (ρ) was calculated. Generalized linear models (QuasiPoisson, overdispersion) with covariates (age, body mass index) were used to compare groups by phylotype. Results. The relative abundance of Bacillota phylotypes is significantly higher in patients with CAD and HFrEF compared to controls (p<0,001). The Bacillota/Bacteroidota dysbiosis ratio was also significantly higher in patients with CAD and HFrEF than in the control group (p<0,001). According to correlation analysis, the relative representation of methanogenic archaea Methanomassiliicoccaceae, Methanomethylophilaceae (ρ=0,49, p<0,0003) directly correlated with ventricular arrhythmias in patients with HFrEF. Conclusion. Patients with CAD and HFrEF differ significantly from the control group in GM composition at the phylotype level. In the study groups, Bacillota/Bacteroidota dysbiosis ratio was also significantly higher than in the control group. The relative abundance of the archaea Methanomassiliicoccaceae and Methanomethylophilaceae is associated with episodes of ventricular tachycardia/ventricular fibrillation.

Super-enhancers (SEs) are large clusters of highly active enhancers that play key regulatory roles in cell identity, development, and disease. While conventional methods classify SEs in a binary fashion—super-enhancer or not—this threshold-based approach can overlook significant intermediate states of enhancer activity. Here, we present a dataset and accompanying framework that facilitate a more nuanced, non-binary examination of SE activation across mouse tissue types (mammary gland, lung tissue, and NMuMG cells) and various experimental conditions (normal, tumor, and drug-treated samples). By consolidating overlapping SE intervals and capturing continuous enhancer activity metrics (e.g., ChIP-seq signal intensities), our dataset reveals gradual transitions between moderate and high enhancer activity levels that are not captured by strictly binary classification. Additionally, the data include extensive functional annotations, linking SE loci to nearby genes and enabling immediate downstream analyses such as clustering and gene ontology enrichment. The flexible approach supports broader investigations of enhancer landscapes, offering a comprehensive platform for understanding how SE activation underpins disease mechanisms, therapeutic response, and developmental processes.

Long non-coding RNAs (lncRNAs) interact with chromatin and recruit epigenetic complexes to specific genomic loci, yet their relationship with super-enhancers (SEs), key regulatory elements frequently reprogrammed in cancer, remains unexplored. We developed an integrative pipeline that combines RNA–chromatin contact data (RNA-Chrom), histone modification–lncRNA expression correlation profiles (HiMoRNA peaks), and super-enhancer annotations (SEdb 3.0) to map lncRNA–SE regulatory axes. Applying this framework to SNHG1 in HCT116 colorectal cancer cells, we identified 21 SNHG1-reactive super-enhancers (Ψ-SEs) among 184 cancer-specific SEs, at which SNHG1 physical contacts co-occur with SNHG1-correlated histone modifications (HiMoRNA peaks), predominantly H3K4me1 (permutation p = 0.001, fold enrichment = 2.03). Comparison with 4145 lncRNAs demonstrated that epigenetic correlations alone do not distinguish SNHG1; instead, the addition of the contact layer is required to delineate the Ψ-SE set. Differential expression (DESeq2) and co-expression analyses in 471 TCGA-COAD tumor samples identified 12 Ψ-SE target genes (including CDC20, PDP1, and TOP1) consistently upregulated in both HCT116 cells and patient tumors and positively correlated with SNHG1, with the co-expression signal robust to tumor purity correction. The proposed Ψ/Ω classification provides a generalizable framework for prioritizing super-enhancers at which lncRNA–chromatin interactions may shape the local epigenetic environment across cancer types.