Positional integration of lung adenocarcinoma susceptibility loci with primary human alveolar epithelial cell epigenomes

Lack of racial and ethnic diversity in lung cancer cell lines contributes to lung cancer health disparities

Lung cancer is the leading cause of cancer death in the United States and worldwide, and a major source of cancer health disparities. Lung cancer cell lines provide key in vitro models for molecular studies of lung cancer development and progression, and for pre-clinical drug testing. To ensure health equity, it is imperative that cell lines representing different lung cancer histological types, carrying different cancer driver genes, and representing different genders, races, and ethnicities should be available. This is particularly relevant for cell lines from Black men, who experience the highest lung cancer mortality in the United States. Here, we undertook a review of the available lung cancer cell lines and their racial and ethnic origin. We noted a marked imbalance in the availability of cell lines from different races and ethnicities. Cell lines from Black patients were strongly underrepresented, and we identified no cell lines from Hispanic/Latin(x) (H/L), American Indian/American Native (AI/AN), or Native Hawaiian or other Pacific Islander (NHOPI) patients. The majority of cell lines were derived from White and Asian patients. Also missing are cell lines representing the cells-of-origin of the major lung cancer histological types, which can be used to model lung cancer development and to study the effects of environmental exposures on lung tissues. To our knowledge, the few available immortalized alveolar epithelial cell lines are all derived from White subjects, and the race and ethnicity of a handful of cell lines derived from bronchial epithelial cells are unknown. The lack of an appropriately diverse collection of lung cancer cell lines and lung cancer cell-of-origin lines severely limits racially and ethnically inclusive lung cancer research. It impedes the ability to develop inclusive models, screen comprehensively for effective compounds, pre-clinically test new drugs, and optimize precision medicine. It thereby hinders the development of therapies that can increase the survival of minority and underserved patients. The noted lack of cell lines from underrepresented groups should constitute a call to action to establish additional cell lines and ensure adequate representation of all population groups in this critical pre-clinical research resource.

Functional characterization of 5p15.33 risk locus in uveal melanoma reveals rs452384 as a functional variant and NKX2.4 as an allele-specific interactor

The TERT/CLPTM1L risk locus on chromosome 5p15.33 is a pleiotropic cancer risk locus in which multiple independent risk alleles have been identified, across well over ten cancer types. We previously conducted a genome-wide association study in uveal melanoma (UM), which uncovered a role for the TERT/CLPTM1L risk locus in this intraocular tumor and identified multiple highly correlated risk alleles. Aiming to unravel the biological mechanisms in UM of this locus, which contains a domain enriched in active chromatin marks and enhancer elements, we demonstrated the allele-specific enhancer activity of this risk region using reporter assays. In UM, we identified the functional variant rs452384, of which the C risk allele is associated with higher gene expression, increased CLPTM1L expression in UM tumors, and a longer telomere length in peripheral blood mononuclear cells. Electrophoretic mobility shift assays and quantitative mass spectrometry identified NKX2.4 as an rs452384-T-specific binding protein, whereas GATA4 preferentially interacted with rs452384-C. Knockdown of NKX2.4 but not GATA4 resulted in increased TERT and CLPTM1L expression. In summary, the UM risk conferred by the 5p locus is at least partly due to rs452384, for which NKX2.4 presents strong differential binding activity and regulates CLPTM1L and TERT expression. Altogether, our work unraveled some of the complex regulatory mechanisms at the 5p15.33 susceptibility region in UM, and this might also shed light on shared mechanisms with other tumor types affected by this susceptibility region.

The landscape of GWAS validation; systematic review identifying 309 validated non-coding variants across 130 human diseases

Background

The remarkable growth of genome-wide association studies (GWAS) has created a critical need to experimentally validate the disease-associated variants, 90% of which involve non-coding variants.

Methods

To determine how the field is addressing this urgent need, we performed a comprehensive literature review identifying 36,676 articles. These were reduced to 1454 articles through a set of filters using natural language processing and ontology-based text-mining. This was followed by manual curation and cross-referencing against the GWAS catalog, yielding a final set of 286 articles.

Results

We identified 309 experimentally validated non-coding GWAS variants, regulating 252 genes across 130 human disease traits. These variants covered a variety of regulatory mechanisms. Interestingly, 70% (215/309) acted through cis-regulatory elements, with the remaining through promoters (22%, 70/309) or non-coding RNAs (8%, 24/309). Several validation approaches were utilized in these studies, including gene expression (n = 272), transcription factor binding (n = 175), reporter assays (n = 171), in vivo models (n = 104), genome editing (n = 96) and chromatin interaction (n = 33).

Conclusions

This review of the literature is the first to systematically evaluate the status and the landscape of experimentation being used to validate non-coding GWAS-identified variants. Our results clearly underscore the multifaceted approach needed for experimental validation, have practical implications on variant prioritization and considerations of target gene nomination. While the field has a long way to go to validate the thousands of GWAS associations, we show that progress is being made and provide exemplars of validation studies covering a wide variety of mechanisms, target genes, and disease areas.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12920-022-01216-w.

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.

MALAT1 Polymorphisms and Lung Cancer Susceptibility in a Chinese Northeast Han Population

Background: LncRNA MALAT1 (metastasis-associated lung adenocarcinoma transcript 1) was competitive endogenous RNA (ceRNA) involved in various molecular processes for metastasis development in lung cancer. Single nucleotide polymorphisms (SNPs) in MALAT1 gene might be predictive markers for lung cancer. In our study, we selected rs619586 and rs3200401 in MALAT1 gene to explore their effects on lung cancer susceptibility. Methods: The case-control study included 444 lung cancer cases and 460 healthy controls. Genotyping was performed by Taqman allelic discrimination method. Logistic regression, Student t-test, and Chi-square test (χ² ) were used to analyze the data. Results: The findings of the study showed that rs3200401 was significantly associated with the risk of non-small cell lung cancer (NSCLC) and lung squamous cell carcinoma (LUSC). Compared with homozygous CC genotype, CT heterozygous genotype decreased risk of NSCLC (P^a = 0.034) and LUSC (P^a = 0.025). In addition, no statistical association was detected between rs619586 and lung cancer susceptibility. The interactions between genes and cigarette smoking were discovered via crossover analysis. However, there were no remarkable gene-environment interactions in additive and multiplicative model. Conclusion: Rs3200401 in lncRNA MALAT1 was associated with the susceptibility of non-small-cell lung cancer and lung squamous cell carcinoma. The gene-environmental (cigarette smoking) interactions were not notable.

Comprehensive epigenomic profiling of human alveolar epithelial differentiation identifies key epigenetic states and transcription factor co-regulatory networks for maintenance of distal lung identity

BACKGROUND

Disruption of alveolar epithelial cell (AEC) differentiation is implicated in distal lung diseases such as chronic obstructive pulmonary disease, idiopathic pulmonary fibrosis, and lung adenocarcinoma that impact morbidity and mortality worldwide. Elucidating underlying disease pathogenesis requires a mechanistic molecular understanding of AEC differentiation. Previous studies have focused on changes of individual transcription factors, and to date no study has comprehensively characterized the dynamic, global epigenomic alterations that facilitate this critical differentiation process in humans.

RESULTS

We comprehensively profiled the epigenomic states of human AECs during type 2 to type 1-like cell differentiation, including the methylome and chromatin functional domains, and integrated this with transcriptome-wide RNA expression data. Enhancer regions were drastically altered during AEC differentiation. Transcription factor binding analysis within enhancer regions revealed diverse interactive networks with enrichment for many transcription factors, including NKX2-1 and FOXA family members, as well as transcription factors with less well characterized roles in AEC differentiation, such as members of the MEF2, TEAD, and AP1 families. Additionally, associations among transcription factors changed during differentiation, implicating a complex network of heterotrimeric complex switching in driving differentiation. Integration of AEC enhancer states with the catalog of enhancer elements in the Roadmap Epigenomics Mapping Consortium and Encyclopedia of DNA Elements (ENCODE) revealed that AECs have similar epigenomic structures to other profiled epithelial cell types, including human mammary epithelial cells (HMECs), with NKX2-1 serving as a distinguishing feature of distal lung differentiation.

CONCLUSIONS

Enhancer regions are hotspots of epigenomic alteration that regulate AEC differentiation. Furthermore, the differentiation process is regulated by dynamic networks of transcription factors acting in concert, rather than individually. These findings provide a roadmap for understanding the relationship between disruption of the epigenetic state during AEC differentiation and development of lung diseases that may be therapeutically amenable.

DCBLD1 is associated with the integrin signaling pathway and has prognostic value in non-small cell lung and invasive breast carcinoma

Germline single nucleotide polymorphisms in the promoter region of the DCBLD1 gene are associated with non-smoking cases of both non-small cell lung carcinoma (NSCLC) and human papillomavirus-negative head and neck cancer. However the clinical relevance and function of DCBLD1 remain unclear. This multicenter retrospective study was designed to evaluate the prognostic value and function of DCBLD1 in the four main solid cancers: NSCLC, invasive breast carcinoma, colorectal adenocarcinoma and prostate adenocarcinoma. We included the following cohorts: GSE81089 NSCLC, METABRIC invasive breast carcinoma, GSE14333 colorectal adenocarcinoma, GSE70770 prostate adenocarcinoma and The Cancer Genome Atlas (TCGA) Firehose Legacy cohorts of all four cancers. DCBLD1 gene expression was associated with a worse overall survival in multivariate analyses for both NSCLC cohorts (TCGA: P = 0.03 and GSE81089: P = 0.04) and both invasive breast carcinoma cohorts (TCGA: P = 0.02 and METABRIC: P < 0.001). Patients with high DCBLD1 expression showed an upregulation of the integrin signaling pathway in comparison to those with low DCBLD1 expression in the TCGA NSCLC cohort (FDR = 5.16 × 10^-14) and TCGA invasive breast carcinoma cohort (FDR = 1.94 × 10^-05).

Functional annotation of lung cancer‒associated genetic variants by cell type‒specific epigenome and long-range chromatin interactome

Functional interpretation of noncoding genetic variants associated with complex human diseases and traits remains a challenge. In an effort to enhance our understanding of common germline variants associated with lung cancer, we categorize regulatory elements based on eight major cell types of human lung tissue. Our results show that 21.68% of lung cancer‒associated risk variants are linked to noncoding regulatory elements, nearly half of which are cell type‒specific. Integrative analysis of high-resolution long-range chromatin interactome maps and single-cell RNA-sequencing data of lung tumors uncovers number of putative target genes of these variants and functionally relevant cell types, which display a potential biological link to cancer susceptibility. The present study greatly expands the scope of functional annotation of lung cancer‒associated genetic risk factors and dictates probable cell types involved in lung carcinogenesis.

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.

Single Nucleotide Polymorphism rs6942067 Is a Risk Factor in Young and in Non-Smoking Patients with HPV Negative Head and Neck Squamous Cell Carcinoma

Genetic factors behind the increasing incidence of human papillomavirus (HPV) negative head and neck squamous cell carcinoma (HNSCC) in young non-smokers are suspected, but have not been identified. Recently, rs6942067, a single nucleotide polymorphism (SNP) located upstream of the DCBLD1 gene, was found associated with non-smoking lung adenocarcinoma. To validate if this SNP is also implicated in HNSCC, participants of The Cancer Genome Atlas HNSCC cohort were investigated for rs6942067 status, associated DCBLD1 expression, and clinical characteristics. Occurrence of the rs6942067 GG genotype is significantly higher in young and in HPV negative non-smoking HNSCC than in other HNSCC. Additionally, rs6942067 GG is associated with higher DCBLD1 expression in HNSCC and patients with high DCBLD1 expression have a worse overall survival at three years, both in univariate and multivariate analysis. Furthermore, high DCBLD1 expression is associated with activation of the integrin signaling pathway and its phosphorylation with EGFR and MET. Collectively, these findings suggest that DCBLD1 plays a critical role in HNSCC and demonstrate an association between rs6942067 and clinical characteristics of young age and HPV negative non-smoking status in HNSCC patients.

The DCBLD receptor family: emerging signaling roles in development, homeostasis and disease

The discoidin, CUB, and LCCL domain-containing (DCBLD) receptor family are composed of the type-I transmembrane proteins DCBLD1 and DCBLD2 (also ESDN and CLCP1). These proteins are highly conserved across vertebrates and possess similar domain structure to that of neuropilins, which act as critical co-receptors in developmental processes. Although DCBLD1 remains largely uncharacterized, the functional and mechanistic roles of DCBLD2 are emerging. This review provides a comprehensive discussion of this presumed receptor family, ranging from structural and signaling aspects to their associations with cancer, physiology, and development.