A single-cell atlas of chromatin accessibility in the human genome.

Current catalogs of regulatory sequences in the human genome are still incomplete and lack cell type resolution. To profile the activity of gene regulatory elements in diverse cell types and tissues in the human body, we applied single-cell chromatin accessibility assays to 30 adult human tissue types from multiple donors. We integrated these datasets with previous single-cell chromatin accessibility data from 15 fetal tissue types to reveal the status of open chromatin for ∼1.2 million candidate cis-regulatory elements (cCREs) in 222 distinct cell types comprised of >1.3 million nuclei. We used these chromatin accessibility maps to delineate cell-type-specificity of fetal and adult human cCREs and to systematically interpret the noncoding variants associated with complex human traits and diseases. This rich resource provides a foundation for the analysis of gene regulatory programs in human cell types across tissues, life stages, and organ systems.

Concurrent inhibition of oncogenic and wild-type RAS-GTP for cancer therapy.

RAS oncogenes (collectively NRAS, HRAS and especially KRAS) are among the most frequently mutated genes in cancer, with common driver mutations occurring at codons 12, 13 and 611. Small molecule inhibitors of the KRAS(G12C) oncoprotein have demonstrated clinical efficacy in patients with multiple cancer types and have led to regulatory approvals for the treatment of non-small cell lung cancer2,3. Nevertheless, KRASG12C mutations account for only around 15% of KRAS-mutated cancers4,5, and there are no approved KRAS inhibitors for the majority of patients with tumours containing other common KRAS mutations. Here we describe RMC-7977, a reversible, tri-complex RAS inhibitor with broad-spectrum activity for the active state of both mutant and wild-type KRAS, NRAS and HRAS variants (a RAS(ON) multi-selective inhibitor). Preclinically, RMC-7977 demonstrated potent activity against RAS-addicted tumours carrying various RAS genotypes, particularly against cancer models with KRAS codon 12 mutations (KRASG12X). Treatment with RMC-7977 led to tumour regression and was well tolerated in diverse RAS-addicted preclinical cancer models. Additionally, RMC-7977 inhibited the growth of KRASG12C cancer models that are resistant to KRAS(G12C) inhibitors owing to restoration of RAS pathway signalling. Thus, RAS(ON) multi-selective inhibitors can target multiple oncogenic and wild-type RAS isoforms and have the potential to treat a wide range of RAS-addicted cancers with high unmet clinical need. A related RAS(ON) multi-selective inhibitor, RMC-6236, is currently under clinical evaluation in patients with KRAS-mutant solid tumours (ClinicalTrials.gov identifier: NCT05379985).

Interpreting non-coding disease-associated human variants using single-cell epigenomics.

Genome-wide association studies (GWAS) have linked hundreds of thousands of sequence variants in the human genome to common traits and diseases. However, translating this knowledge into a mechanistic understanding of disease-relevant biology remains challenging, largely because such variants are predominantly in non-protein-coding sequences that still lack functional annotation at cell-type resolution. Recent advances in single-cell epigenomics assays have enabled the generation of cell type-, subtype- and state-resolved maps of the epigenome in heterogeneous human tissues. These maps have facilitated cell type-specific annotation of candidate cis-regulatory elements and their gene targets in the human genome, enhancing our ability to interpret the genetic basis of common traits and diseases.

Characterizing cis-regulatory elements using single-cell epigenomics.

Cell type-specific gene expression patterns and dynamics during development or in disease are controlled by cis-regulatory elements (CREs), such as promoters and enhancers. Distinct classes of CREs can be characterized by their epigenomic features, including DNA methylation, chromatin accessibility, combinations of histone modifications and conformation of local chromatin. Tremendous progress has been made in cataloguing CREs in the human genome using bulk transcriptomic and epigenomic methods. However, single-cell epigenomic and multi-omic technologies have the potential to provide deeper insight into cell type-specific gene regulatory programmes as well as into how they change during development, in response to environmental cues and through disease pathogenesis. Here, we highlight recent advances in single-cell epigenomic methods and analytical tools and discuss their readiness for human tissue profiling.

Interpreting type 1 diabetes risk with genetics and single-cell epigenomics.

Genetic risk variants that have been identified in genome-wide association studies of complex diseases are primarily non-coding1. Translating these risk variants into mechanistic insights requires detailed maps of gene regulation in disease-relevant cell types2. Here we combined two approaches: a genome-wide association study of type 1 diabetes (T1D) using 520,580 samples, and the identification of candidate cis-regulatory elements (cCREs) in pancreas and peripheral blood mononuclear cells using single-nucleus assay for transposase-accessible chromatin with sequencing (snATAC-seq) of 131,554 nuclei. Risk variants for T1D were enriched in cCREs that were active in T cells and other cell types, including acinar and ductal cells of the exocrine pancreas. Risk variants at multiple T1D signals overlapped with exocrine-specific cCREs that were linked to genes with exocrine-specific expression. At the CFTR locus, the T1D risk variant rs7795896 mapped to a ductal-specific cCRE that regulated CFTR; the risk allele reduced transcription factor binding, enhancer activity and CFTR expression in ductal cells. These findings support a role for the exocrine pancreas in the pathogenesis of T1D and highlight the power of large-scale genome-wide association studies and single-cell epigenomics for understanding the cellular origins of complex disease.

Systematic analysis of binding of transcription factors to noncoding variants.

Many sequence variants have been linked to complex human traits and diseases1, but deciphering their biological functions remains challenging, as most of them reside in noncoding DNA. Here we have systematically assessed the binding of 270 human transcription factors to 95,886 noncoding variants in the human genome using an ultra-high-throughput multiplex protein-DNA binding assay, termed single-nucleotide polymorphism evaluation by systematic evolution of ligands by exponential enrichment (SNP-SELEX). The resulting 828 million measurements of transcription factor-DNA interactions enable estimation of the relative affinity of these transcription factors to each variant in vitro and evaluation of the current methods to predict the effects of noncoding variants on transcription factor binding. We show that the position weight matrices of most transcription factors lack sufficient predictive power, whereas the support vector machine combined with the gapped k-mer representation show much improved performance, when assessed on results from independent SNP-SELEX experiments involving a new set of 61,020 sequence variants. We report highly predictive models for 94 human transcription factors and demonstrate their utility in genome-wide association studies and understanding of the molecular pathways involved in diverse human traits and diseases.

An atlas of gene regulatory elements in adult mouse cerebrum.

The mammalian cerebrum performs high-level sensory perception, motor control and cognitive functions through highly specialized cortical and subcortical structures1. Recent surveys of mouse and human brains with single-cell transcriptomics2-6 and high-throughput imaging technologies7,8 have uncovered hundreds of neural cell types distributed in different brain regions, but the transcriptional regulatory programs that are responsible for the unique identity and function of each cell type remain unknown. Here we probe the accessible chromatin in more than 800,000 individual nuclei from 45 regions that span the adult mouse isocortex, olfactory bulb, hippocampus and cerebral nuclei, and use the resulting data to map the state of 491,818 candidate cis-regulatory DNA elements in 160 distinct cell types. We find high specificity of spatial distribution for not only excitatory neurons, but also most classes of inhibitory neurons and a subset of glial cell types. We characterize the gene regulatory sequences associated with the regional specificity within these cell types. We further link a considerable fraction of the cis-regulatory elements to putative target genes expressed in diverse cerebral cell types and predict transcriptional regulators that are involved in a broad spectrum of molecular and cellular pathways in different neuronal and glial cell populations. Our results provide a foundation for comprehensive analysis of gene regulatory programs of the mammalian brain and assist in the interpretation of noncoding risk variants associated with various neurological diseases and traits in humans.

Leveraging type 1 diabetes human genetic and genomic data in the T1D knowledge portal.

To address the challenge of translating genetic discoveries for type 1 diabetes (T1D) into mechanistic insight, we have developed the T1D Knowledge Portal (T1DKP), an open-access resource for hypothesis development and target discovery in T1D.

Capillary cell-type specialization in the alveolus.

In the mammalian lung, an apparently homogenous mesh of capillary vessels surrounds each alveolus, forming the vast respiratory surface across which oxygen transfers to the blood1. Here we use single-cell analysis to elucidate the cell types, development, renewal and evolution of the alveolar capillary endothelium. We show that alveolar capillaries are mosaics; similar to the epithelium that lines the alveolus, the alveolar endothelium is made up of two intermingled cell types, with complex 'Swiss-cheese'-like morphologies and distinct functions. The first cell type, which we term the 'aerocyte', is specialized for gas exchange and the trafficking of leukocytes, and is unique to the lung. The other cell type, termed gCap ('general' capillary), is specialized to regulate vasomotor tone, and functions as a stem/progenitor cell in capillary homeostasis and repair. The two cell types develop from bipotent progenitors, mature gradually and are affected differently in disease and during ageing. This cell-type specialization is conserved between mouse and human lungs but is not found in alligator or turtle lungs, suggesting it arose during the evolution of the mammalian lung. The discovery of cell type specialization in alveolar capillaries transforms our understanding of the structure, function, regulation and maintenance of the air-blood barrier and gas exchange in health, disease and evolution.

A molecular cell atlas of the human lung from single-cell RNA sequencing.

Although single-cell RNA sequencing studies have begun to provide compendia of cell expression profiles1-9, it has been difficult to systematically identify and localize all molecular cell types in individual organs to create a full molecular cell atlas. Here, using droplet- and plate-based single-cell RNA sequencing of approximately 75,000 human cells across all lung tissue compartments and circulating blood, combined with a multi-pronged cell annotation approach, we create an extensive cell atlas of the human lung. We define the gene expression profiles and anatomical locations of 58 cell populations in the human lung, including 41 out of 45 previously known cell types and 14 previously unknown ones. This comprehensive molecular atlas identifies the biochemical functions of lung cells and the transcription factors and markers for making and monitoring them; defines the cell targets of circulating hormones and predicts local signalling interactions and immune cell homing; and identifies cell types that are directly affected by lung disease genes and respiratory viruses. By comparing human and mouse data, we identified 17 molecular cell types that have been gained or lost during lung evolution and others with substantially altered expression profiles, revealing extensive plasticity of cell types and cell-type-specific gene expression during organ evolution including expression switches between cell types. This atlas provides the molecular foundation for investigating how lung cell identities, functions and interactions are achieved in development and tissue engineering and altered in disease and evolution.

Sex-linked gene traffic underlies the acquisition of sexually dimorphic UV color vision in Heliconius butterflies.

The acquisition of novel sexually dimorphic traits poses an evolutionary puzzle: How do new traits arise and become sex-limited? Recently acquired color vision, sexually dimorphic in animals like primates and butterflies, presents a compelling model for understanding how traits become sex-biased. For example, some Heliconius butterflies uniquely possess UV (ultraviolet) color vision, which correlates with the expression of two differentially tuned UV-sensitive rhodopsins, UVRh1 and UVRh2. To discover how such traits become sexually dimorphic, we studied Heliconius charithonia, which exhibits female-specific UVRh1 expression. We demonstrate that females, but not males, discriminate different UV wavelengths. Through whole-genome shotgun sequencing and assembly of the H. charithonia genome, we discovered that UVRh1 is present on the W chromosome, making it obligately female-specific. By knocking out UVRh1, we show that UVRh1 protein expression is absent in mutant female eye tissue, as in wild-type male eyes. A PCR survey of UVRh1 sex-linkage across the genus shows that species with female-specific UVRh1 expression lack UVRh1 gDNA in males. Thus, acquisition of sex linkage is sufficient to achieve female-specific expression of UVRh1, though this does not preclude other mechanisms, like cis-regulatory evolution from also contributing. Moreover, both this event, and mutations leading to differential UV opsin sensitivity, occurred early in the history of Heliconius. These results suggest a path for acquiring sexual dimorphism distinct from existing mechanistic models. We propose a model where gene traffic to heterosomes (the W or the Y) genetically partitions a trait by sex before a phenotype shifts (spectral tuning of UV sensitivity).

Mapping genetic effects on cell type-specific chromatin accessibility and annotating complex immune trait variants using single nucleus ATAC-seq in peripheral blood.

Gene regulation is highly cell type-specific and understanding the function of non-coding genetic variants associated with complex traits requires molecular phenotyping at cell type resolution. In this study we performed single nucleus ATAC-seq (snATAC-seq) and genotyping in peripheral blood mononuclear cells from 13 individuals. Clustering chromatin accessibility profiles of 96,002 total nuclei identified 17 immune cell types and sub-types. We mapped chromatin accessibility QTLs (caQTLs) in each immune cell type and sub-type using individuals of European ancestry which identified 6,901 caQTLs at FDR < .10 and 4,220 caQTLs at FDR < .05, including those obscured from assays of bulk tissue such as with divergent effects on different cell types. For 3,941 caQTLs we further annotated putative target genes of variant activity using single cell co-accessibility, and caQTL variants were significantly correlated with the accessibility level of linked gene promoters. We fine-mapped loci associated with 16 complex immune traits and identified immune cell caQTLs at 622 candidate causal variants, including those with cell type-specific effects. At the 6q15 locus associated with type 1 diabetes, in line with previous reports, variant rs72928038 was a naïve CD4+ T cell caQTL linked to BACH2 and we validated the allelic effects of this variant on regulatory activity in Jurkat T cells. These results highlight the utility of snATAC-seq for mapping genetic effects on accessible chromatin in specific cell types.

Graded FGF activity patterns distinct cell types within the apical sensory organ of the sea anemone Nematostella vectensis.

Bilaterian animals have evolved complex sensory organs comprised of distinct cell types that function coordinately to sense the environment. Each sensory unit has a defined architecture built from component cell types, including sensory cells, non-sensory support cells, and dedicated sensory neurons. Whether this characteristic cellular composition is present in the sensory organs of non-bilaterian animals is unknown. Here, we interrogate the cell type composition and gene regulatory networks controlling development of the larval apical sensory organ in the sea anemone Nematostella vectensis. Using single cell RNA sequencing and imaging approaches, we reveal two unique cell types in the Nematostella apical sensory organ, GABAergic sensory cells and a putative non-sensory support cell population. Further, we identify the paired-like (PRD) homeodomain gene prd146 as a specific sensory cell marker and show that Prd146+ sensory cells become post-mitotic after gastrulation. Genetic loss of function approaches show that Prd146 is essential for apical sensory organ development. Using a candidate gene knockdown approach, we place prd146 downstream of FGF signaling in the apical sensory organ gene regulatory network. Further, we demonstrate that an aboral FGF activity gradient coordinately regulates the specification of both sensory and support cells. Collectively, these experiments define the genetic basis for apical sensory organ development in a non-bilaterian animal and reveal an unanticipated degree of complexity in a prototypic sensory structure.

Comprehensive analysis of the proximity-dependent nuclear interactome for the oncoprotein NOTCH1 in live cells.

Notch signaling plays a critical role in cell fate decisions in all cell types. Furthermore, gain-of-function mutations in NOTCH1 have been uncovered in many human cancers. Disruption of Notch signaling has recently emerged as an attractive disease treatment strategy. However, the nuclear interaction landscape of the oncoprotein NOTCH1 remains largely unexplored. We therefore employed here a proximity-dependent biotin identification approach to identify in vivo protein associations with the nuclear Notch1 intracellular domain in live cells. We identified a large set of previously reported and unreported proteins that associate with NOTCH1, including general transcription and elongation factors, DNA repair and replication factors, coactivators, corepressors, and components of the NuRD and SWI/SNF chromatin remodeling complexes. We also found that Notch1 intracellular domain associates with protein modifiers and components of other signaling pathways that may influence Notch signal transduction and protein stability such as USP7. We further validated the interaction of NOTCH1 with histone deacetylase 1 or GATAD2B using protein network analysis, proximity-based ligation, in vivo cross-linking and coimmunoprecipitation assays in several Notch-addicted cancer cell lines. Through data mining, we also revealed potential drug targets for the inhibition of Notch signaling. Collectively, these results provide a valuable resource to uncover the mechanisms that fine-tune Notch signaling in tumorigenesis and inform therapeutic targets for Notch-addicted tumors.

Bicistronic expression of nuclear transgenes in Chlamydomonas reinhardtii.

In eukaryotic organisms, proteins are typically translated from monocistronic messenger RNAs containing a single coding sequence (CDS). However, recent long transcript sequencing identified 87 nuclear polycistronic mRNAs in Chlamydomonas reinhardtii natively carrying multiple co-expressed CDSs. In this study, we investigated the dynamics of 22 short intergenic sequences derived from these native polycistronic loci by their application in genetic constructs for synthetic transgene expression. A promising candidate sequence was identified based on the quantification of transformation efficiency and expression strength of a fluorescence reporter protein. Subsequently, the expression of independent proteins from one mRNA was verified by cDNA amplification and protein molecular mass characterization. We demonstrated engineered bicistronic expression in vivo to drive successful co-expression of several terpene synthases with the selection marker aphVIII. Bicistronic transgene design resulted in significantly increased (E)-α-bisabolene production of 7.95 mg L-1 from a single open reading frame, 18.1× fold higher than previous reports. Use of this strategy simplifies screening procedures for identification of high-level expressing transformants, does not require the application of additional fluorescence reporters, and reduces the nucleotide footprint compared to classical monocistronic expression cassettes. Although clear advantages for bicistronic transgene expression were observed, this strategy was found to be limited to the aphVIII marker, and further studies are necessary to gain insights into the underlying mechanism that uniquely permits this co-expression from the algal nuclear genome.

The single-cell opioid responses in the context of HIV (SCORCH) consortium.

Substance use disorders (SUD) and drug addiction are major threats to public health, impacting not only the millions of individuals struggling with SUD, but also surrounding families and communities. One of the seminal challenges in treating and studying addiction in human populations is the high prevalence of co-morbid conditions, including an increased risk of contracting a human immunodeficiency virus (HIV) infection. Of the ~15 million people who inject drugs globally, 17% are persons with HIV. Conversely, HIV is a risk factor for SUD because chronic pain syndromes, often encountered in persons with HIV, can lead to an increased use of opioid pain medications that in turn can increase the risk for opioid addiction. We hypothesize that SUD and HIV exert shared effects on brain cell types, including adaptations related to neuroplasticity, neurodegeneration, and neuroinflammation. Basic research is needed to refine our understanding of these affected cell types and adaptations. Studying the effects of SUD in the context of HIV at the single-cell level represents a compelling strategy to understand the reciprocal interactions among both conditions, made feasible by the availability of large, extensively-phenotyped human brain tissue collections that have been amassed by the Neuro-HIV research community. In addition, sophisticated animal models that have been developed for both conditions provide a means to precisely evaluate specific exposures and stages of disease. We propose that single-cell genomics is a uniquely powerful technology to characterize the effects of SUD and HIV in the brain, integrating data from human cohorts and animal models. We have formed the Single-Cell Opioid Responses in the Context of HIV (SCORCH) consortium to carry out this strategy.

Sequestration from Protease Adaptor Confers Differential Stability to Protease Substrate.

According to the N-end rule, the N-terminal residue of a protein determines its stability. In bacteria, the adaptor ClpS mediates proteolysis by delivering substrates bearing specific N-terminal residues to the protease ClpAP. We now report that the Salmonella adaptor ClpS binds to the N terminus of the regulatory protein PhoP, resulting in PhoP degradation by ClpAP. We establish that the PhoP-activated protein MgtC protects PhoP from degradation by outcompeting ClpS for binding to PhoP. MgtC appears to act exclusively on PhoP, as it did not alter the stability of a different ClpS-dependent ClpAP substrate. Removal of five N-terminal residues rendered PhoP stability independent of both the clpS and mgtC genes. By preserving PhoP protein levels, MgtC enables normal temporal transcription of PhoP-activated genes. The identified mechanism provides a simple means to spare specific substrates from an adaptor-dependent protease.

Chemical Proteomics Strategies for Analyzing Protein Lipidation Reveal the Bacterial O-Mycoloylome.

Protein lipidation dynamically controls protein localization and function within cellular membranes. A unique form of protein O-fatty acylation in Corynebacterium, termed protein O-mycoloylation, involves the attachment of mycolic acids─unusually large and hydrophobic fatty acids─to serine residues of proteins in these organisms' outer mycomembrane. However, as with other forms of protein lipidation, the scope and functional consequences of protein O-mycoloylation are challenging to investigate due to the inherent difficulties of enriching and analyzing lipidated peptides. To facilitate the analysis of protein lipidation and enable the comprehensive profiling and site mapping of protein O-mycoloylation, we developed a chemical proteomics strategy integrating metabolic labeling, click chemistry, cleavable linkers, and a novel liquid chromatography-tandem mass spectrometry (LC-MS/MS) method employing LC separation and complementary fragmentation methods tailored to the analysis of lipophilic, MS-labile O-acylated peptides. Using these tools in the model organism Corynebacterium glutamicum, we identified approximately 30 candidate O-mycoloylated proteins, including porins, mycoloyltransferases, secreted hydrolases, and other proteins with cell envelope-related functions─consistent with a role for O-mycoloylation in targeting proteins to the mycomembrane. Site mapping revealed that many of the proteins contained multiple spatially proximal modification sites, which occurred predominantly at serine residues surrounded by conformationally flexible peptide motifs. Overall, this study (i) discloses the putative protein O-mycoloylome for the first time, (ii) yields new insights into the undercharacterized proteome of the mycomembrane, which is a hallmark of important pathogens (e.g., Corynebacterium diphtheriae, Mycobacterium tuberculosis), and (iii) provides generally applicable chemical strategies for the proteomic analysis of protein lipidation.

Ligand Steric Profile Tunes the Reactivity of Indium Phosphide Clusters.

Indium phosphide quantum dots have become an industrially relevant material for solid-state lighting and wide color gamut displays. The synthesis of indium phosphide quantum dots from indium carboxylates and tris(trimethylsilyl)phosphine (P(SiMe3)3) is understood to proceed through the formation of magic-sized clusters, with In37P20(O2CR)51 being the key isolable intermediate. The reactivity of the In37P20(O2CR)51 cluster is a vital parameter in controlling the conversion to quantum dots. Herein, we report structural perturbations of In37P20(O2CR)51 clusters induced by tuning the steric properties of a series of substituted phenylacetate ligands. This approach allows for control over reactivity with P(SiMe3)3, where meta-substituents enhance the susceptibility to ligand displacement, and para-substituents hinder phosphine diffusion to the core. Thermolysis studies show that with complete cluster dissolution, steric profile can modulate the nucleation period, resulting in a nanocrystal size dependence on ligand steric profile. The enhanced stability from ligand engineering also allows for the isolation and structural characterization by single-crystal X-ray diffraction of a new III-V magic-sized cluster with the formula In26P13(O2CR)39. This intermediate precedes the In37P20(O2CR)51 cluster on the InP QD reaction coordinate. The physical and electronic structure of this cluster are analyzed, providing new insight into previously unrecognized relationships between II-VI and III-V materials and the discrete growth of III-V cluster intermediates.

Velocity dependence of sensory reweighting in human balance control.

The relative contributions of proprioceptive, vestibular, and visual sensory cues to balance control change depending on their availability and reliability. This sensory reweighting is classically supported by nonlinear sway responses to increasing visual surround and/or surface tilt amplitudes. However, recent evidence indicates that visual cues are reweighted based on visual tilt velocity rather than tilt amplitude. Therefore, we designed a study to specifically test the hypothesized velocity dependence of reweighting while expanding on earlier findings for visual reweighting by testing proprioceptive reweighting for standing balance on a tilting surface. Twenty healthy young adults stood with their eyes closed on a toes-up/-down tilting platform. We designed four pseudorandom tilt sequences with either a slow (S) or a fast (F) tilt velocity and different peak-to-peak amplitudes. We used model-based interpretations of measured sway characteristics to estimate the proprioceptive sensory weight (Wprop) within each trial. In addition, root-mean-square values of measured body center of mass sway amplitude (RMS) and velocity (RMSv) were calculated for each tilt sequence. Wprop, RMS, and RMSv values varied depending on the stimulus velocity, exhibiting large effects (all Cohen's d >1.10). In contrast, we observed no significant differences across stimulus amplitudes for Wprop (Cohen's d: 0.02-0.16) and, compared with the differences in velocity, there were much smaller changes in RMS and RMSv values (Cohen's d: 0.05-0.91). These results confirmed the hypothesized velocity, rather than amplitude, dependence of sensory reweighting.NEW & NOTEWORTHY This novel study examined the velocity dependence of sensory reweighting for human balance control using support surface tilt stimuli with independently varied amplitude and velocity. Estimates of the proprioceptive contribution to standing balance, derived from model-based interpretations of sway characteristics, showed greater sensitivity to changes in surface tilt velocity than surface tilt amplitude. These results support a velocity-based mechanism underlying sensory reweighting for human balance control.