An open-source framework for end-to-end analysis of electronic health record data.

With progressive digitalization of healthcare systems worldwide, large-scale collection of electronic health records (EHRs) has become commonplace. However, an extensible framework for comprehensive exploratory analysis that accounts for data heterogeneity is missing. Here we introduce ehrapy, a modular open-source Python framework designed for exploratory analysis of heterogeneous epidemiology and EHR data. ehrapy incorporates a series of analytical steps, from data extraction and quality control to the generation of low-dimensional representations. Complemented by rich statistical modules, ehrapy facilitates associating patients with disease states, differential comparison between patient clusters, survival analysis, trajectory inference, causal inference and more. Leveraging ontologies, ehrapy further enables data sharing and training EHR deep learning models, paving the way for foundational models in biomedical research. We demonstrate ehrapy's features in six distinct examples. We applied ehrapy to stratify patients affected by unspecified pneumonia into finer-grained phenotypes. Furthermore, we reveal biomarkers for significant differences in survival among these groups. Additionally, we quantify medication-class effects of pneumonia medications on length of stay. We further leveraged ehrapy to analyze cardiovascular risks across different data modalities. We reconstructed disease state trajectories in patients with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) based on imaging data. Finally, we conducted a case study to demonstrate how ehrapy can detect and mitigate biases in EHR data. ehrapy, thus, provides a framework that we envision will standardize analysis pipelines on EHR data and serve as a cornerstone for the community.

Sfrp1 inhibits lung fibroblast invasion during transition to injury-induced myofibroblasts.

BACKGROUND: Fibroblast-to-myofibroblast conversion is a major driver of tissue remodelling in organ fibrosis. Distinct lineages of fibroblasts support homeostatic tissue niche functions, yet their specific activation states and phenotypic trajectories during injury and repair have remained unclear. METHODS: We combined spatial transcriptomics, multiplexed immunostainings, longitudinal single-cell RNA-sequencing and genetic lineage tracing to study fibroblast fates during mouse lung regeneration. Our findings were validated in idiopathic pulmonary fibrosis patient tissues in situ as well as in cell differentiation and invasion assays using patient lung fibroblasts. Cell differentiation and invasion assays established a function of SFRP1 in regulating human lung fibroblast invasion in response to transforming growth factor (TGF)ß1. MEASUREMENTS AND MAIN RESULTS: We discovered a transitional fibroblast state characterised by high Sfrp1 expression, derived from both Tcf21-Cre lineage positive and negative cells. Sfrp1 + cells appeared early after injury in peribronchiolar, adventitial and alveolar locations and preceded the emergence of myofibroblasts. We identified lineage-specific paracrine signals and inferred converging transcriptional trajectories towards Sfrp1 + transitional fibroblasts and Cthrc1 + myofibroblasts. TGFß1 downregulated SFRP1 in noninvasive transitional cells and induced their switch to an invasive CTHRC1+ myofibroblast identity. Finally, using loss-of-function studies we showed that SFRP1 modulates TGFß1-induced fibroblast invasion and RHOA pathway activity. CONCLUSIONS: Our study reveals the convergence of spatially and transcriptionally distinct fibroblast lineages into transcriptionally uniform myofibroblasts and identifies SFRP1 as a modulator of TGFß1-driven fibroblast phenotypes in fibrogenesis. These findings are relevant in the context of therapeutic interventions that aim at limiting or reversing fibroblast foci formation.

Ex vivo tissue perturbations coupled to single-cell RNA-seq reveal multilineage cell circuit dynamics in human lung fibrogenesis.

Pulmonary fibrosis develops as a consequence of failed regeneration after injury. Analyzing mechanisms of regeneration and fibrogenesis directly in human tissue has been hampered by the lack of organotypic models and analytical techniques. In this work, we coupled ex vivo cytokine and drug perturbations of human precision-cut lung slices (hPCLS) with single-cell RNA sequencing and induced a multilineage circuit of fibrogenic cell states in hPCLS. We showed that these cell states were highly similar to the in vivo cell circuit in a multicohort lung cell atlas from patients with pulmonary fibrosis. Using micro-CT-staged patient tissues, we characterized the appearance and interaction of myofibroblasts, an ectopic endothelial cell state, and basaloid epithelial cells in the thickened alveolar septum of early-stage lung fibrosis. Induction of these states in the hPCLS model provided evidence that the basaloid cell state was derived from alveolar type 2 cells, whereas the ectopic endothelial cell state emerged from capillary cell plasticity. Cell-cell communication routes in patients were largely conserved in hPCLS, and antifibrotic drug treatments showed highly cell type-specific effects. Our work provides an experimental framework for perturbational single-cell genomics directly in human lung tissue that enables analysis of tissue homeostasis, regeneration, and pathology. We further demonstrate that hPCLS offer an avenue for scalable, high-resolution drug testing to accelerate antifibrotic drug development and translation.

The arginine methyltransferase PRMT7 promotes extravasation of monocytes resulting in tissue injury in COPD.

Extravasation of monocytes into tissue and to the site of injury is a fundamental immunological process, which requires rapid responses via post translational modifications (PTM) of proteins. Protein arginine methyltransferase 7 (PRMT7) is an epigenetic factor that has the capacity to mono-methylate histones on arginine residues. Here we show that in chronic obstructive pulmonary disease (COPD) patients, PRMT7 expression is elevated in the lung tissue and localized to the macrophages. In mouse models of COPD, lung fibrosis and skin injury, reduced expression of PRMT7 associates with decreased recruitment of monocytes to the site of injury and hence less severe symptoms. Mechanistically, activation of NF-κB/RelA in monocytes induces PRMT7 transcription and consequential mono-methylation of histones at the regulatory elements of RAP1A, which leads to increased transcription of this gene that is responsible for adhesion and migration of monocytes. Persistent monocyte-derived macrophage accumulation leads to ALOX5 over-expression and accumulation of its metabolite LTB4, which triggers expression of ACSL4 a ferroptosis promoting gene in lung epithelial cells. Conclusively, inhibition of arginine mono-methylation might offer targeted intervention in monocyte-driven inflammatory conditions that lead to extensive tissue damage if left untreated.

IL-37 regulates allergic inflammation by counterbalancing pro-inflammatory IL-1 and IL-33.

BACKGROUND: Children with asthma have impaired production of interleukin (IL) 37; in mice, IL-37 reduces hallmarks of experimental allergic asthma (EAA). However, it remains unclear how IL-37 exerts its inhibitory properties in asthma. This study aimed to identify the mechanism(s) by which IL-37 controls allergic inflammation. METHODS: IL-37 target cells were identified by single-cell RNA-seq of IL-1R5 and IL-1R8. Airway tissues were isolated by laser-capture microdissection and examined by microarray-based gene expression analysis. Mononuclear cells (MNC) and airway epithelial cells (AECs) were isolated and stimulated with allergen, IL-1ß, or IL-33 together with recombinant human (rh) IL-37. Wild-type, IL-1R1- and IL-33-deficient mice with EAA were treated with rhIL-37. IL-1ß, IL-33, and IL-37 levels were determined in sputum and nasal secretions from adult asthma patients without glucocorticoid therapy. RESULTS: IL-37 target cells included AECs, T cells, and dendritic cells. In mice with EAA, rhIL-37 led to differential expression of >90 genes induced by IL-1ß and IL-33. rhIL-37 reduced production of Th2 cytokines in allergen-activated MNCs from wild-type but not from IL-1R1-deficient mice and inhibited IL-33-induced Th2 cytokine release. Furthermore, rhIL-37 attenuated IL-1ß- and IL-33-induced pro-inflammatory mediator expression in murine AEC cultures. In contrast to wild-type mice, hIL-37 had no effect on EAA in IL-1R1- or IL-33-deficient mice. We also observed that expression/production ratios of both IL-1ß and IL-33 to IL-37 were dramatically increased in asthma patients compared to healthy controls. CONCLUSION: IL-37 downregulates allergic airway inflammation by counterbalancing the disease-amplifying effects of IL-1ß and IL-33.

Single-cell RNA sequencing reveals ex vivo signatures of SARS-CoV-2-reactive T cells through 'reverse phenotyping'.

The in vivo phenotypic profile of T cells reactive to severe acute respiratory syndrome (SARS)-CoV-2 antigens remains poorly understood. Conventional methods to detect antigen-reactive T cells require in vitro antigenic re-stimulation or highly individualized peptide-human leukocyte antigen (pHLA) multimers. Here, we use single-cell RNA sequencing to identify and profile SARS-CoV-2-reactive T cells from Coronavirus Disease 2019 (COVID-19) patients. To do so, we induce transcriptional shifts by antigenic stimulation in vitro and take advantage of natural T cell receptor (TCR) sequences of clonally expanded T cells as barcodes for 'reverse phenotyping'. This allows identification of SARS-CoV-2-reactive TCRs and reveals phenotypic effects introduced by antigen-specific stimulation. We characterize transcriptional signatures of currently and previously activated SARS-CoV-2-reactive T cells, and show correspondence with phenotypes of T cells from the respiratory tract of patients with severe disease in the presence or absence of virus in independent cohorts. Reverse phenotyping is a powerful tool to provide an integrated insight into cellular states of SARS-CoV-2-reactive T cells across tissues and activation states.

An intramembrane aromatic network determines pentameric assembly of Cys-loop receptors.

Cys-loop receptors are pentameric ligand-gated ion channels (pLGICs) that mediate fast synaptic transmission. Here functional pentameric assembly of truncated fragments comprising the ligand-binding N-terminal ectodomains and the first three transmembrane helices, M1-M3, of both the inhibitory glycine receptor (GlyR) alpha1 and the 5HT(3)A receptor subunits was found to be rescued by coexpressing the complementary fourth transmembrane helix, M4. Alanine scanning identified multiple aromatic residues in M1, M3 and M4 as key determinants of GlyR assembly. Homology modeling and molecular dynamics simulations revealed that these residues define an interhelical aromatic network, which we propose determines the geometry of M1-M4 tetrahelical packing such that nascent pLGIC subunits must adopt a closed fivefold symmetry. Because pLGIC ectodomains form random nonstoichiometric oligomers, proper pentameric assembly apparently depends on intersubunit interactions between extracellular domains and intrasubunit interactions between transmembrane segments.

A trimeric quaternary structure is conserved in bacterial and human glutamate transporters.

Neuronal and glial glutamate transporters play a central role in the termination of synaptic transmission and in extracellular glutamate homeostasis in the mammalian central nervous system. They are known to be multimers; however, the number of subunits forming a functional transporter is controversial. We studied the subunit stoichiometry of two distantly related glutamate transporters, the human glial glutamate transporter hEAAT2 and a bacterial glutamate transporter from Escherichia coli, ecgltP. Using blue native polyacrylamide gel electrophoresis, analysis of concatenated transporters, and chemical cross-linking, we demonstrated that human and prokaryotic glutamate transporters expressed in Xenopus laevis oocytes or in mammalian cells are assembled as trimers composed of three identical subunits. In an inducible mammalian cell line expressing hEAAT2 the glutamate uptake currents correlate to the amount of trimeric transporters. Overexpression and purification of ecgltP in E. coli resulted in a homogenous population of trimeric transporters that were functional after reconstitution in lipid vesicles. Our results indicate that an evolutionarily conserved trimeric quaternary structure represents the sole native and functional state of glutamate transporters.