3781 results for "RNA AND sequencing OR transcriptome"

in 19.62 milliseconds.

• Characterizing the Role of the Immune Microenvironment in Multiple Myeloma Progression at a Single Cell Level

  The plasma cell disease multiple myeloma (MM) is highly dependent on its bone marrow microenvironment (BMME). Early alterations within the bone marrow microenvironment that contribute to the progression of multiple myeloma (MM) from its precursor stages could hold the key to identifying novel therapeutic approaches, yet the intrinsic variability in cellular populations between patients together with differences in sample processing and analysis methods have made it difficult to identify consistent changes between datasets. Here, we used single-cell RNA sequencing of BM cells from precursor stages, monoclonal gammopathy of unknown significance (MGUS), smoldering MM (SMM) and newly diagnosed MM and analyzed our data in combination with a previously published dataset that used a similar patient population and sample processing. Despite vast interpatient heterogeneity, some alterations were seen consistently in both datasets. We identified changes in immune cell populations as disease progressed that were characterized by a substantial decrease in memory and naïve CD4 T cells and an increase in CD8+ effector T cells and T-regulatory cells. These alterations were further accompanied by an enrichment of non-clonal memory B cells and an increase of CD14 and CD16 monocytes in MM compared to its precursor stages. These results provide crucial information of the immune changes associated with the progression to clinical MM and can help to develop immune based strategies for patient stratification and early therapeutic intervention. 

  Study phs002756

• Multi-omics bulk and single-cell profiling of epithelioid sarcoma

  Epithelioid sarcoma (EpS) is an aggressive sarcoma, characterized by the loss of SMARCB1 expression. EpS is traditionally classified as distal or proximal according to clinicopathological features, but its molecular characteristics remain largely unknown. To establish an EpS molecular classification and uncover determinants of inter- and intra-patient heterogeneity in EpS, we used multi-omics profiling and integrated the genomic, transcriptional and methylome landscapes with single-cell RNA sequencing on fresh samples as well as spatial transcriptomics. We identified two molecular subtypes of EpS: “distal-like” and “proximal-like”, which were distinct from the histological subtypes. Distal-like tumors include all morphological distal subtypes and a subset of proximal and hybrid EpS. They express a specific molecular single-cell derived epithelial-to-mesenchymal transition signature, which associates with improved patient survival and includes Desmoglein 2 (DSG2), which we identify as a potential routine immunohistochemical biomarker to improve diagnosis. Distal-like EpS also display increased peri-tumoral CD8+ T cell infiltrates and specific tumor-immune cell interactions. Conversely, proximal-like EpS harbor a higher inter-tumoral molecular heterogeneity with some cases resembling other SMARCB1-deficient tumors by DNA methylation profiling, and increased intra-tumoral pro-tumoral macrophage infiltrates. Our study introduces a novel molecular classification of EpS with prognostic value, associated diagnostic biomarkers and distinct tumor immune microenvironments, which extends beyond traditional clinicopathological classification. This paves the way for precision medicine-based therapeutic strategies for patients with EpS.

  Study EGAS50000000973

• Cancer-associated fibroblasts promote drug resistance in ALK-driven lung adenocarcinoma cells by upregulating lipid biosynthesis

  Targeted therapy interventions using tyrosine kinase inhibitors (TKIs) provide encouraging treatment responses in ALK-rearranged lung adenocarcinomas, yet resistances occur almost inevitably. Apart from tumor cell-intrinsic resistance mechanisms, accumulating evidence supports a role of cancer-associated fibroblasts (CAFs) in affecting the therapeutic vulnerability of lung cancer cells. Here, we aimed to investigate underlying molecular networks shaping the therapeutic susceptibility of ALK-driven lung adenocarcinoma cells via tumor microenvironmental cues using three-dimensional (3D) spheroid co-culture settings. We show that CAFs promote therapy resistance of lung tumor cells against ALK inhibition by reducing apoptotic cell death and increasing cell proliferation. Using single-cell RNA-sequencing analysis, we show that genes involved in lipogenesis constitute the major transcriptional difference between TKI-treated homo- and heterotypic lung tumor spheroids. CAF-conditioned medium and CAF-secreted factors HGF and NRG1 were both able to promote resistance of 3D-cultured ALK-rearranged lung tumor cells via AKT signaling, which was accompanied by enhanced de novo lipogenesis and supression of lipid peroxidation. Notably, simultaneous targeting of ALK and SREBP-1 was able to overcome the established CAF-driven lipid metabolic-supportive niche of TKI-resistant lung tumor spheroids. Our findings highlight a crucial role of CAFs in mediating ALK-TKI resistance via lipid metabolic reprogramming and suggest new ways to overcome resistance towards molecular directed drugs by targeting vulnerabilities downstream of oncogenic signaling.

  Study EGAS50000000135

• Relapse CHL study

  Classical Hodgkin Lymphoma (CHL) is characterized by a complex tumor microenvironment (TME) that supports disease progression. While immune cell recruitment by Hodgkin and Reed–Sternberg (HRS) cells is well-documented, the role of non-malignant B cells in relapse remains unclear. Using single-cell RNA sequencing (scRNA-seq) on paired diagnostic and relapsed CHL samples, we identified distinct shifts in B cell populations, particularly an enrichment of naïve B cells and a reduction of memory B cells in early-relapse compared to late-relapse and newly diagnosed CHL. The enrichment of naïve B cells in early relapse biopsies was confirmed in independent validation cohorts using scRNA-seq and immunohistochemistry. Notably, naïve B cells in early-relapse samples exhibited high expression of LGALS9, an immunosuppressive molecule known to interact with T cell immunoglobulin and mucin domain 3 (TIM-3) on regulatory T cells (Tregs). Cell-cell interaction analysis revealed the importance of interactions between LGALS9+ naïve B cells and TIM-3+ Tregs in the early-relapse setting. Spatial analysis by imaging mass cytometry confirmed close proximity of LGALS9+ naïve B cells with TIM-3⁺ CD4⁺ T cells and HRS cells, pointing to their role in shaping an immunosuppressive niche. Our findings highlight a previously unrecognized population of LGALS9⁺ naïve B cells with immunoregulatory potential in early-relapse CHL and provide new insights into the spatial and transcriptional architecture of the relapsed TME in CHL.

  Study EGAS00001008222

• Characterization of HCV-specific CD4 T cells during DAA-Therapy

  Background: Chronic HCV-infection is characterized by a severe impairment of HCV-specific CD4 T cell help that is driven by chronic antigen stimulation. We aimed to study the fate of HCV-specific CD4 T cells after viral elimination. Methods: HCV-specific CD4 T cell responses were longitudinally analyzed using MHC class II tetramer-technology, multicolor flow cytometry and RNA sequencing in a cohort of chronically HCV-infected patients undergoing therapy with direct-acting antivirals. In addition, HCV-specific neutralizing antibodies and CXCL13 levels were analyzed. Results: We observed that the frequency of HCV-specific CD4 T cells increased within two weeks after initiation of DAA therapy. Multicolor flow cytometry revealed a downregulation of exhaustion and activation markers and an upregulation of memory-associated markers. While cells with a Th1 phenotype were the predominant subset at baseline, cells with phenotypic and transcriptional characteristics of follicular T helper cells increasingly shaped the circulating HCV-specific CD4 T cell repertoire, suggesting antigen-independent survival of this subset. These changes were accompanied by a decline of HCV-specific neutralizing antibodies and the germinal center activity. Conclusion: We identified a population of HCV-specific CD4 T cells with a follicular T helper cell signature that is maintained after therapy-induced elimination of persistent infection and may constitute an important target population for vaccination efforts to prevent re-infection and immunotherapeutic approaches for persistent viral infections.

  Study EGAS00001003950

• Probabilistic cell type assignment of single-cell transcriptomic data reveals spatiotemporal microenvironment dynamics in human cancers

  Single-cell RNA sequencing (scRNA-seq) has transformed biomedical research, enabling decomposition of complex tissues into disaggregated, functionally distinct cell types. For many applications, investigators wish to identify cell types with known marker genes. Typically, such cell type assignments are performed through unsupervised clustering followed by manual annotation based on these marker genes, or via "mapping" procedures to existing data. However, the manual interpretation required in the former case scales poorly to large datasets, which are also often prone to batch effects, while existing data for purified cell types must be available for the latter. Furthermore, unsupervised clustering can be error-prone, leading to under- and over- clustering of the cell types of interest. To overcome these issues we present CellAssign, a probabilistic model that leverages prior knowledge of cell type marker genes to annotate scRNA-seq data into pre-defined and de novo cell types. CellAssign automates the process of assigning cells in a highly scalable manner across large datasets while simultaneously controlling for batch and patient effects. We demonstrate the analytical advantages of CellAssign through extensive simulations and exemplify real-world utility to profile the spatial dynamics of high-grade serous ovarian cancer and the temporal dynamics of follicular lymphoma. Our analysis reveals subclonal malignant phenotypes and points towards an evolutionary interplay between immune and cancer cell populations with cancer cells escaping immune recognition.

  Study EGAS00001003452

• Mutational landscape and patterns of clonal evolution in relapsed pediatric acute lymphoblastic leukemia

  Relapse of acute lymphoblastic leukemia (ALL) remains a leading cause of childhood death. Prior studies have shown clonal mutations at relapse often arise from relapse-fated subclones that exist at diagnosis. However, the genomic landscape, evolutionary trajectories and mutational mechanisms driving relapse are incompletely understood. In an analysis of 92 cases of relapsed childhood ALL, incorporating multimodal DNA and RNA sequencing, deep digital mutational tracking and xenografting to formally define clonal structure, we identify 50 significant targets of mutation with distinct patterns of mutational acquisition or enrichment. CREBBP, NOTCH1, and Ras signaling mutations rose from diagnosis subclones, whereas variants in NCOR2, USH2A and NT5C2 were exclusively observed at relapse. Evolutionary modeling and xenografting demonstrated that relapse-fated clones were minor (50%), major (27%) or multiclonal (18%) at diagnosis. Putative second leukemias, including those with lineage shift, were shown to most commonly represent relapse from an ancestral clone rather than a truly independent second primary leukemia. A subset of leukemias prone to repeated relapse exhibited hypermutation driven by at least three distinct mutational processes, resulting in heightened neoepitope burden and potential vulnerability to immunotherapy. Finally, relapse-driving sequence mutations were detected prior to relapse using deep digital PCR at levels comparable to orthogonal approaches to monitor levels of measurable residual disease. These results provide a genomic framework to anticipate and circumvent relapse by earlier detection and targeting of relapse-fated clones.

  Study EGAS00001003975

• Gene Expression in an African American Schizophrenia Dataset

  Schizophrenia is a common and severe psychotic disorder. While some common SNPs and rare copy number variants have been identified as being significantly associated with disease risk, the biological mechanisms remain largely undefined. To identify gene expression abnormalities in schizophrenia, we generated whole-genome gene expression profiles using RNAseq on lymphoblastoid cell lines from schizophrenia cases and controls of African American ancestry. Genes related to immune processes were enriched in the genes differentially expressed in the schizophrenia cases. Approximately half of the differentially expressed genes in this African American dataset were also differentially expressed by affection status in our previously studied European ancestry dataset (PMID: 28418402).

  Study phs002842

• Mapping Cells and Interactions in the Thymus Across Development and Aging (2025-07-28)

  Single cell RNA sequencing on thymic cells from children with syndromic diseases such as Trisomy 21 and inflammatory/autoimmune diseases such as myasthenia gravis. Our group has previously performed atlasing of the healthy human thymus using single cell transcriptomics and obtained an unprecedented understanding of human thymus throughout development. In this study, we intend to survey diseased human thymus and compare this to healthy to obtain a better understanding of thymic dysfunction in human diseases. The diseases we will survey include congenital conditions such as Down's syndrome, CHARGE syndrome, DiGeorge syndrome as well as acquired diseases such as myasthenia gravis and thymic hyperplasia. We intend to perform genetic and genomic study on the samples, including expression analysis. This will not only help us to understand the cells and molecular pathways affected in dysfunctional thymus; but may also lead to generalisable understandings in human genetic diseases and inflammatory/autoimmune diseases. This data is part of a pre-publication release. For information on the proper use of pre-publication data shared by the Wellcome Trust Sanger Institute (including details of any publication moratoria), please see http://www.sanger.ac.uk/datasharing/ . This dataset contains all the data available for this study on 2025-07-28.

  Dataset EGAD00001015660

• Epigenome maps of time-resolved monocyte to macrophage differentiation and innate immune memory

  Innate immune memory is the phenomenon whereby innate immune cells such as monocytes or macrophages undergo functional reprogramming after exposure to microbial components such as LPS. We apply an integrated epigenomic approach to characterize the molecular events involved in LPS-induced tolerance in a time dependent manner. ChIP-seq, RNA-seq, WGBS and ATAC-seq data were generated. This analysis identified epigenetic programs in tolerance and trained macrophages, and the potential transcription factors involved. Experimental set-up Time-course in vitro culture of human monocytes. Two innate immune memory states can be induced in culture through an initial exposure of primary human monocytes to either LPS or BG for 24 hours, followed by removal of stimulus and differentiation to macrophages for an additional 5 days. Cells were collected at baseline (day 0), 1 hour, 4 hour, 24 hour and 6 days.

  Dataset EGAD00001002693

• CEBP-Beta/IL-1-Beta/ TNF-alpha Feedback Loop Drives Drug Resistance to BCL2 and MDM2 Inhibitors in Monocytic Leukemia Cells

  To uncover the mechanisms of MDM2 inhibitor resistance that are mediated by IL-1, we performed RNAseq on three primary AML patient cells treated with idasanutlin in the presence or absence of IL-1β. These three leukemia samples were chosen because they exhibit high blast counts and are protected by IL-1β when treated with idasanutlin, as determined by apoptosis assay. We observed that a majority of the major p53 pathway targets remained upregulated in the presence of IL-1β. We next performed differential expressed gene (DEG) analysis following IL-1β treatment with and without idasanutlin. IL-1β treatment upregulated 3295, 1353, and 680 genes with ≥1 Log2FC, in each of the three patient samples, and a significant majority of these genes (82%, 62.4%, and 47.6%) remained upregulated in the presence of idasanutlin. Pathway analysis showed that IL-1β upregulated genes were associated with NF-κB, Interleukin-1 family signaling pathways, immune and inflammatory associated genes and myeloid maturation markers. To identify downstream targets that contribute to the protective effect of IL-1β during idasanutlin treatment, we performed DEG analysis on IL-1β + idasanutlin treatment compared to idasanutlin treatment alone. Interestingly, IL-1β induces its own expression and other pro-survival cytokines including TNF, CSF2 (GMCSF), IL-6, CSF3 (GCSF), as well as their respective receptors and/or coactivators, including IL1RL1 and multiple TNF receptors. We verified that IL-1β and TNF-α induce the upregulation of their own membrane-bound forms, their corresponding receptors, CSF2RB, and monocytic markers, as determined by flow cytometry immunophenotyping. Furthermore, IL-1β upregulated a number of pro-inflammatory chemokines, including CXCL13, CXCL2, etc. These findings highlight a pleiotropic role of IL-1β. The RNA expression of these three samples in the presence and absence of IL-1β and idasanutlin treatment will be deposited in dbGaP.

  Study phs003479

• Deep single-cell RNA sequencing data for 11138 T cells from tumour, adjacent normal tissue and peripheral blood of treatment-naive CRC patients. The DATA ACCESS AGREEMENT is provided at https://github.com/zhangyybio/single-T-cell-data-access. Applicants can request access to the data by directly downloading it or by sending an email to cancerpku@pku.edu.cn. The process that is used to approve an application includes verifying the institution, participants and research purposes of the application. In general this process will take about two weeks. In principal, any scientific research program complying with the laws and bioethic regulation policies of China will be approved.

  T cells are central players in cancer immunotherapy1, yet some of their fundamental properties such as development and migration within tumours remain elusive. The enormous T cell receptor (TCR) repertoire, required for recognising foreign and self-antigens2,3, could serve as lineage tags to track these T cells in tumours4. Here, we obtained transcriptomes of 11,138 single T cells from 12 colorectal cancer (CRC) patients and developed STARTRAC (Single T-cell Analysis by Rna-seq and Tcr TRACking) indices to quantitatively analyse dynamic relationships among 20 identified T cell subsets with distinct functions and clonalities. While both CD8+ effector and “exhausted” T cells exhibited high clonal expansion, they were independently connected with tumour-resident CD8+ effector memory cells, implicating a TCR-based fate decision. Of the CD4+ T cells, the majority of tumour-infiltrating Tregs showed clonal exclusivity, whereas certain Treg clones were developmentally linked to multiple TH clones. Notably, we identified two IFNG+ TH1-like clusters in tumours, the GZMK+ TEM and CXCL13+ TH1-like clusters, which were associated with distinct IFN-γ-regulating transcription factors, EOMES/RUNX3 and BHLHE40, respectively. Only CXCL13+BHLHE40+ TH1-like cells were preferentially enriched in tumours of microsatellite-instable (MSI) patients, which might explain their favourable response rates to immune-checkpoint blockade. Furthermore, we found IGFLR1 to be highly expressed in both CXCL13+BHLHE40+ TH1-like and CD8+ exhausted T cells and possessed co-stimulatory functions. Our integrated STARTRAC analyses provides a powerful avenue to comprehensively dissect the T cell properties in CRC, which could shed new insights into the dynamic relationships of T cells in other cancers

  Study EGAS00001002791

• Exome Sequencing to Define the Landscape of Plasma Cells in Systemic Light chain Amyloidosis

  Systemic light chain amyloidosis (AL) is characterized by the deposition of immunoglobulin light chains as amyloid fibrils in different organs, where they form toxic protein aggregates. The underlying disease is a plasma cell disorder, but limited whole exome data are available. We report the findings of an exome sequencing study in AL to define a plasma cell signature and compare this to monoclonal gammopathy of undefined significance (MGUS) and myeloma (MM). Twenty-four samples from unselected newly diagnosed untreated AL patients were analysed. CD138+ cells were isolated from bone marrow cells using MACSorting (Miltenyi Biotech, Bisley, UK). Cells were lysed in RLT+ buffer and DNA/RNA extracted using the AllPrep kit (Qiagen, Manchester, UK). Peripheral blood was isolated from patients, white blood cells purified by Ficoll-Pacque and DNA extracted using the QIAamp DNA mini kit (Qiagen). There were a median of 47 acquired nonsynonymous variants (range 4-253) per patient. Thirty percent of patients yielded a mutation in a myeloma driver gene. These mutations were either clonal or sub-clonal. There was evidence of MAPK activation with NRAS mutations, NF-κB activation and DNA repair pathway alterations. Exome data were used to determine the cytogenetic groups of AL samples and identified hyperdiploidy in 33% of cases and t(11;14) in 33% of cases. All translocations were generated via class switch recombination and breakpoints on chromosome 11 fell into the expected regions. A novel translocation, unlikely to be driving disease progression was identified involving RCC1. Finally, further analysis of the APCS gene encoding for the SAP protein did not show any biased genotype in the AL population in comparison to the general population. We conclude that exome sequencing identifies a genetic make-up of AL that is similar to other plasma cell disorders such as MM and MGUS. This not only includes copy number abnormalities and translocations but also nonsynonymous mutations.

  Study EGAS00001001418

• Influence of the Microbiome on Epigenetic Mechanisms in IBD

  Live CD4 T cells were sorted from inflamed and non-inflamed tissue samples of IBD patients or from healthy and IBD blood samples. ATAC-Seq libraries were generated from live CD4 T cells sorted from i) inflamed and non-inflamed tissue samples, ii) healthy and IBD blood samples, or from iii) CD4 T cell subsets polarised from healthy blood samples. After isolating crude nuclei, live CD4+ T cells were treated with Tagment DNA buffer and Tagment DNA Enzyme (Nextera DNA Library Prep Kit, Illumina), and then the DNA was purified by MinElute PCR Purification Kit (Qiagen). Transposed DNA fragments were amplified using specific adapters followed by purification with MinElute PCR Purification Kit (Qiagen). Fragments from 240-360pb were selected in the PippinHT system (Sage Science). The quality of the library and its DNA concentration were assessed by Bioanalyzer instruments (Agilent Technologies) and ultimately submitted for sequencing using Illumina HiSeq 2500 sequencer, V4 chemistry. On the other hand, single cell RNA-Seq libraries were generated exclusively from inflamed and non-inflamed tissue samples of Crohn’s disease patients. Briefly, live CD4 T cells were captured and encapsulated before cDNA amplification using the 10X Genomics Chromium Platform. Samples were prepared as outlined by 10x genomics Single Cell 3’ Reagent Kits v2 user guide. Samples were sequenced on a HiSeq 2500 with the following run parameters: Read 1 – 26 cycles, read 2 – 98 cycles, index 1 – 8 cycles.

  Dataset EGAD00001011066

• Single-Cell Mitochondrial Mutation Lineage Tracing of Non-Dysplastic and Dysplastic Barrett's Esophagus

  This study addresses a critical knowledge gap in understanding Barrett's esophagus (BE), a common metaplastic condition and precursor to esophageal adenocarcinoma, by identifying lineage connections and cell states underlying its development and malignant transformation. We collected biopsies from 13 Barrett's esophagus patients undergoing routine surveillance endoscopy, sampling metaplastic tissue alongside normal squamous esophagus and gastric cardia tissues when clinically appropriate. Our novel approach integrated single-cell RNA sequencing of nearly 180,000 cells with mitochondrial DNA mutation analysis for lineage tracing, coupled with spatial transcriptomics and comprehensive whole-exome sequencing. For one patient with high-grade dysplasia, we performed detailed whole-exome sequencing analysis at 100x coverage on dissociated Barrett's esophagus cells alongside matched normal esophagus and gastric cardia tissues. This genomic profiling revealed critical genetic drivers of malignant transformation, including a truncating mutation in the mutation cluster region of APC (a negative WNT regulator), mutations in the tumor suppressor genes CDKN2A and TP53, and significant copy number alterations including loss of CDKN2A on chromosome 9. These genetic alterations corresponded with the transcriptional changes observed in our single-cell data, particularly WNT pathway dysregulation and the expansion of LGR5-expressing stem cells, providing molecular evidence for how single clones can initiate the dysplastic transformation. Data deposited in dbGaP consist specifically of the whole-exome sequencing raw files, processed variant calls, and copy number analyses from the single patient with high-grade dysplasia and matched normal esophagus and stomach control tissue. Whole exome sequencing was performed on an Illumina NovaSeq 6000. Sequencing data was preprocessed following the Genome Analysis Toolkit's (GATK) Best Practices. Sequencing data was analyzed with Mutect2 (GATK, v4.2.5.0) to identify short somatic mutations including single-nucleotide polymorphisms and insertions and deletions. The Barrett's esophagus sample was run with the two matched normals, as well as the publicly available germline resource somatic-hg38_af-only-gnomad.hg38.vcf.gz (GATK). Somatic variant calls were filtered using FilterMutectCalls (GATK) and loaded into the Integrative Genomics Viewer (v2.12.3) for identification of high-quality variants in known esophageal adenocarcinoma regulators. Copy number alterations in the Barrett's esophagus sample were estimated using CNVkit (v0.9.9). CNV, Mutect2, and SVC analyses will be made available through dbGAP. Copy number alterations in the Barrett's esophagus sample were estimated using CNVkit (v0.9.9) with the default run settings; a pooled reference was generated from the normal esophagus and stomach samples. Copy number results were plotted directly within CNVkit using the function scatter.

  Study phs003949

• Decoding Human Heart Morphogenesis through Single-cell Multi-modal Analyses

  The heart, which is the first organ to develop, is highly dependent on its form to function. However, how diverse cardiac cell types spatially coordinate to create complex morphological structures critical for heart function remains to be elucidated. Here, we show that integration of single cell RNA-sequencing with high-resolution multiplexed error-robust fluorescent in situ hybridization (MERFISH) not only resolves the identity of cardiac cell types developing the human heart but also provides a spatial mapping of individual cells that enables illumination of their organization into cellular communities forming distinct cardiac structures. We discovered that many of these cardiac cell types further specified into subpopulations exclusive to specific communities, supporting their specialization according to cellular ecosystem and anatomic region. In particular, ventricular cardiomyocyte subpopulations displayed an unexpected complex laminar organization across the ventricular wall and formed, with other cell subpopulations, several cellular communities. Interrogating cell-cell interactions within these communities revealed signaling pathways orchestrating the spatial organization of cardiac cell subpopulations during ventricular wall morphogenesis. In vivo conditional genetic mouse models and in vitro human pluripotent stem cell studies confirmed an intricate multicellular PLXN-SEMA crosstalk among specific ventricular cardiomyocyte, fibroblast and endothelial cell subpopulations that directs the compaction of the ventricular wall layers. Thus, these detailed findings into the cellular social interactions and specialization of cardiac cell types constructing and remodeling the human heart offer new insights into structural heart diseases as well as engineering complex multi-cellular tissues for human heart repair. Reprinted from Farah et al. (In press), with permission from Nature.

  Study phs002031

• Non-Coding Autoimmune Risk Variant Defines Role for ICOS in T Peripheral Helper Cell Development

  Our fine-mapping of 76 non-MHC loci associated with risk for rheumatoid arthritis (RA, 11,475 cases, 15,870 controls) has recently identified rs117701653, a non-coding single nucleotide polymorphism (SNP) in the CTLA4/CD28/ICOS locus, as the variant most likely to modulate RA risk within this region, with the minor allele C showing disease protection compared to the major allele A. We found that this SNP exhibits allele-specific protein binding and transcriptional modulation, further supporting its regulatory nature. Here, we aimed to determine whether the functional non-coding variant rs117701653 modulates the expression of genes within the CD28/CTLA4/ICOS locus. To address this question, we employed low-input RNA sequencing using the Illumina SmartSeq2 platform on resting total CD4+ T cells from the blood of 24 healthy subjects (8 AA, 8 AC, 8 CC) recruited based on genotype from the Mass General Brigham biobank. We performed a targeted cis-eQTL mapping analysis for protein coding genes with transcription start sites (TSS) within a 1MB window of SNP rs117701653. These genes included WDR12, NBEAL1, CYP20A1, ABI2, RAPH1, CD28, CTLA4, ICOS, and PARD3B. Surprisingly, allelic variation at rs117701653 correlated strongly with expression of ICOS. This suggests that individuals carrying the protective allele C at rs117701653 express lower levels of ICOS, implicating ICOS as a plausible target gene influenced by this protective allele. The raw sequence data from the 24 healthy individuals are available through dbGAP. Gene expression matrices have been deposited to GEO accession GSE235868.

  Study phs003448

• Case Report: Pre-Clinical Combination Targeting VEGF and PI3K in a Rare, Aggressive Mixed Endometrial Carcinoma

  We report a rare case of a young patient (VENUS 167) initially diagnosed with grade 1 endometrioid endometrial cancer, which, following endocrine treatment, presented with mixed aggressive carcinoma with three distinct histologic patterns: grade 1 endometrioid, large cell neuroendocrine, and undifferentiated carcinoma. The surgical specimen at the time of disease progression was used to establish OPTO.85, a patient-derived organoid (PDO), followed by a corresponding organoid-derived xenograft (ODX). Multi-omic analyses confirmed that OPTO.85 accurately reflected the patient’s tumor characteristics. WES analysis identified oncogenic alterations in PIK3CA, ARID1A, and CTNNB1. Further RNA-seq and ATAC-seq analyses revealed enrichment in VEGF and Wnt signaling pathways, suggesting therapeutic vulnerabilities. A high-throughput drug screen was conducted using ApexBio-Approved and Epigenetic Drug Libraries, along with Kinase Inhibitor and Tool Compound Libraries developed at the Ontario Institute of Cancer Research. The OPTO.85 PDO exhibited sensitivity to PI3K inhibitors and responsiveness to VEGF inhibition. Cediranib demonstrated synergy with BKM120, significantly reducing organoid growth. This combination also showed in vivo efficacy in the ODX model, where dual inhibitors significantly suppressed tumor growth compared to single compounds. This case exemplifies the impact of genomic profiling and patient-derived models in identifying actionable molecular changes in rare cancers with limited therapeutic options and poor prognosis. It highlights that high-throughput sequencing for individual patient tumors and generation of patient-derived models are feasible in endometrial cancer. This preclinical model may assist clinical decision and personalized therapy requiring validation in prospective studies.

  Study EGAS50000001665

• Colorectal cancer organoids expressing BRAF (fusion) genes

  Fusion genes may be oncogenic drivers and potential targets for personalized therapies in a variety of cancer types. The BRAF gene is frequently involved in oncogenic fusions, with fusion frequencies of 0.2-3% throughout different cancers. BRAF can be fused to a wide variety of genes, however, BRAF fusions rarely occur in the exact same gene configuration, making assessment of the fusion relevance and decision for drug treatment challenging. Here, we devised a colorectal cancer organoid-based platform to functionally characterize diverse BRAF fusions, containing various partner genes (AGAP3, DLG1 and TRIM24). We compared these BRAF fusions with respect to cellular behaviour, downstream signaling activation and response to targeted therapies. We found that 5’ partner choice of BRAF affects cellular localization and intracellular signaling capacities of the fusion genes. The DLG1-BRAF fusion gene showed distinct localization to the plasma membrane and exhibited increased levels of MAPK pathway activation under unperturbed conditions. Furthermore, the different BRAF fusions showed varying sensitivities to the targeted inhibition of ERK and BRAF, with the DLG1-BRAF fusion being the most sensitive. Additionally, RNA-sequencing identified distinct subsets of genes affected by the DLG1-BRAF fusion gene. Importantly, all fusion genes conveyed resistance to the clinically relevant EGFR/HER2/HER4-inhibitor afatinib, suggesting that BRAF fusions should be screened alongside other MAPK pathway alterations to identify mCRC patients amenable to cetuximab treatment. In summary, we developed a platform to efficiently assess molecular and cellular effects of fusion genes, and revealed that differential drug responses and distinct gene expression profiles can be induced by different BRAF fusions .

  Study EGAS00001003558

• Healthy_ageing_thymus

  Single cell RNA sequencing of immune and non-immune cells from healthy ageing thymus Thymus is a primary lymphoid organ that creates an environment for the T lymphocyte precursors differentiation into the naive T cells, which involves TCR gene rearrangement, negative and positive selection. Despite thymus critical function for the recognition of the pathogens it starts to atrophy very early in life. Thymic Epithelial Space (TES), where thymocyte education occurs, starts to decline rapidly just after birth, while with the onset of puberty thymus is also progressively replaced with adipose tissue. This leads to a reduced naive T cell output,which decreases an organism's ability to recognize pathogens. In our lab, we have already collected and characterized cell type composition and gene regulatory networks of the fetal tand adult human thymus (Park et al, Science 2020). In the current project we plan to expand the number of paediatric and adult samples in the study with the aim to understand how the process of thymic involution happens. We will use single-cell transcriptomics to zoom on changes that occur in immune and non-immune cell types in different phases of thymic involution with the hope to understand which gene expression changes are causal for the thymic involution process and how naïve T cell production can be increased in elderly. This data is part of a pre-publication release. For information on the proper use of pre-publication data shared by the Wellcome Trust Sanger Institute (including details of any publication moratoria), please see http://www.sanger.ac.uk/datasharing/

  Study EGAS00001004311

• Transcriptomic intra-tumor heterogeneity of colorectal cancer varies depending on tumor location within the colorectum

  Background; Intra-tumor heterogeneity (ITH) of colorectal cancer (CRC) complicates molecular tumor classification, such as transcriptional subtyping. Differences in cellular states, biopsy cell composition, and tumor microenvironment may all lead to ITH. Here we analyze ITH at the transcriptomic and proteomic levels to ascertain whether subtype discordance between multiregional biopsies reflects relevant biological ITH or lack of classifier robustness. Further, we study the impact of tumor location on ITH.MethodsMultiregional biopsies from stage II and III CRC tumors were analyzed by RNA sequencing (41 biopsies, 14 tumors) and multiplex immune protein analysis (89 biopsies, 29 tumors). CRC subtyping were performed using consensus molecular subtypes (CMS), CRC intrinsic subtypes (CRIS), and TUMOR types. ITH-scores and network maps were defined to determine the origin of heterogeneity. A validation cohort was used with one biopsy per tumor (162 tumors).ResultsOverall, inter-tumor transcriptional variation exceeded ITH, and subtyping calls were frequently concordant between multiregional biopsies. Still, some tumors had high ITH and were classified discordantly. Proximal tumors showed ITH related to differences in the microenvironment, while ITH within distal tumors were cancer-cell related. This subtype discordancy reflected actual molecular ITH within the tumors. The relevance of the subtypes was reflected at protein level where several inflammation markers were significantly increased in immune related transcriptional subtypes, which was verified in an independent cohort (Wilcoxon rank sum test; p<0.05).ConclusionOur transcriptomic and proteomic analyses show that the tumor location along the colorectum influence the ITH of CRC, which again influence the concordance of subtyping.

  Study EGAS00001004668

• Transcriptomic Analysis of HIV-Infected Cells

  Defining the attributes of rare CD4 T cells that harbor latent HIV under therapy with anti-retroviral agents has been challenging due to difficulties in identifying and isolating these cells. These cells represent an important barrier to cure, and characterizing them would aid the design and development of novel HIV cure strategies. In this study, we addressed this challenge using a custom microfluidic technology named Focused Interrogation of cells by Nucleic acid Detection and Sequencing (FIND-seq) that isolates HIV-infected cell transcriptomes, based on HIV DNA detection. We applied this method to HIV DNA+ memory CD4 T cells isolated from blood from five persons with HIV who were receiving anti-retroviral agents. After curating samples based on data quality, we proceeded with bioinformatics analysis from three persons. This analysis identified six transcriptomic pathways that were inhibited in HIV DNA+ cells compared to HIV DNA- cells. These included death receptor signaling, necroptosis signaling, and anti-proliferative G-alpha12/13 signaling. Furthermore, we identified two gene clusters, including 60 and 85 genes, that were significantly associated with HIV DNA+ cells identified by co-expression network analysis. These genes included negative regulators of HIV transcription that were higher in HIV DNA+ cells, positive regulators of HIV transcription that were lower in HIV DNA+ cells, and other genes involved in the negative regulation of mRNA translation, additional RNA processing functions, and the regulation of cellular state and fate. We further examined these signatures in transcriptomic data from CD4 T cell subsets sorted by flow cytometry from nine persons with HIV. We found a partial similarity between the smaller cluster of genes and the signature of CCR6- peripheral TFH cells. Raw sequencing data from all eleven human participant samples included in the study are available. These include samples from two participants who were studied only by FIND-seq, three participants who were studied by both FIND-seq and flow cytometry and six participants who were studied only by flow cytometry.

  Study phs003095

• Adult Eosinophilic Esophagitis Registry Atlas

  Eosinophilic esophagitis (EoE) arises in the setting of type 2 cytokine-driven inflammation, is characterized by the infiltration of eosinophils in the esophageal mucosa, and is associated with dysphagia and fibrostenotic complications. Disease development and progression involve complex interactions between epithelial cells, stromal cells and immune cells. Indeed, increased epithelial permeability is consistently found in EoE, potentially leading to increased antigen exposure and inflammation. Persistent inflammation can in turn lead to fibroblast activation and fibrosis, which, together with the mediators produced by eosinophils and mast cells, leads to dysmotility. To decipher the interactions between these cell types across distinct disease states, we obtained biopsies from 7 healthy controls and 15 EoE subjects, including both active disease and clinical remission, to perform single cell RNA-sequencing (scRNA-seq).Each biopsy sample was enzymatically digested within 3 hours of collection with collagenase/dispase to generate a single cell suspension. Suspensions were then loaded on a 10x Chromium controller, and libraries were generated with the 3' Chemistry pipeline (v2 or v3) before sequencing on an Illumina sequencer. The resulting sequence data was demultiplexed and analyzed with the CellRanger pipeline to generate raw fastq files and processed gene-count matrices for each sample. Altogether, our study generated transcriptomic profiles from 421,312 individual cells after quality control and filtering.Our analyses reveal 60 cell subsets and functional alterations in cell states, compositions, and interactions that highlight previously unknown features of EoE. Active disease displayed enrichment of ALOX15+ macrophages, novel PRDM16+ dendritic cells expressing the EoE risk gene ATP10A, and cycling mast cells, with concomitant reduction of TH17 cells. Receptor–ligand expression uncovered eosinophil recruitment programs, increased fibroblast interactions in disease, and IL-9+IL-4+IL-13+TH2 and endothelial cells as potential mast cell interactors. Resolution of inflammation-associated signatures included mast and CD4+ Tissue-resident memory T (TRM) cell contraction and cell type-specific downregulation of eosinophil chemoattractant, growth, and survival factors.

  Study phs003574

• Genomic Basis of Phenotypic Variability of Complex Disorders

  The 520-kbp deletion on chromosome 16p12.1 is associated with multiple neurodevelopmental outcomes, including intellectual disability/developmental delay (ID/DD), schizophrenia, autism, and epilepsy. This variant is inherited in >95% of cases from carrier parents manifesting neuropsychiatric features, while affected children with the deletion were more likely to carry another large Copy Number Variation (CNV) or rare deleterious mutation elsewhere in the genome, or "second-hits", compared to carrier parents. Thus, while the 16p12.1 deletion confers significant risk for a range of disorders, the ultimate phenotypic trajectory is determined by "second-hit" variants in the genetic background. However, outstanding questions remain about the genetic background in these disorders, including the functional effects of "second-hit" variants towards gene expression or developmental phenotypes, the full spectrum of "second-hit" variants, and the mechanisms driving accumulation of these variants.In version 1, we performed deep clinical phenotyping, whole-genome sequencing (WGS), SNP microarray analysis, and RNA-sequencing of lymphoblastoid cell lines of 32 individuals in five large families with the 16p12.1 deletion. We found that specific classes of rare single-nucleotide and structural variants in coding and non-coding regions contributed towards expression changes in affected children based on their inheritance patterns, including a subset of variants that acted synergistically with the deletion to alter expression. We further found that rare variant-mediated expression changes contribute to specific developmental phenotypes in the carrier children, suggesting a potential genetic mechanism for variable expressivity of the deletion.In version 2, we additionally analyzed genetic and phenotypic data from individuals with 16p12.1 deletion and their spouses to assess the role of assortative mating modulating phenotypic outcomes of the 16p12.1 deletion. We identified significant within- and cross- disorder phenotypic correlations between spouses, but observed no significant genetic correlations between spouses, in line with expectations.

  Study phs002450

• RNA-Sequencing of cervical cancers

  Question: Can we identify cervical cancer patients who are at risk for distant metastatic (DM) recurrence following treatment with radiotherapy, concurrent weekly cisplatin and brachytherapy (RTCT). Findings: An immune-based 55 gene risk score was developed using a cohort of 81 patients treated with RTCT that was strongly predictive of DM and cause-specific survival (CSS). The risk score was validated in two independent patient cohorts. A high immune metastatic risk score was associated with a high tumor mutational burden and a ‘cold’, immune-excluded tumor microenvironment at diagnosis. Meaning: The immune gene expression risk score may help to identify patients at risk of DM and potential targets for mitigating this risk.

  Study EGAS50000000087