SEARCH · PubMed Health

Results for “Programming Languages”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 recordsLinked to original sources

PubMed2025

Vcfexpress: flexible, rapid user-expressions to filter and format VCFs.

MOTIVATION: Variant call format (VCF) files are the standard output format for various software tools that identify genetic variation from DNA sequencing experiments. Downstream analyses require the ability to query, filter, and modify them simply and efficiently. Several tools are available to perform these operations from the command line, including BCFTools, vembrane, slivar, and others. RESULTS: Here, we introduce vcfexpress, a new, high-performance toolset for the analysis of VCF files, written in the Rust programming language. It is nearly as fast as BCFTools, but adds functionality to execute user expressions in the lua programming language for precise filtering and reporting of variants from a VCF or BCF file. We demonstrate performance and flexibility by comparing vcfexpress to other tools using the vembrane benchmark. AVAILABILITY AND IMPLEMENTATION: vcfexpress is available under the MIT license at https://github.com/brentp/vcfexpress with code used for the manuscript deposited in https://doi.org/10.5281/zenodo.14756838.

Software↗

PubMed2025

polars-bio-fast, scalable, and out-of-core operations on large genomic interval datasets.

MOTIVATION: Genomic studies very often rely on computationally intensive analyses of relationships between features, which are typically represented as intervals along a 1D coordinate system (such as positions on a chromosome). In this context, the Python programming language is extensively used for manipulating and analyzing data stored in a tabular form of rows and columns, called a DataFrame. Pandas is the most widely used Python DataFrame package and has been criticized for inefficiencies and scalability issues, which its modern alternative-Polars-aims to address with a native backend written in the Rust programming language. RESULTS: polars-bio is a Python library that enables fast, parallel and out-of-core operations on large genomic interval datasets. Its main components are implemented in Rust, using the Apache DataFusion query engine and Apache Arrow for efficient data representation. It is compatible with Polars and Pandas DataFrame formats. In a real-world comparison (107 versus 1.2×106 intervals), our library runs overlap queries 6.5×, nearest queries 15.5×, count_overlaps queries 38×, and coverage queries 15× faster than Bioframe. On equally sized synthetic sets (107 versus 107), the corresponding speedups are 1.6×, 5.5×, 6×, and 6×. In streaming mode, on real and synthetic interval pairs, our implementation uses 90× and 15× less memory for overlap, 4.5× and 6.5× less for nearest, 60× and 12× less for count_overlaps, and 34× and 7× less for coverage than Bioframe. Multi-threaded benchmarks show good scalability characteristics. To the best of our knowledge, polars-bio is the most efficient single-node library for genomic interval DataFrames in Python. AVAILABILITY AND IMPLEMENTATION: polars-bio is an open-source Python package distributed under the Apache License available for major platforms, including Linux, macOS, and Windows in the PyPI registry. The online documentation is https://biodatageeks.org/polars-bio/ and the source code is available on GitHub: https://github.com/biodatageeks/polars-bio and Zenodo: https://doi.org/10.5281/zenodo.16374290. are available at Bioinformatics online.

Software↗

PubMed2025

monarchr: an R package for querying biomedical knowledge graphs.

SUMMARY: Biomedical knowledge graphs (KGs) aggregate and provide a wealth of information, linking genes and their variants, diseases, phenotypes, and much more. While these data are available in raw and API-hosted form, to date, functionality for working with KGs in the R programming language has been limited. We introduce monarchr, a package for querying and manipulating KG data. Support for the expansive Monarch Initiative KG is built in, and monarchr can accommodate any KG in the Knowledge Graph eXchange (KGX) format. This tidy-inspired interface offers researchers an intuitive, iterative approach to querying and visualizing KG data. AVAILABILITY AND IMPLEMENTATION: Source code, documentation, and installation instructions are available at https://github.com/monarch-initiative/monarchr.

Software↗

PubMed2025

ShortCake: an integrated platform for efficient and reproducible single-cell analysis.

SUMMARY: Recent advances in single-cell analysis have introduced new computational challenges. Researchers often need to use multiple analysis tools written in different programming languages while managing version conflicts between related packages within a single workflow. For the research community, minimizing the time spent on environment setup and installation issues is essential. We present ShortCake, a containerized platform that integrates a suite of single-cell analysis tools written in R and Python. ShortCake isolates competing Python tools into separate virtual environments that can be easily accessed within a Jupyter notebook. This enables users to effortlessly transition between various environments, including R, even within a single notebook. Additionally, ShortCake offers multiple "flavors," enabling users to select container images tailored to their specific needs. ShortCake provides a unified environment with fixed versions of various tools, thus streamlining workflows, reducing setup time, and improving reproducibility. AVAILABILITY AND IMPLEMENTATION: The ShortCake image is available on DockerHub (https://hub.docker.com/r/rnakato/shortcake) and Zenodo (DOIs: 10.5281/zenodo.17116765 and 10.5281/zenodo.17118158). The source code is available on GitHub (https://github.com/rnakato/ShortCake).

Single-Cell Analysis↗

PubMed2023

Discussion on the mechanism of Lingguizhugan Decoction in treating hypertension based on network pharmacology and molecular simulation technology.

To explore the mechanism of Lingguizhugan Decoction in treating hypertension based on network pharmacology and molecular simulation. The active ingredients and potential targets were screened by the Systematic Pharmacological Analysis Platform of Traditional Chinese Medicine (TCMSP). Hypertension-related targets were obtained from OMIM and GeneCards databases. Common targets between drug and hypertension were screened in the Venny platform. A protein-protein interaction (PPI) network was constructed in the STRING database using intersection targets. Key targets in PPI network were analyzed by Cytoscape. R language program was used for Gene Ontology (GO) functional annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. Finally, the binding abilities of the main active ingredients to critical targets were verified by molecular simulation. Naringenin, quercetin, kaempferol, and β-sitosterol in Lingguizhugan Decoction, and potential targets such as STAT3, AKT1, TNF, IL6, JUN, PTGS2, MMP9, CASP3, TP53, and MAPK3, were screened out. KEGG Enrichment analysis revealed that the common targets of Lingguizhugan Decoction and hypertension are mainly involved in the lipid and atherosclerosis signaling pathway, AGE-RAGE signaling pathway in diabetic complications, fluid shear stress and atherosclerosis, and IL17 signaling pathway. The molecular simulation results showed that naringenin-MAPK3, quercetin-MMP9, quercetin-PTGS2, and quercetin-TP53 were the top four in the docking scores. Naringenin-MAPK3 and quercetin-MMP9 were stable, with binding free energies of -27.97 ± 1.41 kcal/mol and -21.15 ± 3.17 kcal/mol, respectively. The possible mechanism of Lingguizhugan Decoction in treating hypertension is characterized of multi-component, multi-target, and multi-pathway.Communicated by Ramaswamy H. Sarma.

Network Pharmacology↗

PubMed2025

Realfreq: real-time base modification analysis for nanopore sequencing.

SUMMARY: Nanopore sequencers allow sequencing data to be accessed in real-time. This allows live analysis to be performed, while the sequencing is running, reducing the turnaround time of the results. We introduce realfreq, a framework for obtaining real-time base modification frequencies while a nanopore sequencer is in operation. Realfreq calculates and allows access to the real-time base modification frequency results while the sequencer is running. We demonstrate that the data analysis rate with realfreq on a laptop computer can keep up with the output data rate of a nanopore MinION sequencer, while a desktop computer can keep up with a single PromethION 2 solo flowcell. AVAILABILITY AND IMPLEMENTATION: Realfreq is a free and open-source application implemented in C programming language and shell scripts. The source code and the documentation for realfreq can be found at https://github.com/imsuneth/realfreq. The version used for the manuscript is also available at https://doi.org/10.5281/zenodo.15128668.

Nanopore Sequencing↗

PubMed2026

A novel method for across-chromosome phasing without relative data.

MOTIVATION: Across-chromosome phasing identifies which haplotypes of different chromosomes come from the same parent. This differs from within-chromosome phasing, which uses linkage disequilibrium patterns to determine which alleles were co-inherited within each chromosome but does not match haplotypes across different chromosomes. While across-chromosome phasing can be conducted using genotypes from parents or close relatives, current methods perform poorly for samples of unrelated individuals. Here, we introduce a novel approach for across-chromosome phasing that employs a window-based SNP-similarity metric, eliminating the need for data from close relatives or detection of identical-by-descent haplotypes. RESULTS: Using UK Biobank offspring with both parents genotyped as a gold standard, we evaluated the performance of our method by phasing the offspring without using parental data. In genomic data with no within-chromosome phase errors, our algorithm achieved a mean across-chromosome phasing accuracy of 95%, with 53% of individuals phased perfectly. When data was pre-phased computationally using a standard within-chromosome phasing algorithm, mean accuracy for across-chromosome phasing dropped to 83.1%. Thus, our method is limited primarily by the accuracy of within-chromosome phasing accuracy and can approach near-perfect across-chromosome phasing accuracy as within-chromosome phasing accuracy improves. AVAILABILITY AND IMPLEMENTATION: The implementation was executed within a multi-node computational environment of University of Colorado Boulder Research Computing (Blanca Cluster: https://www.colorado.edu/rc/resources/blanca), employing parallelization techniques in the C programming language. The source code has been made publicly accessible online at https://github.com/emmanuelsapin/AcrossChromosomesPhasing, thereby facilitating reproducibility of the results for researchers with authorized access to the UK Biobank dataset.

Algorithms↗

PubMed2026

An interpretable deep learning framework uncovers features governing CRISPR-Cas9 genome-editing efficiency.

MOTIVATION: CRISPR-Cas9 genome-editing efficiency is strongly influenced by the sequence composition and positional context of single-guide RNAs (sgRNAs). Although numerous deep learning-based models have been developed to predict Cas9 efficiency from sgRNA sequences, most operate as black boxes, offering limited insight into the sequence determinants underlying Cas9 activity. In addition, previous studies often overlook how the positional context of sequence motifs within sgRNAs influences their effects on Cas9 binding or cleavage. RESULTS: We introduce DeepCC9, an interpretable machine learning framework that combines explicit sequence feature extraction with a residual block-based deep architecture to improve interpretability and identify composition- and position-based motifs governing Cas9 genome-editing efficiency. We applied this method to multiple Cas9 variant datasets, achieving superior predictive performance compared with existing methods while enabling direct interpretation of sequence motifs and their positional effects. Our analysis uncovered 74 sequence motifs enriched or depleted at specific positions within sgRNAs and strongly associated with Cas9 efficiency, providing mechanistic insight into sequence features that influence guide performance. Together, these results establish DeepCC9 as a generalizable and interpretable framework for modeling sequence-function relationships and advancing the understanding of the sequence determinants underlying CRISPR-Cas9 genome editing. AVAILABILITY AND IMPLEMENTATION: The authors have implemented their algorithm in the Python programming language (version 3.X), which is accessible using (https://zenodo.org/records/20073890).

Deep Learning↗

PubMed2026

SNPannotator: automated functional annotation of genetic variants and linked proxies.

SUMMARY: Genome-wide association studies (GWASs) have identified thousands of genetic variants associated with complex traits and diseases. However, explaining the mechanisms underlying phenotypic variation remains challenging. Here, we introduce SNPannotator, an automated post-GWAS analysis software package designed to streamline the interpretation of GWAS findings. Our pipeline implements a multi-step process that identifies proxy variants in high linkage disequilibrium (LD) with associated lead variants, then queries comprehensive resources (including Ensembl, the GTEx Portal, the eQTL Catalog, and STRING DB) for genomic position, deleteriousness, regulatory annotations, clinical significance, trait associations, expression (eQTLs) and splicing quantitative trait loci (sQTLs), and functional enrichment analyses and compiles the results into user-friendly reports. This package is implemented in the R programming language and includes auxiliary functions for variant lookup and LD exploration. SNPannotator provides a practical framework for efficiently deriving biologically meaningful insights from GWAS data and for assisting researchers in prioritizing candidate variants for functional validation. AVAILABILITY AND IMPLEMENTATION: The SNPannotator package is available from the Comprehensive R Archive Network (CRAN) at https://cran.r-project.org/web/packages/SNPannotator. The development version and tutorial is available on GitHub (https://github.com/omicslaboratory/SNPannotator). The online version of the package is available at https://omicslab.org/snpannotator.

Software↗

PubMed2025

Accurate detection of tandem repeats exposes ubiquitous reuse of biological sequences.

Tandem repetition is one of the major processes underlying genome evolution and phenotypic diversification. While newly formed tandem repeats are often easy to identify, it is more challenging to detect repeat copies as they diverge over evolutionary timescales. Existing programs for finding tandem repeats return markedly different results, and it is unclear which predictions are more correct and how much room remains for improvement. Here, we introduce DetectRepeats, a new method that uses empirical information about structural repeats to improve the accuracy of repeat detection. We show that DetectRepeats advances the state-of-the-art by finding highly divergent repeats with relatively few false positive detections. We apply DetectRepeats to genomes across the tree of life to discover an enrichment of detectable tandem repeats within different genes, genome regions, and taxa. Furthermore, we use phylogenetic reconciliation to determine that some tandem repeats continue to evolve through intra-repeat unit replacement. In this manner, tandem repeats serve as a renewable genetic resource offering a bountiful source of alternative genetic material. Our work unlocks the confident detection of ancient tandem repeats, opening a doorway to future discoveries. DetectRepeats is part of the DECIPHER package for the R programming language and available via Bioconductor.

Tandem Repeat Sequences↗

PubMed2025

A Systematic Review of Spatial Epidemiological Modeling Approaches Applied During the COVID-19 Pandemic.

BACKGROUND: A wide range of epidemiological modeling approaches have been applied to the SARS-CoV-2 pandemic, which presents an opportunity to assess common approaches applied to specific research questions. Spatial models interrogate how heterogeneities and host movement dynamics influence local and regional patterns of disease, issues that were of great interest for understanding and controlling SARS-CoV-2. OBJECTIVE: Here we present a systematic review of spatial epidemiological modeling approaches of SARS-CoV-2. We describe common themes and highlight unique strategies, providing a foundation for researchers to devise spatial models most appropriate for future pathogens and epidemics. Our review also categorizes the research questions that were addressed with spatial models, highlights parameter estimation techniques, and describes the cyber infrastructure used for model development. METHODS: We conducted a systematic review using Web of Science and a standardized set of keywords, followed by thorough examination of abstracts and full texts to determine which studies met our inclusion criteria. To guide our description and comparisons of models, we developed a Geography, Population, Movement (GPM) framework that conceptualizes the interactions between three distinct subcomponents of any spatial model. The geographic model represents the physical arena in which the model is implemented, the intra-population model describes the transmission and disease processes that occur within distinct spatial units of the geography, and the movement model describes the algorithms that dictate how hosts move among spatial units within the geography. RESULTS: The search identified a total of 193 articles, of which 109 were included in our review. The most abundant intra-population modeling methods were agent-based (47.7%) and compartmental modeling (29.4%) approaches. Movement models ranged in complexity, with the most complex models implementing commuter movement among many points of interest in the geographic arena, which were sometimes parameterized by fine-scale mobility data. Geographic models ranged from describing microcosms, such as single classrooms, all the way up to multi-country models. Of the 63.3% of models studies that specified the programming language used, we detected ten different languages, with Matlab and Python being the most frequent, although only 30.6% of studies provided open-access code for their models. We also described eight specialized software systems that were used to construct agent-based or compartment models of COVID-19. CONCLUSIONS: Our review identified and characterized a variety of spatial modeling strategies and software that were usefully employed to address many relevant epidemiological questions for COVID-19. Future research is needed to quantitatively assess which modeling approaches are most appropriate in specific situations, to answer specific questions, or to apply to certain disease systems. Moreover, future cyberinfrastructure could help to modularize and standardize modeling approaches, which would increase transparency and reproducibility, and which would facilitate a detailed examination of which model attributes relate to model performance in a variety of contexts.

COVID-19↗

PubMed2026

Identification of NR4A2 as a Potential Predictive Biomarker for Atherosclerosis.

INTRODUCTION/OBJECTIVE: Atherosclerosis, a leading cause of death globally, is characterized by the buildup of immune cells and lipids in medium to large-sized arteries. However, its precise mechanism remains unclear. The purpose of this study is to explore innovative and reliable biomarkers as a viable approach for the identification and management of atherosclerosis. METHODS: The atherosclerosis-related datasets GSE100927 and GSE66360 were retrieved from the Gene Expression Omnibus (GEO) database. The Limma package in the R programming language was utilized, applying the criteria of |logFC| > 1 and P < 0.05. Subsequently, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed on the 127 identified DEGs using R. Machine learning techniques were then applied to these data to explore and pinpoint potential biomarkers. The diagnostic potential of these markers was assessed via Receiver Operating Characteristic (ROC) curve analysis. Finally, western blot, real-time quantitative PCR (qRT-PCR), and immunohistochemistry (IHC) were employed to confirm the key biomarkers. RESULTS: Our research indicated that a total of 127 DEGs linked to atherosclerosis were successfully identified. Through the application of machine learning methods, eight critical genes were highlighted. Among these, Nuclear Receptor Subfamily 4 Group A Member-2 (NR4A2) emerged as the most promising marker for further investigation. CIBERSORT analysis revealed that NR4A2 expression levels were significantly correlated with multiple immune cell types, including B cells, plasma cells, and macrophages. Additional validation experiments confirmed that NR4A2 expression was indeed elevated in atherosclerotic plaques, supporting its potential as a biomarker for atherosclerosis. CONCLUSION: Our study identified NR4A2 as a potential immune-related biomarker for the diagnosis and treatment of atherosclerosis.

Atherosclerosis↗

PubMed2026

Weight regain following discontinuation of glucagon-like peptide-1 receptor agonists in adults who are overweight or obese: a systematic review and meta-analysis.

OBJECTIVE: This study aims to explore the effects of glucagon-like peptide-1 receptor agonists (GLP-1RAs) on weight changes and the occurrence of adverse reactions in overweight or obese adults after drug withdrawal. METHODS: Computerized searches were conducted in evidence-based databases such as PubMed, Embase, Cochrane Library and Scopus. The search period was from the establishment of the database to December 2025. Collect randomised controlled trials (RCTs) and controlled trials on GLP-1RAs, including tirzepatide, semaglutide, liraglutide, and dulaglutide, for the treatment of overweight or obese adult patients. The risk of bias in the included studies was assessed using the Cochrane Risk of Bias V2.0 tool provided by the Cochrane Collaboration, and meta-analysis was performed using the R programming language. RESULTS: A total of 699 studies were initially retrieved. Eventually, six studies involving 8,993 patients were included in the quantitative analysis, comprising 5,553 patients in the discontinuation group and 3,440 in the continued treatment group. The results of the meta-analysis showed that, compared with the continued treatment group, the weight difference in the discontinuation group was mean difference (MD) = 17.90%, 95% confidence interval (CI) [14.11-21.69], P < 0.0001. It can be seen that there was a significant rebound in weight after drug withdrawal, and there was statistical heterogeneity among the studies (P = 0.0082). Subgroup analysis further revealed that the weight rebound amplitude after discontinuation of tirzepatide was significantly higher than that of semaglutide. This result suggests that the differences in the mechanism of action of different GLP-1RAs may be the reason for the differences in weight changes after discontinuation. In addition, the percentage difference in body weight between after and before drug withdrawal was MD = 9.11%, 95% CI [7.91-10.30], P < 0.0001, further verifying the trend of weight rebound after drug withdrawal. The summary of adverse reaction reports analyzed and studied indicates that after drug withdrawal, the overall adverse reactions of patients decreased, gastrointestinal adverse reactions decreased, and the incidence of cardiovascular events was not affected by drug withdrawal. CONCLUSION: There is a significant weight rebound phenomenon after discontinuation of GLP-1RAs, and the rebound magnitudes vary among different types of drugs. At the same time, there is a risk of adverse reactions during the use of such drugs.

Humans↗

PubMed2025

AutoPM3: enhancing variant interpretation via LLM-driven PM3 evidence extraction from scientific literature.

MOTIVATION: Rare diseases affect over 300 million people worldwide and are often caused by genetic variants. While variant detection has become cost-effective, interpreting these variants-particularly collecting literature-based evidence like ACMG/AMP PM3-remains complex and time-consuming. RESULTS: We present AutoPM3, a method that automates PM3 evidence extraction from literatures using open-source large language models (LLMs). AutoPM3 combines a Text2SQL-based variant extractor and a retrieval-augmented generation (RAG) module, enhanced by a variant-specific retriever and fine-tuned LLM, to separately process tables and text. We curated PM3-Bench, a dataset of 1027 variant-publication evidence pairs from ClinGen. On openly accessible pairs, AutoPM3 achieved 86.1% accuracy for variant hits and 72.5% recall for in trans variants-outperforming other methods, including those using larger models. We uncovered the effectiveness of AutoPM3's key modules, especially for variant-specific retriever and Text2SQL, through the sequential ablation study. AutoPM3 located evidence in 76 s, demonstrating that open-source LLMs can offer an efficient, cost-effective solution for rare disease diagnosis. AVAILABILITY AND IMPLEMENTATION: AutoPM3 is implemented and freely available under the MIT license at https://github.com/HKU-BAL/AutoPM3.

Genetic Variation↗

PubMed2025

Pre-Meta: priors-augmented retrieval for LLM-based metadata generation.

MOTIVATION: While high-throughput sequencing technologies have dramatically accelerated genomic data generation, the manual processes required for dataset annotation and metadata creation impede the efficient discovery and publication of these resources across disparate public repositories. Large language models (LLMs) have the potential to streamline dataset profiling and discovery. However, their current limitations in generalizing across specialized knowledge domains, particularly in fields such as biomedical genomics, prevent them from fully realizing this potential. This article presents Pre-Meta, an LLM-agnostic and domain-independent data annotation pipeline with an enriched retrieval procedure that leverages related priors-such as pre-generated metadata tags and ontologies-as auxiliary information to improve the accuracy of automated metadata generation. RESULTS: Validated using five selected metadata fields sampled across 1500 papers, the Pre-Meta assisted annotation experiment-without finetuning and prompt optimization-demonstrates a systemic improvement in the annotation task: shown through a 23%, 72%, and 75% accuracy gain from conventional RAG adoptions of GPT-4o mini, Llama 8B, and Mistral 7B respectively. AVAILABILITY AND IMPLEMENTATION: The code, data access, and scripts are available at: https://github.com/SINTEF-SE/LLMDap.

Metadata↗

PubMed2025

Parsing GTF and FASTA files using the eccLib Library.

SUMMARY: Leveraging the Python/C API, eccLib was developed as a high-performance library designed for parsing genomic files and analysing genomic contexts. To the best of the authors' knowledge, it is the fastest Python-based solution available. With eccLib, users can efficiently parse GTF/GFFv3 and FASTA files and utilize the provided methods for additional analysis. AVAILABILITY AND IMPLEMENTATION: This library is implemented in C and distributed under the GPL-3.0 licence. It is compatible with any system that has the Python interpreter (CPython) installed. The use of C enables numerous optimizations at both the implementation and algorithmic levels, which are either unachievable or impractical in Python.

Software↗

PubMed2026

Exploring Implementation of Cantonese Radio Broadcasting as a Mental Health Promotion Initiative for Linguistically Isolated Immigrants.

Digital-first public health efforts often miss linguistically isolated immigrant communities because of structural barriers and stigma. Ethnic legacy media can reach these groups, but keeping programs funded and operational over the long term is frequently difficult. This Practice Note examines a Cantonese-language radio health program in "Los Angeles County" that ran 128 live broadcasts across a full 12-year Chinese zodiac cycle. Instead of a top-down clinical model, the show endured by acting as an informal community navigation hub. We examine two recurring administrative frictions: anonymous off-air calls about urgent economic needs (e.g., hotel job referrals) and on-air audience corrections of mispronounced medical terms. We argue these interactions should be seen not as disruptions but as measurable indicators of structural trust that sustain programs. We offer practical guidance for recruiting undergraduates to expand reach via ethnic print newspapers, using peer-recovery milestones to amplify impact. The note concludes with a pragmatic blueprint for health educators to manage professional boundaries, preserve commercial-clinical separation, and uphold cultural safety in isolated communities.

Chinese American health↗

PubMed2025

Engineered histones reshape chromatin in human cells.

Histone proteins and their variants have been found to play crucial and specialized roles in chromatin organization and the regulation of downstream gene expression; however, the relationship between histone sequence and its effect on chromatin organization remains poorly understood, limiting our functional understanding of sequence variation between distinct subtypes and across evolution and frustrating efforts to rationally design synthetic histones that can be used to engineer specified cell states. Here, we make the first advance towards engineered histone-driven chromatin organization. By expressing libraries of sequence variants of core histones in human cells, we identify variants that dominantly modulate chromatin structure. We further interrogate variants using a combination of imaging, proteomics, and genomics to reveal both cis and trans-acting mechanisms of effect. Functional screening with transcription factor libraries identifies transcriptional programs that are facilitated by engineered histone expression. Double mutation screens combined with protein language models allow us to learn sequence-to-function patterns and design synthetic histone proteins optimized to drive specific chromatin states. This work establishes a foundation for the high-throughput evaluation and engineering of chromatin-associated proteins and positions histones as tunable nodes for understanding and modulating mesoscale chromatin organization.

Journal Article↗