Mapping the cellular biogeography of human bone marrow niches using single-cell transcriptomics and proteomic imaging.

Non-hematopoietic cells are essential contributors to hematopoiesis. However, heterogeneity and spatial organization of these cells in human bone marrow remain largely uncharacterized. We used single-cell RNA sequencing (scRNA-seq) to profile 29,325 non-hematopoietic cells and discovered nine transcriptionally distinct subtypes. We simultaneously profiled 53,417 hematopoietic cells and predicted their interactions with non-hematopoietic subsets. We employed co-detection by indexing (CODEX) to spatially profile over 1.2 million cells. We integrated scRNA-seq and CODEX data to link predicted cellular signaling with spatial proximity. Our analysis revealed a hyperoxygenated arterio-endosteal neighborhood for early myelopoiesis, and an adipocytic localization for early hematopoietic stem and progenitor cells (HSPCs). We used our CODEX atlas to annotate new images and uncovered mesenchymal stromal cell (MSC) expansion and spatial neighborhoods co-enriched for leukemic blasts and MSCs in acute myeloid leukemia (AML) patient samples. This spatially resolved, multiomic atlas of human bone marrow provides a reference for investigation of cellular interactions that drive hematopoiesis.

Efficient recall of Omicron-reactive B cell memory after a third dose of SARS-CoV-2 mRNA vaccine.

We examined antibody and memory B cell responses longitudinally for ∼9-10 months after primary 2-dose SARS-CoV-2 mRNA vaccination and 3 months after a 3rd dose. Antibody decay stabilized between 6 and 9 months, and antibody quality continued to improve for at least 9 months after 2-dose vaccination. Spike- and RBD-specific memory B cells remained durable over time, and 40%-50% of RBD-specific memory B cells simultaneously bound the Alpha, Beta, Delta, and Omicron variants. Omicron-binding memory B cells were efficiently reactivated by a 3rd dose of wild-type vaccine and correlated with the corresponding increase in neutralizing antibody titers. In contrast, pre-3rd dose antibody titers inversely correlated with the fold-change of antibody boosting, suggesting that high levels of circulating antibodies may limit the added protection afforded by repeat short interval boosting. These data provide insight into the quantity and quality of mRNA-vaccine-induced immunity over time through 3 or more antigen exposures.

Seasonal human coronavirus antibodies are boosted upon SARS-CoV-2 infection but not associated with protection.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has rapidly spread within the human population. Although SARS-CoV-2 is a novel coronavirus, most humans had been previously exposed to other antigenically distinct common seasonal human coronaviruses (hCoVs) before the coronavirus disease 2019 (COVID-19) pandemic. Here, we quantified levels of SARS-CoV-2-reactive antibodies and hCoV-reactive antibodies in serum samples collected from 431 humans before the COVID-19 pandemic. We then quantified pre-pandemic antibody levels in serum from a separate cohort of 251 individuals who became PCR-confirmed infected with SARS-CoV-2. Finally, we longitudinally measured hCoV and SARS-CoV-2 antibodies in the serum of hospitalized COVID-19 patients. Our studies indicate that most individuals possessed hCoV-reactive antibodies before the COVID-19 pandemic. We determined that ∼20% of these individuals possessed non-neutralizing antibodies that cross-reacted with SARS-CoV-2 spike and nucleocapsid proteins. These antibodies were not associated with protection against SARS-CoV-2 infections or hospitalizations, but they were boosted upon SARS-CoV-2 infection.

Human epigenetic and transcriptional T cell differentiation atlas for identifying functional T cell-specific enhancers.

The clinical benefit of T cell immunotherapies remains limited by incomplete understanding of T cell differentiation and dysfunction. We generated an epigenetic and transcriptional atlas of T cell differentiation from healthy humans that included exhausted CD8 T cells and applied this resource in three ways. First, we identified modules of gene expression and chromatin accessibility, revealing molecular coordination of differentiation after activation and between central memory and effector memory. Second, we applied this healthy molecular framework to three settings-a neoadjuvant anti-PD1 melanoma trial, a basal cell carcinoma scATAC-seq dataset, and autoimmune disease-associated SNPs-yielding insights into disease-specific biology. Third, we predicted genome-wide cis-regulatory elements and validated this approach for key effector genes using CRISPR interference, providing functional annotation and demonstrating the ability to identify targets for non-coding cellular engineering. These studies define epigenetic and transcriptional regulation of human T cells and illustrate the utility of interrogating disease in the context of a healthy T cell atlas.

Rapid induction of antigen-specific CD4+ T cells is associated with coordinated humoral and cellular immunity to SARS-CoV-2 mRNA vaccination.

SARS-CoV-2 mRNA vaccines have shown remarkable clinical efficacy, but questions remain about the nature and kinetics of T cell priming. We performed longitudinal antigen-specific T cell analyses on healthy SARS-CoV-2-naive and recovered individuals prior to and following mRNA prime and boost vaccination. Vaccination induced rapid antigen-specific CD4+ T cell responses in naive subjects after the first dose, whereas CD8+ T cell responses developed gradually and were variable in magnitude. Vaccine-induced Th1 and Tfh cell responses following the first dose correlated with post-boost CD8+ T cells and neutralizing antibodies, respectively. Integrated analysis revealed coordinated immune responses with distinct trajectories in SARS-CoV-2-naive and recovered individuals. Last, whereas booster vaccination improved T cell responses in SARS-CoV-2-naive subjects, the second dose had little effect in SARS-CoV-2-recovered individuals. These findings highlight the role of rapidly primed CD4+ T cells in coordinating responses to the second vaccine dose in SARS-CoV-2-naive individuals.

Cytomegalovirus Latent Infection is Associated with an Increased Risk of COVID-19-Related Hospitalization.

Some risk factors for severe coronavirus disease 2019 (COVID-19) have been identified, including age, race, and obesity. However, 20%-50% of severe cases occur in the absence of these factors. Cytomegalovirus (CMV) is a herpesvirus that infects about 50% of all individuals worldwide and is among the most significant nongenetic determinants of immune system. We hypothesized that latent CMV infection might influence the severity of COVID-19. Our analyses demonstrate that CMV seropositivity is associated with more than twice the risk of hospitalization due to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. Immune profiling of blood and CMV DNA quantitative polymerase chain reaction in a subset of patients for whom respiratory tract samples were available revealed altered T-cell activation profiles in absence of extensive CMV replication in the upper respiratory tract. These data suggest a potential role for CMV-driven immune perturbations in affecting the outcome of SARS-CoV-2 infection and may have implications for the discrepancies in COVID-19 severity between different human populations.

Genetic predisposition to neuroblastoma mediated by a LMO1 super-enhancer polymorphism.

Neuroblastoma is a paediatric malignancy that typically arises in early childhood, and is derived from the developing sympathetic nervous system. Clinical phenotypes range from localized tumours with excellent outcomes to widely metastatic disease in which long-term survival is approximately 40% despite intensive therapy. A previous genome-wide association study identified common polymorphisms at the LMO1 gene locus that are highly associated with neuroblastoma susceptibility and oncogenic addiction to LMO1 in the tumour cells. Here we investigate the causal DNA variant at this locus and the mechanism by which it leads to neuroblastoma tumorigenesis. We first imputed all possible genotypes across the LMO1 locus and then mapped highly associated single nucleotide polymorphism (SNPs) to areas of chromatin accessibility, evolutionary conservation and transcription factor binding sites. We show that SNP rs2168101 G>T is the most highly associated variant (combined P = 7.47 × 10(-29), odds ratio 0.65, 95% confidence interval 0.60-0.70), and resides in a super-enhancer defined by extensive acetylation of histone H3 lysine 27 within the first intron of LMO1. The ancestral G allele that is associated with tumour formation resides in a conserved GATA transcription factor binding motif. We show that the newly evolved protective TATA allele is associated with decreased total LMO1 expression (P = 0.028) in neuroblastoma primary tumours, and ablates GATA3 binding (P < 0.0001). We demonstrate allelic imbalance favouring the G-containing strand in tumours heterozygous for this SNP, as demonstrated both by RNA sequencing (P < 0.0001) and reporter assays (P = 0.002). These findings indicate that a recently evolved polymorphism within a super-enhancer element in the first intron of LMO1 influences neuroblastoma susceptibility through differential GATA transcription factor binding and direct modulation of LMO1 expression in cis, and this leads to an oncogenic dependency in tumour cells.

Common variants upstream of MLF1 at 3q25 and within CPZ at 4p16 associated with neuroblastoma.

Neuroblastoma is a cancer of the developing sympathetic nervous system that most commonly presents in young children and accounts for approximately 12% of pediatric oncology deaths. Here, we report on a genome-wide association study (GWAS) in a discovery cohort or 2,101 cases and 4,202 controls of European ancestry. We identify two new association signals at 3q25 and 4p16 that replicated robustly in multiple independent cohorts comprising 1,163 cases and 4,396 controls (3q25: rs6441201 combined P = 1.2x10-11, Odds Ratio 1.23, 95% CI:1.16-1.31; 4p16: rs3796727 combined P = 1.26x10-12, Odds Ratio 1.30, 95% CI: 1.21-1.40). The 4p16 signal maps within the carboxypeptidase Z (CPZ) gene. The 3q25 signal resides within the arginine/serine-rich coiled-coil 1 (RSRC1) gene and upstream of the myeloid leukemia factor 1 (MLF1) gene. Increased expression of MLF1 was observed in neuroblastoma cells homozygous for the rs6441201 risk allele (P = 0.02), and significant growth inhibition was observed upon depletion of MLF1 (P < 0.0001) in neuroblastoma cells. Taken together, we show that common DNA variants within CPZ at 4p16 and upstream of MLF1 at 3q25 influence neuroblastoma susceptibility and MLF1 likely plays an important role in neuroblastoma tumorigenesis.

QuantumClone: clonal assessment of functional mutations in cancer based on a genotype-aware method for clonal reconstruction.

Motivation: In cancer, clonal evolution is assessed based on information coming from single nucleotide variants and copy number alterations. Nonetheless, existing methods often fail to accurately combine information from both sources to truthfully reconstruct clonal populations in a given tumor sample or in a set of tumor samples coming from the same patient. Moreover, previously published methods detect clones from a single set of variants. As a result, compromises have to be done between stringent variant filtering [reducing dispersion in variant allele frequency estimates (VAFs)] and using all biologically relevant variants. Results: We present a framework for defining cancer clones using most reliable variants of high depth of coverage and assigning functional mutations to the detected clones. The key element of our framework is QuantumClone, a method for variant clustering into clones based on VAFs, genotypes of corresponding regions and information about tumor purity. We validated QuantumClone and our framework on simulated data. We then applied our framework to whole genome sequencing data for 19 neuroblastoma trios each including constitutional, diagnosis and relapse samples. We confirmed an enrichment of damaging variants within such pathways as MAPK (mitogen-activated protein kinases), neuritogenesis, epithelial-mesenchymal transition, cell survival and DNA repair. Most pathways had more damaging variants in the expanding clones compared to shrinking ones, which can be explained by the increased total number of variants between these two populations. Functional mutational rate varied for ancestral clones and clones shrinking or expanding upon treatment, suggesting changes in clone selection mechanisms at different time points of tumor evolution. Availability and implementation: Source code and binaries of the QuantumClone R package are freely available for download at https://CRAN.R-project.org/package=QuantumClone. Contact: gudrun.schleiermacher@curie.fr or valentina.boeva@inserm.fr. Supplementary information: Supplementary data are available at Bioinformatics online.

MIBG avidity correlates with clinical features, tumor biology, and outcomes in neuroblastoma: A report from the Children's Oncology Group.

BACKGROUND: Prior studies suggest that neuroblastomas that do not accumulate metaiodobenzylguanidine (MIBG) on diagnostic imaging (MIBG non-avid) may have more favorable features compared with MIBG avid tumors. We compared clinical features, biologic features, and clinical outcomes between patients with MIBG nonavid and MIBG avid neuroblastoma. PROCEDURE: Patients had metastatic high- or intermediate-risk neuroblastoma and were treated on Children's Oncology Group protocols A3973 or A3961. Comparisons of clinical and biologic features according to MIBG avidity were made with chi-squared or Fisher exact tests. Event-free (EFS) and overall (OS) survival compared using log-rank tests and modeled using Cox models. RESULTS: Thirty of 343 patients (8.7%) had MIBG nonavid disease. Patients with nonavid tumors were less likely to have adrenal primary tumors (34.5 vs. 57.2%; P = 0.019), bone metastases (36.7 vs. 61.7%; P = 0.008), or positive urine catecholamines (66.7 vs. 91.0%; P < 0.001) compared with patients with MIBG avid tumors. Nonavid tumors were more likely to be MYCN amplified (53.8 vs. 32.6%; P = 0.030) and had lower norepinephrine transporter expression. Patients with MIBG nonavid disease had a 5-year EFS of 50.0% compared with 38.7% for patients with MIBG avid disease (P = 0.028). On multivariate testing in high-risk patients, MIBG avidity was the sole adverse prognostic factor for EFS identified (hazard ratio 1.77; 95% confidence interval 1.04-2.99; P = 0.034). CONCLUSIONS: Patients with MIBG nonavid neuroblastoma have lower rates of adrenal primary tumors, bone metastasis, and catecholamine secretion. Despite being more likely to have MYCN-amplified tumors, these patients have superior outcomes compared with patients with MIBG avid disease.

CODEX: a normalization and copy number variation detection method for whole exome sequencing.

High-throughput sequencing of DNA coding regions has become a common way of assaying genomic variation in the study of human diseases. Copy number variation (CNV) is an important type of genomic variation, but detecting and characterizing CNV from exome sequencing is challenging due to the high level of biases and artifacts. We propose CODEX, a normalization and CNV calling procedure for whole exome sequencing data. The Poisson latent factor model in CODEX includes terms that specifically remove biases due to GC content, exon capture and amplification efficiency, and latent systemic artifacts. CODEX also includes a Poisson likelihood-based recursive segmentation procedure that explicitly models the count-based exome sequencing data. CODEX is compared to existing methods on a population analysis of HapMap samples from the 1000 Genomes Project, and shown to be more accurate on three microarray-based validation data sets. We further evaluate performance on 222 neuroblastoma samples with matched normals and focus on a well-studied rare somatic CNV within the ATRX gene. We show that the cross-sample normalization procedure of CODEX removes more noise than normalizing the tumor against the matched normal and that the segmentation procedure performs well in detecting CNVs with nested structures.

Identification of functionally active, low frequency copy number variants at 15q21.3 and 12q21.31 associated with prostate cancer risk.

Copy number variants (CNVs) are a recently recognized class of human germ line polymorphisms and are associated with a variety of human diseases, including cancer. Because of the strong genetic influence on prostate cancer, we sought to identify functionally active CNVs associated with susceptibility of this cancer type. We queried low-frequency biallelic CNVs from 1,903 men of Caucasian origin enrolled in the Tyrol Prostate Specific Antigen Screening Cohort and discovered two CNVs strongly associated with prostate cancer risk. The first risk locus (P = 7.7 × 10(-4), odds ratio = 2.78) maps to 15q21.3 and overlaps a noncoding enhancer element that contains multiple activator protein 1 (AP-1) transcription factor binding sites. Chromosome conformation capture (Hi-C) data suggested direct cis-interactions with distant genes. The second risk locus (P = 2.6 × 10(-3), odds ratio = 4.8) maps to the α-1,3-mannosyl-glycoprotein 4-ß-N-acetylglucosaminyltransferase C (MGAT4C) gene on 12q21.31. In vitro cell-line assays found this gene to significantly modulate cell proliferation and migration in both benign and cancer prostate cells. Furthermore, MGAT4C was significantly overexpressed in metastatic versus localized prostate cancer. These two risk associations were replicated in an independent PSA-screened cohort of 800 men (15q21.3, combined P = 0.006; 12q21.31, combined P = 0.026). These findings establish noncoding and coding germ line CNVs as significant risk factors for prostate cancer susceptibility and implicate their role in disease development and progression.

Mapping the Cellular Biogeography of Human Bone Marrow Niches Using Single-Cell Transcriptomics and Proteomic Imaging.

The bone marrow is the organ responsible for blood production. Diverse non-hematopoietic cells contribute essentially to hematopoiesis. However, these cells and their spatial organization remain largely uncharacterized as they have been technically challenging to study in humans. Here, we used fresh femoral head samples and performed single-cell RNA sequencing (scRNA-Seq) to profile 29,325 enriched non-hematopoietic bone marrow cells and discover nine transcriptionally distinct subtypes. We next employed CO-detection by inDEXing (CODEX) multiplexed imaging of 18 individuals, including both healthy and acute myeloid leukemia (AML) samples, to spatially profile over one million single cells with a novel 53-antibody panel. We discovered a relatively hyperoxygenated arterio-endosteal niche for early myelopoiesis, and an adipocytic, but not endosteal or perivascular, niche for early hematopoietic stem and progenitor cells. We used our atlas to predict cell type labels in new bone marrow images and used these predictions to uncover mesenchymal stromal cell (MSC) expansion and leukemic blast/MSC-enriched spatial neighborhoods in AML patient samples. Our work represents the first comprehensive, spatially-resolved multiomic atlas of human bone marrow and will serve as a reference for future investigation of cellular interactions that drive hematopoiesis.

A longitudinal single-cell and spatial multiomic atlas of pediatric high-grade glioma.

Pediatric high-grade glioma (pHGG) is an incurable central nervous system malignancy that is a leading cause of pediatric cancer death. While pHGG shares many similarities to adult glioma, it is increasingly recognized as a molecularly distinct, yet highly heterogeneous disease. In this study, we longitudinally profiled a molecularly diverse cohort of 16 pHGG patients before and after standard therapy through single-nucleus RNA and ATAC sequencing, whole-genome sequencing, and CODEX spatial proteomics to capture the evolution of the tumor microenvironment during progression following treatment. We found that the canonical neoplastic cell phenotypes of adult glioblastoma are insufficient to capture the range of tumor cell states in a pediatric cohort and observed differential tumor-myeloid interactions between malignant cell states. We identified key transcriptional regulators of pHGG cell states and did not observe the marked proneural to mesenchymal shift characteristic of adult glioblastoma. We showed that essential neuromodulators and the interferon response are upregulated post-therapy along with an increase in non-neoplastic oligodendrocytes. Through in vitro pharmacological perturbation, we demonstrated novel malignant cell-intrinsic targets. This multiomic atlas of longitudinal pHGG captures the key features of therapy response that support distinction from its adult counterpart and suggests therapeutic strategies which are targeted to pediatric gliomas.

Genetic Predisposition to Neuroblastoma Results from a Regulatory Polymorphism that Promotes the Adrenergic Cell State.

Childhood neuroblastomas exhibit plasticity between an undifferentiated neural crest-like "mesenchymal" cell state and a more differentiated sympathetic "adrenergic" cell state. These cell states are governed by autoregulatory transcriptional loops called core regulatory circuitries (CRCs), which drive the early development of sympathetic neuronal progenitors from migratory neural crest cells during embryogenesis. The adrenergic cell identity of neuroblastoma requires LMO1 as a transcriptional co-factor. Both LMO1 expression levels and the risk of developing neuroblastoma in children are associated with a single nucleotide polymorphism G/T that affects a G ATA motif in the first intron of LMO1. Here we show that wild-type zebrafish with the G ATA genotype develop adrenergic neuroblastoma, while knock-in of the protective T ATA allele at this locus reduces the penetrance of MYCN-driven tumors, which are restricted to the mesenchymal cell state. Whole genome sequencing of childhood neuroblastomas demonstrates that T ATA/ T ATA tumors also exhibit a mesenchymal cell state and are low risk at diagnosis. Thus, conversion of the regulatory G ATA to a T ATA allele in the first intron of LMO1 reduces the neuroblastoma initiation rate by preventing formation of the adrenergic cell state, a mechanism that is conserved over 400 million years of evolution separating zebrafish and humans.

Genetic predisposition to neuroblastoma results from a regulatory polymorphism that promotes the adrenergic cell state.

Childhood neuroblastomas exhibit plasticity between an undifferentiated neural crest-like mesenchymal cell state and a more differentiated sympathetic adrenergic cell state. These cell states are governed by autoregulatory transcriptional loops called core regulatory circuitries (CRCs), which drive the early development of sympathetic neuronal progenitors from migratory neural crest cells during embryogenesis. The adrenergic cell identity of neuroblastoma requires LMO1 as a transcriptional cofactor. Both LMO1 expression levels and the risk of developing neuroblastoma in children are associated with a single nucleotide polymorphism, G/T, that affects a GATA motif in the first intron of LMO1. Here, we showed that WT zebrafish with the GATA genotype developed adrenergic neuroblastoma, while knock-in of the protective TATA allele at this locus reduced the penetrance of MYCN-driven tumors, which were restricted to the mesenchymal cell state. Whole genome sequencing of childhood neuroblastomas demonstrated that TATA/TATA tumors also exhibited a mesenchymal cell state and were low risk at diagnosis. Thus, conversion of the regulatory GATA to a TATA allele in the first intron of LMO1 reduced the neuroblastoma-initiation rate by preventing formation of the adrenergic cell state. This mechanism was conserved over 400 million years of evolution, separating zebrafish and humans.

Human T follicular helper clones seed the germinal center-resident regulatory pool.

The mechanisms by which FOXP3+ T follicular regulatory (Tfr) cells simultaneously steer antibody formation toward microbe or vaccine recognition and away from self-reactivity remain incompletely understood. To explore underappreciated heterogeneity in human Tfr cell development, function, and localization, we used paired TCRVA/TCRVB sequencing to distinguish tonsillar Tfr cells that are clonally related to natural regulatory T cells (nTfr) from those likely induced from T follicular helper (Tfh) cells (iTfr). The proteins iTfr and nTfr cells differentially expressed were used to pinpoint their in situ locations via multiplex microscopy and establish their divergent functional roles. In silico analyses and in vitro tonsil organoid tracking models corroborated the existence of separate Treg-to-nTfr and Tfh-to-iTfr developmental trajectories. Our results identify human iTfr cells as a distinct CD38+, germinal center-resident, Tfh-descended subset that gains suppressive function while retaining the capacity to help B cells, whereas CD38- nTfr cells are elite suppressors primarily localized in follicular mantles. Interventions differentially targeting specific Tfr cell subsets may provide therapeutic opportunities to boost immunity or more precisely treat autoimmune diseases.

IgA deficiency destabilizes homeostasis toward intestinal microbes and increases systemic immune dysregulation.

The ability of most patients with selective immunoglobulin A (IgA) deficiency (SIgAD) to remain apparently healthy has been a persistent clinical conundrum. Compensatory mechanisms, including IgM, have been proposed, yet it remains unclear how secretory IgA and IgM work together in the mucosal system and, on a larger scale, whether the systemic and mucosal anti-commensal responses are redundant or have unique features. To address this gap in knowledge, we developed an integrated host-commensal approach combining microbial flow cytometry and metagenomic sequencing (mFLOW-Seq) to comprehensively define which microbes induce mucosal and systemic antibodies. We coupled this approach with high-dimensional immune profiling to study a cohort of pediatric patients with SIgAD and household control siblings. We found that mucosal and systemic antibody networks cooperate to maintain homeostasis by targeting a common subset of commensal microbes. In IgA-deficiency, we find increased translocation of specific bacterial taxa associated with elevated levels of systemic IgG targeting fecal microbiota. Associated features of immune system dysregulation in IgA-deficient mice and humans included elevated levels of inflammatory cytokines, enhanced follicular CD4 T helper cell frequency and activation, and an altered CD8 T cell activation state. Although SIgAD is clinically defined by the absence of serum IgA, the symptomatology and immune dysregulation were concentrated in the SIgAD participants who were also fecal IgA deficient. These findings reveal that mucosal IgA deficiency leads to aberrant systemic exposures and immune responses to commensal microbes, which increase the likelihood of humoral and cellular immune dysregulation and symptomatic disease in patients with IgA deficiency.

Transcriptional immunogenomic analysis reveals distinct immunological clusters in paediatric nervous system tumours.

BACKGROUND: Cancer immunotherapies including immune checkpoint inhibitors and Chimeric Antigen Receptor (CAR) T-cell therapy have shown variable response rates in paediatric patients highlighting the need to establish robust biomarkers for patient selection. While the tumour microenvironment in adults has been widely studied to delineate determinants of immune response, the immune composition of paediatric solid tumours remains relatively uncharacterized calling for investigations to identify potential immune biomarkers. METHODS: To inform immunotherapy approaches in paediatric cancers with embryonal origin, we performed an immunogenomic analysis of RNA-seq data from 925 treatment-naïve paediatric nervous system tumours (pedNST) spanning 12 cancer types from three publicly available data sets. RESULTS: Within pedNST, we uncovered four broad immune clusters: Paediatric Inflamed (10%), Myeloid Predominant (30%), Immune Neutral (43%) and Immune Desert (17%). We validated these clusters using immunohistochemistry, methylation immune inference and segmentation analysis of tissue images. We report shared biology of these immune clusters within and across cancer types, and characterization of specific immune cell frequencies as well as T- and B-cell repertoires. We found no associations between immune infiltration levels and tumour mutational burden, although molecular cancer entities were enriched within specific immune clusters. CONCLUSIONS: Given the heterogeneity of immune infiltration within pedNST, our findings suggest personalized immunogenomic profiling is needed to guide selection of immunotherapeutic strategies.

Optimizing copy number variation analysis using genome-wide short sequence oligonucleotide arrays.

The detection of copy number variants (CNV) by array-based platforms provides valuable insight into understanding human diversity. However, suboptimal study design and data processing negatively affect CNV assessment. We quantitatively evaluate their impact when short-sequence oligonucleotide arrays are applied (Affymetrix Genome-Wide Human SNP Array 6.0) by evaluating 42 HapMap samples for CNV detection. Several processing and segmentation strategies are implemented, and results are compared to CNV assessment obtained using an oligonucleotide array CGH platform designed to query CNVs at high resolution (Agilent). We quantitatively demonstrate that different reference models (e.g. single versus pooled sample reference) used to detect CNVs are a major source of inter-platform discrepancy (up to 30%) and that CNVs residing within segmental duplication regions (higher reference copy number) are significantly harder to detect (P < 0.0001). After adjusting Affymetrix data to mimic the Agilent experimental design (reference sample effect), we applied several common segmentation approaches and evaluated differential sensitivity and specificity for CNV detection, ranging 39-77% and 86-100% for non-segmental duplication regions, respectively, and 18-55% and 39-77% for segmental duplications. Our results are relevant to any array-based CNV study and provide guidelines to optimize performance based on study-specific objectives.