Search

Supplementary Figure 1. Somatic mutations detected in histiocytic neoplasms. Supplementary Figure 2. Frequencies of known activating kinase mutations in histiocytic neoplasms. Supplementary Figure 3. ARAF mutations and variants of unknown significance detected in BRAFV600E-wildtype, non-Langerhans Cell Histiocytic neoplasms. Supplementary Figure 4. Clinical and histological images of the non-Langerhans cell histiocytosis (non-LCH) lesions from index patients with kinase fusions identified by RNA-seq. Supplementary Figure 5. Gene expression analysis of histiocytic neoplasms by RNA-seq. Supplementary Table 1. Characteristics of patient samples with histiocytic neoplasms used in this study and genomic analysis utilized for each sample. Supplementary Table 2. Whole exome sequencing metrics. Supplementary Table 3. Somatic variants identified by whole exome and whole transcriptome sequencing and validated by droplet-digital PCR and/or custom-capture targeted next-generation sequencing with variant allele frequencies (VAFs). Supplementary Table 4. List of the top 1% of differentially expressed genes across histiocytic samples from mRNA sequencing. Supplementary Table 5. PCR primers with M13F2 and M13R2 tails (blue) for Sanger sequencing used in the targeted sequencing recurrence testing for MAP2K1 exons 2 and 3 and all coding regions of ARAF.

SummaryBackgroundThe emergence of recombinant viruses is a threat to public health. Recombination of viral variants may combine variant-specific features that together catalyze viral escape from treatment or immunity. The selective advantages of recombinant SARS-CoV-2 isolates over their parental lineages remain unknown.MethodsMulti-method amplicon and metagenomic sequencing of a clinical swab and the in vitro grown virus allowed for high-confidence detection of a novel recombinant variant. Mutational, phylogeographic, and structural analyses determined features of the recombinant genome and spike protein. Neutralization assays using infectious as well as pseudotyped viruses and point mutants thereof defined the recombinant’s sensitivity to a panel of monoclonal antibodies and sera from vaccinated and/or convalescent individuals.ResultsA novel Delta-Omicron SARS-CoV-2 recombinant was identified in an unvaccinated, immunosuppressed kidney transplant recipient treated with monoclonal antibody Sotrovimab. The recombination breakpoint is located in the spike N-terminal domain, adjacent to the Sotrovimab quaternary binding site, and results in a 5’-Delta AY.45 and a 3’-Omicron BA.1 mosaic spike protein. Delta and BA.1 are sensitive to Sotrovimab neutralization, whereas the Delta-Omicron recombinant is highly resistant to Sotrovimab, both with and without the RBD resistance mutation E340D.ConclusionsRecombination between circulating SARS-CoV-2 variants can functionally contribute to immune escape. It is critical to validate phenotypes of mosaic viruses and monitor immunosuppressed COVID-19 patients treated with monoclonal antibodies for the selection of recombinant and immune escape variants. (Funded by NYU, the National Institutes of Health, and others)

BackgroundDeveloping genomic resources for a diverse range of species is an important step towards understanding the mechanisms underlying complex traits.Specifically, organisms that exhibit unique, accessible phenotypes-of-interests, allow researchers to address questions that may be ill-suited to traditional model organisms. We sequenced the genome and transcriptome of Alston’s singing mouse (Scotinomys teguina), an emerging model for social cognition and vocal communication. In addition to producing advertisement songs used for mate attraction and male-male competition, these rodents are diurnal, live at high-altitudes, and are obligate insectivores, providing opportunities to explore diverse physiological, ecological, and evolutionary questions.ResultsUsing PromethION, Illumina, and PacBio sequencing, we produced an annotated genome and transcriptome, which were validated using gene expression and functional enrichment analyses. To assess the usefulness of our assemblies, we performed single nuclei sequencing on cells of the orofacial motor cortex, a brain region implicated in song coordination, identifying 12 cell types.ConclusionsThese resources will provide the opportunity to identify the molecular basis of complex traits in singing mice as well as to contribute data that can be used for large-scale comparative analyses.

In 2021, Delta became the predominant SARS-CoV-2 variant worldwide. While vaccines have effectively prevented COVID-19 hospitalization and death, vaccine breakthrough infections increasingly occurred. The precise role of clinical and genomic determinants in Delta infections is not known, and whether they contributed to increased rates of breakthrough infections compared to unvaccinated controls.We studied SARS-CoV-2 variant distribution, dynamics, and adaptive selection over time in relation to vaccine status, phylogenetic relatedness of viruses, full genome mutation profiles, and associated clinical and demographic parameters.We show a steep and near-complete replacement of circulating variants with Delta between May and August 2021 in metropolitan New York. We observed an increase of the Delta sublineage AY.25 (14% in vaccinated, 7% in unvaccinated), its spike mutation S112L, and AY.44 (8% in vaccinated, 2% in unvaccinated) with its nsp12 mutation F192V in breakthroughs. Delta infections were associated with younger age and lower hospitalization rates than Alpha. Delta breakthrough infections increased significantly with time since vaccination, and, after adjusting for confounders, they rose at similar rates as in unvaccinated individuals.We observed a modest adaptation of Delta genomes in breakthrough infections in New York, suggesting an improved genomic framework to support Delta's epidemic growth in times of waning vaccine protection despite limited impact on vaccine escape.The study was supported by NYU institutional funds. The NYULH Genome Technology Center is partially supported by the Cancer Center Support Grant P30CA016087 at the Laura and Isaac Perlmutter Cancer Center.

In 2021, Delta has become the predominant SARS-CoV-2 variant worldwide. While vaccines effectively prevent COVID-19 hospitalization and death, vaccine breakthrough infections increasingly occur. The precise role of clinical and genomic determinants in Delta infections is not known, and whether they contribute to increased rates of breakthrough infections compared to unvaccinated controls. Here, we show a steep and near complete replacement of circulating variants with Delta between May and August 2021 in metropolitan New York. We observed an increase of the Delta sublineage AY.25, its spike mutation S112L, and nsp12 mutation F192V in breakthroughs. Delta infections were associated with younger age and lower hospitalization rates than Alpha. Delta breakthroughs increased significantly with time since vaccination, and, after adjusting for confounders, they rose at similar rates as in unvaccinated individuals. Our data indicate a limited impact of vaccine escape in favor of Delta’s increased epidemic growth in times of waning vaccine protection.

The efficacy of COVID-19 mRNA vaccines is high, but breakthrough infections still occur. We compared the SARS-CoV-2 genomes of 76 breakthrough cases after full vaccination with BNT162b2 (Pfizer/BioNTech), mRNA-1273 (Moderna), or JNJ-78436735 (Janssen) to unvaccinated controls (February-April 2021) in metropolitan New York, including their phylogenetic relationship, distribution of variants, and full spike mutation profiles. Their median age was 48 years; seven required hospitalization and one died. Most breakthrough infections (57/76) occurred with B.1.1.7 (Alpha) or B.1.526 (Iota). Among the 7 hospitalized cases, 4 were infected with B.1.1.7, including 1 death. Both unmatched and matched statistical analyses considering age, sex, vaccine type, and study month as covariates supported the null hypothesis of equal variant distributions between vaccinated and unvaccinated in chi-squared and McNemar tests (p>0.1) highlighting a high vaccine efficacy against B.1.1.7 and B.1.526. There was no clear association among breakthroughs between type of vaccine received and variant. In the vaccinated group, spike mutations in the N-terminal domain and receptor-binding domain that have been associated with immune evasion were overrepresented. The evolving dynamic of SARS-CoV-2 variants requires broad genomic analyses of breakthrough infections to provide real-life information on immune escape mediated by circulating variants and their spike mutations.

Substantial progress has been made in understanding how tumors escape immune surveillance. However, few measures to counteract tumor immune evasion have been developed. Suppression of tumor antigen expression is a common adaptive mechanism that cancers use to evade detection and destruction by the immune system. Epigenetic modifications play a critical role in various aspects of immune invasion, including the regulation of tumor antigen expression. To identify epigenetic regulators of tumor antigen expression, we established a transplantable syngeneic tumor model of immune escape with silenced antigen expression and used this system as a platform for a CRISPR–Cas9 suppressor screen for genes encoding epigenetic modifiers. We found that disruption of the genes encoding either of the chromatin modifiers activating transcription factor 7–interacting protein (Atf7ip) or its interacting partner SET domain bifurcated histone lysine methyltransferase 1 (Setdb1) in tumor cells restored tumor antigen expression. This resulted in augmented tumor immunogenicity concomitant with elevated endogenous retroviral (ERV) antigens and mRNA intron retention. ERV disinhibition was associated with a robust type I interferon response and increased T-cell infiltration, leading to rejection of cells lacking intact Atf7ip or Setdb1. ATF7IP or SETDB1 expression inversely correlated with antigen processing and presentation pathways, interferon signaling, and T-cell infiltration and cytotoxicity in human cancers. Our results provide a rationale for targeting Atf7ip or Setdb1 in cancer immunotherapy.

COVID-19 is a respiratory illness caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and declared by the World Health Organization a global public health emergency. Among the severe outbreaks across South America, Uruguay has become known for curtailing SARS-CoV-2 exceptionally well. To understand the SARS-CoV-2 introductions, local transmissions, and associations with genomic and clinical parameters in Uruguay, we sequenced the viral genomes of 44 outpatients and inpatients in a private healthcare system in its capital, Montevideo, from March to May 2020. We performed a phylogeographic analysis using sequences from our cohort and other studies that indicate a minimum of 23 independent introductions into Uruguay, resulting in five major transmission clusters. Our data suggest that most introductions resulting in chains of transmission originate from other South American countries, with the earliest seeding of the virus in late February 2020, weeks before the borders were closed to all non-citizens and a partial lockdown implemented. Genetic analyses suggest a dominance of S and G clades (G, GH, GR) that make up >90% of the viral strains in our study. In our cohort, lethal outcome of SARS-CoV-2 infection significantly correlated with arterial hypertension, kidney failure, and ICU admission (FDR < 0.01), but not with any mutation in a structural or non-structural protein, such as the spike D614G mutation. Our study contributes genetic, phylodynamic, and clinical correlation data about the exceptionally well-curbed SARS-CoV-2 outbreak in Uruguay, which furthers the understanding of disease patterns and regional aspects of the pandemic in Latin America., SCOPUS: ar.j, info:eu-repo/semantics/published

Effective public response to a pandemic relies upon accurate measurement of the extent and dynamics of an outbreak. Viral genome sequencing has emerged as a powerful approach to link seemingly unrelated cases, and large-scale sequencing surveillance can inform on critical epidemiological parameters. Here, we report the analysis of 236 SARS-CoV2 sequences from cases in the New York City metropolitan area during the initial stages of the 2020 COVID-19 outbreak. The majority of cases throughout the region had no recent travel history or known exposure, and genetically linked cases were spread throughout the region. Comparison to global viral sequences showed that the majority were most related to cases from Europe. Our data are consistent with numerous seed transmissions from multiple sources and a prolonged period of unrecognized community spreading. This work highlights the complementary role of real-time genomic surveillance in addition to traditional epidemiological indicators.

The inheritance of gene expression patterns is dependent on epigenetic regulation, but the establishment and maintenance of epigenetic landscapes during T cell differentiation are incompletely understood. Here we show that two stage-specific Cd4 cis-elements, the previously characterized enhancer E4p and a novel enhancer E4m, coordinately promote Cd4 transcription in mature thymic MHC-II-specific T cells, in part through the canonical Wnt pathway. Specifically, E4p licenses E4m to orchestrate DNA demethylation by TET1 and TET3, which in turn poises the Cd4 locus for transcription in peripheral T cells. Cd4 locus demethylation is important for subsequent Cd4 transcription in activated peripheral T cells wherein these cis-elements become dispensable. By contrast, in developing thymocytes the loss of TET1/3 does not affect Cd4 transcription, highlighting an uncoupled event between transcription and epigenetic modifications. Together our findings reveal an important function for thymic cis-elements in governing gene expression in the periphery via a heritable epigenetic mechanism., The expression of CD4, a critical co-receptor providing T cell help in adaptive immunity, is finely tuned during development. Here the authors show that two enhancer elements, E4p and the newly-defined E4m, coordinate the expression and heritable demethylation of Cd4 in thymocytes but are dispensable for its sustained expression in peripheral T cells.

Atlantic tomcod in the Hudson River Estuary bioaccumulate high hepatic burdens of environmental toxicants. Previously, we demonstrated that Hudson River tomcod developed resistance to TCDD and PCB toxicity probably through strong natural selection during their early life-stages for a variant of the Aryl Hydrocarbon Receptor2 (AHR2). Here, we evaluated the genomic consequences of the resistant genotype by comparing global gene expression in larval tomcod from the Hudson River with expression in larvae from a nearby sensitive population (Shinnecock Bay). We developed an annotated draft tomcod genome to explore the effects of multigenerational exposure to toxicants and a functionally impaired AHR2 on the transcriptome. We used the tomcod genome as a reference in RNA-Seq to compare global gene expression in tomcod larvae from the Hudson River and Shinnecock Bay after experimental exposure of larvae to graded doses of TCDD. We found dramatic differences between offspring from the two populations in the number of genes that were differentially expressed at all doses (0.01, 0.1, and 1 ppb) and even in the vehicle controls. At the two lowest TCDD doses, 250 and 1,141 genes were differentially expressed in Shinnecock Bay larvae compared with 14 and 12, respectively, in Hudson River larvae. At the highest dose (1.0 ppb), 934 genes were differentially expressed in Shinnecock Bay larvae and 173 in Hudson River larvae, but only 28 (16%) of affected genes were shared among both populations. Given the large difference between the two populations in the number and identity of differentially expressed genes, it is likely that the AHR2 pathway interacts directly or indirectly with many genes beyond those known in the AHR2 battery and that other regulatory systems may also respond to TCDD exposure. The effects of chronic multi-generational exposure to environmental toxicants on the genome of Hudson River tomcod are much greater than previously expected.

Introduction In Cameroon, a manifold diversity of HIV strains exists with CRF02_AG and unique recombinant forms (URFs) being the predominant strains. In recent years, a steady increase in URFs and clade F2 viruses has been monitored through partial genome sequencing. There is an information gap in the characterization of emerging URFs along the full genome, which is needed to address the challenges URFs pose towards diagnosis, treatment and HIV‐1 vaccine design. Method Eighteen Cameroonian URFs from samples collected between the years 2000 and 2015 were studied using a newly developed near full genome sequencing (NFGS) protocol based on variable nested RT‐PCRs with a versatile primer set. Near full genomes were characterized for recombination patterns and sequence signatures with possible impact on antiretroviral treatment or Env‐directed immune responses. Third‐generation sequencing (3GS) of near full or half genomes (HGs) gave insight into intra‐patient URF diversity. Results The characterized URFs were composed of a broad variety of subtypes and recombinants including A, F, G, CRF01_AE, CRF02_AG and CRF22_01A1. Phylogenetic analysis unveiled dominant CRF02_AG and F2 recombination patterns. 3GS indicated a high intra‐patient URF diversity with up to four distinct viral sub‐populations present in plasma at the same time. URF pol genomic analysis revealed a number of accessory drug resistance mutations (DRMs) in the ART‐naïve participants. Genotypic env analysis suggests CCR5 usage in 14/18 samples and identified deviations at residues, critical for gp120/gp41 interphase and CD4 binding site broadly neutralizing antibodies in more than half of the studied URFs. V1V2 sites of immune pressure in the human RV144 vaccine study varied in more than a third of URFs. Conclusions This study identified novel mosaic patterns in URFs in Cameroon. In line with the regional predominance of CRF_02AG and the increased prevalence of clade F2, prominent CRF_02AG and F2 background patterns were observed underlying the URFs. In the context of the novel mosaic genomes, the impact of the identified accessory DRMs and Env epitope variations on treatment and immune control remains elusive. The evolving diversity of HIV‐1 URFs in Cameroon requires continuous monitoring to respond to the increasing challenges for diagnosis, antiretroviral treatment and prevention.

Here, we report the draft genome of Streptococcus halitosis sp. nov. strain VT-4, a novel bacterium isolated from the dorsal part of the tongue of a patient with halitosis. The genome comprised 1,880,608 bp with a GC content of 41.0%., Here, we report the draft genome of Streptococcus halitosis sp. nov. strain VT-4, a novel bacterium isolated from the dorsal part of the tongue of a patient with halitosis. The genome comprised 1,880,608 bp with a GC content of 41.0%. There were 1,782 predicted protein-coding genes, including those associated with virulence and antibiotic resistance.

We report here an update to the draft genome sequence of Kluyvera intestini sp. nov. strain GT-16, generated using MinION long-read sequencing technology. The complete genome sequence of the human-derived strain GT-16 measured 5,768,848 bp. An improved high-quality complete genome sequence provides insights into the mobility potential of resistance genes in this species.

We report here an update to the draft genome sequence of Kluyvera intestini sp. nov. strain GT-16, generated using MinION long-read sequencing technology. The complete genome sequence of the human-derived strain GT-16 measured 5,768,848 bp. An improved high-quality complete genome sequence provides insights into the mobility potential of resistance genes in this species.

Near full genome sequencing (NFGS) of HIV-1 is required to assess the genetic composition of HIV-1 strains comprehensively. Population-wide, it enables a determination of the heterogeneity of HIV-1 and the emergence of novel/recombinant strains, while for each individual it constitutes a diagnostic instrument to assist targeted therapeutic measures against viral components. There is still a lack of robust and adaptable techniques for efficient NFGS from miscellaneous HIV-1 subtypes. Using rational primer design, a broad primer set was developed for the amplification and sequencing of diverse HIV-1 group M variants from plasma. Using pure subtypes as well as diverse, unique recombinant forms (URF), variable amplicon approaches were developed for NFGS comprising all functional genes. Twenty-three different genomes composed of subtypes A (A1), B, F (F2), G, CRF01_AE, CRF02_AG, and CRF22_01A1 were successfully determined. The NFGS approach was robust irrespective of viral loads (&ge, 306 copies/mL) and amplification method. Third-generation sequencing (TGS), single genome amplification (SGA), cloning, and bulk sequencing yielded similar outcomes concerning subtype composition and recombinant breakpoint patterns. The introduction of a simple and versatile near full genome amplification, sequencing, and cloning method enables broad application in phylogenetic studies of diverse HIV-1 subtypes and can contribute to personalized HIV therapy and diagnosis.

Histiocytic neoplasms are clonal, hematopoietic disorders characterized by an accumulation of abnormal, monocyte-derived dendritic cells or macrophages in Langerhans cell histiocytosis (LCH) and non-Langerhans cell histiocytosis (non-LCH), respectively. The discovery of BRAFV600E mutations in approximately 50% of these patients provided the first molecular therapeutic target in histiocytosis. However, recurrent driving mutations in the majority of patients with BRAFV600E–wild-type non-LCH are unknown, and recurrent cooperating mutations in non-MAP kinase pathways are undefined for the histiocytic neoplasms. Through combined whole-exome and transcriptome sequencing, we identified recurrent kinase fusions involving BRAF, ALK, and NTRK1, as well as recurrent, activating MAP2K1 and ARAF mutations in patients with BRAFV600E–wild-type non-LCH. In addition to MAP kinase pathway lesions, recurrently altered genes involving diverse cellular pathways were identified. Treatment of patients with MAP2K1- and ARAF-mutated non-LCH using MEK and RAF inhibitors, respectively, resulted in clinical efficacy, demonstrating the importance of detecting and targeting diverse kinase alterations in these disorders. Significance: We provide the first description of kinase fusions in systemic histiocytic neoplasms and activating ARAF and MAP2K1 mutations in non-Langerhans histiocytic neoplasms. Refractory patients with MAP2K1- and ARAF-mutant histiocytoses had clinical responses to MEK inhibition and sorafenib, respectively, highlighting the importance of comprehensive genomic analysis of these disorders. Cancer Discov; 6(2); 154–65. ©2015 AACR. This article is highlighted in the In This Issue feature, p. 109

Histiocytic neoplasms are clonal, hematopoietic disorders characterized by an accumulation of abnormal monocyte-derived dendritic cells or macrophages in Langerhans Cell (LCH) and non-Langerhans (non-LCH) histiocytoses, respectively. The discovery of the BRAF V600E mutation in ~50% of patients with LCH and the non-LCH Erdheim-Chester Disease (ECD) provided the first molecular target in these patients and novel insights into the pathogenesis of these disorders. However, recurrent mutations in the majority of the ~50% of BRAF V600E-wild type patients with non-LCH are unknown. Moreover, recurrent mutations outside of the MAP kinase pathway are undefined throughout histiocytic neoplasms. To address these issues, we performed whole exome sequencing (WES) of frozen biopsies from 24 patients with LCH (n=10) or ECD (n=14) paired with peripheral blood mononuclear cells. 13/24 patients also underwent RNA sequencing (RNA-seq). All mutations in activating kinases were validated by droplet-digital PCR, while targeted-capture next-generation sequencing validated all others. Both adult (n=18; n=2 with LCH) and pediatric cases (n=9; n=8 with LCH) were included. Using combined WES/RNA-seq, activating kinase alterations were identified in 100% of patients. In LCH, 60% and 40% had BRAF V600E and MAP2K1 mutations, respectively. In non-LCH 51%, 14%, 14%, and 7% were BRAFV600E, ARAF, MAP2K1, and NRAS mutant (Fig1A). Overall, a mean of 7 non-synonymous mutations per adult patient was identified (range 1-22) compared with 5 mutations per pediatric patient (range 4-9; p =ns). Mutations affecting diverse cellular processes were found to co-exist with kinase mutations including mutations in epigenetic modifiers and the p38/MAPK pathway. In addition to kinase point mutations, RNA-seq identified recurrent, in-frame kinase fusions-a first for these disorders. All identified fusions were validated using FISH and RT-PCR. This includes novel fusions in BRAF (RNF11-BRAF and CLIP2-BRAF), as well as therapeutically important fusions in ALK (2 separate KIF5B-ALK fusions) and NTRK1 (LMNA-NTRK1;Fig1B). Expression of each fusion in Ba/F3 cells conferred cytokine-independent growth. Importantly, the BRAF fusions were found to be sensitive to MEK inhibition but resistant to vemurafenib while the ALK fusions conferred sensitivity to the ALK inhibitors crizotinib or alectinib. We next interrogated a validation cohort of 37 BRAF V600E-wild type, non-LCH, formalin-fixed, paraffin-embedded tissue samples using targeted mutational profiling for MAP2K1, ARAF, NRAS, KRAS, and PIK3CA. This revealed activating mutations in MAP2K1 (32%; n=12), NRAS (16%; n=6), KRAS (11%; n=4), PIK3CA (8%; n=3), and ARAF (3%; n=1). Three of the investigated non-LCH patients with refractory disease and progressive organ dysfunction were treated with targeted therapies based on the discovery of novel kinase alterations described above. Treatment of 2 refractory MAP2K1- mutant, non-LCH patients with MEK inhibitors (trametinib or cobimetinib) resulted in dramatic clinical improvement (Fig1C). Both patients have been maintained on MEK inhibitor single-agent therapy with a sustained clinical response for >100 days. Further evidence of effective targeted inhibition was found in a refractory ECD patient carrying an ARAF S214A mutation. This patient failed to respond to 3 lines of prior therapies and suffered near blindness due to disease infiltration in the retina and optic nerves. Given a recent report of complete response to sorafenib in a lung cancer patient with an ARAF S214C mutation, we initiated sorafenib. Within 12 weeks, there was improvement in the patientÕs eyesight and decreased infiltrative disease, coinciding with >50% decrease in mutant ARAF DNA in plasma cell-free DNA. Whole exome and transcriptome sequencing identified activating kinase mutations or translocations in all patients with the common downstream effect of activating the MAPK pathway. The preliminary, dramatic, clinical efficacy observed with use of MEK and RAF inhibitors in MAP2K1 - and ARAF-mutated, non-LCH patients further supports the central role of targeting the MAPK pathway in these tumors. The discovery of the discussed mutations and fusions in diverse kinases provides critical new insights into the genetic events central to a spectrum of adult and pediatric histiocytic neoplasms. Figure 1. Figure 1. Disclosures Off Label Use: This abstract describes use of MEK inhibitors (both tremetinib and cobimetinib) as well as sorafenib for MEK1 and ARAF mutant histiocytosis. . Stephens:Foundation Medicine, Inc.: Employment, Equity Ownership. Miller:Foundation Medicine, Inc.: Employment, Equity Ownership. Ross:Foundation Medicine Inc.: Employment. Ali:Foundation Medicine Inc.: Employment. Hyman:Chugai Pharma: Consultancy; Biotherapeutics: Consultancy; Atara: Consultancy, Honoraria.