DNA Copy Number Variations as Markers of Mutagenic Impact

Genome-wide association study of copy number variations with shank traits in a F2 crossbred chicken population

Copy number variations (CNVs) are large-scale changes in the DNA sequence that can affect the genetic structure and phenotype of an organism. The purpose of this study was to investigate the existing CNVs and their associations with the shank diameter (ShD) and shank length (ShL) traits using data from an F2 crossbred chicken population. To carry out the study, 312 chickens were genotyped using the Illumina 60k SNP Beadchip. The shank traits of the birds were measured from day 1 to 12 weeks of age. penncnv and cnvruler tools were used to find copy numbers and regions with copy number changes (CNVR), respectively. The CNVRanger package was used to perform a genome-wide association study between shank traits and CNVs. Gene ontology research in CNVRs was carried out using the david database. In this investigation, 966 CNVs and 606 regions with copy number changes were discovered. The copy number states and variations were randomly distributed along the length of the autosomal chromosomes. Weeks 1-4, 9 and 12 of growth revealed a significant association of copy number variations with shank traits, false discovery rate (FDR-corrected p-value < 0.01), and the majority of CNVs that were statistically significant were found on chromosomes 1-3. These CNV segments are nearby genes such as KCNJ12, FGF6 and MYF5, which are fundamental to growth and development. In addition, gene set analyses revealed terms related to muscle physiology, regulation of cellular processes and potassium channels.

Contribution of Endothelial Dysfunction to Cancer Susceptibility and Progression: A Comprehensive Narrative Review on the Genetic Risk Component

Venous thromboembolism (VTE) is a challenging clinical obstacle in oncological settings, marked by elevated incidence rates and resulting morbidity and mortality. In the context of cancer-associated thrombosis (CAT), endothelial dysfunction (ED) plays a crucial role in promoting a pro-thrombotic environment as endothelial cells lose their ability to regulate blood flow and coagulation. Moreover, emerging research suggests that this disorder may not only contribute to CAT but also impact tumorigenesis itself. Indeed, a dysfunctional endothelium may promote resistance to therapy and favour tumour progression and dissemination. While extensive research has elucidated the multifaceted mechanisms of ED pathogenesis, the genetic component remains a focal point of investigation. This comprehensive narrative review thus delves into the genetic landscape of ED and its potential ramifications on cancer progression. A thorough examination of genetic variants, specifically polymorphisms, within key genes involved in ED pathogenesis, namely eNOS, EDN1, ACE, AGT, F2, SELP, SELE, VWF, ICAM1, and VCAM1, was conducted. Overall, these polymorphisms seem to play a context-dependent role, exerting both oncogenic and tumour suppressor effects depending on the tumour and other environmental factors. In-depth studies are needed to uncover the mechanisms connecting these DNA variations to the pathogenesis of malignant diseases.

Copy Number Variations in Pancreatic Cancer: From Biological Significance to Clinical Utility

Pancreatic ductal adenocarcinoma (PDAC) is the most common type of pancreatic cancer, characterized by high tumor heterogeneity and a poor prognosis. Inter- and intra-tumoral heterogeneity in PDAC is a major obstacle to effective PDAC treatment; therefore, it is highly desirable to explore the tumor heterogeneity and underlying mechanisms for the improvement of PDAC prognosis. Gene copy number variations (CNVs) are increasingly recognized as a common and heritable source of inter-individual variation in genomic sequence. In this review, we outline the origin, main characteristics, and pathological aspects of CNVs. We then describe the occurrence of CNVs in PDAC, including those that have been clearly shown to have a pathogenic role, and further highlight some key examples of their involvement in tumor development and progression. The ability to efficiently identify and analyze CNVs in tumor samples is important to support translational research and foster precision oncology, as copy number variants can be utilized to guide clinical decisions. We provide insights into understanding the CNV landscapes and the role of both somatic and germline CNVs in PDAC, which could lead to significant advances in diagnosis, prognosis, and treatment. Although there has been significant progress in this field, understanding the full contribution of CNVs to the genetic basis of PDAC will require further research, with more accurate CNV assays such as single-cell techniques and larger cohorts than have been performed to date.

Genotype and phenotype data standardization, utilization and integration in the big data era for agricultural sciences

Large-scale genotype and phenotype data have been increasingly generated to identify genetic markers, understand gene function and evolution and facilitate genomic selection. These datasets hold immense value for both current and future studies, as they are vital for crop breeding, yield improvement and overall agricultural sustainability. However, integrating these datasets from heterogeneous sources presents significant challenges and hinders their effective utilization. We established the Genotype-Phenotype Working Group in November 2021 as a part of the AgBioData Consortium (https://www.agbiodata.org) to review current data types and resources that support archiving, analysis and visualization of genotype and phenotype data to understand the needs and challenges of the plant genomic research community. For 2021-22, we identified different types of datasets and examined metadata annotations related to experimental design/methods/sample collection, etc. Furthermore, we thoroughly reviewed publicly funded repositories for raw and processed data as well as secondary databases and knowledgebases that enable the integration of heterogeneous data in the context of the genome browser, pathway networks and tissue-specific gene expression. Based on our survey, we recommend a need for (i) additional infrastructural support for archiving many new data types, (ii) development of community standards for data annotation and formatting, (iii) resources for biocuration and (iv) analysis and visualization tools to connect genotype data with phenotype data to enhance knowledge synthesis and to foster translational research. Although this paper only covers the data and resources relevant to the plant research community, we expect that similar issues and needs are shared by researchers working on animals. Database URL: https://www.agbiodata.org.

Copy number variation of bovine S100A7 as a positional candidate affected body measurements

Beef production is closely related to the national economy and the attention has been paid to the improvement of beef cattle by molecular markers associated. Copy number variations (CNVs) recently have been gained many researches and recognized as an important source of genetic variation. Extensive studies have indicated that CNVs have effects on a large range of economic traits by a wide range of gene copy number alteration. S100A7 is a member of S100 family which is a famous family of Ca2+-binding proteins. S100A7 plays a crucial role in many important phenotypes (progress) including inflammatory diseases, psoriasis, obesity, etc. The aim of our study was to explore the phenotypic effects of CNV located in the S100A7 gene of bovine chromosome 3. We detected S100A7 CNV by qPCR in different cattle breeds, including Qinchuan cattle, Yunling cattle, Xianan cattle and a crossbred group Pinan. The copy number was identified as gain, normal and loss type, our results showed that the gain type was the main type in three types of S100A7 CNV of the whole tested breeds. After CNV detection, association analysis between S100A7 CNV and growth traits was carried out in four cattle breeds. We found significant effects of the CNV on cattle growth traits with differently preferred CNV types such as gain type with better chest depth (p = 0.043) in QC, loss type with better body length (p = 0.008) and rump width (p = 0.014) in YL, normal with better chest girth (p = 0.001), gain with better waist width (p = 0.001) and rump width (p = 0.044) in PN. These results suggested that the S100A7 CNV could affect the phenotypic traits and be used as a promising genetic marker for cattle molecular breeding.

DNA cytosine deamination is associated with recurrent Somatic Copy Number Alterations in stomach adenocarcinoma

Stomach Adenocarcinoma (STAD) is a leading cause of death worldwide. Somatic Copy Number Alterations (SCNAs), which result in Homologous recombination (HR) deficiency in double-strand break repair, are associated with the progression of STAD. However, the landscape of frequent breakpoints of SCNAs (hotspots) and their functional impacts remain poorly understood. In this study, we aimed to explore the frequency and impact of these hotspots in 332 STAD patients and 1,043 cancer cells using data from the Cancer Genome Atlas (TCGA) and Cancer Cell Line Encyclopedia (CCLE). We studied the rates of DSB (Double-Strand Breaks) loci in STAD patients by employing the Non-Homogeneous Poisson Distribution (λ), based on which we identified 145 DSB-hotspots with genes affected. We further verified DNA cytosine deamination as a critical process underlying the burden of DSB in STAD. Finally, we illustrated the clinical impact of the significant biological processes. Our findings highlighted the relationship between DNA cytosine deamination and SCNA in cancer was associated with recurrent Somatic Copy Number Alterations in STAD.

Genomic Landscape of Copy Number Variations and Their Associations with Climatic Variables in the World's Sheep

Sheep show characteristics of phenotypic diversity and adaptation to diverse climatic regions. Previous studies indicated associations between copy number variations (CNVs) and climate-driven adaptive evolution in humans and other domestic animals. Here, we constructed a genomic landscape of CNVs (n = 39,145) in 47 old autochthonous populations genotyped at a set of high-density (600 K) SNPs to detect environment-driven signatures of CNVs using a multivariate regression model. We found 136 deletions and 52 duplications that were significantly (Padj. < 0.05) associated with climatic variables. These climate-mediated selective CNVs are involved in functional candidate genes for heat stress and cold climate adaptation (e.g., B3GNTL1, UBE2L3, and TRAF2), coat and wool-related traits (e.g., TMEM9, STRA6, RASGRP2, and PLA2G3), repairing damaged DNA (e.g., HTT), GTPase activity (e.g., COPG), fast metabolism (e.g., LMF2 and LPIN3), fertility and reproduction (e.g., SLC19A1 and CCDC155), growth-related traits (e.g., ADRM1 and IGFALS), and immune response (e.g., BEGAIN and RNF121) in sheep. In particular, we identified significant (Padj. < 0.05) associations between probes in deleted/duplicated CNVs and solar radiation. Enrichment analysis of the gene sets among all the CNVs revealed significant (Padj. < 0.05) enriched gene ontology terms and pathways related to functions such as nucleotide, protein complex, and GTPase activity. Additionally, we observed overlapping between the CNVs and 140 known sheep QTLs. Our findings imply that CNVs can serve as genomic markers for the selection of sheep adapted to specific climatic conditions.

Identification of Extremely Rare Pathogenic CNVs by Array CGH in Saudi Children with Developmental Delay, Congenital Malformations, and Intellectual Disability

Chromosomal imbalance is implicated in developmental delay (DD), congenital malformations (CM), and intellectual disability (ID), and, thus, precise identification of copy number variations (CNVs) is essential. We therefore aimed to investigate the genetic heterogeneity in Saudi children with DD/CM/ID. High-resolution array comparative genomic hybridization (array CGH) was used to detect disease-associated CNVs in 63 patients. Quantitative PCR was done to confirm the detected CNVs. Giemsa banding-based karyotyping was also performed. Array CGH identified chromosomal abnormalities in 24 patients; distinct pathogenic and/or variants of uncertain significance CNVs were found in 19 patients, and aneuploidy was found in 5 patients including 47,XXY (n = 2), 45,X (n = 2) and a patient with trisomy 18 who carried a balanced Robertsonian translocation. CNVs including 9p24p13, 16p13p11, 18p11 had gains/duplications and CNVs, including 3p23p14, 10q26, 11p15, 11q24q25, 13q21.1q32.1, 16p13.3p11.2, and 20q11.1q13.2, had losses/deletions only, while CNVs including 8q24, 11q12, 15q25q26, 16q21q23, and 22q11q13 were found with both gains or losses in different individuals. In contrast, standard karyotyping detected chromosomal abnormalities in ten patients. The diagnosis rate of array CGH (28%, 18/63 patients) was around two-fold higher than that of conventional karyotyping (15.87%, 10/63 patients). We herein report, for the first time, the extremely rare pathogenic CNVs in Saudi children with DD/CM/ID. The reported prevalence of CNVs in Saudi Arabia adds value to clinical cytogenetics.

Multi-omics profiling of papillary thyroid microcarcinoma reveals different somatic mutations and a unique transcriptomic signature

BACKGROUND

Papillary thyroid microcarcinoma (PTMC) incidence has significantly increased, and some cases still exhibit invasive traits. The entire molecular landscape of PTMC, which can offer hints for the etiology of cancer, is currently absent.

METHODS

We compared our findings with those for PTMC in the TCGA by analyzing the largest study at the current stage of whole exome sequencing and RNA-sequencing data from 64 patients with PTMC. Then, we systematically demonstrated the differences between the two PTMC subtypes based on multi-omics analyses. Additionally, we created a molecular prediction model for the PTMC subtypes and validated them among TCGA patients for individualized integrative assessment.

RESULTS

In addition to the presence of BRAF mutations and RET fusions in the TCGA cohort, we also discovered a new molecular signature named PTMC-inflammatory that implies a potential response to immune intervention, which is enriched with AFP mutations, IGH@-ext fusions, elevated immune-related genes, positive peroxidase antibody, and positive thyroglobulin antibody. Additionally, a molecular prediction model for the PTMC-inflammatory patients was created and validated among TCGA patients, while the prognosis for these patients is poor.

CONCLUSIONS

Our findings comprehensively define the clinical and molecular features of PTMC and may inspire new therapeutic hypotheses.

Molecular heterogeneity and commonalities in pancreatic cancer precursors with gastric and intestinal phenotype

OBJECTIVE

Due to the limited number of modifiable risk factors, secondary prevention strategies based on early diagnosis represent the preferred route to improve the prognosis of pancreatic ductal adenocarcinoma (PDAC). Here, we provide a comparative morphogenetic analysis of PDAC precursors aiming at dissecting the process of carcinogenesis and tackling the heterogeneity of preinvasive lesions.

DESIGN

Targeted and whole-genome low-coverage sequencing, genome-wide methylation and transcriptome analyses were applied on a final collective of 122 morphologically well-characterised low-grade and high-grade PDAC precursors, including intestinal and gastric intraductal papillary mucinous neoplasms (IPMN) and pancreatic intraepithelial neoplasias (PanIN).

RESULTS

Epigenetic regulation of mucin genes determines the phenotype of PDAC precursors. PanIN and gastric IPMN display a ductal molecular profile and numerous similarly regulated pathways, including the Notch pathway, but can be distinguished by recurrent deletions and differential methylation and, in part, by the expression of mucin-like 3. Intestinal IPMN are clearly distinct lesions at the molecular level with a more instable genotype and are possibly related to a different ductal cell compartment.

CONCLUSIONS

PDAC precursors with gastric and intestinal phenotype are heterogeneous in terms of morphology, genetic and epigenetic profile. This heterogeneity is related to a different cell identity and, possibly, to a different aetiology.

Esophageal Squamous Cancer from 4NQO-Induced Mice Model: CNV Alterations

Squamous esophageal carcinoma is a common pathological type of esophageal carcinoma around the world. The prognosis of esophageal carcinoma is usually poor and diagnosed at late stages. Recently, research suggested that genomic instability occurred in esophageal cells during the development of esophageal squamous cell carcinoma (ESCC). Identifying prognostic and specific genomic characteristics, especially at the early hyperplasia stage, is critical. Mice were given 4-nitroquinoline 1-oxide (4NQO) with drinking water to induce esophageal cancer. The immortalized human esophageal epithelial cell line (NE2) was also treated with 4NQO. We performed histologic analyses, immunofluorescence, and immunohistochemical staining to detect DNA damage at different time points. Whole-exome sequencing was accomplished on the esophagus tissues at different pathological stages to detect single-nucleotide variants and copy number variation (CNV) in the genome. Our findings indicate that all mice were tumor-forming, and a series of changes from simple hyperplasia (ESSH) to intraepithelial neoplasia (IEN) to esophageal squamous cell carcinoma (ESCC) was seen at different times. The expression of γ-H2AX increased from ESSH to ESCC. In addition, mutations of the Muc4 gene were detected throughout the pathological stages. Furthermore, CNV burden appeared in the esophageal tissues from the beginning of ESSH and accumulated more in cancer with the deepening of the lesions. This study demonstrates that mutations caused by the early appearance of DNA damage may appear in the early stage of malignant tissue before the emergence of atypia. The detection of CNV and mutations of the Muc4 gene may be used as an ultra-early screening indicator for esophageal cancer.

Comprehensive Multiomic Analysis Identified TUBA1C as a Potential Prognostic Biological Marker of Immune-Related Therapy in Pan-Cancer

TUBA1C is correlated with an unfavourable prognosis and the infiltration of immune cells in several cancers. However, its function as a significant biomarker for the prognosis of immunotherapy in pan-cancer remains unclear. This study aims at assessing the role of TUBA1C in pan-cancer at multiple levels, including mutations, gene expression, methylation, m6A methylation, and immune cell infiltration levels. Data retrieved from major public databases, such as TCGA, GEO, GTEx, GSCA, CancerSEA, HPA, and RNAactDrugs, revealed that TUBA1C expression was high in 33 cancer types. Survival analysis revealed that TUBA1C was a poor prognostic factor for 12 tumour types, and mutations, CNVs, and methylation affected the prognosis of some cancer types. Furthermore, TUBA1C was found to be related to immune-related genes, immune cell infiltration, and the immune microenvironment. In addition, the sensitivity of 10 anticancer drugs was associated with high TUBA1C expression. Therefore, TUBA1C may serve as a viable prognostic biomarker for immunotherapy of pan-cancer.

Molecular characterization, clinical relevance and immune feature of m7G regulator genes across 33 cancer types

Over 170 RNA modifications have been identified after transcriptions, involving in regulation of RNA splicing, processing, translation and decay. Growing evidence has unmasked the crucial role of N6-methyladenosine (m6A) in cancer development and progression, while, as a relative newly found RNA modification, N7-methylguanosine (m7G) is also certified to participate in tumorigenesis via different catalytic machinery from that of m6A. However, system analysis on m7G RNA modification-related regulator genes is lack. In this study, we first investigated the genetic alteration of m7G related regulator genes in 33 cancers, and found mRNA expression levels of most regulator genes were positively correlated with copy number variation (CNV) and negatively correlated with methylation in most cancers. We built a m7G RNA modification model based on the enrichment of the regulator gene scores to evaluate the m7G modification levels in 33 cancers, and investigated the connections of m7G scores to clinical outcomes. Furthermore, we paid close attention to the role of m7G in immunology due to the widely used immune checkpoint blockade therapy. Our results showed the higher m7G scores related to immunosuppression of tumor cells. Further confirmation with phase 3 clinical data with application of anti-PDL1/PDL indicated the impact of m7G modification level on immunotherapy effect. Relevance of m7G regulator genes and drug sensitivity was also evaluated to provide a better treatment choice when treating cancers. In summary, our study uncovered the profile of m7G RNA modification through various cancers, and figured out the connection of m7G modification levels with therapeutical outcomes, providing potential better options of cancer treatment.

Comprehensive Analysis of m6A RNA Methylation Regulators and the Immune Microenvironment to Aid Immunotherapy in Pancreatic Cancer

Pancreatic cancer (PAAD) is one of the most malignant cancers and immune microenvironment has been proved to be involved in pathogenesis of PAAD. m6A modification, related to the expression of m6A regulators, participates in the development of multiple cancers. However, the correlation between m6A regulators and immune microenvironment was largely unknown in PAAD. And because of the small sample size of pancreatic cancer in the TCGA database, it is not enough to draw a convincing conclusion. In the present study, we downloaded seven pancreatic cancer datasets with survival data and removed batch effects among these datasets to be used as the PAAD cohort to analyze the immune landscape of PAAD and the expression pattern of m6A regulators and divided the integrated dataset into cluster 1 and cluster 2 by consensus clustering for m6A regulators. Lower m6A regulators were found to be related to higher immune cell infiltration and a better survival. Moreover, we identified six m6A regulators and constructed the prognostic signature of m6A regulators. Patients with low-risk score had a higher response to immune checkpoint inhibitor and a longer overall survival. To figure out the underlying mechanism, we analyzed the cancer immunity cycle, most altered genes, gene set enrichment analysis (GSEA) and gene set variation analysis (GSVA) in risk subtypes. In summary, the present study proved m6A regulators modulated the PAAD immune microenvironment. And risk scores served as predictive indicator for immunotherapy and played a prognostic role for PAAD patients. Our study provided novel therapeutic targets to improve immunotherapy efficacy.

The Mutation and Expression Level of LRP1B are Associated with Immune Infiltration and Prognosis in Hepatocellular Carcinoma

Purpose

This study aimed to explore the expression level and mutation of LRP1B in hepatocellular carcinoma (HCC) and to analyse the relationship between its prognostic value and immune invasion.

Methods

HCC mutant gene sets were obtained from the Cancer Genome Atlas and International Cancer Genome Consortium databases. The Kaplan–Meier method was used to evaluate the prognostic value of LRP1B expression and mutation load in HCC. The relationships between LRP1B expression level and immune cells and immune marker molecules were analysed by using the TIMER database. The association of LRP1B expression with drug sensitivity was obtained by using CellMiner. Gene set enrichment analysis and co-expression by Spearman correlation analysis were used to explore the internal mechanism of LRP1B in HCC.

Results

Seventeen most commonly mutated genes were screened out, and LRP1B was the only gene associated with HCC prognosis. The copy number variations were significantly correlated with T cell CD8⁺ (P < 0.05). LRP1B expression level was positively correlated with the infiltration degree of macrophage (P < 0.05, R = 0.132), myeloid dendritic cell (P < 0.05, R = 0.093), neutrophil (P < 0.05, R = 0.134) and T cell CD8+ cells (P < 0.05, R = 0.102) and negatively correlated with B cell (P < 0.05, R = −0.014) and T cell CD4+ (P < 0.05, R = −0.075). LRP1B expression level was significantly correlated with immunomarker molecules and drug sensitivity (all P < 0.05). The prediction of lncRNA RUSC1-AS1/hsa-miR-215-5p/LRP1B axis by bioinformatics may be the potential mechanism underlying LRP1B’s effect on HCC prognosis and progression.

Conclusion

LRP1B plays a vital role in HCC prognostic value, which is expected to be a new potential therapeutic target for HCC. LRP1B provides a theoretical basis for the clinical targeted therapy of HCC.

Copy number variation: Characteristics, evolutionary and pathological aspects

Copy number variants (CNVs) were the subject of extensive research in the past years. They are common features of the human genome that play an important role in evolution, contribute to population diversity, development of certain diseases, and influence host–microbiome interactions. CNVs have found application in the molecular diagnosis of many diseases and in non-invasive prenatal care, but their full potential is only emerging. CNVs are expected to have a tremendous impact on screening, diagnosis, prognosis, and monitoring of several disorders, including cancer and cardiovascular disease. Here, we comprehensively review basic definitions of the term CNV, outline mechanisms and factors involved in CNV formation, and discuss their evolutionary and pathological aspects. We suggest a need for better defined distinguishing criteria and boundaries between known types of CNVs.

Catalytically inactive Cas9 impairs DNA replication fork progression to induce focal genomic instability

Catalytically inactive Cas9 (dCas9) has become an increasingly popular tool for targeted gene activation/inactivation, live-cell imaging, and base editing. While dCas9 was reported to induce base substitutions and indels, it has not been associated with structural variations. Here, we show that dCas9 impedes replication fork progression to destabilize tandem repeats in budding yeast. When targeted to the CUP1 array comprising ∼16 repeat units, dCas9 induced its contraction in most cells, especially in the presence of nicotinamide. Replication intermediate analysis demonstrated replication fork stalling in the vicinity of dCas9-bound sites. Genetic analysis indicated that while destabilization is counteracted by the replisome progression complex components Ctf4 and Mrc1 and the accessory helicase Rrm3, it involves single-strand annealing by the recombination proteins Rad52 and Rad59. Although dCas9-mediated replication fork stalling is a potential risk in conventional applications, it may serve as a novel tool for both mechanistic studies and manipulation of genomic instability.

Development and validation of the trans-omics model for pancreatic adenocarcinoma

Aim: To develop a trans-omics-based molecular clinicopathological algorithm for predicting pancreatic adenocarcinoma prognosis, we performed a comprehensive analysis of the expression levels of mRNA, DNA methylation and DNA copy number in The Cancer Genome Atlas dataset. Materials & methods: Based on the least absolute shrinkage and selection operator method - COX regression analysis, a trans-omics-based classifier was established to predict overall survival. Nomogram was constructed by combining the classifier band clinical pathological characterization. Results: Based on trans-omics, we developed a 10-gene-based classifier and a molecular-clinicopathologic nomogram for predicting overall survival with satisfactory accuracy. Conclusion: Trans-omics-based classifier and molecule-clinicopathological nomogram based on the classifier can accurately predict the prognosis of pancreatic adenocarcinoma patients.

TERT-Regulation and Roles in Cancer Formation

Telomerase reverse transcriptase (TERT) is a catalytic subunit of telomerase. Telomerase complex plays a key role in cancer formation by telomere dependent or independent mechanisms. Telomere maintenance mechanisms include complex TERT changes such as gene amplifications, TERT structural variants, TERT promoter germline and somatic mutations, TERT epigenetic changes, and alternative lengthening of telomere. All of them are cancer specific at tissue histotype and at single cell level. TERT expression is regulated in tumors via multiple genetic and epigenetic alterations which affect telomerase activity. Telomerase activity via TERT expression has an impact on telomere length and can be a useful marker in diagnosis and prognosis of various cancers and a new therapy approach. In this review we want to highlight the main roles of TERT in different mechanisms of cancer development and regulation.

The p53-inducible CLDN7 regulates colorectal tumorigenesis and has prognostic significance

Most colorectal cancer (CRC) are characterized by allele loss of the genes located on the short arm of chromosome 17 (17p13.1), including the tumor suppressor p53 gene. Although important, p53 is not the only driver of chromosome 17p loss. In this study, we explored the biological and prognostic significance of genes around p53 on 17p13.1 in CRC. The Cancer Genome Atlas (TCGA) were used to identify differentially expressed genes located between 1000 kb upstream and downstream of p53 gene. The function of CLDN7 was evaluated by both in vitro and in vivo experiments. Quantitative real-time PCR, western blot, and promoter luciferase activity, immunohistochemistry were used to explore the molecular drivers responsible for the development and progression of CRC. The results showed that CLDN7, located between 1000 kb upstream and downstream of p53 gene, were remarkably differentially expressed in tumor and normal tissues. CLDN7 expression also positively associated with p53 level in different stages of the adenoma-carcinoma sequence. Both in vitro and in vivo assays showed that CLDN7 inhibited cell proliferation in p53 wild type CRC cells, but had no effects on p53 mutant CRC cells. Mechanistically, p53 could bind to CLDN7 promoter region and regulate its expression. Clinically, high CLDN7 expression was negatively correlated with tumor size, invasion depth, lymphatic metastasis and AJCC III/IV stage, but was positively associated with favorable prognosis of CRC patients. Collectively, our work uncovers the tumor suppressive function for CLDN7 in a p53-dependent manner, which may mediate colorectal tumorigenesis induced by p53 deletion or mutation.