Faculty Bibliography

Epilepsy is a heterogenous group of disorders defined by recurrent seizure activity due to abnormal synchronized activity of neurons. A growing number of epilepsy cases are believed to be caused by genetic factors and copy number variants (CNV) contribute to up to 5% of epilepsy cases. However, CNVs in epilepsy are usually large deletions or duplications involving multiple neurodevelopmental genes. In patients who underwent seizure focus resection for treatment-resistant epilepsy, whole genome DNA methylation profiling identified 3 main clusters of which one showed strong association with receptor tyrosine kinase (RTK) genes. We identified focal copy number gains involving epidermal growth factor receptor (EGFR) and PDGFRA loci. The dysplastic neurons of cases with amplifications showed marked overexpression of EGFR and PDGFRA, while glial and endothelial cells were negative. Targeted sequencing of regulatory regions and DNA methylation analysis revealed that only enhancer regions of EGFR and gene promoter of PDGFRA were amplified, while coding regions did not show copy number abnormalities or somatic mutations. Somatic focal copy number gains of noncoding regulatory represent a previously unrecognized genetic driver in epilepsy and a mechanism of abnormal activation of RTK genes. Upregulated RTKs provide a potential avenue for therapy in seizure disorders.

PURPOSE/OBJECTIVE:Immunotherapy has recently been shown to improve outcomes for advanced PD-L1-positive triple-negative breast cancer (TNBC) in the Impassion130 trial, leading to FDA approval of the first immune checkpoint inhibitor in combination with taxane chemotherapy. To further develop predictive biomarkers and improve therapeutic efficacy of the combination, interrogation of the tumor immune microenvironment before therapy as well as during each component of treatment is crucial. Here we use single-cell RNA sequencing (scRNA-seq) on tumor biopsies to assess immune cell changes from two patients with advanced TNBC treated in a prospective trial at predefined serial time points, before treatment, on taxane chemotherapy and on chemo-immunotherapy. METHODS:Both patients (one responder and one progressor) received the trial therapy, in cycle 1 nab-paclitaxel given as single agent, in cycle 2 nab-paclitaxel in combination with pembrolizumab. Tumor core biopsies were obtained at baseline, 3 weeks (after cycle 1, chemotherapy alone) and 6 weeks (after cycle 2, chemo-immunotherapy). Single-cell RNA sequencing (scRNA-seq) of both cancer cells and infiltrating immune cells isolated were performed from fresh tumor core biopsy specimens by 10 × chromium sequencing. RESULTS:). In contrast, tumors from the patient with rapid disease progression showed a prevalent and persistent myeloid compartment. CONCLUSIONS:Our study provides a deep cellular analysis of on-treatment changes during chemo-immunotherapy for advanced TNBC, demonstrating not only feasibility of single-cell analyses on serial tumor biopsies but also the heterogeneity of TNBC and differences in on-treatment changes in responder versus progressor.

INTRODUCTION/BACKGROUND:Workers' exposure to metalworking fluid (MWF) has been associated with respiratory disease. As part of a public health investigation of a manufacturing facility, we performed paired environmental and human sampling to evaluate cross-pollination of microbes between environment and host and possible effects on lung pathology present among workers. METHODS:Workplace environmental microbiota was evaluated in air and MWF samples. Human microbiota was evaluated in lung tissue samples from workers with respiratory symptoms found to have lymphocytic bronchiolitis and alveolar ductitis with B-cell follicles and emphysema, lung tissue controls, and in skin, nasal and oral samples from 302 workers from different areas of the facility. In vitro effects of MWF exposure on murine B-cells were assessed. RESULTS:Increased similarity of microbial composition was found between MWF samples and lung tissue samples of case workers compared to controls. Among workers in different locations within the facility, those that worked in machine shop area had skin, nasal and oral microbiota more closely related to the microbiota present in MWF samples. Lung samples from four index cases, and skin and nasal samples from workers in machine shop area were enriched with Pseudomonas, the dominant taxa in MWF. Exposure to used MWF stimulated murine B-cell proliferation in vitro, a hallmark cell subtype found in pathology of index cases. CONCLUSIONS:Evaluation of a manufacturing facility with a cluster of workers with respiratory disease supports cross-pollination of microbes from MWF to humans and suggests the potential for exposure to these microbes to be a health hazard.

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.

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 epi-demiological parameters. Here, we report the analysis of 864 SARS-CoV-2 sequences from cases in the New York City metropolitan area during the COVID-19 outbreak in Spring 2020. The majority of cases had no recent travel history or known exposure, and genetically linked cases were spread throughout the region. Comparison to global viral sequences showed that early transmission was most linked to cases from Europe. Our data are consistent with numerous seeds from multiple sources and a prolonged period of unrecognized community spreading. This work highlights the complementary role of genomic surveillance in addition to traditional epidemiological indicators.

The role of the gut microbiota during coinfection with soil-transmitted helminths (STH) and Plasmodium spp. is poorly understood. We examined peripheral blood and fecal samples from 130 individuals who were either infected with Plasmodium vivax only, coinfected with P. vivax and STH, infected with STH alone, or not infected with either P. vivax or STH. In addition to a complete blood count (CBC) with differential, transcriptional profiling of peripheral blood samples was performed by transcriptome sequencing (RNA-Seq), fecal microbial communities were determined by 16S rRNA gene sequencing, and circulating cytokine levels were measured by bead-based immunoassays. Differences in blood cell counts, including an increased percentage of neutrophils, associated with a transcriptional signature of neutrophil activation, were driven primarily by P. vivax infection. P. vivax infection was also associated with increased levels of interleukin 6 (IL-6), IL-8, and IL-10; these cytokine levels were not affected by STH coinfection. Surprisingly, P. vivax infection was more strongly associated with differences in the microbiota than STH infection. Children infected with only P. vivax exhibited elevated Bacteroides and reduced Prevotella and Clostridiaceae levels, but these differences were not observed in individuals coinfected with STH. We also observed that P. vivax parasitemia was higher in the STH-infected population. When we used machine learning to identify the most important predictors of the P. vivax parasite burden (among P. vivax-infected individuals), bacterial taxa were the strongest predictors of parasitemia. In contrast, circulating transforming growth factor β (TGF-β) was the strongest predictor of the Trichuris trichiura egg burden. This study provides unexpected evidence that the gut microbiota may have a stronger link with P. vivax than with STH infection.IMPORTANCEPlasmodium (malaria) and helminth parasite coinfections are frequent, and both infections can be affected by the host gut microbiota. However, the relationship between coinfection and the gut microbiota is unclear. By performing comprehensive analyses on blood/stool samples from 130 individuals in Colombia, we found that the gut microbiota may have a stronger relationship with the number of P. vivax (malaria) parasites than with the number of helminth parasites infecting a host. Microbiota analysis identified more predictors of the P. vivax parasite burden, whereas analysis of blood samples identified predictors of the helminth parasite burden. These results were unexpected, because we expected each parasite to be associated with greater differences in its biological niche (blood for P. vivax and the intestine for helminths). Instead, we find that bacterial taxa were the strongest predictors of P. vivax parasitemia levels, while circulating TGF-β levels were the strongest predictor of helminth parasite burdens.

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.

The dynamics of viral load (VL) of the 2019 novel coronavirus (severe acute respiratory syndrome coronavirus 2) and its association with different clinical parameters remain poorly characterized in the US patient population. Herein, we investigate associations between VL and parameters, such as severity of symptoms, disposition (admission versus direct discharge), length of hospitalization, admission to the intensive care unit, length of need for oxygen support, and overall survival in a cohort of 205 patients from a tertiary care center in New York City. VL was determined using quantitative PCR and log10 transformed for normalization. Univariate and multivariate regression models were used to test these associations. We found that diagnostic viral load is significantly lower in hospitalized patients than in patients not hospitalized (log10 VL = 3.3 versus 4.0; P = 0.018) after adjusting for age, sex, race, body mass index, and comorbidities. Higher VL was associated with shorter duration of the symptoms in all patients and hospitalized patients only and shorter hospital stay (coefficient = -2.02, -2.61, and -2.18; P < 0.001, P = 0.002, and P = 0.013, respectively). No significant association was noted between VL, admission to intensive care unit, length of oxygen support, and overall survival. Our findings suggest a higher shedding risk in less symptomatic patients, an important consideration for containment strategies in severe acute respiratory syndrome coronavirus 2. Furthermore, we identify a novel association between viral load and history of cancer. Larger studies are warranted to validate our findings.

Splicing alterations are common in disease, such as cancer, where mutations in splicing factor genes are frequently responsible for aberrant splicing. Here we present an alternative mechanism for splicing regulation in T cell acute lymphoblastic leukemia (T-ALL), that involves posttranslational stabilization of the splicing machinery via deubiquitination. We demonstrate there are extensive exon skipping changes in disease affecting proteasomal subunits, cell cycle regulators, and the RNA machinery. We present that the serine/arginine-rich splicing factors (SRSF), controlling exon skipping, are critical for leukemia cell survival. The ubiquitin-specific peptidase 7 (USP7) regulates SRSF6 protein levels via active deubiquitination and USP7 inhibition alters the exon skipping pattern and blocks T-ALL growth. The splicing inhibitor H3B-8800 affects splicing of proteasomal transcripts and proteasome activity and acts synergistically with proteasome inhibitors in inhibiting T-ALL growth. Our study provides the proof-of-principle for regulation of splicing factors via deubiquitination and suggests new therapeutic modalities in T-ALL.