Organoid Modeling of the Tumor Immune Microenvironment.

In vitro cancer cultures, including three-dimensional organoids, typically contain exclusively neoplastic epithelium but require artificial reconstitution to recapitulate the tumor microenvironment (TME). The co-culture of primary tumor epithelia with endogenous, syngeneic tumor-infiltrating lymphocytes (TILs) as a cohesive unit has been particularly elusive. Here, an air-liquid interface (ALI) method propagated patient-derived organoids (PDOs) from >100 human biopsies or mouse tumors in syngeneic immunocompetent hosts as tumor epithelia with native embedded immune cells (T, B, NK, macrophages). Robust droplet-based, single-cell simultaneous determination of gene expression and immune repertoire indicated that PDO TILs accurately preserved the original tumor T cell receptor (TCR) spectrum. Crucially, human and murine PDOs successfully modeled immune checkpoint blockade (ICB) with anti-PD-1- and/or anti-PD-L1 expanding and activating tumor antigen-specific TILs and eliciting tumor cytotoxicity. Organoid-based propagation of primary tumor epithelium en bloc with endogenous immune stroma should enable immuno-oncology investigations within the TME and facilitate personalized immunotherapy testing.

Progenitor identification and SARS-CoV-2 infection in human distal lung organoids.

The distal lung contains terminal bronchioles and alveoli that facilitate gas exchange. Three-dimensional in vitro human distal lung culture systems would strongly facilitate the investigation of pathologies such as interstitial lung disease, cancer and coronavirus disease 2019 (COVID-19) pneumonia caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Here we describe the development of a long-term feeder-free, chemically defined culture system for distal lung progenitors as organoids derived from single adult human alveolar epithelial type II (AT2) or KRT5+ basal cells. AT2 organoids were able to differentiate into AT1 cells, and basal cell organoids developed lumens lined with differentiated club and ciliated cells. Single-cell analysis of KRT5+ cells in basal organoids revealed a distinct population of ITGA6+ITGB4+ mitotic cells, whose offspring further segregated into a TNFRSF12Ahi subfraction that comprised about ten per cent of KRT5+ basal cells. This subpopulation formed clusters within terminal bronchioles and exhibited enriched clonogenic organoid growth activity. We created distal lung organoids with apical-out polarity to present ACE2 on the exposed external surface, facilitating infection of AT2 and basal cultures with SARS-CoV-2 and identifying club cells as a target population. This long-term, feeder-free culture of human distal lung organoids, coupled with single-cell analysis, identifies functional heterogeneity among basal cells and establishes a facile in vitro organoid model of human distal lung infections, including COVID-19-associated pneumonia.

Non-equivalence of Wnt and R-spondin ligands during Lgr5+ intestinal stem-cell self-renewal.

The canonical Wnt/ß-catenin signalling pathway governs diverse developmental, homeostatic and pathological processes. Palmitoylated Wnt ligands engage cell-surface frizzled (FZD) receptors and LRP5 and LRP6 co-receptors, enabling ß-catenin nuclear translocation and TCF/LEF-dependent gene transactivation. Mutations in Wnt downstream signalling components have revealed diverse functions thought to be carried out by Wnt ligands themselves. However, redundancy between the 19 mammalian Wnt proteins and 10 FZD receptors and Wnt hydrophobicity have made it difficult to attribute these functions directly to Wnt ligands. For example, individual mutations in Wnt ligands have not revealed homeostatic phenotypes in the intestinal epithelium-an archetypal canonical, Wnt pathway-dependent, rapidly self-renewing tissue, the regeneration of which is fueled by proliferative crypt Lgr5+ intestinal stem cells (ISCs). R-spondin ligands (RSPO1-RSPO4) engage distinct LGR4-LGR6, RNF43 and ZNRF3 receptor classes, markedly potentiate canonical Wnt/ß-catenin signalling, and induce intestinal organoid growth in vitro and Lgr5+ ISCs in vivo. However, the interchangeability, functional cooperation and relative contributions of Wnt versus RSPO ligands to in vivo canonical Wnt signalling and ISC biology remain unknown. Here we identify the functional roles of Wnt and RSPO ligands in the intestinal crypt stem-cell niche. We show that the default fate of Lgr5+ ISCs is to differentiate, unless both RSPO and Wnt ligands are present. However, gain-of-function studies using RSPO ligands and a new non-lipidated Wnt analogue reveal that these ligands have qualitatively distinct, non-interchangeable roles in ISCs. Wnt proteins are unable to induce Lgr5+ ISC self-renewal, but instead confer a basal competency by maintaining RSPO receptor expression that enables RSPO ligands to actively drive and specify the extent of stem-cell expansion. This functionally non-equivalent yet cooperative interaction between Wnt and RSPO ligands establishes a molecular precedent for regulation of mammalian stem cells by distinct priming and self-renewal factors, with broad implications for precise control of tissue regeneration.

Single-cell RNA-Seq of follicular lymphoma reveals malignant B-cell types and coexpression of T-cell immune checkpoints.

Follicular lymphoma (FL) is a low-grade B-cell malignancy that transforms into a highly aggressive and lethal disease at a rate of 2% per year. Perfect isolation of the malignant B-cell population from a surgical biopsy is a significant challenge, masking important FL biology, such as immune checkpoint coexpression patterns. To resolve the underlying transcriptional networks of follicular B-cell lymphomas, we analyzed the transcriptomes of 34 188 cells derived from 6 primary FL tumors. For each tumor, we identified normal immune subpopulations and malignant B cells, based on gene expression. We used multicolor flow cytometry analysis of the same tumors to confirm our assignments of cellular lineages and validate our predictions of expressed proteins. Comparison of gene expression between matched malignant and normal B cells from the same patient revealed tumor-specific features. Malignant B cells exhibited restricted immunoglobulin (Ig) light chain expression (either Igκ or Igλ), as well the expected upregulation of the BCL2 gene, but also downregulation of the FCER2, CD52, and major histocompatibility complex class II genes. By analyzing thousands of individual cells per patient tumor, we identified the mosaic of malignant B-cell subclones that coexist within a FL and examined the characteristics of tumor-infiltrating T cells. We identified genes coexpressed with immune checkpoint molecules, such as CEBPA and B2M in regulatory T (Treg) cells, providing a better understanding of the gene networks involved in immune regulation. In summary, parallel measurement of single-cell expression in thousands of tumor cells and tumor-infiltrating lymphocytes can be used to obtain a systems-level view of the tumor microenvironment and identify new avenues for therapeutic development.

Registered Dietitian Staffing in Dialysis: Impact on Quality Ratings in Mandate and Nonmandate States.

OBJECTIVE: Texas is the only state to mandate a patient to full-time equivalent (FTE) registered dietitian (RD) ratio (<125 patients:FTE RD) in dialysis. Little research exists about the relationship between patient:FTE RD ratio and quality. Our objectives were to a) examine the effect of the mandate on patient:FTE RD ratios in Texas facilities compared to similar nonmandated facilities and b) examine the association between patient:FTE RD ratio and the patient rating of staff quality or star rating of the facility. DESIGN AND METHODS: Using data from Dialysis Facility Annual Reports, we excluded facilities with less than 125 patients and matched based on region, chain/profit status, and number of patients, pairing each Texas facility with 2 non-Texas facilities. T-tests for difference of means and chi-square tests were performed to compare facility groups characteristics and assess mandate impact. We used correlation between patient:FTE RD ratio and staff quality linearized score, and ANOVA to compare patient:FTE RD ratios between the star rating levels. RESULTS: The patient:FTE RD ratio was higher in non-Texas than Texas facilities (111.84 ± 40.70 vs 90.80 ± 24.02, P < .01). The Pearson correlation between patient:FTE RD ratio and linearized score rating of the dialysis center staff was essentially nonexistent (r = -0.046). We found a trend of higher patient:FTE RD ratios in facilities with lower star ratings for quality and staff but no statistical significance. CONCLUSION: The mandate effectively lowered the patient:FTE RD ratio in Texas dialysis facilities with over 125 patients, indicating that states considering such a mandate may benefit from implementing it. We found no association between patient:FTE RD staffing and quality, either objectively or via patient assessment. Further research should examine whether dialysis facilities are trending toward smaller patient sizes to circumvent staffing mandates and examine the relationship between other professions' staffing ratios and quality.

Promoters recognized by forkhead proteins exist for individual 21U-RNAs.

C. elegans 21U-RNAs are equivalent to the piRNAs discovered in other metazoans and have important roles in gametogenesis and transposon control. The biogenesis and molecular function of 21U-RNAs and piRNAs are poorly understood. Here, we demonstrate that transcription of each 21U-RNA is regulated separately through a conserved upstream DNA motif. We use genomic analysis to show that this motif is associated with low nucleosome occupancy, a characteristic of many promoters that drive expression of protein-coding genes, and that RNA polymerase II is localized to this nucleosome-depleted region. We establish that the most conserved 8-mer sequence in the upstream region of 21U-RNAs, CTGTTTCA, is absolutely required for their individual expression. Furthermore, we demonstrate that the 8-mer is specifically recognized by Forkhead family (FKH) transcription factors and that 21U-RNA expression is diminished in several FKH mutants. Our results suggest that thousands of small noncoding transcription units are regulated by FKH proteins.

Suppression of progenitor differentiation requires the long noncoding RNA ANCR.

Long noncoding RNAs (lncRNAs) regulate diverse processes, yet a potential role for lncRNAs in maintaining the undifferentiated state in somatic tissue progenitor cells remains uncharacterized. We used transcriptome sequencing and tiling arrays to compare lncRNA expression in epidermal progenitor populations versus differentiating cells. We identified ANCR (anti-differentiation ncRNA) as an 855-base-pair lncRNA down-regulated during differentiation. Depleting ANCR in progenitor-containing populations, without any other stimuli, led to rapid differentiation gene induction. In epidermis, ANCR loss abolished the normal exclusion of differentiation from the progenitor-containing compartment. The ANCR lncRNA is thus required to enforce the undifferentiated cell state within epidermis.

Cryopreservation timing is a critical process parameter in a thymic regulatory T-cell therapy manufacturing protocol.

Regulatory T cells (Tregs) are a promising therapy for several immune-mediated conditions but manufacturing a homogeneous and consistent product, especially one that includes cryopreservation, has been challenging. Discarded pediatric thymuses are an excellent source of therapeutic Tregs with advantages including cell quantity, homogeneity and stability. Here we report systematic testing of activation reagents, cell culture media, restimulation timing and cryopreservation to develop a Good Manufacturing Practice (GMP)-compatible method to expand and cryopreserve Tregs. By comparing activation reagents, including soluble antibody tetramers, antibody-conjugated beads and artificial antigen-presenting cells (aAPCs) and different media, we found that the combination of Dynabeads Treg Xpander and ImmunoCult-XF medium preserved FOXP3 expression and suppressive function and resulted in expansion that was comparable with a single stimulation with aAPCs. Cryopreservation tests revealed a critical timing effect: only cells cryopreserved 1-3 days, but not >3 days, after restimulation maintained high viability and FOXP3 expression upon thawing. Restimulation timing was a less critical process parameter than the time between restimulation and cryopreservation. This systematic testing of key variables provides increased certainty regarding methods for in vitro expansion and cryopreservation of Tregs. The ability to cryopreserve expanded Tregs will have broad-ranging applications including enabling centralized manufacturing and long-term storage of cell products.

Control of somatic tissue differentiation by the long non-coding RNA TINCR.

Several of the thousands of human long non-coding RNAs (lncRNAs) have been functionally characterized; however, potential roles for lncRNAs in somatic tissue differentiation remain poorly understood. Here we show that a 3.7-kilobase lncRNA, terminal differentiation-induced ncRNA (TINCR), controls human epidermal differentiation by a post-transcriptional mechanism. TINCR is required for high messenger RNA abundance of key differentiation genes, many of which are mutated in human skin diseases, including FLG, LOR, ALOXE3, ALOX12B, ABCA12, CASP14 and ELOVL3. TINCR-deficient epidermis lacked terminal differentiation ultrastructure, including keratohyalin granules and intact lamellar bodies. Genome-scale RNA interactome analysis revealed that TINCR interacts with a range of differentiation mRNAs. TINCR-mRNA interaction occurs through a 25-nucleotide 'TINCR box' motif that is strongly enriched in interacting mRNAs and required for TINCR binding. A high-throughput screen to analyse TINCR binding capacity to approximately 9,400 human recombinant proteins revealed direct binding of TINCR RNA to the staufen1 (STAU1) protein. STAU1-deficient tissue recapitulated the impaired differentiation seen with TINCR depletion. Loss of UPF1 and UPF2, both of which are required for STAU1-mediated RNA decay, however, did not have differentiation effects. Instead, the TINCR-STAU1 complex seems to mediate stabilization of differentiation mRNAs, such as KRT80. These data identify TINCR as a key lncRNA required for somatic tissue differentiation, which occurs through lncRNA binding to differentiation mRNAs to ensure their expression.

Genome-wide impact of a recently expanded microRNA cluster in mouse.

Variations in microRNA (miRNA) gene and/or target repertoire are likely to be key drivers of phenotypic differences between species. To better understand these changes, we developed a computational method that identifies signatures of species-specific target site gain and loss associated with miRNA acquisition. Interestingly, several of the miRNAs implicated in mouse 3' UTR evolution derive from a single rapidly expanded rodent-specific miRNA cluster. Located in the intron of Sfmbt2, a maternally imprinted polycomb gene, these miRNAs (referred to as the Sfmbt2 cluster) are expressed in both embryonic stem cells and the placenta. One abundant miRNA from the cluster, miR-467a, functionally overlaps with the mir-290-295 cluster in promoting growth and survival of mouse embryonic stem cells. Predicted novel targets of the remaining cluster members are enriched in pathways regulating cell survival. Two relevant species-specific target candidates, Lats2 and Dedd2, were validated in cultured cells. We suggest that the rapid evolution of the Sfmbt2 cluster may be a result of intersex conflict for growth regulation in early mammalian development and could provide a general model for the genomic response to acquisition of miRNAs and similar regulatory factors.

A latent pro-survival function for the mir-290-295 cluster in mouse embryonic stem cells.

MicroRNAs (miRNAs) post-transcriptionally regulate the expression of thousands of distinct mRNAs. While some regulatory interactions help to maintain basal cellular functions, others are likely relevant in more specific settings, such as response to stress. Here we describe such a role for the mir-290-295 cluster, the dominant miRNA cluster in mouse embryonic stem cells (mESCs). Examination of a target list generated from bioinformatic prediction, as well as expression data following miRNA loss, revealed strong enrichment for apoptotic regulators, two of which we validated directly: Caspase 2, the most highly conserved mammalian caspase, and Ei24, a p53 transcriptional target. Consistent with these predictions, mESCs lacking miRNAs were more likely to initiate apoptosis following genotoxic exposure to gamma irradiation or doxorubicin. Knockdown of either candidate partially rescued this pro-apoptotic phenotype, as did transfection of members of the mir-290-295 cluster. These findings were recapitulated in a specific mir-290-295 deletion line, confirming that they reflect miRNA functions at physiological levels. In contrast to the basal regulatory roles previously identified, the pro-survival phenotype shown here may be most relevant to stressful gestations, where pro-oxidant metabolic states induce DNA damage. Similarly, this cluster may mediate chemotherapeutic resistance in a neoplastic context, making it a useful clinical target.

Improving Reliability of Immunological Assays by Defining Minimal Criteria for Cell Fitness.

Human PBMC-based assays are often used as biomarkers for the diagnosis and prognosis of disease, as well as for the prediction and tracking of response to biological therapeutics. However, the development and use of PBMC-based biomarker assays is often limited by poor reproducibility. Complex immunological assays can be further complicated by variation in cell handling before analysis, especially when using cryopreserved cells. Variation in postthaw viability is further increased if PBMC isolation and cryopreservation are done more than a few hours after collection. There is currently a lack of evidence-based standards for the minimal PBMC viability or "fitness" required to ensure the integrity and reproducibility of immune cell-based assays. In this study, we use an "induced fail" approach to examine the effect of thawed human PBMC fitness on four flow cytometry-based assays. We found that cell permeability-based viability stains at the time of thawing did not accurately quantify cell fitness, whereas a combined measurement of metabolic activity and early apoptosis markers did. Investigation of the impact of different types and levels of damage on PBMC-based assays revealed that only when cells were >60-70% live and apoptosis negative did biomarker values cease to be determined by cell fitness rather than the inherent biology of the cells. These data show that, to reproducibly measure immunological biomarkers using cryopreserved PBMCs, minimal acceptable standards for cell fitness should be incorporated into the assay protocol.

CRISPR/Cas9 mutagenesis of the Arabidopsis GROWTH-REGULATING FACTOR (GRF) gene family.

Genome editing in plants typically relies on T-DNA plasmids that are mobilized by Agrobacterium-mediated transformation to deliver the CRISPR/Cas machinery. Here, we introduce a series of CRISPR/Cas9 T-DNA vectors for minimal settings, such as teaching labs. Gene-specific targeting sequences can be inserted as annealed short oligonucleotides in a single straightforward cloning step. Fluorescent markers expressed in mature seeds enable reliable selection of transgenic or transgene-free individuals using a combination of inexpensive LED lamps and colored-glass alternative filters. Testing these tools on the Arabidopsis GROWTH-REGULATING FACTOR (GRF) genes, we were able to create a collection of predicted null mutations in all nine family members with little effort. We then explored the effects of simultaneously targeting two, four and eight GRF genes on the rate of induced mutations at each target locus. In our hands, multiplexing was associated with pronounced disparities: while mutation rates at some loci remained consistently high, mutation rates at other loci dropped dramatically with increasing number of single guide RNA species, thereby preventing a systematic mutagenesis of the family.

Disparities in multimorbidity and mortality among people living with and without HIV across British Columbia's health regions: a population-based cohort study.

OBJECTIVES: Longer survival has increased the likelihood of antiretroviral-treated people living with HIV (PLWH) developing age-associated comorbidities. We compared the burden of multimorbidity and all-cause mortality across HIV status in British Columbia (BC), and assessed the longitudinal effect of multimorbidity on all-cause mortality among PLWH. METHODS: Antiretroviral-treated PLWH aged ≥19 years and 1:4 age-sex-matched HIV-negative individuals from a population-based cohort were followed for ≥1 year during 2001-2012. Diagnoses of seven age-associated comorbidities were identified from provincial administrative databases and grouped into 0, 1, 2, and ≥3 comorbidities. Multimorbidity prevalence and age-standardized mortality rates (ASMRs) in both populations were stratified by BC's health regions. Marginal structural models were used to estimate the effect of multimorbidity on mortality among PLWH, adjusted for time-varying confounders affected by prior multimorbidity. RESULTS: Among 8031 PLWH and 32,124 HIV-negative individuals, 25% versus 11% developed multimorbidity, and 23.53 deaths/1000 person-years (95% confidence interval [95% CI]: 22.02-25.13) versus 3.04 (2.81-3.29) were observed, respectively. PLWH in Northern region had the highest ASMR, but those in South Vancouver Island experienced the greatest difference in mortality compared with HIV-negative individuals. Among PLWH, compared with those with zero comorbidities, adjusted hazard ratios for those with 1, 2, and ≥3 comorbidities were 3.36 (95% CI: 2.86-3.95), 6.92 (5.75-8.33), and 12.87 (10.45-15.85), respectively. CONCLUSION: PLWH across BC's health regions experience excess multimorbidity and associated mortality. We highlight health disparities which are key when planning the distribution of healthcare resources across BC, and provide evidence for improved HIV care models integrating prevention and management of chronic diseases.

RéSUMé: OBJECTIF: Les nouvelles thérapeutiques antirétrovirales (ARV) ont permis une plus longue espérance de vie aux personnes porteuses du VIH. Cependant, le vieillissement augmente la probabilité de développer des comorbidités au sein même de la population des personnes vivant avec le VIH (PVVIH) qui suivent des traitements ARV. On a comparé le fardeau de la multimorbidité et mortalité, toutes causes confondues, lié au statut VIH à travers la Colombie-Britannique. On a aussi évalué l'effet longitudinal de la multimorbidité sur la mortalité, toutes causes confondues, parmi les PVVIH. MéTHODES: L'étude comprit des PVVIH suivant un traitement ARV âgés de ≥19 ans et un groupe témoin séronégatif (4 témoins par PVVIH) comparable en termes d'âge et de sexe, tous provenant d'une cohorte de surveillance continue d'au moins 1 an sur une période allant de 2001 à 2012. Des diagnostics de sept comorbidités liées à l'âge ont été identifiés à partir des bases de données administratives provinciales et regroupés en 0, 1, 2 et ≥3 comorbidités. La prévalence de la multimorbidité et les taux de mortalité normalisés selon l'âge (TMNA) dans les deux populations ont été stratifiés selon les régions sociosanitaires de la Colombie-Britannique. Des modèles structurels marginaux ont été utilisés pour estimer l'effet de la multimorbidité sur la mortalité chez les PVVIH, ajustant pour les facteurs de confusion variables affectés par une multimorbidité antérieure. RéSULTATS: Parmi 8 031 PVVIH et 32 124 personnes séronégatives pour le VIH, 25 % contre 11 % ont développé une multimorbidité et 23,53 décès pour 1 000 personnes-années (intervalle de confiance à 95 % [IC à 95 %]: 22,02­25,13) contre 3,04 (2,81­3,29) ont été observés, respectivement. Les PVVIH de la région septentrionale avaient le TMNA le plus élevé, alors que ceux du sud de l'île de Vancouver ont connu la plus grande différence de mortalité par rapport aux personnes séronégatives. Parmi les PVVIH, les personnes atteintes de 1, 2 et ≥3 comorbidités avaient respectivement 3,36 (IC à 95 %: 2,86­3,95), 6,92 (5,75­8,33) et 12,87 (10,45­15,85) fois plus la probabilité de mourir que les personnes sans comorbidités. CONCLUSION: Les PVVIH des régions sociosanitaires de la Colombie-Britannique connaissent une multimorbidité excessive et la surmortalité s'y associant. Notre étude souligne les disparités-clés en matière de santé qu'il faut prendre en compte lors de la planification de la distribution des ressources de soins de santé à travers la Colombie-Britannique. Elle fournit aussi des preuves pour des modèles de soins du VIH améliorés, y intégrant la prévention et la gestion des maladies chroniques.

Integrated Single-Cell Transcriptomics and Chromatin Accessibility Analysis Reveals Regulators of Mammary Epithelial Cell Identity.

The mammary epithelial cell (MEC) system is a bilayered ductal epithelium of luminal and basal cells, maintained by a lineage of stem and progenitor populations. Here, we used integrated single-cell transcriptomics and chromatin accessibility analysis to reconstruct the cell types of the mouse MEC system and their underlying gene regulatory features in an unbiased manner. We define differentiation states within the secretory type of luminal cells, which forms a continuous spectrum of general luminal progenitor and lactation-committed progenitor cells. By integrating single-cell transcriptomics and chromatin accessibility landscapes, we identify cis- and trans-regulatory elements that are differentially activated in the specific epithelial cell types and our newly defined luminal differentiation states. Our work provides a resource to reveal cis/trans-regulatory elements associated with MEC identity and differentiation that will serve as a reference to determine how the chromatin accessibility landscape changes during breast cancer.

Progenitor identification and SARS-CoV-2 infection in long-term human distal lung organoid cultures.

The distal lung contains terminal bronchioles and alveoli that facilitate gas exchange and is affected by disorders including interstitial lung disease, cancer, and SARS-CoV-2-associated COVID-19 pneumonia. Investigations of these localized pathologies have been hindered by a lack of 3D in vitro human distal lung culture systems. Further, human distal lung stem cell identification has been impaired by quiescence, anatomic divergence from mouse and lack of lineage tracing and clonogenic culture. Here, we developed robust feeder-free, chemically-defined culture of distal human lung progenitors as organoids derived clonally from single adult human alveolar epithelial type II (AT2) or KRT5 + basal cells. AT2 organoids exhibited AT1 transdifferentiation potential, while basal cell organoids progressively developed lumens lined by differentiated club and ciliated cells. Organoids consisting solely of club cells were not observed. Upon single cell RNA-sequencing (scRNA-seq), alveolar organoids were composed of proliferative AT2 cells; however, basal organoid KRT5 + cells contained a distinct ITGA6 + ITGB4 + mitotic population whose proliferation segregated to a TNFRSF12A hi subfraction. Clonogenic organoid growth was markedly enriched within the TNFRSF12A hi subset of FACS-purified ITGA6 + ITGB4 + basal cells from human lung or derivative organoids. In vivo, TNFRSF12A + cells comprised ~10% of KRT5 + basal cells and resided in clusters within terminal bronchioles. To model COVID-19 distal lung disease, we everted the polarity of basal and alveolar organoids to rapidly relocate differentiated club and ciliated cells from the organoid lumen to the exterior surface, thus displaying the SARS-CoV-2 receptor ACE2 on the outwardly-facing apical aspect. Accordingly, basal and AT2 apical-out organoids were infected by SARS-CoV-2, identifying club cells as a novel target population. This long-term, feeder-free organoid culture of human distal lung alveolar and basal stem cells, coupled with single cell analysis, identifies unsuspected basal cell functional heterogeneity and exemplifies progenitor identification within a slowly proliferating human tissue. Further, our studies establish a facile in vitro organoid model for human distal lung infectious diseases including COVID-19-associated pneumonia.

Single-cell multiomic analysis identifies regulatory programs in mixed-phenotype acute leukemia.

Identifying the causes of human diseases requires deconvolution of abnormal molecular phenotypes spanning DNA accessibility, gene expression and protein abundance1-3. We present a single-cell framework that integrates highly multiplexed protein quantification, transcriptome profiling and analysis of chromatin accessibility. Using this approach, we establish a normal epigenetic baseline for healthy blood development, which we then use to deconvolve aberrant molecular features within blood from patients with mixed-phenotype acute leukemia4,5. Despite widespread epigenetic heterogeneity within the patient cohort, we observe common malignant signatures across patients as well as patient-specific regulatory features that are shared across phenotypic compartments of individual patients. Integrative analysis of transcriptomic and chromatin-accessibility maps identified 91,601 putative peak-to-gene linkages and transcription factors that regulate leukemia-specific genes, such as RUNX1-linked regulatory elements proximal to the marker gene CD69. These results demonstrate how integrative, multiomic analysis of single cells within the framework of normal development can reveal both distinct and shared molecular mechanisms of disease from patient samples.

Massively parallel single-cell chromatin landscapes of human immune cell development and intratumoral T cell exhaustion.

Understanding complex tissues requires single-cell deconstruction of gene regulation with precision and scale. Here, we assess the performance of a massively parallel droplet-based method for mapping transposase-accessible chromatin in single cells using sequencing (scATAC-seq). We apply scATAC-seq to obtain chromatin profiles of more than 200,000 single cells in human blood and basal cell carcinoma. In blood, application of scATAC-seq enables marker-free identification of cell type-specific cis- and trans-regulatory elements, mapping of disease-associated enhancer activity and reconstruction of trajectories of cellular differentiation. In basal cell carcinoma, application of scATAC-seq reveals regulatory networks in malignant, stromal and immune cells in the tumor microenvironment. Analysis of scATAC-seq profiles from serial tumor biopsies before and after programmed cell death protein 1 blockade identifies chromatin regulators of therapy-responsive T cell subsets and reveals a shared regulatory program that governs intratumoral CD8+ T cell exhaustion and CD4+ T follicular helper cell development. We anticipate that scATAC-seq will enable the unbiased discovery of gene regulatory factors across diverse biological systems.

Massively parallel digital transcriptional profiling of single cells.

Characterizing the transcriptome of individual cells is fundamental to understanding complex biological systems. We describe a droplet-based system that enables 3' mRNA counting of tens of thousands of single cells per sample. Cell encapsulation, of up to 8 samples at a time, takes place in ∼6 min, with ∼50% cell capture efficiency. To demonstrate the system's technical performance, we collected transcriptome data from ∼250k single cells across 29 samples. We validated the sensitivity of the system and its ability to detect rare populations using cell lines and synthetic RNAs. We profiled 68k peripheral blood mononuclear cells to demonstrate the system's ability to characterize large immune populations. Finally, we used sequence variation in the transcriptome data to determine host and donor chimerism at single-cell resolution from bone marrow mononuclear cells isolated from transplant patients.

Intestinal Enteroendocrine Lineage Cells Possess Homeostatic and Injury-Inducible Stem Cell Activity.

Several cell populations have been reported to possess intestinal stem cell (ISC) activity during homeostasis and injury-induced regeneration. Here, we explored inter-relationships between putative mouse ISC populations by comparative RNA-sequencing (RNA-seq). The transcriptomes of multiple cycling ISC populations closely resembled Lgr5+ ISCs, the most well-defined ISC pool, but Bmi1-GFP+ cells were distinct and enriched for enteroendocrine (EE) markers, including Prox1. Prox1-GFP+ cells exhibited sustained clonogenic growth in vitro, and lineage-tracing of Prox1+ cells revealed long-lived clones during homeostasis and after radiation-induced injury in vivo. Single-cell mRNA-seq revealed two subsets of Prox1-GFP+ cells, one of which resembled mature EE cells while the other displayed low-level EE gene expression but co-expressed tuft cell markers, Lgr5 and Ascl2, reminiscent of label-retaining secretory progenitors. Our data suggest that the EE lineage, including mature EE cells, comprises a reservoir of homeostatic and injury-inducible ISCs, extending our understanding of cellular plasticity and stemness.