DNA copy number changes in human malignant fibrous histiocytomas by array comparative genomic hybridisation

Background

Malignant fibrous histiocytomas (MFHs), or undifferentiated pleomorphic sarcomas, are in general high-grade tumours with extensive chromosomal aberrations. In order to identify recurrent chromosomal regions of gain and loss, as well as novel gene targets of potential importance for MFH development and/or progression, we have analysed DNA copy number changes in 33 MFHs using microarray-based comparative genomic hybridisation (array CGH).

Principal findings

In general, the tumours showed numerous gains and losses of large chromosomal regions. The most frequent minimal recurrent regions of gain were 1p33-p32.3, 1p31.3-p31.2 and 1p21.3 (all gained in 58% of the samples), as well as 1q21.2-q21.3 and 20q13.2 (both 55%). The most frequent minimal recurrent regions of loss were 10q25.3-q26.11, 13q13.3-q14.2 and 13q14.3-q21.1 (all lost in 64% of the samples), as well as 2q36.3-q37.2 (61%), 1q41 (55%) and 16q12.1-q12.2 (52%). Statistical analyses revealed that gain of 1p33-p32.3 and 1p21.3 was significantly associated with better patient survival (P = 0.021 and 0.046, respectively). Comparison with similar array CGH data from 44 leiomyosarcomas identified seven chromosomal regions; 1p36.32-p35.2, 1p21.3-p21.1, 1q32.1-q42.13, 2q14.1-q22.2, 4q33-q34.3, 6p25.1-p21.32 and 7p22.3-p13, which were significantly different in copy number between the MFHs and leiomyosarcomas.

Conclusions

A number of recurrent regions of gain and loss have been identified, some of which were associated with better patient survival. Several specific chromosomal regions with significant differences in copy number between MFHs and leiomyosarcomas were identified, and these aberrations may be used as additional tools for the differential diagnosis of MFHs and leiomyosarcomas.     

Loss of 4q21.23-22.1 Is a Prognostic Marker for Disease Free and Overall Survival in Non-Small Cell Lung Cancer

This study was performed to assess the prognostic relevance of genomic aberrations at chromosome 4q in NSCLC patients. We have previously identified copy number changes at 4q12-q32 to be significantly associated with the early hematogenous dissemination of non-small cell lung cancer (NSCLC), and now aim to narrow down potential hot-spots within this 107 Mb spanning region. Using eight microsatellite markers at position 4q12-35, allelic imbalance (AI) analyses were performed on a preliminary study cohort (n = 86). Positions indicating clinicopathological and prognostic associations in AI analyses were further validated in a larger study cohort using fluorescence in situ hybridization (FISH) in 209 NSCLC patients. Losses at positions 4q21.23 and 4q22.1 were shown to be associated with advanced clinicopathological characteristics as well as with shortened disease free (DFS) and overall survival (OS) (DFS: P = 0.019; OS: P = 0.002). Multivariate analyses identified the losses of 4q21.23-22.1 to be an independent prognostic marker for both DFS and OS in NSCLC (HR 1.64–2.20, all P<0.04), and especially in squamous cell lung cancer (P<0.05). A case report study of a lung cancer patient further revealed a loss of 4q21.23 in disseminated tumor cells (DTCs). Neither gains at the latter positions, nor genomic aberrations at 4q12, 4q31.2 and 4q35.1, indicated a prognostic relevance. In conclusion, our data indicate that loss at 4q21.23-22.1 in NSCLC is of prognostic relevance in NSCLC patients and thus, includes potential new tumor suppressor genes with clinical relevance.     

A Genome-Wide Association Study of Optic Disc Parameters

The optic nerve head is involved in many ophthalmic disorders, including common diseases such as myopia and open-angle glaucoma. Two of the most important parameters are the size of the optic disc area and the vertical cup-disc ratio (VCDR). Both are highly heritable but genetically largely undetermined. We performed a meta-analysis of genome-wide association (GWA) data to identify genetic variants associated with optic disc area and VCDR. The gene discovery included 7,360 unrelated individuals from the population-based Rotterdam Study I and Rotterdam Study II cohorts. These cohorts revealed two genome-wide significant loci for optic disc area, rs1192415 on chromosome 1p22 (p = 6.72×10⁻¹⁹) within 117 kb of the CDC7 gene and rs1900004 on chromosome 10q21.3-q22.1 (p = 2.67×10⁻³³) within 10 kb of the ATOH7 gene. They revealed two genome-wide significant loci for VCDR, rs1063192 on chromosome 9p21 (p = 6.15×10⁻¹¹) in the CDKN2B gene and rs10483727 on chromosome 14q22.3-q23 (p = 2.93×10⁻¹⁰) within 40 kbp of the SIX1 gene. Findings were replicated in two independent Dutch cohorts (Rotterdam Study III and Erasmus Rucphen Family study; N = 3,612), and the TwinsUK cohort (N = 843). Meta-analysis with the replication cohorts confirmed the four loci and revealed a third locus at 16q12.1 associated with optic disc area, and four other loci at 11q13, 13q13, 17q23 (borderline significant), and 22q12.1 for VCDR. ATOH7 was also associated with VCDR independent of optic disc area. Three of the loci were marginally associated with open-angle glaucoma. The protein pathways in which the loci of optic disc area are involved overlap with those identified for VCDR, suggesting a common genetic origin.     

Classification and regression tree and spatial analyses reveal geographic heterogeneity in genome wide linkage study of Indian visceral leishmaniasis

Background

Genome wide linkage studies (GWLS) have provided evidence for loci controlling visceral leishmaniasis on Chromosomes 1p22, 6q27, 22q12 in Sudan and 6q27, 9p21, 17q11-q21 in Brazil. Genome wide studies from the major focus of disease in India have not previously been reported.

Methods and Findings

We undertook a GWLS in India in which a primary ∼10 cM (515 microsatellites) scan was carried out in 58 multicase pedigrees (74 nuclear families; 176 affected, 353 total individuals) and replication sought in 79 pedigrees (102 nuclear families; 218 affected, 473 total individuals). The primary scan provided evidence (≥2 adjacent markers allele-sharing LOD≥0.59; nominal P≤0.05) for linkage on Chromosomes 2, 5, 6, 7, 8, 10, 11, 20 and X, with peaks at 6p25.3-p24.3 and 8p23.1-p21.3 contributed to largely by 31 Hindu families and at Xq21.1-q26.1 by 27 Muslim families. Refined mapping confirmed linkage across all primary scan families at 2q12.2-q14.1 and 11q13.2-q23.3, but only 11q13.2-q23.3 replicated (combined LOD = 1.59; P = 0.0034). Linkage at 6p25.3-p24.3 and 8p23.1-p21.3, and at Xq21.1-q26.1, was confirmed by refined mapping for primary Hindu and Muslim families, respectively, but only Xq21.1-q26.1 replicated across all Muslim families (combined LOD 1.49; P = 0.0045). STRUCTURE and SMARTPCA did not identify population genetic substructure related to religious group. Classification and regression tree, and spatial interpolation, analyses confirm geographical heterogeneity for linkages at 6p25.3-p24.3, 8p23.1-p21.3 and Xq21.1-q26.1, with specific clusters of families contributing LOD scores of 2.13 (P = 0.0009), 1.75 (P = 0.002) and 1.84 (P = 0.001), respectively.

Conclusions

GWLS has identified novel loci that show geographical heterogeneity in their influence on susceptibility to VL in India.     

Expression of the p53 Target Wig-1 Is Associated with HPV Status and Patient Survival in Cervical Carcinoma

The p53 target gene WIG-1 (ZMAT3) is located in chromosomal region 3q26, that is frequently amplified in human tumors, including cervical cancer. We have examined the status of WIG-1 and the encoded Wig-1 protein in cervical carcinoma cell lines and tumor tissue samples. Our analysis of eight cervical cancer lines (Ca Ski, ME-180, MS751, SiHa, SW756, C-4I, C-33A, and HT-3) by spectral karyotype, comparative genomic hybridization and Southern blotting revealed WIG-1 is not the primary target for chromosome 3 gains. However, WIG-1/Wig-1 were readily expressed and WIG-1 mRNA expression was higher in the two HPV-negative cervical cell lines (C33-A, HT-3) than in HPV-positive lines. We then assessed Wig-1 expression by immunohistochemistry in 38 cervical tumor samples. We found higher nuclear Wig-1 expression levels in HPV-negative compared to HPV positive cases (p = 0.002) and in adenocarcinomas as compared to squamous cell lesions (p<0.0001). Cases with moderate nuclear Wig-1 staining and positive cytoplasmic Wig-1 staining showed longer survival than patients with strong nuclear and negative cytoplasmic staining (p = 0.042). Nuclear Wig-1 expression levels were positively associated with age at diagnosis (p = 0.023) and histologic grade (p = 0.034). These results are consistent with a growth-promoting and/or anti-cell death function of nuclear Wig-1 and suggest that Wig-1 expression can serve as a prognostic marker in cervical carcinoma.     

Genomic profiling of advanced-stage oral cancers reveals chromosome 11q alterations as markers of poor clinical outcome

Identifying oral cancer lesions associated with high risk of relapse and predicting clinical outcome remain challenging questions in clinical practice. Genomic alterations may add prognostic information and indicate biological aggressiveness thereby emphasizing the need for genome-wide profiling of oral cancers. High-resolution array comparative genomic hybridization was performed to delineate the genomic alterations in clinically annotated primary gingivo-buccal complex and tongue cancers (n = 60). The specific genomic alterations so identified were evaluated for their potential clinical relevance. Copy-number changes were observed on chromosomal arms with most frequent gains on 3q (60%), 5p (50%), 7p (50%), 8q (73%), 11q13 (47%), 14q11.2 (47%), and 19p13.3 (58%) and losses on 3p14.2 (55%) and 8p (83%). Univariate statistical analysis with correction for multiple testing revealed chromosomal gain of region 11q22.1–q22.2 and losses of 17p13.3 and 11q23–q25 to be associated with loco-regional recurrence (P = 0.004, P = 0.003, and P = 0.0003) and shorter survival (P = 0.009, P = 0.003, and P 0.0001) respectively. The gain of 11q22 and loss of 11q23-q25 were validated by interphase fluorescent in situ hybridization (I-FISH). This study identifies a tractable number of genomic alterations with few underlying genes that may potentially be utilized as biological markers for prognosis and treatment decisions in oral cancers.     

Both Rare and De Novo Copy Number Variants Are Prevalent in Agenesis of the Corpus Callosum but Not in Cerebellar Hypoplasia or Polymicrogyria

Agenesis of the corpus callosum (ACC), cerebellar hypoplasia (CBLH), and polymicrogyria (PMG) are severe congenital brain malformations with largely undiscovered causes. We conducted a large-scale chromosomal copy number variation (CNV) discovery effort in 255 ACC, 220 CBLH, and 147 PMG patients, and 2,349 controls. Compared to controls, significantly more ACC, but unexpectedly not CBLH or PMG patients, had rare genic CNVs over one megabase (p = 1.48×10⁻³; odds ratio [OR] = 3.19; 95% confidence interval [CI] = 1.89–5.39). Rare genic CNVs were those that impacted at least one gene in less than 1% of the combined population of patients and controls. Compared to controls, significantly more ACC but not CBLH or PMG patients had rare CNVs impacting over 20 genes (p = 0.01; OR = 2.95; 95% CI = 1.69–5.18). Independent qPCR confirmation showed that 9.4% of ACC patients had de novo CNVs. These, in comparison to inherited CNVs, preferentially overlapped de novo CNVs previously observed in patients with autism spectrum disorders (p = 3.06×10⁻⁴; OR = 7.55; 95% CI = 2.40–23.72). Interestingly, numerous reports have shown a reduced corpus callosum area in autistic patients, and diminished social and executive function in many ACC patients. We also confirmed and refined previously known CNVs, including significantly narrowing the 8p23.1-p11.1 duplication present in 2% of our current ACC cohort. We found six novel CNVs, each in a single patient, that are likely deleterious: deletions of 1p31.3-p31.1, 1q31.2-q31.3, 5q23.1, and 15q11.2-q13.1; and duplications of 2q11.2-q13 and 11p14.3-p14.2. One ACC patient with microcephaly had a paternally inherited deletion of 16p13.11 that included NDE1. Exome sequencing identified a recessive maternally inherited nonsense mutation in the non-deleted allele of NDE1, revealing the complexity of ACC genetics. This is the first systematic study of CNVs in congenital brain malformations, and shows a much higher prevalence of large gene-rich CNVs in ACC than in CBLH and PMG.     

Selective genomic copy number imbalances and probability of recurrence in early-stage breast cancer

A number of studies of copy number imbalances (CNIs) in breast tumors support associations between individual CNIs and patient outcomes. However, no pattern or signature of CNIs has emerged for clinical use. We determined copy number (CN) gains and losses using high-density molecular inversion probe (MIP) arrays for 971 stage I/II breast tumors and applied a boosting strategy to fit hazards models for CN and recurrence, treating chromosomal segments in a dose-specific fashion (-1 [loss], 0 [no change] and +1 [gain]). The concordance index (C-Index) was used to compare prognostic accuracy between a training (n = 728) and test (n = 243) set and across models. Twelve novel prognostic CNIs were identified: losses at 1p12, 12q13.13, 13q12.3, 22q11, and Xp21, and gains at 2p11.1, 3q13.12, 10p11.21, 10q23.1, 11p15, 14q13.2-q13.3, and 17q21.33. In addition, seven CNIs previously implicated as prognostic markers were selected: losses at 8p22 and 16p11.2 and gains at 10p13, 11q13.5, 12p13, 20q13, and Xq28. For all breast cancers combined, the final full model including 19 CNIs, clinical covariates, and tumor marker-approximated subtypes (estrogen receptor [ER], progesterone receptor, ERBB2 amplification, and Ki67) significantly outperformed a model containing only clinical covariates and tumor subtypes (C-Index _{full model}, train[test]  =  0.72[0.71] ± 0.02 vs. C-Index _{clinical + subtype model}, train[test]  =  0.62[0.62] ± 0.02; p<10⁻⁶). In addition, the full model containing 19 CNIs significantly improved prognostication separately for ER–, HER2+, luminal B, and triple negative tumors over clinical variables alone. In summary, we show that a set of 19 CNIs discriminates risk of recurrence among early-stage breast tumors, independent of ER status. Further, our data suggest the presence of specific CNIs that promote and, in some cases, limit tumor spread.     

Analysis of Molecular Cytogenetic Alteration in Rhabdomyosarcoma by Array Comparative Genomic Hybridization

Rhabdomyosarcoma (RMS) is the most common pediatric soft tissue sarcoma with poor prognosis. The genetic etiology of RMS remains largely unclear underlying its development and progression. To reveal novel genes more precisely and new therapeutic targets associated with RMS, we used high-resolution array comparative genomic hybridization (aCGH) to explore tumor-associated copy number variations (CNVs) and genes in RMS. We confirmed several important genes by quantitative real-time polymerase chain reaction (QRT-PCR). We then performed bioinformatics-based functional enrichment analysis for genes located in the genomic regions with CNVs. In addition, we identified miRNAs located in the corresponding amplification and deletion regions and performed miRNA functional enrichment analysis. aCGH analyses revealed that all RMS showed specific gains and losses. The amplification regions were 12q13.12, 12q13.3, and 12q13.3–q14.1. The deletion regions were 1p21.1, 2q14.1, 5q13.2, 9p12, and 9q12. The recurrent regions with gains were 12q13.3, 12q13.3–q14.1, 12q14.1, and 17q25.1. The recurrent regions with losses were 9p12–p11.2, 10q11.21–q11.22, 14q32.33, 16p11.2, and 22q11.1. The mean mRNA level of GLI1 in RMS was 6.61-fold higher than that in controls (p = 0.0477) by QRT-PCR. Meanwhile, the mean mRNA level of GEFT in RMS samples was 3.92-fold higher than that in controls (p = 0.0354). Bioinformatic analysis showed that genes were enriched in functions such as immunoglobulin domain, induction of apoptosis, and defensin. Proto-oncogene functions were involved in alveolar RMS. miRNAs that located in the amplified regions in RMS tend to be enriched in oncogenic activity (miR-24 and miR-27a). In conclusion, this study identified a number of CNVs in RMS and functional analyses showed enrichment for genes and miRNAs located in these CNVs regions. These findings may potentially help the identification of novel biomarkers and/or drug targets implicated in diagnosis of and targeted therapy for RMS.     

Genetic Variant rs7758229 in 6q26–q27 Is Not Associated with Colorectal Cancer Risk in a Chinese Population

Background

A recent genome-wide association study has identified a new genetic variant rs7758229 in SLC22A3 for colorectal cancer susceptibility in a Japanese population, but it is unknown whether this newly identified variant is associated with colorectal cancer in other populations, including the Chinese population.

Methods

We examined the associations between rs7758229 and colorectal cancer risk among 1,147 cases and 1,203 controls matched by age and sex. Logistic regression model was used to assess the associations.

Results

No significant association was found between rs7758229 and colorectal cancer risk (OR = 0.95, 95%CI  = 0.84–1.09, P = 0.463). Similar results were observed in the stratification of tumor location (OR  = 0.94, 95%CI = 0.80–1.11, P = 0.481 for colon cancer, and OR  = 0.96, 95%CI  = 0.82–1.13, P = 0.621 for rectum cancer).

Conclusions

Our findings did not support an association between rs7758229 in 6q26-q27 and the risk of colorectal cancer in a Chinese population.     

Dual-Color Fluorescence In Situ Hybridization Reveals an Association of Chromosome 8q22 but Not 8p21 Imbalance with High Grade Invasive Breast Carcinoma

We previously reported molecular karyotype analysis of invasive breast tumour core needle biopsies by comparative genomic hybridization (CGH) and fluorescence in situ hybridization (FISH) (Walker et al, Genes Chromosomes Cancer, 2008 May;47(5):405-17). That study identified frequently recurring gains and losses involving chromosome bands 8q22 and 8p21, respectively. Moreover, these data highlighted an association between 8q22 gain and typically aggressive grade 3 tumors. Here we validate and extend our previous investigations through FISH analysis of tumor touch imprints prepared from excised breast tumor specimens. Compared to post-surgical tumor excisions, core needle biopsies are known to be histologically less precise when predicting tumor grade. Therefore investigating these chromosomal aberrations in tumor samples that offer more reliable pathological assessment is likely to give a better overall indication of association. A series of 60 breast tumors were screened for genomic copy number changes at 8q22 and 8p21 by dual-color FISH. Results confirm previous findings that 8p loss (39%) and 8q gain (74%) occur frequently in invasive breast cancer. Both absolute quantification of 8q22 gain across the sample cohort, and a separate relative assessment by 8q22:8p21 copy number ratio, showed that the incidence of 8q22 gain significantly increased with grade (p = 0.004, absolute and p = 0.02, relative). In contrast, no association was found between 8p21 loss and tumor grade. These findings support the notion that 8q22 is a region of interest for invasive breast cancer pathogenesis, potentially harboring one or more genes that, when amplified, precipitate the molecular events that define high tumor grade.     

Loss of Heterozygosity and Copy Number Alterations in Flow-Sorted Bulky Cervical Cancer

Treatment choices for cervical cancer are primarily based on clinical FIGO stage and the post-operative evaluation of prognostic parameters including tumor diameter, parametrial and lymph node involvement, vaso-invasion, infiltration depth, and histological type. The aim of this study was to evaluate genomic changes in bulky cervical tumors and their relation to clinical parameters, using single nucleotide polymorphism (SNP)-analysis.Flow-sorted tumor cells and patient-matched normal cells were extracted from 81 bulky cervical tumors. DNA-index (DI) measurement and whole genome SNP-analysis were performed. Data were analyzed to detect copy number alterations (CNA) and allelic balance state: balanced, imbalanced or pure LOH, and their relation to clinical parameters.The DI varied from 0.92–2.56. Pure LOH was found in ≥40% of samples on chromosome-arms 3p, 4p, 6p, 6q, and 11q, CN gains in >20% on 1q, 3q, 5p, 8q, and 20q, and losses on 2q, 3p, 4p, 11q, and 13q. Over 40% showed gain on 3q. The only significant differences were found between histological types (squamous, adeno and adenosquamous) in the lesser allele intensity ratio (LAIR) (p = 0.035) and in the CNA analysis (p = 0.011). More losses were found on chromosome-arm 2q (FDR = 0.004) in squamous tumors and more gains on 7p, 7q, and 9p in adenosquamous tumors (FDR = 0.006, FDR = 0.004, and FDR = 0.029).Whole genome analysis of bulky cervical cancer shows widespread changes in allelic balance and CN. The overall genetic changes and CNA on specific chromosome-arms differed between histological types. No relation was found with the clinical parameters that currently dictate treatment choice.     

Gene Expression Profiling in Leiomyosarcomas and Undifferentiated Pleomorphic Sarcomas: SRC as a New Diagnostic Marker

Background

Undifferentiated Pleomorphic Sarcoma (UPS) and high-grade Leiomyosarcoma (LMS) are soft tissue tumors with an aggressive clinical behavior, frequently developing local recurrence and distant metastases. Despite several gene expression studies involving soft tissue sarcomas, the potential to identify molecular markers has been limited, mostly due to small sample size, in-group heterogeneity and absence of detailed clinical data.

Materials and Methods

Gene expression profiling was performed for 22 LMS and 22 UPS obtained from untreated patients. To assess the relevance of the gene signature, a meta-analysis was performed using five published studies. Four genes (BAD, MYOCD, SRF and SRC) selected from the gene signature, meta-analysis and functional in silico analysis were further validated by quantitative PCR. In addition, protein-protein interaction analysis was applied to validate the data. SRC protein immunolabeling was assessed in 38 UPS and 52 LMS.

Results

We identified 587 differentially expressed genes between LMS and UPS, of which 193 corroborated with other studies. Cluster analysis of the data failed to discriminate LMS from UPS, although it did reveal a distinct molecular profile for retroperitoneal LMS, which was characterized by the over-expression of smooth muscle-specific genes. Significantly higher levels of expression for BAD, SRC, SRF, and MYOCD were confirmed in LMS when compared with UPS. SRC was the most value discriminator to distinguish both sarcomas and presented the highest number of interaction in the in silico protein-protein analysis. SRC protein labeling showed high specificity and a positive predictive value therefore making it a candidate for use as a diagnostic marker in LMS.

Conclusions

Retroperitoneal LMS presented a unique gene signature. SRC is a putative diagnostic marker to differentiate LMS from UPS.     

Adjusting Breast Cancer Patient Prognosis with Non-HER2-Gene Patterns on Chromosome 17

Background

HER2 and TOP2A gene status are assessed for diagnostic and research purposes in breast cancer with fluorescence in situ hybridization (FISH). However, FISH probes do not target only the annotated gene, while chromosome 17 (chr17) is among the most unstable chromosomes in breast cancer. Here we asked whether the status of specifically targeted genes on chr17 might help in refining prognosis of early high-risk breast cancer patients.

Methods

Copy numbers (CN) for 14 genes on chr17, 4 of which were within and 10 outside the core HER2 amplicon (HER2- and non-HER2-genes, respectively) were assessed with qPCR in 485 paraffin-embedded tumor tissue samples from breast cancer patients treated with adjuvant chemotherapy in the frame of two randomized phase III trials.

Principal Findings

HER2-genes CN strongly correlated to each other (Spearman’s rho >0.6) and were concordant with FISH HER2 status (Kappa 0.6697 for ERBB2 CN). TOP2A CN were not concordant with TOP2A FISH status (Kappa 0.1154). CN hierarchical clustering revealed distinct patterns of gains, losses and complex alterations in HER2- and non-HER2-genes associated with IHC4 breast cancer subtypes. Upon multivariate analysis, non-HER2-gene gains independently predicted for shorter disease-free survival (DFS) and overall survival (OS) in patients with triple-negative cancer, as compared to luminal and HER2-positive tumors (interaction p = 0.007 for DFS and p = 0.011 for OS). Similarly, non-HER2-gene gains were associated with worse prognosis in patients who had undergone breast-conserving surgery as compared to modified radical mastectomy (p = 0.004 for both DFS and OS). Non-HER2-gene losses were unfavorable prognosticators in patients with 1–3 metastatic nodes, as compared to those with 4 or more nodes (p = 0.017 for DFS and p = 0.001 for OS).

Conclusions

TOP2A FISH and qPCR may not identify the same pathology on chr17q. Non-HER2 chr17 CN patterns may further predict outcome in breast cancer patients with known favorable and unfavorable prognosis.     

Germline DNA Copy Number Aberrations Identified as Potential Prognostic Factors for Breast Cancer Recurrence

Breast cancer recurrence (BCR) is a common treatment outcome despite curative-intent primary treatment of non-metastatic breast cancer. Currently used prognostic and predictive factors utilize tumor-based markers, and are not optimal determinants of risk of BCR. Germline-based copy number aberrations (CNAs) have not been evaluated as determinants of predisposition to experience BCR. In this study, we accessed germline DNA from 369 female breast cancer subjects who received curative-intent primary treatment following diagnosis. Of these, 155 experienced BCR and 214 did not, after a median duration of follow up after breast cancer diagnosis of 6.35 years (range = 0.60–21.78) and 8.60 years (range = 3.08–13.57), respectively. Whole genome CNA genotyping was performed on the Affymetrix SNP array 6.0 platform. CNAs were identified using the SNP-Fast Adaptive States Segmentation Technique 2 algorithm implemented in Nexus Copy Number 6.0. Six samples were removed due to poor quality scores, leaving 363 samples for further analysis. We identified 18,561 CNAs with ≥1 kb as a predefined cut-off for observed aberrations. Univariate survival analyses (log-rank tests) identified seven CNAs (two copy number gains and five copy neutral-loss of heterozygosities, CN-LOHs) showing significant differences (P<2.01×10⁻⁵) in recurrence-free survival (RFS) probabilities with and without CNAs.We also observed three additional but distinct CN-LOHs showing significant differences in RFS probabilities (P<2.86×10⁻⁵) when analyses were restricted to stratified cases (luminal A, n = 208) only. After adjusting for tumor stage and grade in multivariate analyses (Cox proportional hazards models), all the CNAs remained strongly associated with the phenotype of BCR. Of these, we confirmed three CNAs at 17q11.2, 11q13.1 and 6q24.1 in representative samples using independent genotyping platforms. Our results suggest further investigations on the potential use of germline DNA variations as prognostic markers in cancer-associated phenotypes.     

Analysis of Gene Expression Profiles of Soft Tissue Sarcoma Using a Combination of Knowledge-Based Filtering with Integration of Multiple Statistics

The diagnosis and treatment of soft tissue sarcomas (STS) have been difficult. Of the diverse histological subtypes, undifferentiated pleomorphic sarcoma (UPS) is particularly difficult to diagnose accurately, and its classification per se is still controversial. Recent advances in genomic technologies provide an excellent way to address such problems. However, it is often difficult, if not impossible, to identify definitive disease-associated genes using genome-wide analysis alone, primarily because of multiple testing problems. In the present study, we analyzed microarray data from 88 STS patients using a combination method that used knowledge-based filtering and a simulation based on the integration of multiple statistics to reduce multiple testing problems. We identified 25 genes, including hypoxia-related genes (e.g., MIF, SCD1, P4HA1, ENO1, and STAT1) and cell cycle- and DNA repair-related genes (e.g., TACC3, PRDX1, PRKDC, and H2AFY). These genes showed significant differential expression among histological subtypes, including UPS, and showed associations with overall survival. STAT1 showed a strong association with overall survival in UPS patients (logrank p = 1.84×10⁻⁶ and adjusted p value 2.99×10⁻³ after the permutation test). According to the literature, the 25 genes selected are useful not only as markers of differential diagnosis but also as prognostic/predictive markers and/or therapeutic targets for STS. Our combination method can identify genes that are potential prognostic/predictive factors and/or therapeutic targets in STS and possibly in other cancers. These disease-associated genes deserve further preclinical and clinical validation.     

Metastasis associated genomic aberrations in stage II rectal cancer

Genomic aberrations of rectal carcinoma, especially DNA copy number changes associated with metastasis were largely unclear. We aim to identify the metastasis associated biomarkers in stage II rectal cancer. Formalin-fixed, paraffin-embedded primary tumor tissues of stage II rectal carcinoma were analyzed by array-based comparative genomic hybridization, and genomic aberrations were identified by Genomic Workbench and SAM software. Copy number changes and mRNA expressions were validated by Real-time PCR in an independent rectal cancer samples. The results showed that the most frequent gains in stage II rectal cancer were at 1q21.2-q23.1, 3p21.31, 11q12.2-q23.3, 12q24.11-q24.31, 12q13.11-q14.1 and losses in 18q11.2-q23, 17q21.33-q22, 13q31.1-q31.3, 21q21.1-q21.3, 8p23.3-p23.1 and 4q22.1-q23. Twenty-two amplifications and five homozygous deletions were also identified. We further found that S100A1 (1q21.3-q23.1), MCM7 (7q22.1) and JUND (19p13.11) were amplified and overexpressed in stage II rectal cancer. Interestingly, the genomic aberrations affected 14 signaling pathways including VEGF signaling pathway and fatty acid metabolism. Most importantly, loss of 13q31.1-q34 and gain of 1q44 were associated with distant metastasis. Our results indicated that these metastasis associated genomic changes may be useful to reveal the pathogenesis of rectal cancer metastasis and identify candidate biomarkers.     

Genetic Predisposition to In Situ and Invasive Lobular Carcinoma of the Breast

Invasive lobular breast cancer (ILC) accounts for 10–15% of all invasive breast carcinomas. It is generally ER positive (ER+) and often associated with lobular carcinoma in situ (LCIS). Genome-wide association studies have identified more than 70 common polymorphisms that predispose to breast cancer, but these studies included predominantly ductal (IDC) carcinomas. To identify novel common polymorphisms that predispose to ILC and LCIS, we pooled data from 6,023 cases (5,622 ILC, 401 pure LCIS) and 34,271 controls from 36 studies genotyped using the iCOGS chip. Six novel SNPs most strongly associated with ILC/LCIS in the pooled analysis were genotyped in a further 516 lobular cases (482 ILC, 36 LCIS) and 1,467 controls. These analyses identified a lobular-specific SNP at 7q34 (rs11977670, OR (95%CI) for ILC = 1.13 (1.09–1.18), P = 6.0×10⁻¹⁰; P-het for ILC vs IDC ER+ tumors = 1.8×10⁻⁴). Of the 75 known breast cancer polymorphisms that were genotyped, 56 were associated with ILC and 15 with LCIS at P<0.05. Two SNPs showed significantly stronger associations for ILC than LCIS (rs2981579/10q26/FGFR2, P-het = 0.04 and rs889312/5q11/MAP3K1, P-het = 0.03); and two showed stronger associations for LCIS than ILC (rs6678914/1q32/LGR6, P-het = 0.001 and rs1752911/6q14, P-het = 0.04). In addition, seven of the 75 known loci showed significant differences between ER+ tumors with IDC and ILC histology, three of these showing stronger associations for ILC (rs11249433/1p11, rs2981579/10q26/FGFR2 and rs10995190/10q21/ZNF365) and four associated only with IDC (5p12/rs10941679; rs2588809/14q24/RAD51L1, rs6472903/8q21 and rs1550623/2q31/CDCA7). In conclusion, we have identified one novel lobular breast cancer specific predisposition polymorphism at 7q34, and shown for the first time that common breast cancer polymorphisms predispose to LCIS. We have shown that many of the ER+ breast cancer predisposition loci also predispose to ILC, although there is some heterogeneity between ER+ lobular and ER+ IDC tumors. These data provide evidence for overlapping, but distinct etiological pathways within ER+ breast cancer between morphological subtypes.