Captive wildlife management survey in Vietnam, 2015-2021.

In Vietnam, breeding and raising a wide range of wildlife species in captive wildlife facilities (CWFs) are common practices but little information on the captive wildlife population is available. We conducted surveys and developed software to create a captive wildlife facilities management (CWFM) system. This database provides up-to-date information on the distribution of CWFs, the number of species, and individuals according to the level of protection outlined by the government and the Convention on International Trade in Endangered Species of Wild Fauna and Flora (CITES) categories. CWFs were located in all provinces and regions, but differed in distribution, number of species and individual animals. The Mekong River Delta region recorded the highest number of CWFs (35.3%) and the highest number of animals (43.1%). In 2021, 95 species belong to the highest level of protection group were being raised at 1824 CWFs; 137 species in 4554 CWFs in CITES appendix II, appendix III, government list IIB; and 139 species in 1499 CWFs belong to the common wildlife. The overall number of CWFs in 50 provinces decreased by a negative compound annual growth rate of -7.2%. However, it is crucial to continue to monitor the changing dynamics to assess the risks of disease transmission from zoonoses originating from wildlife. We recommend periodic compulsory reporting of CWF activities using the CWFM system.

Malignant Progression of an Ancestral Bone Marrow Clone Harboring a CIC-NUTM2A Fusion in Isolated Myeloid Sarcoma.

Myeloid sarcoma is a rare condition consisting of extramedullary myeloid blasts found in association with acute myeloid leukemia or, in the absence of bone marrow involvement. We identified an infant with isolated myeloid sarcoma whose bone marrow was negative for involvement by flow cytometry. Sequencing revealed the fusion oncogene CIC-NUTM2A and identified the sarcoma to be clonally evolved from the bone marrow, which carried the fusion despite the absence of pathology. Murine modeling confirmed the ability of the fusion to transform hematopoietic cells and identified receptor tyrosine kinase (RTK) signaling activation consistent with disruption of the CIC transcriptional repressor. These findings extend the definition of CIC-rearranged malignancies to include hematologic disease, provide insight into the mechanism of oncogenesis, and demonstrate the importance of molecular analysis and tracking of bone marrow involvement over the course of treatment in myeloid sarcoma, including patients that lack flow cytometric evidence of leukemia at diagnosis. IMPLICATIONS: This study illustrates molecular involvement of phenotypically normal bone marrow in myeloid sarcoma, which has significant implications in clinical care. Further, it extends the definition of CIC-rearrangements to include hematologic malignancies and shows evidence of RTK activation that may be exploited therapeutically in cancer(s) driven by these fusions.

The genomic landscape of pediatric acute lymphoblastic leukemia.

Acute lymphoblastic leukemia (ALL) is the most common childhood cancer. Here, using whole-genome, exome and transcriptome sequencing of 2,754 childhood patients with ALL, we find that, despite a generally low mutation burden, ALL cases harbor a median of four putative somatic driver alterations per sample, with 376 putative driver genes identified varying in prevalence across ALL subtypes. Most samples harbor at least one rare gene alteration, including 70 putative cancer driver genes associated with ubiquitination, SUMOylation, noncoding transcripts and other functions. In hyperdiploid B-ALL, chromosomal gains are acquired early and synchronously before ultraviolet-induced mutation. By contrast, ultraviolet-induced mutations precede chromosomal gains in B-ALL cases with intrachromosomal amplification of chromosome 21. We also demonstrate the prognostic significance of genetic alterations within subtypes. Intriguingly, DUX4- and KMT2A-rearranged subtypes separate into CEBPA/FLT3- or NFATC4-expressing subgroups with potential clinical implications. Together, these results deepen understanding of the ALL genomic landscape and associated outcomes.

Management of de novo metastatic hormone-sensitive prostate cancer: A comprehensive report of a single-center experience.

BACKGROUND: Upfront docetaxel or novel hormonal agents (NHA) such as abiraterone and enzalutamide have become the standard of care for metastatic hormone sensitive prostate cancer (mHSPC). We evaluated real-world management of patients treated with these agents at a single center. PATIENTS AND METHODS: Patients with de novo mHSPC treated with upfront docetaxel or an NHA between January 2014 and April 2019 at Mount Sinai Health System were included. We evaluated time to next treatment (TTNT), PSA progression free survival (PFS) and overall survival (OS) after initial treatment with these drugs. Kaplan Meier method and multivariable Cox proportional hazards models were used for analysis. We additionally assessed the prognostic value of post-treatment PSA. RESULTS: We identified 94 de novo mHSPC patients; 52 and 42 treated with upfront docetaxel and NHAs, respectively. NHAs were associated with a median TTNT of 20.7 months compared to 10.1 months with docetaxel (log-rank p = 0.023). We also observed median PSA PFS of 19 months for NHAs and 13.2 months for docetaxel (p = 0.069). However, OS between the two treatment groups was unchanged. Among docetaxel treated patients, TTNT was shorter among those with high metastasis burden (9.63 vs 25.5 months, p = 0.026) which was not observed among NHA treated patients (25.1 vs 20.7 months, p = 0.79). Regardless of treatment, lower post-treatment PSA levels were associated with improved TTNT (58.95 vs. 11.57 vs. 9.4 months for PSA ≤0.2, 0.2-0.4, >0.4ng/ml, respectively; p<0.001). CONCLUSION: Real world data demonstrated a shorter duration of treatment with docetaxel than NHAs, reflecting the time-limited nature of docetaxel regimens compared to the treat-till-progression approach of NHAs. While TTNT was generally longer for NHAs than docetaxel, some docetaxel-treated patients achieved significant periods of time off treatment. In addition, the depth of PSA response following combination treatment may hold prognostic value for mHSPC outcomes.

Targeted Next-Generation Sequencing Reveals Exceptionally High Rates of Molecular Driver Mutations in Never-Smokers With Lung Adenocarcinoma.

INTRODUCTION: Historically, high rates of actionable driver mutations have been reported in never-smokers with lung adenocarcinoma (ADC). In the era of modern, comprehensive cancer mutation sequencing, this relationship necessitates a more detailed analysis. METHODS: All Mount Sinai patients between January 1, 2015, and June 1, 2020, with a diagnosis of ADC of any stage with known smoking status who received genomic testing were included. Most patients were analyzed using the Sema4 hotspot panel or the Oncomine Comprehensive Assay version 3 next-generation sequencing (NGS) panel conducted at Sema4. Patients were considered fully genotyped if they were comprehensively analyzed for alterations in EGFR, KRAS, MET, ALK, RET, ROS1, BRAF, NTRK1-3, and ERBB2, otherwise they were considered partially genotyped. RESULTS: Two hundred and thirty-six never-smokers and 671 smokers met the above criteria. Of the never-smokers, 201 (85%) had a driver mutation with 167 (71%) considered actionable (ie, those with US Food and Drug Administration-approved agents). Among smokers, 439 (65%) had an identified driver mutation with 258 (38%) actionable (P < .0001). When comprehensively sequenced, 95% (70/74) of never-smokers had a driver mutation with 78% (58/74) actionable; whereas, for smokers, 75% (135/180) had a driver with only 47% (74/180) actionable (P < .0001). Within mutations groups, EGFR G719X and KRAS G12Cs were more common to smokers. For stage IV patients harboring EGFR-mutant tumors treated with EGFR-directed therapies, never-smokers had significantly improved OS compared to smokers (hazard ratio = 2.71; P = .025). In multivariable analysis, Asian ancestry and female sex remained significant predictors of (1) OS in stage IV patients and (2) likelihood of harboring a receptor of fusion-based driver. CONCLUSION: Comprehensive NGS revealed driver alterations in 95% of never-smokers, with the majority having an associated therapy available. All efforts should be exhausted to identify or rule out the presence of an actionable driver mutation in all metastatic lung ADC.

Antigen cross-presentation in young tumor-bearing hosts promotes CD8+ T cell terminal differentiation.

The immune system undergoes a progressive functional remodeling with age. Understanding how age bias shapes antitumor immunity is essential in designing effective immunotherapies, especially for pediatric patients. Here, we explore antitumor CD8+ T cell responses generated in young (prepubescent) and adult (presenescent) mice. Using an MHCI-deficient tumor model, we observed that tumor-reactive CD8+ T cells expanded in young tumor-bearing (TB) mice acquired a terminally differentiated phenotype characterized by overexpression of inhibitory receptors and the transcription factor Tox1. Furthermore, tumor-infiltrating CD8+ T cells from young tumors yielded a poor cytokine response compared with CD8+ T cells infiltrating adult tumors. Young migratory dendritic cells (migDCs) from the draining lymph nodes (dLNs), and mononuclear phagocytic cells (MPCs) infiltrating young tumors, were more competent in capturing and cross-presenting tumor antigen, leading to enhanced priming of CD8+ T cells in dLNs and their subsequent terminal differentiation in the tumors. Single-cell transcriptional profiling of tumor-infiltrating MPCs demonstrated that young MPCs are polarized toward an inflammatory, effector phenotype. Consistent with our observations in young versus adult TB mice, analysis of immune infiltrates from pediatric solid tumors showed a correlation between tumor-infiltrating CD8+ T cells with an exhaustion phenotype and the frequency of PD-L1-expressing monocytes/macrophages. Collectively, these data indicate that a young tissue microenvironment contributes to the generation of an immune response skewed toward a less pliable terminal effector state, thus narrowing the window for immunotherapeutic interventions.

Impact of body mass index on the efficacy of aromatase inhibitors in patients with metastatic breast cancer.

PURPOSE: Higher levels of estrogen in obese patients may lead to incomplete inhibition by aromatase inhibitors (AIs). The aim of this study was to determine the impact of body mass index (BMI) on efficacy of AIs in patients with metastatic hormone receptor (HR)-positive breast cancer (BC). METHODS: We performed a retrospective chart review of all female patients with metastatic HR-positive BC on an AI in first- or second-line settings and seen at our academic institution between 2001 and 2020. The primary endpoint was progression-free survival (PFS), defined as the time from start of AI to disease progression or death from any cause. RESULTS: We identified 219 patients who had received an AI in the first- or second-line settings for metastatic HR-positive BC and with documented information on BMI. Of the 219 patients, 56% (123) had a low BMI (defined as < 27 kg/m2) and 44% (96) had a high BMI (≥ 27 kg/m2). The median PFS was 21.9 months (95% CI 14.5 to 28.4) in the low BMI group versus 20.2 months (95% CI 14.3 to 27.5) in the high BMI group (p = 0.73). CONCLUSION: While BMI influences efficacy of AIs in the adjuvant setting, our results suggest that in the metastatic setting, BMI may not impact the efficacy of AIs. This discrepancy could be due to other differences in disease characteristics that make complete aromatase inhibition more important in the adjuvant setting when disease burden is the lowest.

Impact of Peste des petits ruminants for sub-Saharan African farmers: A bioeconomic household production model.

Peste des petits ruminants (PPR), a contagious viral disease affecting small ruminants, has been targeted by the global community for eradication within the next 10 years. However, eradication will require substantial money, human resources, coordination among actors, and individual commitments. The objective was to estimate the cost of PPR at a household level, thereby providing information on the potential economic benefits of PPR eradication and the incentives for small ruminant keepers to actively participate in the PPR Global Eradication Programme. This study focused on four sub-Saharan countries: Ethiopia, Malawi, Mali and Rwanda. Publicly available household level data assembled by FAO were used. A bioeconomic model was built to estimate impacts of PPR for a standardized theoretical area, where each household raised an average small ruminant herd or flock. Model outputs were used to estimate, at a household level, income loss due to a PPR outbreak. We constructed various income scenarios to account for the variability of small ruminant income as a proportion of total household annual income. Household income losses ranged from 2% to 40% of total annual income; percentages varied depending on the income scenario and on the gross annual economic impact of PPR on small ruminant production, which ranged from 20% to 80%, based on results of the bioeconomic model. As expected, the more small ruminant production contributed to household income the greater the impact. Estimates provided herein warn decision makers that, given the heterogeneous impact of PPR on household income, the willingness of small ruminant decision makers to participate in the PPR Global Eradication Programme may vary widely and tailored approaches should be devised and implemented.

Crossroads of highly pathogenic H5N1: overlap between wild and domestic birds in the Black Sea-Mediterranean impacts global transmission.

Understanding transmission dynamics that link wild and domestic animals is a key element of predicting the emergence of infectious disease, an event that has highest likelihood of occurring wherever human livelihoods depend on agriculture and animal trade. Contact between poultry and wild birds is a key driver of the emergence of highly pathogenic avian influenza (HPAI), a process that allows for host switching and accelerated reassortment, diversification, and spread of virus between otherwise unconnected regions. This study addresses questions relevant to the spillover of HPAI at a transmission hotspot: what is the nature of the wild bird-poultry interface in Egypt and adjacent Black Sea-Mediterranean countries and how has this contributed to outbreaks occurring worldwide? Using a spatiotemporal model of infection risk informed by satellite tracking of waterfowl and viral phylogenetics, this study identified ecological conditions that contribute to spillover in this understudied region. Results indicated that multiple ducks (Northern Shoveler and Northern Pintail) hosted segments that shared ancestry with HPAI H5 from both clade 2.2.1 and clade 2.3.4 supporting the role of Anseriformes in linking viral populations in East Asia and Africa over large distances. Quantifying the overlap between wild ducks and H5N1-infected poultry revealed an increasing interface in late winter peaking in early spring when ducks expanded their range before migration, with key differences in the timing of poultry contact risk between local and long-distance migrants.

The genomic profiling and MAMLD1 expression in human and canines with Cushing's disease.

BACKGROUND: Cushing's disease (CD) is defined as hypercortisolemia caused by adrenocorticotropic hormone (ACTH)-secreting pituitary adenomas (corticotroph PA) that afflicts humans and dogs. In order to map common aberrant genomic features of CD between humans and dogs, we performed genomic sequencing and immunostaining on corticotroph PA. METHODS: For inclusion, humans and dog were diagnosed with CD. Whole exome sequencing (WES) was conducted on 6 human corticotroph PA. Transcriptome RNA-Seq was performed on 6 human and 7 dog corticotroph PA. Immunohistochemistry (IHC) was complete on 31 human corticotroph PA. Corticotroph PA were compared with normal tissue and between species analysis were also performed. RESULTS: Eight genes (MAMLD1, MNX1, RASEF, TBX19, BIRC5, TK1, GLDC, FAM131B) were significantly (P < 0.05) overexpressed across human and canine corticotroph PA. IHC revealed MAMLD1 to be positively (3+) expressed in the nucleus of ACTH-secreting tumor cells of human corticotroph PA (22/31, 70.9%), but absent in healthy human pituitary glands. CONCLUSIONS: In this small exploratory cohort, we provide the first preliminary insights into profiling the genomic characterizations of human and dog corticotroph PA with respect to MAMLD1 overexpression, a finding of potential direct impact to CD microadenoma diagnosis. Our study also offers a rationale for potential use of the canine model in development of precision therapeutics.

Genomes for Kids: The Scope of Pathogenic Mutations in Pediatric Cancer Revealed by Comprehensive DNA and RNA Sequencing.

Genomic studies of pediatric cancer have primarily focused on specific tumor types or high-risk disease. Here, we used a three-platform sequencing approach, including whole-genome sequencing (WGS), whole-exome sequencing (WES), and RNA sequencing (RNA-seq), to examine tumor and germline genomes from 309 prospectively identified children with newly diagnosed (85%) or relapsed/refractory (15%) cancers, unselected for tumor type. Eighty-six percent of patients harbored diagnostic (53%), prognostic (57%), therapeutically relevant (25%), and/or cancer-predisposing (18%) variants. Inclusion of WGS enabled detection of activating gene fusions and enhancer hijacks (36% and 8% of tumors, respectively), small intragenic deletions (15% of tumors), and mutational signatures revealing of pathogenic variant effects. Evaluation of paired tumor-normal data revealed relevance to tumor development for 55% of pathogenic germline variants. This study demonstrates the power of a three-platform approach that incorporates WGS to interrogate and interpret the full range of genomic variants across newly diagnosed as well as relapsed/refractory pediatric cancers. SIGNIFICANCE: Pediatric cancers are driven by diverse genomic lesions, and sequencing has proven useful in evaluating high-risk and relapsed/refractory cases. We show that combined WGS, WES, and RNA-seq of tumor and paired normal tissues enables identification and characterization of genetic drivers across the full spectrum of pediatric cancers. This article is highlighted in the In This Issue feature, p. 2945.

Network-based systems pharmacology reveals heterogeneity in LCK and BCL2 signaling and therapeutic sensitivity of T-cell acute lymphoblastic leukemia.

T-cell acute lymphoblastic leukemia (T-ALL) is an aggressive hematological malignancy, and novel therapeutics are much needed. Profiling patient leukemia' drug sensitivities ex vivo, we discovered that 44.4% of childhood and 16.7% of adult T-ALL cases exquisitely respond to dasatinib. Applying network-based systems pharmacology analyses to examine signal circuitry, we identified preTCR-LCK activation as the driver of dasatinib sensitivity, and T-ALL-specific LCK dependency was confirmed in genome-wide CRISPR-Cas9 screens. Dasatinib-sensitive T-ALLs exhibited high BCL-XL and low BCL2 activity and venetoclax resistance. Discordant sensitivity of T-ALL to dasatinib and venetoclax is strongly correlated with T-cell differentiation, particularly with the dynamic shift in LCK vs. BCL2 activation. Finally, single-cell analysis identified leukemia heterogeneity in LCK and BCL2 signaling and T-cell maturation stage, consistent with dasatinib response. In conclusion, our results indicate that developmental arrest in T-ALL drives differential activation of preTCR-LCK and BCL2 signaling in this leukemia, providing unique opportunities for targeted therapy.

Do contrasting patterns of migration movements and disease outbreaks between congeneric waterfowl species reflect differing immunity?

Long-distance migrations influence the dynamics of hostpathogen interactions and understanding the role of migratory waterfowl in the spread of the highly pathogenic avian influenza viruses (HPAIV) is important. While wild geese have been associated with outbreak events, disease ecology of closely related species has not been studied to the same extent. The swan goose (Anser cygnoides) and the bar-headed goose (Anser indicus) are congeneric species with distinctly different HPAIV infection records; the former with few and the latter with numerous records. We compared movements of these species, as well as the more distantly related whooper swan (Cygnus cygnus) through their annual migratory cycle to better understand exposure to HPAIV events and how this compares within and between congeneric and noncongeneric species. In spite of their record of fewer infections, swan geese were more likely to come in contact with disease outbreaks than bar-headed geese. We propose two possible explanations: i) frequent prolonged contact with domestic ducks increases innate immunity in swan geese, and/or ii) the stress of high-elevation migration reduces immunity of bar-headed geese. Continued efforts to improve our understanding of species-level pathogen response is critical to assessing disease transmission risk.

Analysis of Real-World Data to Investigate the Impact of Race and Ethnicity on Response to Programmed Cell Death-1 and Programmed Cell Death-Ligand 1 Inhibitors in Advanced Non-Small Cell Lung Cancers.

BACKGROUND: Racial disparities among clinical trial participants present a challenge to assess whether trial results can be generalized into patients representing diverse races and ethnicities. The objective of this study was to evaluate the impact of race and ethnicity on treatment response in patients with advanced non-small cell lung cancer (aNSCLC) treated with programmed cell death-1 (PD-1) or programmed cell death-ligand 1 (PD-L1) inhibitors through analysis of real-world data (RWD). MATERIALS AND METHODS: A retrospective cohort study of 11,138 patients with lung cancer treated at hospitals within the Mount Sinai Health System was performed. Patients with confirmed aNSCLC who received anti-PD-1/PD-L1 treatment were analyzed for clinical outcomes. Our cohort included 249 patients with aNSCLC who began nivolumab, pembrolizumab, or atezolizumab treatment between November 2014 and December 2018. Time-to-treatment discontinuation (TTD) and overall survival (OS) were the analyzed clinical endpoints. RESULTS: After a median follow-up of 14.8 months, median TTD was 7.8 months (95% confidence interval, 5.4-not estimable [NE]) in 75 African American patients versus 4.6 (2.4-7.2) in 110 White patients (hazard ratio [HR], 0.63). Median OS was not reached (18.4-NE) in African American patients versus 11.6 months (9.7-NE) in White patients (HR, 0.58). Multivariable Cox regression conducted with potential confounders confirmed longer TTD (adjusted HR, 0.65) and OS (adjusted HR, 0.60) in African American versus White patients. Similar real-world response rate (42.6% vs. 43.5%) and disease control rate (59.6% vs. 56.5%) were observed in the African American and White patient populations. Further investigation revealed the African American patient group had lower incidence (14.7%) of putative hyperprogressive diseases (HPD) upon anti-PD-1/PD-L1 treatment than the White patient group (24.5%). CONCLUSION: Analysis of RWD showed longer TTD and OS in African American patients with aNSCLC treated with anti-PD-1/PD-L1 inhibitors. Lower incidence of putative HPD is a possible reason for the favorable outcomes in this patient population. IMPLICATIONS FOR PRACTICE: There is a significant underrepresentation of minority patients in randomized clinical trials, and this study demonstrates that real-world data can be used to investigate the impact of race and ethnicity on treatment response. In retrospective analysis of patients with advanced non-small cell lung cancer treated with programmed cell death-1 or programmed cell death-ligand 1 inhibitors, African American patients had significantly longer time-to-treatment discontinuation and longer overall survival. Analysis of real-world data can yield clinical insights and establish a more complete picture of medical interventions in routine clinical practice.

A composite biomarker of neutrophil-lymphocyte ratio and hemoglobin level correlates with clinical response to PD-1 and PD-L1 inhibitors in advanced non-small cell lung cancers.

BACKGROUND: Immune checkpoint inhibitors (ICIs) have been incorporated into various clinical oncology guidelines for systemic treatment of advanced non-small cell lung cancers (aNSCLC). However, less than 50% (and 20%) of the patients responded to the therapy as a first (or second) line of therapy. PD-L1 immunohistochemistry (IHC) is an extensively studied biomarker of response to ICI, but results from this test have equivocal predictive power. In order to identify other biomarkers that support clinical decision-making around whether to treat with ICIs or not, we performed a retrospective study of patients with aNSCLC who underwent ICI-based therapy in the Mount Sinai Health System between 2014 and 2019. METHODS: We analyzed data from standard laboratory tests performed in patients as a part of the routine clinical workup during treatment, including complete blood counts (CBC) and a comprehensive metabolic panel (CMP), to correlate test results with clinical response and survival. RESULTS: Of 11,138 NSCLC patients identified, 249 had been treated with ICIs. We found associations between high neutrophil-to-lymphocyte ratio (NLR ≥ 5) and poor survival in ICI-treated NSCLC. We further observed that sustained high NLR after initiation of treatment had a more profound impact on survival than baseline NLR, regardless of PD-L1 status. Hazard ratios when comparing patients with NLR ≥ 5 vs. NLR < 5 are 1.7 (p = 0.02), 3.4 (p = 4.2 × 10- 8), and 3.9 (p = 1.4 × 10- 6) at baseline, 2-8 weeks, and 8-14 weeks after treatment start, respectively. Mild anemia, defined as hemoglobin (HGB) less than 12 g/dL was correlated with survival independently of NLR. Finally, we developed a composite NLR and HGB biomarker. Patients with pretreatment NLR ≥ 5 and HGB < 12 g/dL had a median overall survival (OS) of 8.0 months (95% CI 4.5-11.5) compared to the rest of the cohort with a median OS not reached (95% CI 15.9-NE, p = 1.8 × 10- 5), and a hazard ratio of 2.6 (95% CI 1.7-4.1, p = 3.5 × 10- 5). CONCLUSIONS: We developed a novel composite biomarker for ICI-based therapy in NSCLC based on routine CBC tests, which may provide meaningful clinical utility to guide treatment decision. The results suggest that treatment of anemia to elevate HGB before initiation of ICI therapy may improve patient outcomes or the use of alternative non-chemotherapy containing regimens.

Identifying modules of cooperating cancer drivers.

Identifying cooperating modules of driver alterations can provide insights into cancer etiology and advance the development of effective personalized treatments. We present Cancer Rule Set Optimization (CRSO) for inferring the combinations of alterations that cooperate to drive tumor formation in individual patients. Application to 19 TCGA cancer types revealed a mean of 11 core driver combinations per cancer, comprising 2-6 alterations per combination and accounting for a mean of 70% of samples per cancer type. CRSO is distinct from methods based on statistical co-occurrence, which we demonstrate is a suboptimal criterion for investigating driver cooperation. CRSO identified well-studied driver combinations that were not detected by other approaches and nominated novel combinations that correlate with clinical outcomes in multiple cancer types. Novel synergies were identified in NRAS-mutant melanomas that may be therapeutically relevant. Core driver combinations involving NFE2L2 mutations were identified in four cancer types, supporting the therapeutic potential of NRF2 pathway inhibition. CRSO is available at https://github.com/mikekleinsgit/CRSO/.

The acquisition of molecular drivers in pediatric therapy-related myeloid neoplasms.

Pediatric therapy-related myeloid neoplasms (tMN) occur in children after exposure to cytotoxic therapy and have a dismal prognosis. The somatic and germline genomic alterations that drive these myeloid neoplasms in children and how they arise have yet to be comprehensively described. We use whole exome, whole genome, and/or RNA sequencing to characterize the genomic profile of 84 pediatric tMN cases (tMDS: n = 28, tAML: n = 56). Our data show that Ras/MAPK pathway mutations, alterations in RUNX1 or TP53, and KMT2A rearrangements are frequent somatic drivers, and we identify cases with aberrant MECOM expression secondary to enhancer hijacking. Unlike adults with tMN, we find no evidence of pre-existing minor tMN clones (including those with TP53 mutations), but rather the majority of cases are unrelated clones arising as a consequence of cytotoxic therapy. These studies also uncover rare cases of lineage switch disease rather than true secondary neoplasms.

SCISSOR: a framework for identifying structural changes in RNA transcripts.

High-throughput sequencing protocols such as RNA-seq have made it possible to interrogate the sequence, structure and abundance of RNA transcripts at higher resolution than previous microarray and other molecular techniques. While many computational tools have been proposed for identifying mRNA variation through differential splicing/alternative exon usage, challenges in its analysis remain. Here, we propose a framework for unbiased and robust discovery of aberrant RNA transcript structures using short read sequencing data based on shape changes in an RNA-seq coverage profile. Shape changes in selecting sample outliers in RNA-seq, SCISSOR, is a series of procedures for transforming and normalizing base-level RNA sequencing coverage data in a transcript independent manner, followed by a statistical framework for its analysis ( https://github.com/hyochoi/SCISSOR ). The resulting high dimensional object is amenable to unsupervised screening of structural alterations across RNA-seq cohorts with nearly no assumption on the mutational mechanisms underlying abnormalities. This enables SCISSOR to independently recapture known variants such as splice site mutations in tumor suppressor genes as well as novel variants that are previously unrecognized or difficult to identify by any existing methods including recurrent alternate transcription start sites and recurrent complex deletions in 3' UTRs.

Exploration of Coding and Non-coding Variants in Cancer Using GenomePaint.

GenomePaint (https://genomepaint.stjude.cloud/) is an interactive visualization platform for whole-genome, whole-exome, transcriptome, and epigenomic data of tumor samples. Its design captures the inter-relatedness between DNA variations and RNA expression, supporting in-depth exploration of both individual cancer genomes and full cohorts. Regulatory non-coding variants can be inspected and analyzed along with coding variants, and their functional impact further explored by examining 3D genome data from cancer cell lines. Further, GenomePaint correlates mutation and expression patterns with patient outcomes, and supports custom data upload. We used GenomePaint to unveil aberrant splicing that disrupts the RING domain of CREBBP, discover cis activation of the MYC oncogene by duplication of the NOTCH1-MYC enhancer in B-lineage acute lymphoblastic leukemia, and explore the inter- and intra-tumor heterogeneity at EGFR in adult glioblastomas. These examples demonstrate that deep multi-omics exploration of individual cancer genomes enabled by GenomePaint can lead to biological insights for follow-up validation.