Spatial Deconvolution of Cell Types and Cell States at Scale Utilizing TACIT.

Identifying cell types and states remains a time-consuming and error-prone challenge for spatial biology. While deep learning is increasingly used, it is difficult to generalize due to variability at the level of cells, neighborhoods, and niches in health and disease. To address this, we developed TACIT, an unsupervised algorithm for cell annotation using predefined signatures that operates without training data, using unbiased thresholding to distinguish positive cells from background, focusing on relevant markers to identify ambiguous cells in multiomic assays. Using five datasets (5,000,000-cells; 51-cell types) from three niches (brain, intestine, gland), TACIT outperformed existing unsupervised methods in accuracy and scalability. Integration of TACIT-identified cell with a novel Shiny app revealed new phenotypes in two inflammatory gland diseases. Finally, using combined spatial transcriptomics and proteomics, we discover under- and overrepresented immune cell types and states in regions of interest, suggesting multimodality is essential for translating spatial biology to clinical applications.

Inhibition of JAK-STAT pathway corrects salivary gland inflammation and interferon driven immune activation in Sjögren's disease.

OBJECTIVES: Inflammatory cytokines that signal through the Janus kinases-signal transducer and activator of transcription (JAK-STAT) pathway, especially interferons (IFNs), are implicated in Sjögren's disease (SjD). Although inhibition of JAKs is effective in other autoimmune diseases, a systematic investigation of IFN-JAK-STAT signalling and the effect of JAK inhibitor (JAKi) therapy in SjD-affected human tissues has not been fully investigated. METHODS: Human minor salivary glands (MSGs) and peripheral blood mononuclear cells (PBMCs) were investigated using bulk or single-cell (sc) RNA sequencing (RNAseq), immunofluorescence (IF) microscopy and flow cytometry. Ex vivo culture assays on PBMCs and primary salivary gland epithelial cell (pSGEC) lines were performed to model changes in target tissues before and after JAKi. RESULTS: RNAseq and IF showed activated JAK-STAT pathway in SjD MSGs. Elevated IFN-stimulated gene (ISGs) expression associated with clinical variables (eg, focus scores, anti-SSA positivity). scRNAseq of MSGs exhibited cell type-specific upregulation of JAK-STAT and ISGs; PBMCs showed similar trends, including markedly upregulated ISGs in monocytes. Ex vivo studies showed elevated basal pSTAT levels in SjD MSGs and PBMCs that were corrected with JAKi. SjD-derived pSGECs exhibited higher basal ISG expressions and exaggerated responses to IFN-ß, which were normalised by JAKi without cytotoxicity. CONCLUSIONS: SjD patients' tissues exhibit increased expression of ISGs and activation of the JAK-STAT pathway in a cell type-dependent manner. JAKi normalises this aberrant signalling at the tissue level and in PBMCs, suggesting a putative viable therapy for SjD, targeting both glandular and extraglandular symptoms. Predicated on these data, a phase Ib/IIa randomised controlled trial to treat SjD with tofacitinib was initiated.

TScan-II: A genome-scale platform for the de novo identification of CD4+ T cell epitopes.

CD4+ T cells play fundamental roles in orchestrating immune responses and tissue homeostasis. However, our inability to associate peptide human leukocyte antigen class-II (HLA-II) complexes with their cognate T cell receptors (TCRs) in an unbiased manner has hampered our understanding of CD4+ T cell function and role in pathologies. Here, we introduce TScan-II, a highly sensitive genome-scale CD4+ antigen discovery platform. This platform seamlessly integrates the endogenous HLA-II antigen-processing machinery in synthetic antigen-presenting cells and TCR signaling in T cells, enabling the simultaneous screening of multiple HLAs and TCRs. Leveraging genome-scale human, virome, and epitope mutagenesis libraries, TScan-II facilitates de novo antigen discovery and deep exploration of TCR specificity. We demonstrate TScan-II's potential for basic and translational research by identifying a non-canonical antigen for a cancer-reactive CD4+ T cell clone. Additionally, we identified two antigens for clonally expanded CD4+ T cells in Sjögren's disease, which bind distinct HLAs and are expressed in HLA-II-positive ductal cells within affected salivary glands.

Inhibition of JAK-STAT pathway corrects salivary gland inflammation and interferon driven immune activation in Sjögren's Disease.

Objectives: Inflammatory cytokines that signal through the JAK- STAT pathway, especially interferons (IFNs), are implicated in Sjögren's Disease (SjD). Although inhibition of JAKs is effective in other autoimmune diseases, a systematic investigation of IFN-JAK-STAT signaling and effect of JAK inhibitor (JAKi) therapy in SjD-affected human tissues has not been reported. Methods: Human minor salivary glands (MSGs) and peripheral blood mononuclear cells (PBMCs) were investigated using bulk or single cell (sc) RNA sequencing (RNAseq), immunofluorescence microscopy (IF), and flow cytometry. Ex vivo culture assays on PBMCs and primary salivary gland epithelial cell (pSGEC) lines were performed to model changes in target tissues before and after JAKi. Results: RNAseq and IF showed activated JAK-STAT pathway in SjD MSGs. Elevated IFN-stimulated gene (ISGs) expression associated with clinical variables (e.g., focus scores, anti-SSA positivity). scRNAseq of MSGs exhibited cell-type specific upregulation of JAK-STAT and ISGs; PBMCs showed similar trends, including markedly upregulated ISGs in monocytes. Ex vivo studies showed elevated basal pSTAT levels in SjD MSGs and PBMCs that were corrected with JAKi. SjD-derived pSGECs exhibited higher basal ISG expressions and exaggerated responses to IFNß, which were normalized by JAKi without cytotoxicity. Conclusions: SjD patients' tissues exhibit increased expression of ISGs and activation of the JAK-STAT pathway in a cell type-dependent manner. JAKi normalizes this aberrant signaling at the tissue level and in PBMCs, suggesting a putative viable therapy for SjD, targeting both glandular and extraglandular symptoms. Predicated on these data, a Phase Ib/IIa randomized controlled trial to treat SjD with tofacitinib was initiated.

GZMK+CD8+ T cells Target a Specific Acinar Cell Type in Sjögren's Disease.

Sjögren's Disease (SjD) is a systemic autoimmune disease without a clear etiology or effective therapy. Utilizing unbiased single-cell and spatial transcriptomics to analyze human minor salivary glands in health and disease we developed a comprehensive understanding of the cellular landscape of healthy salivary glands and how that landscape changes in SjD patients. We identified novel seromucous acinar cell types and identified a population of PRR4+CST3+WFDC2- seromucous acinar cells that are particularly targeted in SjD. Notably, GZMK+CD8 T cells, enriched in SjD, exhibited a cytotoxic phenotype and were physically associated with immune-engaged epithelial cells in disease. These findings shed light on the immune response's impact on transitioning acinar cells with high levels of secretion and explain the loss of this specific cell population in SjD. This study explores the complex interplay of varied cell types in the salivary glands and their role in the pathology of Sjögren's Disease.

Lysosomal exocytosis of HSP70 stimulates monocytic BMP6 expression in Sjögren's syndrome.

BMP6 is a central cytokine in the induction of Sjögren's syndrome-associated (SS-associated) secretory hypofunction. However, the upstream initiation leading to the production of this cytokine in SS is unknown. In this study, RNA ISH on salivary gland sections taken from patients with SS indicated monocytic lineage cells as a cellular source of BMP6. RNA-Seq data on human salivary glands suggested that TLR4 signaling was an upstream regulator of BMP6, which was confirmed by in vitro cell assays and single-cell transcriptomics of human PBMCs. Further investigation showed that HSP70 was an endogenous natural TLR4 ligand that stimulated BMP6 expression in SS. Release of HSP70 from epithelial cells could be triggered by overexpression of lysosome-associated membrane protein 3 (LAMP3), a protein also associated with SS in several transcriptome studies. In vitro studies supported the idea that HSP70 was released as a result of lysosomal exocytosis initiated by LAMP3 expression, and reverse transcription PCR on RNA from minor salivary glands of patients with SS confirmed a positive correlation between BMP6 and LAMP3 expression. BMP6 expression could be experimentally induced in mice by overexpression of LAMP3, which developed an SS-like phenotype. The newly identified LAMP3/HSP70/BMP6 axis provided an etiological model for SS gland dysfunction and autoimmunity.

Sclerosing Sialadenitis Is Associated With Salivary Gland Hypofunction and a Unique Gene Expression Profile in Sjögren's Syndrome.

Purpose: To develop a novel method to quantify the amount of fibrosis in the salivary gland and to investigate the relationship between fibrosis and specific symptoms associated with Sjögren's syndrome (SS) using this method. Materials and Methods: Paraffin-embedded labial salivary gland (LSG) slides from 20 female SS patients and their clinical and LSG pathology data were obtained from the Sjögren's International Collaborative Clinical Alliance. Relative interstitial fibrosis area (RIFA) in Masson's trichrome-stained LSG sections was quantified from digitally scanned slides and used for correlation analysis. Gene expression levels were assessed by microarray analysis. Core promoter accessibility for RIFA-correlated genes was determined using DNase I hypersensitive sites sequencing analysis. Results: RIFA was significantly correlated with unstimulated whole saliva flow rate in SS patients. Sixteen genes were significantly and positively correlated with RIFA. In a separate analysis, a group of differentially expressed genes was identified by comparing severe and moderate fibrosis groups. This combined set of genes was distinct from differentially expressed genes identified in lung epithelium from idiopathic pulmonary fibrosis patients compared with controls. Single-cell RNA sequencing analysis of salivary glands suggested most of the RIFA-correlated genes are expressed by fibroblasts in the gland and are in a permissive chromatin state. Conclusion: RIFA quantification is a novel method for assessing interstitial fibrosis and the impact of fibrosis on SS symptoms. Loss of gland function may be associated with salivary gland fibrosis, which is likely to be driven by a unique set of genes that are mainly expressed by fibroblasts.

Lysosome-associated membrane protein 3 misexpression in salivary glands induces a Sjögren's syndrome-like phenotype in mice.

OBJECTIVES: Sjögren's syndrome (SS) is an autoimmune sialadenitis with unknown aetiology. Although extensive research implicated an abnormal immune response associated with lymphocytes, an initiating event mediated by salivary gland epithelial cell (SGEC) abnormalities causing activation is poorly characterised. Transcriptome studies have suggested alternations in lysosomal function are associated with SS, but a cause and effect linkage has not been established. In this study, we demonstrated that altered lysosome activity in SGECs by expression of lysosome-associated membrane protein 3 (LAMP3) can initiate an autoimmune response with autoantibody production and salivary dysfunction similar to SS. METHODS: Retroductal cannulation of the submandibular salivary glands with an adeno-associated virus serotype 2 vector encoding LAMP3 was used to establish a model system. Pilocarpine-stimulated salivary flow and the presence of autoantibodies were assessed at several time points post-cannulation. Salivary glands from the mice were evaluated using RNAseq and histologically. RESULTS: Following LAMP3 expression, saliva flow was significantly decreased and serum anti-Ro/SSA and La/SSB antibodies could be detected in the treated mice. Mechanistically, LAMP3 expression increased apoptosis in SGECs and decreased protein expression related to saliva secretion. Analysis of RNAseq data suggested altered lysosomal function in the transduced SGECs, and that the cellular changes can chemoattract immune cells into the salivary glands. Immune cells were activated via toll-like receptors by damage-associated molecular patterns released from LAMP3-expressing SGECs. CONCLUSIONS: These results show a critical role for lysosomal trafficking in the development of SS and establish a causal relationship between LAMP3 misexpression and the development of SS.

SARS-CoV-2 infection of the oral cavity and saliva.

Despite signs of infection-including taste loss, dry mouth and mucosal lesions such as ulcerations, enanthema and macules-the involvement of the oral cavity in coronavirus disease 2019 (COVID-19) is poorly understood. To address this, we generated and analyzed two single-cell RNA sequencing datasets of the human minor salivary glands and gingiva (9 samples, 13,824 cells), identifying 50 cell clusters. Using integrated cell normalization and annotation, we classified 34 unique cell subpopulations between glands and gingiva. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) viral entry factors such as ACE2 and TMPRSS members were broadly enriched in epithelial cells of the glands and oral mucosae. Using orthogonal RNA and protein expression assessments, we confirmed SARS-CoV-2 infection in the glands and mucosae. Saliva from SARS-CoV-2-infected individuals harbored epithelial cells exhibiting ACE2 and TMPRSS expression and sustained SARS-CoV-2 infection. Acellular and cellular salivary fractions from asymptomatic individuals were found to transmit SARS-CoV-2 ex vivo. Matched nasopharyngeal and saliva samples displayed distinct viral shedding dynamics, and salivary viral burden correlated with COVID-19 symptoms, including taste loss. Upon recovery, this asymptomatic cohort exhibited sustained salivary IgG antibodies against SARS-CoV-2. Collectively, these data show that the oral cavity is an important site for SARS-CoV-2 infection and implicate saliva as a potential route of SARS-CoV-2 transmission.

Integrated Single-Cell Atlases Reveal an Oral SARS-CoV-2 Infection and Transmission Axis.

Despite signs of infection, the involvement of the oral cavity in COVID-19 is poorly understood. To address this, single-cell RNA sequencing data-sets were integrated from human minor salivary glands and gingiva to identify 11 epithelial, 7 mesenchymal, and 15 immune cell clusters. Analysis of SARS-CoV-2 viral entry factor expression showed enrichment in epithelia including the ducts and acini of the salivary glands and the suprabasal cells of the mucosae. COVID-19 autopsy tissues confirmed in vivo SARS-CoV-2 infection in the salivary glands and mucosa. Saliva from SARS-CoV-2-infected individuals harbored epithelial cells exhibiting ACE2 expression and SARS-CoV-2 RNA. Matched nasopharyngeal and saliva samples found distinct viral shedding dynamics and viral burden in saliva correlated with COVID-19 symptoms including taste loss. Upon recovery, this cohort exhibited salivary antibodies against SARS-CoV-2 proteins. Collectively, the oral cavity represents a robust site for COVID-19 infection and implicates saliva in viral transmission.

Comparison of Prostate-Specific Membrane Antigen Expression Levels in Human Salivary Glands to Non-Human Primates and Rodents.

Background: Prostate-specific membrane antigen (PSMA) has emerged as a promising target for developing radionuclide therapy (RNT) in prostate cancer; however, accumulation of PSMA-RNT in salivary glands can result in irreversible xerostomia. Methods to prevent PSMA-RNT-related xerostomia could be clinically useful; however, little is known about PSMA expression in salivary glands of preclinical animal models. Using [18F]DCFPyL autoradiography/biodistribution, PSMA expression levels were determined in salivary glands of various preclinical monkey and rodent species and compared with humans. Methods: Binding affinities (Kd) and PSMA levels (Bmax) were determined by in vitro [18F]DCFPyL autoradiography studies. In vivo rodent tissue uptakes (%ID/g) were determined from [18F]DCFPyL biodistributions. Results: [18F]DCFPyL exhibited low nanomolar Kd for submandibular gland (SMG) PSMA across all the species. PSMA levels in human SMG (Bmax = 60.91 nM) were approximately two-fold lower compared with baboon SMG but were two- to three-fold higher than SMG PSMA levels of cynomolgus and rhesus. Rodents had the lowest SMG PSMA levels, with the mouse being 10-fold higher than the rat. In vivo rodent biodistribution studies confirmed these results. Conclusions: SMG of monkeys exhibited comparable PSMA expression to human SMG whereas rodents were lower. However, the results suggest that mice are relatively a better small animal preclinical model than rats for PSMA salivary gland studies.

Identification of RNA aptamer which specifically interacts with PtdIns(3)P.

The phosphinositide PtdIns(3)P plays an important role in autophagy; however, the detailed mechanism of its activity remains unclear. Here, we used a Systematic Evolution of Ligands by EXponential enrichment (SELEX) screening approach to identify an RNA aptamer of 40 nucleotides that specifically recognizes and binds to intracellular lysosomal PtdIns(3)P. Binding occurs in a magnesium concentration- and pH-dependent manner, and consequently inhibits autophagy as determined by LC3II/I conversion, p62 degradation, formation of LC3 puncta, and lysosomal accumulation of Phafin2. These effects in turn inhibited lysosomal acidification, and the subsequent hydrolytic activity of cathepsin D following induction of autophagy. Given the essential role of PtdIns(3)P as a key targeting molecule for autophagy induction, identification of this novel PtdIns(3)P RNA aptamer provides new opportunities for investigating the biological functions and mechanisms of phosphoinositides.

Correction to: ATAC2GRN: optimized ATAC-seq and DNase1-seq pipelines for rapid and accurate genome regulatory network inference.

Following the publication of this article [1], the authors informed us of the following typographical errors in the Results section (the changes are marked in bold).

ATAC2GRN: optimized ATAC-seq and DNase1-seq pipelines for rapid and accurate genome regulatory network inference.

BACKGROUND: Chromatin accessibility profiling assays such as ATAC-seq and DNase1-seq offer the opportunity to rapidly characterize the regulatory state of the genome at a single nucleotide resolution. Optimization of molecular protocols has enabled the molecular biologist to produce next-generation sequencing libraries in several hours, leaving the analysis of sequencing data as the primary obstacle to wide-scale deployment of accessibility profiling assays. To address this obstacle we have developed an optimized and efficient pipeline for the analysis of ATAC-seq and DNase1-seq data. RESULTS: We executed a multi-dimensional grid-search on the NIH Biowulf supercomputing cluster to assess the impact of parameter selection on biological reproducibility and ChIP-seq recovery by analyzing 4560 pipeline configurations. Our analysis improved ChIP-seq recovery by 15% for ATAC-seq and 3% for DNase1-seq and determined that PCR duplicate removal improves biological reproducibility by 36% without significant costs in footprinting transcription factors. Our analyses of down sampled reads identified a point of diminishing returns for increased library sequencing depth, with 95% of the ChIP-seq data of a 200 million read footprinting library recovered by 160 million reads. CONCLUSIONS: We present optimized ATAC-seq and DNase-seq pipelines in both Snakemake and bash formats as well as optimal sequencing depths for ATAC-seq and DNase-seq projects. The optimized ATAC-seq and DNase1-seq analysis pipelines, parameters, and ground-truth ChIP-seq datasets have been made available for deployment and future algorithmic profiling.

Cross-contamination of the human salivary gland HSG cell line with HeLa cells: A STR analysis study.

OBJECTIVES: The human salivary gland (HSG) cell line, labeled as a submandibular ductal cell line, is commonly used as in vitro models to study radiation therapy, Sjögren's syndrome, pleomorphic adenoma, mucocele, epithelial-to-mesenchymal transition, and epigenetics. However, the American Type Culture Collection (ATCC) has recently released a list of cross-contaminated cell lines that included HSG. Despite this notice, some research laboratories still use HSG as a salivary cell model. Therefore, this study examined the authenticity of HSG sampled from three different laboratories. METHODS: DNA was extracted from HSG and additional salivary cell lines (NS-SV-AC, NS-SV-DC, A253, HSY) and submitted for cell line authentication with short tandem repeat (STR) analysis. RESULTS: All HSG samples had STR profiles indicating >80% match with HeLa in both the ATCC and Deutsche Sammlung von Mikroorganismen und Zellkulturen (DSMZ) databases. This confirmed that HSG sampled from three different laboratories and HSY shared a common ancestry (host) with HeLa, whereas NS-SV-AC, NS-SV-DC, and A253 had unique STR profiles. CONCLUSION: Short tandem repeat analysis revealed that HSG was contaminated by the HeLa cell line. Furthermore, because genotyping of the original HSG cell line was not performed during its establishment, it will be difficult to authenticate an uncontaminated sample of HSG.