Spatially resolved multiomics of human cardiac niches.

The function of a cell is defined by its intrinsic characteristics and its niche: the tissue microenvironment in which it dwells. Here we combine single-cell and spatial transcriptomics data to discover cellular niches within eight regions of the human heart. We map cells to microanatomical locations and integrate knowledge-based and unsupervised structural annotations. We also profile the cells of the human cardiac conduction system1. The results revealed their distinctive repertoire of ion channels, G-protein-coupled receptors (GPCRs) and regulatory networks, and implicated FOXP2 in the pacemaker phenotype. We show that the sinoatrial node is compartmentalized, with a core of pacemaker cells, fibroblasts and glial cells supporting glutamatergic signalling. Using a custom CellPhoneDB.org module, we identify trans-synaptic pacemaker cell interactions with glia. We introduce a druggable target prediction tool, drug2cell, which leverages single-cell profiles and drug-target interactions to provide mechanistic insights into the chronotropic effects of drugs, including GLP-1 analogues. In the epicardium, we show enrichment of both IgG+ and IgA+ plasma cells forming immune niches that may contribute to infection defence. Overall, we provide new clarity to cardiac electro-anatomy and immunology, and our suite of computational approaches can be applied to other tissues and organs.

FXR inhibition may protect from SARS-CoV-2 infection by reducing ACE2.

Preventing SARS-CoV-2 infection by modulating viral host receptors, such as angiotensin-converting enzyme 2 (ACE2)1, could represent a new chemoprophylactic approach for COVID-19 that complements vaccination2,3. However, the mechanisms that control the expression of ACE2 remain unclear. Here we show that the farnesoid X receptor (FXR) is a direct regulator of ACE2 transcription in several tissues affected by COVID-19, including the gastrointestinal and respiratory systems. We then use the over-the-counter compound z-guggulsterone and the off-patent drug ursodeoxycholic acid (UDCA) to reduce FXR signalling and downregulate ACE2 in human lung, cholangiocyte and intestinal organoids and in the corresponding tissues in mice and hamsters. We show that the UDCA-mediated downregulation of ACE2 reduces susceptibility to SARS-CoV-2 infection in vitro, in vivo and in human lungs and livers perfused ex situ. Furthermore, we reveal that UDCA reduces the expression of ACE2 in the nasal epithelium in humans. Finally, we identify a correlation between UDCA treatment and positive clinical outcomes after SARS-CoV-2 infection using retrospective registry data, and confirm these findings in an independent validation cohort of recipients of liver transplants. In conclusion, we show that FXR has a role in controlling ACE2 expression and provide evidence that modulation of this pathway could be beneficial for reducing SARS-CoV-2 infection, paving the way for future clinical trials.

Gene-level alignment of single-cell trajectories.

Single-cell data analysis can infer dynamic changes in cell populations, for example across time, space or in response to perturbation, thus deriving pseudotime trajectories. Current approaches comparing trajectories often use dynamic programming but are limited by assumptions such as the existence of a definitive match. Here we describe Genes2Genes, a Bayesian information-theoretic dynamic programming framework for aligning single-cell trajectories. It is able to capture sequential matches and mismatches of individual genes between a reference and query trajectory, highlighting distinct clusters of alignment patterns. Across both real world and simulated datasets, it accurately inferred alignments and demonstrated its utility in disease cell-state trajectory analysis. In a proof-of-concept application, Genes2Genes revealed that T cells differentiated in vitro match an immature in vivo state while lacking expression of genes associated with TNF signaling. This demonstrates that precise trajectory alignment can pinpoint divergence from the in vivo system, thus guiding the optimization of in vitro culture conditions.

Regional Differences in Human Biliary Tissues and Corresponding In Vitro-Derived Organoids.

BACKGROUND AND AIMS: Organoids provide a powerful system to study epithelia in vitro. Recently, this approach was applied successfully to the biliary tree, a series of ductular tissues responsible for the drainage of bile and pancreatic secretions. More precisely, organoids have been derived from ductal tissue located outside (extrahepatic bile ducts; EHBDs) or inside the liver (intrahepatic bile ducts; IHBDs). These organoids share many characteristics, including expression of cholangiocyte markers such as keratin (KRT) 19. However, the relationship between these organoids and their tissues of origin, and to each other, is largely unknown. APPROACH AND RESULTS: Organoids were derived from human gallbladder, common bile duct, pancreatic duct, and IHBDs using culture conditions promoting WNT signaling. The resulting IHBD and EHBD organoids expressed stem/progenitor markers leucine-rich repeat-containing G-protein-coupled receptor 5/prominin 1 and ductal markers KRT19/KRT7. However, RNA sequencing revealed that organoids conserve only a limited number of regional-specific markers corresponding to their location of origin. Of particular interest, down-regulation of biliary markers and up-regulation of cell-cycle genes were observed in organoids. IHBD and EHBD organoids diverged in their response to WNT signaling, and only IHBDs were able to express a low level of hepatocyte markers under differentiation conditions. CONCLUSIONS: Taken together, our results demonstrate that differences exist not only between extrahepatic biliary organoids and their tissue of origin, but also between IHBD and EHBD organoids. This information may help to understand the tissue specificity of cholangiopathies and also to identify targets for therapeutic development.

Snapin promotes HIV-1 transmission from dendritic cells by dampening TLR8 signaling.

HIV-1 traffics through dendritic cells (DCs) en route to establishing a productive infection in T lymphocytes but fails to induce an innate immune response. Within DC endosomes, HIV-1 somehow evades detection by the pattern-recognition receptor (PRR) Toll-like receptor 8 (TLR8). Using a phosphoproteomic approach, we identified a robust and diverse signaling cascade triggered by HIV-1 upon entry into human DCs. A secondary siRNA screen of the identified signaling factors revealed several new mediators of HIV-1 trans-infection of CD4+ T cells in DCs, including the dynein motor protein Snapin. Inhibition of Snapin enhanced localization of HIV-1 with TLR8+ early endosomes, triggered a pro-inflammatory response, and inhibited trans-infection of CD4+ T cells. Snapin inhibited TLR8 signaling in the absence of HIV-1 and is a general regulator of endosomal maturation. Thus, we identify a new mechanism of innate immune sensing by TLR8 in DCs, which is exploited by HIV-1 to promote transmission.

Cell proliferation within small intestinal crypts is the principal driving force for cell migration on villi.

The functional integrity of the intestinal epithelial barrier relies on tight coordination of cell proliferation and migration, with failure to regulate these processes resulting in disease. It is not known whether cell proliferation is sufficient to drive epithelial cell migration during homoeostatic turnover of the epithelium. Nor is it known precisely how villus cell migration is affected when proliferation is perturbed. Some reports suggest that proliferation and migration may not be related while other studies support a direct relationship. We used established cell-tracking methods based on thymine analog cell labeling and developed tailored mathematical models to quantify cell proliferation and migration under normal conditions and when proliferation is reduced and when it is temporarily halted. We found that epithelial cell migration velocities along the villi are coupled to cell proliferation rates within the crypts in all conditions. Furthermore, halting and resuming proliferation results in the synchronized response of cell migration on the villi. We conclude that cell proliferation within the crypt is the primary force that drives cell migration along the villus. This methodology can be applied to interrogate intestinal epithelial dynamics and characterize situations in which processes involved in cell turnover become uncoupled, including pharmacological treatments and disease models.-Parker, A., Maclaren, O. J., Fletcher, A. G., Muraro, D., Kreuzaler, P. A., Byrne, H. M., Maini, P. K., Watson, A. J. M., Pin, C. Cell proliferation within small intestinal crypts is the principal driving force for cell migration on villi.

Integration of hormonal signaling networks and mobile microRNAs is required for vascular patterning in Arabidopsis roots.

As multicellular organisms grow, positional information is continually needed to regulate the pattern in which cells are arranged. In the Arabidopsis root, most cell types are organized in a radially symmetric pattern; however, a symmetry-breaking event generates bisymmetric auxin and cytokinin signaling domains in the stele. Bidirectional cross-talk between the stele and the surrounding tissues involving a mobile transcription factor, SHORT ROOT (SHR), and mobile microRNA species also determines vascular pattern, but it is currently unclear how these signals integrate. We use a multicellular model to determine a minimal set of components necessary for maintaining a stable vascular pattern. Simulations perturbing the signaling network show that, in addition to the mutually inhibitory interaction between auxin and cytokinin, signaling through SHR, microRNA165/6, and PHABULOSA is required to maintain a stable bisymmetric pattern. We have verified this prediction by observing loss of bisymmetry in shr mutants. The model reveals the importance of several features of the network, namely the mutual degradation of microRNA165/6 and PHABULOSA and the existence of an additional negative regulator of cytokinin signaling. These components form a plausible mechanism capable of patterning vascular tissues in the absence of positional inputs provided by the transport of hormones from the shoot.

An integrative analysis of gene expression and molecular interaction data to identify dys-regulated sub-networks in inflammatory bowel disease.

BACKGROUND: Inflammatory bowel disease (IBD) consists of two main disease-subtypes, Crohn's disease (CD) and ulcerative colitis (UC); these subtypes share overlapping genetic and clinical features. Genome-wide microarray data enable unbiased documentation of alterations in gene expression that may be disease-specific. As genetic diseases are believed to be caused by genetic alterations affecting the function of signalling pathways, module-centric optimisation algorithms, whose aim is to identify sub-networks that are dys-regulated in disease, are emerging as promising approaches. RESULTS: In order to account for the topological structure of molecular interaction networks, we developed an optimisation algorithm that integrates databases of known molecular interactions with gene expression data; such integration enables identification of differentially regulated network modules. We verified the performance of our algorithm by testing it on simulated networks; we then applied the same method to study experimental data derived from microarray analysis of CD and UC biopsies and human interactome databases. This analysis allowed the extraction of dys-regulated subnetworks under different experimental conditions (inflamed and uninflamed tissues in CD and UC). Optimisation was performed to highlight differentially expressed network modules that may be common or specific to the disease subtype. CONCLUSIONS: We show that the selected subnetworks include genes and pathways of known relevance for IBD; in particular, the solutions found highlight cross-talk among enriched pathways, mainly the JAK/STAT signalling pathway and the EGF receptor signalling pathway. In addition, integration of gene expression with molecular interaction data highlights nodes that, although not being differentially expressed, interact with differentially expressed nodes and are part of pathways that are relevant to IBD. The method proposed here may help identifying dys-regulated sub-networks that are common in different diseases and sub-networks whose dys-regulation is specific to a particular disease.

The role of auxin and cytokinin signalling in specifying the root architecture of Arabidopsis thaliana.

Auxin and cytokinin are key hormonal signals that control the cellular architecture of the primary root and the initiation of new lateral root organs in the plant Arabidopsis thaliana. Both developmental processes are regulated by cross-talk between these hormones and their signalling pathways. In this paper, sub-cellular and multi-cellular mathematical models are developed to investigate how interactions between auxin and cytokinin influence the size and location of regions of division and differentiation within the primary root, and describe how their cross-regulation may cause periodic branching of lateral roots. We show how their joint activity may influence tissue-specific oscillations in gene expression, as shown in Moreno-Risueno et al. (2010) and commented upon in Traas and Vernoux (2010), and we propose mechanisms that may generate synchronisation of such periodic behaviours inside a cell and with its neighbours. Using a multi-cellular model, we also analyse the roles of cytokinin and auxin in specifying the three main regions of the primary root (elongation, transition and division zones), our simulation results being in good agreement with independent experimental observations. We then use our model to generate testable predictions concerning the effect of varying the concentrations of the auxin efflux transporters on the sizes of the different root regions. In particular, we predict that over-expression of the transporters will generate a longer root with a longer elongation zone and a smaller division zone than that of a wild type root. This root will contain fewer cells than its wild type counterpart. We conclude that our model can provide a useful tool for investigating the response of cell division and elongation to perturbations in hormonal signalling.

Single-cell atlas of human liver development reveals pathways directing hepatic cell fates.

The liver has been studied extensively due to the broad number of diseases affecting its vital functions. However, therapeutic advances have been hampered by the lack of knowledge concerning human hepatic development. Here, we addressed this limitation by describing the developmental trajectories of different cell types that make up the human liver at single-cell resolution. These transcriptomic analyses revealed that sequential cell-to-cell interactions direct functional maturation of hepatocytes, with non-parenchymal cells playing essential roles during organogenesis. We utilized this information to derive bipotential hepatoblast organoids and then exploited this model system to validate the importance of signalling pathways in hepatocyte and cholangiocyte specification. Further insights into hepatic maturation also enabled the identification of stage-specific transcription factors to improve the functionality of hepatocyte-like cells generated from human pluripotent stem cells. Thus, our study establishes a platform to investigate the basic mechanisms directing human liver development and to produce cell types for clinical applications.