Functional characterization of a chr13q22.1 pancreatic cancer risk locus reveals long-range interaction and allele-specific effects on DIS3 expression

Interpretation of the role of germline and somatic non-coding mutations in cancer: expression and chromatin conformation informed analysis

BACKGROUND

There has been extensive scrutiny of cancer driving mutations within the exome (especially amino acid altering mutations) as these are more likely to have a clear impact on protein functions, and thus on cell biology. However, this has come at the neglect of systematic identification of regulatory (non-coding) variants, which have recently been identified as putative somatic drivers and key germline risk factors for cancer development. Comprehensive understanding of non-coding mutations requires understanding their role in the disruption of regulatory elements, which then disrupt key biological functions such as gene expression.

MAIN BODY

We describe how advancements in sequencing technologies have led to the identification of a large number of non-coding mutations with uncharacterized biological significance. We summarize the strategies that have been developed to interpret and prioritize the biological mechanisms impacted by non-coding mutations, focusing on recent annotation of cancer non-coding variants utilizing chromatin states, eQTLs, and chromatin conformation data.

CONCLUSION

We believe that a better understanding of how to apply different regulatory data types into the study of non-coding mutations will enhance the discovery of novel mechanisms driving cancer.

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 < 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.

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.

Susceptibility loci for pancreatic cancer in the Brazilian population

BACKGROUND

Pancreatic adenocarcinoma (PA) is a very aggressive cancer and has one of the poorest prognoses. Usually, the diagnosis is late and resistant to conventional treatment. Environmental and genetic factors contribute to the etiology, such as tobacco and alcohol consumption, chronic pancreatitis, diabetes and obesity. Somatic mutation in pancreatic cancer cells are known and SNP profile by GWAS could access novel genetic risk factors for this disease in different population context. Here we describe a SNP panel for Brazilian pancreatic cancer, together with clinical and epidemiological data.

METHODS

78 pancreatic adenocarcinoma and 256 non-pancreatic cancer subjects had 25 SNPs genotyped by real-time PCR. Unconditional logistic regression methods were used to assess the main effects on PA risk, using allelic, co-dominant and dominant inheritance models.

RESULTS

9 SNPs were nominally associated with pancreatic adenocarcinoma risk, with 5 of the minor alleles conferring protective effect while 4 related as risk factor. In epidemiological and clinical data, tobacco smoking, diabetes and pancreatitis history were significantly related to pancreatic adenocarcinoma risk. Polygenic risk scores computed using the SNPs in the study showed strong associations with PA risk.

CONCLUSION

We could assess for the first time some SNPs related with PA in Brazilian populations, a result that could be used for genetic screening in risk population such as familial pancreatic cancer, smokers, alcohol users and diabetes patients.

Identification of RNA-binding protein SNRPA1 for prognosis in prostate cancer

Prostate cancer is one of the deadliest cancers in men. RNA-binding proteins play a critical role in human cancers; however, whether they have a significant effect on the prognosis of prostate cancer has yet to be elucidated. In the present study, we performed a comprehensive analysis of RNA sequencing and clinical data from the Cancer Genome Atlas dataset and obtained differentially expressed RNA-binding proteins between prostate cancer and benign tissues. We constructed a protein-protein interaction network and Cox regression analyses were conducted to identify prognostic hub RNA-binding proteins. SNRPA1 was associated with the highest risk of poor prognosis and was therefore selected for further analysis. SNRPA1 expression was positively correlated with Gleason score and pathological TNM stage in prostate cancer patients. Furthermore, the expression profile of SNRPA1 was validated using the Oncomine, Human Protein Atlas, and Cancer Cell Line Encyclopedia databases. Meanwhile, the prognostic profile of SNRPA1 was successfully verified in GSE70769. Additionally, the results of molecular experiments revealed the proliferative role of SNRPA1 in prostate cancer cells. In summary, our findings evidenced a relationship between RNA-binding proteins and prostate cancer and indicated the prognostic significance of SNRPA1 in prostate cancer.

Genome-wide association study of prevalent and persistent cervical high-risk human papillomavirus (HPV) infection

Background

Genetic factors may influence the susceptibility to high-risk (hr) human papillomavirus (HPV) infection and persistence. We conducted the first genome-wide association study (GWAS) to identify variants associated with cervical hrHPV infection and persistence.

Methods

Participants were 517 Nigerian women evaluated at baseline and 6 months follow-up visits for HPV. HPV was characterized using SPF₁₀/LiPA₂₅. hrHPV infection was positive if at least one carcinogenic HPV genotype was detected in a sample provided at the baseline visit and persistent if at least one carcinogenic HPV genotype was detected in each of the samples provided at the baseline and follow-up visits. Genotyping was done using the Illumina Multi-Ethnic Genotyping Array (MEGA) and imputation was done using the African Genome Resources Haplotype Reference Panel. Association analysis was done for hrHPV infection (125 cases/392 controls) and for persistent hrHPV infection (51 cases/355 controls) under additive genetic models adjusted for age, HIV status and the first principal component (PC) of the genotypes.

Results

The mean (±SD) age of the study participants was 38 (±8) years, 48% were HIV negative, 24% were hrHPV positive and 10% had persistent hrHPV infections. No single variant reached genome-wide significance (p < 5 X 10^− 8). The top three variants associated with hrHPV infections were intronic variants clustered in KLF12 (all OR: 7.06, p = 1.43 × 10^− 6). The top variants associated with cervical hrHPV persistence were in DAP (OR: 6.86, p = 7.15 × 10^− 8), NR5A2 (OR: 3.65, p = 2.03 × 10^− 7) and MIR365–2 (OR: 7.71, p = 2.63 × 10^− 7) gene regions.

Conclusions

This exploratory GWAS yielded suggestive candidate risk loci for cervical hrHPV infection and persistence. The identified loci have biological annotation and functional data supporting their role in hrHPV infection and persistence. Given our limited sample size, larger discovery and replication studies are warranted to further characterize the reported associations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12881-020-01156-1.

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.

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.

Emerging epigenomic landscapes of pancreatic cancer in the era of precision medicine

Genetic studies have advanced our understanding of pancreatic cancer at a mechanistic and translational level. Genetic concepts and tools are increasingly starting to be applied to clinical practice, in particular for precision medicine efforts. However, epigenomics is rapidly emerging as a promising conceptual and methodological paradigm for advancing the knowledge of this disease. More importantly, recent studies have uncovered potentially actionable pathways, which support the prediction that future trials for pancreatic cancer will involve the vigorous testing of epigenomic therapeutics. Thus, epigenomics promises to generate a significant amount of new knowledge of both biological and medical importance.

In pancreatic cancer, the epigenomic landscape can strongly impact the disease phenotype. Here, the authors discuss recent advances in our understanding of pancreatic cancer epigenomics, and how this knowledge can integrate with precision medicine approaches in this lethal disease.

Genetic variants in SLC22A3 contribute to the susceptibility to colorectal cancer

Previous a genome-wide association study (GWAS) of colorectal cancer in Japanese population has identified a risk region at the chromosome 6q26-q27 associated with colorectal cancer risk. However, the causal gene at this locus remained unclear. In our study, we enrolled a total of 14 candidate functional single nucleotide polymorphisms (SNPs) at 6q26-q27 (318 kb), and then genotyped them by TaqMan method in a Chinese population including 1,147 colorectal cancer cases and 1,203 controls. Among that, 5 SNPs were identified statistical association with colorectal cancer risk by logistic regression analysis. Of which, SNP rs420038 G > A in SLC22A3 was related to decreased risk of colorectal cancer (adjusted odds ratio (OR) = 0.79, 95% confidence interval (CI) = 0.67-0.94, p = 0.007), and also associated with lower expression of SLC22A3 (p = 0.040) using expression quantitative trait loci (eQTL) analysis. Moreover, by the luciferase assays, we found that compared to the G allele of rs420038, the A allele could suppress the activity of the promoter in SLC22A3. Furthermore, the expression of SLC22A3 was significantly higher in colorectal cancer tissues than that in paired normal tissues (p < 0.001). Meanwhile, the phenotypes of proliferation, migration, invasion, cell cycle and apoptosis of colorectal cancer cell were significantly affected by SLC22A3 in vitro. Our results revealed a novel susceptible locus, rs420038 in SLC22A3, which may be involved in colorectal cancer development and progression.

Novel insights into chromosomal conformations in cancer

Exploring gene function is critical for understanding the complexity of life. DNA sequences and the three-dimensional organization of chromatin (chromosomal interactions) are considered enigmatic factors underlying gene function, and interactions between two distant fragments can regulate transactivation activity via mediator proteins. Thus, a series of chromosome conformation capture techniques have been developed, including chromosome conformation capture (3C), circular chromosome conformation capture (4C), chromosome conformation capture carbon copy (5C), and high-resolution chromosome conformation capture (Hi-C). The application of these techniques has expanded to various fields, but cancer remains one of the major topics. Interactions mediated by proteins or long noncoding RNAs (lncRNAs) are typically found using 4C-sequencing and chromatin interaction analysis by paired-end tag sequencing (ChIA-PET). Currently, Hi-C is used to identify chromatin loops between cancer risk-associated single-nucleotide polymorphisms (SNPs) found by genome-wide association studies (GWAS) and their target genes. Chromosomal conformations are responsible for altered gene regulation through several typical mechanisms and contribute to the biological behavior and malignancy of different tumors, particularly prostate cancer, breast cancer and hematologic neoplasms. Moreover, different subtypes may exhibit different 3D-chromosomal conformations. Thus, C-tech can be used to help diagnose cancer subtypes and alleviate cancer progression by destroying specific chromosomal conformations. Here, we review the fundamentals and improvements in chromosome conformation capture techniques and their clinical applications in cancer to provide insight for future research.

Whole-exome analysis in osteosarcoma to identify a personalized therapy

Osteosarcoma is the most common pediatric primary non-hematopoietic bone tumor. Survival of these young patients is related to the response to chemotherapy and development of metastases. Despite many advances in cancer research, chemotherapy regimens for osteosarcoma are still based on non-selective cytotoxic drugs. It is essential to investigate new specific molecular therapies for osteosarcoma to increase the survival rate of these patients. We performed exomic sequence analyses of 8 diagnostic biopsies of patients with conventional high grade osteosarcoma to advance our understanding of their genetic underpinnings and to correlate the genetic alteration with the clinical and pathological features of each patient to identify a personalized therapy. We identified 18,275 somatic variations in 8,247 genes and we found three mutated genes in 7/8 (87%) samples (KIF1B, NEB and KMT2C). KMT2C showed the highest number of variations; it is an important component of a histone H3 lysine 4 methyltransferase complex and it is one of the histone modifiers previously implicated in carcinogenesis, never studied in osteosarcoma. Moreover, we found a group of 15 genes that showed variations only in patients that did not respond to therapy and developed metastasis and some of these genes are involved in carcinogenesis and tumor progression in other tumors. These data could offer the opportunity to get a key molecular target to identify possible new strategies for early diagnosis and new therapeutic approaches for osteosarcoma and to provide a tailored treatment for each patient based on their genetic profile.

Targeting RNA for processing or destruction by the eukaryotic RNA exosome and its cofactors

In this review, Zinder and Lima highlight recent advances that have illuminated roles for the RNA exosome and its cofactors in specific biological pathways, alongside studies that attempted to dissect these activities through structural and biochemical characterization of nuclear and cytoplasmic RNA exosome complexes.

The eukaryotic RNA exosome is an essential and conserved protein complex that can degrade or process RNA substrates in the 3′-to-5′ direction. Since its discovery nearly two decades ago, studies have focused on determining how the exosome, along with associated cofactors, achieves the demanding task of targeting particular RNAs for degradation and/or processing in both the nucleus and cytoplasm. In this review, we highlight recent advances that have illuminated roles for the RNA exosome and its cofactors in specific biological pathways, alongside studies that attempted to dissect these activities through structural and biochemical characterization of nuclear and cytoplasmic RNA exosome complexes.