Transcriptional profile of the rat cardiovascular system at single cell resolution

Background

Despite the critical role of the cardiovascular system, our understanding of its cellular and transcriptional diversity remains limited. We therefore sought to characterize the cellular composition, phenotypes, molecular pathways, and communication networks between cell types at the tissue and sub-tissue level across the cardiovascular system of the healthy Wistar rat, an important model in preclinical cardiovascular research. We obtained high quality tissue samples under controlled conditions that reveal a level of cellular detail so far inaccessible in human studies.

Methods and Results

We performed single nucleus RNA-sequencing in 78 samples in 10 distinct regions including the four chambers of the heart, ventricular septum, sinoatrial node, atrioventricular node, aorta, pulmonary artery, and pulmonary veins (PV), which produced an aggregate map of 505,835 nuclei. We identified 26 distinct cell types and additional subtypes, including a number of rare cell types such as PV cardiomyocytes and non-myelinating Schwann cells (NMSCs), and unique groups of vascular smooth muscle cells (VSMCs), endothelial cells (ECs) and fibroblasts (FBs), which gave rise to a detailed cell type distribution across tissues. We demonstrated differences in the cellular composition across different cardiac regions and tissue-specific differences in transcription for each cell type, highlighting the molecular diversity and complex tissue architecture of the cardiovascular system. Specifically, we observed great transcriptional heterogeneities among ECs and FBs. Importantly, several cell subtypes had a unique regional localization such as a subtype of VSMCs enriched in the large vasculature. We found the cellular makeup of PV tissue is closer to heart tissue than to the large arteries. We further explored the ligand-receptor repertoire across cell clusters and tissues, and observed tissue-enriched cellular communication networks, including heightened Nppa - Npr1 / 2 / 3 signaling in the sinoatrial node.

Conclusions

Through a large single nucleus sequencing effort encompassing over 500,000 nuclei, we broadened our understanding of cellular transcription in the healthy cardiovascular system. The existence of tissue-restricted cellular phenotypes suggests regional regulation of cardiovascular physiology. The overall conservation in gene expression and molecular pathways across rat and human cell types, together with our detailed transcriptional characterization of each cell type, offers the potential to identify novel therapeutic targets and improve preclinical models of cardiovascular disease.

Single-nucleus profiling of human dilated and hypertrophic cardiomyopathy

Heart failure encompasses a heterogeneous set of clinical features that converge on impaired cardiac contractile function1,2 and presents a growing public health concern. Previous work has highlighted changes in both transcription and protein expression in failing hearts3,4, but may overlook molecular changes in less prevalent cell types. Here we identify extensive molecular alterations in failing hearts at single-cell resolution by performing single-nucleus RNA sequencing of nearly 600,000 nuclei in left ventricle samples from 11 hearts with dilated cardiomyopathy and 15 hearts with hypertrophic cardiomyopathy as well as 16 non-failing hearts. The transcriptional profiles of dilated or hypertrophic cardiomyopathy hearts broadly converged at the tissue and cell-type level. Further, a subset of hearts from patients with cardiomyopathy harbour a unique population of activated fibroblasts that is almost entirely absent from non-failing samples. We performed a CRISPR-knockout screen in primary human cardiac fibroblasts to evaluate this fibrotic cell state transition; knockout of genes associated with fibroblast transition resulted in a reduction of myofibroblast cell-state transition upon TGFβ1 stimulation for a subset of genes. Our results provide insights into the transcriptional diversity of the human heart in health and disease as well as new potential therapeutic targets and biomarkers for heart failure.

MACA: marker-based automatic cell-type annotation for single-cell expression data

SUMMARY

Accurately identifying cell types is a critical step in single-cell sequencing analyses. Here, we present marker-based automatic cell-type annotation (MACA), a new tool for annotating single-cell transcriptomics datasets. We developed MACA by testing four cell-type scoring methods with two public cell-marker databases as reference in six single-cell studies. MACA compares favorably to four existing marker-based cell-type annotation methods in terms of accuracy and speed. We show that MACA can annotate a large single-nuclei RNA-seq study in minutes on human hearts with ∼290K cells. MACA scales easily to large datasets and can broadly help experts to annotate cell types in single-cell transcriptomics datasets, and we envision MACA provides a new opportunity for integration and standardization of cell-type annotation across multiple datasets.

AVAILABILITY AND IMPLEMENTATION

MACA is written in python and released under GNU General Public License v3.0. The source code is available at https://github.com/ImXman/MACA.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Deep learning enables genetic analysis of the human thoracic aorta

Enlargement or aneurysm of the aorta predisposes to dissection, an important cause of sudden death. We trained a deep learning model to evaluate the dimensions of the ascending and descending thoracic aorta in 4.6 million cardiac magnetic resonance images from the UK Biobank. We then conducted genome-wide association studies in 39,688 individuals, identifying 82 loci associated with ascending and 47 with descending thoracic aortic diameter, of which 14 loci overlapped. Transcriptome-wide analyses, rare-variant burden tests and human aortic single nucleus RNA sequencing prioritized genes including SVIL, which was strongly associated with descending aortic diameter. A polygenic score for ascending aortic diameter was associated with thoracic aortic aneurysm in 385,621 UK Biobank participants (hazard ratio = 1.43 per s.d., confidence interval 1.32-1.54, P = 3.3 × 10-20). Our results illustrate the potential for rapidly defining quantitative traits with deep learning, an approach that can be broadly applied to biomedical images.

Myocyte-Specific Upregulation of ACE2 in Cardiovascular Disease: Implications for SARS-CoV-2-Mediated Myocarditis

Supplemental Digital Content is available in the text.

Transcriptional and Cellular Diversity of the Human Heart

BACKGROUND

The human heart requires a complex ensemble of specialized cell types to perform its essential function. A greater knowledge of the intricate cellular milieu of the heart is critical to increase our understanding of cardiac homeostasis and pathology. As recent advances in low-input RNA sequencing have allowed definitions of cellular transcriptomes at single-cell resolution at scale, we have applied these approaches to assess the cellular and transcriptional diversity of the nonfailing human heart.

METHODS

Microfluidic encapsulation and barcoding was used to perform single nuclear RNA sequencing with samples from 7 human donors, selected for their absence of overt cardiac disease. Individual nuclear transcriptomes were then clustered based on transcriptional profiles of highly variable genes. These clusters were used as the basis for between-chamber and between-sex differential gene expression analyses and intersection with genetic and pharmacologic data.

RESULTS

We sequenced the transcriptomes of 287 269 single cardiac nuclei, revealing 9 major cell types and 20 subclusters of cell types within the human heart. Cellular subclasses include 2 distinct groups of resident macrophages, 4 endothelial subtypes, and 2 fibroblast subsets. Comparisons of cellular transcriptomes by cardiac chamber or sex reveal diversity not only in cardiomyocyte transcriptional programs but also in subtypes involved in extracellular matrix remodeling and vascularization. Using genetic association data, we identified strong enrichment for the role of cell subtypes in cardiac traits and diseases. Intersection of our data set with genes on cardiac clinical testing panels and the druggable genome reveals striking patterns of cellular specificity.

CONCLUSIONS

Using large-scale single nuclei RNA sequencing, we defined the transcriptional and cellular diversity in the normal human heart. Our identification of discrete cell subtypes and differentially expressed genes within the heart will ultimately facilitate the development of new therapeutics for cardiovascular diseases.

Myocyte Specific Upregulation of ACE2 in Cardiovascular Disease: Implications for SARS-CoV-2 mediated myocarditis

Coronavirus disease 2019 (COVID-19) is a global pandemic caused by a novel severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2). SARS-CoV-2 infection of host cells occurs predominantly via binding of the viral surface spike protein to the human angiotensin-converting enzyme 2 (ACE2) receptor. Hypertension and pre-existing cardiovascular disease are risk factors for morbidity from COVID-19, and it remains uncertain whether the use of angiotensin converting enzyme inhibitors (ACEi) or angiotensin receptor blockers (ARB) impacts infection and disease. Here, we aim to shed light on this question by assessing ACE2 expression in normal and diseased human myocardial samples profiled by bulk and single nucleus RNA-seq.

Conformational Adaption May Explain the Slow Dissociation Kinetics of Roniciclib (BAY 1000394), a Type I CDK Inhibitor with Kinetic Selectivity for CDK2 and CDK9

Roniciclib (BAY 1000394) is a type I pan-CDK (cyclin-dependent kinase) inhibitor which has revealed potent efficacy in xenograft cancer models. Here, we show that roniciclib displays prolonged residence times on CDK2 and CDK9, whereas residence times on other CDKs are transient, thus giving rise to a kinetic selectivity of roniciclib. Surprisingly, variation of the substituent at the 5-position of the pyrimidine scaffold results in changes of up to 3 orders of magnitude of the drug-target residence time. CDK2 X-ray cocrystal structures have revealed a DFG-loop adaption for the 5-(trifluoromethyl) substituent, while for hydrogen and bromo substituents the DFG loop remains in its characteristic type I inhibitor position. In tumor cells, the prolonged residence times of roniciclib on CDK2 and CDK9 are reflected in a sustained inhibitory effect on retinoblastoma protein (RB) phosphorylation, indicating that the target residence time on CDK2 may contribute to sustained target engagement and antitumor efficacy.

Functional and Structural Characterization of Bub3·BubR1 Interactions Required for Spindle Assembly Checkpoint Signaling in Human Cells

The spindle assembly checkpoint (SAC) is an essential safeguarding mechanism devised to ensure equal chromosome distribution in daughter cells upon mitosis. The proteins Bub3 and BubR1 are key components of the mitotic checkpoint complex, an essential part of the molecular machinery on which the SAC relies. In the present work we have performed a detailed functional and biochemical characterization of the interaction between human Bub3 and BubR1 in cells and in vitro Our results demonstrate that genetic knockdown of Bub3 abrogates the SAC, promotes apoptosis, and inhibits the proliferation of human cancer cells. We also show that the integrity of the human mitotic checkpoint complex depends on the specific recognition between BubR1 and Bub3, for which the BubR1 Gle2 binding sequence motif is essential. This 1:1 binding event is high affinity, enthalpy-driven and with slow dissociation kinetics. The affinity, kinetics, and thermodynamic parameters of the interaction are differentially modulated by small regions in the N and C termini of the Gle2 binding domain sequence, suggesting the existence of "hotspots" for this protein-protein interaction. Furthermore, we show that specific disruption of endogenous BubR1·Bub3 complexes in human cancer cells phenocopies the effects observed in gene targeting experiments. Our work enhances the current understanding of key members of the SAC and paves the road for the pursuit of novel targeted cancer therapies based on SAC inhibition.

Discovery of a Potent Class I Protein Arginine Methyltransferase Fragment Inhibitor

Protein methyltransferases (PMTs) are a promising target class in oncology and other disease areas. They are composed of SET domain methyltransferases and structurally unrelated Rossman-fold enzymes that include protein arginine methyltransferases (PRMTs). In the absence of a well-defined medicinal chemistry tool-kit focused on PMTs, most current inhibitors were identified by screening large and diverse libraries of leadlike molecules. So far, no successful fragment-based approach was reported against this target class. Here, by deconstructing potent PRMT inhibitors, we find that chemical moieties occupying the substrate arginine-binding site can act as efficient fragment inhibitors. Screening a fragment library against PRMT6 produced numerous hits, including a 300 nM inhibitor (ligand efficiency of 0.56) that decreased global histone 3 arginine 2 methylation in cells, and can serve as a warhead for the development of PRMT chemical probes.

An autoinhibited state in the structure of Thermotoga maritima NusG

NusG is a conserved regulatory protein interacting with RNA polymerase (RNAP) and other proteins to form multicomponent complexes that modulate transcription. The crystal structure of Thermotoga maritima NusG (TmNusG) shows a three-domain architecture, comprising well-conserved amino-terminal (NTD) and carboxy-terminal (CTD) domains with an additional, species-specific domain inserted into the NTD. NTD and CTD directly contact each other, occluding a surface of the NTD for binding to RNAP and a surface on the CTD interacting either with transcription termination factor Rho or transcription antitermination factor NusE. NMR spectroscopy confirmed the intramolecular NTD-CTD interaction up to the optimal growth temperature of Thermotoga maritima. The domain interaction involves a dynamic equilibrium between open and closed states and contributes significantly to the overall fold stability of the protein. Wild-type TmNusG and deletion variants could not replace endogenous Escherichia coli NusG, suggesting that the NTD-CTD interaction of TmNusG represents an autoinhibited state.

The crystal structure of PPIL1 bound to cyclosporine A suggests a binding mode for a linear epitope of the SKIP protein

Background

The removal of introns from pre-mRNA is carried out by a large macromolecular machine called the spliceosome. The peptidyl-prolyl cis/trans isomerase PPIL1 is a component of the human spliceosome and binds to the spliceosomal SKIP protein via a binding site distinct from its active site.

Principal Findings

Here, we have studied the PPIL1 protein and its interaction with SKIP biochemically and by X-ray crystallography. A minimal linear binding epitope derived from the SKIP protein could be determined using a peptide array. A 36-residue region of SKIP centred on an eight-residue epitope suffices to bind PPIL1 in pull-down experiments. The crystal structure of PPIL1 in complex with the inhibitor cyclosporine A (CsA) was obtained at a resolution of 1.15 Å and exhibited two bound Cd²⁺ ions that enabled SAD phasing. PPIL1 residues that have previously been implicated in binding of SKIP are involved in the coordination of Cd²⁺ ions in the present crystal structure. Employing the present crystal structure, the determined minimal binding epitope and previously published NMR data [1], a molecular docking study was performed. In the docked model of the PPIL1·SKIP interaction, a proline residue of SKIP is buried in a hydrophobic pocket of PPIL1. This hydrophobic contact is encircled by several hydrogen bonds between the SKIP peptide and PPIL1.

Conclusion

We characterized a short, linear epitope of SKIP that is sufficient to bind the PPIL1 protein. Our data indicate that this SKIP peptide could function in recruiting PPIL1 into the core of the spliceosome. We present a molecular model for the binding mode of SKIP to PPIL1 which emphasizes the versatility of cyclophilin-type PPIases to engage in additional interactions with other proteins apart from active site contacts despite their limited surface area.

Structural and biochemical studies of botulinum neurotoxin serotype C1 light chain protease: implications for dual substrate specificity

Clostridial neurotoxins are the causative agents of the neuroparalytic disease botulism and tetanus. They block neurotransmitter release through specific proteolysis of one of the three soluble N-ethylmaleimide-sensitive-factor attachment protein receptors (SNAREs) SNAP-25, syntaxin, and synaptobrevin, which constitute part of the synaptic vesicle fusion machinery. The catalytic component of the clostridial neurotoxins is their light chain (LC), a Zn2+ endopeptidase. There are seven structurally and functionally related botulinum neurotoxins (BoNTs), termed serotype A to G, and tetanus neurotoxin (TeNT). Each of them exhibits unique specificity for their target SNAREs and peptide bond(s) they cleave. The mechanisms of action for substrate recognition and target cleavage are largely unknown. Here, we report structural and biochemical studies of BoNT/C1-LC, which is unique among BoNTs in that it exhibits dual specificity toward both syntaxin and SNAP-25. A distinct pocket (S1') near the active site likely achieves the correct register for the cleavage site by only allowing Ala as the P1' residue for both SNAP-25 and syntaxin. Mutations of this SNAP-25 residue dramatically reduce enzymatic activity. The remote alpha-exosite that was previously identified in the complex of BoNT/A-LC and SNAP-25 is structurally conserved in BoNT/C1. However, mutagenesis experiments show that the alpha-exosite of BoNT/C1 plays a less stringent role in substrate discrimination in comparison to that of BoNT/A, which could account for its dual substrate specificity.

Inhibition of metalloprotease botulinum serotype A from a pseudo-peptide binding mode to a small molecule that is active in primary neurons

An efficient research strategy integrating empirically guided, structure-based modeling and chemoinformatics was used to discover potent small molecule inhibitors of the botulinum neurotoxin serotype A light chain. First, a modeled binding mode for inhibitor 2-mercapto-3-phenylpropionyl-RATKML (K(i) = 330 nM) was generated, and required the use of a molecular dynamic conformer of the enzyme displaying the reorientation of surface loops bordering the substrate binding cleft. These flexible loops are conformationally variable in x-ray crystal structures, and the model predicted that they were pivotal for providing complementary binding surfaces and solvent shielding for the pseudo-peptide. The docked conformation of 2-mercapto-3-phenylpropionyl-RATKML was then used to refine our pharmacophore for botulinum serotype A light chain inhibition. Data base search queries derived from the pharmacophore were employed to mine small molecule (non-peptidic) inhibitors from the National Cancer Institute's Open Repository. Four of the inhibitors possess K(i) values ranging from 3.0 to 10.0 microM. Of these, NSC 240898 is a promising lead for therapeutic development, as it readily enters neurons, exhibits no neuronal toxicity, and elicits dose-dependent protection of synaptosomal-associated protein (of 25 kDa) in a primary culture of embryonic chicken neurons. Isothermal titration calorimetry showed that the interaction between NSC 240898 and the botulinum A light chain is largely entropy-driven, and occurs with a 1:1 stoichiometry and a dissociation constant of 4.6 microM.