Functional characterization of a multi-cancer risk locus on chr5p15.33 reveals regulation of TERT by ZNF148

E4F1 and ZNF148 are transcriptional activators of the -57A > C and wild-type TERT promoter

Point mutations within the TERT promoter are the most recurrent somatic noncoding mutations identified across different cancer types, including glioblastoma, melanoma, hepatocellular carcinoma, and bladder cancer. They are most abundant at -146C > T and -124C > T, and rarer at -57A > C, with the latter originally described as a familial case, but subsequently shown also to occur somatically. All three mutations create de novo E26-specific (ETS) binding sites and result in activation of the TERT gene, allowing cancer cells to achieve replicative immortality. Here, we used a systematic proteomics screen to identify transcription factors preferentially binding to the -146C > T, -124C > T, and -57A > C mutations. Although we confirmed binding of multiple ETS factors to the mutant -146C > T and -124C > T sequences, we identified E4F1 as a -57A > C-specific binder and ZNF148 as a TERT wild-type (WT) promoter binder that showed reduced interaction with the -124C > T allele. Both proteins are activating transcription factors that bind specifically to the -57A > C and WT (at position 124) TERT promoter sequence in corresponding cell lines, and up-regulate TERT transcription and telomerase activity. Our work describes new regulators of TERT gene expression with possible roles in cancer.

Human-genome single nucleotide polymorphisms affecting transcription factor binding and their role in pathogenesis

Single nucleotide polymorphisms (SNPs) are the most common type of variation in the human genome. The vast majority of SNPs identified in the human genome do not have any effect on the phenotype; however, some can lead to changes in the function of a gene or the level of its expression. Most SNPs associated with certain traits or pathologies are mapped to regulatory regions of the genome and affect gene expression by changing transcription factor binding sites. In recent decades, substantial effort has been invested in searching for such regulatory SNPs (rSNPs) and understanding the mechanisms by which they lead to phenotypic differences, primarily to individual differences in susceptibility to diseases and in sensitivity to drugs. The development of the NGS (next-generation sequencing) technology has contributed not only to the identification of a huge number of SNPs and to the search for their association (genome-wide association studies, GWASs) with certain diseases or phenotypic manifestations, but also to the development of more productive approaches to their functional annotation. It should be noted that the presence of an association does not allow one to identify a functional, truly disease-associated DNA sequence variant among multiple marker SNPs that are detected due to linkage disequilibrium. Moreover, determination of associations of genetic variants with a disease does not provide information about the functionality of these variants, which is necessary to elucidate the molecular mechanisms of the development of pathology and to design effective methods for its treatment and prevention. In this regard, the functional analysis of SNPs annotated in the GWAS catalog, both at the genome-wide level and at the level of individual SNPs, became especially relevant in recent years. A genome-wide search for potential rSNPs is possible without any prior knowledge of their association with a trait. Thus, mapping expression quantitative trait loci (eQTLs) makes it possible to identify an SNP for which - among transcriptomes of homozygotes and heterozygotes for its various alleles - there are differences in the expression level of certain genes, which can be located at various distances from the SNP. To predict rSNPs, approaches based on searches for allele-specific events in RNA-seq, ChIP-seq, DNase-seq, ATAC-seq, MPRA, and other data are also used. Nonetheless, for a more complete functional annotation of such rSNPs, it is necessary to establish their association with a trait, in particular, with a predisposition to a certain pathology or sensitivity to drugs. Thus, approaches to finding SNPs important for the development of a trait can be categorized into two groups: (1) starting from data on an association of SNPs with a certain trait, (2) starting from the determination of allele-specific changes at the molecular level (in a transcriptome or regulome). Only comprehensive use of strategically different approaches can considerably enrich our knowledge about the role of genetic determinants in the molecular mechanisms of trait formation, including predisposition to multifactorial diseases.

N -acetyltransferase 2 haplotype modifies risks for both dyslipidemia and urinary bladder cancer

A novel haplotype in N -acetyltransferase 2 ( NAT2 ) composed of seven non-coding variants (rs1495741, rs4921913, rs4921914, rs4921915, rs146812806, rs35246381, and rs35570672) has been linked to dyslipidemia by multiple, independent genome-wide association studies. The haplotype is located approximately 14 kb downstream of NAT2-coding region (ch8:18,272,377-18,272,881; GRCh38/hg38) and represents a non-coding, intergenic haplotype. Interestingly, the same dyslipidemia NAT2 haplotype is also linked to urinary bladder cancer risk. Dyslipidemia risk alleles are associated with rapid acetylator phenotype, whereas bladder cancer risk alleles are associated with slow acetylator, suggesting that the level of systemic NAT2 activity modifies the risk of these pathologies. We speculate that rs1495741 (and its associated haplotype) belongs to a distal regulatory element of human NAT2 gene (e.g., enhancer or silencer), and the genetic variation at the newly discovered haplotype results in a differential level of NAT2 gene expression. Understanding how this NAT2 haplotype contributes to not only urinary bladder cancer but also to dyslipidemia will ultimately help devise strategies to identify and protect susceptible individuals.

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.

Pancreatic Tumorigenesis: Precursors, Genetic Risk Factors and Screening

Pancreatic cancer (PC) is a highly malignant and aggressive tumor. Despite medical advancement, the silent nature of PC results in only 20% of all cases considered resectable at the time of diagnosis. It is projected to become the second leading cause in 2030. Most pancreatic cancer cases are diagnosed in the advanced stages. Such cases are typically unresectable and are associated with a 5-year survival of less than 10%. Although there is no guideline consensus regarding recommendations for screening for pancreatic cancer, early detection has been associated with better outcomes. In addition to continued utilization of imaging and conventional tumor markers, clinicians should be aware of novel testing modalities that may be effective for early detection of pancreatic cancer in individuals with high-risk factors. The pathogenesis of PC is not well understood; however, various modifiable and non-modifiable factors have been implicated in pancreatic oncogenesis. PC detection in the earlier stages is associated with better outcomes; nevertheless, most oncological societies do not recommend universal screening as it may result in a high false-positive rate. Therefore, targeted screening for high-risk individuals represents a reasonable option. In this review, we aimed to summarize the pathogenesis, genetic risk factors, high-risk population, and screening modalities for PC.

Functional studies of lung cancer GWAS beyond association

Fourteen years after the first genome-wide association study (GWAS) of lung cancer was published, approximately 45 genomic loci have now been significantly associated with lung cancer risk. While functional characterization was performed for several of these loci, a comprehensive summary of the current molecular understanding of lung cancer risk has been lacking. Further, many novel computational and experimental tools now became available to accelerate the functional assessment of disease-associated variants, moving beyond locus-by-locus approaches. In this review, we first highlight the heterogeneity of lung cancer GWAS findings across histological subtypes, ancestries and smoking status, which poses unique challenges to follow-up studies. We then summarize the published lung cancer post-GWAS studies for each risk-associated locus to assess the current understanding of biological mechanisms beyond the initial statistical association. We further summarize strategies for GWAS functional follow-up studies considering cutting-edge functional genomics tools and providing a catalog of available resources relevant to lung cancer. Overall, we aim to highlight the importance of integrating computational and experimental approaches to draw biological insights from the lung cancer GWAS results beyond association.

Targeted Long-Read Bisulfite Sequencing Identifies Differences in the TERT Promoter Methylation Profiles between TERT Wild-Type and TERT Mutant Cancer Cells

Background: TERT promoter methylation, located several hundred base pairs upstream of the transcriptional start site, is cancer specific and correlates with increased TERT mRNA expression and poorer patient outcome. Promoter methylation, however, is not mutually exclusive to TERT activating genetic alterations, as predicted for functionally redundant mechanisms. To annotate the altered patterns of TERT promoter methylation and their relationship with gene expression, we applied a Pacific Biosciences-based, long-read, bisulfite-sequencing technology and compared the differences in the methylation marks between wild-type and mutant cancers in an allele-specific manner. Results: We cataloged TERT genetic alterations (i.e., promoter point mutations or structural variations), allele-specific promoter methylation patterns, and allele-specific expression levels in a cohort of 54 cancer cell lines. In heterozygous mutant cell lines, the mutant alleles were significantly less methylated than their silent, mutation-free alleles (p < 0.05). In wild-type cell lines, by contrast, both epialleles were equally methylated to high levels at the TERT distal promoter, but differentially methylated in the proximal regions. ChIP analysis showed that epialleles with the hypomethylated proximal and core promoter were enriched in the active histone mark H3K4me2/3, whereas epialleles that were methylated in those regions were enriched in the repressive histone mark H3K27me3. Decitabine therapy induced biallelic expression in the wild-type cancer cells, whereas the mutant cell lines were unaffected. Conclusions: Long-read bisulfite sequencing analysis revealed differences in the methylation profiles and responses to demethylating agents between TERT wild-type and genetically altered cancer cell lines. The causal relation between TERT promoter methylation and gene expression remains to be established.

Transmembrane Protein 170B is a Prognostic Biomarker and Associated With Immune Infiltrates in Pancreatic Adenocarcinoma

Background: Pancreatic adenocarcinoma (PAAD) is among the most common types of cancer with a poor prognosis. Transmembrane protein 170B (TMEM170B) has been reported to suppress breast cancer proliferation, metastasis, and tumorigenesis and is related to prognosis. However, its role in PAAD and the underlying molecular mechanisms are yet to be investigated.

Patients and methods: We performed a comprehensive analysis of RNA sequencing data obtained from the Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) databases to determine TMEM170B expression. Immunostaining and real-time polymerase chain reaction (RT-PCR) were done to determine TMEM170B expression in human pancreatic cancer cell lines and tissue specimens. Furthermore, the correlation of TMEM170B with clinicopathological features and PAAD prognosis was investigated, and the mechanisms were explored through enrichment analysis and immune cell infiltration analysis.

Results: TCGA and GEO dataset analysis revealed that TMEM170B expression in PAAD tissue samples was significantly lower than that in non-tumorous tissues, which was further confirmed by immunohistochemistry and RT-PCR. Low TMEM170B expression was associated with poor differentiation (p = 0.014). Multivariate analysis identified that TMEM170B is an independent indicator for overall survival [hazard ratio (HR) = 0.116, 95% confidence interval (CI) = 0.014–0.995; p = 0.049] and disease-free survival (HR = 0.19, 95% CI = 0.04–0.910; p = 0.038) in patients with PAAD. Additionally, TMEM170B was involved in immune-related gene sets, including those related to chemokine signaling pathways and innate and adaptive immunity. High TMEM170B expression was linked to antitumor immune microenvironment with a high infiltration of B cells, T cells, dendritic cells, monocytes, M1 macrophages, neutrophil, and natural killer cells and a low infiltration of Tregs and myeloid-derived suppressor cells (all p &lt; 0.05).

Plain Language Summary: There is an urgent need to identify clinical prognostic biomarkers and targeted drugs for pancreatic cancer treatment. In this study, the expression status and prognostic value of transmembrane protein 170B (TMEM170B) in pancreatic adenocarcinoma were elucidated. Furthermore, TMEM170B, as a tumor suppressor gene, induced antitumor immune effects, including increased tumor infiltration of immune effector cells and reduced levels of inhibitory immune molecules and regulatory cells. Therefore, TMEM170B could be regarded as a novel target in preventing the progression of pancreatic cancer.

Conclusion: The findings suggest that low TMEM170B expression is remarkably correlated with poor PAAD prognosis, which might provide a therapeutic target for PAAD.

High throughput screen for the improvement of inducible promoters for tumor microenvironment cues

Cancer immunotherapies are highly potent and are gaining wide clinical usage. However, severe side effects require focusing effector immune cell activities on the tumor microenvironment (TME). We recently developed a chimeric antigen receptor tumor-induced vector (CARTIV), a synthetic promoter activated by TME factors. To improve CARTIV functions including background, activation levels, and synergism, we screened a library of promoters with variations in key positions. Here, we present a screening method involving turning ON/OFF stimulating TNFα and IFNγ cytokines, followed by sequential cell sorting. Sequencing of enriched promoters identified seventeen candidates, which were cloned and whose activities were then validated, leading to the identification of two CARTIVs with lower background and higher induction. We further combined a third hypoxia element with the two-factor CARTIV, demonstrating additional modular improvement. Our study presents a method of fine-tuning synthetic promoters for desired immunotherapy needs.

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.

Telomere length is not a main factor for the development of islet autoimmunity and type 1 diabetes in the TEDDY study

The Environmental Determinants of Diabetes in the Young (TEDDY) study enrolled 8676 children, 3-4 months of age, born with HLA-susceptibility genotypes for islet autoimmunity (IA) and type 1 diabetes (T1D). Whole-genome sequencing (WGS) was performed in 1119 children in a nested case-control study design. Telomere length was estimated from WGS data using five tools: Computel, Telseq, Telomerecat, qMotif and Motif_counter. The estimated median telomere length was 5.10 kb (IQR 4.52-5.68 kb) using Computel. The age when the blood sample was drawn had a significant negative correlation with telomere length (P = 0.003). European children, particularly those from Finland (P = 0.041) and from Sweden (P = 0.001), had shorter telomeres than children from the U.S.A. Paternal age (P = 0.019) was positively associated with telomere length. First-degree relative status, presence of gestational diabetes in the mother, and maternal age did not have a significant impact on estimated telomere length. HLA-DR4/4 or HLA-DR4/X children had significantly longer telomeres compared to children with HLA-DR3/3 or HLA-DR3/9 haplogenotypes (P = 0.008). Estimated telomere length was not significantly different with respect to any IA (P = 0.377), IAA-first (P = 0.248), GADA-first (P = 0.248) or T1D (P = 0.861). These results suggest that telomere length has no major impact on the risk for IA, the first step to develop T1D. Nevertheless, telomere length was shorter in the T1D high prevalence populations, Finland and Sweden.

Common variants in breast cancer risk loci predispose to distinct tumor subtypes

BACKGROUND

Genome-wide association studies (GWAS) have identified multiple common breast cancer susceptibility variants. Many of these variants have differential associations by estrogen receptor (ER) status, but how these variants relate with other tumor features and intrinsic molecular subtypes is unclear.

METHODS

Among 106,571 invasive breast cancer cases and 95,762 controls of European ancestry with data on 173 breast cancer variants identified in previous GWAS, we used novel two-stage polytomous logistic regression models to evaluate variants in relation to multiple tumor features (ER, progesterone receptor (PR), human epidermal growth factor receptor 2 (HER2) and grade) adjusting for each other, and to intrinsic-like subtypes.

RESULTS

Eighty-five of 173 variants were associated with at least one tumor feature (false discovery rate < 5%), most commonly ER and grade, followed by PR and HER2. Models for intrinsic-like subtypes found nearly all of these variants (83 of 85) associated at p < 0.05 with risk for at least one luminal-like subtype, and approximately half (41 of 85) of the variants were associated with risk of at least one non-luminal subtype, including 32 variants associated with triple-negative (TN) disease. Ten variants were associated with risk of all subtypes in different magnitude. Five variants were associated with risk of luminal A-like and TN subtypes in opposite directions.

CONCLUSION

This report demonstrates a high level of complexity in the etiology heterogeneity of breast cancer susceptibility variants and can inform investigations of subtype-specific risk prediction.

Dysplastic nevus syndrome and pancreatic cancer: A case report

Multiple primary cancers may occur in the same patient, with a prevalence that follows an ascendant trend. Their development is dictated by a complex interplay between a variety of factors, both patient-dependent and external. The case of a 38-year-old female patient diagnosed and treated for pancreatic cancer (PC) is presented in whom the digital dermoscopic monitoring of melanocytic nevi revealed a marked change of two nevi that acquired rapidly highly atypical features. They were surgically excised and the histopathological examination revealed two completely excised dysplastic compound nevi. Clinicians should be aware of the strong association between dysplastic nevus syndrome and PC, a malignancy associated with an extremely poor prognosis. Familial atypical multiple mole melanoma syndrome (FAMMM) predisposes to the development of melanoma, pancreatic cancer and other neoplasms. The common genetic background of PC and hereditary melanoma is discussed and the importance of regular skin checkup and screening for PC in these patients is underlined.

A 584 bp deletion in CTRB2 inhibits chymotrypsin B2 activity and secretion and confers risk of pancreatic cancer

Genome-wide association studies (GWASs) have discovered 20 risk loci in the human genome where germline variants associate with risk of pancreatic ductal adenocarcinoma (PDAC) in populations of European ancestry. Here, we fine-mapped one such locus on chr16q23.1 (rs72802365, p = 2.51 × 10-17, OR = 1.36, 95% CI = 1.31-1.40) and identified colocalization (PP = 0.87) with aberrant exon 5-7 CTRB2 splicing in pancreatic tissues (pGTEx = 1.40 × 10-69, βGTEx = 1.99; pLTG = 1.02 × 10-30, βLTG = 1.99). Imputation of a 584 bp structural variant overlapping exon 6 of CTRB2 into the GWAS datasets resulted in a highly significant association with pancreatic cancer risk (p = 2.83 × 10-16, OR = 1.36, 95% CI = 1.31-1.42), indicating that it may underlie this signal. Exon skipping attributable to the deletion (risk) allele introduces a premature stop codon in exon 7 of CTRB2, yielding a truncated chymotrypsinogen B2 protein that lacks chymotrypsin activity, is poorly secreted, and accumulates intracellularly in the endoplasmic reticulum (ER). We propose that intracellular accumulation of a nonfunctional chymotrypsinogen B2 protein leads to ER stress and pancreatic inflammation, which may explain the increased pancreatic cancer risk in carriers of CTRB2 exon 6 deletion alleles.

Large-scale cross-cancer fine-mapping of the 5p15.33 region reveals multiple independent signals

Genome-wide association studies (GWASs) have identified thousands of cancer risk loci revealing many risk regions shared across multiple cancers. Characterizing the cross-cancer shared genetic basis can increase our understanding of global mechanisms of cancer development. In this study, we collected GWAS summary statistics based on up to 375,468 cancer cases and 530,521 controls for fourteen types of cancer, including breast (overall, estrogen receptor [ER]-positive, and ER-negative), colorectal, endometrial, esophageal, glioma, head/neck, lung, melanoma, ovarian, pancreatic, prostate, and renal cancer, to characterize the shared genetic basis of cancer risk. We identified thirteen pairs of cancers with statistically significant local genetic correlations across eight distinct genomic regions. Specifically, the 5p15.33 region, harboring the TERT and CLPTM1L genes, showed statistically significant local genetic correlations for multiple cancer pairs. We conducted a cross-cancer fine-mapping of the 5p15.33 region based on eight cancers that showed genome-wide significant associations in this region (ER-negative breast, colorectal, glioma, lung, melanoma, ovarian, pancreatic, and prostate cancer). We used an iterative analysis pipeline implementing a subset-based meta-analysis approach based on cancer-specific conditional analyses and identified ten independent cross-cancer associations within this region. For each signal, we conducted cross-cancer fine-mapping to prioritize the most plausible causal variants. Our findings provide a more in-depth understanding of the shared inherited basis across human cancers and expand our knowledge of the 5p15.33 region in carcinogenesis.

A Panel of rSNPs Demonstrating Allelic Asymmetry in Both ChIP-seq and RNA-seq Data and the Search for Their Phenotypic Outcomes through Analysis of DEGs

Currently, the detection of the allele asymmetry of gene expression from RNA-seq data or the transcription factor binding from ChIP-seq data is one of the approaches used to identify the functional genetic variants that can affect gene expression (regulatory SNPs or rSNPs). In this study, we searched for rSNPs using the data for human pulmonary arterial endothelial cells (PAECs) available from the Sequence Read Archive (SRA). Allele-asymmetric binding and expression events are analyzed in paired ChIP-seq data for H3K4me3 mark and RNA-seq data obtained for 19 individuals. Two statistical approaches, weighted z-scores and predicted probabilities, were used to improve the efficiency of finding rSNPs. In total, we identified 14,266 rSNPs associated with both allele-specific binding and expression. Among them, 645 rSNPs were associated with GWAS phenotypes; 4746 rSNPs were reported as eQTLs by GTEx, and 11,536 rSNPs were located in 374 candidate transcription factor binding motifs. Additionally, we searched for the rSNPs associated with gene expression using an SRA RNA-seq dataset for 281 clinically annotated human postmortem brain samples and detected eQTLs for 2505 rSNPs. Based on these results, we conducted Gene Ontology (GO), Disease Ontology (DO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses and constructed the protein-protein interaction networks to represent the top-ranked biological processes with a possible contribution to the phenotypic outcome.

Regulatory SNPs: Altered Transcription Factor Binding Sites Implicated in Complex Traits and Diseases

The vast majority of the genetic variants (mainly SNPs) associated with various human traits and diseases map to a noncoding part of the genome and are enriched in its regulatory compartment, suggesting that many causal variants may affect gene expression. The leading mechanism of action of these SNPs consists in the alterations in the transcription factor binding via creation or disruption of transcription factor binding sites (TFBSs) or some change in the affinity of these regulatory proteins to their cognate sites. In this review, we first focus on the history of the discovery of regulatory SNPs (rSNPs) and systematized description of the existing methodical approaches to their study. Then, we brief the recent comprehensive examples of rSNPs studied from the discovery of the changes in the TFBS sequence as a result of a nucleotide substitution to identification of its effect on the target gene expression and, eventually, to phenotype. We also describe state-of-the-art genome-wide approaches to identification of regulatory variants, including both making molecular sense of genome-wide association studies (GWAS) and the alternative approaches the primary goal of which is to determine the functionality of genetic variants. Among these approaches, special attention is paid to expression quantitative trait loci (eQTLs) analysis and the search for allele-specific events in RNA-seq (ASE events) as well as in ChIP-seq, DNase-seq, and ATAC-seq (ASB events) data.

SMDB: pivotal somatic sequence alterations reprogramming regulatory cascades

Binding motifs for transcription factors, RNA-binding proteins, microRNAs (miRNAs), etc. are vital for proper gene transcription and translation regulation. Sequence alteration mechanisms including single nucleotide mutations, insertion, deletion, RNA editing and single nucleotide polymorphism can lead to gains and losses of binding motifs; such consequentially emerged or vanished binding motifs are termed ‘somatic motifs’ by us. Somatic motifs have been studied sporadically but have never been curated into a comprehensive resource. By analyzing various types of sequence altering data from large consortiums, we successfully identified millions of somatic motifs, including those for important transcription factors, RNA-binding proteins, miRNA seeds and miRNA–mRNA 3′-UTR target motifs. While a few of these somatic motifs have been well studied, our results contain many novel somatic motifs that occur at high frequency and are thus likely to cause important biological repercussions. Genes targeted by these altered motifs are excellent candidates for further mechanism studies. Here, we present the first database that hosts millions of somatic motifs ascribed to a variety of sequence alteration mechanisms.

A Transcriptome-Wide Association Study Identifies Novel Candidate Susceptibility Genes for Pancreatic Cancer

BACKGROUND

Although 20 pancreatic cancer susceptibility loci have been identified through genome-wide association studies in individuals of European ancestry, much of its heritability remains unexplained and the genes responsible largely unknown.

METHODS

To discover novel pancreatic cancer risk loci and possible causal genes, we performed a pancreatic cancer transcriptome-wide association study in Europeans using three approaches: FUSION, MetaXcan, and Summary-MulTiXcan. We integrated genome-wide association studies summary statistics from 9040 pancreatic cancer cases and 12 496 controls, with gene expression prediction models built using transcriptome data from histologically normal pancreatic tissue samples (NCI Laboratory of Translational Genomics [n = 95] and Genotype-Tissue Expression v7 [n = 174] datasets) and data from 48 different tissues (Genotype-Tissue Expression v7, n = 74-421 samples).

RESULTS

We identified 25 genes whose genetically predicted expression was statistically significantly associated with pancreatic cancer risk (false discovery rate < .05), including 14 candidate genes at 11 novel loci (1p36.12: CELA3B; 9q31.1: SMC2, SMC2-AS1; 10q23.31: RP11-80H5.9; 12q13.13: SMUG1; 14q32.33: BTBD6; 15q23: HEXA; 15q26.1: RCCD1; 17q12: PNMT, CDK12, PGAP3; 17q22: SUPT4H1; 18q11.22: RP11-888D10.3; and 19p13.11: PGPEP1) and 11 at six known risk loci (5p15.33: TERT, CLPTM1L, ZDHHC11B; 7p14.1: INHBA; 9q34.2: ABO; 13q12.2: PDX1; 13q22.1: KLF5; and 16q23.1: WDR59, CFDP1, BCAR1, TMEM170A). The association for 12 of these genes (CELA3B, SMC2, and PNMT at novel risk loci and TERT, CLPTM1L, INHBA, ABO, PDX1, KLF5, WDR59, CFDP1, and BCAR1 at known loci) remained statistically significant after Bonferroni correction.

CONCLUSIONS

By integrating gene expression and genotype data, we identified novel pancreatic cancer risk loci and candidate functional genes that warrant further investigation.

Massively parallel reporter assays of melanoma risk variants identify MX2 as a gene promoting melanoma

Genome-wide association studies (GWAS) have identified ~20 melanoma susceptibility loci, most of which are not functionally characterized. Here we report an approach integrating massively-parallel reporter assays (MPRA) with cell-type-specific epigenome and expression quantitative trait loci (eQTL) to identify susceptibility genes/variants from multiple GWAS loci. From 832 high-LD variants, we identify 39 candidate functional variants from 14 loci displaying allelic transcriptional activity, a subset of which corroborates four colocalizing melanocyte cis-eQTL genes. Among these, we further characterize the locus encompassing the HIV-1 restriction gene, MX2 (Chr21q22.3), and validate a functional intronic variant, rs398206. rs398206 mediates the binding of the transcription factor, YY1, to increase MX2 levels, consistent with the cis-eQTL of MX2 in primary human melanocytes. Melanocyte-specific expression of human MX2 in a zebrafish model demonstrates accelerated melanoma formation in a BRAF^V600E background. Our integrative approach streamlines GWAS follow-up studies and highlights a pleiotropic function of MX2 in melanoma susceptibility.

There are more than 20 known melanoma susceptibility genes. Here, using a massively parallel reporter assay, the authors identify risk-associated variants that alter gene transcription, and demonstrate that expression of one such gene, MX2, leads to the promotion of melanoma in a zebrafish model.

Validation of genome-wide association study-identified single nucleotide polymorphisms in a case-control study of pancreatic cancer from Taiwan

Background

Due to differences in genetic background, it is unclear whether the genetic loci identified by the previous genome-wide association studies (GWAS) of pancreatic cancer also play significant roles in the development of pancreatic cancer among the Taiwanese population.

Methods

This study aimed to validate the 25 pancreatic cancer GWAS-identified single nucleotide polymorphisms (SNPs) in a case-control study (278 cases and 658 controls) of pancreatic cancer conducted in Taiwan. Statistical analyses were conducted to determine the associations between the GWAS-identified SNPs and pancreatic cancer risk. Gene-environment interaction analysis was conducted to evaluate the interactions between SNPs and environmental factors on pancreatic cancer risk.

Results

Among the 25 GWAS-identified SNPs, 7 (rs2816938 (~ 11 kb upstream of NR5A2), rs10094872 (~ 28 kb upstream of MYC), rs9581943 (200 bp upstream of PDX1) and 4 chromosome 13q22.1 SNPs: rs4885093, rs9573163, rs9543325, rs9573166) showed a statistically significant association with pancreatic cancer risk in the current study. Additional analyses showed two significant gene-environment interactions (between poor oral hygiene and NR5A2 rs2816938 and between obesity and PDX1 rs9581943) on the risk of pancreatic cancer.

Conclusions

The current study confirmed the associations between 7 of the 25 GWAS-identified SNPs and pancreatic risk among the Taiwanese population.

Furthermore, pancreatic cancer was jointly influenced by lifestyle and medical factors, genetic polymorphisms, and gene-environment interaction. Additional GWAS is needed to determine the genetic polymorphisms that are more relevant to the pancreatic cancer cases occurring in Taiwan.

Telomeres and Telomere Length: A General Overview

Telomeres are highly conserved tandem nucleotide repeats that include proximal double-stranded and distal single-stranded regions that in complex with shelterin proteins afford protection at chromosomal ends to maintain genomic integrity. Due to the inherent limitations of DNA replication and telomerase suppression in most somatic cells, telomeres undergo age-dependent incremental attrition. Short or dysfunctional telomeres are recognized as DNA double-stranded breaks, triggering cells to undergo replicative senescence. Telomere shortening, therefore, acts as a counting mechanism that drives replicative senescence by limiting the mitotic potential of cells. Telomere length, a complex hereditary trait, is associated with aging and age-related diseases. Epidemiological data, in general, support an association with varying magnitudes between constitutive telomere length and several disorders, including cancers. Telomere attrition is also influenced by oxidative damage and replicative stress caused by genetic, epigenetic, and environmental factors. Several single nucleotide polymorphisms at different loci, identified through genome-wide association studies, influence inter-individual variation in telomere length. In addition to genetic factors, environmental factors also influence telomere length during growth and development. Telomeres hold potential as biomarkers that reflect the genetic predisposition together with the impact of environmental conditions and as targets for anti-cancer therapies.

ZBP-89 negatively regulates self-renewal of liver cancer stem cells via suppression of Notch1 signaling pathway

Liver cancer stem cells (LCSCs) initiate hepatocellular carcinoma (HCC) and contribute to its recurrence and treatment resistance. Studies have suggested ZBP-89 as a candidate tumor suppressor in HCC. We explored the role of ZBP-89 in the regulation of LCSCs. This study was performed in liver tissue samples from 104 HCC patients, 2 cell lines and mouse tumor models. We demonstrated that ZBP-89 was weakly expressed in LCSCs. Patients with high expression of LCSC markers displayed reduced survivals and higher recurrence rates after curative surgical operation. The expression of ZBP-89 was predictive for decreased recurrence. LCSC markers were negatively correlated with ZBP-89 in HCC tissues and in enriched liver tumor spheres. The exogenous expression of ZBP-89 attenuated the tumor-sphere formation and secondary colony formation capabilities of LCSCs in vitro and tumorigenicity in vivo. Furthermore, the negative effect of ZBP-89 on cancer stemness was Notch1-dependent. Localized with Notch1 intracellular domain (NICD1) in the nucleus, ZBP-89 repressed the Notch1 signaling pathway by competitive binding to NICD1 with MAML1. Collectively, ZBP-89 negatively regulates HCC stemness via inhibiting the Notch1 signaling.

PAQMAN: Protein-nucleic acid affinity quantification by MAss spectrometry in nuclear extracts

In recent years, various mass spectrometry-based approaches have been developed to determine global protein-DNA binding specificities using DNA affinity purifications from crude nuclear extracts. However, these assays are semi-quantitative and do not provide information about interaction affinities. We recently developed a technology that we call Protein-nucleic acid Affinity Quantification by MAss spectrometry in Nuclear extracts or PAQMAN, that can be used to determine apparent affinities between multiple nuclear proteins and a nucleic acid sequence of interest in one experiment. In PAQMAN, a series of affinity purifications with increasing bait concentrations and fixed amounts of crude nuclear extracts are combined with isobaric stable isotope labeling and quantitative mass spectrometry to generate Hill-like K_d curves for dozens of proteins in a single experiment. Here, we apply PAQMAN to determine apparent affinities for a genetic variant, rs36115365-C, which regulates TERT expression and is associated with an increased risk to develop various malignancies. Furthermore, we describe a detailed protocol for this method including important quality checks.

Telomere length and associated factors in older adults with hypertension

Telomeres and telomerase play important roles in the occurrence and development of hypertension. This review was performed to clarify the factors that influence telomere length and telomerase activity in older patients and elucidate the association of these factors with hypertension. A PubMed search and critical review of studies assessing the risk factors underlying the association of hypertension with telomere length and telomerase activity was performed. Telomere length and telomerase activity were found to be associated with hypertension. The factors that influence telomere length and telomerase activity in older patients with hypertension include genetics, demographics, social and environmental factors, chronic disease, psychological factors, and antihypertensive drug treatment. A better understanding of the molecular mechanisms underlying the association of hypertension with telomere length and telomerase activity may help to reduce the incidence of hypertension.

Lung Cancer Risk in Never-Smokers of European Descent is Associated With Genetic Variation in the 5p15.33 TERT-CLPTM1Ll Region

INTRODUCTION

Inherited susceptibility to lung cancer risk in never-smokers is poorly understood. The major reason for this gap in knowledge is that this disease is relatively uncommon (except in Asians), making it difficult to assemble an adequate study sample. In this study we conducted a genome-wide association study on the largest, to date, set of European-descent never-smokers with lung cancer.

METHODS

We conducted a two-phase (discovery and replication) genome-wide association study in never-smokers of European descent. We further augmented the sample by performing a meta-analysis with never-smokers from the recent OncoArray study, which resulted in a total of 3636 cases and 6295 controls. We also compare our findings with those in smokers with lung cancer.

RESULTS

We detected three genome-wide statistically significant single nucleotide polymorphisms rs31490 (odds ratio [OR]: 0.769, 95% confidence interval [CI]: 0.722-0.820; p value 5.31 × 10-16), rs380286 (OR: 0.770, 95% CI: 0.723-0.820; p value 4.32 × 10-16), and rs4975616 (OR: 0.778, 95% CI: 0.730-0.829; p value 1.04 × 10-14). All three mapped to Chromosome 5 CLPTM1L-TERT region, previously shown to be associated with lung cancer risk in smokers and in never-smoker Asian women, and risk of other cancers including breast, ovarian, colorectal, and prostate.

CONCLUSIONS

We found that genetic susceptibility to lung cancer in never-smokers is associated to genetic variants with pan-cancer risk effects. The comparison with smokers shows that top variants previously shown to be associated with lung cancer risk only confer risk in the presence of tobacco exposure, underscoring the importance of gene-environment interactions in the etiology of this disease.

Assessing a single SNP located at TERT/CLPTM1L multi-cancer risk region as a genetic modifier for risk of pancreatic cancer and melanoma in Dutch CDKN2A mutation carriers

Carriers of pathogenic variants in CDKN2A have a 70% life-time risk of developing melanoma and 15–20% risk of developing pancreatic cancer (PC). In the Netherlands, a 19-bp deletion in exon 2 of CDKN2A (p16-Leiden mutation) accounts for most hereditary melanoma cases. Clinical experience suggests variability in occurrence of melanoma and PC in p16-Leiden families. Thereby, the risk of developing cancer could be modified by both environmental and genetic contributors, suggesting that identification of genetic modifiers could improve patients’ surveillance. In a recent genome-wide association study (GWAS), rs36115365-C was found to significantly modify risk of PC and melanoma in the European population. This SNP is located on chr5p15.33 and has allele-specific regulatory activities on TERT expression. Herein, we investigated the modifying capacities of rs36115365-C on PC and melanoma in a cohort of 283 p16-Leiden carriers including 29 diagnosed with PC, 171 diagnosed with melanoma, 21 diagnosed with both PC and melanoma and 62 with neither PC nor melanoma. In contrast to previously reported findings, we did not find a significant association of PC risk with risk variant presence as determined by Generalized Estimating Equations (GEE) modelling. Interestingly, carrier-ship of the risk variant had a significant protective effect for melanoma (OR − 0.703 [95% CI − 1.201 to − 0.205], p = 0.006); however, the observed association was no longer significant after exclusion of probands to assess possible influence of ascertainment. Collectively, genetic modifiers for the prediction of PC and melanoma risk in p16-Leiden carriers remain to be determined.

Telomere length, telomerase reverse transcriptase promoter mutations, and melanoma risk

Telomere repeats at chromosomal ends, critical for genomic integrity, undergo age-dependent attrition and telomere length has been associated with different disorders including cancers. In this study, based on 1469 patients and 1158 healthy controls, we show a statistically significant (P = 6 × 10^-10 ) association between increased telomere length and melanoma risk. Mendelian randomization, using 5 telomere length-associated polymorphisms, ruled out confounding factors or reverse causality and showed association between increased telomere length and melanoma risk with odds ratio of 2.66 (95% confidence interval: 2.07-3.25). Age-dependent telomere attrition was faster in melanoma cases than controls (P = .01). The carriers of a highly penetrant germline -57A>C TERT promoter mutation, in a previously reported melanoma family, had longer telomeres than the noncarriers. The mutation causes increased TERT and telomerase levels through creation of a binding motif for E-twenty six (ETS) transcription factors and the carriers develop melanoma with an early age of onset and rapid progression to metastasis. In analogy, we hypothesize that increased telomere length in melanoma patients reflects stochastic increased telomerase levels due to common genetic variation. Paradoxically, we observed shorter telomeres (P = 1 × 10^-5 ) in primary tumors from unrelated melanoma patients with (121) than without (170) somatic TERT promoter mutations that similar to the germline mutation, also create binding motifs for ETS transcription factors. However, the age-dependent telomere attrition was faster in tumors with the TERT promoter mutations than in those without such mutations. Besides a robust association between increased telomere length and risk, our data show a perturbed telomere homeostasis in melanoma.

Telomere length and survival in primary cutaneous melanoma patients

Telomere repeats at chromosomal ends, critical to genomic integrity, undergo age-dependent attrition. Telomere length, a polygenic trait, has been associated with risk of several disorders including cancers. In contrast to association of long telomeres with increased risk of several cancers, including melanoma, emerging reports suggest that short telomeres predict poor survival in patients with different cancers. In this study based on 1019 stage I and II cutaneous melanoma patients, we show an association between the patients with short telomeres and poor melanoma-specific survival (HR 2.05, 95% CI 1.33-3.16) compared to patients with long telomeres. Due to inverse correlation between age and telomere length (r -0.19, P < 0.0001), we stratified the patients into quantiles based on age at diagnosis and also carried out age-matched analysis. The effect of short telomeres on survival was determined by using multivariate Cox regression that included composite genetic risk score computed from genotyping of the patients for telomere-length associated polymorphisms. The effect of decreased telomere length on poor melanoma-specific survival was particularly strong in patients within the age quantile below 30 years (HR 3.82, 95% CI 1.10-13.30) and between 30-40 years (HR 2.69, 95% CI 1.03-7.03). Our study shows that in contrast to increased melanoma risk associated with increased telomere length, decreased telomere length predicts poor survival in melanoma subgroups.

Peritumoral overexpression of ZBP-89 is associated with unfavorable disease-free survival rates in patients with hepatocellular carcinoma following hepatectomy

Previous studies have revealed that the peritumoral environment has a profound influence on tumor initiation and progression. Zinc-binding protein-89 (ZBP-89) has been observed to be involved with tumor development, recurrence, and metastasis. High intratumoral expression of ZBP-89 has been associated with improved prognosis in several tumor types. However, the prognostic values of peritumoral expression of ZBP-89 remain to be elucidated in patients with hepatocellular carcinoma (HCC) following curative resection. In the present study, peritumoral ZBP-89 expression was examined using immunohistochemistry in 102 HCC patients who had received curative hepatectomy. Expression of ZBP-89 protein was positive in 66.3% of the peritumoral samples from 102 HCC patients. HCC patients with high peritumoral ZBP-89 expression exhibited significantly shorter disease-free survival (DFS) times (P=0.012) than those patients with low peritumoral ZBP-89 expression. Additionally, high ZBP-89 expression in peritumoral HCC tissue was positively associated with the presence of liver cirrhosis. Univariate and multivariate Cox proportional hazard regression analyses demonstrated that albumin levels ≤35 g/l, multiple tumors, tumor sizes ≥5 cm, and macroscopic vascular invasion may serve as independent prognostic factors for overall survival (OS) [hazard ratio (HR)=2.031; P=0.014] in patients with HCC. The multivariate Cox regression model identified that high ZBP-89 expression, multiple tumors and macroscopic vascular invasion were independent prognostic factors for shorter DFS durations. High expression of ZBP-89 in peritumoral HCC tissues was associated with a shorter DFS in HCC patients following curative hepatectomy. Additionally, high ZBP-89 expression in peritumoral HCC tissue was positively associated with the presence of liver cirrhosis in HCC patients, indicating that cirrhosis accompanied by high ZBP-89 expression may be a contributing factor to the poor prognosis of patients with HCC. Therefore, peritumoral ZBP-89 expression may be a good prognostic marker to predict DFS time in HCC patients following curative hepatectomy and may provide novel insights into the molecular mechanisms of HCC initiation.

Beginning at the ends: telomeres and human disease

Studies of rare and common illnesses have led to remarkable progress in the understanding of the role of telomeres (nucleoprotein complexes at chromosome ends essential for chromosomal integrity) in human disease. Telomere biology disorders encompass a growing spectrum of conditions caused by rare pathogenic germline variants in genes encoding essential aspects of telomere function. Dyskeratosis congenita, a disorder at the severe end of this spectrum, typically presents in childhood with the classic triad of abnormal skin pigmentation, nail dystrophy, and oral leukoplakia, accompanied by a very high risk of bone marrow failure, cancer, pulmonary fibrosis, and other medical problems. In contrast, the less severe end of the telomere biology disorder spectrum consists of middle-age or older adults with just one feature typically seen in dyskeratosis congenita, such as pulmonary fibrosis or bone marrow failure. In the common disease realm, large-scale molecular epidemiology studies have discovered novel associations between illnesses, such as cancer, heart disease, and mental health, and both telomere length and common genetic variants in telomere biology genes. This review highlights recent findings of telomere biology in human disease from both the rare and common disease perspectives. Multi-disciplinary collaborations between clinicians, basic scientists, and epidemiologist are essential as we seek to incorporate new telomere biology discoveries to improve health outcomes.