TENET 2.0: Identification of key transcriptional regulators and enhancers in lung adenocarcinoma.

Lung cancer is the leading cause of cancer-related death and lung adenocarcinoma is its most common subtype. Although genetic alterations have been identified as drivers in subsets of lung adenocarcinoma, they do not fully explain tumor development. Epigenetic alterations have been implicated in the pathogenesis of tumors. To identify epigenetic alterations driving lung adenocarcinoma, we used an improved version of the Tracing Enhancer Networks using Epigenetic Traits method (TENET 2.0) in primary normal lung and lung adenocarcinoma cells. We found over 32,000 enhancers that appear differentially activated between normal lung and lung adenocarcinoma. Among the identified transcriptional regulators inactivated in lung adenocarcinoma vs. normal lung, NKX2-1 was linked to a large number of silenced enhancers. Among the activated transcriptional regulators identified, CENPA, FOXM1, and MYBL2 were linked to numerous cancer-specific enhancers. High expression of CENPA, FOXM1, and MYBL2 is particularly observed in a subgroup of lung adenocarcinomas and is associated with poor patient survival. Notably, CENPA, FOXM1, and MYBL2 are also key regulators of cancer-specific enhancers in breast adenocarcinoma of the basal subtype, but they are associated with distinct sets of activated enhancers. We identified individual lung adenocarcinoma enhancers linked to CENPA, FOXM1, or MYBL2 that were associated with poor patient survival. Knockdown experiments of FOXM1 and MYBL2 suggest that these factors regulate genes involved in controlling cell cycle progression and cell division. For example, we found that expression of TK1, a potential target gene of a MYBL2-linked enhancer, is associated with poor patient survival. Identification and characterization of key transcriptional regulators and associated enhancers in lung adenocarcinoma provides important insights into the deregulation of lung adenocarcinoma epigenomes, highlighting novel potential targets for clinical intervention.

Genome-wide integration of microRNA and transcriptomic profiles of differentiating human alveolar epithelial cells.

The alveolar epithelium is comprised of two cell types, alveolar epithelial type 1 (AT1) and type 2 (AT2) cells, the latter being capable of self-renewal and transdifferentiation into AT1 cells for normal maintenance and restoration of epithelial integrity following injury. MicroRNAs (miRNAs) are critical regulators of several biological processes, including cell differentiation; however, their role in establishment/maintenance of cellular identity in adult alveolar epithelium is not well understood. To investigate this question, we performed genome-wide analysis of sequential changes in miRNA and gene expression profiles using a well-established model in which human AT2 (hAT2) cells transdifferentiate into AT1-like cells over time in culture that recapitulates many aspects of transdifferentiation in vivo. We defined three phases of miRNA expression during the transdifferentiation process as "early," "late," and "consistently" changed, which were further subclassified as up- or downregulated. miRNAs with altered expression at all time points during transdifferentiation were the largest subgroup, suggesting the need for consistent regulation of signaling pathways to mediate this process. Target prediction analysis and integration with previously published gene expression data identified glucocorticoid signaling as the top pathway regulated by miRNAs. Serum/glucocorticoid-regulated kinase 1 (SGK1) emerged as a central regulatory factor, whose downregulation correlated temporally with gain of hsa-miR-424 and hsa-miR-503 expression. Functional validation demonstrated specific targeting of these miRNAs to the 3'-untranslated region of SGK1. These data demonstrate the time-related contribution of miRNAs to the alveolar transdifferentiation process and suggest that inhibition of glucocorticoid signaling is necessary to achieve the AT1-like cell phenotype.

Epigenome-wide analysis of DNA methylation in lung tissue shows concordance with blood studies and identifies tobacco smoke-inducible enhancers.

Smoking-associated DNA hypomethylation has been observed in blood cells and linked to lung cancer risk. However, its cause and mechanistic relationship to lung cancer remain unclear. We studied the association between tobacco smoking and epigenome-wide methylation in non-tumor lung (NTL) tissue from 237 lung cancer cases in the Environment And Genetics in Lung cancer Etiology study, using the Infinium HumanMethylation450 BeadChip. We identified seven smoking-associated hypomethylated CpGs (P < 1.0 × 10-7), which were replicated in NTL data from The Cancer Genome Atlas. Five of these loci were previously reported as hypomethylated in smokers' blood, suggesting that blood-based biomarkers can reflect changes in the target tissue for these loci. Four CpGs border sequences carrying aryl hydrocarbon receptor binding sites and enhancer-specific histone modifications in primary alveolar epithelium and A549 lung adenocarcinoma cells. A549 cell exposure to cigarette smoke condensate increased these enhancer marks significantly and stimulated expression of predicted target xenobiotic response-related genes AHRR (P = 1.13 × 10-62) and CYP1B1 (P < 2.49 × 10-61). Expression of both genes was linked to smoking-related transversion mutations in lung tumors. Thus, smoking-associated hypomethylation may be a consequence of enhancer activation, revealing environmentally-induced regulatory elements implicated in lung carcinogenesis.

CLDN18.1 attenuates malignancy and related signaling pathways of lung adenocarcinoma in vivo and in vitro.

Claudins are a family of transmembrane proteins integral to the structure and function of tight junctions (TJ). Disruption of TJ and alterations in claudin expression are important features of invasive and metastatic cancer cells. Expression of CLDN18.1, the lung-specific isoform of CLDN18, is markedly decreased in lung adenocarcinoma (LuAd). Furthermore, we recently observed that aged Cldn18 -/- mice have increased propensity to develop LuAd. We now demonstrate that CLDN18.1 expression correlates inversely with promoter methylation and with LuAd patient mortality. In addition, when restored in LuAd cells that have lost expression, CLDN18.1 markedly attenuates malignant properties including xenograft tumor growth in vivo as well as cell proliferation, migration, invasion and anchorage-independent colony formation in vitro. Based on high throughput analyses of Cldn18 -/- murine lung alveolar epithelial type II cells, as well as CLDN18.1-repleted human LuAd cells, we hypothesized and subsequently confirmed by Western analysis that CLDN18.1 inhibits insulin-like growth factor-1 receptor (IGF-1R) and AKT phosphorylation. Consistent with recent data in Cldn18 -/- knockout mice, expression of CLDN18.1 in human LuAd cells also decreased expression of transcriptional co-activator with PDZ-binding motif (TAZ) and Yes-associated protein (YAP) and their target genes, contributing to its tumor suppressor activity. Moreover, analysis of LuAd cells in which YAP and/or TAZ are silenced with siRNA suggests that inhibition of TAZ, and possibly YAP, is also involved in CLDN18.1-mediated AKT inactivation. Taken together, these data indicate a tumor suppressor role for CLDN18.1 in LuAd mediated by a regulatory network that encompasses YAP/TAZ, IGF-1R and AKT signaling.

Cross-Species Transcriptome Profiling Identifies New Alveolar Epithelial Type I Cell-Specific Genes.

Diseases involving the distal lung alveolar epithelium include chronic obstructive pulmonary disease, idiopathic pulmonary fibrosis, and lung adenocarcinoma. Accurate labeling of specific cell types is critical for determining the contribution of each to the pathogenesis of these diseases. The distal lung alveolar epithelium is composed of two cell types, alveolar epithelial type 1 (AT1) and type 2 (AT2) cells. Although cell type-specific markers, most prominently surfactant protein C, have allowed detailed lineage tracing studies of AT2 cell differentiation and the cells' roles in disease, studies of AT1 cells have been hampered by a lack of genes with expression unique to AT1 cells. In this study, we performed genome-wide expression profiling of multiple rat organs together with purified rat AT2, AT1, and in vitro differentiated AT1-like cells, resulting in the identification of 54 candidate AT1 cell markers. Cross-referencing with genes up-regulated in human in vitro differentiated AT1-like cells narrowed the potential list to 18 candidate genes. Testing the top four candidate genes at RNA and protein levels revealed GRAM domain 2 (GRAMD2), a protein of unknown function, as highly specific to AT1 cells. RNA sequencing (RNAseq) confirmed that GRAMD2 is transcriptionally silent in human AT2 cells. Immunofluorescence verified that GRAMD2 expression is restricted to the plasma membrane of AT1 cells and is not expressed in bronchial epithelial cells, whereas reverse transcription-polymerase chain reaction confirmed that it is not expressed in endothelial cells. Using GRAMD2 as a new AT1 cell-specific gene will enhance AT1 cell isolation, the investigation of alveolar epithelial cell differentiation potential, and the contribution of AT1 cells to distal lung diseases.

The 78-kD Glucose-Regulated Protein Regulates Endoplasmic Reticulum Homeostasis and Distal Epithelial Cell Survival during Lung Development.

Bronchopulmonary dysplasia (BPD), a chronic lung disease of prematurity, has been linked to endoplasmic reticulum (ER) stress. To investigate a causal role for ER stress in BPD pathogenesis, we generated conditional knockout (KO) mice (cGrp78(f/f)) with lung epithelial cell-specific KO of Grp78, a gene encoding the ER chaperone 78-kD glucose-regulated protein (GRP78), a master regulator of ER homeostasis and the unfolded protein response (UPR). Lung epithelial-specific Grp78 KO disrupted lung morphogenesis, causing developmental arrest, increased alveolar epithelial type II cell apoptosis, and decreased surfactant protein and type I cell marker expression in perinatal lungs. cGrp78(f/f) pups died immediately after birth, likely owing to respiratory distress. Importantly, Grp78 KO triggered UPR activation with marked induction of the proapoptotic transcription factor CCAAT/enhancer-binding proteins (C/EBP) homologous protein (CHOP). Increased expression of genes involved in oxidative stress and cell death and decreased expression of genes encoding antioxidant enzymes suggest a role for oxidative stress in alveolar epithelial cell (AEC) apoptosis. Increased Smad3 phosphorylation and expression of transforming growth factor-ß/Smad3 targets Cdkn1a (encoding p21) and Gadd45a suggest that interactions among the apoptotic arm of the UPR, oxidative stress, and transforming growth factor-ß/Smad signaling pathways contribute to Grp78 KO-induced AEC apoptosis and developmental arrest. Chemical chaperone Tauroursodeoxycholic acid reduced UPR activation and apoptosis in cGrp78(f/f) lungs cultured ex vivo, confirming a role for ER stress in observed AEC abnormalities. These results demonstrate a key role for GRP78 in AEC survival and gene expression during lung development through modulation of ER stress, and suggest the UPR as a potential therapeutic target in BPD.

Combinations of differentiation markers distinguish subpopulations of alveolar epithelial cells in adult lung.

Distal lung epithelium is maintained by proliferation of alveolar type II (AT2) cells and, for some daughter AT2 cells, transdifferentiation into alveolar type I (AT1) cells. We investigated if subpopulations of alveolar epithelial cells (AEC) exist that represent various stages in transdifferentiation from AT2 to AT1 cell phenotypes in normal adult lung and if they can be identified using combinations of cell-specific markers. Immunofluorescence microscopy showed that, in distal rat and mouse lungs, ∼ 20-30% of NKX2.1(+) (or thyroid transcription factor 1(+)) cells did not colocalize with pro-surfactant protein C (pro-SP-C), a highly specific AT2 cell marker. In distal rat lung, NKX2.1(+) cells coexpressed either pro-SP-C or the AT1 cell marker homeodomain only protein x (HOPX). Not all HOPX(+) cells colocalize with the AT1 cell marker aquaporin 5 (AQP5), and some AQP5(+) cells were NKX2.1(+). HOPX was expressed earlier than AQP5 during transdifferentiation in rat AEC primary culture, with robust expression of both by day 7. We speculate that NKX2.1 and pro-SP-C colocalize in AT2 cells, NKX2.1 and HOPX or AQP5 colocalize in intermediate or transitional cells, and HOPX and AQP5 are expressed without NKX2.1 in AT1 cells. These findings suggest marked heterogeneity among cells previously identified as exclusively AT1 or AT2 cells, implying the presence of subpopulations of intermediate or transitional AEC in normal adult lung.

Integrated transcriptomic and epigenomic analysis of primary human lung epithelial cell differentiation.

Elucidation of the epigenetic basis for cell-type specific gene regulation is key to gaining a full understanding of how the distinct phenotypes of differentiated cells are achieved and maintained. Here we examined how epigenetic changes are integrated with transcriptional activation to determine cell phenotype during differentiation. We performed epigenomic profiling in conjunction with transcriptomic profiling using in vitro differentiation of human primary alveolar epithelial cells (AEC). This model recapitulates an in vivo process in which AEC transition from one differentiated cell type to another during regeneration following lung injury. Interrogation of histone marks over time revealed enrichment of specific transcription factor binding motifs within regions of changing chromatin structure. Cross-referencing of these motifs with pathways showing transcriptional changes revealed known regulatory pathways of distal alveolar differentiation, such as the WNT and transforming growth factor beta (TGFB) pathways, and putative novel regulators of adult AEC differentiation including hepatocyte nuclear factor 4 alpha (HNF4A), and the retinoid X receptor (RXR) signaling pathways. Inhibition of the RXR pathway confirmed its functional relevance for alveolar differentiation. Our incorporation of epigenetic data allowed specific identification of transcription factors that are potential direct upstream regulators of the differentiation process, demonstrating the power of this approach. Integration of epigenomic data with transcriptomic profiling has broad application for the identification of regulatory pathways in other models of differentiation.

Knockout mice reveal key roles for claudin 18 in alveolar barrier properties and fluid homeostasis.

Claudin proteins are major constituents of epithelial and endothelial tight junctions (TJs) that regulate paracellular permeability to ions and solutes. Claudin 18, a member of the large claudin family, is highly expressed in lung alveolar epithelium. To elucidate the role of claudin 18 in alveolar epithelial barrier function, we generated claudin 18 knockout (C18 KO) mice. C18 KO mice exhibited increased solute permeability and alveolar fluid clearance (AFC) compared with wild-type control mice. Increased AFC in C18 KO mice was associated with increased ß-adrenergic receptor signaling together with activation of cystic fibrosis transmembrane conductance regulator, higher epithelial sodium channel, and Na-K-ATPase (Na pump) activity and increased Na-K-ATPase ß1 subunit expression. Consistent with in vivo findings, C18 KO alveolar epithelial cell (AEC) monolayers exhibited lower transepithelial electrical resistance and increased solute and ion permeability with unchanged ion selectivity. Claudin 3 and claudin 4 expression was markedly increased in C18 KO mice, whereas claudin 5 expression was unchanged and occludin significantly decreased. Microarray analysis revealed changes in cytoskeleton-associated gene expression in C18 KO mice, consistent with observed F-actin cytoskeletal rearrangement in AEC monolayers. These findings demonstrate a crucial nonredundant role for claudin 18 in the regulation of alveolar epithelial TJ composition and permeability properties. Increased AFC in C18 KO mice identifies a role for claudin 18 in alveolar fluid homeostasis beyond its direct contributions to barrier properties that may, at least in part, compensate for increased permeability.

Claudin 4 knockout mice: normal physiological phenotype with increased susceptibility to lung injury.

Claudins are tight junction proteins that regulate paracellular ion permeability of epithelium and endothelium. Claudin 4 has been reported to function as a paracellular sodium barrier and is one of three major claudins expressed in lung alveolar epithelial cells (AEC). To directly assess the role of claudin 4 in regulation of alveolar epithelial barrier function and fluid homeostasis in vivo, we generated claudin 4 knockout (Cldn4 KO) mice. Unexpectedly, Cldn4 KO mice exhibited normal physiological phenotype although increased permeability to 5-carboxyfluorescein and decreased alveolar fluid clearance were noted. Cldn4 KO AEC monolayers exhibited unchanged ion permeability, higher solute permeability, and lower short-circuit current compared with monolayers from wild-type mice. Claudin 3 and 18 expression was similar between wild-type and Cldn4 KO alveolar epithelial type II cells. In response to either ventilator-induced lung injury or hyperoxia, claudin 4 expression was markedly upregulated in wild-type mice, whereas Cldn4 KO mice showed greater degrees of lung injury. RNA sequencing, in conjunction with differential expression and upstream analysis after ventilator-induced lung injury, suggested Egr1, Tnf, and Il1b as potential mediators of increased lung injury in Cldn4 KO mice. These results demonstrate that claudin 4 has little effect on normal lung physiology but may function to protect against acute lung injury.

Cell of origin alters myeloid-mediated immunosuppression in lung adenocarcinoma.

Solid carcinomas are often highly heterogenous cancers, arising from multiple epithelial cells of origin. Yet, how the cell of origin influences the response of the tumor microenvironment is poorly understood. Lung adenocarcinoma (LUAD) arises in the distal alveolar epithelium which is populated primarily by alveolar epithelial type I (AT1) and type II (AT2) cells. It has been previously reported that Gramd2 + AT1 cells can give rise to a histologically-defined LUAD that is distinct in pathology and transcriptomic identity from that arising from Sftpc + AT2 cells1,2. To determine how cells of origin influence the tumor immune microenvironment (TIME) landscape, we comprehensively characterized transcriptomic, molecular, and cellular states within the TIME of Gramd2 + AT1 and Sftpc + AT2-derived LUAD using KRASG12D oncogenic driver mouse models. Myeloid cells within the Gramd2 + AT1-derived LUAD TIME were increased, specifically, immunoreactive monocytes and tumor associated macrophages (TAMs). In contrast, the Sftpc + AT2 LUAD TIME was enriched for Arginase-1+ myeloid derived suppressor cells (MDSC) and TAMs expressing profiles suggestive of immunosuppressive function. Validation of immune infiltration was performed using flow cytometry, and intercellular interaction analysis between the cells of origin and major myeloid cell populations indicated that cell-type specific markers SFTPD in AT2 cells and CAV1 in AT1 cells mediated unique interactions with myeloid cells of the differential immunosuppressive states within each cell of origin mouse model. Taken together, Gramd2 + AT1-derived LUAD presents with an anti-tumor, immunoreactive TIME, while the TIME of Sftpc + AT2-derived LUAD has hallmarks of immunosuppression. This study suggests that LUAD cell of origin influences the composition and suppression status of the TIME landscape and may hold critical implications for patient response to immunotherapy.

Culture Conditions Differentially Regulate the Inflammatory Niche and Cellular Phenotype of Tracheo-Bronchial Basal Stem Cells.

Human bronchial epithelial cells (HBECs) derived from the tracheo-bronchial regions of human airways provide an excellent in vitro model for studying pathological mechanisms and evaluating therapeutics in human airway cells. This cell population comprises a mixed population of basal cells (BCs), the predominant stem cell in airways capable of both self-renewal and functional differentiation. Despite their potential for regenerative medicine, BCs exhibit significant phenotypic variability in culture. To investigate how culture conditions influence BC phenotype and function, we expanded three independent BC isolates in three media, airway epithelial cell growth medium (AECGM), dual-SMAD inhibitor (DSI)-enriched AECGM, and Pneumacult Ex plus (PEx+). Extensive RNA sequencing, immune assays and electrical measurements revealed that PEx+ media significantly drove cell proliferation and a broad pro-inflammatory phenotype in BCs. In contrast, BCs expanded in AECGM, displayed increased expression of structural and extracellular matrix components at high passage. Whereas culture in AECGM increased expression of some cytokines at high passage, DSI suppressed inflammation altogether thus implicating TGF-ß in BC inflammatory processes. Differentiation capacity declined with time in culture irrespective of expansion media except for PLUNC expressing secretory cells that were elevated at high passage in AECGM and PEx+ suggestive of an immune modulatory role of PLUNC in BCs. These findings underscore the profound impact of media conditions on inflammatory niche and function of in vitro expanded BCs. The broad pro-inflammatory phenotype driven by PEx+ media, in particular, should be considered in the development of cell-based models for airway diseases and therapeutic application.

PEACOCK: a machine learning approach to assess the validity of cell type-specific enhancer-gene regulatory relationships.

The vast majority of disease-associated variants identified in genome-wide association studies map to enhancers, powerful regulatory elements which orchestrate the recruitment of transcriptional complexes to their target genes' promoters to upregulate transcription in a cell type- and timing-dependent manner. These variants have implicated thousands of enhancers in many common genetic diseases, including nearly all cancers. However, the etiology of most of these diseases remains unknown because the regulatory target genes of the vast majority of enhancers are unknown. Thus, identifying the target genes of as many enhancers as possible is crucial for learning how enhancer regulatory activities function and contribute to disease. Based on experimental results curated from scientific publications coupled with machine learning methods, we developed a cell type-specific score predictive of an enhancer targeting a gene. We computed the score genome-wide for every possible cis enhancer-gene pair and validated its predictive ability in four widely used cell lines. Using a pooled final model trained across multiple cell types, all possible gene-enhancer regulatory links in cis (~17 M) were scored and added to the publicly available PEREGRINE database ( www.peregrineproj.org ). These scores provide a quantitative framework for the enhancer-gene regulatory prediction that can be incorporated into downstream statistical analyses.

GRAMD2 + alveolar type I cell plasticity facilitates cell state transitions in organoid culture.

Alveolar epithelial regeneration is critical for normal lung function and becomes dysregulated in disease. While alveolar type 2 (AT2) and club cells are known distal lung epithelial progenitors, determining if alveolar epithelial type 1 (AT1) cells also contribute to alveolar regeneration has been hampered by lack of highly specific mouse models labeling AT1 cells. To address this, the Gramd2 CreERT2 transgenic strain was generated and crossed to Rosa mTmG mice. Extensive cellular characterization, including distal lung immunofluorescence and cytospin staining, confirmed that GRAMD2 + AT1 cells are highly enriched for green fluorescent protein (GFP). Interestingly, Gramd2 CreERT2 GFP + cells were able to form organoids in organoid co-culture with Mlg fibroblasts. Temporal scRNAseq revealed that Gramd2 + AT1 cells transition through numerous intermediate lung epithelial cell states including basal, secretory and AT2 cell in organoids while acquiring proliferative capacity. Our results indicate that Gramd2 + AT1 cells are highly plastic suggesting they may contribute to alveolar regeneration.

Alveolar type I cells can give rise to KRAS-induced lung adenocarcinoma.

Lung adenocarcinoma (LUAD) is the most prevalent subtype of lung cancer and presents clinically with a high degree of biological heterogeneity and distinct clinical outcomes. The current paradigm of LUAD etiology posits alveolar epithelial type II (AT2) cells as the primary cell of origin, while the role of AT1 cells in LUAD oncogenesis remains unknown. Here, we examine oncogenic transformation in mouse Gram-domain containing 2 (Gramd2)+ AT1 cells via oncogenic KRASG12D. Activation of KRASG12D in AT1 cells induces multifocal LUAD, primarily of papillary histology. Furthermore, KRT8+ intermediate cell states were observed in both AT2- and AT1-derived LUAD, but SCGB3A2+, another intermediate cell marker, was primarily associated with AT1 cells, suggesting different mechanisms of tumor evolution. Collectively, our study reveals that Gramd2+ AT1 cells can serve as a cell of origin for LUAD and suggests that distinct subtypes of LUAD based on cell of origin be considered in the development of therapeutics.

Development of human alveolar epithelial cell models to study distal lung biology and disease.

Many acute and chronic diseases affect the distal lung alveoli. Alveolar epithelial cell (AEC) lines are needed to better model these diseases. We used de-identified human remnant transplant lungs to develop a method to establish AEC lines. The lines grow well in 2-dimensional (2D) culture as epithelial monolayers expressing lung progenitor markers. In 3-dimensional (3D) culture with fibroblasts, Matrigel, and specific media conditions, the cells form alveolar-like organoids expressing mature AEC markers including aquaporin 5 (AQP5), G-protein-coupled receptor class C group 5 member A (GPRC5A), and surface marker HTII280. Single-cell RNA sequencing of an AEC line in 2D versus 3D culture revealed increased cellular heterogeneity and induction of cytokine and lipoprotein signaling in 3D organoids. Our approach yields lung progenitor lines that retain the ability to differentiate along the alveolar cell lineage despite long-term expansion and provides a valuable system to model and study the distal lung in vitro.

Indole-3-carbinol downregulation of telomerase gene expression requires the inhibition of estrogen receptor-alpha and Sp1 transcription factor interactions within the hTERT promoter and mediates the G1 cell cycle arrest of human breast cancer cells.

Indole-3-carbinol (I3C), a naturally occurring hydrolysis product of glucobrassicin from cruciferous vegetables such as broccoli, cabbage and Brussels sprouts, is an anticancer phytochemical that triggers complementary sets of antiproliferative pathways to induce a cell cycle arrest of estrogen-responsive MCF7 breast cancer cells. I3C strongly downregulated transcript expression of the catalytic subunit of the human telomerase (hTERT) gene, which correlated with the dose-dependent indole-mediated G(1) cell cycle arrest without altering the transcript levels of the RNA template (hTR) for telomerase elongation. Exogenous expression of hTERT driven by a constitutive promoter prevented the I3C-induced cell cycle arrest and rescued the I3C inhibition of telomerase enzymatic activity and activation of cellular senescence. Time course studies showed that I3C downregulated expression of estrogen receptor-alpha (ERα) and cyclin-dependent kinase-6 transcripts levels (which is regulated through the Sp1 transcription factor) prior to the downregulation of hTERT suggesting a mechanistic link. Chromatin immunoprecipitation assays demonstrated that I3C disrupted endogenous interactions of both ERα and Sp1 with an estrogen response element-Sp1 composite element within the hTERT promoter. I3C inhibited 17ß-estradiol stimulated hTERT expression and stimulated the production of threonine-phosphorylated Sp1, which inhibits Sp1-DNA interactions. Exogenous expression of both ERα and Sp1, but not either alone, in MCF7 cells blocked the I3C-mediated downregulation of hTERT expression. These results demonstrate that I3C disrupts the combined ERα- and Sp1-driven transcription of hTERT gene expression, which plays a significant role in the I3C-induced cell cycle arrest of human breast cancer cells.

The Sulfotransferase SULT1C2 Is Epigenetically Activated and Transcriptionally Induced by Tobacco Exposure and Is Associated with Patient Outcome in Lung Adenocarcinoma.

Lung cancer is the leading cause of cancer-related death. Tobacco exposure is associated with 80-90% of lung cancer cases. The SULT1C2 sulfotransferase modifies xenobiotic compounds to enhance secretion but can also render these compounds carcinogenic. To determine if SULT1C2 contributes to tobacco-related carcinogenesis in the lung, we analyzed the expression and epigenetic state of SULT1C2 in human lung adenocarcinoma (LUAD) samples and in LUAD cell lines exposed to cigarette smoke condensate (CSC). SULT1C2 expression was significantly positively correlated to overall LUAD patient survival in smokers, was elevated in LUAD tumors compared to adjacent non-tumor lung, and was significantly correlated with levels of patient exposure to tobacco smoke. SULT1C2 promoter DNA methylation was inversely correlated with expression in LUAD, and hypomethylation of the SULT1C2 promoter was observed in Asian patients, as compared to Caucasians. In vitro analysis of LUAD cell lines indicates that CSC stimulates expression of SULT1C2 in a dose-dependent and cell-line-specific manner. In vitro methylation of the SULT1C2 promoter significantly decreased transcriptional activity of a reporter plasmid, and SULT1C2 expression was activated by the DNA demethylating agent 5-Aza-2'-deoxycytidine in a cell line in which the SULT1C2 promoter was hypermethylated. An aryl hydrocarbon receptor (AHR) binding site was detected spanning critical methylation sites upstream of SULT1C2. CSC exposure significantly increased AHR binding to this predicted binding site in the SULT1C2 promoter in multiple lung cell lines. Our data suggest that CSC exposure leads to activation of the AHR transcription factor, increased binding to the SULT1C2 promoter, and upregulation of SULT1C2 expression and that this process is inhibited by DNA methylation at the SULT1C2 locus. Additionally, our results suggest that the level of SULT1C2 promoter methylation and gene expression in normal lung varies depending on the race of the patient, which could in part reflect the molecular mechanisms of racial disparities seen in lung cellular responses to cigarette smoke exposure.

PEREGRINE: A genome-wide prediction of enhancer to gene relationships supported by experimental evidence.

Enhancers are powerful and versatile agents of cell-type specific gene regulation, which are thought to play key roles in human disease. Enhancers are short DNA elements that function primarily as clusters of transcription factor binding sites that are spatially coordinated to regulate expression of one or more specific target genes. These regulatory connections between enhancers and target genes can therefore be characterized as enhancer-gene links that can affect development, disease, and homeostatic cellular processes. Despite their implication in disease and the establishment of cell identity during development, most enhancer-gene links remain unknown. Here we introduce a new, publicly accessible database of predicted enhancer-gene links, PEREGRINE. The PEREGRINE human enhancer-gene links interactive web interface incorporates publicly available experimental data from ChIA-PET, eQTL, and Hi-C assays across 78 cell and tissue types to link 449,627 enhancers to 17,643 protein-coding genes. These enhancer-gene links are made available through the new Enhancer module of the PANTHER database and website where the user may easily access the evidence for each enhancer-gene link, as well as query by target gene and enhancer location.

Integrated Single-Cell RNA-Sequencing Analysis of Aquaporin 5-Expressing Mouse Lung Epithelial Cells Identifies GPRC5A as a Novel Validated Type I Cell Surface Marker.

Molecular and functional characterization of alveolar epithelial type I (AT1) cells has been challenging due to difficulty in isolating sufficient numbers of viable cells. Here we performed single-cell RNA-sequencing (scRNA-seq) of tdTomato+ cells from lungs of AT1 cell-specific Aqp5-Cre-IRES-DsRed (ACID);R26tdTomato reporter mice. Following enzymatic digestion, CD31-CD45-E-cadherin+tdTomato+ cells were subjected to fluorescence-activated cell sorting (FACS) followed by scRNA-seq. Cell identity was confirmed by immunofluorescence using cell type-specific antibodies. After quality control, 92 cells were analyzed. Most cells expressed 'conventional' AT1 cell markers (Aqp5, Pdpn, Hopx, Ager), with heterogeneous expression within this population. The remaining cells expressed AT2, club, basal or ciliated cell markers. Integration with public datasets identified three robust AT1 cell- and lung-enriched genes, Ager, Rtkn2 and Gprc5a, that were conserved across species. GPRC5A co-localized with HOPX and was not expressed in AT2 or airway cells in mouse, rat and human lung. GPRC5A co-localized with AQP5 but not pro-SPC or CC10 in mouse lung epithelial cell cytospins. We enriched mouse AT1 cells to perform molecular phenotyping using scRNA-seq. Further characterization of putative AT1 cell-enriched genes revealed GPRC5A as a conserved AT1 cell surface marker that may be useful for AT1 cell isolation.