Genome-wide scan identifies a copy number variable region at 3q26 that regulates PPM1L in APC mutation-negative familial colorectal cancer patients

Identification of Key MicroRNAs and Genes between Colorectal Adenoma and Colorectal Cancer via Deep Learning on GEO Databases and Bioinformatics

Background

Deep learning techniques are gaining momentum in medical research. Colorectal adenoma (CRA) is a precancerous lesion that may develop into colorectal cancer (CRC) and its etiology and pathogenesis are unclear. This study aims to identify transcriptome differences between CRA and CRC via deep learning on Gene Expression Omnibus (GEO) databases and bioinformatics in the Chinese population.

Methods

In this study, three microarray datasets from the GEO database were used to identify the differentially expressed genes (DEGs) and differentially expressed miRNAs (DEMs) in CRA and CRC. The FunRich software was performed to predict the targeted mRNAs of DEMs. The targeted mRNAs were overlapped with DEGs to determine the key DEGs. Molecular mechanisms of CRA and CRC were evaluated using enrichment analysis. Cytoscape was used to construct protein-protein interaction (PPI) and miRNA-mRNA regulatory networks. We analyzed the expression of key DEMs and DEGs, their prognosis, and correlation with immune infiltration based on the Kaplan-Meier plotter, UALCAN, and TIMER databases.

Results

A total of 38 DEGs are obtained after the intersection, including 11 upregulated genes and 27 downregulated genes. The DEGs were involved in the pathways, including epithelial-to-mesenchymal transition, sphingolipid metabolism, and intrinsic pathway for apoptosis. The expression of has-miR-34c (P = 0.036), hsa-miR-320a (P = 0.045), and has-miR-338 (P = 0.0063) was correlated with the prognosis of CRC patients. The expression levels of BCL2, PPM1L, ARHGAP44, and PRKACB in CRC tissues were significantly lower than normal tissues (P < 0.001), while the expression levels of TPD52L2 and WNK4 in CRC tissues were significantly higher than normal tissues (P < 0.01). These key genes are significantly associated with the immune infiltration of CRC.

Conclusion

This preliminary study will help identify patients with CRA and early CRC and establish prevention and monitoring strategies to reduce the incidence of CRC.

circFAM120B functions as a tumor suppressor in esophageal squamous cell carcinoma via the miR-661/PPM1L axis and the PKR/p38 MAPK/EMT pathway

Extensive changes of circRNA expression underscore their essential contributions to multiple hallmarks of cancers; however, their functions and mechanisms of action in esophageal squamous cell carcinoma (ESCC) remain undetermined. Here, we adopted a three-stage approach by first screening for significantly differentially expressed circRNAs in ESCC and performing an external validation study, followed by the functional analyses. The properties of circRNAs were evaluated using Sanger sequencing, RNase R digestion, actinomycin D treatment, subcellular localization analysis, and fluorescence in situ hybridization. Target transcripts were predicted using online tools and verified by dual-luciferase, RNA immunoprecipitation, qRT-PCR, and western blot. Biotin-labeled RNA-protein pull-down, mass spectrometry, and RNA immunoprecipitation were employed to identify proteins interacting with circRNAs. Gain- and loss-of-function experiments were performed to uncover the roles of circRNAs, their target genes, and binding proteins in the proliferation, metastasis, and invasion. We observed that circFAM120B (hsa_circ_0001666) was frequently downregulated in cancer tissues and patient plasma, and its expression level was related to overall survival in ESCC patients. Overexpression of circFAM120B inhibited the proliferation, metastasis, and invasion of ESCC while silencing it enhanced malignant phenotypes. Mechanistically, circFAM120B was predominantly located in the cytoplasm, guarantying its sponging for miR-661 to restore the expression of PPM1L, a tumor suppressor. We observed that circFAM120B could reduce the stability of RNA-dependent protein kinase (PKR) by promoting its ubiquitination-dependent degradation and subsequently regulating the p38 MAPK signaling pathway, resulting in the repression of EMTs in ESCC cells. Our findings suggest that circFAM120B is a promising biomarker of ESCC, which acts as a tumor suppressor via the circFAM120B/miR-661/PPM1L axis and PKR/p38 MAPK/EMT pathway, supporting its significance as a candidate therapeutic target.

Comprehensive analysis of transcriptome data for identifying biomarkers and therapeutic targets in head and neck squamous cell carcinoma

Background

Head and neck squamous cell carcinoma (HNSCC) is one of the most common malignancy worldwide. Accumulating evidences have highlighted the importance of transcriptome data during HNSCC tumorigenesis. The aim of this study was to identify significant genes as effective biomarkers for HNSCC and constructed miRNA-mRNA regulatory network for a more comprehensive understanding of the underlying molecular mechanisms.

Methods

A total of four independent microarrays conducted on HNSCC samples were downloaded from the Gene Expression Omnibus (GEO) and analyzed through R software. FunRich was applied to predict potential transcription factors and targeted genes of miRNAs. Protein-protein interaction (PPI) network and miRNA-mRNA regulatory network were constructed in Cytoscape. Additionally, the database for annotation, visualization, and integrated discovery (DAVID) was utilized to perform GO and KEGG pathway enrichment analyses. Validation of gene expression levels was conducted by online databases and qPCR experiments.

Results

A total of 35 and 193 differentially expressed miRNAs (DEMs) and mRNAs (DEGs) were screened out by the limma package in R. The interactive network of the overlapping DEGs presented three significant modules and ten hub genes (FN1, MMP3, SPP1, STAT1, LOX, CXCL5, CXCL11, ISG15, IFIT3, and RSAD2). Predicted target genes of DEMs were visualized in Cytoscape and six miRNA-mRNA regulatory pairs were identified. Further validation demonstrated the upregulation of SLC16A1 and COL4A1 in HNSCC.

Conclusions

We performed an integrated and comprehensive bioinformatics analysis of miRNAs and mRNAs in HNSCC, contributing to explore the underlying regulatory mechanisms and to identify genetic biomarkers and therapeutic targets for HNSCC.

A yeast phenomic model for the influence of Warburg metabolism on genetic buffering of doxorubicin

BACKGROUND

The influence of the Warburg phenomenon on chemotherapy response is unknown. Saccharomyces cerevisiae mimics the Warburg effect, repressing respiration in the presence of adequate glucose. Yeast phenomic experiments were conducted to assess potential influences of Warburg metabolism on gene-drug interaction underlying the cellular response to doxorubicin. Homologous genes from yeast phenomic and cancer pharmacogenomics data were analyzed to infer evolutionary conservation of gene-drug interaction and predict therapeutic relevance.

METHODS

Cell proliferation phenotypes (CPPs) of the yeast gene knockout/knockdown library were measured by quantitative high-throughput cell array phenotyping (Q-HTCP), treating with escalating doxorubicin concentrations under conditions of respiratory or glycolytic metabolism. Doxorubicin-gene interaction was quantified by departure of CPPs observed for the doxorubicin-treated mutant strain from that expected based on an interaction model. Recursive expectation-maximization clustering (REMc) and Gene Ontology (GO)-based analyses of interactions identified functional biological modules that differentially buffer or promote doxorubicin cytotoxicity with respect to Warburg metabolism. Yeast phenomic and cancer pharmacogenomics data were integrated to predict differential gene expression causally influencing doxorubicin anti-tumor efficacy.

RESULTS

Yeast compromised for genes functioning in chromatin organization, and several other cellular processes are more resistant to doxorubicin under glycolytic conditions. Thus, the Warburg transition appears to alleviate requirements for cellular functions that buffer doxorubicin cytotoxicity in a respiratory context. We analyzed human homologs of yeast genes exhibiting gene-doxorubicin interaction in cancer pharmacogenomics data to predict causality for differential gene expression associated with doxorubicin cytotoxicity in cancer cells. This analysis suggested conserved cellular responses to doxorubicin due to influences of homologous recombination, sphingolipid homeostasis, telomere tethering at nuclear periphery, actin cortical patch localization, and other gene functions.

CONCLUSIONS

Warburg status alters the genetic network required for yeast to buffer doxorubicin toxicity. Integration of yeast phenomic and cancer pharmacogenomics data suggests evolutionary conservation of gene-drug interaction networks and provides a new experimental approach to model their influence on chemotherapy response. Thus, yeast phenomic models could aid the development of precision oncology algorithms to predict efficacious cytotoxic drugs for cancer, based on genetic and metabolic profiles of individual tumors.

A Humanized Yeast Phenomic Model of Deoxycytidine Kinase to Predict Genetic Buffering of Nucleoside Analog Cytotoxicity

Knowledge about synthetic lethality can be applied to enhance the efficacy of anticancer therapies in individual patients harboring genetic alterations in their cancer that specifically render it vulnerable. We investigated the potential for high-resolution phenomic analysis in yeast to predict such genetic vulnerabilities by systematic, comprehensive, and quantitative assessment of drug–gene interaction for gemcitabine and cytarabine, substrates of deoxycytidine kinase that have similar molecular structures yet distinct antitumor efficacy. Human deoxycytidine kinase (dCK) was conditionally expressed in the Saccharomyces cerevisiae genomic library of knockout and knockdown (YKO/KD) strains, to globally and quantitatively characterize differential drug–gene interaction for gemcitabine and cytarabine. Pathway enrichment analysis revealed that autophagy, histone modification, chromatin remodeling, and apoptosis-related processes influence gemcitabine specifically, while drug–gene interaction specific to cytarabine was less enriched in gene ontology. Processes having influence over both drugs were DNA repair and integrity checkpoints and vesicle transport and fusion. Non-gene ontology (GO)-enriched genes were also informative. Yeast phenomic and cancer cell line pharmacogenomics data were integrated to identify yeast–human homologs with correlated differential gene expression and drug efficacy, thus providing a unique resource to predict whether differential gene expression observed in cancer genetic profiles are causal in tumor-specific responses to cytotoxic agents.

Genomic profiling of multiple sequentially acquired tumor metastatic sites from an "exceptional responder" lung adenocarcinoma patient reveals extensive genomic heterogeneity and novel somatic variants driving treatment response

We used next-generation sequencing to identify somatic alterations in multiple metastatic sites from an “exceptional responder” lung adenocarcinoma patient during his 7-yr course of ERBB2-directed therapies. The degree of heterogeneity was unprecedented, with ∼1% similarity between somatic alterations of the lung and lymph nodes. One novel translocation, PLAG1-ACTA2, present in both sites, up-regulated ACTA2 expression. ERBB2, the predominant driver oncogene, was amplified in both sites, more pronounced in the lung, and harbored an L869R mutation in the lymph node. Functional studies showed increased proliferation, migration, metastasis, and resistance to ERBB2-directed therapy because of L869R mutation and increased migration because of ACTA2 overexpression. Within the lung, a nonfunctional CDK12, due to a novel G879V mutation, correlated with down-regulation of DNA damage response genes, causing genomic instability, and sensitivity to chemotherapy. We propose a model whereby a subclone metastasized early from the primary site and evolved independently in lymph nodes.

Germline DNA copy number variations as potential prognostic markers for non-muscle invasive bladder cancer progression

Accumulating evidence has suggested that germline DNA copy number variations (CNVs) affect various disorders, including human malignancies. However, the significance of CNVs in non-muscle invasive bladder cancer (NMIBC) remains unclear. The purpose of the present study was to identify the role of CNVs in NMIBC. Array comparative genomic hybridization (CGH) analysis was performed to search for candidate CNVs associated with NMIBC susceptibility. Quantitative polymerase chain reaction was carried out to evaluate CNVs associated with patient outcome in 189 NMIBC cases. In total, 11 CNVs were associated with NMIBC risk in array CGH analysis. Out of the 189 CNVs examined, family with sequence similarity 81 member A (FAM81A) and proprotein convertase subtilisin/kexin type 6 (PCSK6) CNVs exhibited a significant association with recurrence and disease progression in NMIBC. PCSK6 has been reported to regulate proliferation and tumor progression in breast and prostate malignancies. Notably, patients with pT1 stage had significantly lower PCSK6 relative copy number than those with pTa (P=0.0196). In multivariate analyses, PCSK6 copy number was an independent prognostic factor for progression-free survival (P=0.0456; risk ratio, 2.17; 95% confidence interval, 1.02-4.82). These data suggest that PCSK6 CNV is a potential new tumor marker for estimating disease progression in NMIBC.

Chromosome 19q13 disruption alters expressions of CYP2A7, MIA and MIA-RAB4B lncRNA and contributes to FAP-like phenotype in APC mutation-negative familial colorectal cancer patients

Familial adenomatous polyposis (FAP) is an autosomal-dominantly inherited form of colorectal cancer (CRC) caused by mutation in the adenomatous polyposis coli (APC) gene. Our ability to exhaustively screen for APC mutations identify microsatellite-stable and APC-mutation negative familial CRC patients, enabling us to search for novel genes. We performed genome-wide scan on two affected siblings of one family and 88 ethnicity- and gender-matched healthy controls to identify deletions shared by the siblings. Combined loss of heterozygosity, copy number and allelic-specific copy number analysis uncovered 5 shared deletions. Long-range polymerase chain reaction (PCR) confirmed chromosome 19q13 deletion, which was subsequently found in one other family. The 32 kb deleted region harbors the CYP2A7 gene and was enriched with enhancer, repressor and insulator sites. The wildtype allele was lost in the polyps of the proband. Further, real-time RT-PCR assays showed that expressions of MIA and MIA-RAB4B located 35 kb upstream of the deletion, were up-regulated in the polyps compared to the matched mucosa of the proband. MIA-RAB4B, the read-through long non-coding RNA (lncRNA), RAB4B, PIM2 and TAOK1 share common binding site of a microRNA, miR-24, in their 3'UTRs. PIM2 and TAOK1, two target oncogenes of miR-24, were co-ordinately up-regulated with MIA-RAB4B in the polyps, suggesting that MIA-RAB4B could function as competitive endogenous RNA to titrate miR-24 away from its other targets. The data suggest that the 19.13 deletion disrupted chromatin boundary, leading to altered expression of several genes and lncRNA, could contribute to colorectal cancer via novel genetic and epigenetic mechanisms.

Cyclin E Deregulation and Genomic Instability

Precise replication of genetic material and its equal distribution to daughter cells are essential to maintain genome stability. In eukaryotes, chromosome replication and segregation are temporally uncoupled, occurring in distinct intervals of the cell cycle, S and M phases, respectively. Cyclin E accumulates at the G1/S transition, where it promotes S phase entry and progression by binding to and activating CDK2. Several lines of evidence from different models indicate that cyclin E/CDK2 deregulation causes replication stress in S phase and chromosome segregation errors in M phase, leading to genomic instability and cancer. In this chapter, we will discuss the main findings that link cyclin E/CDK2 deregulation to genomic instability and the molecular mechanisms by which cyclin E/CDK2 induces replication stress and chromosome aberrations during carcinogenesis.

Genomics of Colorectal Cancer in African Americans

Genome-wide studies are increasingly becoming a must, especially for complex diseases such as cancer where multiple genes and diverse molecular mechanisms are known to be involved in genes' function alteration. In this review, we report our latest genomic and epigenomic findings in African-American colorectal cancer patients. This population suffers a higher burden of the disease and most investigators in this field are looking for the underlying genetic and epigenetic targets that might be responsible for this disparity. We here report genome-wide copy number variations, single nucleotide mutations and DNA methylation findings that might be specific to this population.

Somatic gene copy number alterations in colorectal cancer: new quest for cancer drivers and biomarkers

Colorectal cancer (CRC) results from the accumulation of genetic alterations, and somatic copy number alterations (CNAs) are crucial for the development of CRC. Genome-wide survey of CNAs provides opportunities for identifying cancer driver genes in an unbiased manner. The detection of aberrant CNAs may provide novel markers for the early diagnosis and personalized treatment of CRC. A major challenge in array-based profiling of CNAs is to distinguish the alterations that play causative roles from the random alterations that accumulate during colorectal carcinogenesis. In this view, we systematically discuss the frequent CNAs in CRC, focusing on functional genes that have potential diagnostic, prognostic and therapeutic significance.

The association of copy number variation and percent mammographic density

Background

Percent mammographic density (PD) estimates the proportion of stromal, fat, and epithelial breast tissues on the mammogram image. Adjusted for age and body mass index (BMI), PD is one of the strongest risk factors for breast cancer [1]. Inherited factors are hypothesized to explain between 30 and 60% of the variance in this trait [2–5]. However, previously identified common genetic variants account for less than 6% of the variance in PD, leaving much of the genetic contribution to this trait unexplained. We performed the first study to examine whether germline copy number variation (CNV) are associated with PD. Two genome-wide association studies (GWAS) of percent density conducted on the Illumina 660W-Quad were used to identify and replicate the association between candidate CNVs and PD: the Minnesota Breast Cancer Family Study (MBCFS) and controls from the Mayo Venous Thromboembolism (Mayo VTE) Case–Control Study, with 585 and 328 women, respectively. Linear models were utilized to examine the association of each probe with PD, adjusted for age, menopausal status and BMI. Segmentation was subsequently performed on the probe-level test statistics to identify candidate CNV regions that were associated with PD.

Results

Sixty-one probes from five chromosomal regions [3q26.1 (2 regions), 8q24.22, 11p15.3, and 17q22] were significantly associated with PD in MBCFS (p-values <0.0001). A CNV at 3q26.1 showed the greatest evidence for association with PD; a region without any known SNPs. Conversely, the CNV at 17q22 was largely due to the association between SNPs and PD in the region. SNPs in the 8q24.22 region have been shown to be associated with risk of many cancers; however, SNPs in this region were not responsible for the observed CNV association. While we were unable to replicate the associations with PD, two of the five CNVs (3q26.1 and 11p15.3) were also observed in the Mayo VTE controls.

Conclusions

CNVs may help to explain some of the variability in PD that is currently unexplained by SNPs. While we were able to replicate the existence of two CNVs across the two GWAS studies, we were unable to replicate the associations with PD. Even so, the proximity of the identified CNV regions to loci known to be associated with breast cancer risk suggests further investigation and potentially shared genetic mechanisms underlying the PD and breast cancer association.

Electronic supplementary material

The online version of this article (doi:10.1186/s13104-015-1212-y) contains supplementary material, which is available to authorized users.

Identification of rare germline copy number variations over-represented in five human cancer types

BACKGROUND

Copy number variations (CNVs) are increasingly recognized as significant disease susceptibility markers in many complex disorders including cancer. The availability of a large number of chromosomal copy number profiles in both malignant and normal tissues in cancer patients presents an opportunity to characterize not only somatic alterations but also germline CNVs, which may confer increased risk for cancer.

RESULTS

We explored the germline CNVs in five cancer cohorts from the Cancer Genome Atlas (TCGA) consisting of 351 brain, 336 breast, 342 colorectal, 370 renal, and 314 ovarian cancers, genotyped on Affymetrix SNP6.0 arrays. Comparing these to ~3000 normal controls from another study, our case-control association study revealed 39 genomic loci (9 brain, 3 breast, 4 colorectal, 11 renal, and 12 ovarian cancers) as potential candidates of tumor susceptibility loci. Many of these loci are new and in some cases are associated with a substantial increase in disease risk. The majority of the observed loci do not overlap with coding sequences; however, several observed genomic loci overlap with known cancer genes including RET in brain cancers, ERBB2 in renal cell carcinomas, and DCC in ovarian cancers, all of which have not been previously associated with germline changes in cancer.

CONCLUSIONS

This large-scale genome-wide association study for CNVs across multiple cancer types identified several novel rare germline CNVs as cancer predisposing genomic loci. These loci can potentially serve as clinically useful markers conferring increased cancer risk.

Continuing difficulties in interpreting CNV data: lessons from a genome-wide CNV association study of Australian HNPCC/lynch syndrome patients

Background

Hereditary non-polyposis colorectal cancer (HNPCC)/Lynch syndrome (LS) is a cancer syndrome characterised by early-onset epithelial cancers, especially colorectal cancer (CRC) and endometrial cancer. The aim of the current study was to use SNP-array technology to identify genomic aberrations which could contribute to the increased risk of cancer in HNPCC/LS patients.

Methods

Individuals diagnosed with HNPCC/LS (100) and healthy controls (384) were genotyped using the Illumina Human610-Quad SNP-arrays. Copy number variation (CNV) calling and association analyses were performed using Nexus software, with significant results validated using QuantiSNP. TaqMan Copy-Number assays were used for verification of CNVs showing significant association with HNPCC/LS identified by both software programs.

Results

We detected copy number (CN) gains associated with HNPCC/LS status on chromosome 7q11.21 (28% cases and 0% controls, Nexus; p = 3.60E-20 and QuantiSNP; p < 1.00E-16) and 16p11.2 (46% in cases, while a CN loss was observed in 23% of controls, Nexus; p = 4.93E-21 and QuantiSNP; p = 5.00E-06) via in silico analyses. TaqMan Copy-Number assay was used for validation of CNVs showing significant association with HNPCC/LS. In addition, CNV burden (total CNV length, average CNV length and number of observed CNV events) was significantly greater in cases compared to controls.

Conclusion

A greater CNV burden was identified in HNPCC/LS cases compared to controls supporting the notion of higher genomic instability in these patients. One intergenic locus on chromosome 7q11.21 is possibly associated with HNPCC/LS and deserves further investigation. The results from this study highlight the complexities of fluorescent based CNV analyses. The inefficiency of both CNV detection methods to reproducibly detect observed CNVs demonstrates the need for sequence data to be considered alongside intensity data to avoid false positive results.

Genome wide analysis of chromosomal alterations in oral squamous cell carcinomas revealed over expression of MGAM and ADAM9

Despite the advances in diagnosis and treatment of oral squamous cell carcinoma (OSCC), mortality and morbidity rates have not improved over the past decade. A major drawback in diagnosis and treatment of OSCC is the lack of knowledge relating to how genetic instability in oral cancer genomes affects oral carcinogenesis. Hence, the key aim of this study was to identify copy number alterations (CNAs) that may be cancer associated in OSCC using high-resolution array comparative genomic hybridization (aCGH). To our knowledge this is the first study to use ultra-high density aCGH microarrays to profile a large number of OSCC genomes (n = 46). The most frequently amplified CNAs were located on chromosome 11q11(52%), 2p22.3(52%), 1q21.3-q22(54%), 6p21.32(59%), 20p13(61%), 7q34(52% and 72%),8p11.23-p11.22(80%), 8q11.1-q24.4(54%), 9q13-q34.3(54%), 11q23.3-q25(57%); 14q21.3-q31.1(54%); 14q31.3-q32.33(57%), 20p13-p12.3(54%) and 20q11.21-q13.33(52%). The most frequently deleted chromosome region was located on 3q26.1 (54%). In order to verify the CNAs from aCGH using quantitative polymerase chain reaction (qPCR), the three top most amplified regions and their associated genes, namely ADAM5P (8p11.23-p11.22), MGAM (7q34) and SIRPB1 (20p13.1), were selected in this study. The ADAM5P locus was found to be amplified in 39 samples and deleted in one; MGAM (24 amplifications and 3 deletions); and SIRPB1 (12 amplifications, others undetermined). On the basis of putative cancer-related annotations, two genes, namely ADAM metallopeptidase domain 9 (ADAM9) and maltase-glucoamylase alpha-glucosidase (MGAM), that mapped to CNA regions were selected for further evaluation of their mRNA expression using reverse transcriptase qPCR. The over-expression of MGAM was confirmed with a 6.6 fold increase in expression at the mRNA level whereas the fold change in ADAM9 demonstrated a 1.6 fold increase. This study has identified significant regions in the OSCC genome that were amplified and resulted in consequent over-expression of the MGAM and ADAM9 genes that may be utilized as biological markers for OSCC.

Germline copy number variations associated with breast cancer susceptibility in a Japanese population

Although copy number variations (CNVs) are expected to affect various diseases, little is known about the association between CNVs and breast cancer susceptibility. Therefore, we investigated this relation. Array comparative genomic hybridization was performed to search for candidate CNVs related to breast cancer susceptibility. Subsequent quantitative real-time polymerase chain reaction was carried out for confirmation. We found seven CNV markers associated with breast cancer risk. The means of the relative copy numbers of patients with a history of breast cancer and women in the control group were 0.8 and 1.8 for Hs06535529_cn on 1p36.12 (P < 0.0001), 2.9 and 2.2 for Hs03103056_cn on 3q26.1 (P < 0.0001), 1.2 and 1.8 for Hs03899300_cn on 15q26.3 (P < 0.0001), 1.0 and 1.5 for Hs03908783_cn on 15q26.3 (P < 0.0001), and 1.1 and 1.7 for Hs03898338_cn on 15q26.3 (P < 0.0001), respectively. Interestingly, nine or more copies of Hs04093415_cn on 22q12.3 were found only in 8/193 (4.1 %) patients with a history of breast cancer and in none of the controls (P = 0.0081). Similarly, 12 or more copies of Hs040908898_cn on 22q12.3 were found only in 7/193 (3.6 %) patients with a history of breast cancer and in none of the controls (P = 0.016). A combination of two CNVs resulted in 80.3 % sensitivity, 80.6 % specificity, 82.4 % positive predictive value, and 78.3 % negative predictive value for the prediction of breast cancer susceptibility. These findings may lead to a new means of risk assessment for breast cancer. Confirmatory studies using independent data sets are needed to support our findings.

Deletion polymorphism at chromosome 3q26.1 and oral squamous cell carcinoma

Several recent studies have investigated DNA instability in malignancies including deletions and duplications of part of the chromosome using array-based comparative genomic hybridization (CGH) analysis. Using the same approach on oral squamous cell carcinoma (OSCC) tissue samples, we found a frequent deletion at chromosome 3q26.1 in OSCC patients; this polymorphism showed a gene frequency of 0.293-0.368 in healthy volunteers (n=60) and 0.129-0.195 in OSCC patients (n=54). Detailed analysis around the polymorphic region revealed the deletion breakage point. A significant association of gene frequency for the deletion polymorphism between healthy volunteers and patients implicated genetic factors related to this polymorphism in the development of OSCC. Currently, no gene is predicted to lie within the 3,606-kbp region around the polymorphism. Thus, although a single-gene model could not explain the occurrence of OSCC, we believe that examining this polymorphism could be useful in identifying risk factors for OSCC.

Genomic aberrations in an African American colorectal cancer cohort reveals a MSI-specific profile and chromosome X amplification in male patients

Objective

DNA aberrations that cause colorectal cancer (CRC) occur in multiple steps that involve microsatellite instability (MSI) and chromosomal instability (CIN). Herein, we studied CRCs from AA patients for their CIN and MSI status.

Experimental Design

Array CGH was performed on 30 AA colon tumors. The MSI status was established. The CGH data from AA were compared to published lists of 41 TSG and oncogenes in Caucasians and 68 cancer genes, proposed via systematic sequencing for somatic mutations in colon and breast tumors. The patient-by-patient CGH profiles were organized into a maximum parsimony cladogram to give insights into the tumors' aberrations lineage.

Results

The CGH analysis revealed that CIN was independent of age, gender, stage or location. However, both the number and nature of aberrations seem to depend on the MSI status. MSI-H tumors clustered together in the cladogram. The chromosomes with the highest rates of CGH aberrations were 3, 5, 7, 8, 20 and X. Chromosome X was primarily amplified in male patients. A comparison with Caucasians revealed an overall similar aberration profile with few exceptions for the following genes; THRB, RAF1, LPL, DCC, XIST, PCNT, STS and genes on the 20q12-q13 cytoband. Among the 68 CAN genes, all showed some level of alteration in our cohort.

Conclusion

Chromosome X amplification in male patients with CRC merits follow-up. The observed CIN may play a distinctive role in CRC in AAs. The clustering of MSI-H tumors in global CGH data analysis suggests that chromosomal aberrations are not random.

Foxa1 and Foxa2 are essential for sexual dimorphism in liver cancer

Hepatocellular carcinoma (HCC) is sexually dimorphic in both rodents and humans, with significantly higher incidence in males, an effect that is dependent on sex hormones. The molecular mechanisms by which estrogens prevent and androgens promote liver cancer remain unclear. Here, we discover that sexually dimorphic HCC is completely reversed in Foxa1- and Foxa2-deficient mice after diethylnitrosamine-induced hepatocarcinogenesis. Coregulation of target genes by Foxa1/a2 and either the estrogen receptor (ERα) or the androgen receptor (AR) was increased during hepatocarcinogenesis in normal female or male mice, respectively, but was lost in Foxa1/2-deficient mice. Thus, both estrogen-dependent resistance to and androgen-mediated facilitation of HCC depend on Foxa1/2. Strikingly, single nucleotide polymorphisms at FOXA2 binding sites reduce binding of both FOXA2 and ERα to their targets in human liver and correlate with HCC development in women. Thus, Foxa factors and their targets are central for the sexual dimorphism of HCC.

Germline DNA copy number variation in familial and early-onset breast cancer

Introduction

Genetic factors predisposing individuals to cancer remain elusive in the majority of patients with a familial or clinical history suggestive of hereditary breast cancer. Germline DNA copy number variation (CNV) has recently been implicated in predisposition to cancers such as neuroblastomas as well as prostate and colorectal cancer. We evaluated the role of germline CNVs in breast cancer susceptibility, in particular those with low population frequencies (rare CNVs), which are more likely to cause disease."

Methods

Using whole-genome comparative genomic hybridization on microarrays, we screened a cohort of women fulfilling criteria for hereditary breast cancer who did not carry BRCA1/BRCA2 mutations.

Results

The median numbers of total and rare CNVs per genome were not different between controls and patients. A total of 26 rare germline CNVs were identified in 68 cancer patients, however, a proportion that was significantly different (P = 0.0311) from the control group (23 rare CNVs in 100 individuals). Several of the genes affected by CNV in patients and controls had already been implicated in cancer.

Conclusions

This study is the first to explore the contribution of germline CNVs to BRCA1/2-negative familial and early-onset breast cancer. The data suggest that rare CNVs may contribute to cancer predisposition in this small cohort of patients, and this trend needs to be confirmed in larger population samples.

Mapping of genetic abnormalities of primary tumours from metastatic CRC by high-resolution SNP arrays

Background

For years, the genetics of metastatic colorectal cancer (CRC) have been studied using a variety of techniques. However, most of the approaches employed so far have a relatively limited resolution which hampers detailed characterization of the common recurrent chromosomal breakpoints as well as the identification of small regions carrying genetic changes and the genes involved in them.

Methodology/Principal Findings

Here we applied 500K SNP arrays to map the most common chromosomal lesions present at diagnosis in a series of 23 primary tumours from sporadic CRC patients who had developed liver metastasis. Overall our results confirm that the genetic profile of metastatic CRC is defined by imbalanced gains of chromosomes 7, 8q, 11q, 13q, 20q and X together with losses of the 1p, 8p, 17p and 18q chromosome regions. In addition, SNP-array studies allowed the identification of small (<1.3 Mb) and extensive/large (>1.5 Mb) altered DNA sequences, many of which contain cancer genes known to be involved in CRC and the metastatic process. Detailed characterization of the breakpoint regions for the altered chromosomes showed four recurrent breakpoints at chromosomes 1p12, 8p12, 17p11.2 and 20p12.1; interestingly, the most frequently observed recurrent chromosomal breakpoint was localized at 17p11.2 and systematically targeted the FAM27L gene, whose role in CRC deserves further investigations.

Conclusions/Significance

In summary, in the present study we provide a detailed map of the genetic abnormalities of primary tumours from metastatic CRC patients, which confirm and extend on previous observations as regards the identification of genes potentially involved in development of CRC and the metastatic process.