Stromal Signals Dominate Gene Expression Signature Scores That Aim to Describe Cancer Cell-intrinsic Stemness or Mesenchymality Characteristics.

Epithelial-to-mesenchymal transition (EMT) in cancer cells confers migratory abilities, a crucial aspect in the metastasis of tumors that frequently leads to death. In multiple studies, authors proposed gene expression signatures for EMT, stemness, or mesenchymality of tumors based on bulk tumor expression profiling. However, recent studies suggested that noncancerous cells from the microenvironment or macroenvironment heavily influence such signature profiles. Here, we strengthen these findings by investigating 11 published and frequently referenced gene expression signatures that were proposed to describe EMT-related (EMT, mesenchymal, or stemness) characteristics in various cancer types. By analyses of bulk, single-cell, and pseudobulk expression data, we show that the cell type composition of a tumor sample frequently dominates scores of these EMT-related signatures. A comprehensive, integrated analysis of bulk RNA sequencing (RNA-seq) and single-cell RNA-seq data shows that stromal cells, most often fibroblasts, are the main drivers of EMT-related signature scores. We call attention to the risk of false conclusions about tumor properties when interpreting EMT-related signatures, especially in a clinical setting: high patient scores of EMT-related signatures or calls of "stemness subtypes" often result from low cancer cell content in tumor biopsies rather than cancer cell-specific stemness or mesenchymal/EMT characteristics. SIGNIFICANCE: Cancer self-renewal and migratory abilities are often characterized via gene module expression profiles, also called EMT or stemness gene expression signatures. Using published clinical tumor samples, cancer cell lines, and single cancer cells, we highlight the dominating influence of noncancer cells in low cancer cell content biopsies on their scores. We caution on their application for low cancer cell content clinical cancer samples with the intent to assign such characteristics or subtypes.

The effector program of human CD8 T cells supports tissue remodeling.

CD8 T lymphocytes are classically viewed as cytotoxic T cells. Whether human CD8 T cells can, in parallel, induce a tissue regeneration program is poorly understood. Here, antigen-specific assay systems revealed that human CD8 T cells not only mediated cytotoxicity but also promoted tissue remodeling. Activated CD8 T cells could produce the epidermal growth factor receptor (EGFR)-ligand amphiregulin (AREG) and sensitize epithelial cells for enhanced regeneration potential. Blocking the EGFR or the effector cytokines IFN-γ and TNF could inhibit tissue remodeling. This regenerative program enhanced tumor spheroid and stem cell-mediated organoid growth. Using single-cell gene expression analysis, we identified an AREG+, tissue-resident CD8 T cell population in skin and adipose tissue from patients undergoing abdominal wall or abdominoplasty surgery. These tissue-resident CD8 T cells showed a strong TCR clonal relation to blood PD1+TIGIT+ CD8 T cells with tissue remodeling abilities. These findings may help to understand the complex CD8 biology in tumors and could become relevant for the design of therapeutic T cell products.

Explainable artificial intelligence for omics data: a systematic mapping study.

Researchers increasingly turn to explainable artificial intelligence (XAI) to analyze omics data and gain insights into the underlying biological processes. Yet, given the interdisciplinary nature of the field, many findings have only been shared in their respective research community. An overview of XAI for omics data is needed to highlight promising approaches and help detect common issues. Toward this end, we conducted a systematic mapping study. To identify relevant literature, we queried Scopus, PubMed, Web of Science, BioRxiv, MedRxiv and arXiv. Based on keywording, we developed a coding scheme with 10 facets regarding the studies' AI methods, explainability methods and omics data. Our mapping study resulted in 405 included papers published between 2010 and 2023. The inspected papers analyze DNA-based (mostly genomic), transcriptomic, proteomic or metabolomic data by means of neural networks, tree-based methods, statistical methods and further AI methods. The preferred post-hoc explainability methods are feature relevance (n = 166) and visual explanation (n = 52), while papers using interpretable approaches often resort to the use of transparent models (n = 83) or architecture modifications (n = 72). With many research gaps still apparent for XAI for omics data, we deduced eight research directions and discuss their potential for the field. We also provide exemplary research questions for each direction. Many problems with the adoption of XAI for omics data in clinical practice are yet to be resolved. This systematic mapping study outlines extant research on the topic and provides research directions for researchers and practitioners.

Using single-cell chromatin accessibility sequencing to characterize CD4+ T cells from murine tissues.

The Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq) is a cutting-edge technology that enables researchers to assess genome-wide chromatin accessibility and to characterize cell type specific gene-regulatory programs. Recent technological progress allows for using this technology also on the single-cell level. In this article, we describe the whole value chain from the isolation of T cells from murine tissues to a complete bioinformatic analysis workflow. We start with methods for isolating scATAC-seq-ready CD4+ T cells from murine tissues such as visceral adipose tissue, skin, colon, and secondary lymphoid tissues such as the spleen. We describe the preparation of nuclei and quality control parameters during library preparation. Based on publicly available sequencing data that was generated using these protocols, we describe a step-by-step bioinformatic analysis pipeline for data pre-processing and downstream analysis. Our analysis workflow will follow the R-based bioinformatics framework ArchR, which is currently well established for scATAC-seq datasets. All in all, this work serves as a one-stop shop for generating and analyzing chromatin accessibility landscapes in T cells.

Chronic chromosome instability induced by Plk1 results in immune suppression in breast cancer.

Chromosome instability (CIN) contributes to resistance to therapies and tumor evolution. Although natural killer (NK) cells can eliminate cells with complex karyotypes, high-CIN human tumors have an immunosuppressive phenotype. To understand which CIN-associated molecular features alter immune recognition during tumor evolution, we overexpress Polo-like kinase 1 (Plk1) in a Her2+ breast cancer model. These high-CIN tumors activate a senescence-associated secretory phenotype (SASP), upregulate PD-L1 and CD206, and induce non-cell-autonomous nuclear factor κB (NF-κß) signaling, facilitating immune evasion. Single-cell RNA sequencing from pre-neoplastic mammary glands unveiled the presence of Arg1+ macrophages, NK cells with reduced effector functions, and increased resting regulatory T cell infiltration. We further show that high PLK1-expressing human breast tumors display gene expression patterns associated with SASP, NF-κß signaling, and immune suppression. These findings underscore the need to understand the immune landscape in CIN tumors to identify more effective therapies, potentially combining immune checkpoint or NF-κß inhibitors with current treatments.

DNA-methylation signature accurately differentiates pancreatic cancer from chronic pancreatitis in tissue and plasma.

OBJECTIVE: Pancreatic ductal adenocarcinoma (PDAC) is a lethal malignancy. Differentiation from chronic pancreatitis (CP) is currently inaccurate in about one-third of cases. Misdiagnoses in both directions, however, have severe consequences for patients. We set out to identify molecular markers for a clear distinction between PDAC and CP. DESIGN: Genome-wide variations of DNA-methylation, messenger RNA and microRNA level as well as combinations thereof were analysed in 345 tissue samples for marker identification. To improve diagnostic performance, we established a random-forest machine-learning approach. Results were validated on another 48 samples and further corroborated in 16 liquid biopsy samples. RESULTS: Machine-learning succeeded in defining markers to differentiate between patients with PDAC and CP, while low-dimensional embedding and cluster analysis failed to do so. DNA-methylation yielded the best diagnostic accuracy by far, dwarfing the importance of transcript levels. Identified changes were confirmed with data taken from public repositories and validated in independent sample sets. A signature of six DNA-methylation sites in a CpG-island of the protein kinase C beta type gene achieved a validated diagnostic accuracy of 100% in tissue and in circulating free DNA isolated from patient plasma. CONCLUSION: The success of machine-learning to identify an effective marker signature documents the power of this approach. The high diagnostic accuracy of discriminating PDAC from CP could have tremendous consequences for treatment success, once the result from still a limited number of liquid biopsy samples would be confirmed in a larger cohort of patients with suspected pancreatic cancer.

Persister cell phenotypes contribute to poor patient outcomes after neoadjuvant chemotherapy in PDAC.

Neoadjuvant chemotherapy can improve the survival of individuals with borderline and unresectable pancreatic ductal adenocarcinoma; however, heterogeneous responses to chemotherapy remain a significant clinical challenge. Here, we performed RNA sequencing (n = 97) and multiplexed immunofluorescence (n = 122) on chemo-naive and postchemotherapy (post-CTX) resected patient samples (chemoradiotherapy excluded) to define the impact of neoadjuvant chemotherapy. Transcriptome analysis combined with high-resolution mapping of whole-tissue sections identified GATA6 (classical), KRT17 (basal-like) and cytochrome P450 3A (CYP3A) coexpressing cells that were preferentially enriched in post-CTX resected samples. The persistence of GATA6hi and KRT17hi cells post-CTX was significantly associated with poor survival after mFOLFIRINOX (mFFX), but not gemcitabine (GEM), treatment. Analysis of organoid models derived from chemo-naive and post-CTX samples demonstrated that CYP3A expression is a predictor of chemotherapy response and that CYP3A-expressing drug detoxification pathways can metabolize the prodrug irinotecan, a constituent of mFFX. These findings identify CYP3A-expressing drug-tolerant cell phenotypes in residual disease that may ultimately inform adjuvant treatment selection.

Immune cell type signature discovery and random forest classification for analysis of single cell gene expression datasets.

Background: Robust immune cell gene expression signatures are central to the analysis of single cell studies. Nearly all known sets of immune cell signatures have been derived by making use of only single gene expression datasets. Utilizing the power of multiple integrated datasets could lead to high-quality immune cell signatures which could be used as superior inputs to machine learning-based cell type classification approaches. Results: We established a novel workflow for the discovery of immune cell type signatures based primarily on gene-versus-gene expression similarity. It leverages multiple datasets, here seven single cell expression datasets from six different cancer types and resulted in eleven immune cell type-specific gene expression signatures. We used these to train random forest classifiers for immune cell type assignment for single-cell RNA-seq datasets. We obtained similar or better prediction results compared to commonly used methods for cell type assignment in independent benchmarking datasets. Our gene signature set yields higher prediction scores than other published immune cell type gene sets in random forest-based cell type classification. We further demonstrate how our approach helps to avoid bias in downstream statistical analyses by re-analysis of a published IFN stimulation experiment. Discussion and conclusion: We demonstrated the quality of our immune cell signatures and their strong performance in a random forest-based cell typing approach. We argue that classifying cells based on our comparably slim sets of genes accompanied by a random forest-based approach not only matches or outperforms widely used published approaches. It also facilitates unbiased downstream statistical analyses of differential gene expression between cell types for significantly more genes compared to previous cell classification algorithms.

Single-cell transcriptomics and data analyses for prokaryotes-Past, present and future concepts.

Transcriptomics, or more specifically mRNA sequencing, is a powerful tool to study gene expression at the single-cell level (scRNA-seq) which enables new insights into a plethora of biological processes. While methods for single-cell RNA-seq in eukaryotes are well established, application to prokaryotes is still challenging. Reasons for that are rigid and diverse cell wall structures hampering lysis, the lack of polyadenylated transcripts impeding mRNA enrichment, and minute amounts of RNA requiring amplification steps before sequencing. Despite those obstacles, several promising scRNA-seq approaches for bacteria have been published recently, albeit difficulties in the experimental workflow and data processing and analysis remain. In particular, bias is often introduced by amplification which makes it difficult to distinguish between technical noise and biological variation. Future optimization of experimental procedures and data analysis algorithms are needed for the improvement of scRNA-seq but also to aid in the emergence of prokaryotic single-cell multi-omics. to help address 21st century challenges in the biotechnology and health sector.

Integrated genetic analyses of immunodeficiency-associated Epstein-Barr virus- (EBV) positive primary CNS lymphomas.

Immunodeficiency-associated primary CNS lymphoma (PCNSL) represents a distinct clinicopathological entity, which is typically Epstein-Barr virus-positive (EBV+) and carries an inferior prognosis. Genetic alterations that characterize EBV-related CNS lymphomagenesis remain unclear precluding molecular classification and targeted therapies. In this study, a comprehensive genetic analysis of 22 EBV+ PCNSL, therefore, integrated clinical and pathological information with exome and RNA sequencing (RNASeq) data. EBV+ PCNSL with germline controls carried a median of 55 protein-coding single nucleotide variants (SNVs; range 24-217) and 2 insertions/deletions (range 0-22). Genetic landscape was largely shaped by aberrant somatic hypermutation with a median of 41.01% (range 31.79-53.49%) of SNVs mapping to its target motifs. Tumors lacked established SNVs (MYD88, CD79B, PIM1) and copy number variants (CDKN2A, HLA loss) driving EBV- PCNSL. Instead, EBV+ PCNSL were characterized by SOCS1 mutations (26%), predicted to disinhibit JAK/STAT signaling, and mutually exclusive gain-of-function NOTCH pathway SNVs (26%). Copy number gains were enriched on 11q23.3, a locus directly targeted for chromosomal aberrations by EBV, that includes SIK3 known to protect from cytotoxic T-cell responses. Losses covered 5q31.2 (STING), critical for sensing viral DNA, and 17q11 (NF1). Unsupervised clustering of RNASeq data revealed two distinct transcriptional groups, that shared strong expression of CD70 and IL1R2, previously linked to tolerogenic tumor microenvironments. Correspondingly, deconvolution of bulk RNASeq data revealed elevated M2-macrophage, T-regulatory cell, mast cell and monocyte fractions in EBV+ PCNSL. In addition to novel insights into the pathobiology of EBV+ PCNSL, the data provide the rationale for the exploration of targeted therapies including JAK-, NOTCH- and CD70-directed approaches.

Spatially Resolved Multi-Omics Single-Cell Analyses Inform Mechanisms of Immune Dysfunction in Pancreatic Cancer.

BACKGROUND & AIMS: As pancreatic ductal adenocarcinoma (PDAC) continues to be recalcitrant to therapeutic interventions, including poor response to immunotherapy, albeit effective in other solid malignancies, a more nuanced understanding of the immune microenvironment in PDAC is urgently needed. We aimed to unveil a detailed view of the immune micromilieu in PDAC using a spatially resolved multimodal single-cell approach. METHODS: We applied single-cell RNA sequencing, spatial transcriptomics, multiplex immunohistochemistry, and mass cytometry to profile the immune compartment in treatment-naïve PDAC tumors and matched adjacent normal pancreatic tissue, as well as in the systemic circulation. We determined prognostic associations of immune signatures and performed a meta-analysis of the immune microenvironment in PDAC and lung adenocarcinoma on single-cell level. RESULTS: We provided a spatially resolved fine map of the immune landscape in PDAC. We substantiated the exhausted phenotype of CD8 T cells and immunosuppressive features of myeloid cells, and highlighted immune subsets with potentially underappreciated roles in PDAC that diverged from immune populations within adjacent normal areas, particularly CD4 T cell subsets and natural killer T cells that are terminally exhausted and acquire a regulatory phenotype. Differential analysis of immune phenotypes in PDAC and lung adenocarcinoma revealed the presence of extraordinarily immunosuppressive subtypes in PDAC, along with a distinctive immune checkpoint composition. CONCLUSIONS: Our study sheds light on the multilayered immune dysfunction in PDAC and presents a holistic view of the immune landscape in PDAC and lung adenocarcinoma, providing a comprehensive resource for functional studies and the exploration of therapeutically actionable targets in PDAC.

High-resolution transcriptomic and epigenetic profiling identifies novel regulators of COPD.

Patients with chronic obstructive pulmonary disease (COPD) are still waiting for curative treatments. Considering its environmental cause, we hypothesized that COPD will be associated with altered epigenetic signaling in lung cells. We generated genome-wide DNA methylation maps at single CpG resolution of primary human lung fibroblasts (HLFs) across COPD stages. We show that the epigenetic landscape is changed early in COPD, with DNA methylation changes occurring predominantly in regulatory regions. RNA sequencing of matched fibroblasts demonstrated dysregulation of genes involved in proliferation, DNA repair, and extracellular matrix organization. Data integration identified 110 candidate regulators of disease phenotypes that were linked to fibroblast repair processes using phenotypic screens. Our study provides high-resolution multi-omic maps of HLFs across COPD stages. We reveal novel transcriptomic and epigenetic signatures associated with COPD onset and progression and identify new candidate regulators involved in the pathogenesis of chronic lung diseases. The presence of various epigenetic factors among the candidates demonstrates that epigenetic regulation in COPD is an exciting research field that holds promise for novel therapeutic avenues for patients.

Corrigendum: Clinical impact of molecular subtyping of pancreatic cancer.

[This corrects the article DOI: 10.3389/fcell.2021.743908.].

Neuroblastoma arises in early fetal development and its evolutionary duration predicts outcome.

Neuroblastoma, the most frequent solid tumor in infants, shows very diverse outcomes from spontaneous regression to fatal disease. When these different tumors originate and how they evolve are not known. Here we quantify the somatic evolution of neuroblastoma by deep whole-genome sequencing, molecular clock analysis and population-genetic modeling in a comprehensive cohort covering all subtypes. We find that tumors across the entire clinical spectrum begin to develop via aberrant mitoses as early as the first trimester of pregnancy. Neuroblastomas with favorable prognosis expand clonally after short evolution, whereas aggressive neuroblastomas show prolonged evolution during which they acquire telomere maintenance mechanisms. The initial aneuploidization events condition subsequent evolution, with aggressive neuroblastoma exhibiting early genomic instability. We find in the discovery cohort (n = 100), and validate in an independent cohort (n = 86), that the duration of evolution is an accurate predictor of outcome. Thus, insight into neuroblastoma evolution may prospectively guide treatment decisions.

Integrated single-cell profiling dissects cell-state-specific enhancer landscapes of human tumor-infiltrating CD8+ T cells.

Despite extensive studies on the chromatin landscape of exhausted T cells, the transcriptional wiring underlying the heterogeneous functional and dysfunctional states of human tumor-infiltrating lymphocytes (TILs) is incompletely understood. Here, we identify gene-regulatory landscapes in a wide breadth of functional and dysfunctional CD8+ TIL states covering four cancer entities using single-cell chromatin profiling. We map enhancer-promoter interactions in human TILs by integrating single-cell chromatin accessibility with single-cell RNA-seq data from tumor-entity-matching samples and prioritize cell-state-specific genes by super-enhancer analysis. Besides revealing entity-specific chromatin remodeling in exhausted TILs, our analyses identify a common chromatin trajectory to TIL dysfunction and determine key enhancers, transcriptional regulators, and deregulated genes involved in this process. Finally, we validate enhancer regulation at immunotherapeutically relevant loci by targeting non-coding regulatory elements with potent CRISPR activators and repressors. In summary, our study provides a framework for understanding and manipulating cell-state-specific gene-regulatory cues from human tumor-infiltrating lymphocytes.

Feasibility and outcome of reproducible clinical interpretation of high-dimensional molecular data: a comparison of two molecular tumor boards.

BACKGROUND: Structured and harmonized implementation of molecular tumor boards (MTB) for the clinical interpretation of molecular data presents a current challenge for precision oncology. Heterogeneity in the interpretation of molecular data was shown for patients even with a limited number of molecular alterations. Integration of high-dimensional molecular data, including RNA- (RNA-Seq) and whole-exome sequencing (WES), is expected to further complicate clinical application. To analyze challenges for MTB harmonization based on complex molecular datasets, we retrospectively compared clinical interpretation of WES and RNA-Seq data by two independent molecular tumor boards. METHODS: High-dimensional molecular cancer profiling including WES and RNA-Seq was performed for patients with advanced solid tumors, no available standard therapy, ECOG performance status of 0-1, and available fresh-frozen tissue within the DKTK-MASTER Program from 2016 to 2018. Identical molecular profiling data of 40 patients were independently discussed by two molecular tumor boards (MTB) after prior annotation by specialized physicians, following independent, but similar workflows. Identified biomarkers and resulting treatment options were compared between the MTBs and patients were followed up clinically. RESULTS: A median of 309 molecular aberrations from WES and RNA-Seq (n = 38) and 82 molecular aberrations from WES only (n = 3) were considered for clinical interpretation for 40 patients (one patient sequenced twice). A median of 3 and 2 targeted treatment options were identified per patient, respectively. Most treatment options were identified for receptor tyrosine kinase, PARP, and mTOR inhibitors, as well as immunotherapy. The mean overlap coefficient between both MTB was 66%. Highest agreement rates were observed with the interpretation of single nucleotide variants, clinical evidence levels 1 and 2, and monotherapy whereas the interpretation of gene expression changes, preclinical evidence levels 3 and 4, and combination therapy yielded lower agreement rates. Patients receiving treatment following concordant MTB recommendations had significantly longer overall survival than patients receiving treatment following discrepant recommendations or physician's choice. CONCLUSIONS: Reproducible clinical interpretation of high-dimensional molecular data is feasible and agreement rates are encouraging, when compared to previous reports. The interpretation of molecular aberrations beyond single nucleotide variants and preclinically validated biomarkers as well as combination therapies were identified as additional difficulties for ongoing harmonization efforts.

Bcl-xL as prognostic marker and potential therapeutic target in cholangiocarcinoma.

Intrahepatic, perihilar, and distal cholangiocarcinoma (iCCA, pCCA, dCCA) are highly malignant tumours with increasing mortality rates due to therapy resistances. Among the mechanisms mediating resistance, overexpression of anti-apoptotic Bcl-2 proteins (Bcl-2, Bcl-xL , Mcl-1) is particularly important. In this study, we investigated whether antiapoptotic protein patterns are prognostically relevant and potential therapeutic targets in CCA. Bcl-2 proteins were analysed in a pan-cancer cohort from the NCT/DKFZ/DKTK MASTER registry trial (n = 1140, CCA n = 72) via RNA-sequencing and transcriptome-based protein activity interference revealing high ranks of CCA for Bcl-xL and Mcl-1. Expression of Bcl-xL , Mcl-1, and Bcl-2 was assessed in human CCA tissue and cell lines compared with cholangiocytes by immunohistochemistry, immunoblotting, and quantitative-RT-PCR. Immunohistochemistry confirmed the upregulation of Bcl-xL and Mcl-1 in iCCA tissues. Cell death of CCA cell lines upon treatment with specific small molecule inhibitors of Bcl-xL (Wehi-539), of Mcl-1 (S63845), and Bcl-2 (ABT-199), either alone, in combination with each other or together with chemotherapeutics was assessed by flow cytometry. Targeting Bcl-xL induced cell death and augmented the effect of chemotherapy in CCA cells. Combined inhibition of Bcl-xL and Mcl-1 led to a synergistic increase in cell death in CCA cell lines. Correlation between Bcl-2 protein expression and survival was analysed within three independent patient cohorts from cancer centers in Germany comprising 656 CCA cases indicating a prognostic value of Bcl-xL in CCA depending on the CCA subtype. Collectively, these observations identify Bcl-xL as a key protein in cell death resistance of CCA and may pave the way for clinical application.

The antileukemic activity of decitabine upon PML/RARA-negative AML blasts is supported by all-trans retinoic acid: in vitro and in vivo evidence for cooperation.

The prognosis of AML patients with adverse genetics, such as a complex, monosomal karyotype and TP53 lesions, is still dismal even with standard chemotherapy. DNA-hypomethylating agent monotherapy induces an encouraging response rate in these patients. When combined with decitabine (DAC), all-trans retinoic acid (ATRA) resulted in an improved response rate and longer overall survival in a randomized phase II trial (DECIDER; NCT00867672). The molecular mechanisms governing this in vivo synergism are unclear. We now demonstrate cooperative antileukemic effects of DAC and ATRA on AML cell lines U937 and MOLM-13. By RNA-sequencing, derepression of >1200 commonly regulated transcripts following the dual treatment was observed. Overall chromatin accessibility (interrogated by ATAC-seq) and, in particular, at motifs of retinoic acid response elements were affected by both single-agent DAC and ATRA, and enhanced by the dual treatment. Cooperativity regarding transcriptional induction and chromatin remodeling was demonstrated by interrogating the HIC1, CYP26A1, GBP4, and LYZ genes, in vivo gene derepression by expression studies on peripheral blood blasts from AML patients receiving DAC + ATRA. The two drugs also cooperated in derepression of transposable elements, more effectively in U937 (mutated TP53) than MOLM-13 (intact TP53), resulting in a "viral mimicry" response. In conclusion, we demonstrate that in vitro and in vivo, the antileukemic and gene-derepressive epigenetic activity of DAC is enhanced by ATRA.

Transcriptomic diversity in human medullary thymic epithelial cells.

The induction of central T cell tolerance in the thymus depends on the presentation of peripheral self-epitopes by medullary thymic epithelial cells (mTECs). This promiscuous gene expression (pGE) drives mTEC transcriptomic diversity, with non-canonical transcript initiation, alternative splicing, and expression of endogenous retroelements (EREs) representing important but incompletely understood contributors. Here we map the expression of genome-wide transcripts in immature and mature human mTECs using high-throughput 5' cap and RNA sequencing. Both mTEC populations show high splicing entropy, potentially driven by the expression of peripheral splicing factors. During mTEC maturation, rates of global transcript mis-initiation increase and EREs enriched in long terminal repeat retrotransposons are up-regulated, the latter often found in proximity to differentially expressed genes. As a resource, we provide an interactive public interface for exploring mTEC transcriptomic diversity. Our findings therefore help construct a map of transcriptomic diversity in the healthy human thymus and may ultimately facilitate the identification of those epitopes which contribute to autoimmunity and immune recognition of tumor antigens.