Widespread Transcriptional Scanning in the Testis Modulates Gene Evolution Rates.

The testis expresses the largest number of genes of any mammalian organ, a finding that has long puzzled molecular biologists. Our single-cell transcriptomic data of human and mouse spermatogenesis provide evidence that this widespread transcription maintains DNA sequence integrity in the male germline by correcting DNA damage through a mechanism we term transcriptional scanning. We find that genes expressed during spermatogenesis display lower mutation rates on the transcribed strand and have low diversity in the population. Moreover, this effect is fine-tuned by the level of gene expression during spermatogenesis. The unexpressed genes, which in our model do not benefit from transcriptional scanning, diverge faster over evolutionary timescales and are enriched for sensory and immune-defense functions. Collectively, we propose that transcriptional scanning shapes germline mutation signatures and modulates mutation rates in a gene-specific manner, maintaining DNA sequence integrity for the bulk of genes but allowing for faster evolution in a specific subset.

Phylogenetic debugging of a complete human biosynthetic pathway transplanted into yeast.

Cross-species pathway transplantation enables insight into a biological process not possible through traditional approaches. We replaced the enzymes catalyzing the entire Saccharomyces cerevisiae adenine de novo biosynthesis pathway with the human pathway. While the 'humanized' yeast grew in the absence of adenine, it did so poorly. Dissection of the phenotype revealed that PPAT, the human ortholog of ADE4, showed only partial function whereas all other genes complemented fully. Suppressor analysis revealed other pathways that play a role in adenine de-novo pathway regulation. Phylogenetic analysis pointed to adaptations of enzyme regulation to endogenous metabolite level 'setpoints' in diverse organisms. Using DNA shuffling, we isolated specific amino acids combinations that stabilize the human protein in yeast. Thus, using adenine de novo biosynthesis as a proof of concept, we suggest that the engineering methods used in this study as well as the debugging strategies can be utilized to transplant metabolic pathway from any origin into yeast.

Desmosome mutations impact the tumor microenvironment to promote melanoma proliferation.

Desmosomes are transmembrane protein complexes that contribute to cell-cell adhesion in epithelia and other tissues. Here, we report the discovery of frequent genetic alterations in the desmosome in human cancers, with the strongest signal seen in cutaneous melanoma where desmosomes are mutated in over 70% of cases. In primary but not metastatic melanoma biopsies, the burden of coding mutations on desmosome genes associates with a strong reduction in desmosome gene expression. Analysis by spatial transcriptomics suggests that these expression decreases occur in keratinocytes in the microenvironment rather than in primary melanoma tumor cells. In further support of a microenvironmental origin, we find that loss-of-function knockdowns of the desmosome in keratinocytes yield markedly increased proliferation of adjacent melanocytes in keratinocyte/melanocyte co-cultures. Thus, gradual accumulation of desmosome mutations in neighboring cells may prime melanocytes for neoplastic transformation.

GABA Regulates Electrical Activity and Tumor Initiation in Melanoma.

Oncogenes can initiate tumors only in certain cellular contexts, which is referred to as oncogenic competence. In melanoma, whether cells in the microenvironment can endow such competence remains unclear. Using a combination of zebrafish transgenesis coupled with human tissues, we demonstrate that GABAergic signaling between keratinocytes and melanocytes promotes melanoma initiation by BRAFV600E. GABA is synthesized in melanoma cells, which then acts on GABA-A receptors in keratinocytes. Electron microscopy demonstrates specialized cell-cell junctions between keratinocytes and melanoma cells, and multielectrode array analysis shows that GABA acts to inhibit electrical activity in melanoma/keratinocyte cocultures. Genetic and pharmacologic perturbation of GABA synthesis abrogates melanoma initiation in vivo. These data suggest that GABAergic signaling across the skin microenvironment regulates the ability of oncogenes to initiate melanoma. SIGNIFICANCE: This study shows evidence of GABA-mediated regulation of electrical activity between melanoma cells and keratinocytes, providing a new mechanism by which the microenvironment promotes tumor initiation. This provides insights into the role of the skin microenvironment in early melanomas while identifying GABA as a potential therapeutic target in melanoma. See related commentary by Ceol, p. 2128. This article is featured in Selected Articles from This Issue, p. 2109.

Cancer cell states recur across tumor types and form specific interactions with the tumor microenvironment.

Transcriptional heterogeneity among malignant cells of a tumor has been studied in individual cancer types and shown to be organized into cancer cell states; however, it remains unclear to what extent these states span tumor types, constituting general features of cancer. Here, we perform a pan-cancer single-cell RNA-sequencing analysis across 15 cancer types and identify a catalog of gene modules whose expression defines recurrent cancer cell states including 'stress', 'interferon response', 'epithelial-mesenchymal transition', 'metal response', 'basal' and 'ciliated'. Spatial transcriptomic analysis linked the interferon response in cancer cells to T cells and macrophages in the tumor microenvironment. Using mouse models, we further found that induction of the interferon response module varies by tumor location and is diminished upon elimination of lymphocytes. Our work provides a framework for studying how cancer cell states interact with the tumor microenvironment to form organized systems capable of immune evasion, drug resistance and metastasis.

Cooperation between melanoma cell states promotes metastasis through heterotypic cluster formation.

Melanomas can have multiple coexisting cell states, including proliferative (PRO) versus invasive (INV) subpopulations that represent a "go or grow" trade-off; however, how these populations interact is poorly understood. Using a combination of zebrafish modeling and analysis of patient samples, we show that INV and PRO cells form spatially structured heterotypic clusters and cooperate in the seeding of metastasis, maintaining cell state heterogeneity. INV cells adhere tightly to each other and form clusters with a rim of PRO cells. Intravital imaging demonstrated cooperation in which INV cells facilitate dissemination of less metastatic PRO cells. We identified the TFAP2 neural crest transcription factor as a master regulator of clustering and PRO/INV states. Isolation of clusters from patients with metastatic melanoma revealed a subset with heterotypic PRO-INV clusters. Our data suggest a framework for the co-existence of these two divergent cell populations, in which heterotypic clusters promote metastasis via cell-cell cooperation.

Molecular Insights into Bacteriophage Evolution toward Its Host.

Bacteriophages (phages), viruses that infect bacteria, are considered to be highly host-specific. To add to the knowledge about the evolution and development of bacteriophage speciation toward its host, we conducted a 21-day experiment with the broad host-range bacteriophage Aquamicrobium phage P14. We incubated the phage, which was previously isolated and enriched with the Alphaproteobacteria Aquamicrobium H14, with the Betaproteobacteria Alcaligenaceae H5. During the experiment, we observed an increase in the phage's predation efficacy towards Alcaligenaceae H5. Furthermore, genome analysis and the comparison of the bacteriophage's whole genome indicated that rather than being scattered evenly along the genome, mutations occur in specific regions. In total, 67% of the mutations with a frequency higher than 30% were located in genes that encode tail proteins, which are essential for host recognition and attachment. As control, we incubated the phage with the Alphaproteobacteria Aquamicrobium H8. In both experiments, most of the mutations appeared in the gene encoding the tail fiber protein. However, mutations in the gene encoding the tail tubular protein B were only observed when the phage was incubated with Alcaligenaceae H5. This highlights the phage's tail as a key player in its adaptation to different hosts. We conclude that mutations in the phage's genome were mainly located in tail-related regions. Further investigation is needed to fully characterize the adaptation mechanisms of the Aquamicrobium phage P14.

Author Correction: Integrating microarray-based spatial transcriptomics and single-cell RNA-seq reveals tissue architecture in pancreatic ductal adenocarcinomas.

A Correction to this paper has been published: https://doi.org/10.1038/s41587-019-0392-8.

Integrating microarray-based spatial transcriptomics and single-cell RNA-seq reveals tissue architecture in pancreatic ductal adenocarcinomas.

Single-cell RNA sequencing (scRNA-seq) enables the systematic identification of cell populations in a tissue, but characterizing their spatial organization remains challenging. We combine a microarray-based spatial transcriptomics method that reveals spatial patterns of gene expression using an array of spots, each capturing the transcriptomes of multiple adjacent cells, with scRNA-Seq generated from the same sample. To annotate the precise cellular composition of distinct tissue regions, we introduce a method for multimodal intersection analysis. Applying multimodal intersection analysis to primary pancreatic tumors, we find that subpopulations of ductal cells, macrophages, dendritic cells and cancer cells have spatially restricted enrichments, as well as distinct coenrichments with other cell types. Furthermore, we identify colocalization of inflammatory fibroblasts and cancer cells expressing a stress-response gene module. Our approach for mapping the architecture of scRNA-seq-defined subpopulations can be applied to reveal the interactions inherent to complex tissues.

The Stress-Like Cancer Cell State Is a Consistent Component of Tumorigenesis.

Transcriptional profiling of tumors has revealed a stress-like state among the cancer cells with the concerted expression of genes such as fos, jun, and heat-shock proteins, though this has been controversial given possible dissociation-effects associated with single-cell RNA sequencing. Here, we validate the existence of this state using a combination of zebrafish melanoma modeling, spatial transcriptomics, and human samples. We found that the stress-like subpopulation of cancer cells is present from the early stages of tumorigenesis. Comparing with previously reported single-cell RNA sequencing datasets from diverse cancer types, including triple-negative breast cancer, oligodendroglioma, and pancreatic adenocarcinoma, indicated the conservation of this state during tumorigenesis. We also provide evidence that this state has higher tumor-seeding capabilities and that its induction leads to increased growth under both MEK and BRAF inhibitors. Collectively, our study supports the stress-like cells as a cancer cell state expressing a coherent set of genes and exhibiting drug-resistance properties.

New skin for the old RNA-Seq ceremony: the age of single-cell multi-omics.

New methods for simultaneously quantifying protein and gene expression at the single-cell level have the power to identify cell types and to classify cell populations.

Genome Analysis of a Novel Broad Host Range Proteobacteria Phage Isolated from a Bioreactor Treating Industrial Wastewater.

Bacteriophages are viruses that infect bacteria, and consequently they have a major impact on the development of a microbial population. In this study, the genome of a novel broad host range bacteriophage, Aquamicrobium phage P14, isolated from a wastewater treatment plant, was analyzed. The Aquamicrobium phage P14 was found to infect members of different Proteobacteria classes (Alphaproteobacteria and Betaproteobacteria). This phage contains a 40,551 bp long genome and 60% of its genes had blastx hits. Furthermore, the bacteriophage was found to share more than 50% of its genes with several podoviruses and has the same gene order as other polyvalent bacteriophages. The results obtained in this study led to the conclusion that indeed general features of the genome of the Aquamicrobium phage P14 are shared with other broad host range bacteriophages, however further analysis of the genome is needed in order to identify the specific mechanisms which enable the bacteriophage to infect both Alphaproteobacteria and Betaproteobacteria.

A Single-Cell Transcriptomic Map of the Human and Mouse Pancreas Reveals Inter- and Intra-cell Population Structure.

Although the function of the mammalian pancreas hinges on complex interactions of distinct cell types, gene expression profiles have primarily been described with bulk mixtures. Here we implemented a droplet-based, single-cell RNA-seq method to determine the transcriptomes of over 12,000 individual pancreatic cells from four human donors and two mouse strains. Cells could be divided into 15 clusters that matched previously characterized cell types: all endocrine cell types, including rare epsilon-cells; exocrine cell types; vascular cells; Schwann cells; quiescent and activated stellate cells; and four types of immune cells. We detected subpopulations of ductal cells with distinct expression profiles and validated their existence with immuno-histochemistry stains. Moreover, among human beta- cells, we detected heterogeneity in the regulation of genes relating to functional maturation and levels of ER stress. Finally, we deconvolved bulk gene expression samples using the single-cell data to detect disease-associated differential expression. Our dataset provides a resource for the discovery of novel cell type-specific transcription factors, signaling receptors, and medically relevant genes.