Molecular, clinical and therapeutic determinants of outcome in NPM1 mutated AML

Although NPM1-mutated acute myeloid leukemia (AML) carries a generally favorable prognosis, many patients still relapse and die. Previous studies identified several molecular and clinical features associated with poor outcome, however only FLT3-ITD mutation and adverse karyotype are currently used for risk stratification due to inconsistent results and uncertainty around how other factors should influence treatment, particularly given the strong prognostic impact of post-induction measurable residual disease (MRD). Here we analyzed a large group of patients with NPM1mut AML enrolled in prospective trials (NCRI AML17 and AML19, n=1357) to delineate the impact of baseline molecular and clinical features, post induction MRD status and treatment intensity on outcome. FLT3-ITD (HR 1.28, 95%CI 1.01-1.63), DNMT3A (HR 1.65, 95%CI 1.32-2.05), WT1 (HR 1.74, 95%CI 1272-2.38) and non-ABD NPM1 mutations (HR 1.64, 95%CI 1.22-2.21) were independently associated with poorer overall survival (OS). These factors were also strongly associated with MRD positivity. For patients achieving MRD negativity, these mutations (except FLT3-ITD) were associated with an increased cumulative incidence of relapse (CIR) and poorer OS. However, apart from the few patients with adverse cytogenetics, we could not identify any group of MRD negative patients with a CIR >40% or with benefit from allograft in first remission. Intensified chemotherapy with the FLAG-Ida regimen was associated with improved outcomes in all subgroups, with greater benefits observed in the highest risk molecular subgroups.

Donor whole blood DNA methylation is not a strong predictor of acute graft versus host disease in unrelated donor allogeneic haematopoietic cell transplantation

Allogeneic hematopoietic cell transplantation (HCT) is used to treat many blood-based disorders and malignancies, however it can also result in serious adverse events, such as the development of acute graft-versus-host disease (aGVHD). This study aimed to develop a donor-specific epigenetic classifier to reduce incidence of aGVHD by improving donor selection. Genome-wide DNA methylation was assessed in a discovery cohort of 288 HCT donors selected based on recipient aGVHD outcome; this cohort consisted of 144 cases with aGVHD grades III-IV and 144 controls with no aGVHD. We applied a machine learning algorithm to identify CpG sites predictive of aGVHD. Receiver operating characteristic (ROC) curve analysis of these sites resulted in a classifier with an encouraging area under the ROC curve (AUC) of 0.91. To test this classifier, we used an independent validation cohort (n = 288) selected using the same criteria as the discovery cohort. Attempts to validate the classifier failed with the AUC falling to 0.51. These results indicate that donor DNA methylation may not be a suitable predictor of aGVHD in an HCT setting involving unrelated donors, despite the initial promising results in the discovery cohort. Our work highlights the importance of independent validation of machine learning classifiers, particularly when developing classifiers intended for clinical use.

Single-cell analysis of psoriasis resolution demonstrates an inflammatory fibroblast state targeted by IL-23 blockade

Biologic therapies targeting the IL-23/IL-17 axis have transformed the treatment of psoriasis. However, the early mechanisms of action of these drugs remain poorly understood. Here, we perform longitudinal single-cell RNA-sequencing in affected individuals receiving IL-23 inhibitor therapy. By profiling skin at baseline, day 3 and day 14 of treatment, we demonstrate that IL-23 blockade causes marked gene expression shifts, with fibroblast and myeloid populations displaying the most extensive changes at day 3. We also identify a transient WNT5A+/IL24+ fibroblast state, which is only detectable in lesional skin. In-silico and in-vitro studies indicate that signals stemming from these WNT5A+/IL24+ fibroblasts upregulate multiple inflammatory genes in keratinocytes. Importantly, the abundance of WNT5A+/IL24+ fibroblasts is significantly reduced after treatment. This observation is validated in-silico, by deconvolution of multiple transcriptomic datasets, and experimentally, by RNA in-situ hybridization. These findings demonstrate that the evolution of inflammatory fibroblast states is a key feature of resolving psoriasis skin.

Repurposing Metformin for periodontal disease management as a form of oral-systemic preventive medicine

BACKGROUND

Despite the improvements in treatment over the last decades, periodontal disease (PD) affects millions of people around the world and the only treatment available is based on controlling microbial load. Diabetes is known to increase the risk of PD establishment and progression, and recently, glucose metabolism modulation by pharmaceutical or dietarian means has been emphasised as a significant modulator of non-communicable disease development.

METHODS

The impact of pharmaceutically controlling glucose metabolism in non-diabetic animals and humans (REBEC, UTN code: U1111-1276-1942) was investigated by repurposing Metformin, as a mean to manage periodontal disease and its associated systemic risk factors.

RESULTS

We found that glucose metabolism control via use of Metformin aimed at PD management resulted in significant prevention of bone loss during induced periodontal disease and age-related bone loss in vivo. Metformin also influenced the bacterial species present in the oral environment and impacted the metabolic epithelial and stromal responses to bacterial dysbiosis at a single cell level. Systemically, Metformin controlled blood glucose levels and age-related weight gain when used long-term. Translationally, our pilot randomized control trial indicated that systemic Metformin was safe to use in non-diabetic patients and affected the periodontal tissues. During the medication window, patients showed stable levels of systemic blood glucose, lower circulating hsCRP and lower insulin levels after periodontal treatment when compared to placebo. Finally, patients treated with Metformin had improved periodontal parameters when compared to placebo treated patients.

CONCLUSION

This is the first study to demonstrate that systemic interventions using Metformin in non-diabetic individuals aimed at PD prevention have oral-systemic effects constituting a possible novel form of preventive medicine for oral-systemic disease management.

B cell profiles, antibody repertoire and reactivity reveal dysregulated responses with autoimmune features in melanoma

B cells are known to contribute to the anti-tumor immune response, especially in immunogenic tumors such as melanoma, yet humoral immunity has not been characterized in these cancers to detail. Here we show comprehensive phenotyping in samples of circulating and tumor-resident B cells as well as serum antibodies in melanoma patients. Memory B cells are enriched in tumors compared to blood in paired samples and feature distinct antibody repertoires, linked to specific isotypes. Tumor-associated B cells undergo clonal expansion, class switch recombination, somatic hypermutation and receptor revision. Compared with blood, tumor-associated B cells produce antibodies with proportionally higher levels of unproductive sequences and distinct complementarity determining region 3 properties. The observed features are signs of affinity maturation and polyreactivity and suggest an active and aberrant autoimmune-like reaction in the tumor microenvironment. Consistent with this, tumor-derived antibodies are polyreactive and characterized by autoantigen recognition. Serum antibodies show reactivity to antigens attributed to autoimmune diseases and cancer, and their levels are higher in patients with active disease compared to post-resection state. Our findings thus reveal B cell lineage dysregulation with distinct antibody repertoire and specificity, alongside clonally-expanded tumor-infiltrating B cells with autoimmune-like features, shaping the humoral immune response in melanoma.

Psoriatic and rheumatoid arthritis joints differ in the composition of CD8+ tissue-resident memory T cell subsets

CD69+CD103+ tissue-resident memory T (TRM) cells are important drivers of inflammation. To decipher their role in inflammatory arthritis, we apply single-cell, high-dimensional profiling to T cells from the joints of patients with psoriatic arthritis (PsA) or rheumatoid arthritis (RA). We identify three groups of synovial CD8+CD69+CD103+ TRM cells: cytotoxic and regulatory T (Treg)-like TRM cells are present in both PsA and RA, while CD161+CCR6+ type 17-like TRM cells with a pro-inflammatory cytokine profile (IL-17A+TNFα+IFNγ+) are specifically enriched in PsA. In contrast, only one population of CD4+CD69+CD103+ TRM cells is detected and at similarly low frequencies in both diseases. Type 17-like CD8+ TRM cells have a distinct transcriptomic signature and a polyclonal, but distinct, TCR repertoire. Type 17-like cells are also enriched in CD8+CD103- T cells in PsA compared with RA. These findings illustrate differences in the immunopathology of PsA and RA, with a particular enrichment for type 17 CD8+ T cells in the PsA joint.

Abstract 6012: Comparison of pre-analytical methods for DNA methylation analysis of cfDNA

Liquid biopsies in the form of cell-free circulating tumor DNA (ctDNA) isolated from blood plasma or serum have shown great promise to be used as minimally-invasive markers in cancer patients providing an insight into the composition of both primary tumor and metastases, and an opportunity to assess clonal dynamics throughout the course of disease. CtDNA carries not only tumor- specific changes in its sequence but also distinctive epigenetic marks, namely DNA methylation patterns similar to the patterns in the primary tumor. Deregulation of tissue-specific DNA methylation patterns is an early event in carcinogenesis holding great promise to be exploited for early detection of cancer. However, in contrast to mutation detection, targeted DNA methylation sequencing of cfDNA remains a challenge due to low concentration and heavy fragmentation. The limiting step in detection of any analyte is the recovery of the DNA molecules. The pre-analytical steps in the DNA methylation detection pipeline include extraction of cfDNA molecules and their conversion using sodium bisulfite. Bisulfite conversion is a harsh chemical reaction that induces further fragmentation and degradation of already scarce cfDNA. In this study we compared commercial methods for cfDNA extraction and for bisulfite conversion to evaluate the yield, degradation and procedural loss. We compared five commercial kits for cfDNA extraction from blood plasma, three column-based (Qiagen’s QIAamp Circulating nucleic acids kit, Norgen’s Plasma/Serum Cell-Free Circulating DNA Purification and Macherey-Nagel’s NucleoSpin® Plasma XS) and two based on magnetic beads cleanup (Bioo Scientific’s Next Prep magnetic Cell-Free Circulating DNA isolation kitand Applied Biosystems’ Magmax Cell-Free Circulating DNA isolation kit. The kits used for bisulfite conversion were Zymo EZ DNA Methylation Lightning with column based purification and with magnetic beads purification, Zymo EZ DNA Methylation Direct with column based purification and with magnetic bead purification, Qiagen Epitect Fast DNA Bisulfite and Diagenode Premium Bisulfite kits. Same samples were processed by each of the kits, Sample degradation was determined using Agilent Bioanalyzer dsDNA kit. Yield was evaluated by Picogreen, absolute qPCR for single-copy gene and for contamination of gDNA by a B-cell specific assay, while procedural loss was determined based on the recovery of a synthetic spike-in.

Citation Format: Miljana Tanic, Oluwadunni Emiloju, Pawan Dhami, Stephan Beck. Comparison of pre-analytical methods for DNA methylation analysis of cfDNA. [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 6012.

Loss of the crumbs cell polarity complex disrupts epigenetic transcriptional control and cell cycle progression in the developing retina

The crumbs cell polarity complex plays a crucial role in apical-basal epithelial polarity, cellular adhesion, and morphogenesis. Homozygous variants in human CRB1 result in autosomal recessive Leber congenital amaurosis (LCA) and retinitis pigmentosa (RP), with no established genotype-phenotype correlation. The associated protein complexes have key functions in developmental pathways; however, the underlying disease mechanism remains unclear. Using the oko meduzym289/m289 (crb2a-/- ) zebrafish, we performed integrative transcriptomic (RNA-seq data) and methylomic [reduced representation bisulphite sequencing (RRBS)] analysis of whole retina to identify dysregulated genes and pathways. Delayed retinal cell specification was identified in both the crb2a-/- zebrafish and CRB1 patient-derived retinal organoids, highlighting the dysfunction of cell cycle modulation and epigenetic transcriptional control. Differential DNA methylation analysis revealed novel hypermethylated pathways involving biological adhesion, Hippo, and transforming growth factor β (TGFβ) signalling. By integrating gene expression with DNA methylation using functional epigenetic modules (FEM), we identified six key modules involving cell cycle control and disturbance of TGFβ, bone morphogenetic protein (BMP), Hippo, and SMAD protein signal transduction pathways, revealing significant interactome hotspots relevant to crb2a function and confirming the epigenetic control of gene regulation in early retinal development, which points to a novel mechanism underlying CRB1-retinopathies. © 2023 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Spatially resolved transcriptomics reveals pro-inflammatory fibroblast involved in lymphocyte recruitment through CXCL8 and CXCL10

The interplay among different cells in a tissue is essential for maintaining homeostasis. Although disease states have been traditionally attributed to individual cell types, increasing evidence and new therapeutic options have demonstrated the primary role of multicellular functions to understand health and disease, opening new avenues to understand pathogenesis and develop new treatment strategies. We recently described the cellular composition and dynamics of the human oral mucosa; however, the spatial arrangement of cells is needed to better understand a morphologically complex tissue. Here, we link single-cell RNA sequencing, spatial transcriptomics, and high-resolution multiplex fluorescence in situ hybridisation to characterise human oral mucosa in health and oral chronic inflammatory disease. We deconvolved expression for resolution enhancement of spatial transcriptomic data and defined highly specialised epithelial and stromal compartments describing location-specific immune programs. Furthermore, we spatially mapped a rare pathogenic fibroblast population localised in a highly immunogenic region, responsible for lymphocyte recruitment through CXCL8 and CXCL10 and with a possible role in pathological angiogenesis through ALOX5AP. Collectively, our study provides a comprehensive reference for the study of oral chronic disease pathogenesis.

Single-cell transcriptomic and spatial landscapes of the developing human pancreas

Current differentiation protocols have not been successful in reproducibly generating fully functional human beta cells in vitro, partly due to incomplete understanding of human pancreas development. Here, we present detailed transcriptomic analysis of the various cell types of the developing human pancreas, including their spatial gene patterns. We integrated single-cell RNA sequencing with spatial transcriptomics at multiple developmental time points and revealed distinct temporal-spatial gene cascades. Cell trajectory inference identified endocrine progenitor populations and branch-specific genes as the progenitors differentiate toward alpha or beta cells. Spatial differentiation trajectories indicated that Schwann cells are spatially co-located with endocrine progenitors, and cell-cell connectivity analysis predicted that they may interact via L1CAM-EPHB2 signaling. Our integrated approach enabled us to identify heterogeneity and multiple lineage dynamics within the mesenchyme, showing that it contributed to the exocrine acinar cell state. Finally, we have generated an interactive web resource for investigating human pancreas development for the research community.

Abstract 5710: Identification of convergent gene repression mechanisms through integrative genomic and DNA methylation analysis in TRACERx

Introduction: Lung TRACERx is a prominent study employing multi-region and longitudinal multi-omics sequencing to unravel the evolutionary trajectories of lung cancer. Aberrant DNA methylation patterns have been widely described in nearly all human cancers, yet their interplay with DNA mutations in lung cancer is not well understood. Incorporating the contribution of epigenetic modifications to cancer evolution trajectories within TRACERx could improve our understanding of the intricate relationship between genetic and epigenetic changes in non-small cell lung cancer (NSCLC) evolution.

Methods: Multi-region sampling from 38 TRACERx patients including 112 tumor regions and 37 matched normal adjacent tissue samples was performed. Reduced representation bisulfite sequencing (RRBS) was performed to assess DNA methylation and the CAMDAC (Larose-Cadieux et al, 2020) was applied to estimate purified tumor methylation rates and correct for copy number changes. Whole exome sequencing and somatic copy number alterations (SCNAs) were inferred using the ASCAT tool (Van Loo et al, 2010) and Methsig (Pan et al, 2021) was performed to discover new methylation driver genes.

Results: Using multi-region sequencing, we identified ubiquitous hypermethylation of 29 known driver genes in both lung adenocarcinoma (LUAD) and squamous cell lung cancer (LUSC), together with an additional 9 and 27 genes exclusive to LUSC and LUAD, respectively. We also identified 13 and 7 driver genes non-ubiquitously hypermethylated exclusively in LUSC and LUAD, respectively. Using a differential methylation based approach, we describe a method to determine the extent of intra-tumor methylation heterogeneity akin to established ITH scores based on genomics data. In addition, we report the identification of novel subtype-specific methylation driver genes enriched in HOX family members which are related to cancer progression. Through integration of DNA methylation and genomic sequencing data, we identify parallel mechanisms contributing towards ubiquitous tumor suppressor gene alterations. At the patient level, multiple driver genes such as NSD1, GATA3 and MGA were subject to repression by both copy number loss and DNA hypermethylation. Finally, we describe dosage-compensation of genes such as the Notch ligands JAG2 and DLK1 that are proximal to amplified oncogenes and hypermethylated during tumor evolution.

Conclusion: We describe the contribution of DNA methylation and genomic alterations to altering the landscape of NSCLC. Leveraging DNA methylation, we can determine the extent of convergent repression mechanisms in different regions of the same tumor, assess DNA methylation heterogeneity, and discover DNA methylation-based driver genes in NSCLC.

Citation Format: Francisco Gimeno-Valiente, Carla Castignani, Elizabeth Larose-Cadieux, Kezhong Chen, Nana Mensah, Olga Chervova, Thomas Watkins, Pawan Dhami, Heli Vaikkinen, Andrew Feber, TRACERx Consortium, Jonas Demeulemeester, Miljana Tanic, Stephan Beck, Peter Van Loo, Charles Swanton, Nnennaya Kanu. Identification of convergent gene repression mechanisms through integrative genomic and DNA methylation analysis in TRACERx [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 5710.

Systematic Evaluation of the Immune Environment of Small Intestinal Neuroendocrine Tumors

PURPOSE

The immune tumor microenvironment and the potential therapeutic opportunities for immunotherapy in small intestinal neuroendocrine tumors (siNET) have not been fully defined.

EXPERIMENTAL DESIGN

Herein, we studied 40 patients with primary and synchronous metastatic siNETs, and matched blood and normal tissue obtained during surgery. We interrogated the immune checkpoint landscape using multi-parametric flow cytometry. In addition, matched FFPE tissue was obtained for multi-parametric IHC to determine the relative abundance and distribution of T-cell infiltrate. Tumor mutational burden (TMB) was also assessed and correlated with immune infiltration.

RESULTS

Effector tumor-infiltrating lymphocytes (TIL) had a higher expression of PD-1 in the tumor microenvironment compared with the periphery. In addition, CD8+ TILs had a significantly higher co-expression of PD-1/ICOS and PD-1/CTLA-4 (cytotoxic T lymphocyte antigen-4) and higher levels of PD-1 expression compared with normal tissue. IHC revealed that the majority of cases have ≤10% intra-tumoral T cells but a higher number of peri-tumoral T cells, demonstrating an "exclusion" phenotype. Finally, we confirmed that siNETs have a low TMB compared with other tumor types in the TCGA database but did not find a correlation between TMB and CD8/Treg ratio.

CONCLUSIONS

Taken together, these results suggest that a combination therapy approach will be required to enhance the immune response, using PD-1 as a checkpoint immunomodulator backbone in combination with other checkpoint targeting molecules (CTLA-4 or ICOS), or with drugs targeting other pathways to recruit "excluded" T cells into the tumor microenvironment to treat patients with siNETs.

Comparison and imputation-aided integration of five commercial platforms for targeted DNA methylome analysis

Targeted bisulfite sequencing (TBS) has become the method of choice for the cost-effective, targeted analysis of the human methylome at base-pair resolution. In this study, we benchmarked five commercially available TBS platforms-three hybridization capture-based (Agilent, Roche and Illumina) and two reduced-representation-based (Diagenode and NuGen)-across 11 samples. Two samples were also compared with whole-genome DNA methylation sequencing with the Illumina and Oxford Nanopore platforms. We assessed workflow complexity, on/off-target performance, coverage, accuracy and reproducibility. Although all platforms produced robust and reproducible data, major differences in the number and identity of the CpG sites covered make it difficult to compare datasets generated on different platforms. To overcome this limitation, we applied imputation and show that it improves interoperability from an average of 10.35% (0.8 million) to 97% (7.6 million) common CpG sites. Our study provides guidance on which TBS platform to use for different methylome features and offers an imputation-based harmonization solution that allows comparative, integrative analysis.

Signatures of copy number alterations in human cancer

Gains and losses of DNA are prevalent in cancer and emerge as a consequence of inter-related processes of replication stress, mitotic errors, spindle multipolarity and breakage-fusion-bridge cycles, among others, which may lead to chromosomal instability and aneuploidy1,2. These copy number alterations contribute to cancer initiation, progression and therapeutic resistance3-5. Here we present a conceptual framework to examine the patterns of copy number alterations in human cancer that is widely applicable to diverse data types, including whole-genome sequencing, whole-exome sequencing, reduced representation bisulfite sequencing, single-cell DNA sequencing and SNP6 microarray data. Deploying this framework to 9,873 cancers representing 33 human cancer types from The Cancer Genome Atlas6 revealed a set of 21 copy number signatures that explain the copy number patterns of 97% of samples. Seventeen copy number signatures were attributed to biological phenomena of whole-genome doubling, aneuploidy, loss of heterozygosity, homologous recombination deficiency, chromothripsis and haploidization. The aetiologies of four copy number signatures remain unexplained. Some cancer types harbour amplicon signatures associated with extrachromosomal DNA, disease-specific survival and proto-oncogene gains such as MDM2. In contrast to base-scale mutational signatures, no copy number signature was associated with many known exogenous cancer risk factors. Our results synthesize the global landscape of copy number alterations in human cancer by revealing a diversity of mutational processes that give rise to these alterations.

Two subsets of human marginal zone B cells resolved by global analysis of lymphoid tissues and blood

B cells generate antibodies that are essential for immune protection, but their subgroups are poorly defined. Here, we perform undirected deep profiling of B cells in matched human lymphoid tissues from deceased transplant organ donors and blood. In addition to identifying unanticipated features of tissue-based B cell differentiation, we resolve two subsets of marginal zone B (MZB) cells differing in cell surface and transcriptomic profiles, clonal relationships to other subsets, enrichment of genes in the NOTCH pathway, distribution bias within splenic marginal zone microenvironment, and immunoglobulin repertoire diversity and hypermutation frequency. Each subset is present in spleen, gut-associated lymphoid tissue, mesenteric lymph nodes, and blood. MZB cells and the lineage from which they are derived are depleted in lupus nephritis. Here, we show that this depletion is of only one MZB subset. The other remains unchanged as a proportion of total B cells compared with health. Thus, it is important to factor MZB cell heterogeneity into studies of human B cell responses and pathology.

Whole-genome sequencing of single circulating tumor cells from neuroendocrine neoplasms

Single-cell profiling of circulating tumor cells (CTCs) as part of a minimally invasive liquid biopsy presents an opportunity to characterize and monitor tumor heterogeneity and evolution in individual patients. In this study, we aimed to compare single-cell copy number variation (CNV) data with tissue and define the degree of intra- and inter-patient genomic heterogeneity. We performed next-generation sequencing (NGS) whole-genome CNV analysis of 125 single CTCs derived from seven patients with neuroendocrine neoplasms (NEN) alongside matched white blood cells (WBC), formalin-fixed paraffin-embedded (FFPE), and fresh frozen (FF) samples. CTC CNV profiling demonstrated recurrent chromosomal alterations in previously reported NEN copy number hotspots, including the prognostically relevant loss of chromosome 18. Unsupervised hierarchical clustering revealed CTCs with distinct clonal lineages as well as significant intra- and inter-patient genomic heterogeneity, including subclonal alterations not detectable by bulk analysis and previously unreported in NEN. Notably, we also demonstrated the presence of genomically distinct CTCs according to the enrichment strategy utilized (EpCAM-dependent vs size-based). This work has significant implications for the identification of therapeutic targets, tracking of evolutionary change, and the implementation of CTC-biomarkers in cancer.

Human marginal zone B cell development from early T2 progenitors

B cells emerge from the bone marrow as transitional (TS) B cells that differentiate through T1, T2, and T3 stages to become naive B cells. We have identified a bifurcation of human B cell maturation from the T1 stage forming IgMhi and IgMlo developmental trajectories. IgMhi T2 cells have higher expression of α4β7 integrin and lower expression of IL-4 receptor (IL4R) compared with the IgMlo branch and are selectively recruited into gut-associated lymphoid tissue. IgMhi T2 cells also share transcriptomic features with marginal zone B cells (MZBs). Lineage progression from T1 cells to MZBs via an IgMhi trajectory is identified by pseudotime analysis of scRNA-sequencing data. Reduced frequency of IgMhi gut-homing T2 cells is observed in severe SLE and is associated with reduction of MZBs and their putative IgMhi precursors. The collapse of the gut-associated MZB maturational axis in severe SLE affirms its existence in health.

Undifferentiated Sarcomas Develop through Distinct Evolutionary Pathways

Undifferentiated sarcomas (USARCs) of adults are diverse, rare, and aggressive soft tissue cancers. Recent sequencing efforts have confirmed that USARCs exhibit one of the highest burdens of structural aberrations across human cancer. Here, we sought to unravel the molecular basis of the structural complexity in USARCs by integrating DNA sequencing, ploidy analysis, gene expression, and methylation profiling. We identified whole genome duplication as a prevalent and pernicious force in USARC tumorigenesis. Using mathematical deconvolution strategies to unravel the complex copy-number profiles and mutational timing models we infer distinct evolutionary pathways of these rare cancers. In addition, 15% of tumors exhibited raised mutational burdens that correlated with gene expression signatures of immune infiltration, and good prognosis.

Selfish mutations dysregulating RAS-MAPK signaling are pervasive in aged human testes

Mosaic mutations present in the germline have important implications for reproductive risk and disease transmission. We previously demonstrated a phenomenon occurring in the male germline, whereby specific mutations arising spontaneously in stem cells (spermatogonia) lead to clonal expansion, resulting in elevated mutation levels in sperm over time. This process, termed "selfish spermatogonial selection," explains the high spontaneous birth prevalence and strong paternal age-effect of disorders such as achondroplasia and Apert, Noonan and Costello syndromes, with direct experimental evidence currently available for specific positions of six genes (FGFR2, FGFR3, RET, PTPN11, HRAS, and KRAS). We present a discovery screen to identify novel mutations and genes showing evidence of positive selection in the male germline, by performing massively parallel simplex PCR using RainDance technology to interrogate mutational hotspots in 67 genes (51.5 kb in total) in 276 biopsies of testes from five men (median age, 83 yr). Following ultradeep sequencing (about 16,000×), development of a low-frequency variant prioritization strategy, and targeted validation, we identified 61 distinct variants present at frequencies as low as 0.06%, including 54 variants not previously directly associated with selfish selection. The majority (80%) of variants identified have previously been implicated in developmental disorders and/or oncogenesis and include mutations in six newly associated genes (BRAF, CBL, MAP2K1, MAP2K2, RAF1, and SOS1), all of which encode components of the RAS-MAPK pathway and activate signaling. Our findings extend the link between mutations dysregulating the RAS-MAPK pathway and selfish selection, and show that the aging male germline is a repository for such deleterious mutations.

UroMark-a urinary biomarker assay for the detection of bladder cancer

Background

Bladder cancer (BC) is one of the most common cancers in the western world and ranks as the most expensive to manage, due to the need for cystoscopic examination. BC shows frequent changes in DNA methylation, and several studies have shown the potential utility of urinary biomarkers by detecting epigenetic alterations in voided urine. The aim of this study is to develop a targeted bisulfite next-generation sequencing assay to diagnose BC from urine with high sensitivity and specificity.

Results

We defined a 150 CpG loci biomarker panel from a cohort of 86 muscle-invasive bladder cancers and 30 normal urothelium. Based on this panel, we developed the UroMark assay, a next-generation bisulphite sequencing assay and analysis pipeline for the detection of bladder cancer from urinary sediment DNA. The 150 loci UroMark assay was validated in an independent cohort (n = 274, non-cancer (n = 167) and bladder cancer (n = 107)) voided urine samples with an AUC of 97%. The UroMark classifier sensitivity of 98%, specificity of 97% and NPV of 97% for the detection of primary BC was compared to non-BC urine.

Conclusions

Epigenetic urinary biomarkers for detection of BC have the potential to revolutionise the management of this disease. In this proof of concept study, we show the development and utility of a novel high-throughput, next-generation sequencing-based biomarker for the detection of BC-specific epigenetic alterations in urine.

Electronic supplementary material

The online version of this article (doi:10.1186/s13148-016-0303-5) contains supplementary material, which is available to authorized users.

Identifying gene regulatory elements by genomic microarray mapping of DNaseI hypersensitive sites

The identification of cis-regulatory elements is central to understanding gene transcription. Hypersensitivity of cis-regulatory elements to digestion with DNaseI remains the gold-standard approach to locating such elements. Traditional methods used to identify DNaseI hypersensitive sites are cumbersome and can only be applied to short stretches of DNA at defined locations. Here we report the development of a novel genomic array-based approach to DNaseI hypersensitive site mapping (ADHM) that permits precise, large-scale identification of such sites from as few as 5 million cells. Using ADHM we identified all previously recognized hematopoietic regulatory elements across 200 kb of the mouse T-cell acute lymphocytic leukemia-1 (Tal1) locus, and, in addition, identified two novel elements within the locus, which show transcriptional regulatory activity. We further validated the ADHM protocol by mapping the DNaseI hypersensitive sites across 250 kb of the human TAL1 locus in CD34+ primary stem/progenitor cells and K562 cells and by mapping the previously known DNaseI hypersensitive sites across 240 kb of the human alpha-globin locus in K562 cells. ADHM provides a powerful approach to identifying DNaseI hypersensitive sites across large genomic regions.

The DNA sequence and biological annotation of human chromosome 1

The reference sequence for each human chromosome provides the framework for understanding genome function, variation and evolution. Here we report the finished sequence and biological annotation of human chromosome 1. Chromosome 1 is gene-dense, with 3,141 genes and 991 pseudogenes, and many coding sequences overlap. Rearrangements and mutations of chromosome 1 are prevalent in cancer and many other diseases. Patterns of sequence variation reveal signals of recent selection in specific genes that may contribute to human fitness, and also in regions where no function is evident. Fine-scale recombination occurs in hotspots of varying intensity along the sequence, and is enriched near genes. These and other studies of human biology and disease encoded within chromosome 1 are made possible with the highly accurate annotated sequence, as part of the completed set of chromosome sequences that comprise the reference human genome.

DNA sequence and analysis of human chromosome 9

Chromosome 9 is highly structurally polymorphic. It contains the largest autosomal block of heterochromatin, which is heteromorphic in 6-8% of humans, whereas pericentric inversions occur in more than 1% of the population. The finished euchromatic sequence of chromosome 9 comprises 109,044,351 base pairs and represents >99.6% of the region. Analysis of the sequence reveals many intra- and interchromosomal duplications, including segmental duplications adjacent to both the centromere and the large heterochromatic block. We have annotated 1,149 genes, including genes implicated in male-to-female sex reversal, cancer and neurodegenerative disease, and 426 pseudogenes. The chromosome contains the largest interferon gene cluster in the human genome. There is also a region of exceptionally high gene and G + C content including genes paralogous to those in the major histocompatibility complex. We have also detected recently duplicated genes that exhibit different rates of sequence divergence, presumably reflecting natural selection.

The DNA sequence and comparative analysis of human chromosome 10

The finished sequence of human chromosome 10 comprises a total of 131,666,441 base pairs. It represents 99.4% of the euchromatic DNA and includes one megabase of heterochromatic sequence within the pericentromeric region of the short and long arm of the chromosome. Sequence annotation revealed 1,357 genes, of which 816 are protein coding, and 430 are pseudogenes. We observed widespread occurrence of overlapping coding genes (either strand) and identified 67 antisense transcripts. Our analysis suggests that both inter- and intrachromosomal segmental duplications have impacted on the gene count on chromosome 10. Multispecies comparative analysis indicated that we can readily annotate the protein-coding genes with current resources. We estimate that over 95% of all coding exons were identified in this study. Assessment of single base changes between the human chromosome 10 and chimpanzee sequence revealed nonsense mutations in only 21 coding genes with respect to the human sequence.

The DNA sequence and analysis of human chromosome 13

Chromosome 13 is the largest acrocentric human chromosome. It carries genes involved in cancer including the breast cancer type 2 (BRCA2) and retinoblastoma (RB1) genes, is frequently rearranged in B-cell chronic lymphocytic leukaemia, and contains the DAOA locus associated with bipolar disorder and schizophrenia. We describe completion and analysis of 95.5 megabases (Mb) of sequence from chromosome 13, which contains 633 genes and 296 pseudogenes. We estimate that more than 95.4% of the protein-coding genes of this chromosome have been identified, on the basis of comparison with other vertebrate genome sequences. Additionally, 105 putative non-coding RNA genes were found. Chromosome 13 has one of the lowest gene densities (6.5 genes per Mb) among human chromosomes, and contains a central region of 38 Mb where the gene density drops to only 3.1 genes per Mb.

The DNA sequence and analysis of human chromosome 6

Chromosome 6 is a metacentric chromosome that constitutes about 6% of the human genome. The finished sequence comprises 166,880,988 base pairs, representing the largest chromosome sequenced so far. The entire sequence has been subjected to high-quality manual annotation, resulting in the evidence-supported identification of 1,557 genes and 633 pseudogenes. Here we report that at least 96% of the protein-coding genes have been identified, as assessed by multi-species comparative sequence analysis, and provide evidence for the presence of further, otherwise unsupported exons/genes. Among these are genes directly implicated in cancer, schizophrenia, autoimmunity and many other diseases. Chromosome 6 harbours the largest transfer RNA gene cluster in the genome; we show that this cluster co-localizes with a region of high transcriptional activity. Within the essential immune loci of the major histocompatibility complex, we find HLA-B to be the most polymorphic gene on chromosome 6 and in the human genome.

The DNA sequence and comparative analysis of human chromosome 20

The finished sequence of human chromosome 20 comprises 59,187,298 base pairs (bp) and represents 99.4% of the euchromatic DNA. A single contig of 26 megabases (Mb) spans the entire short arm, and five contigs separated by gaps totalling 320 kb span the long arm of this metacentric chromosome. An additional 234,339 bp of sequence has been determined within the pericentromeric region of the long arm. We annotated 727 genes and 168 pseudogenes in the sequence. About 64% of these genes have a 5' and a 3' untranslated region and a complete open reading frame. Comparative analysis of the sequence of chromosome 20 to whole-genome shotgun-sequence data of two other vertebrates, the mouse Mus musculus and the puffer fish Tetraodon nigroviridis, provides an independent measure of the efficiency of gene annotation, and indicates that this analysis may account for more than 95% of all coding exons and almost all genes.

Integration of cytogenetic landmarks into the draft sequence of the human genome

We have placed 7,600 cytogenetically defined landmarks on the draft sequence of the human genome to help with the characterization of genes altered by gross chromosomal aberrations that cause human disease. The landmarks are large-insert clones mapped to chromosome bands by fluorescence in situ hybridization. Each clone contains a sequence tag that is positioned on the genomic sequence. This genome-wide set of sequence-anchored clones allows structural and functional analyses of the genome. This resource represents the first comprehensive integration of cytogenetic, radiation hybrid, linkage and sequence maps of the human genome; provides an independent validation of the sequence map and framework for contig order and orientation; surveys the genome for large-scale duplications, which are likely to require special attention during sequence assembly; and allows a stringent assessment of sequence differences between the dark and light bands of chromosomes. It also provides insight into large-scale chromatin structure and the evolution of chromosomes and gene families and will accelerate our understanding of the molecular bases of human disease and cancer.

High-resolution landmark framework for the sequence-ready mapping of Xq23-q26.1

We have established a landmark framework map over 20-25 Mb of the long arm of the human X chromosome using yeast artificial chromosome (YAC) clones. The map has approximately one landmark per 45 kb of DNA and stretches from DXS7531 in proximal Xq23 to DXS895 in proximal Xq26, connecting to published framework maps on its proximal and distal sides. There are three gaps in the framework map resulting from the failure to obtain clone coverage from the YAC resources available. Estimates of the maximum sizes of these gaps have been obtained. The four YAC contigs have been positioned and oriented using somatic-cell hybrids and fluorescence in situ hybridization, and the largest is estimated to cover approximately 15 Mb of DNA. The framework map is being used to assemble a sequence-ready map in large-insert bacterial clones, as part of an international effort to complete the sequence of the X chromosome. PAC and BAC contigs currently cover 18 Mb of the region, and from these, 12 Mb of finished sequence is available.