SNP array profiling of childhood adrenocortical tumors reveals distinct pathways of tumorigenesis and highlights candidate driver genes

Childhood adrenocortical tumors (ACT) are rare malignancies, except in southern Brazil, where a higher incidence rate is associated to a high frequency of the founder R337H TP53 mutation. To date, copy number alterations in these tumors have only been analyzed by low-resolution comparative genomic hybridization.

Bayesian Hidden Markov Modeling of Array CGH Data

Genomic alterations have been linked to the development and progression of cancer. The technique of Comparative Genomic Hybridization (CGH) yields data consisting of fluorescence intensity ratios of test and reference DNA samples. The intensity ratios provide information about the number of copies in DNA. Practical issues such as the contamination of tumor cells in tissue specimens and normalization errors necessitate the use of statistics for learning about the genomic alterations from array-CGH data. As increasing amounts of array CGH data become available, there is a growing need for automated algorithms for characterizing genomic profiles. Specifically, there is a need for algorithms that can identify gains and losses in the number of copies based on statistical considerations, rather than merely detect trends in the data. We adopt a Bayesian approach, relying on the hidden Markov model to account for the inherent dependence in the intensity ratios. Posterior inferences are made about gains and losses in copy number. Localized amplifications (associated with oncogene mutations) and deletions (associated with mutations of tumor suppressors) are identified using posterior probabilities. Global trends such as extended regions of altered copy number are detected. Since the posterior distribution is analytically intractable, we implement a Metropolis-within-Gibbs algorithm for efficient simulation-based inference. Publicly available data on pancreatic adenocarcinoma, glioblastoma multiforme and breast cancer are analyzed, and comparisons are made with some widely-used algorithms to illustrate the reliability and success of the technique.

A de novo 1.1-1.6 Mb subtelomeric deletion of chromosome 20q13.33 in a patient with learning difficulties but without obvious dysmorphic features

We report on a de novo submicroscopic deletion of 20q13.33 identified by subtelomeric fluorescence in situ hybridization (FISH) in a 4-year-old girl with learning difficulties, hyperlaxity and strabismus, but without obvious dysmorphic features. Further investigations by array-based comparative genomic hybridization (array-CGH) and FISH analysis allowed us to delineate the smallest reported subterminal deletion of chromosome 20q, spanning a 1.1-1.6 Mb with a breakpoint localized between BAC RP5-887L7 and RP11-261N11. The genes CHRNA4 and KCNQ2 implicated in autosomal dominant epilepsy are included in the deletion interval. Subterminal 20q deletions as found in the present patient have, to our knowledge, only been reported in three patients. We review the clinical and behavioral phenotype of such "pure" subterminal 20q deletions.

Losses of chromosome arms 4q, 8p, 13q and gain of 8q are correlated with increasing chromosomal instability in hepatocellular carcinoma

OBJECTIVE: Chromosomal instability is a key feature in hepatocellular carcinoma (HCC). Array comparative genomic hybridization (aCGH) revealed recurring structural aberrations, whereas fluorescence in situ hybridization (FISH) indicated an increasing number of numerical aberrations in dedifferentiating HCC. Therefore, we examined whether there was a correlation between structural and numerical aberrations of chromosomal instability in HCC. METHODS AND RESULTS: 27 HCC (5 well, 10 moderately, 12 lower differentiated) already cytogenetically characterized by aCGH were analyzed. FISH analysis using probes for chromosomes 1, 3, 7, 8 and 17 revealed 1.46-4.24 signals/nucleus, which correlated with the histological grade (well vs. moderately,p < 0.0003; moderately vs. lower, p < 0.004). The number of chromosomes to each other was stable with exceptions only seen for chromosome 8. Loss of 4q and 13q, respectively, were correlated with the number of aberrations detected by aCGH (p < 0.001, p < 0.005; Mann-Whitney test). Loss of 4q and gain of 8q were correlated with an increasing number of numerical aberrations detected by FISH (p < 0.020, p < 0.031). Loss of 8p was correlated with the number of structural imbalances seen in aCGH (p < 0.048), but not with the number of numerical changes seen in FISH. CONCLUSION: We found that losses of 4q, 8p and 13q were closely correlated with an increasing number of aberrations detected by aCGH, whereas a loss of 4q and a gain of 8q were also observed in the context of polyploidization, the cytogenetic correlate of morphological dedifferentiation.

Loss of DAXX and ATRX are associated with chromosome instability and reduced survival of patients with pancreatic neuroendocrine tumors

BACKGROUND & AIMS Sporadic pancreatic neuroendocrine tumors (pNETs) are rare and genetically heterogeneous. Chromosome instability (CIN) has been detected in pNETs from patients with poor outcomes, but no specific genetic factors have been associated with CIN. Mutations in death domain-associated protein gene (DAXX) or ATR-X gene (ATRX) (which both encode proteins involved in chromatin remodeling) have been detected in 40% of pNETs, in association with activation of alternative lengthening of telomeres. We investigated whether loss of DAXX or ATRX, and consequent alternative lengthening of telomeres, are related to CIN in pNETs. We also assessed whether loss of DAXX or ATRX is associated with specific phenotypes of pNETs. METHODS We collected well-differentiated primary pNET samples from 142 patients at the University Hospital Zurich and from 101 patients at the University Hospital Bern (both located in Switzerland). Clinical follow-up data were obtained for 149 patients from general practitioners and tumor registries. The tumors were reclassified into 3 groups according to the 2010 World Health Organization classification. Samples were analyzed by immunohistochemistry and telomeric fluorescence in situ hybridization. We correlated loss of DAXX, or ATRX, expression, and activation of alternative lengthening of telomeres with data from comparative genomic hybridization array studies, as well as with clinical and pathological features of the tumors and relapse and survival data. RESULTS Loss of DAXX or ATRX protein and alternative lengthening of telomeres were associated with CIN in pNETs. Furthermore, loss of DAXX or ATRX correlated with tumor stage and metastasis, reduced time of relapse-free survival, and decreased time of tumor-associated survival. CONCLUSIONS Loss of DAXX or ATRX is associated with CIN in pNETs and shorter survival times of patients. These results support the hypothesis that DAXX- and ATRX-negative tumors are a more aggressive subtype of pNET, and could lead to identification of strategies to target CIN in pancreatic tumors.

Congenital aural atresia associated with agenesis of internal carotid artery in a girl with a FOXI3 deletion.

We report on the molecular characterization of a microdeletion of approximately 2.5 Mb at 2p11.2 in a female baby with left congenital aural atresia, microtia, and ipsilateral internal carotid artery agenesis. The deletion was characterized by fluorescence in situ hybridization, array comparative genomic hybridization, and whole genome re-sequencing. Among the genes present in the deleted region, we focused our attention on the FOXI3 gene. Foxi3 is a member of the Foxi class of Forkhead transcription factors. In mouse, chicken and zebrafish Foxi3 homologues are expressed in the ectoderm and endoderm giving rise to elements of the jaw as well as external, middle and inner ear. Homozygous Foxi3-/- mice have recently been generated and show a complete absence of the inner, middle, and external ears as well as severe defects in the jaw and palate. Recently, a 7-bp duplication within exon 1 of FOXI3 that produces a frameshift and a premature stop codon was found in hairless dogs. Mild malformations of the outer auditory canal (closed ear canal) and ear lobe have also been noted in a fraction of FOXI3 heterozygote Peruvian hairless dogs. Based on the phenotypes of Foxi3 mutant animals, we propose that FOXI3 may be responsible for the phenotypic features of our patient. Further characterization of the genomic region and the analysis of similar patients may help to demonstrate this point. © 2015 Wiley Periodicals, Inc.

The search for the primary tumor in metastasized gastroenteropancreatic neuroendocrine neoplasm.

Gastroenteropancreatic neuroendocrine tumors (NETs) often present as liver metastasis from a carcinoma of unknown primary. We recently showed that primary NETs from the pancreas, small intestine and stomach as well as their respective liver metastases differ from each other by the expression profile of the three genes CD302, PPWD1 and ABHB14B. The gene and protein expression of CD302, PPWD1, and ABHB14B was studied in abdominal NET metastases to identify the site of the respective primary tumors. Cryopreserved tissue from NET metastases collected in different institutions (group A: 29, group B: 50, group C: 132 specimens) were examined by comparative genomic hybridization (Agilent 105 K), gene expression analysis (Agilent 44 K) (groups A and B) and immunohistochemistry (group C). The data were blindly evaluated, i.e. without knowing the site of the primary. Gene expression analysis correctly revealed the primary in the ileum in 94 % of the cases of group A and in 58 % of group B. A pancreatic primary was predicted in 83 % (group A) and 20 % (group B), respectively. The combined sensitivity of group A and B was 75 % for ileal NETs and 38 % for pancreatic NETs. Immunohistochemical analysis of group C revealed an overall sensitivity of 80 %. Gene and protein expression analysis of CD302 and PPWD1 in NET metastases correctly identifies the primary in the pancreas or the ileum in 80 % of the cases, provided that the tissue is well preserved. Immunohistochemical profiling revealed CD302 as the best marker for ileal and PPWD1 for pancreatic detection.

Expression cloning of a candidate prostate cancer oncogene

Genetic analysis, both karyotyping and comparative genomic hybridization, of prostate cancer cell lines and specimens have revealed multiple areas of concordant increases in DNA content. An increase of DNA in specific regions of the genome in cancer is often associated with the amplification of oncogenes. Based on these observations we have hypothesized that oncogenes are involved in the initiation or progression of prostate cancer. An expression cloning approach was utilized to identify candidate oncogenes in prostate cancer. ^ A full-length, unidirectional cDNA expression library was constructed from DU145 prostate cancer cells. The cDNA library was screened using CP12, a rat prostate epithelial cell line. In soft agarose assays, CP12 (parental or vector transfected) do not form colonies. However, upon the introduction of a number of known oncogenes CP12 becomes anchorage independent in soft agarose. Based on this in-vitro phenotypic shift, a DU145 cDNA library was stably transfected into CP12, and selected for anchorage independence. Two hundred fifty nine anchorage independent clones were isolated. Some colonies contained more than one insert, bringing the candidate oncogene pool to approximately 400. Seven inserts were sequenced at random. Using the sequences obtained, GenBank was screened, and matches were found with p53, PARG1, a mitochondrial ATPase, RNF6, and three unknown genes that mapped to Unigene clusters. As the pool of cDNA inserts appeared promising, overexpressed genes were further selected. From 259 clones, 17 clones were overexpressed more than 6-fold in DU145 compared to Normal Prostate. From the 17 clones, 12 cDNA inserts were strongly expressed in DU145 and were isolated for sequencing. ^ Two of the sequences, 1G6 and 3E9, were identical. Expression of 1G6/2G9/3E9 was tested by RT-PCR. 1G6/2G9/3E9 was not expressed in normal prostate, but was expressed in all prostate cancer cell lines tested as well as six prostate cancer samples. When retransfected into CP12, 1G6/2G9/3E9 induced the formation of foci and anchorage independent colonies. Thus, functional and expression data suggest that 1G6/2G9/3E9 may be a prostate cancer oncogene. ^

Physical and functional mapping of a novel tumor suppressor locus for prostate cancer within chromosome 10p12.31-q11

Prostate cancer remains the second leading cause of male cancer deaths in the United States, yet the molecular mechanisms underlying this disease remain largely unknown. Cytogenetic and molecular analyses of prostate tumors suggest a consistent association with the loss of chromosome 10. Previously, we have defined a novel tumor suppressor locus PAC-1 within chromosome 10pter-q11. Introduction of the short arm of chromosome 10 into a prostatic adenocarcinoma cell line PC-3H resulted in dramatic tumor suppression and restoration of a programmed cell death pathway. Using a combined approach of comparative genomic hybridization and microsatellite analysis of PC-3H, I have identified a region of hemizygosity within 10p12-p15. This region has been shown to be involved in frequent loss of heterozygosity in gliomas and melanoma. To functionally dissect the region within chromosome 10p containing PAC-1, we developed a strategy of serial microcell fusion, a technique that allows the transfer of defined fragments of chromosome 10p into PC-3H. Serial microcell fusion was used to transfer defined 10p fragments into a mouse A9 fibrosarcoma cell line. Once characterized by FISH and microsatellite analyses, the 10p fragments were subsequently transferred into PC-3H to generate a panel of microcell hybrid clones containing overlapping deletions of chromosome 10p. In vivo and microsatellite analyses of these PC hybrids identified a small chromosome 10p fragment (an estimated 31 Mb in size inclusive of the centromere) that when transferred into the PC-3H background, resulted in significant tumor suppression and limited a region of functional tumor suppressor activity to chromosome 10p12.31-q11. This region coincides with a region of LOH demonstrated in prostate cancer. These studies demonstrate the utility of this approach as a powerful tool to limit regions of functional tumor suppressor activity. Furthermore, these data used in conjunction with data generated by the Human Genome Project lent a focused approach to identify candidate tumor suppressor genes involved in prostate cancer. ^

High-resolution microarray analysis of chromosome 20q in human colon cancer metastasis model systems

Amplification of human chromosome 20q DNA is the most frequently occurring chromosomal abnormality detected in sporadic colorectal carcinomas and shows significant correlation with liver metastases. Through comprehensive high-resolution microarray comparative genomic hybridization and microarray gene expression profiling, we have characterized chromosome 20q amplicon genes associated with human colorectal cancer metastasis in two in vitro metastasis model systems. The results revealed increasing complexity of the 20q genomic profile from the primary tumor-derived cell lines to the lymph node and liver metastasis derived cell lines. Expression analysis of chromosome 20q revealed a subset of over expressed genes residing within the regions of genomic copy number gain in all the tumor cell lines, suggesting these are Chromosome 20q copy number responsive genes. Bases on their preferential expression levels in the model system cell lines and known biological function, four of the over expressed genes mapping to the common intervals of genomic copy gain were considered the most promising candidate colorectal metastasis-associated genes. Validation of genomic copy number and expression array data was carried out on these genes, with one gene, DNMT3B, standing out as expressed at a relatively higher levels in the metastasis-derived cell lines compared with their primary-derived counterparts in both the models systems analyzed. The data provide evidence for the role of chromosome 20q genes with low copy gain and elevated expression in the clonal evolution of metastatic cells and suggests that such genes may serve as early biomarkers of metastatic potential. The data also support the utility of the combined microarray comparative genomic hybridization and expression array analysis for identifying copy number responsive genes in areas of low DNA copy gain in cancer cells. ^

Chromosomal aberrations in PARP−/− mice: Genome stabilization in immortalized cells by reintroduction of poly(ADP-ribose) polymerase cDNA

Depletion of poly(ADP-ribose) polymerase (PARP) increases the frequency of recombination, gene amplification, sister chromatid exchanges, and micronuclei formation in cells exposed to genotoxic agents, implicating PARP in the maintenance of genomic stability. Flow cytometric analysis now has revealed an unstable tetraploid population in immortalized fibroblasts derived from PARP−/− mice. Comparative genomic hybridization detected partial chromosomal gains in 4C5-ter, 5F-ter, and 14A1-C1 in PARP−/−mice and immortalized PARP−/−fibroblasts. Neither the chromosomal gains nor the tetraploid population were apparent in PARP−/− cells stably transfected with PARP cDNA [PARP−/−(+PARP)], indicating negative selection of cells with these genetic aberrations after reintroduction of PARP cDNA. Although the tumor suppressor p53 was not detectable in PARP−/− cells, p53 expression was partially restored in PARP−/− (+PARP) cells. Loss of 14D3-ter that encompasses the tumor suppressor gene Rb-1 in PARP−/− mice was associated with a reduction in retinoblastoma(Rb) expression; increased expression of the oncogene Jun was correlated with a gain in 4C5-ter that harbors this oncogene. These results further implicate PARP in the maintenance of genomic stability and suggest that altered expression of p53, Rb, and Jun, as well as undoubtedly many other proteins may be a result of genomic instability associated with PARP deficiency.

Spectral karyotyping suggests additional subsets of colorectal cancers characterized by pattern of chromosome rearrangement

The abundant chromosome abnormalities in most carcinomas are probably a reflection of genomic instability present in the tumor, so the pattern and variability of chromosome abnormalities will reflect the mechanism of instability combined with the effects of selection. Chromosome rearrangement was investigated in 17 colorectal carcinoma-derived cell lines. Comparative genomic hybridization showed that the chromosome changes were representative of those found in primary tumors. Spectral karyotyping (SKY) showed that translocations were very varied and mostly unbalanced, with no translocation occurring in more than three lines. At least three karyotype patterns could be distinguished. Some lines had few chromosome abnormalities: they all showed microsatellite instability, the replication error (RER)+ phenotype. Most lines had many chromosome abnormalities: at least seven showed a surprisingly consistent pattern, characterized by multiple unbalanced translocations and intermetaphase variation, with chromosome numbers around triploid, 6–16 structural aberrations, and similarities in gains and losses. Almost all of these were RER−, but one, LS411, was RER+. The line HCA7 showed a novel pattern, suggesting a third kind of genomic instability: multiple reciprocal translocations, with little numerical change or variability. This line was also RER+. The coexistence in one tumor of two kinds of genomic instability is to be expected if the underlying defects are selected for in tumor evolution.

Quantitative analysis of chromosomal CGH in human breast tumors associates copy number abnormalities with p53 status and patient survival

We present a general method for rigorously identifying correlations between variations in large-scale molecular profiles and outcomes and apply it to chromosomal comparative genomic hybridization data from a set of 52 breast tumors. We identify two loci where copy number abnormalities are correlated with poor survival outcome (gain at 8q24 and loss at 9q13). We also identify a relationship between abnormalities at two loci and the mutational status of p53. Gain at 8q24 and loss at 5q15-5q21 are linked with mutant p53. The 9q and 5q losses suggest the possibility of gene products involved in breast cancer progression. The analytical techniques are general and also are applicable to the analysis of array-based expression data.

Gain of chromosome 3q defines the transition from severe dysplasia to invasive carcinoma of the uterine cervix.

We have chosen tumors of the uterine cervix as a model system to identify chromosomal aberrations that occur during carcinogenesis. A phenotype/genotype correlation was established in defined regions of archived, formalin-fixed, and hematoxylin/eosin-stained tissue sections that were dissected from normal cervical epithelium (n = 3), from mild (n = 4), moderate (n = 6), and severe dysplasias/carcinomas in situ (CIS) (n = 13), and from invasive carcinomas (n = 10) and investigated by comparative genomic hybridization. The same tissues were analyzed for DNA ploidy, proliferative activity, and the presence of human papillomavirus (HPV) sequences. The results show that an increase in proliferative activity and tetraploidization had occurred already in mildly dysplastic lesions. No recurrent chromosomal aberrations were observed in DNA extracted from normal epithelium or from mild and moderate dysplasias, indicating that the tetraploidization precedes the loss or gain of specific chromosomes. A gain of chromosome 3q became visible in one of the severe dysplasias/CIS. Notably, chromosome 3q was overrepresented in 90% of the carcinomas and was also found to have undergone a high-level copy-number increase (amplification). We therefore conclude that the gain of chromosome 3q that occurs in HPV16-infected, aneuploid cells represents a pivotal genetic aberration at the transition from severe dysplasia/CIS to invasive cervical carcinoma.

Caracterização de rearranjos cromossômicos citogeneticamente equilibrados associados a quadros clínicos

Este estudo teve como objetivos (a) identificar mecanismos pelos quais rearranjos cromossômicos citogeneticamente equilibrados possam estar associados de maneira causal a determinados quadros clínicos e (b) contribuir para a compreensão dos mecanismos de formação desses rearranjos. Para isso, foram estudados 45 rearranjos cromossômicos citogeneticamente equilibrados (29 translocações, 10 inversões e seis rearranjos complexos), detectados em pacientes que apresentavam malformações congênitas, comprometimento do desenvolvimento neuropsicomotor ou déficit intelectual. Foram 31 rearranjos cromossômicos esporádicos, três familiais que segregavam com o quadro clínico e mais 11 rearranjos cromossômicos herdados de genitores fenotipicamente normais. Inicialmente os pontos de quebra desses rearranjos foram mapeados por hibridação in situ fluorescente (FISH). A busca por microdeleções e duplicações genômicas foi realizada por a-CGH. A investigação dos pontos de quebra prosseguiu com a aplicação da técnica de Mate-Pair Sequencing (MPS), que permite localizar as quebras em segmentos de 100 pb - 1 kb, na maioria dos casos. Para obter os segmentos de junção das quebras no nível de pares de bases, os segmentos delimitados por MPS foram sequenciados pelo método de Sanger. A análise por aCGH revelou microdeleções ou microduplicações localizadas nos cromossomos rearranjados, em 12 dos 45 pacientes investigados (27%). A análise de 27 rearranjos por MPS permitiu a caracterização dos pontos de junção das quebras. MPS expandiu o número de pontos de quebra, detectados por análise do cariótipo ou aCGH, de 114 para 156 (em resolução < 2kb, na maioria dos casos). O número de pontos de quebra/rearranjo variou de 2 a 20. Os 156 pontos de quebra resultaram em 86 variantes estruturais equilibradas e outras 32 variantes não equilibradas. Perdas e ganhos de segmentos submiscroscópicos nos cromossomos rearranjados constituíram a principal causa ou, provavelmente, contribuíram para o quadro clínico de 12 dos 45 pacientes. Em cinco desses 12 rearranjos foram detectadas por MPS a interrupção de genes já relacionados à doença, ou provável alteração de sua região reguladora, contribundo para o quadro clínico. Em quatro dos 33 rearranjos não associados a perdas ou ganhos de segmentos, a análise por MPS revelou a interrupção de genes que já foram anteriormente relacionados a doenças, explicando-se, assim, as características clínicas dos portadores; outro rearranjo pode ter levando alteração da expressão gênica de gene sensível a dosagem e ao quadro clínico. Um rearranjo cromossômico familial, identificado na análise após bandamento G como uma translocação equilibrada, t(2;22)(p14;q12), segregava com quadro de atraso do desenvolvimento neuropsicomotor e dificuldade de aprendizado associados a dismorfismos. A combinação das análises por FISH, aCGH e MPS revelou que se tratava, na verdade, de rearranjo complexo entre os cromossomos 2, 5 e 22, incluindo 10 quebras. A segregação de diferentes desequilíbrios submicroscópicos em indivíduos afetados e clinicamente normais permitiu a compreensão da variabilidade clínica observada na família. Rearranjos equilibrados detectados em indivíduos afetados, mas herdados de genitores clinicamente normais, são, em geral, considerados como não tendo relação com o quadro clínico, apesar da possibilidade de desequilíbrios cromossômicos gerados por permuta desigual na meiose do genitor portador do rearranjo. Neste trabalho, a investigação de 11 desses rearranjos por aCGH não revelou perdas ou ganhos de segmentos nos cromossomos rearranjados. No entanto, a análise por aCGH da portadora de um desses rearranjos - inv(12)mat - revelou deleção de 8,7 Mb no cromossomo 8, como causa de seu fenótipo clínico. Essa deleção estava relacionada com outro rearranjo equilibrado também presente em sua mãe, independente da inversão. Para compreender os mecanismos de formação de rearranjos citogeneticamente equilibrados, investigamos os segmentos de junção no nível de pares de base. A análise por MPS que levou, na maioria dos casos, ao mapeamento dos pontos de quebras em segmentos <1kb permitiu o sequenciamento pelo método de Sanger de 51 segmentos de junções de 17 rearranjos. A ocorrência de blunt fusions ou inserções e deleções <10 pb, e a ausência de homologia ou a presença de micro homologia de 2 pb a 4 pb de extensão indicaram o mecanismo de junção de extremidades não homólogas (non-homologous end joinging; NHEJ), na maioria das 51 junções caracterizadas. As características de três dos quatro rearranjos mais complexos, com 17-20 quebras, indicaram sua formação pelo mecanismo de chromothripsis. Este estudo mostra a importância da análise genômica de variações de número de cópias por microarray, juntamente com o mapeamento dos pontos de quebra por MPS, para determinar a estrutura de rearranjos cromossômicos citogeneticamente equilibrados e seu impacto clínico. O mapeamento dos segmentos de junção por MPS, permitindo o sequenciamento pelo método de Sanger, foi essencial para a compreensão de mecanismos de formação desses rearranjos