Epigenetic Regulation of DLK1-DIO3 Region in Thyroid Carcinoma.

Non-coding RNAs (ncRNAs) have emerged as pivotal regulators in cellular biology, dispelling their former perception as 'junk transcripts'. Notably, the DLK1-DIO3 region harbors numerous ncRNAs, including long non-coding RNAs (lncRNAs) and over 50 microRNA genes. While papillary thyroid cancer showcases a pervasive decrease in DLK1-DIO3-derived ncRNA expression, the precise mechanisms driving this alteration remain elusive. We hypothesized that epigenetic alterations underlie shifts in ncRNA expression during thyroid cancer initiation and progression. This study aimed to elucidate the epigenetic mechanisms governing DLK1-DIO3 region expression in this malignancy. We have combined the analysis of DNA methylation by bisulfite sequencing together with that of histone modifications through ChIP-qPCR to gain insights into the epigenetic contribution to thyroid cancer in cell lines representing malignancies with different genetic backgrounds. Our findings characterize the region's epigenetic signature in thyroid cancer, uncovering distinctive DNA methylation patterns, particularly within CpG islands on the lncRNA MEG3-DMR, which potentially account for its downregulation in tumors. Pharmacological intervention targeting DNA methylation combined with histone deacetylation restored ncRNA expression. These results contribute to the understanding of the epigenetic mechanisms controlling the DLK1-DIO3 region in thyroid cancer, highlighting the combined role of DNA methylation and histone marks in regulating the locus' expression.

Cancer cell plasticity defines response to immunotherapy in cutaneous squamous cell carcinoma.

Immune checkpoint blockade (ICB) approaches have changed the therapeutic landscape for many tumor types. However, half of cutaneous squamous cell carcinoma (cSCC) patients remain unresponsive or develop resistance. Here, we show that, during cSCC progression in male mice, cancer cells acquire epithelial/mesenchymal plasticity and change their immune checkpoint (IC) ligand profile according to their features, dictating the IC pathways involved in immune evasion. Epithelial cancer cells, through the PD-1/PD-L1 pathway, and mesenchymal cancer cells, through the CTLA-4/CD80 and TIGIT/CD155 pathways, differentially block antitumor immune responses and determine the response to ICB therapies. Accordingly, the anti-PD-L1/TIGIT combination is the most effective strategy for blocking the growth of cSCCs that contain both epithelial and mesenchymal cancer cells. The expression of E-cadherin/Vimentin/CD80/CD155 proteins in cSCC, HNSCC and melanoma patient samples predicts response to anti-PD-1/PD-L1 therapy. Collectively, our findings indicate that the selection of ICB therapies should take into account the epithelial/mesenchymal features of cancer cells.

IκBα controls dormancy in hematopoietic stem cells via retinoic acid during embryonic development.

Recent findings suggest that Hematopoietic Stem Cells (HSC) and progenitors arise simultaneously and independently of each other already in the embryonic aorta-gonad mesonephros region, but it is still unknown how their different features are established. Here, we uncover IκBα (Nfkbia, the inhibitor of NF-κB) as a critical regulator of HSC proliferation throughout development. IκBα balances retinoic acid signaling levels together with the epigenetic silencer, PRC2, specifically in HSCs. Loss of IκBα decreases proliferation of HSC and induces a dormancy related gene expression signature instead. Also, IκBα deficient HSCs respond with superior activation to in vitro culture and in serial transplantation. At the molecular level, chromatin regions harboring binding motifs for retinoic acid signaling are hypo-methylated for the PRC2 dependent H3K27me3 mark in IκBα deficient HSCs. Overall, we show that the proliferation index in the developing HSCs is regulated by a IκBα-PRC2 axis, which controls retinoic acid signaling.

DNA methylation profiling of myelodysplastic syndromes and clinical response to azacitidine: A multicentre retrospective study.

Real-world data have revealed that a substantial portion of patients with myelodysplastic syndromes (MDS) does not respond to epigenetic therapy with hypomethylating agents (HMAs). The cellular and molecular reasons for this resistance to the demethylating agent and biomarkers that would be able to predict the treatment refractoriness are largely unknown. In this study, we shed light on this enigma by characterizing the epigenomic profiles of patients with MDS treated with azacitidine. Our approach provides a comprehensive view of the evolving DNA methylation architecture of the disease and holds great potential for advancing our understanding of MDS treatment responses to HMAs.

Single-cell Multiomics Analysis of Myelodysplastic Syndromes and Clinical Response to Hypomethylating Therapy.

Alterations in epigenetic marks, such as DNA methylation, represent a hallmark of cancer that has been successfully exploited for therapy in myeloid malignancies. Hypomethylating agents (HMA), such as azacitidine, have become standard-of-care therapy to treat myelodysplastic syndromes (MDS), myeloid neoplasms that can evolve into acute myeloid leukemia. However, our capacity to identify who will respond to HMAs, and the duration of response, remains limited. To shed light on this question, we have leveraged the unprecedented analytic power of single-cell technologies to simultaneously map the genome and immunoproteome of MDS samples throughout clinical evolution. We were able to chart the architecture and evolution of molecular clones in precious paired bone marrow MDS samples at diagnosis and posttreatment to show that a combined imbalance of specific cell lineages with diverse mutational profiles is associated with the clinical response of patients with MDS to hypomethylating therapy. SIGNIFICANCE: MDS are myeloid clonal hemopathies with a low 5-year survival rate, and approximately half of the cases do not respond to standard HMA therapy. Our innovative single-cell multiomics approach offers valuable biological insights and potential biomarkers associated with the demethylating agent efficacy. It also identifies vulnerabilities that can be targeted using personalized combinations of small drugs and antibodies.

Cis inhibition of NOTCH1 through JAGGED1 sustains embryonic hematopoietic stem cell fate.

Hematopoietic stem cells (HSCs) develop from the hemogenic endothelium (HE) in the aorta- gonads-and mesonephros (AGM) region and reside within Intra-aortic hematopoietic clusters (IAHC) along with hematopoietic progenitors (HPC). The signalling mechanisms that distinguish HSCs from HPCs are unknown. Notch signaling is essential for arterial specification, IAHC formation and HSC activity, but current studies on how Notch segregates these different fates are inconsistent. We now demonstrate that Notch activity is highest in a subset of, GFI1 + , HSC-primed HE cells, and is gradually lost with HSC maturation. We uncover that the HSC phenotype is maintained due to increasing levels of NOTCH1 and JAG1 interactions on the surface of the same cell (cis) that renders the NOTCH1 receptor from being activated. Forced activation of the NOTCH1 receptor in IAHC activates a hematopoietic differentiation program. Our results indicate that NOTCH1-JAG1 cis-inhibition preserves the HSC phenotype in the hematopoietic clusters of the embryonic aorta.

ERBB4-Mediated Signaling Is a Mediator of Resistance to PI3K and BTK Inhibitors in B-cell Lymphoid Neoplasms.

BTK and PI3K inhibitors are among the drugs approved for the treatment of patients with lymphoid neoplasms. Although active, their ability to lead to long-lasting complete remission is rather limited, especially in the lymphoma setting. This indicates that tumor cells often develop resistance to the drugs. We started from a marginal zone lymphoma cell line, Karpas-1718, kept under prolonged exposure to the PI3Kδ inhibitor idelalisib until acquisition of resistance, or with no drug. Cells underwent transcriptome, miRNA and methylation profiling, whole-exome sequencing, and pharmacologic screening, which led to the identification of the overexpression of ERBB4 and its ligands HBEGF and NRG2 in the resistant cells. Cellular and genetic experiments demonstrated the involvement of this axis in blocking the antitumor activity of various BTK/PI3K inhibitors, currently used in the clinical setting. Addition of recombinant HBEGF induced resistance to BTK/PI3K inhibitors in parental cells and in additional lymphoma models. Combination with the ERBB inhibitor lapatinib was beneficial in resistant cells and in other lymphoma models already expressing the identified resistance factors. An epigenetic reprogramming sustained the expression of the resistance-related factors, and pretreatment with demethylating agents or EZH2 inhibitors overcame the resistance. Resistance factors were also shown to be expressed in clinical specimens. In conclusion, we showed that the overexpression of ERBB4 and its ligands represents a novel mechanism of resistance for lymphoma cells to bypass the antitumor activity of BTK and PI3K inhibitors and that targeted pharmacologic interventions can restore sensitivity to the small molecules.

Gastric metaplasia as a precursor of nonconventional dysplasia in inflammatory bowel disease.

Gastric metaplasia in colonic mucosa with inflammatory bowel disease (IBD) develops as an adaptation mechanism. The association between gastric metaplasia and nonconventional and/or conventional dysplasia as precursors of colitis-associated colorectal cancer is unknown. To address this question, we retrospectively reviewed a series of 33 IBD colectomies to identify gastric metaplasia in 76 precursor lesions. We obtained 61 nonconventional and 15 conventional dysplasias. Among nonconventional dysplasia, 31 (50.8 %) were low-grade (LGD), 4 (6.5 %) were high-grade (HGD), 9 (14.8 %) had both LGD and HGD, and 17 (27.9 %) had no dysplasia (ND), while 14 (93 %) conventional dysplasias had LGD, and 1 (7 %) had LGD and HGD. Gastric metaplasia was assessed by concomitant immunoexpression of MUC5AC and loss of CDX2 staining. Expression of a p53-mut pattern was considered as a surrogate for gene mutation, and complete loss of MLH1 staining as presence of MLH1 hypermethylation. In nonconventional dysplasia, MUC5AC immunoexpression decreased as the degree of dysplasia increased, being 78 % in LGD and 39 % in HGD (p = 0.006). CDX2 was lost in epithelial glands with high expression of MUC5AC (p < 0.001). The p53-mut pattern was observed in 77 % HGD, 45 % LGD, and in 6 % with ND (p < 0.001). Neither nonconventional nor conventional dysplasia showed complete loss of MLH1 staining. Gastric metaplasia was also present in mucosa adjacent to nonconventional dysplasia with chronic changes or active inflammation. Our results show that gastric metaplasia appears in IBD-inflamed colon mucosa, it is the substrate of most nonconventional dysplasia and occurs prior to p53 alterations.

Serum DNA methylome of the colorectal cancer serrated pathway enables non-invasive detection.

The clinical relevance of the colorectal cancer serrated pathway is evident, but the screening of serrated lesions remains challenging. We aimed to characterize the serum methylome of the serrated pathway and to evaluate circulating cell-free DNA (cfDNA) methylomes as a potential source of biomarkers for the non-invasive detection of serrated lesions. We collected serum samples from individuals with serrated adenocarcinoma (SAC), traditional serrated adenomas, sessile serrated lesions, hyperplastic polyps and individuals with no colorectal findings. First, we quantified cfDNA methylation with the MethylationEPIC array. Then, we compared the methylation profiles with tissue and serum datasets. Finally, we evaluated the utility of serum cfDNA methylation biomarkers. We identified a differential methylation profile able to distinguish high-risk serrated lesions from no serrated neoplasia, showing concordance with tissue methylation from SAC and sessile serrated lesions. Serum methylation profiles are pathway-specific, clearly separating serrated lesions from conventional adenomas. The combination of ninjurin 2 (NINJ2) and glutamate-rich 1 (ERICH1) methylation discriminated high-risk serrated lesions and SAC with 91.4% sensitivity (64.4% specificity), while zinc finger protein 718 (ZNF718) methylation reported 100% sensitivity for the detection of SAC (96% specificity). This is the first study exploring the serum methylome of serrated lesions. Differential methylation of cfDNA can be used for the non-invasive detection of colorectal serrated lesions.

Chromatin insulation orchestrates matrix metalloproteinase gene cluster expression reprogramming in aggressive breast cancer tumors.

BACKGROUND: Triple-negative breast cancer (TNBC) is an aggressive subtype that exhibits a high incidence of distant metastases and lacks targeted therapeutic options. Here we explored how the epigenome contributes to matrix metalloprotease (MMP) dysregulation impacting tumor invasion, which is the first step of the metastatic process. METHODS: We combined RNA expression and chromatin interaction data to identify insulator elements potentially associated with MMP gene expression and invasion. We employed CRISPR/Cas9 to disrupt the CCCTC-Binding Factor (CTCF) binding site on an insulator element downstream of the MMP8 gene (IE8) in two TNBC cellular models. We characterized these models by combining Hi-C, ATAC-seq, and RNA-seq with functional experiments to determine invasive ability. The potential of our findings to predict the progression of ductal carcinoma in situ (DCIS), was tested in data from clinical specimens. RESULTS: We explored the clinical relevance of an insulator element located within the Chr11q22.2 locus, downstream of the MMP8 gene (IE8). This regulatory element resulted in a topologically associating domain (TAD) boundary that isolated nine MMP genes into two anti-correlated expression clusters. This expression pattern was associated with worse relapse-free (HR = 1.57 [1.06 - 2.33]; p = 0.023) and overall (HR = 2.65 [1.31 - 5.37], p = 0.005) survival of TNBC patients. After CRISPR/Cas9-mediated disruption of IE8, cancer cells showed a switch in the MMP expression signature, specifically downregulating the pro-invasive MMP1 gene and upregulating the antitumorigenic MMP8 gene, resulting in reduced invasive ability and collagen degradation. We observed that the MMP expression pattern predicts DCIS that eventually progresses into invasive ductal carcinomas (AUC = 0.77, p < 0.01). CONCLUSION: Our study demonstrates how the activation of an IE near the MMP8 gene determines the regional transcriptional regulation of MMP genes with opposing functional activity, ultimately influencing the invasive properties of aggressive forms of breast cancer.

3-D chromatin conformation, accessibility, and gene expression profiling of triple-negative breast cancer.

OBJECTIVES: Triple-negative breast cancer (TNBC) is a highly aggressive breast cancer subtype with limited treatment options. Unlike other breast cancer subtypes, the scarcity of specific therapies and greater frequencies of distant metastases contribute to its aggressiveness. We aimed to find epigenetic changes that aid in the understanding of the dissemination process of these cancers. DATA DESCRIPTION: Using CRISPR/Cas9, our experimental approach led us to identify and disrupt an insulator element, IE8, whose activity seemed relevant for cell invasion. The experiments were performed in two well-established TNBC cellular models, the MDA-MB-231 and the MDA-MB-436. To gain insights into the underlying molecular mechanisms of TNBC invasion ability, we generated and characterized high-resolution chromatin interaction (Hi-C) and chromatin accessibility (ATAC-seq) maps in both cell models and complemented these datasets with gene expression profiling (RNA-seq) in MDA-MB-231, the cell line that showed more significant changes in chromatin accessibility. Altogether, our data provide a comprehensive resource for understanding the spatial organization of the genome in TNBC cells, which may contribute to accelerating the discovery of TNBC-specific alterations triggering advances for this devastating disease.

Formaldehyde regulates S-adenosylmethionine biosynthesis and one-carbon metabolism.

One-carbon metabolism is an essential branch of cellular metabolism that intersects with epigenetic regulation. In this work, we show how formaldehyde (FA), a one-carbon unit derived from both endogenous sources and environmental exposure, regulates one-carbon metabolism by inhibiting the biosynthesis of S-adenosylmethionine (SAM), the major methyl donor in cells. FA reacts with privileged, hyperreactive cysteine sites in the proteome, including Cys120 in S-adenosylmethionine synthase isoform type-1 (MAT1A). FA exposure inhibited MAT1A activity and decreased SAM production with MAT-isoform specificity. A genetic mouse model of chronic FA overload showed a decrease n SAM and in methylation on selected histones and genes. Epigenetic and transcriptional regulation of Mat1a and related genes function as compensatory mechanisms for FA-dependent SAM depletion, revealing a biochemical feedback cycle between FA and SAM one-carbon units.

Serum methylation of GALNT9, UPF3A, WARS, and LDB2 as noninvasive biomarkers for the early detection of colorectal cancer and advanced adenomas.

BACKGROUND: Early detection has proven to be the most effective strategy to reduce the incidence and mortality of colorectal cancer (CRC). Nevertheless, most current screening programs suffer from low participation rates. A blood test may improve both the adherence to screening and the selection to colonoscopy. In this study, we conducted a serum-based discovery and validation of cfDNA methylation biomarkers for CRC screening in a multicenter cohort of 433 serum samples including healthy controls, benign pathologies, advanced adenomas (AA), and CRC. RESULTS: First, we performed an epigenome-wide methylation analysis with the MethylationEPIC array using a sample pooling approach, followed by a robust prioritization of candidate biomarkers for the detection of advanced neoplasia (AN: AA and CRC). Then, candidate biomarkers were validated by pyrosequencing in independent individual cfDNA samples. We report GALNT9, UPF3A, WARS, and LDB2 as new noninvasive biomarkers for the early detection of AN. The combination of GALNT9/UPF3A by logistic regression discriminated AN with 78.8% sensitivity and 100% specificity, outperforming the commonly used fecal immunochemical test and the methylated SEPT9 blood test. CONCLUSIONS: Overall, this study highlights the utility of cfDNA methylation for CRC screening. Our results suggest that the combination methylated GALNT9/UPF3A has the potential to serve as a highly specific and sensitive blood-based test for screening and early detection of CRC.

Culture expansion of CAR T cells results in aberrant DNA methylation that is associated with adverse clinical outcome.

Chimeric antigen receptor (CAR) T cells provide new perspectives for treatment of hematological malignancies. Manufacturing of these cellular products includes culture expansion procedures, which may affect cellular integrity and therapeutic outcome. In this study, we investigated culture-associated epigenetic changes in CAR T cells and found continuous gain of DNAm, particularly within genes that are relevant for T cell function. Hypermethylation in many genes, such as TCF7, RUNX1, and TOX, was reflected by transcriptional downregulation. 332 CG dinucleotides (CpGs) showed an almost linear gain in methylation with cell culture time, albeit neighboring CpGs were not coherently regulated on the same DNA strands. An epigenetic signature based on 14 of these culture-associated CpGs predicted cell culture time across various culture conditions. Notably, even in CAR T cell products of similar culture time higher DNAm levels at these CpGs were associated with significantly reduced long-term survival post transfusion. Our data demonstrate that cell culture expansion of CAR T cells evokes DNA hypermethylation at specific sites in the genome and the signature may also reflect loss of potential in CAR T cell products. Hence, reduced cultivation periods are beneficial to avoid dysfunctional methylation programs that seem to be associated with worse therapeutic outcome.

Cancer metastasis under the magnifying glass of epigenetics and epitranscriptomics.

Most of the cancer-associated mortality and morbidity can be attributed to metastasis. The role of epigenetic and epitranscriptomic alterations in cancer origin and progression has been extensively demonstrated during the last years. Both regulations share similar mechanisms driven by DNA or RNA modifiers, namely writers, readers, and erasers; enzymes responsible of respectively introducing, recognizing, or removing the epigenetic or epitranscriptomic modifications. Epigenetic regulation is achieved by DNA methylation, histone modifications, non-coding RNAs, chromatin accessibility, and enhancer reprogramming. In parallel, regulation at RNA level, named epitranscriptomic, is driven by a wide diversity of chemical modifications in mostly all RNA molecules. These two-layer regulatory mechanisms are finely controlled in normal tissue, and dysregulations are associated with every hallmark of human cancer. In this review, we provide an overview of the current state of knowledge regarding epigenetic and epitranscriptomic alterations governing tumor metastasis, and compare pathways regulated at DNA or RNA levels to shed light on a possible epi-crosstalk in cancer metastasis. A deeper understanding on these mechanisms could have important clinical implications for the prevention of advanced malignancies and the management of the disseminated diseases. Additionally, as these epi-alterations can potentially be reversed by small molecules or inhibitors against epi-modifiers, novel therapeutic alternatives could be envisioned.

G protein-coupled receptor 183 mediates the sensitization of Burkitt lymphoma tumors to CD47 immune checkpoint blockade by anti-CD20/PI3Kδi dual therapy.

Background: Immunotherapy-based regimens have considerably improved the survival rate of B-cell non-Hodgkin lymphoma (B-NHL) patients in the last decades; however, most disease subtypes remain almost incurable. TG-1801, a bispecific antibody that targets CD47 selectively on CD19+ B-cells, is under clinical evaluation in relapsed/refractory (R/R) B-NHL patients either as a single-agent or in combination with ublituximab, a new generation CD20 antibody. Methods: A set of eight B-NHL cell lines and primary samples were cultured in vitro in the presence of bone marrow-derived stromal cells, M2-polarized primary macrophages, and primary circulating PBMCs as a source of effector cells. Cell response to TG-1801 alone or combined with the U2 regimen associating ublituximab to the PI3Kδ inhibitor umbralisib, was analyzed by proliferation assay, western blot, transcriptomic analysis (qPCR array and RNA sequencing followed by gene set enrichment analysis) and/or quantification of antibody-dependent cell death (ADCC) and antibody-dependent cell phagocytosis (ADCP). CRISPR-Cas9 gene edition was used to selectively abrogate GPR183 gene expression in B-NHL cells. In vivo, drug efficacy was determined in immunodeficient (NSG mice) or immune-competent (chicken embryo chorioallantoic membrane (CAM)) B-NHL xenograft models. Results: Using a panel of B-NHL co-cultures, we show that TG-1801, by disrupting the CD47-SIRPα axis, potentiates anti-CD20-mediated ADCC and ADCP. This led to a remarkable and durable antitumor effect of the triplet therapy composed by TG-1801 and U2 regimen, in vitro, as well as in mice and CAM xenograft models of B-NHL. Transcriptomic analysis also uncovered the upregulation of the G protein-coupled and inflammatory receptor, GPR183, as a crucial event associated with the efficacy of the triplet combination. Genetic depletion and pharmacological inhibition of GPR183 impaired ADCP initiation, cytoskeleton remodeling and cell migration in 2D and 3D spheroid B-NHL co-cultures, and disrupted macrophage-mediated control of tumor growth in B-NHL CAM xenografts. Conclusions: Altogether, our results support a crucial role for GPR183 in the recognition and elimination of malignant B cells upon concomitant targeting of CD20, CD47 and PI3Kδ, and warrant further clinical evaluation of this triplet regimen in B-NHL.

Alterations in p53, Microsatellite Stability and Lack of MUC5AC Expression as Molecular Features of Colorectal Carcinoma Associated with Inflammatory Bowel Disease.

Colitis-associated colorectal carcinoma (CAC) occurs in inflammatory bowel disease (IBD) because of the "chronic inflammation-dysplasia-cancer" carcinogenesis pathway characterized by p53 alterations in the early stages. Recently, gastric metaplasia (GM) has been described as the initial event of the serrated colorectal cancer (CRC) process, resulting from chronic stress on the colon mucosa. The aim of the study is to characterize CAC analyzing p53 alterations and microsatellite instability (MSI) to explore their relationship with GM using a series of CRC and the adjacent intestinal mucosa. Immunohistochemistry was performed to assess p53 alterations, MSI and MUC5AC expression as a surrogate for GM. The p53 mut-pattern was found in more than half of the CAC, most frequently stable (MSS) and MUC5AC negative. Only six tumors were unstable (MSI-H), being with p53 wt-pattern (p = 0.010) and MUC5AC positive (p = 0.005). MUC5AC staining was more frequently observed in intestinal mucosa, inflamed or with chronic changes, than in CAC, especially in those with p53 wt-pattern and MSS. Based on our results, we conclude that, as in the serrated pathway of CRC, in IBD GM occurs in inflamed mucosa, persists in those with chronic changes and disappears with the acquisition of p53 mutations.

A Large-Scale Exome-Wide Association Study Identifies Novel Germline Mutations in Lung Cancer.

Rationale: Genome-wide association studies have identified common variants of lung cancer. However, the contribution of rare exome-wide variants, especially protein-coding variants, to cancers remains largely unexplored. Objectives: To evaluate the role of human exomes in genetic predisposition to lung cancer. Methods: We performed exome-wide association studies to detect the association of exomes with lung cancer in 30,312 patients and 652,902 control subjects. A scalable and accurate implementation of a generalized mixed model was used to detect the association signals for loss-of-function, missense, and synonymous variants and gene-level sets. Furthermore, we performed association and Bayesian colocalization analyses to evaluate their relationships with intermediate exposures. Measurements and Main Results: We systematically analyzed 216,739 single-nucleotide variants in the human exome. The loss-of-function variants exhibited the most notable effects on lung cancer risk. We identified four novel variants, including two missense variants (rs202197044TET3 [Pmeta (P values of meta-analysis) = 3.60 × 10-8] and rs202187871POT1 [Pmeta = 2.21 × 10-8]) and two synonymous variants (rs7447927TMEM173 [Pmeta = 1.32 × 10-9] and rs140624366ATRN [Pmeta = 2.97 × 10-9]). rs202197044TET3 was significantly associated with emphysema (odds ratio, 3.55; Pfdr = 0.015), whereas rs7447927POT1 was strongly associated with telomere length (ß = 1.08; Pfdr (FDR corrected P value) = 3.76 × 10-53). Functional evidence of expression of quantitative trait loci, splicing quantitative trait loci, and isoform expression was found for the four novel genes. Gene-level association tests identified several novel genes, including POT1 (protection of telomeres 1), RTEL1, BSG, and ZNF232. Conclusions: Our findings provide insights into the genetic architecture of human exomes and their role in lung cancer predisposition.