SARS-CoV-2 mRNA vaccination-induced immunological memory in human nonlymphoid and lymphoid tissues.

Tissue-resident lymphocytes provide organ-adapted protection against invading pathogens. Whereas their biology has been examined in great detail in various infection models, their generation and functionality in response to vaccination have not been comprehensively analyzed in humans. We therefore studied SARS-CoV-2 mRNA vaccine-specific T cells in surgery specimens of kidney, liver, lung, bone marrow, and spleen compared with paired blood samples from largely virus-naive individuals. As opposed to lymphoid tissues, nonlymphoid organs harbored significantly elevated frequencies of spike-specific CD4+ T cells compared with blood showing hallmarks of tissue residency and an expanded memory pool. Organ-derived CD4+ T cells further exhibited increased polyfunctionality over those detected in blood. Single-cell RNA-Seq together with T cell receptor repertoire analysis indicated that the clonotype rather than organ origin is a major determinant of transcriptomic state in vaccine-specific CD4+ T cells. In summary, our data demonstrate that SARS-CoV-2 vaccination entails acquisition of tissue memory and residency features in organs distant from the inoculation site, thereby contributing to our understanding of how local tissue protection might be accomplished.

90K/LGALS3BP expression is upregulated in COVID-19 but may not restrict SARS-CoV-2 infection.

Glycoprotein 90K, encoded by the interferon-stimulated gene LGALS3BP, displays broad antiviral activity. It reduces HIV-1 infectivity by interfering with Env maturation and virion incorporation, and increases survival of Influenza A virus-infected mice via antiviral innate immune signaling. Its antiviral potential in SARS-CoV-2 infection remains largely unknown. Here, we analyzed the expression of 90K/LGALS3BP in 44 hospitalized COVID-19 patients at multiple levels. We quantified 90K protein concentrations in serum and PBMCs as well as LGALS3BP mRNA levels. Complementary, we analyzed two single cell RNA-sequencing datasets for expression of LGALS3BP in respiratory specimens and PBMCs from COVID-19 patients. Finally, we analyzed the potential of 90K to interfere with SARS-CoV-2 infection of HEK293T/ACE2, Calu-3 and Caco-2 cells using authentic virus. 90K protein serum concentrations were significantly elevated in COVID-19 patients compared to uninfected sex- and age-matched controls. Furthermore, PBMC-associated concentrations of 90K protein were overall reduced by SARS-CoV-2 infection in vivo, suggesting enhanced secretion into the extracellular space. Mining of published PBMC scRNA-seq datasets uncovered monocyte-specific induction of LGALS3BP mRNA expression in COVID-19 patients. In functional assays, neither 90K overexpression in susceptible cell lines nor exogenous addition of purified 90K consistently inhibited SARS-CoV-2 infection. Our data suggests that 90K/LGALS3BP contributes to the global type I IFN response during SARS-CoV-2 infection in vivo without displaying detectable antiviral properties in vitro.

Patient-specific modeling of stroma-mediated chemoresistance of pancreatic cancer using a three-dimensional organoid-fibroblast co-culture system.

BACKGROUND: Cancer-associated fibroblasts (CAFs) are considered to play a fundamental role in pancreatic ductal adenocarcinoma (PDAC) progression and chemoresistance. Patient-derived organoids have demonstrated great potential as tumor avatars for drug response prediction in PDAC, yet they disregard the influence of stromal components on chemosensitivity. METHODS: We established direct three-dimensional (3D) co-cultures of primary PDAC organoids and patient-matched CAFs to investigate the effect of the fibroblastic compartment on sensitivity to gemcitabine, 5-fluorouracil and paclitaxel treatments using an image-based drug assay. Single-cell RNA sequencing was performed for three organoid/CAF pairs in mono- and co-culture to uncover transcriptional changes induced by tumor-stroma interaction. RESULTS: Upon co-culture with CAFs, we observed increased proliferation and reduced chemotherapy-induced cell death of PDAC organoids. Single-cell RNA sequencing data evidenced induction of a pro-inflammatory phenotype in CAFs in co-cultures. Organoids showed increased expression of genes associated with epithelial-to-mesenchymal transition (EMT) in co-cultures and several potential receptor-ligand interactions related to EMT were identified, supporting a key role of CAF-driven induction of EMT in PDAC chemoresistance. CONCLUSIONS: Our results demonstrate the potential of personalized PDAC co-cultures models not only for drug response profiling but also for unraveling the molecular mechanisms involved in the chemoresistance-supporting role of the tumor stroma.

ELF5 is a potential respiratory epithelial cell-specific risk gene for severe COVID-19.

Despite two years of intense global research activity, host genetic factors that predispose to a poorer prognosis of COVID-19 infection remain poorly understood. Here, we prioritise eight robust (e.g., ELF5) or suggestive but unreported (e.g., RAB2A) candidate protein mediators of COVID-19 outcomes by integrating results from the COVID-19 Host Genetics Initiative with population-based plasma proteomics using statistical colocalisation. The transcription factor ELF5 (ELF5) shows robust and directionally consistent associations across different outcome definitions, including a >4-fold higher risk (odds ratio: 4.88; 95%-CI: 2.47-9.63; p-value < 5.0 × 10-6) for severe COVID-19 per 1 s.d. higher genetically predicted plasma ELF5. We show that ELF5 is specifically expressed in epithelial cells of the respiratory system, such as secretory and alveolar type 2 cells, using single-cell RNA sequencing and immunohistochemistry. These cells are also likely targets of SARS-CoV-2 by colocalisation with key host factors, including ACE2 and TMPRSS2. In summary, large-scale human genetic studies together with gene expression at single-cell resolution highlight ELF5 as a risk gene for severe COVID-19, supporting a role of epithelial cells of the respiratory system in the adverse host response to SARS-CoV-2.

Intratumoral Heterogeneity and Immune Modulation in Lung Adenocarcinoma in Female Smokers and Never Smokers.

SIGNIFICANCE: Single-cell analysis of healthy lung tissue and lung cancer reveals distinct tumor cell populations, including cells with differential immune modulating capacity between smokers and never smokers, which could guide future therapeutic strategies.

Single-cell analysis of patient-derived PDAC organoids reveals cell state heterogeneity and a conserved developmental hierarchy.

Pancreatic ductal adenocarcinoma (PDAC) is projected to be the second leading cause of cancer mortality by 2030. Bulk transcriptomic analyses have distinguished 'classical' from 'basal-like' tumors with more aggressive clinical behavior. We derive PDAC organoids from 18 primary tumors and two matched liver metastases, and show that 'classical' and 'basal-like' cells coexist in individual organoids. By single-cell transcriptome analysis of PDAC organoids and primary PDAC, we identify distinct tumor cell states shared across patients, including a cycling progenitor cell state and a differentiated secretory state. Cell states are connected by a differentiation hierarchy, with 'classical' cells concentrated at the endpoint. In an imaging-based drug screen, expression of 'classical' subtype genes correlates with better drug response. Our results thus uncover a functional hierarchy of PDAC cell states linked to transcriptional tumor subtypes, and support the use of PDAC organoids as a clinically relevant model for in vitro studies of tumor heterogeneity.

Functional States in Tumor-Initiating Cell Differentiation in Human Colorectal Cancer.

Intra-tumor heterogeneity of tumor-initiating cell (TIC) activity drives colorectal cancer (CRC) progression and therapy resistance. Here, we used single-cell RNA-sequencing of patient-derived CRC models to decipher distinct cell subpopulations based on their transcriptional profiles. Cell type-specific expression modules of stem-like, transit amplifying-like, and differentiated CRC cells resemble differentiation states of normal intestinal epithelial cells. Strikingly, identified subpopulations differ in proliferative activity and metabolic state. In summary, we here show at single-cell resolution that transcriptional heterogeneity identifies functional states during TIC differentiation. Furthermore, identified expression signatures are linked to patient prognosis. Targeting transcriptional states associated to cancer cell differentiation might unravel novel vulnerabilities in human CRC.

Single-Nucleus and In Situ RNA-Sequencing Reveal Cell Topographies in the Human Pancreas.

BACKGROUND & AIMS: Molecular evidence of cellular heterogeneity in the human exocrine pancreas has not been yet established because of the local concentration and cascade of hydrolytic enzymes that can rapidly degrade cells and RNA upon pancreatic resection. We sought to better understand the heterogeneity and cellular composition of the pancreas in neonates and adults in healthy and diseased conditions using single-cell sequencing approaches. METHODS: We innovated single-nucleus RNA-sequencing protocols and profiled more than 120,000 cells from pancreata of adult and neonatal human donors. We validated the single-nucleus findings using RNA fluorescence in situ hybridization, in situ sequencing, and computational approaches. RESULTS: We created the first comprehensive atlas of human pancreas cells including epithelial and nonepithelial constituents, and uncovered 3 distinct acinar cell types, with possible implications for homeostatic and inflammatory processes of the pancreas. The comparison with neonatal single-nucleus sequencing data showed a different cellular composition of the endocrine tissue, highlighting the tissue dynamics occurring during development. By applying spatial cartography, involving cell proximity mapping through in situ sequencing, we found evidence of specific cell type neighborhoods, dynamic topographies in the endocrine and exocrine pancreas, and principles of morphologic organization of the organ. Furthermore, similar analyses in chronic pancreatitis biopsy samples showed the presence of acinar-REG+ cells, a reciprocal association between macrophages and activated stellate cells, and a new potential role of tuft cells in this disease. CONCLUSIONS: Our human pancreas cell atlas can be interrogated to understand pancreatic cell biology and provides a crucial reference set for comparisons with diseased tissue samples to map the cellular foundations of pancreatic diseases.

COVID-19 severity correlates with airway epithelium-immune cell interactions identified by single-cell analysis.

To investigate the immune response and mechanisms associated with severe coronavirus disease 2019 (COVID-19), we performed single-cell RNA sequencing on nasopharyngeal and bronchial samples from 19 clinically well-characterized patients with moderate or critical disease and from five healthy controls. We identified airway epithelial cell types and states vulnerable to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. In patients with COVID-19, epithelial cells showed an average three-fold increase in expression of the SARS-CoV-2 entry receptor ACE2, which correlated with interferon signals by immune cells. Compared to moderate cases, critical cases exhibited stronger interactions between epithelial and immune cells, as indicated by ligand-receptor expression profiles, and activated immune cells, including inflammatory macrophages expressing CCL2, CCL3, CCL20, CXCL1, CXCL3, CXCL10, IL8, IL1B and TNF. The transcriptional differences in critical cases compared to moderate cases likely contribute to clinical observations of heightened inflammatory tissue damage, lung injury and respiratory failure. Our data suggest that pharmacologic inhibition of the CCR1 and/or CCR5 pathways might suppress immune hyperactivation in critical COVID-19.

Modeling glioblastoma invasion using human brain organoids and single-cell transcriptomics.

BACKGROUND: Glioblastoma (GBM) consists of devastating neoplasms with high invasive capacity, which have been difficult to study in vitro in a human-derived model system. Therapeutic progress is also limited by cellular heterogeneity within and between tumors, among other factors such as therapy resistance. To address these challenges, we present an experimental model using human cerebral organoids as a scaffold for patient-derived GBM cell invasion. METHODS: This study combined tissue clearing and confocal microscopy with single-cell RNA sequencing of GBM cells before and after co-culture with organoid cells. RESULTS: We show that tumor cells within organoids extend a network of long microtubes, recapitulating the in vivo behavior of GBM. Transcriptional changes implicated in the invasion process are coherent across patient samples, indicating that GBM cells reactively upregulate genes required for their dispersion. Potential interactions between GBM and organoid cells identified by an in silico receptor-ligand pairing screen suggest functional therapeutic targets. CONCLUSIONS: Taken together, our model has proven useful for studying GBM invasion and transcriptional heterogeneity in vitro, with applications for both pharmacological screens and patient-specific treatment selection on a time scale amenable to clinical practice.

SARS-CoV-2 receptor ACE2 and TMPRSS2 are primarily expressed in bronchial transient secretory cells.

The SARS-CoV-2 pandemic affecting the human respiratory system severely challenges public health and urgently demands for increasing our understanding of COVID-19 pathogenesis, especially host factors facilitating virus infection and replication. SARS-CoV-2 was reported to enter cells via binding to ACE2, followed by its priming by TMPRSS2. Here, we investigate ACE2 and TMPRSS2 expression levels and their distribution across cell types in lung tissue (twelve donors, 39,778 cells) and in cells derived from subsegmental bronchial branches (four donors, 17,521 cells) by single nuclei and single cell RNA sequencing, respectively. While TMPRSS2 is strongly expressed in both tissues, in the subsegmental bronchial branches ACE2 is predominantly expressed in a transient secretory cell type. Interestingly, these transiently differentiating cells show an enrichment for pathways related to RHO GTPase function and viral processes suggesting increased vulnerability for SARS-CoV-2 infection. Our data provide a rich resource for future investigations of COVID-19 infection and pathogenesis.

Fit-for-purpose treatment goals for produced waters in shale oil and gas fields.

Hydraulic fracturing (HF), or "fracking," is the driving force behind the "shale gas revolution," completely transforming the United States energy industry over the last two decades. HF requires that 4-6 million gallons per well (15,000-23,000 m3/well) of water be pumped underground to stimulate the release of entrapped hydrocarbons from unconventional (i.e., shale or carbonate) formations. Estimated U.S. produced water volumes exceed 150 billion gallons/year across the industry from unconventional wells alone and are projected to grow for at least another two decades. Concerns over the environmental impact from accidental or incidental release of produced water from HF wells ("U-PW"), along with evolving regulatory and economic drivers, has spurred great interest in technological innovation to enhance U-PW recycling and reuse. In this review, we analyze U-PW quantity and composition based on the latest U.S. Geographical Survey data, identify key contamination metrics useful in tracking water quality improvement in the context of HF operations, and suggest "fit-for-purpose treatment" to enhance cost-effective regulatory compliance, water recovery/reuse, and resource valorization. Drawing on industrial practice and technoeconomic constraints, we further assess the challenges associated with U-PW treatment for onshore U.S. operations. Presented are opportunities for targeted end-uses of treated U-PW. We highlight emerging technologies that may enhance cost-effective U-PW management as HF activities grow and evolve in the coming decades.

In Situ Methods for Identifying Reactive Surface Intermediates during Hydrogenolysis Reactions: C-O Bond Cleavage on Nanoparticles of Nickel and Nickel Phosphides.

Identifying individual reactive intermediates within the "zoo" of organometallic species that form on catalytic surfaces during reactions is a long-standing challenge in heterogeneous catalysis. Here, we identify distinct reactive intermediates, all of which exist at low coverages, that lead to distinguishable reaction pathways during the hydrogenolysis of 2-methyltetrahydrofuran (MTHF) on Ni, Ni12P5, and Ni2P catalysts by combining advanced spectroscopic methods with quantum chemical calculations. Each of these reactive complexes cleaves specific C-O bonds, gives rise to unique products, and exhibits different apparent activation barriers for ring opening. The spectral features of the reactive intermediates are extracted by collecting in situ infrared spectra while sinusoidally modulating the H2 pressure during MTHF hydrogenolysis and applying phase-sensitive detection (PSD), which suppresses the features of inactive surface species. The combined spectra of all reactive species are deconvoluted using singular-value decomposition techniques that yield spectra and changes in surface coverage for each set of kinetically differentiable species. These deconvoluted spectra are consistent with predicted spectral features for the reactive surface intermediates implicated by detailed kinetic measurements and DFT calculations. Notably, these methods give direct evidence for several anticipated differences in the coordination and composition of reactive MTHF-derived species. The compositions of the most abundant reactive intermediate (MARI) on Ni, Ni12P5, and Ni2P nanoparticles during the C-O bond rupture of MTHF are identical; however, MARI changes orientation from Ni3(µ3-C5H10O) to Ni3(η5-C5H10O) (i.e., lies more parallel with the catalyst surface) with increasing phosphorus content. The shift in binding configuration with phosphorus content suggests that the decrease in steric hindrance to rupture the 3C-O bond is the fundamental cause of increased selectivity toward 3C-O bond rupture. Previous kinetic measurements and DFT calculations indicate that C-O bond rupture occurs on Ni ensembles on Ni, Ni12P5, and Ni2P catalysts; however, the addition of more electronegative phosphorus atoms that withdraw a small charge from Ni ensembles results in differences in the binding configuration, activation enthalpy, and selectivity. The results from this in situ spectroscopic methodology support previous proposals that the manipulation of the electronic structure of metal ensembles by the introduction of phosphorus provides strategies for designing catalysts for the selective cleavage of hindered C-X bonds and demonstrate the utility of this approach in identifying individual reactive species within the zoo.

Pheno-seq - linking visual features and gene expression in 3D cell culture systems.

Patient-derived 3D cell culture systems are currently advancing cancer research since they potentiate the molecular analysis of tissue-like properties and drug response under well-defined conditions. However, our understanding of the relationship between the heterogeneity of morphological phenotypes and the underlying transcriptome is still limited. To address this issue, we here introduce "pheno-seq" to directly link visual features of 3D cell culture systems with profiling their transcriptome. As prototypic applications breast and colorectal cancer (CRC) spheroids were analyzed by pheno-seq. We identified characteristic gene expression signatures of epithelial-to-mesenchymal transition that are associated with invasive growth behavior of clonal breast cancer spheroids. Furthermore, we linked long-term proliferative capacity in a patient-derived model of CRC to a lowly abundant PROX1-positive cancer stem cell subtype. We anticipate that the ability to integrate transcriptome analysis and morphological patho-phenotypes of cancer cells will provide novel insight on the molecular origins of intratumor heterogeneity.

Unraveling mitotic protein networks by 3D multiplexed epitope drug screening.

Three-dimensional protein localization intricately determines the functional coordination of cellular processes. The complex spatial context of protein landscape has been assessed by multiplexed immunofluorescent staining or mass spectrometry, applied to 2D cell culture with limited physiological relevance or tissue sections. Here, we present 3D SPECS, an automated technology for 3D Spatial characterization of Protein Expression Changes by microscopic Screening. This workflow comprises iterative antibody staining, high-content 3D imaging, and machine learning for detection of mitoses. This is followed by mapping of spatial protein localization into a spherical, cellular coordinate system, a basis for model-based prediction of spatially resolved affinities of proteins. As a proof-of-concept, we mapped twelve epitopes in 3D-cultured spheroids and investigated the network effects of twelve mitotic cancer drugs. Our approach reveals novel insights into spindle fragility and chromatin stress, and predicts unknown interactions between proteins in specific mitotic pathways. 3D SPECS's ability to map potential drug targets by multiplexed immunofluorescence in 3D cell culture combined with our automated high-content assay will inspire future functional protein expression and drug assays.

Screening drug effects in patient-derived cancer cells links organoid responses to genome alterations.

Cancer drug screening in patient-derived cells holds great promise for personalized oncology and drug discovery but lacks standardization. Whether cells are cultured as conventional monolayer or advanced, matrix-dependent organoid cultures influences drug effects and thereby drug selection and clinical success. To precisely compare drug profiles in differently cultured primary cells, we developed DeathPro, an automated microscopy-based assay to resolve drug-induced cell death and proliferation inhibition. Using DeathPro, we screened cells from ovarian cancer patients in monolayer or organoid culture with clinically relevant drugs. Drug-induced growth arrest and efficacy of cytostatic drugs differed between the two culture systems. Interestingly, drug effects in organoids were more diverse and had lower therapeutic potential. Genomic analysis revealed novel links between drug sensitivity and DNA repair deficiency in organoids that were undetectable in monolayers. Thus, our results highlight the dependency of cytostatic drugs and pharmacogenomic associations on culture systems, and guide culture selection for drug tests.

Islamic Medicine in the Middle Ages.

The Islamic culture flourished between the 9th and 13th centuries. Scholars from this era made significant contributions in mathematics, science and medicine. Caliphs and physicians built hospitals that provided universal care and the foundation for medical education. Physician-scientists made significant advances in medical care, surgery and pharmacology. Notable authorities include al-Razi (865-925 CE) who wrote the Kitab al-Hawi fi al-tibb (The Comprehensive Book on Medicine), a 23-volume textbook that provided the main medical curriculum for European schools into the 14th century. Ibn Sina (980-1037 CE), an extraordinary Persian polymath, wrote al Qanun fi al-Tibb (The Canon of Medicine), an encyclopedic treatment of medicine that combined his own observations with medical information from Galen and philosophy from Aristotle. Mansur (1380-1422 CE) wrote the first color illustrated book on anatomy. Other important physicians compiled information on the use of medication from plants, advanced surgical techniques, including cataract extraction and studied physiology, including the pulmonary circulation. These books and ideas provided the basis for medical care in Europe during its recovery from the Dark Ages.

High-content siRNA screen reveals global ENaC regulators and potential cystic fibrosis therapy targets.

Dysfunction of ENaC, the epithelial sodium channel that regulates salt and water reabsorption in epithelia, causes several human diseases, including cystic fibrosis (CF). To develop a global understanding of molecular regulators of ENaC traffic/function and to identify of candidate CF drug targets, we performed a large-scale screen combining high-content live-cell microscopy and siRNAs in human airway epithelial cells. Screening over 6,000 genes identified over 1,500 candidates, evenly divided between channel inhibitors and activators. Genes in the phosphatidylinositol pathway were enriched on the primary candidate list, and these, along with other ENaC activators, were examined further with secondary siRNA validation. Subsequent detailed investigation revealed ciliary neurotrophic factor receptor (CNTFR) as an ENaC modulator and showed that inhibition of (diacylglycerol kinase, iota) DGKι, a protein involved in PiP2 metabolism, downgrades ENaC activity, leading to normalization of both Na+ and fluid absorption in CF airways to non-CF levels in primary human lung cells from CF patients.

Functional genomics assays to study CFTR traffic and ENaC function.

As several genomes have been sequenced, post-genomic approaches like transcriptomics and proteomics, identifying gene products differentially expressed in association with a given pathology, have held promise both of understanding the pathways associated with the respective disease and as a fast track to therapy. Notwithstanding, these approaches cannot distinguish genes and proteins with mere secondary pathological association from those primarily involved in the basic defect(s). New global strategies and tools identifying gene products responsible for the basic cellular defect(s) in CF pathophysiology currently being performed are presented here. These include high-content screens to determine proteins affecting function and trafficking of CFTR and ENaC.

Automatic identification of subcellular phenotypes on human cell arrays.

Light microscopic analysis of cell morphology provides a high-content readout of cell function and protein localization. Cell arrays and microwell transfection assays on cultured cells have made cell phenotype analysis accessible to high-throughput experiments. Both the localization of each protein in the proteome and the effect of RNAi knock-down of individual genes on cell morphology can be assayed by manual inspection of microscopic images. However, the use of morphological readouts for functional genomics requires fast and automatic identification of complex cellular phenotypes. Here, we present a fully automated platform for high-throughput cell phenotype screening combining human live cell arrays, screening microscopy, and machine-learning-based classification methods. Efficiency of this platform is demonstrated by classification of eleven subcellular patterns marked by GFP-tagged proteins. Our classification method can be adapted to virtually any microscopic assay based on cell morphology, opening a wide range of applications including large-scale RNAi screening in human cells.