Chromatin accessibility profiling by ATAC-seq

Dual-role transcription factors stabilize intermediate expression levels

Precise control of gene expression levels is essential for normal cell functions, yet how they are defined and tightly maintained, particularly at intermediate levels, remains elusive. Here, using a series of newly developed sequencing, imaging, and functional assays, we uncover a class of transcription factors with dual roles as activators and repressors, referred to as condensate-forming level-regulating dual-action transcription factors (TFs). They reduce high expression but increase low expression to achieve stable intermediate levels. Dual-action TFs directly exert activating and repressing functions via condensate-forming domains that compartmentalize core transcriptional unit selectively. Clinically relevant mutations in these domains, which are linked to a range of developmental disorders, impair condensate selectivity and dual-action TF activity. These results collectively address a fundamental question in expression regulation and demonstrate the potential of level-regulating dual-action TFs as powerful effectors for engineering controlled expression levels.

Single-cell technologies in psoriasis

Psoriasis is a chronic and recurrent inflammatory skin disorder. The primary manifestation of psoriasis arises from disturbances in the cutaneous immune microenvironment, but the specific functions of the cellular components within this microenvironment remain unknown. Recent advancements in single-cell technologies have enabled the detection of multi-omics at the level of individual cells, including single-cell transcriptome, proteome, and metabolome, which have been successfully applied in studying autoimmune diseases, and other pathologies. These techniques allow the identification of heterogeneous cell clusters and their varying contributions to disease development. Considering the immunological traits of psoriasis, an in-depth exploration of immune cells and their interactions with cutaneous parenchymal cells can markedly advance our comprehension of the mechanisms underlying the onset and recurrence of psoriasis. In this comprehensive review, we present an overview of recent applications of single-cell technologies in psoriasis, aiming to improve our understanding of the underlying mechanisms of this disorder.

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.

Primary microglia cell cultures in translational research: Strengths and limitations

Microglia are the resident macrophages in the central nervous system, accounting for 10-15% of the cell mass in the brain. Next to their physiological role in development, monitoring neuronal function and the maintenance of homeostasis, microglia are crucial in the brain's immune defense. Brain injury and chronic neurological disorders are associated with neuroinflammation, in which microglia activation is a central element. Microglia acquire a wide spectrum of activation states in the diseased or injured brain, some of which are neurotoxic. The investigation of microglia (patho)physiology and therapeutic interventions targeting neuroinflammation is a substantial challenge. In addition to in vivo approaches, the application of in vitro model systems has gained significant ground and is essential to complement in vivo work. Primary microglia cultures have proved to be a useful tool. Microglia cultures have offered the opportunity to explore the mechanistic, molecular elements of microglia activation, the microglia secretome, and the efficacy of therapeutic treatments against neuroinflammation. As all model systems, primary microglia cultures have distinct strengths and limitations to be weighed when experiments are designed and when data are interpreted. Here, we set out to provide a succinct overview of the advantages and pitfalls of the use of microglia cultures, which instructs the refinement and further development of this technique to remain useful in the toolbox of microglia researchers. Since there is no conclusive therapy to combat neurotoxicity linked to neuroinflammation in acute brain injury or neurodegenerative disorders, these research tools remain essential to explore therapeutic opportunities.

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.

Untangling the gordian knot: The intertwining interactions between developmental hormone signaling and epigenetic mechanisms in insects

Hormones control developmental and physiological processes, often by regulating the expression of multiple genes simultaneously or sequentially. Crosstalk between hormones and epigenetics is pivotal to dynamically coordinate this process. Hormonal signals can guide the addition and removal of epigenetic marks, steering gene expression. Conversely, DNA methylation, histone modifications and non-coding RNAs can modulate regional chromatin structure and accessibility and regulate the expression of numerous (hormone-related) genes. Here, we provide a review of the interplay between the classical insect hormones, ecdysteroids and juvenile hormones, and epigenetics. We summarize the mode-of-action and roles of these hormones in post-embryonic development, and provide a general overview of epigenetic mechanisms. We then highlight recent advances on the interactions between these hormonal pathways and epigenetics, and their involvement in development. Furthermore, we give an overview of several 'omics techniques employed in the field. Finally, we discuss which questions remain unanswered and possible avenues for future research.

Integrated analysis of ATAC-seq and transcriptomic reveals the ScDof3-ScproC molecular module regulating the cold acclimation capacity of potato

Low temperature severely affects the geographical distribution and production of potato, which may incur cold damage in early spring or winter. Cultivated potatoes, mainly derived from Solanum tuberosum, are sensitive to freezing stress, but wild species of potato such as S. commersonii exhibit both constitutive freezing tolerance and/or cold acclimation tolerance. Hence, such wild species could assist in cold hardiness breeding. Yet the key transcription factors and their downstream functional genes that confer freezing tolerance are far from clear, hindering the breeding process. Here, we used ATAC-seq (Assay for Transposase-Accessible Chromatin with high-throughput sequencing) alongside RNA-seq to investigate the variation in chromatin accessibility and patterns of gene expression in freezing-tolerant CMM5 (S. commersonii), before and after its cold treatment. Our results suggest that after exposure to cold, transcription factors including Dof3, ABF2, PIF4, and MYB4 were predicted to further control the genes active in the synthetic/metabolic pathways of plant hormones, namely abscisic acid, polyamine, and reductive glutathione (among others). This suggests these transcription factors could regulate freezing tolerance of CMM5 leaves. In particular, ScDof3 was proven to regulate the expression of ScproC (pyrroline-5-carboxylate reductase, P5CR) according to dual-LUC assays. Overexpressing ScDof3 in Nicotiana benthamiana leaves led to an increase in both the proline content and expression level of NbproC (homolog of ScproC). These results demonstrate the ScDof3-ScproC module regulates the proline content and thus promotes freezing tolerance in potato. Our research provides valuable genetic resources to further study the molecular mechanisms underpinning cold tolerance in potato.

Integrative ATAC-seq and RNA-seq analysis of myogenic differentiation of ovine skeletal muscle satellite cell

Skeletal muscle satellite cells (SMSCs) play an important role in regulating muscle growth and regeneration. Chromatin accessibility allows physical interactions that synergistically regulate gene expression through enhancers, promoters, insulators, and chromatin binding factors. However, the chromatin accessibility altas and its regulatory role in ovine myoblast differentiation is still unclear. Therefore, ATAC-seq and RNA-seq analysis were performed on ovine SMSCs at the proliferation stage (SCG) and differentiation stage (SCD). 17,460 DARs (differential accessibility regions) and 3732 DEGs (differentially expressed genes) were identified. Based on joint analysis of ATAC-seq and RNA-seq, we revealed that PI3K-Akt, TGF-β and other signaling pathways regulated SMSCs differentiation. We identified two novel candidate genes, FZD5 and MAP2K6, which may affect the proliferation and differentiation of SMSCs. Our data identify potential cis regulatory elements of ovine SMSCs. This study can provide a reference for exploring the mechanisms of the differentiation and regeneration of SMSCs in the future.

PBX/Knotted 1 homeobox-2 (PKNOX2) is a novel regulator of myocardial fibrosis

Much effort has been made to uncover the cellular heterogeneities of human hearts by single-nucleus RNA sequencing. However, the cardiac transcriptional regulation networks have not been systematically described because of the limitations in detecting transcription factors. In this study, we optimized a pipeline for isolating nuclei and conducting single-nucleus RNA sequencing targeted to detect a higher number of cell signal genes and an optimal number of transcription factors. With this unbiased protocol, we characterized the cellular composition of healthy human hearts and investigated the transcriptional regulation networks involved in determining the cellular identities and functions of the main cardiac cell subtypes. Particularly in fibroblasts, a novel regulator, PKNOX2, was identified as being associated with physiological fibroblast activation in healthy hearts. To validate the roles of these transcription factors in maintaining homeostasis, we used single-nucleus RNA-sequencing analysis of transplanted failing hearts focusing on fibroblast remodelling. The trajectory analysis suggested that PKNOX2 was abnormally decreased from fibroblast activation to pathological myofibroblast formation. Both gain- and loss-of-function in vitro experiments demonstrated the inhibitory role of PKNOX2 in pathological fibrosis remodelling. Moreover, fibroblast-specific overexpression and knockout of PKNOX2 in a heart failure mouse model induced by transverse aortic constriction surgery significantly improved and aggravated myocardial fibrosis, respectively. In summary, this study established a high-quality pipeline for single-nucleus RNA-sequencing analysis of heart muscle. With this optimized protocol, we described the transcriptional regulation networks of the main cardiac cell subtypes and identified PKNOX2 as a novel regulator in suppressing fibrosis and a potential therapeutic target for future translational studies.

Gfi1 controls the formation of effector CD8 T cells during chronic infection and cancer

During chronic infections and tumor progression, CD8 T cells gradually lose their effector functions and become exhausted. These exhausted CD8 T cells are heterogeneous and comprised of different subsets, including self-renewing progenitors that give rise to Ly108 - CX3CR1 + effector-like cells. Generation of these effector-like cells is essential for the control of chronic infections and tumors, albeit limited. However, the precise cues and mechanisms directing the formation and maintenance of exhausted effector-like are incompletely understood. Using genetic mouse models challenged with LCMV Clone 13 or syngeneic tumors, we show that the expression of a transcriptional repressor, growth factor independent 1 (Gfi1) is dynamically regulated in exhausted CD8 T cells, which in turn regulates the formation of exhausted effector-like cells. Gfi1 deletion in T cells dysregulates the chromatin accessibility and transcriptomic programs associated with the differentiation of LCMV Clone 13-specific CD8 T cell exhaustion, preventing the formation of effector-like and terminally exhausted cells while maintaining progenitors and a newly identified Ly108 + CX3CR1 + state. These Ly108 + CX3CR1 + cells have a distinct chromatin profile and may represent an alternative target for therapeutic interventions to combat chronic infections and cancer. In sum, we show that Gfi1 is a critical regulator of the formation of exhausted effector-like cells.

Inference of genomic landscapes using ordered Hidden Markov Models with emission densities (oHMMed)

BACKGROUND

Genomes are inherently inhomogeneous, with features such as base composition, recombination, gene density, and gene expression varying along chromosomes. Evolutionary, biological, and biomedical analyses aim to quantify this variation, account for it during inference procedures, and ultimately determine the causal processes behind it. Since sequential observations along chromosomes are not independent, it is unsurprising that autocorrelation patterns have been observed e.g., in human base composition. In this article, we develop a class of Hidden Markov Models (HMMs) called oHMMed (ordered HMM with emission densities, the corresponding R package of the same name is available on CRAN): They identify the number of comparably homogeneous regions within autocorrelated observed sequences. These are modelled as discrete hidden states; the observed data points are realisations of continuous probability distributions with state-specific means that enable ordering of these distributions. The observed sequence is labelled according to the hidden states, permitting only neighbouring states that are also neighbours within the ordering of their associated distributions. The parameters that characterise these state-specific distributions are inferred.

RESULTS

We apply our oHMMed algorithms to the proportion of G and C bases (modelled as a mixture of normal distributions) and the number of genes (modelled as a mixture of poisson-gamma distributions) in windows along the human, mouse, and fruit fly genomes. This results in a partitioning of the genomes into regions by statistically distinguishable averages of these features, and in a characterisation of their continuous patterns of variation. In regard to the genomic G and C proportion, this latter result distinguishes oHMMed from segmentation algorithms based in isochore or compositional domain theory. We further use oHMMed to conduct a detailed analysis of variation of chromatin accessibility (ATAC-seq) and epigenetic markers H3K27ac and H3K27me3 (modelled as a mixture of poisson-gamma distributions) along the human chromosome 1 and their correlations.

CONCLUSIONS

Our algorithms provide a biologically assumption free approach to characterising genomic landscapes shaped by continuous, autocorrelated patterns of variation. Despite this, the resulting genome segmentation enables extraction of compositionally distinct regions for further downstream analyses.

HYENA detects oncogenes activated by distal enhancers in cancer

Somatic structural variations (SVs) in cancer can shuffle DNA content in the genome, relocate regulatory elements, and alter genome organization. Enhancer hijacking occurs when SVs relocate distal enhancers to activate proto-oncogenes. However, most enhancer hijacking studies have only focused on protein-coding genes. Here, we develop a computational algorithm "HYENA" to identify candidate oncogenes (both protein-coding and non-coding) activated by enhancer hijacking based on tumor whole-genome and transcriptome sequencing data. HYENA detects genes whose elevated expression is associated with somatic SVs by using a rank-based regression model. We systematically analyze 1,146 tumors across 25 types of adult tumors and identify a total of 108 candidate oncogenes including many non-coding genes. A long non-coding RNA TOB1-AS1 is activated by various types of SVs in 10% of pancreatic cancers through altered 3-dimensional genome structure. We find that high expression of TOB1-AS1 can promote cell invasion and metastasis. Our study highlights the contribution of genetic alterations in non-coding regions to tumorigenesis and tumor progression.

CUT&Tag Applied to Zebrafish Adult Tail Fins Reveals a Return of Embryonic H3K4me3 Patterns During Regeneration

Regenerative potential is governed by a complex process of transcriptional reprogramming, involving chromatin reorganization and dynamics in transcription factor binding patterns throughout the genome. The degree to which chromatin and epigenetic changes contribute to this process remains partially understood. Here we provide a modified CUT&Tag protocol suitable for improved characterization and interrogation of epigenetic changes during adult fin regeneration in zebrafish. Our protocol generates data that recapitulates results from previously published ChIP-Seq methods, requires far fewer cells as input, and significantly improves signal to noise ratios. We deliver high-resolution enrichment maps for H3K4me3 of uninjured and regenerating fin tissues. During regeneration, we find that H3K4me3 levels increase over gene promoters which become transcriptionally active and genes which lose H3K4me3 become silenced. Interestingly, these epigenetic reprogramming events recapitulate the H3K4me3 patterns observed in developing fin folds of 24-hour old zebrafish embryos. Our results indicate that changes in genomic H3K4me3 patterns during fin regeneration occur in a manner consistent with reactivation of developmental programs, demonstrating CUT&Tag to be an effective tool for profiling chromatin landscapes in regenerating tissues.

Strategies for non-viral vectors targeting organs beyond the liver

In recent years, nanoparticles have evolved to a clinical modality to deliver diverse nucleic acids. Rising interest in nanomedicines comes from proven safety and efficacy profiles established by continuous efforts to optimize physicochemical properties and endosomal escape. However, despite their transformative impact on the pharmaceutical industry, the clinical use of non-viral nucleic acid delivery is limited to hepatic diseases and vaccines due to liver accumulation. Overcoming liver tropism of nanoparticles is vital to meet clinical needs in other organs. Understanding the anatomical structure and physiological features of various organs would help to identify potential strategies for fine-tuning nanoparticle characteristics. In this Review, we discuss the source of liver tropism of non-viral vectors, present a brief overview of biological structure, processes and barriers in select organs, highlight approaches available to reach non-liver targets, and discuss techniques to accelerate the discovery of non-hepatic therapies.

Genetic variation across and within individuals

Germline variation and somatic mutation are intricately connected and together shape human traits and disease risks. Germline variants are present from conception, but they vary between individuals and accumulate over generations. By contrast, somatic mutations accumulate throughout life in a mosaic manner within an individual due to intrinsic and extrinsic sources of mutations and selection pressures acting on cells. Recent advancements, such as improved detection methods and increased resources for association studies, have drastically expanded our ability to investigate germline and somatic genetic variation and compare underlying mutational processes. A better understanding of the similarities and differences in the types, rates and patterns of germline and somatic variants, as well as their interplay, will help elucidate the mechanisms underlying their distinct yet interlinked roles in human health and biology.

Construction of single-cell cross-species chromatin accessibility landscapes with combinatorial-hybridization-based ATAC-seq

Despite recent advances in single-cell genomics, the lack of maps for single-cell candidate cis-regulatory elements (cCREs) in non-mammal species has limited our exploration of conserved regulatory programs across vertebrates and invertebrates. Here, we developed a combinatorial-hybridization-based method for single-cell assay for transposase-accessible chromatin using sequencing (scATAC-seq) named CH-ATAC-seq, enabling the construction of single-cell accessible chromatin landscapes for zebrafish, Drosophila, and earthworms (Eisenia andrei). By integrating scATAC censuses of humans, monkeys, and mice, we systematically identified 152 distinct main cell types and around 0.8 million cell-type-specific cCREs. Our analysis provided insights into the conservation of neural, muscle, and immune lineages across species, while epithelial cells exhibited a higher organ-origin heterogeneity. Additionally, a large-scale gene regulatory network (GRN) was constructed in four vertebrates by integrating scRNA-seq censuses. Overall, our study provides a valuable resource for comparative epigenomics, identifying the evolutionary conservation and divergence of gene regulation across different species.

Cryopreservation of cerebrospinal fluid cells preserves the transcriptional landscape for single-cell analysis

Cerebrospinal fluid (CSF) matrix biomarkers have become increasingly valuable surrogate markers of neuropsychiatric diseases in research and clinical practice. In contrast, CSF cells have been rarely investigated due to their relative scarcity and fragility, and lack of common collection and cryopreservation protocols, with limited exceptions for neurooncology and primary immune-based diseases like multiple sclerosis. the advent of a microfluidics-based multi-omics approach to studying individual cells has allowed for the study of cellular phenotyping, intracellular dynamics, and intercellular relationships that provide multidimensionality unable to be obtained through acellular fluid-phase analyses. challenges to cell-based research include site-to-site differences in handling, storage, and thawing methods, which can lead to inaccuracy and inter-assay variability. In the present study, we performed single-cell RNA sequencing (10x Genomics) on fresh or previously cryopreserved human CSF samples from three alternative cryopreservation methods: Fetal Bovine Serum with Dimethyl sulfoxide (FBS/DMSO), FBS/DMSO after a DNase step (a step often included in epigenetic studies), and cryopreservation using commercially available Recovery© media. In comparing relative differences between fresh and cryopreserved samples, we found little effect of the cryopreservation method on being able to resolve donor-linked cell type proportions, markers of cellular stress, and overall gene expression at the single-cell level, whereas donor-specific differences were readily discernable. We further demonstrate the compatibility of fresh and cryopreserved CSF immune cell sequencing using biologically relevant sexually dimorphic gene expression differences by donor. Our findings support the utility and interchangeability of FBS/DMSO and Recovery cryopreservation with fresh sample analysis, providing a methodological grounding that will enable researchers to further expand our understanding of the CSF immune cell contributions to neurological and psychiatric disease.

Hi-C, a chromatin 3D structure technique advancing the functional genomics of immune cells

The functional performance of immune cells relies on a complex transcriptional regulatory network. The three-dimensional structure of chromatin can affect chromatin status and gene expression patterns, and plays an important regulatory role in gene transcription. Currently available techniques for studying chromatin spatial structure include chromatin conformation capture techniques and their derivatives, chromatin accessibility sequencing techniques, and others. Additionally, the recently emerged deep learning technology can be utilized as a tool to enhance the analysis of data. In this review, we elucidate the definition and significance of the three-dimensional chromatin structure, summarize the technologies available for studying it, and describe the research progress on the chromatin spatial structure of dendritic cells, macrophages, T cells, B cells, and neutrophils.

PRMT1 promotes pancreatic cancer development and resistance to chemotherapy

Pancreatic ductal adenocarcinoma (PDAC) remains one of the most lethal types of cancer, and novel treatment regimens are direly needed. Epigenetic regulation contributes to the development of various cancer types, but its role in the development of and potential as a therapeutic target for PDAC remains underexplored. Here, we show that PRMT1 is highly expressed in murine and human pancreatic cancer and is essential for cancer cell proliferation and tumorigenesis. Deletion of PRMT1 delays pancreatic cancer development in a KRAS-dependent mouse model, and multi-omics analyses reveal that PRMT1 depletion leads to global changes in chromatin accessibility and transcription, resulting in reduced glycolysis and a decrease in tumorigenic capacity. Pharmacological inhibition of PRMT1 in combination with gemcitabine has a synergistic effect on pancreatic tumor growth in vitro and in vivo. Collectively, our findings implicate PRMT1 as a key regulator of pancreatic cancer development and a promising target for combination therapy.

Arresting the bad seed: HDAC3 regulates proliferation of different microglia after ischemic stroke

The accumulation of self-renewed polarized microglia in the penumbra is a critical neuroinflammatory process after ischemic stroke, leading to secondary demyelination and neuronal loss. Although known to regulate tumor cell proliferation and neuroinflammation, HDAC3's role in microgliosis and microglial polarization remains unclear. We demonstrated that microglial HDAC3 knockout (HDAC3-miKO) ameliorated poststroke long-term functional and histological outcomes. RNA-seq analysis revealed mitosis as the primary process affected in HDAC3-deficent microglia following stroke. Notably, HDAC3-miKO specifically inhibited proliferation of proinflammatory microglia without affecting anti-inflammatory microglia, preventing microglial transition to a proinflammatory state. Moreover, ATAC-seq showed that HDAC3-miKO induced closing of accessible regions enriched with PU.1 motifs. Overexpressing microglial PU.1 via an AAV approach reversed HDAC3-miKO-induced proliferation inhibition and protective effects on ischemic stroke, indicating PU.1 as a downstream molecule that mediates HDAC3's effects on stroke. These findings uncovered that HDAC3/PU.1 axis, which mediated differential proliferation-related reprogramming in different microglia populations, drove poststroke inflammatory state transition, and contributed to pathophysiology of ischemic stroke.

Application of Single Cell Gene Expression Technologies to Neurotoxicology

Neurotoxicological research faces the challenge of linking biological changes resulting from exposures to neuronal function. An additional challenge is understanding cell-type specific differences and selective vulnerabilities of distinct neuronal populations to toxic insults. Single cell RNA-sequencing (scRNA-seq) allows for measurement of the transcriptome of individual cells. This makes it a valuable tool for validating and characterizing cell types present in multicell type samples in complex tissue or cell culture models, but also for understanding how different cell types respond to toxic insults. Pathway analysis of differentially expressed genes can provide in depth insights into underlying cell type-specific mechanisms of neurotoxicity. Toxicological data often has to be translated to outcomes for human health which requires an understanding of inter-species differences. Transcriptomic data aids in understanding these differences, including understanding developmental timelines of different species. We believe that scRNA-seq holds exciting promises for future neurotoxicological research.

Nano-CUT&Tag for multimodal chromatin profiling at single-cell resolution

The ability to comprehensively analyze the chromatin state with single-cell resolution is crucial for understanding gene regulatory principles in heterogenous tissues or during development. Recently, we developed a nanobody-based single-cell CUT&Tag (nano-CT) protocol to simultaneously profile three epigenetic modalities-two histone marks and open chromatin state-from the same single cell. Nano-CT implements a new set of secondary nanobody-Tn5 fusion proteins to direct barcoded tagmentation by Tn5 transposase to genomic targets labeled by primary antibodies raised in different species. Such nanobody-Tn5 fusion proteins are currently not commercially available, and their in-house production and purification can be completed in 3-4 d by following our detailed protocol. The single-cell indexing in nano-CT is performed on a commercially available platform, making it widely accessible to the community. In comparison to other multimodal methods, nano-CT stands out in data complexity, low sample requirements and the flexibility to choose two of the three modalities. In addition, nano-CT works efficiently with fresh brain samples, generating multimodal epigenomic profiles for thousands of brain cells at single-cell resolution. The nano-CT protocol can be completed in just 3 d by users with basic skills in standard molecular biology and bioinformatics, although previous experience with single-cell assay for transposase-accessible chromatin using sequencing (scATAC-seq) is beneficial for more in-depth data analysis. As a multimodal assay, nano-CT holds immense potential to reveal interactions of various chromatin modalities, to explore epigenetic heterogeneity and to increase our understanding of the role and interplay that chromatin dynamics has in cellular development.

Moving toward precision medicine to predict drug sensitivity in patients with metastatic breast cancer

Tumor heterogeneity represents a major challenge in breast cancer, being associated with disease progression and treatment resistance. Precision medicine has been extensively applied to dissect tumor heterogeneity and, through a deeper molecular understanding of the disease, to personalize therapeutic strategies. In the last years, technological advances have widely improved the understanding of breast cancer biology and several trials have been developed to translate these new insights into clinical practice, with the ultimate aim of improving patients' outcomes. In the era of molecular oncology, genomics analyses and other methodologies are shaping a new treatment algorithm in breast cancer care. In this manuscript, we review the main steps of precision medicine to predict drug sensitivity in breast cancer from a translational point of view. Genomic developments and their clinical implications are discussed, along with technological advancements that could broaden precision medicine applications. Current achievements are put into perspective to provide an overview of the state-of-art of breast cancer precision oncology as well as to identify future research directions.

Integrated ATAC-seq and RNA-seq data analysis identifies transcription factors related to rice stripe virus infection in Oryza sativa

Animal studies have shown that virus infection causes changes in host chromatin accessibility, but little is known about changes in chromatin accessibility of plants infected by viruses and its potential impact. Here, rice infected by rice stripe virus (RSV) was used to investigate virus-induced changes in chromatin accessibility. Our analysis identified a total of 6462 open- and 3587 closed-differentially accessible chromatin regions (DACRs) in rice under RSV infection by ATAC-seq. Additionally, by integrating ATAC-seq and RNA-seq, 349 up-regulated genes in open-DACRs and 126 down-regulated genes in closed-DACRs were identified, of which 34 transcription factors (TFs) were further identified by search of upstream motifs. Transcription levels of eight of these TFs were validated by reverse transcription-PCR. Importantly, four of these TFs (OsWRKY77, OsWRKY28, OsZFP12 and OsERF91) interacted with RSV proteins and are therefore predicted to play important roles in RSV infection. This is the first application of ATAC-seq and RNA-seq techniques to analyse changes in rice chromatin accessibility caused by RSV infection. Integrating ATAC-seq and RNA-seq provides a new approach to select candidate TFs in response to virus infection.

Improved metagenome assemblies through selective enrichment of bacterial genomic DNA from eukaryotic host genomic DNA using ATAC-seq

Genomics can be used to study the complex relationships between hosts and their microbiota. Many bacteria cannot be cultured in the laboratory, making it difficult to obtain adequate amounts of bacterial DNA and to limit host DNA contamination for the construction of metagenome-assembled genomes (MAGs). For example, Wolbachia is a genus of exclusively obligate intracellular bacteria that live in a wide range of arthropods and some nematodes. While Wolbachia endosymbionts are frequently described as facultative reproductive parasites in arthropods, the bacteria are obligate mutualistic endosymbionts of filarial worms. Here, we achieve 50-fold enrichment of bacterial sequences using ATAC-seq (Assay for Transposase-Accessible Chromatin using sequencing) with Brugia malayi nematodes, containing Wolbachia (wBm). ATAC-seq uses the Tn5 transposase to cut and attach Illumina sequencing adapters to accessible DNA lacking histones, typically thought to be open chromatin. Bacterial and mitochondrial DNA in the lysates are also cut preferentially since they lack histones, leading to the enrichment of these sequences. The benefits of this include minimal tissue input (<1 mg of tissue), a quick protocol (<4 h), low sequencing costs, less bias, correct assembly of lateral gene transfers and no prior sequence knowledge required. We assembled the wBm genome with as few as 1 million Illumina short paired-end reads with >97% coverage of the published genome, compared to only 12% coverage with the standard gDNA libraries. We found significant bacterial sequence enrichment that facilitated genome assembly in previously published ATAC-seq data sets from human cells infected with Mycobacterium tuberculosis and C. elegans contaminated with their food source, the OP50 strain of E. coli. These results demonstrate the feasibility and benefits of using ATAC-seq to easily obtain bacterial genomes to aid in symbiosis, infectious disease, and microbiome research.

Quantitative transcriptomic and epigenomic data analysis: a primer

The advent of microarray and second generation sequencing technology has revolutionized the field of molecular biology, allowing researchers to quantitatively assess transcriptomic and epigenomic features in a comprehensive and cost-efficient manner. Moreover, technical advancements have pushed the resolution of these sequencing techniques to the single cell level. As a result, the bottleneck of molecular biology research has shifted from the bench to the subsequent omics data analysis. Even though most methodologies share the same general strategy, state-of-the-art literature typically focuses on data type specific approaches and already assumes expert knowledge. Here, however, we aim at providing conceptual insight in the principles of genome-wide quantitative transcriptomic and epigenomic (including open chromatin assay) data analysis by describing a generic workflow. By starting from a general framework and its assumptions, the need for alternative or additional data-analytical solutions when working with specific data types becomes clear, and are hence introduced. Thus, we aim to enable readers with basic omics expertise to deepen their conceptual and statistical understanding of general strategies and pitfalls in omics data analysis and to facilitate subsequent progression to more specialized literature.

Predicting the unpredictable: the regulatory nature and promiscuity of herbicide cross resistance

The emergence of herbicide-resistant weeds is a significant threat to modern agriculture. Cross resistance, a phenomenon where resistance to one herbicide confers resistance to another, is a particular concern owing to its unpredictability. Nontarget-site (NTS) cross resistance is especially challenging to predict, as it arises from genes that encode enzymes that do not directly involve the herbicide target site and can affect multiple herbicides. Recent advancements in genomic and structural biology techniques could provide new venues for predicting NTS resistance in weed species. In this review, we present an overview of the latest approaches that could be used. We discuss the use of genomic and epigenomics techniques such as ATAC-seq and DAP-seq to identify transcription factors and cis-regulatory elements associated with resistance traits. Enzyme/protein structure prediction and docking analysis are discussed as an initial step for predicting herbicide binding affinities with key enzymes to identify candidates for subsequent in vitro validation. We also provide example analyses that can be deployed toward elucidating cross resistance and its regulatory patterns. Ultimately, our review provides important insights into the latest scientific advancements and potential directions for predicting and managing herbicide cross resistance in weeds. © 2023 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

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.

Multi-omics integrative analysis revealed characteristic changes in blood cell immunity and amino acid metabolism in a silkworm model of hyperproteinemia

Hyperproteinemia is a serious metabolic disease of both humans and animals characterized by an abnormally high plasma protein concentration (HPPC). Although hyperproteinemia can cause an imbalance in blood cell homeostasis, the functional changes to blood cells remain unclear. Here, a HPPC silkworm model was used to assess changes to the chromatin accessibility and transcript levels of genes related to blood cell metabolism and immune function. The results showed that HPPC enhanced phagocytosis of blood cells, increased chromatin accessibility and transcript levels of genes involved in cell phagocytosis, proliferation, stress, and programmed death, while genes associated with aromatic amino acid metabolism, and antibacterial peptide synthesis were inhibited in blood cells. Further analysis of the chromatin accessibility of the promoter region found that the high chromatin accessibility of genes sensitive to HPPC, was related to histone modifications, including tri-methylation of lysine residue 4 of histone H3 and acetylation of lysine residue 27 of histone H3. Changes to the chromatin accessibility and transcript levels of genes related to immune function and amino acid metabolism in the blood cells of the HPPC silkworm model provided useful references for future studies of the mechanisms underlying epigenomic regulation mediated by hyperproteinemia.

TRPS1 modulates chromatin accessibility to regulate estrogen receptor alpha (ER) binding and ER target gene expression in luminal breast cancer cells

Common genetic variants in the repressive GATA-family transcription factor (TF) TRPS1 locus are associated with breast cancer risk, and luminal breast cancer cell lines are particularly sensitive to TRPS1 knockout. We introduced an inducible degron tag into the native TRPS1 locus within a luminal breast cancer cell line to identify the direct targets of TRPS1 and determine how TRPS1 mechanistically regulates gene expression. We acutely deplete over 80 percent of TRPS1 from chromatin within 30 minutes of inducing degradation. We find that TRPS1 regulates transcription of hundreds of genes, including those related to estrogen signaling. TRPS1 directly regulates chromatin structure, which causes estrogen receptor alpha (ER) to redistribute in the genome. ER redistribution leads to both repression and activation of dozens of ER target genes. Downstream from these primary effects, TRPS1 depletion represses cell cycle-related gene sets and reduces cell doubling rate. Finally, we show that high TRPS1 activity, calculated using a gene expression signature defined by primary TRPS1-regulated genes, is associated with worse breast cancer patient prognosis. Taken together, these data suggest a model in which TRPS1 modulates the genomic distribution of ER, both activating and repressing transcription of genes related to cancer cell fitness.

Application of Single-Cell Assay for Transposase-Accessible Chromatin with High Throughput Sequencing in Plant Science: Advances, Technical Challenges, and Prospects

The Single-cell Assay for Transposase-Accessible Chromatin with high throughput sequencing (scATAC-seq) has gained increasing popularity in recent years, allowing for chromatin accessibility to be deciphered and gene regulatory networks (GRNs) to be inferred at single-cell resolution. This cutting-edge technology now enables the genome-wide profiling of chromatin accessibility at the cellular level and the capturing of cell-type-specific cis-regulatory elements (CREs) that are masked by cellular heterogeneity in bulk assays. Additionally, it can also facilitate the identification of rare and new cell types based on differences in chromatin accessibility and the charting of cellular developmental trajectories within lineage-related cell clusters. Due to technical challenges and limitations, the data generated from scATAC-seq exhibit unique features, often characterized by high sparsity and noise, even within the same cell type. To address these challenges, various bioinformatic tools have been developed. Furthermore, the application of scATAC-seq in plant science is still in its infancy, with most research focusing on root tissues and model plant species. In this review, we provide an overview of recent progress in scATAC-seq and its application across various fields. We first conduct scATAC-seq in plant science. Next, we highlight the current challenges of scATAC-seq in plant science and major strategies for cell type annotation. Finally, we outline several future directions to exploit scATAC-seq technologies to address critical challenges in plant science, ranging from plant ENCODE(The Encyclopedia of DNA Elements) project construction to GRN inference, to deepen our understanding of the roles of CREs in plant biology.

Gestational diabetes in mice induces hematopoietic memory that affects the long-term health of the offspring

Gestational diabetes is a common medical complication of pregnancy that is associated with adverse perinatal outcomes and an increased risk of metabolic diseases and atherosclerosis in adult offspring. The mechanisms responsible for this delayed pathological transmission remain unknown. In mouse models, we found that the development of atherosclerosis in adult offspring born to diabetic pregnancy can be in part linked to hematopoietic alterations. Although they do not show any gross metabolic disruptions, the adult offspring maintain hematopoietic features associated with diabetes, indicating the acquisition of a lasting diabetic hematopoietic memory. We show that the induction of this hematopoietic memory during gestation relies on the activity of the advanced glycation end product receptor (AGER) and the nucleotide binding and oligomerization domain-like receptor family pyrin domain-containing 3 (NLRP3) inflammasome, which lead to increased placental inflammation. In adult offspring, we find that this memory is associated with DNA methyltransferase 1 (DNMT1) upregulation and epigenetic changes in hematopoietic progenitors. Together, our results demonstrate that the hematopoietic system can acquire a lasting memory of gestational diabetes and that this memory constitutes a pathway connecting gestational health to adult pathologies.

MethMotif.Org 2024: a database integrating context-specific transcription factor-binding motifs with DNA methylation patterns

MethMotif (https://methmotif.org) is a publicly available database that provides a comprehensive repository of transcription factor (TF)-binding profiles, enriched with DNA methylation patterns. In this release, we have enhanced the platform, expanding our initial collection to over 700 position weight matrices (PWM), all of which include DNA methylation profiles. One of the key advancements in this release is the segregation of TF-binding motifs based on their cofactors and DNA methylation status. We have previously demonstrated that gene ontology (GO) enriched terms associated with TF target genes may differ based on their association with alternative cofactors and DNA methylation status. MethMotif provides precomputed GO annotations for each human TF of interest, as well as for TF-co-TF complexes, enabling a comprehensive analysis of TF functions in the context of their co-factors. Additionally, MethMotif has been updated to encompass data for two new species, Mus musculus and Arabidopsis thaliana, widening its applicability to a broader community. MethMotif stands out as the first and only TF-binding motifs database to incorporate context-specific PWM coupled with epigenetic information, thereby enlightening context-specific TF functions. This enhancement allows the community to explore and gain deeper insights into the regulatory mechanisms governing transcriptional processes.

MCT1-governed pyruvate metabolism is essential for antibody class-switch recombination through H3K27 acetylation

Monocarboxylate transporter 1 (MCT1) exhibits essential roles in cellular metabolism and energy supply. Although MCT1 is highly expressed in activated B cells, it is not clear how MCT1-governed monocarboxylates transportation is functionally coupled to antibody production during the glucose metabolism. Here, we report that B cell-lineage deficiency of MCT1 significantly influences the class-switch recombination (CSR), rendering impaired IgG antibody responses in Mct1f/fMb1Cre mice after immunization. Metabolic flux reveals that glucose metabolism is significantly reprogrammed from glycolysis to oxidative phosphorylation in Mct1-deficient B cells upon activation. Consistently, activation-induced cytidine deaminase (AID), is severely suppressed in Mct1-deficient B cells due to the decreased level of pyruvate metabolite. Mechanistically, MCT1 is required to maintain the optimal concentration of pyruvate to secure the sufficient acetylation of H3K27 for the elevated transcription of AID in activated B cells. Clinically, we found that MCT1 expression levels are significantly upregulated in systemic lupus erythematosus patients, and Mct1 deficiency can alleviate the symptoms of bm12-induced murine lupus model. Collectively, these results demonstrate that MCT1-mediated pyruvate metabolism is required for IgG antibody CSR through an epigenetic dependent AID transcription, revealing MCT1 as a potential target for vaccine development and SLE disease treatment.

Cis-Regulatory Elements in Mammals

In cis-regulatory elements, enhancers and promoters with complex molecular interactions are used to coordinate gene transcription through physical proximity and chemical modifications. These processes subsequently influence the phenotypic characteristics of an organism. An in-depth exploration of enhancers and promoters can substantially enhance our understanding of gene regulatory networks, shedding new light on mammalian development, evolution and disease pathways. In this review, we provide a comprehensive overview of the intrinsic structural attributes, detection methodologies as well as the operational mechanisms of enhancers and promoters, coupled with the relevant novel and innovative investigative techniques used to explore their actions. We further elucidated the state-of-the-art research on the roles of enhancers and promoters in the realms of mammalian development, evolution and disease, and we conclude with forward-looking insights into prospective research avenues.

Integration of ATAC-Seq and RNA-Seq Analysis to Identify Key Genes in the Longissimus Dorsi Muscle Development of the Tianzhu White Yak

During the postnatal stages, skeletal muscle development undergoes a series of meticulously regulated alterations in gene expression. However, limited studies have employed chromatin accessibility to unravel the underlying molecular mechanisms governing muscle development in yak species. Therefore, we conducted an analysis of both gene expression levels and chromatin accessibility to comprehensively characterize the dynamic genome-wide chromatin accessibility during muscle growth and development in the Tianzhu white yak, thereby elucidating the features of accessible chromatin regions throughout this process. Initially, we compared the differences in chromatin accessibility between two groups and observed that calves exhibited higher levels of chromatin accessibility compared to adult cattle, particularly within ±2 kb of the transcription start site (TSS). In order to investigate the correlation between alterations in chromatin accessible regions and variations in gene expression levels, we employed a combination of ATAC-seq and RNA-seq techniques, leading to the identification of 18 central transcriptional factors (TFs) and 110 key genes with significant effects. Through further analysis, we successfully identified several TFs, including Sp1, YY1, MyoG, MEF2A and MEF2C, as well as a number of candidate genes (ANKRD2, ANKRD1, BTG2 and LMOD3) which may be closely associated with muscle growth and development. Moreover, we constructed an interactive network program encompassing hub TFs and key genes related to muscle growth and development. This innovative approach provided valuable insights into the molecular mechanism underlying skeletal muscle development in the postnatal stages of Tianzhu white yaks while also establishing a solid theoretical foundation for future research on yak muscle development.

Dynamics of chromatin accessibility and genome wide control of desiccation tolerance in the resurrection plant Haberlea rhodopensis

BACKGROUND

Drought is one of the main consequences of global climate change and this problem is expected to intensify in the future. Resurrection plants evolved the ability to withstand the negative impact of long periods of almost complete desiccation and to recover at rewatering. In this respect, many physiological, transcriptomic, proteomic and genomic investigations have been performed in recent years, however, few epigenetic control studies have been performed on these valuable desiccation-tolerant plants so far.

RESULTS

In the present study, for the first time for resurrection plants we provide evidences about the differential chromatin accessibility of Haberlea rhodopensis during desiccation stress by ATAC-seq (Assay for Transposase Accessible Chromatin with high-throughput sequencing). Based on gene similarity between species, we used the available genome of the closely related resurrection plant Dorcoceras hygrometricum to identify approximately nine hundred transposase hypersensitive sites (THSs) in H. rhodopensis. The majority of them corresponds to proximal and distal regulatory elements of different genes involved in photosynthesis, carbon metabolism, synthesis of secondary metabolites, cell signalling and transcriptional regulation, cell growth, cell wall, stomata conditioning, chaperons, oxidative stress, autophagy and others. Various types of binding motifs recognized by several families of transcription factors have been enriched from the THSs found in different stages of drought. Further, we used the previously published RNA-seq data from H. rhodopensis to evaluate the expression of transcription factors putatively interacting with the enriched motifs, and the potential correlation between the identified THS and the expression of their corresponding genes.

CONCLUSIONS

These results provide a blueprint for investigating the epigenetic regulation of desiccation tolerance in resurrection plant H. rhodopensis and comparative genomics between resurrection and non-resurrection species with available genome information.

Transcriptional changes during isoproterenol-induced cardiac fibrosis in mice

Introduction: β-adrenergic stimulation using β-agonists such as isoproterenol has been routinely used to induce cardiac fibrosis in experimental animal models. Although transcriptome changes in surgical models of cardiac fibrosis such as transverse aortic constriction (TAC) and coronary artery ligation (CAL) are well-studied, transcriptional changes during isoproterenol-induced cardiac fibrosis are not well-explored. Methods: Cardiac fibrosis was induced in male C57BL6 mice by administration of isoproterenol for 4, 8, or 11 days at 50 mg/kg/day dose. Temporal changes in gene expression were studied by RNA sequencing. Results and discussion: We observed a significant alteration in the transcriptome profile across the different experimental groups compared to the saline group. Isoproterenol treatment caused upregulation of genes associated with ECM organization, cell-cell contact, three-dimensional structure, and cell growth, while genes associated with fatty acid oxidation, sarcoplasmic reticulum calcium ion transport, and cardiac muscle contraction are downregulated. A number of known long non-coding RNAs (lncRNAs) and putative novel lncRNAs exhibited differential regulation. In conclusion, our study shows that isoproterenol administration leads to the dysregulation of genes relevant to ECM deposition and cardiac contraction, and serves as an excellent alternate model to the surgical models of heart failure.

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.

Integrating omics databases for enhanced crop breeding

Crop plant breeding involves selecting and developing new plant varieties with desirable traits such as increased yield, improved disease resistance, and enhanced nutritional value. With the development of high-throughput technologies, such as genomics, transcriptomics, and metabolomics, crop breeding has entered a new era. However, to effectively use these technologies, integration of multi-omics data from different databases is required. Integration of omics data provides a comprehensive understanding of the biological processes underlying plant traits and their interactions. This review highlights the importance of integrating omics databases in crop plant breeding, discusses available omics data and databases, describes integration challenges, and highlights recent developments and potential benefits. Taken together, the integration of omics databases is a critical step towards enhancing crop plant breeding and improving global food security.

Improved quality metrics for association and reproducibility in chromatin accessibility data using mutual information

BACKGROUND

Correlation metrics are widely utilized in genomics analysis and often implemented with little regard to assumptions of normality, homoscedasticity, and independence of values. This is especially true when comparing values between replicated sequencing experiments that probe chromatin accessibility, such as assays for transposase-accessible chromatin via sequencing (ATAC-seq). Such data can possess several regions across the human genome with little to no sequencing depth and are thus non-normal with a large portion of zero values. Despite distributed use in the epigenomics field, few studies have evaluated and benchmarked how correlation and association statistics behave across ATAC-seq experiments with known differences or the effects of removing specific outliers from the data. Here, we developed a computational simulation of ATAC-seq data to elucidate the behavior of correlation statistics and to compare their accuracy under set conditions of reproducibility.

RESULTS

Using these simulations, we monitored the behavior of several correlation statistics, including the Pearson's R and Spearman's [Formula: see text] coefficients as well as Kendall's [Formula: see text] and Top-Down correlation. We also test the behavior of association measures, including the coefficient of determination R[Formula: see text], Kendall's W, and normalized mutual information. Our experiments reveal an insensitivity of most statistics, including Spearman's [Formula: see text], Kendall's [Formula: see text], and Kendall's W, to increasing differences between simulated ATAC-seq replicates. The removal of co-zeros (regions lacking mapped sequenced reads) between simulated experiments greatly improves the estimates of correlation and association. After removing co-zeros, the R[Formula: see text] coefficient and normalized mutual information display the best performance, having a closer one-to-one relationship with the known portion of shared, enhanced loci between simulated replicates. When comparing values between experimental ATAC-seq data using a random forest model, mutual information best predicts ATAC-seq replicate relationships.

CONCLUSIONS

Collectively, this study demonstrates how measures of correlation and association can behave in epigenomics experiments. We provide improved strategies for quantifying relationships in these increasingly prevalent and important chromatin accessibility assays.

Global identification of SWI/SNF targets reveals compensation by EP400

Mammalian SWI/SNF chromatin remodeling complexes move and evict nucleosomes at gene promoters and enhancers to modulate DNA access. Although SWI/SNF subunits are commonly mutated in disease, therapeutic options are limited by our inability to predict SWI/SNF gene targets and conflicting studies on functional significance. Here, we leverage a fast-acting inhibitor of SWI/SNF remodeling to elucidate direct targets and effects of SWI/SNF. Blocking SWI/SNF activity causes a rapid and global loss of chromatin accessibility and transcription. Whereas repression persists at most enhancers, we uncover a compensatory role for the EP400/TIP60 remodeler, which reestablishes accessibility at most promoters during prolonged loss of SWI/SNF. Indeed, we observe synthetic lethality between EP400 and SWI/SNF in cancer cell lines and human cancer patient data. Our data define a set of molecular genomic features that accurately predict gene sensitivity to SWI/SNF inhibition in diverse cancer cell lines, thereby improving the therapeutic potential of SWI/SNF inhibitors.

Vaginal microbes alter epithelial transcriptomic and epigenomic modifications providing insight into the molecular mechanisms for susceptibility to adverse reproductive outcomes

The cervicovaginal microbiome is highly associated with women's health with microbial communities dominated by Lactobacillus spp. being considered optimal. Conversely, a lack of lactobacilli and a high abundance of strict and facultative anaerobes including Gardnerella vaginalis , have been associated with adverse reproductive outcomes. However, the molecular pathways modulated by microbe interactions with the cervicovaginal epithelia remain unclear. Using RNA-sequencing, we characterize the in vitro cervicovaginal epithelial transcriptional response to different vaginal bacteria and their culture supernatants. We showed that G. vaginalis upregulated genes were associated with an activated innate immune response including anti-microbial peptides and inflammasome pathways, represented by NLRP3-mediated increases in caspase-1, IL-1β and cell death. Cervicovaginal epithelial cells exposed to L. crispatus showed limited transcriptomic changes, while exposure to L. crispatus culture supernatants resulted in a shift in the epigenomic landscape of cervical epithelial cells. ATAC-sequencing confirmed epigenetic changes with reduced chromatin accessibility. This study reveals new insight into host-microbe interactions in the lower reproductive tract and suggest potential therapeutic strategies leveraging the vaginal microbiome to improve reproductive health.

Next-Generation sequencing transforming clinical practice and precision medicine

Next-generation sequencing (NGS) has revolutionized the field of genomics and is rapidly transforming clinical diagnosis and precision medicine. This advanced sequencing technology enables the rapid and cost-effective analysis of large-scale genomic data, allowing comprehensive exploration of the genetic landscape of diseases. In clinical diagnosis, NGS has proven to be a powerful tool for identifying disease-causing variants, enabling accurate and early detection of genetic disorders. Additionally, NGS facilitates the identification of novel disease-associated genes and variants, aiding in the development of targeted therapies and personalized treatment strategies. NGS greatly benefits precision medicine by enhancing our understanding of disease mechanisms and enabling the identification of specific molecular markers for disease subtypes, thus enabling tailored medical interventions based on individual characteristics. Furthermore, NGS contributes to the development of non-invasive diagnostic approaches, such as liquid biopsies, which can monitor disease progression and treatment response. The potential of NGS in clinical diagnosis and precision medicine is vast, yet challenges persist in data analysis, interpretation, and protocol standardization. This review highlights NGS applications in disease diagnosis, prognosis, and personalized treatment strategies, while also addressing challenges and future prospects in fully harnessing genomic potential within clinical practice.

Recent methodological advances towards single-cell proteomics

Studying the central dogma at the single-cell level has gained increasing attention to reveal hidden cell lineages and functions that cannot be studied using traditional bulk analyses. Nonetheless, most single-cell studies exploiting genomic and transcriptomic levels fail to address information on proteins that are central to many important biological processes. Single-cell proteomics enables understanding of the functional status of individual cells and is particularly crucial when the specimen is composed of heterogeneous entities of cells. With the growing importance of this field, significant methodological advancements have emerged recently. These include miniaturized and automated sample preparation, multi-omics analyses, and combined analyses of multiple techniques such as mass spectrometry and microscopy. Moreover, artificial intelligence and single-molecule detection technologies have advanced throughput and improved sensitivity limitations, respectively, over conventional methods. In this review, we summarize cutting-edge methodologies for single-cell proteomics and relevant emerging technologies that have been reported in the last 5 years, and provide an outlook on this research field.

Regulation of CTCF loop formation during pancreatic cell differentiation

Transcription reprogramming during cell differentiation involves targeting enhancers to genes responsible for establishment of cell fates. To understand the contribution of CTCF-mediated chromatin organization to cell lineage commitment, we analyzed 3D chromatin architecture during the differentiation of human embryonic stem cells into pancreatic islet organoids. We find that CTCF loops are formed and disassembled at different stages of the differentiation process by either recruitment of CTCF to new anchor sites or use of pre-existing sites not previously involved in loop formation. Recruitment of CTCF to new sites in the genome involves demethylation of H3K9me3 to H3K9me2, demethylation of DNA, recruitment of pioneer factors, and positioning of nucleosomes flanking the new CTCF sites. Existing CTCF sites not involved in loop formation become functional loop anchors via the establishment of new cohesin loading sites containing NIPBL and YY1 at sites between the new anchors. In both cases, formation of new CTCF loops leads to strengthening of enhancer promoter interactions and increased transcription of genes adjacent to loop anchors. These results suggest an important role for CTCF and cohesin in controlling gene expression during cell differentiation.

Chromatin Accessibility Landscape of Human Triple-negative Breast Cancer Cell Lines Reveals Variation by Patient Donor Ancestry

African American (AA) women have an excessive risk of developing triple-negative breast cancer (TNBC). We employed Assay for Transposase-Accessible Chromatin using sequencing to characterize differences in chromatin accessibility between nine commonly used TNBC cell lines derived from patients of European and African ancestry. Principal component and chromosome mapping analyses of accessibility peaks with the most variance revealed separation of chromatin profiles by patient group. Motif enrichment and footprinting analyses of disparate open chromatin regions revealed differences in transcription factor activity, identifying 79 with ancestry-associated binding patterns (FDR < 0.01). AA TNBC cell lines exhibited increased accessibility for 62 transcription factors associated with epithelial-to-mesenchymal transition, cancer stemness/chemotherapeutic resistance, proliferation, and aberrant p53 regulation, as well as KAISO, which has been previously linked to aggressive tumor characteristics in AA patients with cancer. Differential Assay for Transposase-Accessible Chromatin signal analysis identified 1,596 genes located within promoters of differentially open chromatin regions in AA-derived TNBC, identifying DNA methyltransferase 1 as the top upregulated gene associated with African ancestry. Pathway analyses with these genes revealed enrichment in several pathways, including hypoxia. Culturing cells under hypoxia showed ancestry-specific stress responses that led to the identification of a core set of AA-associated transcription factors, which included members of the Kruppel-like factor and Sp subfamilies, as well as KAISO, and identified ZDHHC1, a gene previously implicated in immunity and STING activation, as the top upregulated AA-specific gene under hypoxia. Together, these data reveal a differential chromatin landscape in TNBC associated with donor ancestry. The open chromatin structure of AA TNBC may contribute to a more lethal disease.

SIGNIFICANCE

We identify an ancestry-associated open chromatin landscape and related transcription factors that may contribute to aggressive TNBC in AA women. Furthermore, this study advocates for the inclusion of diversely sourced cell lines in experimental in vitro studies to advance health equity at all levels of scientific research.

Location, location, location: mapping the lymphoma tumor microenvironment using spatial transcriptomics

Lymphomas are a heterogenous group of lymphoid neoplasms with a wide variety of clinical presentations. Response to treatment and prognosis differs both between and within lymphoma subtypes. Improved molecular and genetic profiling has increased our understanding of the factors which drive these clinical dynamics. Immune and non-immune cells within the lymphoma tumor microenvironment (TME) can both play a key role in antitumor immune responses and conversely also support lymphoma growth and survival. A deeper understanding of the lymphoma TME would identify key lymphoma and immune cell interactions which could be disrupted for therapeutic benefit. Single cell RNA sequencing studies have provided a more comprehensive description of the TME, however these studies are limited in that they lack spatial context. Spatial transcriptomics provides a comprehensive analysis of gene expression within tissue and is an attractive technique in lymphoma to both disentangle the complex interactions between lymphoma and TME cells and improve understanding of how lymphoma cells evade the host immune response. This article summarizes current spatial transcriptomic technologies and their use in lymphoma research to date. The resulting data has already enriched our knowledge of the mechanisms and clinical impact of an immunosuppressive TME in lymphoma and the accrual of further studies will provide a fundamental step in the march towards personalized medicine.

Epigenetic programming in the ovarian reserve

The ovarian reserve defines female reproductive lifespan, which in humans spans decades. The ovarian reserve consists of oocytes residing in primordial follicles arrested in meiotic prophase I and is maintained independent of DNA replication and cell proliferation, thereby lacking stem cell-based maintenance. Largely unknown is how cellular states of the ovarian reserve are established and maintained for decades. Our recent study revealed that a distinct chromatin state is established during ovarian reserve formation in mice, uncovering a novel window of epigenetic programming in female germline development. We showed that an epigenetic regulator, Polycomb Repressive Complex 1 (PRC1), establishes a repressive chromatin state in perinatal mouse oocytes that is essential for prophase I-arrested oocytes to form the ovarian reserve. Here we discuss the biological roles and mechanisms underlying epigenetic programming in ovarian reserve formation, highlighting current knowledge gaps and emerging research areas in female reproductive biology.

Nucleosome repositioning in chronic lymphocytic leukemia

The location of nucleosomes in the human genome determines the primary chromatin structure and regulates access to regulatory regions. However, genome-wide information on deregulated nucleosome occupancy and its implications in primary cancer cells is scarce. Here, we conducted a genome-wide comparison of high-resolution nucleosome maps in peripheral blood B cells from patients with chronic lymphocytic leukemia (CLL) and healthy individuals at single-base-pair resolution. Our investigation uncovered significant changes of nucleosome positioning in CLL. Globally, the spacing between nucleosomes-the nucleosome repeat length (NRL)-is shortened in CLL. This effect is stronger in the more aggressive IGHV-unmutated CLL subtype than in the IGHV-mutated CLL subtype. Changes in nucleosome occupancy at specific sites are linked to active chromatin remodeling and reduced DNA methylation. Nucleosomes lost or gained in CLL marks differential binding of 3D chromatin organizers such as CTCF as well as immune response-related transcription factors and delineated mechanisms of epigenetic deregulation. The principal component analysis of nucleosome occupancy in cancer-specific regions allowed the classification of samples between cancer subtypes and normal controls. Furthermore, patients could be better assigned to CLL subtypes according to differential nucleosome occupancy than based on DNA methylation or gene expression. Thus, nucleosome positioning constitutes a novel readout to dissect molecular mechanisms of disease progression and to stratify patients. Furthermore, we anticipate that the global nucleosome repositioning detected in our study, such as changes in the NRL, can be exploited for liquid biopsy applications based on cell-free DNA to stratify patients and monitor disease progression.