Bibliometric review of ATAC-Seq and its application in gene expression

ProOvErlap: Assessing feature proximity/overlap and testing statistical significance from genomic intervals

Feature overlap is a critical concept in bioinformatics and occurs when two genomic intervals, usually represented as BED files, are located in the same genomic regions. Instead, feature proximity refers to the spatial proximity of genomic elements. For example, promoters typically overlap or are close to the genes they regulate. Overlap and proximity are also important in epigenetic studies. Here, the overlap of regions enriched for specific epigenetic modifications or accessible chromatin can elucidate complex molecular phenotypes. Consequently, the ability to analyze and interpret feature overlap and proximity is essential for understanding the biological processes that contribute to a given phenotype. To address this need, we present a computational method capable of analyzing data represented in the BED format. This method aims to quantitatively assess the degree of proximity or overlap between genomic features and to determine the statistical significance of these events in the context of a non-parametric randomization test. The aim is to understand whether the observed state differs from what would be expected by chance. The method is designed to be easy to use, requiring only a single command line to run, allowing straightforward overlap and proximity analysis. It also provides clear visualizations and publication-quality figures. In conclusion, this study highlights the importance of feature overlap and proximity in epigenetic studies and presents a method to improve the systematic assessment and interpretation of these features. A new resource for identifying biologically significant interactions between genomic features in both healthy and disease states.

Genome-wide chromatin accessibility and transcriptome analysis reveal the up-regulation of immunosuppressive genes in macrophages under simulated microgravity

Long-term space missions are of growing research interest because of the space exploration. However, plenty of works focused on the impaired immune response, less attention has been paid to the activation of immunosuppressive or anti-inflammatory function. The molecular mechanism of immune disorder induced by microgravity still needs investigation. Here, we used a random positioning machine to generate a simulated microgravity environment and evaluated its effects on mouse RAW 264.7 macrophage cell line. We used ATAC-seq and RNA-seq for revealing the mechanism at chromatin level and gene level. From ATAC-seq, we obtained an average of 75,700,675 paired-end clean reads for each library and the mapping rates averaged at 96.8 %. The number of differential accessible regions were 510 for increased peaks, 638 for decreased peaks. From RNA-seq, we obtained 278 differentially expressed genes, of which 104 were down-regulated and 174 were up-regulated genes. Through ATAC-seq and RNA-seq multi-omics analysis, we identified a group of 17 genes. Then we chose 6 up-regulated genes (CD83, CEBPD, CXCR5, DUSP6, SEMA4B, TNFRSF22) that related to immunosuppressive function for further confirmation. The qRT-PCR results were consistent with sequencing results, which indicated that simulated microgravity leads to the up-regulated expression of immunosuppressive genes of macrophages. Taken together, our results offered novel insights for understanding the brief principles and mechanisms of simulated microgravity induced immune dysfunction to macrophage.

Epigenetics factors in schizophrenia: future directions for etiologic and therapeutic study approaches

Schizophrenia is a complex, heterogeneous, and highly disabling severe mental disorder whose pathogenesis has not yet been fully elucidated. Epigenetics, as a bridge between genetic and environmental factors, plays an important role in the pathophysiology of schizophrenia. Over the past decade, epigenetic-wide association studies have rapidly become an important branch of psychiatric research, especially in deciphering the molecular mechanisms of schizophrenia. This review systematically analyzes recent advances in epigenome-wide association studies (EWAS) of schizophrenia, focusing on technological developments. We synthesize findings from large-scale EWAS alongside emerging evidence on DNA methylation patterns, histone modifications, and regulatory networks, emphasizing their roles in disease mechanisms and treatment responses. In addition, this review provides a prospective outlook, evaluating the impact that technological developments may have on future studies of schizophrenia. With the continuous advancement of high-throughput sequencing technology and the increasing maturity of big data analysis methods, epigenetics is expected to have a significant impact on the early diagnosis, prognosis assessment and even personalized treatment of schizophrenia.

An unusual transfusion-dependent hemoglobin H disease caused by a novel complex inverted duplication involving the α-globin regulatory elements and α-thalassemia--SEA deletion

Four multi-species conserved sequences (MCSs) are important enhancers which affect α-globin expression. Deletions of MCS can cause α-thalassemia. So far, duplication of MCS has never been reported to account for thalassemia. In this study, an unusual transfusion-dependent case of hemoglobin H disease was identified by whole-genome sequencing, optical genome mapping and longer PCR with special primers, which was caused by a familial 96,620-bp inverted duplication (from MCS-R1 to MCS-R4), inserted between chr16:199348 and 199349 (GRCh37/hg19) within MCSs. The duplication segment included an inverted repeat sequence from chr16:102712 to176193 and one direct repeat sequence from chr16:176208 to 199348. The associated α-thalassemia trait was confirmed to result from disrupted topological chromatin domains using ATAC-seq and the dual‑luciferase reporter assay system. This case presents a new mechanism of α-thalassemia, and may aid our understanding of the effects of enhancers on gene expression and the differential contribution of the four enhancer elements in the human a-globin locus.

Integration of ATAC-seq and RNA-seq reveals the dynamics of chromatin accessibility and gene expression in zoysiagrass response to drought

KEY MESSAGE

The 'X4' accession of zoysiagrass demonstrated superior drought tolerance compared to other accessions. Integration analysis of transcriptomics and epigenomics revealed a positive correlation between ATAC-seq peak intensity and gene expression levels. Six motifs involved in regulating drought responses were identified, which are similar to the domains of the ERF and C2H2 transcription factor families. Heterologous expression of Zja11G000860 in yeast enhanced tolerance to drought stress, allowing robust growth even at high PEG6000 concentrations. Zoysiagrass is renowned for its drought tolerance and serves as an exceptional domestic turfgrass in China. However, the changes in chromatin accessibility during drought in zoysiagrass are not well understood. We conducted a preliminary evaluation of the phenotypic changes in drought tolerance for six zoysiagrass cultivars, taking into account their growth characteristics and physiological traits under drought conditions. Additionally, we utilized an integrated multi-omics strategy, encompassing RNA sequencing (RNA-seq), Assay for Transposase Accessible Chromatin using high-throughput sequencing (ATAC-seq), and reverse transcription quantitative PCR (RT-qPCR) verification experiments, to gain deeper understanding of the chromatin accessibility patterns linked to gene expression in response to drought stress in zoysiagrass. Preliminary analysis of the trends in relative water content and proline content suggested that the variety 'X4' exhibited superior drought tolerance compared to the other five accessions. The KEGG pathway enrichment analysis revealed that zoysiagrass responded to environmental stress by regulating stress response and antioxidant defense pathways. Notably, the expression levels of genes Zja03G031540 and Zja11G000860 were significantly increased in the 'X4' zoysiagrass genotype, which exhibited improved drought tolerance, compared to the 'X1' zoysiagrass genotype with reduced drought tolerance. This study suggested that 63 high-confidence genes are related to drought stress, including Zja03G031540 and Zja11G000860. Additionally, six motifs regulating drought responses were unearthed. Furthermore, the heterologous expression of Zja11G000860 in yeast enhanced tolerance to drought stress. The study discovered a positive correlation between ATAC-seq peak intensity and gene expression levels. The expression of high-confidence genes was linked to zoysiagrass resistance evaluation and phenotypic traits, implying that these genes are involved in responding to external drought stress. This study combined ATAC-seq and RNA-seq technologies for the first time to identify drought-related gene expression in zoysiagrass, elucidating the grass adaptation to environmental stress and the regulatory mechanisms underlying stress responses, and laying the groundwork for zoysiagrass improvement and breeding.

Integrative single-cell RNA-seq and ATAC-seq analysis of the evolutionary trajectory features of adipose-derived stem cells induced into astrocytes

This study employs single-cell RNA sequencing (scRNA-seq) and assay for transposase-accessible chromatin with high-throughput sequencing technologies (scATAC-seq) to perform joint sequencing on cells at various time points during the induction of adipose-derived stem cells (ADSCs) into astrocytes. We applied bioinformatics approaches to investigate the differentiation trajectories of ADSCs during their induced differentiation into astrocytes. Pseudotemporal analysis was used to infer differentiation trajectories. Additionally, we assessed chromatin accessibility patterns during the differentiation process. Key transcription factors driving the differentiation of ADSCs into astrocytes were identified using motif and footprint methods. Our analysis revealed significant shifts in gene expression during the induction process, with astrocyte-related genes upregulated and stem cell-related genes downregulated. ADSCs first differentiated into neural stem cell-like cells with high plasticity, which further matured into astrocytes via two distinct pathways. Marked changes in chromatin accessibility were observed during ADSC-induced differentiation, affecting transcription regulation and cell function. Transcription factors analysis identified NFIA/B/C/X and CEBPA/B/D as key regulators in ADSCs differentiation into astrocytes. We observed a correlation between chromatin accessibility and gene expression, with ADSCs exhibiting broad chromatin accessibility prior to lineage commitment, where chromatin opening precedes transcription initiation. In summary, we found that ADSCs first enter a neural stem cell-like state before differentiating into astrocytes. ADSCs also display extensive chromatin accessibility prior to astrocyte differentiation, although transcription has not yet been initiated. These findings offer a theoretical framework for understanding the molecular mechanisms underlying this process.

Advances on the Role of Lung Macrophages in the Pathogenesis of Chronic Obstructive Pulmonary Disease in the Era of Single-Cell Genomics

Chronic Obstructive Pulmonary Disease (COPD) is a heterogeneous respiratory disorder characterized by persistent airflow limitation. The diverse pathogenic mechanisms underlying COPD progression remain incompletely understood. Macrophages, serving as the most representative immune cells in the respiratory tract, constitute the first line of innate immune defense and maintain pulmonary immunological homeostasis. Recent advances have provided deeper insights into the phenotypic and functional alterations of pulmonary macrophages and their role in COPD pathogenesis. Notably, the advent of single-cell RNA sequencing has revolutionized our understanding of macrophage molecular heterogeneity in COPD. Herein, we review principal investigations concerning the sophisticated mechanisms through which pulmonary macrophages influence COPD, encompassing inflammatory mediator production, protease/antiprotease release, and phagocytic activity. Additionally, we synthesize findings from available literature regarding all identified pulmonary macrophage sub-populations in COPD, thereby advancing our comprehension of macrophage heterogeneity's significance in the complex pathophysiological mechanisms of COPD.

Robust estimation of cancer and immune cell-type proportions from bulk tumor ATAC-Seq data

Assay for Transposase-Accessible Chromatin sequencing (ATAC-Seq) is a widely used technique to explore gene regulatory mechanisms. For most ATAC-Seq data from healthy and diseased tissues such as tumors, chromatin accessibility measurement represents a mixed signal from multiple cell types. In this work, we derive reliable chromatin accessibility marker peaks and reference profiles for most non-malignant cell types frequently observed in the microenvironment of human tumors. We then integrate these data into the EPIC deconvolution framework (Racle et al., 2017) to quantify cell-type heterogeneity in bulk ATAC-Seq data. Our EPIC-ATAC tool accurately predicts non-malignant and malignant cell fractions in tumor samples. When applied to a human breast cancer cohort, EPIC-ATAC accurately infers the immune contexture of the main breast cancer subtypes.

Integrative Analysis of ATAC-Seq and RNA-Seq through Machine Learning Identifies 10 Signature Genes for Breast Cancer Intrinsic Subtypes

Breast cancer is a heterogeneous disease composed of various biologically distinct subtypes, each characterized by unique molecular features. Its formation and progression involve a complex, multistep process that includes the accumulation of numerous genetic and epigenetic alterations. Although integrating RNA-seq transcriptome data with ATAC-seq epigenetic information provides a more comprehensive understanding of gene regulation and its impact across different conditions, no classification model has yet been developed for breast cancer intrinsic subtypes based on such integrative analyses. In this study, we employed machine learning algorithms to predict intrinsic subtypes through the integrative analysis of ATAC-seq and RNA-seq data. We identified 10 signature genes (CDH3, ERBB2, TYMS, GREB1, OSR1, MYBL2, FAM83D, ESR1, FOXC1, and NAT1) using recursive feature elimination with cross-validation (RFECV) and a support vector machine (SVM) based on SHAP (SHapley Additive exPlanations) feature importance. Furthermore, we found that these genes were primarily associated with immune responses, hormone signaling, cancer progression, and cellular proliferation.

Novel prognostic alternative splicing events in colorectal Cancer: Impact on immune infiltration and therapy response

OBJECTIVE

This study aims to comprehensively analyze alternative splicing (AS) features in colorectal cancer (CRC) using integrative multi-omics and to elucidate their relationship with the CRC immune microenvironment.

METHODS

Transcriptomic data, clinical information, and Percent Spliced In (PSI) values of AS events for CRC patients were obtained from The Cancer Genome Atlas (TCGA) and TCGA SpliceSeq databases. Differentially expressed AS events were identified. Univariate Cox analysis was used to pinpoint prognosis-related AS events. A prognostic risk model was developed and validated using multivariate Cox analysis, patient survival analysis, and the area under the receiver operating characteristic (ROC) curve (AUC). Gene Set Enrichment Analysis (GSEA), immune infiltration, immunotherapy, chemotherapy sensitivity analyses, and regulatory relationships between AS events and splicing factors (SFs) were conducted. Single-cell sequencing was used to study the distribution of key factors. siRNA and overexpression vectors were utilized to silence/overexpress BCAS1 in CRC cells and evaluate their effects on cell growth, migration, and invasion. Furthermore, the interaction between BCAS1 and ANO7 pre-mRNA was investigated using RIP-PCR.

RESULTS

82 prognosis-related AS events were identified in CRC patients. A 15-AS prognostic model was constructed, which correlated with immune cell infiltration and showed differences in immunotherapy and chemotherapy sensitivity. BCAS1 was identified as a potential regulator of the ANO7|58341|AT splicing event in CRC. Single-cell sequencing analysis revealed the distribution of BCAS1 and ANO7 in cancer stem cells. In vitro experiments demonstrated that overexpression of BCAS1 and silencing of ANO7 inhibit the proliferation, migration, and invasion of CRC cells. Moreover, BCAS1 suppresses the progression of CRC by modulating ANO7 alternative splicing.

CONCLUSION

This study provides new insights into the role of alternative splicing in colorectal cancer, particularly the potential regulatory action of BCAS1 on the ANO7|58341|AT splicing event. It also identifies the impact of alternative splicing on the tumor microenvironment and potential implications for immunotherapy, highlighting its relevance for the in-depth study and treatment of CRC.

Deconvoluting multi-person biological mixtures and accurate characterization and identification of separated contributors using non-targeted single-cell DNA sequencing

The genetic characterization and identification of individuals who contributed to biological mixtures are complex and mostly unresolved tasks. These tasks are relevant in various fields, particularly in forensic investigations, which frequently encounters crime scene stains generated by more than one person. Currently, forensic mixture deconvolution is mostly performed subsequent to forensic DNA profiling at the level of the mixed DNA profiles, which comes with several limitations. Some previous studies attempted at separating single cells prior to forensic DNA profiling. However, these approaches are biased at selection of the cells and, due to their targeted DNA analysis on low template DNA, provide incomplete and unreliable forensic DNA profiles. We recently demonstrated the feasibility of performing mixture deconvolution prior to forensic DNA profiling through the utilization of a non-targeted single-cell transcriptome sequencing (scRNA-seq). In addition to individual-specific mixture deconvolution, this approach also allowed accurate characterisation of biological sex, biogeographic ancestry and individual identification of the separated mixture contributors. However, RNA has the forensic disadvantage of being prone to degradation, and sequencing RNA - focussing on coding regions - limits the number of single nucleotide polymorphisms (SNPs) utilized for genetic mixture deconvolution, characterization, and identification. These limitations can be overcome by performing single-cell sequencing on the level of DNA instead of RNA. Here, for the first time, we applied non-targeted single-cell DNA sequencing (scDNA-seq) by applying the scATAC-seq (Assay for Transposase-Accessible Chromatin with sequencing) technique to address the challenges of mixture deconvolution in the forensic context. We demonstrated that scATAC-seq, together with our recently developed De-goulash data analysis pipeline, is capable of deconvoluting complex blood mixtures of five individuals from both sexes with varying biogeographic ancestries. We further showed that our approach achieved correct genetic characterization of the biological sex and the biogeographic ancestry of each of the separated mixture contributors and established their identity. Furthermore, by analysing in-silico generated scATAC-seq data mixtures, we demonstrated successful individual-specific mixture deconvolution of i) highly complex mixtures of 11 individuals, ii) balanced mixtures containing as few as 20 cells (10 per each individual), and iii) imbalanced mixtures with a ratio as low as 1:80. Overall, our proof-of-principle study demonstrates the general feasibility of scDNA-seq in general, and scATAC-seq in particular, for mixture deconvolution, genetic characterization and individual identification of the separated mixture contributors. Furthermore, it shows that compared to scRNA-seq, scDNA-seq detects more SNPs from fewer cells, providing higher sensitivity, that is valuable in forensic genetics.

Next-generation sequencing strategies in venous thromboembolism: in whom and for what purpose?

This invited review follows the oral presentation "To Sequence or Not to Sequence, That Is Not the Question; But 'When, Who, Which and What For?' Is" given during the State of the Art session "Translational Genomics in Thrombosis: From OMICs to Clinics" of the International Society on Thrombosis and Haemostasis 2023 Congress. Emphasizing the power of next-generation sequencing technologies and the diverse strategies associated with DNA variant analysis, this review highlights the unresolved questions and challenges in their implementation both for the clinical diagnosis of venous thromboembolism and in translational research.

Cactus: A user-friendly and reproducible ATAC-Seq and mRNA-Seq analysis pipeline for data preprocessing, differential analysis, and enrichment analysis

The ever decreasing cost of Next-Generation Sequencing coupled with the emergence of efficient and reproducible analysis pipelines has rendered genomic methods more accessible. However, downstream analyses are basic or missing in most workflows, creating a significant barrier for non-bioinformaticians. To help close this gap, we developed Cactus, an end-to-end pipeline for analyzing ATAC-Seq and mRNA-Seq data, either separately or jointly. Its Nextflow-, container-, and virtual environment-based architecture ensures efficient and reproducible analyses. Cactus preprocesses raw reads, conducts differential analyses between conditions, and performs enrichment analyses in various databases, including DNA-binding motifs, ChIP-Seq binding sites, chromatin states, and ontologies. We demonstrate the utility of Cactus in a multi-modal and multi-species case study as well as by showcasing its unique capabilities as compared to other ATAC-Seq pipelines. In conclusion, Cactus can assist researchers in gaining comprehensive insights from chromatin accessibility and gene expression data in a quick, user-friendly, and reproducible manner.

Androgen signaling inhibits de novo lipogenesis to alleviate lipid deposition in zebrafish

Testosterone is closely associated with lipid metabolism and known to affect body fat composition and muscle mass in males. However, the mechanisms by which testosterone acts on lipid metabolism are not yet fully understood, especially in teleosts. In this study, cyp17a1-/- zebrafish ( Danio rerio) exhibited excessive visceral adipose tissue (VAT), lipid content, and up-regulated expression and activity of hepatic de novo lipogenesis (DNL) enzymes. The assay for transposase accessible chromatin with sequencing (ATAC-seq) results demonstrated that chromatin accessibility of DNL genes was increased in cyp17a1-/- fish compared to cyp17a1+/+ male fish, including stearoyl-CoA desaturase ( scd) and fatty acid synthase ( fasn). Androgen response element (ARE) motifs in the androgen signaling pathway were significantly enriched in cyp17a1+/+ male fish but not in cyp17a1-/- fish. Both androgen receptor ( ar)-/- and wild-type (WT) zebrafish administered with Ar antagonist flutamide displayed excessive visceral adipose tissue, lipid content, and up-regulated expression and activity of hepatic de novo lipogenesis enzymes. The Ar agonist BMS-564929 reduced the content of VAT and lipid content, and down-regulated acetyl-CoA carboxylase a ( acaca), fasn, and scd expression. Mechanistically, the rescue effect of testosterone on cyp17a1-/- fish in terms of phenotypes was abolished when ar was additionally depleted. Collectively, these findings reveal that testosterone inhibits lipid deposition by down-regulating DNL genes via Ar in zebrafish, thus expanding our understanding of the relationship between testosterone and lipid metabolism in teleosts.

Dynamics of Chromatin Opening across Larval Development in the Urochordate Ascidian Ciona savignyi

Ascidian larvae undergo tail elongation and notochord lumenogenesis, making them an ideal model for investigating tissue morphogenesis in embryogenesis. The cellular and mechanical mechanisms of these processes have been studied; however, the underlying molecular regulatory mechanism remains to be elucidated. In this study, assays for transposase-accessible chromatin using sequencing (ATAC-seq) and RNA sequencing (RNA-seq) were applied to investigate potential regulators of the development of ascidian Ciona savignyi larvae. Our results revealed 351 and 138 differentially accessible region genes through comparisons of ATAC-seq data between stages 21 and 24 and between stages 24 and 25, respectively. A joint analysis of RNA-seq and ATAC-seq data revealed a correlation between chromatin accessibility and gene transcription. We further verified the tissue expression patterns of 12 different genes. Among them, Cs-matrix metalloproteinase 24 (MMP24) and Cs-krüppel-like factor 5 (KLF5) were highly expressed in notochord cells. Functional assay results demonstrated that both genes are necessary for notochord lumen formation and expansion. Finally, we performed motif enrichment analysis of the differentially accessible regions in different tailbud stages and summarized the potential roles of these motif-bearing transcription factors in larval development. Overall, our study found a correlation between gene expression and chromatin accessibility and provided a vital resource for understanding the mechanisms of the development of ascidian embryos.

Comprehensive analysis of PSMD family members and validation of PSMD9 as a potential therapeutic target in human glioblastoma

AIMS

PSMD family members, as important components of the 26S proteasome, are well known to be involved in protein degradation. However, their role in glioblastoma (GBM) has not been rigorously investigated. We aimed to perform systematic analysis of the expression signature, prognostic significance and functions of PSMD family genes in GBM to reveal potential prognostic markers and new therapeutic targets among PSMD family members.

METHODS

In this study, we systemically analyzed PSMD family members in terms of their expression profiles, prognostic implications, DNA methylation levels, and genetic alterations; the relationships between their expression levels and immune infiltration and drug sensitivity; and their potential functional enrichment in GBM through bioinformatics assessment. Moreover, in vitro and in vivo experiments were used to validate the biological functions of PSMD9 and its targeted therapeutic effect in GBM.

RESULTS

The mRNA levels of PSMD5/8/9/10/11/13/14 were higher in GBM than in normal brain tissues, and the mRNA levels of PSMD1/4/5/8/9/11/12 were higher in high-grade glioma (WHO grade III & IV) than in low-grade glioma (WHO grade II). High mRNA expression of PSMD2/6/8/9/12/13/14 and low mRNA expression of PSMD7 were associated with poor overall survival (OS). Multivariate Cox regression analysis identified PSMD2/5/6/8/9/10/11/12 as independent prognostic factors for OS prediction. In addition, the protein-protein interaction network and gene set enrichment analysis results suggested that PSMD family members and their interacting molecules were involved in the regulation of the cell cycle, cell invasion and migration, and other biological processes in GBM. In addition, knockdown of PSMD9 inhibited cell proliferation, invasion and migration and induced G2/M cell cycle arrest in LN229 and A172 GBM cells. Moreover, PSMD9 promoted the malignant progression of GBM in vivo. GBM cell lines with high PSMD9 expression were more resistant to panobinostat, a potent deacetylase inhibitor, than those with low PSMD9 expression. In vitro and in vivo experiments further validated that PSMD9 overexpression rescued the GBM inhibitory effect of panobinostat.

CONCLUSION

This study provides new insights into the value of the PSMD family in human GBM diagnosis and prognosis evaluation, and we further identified PSMD9 as a potential therapeutic target. These findings may lead to the development of effective therapeutic strategies for GBM.

Dynamic chromatin accessibility landscapes of osteoblast differentiation and mineralization

Bone acts as a self-healing organ, which undergoes continuous regeneration process that is tightly regulated by the cooperation of osteoclasts with the capability of bone resorption and osteoblasts with the capability of bone formation. Generally, bone marrow derived mesenchymal stem cells (BMSCs) differentiated to final osteoblasts have been considered as critical role in bone remodeling. In this regard, several transcription factors (TFs) whose binding sites are initially hidden deep within accessible chromatin that participate in modulating osteoblast differentiation and bone matrix mineralization. Then, it is necessary to explore further the dynamic changes about the epigenetic transcription machinery during osteoblastogenesis. Here, we performed the chromatin accessibility and transcriptomic landscape of osteoblast differentiation and mineralization by using transposase-accessible chromatin sequencing (ATAC-seq) and RNA sequencing (RNA-Seq). Our data found that global chromatin accessibility during osteoblastogenesis was extensively improved. Above this, it is shown that key target genes including Col6a3, Serpina3n, Ms4a4d, Lyz2, Phf11b and Grin3a were enriched in differential loci RNA-seq and ATAC-Seq peaks with continuous changed tendency during osteoblasts differentiation and mineralization. In addition, Analysis of Motif Enrichment (AME) was used to elucidate TFs which modulated these target genes. In this study, it was shown for the first time that these important TFs including MEF2A, PRRX1, Shox2 and HOXB13 could alter promoter accessibility of target genes during osteoblastogenesis. This helps us understand how TF binding motif accessibility influences osteoblast differentiation. In addition, it also suggests that modulating the chromatin accessibility of osteogenesis could be developed as the promising strategies to regulate bone regeneration.

3D genome structural variations play important roles in regulating seed oil content of Brassica napus

Dissecting the complex regulatory mechanism of seed oil content (SOC) is one of the main research goals in Brassica napus. Increasing evidence suggests that genome architecture is linked to multiple biological functions. However, the effect of genome architecture on SOC regulation remains unclear. Here, we used high-throughput chromatin conformation capture to characterize differences in the three-dimensional (3D) landscape of genome architecture of seeds from two B. napus lines, N53-2 (with high SOC) and Ken-C8 (with low SOC). Bioinformatics analysis demonstrated that differentially accessible regions and differentially expressed genes between N53-2 and Ken-C8 were preferentially enriched in regions with quantitative trait loci (QTLs)/associated genomic regions (AGRs) for SOC. A multi-omics analysis demonstrated that expression of SOC-related genes was tightly correlated with genome structural variations in QTLs/AGRs of B. napus. The candidate gene BnaA09g48250D, which showed structural variation in a QTL/AGR on chrA09, was identified by fine-mapping of a KN double-haploid population derived from hybridization of N53-2 and Ken-C8. Overexpression and knockout of BnaA09g48250D led to significant increases and decreases in SOC, respectively, in the transgenic lines. Taken together, our results reveal the 3D genome architecture of B. napus seeds and the roles of genome structural variations in SOC regulation, enriching our understanding of the molecular mechanisms of SOC regulation from the perspective of spatial chromatin structure.

Dissecting the Immune System through Gene Regulation

The immune system plays a dual role in human health, functioning both as a protector against pathogens and, at times, as a contributor to disease. This feature emphasizes the importance to uncover the underlying causes of its malfunctions, necessitating an in-depth analysis in both pathological and physiological conditions to better understand the immune system and immune disorders. Recent advances in scientific technology have enabled extensive investigations into gene regulation, a crucial mechanism governing cellular functionality. Studying gene regulatory mechanisms within the immune system is a promising avenue for enhancing our understanding of immune cells and the immune system as a whole. The gene regulatory mechanisms, revealed through various methodologies, and their implications in the field of immunology are discussed in this chapter.

Analysis of Chromatin Accessibility Changes Induced by BMMC Recognition of Foot-and-Mouth Disease Virus-like Particles through ATAC-seq

Mast cells can recognize foot-and-mouth disease virus-like particles (FMDV-VLPs) via mannose receptors (MRs) to produce differentially expressed cytokines. The regulatory role of chromatin accessibility in this process is unclear. Bone marrow-derived mast cells (BMMCs) were cultured, and an assay of transposase-accessible chromatin sequencing (ATAC-seq) was applied to demonstrate the regulation of chromatin accessibility in response to the BMMCs' recognition of FMDV-VLPs. A pathway enrichment analysis showed that peaks associated with the nuclear factor-κB (NF-κB), mitogen-activated protein kinase (MAPK), phosphatidylinositol 3 kinase-protein kinase B (PI3K-Akt), and other signaling pathways, especially the NF-κB pathway, were involved in the BMMCs' recognition of VLPs. Moreover, transcription factors including SP1, NRF1, AP1, GATA3, microphthalmia-associated transcription factor (MITF), and NF-κB-p65 may bind to the motifs with altered chromatin accessibility to regulate gene transcription. Furthermore, the expression of NF-κB, interleukin (IL)-9, tumor necrosis factor (TNF)-α, and interferon (IFN)-γ in the BMMCs of the VLP group increased compared with that of the BMMCs in the control group, whereas the expression of IL-10 did not differ significantly between groups. After inhibiting the MRs, the expression of NF-κB, IL-9, TNF-α, and IFN-γ decreased significantly, whereas the expression of IL-10 increased. The expression of MAPK and IL-6 showed no significant change after MR inhibition. This study demonstrated that MRs expressed on BMMCs can affect the NF-κB pathway by changing chromatin accessibility to regulate the transcription of specific cytokines, ultimately leading to the differential expression of cytokines. These data provide a theoretical basis and new ideas for the development of a novel vaccine for FMD.

Prognostic hub gene CBX2 drives a cancer stem cell-like phenotype in HCC revealed by multi-omics and multi-cohorts

Hepatocellular carcinoma (HCC) is a malignant tumor with a high prevalence and fatality rate. CBX2 has been demonstrated to impact the development and advancement of various cancers, albeit it has received limited attention in relation to HCC. In this study, CBX2 and CEP55 were screened out with the refined triple regulatory networks constructed by total RNA-seq datasets (TCGA-LIHC, GSE140845) and a robust prognostic model. Aberrantly higher expression levels of CBX2 and CEP55 in HCC may be caused by CNV alterations, promoter hypo-methylation, open chromatin accessibility, and greater active marks such as H3K4me3, H3K4me1, and H3K27ac. Functionally, CBX2, which was highly correlated with CD44, shaped a cancer stem cell-like phenotype by positively regulating cell-cycle progression, proliferation, invasion, metastasis, wound healing, and radiation resistance, revealed by combining bulk RNA-seq and scRNA-seq datasets. CBX2 knockdown validated its role in affecting the cell cycle. Importantly, we revealed CBX2 could activate gene by cooperating with co-regulators or not rather than a recognizer of the repressive mark H3K27me3. For instance, we uncovered CBX2 bound to promoter of CTNNB1 and CEP55 to augment their expressions. CBX2 showed a highly positive correlation with CEP55 at pan-cancer level. In addition, CBX2 and CEP55 may enhance extracellular matrix reprograming via cancer-associated fibroblast. Surprisingly, patients with high expression of CBX2 or CEP55 exhibited a higher response to immunotherapy, indicating that CBX2 and CEP55 may be promising therapeutic targets for HCC patients.

ATF4 renders human T-cell acute lymphoblastic leukemia cell resistance to FGFR1 inhibitors through amino acid metabolic reprogramming

Abnormalities of FGFR1 have been reported in multiple malignancies, suggesting FGFR1 as a potential target for precision treatment, but drug resistance remains a formidable obstacle. In this study, we explored whether FGFR1 acted a therapeutic target in human T-cell acute lymphoblastic leukemia (T-ALL) and the molecular mechanisms underlying T-ALL cell resistance to FGFR1 inhibitors. We showed that FGFR1 was significantly upregulated in human T-ALL and inversely correlated with the prognosis of patients. Knockdown of FGFR1 suppressed T-ALL growth and progression both in vitro and in vivo. However, the T-ALL cells were resistant to FGFR1 inhibitors AZD4547 and PD-166866 even though FGFR1 signaling was specifically inhibited in the early stage. Mechanistically, we found that FGFR1 inhibitors markedly increased the expression of ATF4, which was a major initiator for T-ALL resistance to FGFR1 inhibitors. We further revealed that FGFR1 inhibitors induced expression of ATF4 through enhancing chromatin accessibility combined with translational activation via the GCN2-eIF2α pathway. Subsequently, ATF4 remodeled the amino acid metabolism by stimulating the expression of multiple metabolic genes ASNS, ASS1, PHGDH and SLC1A5, maintaining the activation of mTORC1, which contributed to the drug resistance in T-ALL cells. Targeting FGFR1 and mTOR exhibited synergistically anti-leukemic efficacy. These results reveal that FGFR1 is a potential therapeutic target in human T-ALL, and ATF4-mediated amino acid metabolic reprogramming contributes to the FGFR1 inhibitor resistance. Synergistically inhibiting FGFR1 and mTOR can overcome this obstacle in T-ALL therapy.

An Integrative ATAC-Seq and RNA-Seq Analysis of the Endometrial Tissues of Meishan and Duroc Pigs

Meishan pigs are a well-known indigenous pig breed in China characterized by a high fertility. Notably, the number of endometrial grands is significantly higher in Meishan pigs than Duroc pigs. The characteristics of the endometrial tissue are related to litter size. Therefore, we used the assay for transposase-accessible chromatin with sequencing (ATAC-seq) and RNA-sequencing (RNA-seq) to analyze the mechanisms underlying the differences in fecundity between the breeds. We detected the key transcription factors, including Double homeobox (Dux), Ladybird-like homeobox gene 2 (LBX2), and LIM homeobox 8 (Lhx8), with potentially pivotal roles in the regulation of the genes related to endometrial development. We identified the differentially expressed genes between the breeds, including SOX17, ANXA4, DLX3, DMRT1, FLNB, IRF6, CBFA2T2, TFCP2L1, EFNA5, SLIT2, and CYFIP2, with roles in epithelial cell differentiation, fertility, and ovulation. Interestingly, ANXA4, CBFA2T2, and TFCP2L1, which were upregulated in the Meishan pigs in the RNA-seq analysis, were identified again by the integration of the ATAC-seq and RNA-seq data. Moreover, we identified genes in the cancer or immune pathways, FoxO signaling, Wnt signaling, and phospholipase D signaling pathways. These ATAC-seq and RNA-seq analyses revealed the accessible chromatin and potential mechanisms underlying the differences in the endometrial tissues between the two types of pigs.

Factors and Methods for the Detection of Gene Expression Regulation

Gene-expression regulation involves multiple processes and a range of regulatory factors. In this review, we describe the key factors that regulate gene expression, including transcription factors (TFs), chromatin accessibility, histone modifications, DNA methylation, and RNA modifications. In addition, we also describe methods that can be used to detect these regulatory factors.

Sample Preparation and Integrative Data Analysis of a Droplet-based Single-Cell ATAC-sequencing Using Murine Thymic Epithelial Cells

Accessible chromatin regions modulate gene expression by acting as cis-regulatory elements. Understanding the epigenetic landscape by mapping accessible regions of DNA is therefore imperative to decipher mechanisms of gene regulation under specific biological contexts of interest. The assay for transposase-accessible chromatin sequencing (ATAC-seq) has been widely used to detect accessible chromatin and the recent introduction of single-cell technology has increased resolution to the single-cell level. In a recent study, we used droplet-based, single-cell ATAC-seq technology (scATAC-seq) to reveal the epigenetic profile of the transit-amplifying subset of thymic epithelial cells (TECs), which was identified previously using single-cell RNA-sequencing technology (scRNA-seq). This protocol allows the preparation of nuclei from TECs in order to perform droplet-based scATAC-seq and its integrative analysis with scRNA-seq data obtained from the same cell population. Integrative analysis has the advantage of identifying cell types in scATAC-seq data based on cell cluster annotations in scRNA-seq analysis.

From single- to multi-omics: future research trends in medicinal plants

Medicinal plants are the main source of natural metabolites with specialised pharmacological activities and have been widely examined by plant researchers. Numerous omics studies of medicinal plants have been performed to identify molecular markers of species and functional genes controlling key biological traits, as well as to understand biosynthetic pathways of bioactive metabolites and the regulatory mechanisms of environmental responses. Omics technologies have been widely applied to medicinal plants, including as taxonomics, transcriptomics, metabolomics, proteomics, genomics, pangenomics, epigenomics and mutagenomics. However, because of the complex biological regulation network, single omics usually fail to explain the specific biological phenomena. In recent years, reports of integrated multi-omics studies of medicinal plants have increased. Until now, there have few assessments of recent developments and upcoming trends in omics studies of medicinal plants. We highlight recent developments in omics research of medicinal plants, summarise the typical bioinformatics resources available for analysing omics datasets, and discuss related future directions and challenges. This information facilitates further studies of medicinal plants, refinement of current approaches and leads to new ideas.

Integrated analyses of the methylome and transcriptome to unravel sex differences in the perirenal fat from suckling lambs

In sheep, differences were observed regarding fat accumulation and fatty acid (FA) composition between males and females, which may impact the quality and organoleptic characteristics of the meat. The integration of different omics technologies is a relevant approach for investigating biological and genetic mechanisms associated with complex traits. Here, the perirenal tissue of six male and six female Assaf suckling lambs was evaluated using RNA sequencing and whole-genome bisulfite sequencing (WGBS). A multiomic discriminant analysis using multiblock (s)PLS-DA allowed the identification of 314 genes and 627 differentially methylated regions (within these genes), which perfectly discriminate between males and females. These candidate genes overlapped with previously reported QTLs for carcass fat volume and percentage of different FAs in milk and meat from sheep. Additionally, differentially coexpressed (DcoExp) modules of genes between males (nine) and females (three) were identified that harbour 22 of these selected genes. Interestingly, these DcoExp were significantly correlated with fat percentage in different deposits (renal, pelvic, subcutaneous and intramuscular) and were associated with relevant biological processes for adipogenesis, adipocyte differentiation, fat volume and FA composition. Consequently, these genes may potentially impact adiposity and meat quality traits in a sex-specific manner, such as juiciness, tenderness and flavour.

Integration of RNA-seq and ATAC-seq identifies muscle-regulated hub genes in cattle

As the main product of livestock, muscle itself plays an irreplaceable role in maintaining animal body movement and regulating metabolism. Therefore, it is of great significance to explore its growth, development and regeneration to improve the meat yield and quality of livestock. In this study, we attempted to use RNA-seq and ATAC-seq techniques to identify differentially expressed genes (DEGs) specifically expressed in bovine skeletal muscle as potential candidates for studying the regulatory mechanisms of muscle development. Microarray data from 8 tissue samples were selected from the GEO database for analysis. First, we obtained gene modules related to each tissue through WGCNA analysis. Through Gene Ontology (GO) functional annotation, the module of lightyellow (ME_lightyellow) was closely related to muscle development, and 213 hub genes were screened as follow-up research targets. Further, the difference analysis showed that, except for PREB, all other candidate hub genes were up-regulated (muscle group vs. other-group). ATAC-seq analysis showed that muscle-specific accessible chromatin regions were mainly located in promoter of genes related to muscle structure development (GO:0061061), muscle cell development (GO:0055001) and muscle system process (GO:0003012), which were involved in cAMP, CGMP-PKG, MAPK, and other signaling pathways. Next, we integrated the results of RNA-seq and ATAC-seq analysis, and 54 of the 212 candidate hub genes were identified as key regulatory genes in skeletal muscle development. Finally, through motif analysis, 22 of the 54 key genes were found to be potential target genes of transcription factor MEF2C. Including CAPN3, ACTN2, MB, MYOM3, SRL, CKM, ALPK3, MAP3K20, UBE2G1, NEURL2, CAND2, DOT1L, HRC, MAMSTR, FSD2, LRRC2, LSMEM1, SLC29A2, FHL3, KLHL41, ATXN7L2, and PDRG1. This provides a potential reference for studying the molecular mechanism of skeletal muscle development in mammals.