LKB1 drives stasis and C/EBP-mediated reprogramming to an alveolar type II fate in lung cancer.

LKB1 is among the most frequently altered tumor suppressors in lung adenocarcinoma. Inactivation of Lkb1 accelerates the growth and progression of oncogenic KRAS-driven lung tumors in mouse models. However, the molecular mechanisms by which LKB1 constrains lung tumorigenesis and whether the cancer state that stems from Lkb1 deficiency can be reverted remains unknown. To identify the processes governed by LKB1 in vivo, we generated an allele which enables Lkb1 inactivation at tumor initiation and subsequent Lkb1 restoration in established tumors. Restoration of Lkb1 in oncogenic KRAS-driven lung tumors suppressed proliferation and led to tumor stasis. Lkb1 restoration activated targets of C/EBP transcription factors and drove neoplastic cells from a progenitor-like state to a less proliferative alveolar type II cell-like state. We show that C/EBP transcription factors govern a subset of genes that are induced by LKB1 and depend upon NKX2-1. We also demonstrate that a defining factor of the alveolar type II lineage, C/EBPα, constrains oncogenic KRAS-driven lung tumor growth in vivo. Thus, this key tumor suppressor regulates lineage-specific transcription factors, thereby constraining lung tumor development through enforced differentiation.

LKB1 inactivation modulates chromatin accessibility to drive metastatic progression.

Metastasis is the leading cause of cancer-related deaths and enables cancer cells to compromise organ function by expanding in secondary sites. Since primary tumours and metastases often share the same constellation of driver mutations, the mechanisms that drive their distinct phenotypes are unclear. Here we show that inactivation of the frequently mutated tumour suppressor gene LKB1 (encoding liver kinase B1) has evolving effects throughout the progression of lung cancer, which leads to the differential epigenetic re-programming of early-stage primary tumours compared with late-stage metastases. By integrating genome-scale CRISPR-Cas9 screening with bulk and single-cell multi-omic analyses, we unexpectedly identify LKB1 as a master regulator of chromatin accessibility in lung adenocarcinoma primary tumours. Using an in vivo model of metastatic progression, we further show that loss of LKB1 activates the early endoderm transcription factor SOX17 in metastases and a metastatic-like sub-population of cancer cells within primary tumours. The expression of SOX17 is necessary and sufficient to drive a second wave of epigenetic changes in LKB1-deficient cells that enhances metastatic ability. Overall, our study demonstrates how the downstream effects of an individual driver mutation can change throughout cancer development, with implications for stage-specific therapeutic resistance mechanisms and the gene regulatory underpinnings of metastatic evolution.

High-throughput single-cell chromatin accessibility CRISPR screens enable unbiased identification of regulatory networks in cancer.

Chromatin accessibility profiling can identify putative regulatory regions genome wide; however, pooled single-cell methods for assessing the effects of regulatory perturbations on accessibility are limited. Here, we report a modified droplet-based single-cell ATAC-seq protocol for perturbing and evaluating dynamic single-cell epigenetic states. This method (Spear-ATAC) enables simultaneous read-out of chromatin accessibility profiles and integrated sgRNA spacer sequences from thousands of individual cells at once. Spear-ATAC profiling of 104,592 cells representing 414 sgRNA knock-down populations reveals the temporal dynamics of epigenetic responses to regulatory perturbations in cancer cells and the associations between transcription factor binding profiles.

ArchR is a scalable software package for integrative single-cell chromatin accessibility analysis.

The advent of single-cell chromatin accessibility profiling has accelerated the ability to map gene regulatory landscapes but has outpaced the development of scalable software to rapidly extract biological meaning from these data. Here we present a software suite for single-cell analysis of regulatory chromatin in R (ArchR; https://www.archrproject.com/ ) that enables fast and comprehensive analysis of single-cell chromatin accessibility data. ArchR provides an intuitive, user-focused interface for complex single-cell analyses, including doublet removal, single-cell clustering and cell type identification, unified peak set generation, cellular trajectory identification, DNA element-to-gene linkage, transcription factor footprinting, mRNA expression level prediction from chromatin accessibility and multi-omic integration with single-cell RNA sequencing (scRNA-seq). Enabling the analysis of over 1.2 million single cells within 8 h on a standard Unix laptop, ArchR is a comprehensive software suite for end-to-end analysis of single-cell chromatin accessibility that will accelerate the understanding of gene regulation at the resolution of individual cells.

Finding needles in a haystack: dissecting tumor heterogeneity with single-cell transcriptomic and chromatin accessibility profiling.

Tumor evolution often results in a wealth of heterogeneous cancer cell types within a single tumor - heterogeneity that can include epigenetic and gene expression changes that are impossible to identify from histological features alone. The invasion of cancer cells into nearby healthy tissue, accompanied by the infiltration of responding immune cells, results in an even more complex architecture of tumor and non-tumor cells. However, bulk genomics-based methods can only assay the aggregate transcriptomic and epigenetic profiles across all of this rich cellular diversity. Such bulk averaging hides small subpopulations of tumor cells with unique phenotypes that might result in therapeutic resistance or metastatic progression. The advent of single-cell-based genomics assays for measuring transcription and chromatin accessibility - particularly scRNA-seq and scATAC-seq - has enabled the dissection of cell-types within tumors at a scale and resolution capable of unraveling the epigenetic and gene expression programs of rare and unique cellular subpopulations. This Review focuses on recent advances in scRNA-seq and scATAC-seq technologies and their application to cancer biology in the context of furthering our understanding of tumor heterogeneity.

Zmat3 Is a Key Splicing Regulator in the p53 Tumor Suppression Program.

Although TP53 is the most commonly mutated gene in human cancers, the p53-dependent transcriptional programs mediating tumor suppression remain incompletely understood. Here, to uncover critical components downstream of p53 in tumor suppression, we perform unbiased RNAi and CRISPR-Cas9-based genetic screens in vivo. These screens converge upon the p53-inducible gene Zmat3, encoding an RNA-binding protein, and we demonstrate that ZMAT3 is an important tumor suppressor downstream of p53 in mouse KrasG12D-driven lung and liver cancers and human carcinomas. Integrative analysis of the ZMAT3 RNA-binding landscape and transcriptomic profiling reveals that ZMAT3 directly modulates exon inclusion in transcripts encoding proteins of diverse functions, including the p53 inhibitors MDM4 and MDM2, splicing regulators, and components of varied cellular processes. Interestingly, these exons are enriched in NMD signals, and, accordingly, ZMAT3 broadly affects target transcript stability. Collectively, these studies reveal ZMAT3 as a novel RNA-splicing and homeostasis regulator and a key component of p53-mediated tumor suppression.

CRISPR screens in cancer spheroids identify 3D growth-specific vulnerabilities.

Cancer genomics studies have identified thousands of putative cancer driver genes1. Development of high-throughput and accurate models to define the functions of these genes is a major challenge. Here we devised a scalable cancer-spheroid model and performed genome-wide CRISPR screens in 2D monolayers and 3D lung-cancer spheroids. CRISPR phenotypes in 3D more accurately recapitulated those of in vivo tumours, and genes with differential sensitivities between 2D and 3D conditions were highly enriched for genes that are mutated in lung cancers. These analyses also revealed drivers that are essential for cancer growth in 3D and in vivo, but not in 2D. Notably, we found that carboxypeptidase D is responsible for removal of a C-terminal RKRR motif2 from the α-chain of the insulin-like growth factor 1 receptor that is critical for receptor activity. Carboxypeptidase D expression correlates with patient outcomes in patients with lung cancer, and loss of carboxypeptidase D reduced tumour growth. Our results reveal key differences between 2D and 3D cancer models, and establish a generalizable strategy for performing CRISPR screens in spheroids to reveal cancer vulnerabilities.

Massively parallel single-cell chromatin landscapes of human immune cell development and intratumoral T cell exhaustion.

Understanding complex tissues requires single-cell deconstruction of gene regulation with precision and scale. Here, we assess the performance of a massively parallel droplet-based method for mapping transposase-accessible chromatin in single cells using sequencing (scATAC-seq). We apply scATAC-seq to obtain chromatin profiles of more than 200,000 single cells in human blood and basal cell carcinoma. In blood, application of scATAC-seq enables marker-free identification of cell type-specific cis- and trans-regulatory elements, mapping of disease-associated enhancer activity and reconstruction of trajectories of cellular differentiation. In basal cell carcinoma, application of scATAC-seq reveals regulatory networks in malignant, stromal and immune cells in the tumor microenvironment. Analysis of scATAC-seq profiles from serial tumor biopsies before and after programmed cell death protein 1 blockade identifies chromatin regulators of therapy-responsive T cell subsets and reveals a shared regulatory program that governs intratumoral CD8+ T cell exhaustion and CD4+ T follicular helper cell development. We anticipate that scATAC-seq will enable the unbiased discovery of gene regulatory factors across diverse biological systems.

Integrin-ß4 identifies cancer stem cell-enriched populations of partially mesenchymal carcinoma cells.

Neoplastic cells within individual carcinomas often exhibit considerable phenotypic heterogeneity in their epithelial versus mesenchymal-like cell states. Because carcinoma cells with mesenchymal features are often more resistant to therapy and may serve as a source of relapse, we sought to determine whether such cells could be further stratified into functionally distinct subtypes. Indeed, we find that a basal epithelial marker, integrin-ß4 (ITGB4), can be used to enable stratification of mesenchymal-like triple-negative breast cancer (TNBC) cells that differ from one another in their relative tumorigenic abilities. Notably, we demonstrate that ITGB4+ cancer stem cell (CSC)-enriched mesenchymal cells reside in an intermediate epithelial/mesenchymal phenotypic state. Among patients with TNBC who received chemotherapy, elevated ITGB4 expression was associated with a worse 5-year probability of relapse-free survival. Mechanistically, we find that the ZEB1 (zinc finger E-box binding homeobox 1) transcription factor activity in highly mesenchymal SUM159 TNBC cells can repress expression of the epithelial transcription factor TAp63α (tumor protein 63 isoform 1), a protein that promotes ITGB4 expression. In addition, we demonstrate that ZEB1 and ITGB4 are important in modulating the histopathological phenotypes of tumors derived from mesenchymal TNBC cells. Hence, mesenchymal carcinoma cell populations are internally heterogeneous, and ITGB4 is a mechanistically driven prognostic biomarker that can be used to identify the more aggressive subtypes of mesenchymal carcinoma cells in TNBC. The ability to rapidly isolate and mechanistically interrogate the CSC-enriched, partially mesenchymal carcinoma cells should further enable identification of novel therapeutic opportunities to improve the prognosis for high-risk patients with TNBC.

Detection and Identification of Salmonella enterica, Escherichia coli, and Shigella spp. via PCR-electrospray ionization mass spectrometry: isolate testing and analysis of food samples.

An assay to identify the common food-borne pathogens Salmonella, Escherichia coli, Shigella, and Listeria monocytogenes was developed in collaboration with Ibis Biosciences (a division of Abbott Molecular) for the Plex-ID biosensor system, a platform that uses electrospray ionization mass spectroscopy (ESI-MS) to detect the base composition of short PCR amplicons. The new food-borne pathogen (FBP) plate has been experimentally designed using four gene segments for a total of eight amplicon targets. Initial work built a DNA base count database that contains more than 140 Salmonella enterica, 139 E. coli, 11 Shigella, and 36 Listeria patterns and 18 other Enterobacteriaceae organisms. This assay was tested to determine the scope of the assay's ability to detect and differentiate the enteric pathogens and to improve the reference database associated with the assay. More than 800 bacterial isolates of S. enterica, E. coli, and Shigella species were analyzed. Overall, 100% of S. enterica, 99% of E. coli, and 73% of Shigella spp. were detected using this assay. The assay was also able to identify 30% of the S. enterica serovars to the serovar level. To further characterize the assay, spiked food matrices and food samples collected during regulatory field work were also studied. While analysis of preenrichment media was inconsistent, identification of S. enterica from selective enrichment media resulted in serovar-level identifications for 8 of 10 regulatory samples. The results of this study suggest that this high-throughput method may be useful in clinical and regulatory laboratories testing for these pathogens.

Investigation of the Reactivity of Oligodeoxynucleotides with Glyoxal and KMnO(4) Chemical Probes by Electrospray Ionization Mass Spectrometry.

The reactions of two well-known chemical probes, glyoxal and potassium permanganate (KMnO(4)), with oligodeoxynucleotides were monitored by electrospray ionization (ESI) mass spectrometry to evaluate the influence of the sequence of DNA, its secondary structure, and interactions with associated ligands on the reactivity of the two probes. Glyoxal, a guanine-reactive probe, incorporated a mass shift of 58 Da, and potassium permanganate (KMnO(4)) is a thymine-reactive probe that resulted in a mass shift of 34 Da. The reactions depended on the accessibility of the nucleobases, and the peak abundances of the adducts in the ESI-mass spectra were used to quantify the extent of the chemical probe reactions. In this study, both mixed-base sequences were studied as well as control sequences in which one reactive site was located at the terminus or center of the oligodeoxynucleotide while the surrounding bases were a second, different nucleobase. In addition, the reactions of the chemical probes with non-covalent complexes formed between DNA and either actinomycin D or ethidium bromide, both known to interact with single strand DNA, were evaluated.

Characterization of aziridinylbenzoquinone DNA cross-links by liquid chromatography-infrared multiphoton dissociation-mass spectrometry.

DNA cross-linking was evaluated by liquid chromatography-tandem mass spectrometry to determine the relative cross-linking abilities of two aziridinylbenzoquinones. Reactivities of RH1 (2,5-diaziridinyl-3-[hydroxymethyl]-6-methyl-1,4-benzoquinone), a clinically studied antitumor cross-linking agent, and an analogue containing a phenyl group (2,5-diaziridinyl-3-[hydroxymethyl]-6-phenyl-1,4-benzoquinone, PhRH1) rather than a methyl group were compared. The bulky phenyl substituent was added to determine the impact of steric hindrance on the formation of cross-links within a double helical structure. Cross-links formed by RH1 and PhRH1 were observed at 5'-dGNC sites as well as 5'-dGAAC/dGTTC sites. RH1 was more effective at forming cross-links than PhRH1 for a variety of duplexes. Infrared multiphoton dissociation (IRMPD) and collision-induced dissociation results confirmed the presence and the location of the cross-links within the duplexes, and IRMPD was used to identify the dissociation pathways of the cross-linked duplexes.

Examination of the effect of the annealing cation on higher order structures containing guanine or isoguanine repeats.

Isoguanine (2-oxo-6-amino-guanine), a natural but non-standard base, exhibits unique self-association properties compared to its isomer, guanine, and results in formation of different higher order DNA structures. In this work, the higher order structures formed by oligonucleotides containing guanine repeats or isoguanine repeats after annealing in solutions containing various cations are evaluated by electrospray ionization mass spectrometry (ESI-MS) and circular dichroism (CD) spectroscopy. The guanine-containing strand (G9) consistently formed quadruplexes upon annealing, whereas the isoguanine strand (Ig9) formed both pentaplexes and quadruplexes depending on the annealing cation. Quadruplex formation with G9 showed some dependence on the identity of the cation present during annealing with high relative quadruplex formation detected with six of ten cations. Analogous annealing experiments with Ig9 resulted in complex formation with all ten cations, and the majority of the resulting complexes were pentaplexes. CD results indicated most of the original complexes survived the desalting process necessary for ESI-MS analysis. In addition, several complexes, especially the pentaplexes, were found to be capable of cation exchange with ammonium ions. Ab initio calculations were conducted for isoguanine tetrads and pentads coordinated with all ten cations to predict the most energetically stable structures of the complexes in the gas phase. The observed preference of forming quadruplexes versus pentaplexes as a function of the coordinated cation can be interpreted by the calculated reaction energies of both the tetrads and pentads in combination with the distortion energies of tetrads.

Chemogenomic approaches to elucidation of gene function and genetic pathways.

The approximately 6,000 strains in the yeast deletion collection can be studied in a single culture by using a microarray to detect the 20 bp DNA "barcodes" or "tags" contained in each strain. Barcode intensities measured by microarray are compared across time-points or across conditions to analyze the relative fitness of each strain. The development of this pooled fitness assay has greatly facilitated the functional annotation of the yeast genome by making genome-wide gene-deletion studies faster and easier, and has led to the development of high throughput methods for studying drug action in yeast. Pooled screens can be used for identifying gene functions, measuring the functional relatedness of gene pairs to group genes into pathways, identifying drug targets, and determining a drug's mechanism of action. This process involves five main steps: preparing aliquots of pooled cells, pooled growth, isolation of genomic DNA and PCR amplification of the barcodes, array hybridization, and data analysis. In addition to yeast fitness applications, the general method of studying pooled samples with barcode arrays can also be adapted for use with other types of samples, such as mutant collections in other organisms, siRNA vectors, and molecular inversion probes.

Evaluation of binding selectivities and affinities of platinum-based quadruplex interactive complexes by electrospray ionization mass spectrometry.

The quadruplex binding affinities and selectivities of two large pi-surface Pt(II) phenanthroimidazole complexes, as well as a smaller pi-surface platinum bipyridine complex and a larger Ru(II) complex, were evaluated by electrospray ionization mass spectrometry. Circular dichroism (CD) spectroscopy was used to determine the structures of various quadruplexes and to study the thermal denaturation of the quadruplexes in the absence and presence of the metal complexes. In addition, chemical probe reactions with glyoxal were used to monitor the changes in the quadruplex conformation because of association with the complexes. The platinum phenanthroimidazole complexes show increased affinity for several of the quadruplexes with elongated loops between guanine repeats. Quadruplexes with shorter loops exhibited insubstantial binding to the transition metal complexes. Similarly binding to duplex and single strand oligonucleotides was low overall. Although the ruthenium-based metal complex showed somewhat enhanced quadruplex binding, the Pt(II) complexes had higher quadruplex affinities and selectivities that are attributed to their square planar geometries. The chemical probe reactions using glyoxal indicated increased reactivity when the platinum phenanthroimidazole complexes were bound to the quadruplexes, thus suggesting a conformational change that alters guanine accessibility.

ESI-MS characterization of a novel pyrrole-inosine nucleoside that interacts with guanine bases.

Based on binding studies undertaken by electrospray ionization-mass spectrometry, a synthetic pyrrole-inosine nucleoside, 1, capable of forming an extended three-point Hoogsteen-type hydrogen-bonding interaction with guanine, is shown to form specific complexes with two different quadruplex DNA structures [dTG(4)T](4) and d(T(2)G(4))(4) as well as guanine-rich duplex DNA. The binding interactions of two other analogs were evaluated in order to unravel the structural features that contribute to specific DNA recognition. The importance of the Hoogsteen interactions was confirmed through the absence of specific binding when the pyrrole NH hydrogen-bonding site was blocked or removed. While 2, with a large blocking group, was not found to interact with virtually any form of DNA, 3, with the pyrrole functionality missing, was found to interact non-specifically with several types of DNA. The specific binding of 1 to guanine-rich DNA emphasizes the necessity of careful ligand design for specific sequence recognition.

The chemical genomic portrait of yeast: uncovering a phenotype for all genes.

Genetics aims to understand the relation between genotype and phenotype. However, because complete deletion of most yeast genes ( approximately 80%) has no obvious phenotypic consequence in rich medium, it is difficult to study their functions. To uncover phenotypes for this nonessential fraction of the genome, we performed 1144 chemical genomic assays on the yeast whole-genome heterozygous and homozygous deletion collections and quantified the growth fitness of each deletion strain in the presence of chemical or environmental stress conditions. We found that 97% of gene deletions exhibited a measurable growth phenotype, suggesting that nearly all genes are essential for optimal growth in at least one condition.

Genome-wide analysis of barcoded Saccharomyces cerevisiae gene-deletion mutants in pooled cultures.

The availability of a near-complete (96%) collection of gene-deletion mutants in Saccharomyces cerevisiae greatly facilitates the systematic analyses of gene function in yeast. The unique 20 bp DNA 'barcodes' or 'tags' in each deletion strain enable the individual fitness of thousands of deletion mutants to be resolved from a single pooled culture. Here, we present protocols for the study of pooled cultures of tagged yeast deletion mutants with a tag microarray. This process involves five main steps: pooled growth, isolation of genomic DNA, PCR amplification of the barcodes, array hybridization and data analysis. Pooled deletion screening can be used to study gene function, uncover a compound's mode of action and identify drug targets. In addition to these applications, the general method of studying pooled samples with barcode arrays can also be adapted for use with other types of samples, such as mutant collections in other organisms, short interfering RNA vectors and molecular inversion probes.

Head tilt-translation combinations distinguished at the level of neurons.

Angular and linear accelerations of the head occur throughout everyday life, whether from external forces such as in a vehicle or from volitional head movements. The relative timing of the angular and linear components of motion differs depending on the movement. The inner ear detects the angular and linear components with its semicircular canals and otolith organs, respectively, and secondary neurons in the vestibular nuclei receive input from these vestibular organs. Many secondary neurons receive both angular and linear input. Linear information alone does not distinguish between translational linear acceleration and angular tilt, with its gravity-induced change in the linear acceleration vector. Instead, motions are thought to be distinguished by use of both angular and linear information. However, for combined motions, composed of angular tilt and linear translation, the infinite range of possible relative timing of the angular and linear components gives an infinite set of motions among which to distinguish the various types of movement. The present research focuses on motions consisting of angular tilt and horizontal translation, both sinusoidal, where the relative timing, i.e. phase, of the tilt and translation can take any value in the range -180 degrees to 180 degrees . The results show how hypothetical neurons receiving convergent input can distinguish tilt from translation, and that each of these neurons has a preferred combined motion, to which the neuron responds maximally. Also shown are the values of angular and linear response amplitudes and phases that can cause a neuron to be tilt-only or translation-only. Such neurons turn out to be sufficient for distinguishing between combined motions, with all of the possible relative angular-linear phases. Combinations of other neurons, as well, are shown to distinguish motions. Relative response phases and in-phase firing-rate modulation are the key to identifying specific motions from within this infinite set of combined motions.