Expressed Barcoding Enables High-Resolution Tracking of the Evolution of Drug Tolerance.

For a majority of patients with non-small cell lung cancer with EGFR mutations, treatment with EGFR inhibitors (EGFRi) induces a clinical response. Despite this initial reduction in tumor size, residual disease persists that leads to disease relapse. Elucidating the preexisting biological differences between sensitive cells and surviving drug-tolerant persister cells and deciphering how drug-tolerant cells evolve in response to treatment could help identify strategies to improve the efficacy of EGFRi. In this study, we tracked the origins and clonal evolution of drug-tolerant cells at a high resolution by using an expressed barcoding system coupled with single-cell RNA sequencing. This platform enabled longitudinal profiling of gene expression and drug sensitivity in response to EGFRi across a large number of clones. Drug-tolerant cells had higher expression of key survival pathways such as YAP and EMT at baseline and could also differentially adapt their gene expression following EGFRi treatment compared with sensitive cells. In addition, drug combinations targeting common downstream components (MAPK) or orthogonal factors (chemotherapy) showed greater efficacy than EGFRi alone, which is attributable to broader targeting of the heterogeneous EGFRi-tolerance mechanisms present in tumors. Overall, this approach facilitates thorough examination of clonal evolution in response to therapy that could inform the development of improved diagnostic approaches and treatment strategies for targeting drug-tolerant cells. SIGNIFICANCE: The evolution and heterogeneity of EGFR inhibitor tolerance are identified in a large number of clones at enhanced cellular and temporal resolution using an expressed barcode technology coupled with single-cell RNA sequencing.

SHP2 blockade enhances anti-tumor immunity via tumor cell intrinsic and extrinsic mechanisms.

SHP2 is a ubiquitous tyrosine phosphatase involved in regulating both tumor and immune cell signaling. In this study, we discovered a novel immune modulatory function of SHP2. Targeting this protein with allosteric SHP2 inhibitors promoted anti-tumor immunity, including enhancing T cell cytotoxic function and immune-mediated tumor regression. Knockout of SHP2 using CRISPR/Cas9 gene editing showed that targeting SHP2 in cancer cells contributes to this immune response. Inhibition of SHP2 activity augmented tumor intrinsic IFNγ signaling resulting in enhanced chemoattractant cytokine release and cytotoxic T cell recruitment, as well as increased expression of MHC Class I and PD-L1 on the cancer cell surface. Furthermore, SHP2 inhibition diminished the differentiation and inhibitory function of immune suppressive myeloid cells in the tumor microenvironment. SHP2 inhibition enhanced responses to anti-PD-1 blockade in syngeneic mouse models. Overall, our study reveals novel functions of SHP2 in tumor immunity and proposes that targeting SHP2 is a promising strategy for cancer immunotherapy.

TGFß-blockade uncovers stromal plasticity in tumors by revealing the existence of a subset of interferon-licensed fibroblasts.

Despite the increasing interest in targeting stromal elements of the tumor microenvironment, we still face tremendous challenges in developing adequate therapeutics to modify the tumor stromal landscape. A major obstacle to this is our poor understanding of the phenotypic and functional heterogeneity of stromal cells in tumors. Herein, we perform an unbiased interrogation of tumor mesenchymal cells, delineating the co-existence of distinct subsets of cancer-associated fibroblasts (CAFs) in the microenvironment of murine carcinomas, each endowed with unique phenotypic features and functions. Furthermore, our study shows that neutralization of TGFß in vivo leads to remodeling of CAF dynamics, greatly reducing the frequency and activity of the myofibroblast subset, while promoting the formation of a fibroblast population characterized by strong response to interferon and heightened immunomodulatory properties. These changes correlate with the development of productive anti-tumor immunity and greater efficacy of PD1 immunotherapy. Along with providing the scientific rationale for the evaluation of TGFß and PD1 co-blockade in the clinical setting, this study also supports the concept of plasticity of the stromal cell landscape in tumors, laying the foundation for future investigations aimed at defining pathways and molecules to program CAF composition for cancer therapy.

Target (MexB)- and Efflux-Based Mechanisms Decreasing the Effectiveness of the Efflux Pump Inhibitor D13-9001 in Pseudomonas aeruginosa PAO1: Uncovering a New Role for MexMN-OprM in Efflux of ß-Lactams and a Novel Regulatory Circuit (MmnRS) Controlling MexMN Expression.

Efflux pumps contribute to antibiotic resistance in Gram-negative pathogens. Correspondingly, efflux pump inhibitors (EPIs) may reverse this resistance. D13-9001 specifically inhibits MexAB-OprM in Pseudomonas aeruginosa Mutants with decreased susceptibility to MexAB-OprM inhibition by D13-9001 were identified, and these fell into two categories: those with alterations in the target MexB (F628L and ΔV177) and those with an alteration in a putative sensor kinase of unknown function, PA1438 (L172P). The alterations in MexB were consistent with reported structural studies of the D13-9001 interaction with MexB. The PA1438L172P alteration mediated a >150-fold upregulation of MexMN pump gene expression and a >50-fold upregulation of PA1438 and the neighboring response regulator gene, PA1437. We propose that these be renamed mmnR and mmnS for MexMN regulator and MexMN sensor, respectively. MexMN was shown to partner with the outer membrane channel protein OprM and to pump several ß-lactams, monobactams, and tazobactam. Upregulated MexMN functionally replaced MexAB-OprM to efflux these compounds but was insusceptible to inhibition by D13-9001. MmnSL172P also mediated a decrease in susceptibility to imipenem and biapenem that was independent of MexMN-OprM. Expression of oprD, encoding the uptake channel for these compounds, was downregulated, suggesting that this channel is also part of the MmnSR regulon. Transcriptome sequencing (RNA-seq) of cells encoding MmnSL172P revealed, among other things, an interrelationship between the regulation of mexMN and genes involved in heavy metal resistance.

Short-term Low-Dose mTORC1 Inhibition in Aged Rats Counter-Regulates Age-Related Gene Changes and Blocks Age-Related Kidney Pathology.

Rapalogs, inhibitors of mTORC1 (mammalian target of rapamycin complex 1), increase life span and delay age-related phenotypes in many species. However, the molecular mechanisms have not been fully elucidated. We determined gene expression changes comparing 6- and 24-month-old rats in the kidney, liver, and skeletal muscle, and asked which of these changes were counter-regulated by a clinically-translatable (short-term and low-concentration) treatment, with a rapalog (RAD001). Surprisingly, RAD001 had a more pronounced effect on the kidney under this regimen in comparison to the liver or skeletal muscle. Histologic evaluation of kidneys revealed that the severity of chronic progressive nephropathy lesions was lower in kidneys from 24-month-old rats treated with RAD001 compared with vehicle. In addition to other gene expression changes, c-Myc, which has been shown to regulate aging, was induced by aging in the kidney and counter-regulated by RAD001. RAD001 caused a decrease in c-Myc protein, which could be rescued by a proteasome inhibitor. These findings point to settings for use of mTORC1 inhibitors to treat age-related disorders, and highlight c-Myc regulation as one of the potential mechanisms by which mTORC1 inhibition is perturbing age-related phenotypes.

Genetic Determinants of Drug Resistance in Mycobacterium tuberculosis and Their Diagnostic Value.

RATIONALE: The development of molecular diagnostics that detect both the presence of Mycobacterium tuberculosis in clinical samples and drug resistance-conferring mutations promises to revolutionize patient care and interrupt transmission by ensuring early diagnosis. However, these tools require the identification of genetic determinants of resistance to the full range of antituberculosis drugs. OBJECTIVES: To determine the optimal molecular approach needed, we sought to create a comprehensive catalog of resistance mutations and assess their sensitivity and specificity in diagnosing drug resistance. METHODS: We developed and validated molecular inversion probes for DNA capture and deep sequencing of 28 drug-resistance loci in M. tuberculosis. We used the probes for targeted sequencing of a geographically diverse set of 1,397 clinical M. tuberculosis isolates with known drug resistance phenotypes. We identified a minimal set of mutations to predict resistance to first- and second-line antituberculosis drugs and validated our predictions in an independent dataset. We constructed and piloted a web-based database that provides public access to the sequence data and prediction tool. MEASUREMENTS AND MAIN RESULTS: The predicted resistance to rifampicin and isoniazid exceeded 90% sensitivity and specificity but was lower for other drugs. The number of mutations needed to diagnose resistance is large, and for the 13 drugs studied it was 238 across 18 genetic loci. CONCLUSIONS: These data suggest that a comprehensive M. tuberculosis drug resistance diagnostic will need to allow for a high dimension of mutation detection. They also support the hypothesis that currently unknown genetic determinants, potentially discoverable by whole-genome sequencing, encode resistance to second-line tuberculosis drugs.

Complete Genome Sequence of Lactobacillus acidophilus FSI4, Isolated from Yogurt.

A new Lactobacillus acidophilus strain, FSI4, isolated from yogurt, was isolated and sequenced in our laboratory. Our data, although supportive of previous conclusions regarding the remarkable stability of L. acidophilus species, indicate accumulating mutations in commercial L. acidophilus strains that warrant further study of the effect of damaged genes on the competitiveness of these bacteria in gut microbiota.

Polymorphisms in the gene that encodes the iron transport protein ferroportin 1 influence susceptibility to tuberculosis.

BACKGROUND: We studied the association between iron intake and polymorphisms in the iron transporter gene SLC40A1 and the risk of tuberculosis. METHODS: We compared iron intake, the frequency of SLC40A1 mutations, and interactions among these variables among 98 tuberculosis patients and 125 controls in Kwazulu-Natal, South Africa. RESULTS: Four SLC40A1 single-nucleotide polymorphisms (SNPs) were associated with an increased risk of tuberculosis and 1 SNP with reduced risk. We also found a gene-environment interaction for 4 nonexonic SNPs and iron intake. CONCLUSIONS: This pilot study demonstrated an association between polymorphisms in SLC40A1 and tuberculosis and provided evidence of an interaction between dietary iron and SLC40A1.

A rare penetrant mutation in CFH confers high risk of age-related macular degeneration.

Two common variants in the gene encoding complement factor H (CFH), the Y402H substitution (rs1061170, c.1204C>T)(1-4) and the intronic rs1410996 SNP(5,6), explain 17% of age-related macular degeneration (AMD) liability. However, proof for the involvement of CFH, as opposed to a neighboring transcript, and knowledge of the potential mechanism of susceptibility alleles are lacking. Assuming that rare functional variants might provide mechanistic insights, we used genotype data and high-throughput sequencing to discover a rare, high-risk CFH haplotype with a c.3628C>T mutation that resulted in an R1210C substitution. This allele has been implicated previously in atypical hemolytic uremic syndrome, and it abrogates C-terminal ligand binding(7,8). Genotyping R1210C in 2,423 AMD cases and 1,122 controls demonstrated high penetrance (present in 40 cases versus 1 control, P = 7.0 × 10(-6)) and an association with a 6-year-earlier onset of disease (P = 2.3 × 10(-6)). This result suggests that loss-of-function alleles at CFH are likely to drive AMD risk. This finding represents one of the first instances in which a common complex disease variant has led to the discovery of a rare penetrant mutation.

Use of whole genome sequencing to estimate the mutation rate of Mycobacterium tuberculosis during latent infection.

Tuberculosis poses a global health emergency, which has been compounded by the emergence of drug-resistant Mycobacterium tuberculosis (Mtb) strains. We used whole-genome sequencing to compare the accumulation of mutations in Mtb isolated from cynomolgus macaques with active, latent or reactivated disease. We sequenced 33 Mtb isolates from nine macaques with an average genome coverage of 93% and an average read depth of 117×. Based on the distribution of SNPs observed, we calculated the mutation rates for these disease states. We found a similar mutation rate during latency as during active disease or in a logarithmically growing culture over the same period of time. The pattern of polymorphisms suggests that the mutational burden in vivo is because of oxidative DNA damage. We show that Mtb continues to acquire mutations during disease latency, which may explain why isoniazid monotherapy for latent tuberculosis is a risk factor for the emergence of isoniazid resistance.

Evaluation of second-generation sequencing of 19 dilated cardiomyopathy genes for clinical applications.

Medical sequencing for diseases with locus and allelic heterogeneities has been limited by the high cost and low throughput of traditional sequencing technologies. "Second-generation" sequencing (SGS) technologies allow the parallel processing of a large number of genes and, therefore, offer great promise for medical sequencing; however, their use in clinical laboratories is still in its infancy. Our laboratory offers clinical resequencing for dilated cardiomyopathy (DCM) using an array-based platform that interrogates 19 of more than 30 genes known to cause DCM. We explored both the feasibility and cost effectiveness of using PCR amplification followed by SGS technology for sequencing these 19 genes in a set of five samples enriched for known sequence alterations (109 unique substitutions and 27 insertions and deletions). While the analytical sensitivity for substitutions was comparable to that of the DCM array (98%), SGS technology performed better than the DCM array for insertions and deletions (90.6% versus 58%). Overall, SGS performed substantially better than did the current array-based testing platform; however, the operational cost and projected turnaround time do not meet our current standards. Therefore, efficient capture methods and/or sample pooling strategies that shorten the turnaround time and decrease reagent and labor costs are needed before implementing this platform into routine clinical applications.

Filter-based hybridization capture of subgenomes enables resequencing and copy-number detection.

To exploit contemporary sequencing technologies for targeted genetic analyses, we developed a hybridization enrichment strategy for DNA capture that uses PCR products as subgenomic traps. We applied this strategy to 115 kilobases of the human genome encompassing 47 genes implicated in cardiovascular disease. Massively parallel sequencing of captured subgenomic libraries interrogated 99.8% of targeted nucleotides >or=20 times ( approximately 40,000-fold enrichment), enabling sensitive and specific detection of sequence variation and copy-number variation.

SOS1 mutations are rare in human malignancies: implications for Noonan Syndrome patients.

Germ line gain-of-function mutations in several members of the RAS/ERK pathway, including PTPN11, KRAS, and RAF1, cause the autosomal dominant genetic disorder Noonan Syndrome (NS). NS patients are at increased risk of leukemia/myeloproliferative disease and possibly some solid tumors, such as neuroblastoma. Recently, SOS1 gain of function mutations have also been shown to cause NS. Somatic PTPN11, KRAS, and RAF1 mutations occur (although at different frequencies) in a variety of sporadic neoplasms, but whether SOS1 mutations are associated with human cancer has not been evaluated. We sequenced DNA from a total of 810 primary malignancies, including pancreatic, lung, breast, and colon carcinomas, and acute myelogenous leukemia, as well as several neuroblastoma cell lines. From this large, diverse series, missense SOS1 mutations were identified in a single pancreatic tumor, one lung adenocarcinoma, and a T-cell acute lymphoblastic leukemia cell line. Our findings suggest that SOS1 is not a significant human oncogene in most cancers. Furthermore, NS patients with SOS1 mutations may not be at increased risk of developing cancer.

Highly multiplexed molecular inversion probe genotyping: over 10,000 targeted SNPs genotyped in a single tube assay.

Large-scale genetic studies are highly dependent on efficient and scalable multiplex SNP assays. In this study, we report the development of Molecular Inversion Probe technology with four-color, single array detection, applied to large-scale genotyping of up to 12,000 SNPs per reaction. While generating 38,429 SNP assays using this technology in a population of 30 trios from the Centre d'Etude Polymorphisme Humain family panel as part of the International HapMap project, we established SNP conversion rates of approximately 90% with concordance rates >99.6% and completeness levels >98% for assays multiplexed up to 12,000plex levels. Furthermore, these individual metrics can be "traded off" and, by sacrificing a small fraction of the conversion rate, the accuracy can be increased to very high levels. No loss of performance is seen when scaling from 6,000plex to 12,000plex assays, strongly validating the ability of the technology to suppress cross-reactivity at high multiplex levels. The results of this study demonstrate the suitability of this technology for comprehensive association studies that use targeted SNPs in indirect linkage disequilibrium studies or that directly screen for causative mutations.