Cell-type-specific and disease-associated expression quantitative trait loci in the human lung.

Common genetic variants confer substantial risk for chronic lung diseases, including pulmonary fibrosis. Defining the genetic control of gene expression in a cell-type-specific and context-dependent manner is critical for understanding the mechanisms through which genetic variation influences complex traits and disease pathobiology. To this end, we performed single-cell RNA sequencing of lung tissue from 66 individuals with pulmonary fibrosis and 48 unaffected donors. Using a pseudobulk approach, we mapped expression quantitative trait loci (eQTLs) across 38 cell types, observing both shared and cell-type-specific regulatory effects. Furthermore, we identified disease interaction eQTLs and demonstrated that this class of associations is more likely to be cell-type-specific and linked to cellular dysregulation in pulmonary fibrosis. Finally, we connected lung disease risk variants to their regulatory targets in disease-relevant cell types. These results indicate that cellular context determines the impact of genetic variation on gene expression and implicates context-specific eQTLs as key regulators of lung homeostasis and disease.

Cell type-specific and disease-associated eQTL in the human lung.

Common genetic variants confer substantial risk for chronic lung diseases, including pulmonary fibrosis (PF). Defining the genetic control of gene expression in a cell-type-specific and context-dependent manner is critical for understanding the mechanisms through which genetic variation influences complex traits and disease pathobiology. To this end, we performed single-cell RNA-sequencing of lung tissue from 67 PF and 49 unaffected donors. Employing a pseudo-bulk approach, we mapped expression quantitative trait loci (eQTL) across 38 cell types, observing both shared and cell type-specific regulatory effects. Further, we identified disease-interaction eQTL and demonstrated that this class of associations is more likely to be cell-type specific and linked to cellular dysregulation in PF. Finally, we connected PF risk variants to their regulatory targets in disease-relevant cell types. These results indicate that cellular context determines the impact of genetic variation on gene expression, and implicates context-specific eQTL as key regulators of lung homeostasis and disease.

Neutrophil trafficking to the site of infection requires Cpt1a-dependent fatty acid ß-oxidation.

Cellular metabolism influences immune cell function, with mitochondrial fatty acid ß-oxidation and oxidative phosphorylation required for multiple immune cell phenotypes. Carnitine palmitoyltransferase 1a (Cpt1a) is considered the rate-limiting enzyme for mitochondrial metabolism of long-chain fatty acids, and Cpt1a deficiency is associated with infant mortality and infection risk. This study was undertaken to test the hypothesis that impairment in Cpt1a-dependent fatty acid oxidation results in increased susceptibility to infection. Screening the Cpt1a gene for common variants predicted to affect protein function revealed allele rs2229738_T, which was associated with pneumonia risk in a targeted human phenome association study. Pharmacologic inhibition of Cpt1a increases mortality and impairs control of the infection in a murine model of bacterial pneumonia. Susceptibility to pneumonia is associated with blunted neutrophilic responses in mice and humans that result from impaired neutrophil trafficking to the site of infection. Chemotaxis responsible for neutrophil trafficking requires Cpt1a-dependent mitochondrial fatty acid oxidation for amplification of chemoattractant signals. These findings identify Cpt1a as a potential host determinant of infection susceptibility and demonstrate a requirement for mitochondrial fatty acid oxidation in neutrophil biology.

Increased Development of Th1, Th17, and Th1.17 Cells Under T1 Polarizing Conditions in Juvenile Idiopathic Arthritis.

In juvenile idiopathic arthritis (JIA) inflammatory T cells and their produced cytokines are drug targets and play a role in disease pathogenesis. Despite their clinical importance, the sources and types of inflammatory T cells involved remain unclear. T cells respond to polarizing factors to initiate types of immunity to fight infections, which include immunity types 1 (T1), 2 (T2), and 3 (T17). Polarizing factors drive CD4+ T cells towards T helper (Th) cell subtypes and CD8+ T cells towards cytotoxic T cell (Tc) subtypes. T1 and T17 polarization are associated with autoimmunity and production of the cytokines IFNγ and IL-17 respectively. We show that JIA and child healthy control (HC) peripheral blood mononuclear cells are remarkably similar, with the same frequencies of CD4+ and CD8+ naïve and memory T cell subsets, T cell proliferation, and CD4+ and CD8+ T cell subsets upon T1, T2, and T17 polarization. Yet, under T1 polarizing conditions JIA cells produced increased IFNγ and inappropriately produced IL-17. Under T17 polarizing conditions JIA T cells produced increased IL-17. Gene expression of IFNγ, IL-17, Tbet, and RORγT by quantitative PCR and RNA sequencing revealed activation of immune responses and inappropriate activation of IL-17 signaling pathways in JIA polarized T1 cells. The polarized JIA T1 cells were comprised of Th and Tc cells, with Th cells producing IFNγ (Th1), IL-17 (Th17), and both IFNγ-IL-17 (Th1.17) and Tc cells producing IFNγ (Tc1). The JIA polarized CD4+ T1 cells expressed both Tbet and RORγT, with higher expression of the transcription factors associated with higher frequency of IL-17 producing cells. T1 polarized naïve CD4+ cells from JIA also produced more IFNγ and more IL-17 than HC. We show that in JIA T1 polarization inappropriately generates Th1, Th17, and Th1.17 cells. Our data provides a tool for studying the development of heterogeneous inflammatory T cells in JIA under T1 polarizing conditions and for identifying pathogenic immune cells that are important as drug targets and diagnostic markers.

Synergistic action of WDR5 and HDM2 inhibitors in SMARCB1-deficient cancer cells.

Rhabdoid tumors (RT) are rare and deadly pediatric cancers driven by loss of SMARCB1, which encodes the SNF5 component of the SWI/SNF chromatin remodeler. Loss of SMARCB1 is associated with a complex set of phenotypic changes including vulnerability to inhibitors of protein synthesis and of the p53 ubiquitin-ligase HDM2. Recently, we discovered small molecule inhibitors of the 'WIN' site of WDR5, which in MLL-rearranged leukemia cells decrease the expression of a set of genes linked to protein synthesis, inducing a translational choke and causing p53-dependent inhibition of proliferation. Here, we characterize how WIN site inhibitors act in RT cells. As in leukemia cells, WIN site inhibition in RT cells causes the comprehensive displacement of WDR5 from chromatin, resulting in a decrease in protein synthesis gene expression. Unlike leukemia cells, however, the growth response of RT cells to WIN site blockade is independent of p53. Exploiting this observation, we demonstrate that WIN site inhibitor synergizes with an HDM2 antagonist to induce p53 and block RT cell proliferation in vitro. These data reveal a p53-independent action of WIN site inhibitors and forecast that future strategies to treat RT could be based on dual WDR5/HDM2 inhibition.

WIN site inhibition disrupts a subset of WDR5 function.

WDR5 nucleates the assembly of histone-modifying complexes and acts outside this context in a range of chromatin-centric processes. WDR5 is also a prominent target for pharmacological inhibition in cancer. Small-molecule degraders of WDR5 have been described, but most drug discovery efforts center on blocking the WIN site of WDR5, an arginine binding cavity that engages MLL/SET enzymes that deposit histone H3 lysine 4 methylation (H3K4me). Therapeutic application of WIN site inhibitors is complicated by the disparate functions of WDR5, but is generally guided by two assumptions-that WIN site inhibitors disable all functions of WDR5, and that changes in H3K4me drive the transcriptional response of cancer cells to WIN site blockade. Here, we test these assumptions by comparing the impact of WIN site inhibition versus WDR5 degradation on H3K4me and transcriptional processes. We show that WIN site inhibition disables only a specific subset of WDR5 activity, and that H3K4me changes induced by WDR5 depletion do not explain accompanying transcriptional responses. These data recast WIN site inhibitors as selective loss-of-function agents, contradict H3K4me as a relevant mechanism of action for WDR5 inhibitors, and indicate distinct clinical applications of WIN site inhibitors and WDR5 degraders.

Multiple interactions of the oncoprotein transcription factor MYC with the SWI/SNF chromatin remodeler.

The SNF5 subunit of the SWI/SNF chromatin remodeling complex has been shown to act as a tumor suppressor through multiple mechanisms, including impairing the ability of the oncoprotein transcription factor MYC to bind chromatin. Beyond SNF5, however, it is unknown to what extent MYC can access additional SWI/SNF subunits or how these interactions affect the ability of MYC to drive transcription, particularly in SNF5-null cancers. Here, we report that MYC interacts with multiple SWI/SNF components independent of SNF5. We show that MYC binds the pan-SWI/SNF subunit BAF155 through the BAF155 SWIRM domain, an interaction that is inhibited by the presence of SNF5. In SNF5-null cells, MYC binds with remaining SWI/SNF components to essential genes, although for a purpose that is distinct from chromatin remodeling. Analysis of MYC-SWI/SNF target genes in SNF5-null cells reveals that they are associated with core biological functions of MYC linked to protein synthesis. These data reveal that MYC can bind SWI/SNF in an SNF5-independent manner and that SNF5 modulates access of MYC to core SWI/SNF complexes. This work provides a framework in which to interrogate the influence of SWI/SNF on MYC function in cancers in which SWI/SNF or MYC are altered.

Coordinated interactions between endothelial cells and macrophages in the islet microenvironment promote ß cell regeneration.

Endogenous ß cell regeneration could alleviate diabetes, but proliferative stimuli within the islet microenvironment are incompletely understood. We previously found that ß cell recovery following hypervascularization-induced ß cell loss involves interactions with endothelial cells (ECs) and macrophages (MΦs). Here we show that proliferative ECs modulate MΦ infiltration and phenotype during ß cell loss, and recruited MΦs are essential for ß cell recovery. Furthermore, VEGFR2 inactivation in quiescent ECs accelerates islet vascular regression during ß cell recovery and leads to increased ß cell proliferation without changes in MΦ phenotype or number. Transcriptome analysis of ß cells, ECs, and MΦs reveals that ß cell proliferation coincides with elevated expression of extracellular matrix remodeling molecules and growth factors likely driving activation of proliferative signaling pathways in ß cells. Collectively, these findings suggest a new ß cell regeneration paradigm whereby coordinated interactions between intra-islet MΦs, ECs, and extracellular matrix mediate ß cell self-renewal.

Combinatorial Transcriptional Profiling of Mouse and Human Enteric Neurons Identifies Shared and Disparate Subtypes In Situ.

BACKGROUND & AIMS: The enteric nervous system (ENS) coordinates essential intestinal functions through the concerted action of diverse enteric neurons (ENs). However, integrated molecular knowledge of EN subtypes is lacking. To compare human and mouse ENs, we transcriptionally profiled healthy ENS from adult humans and mice. We aimed to identify transcripts marking discrete neuron subtypes and visualize conserved EN subtypes for humans and mice in multiple bowel regions. METHODS: Human myenteric ganglia and adjacent smooth muscle were isolated by laser-capture microdissection for RNA sequencing. Ganglia-specific transcriptional profiles were identified by computationally subtracting muscle gene signatures. Nuclei from mouse myenteric neurons were isolated and subjected to single-nucleus RNA sequencing, totaling more than 4 billion reads and 25,208 neurons. Neuronal subtypes were defined using mouse single-nucleus RNA sequencing data. Comparative informatics between human and mouse data sets identified shared EN subtype markers, which were visualized in situ using hybridization chain reaction. RESULTS: Several EN subtypes in the duodenum, ileum, and colon are conserved between humans and mice based on orthologous gene expression. However, some EN subtype-specific genes from mice are expressed in completely distinct morphologically defined subtypes in humans. In mice, we identified several neuronal subtypes that stably express gene modules across all intestinal segments, with graded, regional expression of 1 or more marker genes. CONCLUSIONS: Our combined transcriptional profiling of human myenteric ganglia and mouse EN provides a rich foundation for developing novel intestinal therapeutics. There is congruency among some EN subtypes, but we note multiple species differences that should be carefully considered when relating findings from mouse ENS research to human gastrointestinal studies.

WDR5 is a conserved regulator of protein synthesis gene expression.

WDR5 is a highly-conserved nuclear protein that performs multiple scaffolding functions in the context of chromatin. WDR5 is also a promising target for pharmacological inhibition in cancer, with small molecule inhibitors of an arginine-binding pocket of WDR5 (the 'WIN' site) showing efficacy against a range of cancer cell lines in vitro. Efforts to understand WDR5, or establish the mechanism of action of WIN site inhibitors, however, are stymied by its many functions in the nucleus, and a lack of knowledge of the conserved gene networks-if any-that are under its control. Here, we have performed comparative genomic analyses to identify the conserved sites of WDR5 binding to chromatin, and the conserved genes regulated by WDR5, across a diverse panel of cancer cell lines. We show that a specific cohort of protein synthesis genes (PSGs) are invariantly bound by WDR5, demonstrate that the WIN site anchors WDR5 to chromatin at these sites, and establish that PSGs are bona fide, acute, and persistent targets of WIN site blockade. Together, these data reveal that WDR5 plays a predominant transcriptional role in biomass accumulation and provide further evidence that WIN site inhibitors act to repress gene networks linked to protein synthesis homeostasis.

The LRRC8 volume-regulated anion channel inhibitor, DCPIB, inhibits mitochondrial respiration independently of the channel.

There has been a resurgence of interest in the volume-regulated anion channel (VRAC) since the recent cloning of the LRRC8A-E gene family that encodes VRAC. The channel is a heteromer comprised of LRRC8A and at least one other family member; disruption of LRRC8A expression abolishes VRAC activity. The best-in-class VRAC inhibitor, DCPIB, suffers from off-target activity toward several different channels and transporters. Considering that some anion channel inhibitors also suppress mitochondrial respiration, we systematically explored whether DCPIB inhibits respiration in wild type (WT) and LRRC8A-knockout HAP-1 and HEK-293 cells. Knockout of LRRC8A had no apparent effects on cell morphology, proliferation rate, mitochondrial content, or expression of several mitochondrial genes in HAP-1 cells. Addition of 10 µM DCPIB, a concentration typically used to inhibit VRAC, suppressed basal and ATP-linked respiration in part through uncoupling the inner mitochondrial membrane (IMM) proton gradient and membrane potential. Additionally, DCPIB inhibits the activity of complex I, II, and III of the electron transport chain (ETC). Surprisingly, the effects of DCPIB on mitochondrial function are also observed in HAP-1 and HEK-293 cells which lack LRRC8A expression. Finally, we demonstrate that DCPIB activates ATP-inhibitable potassium channels comprised of heterologously expressed Kir6.2 and SUR1 subunits. These data indicate that DCPIB suppresses mitochondrial respiration and ATP production by dissipating the mitochondrial membrane potential and inhibiting complexes I-III of the ETC. They further justify the need for the development of sharper pharmacological tools for evaluating the integrative physiology and therapeutic potential of VRAC in human diseases.

Novel kidney dissociation protocol and image-based flow cytometry facilitate improved analysis of injured proximal tubules.

Flow cytometry studies on injured kidney tubules are complicated by the low yield of nucleated single cells. Furthermore, cell-specific responses such as cell cycle dynamics in vivo have conventionally relied on indirect immunohistochemistry and proximal tubule markers that may be downregulated in injury. Here, we report a new tissue dissociation protocol for the kidney with an early fixation step that greatly enhances the yield of single cells. Genetic labeling of the proximal tubule with either mT/mG "tomato" or R26Fucci2aR (Fucci) cell cycle reporter mice allows us to follow proximal tubule-specific changes in cell cycle after renal injury. Image-based flow cytometry (FlowSight) enables gating of the cell cycle and concurrent visualization of the cells with bright field and fluorescence. We used the Fucci mouse in conjunction with FlowSight to identify a discrete polyploid population in proximal tubules after aristolochic acid injury. The tissue dissociation protocol in conjunction with genetic labeling and image-based flow cytometry is a tool that can improve our understanding of any discrete cell population after injury.

Ectonucleoside Triphosphate Diphosphohydrolase-3 Antibody Targets Adult Human Pancreatic ß Cells for In Vitro and In Vivo Analysis.

Identification of cell-surface markers specific to human pancreatic ß cells would allow in vivo analysis and imaging. Here we introduce a biomarker, ectonucleoside triphosphate diphosphohydrolase-3 (NTPDase3), that is expressed on the cell surface of essentially all adult human ß cells, including those from individuals with type 1 or type 2 diabetes. NTPDase3 is expressed dynamically during postnatal human pancreas development, appearing first in acinar cells at birth, but several months later its expression declines in acinar cells while concurrently emerging in islet ß cells. Given its specificity and membrane localization, we utilized an NTPDase3 antibody for purification of live human ß cells as confirmed by transcriptional profiling, and, in addition, for in vivo imaging of transplanted human ß cells. Thus, NTPDase3 is a cell-surface biomarker of adult human ß cells, and the antibody directed to this protein should be a useful new reagent for ß cell sorting, in vivo imaging, and targeting.

Alternative splicing of ALCAM enables tunable regulation of cell-cell adhesion through differential proteolysis.

While many adhesion receptors are known to influence tumor progression, the mechanisms by which they dynamically regulate cell-cell adhesion remain elusive. We previously identified Activated Leukocyte Cell Adhesion Molecule (ALCAM) as a clinically relevant driver of metastasis and hypothesized that a tunable mechanism of ectodomain shedding regulates its contribution to dissemination. To test this hypothesis, we examined an under-explored ALCAM splice variant (ALCAM-Iso2) and demonstrated that loss of the membrane-proximal region of ALCAM (exon 13) increased metastasis four-fold. Mechanistic studies identified a novel MMP14-dependent membrane distal cleavage site in ALCAM-Iso2, which mediated a ten-fold increase in shedding, thereby decreasing cellular cohesion. Importantly, the loss of cohesion is not limited to the cell capable of shedding because the released extracellular domain diminished cohesion of non-shedding cells through disruption of ALCAM-ALCAM interactions. ALCAM-Iso2-dominated expression in bladder cancer tissue, compared to normal bladder, further emphasizes that ALCAM alternative splicing may contribute to clinical disease progression. The requirement for both the loss of exon 13 and the gain of metalloprotease activity suggests that ALCAM shedding and concomitant regulation of tumor cell adhesion is a locally tunable process.

Obesity Alters B Cell and Macrophage Populations in Brown Adipose Tissue.

OBJECTIVE: The prevalence of obesity continues to rise, and it is understood that regulation of white adipose tissue (WAT) function is important to systemic metabolic homeostasis. Immune cells play a central role in the maintenance of WAT, and their compositions change in number and inflammatory phenotype with the progression of obesity. Because of its energy-burning capabilities, brown adipose tissue (BAT) has become a focus of obesity research. Although novel studies have focused on the function of brown adipocytes in thermogenesis, the tissue as a whole has not been immunologically characterized. METHODS: BAT immune cell populations were analyzed by flow cytometry and immunohistochemistry in mice with diet-induced obesity (3, 8, or 16 weeks of diet) and in aged mice (1, 6-7, and 10-15 months). RESULTS: The data confirmed the presence of macrophages and eosinophils, as previously reported, and showed that 20% to 30% of the immune cells in BAT were B cells. The number of B cells and eosinophils increased with diet-induced obesity, whereas macrophages decreased. There was no change in number of any immune cell quantified with age. CONCLUSIONS: These studies reveal a novel finding of B220 + B cells in BAT and show that BAT immune cell populations change in response to diet-induced obesity.

Chronic HIV-1 Infection Impairs Superantigen-Induced Activation of Peripheral CD4+CXCR5+PD-1+ Cells, With Relative Preservation of Recall Antigen-Specific Responses.

Peripheral CD4+CXCR5+PD-1+ T cells are a putative circulating counterpart to germinal center T follicular helper (TFH) cells. They show both phenotypic and functional similarities to TFH cells, which provide necessary help for the differentiation of B cells to antibody-secreting plasmablasts. In this study, we evaluated the frequency, phenotypes, and responses of peripheral TFH-like (pTFH) cells to superantigen and recall antigen stimulation in 10 healthy and 34 chronically infected treatment-naive HIV-1+ individuals. There was no difference in the frequency of pTFH cells between HIV+ and HIV- individuals. Surface expression of ICOS, but not CD40L, was higher on pTFH cells at baseline in HIV+ individuals. Compared with HIV- individuals, pTFH cells from HIV+ individuals had decreased maximal expression of ICOS and CD40L in response to in vitro superantigen stimulation. This decreased response did not correlate with viral control, CD4 T-cell count, duration of infection, or the degree of neutralizing antibody breadth. Despite a decreased maximal response, pTFH responses to HIV Gag and tetanus toxoid recall antigens were preserved.

Distributions of Cells and Neurons across the Cortical Sheet in Old World Macaques.

According to previous research, cell and neuron densities vary across neocortex in a similar manner across primate taxa. Here, we provide a more extensive examination of this effect in macaque monkeys. We separated neocortex from the underlying white matter in 4 macaque monkey hemispheres (1 Macaca nemestrina, 2 Macaca radiata, and 1 Macaca mulatta), manually flattened the neocortex, and divided it into smaller tissue pieces for analysis. The number of cells and neurons were determined for each piece across the cortical sheet using flow cytometry. Primary visual cortex had the most densely packed neurons and primary motor cortex had the least densely packed neurons. With respect to differences in brain size between cases, there was little variability in the total cell and neuron numbers within specific areas, and overall trends were similar to what has been previously described in Old World baboons and other primates. The average hemispheric total cell number per hemisphere ranged from 2.9 to 3.7 billion, while the average total neuron number ranged from 1.3 to 1.7 billion neurons. The visual cortex neuron densities were predictably higher, ranging from 18.2 to 34.7 million neurons/cm2 in macaques, in comparison to a range of 9.3-17.7 million neurons/cm2 across cortex as a whole. The results support other evidence that neuron surface densities vary across the cortical sheet in a predictable pattern within and across primate taxa.

Cortical cell and neuron density estimates in one chimpanzee hemisphere.

The density of cells and neurons in the neocortex of many mammals varies across cortical areas and regions. This variability is, perhaps, most pronounced in primates. Nonuniformity in the composition of cortex suggests regions of the cortex have different specializations. Specifically, regions with densely packed neurons contain smaller neurons that are activated by relatively few inputs, thereby preserving information, whereas regions that are less densely packed have larger neurons that have more integrative functions. Here we present the numbers of cells and neurons for 742 discrete locations across the neocortex in a chimpanzee. Using isotropic fractionation and flow fractionation methods for cell and neuron counts, we estimate that neocortex of one hemisphere contains 9.5 billion cells and 3.7 billion neurons. Primary visual cortex occupies 35 cm(2) of surface, 10% of the total, and contains 737 million densely packed neurons, 20% of the total neurons contained within the hemisphere. Other areas of high neuron packing include secondary visual areas, somatosensory cortex, and prefrontal granular cortex. Areas of low levels of neuron packing density include motor and premotor cortex. These values reflect those obtained from more limited samples of cortex in humans and other primates.

Multiparameter analysis of stimulated human peripheral blood mononuclear cells: A comparison of mass and fluorescence cytometry.

Mass and fluorescence cytometry are quantitative single cell flow cytometry approaches that are powerful tools for characterizing diverse tissues and cellular systems. Here mass cytometry was directly compared with fluorescence cytometry by studying phenotypes of healthy human peripheral blood mononuclear cells (PBMC) in the context of superantigen stimulation. One mass cytometry panel and five fluorescence cytometry panels were used to measure 20 well-established lymphocyte markers of memory and activation. Comparable frequencies of both common and rare cell subpopulations were observed with fluorescence and mass cytometry using biaxial gating. The unsupervised high-dimensional analysis tool viSNE was then used to analyze data sets generated from both mass and fluorescence cytometry. viSNE analysis effectively characterized PBMC using eight features per cell and identified similar frequencies of activated CD4+ T cells with both technologies. These results suggest combinations of unsupervised analysis programs and extended multiparameter cytometry will be indispensable tools for detecting perturbations in protein expression in both health and disease.