Transcript map of the 8p23 putative tumor suppressor region

Tumor suppressor role of the complement inhibitor CSMD1 and its role in TNF-induced neuroinflammation in gliomas

BACKGROUND

The complement inhibitor CSMD1 acts as a tumor suppressor in various types of solid cancers. Despite its high level of expression in the brain, its function in gliomas, malignant brain tumors originating from glial cells, has not been investigated.

METHODS

Three cohorts of glioma patients comprising 1500 patients were analyzed in our study along with their clinical data. H4, U-118 and U-87 cell lines were used to investigate the tumor suppressor function of CSMD1 in gliomas. PDGFB-induced brain tumor model was utilized for the validation of in vitro data.

RESULTS

The downregulation of CSMD1 expression correlated with reduced overall and disease-free survival, elevated tumor grade, wild-type IDH genotype, and intact 1p/19q status. Moreover, enhanced activity was noted in the neuroinflammation pathway. Importantly, ectopic expression of CSMD1 in glioma cell lines led to decreased aggressiveness in vitro. Mechanically, CSMD1 obstructed the TNF-induced NF-kB and STAT3 signaling pathways, effectively suppressing the secretion of IL-6 and IL-8. There was also reduced survival in PDGFB-induced brain tumors in mice when Csmd1 was downregulated.

CONCLUSIONS

Our study has identified CSMD1 as a tumor suppressor in gliomas and elucidated its role in TNF-induced neuroinflammation, contributing to a deeper understanding of glioma pathogenesis.

The Diverse Role of CUB and Sushi Multiple Domains 1 (CSMD1) in Human Diseases

CUB and Sushi Multiple Domains 1 (CSMD1), a tumour suppressor gene, encodes a large membrane-bound protein including a single transmembrane domain. This transmembrane region has a potential tyrosine phosphorylation site, suggesting that CSMD1 is involved in controlling cellular functions. Although the specific mechanisms of action for CSMD1 have not yet been uncovered, it has been linked to a number of processes including development, complement control, neurodevelopment, and cancer progression. In this review, we summarise CSMD1 functions in the cellular processes involved in the complement system, metastasis, and Epithelial mesenchymal transition (EMT) and also in the diseases schizophrenia, Parkinson's disease, and cancer. Clarifying the association between CSMD1 and the aforementioned diseases will contribute to the development of new diagnosis and treatment methods for these diseases. Recent studies in certain cancer types, e.g., gastric cancer, oesophageal cancer, and head and neck squamous cell carcinomas, have indicated the involvement of CSMD1 in response to immunotherapy.

Exploring the Genetic Association between Obesity and Serum Lipid Levels Using Bivariate Methods

It is crucial to understand the genetic mechanisms and biological pathways underlying the relationship between obesity and serum lipid levels. Structural equation models (SEMs) were constructed to calculate heritability for body mass index (BMI), total cholesterol (TC), triglyceride (TG), high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), and the genetic connections between BMI and the four classes of lipids using 1197 pairs of twins from the Chinese National Twin Registry (CNTR). Bivariate genomewide association studies (GWAS) were performed to identify genetic variants associated with BMI and lipids using the records of 457 individuals, and the results were further validated in 289 individuals. The genetic background affecting BMI may differ by gender, and the heritability of males and females was 71% (95% CI [.66, .75]) and 39% (95% CI [.15, .71]) respectively. BMI was positively correlated with TC, TG and LDL-C in phenotypic and genetic correlation, while negatively correlated with HDL-C. There were gender differences in the correlation between BMI and lipids. Bivariate GWAS analysis and validation stage found 7 genes (LOC105378740, LINC02506, CSMD1, MELK, FAM81A, ERAL1 and MIR144) that were possibly related to BMI and lipid levels. The significant biological pathways were the regulation of cholesterol reverse transport and the regulation of high-density lipoprotein particle clearance (p < .001). BMI and blood lipid levels were affected by genetic factors, and they were genetically correlated. There might be gender differences in their genetic correlation. Bivariate GWAS analysis found MIR144 gene and its related biological pathways may influence obesity and lipid levels.

CSMD1 Shows Complex Patterns of Somatic Copy Number Alterations and Expressions of mRNAs and Target Micro RNAs in Esophageal Squamous Cell Carcinoma

Simple Summary

Human Cub and Sushi Multiple Domains 1 (CSMD1) is a novel candidate tumor-suppressor gene. We investigated CSMD1 in esophageal squamous cell carcinoma (ESCC) by performing an integrated analysis of somatic DNA alterations (i.e., copy number alteration, allelic imbalance, and loss of heterozygosity) with RNA expressions (mRNA and target miRNAs) on specimens from the same ESCC patients, using data from SNP, miRNA, and RT-PCR arrays. Our results indicate that the CSMD1 gene may play a role in the development of ESCC through complex patterns involving somatic alterations and mRNA expression. Furthermore, somatic copy number alterations in SNPs located in non-coding regions of CSMD1 appear to influence expression of both this gene and its target miRNAs.

Abstract

Background: Human Cub and Sushi Multiple Domains 1 (CSMD1) is a novel candidate tumor-suppressor gene that codes for multiple domains, including complement regulatory and adhesion proteins, and has recently been shown to have alterations in multiple cancers. We investigated CSMD1 in esophageal squamous cell carcinoma (ESCC) by performing an integrated analysis on somatic copy number alterations (CNAs), including copy-number gain or loss, allelic imbalance (AI), loss of heterozygosity (LOH), and the expressions of mRNA and its target miRNAs on specimens from the same patients with ESCC. Results: (i) Two-thirds of ESCC patients had all three types of alterations studied—somatic DNA alterations in 70%, and abnormal expressions of CSMD1 RNA in 69% and in target miRNAs in 66%; patterns among these alterations were complex. (ii) In total, 97% of 888 CSMD1 SNPs studied showed somatic DNA alterations, with most located near exons 4–11, 24–25, 39–40, 55–56, and 69–70. (iii) In total, 68% of SNPs with a CNA were correlated with expression of CSMD1. (iv) A total of 33 correlations between non-coding SNPs and expression of CSMD1 target miRs were found. Conclusions: Our results indicate that the CSMD1 gene may play a role in ESCC through complex patterns of DNA alterations and RNA and miRNA expressions. Alterations in some somatic SNPs in non-coding regions of CSMD1 appear to influence expression of this gene and its target miRNAs.

CSMD1 suppresses cancer progression by inhibiting proliferation, epithelial-mesenchymal transition, chemotherapy-resistance and inducing immunosuppression in esophageal squamous cell carcinoma

Human CUB and Sushi multiple domains (CSMD1) is considered a crucial role in cancer progression, but the specific function in esophageal squamous cell carcinoma (ESCC) is not clear. Understanding the role of CSMD1 in ESCC progression may lead to a novel strategy for ESCC treatment. Here, we found that both CSMD1 mRNA and protein levels were downregulated in ESCC tissues. Reduced CSMD1 expression was correlated with a poor prognosis in ESCC patients. CSMD1 expression inhibited proliferation, migration and invasion in ESCC cell lines in vitro. CSMD1 deficiency in established xenografted tumors increases tumor size and weight. We further found that CSMD1-overexpression cells are more sensitive to chemotherapy. Moreover, we addressed the role of CSMD1 in the CD8⁺ T cell immune response. An in vitro killing assay showed that the cytotoxicity of CD8⁺ T cells was inhibited in CSMD1-overexpression tumor cells. In vivo, in CSMD1 deficiency tumor-bearing mice activation and expansion of CD8⁺ T cells were increased. Further investigation showed that CSMD1 expression on tumor cells was positively correlated with CD8⁺ T cells infiltration and cytokines secretion. These findings highlight that CSMD1 is a tumor suppressor gene in ESCC patients and a positive regulator of CD8⁺ T cells expansion and activation, and could increase cytokines secretion, indicating that tumor cell-associated CSMD1 might be a target for ESCC.

New Insight into the human genetic diversity in North African populations by genotyping of SNPs in DRD3 , CSMD1 and NRG1 genes

Background

The single nucleotide polymorphisms (SNPs) of the dopamine D3 receptor (DRD3), the CUB and sushi multiple domains 1 (CSMD1) and the neuregulin 1 (NRG1) genes were used to study the genetic diversity and affinity among North African populations and to examine their genetic relationships in worldwide populations.

Methods

The rs3773678, rs3732783 and rs6280 SNPs of the DRD3 gene located on chromosome 3, the rs10108270 SNP of the CSMD1 gene and the rs383632, rs385396 and rs1462906 SNPs of the NRG1 gene located on chromosome 8 were analysed in 366 individuals from seven North African populations (Libya, Kairouan, Mehdia, Sousse, Kesra, Smar and Kerkennah).

Results

The low values of F_ST indicated that only 0.27%–1.65% of the genetic variability was due to the differences between the populations. The Kairouan population has the lowest average heterozygosity among the North African populations. Haplotypes composed of the ancestral alleles ACC and ACAT were more frequent in the Kairouan population than in other North African populations. The PCA and the haplotypic analysis showed that the genetic structure of populations in North Africa was closer to that of Europeans, Admixed Americans, South Asians and East Asians. However, analysis of the rs3732783 and rs6280 SNPs revealed that the CT microhaplotype was specific to the North African population.

Conclusions

The Kairouan population exhibited a relatively low rate of genetic variability. The North African population has undergone significant gene flow but also evolutionary forces that have made it genetically distinct from other populations.

Association of CSMD1 with Tumor Mutation Burden and Other Clinical Outcomes in Gastric Cancer

Background

Immunotherapy is considered as a powerful and promising clinical approach for the treatment of gastric cancer (GC). However, it is still challenging to precisely screen patients who potentially benefit from immune checkpoint therapy (ICT). Identification of potential biomarkers for selecting patients sensitive to immunotherapy was urgently needed.

Methods

Public sequence data and corresponding clinical data were used to explore the potential biomarkers for immunotherapy.

Results

We found that CSMD1 is the most frequently mutated gene and its mutation is highly correlated with prognosis in gastric cancer patients. Interestingly, patients with mutated CSMD1 exhibit a high mutation burden and upregulated PDL1 expression. The ratio of microsatellite instability (MSI) in the CSMD1 mutation cohort was higher than that in the cohort without CSMD1 mutation. Furthermore, patients with CSMD1 mutation have been found to possess a higher number of activated CD4+ T cells and neoantigens.

Conclusion

CSMD1 mutation may act as a novel biomarker for assessing the survival and immune therapy response in patients with gastric cancer.

CSMD1 Mutation Related to Immunity Can Be Used as a Marker to Evaluate the Clinical Therapeutic Effect and Prognosis of Patients with Esophageal Cancer

INTRODUCTION

As a highly aggressive tumor with a poor prognosis, esophageal cancer (ESCA)'s relationship with gene mutations is unclear. Therefore, we tried to explore the role of gene mutation in ESCA progression and its relationship with immune response, clinical treatment, and prognosis.

METHODS

In addition to copy number variation (CNV) situations of common genes obtained from 2 public databases, the relationship between mutations and prognosis/tumor mutational burden (TMB) was also analyzed. Kaplan-Meier survival and Cox regression analysis were used to identify the CSMD1 mutation status as an independent predictor of prognosis. We also enriched related functions and pathways. Next, the relationship between 22 immune cells and CSMD1 mutation status was analyzed. In addition to the differences in the expression levels of immune checkpoint inhibitors (ICIs)-related genes between the high TMB and low TMB groups, the differences in the expression levels of ICIs/m6a/multi-drug resistance-related genes and the sensitivity of three chemotherapeutic drugs between CSMD1 mutant and the wild group were also compared. In addition to differences and prognostic analysis of CSMD1 expression, the correlation analysis between the expression of these genes/immune cells and the expression of CSMD1 was also performed. Finally, a nomogram that could efficiently and conveniently predict the survival probability of ESCA patients was constructed and verified.

RESULTS

We obtained 17 frequently mutated genes distribution. Mutation and loss of CSMD1 are frequent in ESCA. Only CSMD1 mutation can be used as an independent predictor of poor prognosis. Patients in the high TMB group have a lower survival probability. Wild CSMD1 may be involved in immune-related pathways. More helper T cells and fewer resting state dendritic cells were found in the CSMD1 mutant group. The PD-1 expression in the high TMB group showed higher. Paclitaxel sensitivity and ABCC1 expression were higher in the wild CSMD1 group. Most cancers show differential expression of CSMD1. Except for the prognosis of ESCA, the expression of CSMD1 is related to immune cell content and the expression of ICIs/m6a/multi-drug resistance related genes.

DISCUSSION

CSMD1 mutation could be used as an immune-related biomarker to predict prognosis and treatment effect of paclitaxel. Mutation and loss of CSMD1 may promote the progression of ESCA.

Insilico Functional Analysis of Genome-Wide Dataset From 17,000 Individuals Identifies Candidate Malaria Resistance Genes Enriched in Malaria Pathogenic Pathways

Recent genome-wide association studies (GWASs) of severe malaria have identified several association variants. However, much about the underlying biological functions are yet to be discovered. Here, we systematically predicted plausible candidate genes and pathways from functional analysis of severe malaria resistance GWAS summary statistics (N = 17,000) meta-analysed across 11 populations in malaria endemic regions. We applied positional mapping, expression quantitative trait locus (eQTL), chromatin interaction mapping, and gene-based association analyses to identify candidate severe malaria resistance genes. We further applied rare variant analysis to raw GWAS datasets (N = 11,000) of three malaria endemic populations including Kenya, Malawi, and Gambia and performed various population genetic structures of the identified genes in the three populations and global populations. We performed network and pathway analyses to investigate their shared biological functions. Our functional mapping analysis identified 57 genes located in the known malaria genomic loci, while our gene-based GWAS analysis identified additional 125 genes across the genome. The identified genes were significantly enriched in malaria pathogenic pathways including multiple overlapping pathways in erythrocyte-related functions, blood coagulations, ion channels, adhesion molecules, membrane signalling elements, and neuronal systems. Our population genetic analysis revealed that the minor allele frequencies (MAF) of the single nucleotide polymorphisms (SNPs) residing in the identified genes are generally higher in the three malaria endemic populations compared to global populations. Overall, our results suggest that severe malaria resistance trait is attributed to multiple genes, highlighting the possibility of harnessing new malaria therapeutics that can simultaneously target multiple malaria protective host molecular pathways.

Germline variants in DNA repair genes are associated with young-onset head and neck cancer

The genetic predisposition to head and neck carcinomas (HNSCC) and how the known risk factors (papillomavirus infection, alcohol, and tobacco consumption) contribute to the early-onset disease are barely explored. Although HNSCC at early onset is rare, its frequency is increasing in recent years. Germline and somatic variants were assessed to build a comprehensive genetic influence pattern in HNSCC predisposition and patient outcome. Whole-exome sequencing was performed in 45 oral and oropharynx carcinomas paired with normal samples of young adults (≤49 years). We found FANCG, CDKN2A, and TPP germline variants previously associated with HNSCC risk. At least one germline variant in DNA repair pathway genes was detected in 67% of cases. Germline and somatic variants (including copy number variations) in FAT1 gene were identified in 9 patients (20%) and 12 tumors (30%), respectively. Somatic variants were found in HNSCC associated genes, such as TP53, CDKN2A, and PIK3CA. To date, 55 of 521 cases from the large cohort of TCGA presented < 49 years old. A comparison between the somatic alterations of TCGA-HNSCC at early onset and our dataset revealed strong similarities. Protein-protein interaction analysis between somatic and germline altered genes revealed a central role of TP53. Altogether, germline alterations in DNA repair genes potentially contribute to an increased risk of developing HNSCC at early-onset, while FAT1 could impact the prognosis.

Complement inhibitor CSMD1 modulates epidermal growth factor receptor oncogenic signaling and sensitizes breast cancer cells to chemotherapy

BACKGROUND

Human CUB and Sushi multiple domains 1 (CSMD1) is a large membrane-bound tumor suppressor in breast cancer. The current study aimed to elucidate the molecular mechanism underlying the effect of CSMD1 in highly invasive triple negative breast cancer (TNBC).

METHODS

We examined the antitumor action of CSMD1 in three TNBC cell lines overexpressing CSMD1, MDA-MB-231, BT-20 and MDA-MB-486, in vitro using scanning electron microscopy, proteome array, qRT-PCR, immunoblotting, proximity ligation assay, ELISA, co-immunoprecipitation, immunofluorescence, tumorsphere formation assays and flow cytometric analysis. The mRNA expression pattern and clinical relevance of CSMD1 were evaluated in 3520 breast cancers from a modern population-based cohort.

RESULTS

CSMD1-expressing cells had distinct morphology, with reduced deposition of extracellular matrix components. We found altered expression of several cancer-related molecules, as well as diminished expression of signaling receptors including Epidermal Growth Factor Receptor (EGFR), in CSMD1-expressing cells compared to control cells. A direct interaction of CSMD1 and EGFR was identified, with the EGF-EGFR induced signaling cascade impeded in the presence of CSMD1. Accordingly, we detected increased  ubiquitination levels of EGFR upon activation in CSMD1-expressing cells, as well as increased degradation kinetics and chemosensitivity. Accordingly, CSMD1 expression rendered tumorspheres pretreated with gefitinib more sensitive to chemotherapy. In addition, higher mRNA levels of CSMD1 tend to be associated with better outcome of triple negative breast cancer patients treated with chemotherapy.

CONCLUSIONS

Our results indicate that CSMD1 cross-talks with the EGFR endosomal trafficking cascade in a way that renders highly invasive breast cancer cells sensitive to chemotherapy. Our study unravels one possible underlying molecular mechanism of CSMD1 tumor suppressor function and may provide novel avenues for design of better treatment.

Inhibition of CUB and sushi multiple domains 1 (CSMD1) expression by miRNA-190a-3p enhances hypertrophic scar-derived fibroblast migration in vitro

Background

Hypertrophic scar (HTS) is a fibroproliferative skin disorder characterized by excessive cell proliferation, migration, and extracellular matrix (ECM) deposition. The CUB and Sushi multiple domains 1 (CSMD1) has previously been identified as the key regulatory gene of hypertrophic scar by a large sample GWAS study. However, further research has not yet been conducted to verify this finding in other HTS patients and to determine the underlying mechanism.

Results

In this study, we verified that CSMD1 was downregulated in both HTS tissue and HTS-derived fibroblasts. The knockdown of CSMD1 resulted in enhanced migration and fibronectin1 (FN1) secretion in fibroblasts in vitro. In addition, the upstream and downstream regulatory mechanisms of CSMD1 were also investigated through microRNA (miRNA) databases screening and RNA-sequencing (RNA-seq) respectively. The screening of four common microRNA (miRNA) databases suggested that miR-190a-3p binds to the CSMD1 and may regulate its expression. We confirmed that miR-190a-3p directly targeted the CSMD1–3′-UTR using luciferase reporter assays. Furthermore, the overexpression of miR-190a-3p showed promotion of migratory activity and FN1 secretion in fibroblasts, resembling the effect of CSMD1 knockdown; whereas the knockdown of miR-190a-3p exerted the opposite effect. Finally, transcriptomic analysis showed activation of Janus kinase-signal transducer and activator of transcription (JAK/STAT) signaling pathway in the CSMD1 knockdown fibroblasts.

Conclusions

This study has validated the conclusions of the previous GWAS study conducted in Chinese population. In vitro experiments have provided further evidence on the function of CSMD1 in the development of HTS, and have also revealed the underlying upstream and downstream regulating mechanisms. Additionally, the JAK/STAT signaling pathway identified using RNA-seq might provide a potential treatment approach, especially for HTS.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07920-8.

CUB and Sushi Multiple Domains (CSMD1) Gene Polymorphisms and Susceptibilities to Idiopathic Parkinson's Disease in Northern Chinese Han Population: A Case-Control Study

Evidence has shown that the CUB and Sushi Multiple Domains (CSMD1) gene is an inhibitor of the complement activation pathway and is also involved in central nervous system inflammation. Previous studies have revealed that the CSMD1 gene is related to familial Parkinson's disease. This study aimed to investigate the relationship between CSMD1 gene and susceptibility to Parkinson's disease in population of northern China. A case-control study was performed on 423 Parkinson's disease patients and 465 healthy controls matched for age and sex. DNA from enrolled subjects were extracted from the peripheral blood, and single nucleotide polymorphisms (SNPs) rs12681349 (C＞T), rs10503253 (C＞A), and rs1983474 (T＞G) within CSMD1 gene were genotyped using polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP). Genotype frequency of rs10503253 (CA versus CC : OR = 1.554, 95% CI = 1.169–2.066, p=0.002) and rs1983474 (GG versus TT : OR = 0.599, 95% CI = 0.401–0.895, p=0.012) was significantly different between PD cases and controls, but not for rs12681349. Comprehensive and subgroup analysis indicated that rs10503252 showed significant statistical differences in the dominant model (AA + CA versus CC : OR = 0.677, 95% CI = 0.517–0.886, p=0.004), late-onset cohort (CA versus CC : OR = 1.570, 95% CI = 1.159–2.126, p=0.004), and the female cohort (CA versus CC : OR = 0.687, 95% CI = 0.497–0.952, p=0.023), compared with the matched control group. The difference of recessive model of rs1983474 (GG versus TT + TG : OR = 1.837, 95% CI = 1.287–2.620, p=0.001) was significant in Parkinson's disease. According to the subgroup analysis, results indicated that late-onset cohort (GG versus TT : OR = 0.643, 95% CI = 0.420–0.985, p=0.042), male cohort (TG versus TT : OR = 2.160, 95% CI = 1.162–4.016, p=0.015), and female group (GG versus TT : OR = 0.418, 95% CI = 0.234–0.746, p=0.003) of rs1983474 were significantly associated with Parkinson's disease susceptibility. In both genotype and subgroup analysis, we failed to find any relationship between rs12681349 polymorphism and Parkinson's disease risk. Our results indicate that the rs10503253 and rs1983474 gene polymorphism may be associated with idiopathic Parkinson's disease susceptibility in Chinese population. Nevertheless, these conclusions need to be further verified by more studies.

The Role of Csmd1 during Mammary Gland Development

The Cub Sushi Multiple Domains-1 (CSMD1) protein is a tumour suppressor which has been shown to play a role in regulating human mammary duct development in vitro. CSMD1 knockdown in vitro demonstrated increased cell proliferation, invasion and motility. However, the role of Csmd1 in vivo is poorly characterised when it comes to ductal development and is therefore an area which warrants further exploration. In this study a Csmd1 knockout (KO) mouse model was used to identify the role of Csmd1 in regulating mammary gland development during puberty. Changes in duct development and protein expression patterns were analysed by immunohistochemistry. This study identified increased ductal development during the early stages of puberty in the KO mice, characterised by increased ductal area and terminal end bud number at 6 weeks. Furthermore, increased expression of various proteins (Stat1, Fak, Akt, Slug/Snail and Progesterone receptor) was shown at 4 weeks in the KO mice, followed by lower expression levels from 6 weeks in the KO mice compared to the wild type mice. This study identifies a novel role for Csmd1 in mammary gland development, with Csmd1 KO causing significantly more rapid mammary gland development, suggesting an earlier adult mammary gland formation.

Analyses of multi-omics differences between patients with high and low PD1/PDL1 expression in lung squamous cell carcinoma

BACKGROUND

Immunotherapy has achieved excellent results in patients with lung squamous cell carcinoma. However, in which population it can exert the greatest effect is still unknown. Some studies have suggested that its effect is related to the expression level of PD1. Analyzing the relationship between PD1 expression level and genetic differences in lung squamous cell carcinoma patients will be helpful in understanding the underlying causes of this immunotherapy effect and provide a reference for clinical practice.

METHODS

In this study, we used RNA-seq, miRNA-seq, methylation array, mutation profiles, and copy number variation data from the TCGA database and RNA-seq data from the GEO database to analyze the distinctive genomic patterns associated with PD1 and PDL1 expression. RNA-seq data from 44 LUSC patients who underwent surgery at Zhongshan Hospital were also included in the study.

RESULTS

After grouping LUSC patients according to the expression levels of PD1 and PDL1, we found no significant difference in survival between the two groups. However, 178 genes, including IL-21, KLRC3, and KLRC4, were significantly upregulated in both the TCGA and GEO databases in the high expression group, and there was a precise correlation between gene expression and epigenetic changes in the two groups. At the same time, the overall level of somatic mutations was not significantly different between the two groups. It is worth noting that the gene enrichment results showed that the differential pathways were mainly enriched in immune regulation, especially T cell-related immune activities.

CONCLUSION

We found that the differences in gene expression between the two groups were related to immunity, which may affect the effectiveness of immunotherapy. We hope our results can provide a reference for further research and help in finding other targets to improve the efficacy of immunotherapy.

Analysis of Copy-Number Variations and Feline Mammary Carcinoma Survival

Feline mammary carcinomas (FMCs) are highly malignant. As the disease-free survival (DFS) and overall survival (OS) are short, prognostication is crucial. Copy-number variations (CNVs) analysis by next-generation sequencing serves to identify critical cancer-related genomic regions. Thirty-three female cats with FMCs were followed during two years after surgery. Tumours represented tubulopapillary and solid carcinomas encompassing six molecular subtypes. Regardless of the histopathological diagnosis, molecular subtypes showed important differences in survival. Luminal A tumours exhibited the highest DFS (p = 0.002) and cancer-specific OS (p = 0.001), and the lowest amount of CNVs (p = 0.0001). In contrast, basal-like triple-negative FMCs had the worst outcome (DFS, p < 0.0001; and OS, p < 0.00001) and were the most aberrant (p = 0.05). In the multivariate analysis, copy-number losses (CNLs) in chromosome B1 (1-23 Mb) harbouring several tumour-repressors (e.g. CSMD1, MTUS1, MSR1, DBC2, and TUSC3) negatively influenced DFS. Whereas, copy-number gains (CNGs) in B4 (1-29 Mb) and F2 (64-82.3 Mb) comprising epithelial to mesenchymal transition genes and metastasis-promoting transcription factors (e.g. GATA3, VIM, ZEB1, and MYC) negatively influenced DFS and cancer-specific OS. These data evidence an association between specific CNVs in chromosomes B1, B4 and F2, and poor prognosis in FMCs.

Clinical Significance of CUB and Sushi Multiple Domains 1 Inactivation in Head and Neck Squamous Cell Carcinoma

Although the genetic alteration of CUB and Sushi multiple domains 1 (CSMD1) is known to be associated with poor prognosis in several cancers, there is a lack of clinical relevance in head and neck cancer. The aim of this study was to offer insight into the clinical significance of CSMD1, utilizing a multimodal approach that leverages publicly available independent genome-wide expression datasets. CSMD1-related genes were found and analyzed to examine the clinical significance of CSMD1 inactivation in the HNSCC cohort of publicly available databases. We analyzed the frequency of somatic mutations, clinicopathologic characteristics, association with immunotherapy-related gene signatures, and the pathways of gene signatures. We found 363 CSMD1-related genes. The prognosis of the CSMD1-inactivated subgroup was poor. FBXW7, HLA-A, MED1, NOTCH2, NOTCH3, and TP53 had higher mutation rates in the CSMD1-inactivated subgroups. The Interferon-gamma score and immune signature score were elevated in CSMD1-inactivated subgroups. We identified several CSMD1-related pathways, such as the phosphatidylinositol signaling system and inositol phosphate metabolism. Our study using three large and independent datasets suggests that CSMD1-related gene signatures are associated with the prognosis of HNSCC patients.

Cloud-based interactive analytics for terabytes of genomic variants data

Motivation

Large scale genomic sequencing is now widely used to decipher questions in diverse realms such as biological function, human diseases, evolution, ecosystems, and agriculture. With the quantity and diversity these data harbor, a robust and scalable data handling and analysis solution is desired.

Results

We present interactive analytics using a cloud-based columnar database built on Dremel to perform information compression, comprehensive quality controls, and biological information retrieval in large volumes of genomic data. We demonstrate such Big Data computing paradigms can provide orders of magnitude faster turnaround for common genomic analyses, transforming long-running batch jobs submitted via a Linux shell into questions that can be asked from a web browser in seconds. Using this method, we assessed a study population of 475 deeply sequenced human genomes for genomic call rate, genotype and allele frequency distribution, variant density across the genome, and pharmacogenomic information.

Availability and implementation

Our analysis framework is implemented in Google Cloud Platform and BigQuery. Codes are available at https://github.com/StanfordBioinformatics/mvp_aaa_codelabs.

Contact

cuiping@stanford.edu or ptsao@stanford.edu.

Supplementary information

Supplementary data are available at Bioinformatics online.

Differences in cancer gene copy number alterations between Epstein-Barr virus-positive and Epstein-Barr virus-negative nasopharyngeal carcinoma

BACKGROUND

Nasopharyngeal carcinoma (NPC) treatment is mainly based on clinical staging. We hypothesize that better understanding of the molecular heterogeneity of NPC can aid in better treatment decisions. Therefore, the purpose of this study was to present our exploration of cancer gene copy-number alterations (CNAs) of Epstein-Barr virus (EBV)-positive and EBV-negative NPC.

METHODS

Multiplex ligation-dependent probe amplification was applied to detect CNAs of 36 cancer genes (n = 103). Correlation between CNAs, clinicopathological features, and survival were examined.

RESULTS

The CNAs occurred significantly more in EBV-negative NPC, with PIK3CA and MCCC1 (P < .001) gain/amplification occurring more frequently. Gain/amplification of cyclin-L1 (CCNL1) and PTK2 (P < .001) predict worse disease-free survival (DFS) in EBV-positive NPC.

CONCLUSION

The EBV-positive and EBV-negative NPC show some similarities in cancer gene CNAs suggesting a common pathogenic route but also important differences possibly indicating divergence in oncogenesis. Copy number gain/amplification of CCNL1 and PTK2 are possibly good predictors of survival in EBV-positive NPC.

Novel potential inhibitors of complement system and their roles in complement regulation and beyond

The complement system resembles a double-edged sword since its activation can either benefit or harm the host. Thus, regulation of this system is of utmost importance and performed by several circulating and membrane-bound complement inhibitors. The pool of well-established regulators has recently been enriched with proteins that either share structural homology to known complement inhibitors such as Sushi domain-containing (SUSD) protein family and Human CUB and Sushi multiple domains (CSMD) families or extracellular matrix (ECM) macromolecules that interact with and modulate complement activity. In this review, we summarize the current knowledge about newly discovered complement inhibitors and discuss their implications in complement regulation, as well as in processes beyond complement regulation such cancer development. Understanding the behavior of these proteins will introduce new mechanisms of complement regulation and may provide new avenues in the development of novel therapies.

Complement inhibitor CSMD1 acts as tumor suppressor in human breast cancer

Human CUB and Sushi multiple domains 1 (CSMD1) is a membrane-bound complement inhibitor suggested to act as a putative tumor suppressor gene, since allelic loss of this region encompassing 8p23 including CSMD1 characterizes various malignancies. Here, we assessed the role of CSMD1 as a tumor suppressor gene in the development of breast cancer in vitro and in vivo. We found that human breast tumor tissues expressed CSMD1 at lower levels compared to that in normal mammary tissues. The decreased expression of CSMD1 was linked to a shorter overall survival of breast cancer patients. We also revealed that expression of CSMD1 in human breast cancer cells BT-20 and MDA-MB-231 significantly inhibited their malignant phenotypes, including migration, adhesion and invasion. Conversely, stable silencing of CSMD1 expression in T47D cells enhanced cancer cell migratory, adherent and clonogenic abilities. Moreover, expression of CSMD1 in the highly invasive MDA-MB-231 cells diminished their signaling potential as well as their stem cell-like properties as assessed by measurement of aldehyde dehydrogenase activity. In a xenograft model, expression of CSMD1 blocked the ability of cancer cells to metastasize to secondary sites in vivo, likely via inhibiting local invasion but not the extravasation into distant tissues. Taken together, these findings demonstrate the role of CSMD1 as a tumor suppressor gene in breast cancer.

The role of complement inhibitors beyond controlling inflammation

The complement system is an arm of innate immunity that aids in the removal of pathogens and dying cells. Due to its harmful, pro-inflammatory potential, complement is controlled by several soluble and membrane-bound inhibitors. This family of complement regulators has been recently extended by the discovery of several new members, and it is becoming apparent that these proteins harbour additional functions. In this review, the current state of knowledge of the physiological functions of four complement regulators will be described: cartilage oligomeric matrix protein (COMP), CUB and sushi multiple domains 1 (CSMD1), sushi domain-containing protein 4 (SUSD4) and CD59. Complement activation is involved in both the development of and defence against cancer. COMP expression is pro-oncogenic, whereas CSMD1 and SUSD4 act as tumour suppressors. These effects may be related in part to the complex influence of complement on cancer but also depend on unrelated functions such as the protection of cells from endoplasmic reticulum stress conveyed by intracellular COMP. CD59 is the main inhibitor of the membrane attack complex, and its deficiency leads to complement attack on erythrocytes and severe haemolytic anaemia, which is now amenable to treatment with an inhibitor of C5 cleavage. Unexpectedly, the intracellular pool of CD59 is crucial for insulin secretion from pancreatic β-cells. This finding is one of several relating to the intracellular functions of complement proteins, which until recently were only considered to be present in the extracellular space. Understanding the alternative functions of complement inhibitors may unravel unexpected links between complement and other physiological systems, but is also important for better design of therapeutic complement inhibition.

Loss of CSMD1 expression disrupts mammary duct formation while enhancing proliferation, migration and invasion

The CUB and sushi multiple domains 1 (CSMD1) gene maps to chromosome 8p23, a region deleted in many cancers. Loss of CSMD1 expression is associated with poor prognosis in breast cancer suggesting that it acts as a tumour suppressor in this cancer. However, the function of CSMD1 is largely unknown. Herein, we investigated CSMD1 functions in cell line models. CSMD1 expression was suppressed in MCF10A and LNCaP cells using short hairpin RNA. Functional assays were performed focusing on the 'normal' MCF10A cell line. Suppression of CSMD1 significantly increased the proliferation, cell migration and invasiveness of MCF10A cells compared to shcontrols. shCSMD1 cells also showed significantly reduced adhesion to Matrigel and fibronectin. In a three-dimensional Matrigel model of MCF10A cells, reduced CSMD1 expression resulted in the development of larger and more poorly differentiated breast acini-like structures that displayed impaired lumen formation. Loss of CSMD1 expression disrupts a model of mammary duct formation while enhancing proliferation, migration and invasion. Our data suggest that CSMD1 is involved in the suppression of a transformed phenotype.

Copy number variation in archival melanoma biopsies versus benign melanocytic lesions

BACKGROUND

Skin melanocytes can give rise to benign and malignant neoplasms. Discrimination of an early melanoma from an unusual/atypical benign nevus can represent a significant challenge. However, previous studies have shown that in contrast to benign nevi, melanoma demonstrates pervasive chromosomal aberrations.

OBJECTIVE

This substantial difference between melanoma and benign nevi can be exploited to discriminate between melanoma and benign nevi.

METHODS

Array-comparative genomic hybridization (aCGH) is an approach that can be used on DNA extracted from formalin-fixed paraffin-embedded (FFPE) tissues to assess the entire genome for the presence of changes in DNA copy number. In this study, high resolution, genome-wide single-nucleotide polymorphism (SNP) arrays were utilized to perform comprehensive and detailed analyses of recurrent copy number aberrations in 41 melanoma samples in comparison with 21 benign nevi.

RESULTS

We found statistically significant copy number gains and losses within melanoma samples. Some of the identified aberrations are previously implicated in melanoma. Moreover, novel regions of copy number alterations were identified, revealing new candidate genes potentially involved in melanoma pathogenesis.

CONCLUSIONS

Taken together, these findings can help improve melanoma diagnosis and introduce novel melanoma therapeutic targets.

Genome-Wide Association Studies of Multiple Keratinocyte Cancers

There is strong evidence for a role of environmental risk factors involved in susceptibility to develop multiple keratinocyte cancers (mKCs), but whether genes are also involved in mKCs susceptibility has not been thoroughly investigated. We investigated whether single nucleotide polymorphisms (SNPs) are associated with susceptibility for mKCs. A genome-wide association study (GWAS) of 1,666 cases with mKCs and 1,950 cases with single KC (sKCs; controls) from Harvard cohorts (the Nurses' Health Study [NHS], NHS II, and the Health Professionals Follow-Up Study) and the Framingham Heart Study was carried-out using over 8 million SNPs (stage-1). We sought to replicate the most significant statistical associations (p-value≤ 5.5x10-6) in an independent cohort of 574 mKCs and 872 sKCs from the Rotterdam Study. In the discovery stage, 40 SNPs with suggestive associations (p-value ≤5.5x10-6) were identified, with eight independent SNPs tagging all 40 SNPs. The most significant SNP was located at chromosome 9 (rs7468390; p-value = 3.92x10-7). In stage-2, none of these SNPs replicated and only two of them were associated with mKCs in the same direction in the combined meta-analysis. We tested the associations for 19 previously reported basal cell carcinoma-related SNPs (candidate gene association analysis), and found that rs1805007 (MC1R locus) was significantly associated with risk of mKCs (p-value = 2.80x10-4). Although the suggestive SNPs with susceptibility for mKCs were not replicated, we found that previously identified BCC variants may also be associated with mKC, which the most significant association (rs1805007) located at the MC1R gene.

miR-10b exerts oncogenic activity in human hepatocellular carcinoma cells by targeting expression of CUB and sushi multiple domains 1 (CSMD1)

BACKGROUND

Hepatocellular carcinoma (HCC) is a lethal disease, while the precise underlying molecular mechanisms of HCC pathogenesis remain to be defined. MicroRNA (miRNA), a class of non-coding small RNAs, can post-transcriptionally regulate gene expression. Altered miRNA expression has been reported in HCCs. This study assessed expression and the oncogenic activity of miRNA-10b (miR-10b) in HCC.

METHODS

Forty-five paired human HCC and adjacent non-tumor tissues were collected for qRT-PCR and immunohistochemistry analysis of miR-10b and CUB and Sushi multiple domains 1 (CSMD1), respectively. We analyzed the clinicopathological data from these patients to further determine if there was an association between miR-10b and CSMD1. HCC cell lines were used to assess the effects of miR-10b mimics or inhibitors on cell viability, migration, invasion, cell cycle distribution, and colony formation. Luciferase assay was used to assess miR-10b binding to the 3'-untranslated region (3'-UTR) of CSMD1.

RESULTS

miR-10b was highly expressed in HCC tissues compared to normal tissues. In vitro, overexpression of miR-10b enhanced HCC cell viability, migration, and invasion; whereas, downregulation of miR-10b expression suppressed these properties in HCC cells. Injection of miR-10b mimics into tumor cell xenografts also promoted xenograft growth in nude mice. Bioinformatics and luciferase reporter assay demonstrated that CSMD1 was the target gene of miR-10b. Immunocytochemical, immunohistochemical, and qRT-PCR data indicated that miR-10b decreased CSMD1 expression in HCC cells.

CONCLUSIONS

We showed that miR-10b is overexpressed in HCC tissues and miR-10b mimics promoted HCC cell viability and invasion via targeting CSMD1 expression. Our findings suggest that miR-10b acts as an oncogene by targeting the tumor suppressor gene, CSMD1, in HCC.

Genome-wide Association Study of Postburn Scarring Identifies a Novel Protective Variant

OBJECTIVE

To identify genetic variants associated with the severity of postburn hypertrophic scarring (HTS) using a genome-wide approach.

BACKGROUND

Risk of severe postburn HTS is known to depend on race, but the genetic determinants of HTS are unknown.

METHODS

We conducted a genome-wide association study (GWAS) in a prospective cohort of adults admitted with deep-partial-thickness burns from 2007 through 2014. Scar severity was assessed over time using the Vancouver Scar Scale (VSS), and DNA was genotyped with a >500,000-marker array. We performed association testing of single-nucleotide polymorphisms (SNPs) with minor allele frequency (MAF) >0.01 using linear regression of VSS height score on genotype adjusted for patient and injury characteristics as well as population genetic structure. Array-wide significance was based on Bonferroni correction for multiple testing.

RESULTS

Of 538 patients (median age 40 years, median burn size 6.0% of body surface area), 71% were men and 76% were White. The mean VSS height score was 1.2 (range: 0-3). Of 289,639 SNPs tested, a variant in the CUB and Sushi multiple domains 1 (CSMD1) gene (rs11136645; MAF = 0.49), was significantly associated with decreased scar height (regression coefficient = -0.23, P = 7.9 × 10).

CONCLUSIONS

In the first published GWAS of HTS, we report that a common intronic variant in the CSMD1 gene is associated with reduced severity of postburn HTS. If this association is confirmed in an independent cohort, investigating the potential role of CSMD1 in wound healing may elucidate HTS pathophysiology.

Genetic imbalances detected by multiplex ligation-dependent probe amplification in a cohort of patients with oral squamous cell carcinoma-the first step towards clinical personalized medicine

Oral tumors are a growing health problem worldwide; thus, it is mandatory to establish genetic markers in order to improve diagnosis and early detection of tumors, control relapses and, ultimately, delineate individualized therapies. This study was the first to evaluate and discuss the clinical applicability of a multiplex ligation-dependent probe amplification (MLPA) probe panel directed to head and neck cancer. Thirty primary oral squamous cell tumors were analyzed using the P428 MLPA probe panel. We detected genetic imbalances in 26 patients and observed a consistent pattern of distribution of genetic alterations in terms of losses and gains for some chromosomes, particularly for chromosomes 3, 8, and 11. Regarding the latter, some specific genes were highlighted due to frequent losses of genetic material--RARB, FHIT, CSMD1, GATA4, and MTUS1--and others due to gains--MCCC1, MYC, WISP1, PTK2, CCND1, FGF4, FADD, and CTTN. We also verified that the gains of MYC and WISP1 genes seem to suggest higher propensity of tumors localized in the floor of the mouth. This study proved the value of this MLPA probe panel for a first-tier analysis of oral tumors. The probemix was developed to include target regions that have been already shown to be of diagnostic/prognostic relevance for oral tumors. Furthermore, this study emphasized several of those specific genetic targets, suggesting its importance to oral tumor development, to predict patients' outcomes, and also to guide the development of novel molecular therapies.

CSMD1 exhibits antitumor activity in A375 melanoma cells through activation of the Smad pathway

In this work, we studied the effects of CUB and Sushi multiple domains 1 gene (CSMD1) expression in A375 melanoma cells in vivo and in vitro. CSDM1 expression decreased proliferation and migration, and increased apoptosis and G(1) arrest in A375 cells in vitro. Expression of CSDM1 in established xenografted tumors decreased tumor size and weight, and decreased the density of intratumor microvessels. The survival rate of mice with tumors expressing CSMD1 was significantly higher than mice with tumors that did not express CSDM1. These results confirm the role of CSDM1 as a tumor suppressor gene in melanoma cells. Furthermore, we found that CSMD1 can interact with Smad3, activate Smad1, Smad2, and Smad3, and increase the expression of Smad4. These results might prove helpful for the development of novel therapies for melanoma treatment.

Neuropsychological deficits in mice depleted of the schizophrenia susceptibility gene CSMD1

Recent meta-analyses of schizophrenia genome-wide association studies (GWASs) have identified the CUB and SUSHI multiple domains 1 (CSMD1) gene as a statistically strong risk factor. CSMD1 is a complement control-related protein suggested to inhibit the classical complement pathway, being expressed in developing neurons. However, expression of CSMD1 is largely uncharacterized and relevance for behavioral phenotypes is not previously demonstrated. Here, we assess neuropsychological behaviors of a Csmd1 knockout (KO) mouse in a selection of standard behavioral tests. Deregulation of neuropsychological responses were observed in both the open field and the elevated plus maze tests, in which KO mice spent 55% and 33% less time than WT littermate mice in open areas, respectively. Altered behaviors were also observed in tail suspension and to higher acoustic stimuli, for which Csmd1 KO mice showed helplessness and moderate increase in startle amplitude, respectively. Furthermore, Csmd1 KO mice also displayed increased weight-gain and glucose tolerance, similar to a major phenotype of the metabolic syndrome that also has been associated to the human CSMD1 locus. Consistent with a role in the control of behaviors, Csmd1 was found highly expressed in the central nervous system (CNS), and with some expression in visceral fat and ovary, under tissue-specific control by a novel promoter-associated lncRNA. In summary, disruption of Csmd1 induces behaviors reminiscent of blunted emotional responses, anxiety and depression. These observations suggest an influence of the CSMD1 schizophrenia susceptibility gene on psychopathology and endophenotypes of the negative symptom spectra.

The novel complement inhibitor human CUB and Sushi multiple domains 1 (CSMD1) protein promotes factor I-mediated degradation of C4b and C3b and inhibits the membrane attack complex assembly

CUB and Sushi multiple domains 1 (CSMD1) is a transmembrane protein containing 15 consecutive complement control protein (CCP) domains, which are characteristic for complement inhibitors. We expressed a membrane-bound fragment of human CSMD1 composed of the 15 C-terminal CCP domains and demonstrated that it inhibits deposition of C3b by the classical pathway on the surface of Chinese hamster ovary cells by 70% at 6% serum and of C9 (component of membrane attack complex) by 90% at 1.25% serum. Furthermore, this fragment of CSMD1 served as a cofactor to factor I-mediated degradation of C3b. In all functional assays performed, well-characterized complement inhibitors were used as positive controls, whereas Coxsackie adenovirus receptor, a protein with no effect on complement, was a negative control. Moreover, attenuation of expression in human T47 breast cancer cells that express endogenous CSMD1 significantly increased C3b deposition on these cells by 45% at 8% serum compared with that for the controls. Furthermore, by expressing a soluble 17-21 CCP fragment of CSMD1, we found that CSMD1 inhibits complement by promoting factor I-mediated C4b/C3b degradation and inhibition of MAC assembly at the level of C7. Our results revealed a novel complement inhibitor for the classical and lectin pathways.

Somatic mutations, allele loss, and DNA methylation of the Cub and Sushi Multiple Domains 1 (CSMD1) gene reveals association with early age of diagnosis in colorectal cancer patients

Background

The Cub and Sushi Multiple Domains 1 (CSMD1) gene, located on the short arm of chromosome 8, codes for a type I transmembrane protein whose function is currently unknown. CSMD1 expression is frequently lost in many epithelial cancers. Our goal was to characterize the relationships between CSMD1 somatic mutations, allele imbalance, DNA methylation, and the clinical characteristics in colorectal cancer patients.

Methods

We sequenced the CSMD1 coding regions in 54 colorectal tumors using the 454FLX pyrosequencing platform to interrogate 72 amplicons covering the entire coding sequence. We used heterozygous SNP allele ratios at multiple CSMD1 loci to determine allelic balance and infer loss of heterozygosity. Finally, we performed methylation-specific PCR on 76 colorectal tumors to determine DNA methylation status for CSMD1 and known methylation targets ALX4, RUNX3, NEUROG1, and CDKN2A.

Results

Using 454FLX sequencing and confirming with Sanger sequencing, 16 CSMD1 somatic mutations were identified in 6 of the 54 colorectal tumors (11%). The nonsynonymous to synonymous mutation ratio of the 16 somatic mutations was 15∶1, a ratio significantly higher than the expected 2∶1 ratio (p = 0.014). This ratio indicates a presence of positive selection for mutations in the CSMD1 protein sequence. CSMD1 allelic imbalance was present in 19 of 37 informative cases (56%). Patients with allelic imbalance and CSMD1 mutations were significantly younger (average age, 41 years) than those without somatic mutations (average age, 68 years). The majority of tumors were methylated at one or more CpG loci within the CSMD1 coding sequence, and CSMD1 methylation significantly correlated with two known methylation targets ALX4 and RUNX3. C:G>T:A substitutions were significantly overrepresented (47%), suggesting extensive cytosine methylation predisposing to somatic mutations.

Conclusions

Deep amplicon sequencing and methylation-specific PCR reveal that CSMD1 alterations can correlate with earlier clinical presentation in colorectal tumors, thus further implicating CSMD1 as a tumor suppressor gene.

Assessment of behaviors modeling aspects of schizophrenia in Csmd1 mutant mice

Schizophrenia is a debilitating psychotic disorder that affects up to 1.5% of the population worldwide. Two recent studies in humans identified genome-wide significant associations between schizophrenia and single-nucleotide polymorphisms (SNPs) in an intron of CSMD1. The effect of deleting CSMD1 on mouse behavior is unknown. The present study utilized mice with a mutant Csmd1 allele in which the first exon had been ablated (KO mice). All Csmd1 transcripts that included the first exon were absent in the brains of KO mice, but there was persistent expression of at least one other transcript that does not include the first exon. Wild type (WT), heterozygous (HET), and KO mice were assessed using several well-established behavioral paradigms that model aspects of schizophrenia. Csmd1 KO mice did not differ from wild-type littermates for sensorimotor gating (measured as prepulse inhibition), social interaction, anhedonia (measured by sucrose preference), or sensitivity to the locomotor stimulant effects of the dopaminergic agent d-amphetamine. These data demonstrate that loss of Csmd1 transcripts that include the first exon does not alter multiple well-established behaviors that model aspects of schizophrenia. The SNP most strongly associated with schizophrenia in humans is between exons 3 and 4; therefore, ablation of exon 1 appeared to be a logical animal model. Nevertheless, future studies should consider alternative mouse models including gain-of-function mutations, and loss-of-function mutations that target alternative transcripts of Csmd1.

Molecular analysis of desmoid tumors with a high-density single-nucleotide polymorphism array identifies new molecular candidate lesions

BACKGROUND

Desmoid tumors are neoplastic proliferations of connective tissues. The mutation status of the gene coding for catenin (cadherin-associated protein) beta 1 (CTNNB1) and trisomy 8 on the chromosomal level have been described to have prognostic relevance.

PATIENTS AND METHODS

In order to elucidate new molecular mechanisms underlying these tumors, we carried out a molecular analysis with a genome-wide human high-density single-nucleotide polymorphism (SNP) array, in 9 patients.

RESULTS

Single samples showed numerical aberrations on chromosomes (Chrs) 20 and 6 with either trisomy 20 or monosomy 6. No trisomy 8 could be detected. Recurrent heterozygous deletions were found in Chr 5q (including the APC gene locus, n = 3) and Chr 8p23 (n = 4, containing coding regions for the potential tumor suppressor gene CSMD1). This novel deletion in 8p23 showed an association with local recurrence. In addition, structural chromosomal changes (gain of Chrs 8 and 20) were found in a minority of cases.

CONCLUSION

The genomic alteration affecting the candidate gene CSMD1 could be important in the development of desmoid tumors.

Comparative analysis of algorithms for integration of copy number and expression data

Chromosomal instability is a hallmark of cancer, and genes that display abnormal expression in aberrant chromosomal regions are likely to be key players in tumor progression. Identifying such driver genes reliably requires computational methods that can integrate genome-scale data from several sources. We compared the performance of ten algorithms that integrate copy-number and transcriptomics data from 15 head and neck squamous cell carcinoma cell lines, 129 lung squamous cell carcinoma primary tumors and simulated data. Our results revealed clear differences between the methods in terms of sensitivity and specificity as well as in their performance in small and large sample sizes. Results of the comparison are available at http://csbi.ltdk.helsinki.fi/cn2gealgo/.

Integrated genomic analysis of relapsed childhood acute lymphoblastic leukemia reveals therapeutic strategies

Despite an increase in survival for children with acute lymphoblastic leukemia (ALL), the outcome after relapse is poor. To understand the genetic events that contribute to relapse and chemoresistance and identify novel targets of therapy, 3 high-throughput assays were used to identify genetic and epigenetic changes at relapse. Using matched diagnosis/relapse bone marrow samples from children with relapsed B-precursor ALL, we evaluated gene expression, copy number abnormalities (CNAs), and DNA methylation. Gene expression analysis revealed a signature of differentially expressed genes from diagnosis to relapse that is different for early (< 36 months) and late (≥ 36 months) relapse. CNA analysis discovered CNAs that were shared at diagnosis and relapse and others that were new lesions acquired at relapse. DNA methylation analysis found increased promoter methylation at relapse. There were many genetic alterations that evolved from diagnosis to relapse, and in some cases these genes had previously been associated with chemoresistance. Integration of the results from all 3 platforms identified genes of potential interest, including CDKN2A, COL6A2, PTPRO, and CSMD1. Although our results indicate that a diversity of genetic changes are seen at relapse, integration of gene expression, CNA, and methylation data suggest a possible convergence on the WNT and mitogen-activated protein kinase pathways.

The molecular biology of head and neck cancer

Head and neck squamous cell carcinomas (HNSCCs) are caused by tobacco and alcohol consumption and by infection with high-risk types of human papillomavirus (HPV). Tumours often develop within preneoplastic fields of genetically altered cells. The persistence of these fields after treatment presents a major challenge, because it might lead to local recurrences and second primary tumours that are responsible for a large proportion of deaths. Aberrant signalling pathways have been identified in HNSCCs and inhibition of epidermal growth factor receptor (EGFR) has proved a successful therapeutic strategy. In this Review, we discuss the recent literature on tumour heterogeneity, field cancerization, molecular pathogenesis and the underlying causative cancer genes that can be exploited for novel and personalized treatments of patients with HNSCC.

The molecular biology of head and neck cancer

Head and neck squamous cell carcinomas (HNSCCs) are caused by tobacco and alcohol consumption and by infection with high-risk types of human papillomavirus (HPV). Tumours often develop within preneoplastic fields of genetically altered cells. The persistence of these fields after treatment presents a major challenge, because it might lead to local recurrences and second primary tumours that are responsible for a large proportion of deaths. Aberrant signalling pathways have been identified in HNSCCs and inhibition of epidermal growth factor receptor (EGFR) has proved a successful therapeutic strategy. In this Review, we discuss the recent literature on tumour heterogeneity, field cancerization, molecular pathogenesis and the underlying causative cancer genes that can be exploited for novel and personalized treatments of patients with HNSCC.

Genetic variations in ATP2B1, CSK, ARSG and CSMD1 loci are related to blood pressure and/or hypertension in two Korean cohorts

Blood pressure, one of the important vital signs, is affected by multiple genetic and environmental factors. Recently, several genome-wide association (GWA) studies have successfully identified genetic factors that influence blood pressure and hypertension risk. In this study, we report results of the Korean Association REsource (KARE, 8842 subjects) GWA study on blood pressure and hypertension risk. In all, 10 single-nucleotide polymorphisms (SNPs) that showed significant association with hypertension were further analysed for replication associations in the Health2 project (7861 subjects). Among these 10 SNPs, 3 were replicated in the Health2 cohort for an association with systolic or diastolic blood pressure. The most significant SNP (rs17249754 located in ATPase, Ca(++) transporting, plasma membrane 1 (ATP2B1)) has been previously reported, and the other two SNPs are rs1378942 in the c-src tyrosine kinase (CSK) gene and rs12945290 in the arylsulphatase G (ARSG) gene. An additional hypertension case-control study confirmed that rs17249754 (in ATP2B1) increases hypertension risk in both the KARE and Health2 (meta-analysis, P-value=4.25 x 10(-9)) cohorts. One more SNP, rs995322, located in the CUB and Sushi multiple domains 1 (CSMD1), is also associated with increased risk of hypertension (meta-analysis, P-value=1.00 x 10(-4)). Despite the difficulty of obtaining replication results for a complex trait genetic association between blood pressure and hypertension, we were able to identify consistent genetic factors in both the Korean cohorts in ATP2B1, CSK, ARSG and CSMD1 genes.

Comparative analyses of chromosome alterations in soft-tissue metastases within and across patients with castration-resistant prostate cancer

Androgen deprivation is the mainstay of therapy for progressive prostate cancer. Despite initial and dramatic tumor inhibition, most men eventually fail therapy and die of metastatic castration-resistant (CR) disease. Here, we characterize the profound degree of genomic alteration found in CR tumors using array comparative genomic hybridization (array CGH), gene expression arrays, and fluorescence in situ hybridization (FISH). Bycluster analysis, we show that the similarity of the genomic profiles from primary and metastatic tumors is driven by the patient. Using data adjusted for this similarity, we identify numerous high-frequency alterations in the CR tumors, such as 8p loss and chromosome 7 and 8q gain. By integrating array CGH and expression array data, we reveal genes whose correlated values suggest they are relevant to prostate cancer biology. We find alterations that are significantly associated with the metastases of specific organ sites, and others with CR tumors versus the tumors of patients with localized prostate cancer not treated with androgen deprivation. Within the high-frequency sites of loss in CR metastases, we find an overrepresentation of genes involved in cellular lipid metabolism, including PTEN. Finally, using FISH, we verify the presence of a gene fusion between TMPRSS2 and ERG suggested by chromosome 21 deletions detected by array CGH. We find the fusion in 54% of our CR tumors, and 81% of the fusion-positive tumors contain cells with multiple copies of the fusion. Our investigation lays the foundation for a better understanding of and possible therapeutic targets for CR disease, the poorly responsive and final stage of prostate cancer.

Loss of CSMD1 expression is associated with high tumour grade and poor survival in invasive ductal breast carcinoma

CUB and SUSHI multiple domain protein 1 (CSMD1) is a candidate tumour suppressor gene that maps to chromosome 8p23, a region deleted in many tumour types including 50% of breast cancers. CSMD1 has homologies to proteins implicated in carcinogenesis. We aimed to study the expression pattern of the CSMD1 protein and evaluate its prognostic importance in invasive ductal carcinoma (IDC). An anti-CSMD1 antibody was developed and validated. The expression pattern of CSMD1 in normal breast and IDC samples was investigated by immunohistochemistry in 275 patients. Univariate and multivariate Cox regression analyses were performed. In normal breast duct epithelial cells, luminal, membranous and cytoplasmic CSMD1 staining was identified. Reduced expression of CSMD1 was detected in 79/275 (28.7%) of IDC cases. Low CSMD1 expression was significantly associated with high tumour grade (P = 0.003). CSMD1 expression was associated with overall survival (OS; HR = 0.607, 95%CI: 0.4-0.91, P = 0.018) but not with disease-free survival (DFS; HR = 0.81, 95%CI: 0.46-1.43, P = 0.48). Multivariate analysis showed that CSMD1, together with Nottingham Prognostic Index, was considered an independent predictor of OS (HR = 0.607, 95%CI: 0.4-0.91, P = 0.018) but not DFS (HR = 0.84, 95%CI: 0.46-1.5, P = 0.573). Reduction of CSMD1 expression was significantly associated with high tumour grade and decreased OS. Therefore, our results support the idea that CSMD1 is a tumour suppressor gene and suggest its possible use as a new prognostic biomarker. The membrane expression pattern of CSMD1 suggests that it may be a receptor or co-receptor involved in the process of signal transduction.

Patterns of gene expression and copy-number alterations in von-hippel lindau disease-associated and sporadic clear cell carcinoma of the kidney

Recent insights into the role of the von-Hippel Lindau (VHL) tumor suppressor gene in hereditary and sporadic clear-cell renal cell carcinoma (ccRCC) have led to new treatments for patients with metastatic ccRCC, although virtually all patients eventually succumb to the disease. We performed an integrated, genome-wide analysis of copy-number changes and gene expression profiles in 90 tumors, including both sporadic and VHL disease-associated tumors, in hopes of identifying new therapeutic targets in ccRCC. We identified 14 regions of nonrandom copy-number change, including 7 regions of amplification (1q, 2q, 5q, 7q, 8q, 12p, and 20q) and 7 regions of deletion (1p, 3p, 4q, 6q, 8p, 9p, and 14q). An analysis aimed at identifying the relevant genes revealed VHL as one of three genes in the 3p deletion peak, CDKN2A and CDKN2B as the only genes in the 9p deletion peak, and MYC as the only gene in the 8q amplification peak. An integrated analysis to identify genes in amplification peaks that are consistently overexpressed among amplified samples confirmed MYC as a potential target of 8q amplification and identified candidate oncogenes in the other regions. A comparison of genomic profiles revealed that VHL disease-associated tumors are similar to a subgroup of sporadic tumors and thus more homogeneous overall. Sporadic tumors without evidence of biallelic VHL inactivation fell into two groups: one group with genomic profiles highly dissimilar to the majority of ccRCC and a second group with genomic profiles that are much more similar to tumors with biallelic inactivation of VHL.

Association of chromosome 8q variants with prostate cancer risk in Caucasian and Hispanic men

Genotyping of a 615 kb region within 8q24 with 49 haplotype-tagged single-nucleotide polymorphisms (SNPs) in 2109 samples (797 cases and 1312 controls) of two ethnic/racial groups found SNPs that are significantly associated with the risk for prostate cancer (PCa). The highest significance in Caucasian men was found for rs6983267; the AA genotype reduced the risk for PCa [odds ratio (OR) = 0.48, 95% confidence interval (CI) = 0.35-0.65, P = 2.74 x 10(-6)]. This SNP also had a significant independent effect from other SNPs in the region in this group. In Hispanic men, rs7837328 and rs921146 showed independent effects (OR = 2.55, 95% CI = 1.51-4.31, P = 4.33 x 10(-4), OR = 2.09, 95% CI = 1.40-3.12, P = 3.13 x 10(-4), respectively). Significant synergist effects for increasing numbers of high-risk alleles were found in both ethnicities. Haplotype analysis revealed major haplotypes, containing the non-risk alleles, conferred protection against PCa. We found high linkage disequilibrium between significant SNPs within the region and SNPs within the CUB and Sushi Multiple Domains 1 gene (CSMD1), on the short arm of chromosome 8 in both ethnicities. These data suggest that multiple interacting SNPs within 8q24, as well as different regions on chromosome 8 far beyond this 8q24 candidate region, may confer increased risk of PCa. This is the first report to investigate the involvement of 8q24 variants in the susceptibility for PCa in Hispanic men.

Analysis of DNA copy number alterations in ovarian serous tumors identifies new molecular genetic changes in low-grade and high-grade carcinomas

Ovarian serous carcinoma, the most common and lethal type of ovarian cancer, is thought to develop from two distinct molecular pathways. High-grade (HG) serous carcinomas contain frequent TP53 mutations, whereas low-grade (LG) carcinomas arise from serous borderline tumors (SBT) and harbor mutations in KRAS/BRAF/ERBB2 pathway. However, the molecular alterations involved in the progression from SBT to LG carcinoma remain unknown. In addition, the extent of deletion of tumor suppressors in ovarian serous carcinomas has not been well studied. To further address these two issues, we assessed DNA copy number changes among affinity-purified tumor cells from 37 ovarian serous neoplasms including SBT, LG, and HG tumors using high-density 250K single nucleotide polymorphism arrays. Chromosomal instability index as measured by changes in DNA copy number was significantly higher in HG than in LG serous carcinomas. Hemizygous ch1p36 deletion was common in LG serous carcinomas but was rarely seen in SBT. This region contains several candidate tumor suppressors including miR-34a. In contrast, in HG serous carcinomas, significant numbers of amplifications and deletions, including homozygous deletions, were identified. Among homozygous deletions, loci containing Rb1, CDKN2A/B, CSMD1, and DOCK4 were most common, being present in 10.6%, 6.4%, 6.4%, and 4.3%, respectively, in independent 47 affinity-purified HG serous carcinomas. Except for the CDKN2A/B region, these homozygous deletions were not present in either SBT or LG tumors. Our study provides a genome-wide homozygous deletion profile in HG serous carcinomas, which can serve as a molecular foundation to study tumor suppressors in ovarian cancer.

Allelic imbalances and homozygous deletion on 8p23.2 for stepwise progression of hepatocarcinogenesis

Early hepatocellular carcinoma (eHCC) originates from the hepatocytes of chronic liver disease and develops into classical hepatocellular carcinoma (HCC). To identify sequential genetic changes in multistep hepatocarcinogenesis, we analyzed molecular karyotypes using oligonucleotide genotyping 50K arrays. First, 1q21.3-44 gain and loss of heterozygosity (LOH) on 1p36.21-36.32 and 17p13.1-13.3 were frequently observed in eHCC, but not in chronic liver diseases, suggesting that such chromosomal aberrations are early, possibly causative events in liver cancer. Next, we detected 25 chromosomal loci associated with liver cancer progression in five HCCs with nodule-in-nodule appearance, in which the inner nodule develops within eHCC lesion. Using these chromosomal regions as independent variables, decision tree analysis was applied on 14 early and 25 overt HCCs, and extracted combination of chromosomal gains on 5q11.1-35.3 and 8q11.1-24.3 and LOH on 4q11-34.3 and 8p11.21-23.3 as distinctive attributes, which can classify early and overt HCCs recursively. In these four altered regions identified as late events of hepatocarcinogenesis, two tumors in 32 overt HCCs analyzed in the present study and one in a set of independent samples of 36 overt HCCs in our previous study harbored a homozygous deletion near the CSMD1 locus on 8p23.2. CSMD1 messenger RNA expression was decreased in HCC without 8p23.2 deletion, possibly due to hypermethylation of the CpG islands in its promoter region.

CONCLUSION

1q gain and 1p and 17p LOH are early molecular events, whereas gains in 5q and 8q and LOH on 4q and 8p only occur in advanced HCC, and inactivation of the putative suppressor gene, CSMD1, may be the key event in progression of liver cancer.

High-resolution copy number and gene expression microarray analyses of head and neck squamous cell carcinoma cell lines of tongue and larynx

Gene amplifications and deletions are frequent in head and neck squamous cell carcinomas (SCC) but the association of these alterations with gene expression is mostly unknown. Here, we characterized genome-wide copy number and gene expression changes on microarrays for 18 oral tongue SCC (OTSCC) cell lines. We identified a number of altered regions including nine high-level amplifications such as 6q12-q14 (CD109, MYO6), 9p24 (JAK2, CD274, SLC1A1, RLN1), 11p12-p13 (TRAF6, COMMD9, TRIM44, FJX1, CD44, PDHX, APIP), 11q13 (FADD, PPFIA1, CTTN), and 14q24 (ABCD4, HBLD1, LTBP2, ZNF410, COQ6, ACYP1, JDP2) where 9% to 64% of genes showed overexpression. Across the whole genome, 26% of the amplified genes had associated overexpression in OTSCC. Furthermore, our data implicated that OTSCC cell lines harbored similar genomic alterations as laryngeal SCC cell lines We have previously analyzed, suggesting that despite differences in clinicopathological features there are no marked differences in molecular genetic alterations of these two HNSCC sites. To identify genes whose expression was associated with copy number increase in head and neck SCC, a statistical analysis for oral tongue and laryngeal SCC cell line data were performed. We pinpointed 1,192 genes that had a statistically significant association between copy number and gene expression. These results suggest that genomic alterations with associated gene expression changes play an important role in the malignant behavior of head and neck SCC. The identified genes provide a basis for further functional validation and may lead to the identification of novel candidates for targeted therapies. This article contains Supplementary Material available at http://www.interscience.wiley.com/jpages/1045-2257/suppmat.