Regulatory Architecture of the LßT2 Gonadotrope Cell Underlying the Response to Gonadotropin-Releasing Hormone.

The LßT2 mouse pituitary cell line has many characteristics of a mature gonadotrope and is a widely used model system for studying the developmental processes and the response to gonadotropin-releasing hormone (GnRH). The global epigenetic landscape, which contributes to cell-specific gene regulatory mechanisms, and the single-cell transcriptome response variation of LßT2 cells have not been previously investigated. Here, we integrate the transcriptome and genome-wide chromatin accessibility state of LßT2 cells during GnRH stimulation. In addition, we examine cell-to-cell variability in the transcriptional response to GnRH using Gel bead-in-Emulsion Drop-seq technology. Analysis of a bulk RNA-seq data set obtained 45 min after exposure to either GnRH or vehicle identified 112 transcripts that were regulated >4-fold by GnRH (FDR < 0.05). The top regulated transcripts constitute, as determined by Bayesian massive public data integration analysis, a human pituitary-relevant coordinated gene program. Chromatin accessibility [assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq)] data sets generated from GnRH-treated LßT2 cells identified more than 58,000 open chromatin regions, some containing notches consistent with bound transcription factor footprints. The study of the most prominent open regions showed that 75% were in transcriptionally active promoters or introns, supporting their involvement in active transcription. Lhb, Cga, and Egr1 showed significantly open chromatin over their promoters. While Fshb was closed over its promoter, several discrete significantly open regions were found at -40 to -90 kb, which may represent novel upstream enhancers. Chromatin accessibility determined by ATAC-seq was associated with high levels of gene expression determined by RNA-seq. We obtained high-quality single-cell Gel bead-in-Emulsion Drop-seq transcriptome data, with an average of >4,000 expressed genes/cell, from 1,992 vehicle- and 1,889 GnRH-treated cells. While the individual cell expression patterns showed high cell-to-cell variation, representing both biological and measurement variation, the average expression patterns correlated well with bulk RNA-seq data. Computational assignment of each cell to its precise cell cycle phase showed that the response to GnRH was unaffected by cell cycle. To our knowledge, this study represents the first genome-wide epigenetic and single-cell transcriptomic characterization of this important gonadotrope model. The data have been deposited publicly and should provide a resource for hypothesis generation and further study.

Erratum: Lkb1 inactivation drives lung cancer lineage switching governed by Polycomb Repressive Complex 2.

This corrects the article DOI: 10.1038/ncomms14922.

Lkb1 inactivation drives lung cancer lineage switching governed by Polycomb Repressive Complex 2.

Adenosquamous lung tumours, which are extremely poor prognosis, may result from cellular plasticity. Here, we demonstrate lineage switching of KRAS+ lung adenocarcinomas (ADC) to squamous cell carcinoma (SCC) through deletion of Lkb1 (Stk11) in autochthonous and transplant models. Chromatin analysis reveals loss of H3K27me3 and gain of H3K27ac and H3K4me3 at squamous lineage genes, including Sox2, ΔNp63 and Ngfr. SCC lesions have higher levels of the H3K27 methyltransferase EZH2 than the ADC lesions, but there is a clear lack of the essential Polycomb Repressive Complex 2 (PRC2) subunit EED in the SCC lesions. The pattern of high EZH2, but low H3K27me3 mark, is also prevalent in human lung SCC and SCC regions within ADSCC tumours. Using FACS-isolated populations, we demonstrate that bronchioalveolar stem cells and club cells are the likely cells-of-origin for SCC transitioned tumours. These findings shed light on the epigenetics and cellular origins of lineage-specific lung tumours.