3781 results for "RNA AND sequencing OR transcriptome"

in 44.69 milliseconds.

• Ongoing Replication Stress Tolerance and Clonal T Cell Responses Distinguish Liver and Lung Recurrence and Patient Outcomes in Pancreatic Ductal Adenocarcinoma

  Patients with metastatic pancreatic ductal adenocarcinoma (PDAC) survive longer when disease spreads to the lung but not to the liver. We generated overlapping, multi-omic datasets to identify molecular and cellular features that distinguish patients whose disease develops liver metastasis (liver cohort) from those whose disease progression results in lung metastasis without liver metastases (lung cohort). Lung cohort patients generally survived longer than liver cohort patients, independent of tumor subtype. We developed a pORG gene signature that distinguishes primary tumors in the liver and lung cohorts. We identified ongoing replication stress (RS) response pathways in high pORG/liver cohort tumors, while low pORG/lung cohort tumors had greater densities of lymphocytes and shared T cell clonal responses. Our study demonstrates that liver-avid PDAC is associated with tolerance to ongoing RS, limited tumor immunity, and less favorable outcomes; whereas low RS, lung-avid/liver-averse tumors are associated with active tumor immunity that may account for favorable outcomes. As expected, we found high frequencies of KRAS, TP53, CDKN2A, and SMAD4 gene alterations in our tumor samples. High pORG activity (GSVA scores) in our primary tumors appear to be positively correlated with alteration frequencies in both TP53 and CDKN2A. We did not see a correlation with alterations in KRAS as almost all the tumor samples have KRAS alterations. From a de-identified dataset of 1,873 patients diagnosed with and/or treated for PDAC at our institution between 2004 and 2020, we identified 422 patients for which we had specimens with sequencing data (N=374) and/or specific evidence of disease metastasis site(s) from the OHSU cancer registry or disease-relevant computed tomography (CT) scans to allow cohort classification. Note that our study includes RNA-Seq, DNA-Seq panel, and TCR-Seq (T-cell Receptor Sequencing) data, but that this dbGaP submission just includes the patients/samples with RNA-Seq and/or DNA-Seq panel data. Therefore, this submission includes 290 samples from 278 patients. TCR sequence data is available on the Adaptive Biotechnologies platform. Clinical course timepoints, patient demographics, stage, grade, nodal involvement, resection margins, and angiolymphatic invasion were provided as deidentified data by the OHSU cancer registry with quality control data verification in a subset by pathologists (BB and TM). We reviewed all available computed tomography (CT) scans for all patients with primary tumor resection dates recorded by the cancer registrar, with tumor samples analyzed by RNA-seq, DNA-seq, or TCR-seq, and/or with additional information indicating metastatic spread (e.g., metastatic samples received for related studies). We abstracted the site of all lesions proven to be metastatic by biopsy and/or that clearly increased in size during progression or decreased in size during treatment as long as a radiologist described the lesion as “likely”, “suspicious for”, “concerning for”, or “favor” metastasis. Clinical imaging was reviewed by a radiologist (AG) to validate patient assignments to the liver, lung, and neither liver nor lung (other recurrence site) cohorts. Time to recurrence was calculated from the earliest of either the recurrence date provided by the OHSU cancer registry, or the date of earliest lesion abstracted from CT reports. The subject attributes reported in this submission were up to date at the time of our data freeze, which occurred July 2021.

  Study phs003597

• B cell receptor silencing reveals origin and dependencies of high-grade B cell lymphomas with MYC and BCL2 rearrangements

  High-grade B cell lymphomas with MYC and BCL2 rearrangements (HGBCL-DH-BCL2) represent highly aggressive Germinal center (GC) derived B cell malignancies with dismal prognosis and limited therapeutic options. In the present study, we investigated the expression of BCR complex components, including immunoglobulin (IG) heavy (H) and light (L) chains, and CD79A and CD79B signaling subunits in HGBCL-DH-BCL2, scoring 65% of cases with undetectable IGH immunoreactivity (IGHUND). Transcriptomics (bulk RNA-seq) and genetic (whole exome/genome sequencing) analyses indicated that IGHUND and IGH+ HGBCL-DH-BCL2 lymphomas represent distinct biological entities. Specifically, IGHUND HGBCL-DH-BCL2 exhibit a GC DZ-like molecular program, and select subset-specific gene mutations pointing to distinct evolutionary trajectories respect to the IGH+ counterparts. IG profiling revealed that IGHUND HGBCL-DH-BCL2 preserved potentially productive, hypermutated IGH variable (V) gene rearrangements, indicating a selective advantage conferred to the malignant B cells by residual IGH molecules. Accordingly, disruption of the productively rearranged IGHG1 gene in BCRnull HGBCL-DH-BCL2 cell line model confirmed the dependency on IGH expression for optimal fitness. Using transcriptomics and in situ RNA-scope analysis, we reported that IGHUND HGBCL-DH-BCL2 originated from IG-class switched GC B cells, while IGH-expressing tumors were selected for continuous IGM expression. Moreover, while IGH+ HGBCL-DH-BCL2 favor IGM/IG-Kappa expression, IGHUND counterparts complete IG-Lambda light chain rearrangements. Notably, BCR silencing, caused by accelerated IGH turnover and reduced IGH expression, precedes HGBCL-DH-BCL2 onset, inducing RAG1/2-dependent IG light chain editing and eventually facilitating t(8;22)/IGL::MYC translocations. Finally, IGH silencing observed in HGBCL-DH-BCL2 cell models, which leads to intracellular retention/silencing of CD79B, results in reduced sensitivity to the CD79B-targeting antibody-drug conjugate Polatuzumab-vedotin. Collectively, HGBCL-DH-BCL2 originates from isotype-switched t(14;18)+ GC B cells that undergo IG light chain editing fuelling intra-clonal diversification, BCR silencing, and the acquisition of t(8;22) translocation, while preserving IGH dependency.

  Study EGAS50000001047

• Identifying transcriptional programs underlying anti-EGFR small molecule response and resistance with TraCe-seq

  Genetic and non-genetic heterogeneity within cancer cell populations represents a major challenge to anti-cancer therapies. We currently lack robust methods to determine how pre-existing and adaptive features affect cellular responses to therapies. Here, by conducting clonal fitness mapping and transcriptional characterization using expressed barcodes and single-cell RNA-sequencing, we have developed TraCe-seq, a method that captures at clonal resolution the origin, fate, and differential early adaptive transcriptional programs of cells in a complex population in response to distinct treatments. We used TraCe-seq to benchmark how next-generation dual EGFR inhibitors-degraders compare to standard EGFR kinase inhibitors in EGFR-mutant lung cancer cells. To QC the TraCe-seq strategy, single-cell RNA-seq libraries were generated from a variety of human cancer cell lines transduced with the TraCe-seq library to validate the TraCe-seq strategy. Specifically, 5 different cell lines (PC9, MCF-10A, MDA-MB-231, NCI-H358, and NCI-H1373) were each transduced with a unique TraCe-seq barcode. The transduced cells were selected with puromycin only, dissociated to single cell suspensions, and then mixed together. The complex mixture of the 5 cell lines was profiled by 10X scRNA-seq. Furthermore, transduced NCI-H1373 cells were sorted by FACS to enrich for the top 50% of eGFP positive cells, and sorted cells were cultured briefly and used to construct scRNA-seq libraries and profiled by 10x scRNA-seq. To carry out the full TraCe-seq experiment, ~600 PC9 cells carrying unique TraCe-seq barcodes were expanded over 12 doublings to establish the barcoded population. A subset of the barcoded PC9 population was used to generate scRNA-seq libraries and profiled by 10x scRNA-seq prior to treatment to establish a baseline transcription profile for each barcoded clone. The rest of the cells were then treated for four days with 1 µM erlotinib, 1 µM GNE-069, or 1 µM GNE-104 respectively. scRNA-seq libraries were then generated form the treated cells and profiled by 10x scRNA-seq.

  Study EGAS00001005405

• Mapping the epigenomic landscape of human monocytes following innate immune activation reveals context-specific mechanisms driving endotoxin tolerance

  Background Monocytes are key mediators of innate immunity to infection, undergoing profound and dynamic changes in epigenetic state and immune function which are broadly protective but may be dysregulated in disease. Here, we aimed to advance understanding of epigenetic regulation following innate immune activation, acutely and in endotoxin tolerant states. Methods We exposed human primary monocytes from healthy donors (n=6) to interferon-γ or differing combinations of endotoxin (lipopolysaccharide), including acute response (2hr) and two models of endotoxin tolerance: repeated stimulations (6+6hr) and prolonged exposure to endotoxin (24hr). Another subset of monocytes was left untreated (naïve). We identified context-specific regulatory elements based on epigenetic signatures for chromatin accessibility (ATAC-seq) and regulatory non-coding RNAs from total RNA sequencing. Results We present an atlas of differential gene expression for endotoxin and interferon response, identifying widespread context specific changes. Across assayed states, only 24-29% of genes showing differential exon usage are also differential at the gene level. Overall, 19.9% (6884 of 34616) of repeatedly observed ATAC peaks were differential in at least one condition, the majority upregulated on stimulation and located in distal regions (64.1% vs 45.9% of non-differential peaks) within which sequences were less conserved than non-differential peaks. We identified enhancer-derived RNA signatures specific to different monocyte states that correlated with chromatin accessibility changes. The endotoxin tolerance models showed distinct chromatin accessibility and transcriptomic signatures, with integrated analysis identifying genes and pathways involved in the inflammatory response, detoxification, metabolism and wound healing. We leveraged eQTL mapping for the same monocyte activation states to link potential enhancers with specific genes, identifying 1,946 unique differential ATAC peaks with 1,340 expression associated genes. We further use this to inform understanding of reported GWAS, for example involving FCHO1 and coronary artery disease. Conclusion This study reports context-specific regulatory elements based on transcriptomic profiling and epigenetic signatures for enhancer-derived RNAs and chromatin accessibility in immune tolerant monocyte states, and demonstrates the informativeness of linking such elements and eQTL to inform future mechanistic studies aimed at defining therapeutic targets of immunosuppression and diseases.

  Study EGAS00001007362

• Genomic Profiling of Peripheral T-cell Lymphoma

  Peripheral T-cell lymphomas (PTCLs) are malignant and highly aggressive hematologic tumors arising from mature post-thymic T cells. The diagnosis of PTCL includes diverse lymphoma subgroups, altogether accounting for about 15% of all non-Hodgkin lymphomas. Despite recent advances in uncovering the mutational landscape of PTCLs, the prognoses of most PTCLs remain poor. To overcome the limits of exome sequencing and identify new genomic abnormalities in PTCL, we performed RNAseq and whole genome sequencing. From our analyses, we identified a recurrent in-frame deletion (VAV1 Δ778-786) generated by a focal deletion-driven alternative splicing mechanism as well as novel VAV1 gene fusions (VAV1-THAP4, VAV1-MYO1F, and VAV1-S100A7). These genetic lesions result in increased activation of VAV1 catalytic-dependent (MAPK, JNK) and non-catalytic-dependent (nuclear factor of activated T cells, NFAT) VAV1 effector pathways. These results support an oncogenic role for VAV1 signaling in the pathogenesis of PTCL (Abate et al., 2017). Recently, we have also identified a novel chimeric oncogene, FYN-TRAF3IP2, in the RNA-seq and whole genome sequencing analyses. FYN-TRAF3IP2 leads to aberrant NF-κβ activity downstream of T-cell receptor signaling and induces NF-κβ-driven T-cell transformation when expressed in hematopoietic progenitors in mice. This novel gene fusion supports the oncogenic driver role of NF-κβ signaling in the pathogenesis of PTCL (Kim et al., in review).

  Study phs001962

• Whole Tumor spatial heterogeneity in metastatic melanoma

  Heterogeneous inter-tumoral responses, sustained periods of apparent clinical benefit despite lack of objective response to various therapies, and even spontaneous remission are well known within a subpopulation of advanced melanoma patients. The molecular and cellular dynamics facilitating long-term survival remain poorly defined, particularly in the current era ofexposure to multiple potentially active therapies. We studied an exceptional case of long-term survival in a patient with non-responding metastatic melanoma in order to characterize the clonal and microenvironmental factors.We performed immunogenomic analysis of 3 metachronous tumors, and of 67 intratumor sub-regions of one metastatic lesion, usingwhole exome sequencing, RNA-sequencing, immunohistochemistry, and T cell receptor sequencing. Longitudinal analysesidentified mutations in severalknownresistance genes affecting distinctmetastases. Intra-tumoral immune signatures revealed spatially-distinct pockets of immune activation and suppression throughout the metastasis. A specific T cell VBeta CDR3 rearrangement was identified as both dominant and recurrent not only across multiple spatial points within a single tumor mass, but also across metachronous tumors spanning the patient’s disease course. Immunophenotyping of this T cell population with single-cell sequencingsuggestedrepeatedT-cell priming events leading to the persistence of both activated and exhausted T cells bearing this same TCR-Beta at any given time.Our findings highlight genomic and immune heterogeneityin metastatic disease, extending to the intra-tumoral leveland provide a rationale for integrative spatiotemporal analysis of metastatic lesions in the context of immune checkpoint blockade.

  Study EGAS00001003292

• Differential Presence of Exons in Cell-Free DNA Reveals Different Patterns in Colorectal Cancer Between Metastatic and Non-Metastatic Patients

  Tumor genotyping during disease follow-up is gaining an increasing interest for monitoring and clinical management in colorectal cancer (CRC). Growing evidences suggest that single biopsies are inefficient for molecular profiling, due to clonal evolution and intratumoral heterogeneity of primary tumors and/or metastasis. Circulating free DNA (cfDNA) has emerged as an alternative source of genetic material that seems to be representative of the continuously changing tumor molecular features. Thus, the so-called 'liquid biopsy' also offers the additional advantages of accessibility and non-invasiveness for the patient. Novel strategies based on next-generation sequencing (NGS) of cfDNA are being currently developed, with most studies focusing on targeted deep sequencing panels of potential clinically actionable genes. Here we describe an alternative approach, consisting of exome-sequencing in cfDNA at a relatively shallower depth to provide a more general overview of the circulating genome. In this study, two groups of colorectal cancer patients, with either disseminated (M) or localized disease (N), were compared using whole-exome sequencing and a RNA-seq bioinformatics analysis pipeline. Thus, a set of 379 exons present at different levels in cfDNA from both groups were identified, giving rise to a completely new concept in this field, termed as 'Differential Presence of Exons'. Differentially present exons (DPEs) were subsequently used for M and N patients clustering and classification, allowing for the design of a predictive algorithm with encouraging preliminary results in a small subset of patients that could not be initially classified attending to the selection criteria (U group, unassignable).

  Study EGAS00001002687

• Determination of cell specific regulatory enhancers in hematopoetic models

  The BLUEPRINT project is a large-scale project investigating epigenetic mechanisms involved in blood formation, in health and disease. The human variation workpackage (WP10, led by NS) of the project seeks to characterize the effect of common sequence variation on the epigenome status of a cell. To do this, the project will use highly purified blood cells to minimise "experimental noise" and therefore enhance the power to discover modest effects. Two peripheral blood cell types, the CD14+CD16- monocyte (an important central orchestrator of adaptive immunity and a bridge between innate and adaptive immunity) and the CD65+CD9- neutrophilic granulocyte (the frontline cell for innate immunity) have been selected for this purpose. The two types of cells will be obtained at high purity from adult blood (AB) of 200 healthy males and females, respectively. Cells will be purified by using already validated and fully operational protocols that are based on density gradient centrifugation of the buffy coat obtained from whole blood, followed by magnetic bead-based purification using monoclonal antibodies against Cluster of Differentiation (CD) lineage-specific cell surface markers. Units of 475 ml of AB will be obtained from consenting volunteers of the Cambridge BioResource (CBR), a panel of 10,000 healthy volunteers local to Cambridge who have already consented to participate in biomedical research and of whom biological samples (DNA, plasma, serum) and lifestyle data have been deposited in a repository and database, respectively. We are requesting funding from the Human Diversity project to sequence the genomes of the 200 CBR volunteers at low pass (6x coverage). Nuclei, DNA and RNA will be recovered from the purified cells and made available for RNA-seq, DNA-seq and ChIP-seq and genomic DNA for entire genome sequencing will be recovered from the DNA repository.

  Dataset EGAD00001002205

• RNA sequencing of a glioblastoma PDX cohort

  Patient-derived xenograft (PDX) models represent an essential tool for the preclinical evaluation of novel targeted therapies and their combinations. In this project, we molecularly characterized a comprehensive panel of glioma PDX models with well conserved molecular and pathological features over multiple passages. PDX samples were subjected to next generation sequencing (NGS). Characterization of cancer-relevant features including driver mutations or cellular processes was performed using mutational and gene expression data in order to identify potential biomarker or treatment targets in glioblastoma. In summary, we reported a newly established and molecularly characterized panel of glioma PDX models with high relevance for translational preclinical research.

  Study EGAS00001007119

• Evolution and clinical impact of genetic epistasis within EGFR-mutant lung cancers: multi-timepoint exome sequencing of a single patient's disease

  This dataset is a time-series of EGFR-mutant NSCLC clinical specimens from an individual patient profiled using tumor-based whole exome sequencing and the data is in BAM format. DNA was extracted from FFPE for primary tumor and frozen tumor tissue samples and matched non-tumor tissue using the Qiagen Allprep DNA/RNA Mini Kit.  The library preparation protocol was based on the Agilent SureSelect Library Prep and Capture System. DNA was resuspended in a low TE buffer and sheared (Duty Cycle 5%; Intensity 175; Cycles/Burst: 200; Time: 300s, Corvaris S2 Utrasonicator).  Bar-coded exome libraries were prepared using the Agilent Sure Select V5 library kit per manfucaturer’s specifications. The libraries were run on the HiSeq2500. Raw paired end reads (100bp) in FastQ format generated by the Illumina pipeline were aligned to the full hg19 genomic assembly obtained from USCS, gencode 14, using bwa version 0.7.12. Picard tools version 1.117 was used to sort, remove duplicate reads and generate QC statistics. Tumor DNA was sequenced to median depth of 303X (range 114.39-383.41) and the matched germline DNA to average depth of 231.65.

  Dataset EGAD00001003769

• Single-cell level characterization of B cell depletion and repopulation following rituximab in systemic lupus erythematosus

  Rituximab, a CD20+ B cell depletion therapy, is frequently used in the treatment of systemic lupus erythematosus (SLE). However, variability in patient response highlights the need for a deeper understanding of the underlying immune cell dynamics of B cell depletion and repopulation. In this study, we conducted longitudinal single-cell profiling of nine SLE patients treated with rituximab from pretreatment to up to 15 months post-treatment. These were compared to eight healthy controls. We profiled PBMCs via 10X Genomics single-cell RNA, surface protein (CITE-seq), B cell receptor (BCR), and T cell receptor (TCR) sequencing and sequenced bulk BCR repertoires in parallel. For single cell sequencing, 10 pools were created with an equal number of cells from 4 samples each. Samples collected at different timepoints from the same patient were distributed across different pools to allow demultiplexing of individuals using genotypes captured from the scRNA-seq data. Libraries were pooled using a ratio of GEX:CITE:TCR:BCR=9:2:1:1 and were sequenced across two lanes of the NovaSeq 6000. For bulk BCR sequencing, individual samples were demultiplexed using primer barcodes and paired-end BCR amplicon reads were merged prior to submission; data are provided as unmapped single-end reads.

  Dataset EGAD00001015817

• Analysis of Donor Pancreata Defines the Transcriptomic Signature and Microenvironment of Early Neoplastic Pancreatic Lesions

  The adult healthy human pancreas has been poorly studied given lack of indication to obtain tissue from the pancreas in the absence of disease and rapid postmortem degradation. We obtained pancreata from brain dead donors thus avoiding any warm ischemia time. We found that neoplastic lesions occur frequently in healthy organs, in donors as young as in their 3rd decade of life. Using the 10X Genomics Platform, we performed single cell RNA sequencing on 7 donor pancreata, 6 of which we separately sequenced tissue from the pancreas head and tail, for a total of 11 sequencing runs. Of note, because one of our runs (containing sequencing from pancreas head and tail) did not pass qc, our final analysis in our study utilized 6 donor pancreata, 5 of which contained pancreas head and tail sequenced separately. This is a robust single-cell dataset of healthy, nonpathologic pancreas tissue and includes acinar, ductal, and non-epithelial cells (myeloid cells, fibroblasts, endothelial cells, lymphocytes).Briefly, the donor pancreas was dissected by a clinically-trained pancreatobiliary surgeon. The organ was dissected into head, body, and tail, with portions of each placed into DMEM media with 1% BSA/10µM Y27632. Tissue was minced into 1mm3 pieces, then digested with 1 mg/mL collagenase P for 20-30 min at 37°C with gentle agitation. Digested tissue was rinsed three times with DMEM/1% BSA/10µM Y27632, then filtered through a 40µm mesh. Resulting cells were submitted to the University of Michigan Advanced Genomics Core for single cell sequencing using the 10x Genomics Platform. FASTQ files of all sequencing runs along with clinical annotation of donors is available through this dbGaP submission.

  Study phs003229

• Exploratory novel biomarker and resistance mechanism of milademetan, an MDM2 inhibitor, in MDM2 amplified intimal sarcoma from an open-label phase 1b/2 trial ���NCCH1806/MK004���

  Amplified Murine double minute 2 (MDM2) is found in > 70% of intimal sarcoma, known as one of the ultra-rare sarcomas. Milademetan (DS3032, RAIN-32) is a novel, specific, small-molecule MDM2 inhibitor that disrupts MDM2 and the tumor suppressor protein p53 interactions in tumor cells. We conducted a phase 1b/2 trial (Trial registration No: JMA-IIA00402) in patients with amplified MDM2 wild-type TP53 intimal sarcoma as a sub-study under the nationwide large registry for rare cancers in Japan (MASTERKEY Project). Eleven patients were enrolled, and ten were included in the efficacy analysis. Two (20%) patients had durable responses for > 15 months. Milademetan provided clinical benefits in patients with amplified MDM2 intimal sarcoma. Predictive biomarkers other than amplified MDM2 and acquired resistance mechanisms for milademetan are unknown. Whole-exome and RNA sequencing analyses of pre-treatment tissue samples were conducted to identify determinants of response. Genomic alterations were analyzed for 10 patients, and gene expression was analyzed for 9 patients using their pre-treatment tissue samples. Targeted sequencing of cell-free DNA (cfDNA) samples (liquid biopsy) was also conducted sequentially at three points [before treatment with milademetan (baseline), at Cycle 2 Day1, and at the time of disease progression] to identify determinants of response and resistance. From whole-exome and RNA sequencing analyses of pre-treatment tissue samples, we could not find any molecular pathways associated with the anti-tumor activity of milademetan. Focusing on 8 genes (CDK4, CDKN2A, CDKN2B, EGFR, ERBB3, MDM2, PDGFRA, TP53) known to be frequently affected in intimal sarcoma and 10 genes (AKT1, ATM, BBC3, CDKN1A, CDKN1C, CHEK2, MDM4, PMAIP1, PPM1D, TWIST1) reported to be associated with MDM2 inhibitor responses; we found that anti-tumor activity correlated with amplified TWIST1 (p-value = 0.028) and negatively with CDKN2A loss (p-value = 0.071). Eight of the 10 patients had their cfDNA collected sequentially at baseline, at Cycle 2 Day1 and disease progression; however, one did not consent to the exploratory analysis study using cfDNA, and one had cfDNA collected at baseline but not at disease progression due to ongoing treatment. Of the eight patients, TP53 mutations in cfDNA were detected in one and five patients at baseline and disease progression, respectively. The cfDNA allele frequency of TP53 mutations increased with disease progression. CDKN2A loss and amplified TWIST1 could be associated with the anti-tumor activity of milademetan in patients with MDM2 amplified intimal sarcoma. Acquired TP53 mutations were detected in sequential liquid biopsies as loss-of-function mutations, and these TP53 mutations might compromise the anti-tumor activity.

  Study JGAS000619

• Single cell sequencing of breast cancer T cells reveals a tissue-resident memory subset associated with improved prognosis

  The quantity of tumor infiltrating lymphocytes (TILs) in breast cancer (BC) is a robust prognostic factor for improved patient survival particularly in triple negative and HER2-overexpressing BC subtypes. Whilst T cells are the predominant TIL population, the relationship between quantitative and qualitative differences in T cell subpopulations and patient prognosis remains unknown. We performed single-cell RNA sequencing (scRNA-seq) of 6311 T cells isolated from human BCs, and show for the first time that significant heterogeneity exists in the infiltrating T cell population. We demonstrated that high TIL BCs contained significant numbers of CD8+ T cells with features of tissue-resident memory T (TRM) cell differentiation and that these CD8+ TRM cells expressed high levels of immune checkpoint molecules and effector proteins. A CD8+ TRM gene signature developed from the scRNA-seq data was significantly associated with improved patient survival in early-stage triple negative breast cancer (TNBC) and provided better prognostication compared with CD8 expression alone. Our data suggests that CD8+ TRM contribute to BC immunosurveillance and are the key targets of modulation through immune checkpoint inhibition. Further understanding of the development, maintenance and regulation of TRM cells will be crucial for successful immunotherapeutic development in BC.

  Study EGAS00001002845

• Whole-Genome Sequencing Suggests Schizophrenia Risk Mechanisms in Humans with 22q11.2 Deletion Syndrome

  Chromosome 22q11.2 microdeletions impart a high but incomplete risk for schizophrenia. Possible mechanisms include genome-wide effects of DGCR8 haploinsufficiency. In a proof-of-principle study to assess the power of this model, we used high-quality, whole-genome sequencing of nine individuals with 22q11.2 deletions and extreme phenotypes (schizophrenia, or no psychotic disorder at age >50 years). The schizophrenia group had a greater burden of rare, damaging variants impacting protein-coding neurofunctional genes, including genes involved in neuron projection (nominal P = 0.02, joint burden of three variant types). Variants in the intact 22q11.2 region were not major contributors. Restricting to genes affected by a DGCR8 mechanism tended to amplify between-group differences. Damaging variants in highly conserved long intergenic noncoding RNA genes also were enriched in the schizophrenia group (nominal P = 0.04). The findings support the 22q11.2 deletion model as a threshold-lowering first hit for schizophrenia risk. If applied to a larger and thus better-powered cohort, this appears to be a promising approach to identify genome-wide rare variants in coding and noncoding sequence that perturb gene networks relevant to idiopathic schizophrenia. Similarly designed studies exploiting genetic models may prove useful to help delineate the genetic architecture of other complex phenotypes.

  Study EGAS00001002344

• Single-Cell ATAC and RNA Sequencing of Human Breast Cancer Reveals Salient Cancer-Specific Enhancers

  Matched single-cell chromatin accessibility sequencing (scATAC-seq) and single-cell transcriptomic sequencing (scRNA-seq) were performed for 12 human breast tumor samples, 4 normal human mammary reduction samples, and 4 human breast cancer cell lines (HCC1143, SUM149PT, MCF7, and T47D). Both scATAC-seq and scRNA-seq assays resolved the underlying cellular heterogeneity in chromatin accessibility and transcriptional output for each patient sample. In this work, we quantitatively linked variation in chromatin accessibility to gene expression across malignant and non-malignant cell types from patient samples. We then showed that our findings were recapitulated in the breast cancer cell line samples. These data will serve as an important resource for the single-cell genomics and breast cancer research communities.

  Study phs003253

• Immune Responses in Checkpoint Myocarditis Across Heart, Blood, and Tumor

  Immune checkpoint inhibitors (ICIs) are widely used in anti-cancer therapies, but they can cause morbid and potentially fatal immune-related adverse events (irAEs), such as ICI-related myocarditis (irMyocarditis). The pathogenesis of irMyocarditis and its relationship to anti-tumor immunity remain poorly understood. We sought to define immune responses in heart, tumor, and blood in irMyocarditis patients and controls by leveraging single‐cell RNA sequencing (scRNA‐seq) coupled with T-cell Receptor Sequencing (TCR-Seq). Our analysis demonstrated increased frequencies of cytotoxic T cells, inflammatory mononuclear phagocytes (MNPs), conventional dendritic cells (cDCs), and inflammatory fibroblasts in irMyocarditis heart tissue. Additionally, we revealed decreased frequencies of plasmacytoid dendritic cells, cDCs, and B lineage cells but an increased frequency of MNPs in the blood of irMyocarditis patients. Raw data for heart and blood can be found in dbGaP. All processed data and raw TCRseq data from heart tissue and tumor samples for this study can be found in GEO (Accession number GSE228597).Subject IDs with the prefix "SIC" were collected by the Severe Immunotherapy Complications (SIC) service at MGH due to suspicion of immune-related adverse events. Subject IDs with the prefix "donor" were collected by the MGH Melanoma Biobank at scheduled timepoints; this biobank systematically collects samples from patients at MGH at clinically relevant timepoints.

  Study phs003413

• A Perioperative Study of Safusidenib in Patients with IDH1-Mutated Glioma

  10 participants with a radiological or pathological diagnosis of WHO Grade 2 lower-grade glioma (LGG) were recruited at a single center starting in December, 2022. Eligibility criteria permitted participants who had undergone prior surgery for WHO Grade 2 LGG, but excluded those with a history of prior radiation or chemotherapy.In Part A of the study, participants underwent an initial craniotomy with tumor biopsy and a lumbar puncture for cerebrospinal fluid (CSF) collection at study entry. Following sample collection, participants received safusidenib at the recommended phase 2 doses of 250 mg twice daily (b.i.d.) for a period of 22–36 days. This was followed by a second craniotomy for maximal safe tumor resection and a repeat lumbar puncture within 5 hours of the final pre-resection dose of safusidenib. Tumor, blood, and CSF samples were collected at both surgical time points for comprehensive translational research.Molecular analyses were performed on samples collected at both time points. These included pharmacokinetics (PK), targeted metabolomics, whole-genome sequencing (WGS) of tumor tissue, single-nucleus RNA sequencing (snRNA-seq), and DNA methylation profiling on all patients. Targeted metabolomics was conducted on samples from 9 participants, while spatial transcriptomics and spatial metabolomics were performed on tumor tissue from 3 participants each. Histone methylation profiling was conducted on samples from 8 participants.

  Study phs003976

• Integrated analysis of relapsed B-cell precursor Acute Lymphoblastic Leukemia identifies subtype-specific cytokine and metabolic signatures

  Recent efforts reclassified B-Cell Precursor Acute Lymphoblastic Leukemia (BCP-ALL) into more refined subtypes. Nevertheless, outcomes of relapsed BCP-ALL remain unsatisfactory, particularly in adult patients where the molecular basis of relapse is still poorly understood. To elucidate the evolution of relapse in BCP-ALL, we established a comprehensive multi-omics dataset including DNA-sequencing, RNA-sequencing, DNA methylation array and proteome MASS-spec data from matched diagnosis and relapse samples of BCP-ALL patients (n = 50) including the subtypes DUX4, Ph-like and two aneuploid subtypes. Relapse-specific alterations were enriched for chromatin modifiers, nucleotide and steroid metabolism including the novel candidates FPGS, AGBL and ZNF483. The proteome expression analysis unraveled deregulation of metabolic pathways at relapse including the key proteins G6PD, TKT, GPI and PGD. Moreover, we identified a novel relapse-specific gene signature specific for DUX4 BCP-ALL patients highlighting chemotaxis and cytokine environment as a possible driver event at relapse. This study presents novel insights at distinct molecular levels of relapsed BCP-ALL based on a comprehensive multi-omics integrated data set including a valuable proteomics data set. The relapse specific aberrations reveal metabolic signatures on genomic and proteomic levels in BCP-ALL relapse. Furthermore, the chemokine expression signature in DUX4 relapse underscores the distinct status of DUX4-fusion BCP-ALL.

  Study EGAS00001002856

• Genetic characterization of B-cell prolymphocytic leukemia: a hierarchical prognostic model involving MYC and TP53 abnormalities - RNA-seq

  B-cell prolymphocytic leukemia (B-PLL) is a rare disease, whose molecular pathology is largely unknown. We report here the cytogenetic and molecular findings in a large series of 34 B-PLL. Karyotype (K) was complex (≥3 abnormalities) in 73%, and highly complex (HCK≥5) in 45%. The most frequent chromosomal aberrations were: translocation targeting the MYC gene [t(MYC)] (62%), 17p deletion including TP53 gene (38%), trisomy 18/18q (30%), 13q14 deletion (29%), trisomy 3 (24%), trisomy 12 (24%) and 8p deletion (23%). Whole-Exome Sequencing performed in 16 patients revealed recurrent mutations in TP53 (6/16, 38%), MYD88 (n=4), BCOR (n=4), MYC (n=3), SF3B1 (n=3), FAT1 (n=3), SETD2 (n=2), CHD2 (n=2), CXCR4 (n=2) and BCLAF1 (n=2). The main group of patients (21/34, 62%) had a t(MYC) associated with a higher percentage of prolymphocytes (p=0.03), CD38 expression (p<0.001), lower K complexity (p=0.0004), mutations in MYC and in genes involved in RNA metabolism and chromatin remodeling. Principal component analysis of gene expression data showed that patients with t(MYC) clustered together. A second group with MYC gain (5/34, 15%), was associated with HCK≥5 (p=0,01) and trisomy 3 (p=0,008). Altogether, 26/34 patients (76%) had a MYC activation, translocation or gain, that were mutually exclusive. We identified 3 distinct cytogenetic prognostic groups (p=0.0006): lower risk: absence of MYC activation (median not reached); intermediate risk: MYC activation without del17p (125 months); high risk: MYC activation + del17p (11 months). Our results show that cytogenetic analysis is a useful diagnostic tool in B-PLL that improves prognostic stratification.

  Study EGAS00001003274

• Genomic Analysis of Prostate Tumor Heterogeneity in Metastasis

  To investigate the molecular changes, mutations, and signaling pathways that characterize the initial steps of prostate cancer metastasis to the lymph nodes, we analyzed 19 discrete lymph node (LN) metastases, 97 primary tumor foci, 39 benign lymph nodes, and 10 normal adjacent prostate tissue samples from 43 patients enrolled in a clinical trial (NCT01808222) of extended lymphadenectomy and metastasis excision guided by [18F]-fluciclovine PET imaging for high risk, newly diagnosed prostate cancer conducted at Emory University. These patients fell within high or very high-risk groups (T3a, Gleason score 8-10, or PSA greater than 20 ng/ml). Of these patients, 16 had LN metastases at least 4 mm in diameter, and 20 patients had no metastases. For metastatic patients, we performed RNA and Whole Exome Sequencing (WES) on three primary tumor foci, one benign LN, and all LN metastases greater than 4mm in diameter. For non-metastatic patients, we sequenced two primary tumor foci and one benign LN. In total, RNAseq was performed on 165 tissue samples from 43 patients, and WES was performed on 137 samples from 35 patients.

  Study phs003404

• Harnessing Epigenetic Regulators to improve HSC-based lentiviral gene therapy

  Curative benefits of autologous and allogeneic transplantation of hematopoietic stem cells (HSCs) have been proven for various diseases. However, the low number of true HSCs that can be collected from patients and subsequently in vitro maintenance and expansion of true HSCs for genetic correction remain challenging. Addressing this issue, we here focused on optimizing culture conditions to improve the ex vivo expansion of true HSCs for gene therapy purposes. In particular, we explore the use of epigenetic regulators to enhance the effectiveness of HSC-based lentiviral (LV) gene therapy. The HDAC inhibitor Quisinostat and the bromodomain inhibitor CPI203 each promote ex vivo expansion of functional HSCs, as validated by xenotransplantation assays and single cell RNA-sequencing analysis. We confirmed the stealth effect of LV transduction on the loss of HSC numbers in commonly used culture protocols, while addition of Quisinostat or CPI203 improved expansion of HSCs in transduction protocols. Of note, we demonstrated that addition of Quisinostat improved LV transduction efficiency of HSCs and early progenitors. Our suggested culture conditions highlight the potential therapeutic effect of epigenetic regulators in hematopoietic stem cell biology and their clinical applications to advance HSC-based gene correction.

  Study EGAS50000000175

• Programs, Origins, and Immunomodulatory Functions of Myeloid Cells in Gliomas

  Gliomas are challenging malignancies characterized by an immunosuppressive microenvironment filled with myeloid cells, whose roles in immune modulation are not fully understood. In this study, we utilized single-cell RNA sequencing (scRNA-seq) data from 85 primary human glioma tumors, related blood samples, and brain tumor organoid models. We found that nearly all glioma-associated myeloid cells express at least one of four distinct immunomodulatory activity programs. These programs include:1. Scavenger Immunosuppressive Program2. Complement Immunosuppressive Program3. Microglial Inflammatory Program4. Systemic Inflammatory ProgramAll four programs are found in both lower-grade and high-grade gliomas and are expressed across various myeloid cell types, whether derived from blood or resident origins. We validated the remarkable plasticity of myeloid cells by showing that the introduction of peripheral blood monocytes into glioma organoid models can induce the emergence of identities resembling microglia, macrophages, and dendritic cells, as well as all four immunomodulatory programs. Our study provides a valuable resource and a new framework for understanding the role of immunomodulatory myeloid cells in gliomas, laying the groundwork for developing effective immunotherapy strategies for glioma patients.

  Study phs003756

• Defining Cutaneous Gene Expression Signatures in Juvenile Dermatomyositis

  Skin inflammation in juvenile dermatomyositis (JDM) can signal disease onset or flare up and prevent complete disease remission. Our objective was to evaluate JDM cutaneous expression signatures through tape stripping in both lesional and non-lesional skin and as compared to paired blood. A total of 28 JDM and 20 control (CTL) patients participated in the study. All patients had tape strip sampling of non-lesional (NL) skin at enrollment visit. Within our JDM cohort, paired lesional (L) skin was additionally sampled if rash was present. 16/28 of JDM patients had at least one additional skin sample associated with a follow-up (FU) visit (total number of follow-up samples included = 26 NL and 4 L). We identified a common innate immune signature in both lesional and non-lesional JDM skin. Using unsupervised hierarchical clustering of skin expression data, we distinguished a JDM molecular subgroup characterized by a mitochondrial dysfunction signature. In comparison to blood expression data, we determined that skin may more effectively highlight the JDM subgroup characterized by a mitochondrial dysfunction expression signature. Data available include bulk RNA-sequencing data from tape strips and blood of the JDM and control patients from this study.

  Study phs003884

• Single cell atlas of large B-cell lymphoma

  Large B-cell lymphomas (LBCL) are clinically and biologically heterogeneous lymphoid malignancies, shaped by complex tumor and microenvironments that play a critical role in disease progression. There has been no direct and comprehensive single-cell characterization of tumor microenvironment in LBCL that simultaneously captures lymphoid, myeloid, and stromal populations across a sufficiently large cohort to reveal reproducible and biologically meaningful patterns. Here, we performed single-nucleus multiome profiling, coupled with bulk DNA and RNA sequencing, on 232 tumor and control biopsies to comprehensively map cell subpopulations of LBCL. Through this integrative approach, we identified Lymphoma Microenvironment Archetype Profiles (LymphoMAPs)—distinct, reproducible patterns defined by the co-occurrence of functionally specialized cell subsets. We further characterized the cell-cell communication pathways that shape these transcriptional phenotypes within archetypes. Importantly, we discovered associations between specific LymphoMAPs and clinical and molecular features of LBCL, including response to CD19 CAR T-cell therapy. Our study provides a valuable resource of lymphoid, myeloid, and non-hematopoietic cell states in LBCL and offers a framework for understanding how the abundance and functional state of cells within the tumor microenvironment may drive inter-patient variability in therapeutic outcome.

  Study EGAS50000001022

• Omic Studies in Children's Respiratory and Environmental Workgroup (CREW) Cohorts

  Children's Respiratory Research Workgroup (CREW) is an NIH-funded project consisting of multiple individual U.S. cohorts and three scientific centers working together to identify phenotypes and causes of childhood asthma. CREW will include data from a large number of children and their families, with broad diversity in terms of ethnicity, family characteristics, neighborhoods and geographic locations. One of the primary goals of CREW is to put together sets of data and samples of participating cohorts to identify phenotypes of childhood asthma (i.e. specific subtypes of asthma that can be distinguished by clinical features such as natural history, triggers, exacerbation frequency, concurrent allergies, lung function, sex, etc). Omics studies will include genotyping, methylation, and RNA sequencing. By integrating the DNA genotyping with the expression, immunologic and clinical data, we will identify eQTL, inflammation QTL, and lung function QTL, as well as QTL associated with other biomarkers. Analyses will identify genetic, epigenetic, transcriptional and protein/metabolite profiles that are associated with specific exposures or phenotypes. As we obtain mechanistic insights about personal and early life risk factors, we will connect asthma phenotypes with underlying causes and pathogenic mechanisms to define endotypes of childhood asthma.

  Study phs004036

• Biallelic loss of function PAX4 variants are a Novel Cause of Transient Neonatal Diabetes

  The identification of genes in which loss of function (LoF) variants cause neonatal diabetes mellitus (NDM) can provide unique insights into the development and function of human pancreatic beta cells. We describe the identification of biallelic PAX4 LoF variants in 2 unrelated individuals with transient NDM: one due to a homozygous p.(Arg126*) PAX4 stop gain variant and a second due to a homozygous c.-352_104del 2.65kb deletion affecting the first 4 exons and promoter of the PAX4 Matched Annotation from NCBI and EMBL-EBI select transcript. We confirmed reduced transcript expression due to nonsense mediated decay for the e p.(Arg126*) variant in CRISPR-edited iPSC-derived pancreatic endoderm cells. Analysis of PAX4 binding through integration of CUT&RUN and RNA-sequencing in PAX4 depleted endocrine cells identified genes directly regulated by PAX4 involved in both pancreatic islet development and glucose-sensitive insulin secretion. The remitting and relapsing diabetes in both probands demonstrates that in contrast to mice PAX4 is not essential for the development of pancreatic beta-cells and supports differences in the requirements of functional beta-cell mass across the life-course.

  Study EGAS50000000857

• Germline mutations in the transcription factor IKZF5 cause thrombocytopenia.

  To identify novel causes of hereditary thrombocytopenia, we performed a genetic association analysis of whole-genome sequencing data from 13,037 individuals enrolled in the NIHR BioResource, including 233 cases with isolated thrombocytopenia. We found an association between rare variants in the transcription factor encoding gene IKZF5 and thrombocytopenia. We report five causal missense variants in or near IKZF5 zinc fingers (Znfs), of which two occurred de novo and three co-segregated in three pedigrees. A canonical DNA-Znf binding model predicts that three of the variants alter DNA recognition. Expression studies showed that chromatin binding was disrupted in mutant compared to wild-type (WT) IKZF5 and electron microscopy revealed a reduced quantity of alpha granules in normally sized platelets. Proplatelet formation (PPF) was reduced in megakaryocytes (MKs) from seven cases relative to six controls. Comparison of RNA-seq data from platelets, monocytes, neutrophils and CD4 T-cells from three cases and 14 healthy controls showed 1,194 differentially expressed genes in platelets but only four DEGs in each of the other blood cell types. In conclusion, IKZF5 is a novel transcriptional regulator of megakaryopoiesis and the eighth transcription factor associated with dominant thrombocytopenia in humans.

  Study EGAS00001003736

• Oncoprint GSCCs

  Glioblastoma (GBM) remains the most malignant primary brain tumor, with a median survival rarely exceeding 2 years. Tumor heterogeneity and an immunosuppressive microenvironment are key factors contributing to the poor response rates of current therapeutic approaches. GBM-associated macrophages (GAMs) often exhibit immunosuppressive features that promote tumor progression. However, their dynamic interactions with GBM tumor cells remain poorly understood. Here, we used patient-derived GBM stem cell cultures and combined single-cell RNA sequencing of GAM-GBM co-cultures and real-time in vivo monitoring of GAM-GBM interactions in orthotopic zebrafish xenograft models to provide insight into the cellular, molecular, and spatial heterogeneity. Our analyses revealed substantial heterogeneity across GBM patients in GBM-induced GAM polarization and the ability to attract and activate GAMs – features that correlated with patient survival. Differential gene expression analysis, immunohistochemistry on original tumor samples, and knock-out experiments in zebrafish subsequently identified LGALS1 as a primary regulator of immunosuppression. Overall, our work highlights that GAM-GBM interactions can be studied in a clinically relevant way using co-cultures and avatar models, while offering new opportunities to identify promising immune-modulating targets.

  Study EGAS00001007481

• Biased allelic expression in human primary fibroblast single cells.

  The study of gene expression in mammalian single cells using genomic technologies now provides the possibility to investigate the patterns of allelic gene expression. We have used single-cell RNA sequencing to detect the allele-specific mRNA level in 203 single human primary fibroblast cells over 133,633 unique heterozygous single nucleotide variants (hetSNVs). We have observed that at the snapshot of analyses each cell contains mostly transcripts from one allele from the majority of genes; indeed 76.4% of the hetSNVs display stochastic monoallelic expression in single cells. Remarkably, adjacent hetSNVs exhibit haplotype consistent allelic ratio; in contrast distant sites located in two different genes are independent of the haplotype structure. Moreover, the allele-specific expression in single cells correlated with the abundance of the cellular transcript. We observed that genes expressing both alleles in majority of the single cells at a given time point are rare and enriched in highly expressed genes. The relative abundance of each allele in a cell is controlled by some regulatory mechanisms since we observed related single-cell allelic profiles according to genes. Overall, these results have direct implications in cellular phenotypic variability.

  Study EGAS00001001009

• Germline biallelic mutation affecting the transcription factor Helios causes pleiotropic defects of immunity

  We studied a patient with recurrent respiratory infections and hypogammaglobulinemia and identified a germline homozygous missense mutation in IKZF2 encoding Helios (p.Ile325Val). We show that HeliosI325V retains DNA-binding and dimerization properties, but loses interaction with several partners, including epigenetic remodelers HDAC1, HDAC3 and the ATAC complex. Single-cell RNA-sequencing of peripheral blood mononuclear cells revealed gene expression signatures indicative of a shift towards pro-inflammatory, effector-like status in the patient’s T cells. We observed an upregulation of the anti-inflammatory gene, TSC22D3, encoding the glucocorticoid-induced leucine zipper (GILZ), whose expression was reported to increase upon IL-2 deprivation. We validated this finding in primary T cells, where we observed a pronounced reduction in IL-2 production upon their stimulation, together with a defect in their ability to proliferate, rendering them more anergic. Collectively, we identify a novel germline-encoded inborn error of immunity and define a role for Helios in conventional T cells, whereby interactions with specific binding partners is necessary to mediate the transcriptional programs that enable T-cell homeostasis in health and disease.

  Study EGAS00001005675

• Molecular_risk_stratification_in_patients_with_T1_colorectal_cancer_WES

  The introduction of bowel cancer screening has led to a significant increase in the proportion of patients being diagnosed with asymptomatic, early-stage colorectal cancer (CRC). Although the majority of these patients are successfully treated with surgery alone, a small proportion of patients have ‘born-to-be-bad’ aggressive lesions with early dissemination leading to distant metastases. Current standard of care histological assessment is unable to distinguish between these aggressive versus non-aggressive early lesions which is essential to provide appropriate clinical management decisions. This study aims to carry out molecular and histological profiling of approximately 300 T1 CRCs in order to develop a molecular stratifier based on the risk of relapse in early-invasive CRC. This novel T1 cohort will represent the world’s largest molecularly characterised T1 cohort of samples, with digital pathology assessment alongside whole exome sequencing, copy number variation analysis and 3’ RNA-seq. This data will be used to generate a robust panel of molecular and/or histological markers applicable to formalin-fixed paraffin embedded (FFPE) archival tissue which discriminates between T1 lesions based on risk of relapse, which will ultimately be used to inform clinical management of CRC at the earliest stages of the disease.

  Study EGAS00001005734

• Identification of Germline Monoallelic Mutations in IKZF2 in Patients with Immune Dysregulation.

  Helios, encoded by IKZF2, is a member of the Ikaros family of transcription factors with pivotal roles in T-follicular helper, NK- and T-regulatory cell physiology. Somatic IKZF2 mutations are frequently found in lymphoid malignancies. Although germline mutations in IKZF1 and IKZF3, encoding Ikaros and Aiolos, have recently been identified in patients with phenotypically similar immunodeficiency syndromes, the effect of germline mutations in IKZF2 on human hematopoiesis and immunity remains enigmatic. We identified germline IKZF2 mutations (one nonsense (p.R291X)- and 4 distinct missense variants) in six patients with systemic lupus erythematosus, immune thrombocytopenia or EBV-associated hemophagocytic lymphohistiocytosis. Patients exhibited hypogammaglobulinemia, decreased number of T-follicular helper and NKcells. Single-cell RNA sequencing of PBMCs from the patient carrying the R291X variant revealed upregulation of pro-inflammatory genes associated with T-cell receptor activation and T-cell exhaustion. Functional assays revealed the inability of HeliosR291X to homodimerize and bind target DNA as dimers. Moreover, proteomic analysis by proximity-dependent Biotin Identification revealed aberrant interaction of 3/5 Helios mutants with core components of the NuRD complex conveying HELIOS mediated epigenetic and transcriptional dysregulation.

  Study EGAS00001005874

• Patient-tailored design for selective co-inhibition of leukemic cell subpopulations

  The extensive primary and secondary drug resistance in many cancer types requires rational approaches to design personalized and selective combinatorial therapies that do not only show synergistic effect in overall cancer cell killing but also result in minimal toxic side effects on non-malignant cells. To address the combinatorial explosion in the number of relevant combinations, we implemented a machine learning approach that prioritizes patient-customized drug combinations with a desired synergy-efficacy-toxicity balance by combining single-cell RNA-sequencing with ex vivo single-agent testing in scarce patient-derived primary cells. When applied to two diagnostic and two refractory AML patient cases, each with a different genetic background, our integrated approach accurately predicted patient-specific combinations that were shown to result not only in synergistic cancer cell co-inhibition but were also capable of targeting specific AML cell subpopulations that emerge in differing stages of disease pathogenesis or treatment regimens. Our data-driven approach provides an unbiased means for systematic identification of personalized combinatorial regimens that selectively co-inhibit leukemic cells while avoiding inhibition of non-malignant cells, and highlight the relevance of considering cell heterogeneity for personalized cancer therapy.

  Study EGAS00001004614

• Targeted therapy of advanced parathyroid carcinoma guided by genomic and transcriptomic profiling (hipo_021)

  Parathyroid carcinoma (PC) is an orphan malignancy with a high risk of recurrence after surgery. Tumor-directed systemic treatments for PC are not established. We used whole-genome and RNA sequencing in four patients with advanced PC to identify molecular alterations that could guide clinical management. In two cases, the genomic and transcriptomic profiles provided targets for experimental therapies that resulted in biochemical response and prolonged disease stabilization: (i) immune checkpoint inhibition with pembrolizumab based on high tumor mutational burden and a single-base substitution signature associated with APOBEC (apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like) overactivation; (ii) multi-receptor tyrosine kinase inhibition with lenvatinib due to overexpression of FGFR1 and RET and, later in the course of disease, PARP inhibition with olaparib prompted by signs of defective homologous recombination DNA repair. In addition, our data provided new insights into the molecular landscape of PC with respect to the genome-wide footprints of specific mutational processes and pathogenic germline alterations. These data underscore the potential of comprehensive molecular analyses to improve care for patients with ultra-rare cancers based on insight into disease biology.

  Study EGAS00001006747

• The impact of BNT162b2 mRNA vaccine against SARS-CoV-2 on adaptive and innate immune responses

  The mRNA-based BNT162b2 protects against severe disease and mortality caused by SARS-CoV-2 through induction of specific antibody and T-cell responses. Much less is known about its broad effects on immune responses against other pathogens. In the present study, we investigated the specific adaptive immune responses induced by BNT162b2 vaccination against various SARS-CoV-2 variants, as well as its effects on the responsiveness of human immune cells upon stimulation with heterologous viral, bacterial, and fungal pathogens. BNT162b2 vaccination induced effective humoral and cellular immunity against SARS-CoV-2 that started to wane after six months. We also observed long-term transcriptional changes in immune cells after vaccination, as assessed by RNA sequencing. Additionally, vaccination with BNT162b2 modulated innate immune responses as measured by the production of inflammatory cytokines when stimulated with various microbial stimuli other than SARS-CoV-2, including higher IL-1/IL-6 release and decreased production of IFN-α. Altogether, these data expand our knowledge regarding the overall immunological effects of this new class of vaccines and underline the need of additional studies to elucidate their effects on both innate and adaptive immune responses.

  Study EGAS00001006818

• Single-cell profiling of co-cultures of GSCCs and macrophages

  Glioblastoma (GBM) remains the most malignant primary brain tumor, with a median survival rarely exceeding 2 years. Tumor heterogeneity and an immunosuppressive microenvironment are key factors contributing to the poor response rates of current therapeutic approaches. GBM-associated macrophages (GAMs) often exhibit immunosuppressive features that promote tumor progression. However, their dynamic interactions with GBM tumor cells remain poorly understood. Here, we used patient-derived GBM stem cell cultures and combined single-cell RNA sequencing of GAM-GBM co-cultures and real-time in vivo monitoring of GAM-GBM interactions in orthotopic zebrafish xenograft models to provide insight into the cellular, molecular, and spatial heterogeneity. Our analyses revealed substantial heterogeneity across GBM patients in GBM-induced GAM polarization and the ability to attract and activate GAMs – features that correlated with patient survival. Differential gene expression analysis, immunohistochemistry on original tumor samples, and knock-out experiments in zebrafish subsequently identified LGALS1 as a primary regulator of immunosuppression. Overall, our work highlights that GAM-GBM interactions can be studied in a clinically relevant way using co-cultures and avatar models, while offering new opportunities to identify promising immune-modulating targets.

  Study EGAS00001007482

• SINGLE-CELL RNA SEQUENCING Single-cell RNA sequencing was performed on 13 ‘mild-moderate’ and 10 ‘critical’ COVID19 PBMC samples

  Summary Epidemiological and clinical reports have indicated that SARS-CoV-2 virulence hinges upon the triggering of an aberrant host immune response, more so than on direct virus-induced cellular damage. To elucidate the immunopathology underlying COVID-19 severity, we performed cytokine and multiplex immune profiling in mild-moderate and critically ill COVID-19 patients. Hypercytokinemia in COVID-19 differed from the IFN-γ-driven cytokine storm in macrophage activation syndrome, and was more pronounced in critical versus mild-moderate COVID-19. Systems modelling of cytokine levels paired with deep-immune profiling showed that classical monocytes drive this hyper-inflammatory phenotype and that a reduction in T-lymphocytes correlates with disease severity, with CD8+ cells being disproportionately affected. Expression of antigen presenting machinery was also reduced in critical disease. Furthermore, we found that neutrophils contributed to disease severity and local tissue damage by amplification of hypercytokinemia and the formation of neutrophil extracellular traps. Together these findings suggest a myeloid-driven immunopathology, in which hyperactivated neutrophils and an ineffective adaptive immune system act as mediators of COVID-19 disease severity.

  Study EGAS00001005039

• Fecal microbiome of T2D patients undergoing semaglutide or empagliflozin treatment

  The gut microbiome has been shown to be affected by the use of many human-targeted medications, and the interaction can be bidirectional. This has been clearly demonstrated for type 2 diabetes medications that have been in clinical use for several decades. However, the bidirectional effects of novel type 2 diabetes drugs semaglutide, empagliflozin, and the gut microbiome have yet to be clearly described. We investigated the effect of semaglutide and empagliflozin initiation on the gut microbiome of type 2 diabetes patients. In addition, we analyzed whether the pre-treatment gut microbiome can predict the treatment efficacy. Twenty subjects (10 women and 10 men) were enrolled in the study between November 2019 and January 2023 from the University of Tartu Hospital, Tartu, Estonia. All participants had been diagnosed with type 2 diabetes. The study group taking semaglutide included 6 women and 4 men, and the group prescribed empagliflozin included 4 women and 6 men. Gut microbiome fecal samples donated at four timepoints (Baseline, Month 1, Month 3; Month 12) were studied using 16S ribosomal RNA gene sequencing and analysis.

  Dataset EGAD50000000756

• Access to "A Spatially Resolved Single-Cell Atlas of the Human Fetal Olfactory System"

  The human nasal region arises from neural crest and placodal lineages, yet its early development remains poorly understood owing to limited fetal tissue access and structural complexity. Here we present an integrated single-nucleus and spatial transcriptomic atlas of the human fetal nasal region, generated from 10 fetuses (males and females) between 7 and 12 post-conceptional weeks. Single-nucleus RNA sequencing resolved 32 distinct cell types, while integration with multiplexed error-robust fluorescence in situ hybridization (MERFISH) enabled spatial and temporal mapping of gene expression dynamics across the olfactory epithelium (OE) and adjacent tissues. We identify novel markers of olfactory sensory neuron differentiation and pathways governing epithelial patterning and OE morphogenesis. Notably, spatially organized expression of 169 olfactory receptor genes reveals molecular evidence for the “one neuron-one receptor” principle during the first trimester. Together, this work establishes the first molecular and spatial framework of early human olfactory development and provides a foundational resource for studies of sensory neurogenesis and congenital disorders.

  Dac EGAC50000000688

• Single-cell and bulk RNA-sequencing of nivolumab-treated glioblastoma

  Glioblastoma is a highly aggressive type of brain tumour for which there is no curative treatment available. Immunotherapies have shown limited responses in unselected glioblastoma patients, and it is currently not well understood what separates the few responders from the non-responders at the cellular and molecular level. Here we investigate the phenotypic and transcriptional evolutionary dynamics at single-cell resolution during nivolumab immune checkpoint treatment of eight glioblastoma patients enrolled in a translational phase I/II clinical trial.

  Study EGAS00001007110

• Genetics of gene expression in human macrophage response to Salmonella

  Understanding how cells sense and respond to their environment, and how these responses are modulated by genetic variation, are fundamental biological problems, particularly for understanding how pathogenic organisms invade and manipulate the cells of the human immune system. Macrophages recognize and respond to many important human pathogens including HIV-1, Mycobacteria tuberculosis and Salmonella. This study will focus on the cellular response of human macrophages to Salmonella infection and how this response is modulated by the genetic bacground of the individual as well as additional pro-inflammatory stimulus (interferon-gamma priming). We will acquire 100 human induced pluripotent stem cell lines from the HipSci project, differentiate the cells in vitro into macrophages and expose them to four environmental conditions: (i) no stimulation, (ii) interferon-gamma (18h), (iii) Salmonella typhimurium SL1344 (5h), (iv) interferon-gamma (18h) + Salmonella (5h).Subsequently, we will isolate RNA from the samples for sequencing.

  Dataset EGAD00001003204

• Cellular Profiling Identifies Targetable T Cell Phenotypes in Lymphocytic Variant Hypereosinophilic Syndrome

  Lymphocytic variant hypereosinophilic syndrome (LHES) is driven by aberrant Th2‑skewed lymphocyte clones that promote tissue eosinophilia. We performed multimodal single‑cell RNA‑sequencing (scRNA‑seq), cell‑surface protein (ADT) profiling, and paired αβ‑TCR sequencing on peripheral‑blood mononuclear cells (PBMCs) from two clinically characterized LHES patients and one idiopathic HES control. One LHES patient was sequenced at two timepoints, pre-JAK inhibitor treatment and post-JAK inhibitor treatment. Single‑cell regulatory‑network inference (SCENIC), CellChat ligand–receptor modelling, and scRNA variant calling were integrated with flow‑cytometric validation and 5,000× amplicon sequencing. Two clonally expanded CD3‑CD4⁺ central‑/effector‑memory T‑cell clusters were observed exclusively in LHES. These cells showed Th2‑associated regulons (GATA3, BATF, MAF), high expression of CCR4, CD52, PTGDR2 and IL13, and down‑regulation of CD247 (CD3ζ), a mechanism experimentally confirmed to abolish surface CD3/TCR. Autocrine IL‑7 signaling, predicted by CellChat, provides an antigen‑independent proliferative cue. Therapeutically actionable targets (CCR4, CD52, IL‑5, IL‑7, TNF) were prioritized; one refractory patient achieved complete clinical and cytologic remission after mogamulizumab (anti‑CCR4) therapy. Longitudinal sequencing revealed clonal evolution with acquisition of a STAT3 p.G618R mutation under JAK‑inhibition, highlighting a resistance mechanism. Raw scRNA‑seq (5′ cDNA, VDJ, ADT) data, and sample‑level metadata will be distributed through dbGaP. These resources provide the first high‑resolution immune atlas of LHES and a framework for rational, precision immunotherapy in eosinophilic disorders.

  Study phs004041

• 13 WGS and 8 TCS data of ENKL patients treated with pembrolizumab.

  Whole-genome sequencing (WGS) was performed for 13 pairs of tumor-normal samples from patients diagnosed with NKTL. Genomic DNA from tumor tissue was extracted with QIAamp DNA Mini Kit. The DNA for the matching normal was obtained from blood or buccal swabs and purified by Blood and Cell Culture DNA Mini kit or E.Z.N.A. Tissue DNA Kit (Omega Bio-tek) according to manufacturer’s instructions. The quantity and quality were assessed by Quant-iT PicoGreen dsDNA Assay Kit (Invitrogen) and agarose gel electrophoresis. All sequencing libraries were prepared using TruSeq Nano DNA Library Prep Kit (Illumina). Paired-end sequencing was performed on Illumina HiSeq 2000 or HiSeq X Ten as 2x101 bp or 2x151 bp, respectively. 8 NKTL FFPE specimens were screened for somatic mutations using deep targeted capture sequencing (TCS). FFPE rolls or slides were extracted using QIAamp DNA FFPE Tissue kit (QIAGEN). The FFPE genomic DNA was treated with NEBNext FFPE DNA Repair Mix and assessed by Quant-it PicoGreen dsDNA Assay Kit (Invitrogen). The library was generated from 10-200 ng DNA with SureSelectXT Low Input Target Enrichment System for Illumina Paired-End Sequencing Library (Agilent Technologies) according to manufacturer’s instructions. RNA based probe was designed with SureDesign (Agilent Technologies) to target-capture 140 genes. Next, the captured libraries were pooled in equimolar concentration and sequenced on Illumina Novaseq 6000 platform with SP or S1 chip. Reads aligning to 40 selected genes were isolated post-alignment for this submission. Prefix used in filenames: T - Tumor samples N - Matched-Normal samples

  Dataset EGAD00001004140

• Using Genomics to Reduce Breast Cancer Disparities in the African Diaspora

  Differences in breast cancer incidence and mortality rates between North American Caucasian and African American women are well-described and transcend socioeconomic issues. Black women are diagnosed with breast cancer at a younger median age; have more clinically aggressive disease and stage-for-stage; and have higher mortality rates than age-matched Caucasian women. Black women in West Africa, the origin of the slave trade in the US in the 19th century and thus the founder population for most African Americans, have even higher rates of early-onset, poor-prognosis breast cancer than African American women. Racial difference in the distribution of intrinsic molecular subtypes has been well characterized in the US and throughout the African Diaspora as well. Despite the large efforts on characterizing racial/ethnic differences, however, the reasons women of African ancestry are disproportionately affected by breast cancer incidence and mortality remain poorly understood - largely due to paucity of data on inherent genomic differences that contribute to the disparities in incidence and progression of breast cancer across populations. West Africa Breast Cancer Study (WABCS) is an initiative that aims to comprehensively understand the genetic architecture of breast cancer in West Africans, the founder population of a large proportion of black women in the United States. The objective of the study was to provide a better understanding of the molecular genetic factors that influence prognosis in Nigerian breast cancer patients, and determine which of these alterations may be amenable to available therapy. To that end, we examined the molecular features of breast cancers of indigenous African women using a combination of whole-genome, whole-exome, and transcriptome sequencing (WGS, WES, and RNA-seq) on 194 tumors from Nigerian patients. The goal of this project was to obtain answers to two related research questions using an unscreened population without genetic admixture in Nigeria: 1) why are women of African ancestry more likely to develop aggressive young onset breast cancer? 2) What are the associated genomic and non-genomic risk factors? We hypothesize that the genomic determinants of breast cancer molecular subtypes in women of African ancestry are also molecular drivers of tumor progression and represent targets for interventions to improve clinical outcomes and close the mortality gap. By identifying causal links between genetic variants that promote aggressive tumor progression in Nigerian women in comparison to women from different population found in TCGA and ICGC, the present dataset will have significant public health impact on millions of women in the African Diaspora. The potential to identify novel pathways for interventions to reduce the increasing mortality gap between women of African and European ancestry is huge

  Study phs001687

• RNA-seq analysis of transcriptome variation with human ESC subclones

  Many studies over the past 10 years, culminating in the recent report of the International Stem Cell Initiative (ISCI, 2011) have shown that hPSC acquire genetic and epigenetic changes during their time in culture. Many of the genetic changes are non-random and recurrent, probably because they provide a selective growth advantage to the undifferentiated cells. Some are shared by embryonal carcinoma cells, the malignant counterparts of ES cells. The origins of these growth advantages are poorly understood, but may come from altered cell cycle dynamics, resistance to apoptosis or altered patterns of differentiation. Less is known about the nature and consequences of epigenetic changes, but it is likely that these similarly affect hPSC behaviour; e.g., enhanced expression of DLK1, an imprinted gene, is associated with altered hPSC growth (Enver et al 2005). Inevitably, these genetic and epigenetic changes will impact on our ability to use hPSC for regenerative medicine, either because malignant transformation of the undifferentiated cells or their differentiated derivatives to be used for transplantation compromises safety, or because they impede the function of those differentiated derivatives, or because they affect the efficiency with which the undifferentiated cells can be expanded and differentiated into desired cell types. Focusing initially upon the existing clinical grade hESC lines, later moving to iPSC, we will Consolidate and extend knowledge of the rate, type and functional impact of the genetic variations that occur during hPSC culture. We will use whole genome and exome sequencing as well as SNP arrays, together with clonal analysis and other cytogenetics techniques. Common changes will be compared with those found in the normal human population, at low frequency in the original cell population or observed during iPSC generation in the HIPSCI project currently based at the WTSI. These studies will provide a better understanding of the range of genetic changes that occur in hPSC beyond the CNVs already identified. In conjunction with cancer genome resources and expertise at WTSI, bioinformatic analyses of these hPSC data will allow us to assess potential impact on hPSC behaviour pertinent to applications in regenerative medicine, notably the likelihood that specific changes arising in undifferentiated PSC cultures may be associated with potential malignant transformation of differentiated progeny 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/ .

  Dataset EGAD00001004266

• Single-cell RNA-sequencing of CD19 Chimeric Antigen Receptor (CAR) T-cells (CD19 CAR T-cells) from a Phase I clinical study in paediatric ALL: CARPALL

  Chimeric antigen receptor (CAR)-modified T-cells have become established as an effective treatment of haematological cancers. In the context of relapsed and refractory childhood pre-B cell acute lymphoblastic leukaemia (B ALL), CD19 targeting CAR T-cells often induce durable remissions. Previously, we generated a novel low-affinity CAR incorporating a CD19-specific single-chain variable fragment (scFV) called CAT, displaying a faster off-rate of interaction than the FMC63 CD19 binder used in prior clinical studies. Here, we systematically analysed CD19 CAR T-cells of ten children with relapsed or refractory B ALL enrolled in the CARPALL trial (NCT02443831). To characterize persisting CD19 CAR T-cells, we performed high throughput single-cell gene expression and T-cell receptor (TCR) sequencing of infusion products and serial blood and bone marrow samples up to five years post-infusion. We isolated CAR T-cells from peripheral blood or bone marrow by flow cytometry for CD3 and CAR expression, prior to single cell sequencing (Chromium 10X) platform.

  Dataset EGAD00001010018

• Whole genome sequencing data from tumor and normal samples

  Whole genome sequencing FASTQ files from tumor and normal samples from 202 patients. All experiments reported in this article were performed in the Clinical Laboratory Improvement Amendments (CLIA)-certified and College of American Pathologists (CAP)-accredited laboratories at Personalis Inc., as guided by the Association for Molecular Pathology (AMP) and CAP’s joint recommendations52. Preferred input material consisting of FFPE slides (enabling H&E and macrodissection which is part of the standard NeXT Personal process) was not available for the majority (94%) of tumor samples. Most tumor samples (72%) were processed from FFPE curls, with the remaining processed from pre-extracted DNA and fresh frozen material. Where FFPE material on slide was available, tumor sections were macrodissected to increase tumor content, with a minimum allowable cellularity threshold of 20% determined by pathological review. The tumor cellularity threshold of 20% was chosen based on prior work demonstrating that no significant correlation is observed between tumor purity and assay limit of detection (LOD) at this threshold31. Genomic DNA was isolated from tumor and normal samples and WGS libraries prepared as previously described31. Genomic DNA was isolated from matched tumor and normal samples using the Qiagen AllPrep DNA/RNA FFPE Tissue Kit or the QIAamp DNA Mini Kit (QIAGEN, Germantown, MD, USA) using internally optimized workflows. WGS sequencing libraries were prepared with 120-500 ng of acoustically sheared genomic DNA (Covaris LLC, Woburn, MA, USA) using the KAPA HyperPrep Kit (Roche Sequencing Solutions, Pleasanton, CA, USA) and customized methods. Libraries were cleaned-up using AMPure XP beads and then quantified using the KAPA Library Quantification Kit (Roche Sequencing Solutions, Pleasanton, CA, USA), before being sequenced to 30X depth of coverage using a NovaSeq 6000 instrument (Illumina, San Diego, CA, USA).

  Dataset EGAD50000001909

• Genomic Analysis of Extra-Nodal Natural Killer/T-Cell Lymphoma (ENKTCL)

  In this study, we performed extensive genetic profiling on 209 extra-nodal natural killer/T-cell lymphoma (ENKTCL) nasal-type cases, including 174 Chinese from four hospitals in China, and 31 Hispanic and 4 Caucasian patients from City of Hope and Eberhard Karls University of Tübingen. All of the cases were reviewed by experienced hematopathologists to confirm diagnosis. Written informed consents from patients were reviewed and approved by the institutional review board of each institution. We deposited targeted exon sequencing data of a selected 333-gene panel on the ENKTCL cases in dbGaP. DNA was extracted from formalin-fixed paraffin-embedded (FFPE) tissues with an AllPrep DNA/RNA FFPE Kit. A custom platform (Agilent SureSelect) was used to capture target exons from genomic DNA of ENKTCL tumors. DNA was fragmented to 200bp using a Covaris S220 Focused-ultrasonicator and was end-repaired and ligated to Illumina adaptors with a KAPA Hyper Prep Kit. Ligated products were purified and amplified with a 7-cycle PCR. Exome capture from the PCR products was performed using a SureSelect Target Enrichment Kit. The captured exomes were further amplified with an 8-cycle PCR and purified for cluster generation using the cBot cluster generation system with a HiSeq PE Cluster Kit V4. 100bp paired-end reads were generated from an Illumina HiSeq 2500 system with a HiSeq SBS Kit V4.

  Study phs002925

• The Haplotype-Resolved Genome and Epigenome of the Aneuploid HeLa Cancer Cell Line

  The genomic architecture of HeLa remains largely unexplored beyond its karyotype, partly because like many cancers, its extensive aneuploidy renders such analyses challenging. We carried out haplotype-resolved whole-genome sequencing of the HeLa CCL-2 strain, examined point- and indel-mutation variations, mapped copy-number variations and loss of heterozygosity regions, and phased variants across full chromosome arms. We also investigated variation and copy-number profiles for HeLa S3 and eight additional strains. We find that HeLa is relatively stable in terms of point variation, with few new mutations accumulating after early passaging. Haplotype resolution facilitated reconstruction of an amplified, highly rearranged region of chromosome 8q24.21 at which integration of the human papilloma virus type 18 (HPV-18) genome occurred and that is likely to be the event that initiated tumorigenesis. We combined these maps with RNA-seq and ENCODE Project data sets to phase the HeLa epigenome. This revealed strong, haplotype specific activation of the proto-oncogene MYC by the integrated HPV-18 genome approximately 500 kilobases upstream, and enabled global analyses of the relationship between gene dosage and expression. These data provide an extensively phased, high-quality reference genome for past and future experiments relying on HeLa, and demonstrate the value of haplotype resolution for characterizing cancer genomes and epigenomes. Reprinted from Adey et. al. 2013, Nature, with permission from Nature Publishing Group.

  Study phs000642