Comparison of the Agilent, ROMA/NimbleGen and Illumina platforms for classification of copy number alterations in human breast tumors

<p>Abstract</p> <p>Background</p> <p>Microarray Comparative Genomic Hybridization (array CGH) provides a means to examine DNA copy number aberrations. Various platforms, brands and underlying technologies are available, facing the user with many choices regarding platform sensitivity and number, localization, and density distribution of probes.</p> <p>Results</p> <p>We evaluate three different platforms presenting different nature and arrangement of the probes: The Agilent Human Genome CGH Microarray 44 k, the ROMA/NimbleGen Representational Oligonucleotide Microarray 82 k, and the Illumina Human-1 Genotyping 109 k BeadChip, with Agilent being gene oriented, ROMA/NimbleGen being genome oriented, and Illumina being genotyping oriented. We investigated copy number changes in 20 human breast tumor samples representing different gene expression subclasses, using a suite of graphical and statistical methods designed to work across platforms. Despite substantial differences in the composition and spatial distribution of probes, the comparison revealed high overall concordance. Notably however, some short amplifications and deletions of potential biological importance were not detected by all platforms. Both correlation and cluster analysis indicate a somewhat higher similarity between ROMA/NimbleGen and Illumina than between Agilent and the other two platforms. The programs developed for the analysis are available from <url>http://www.ifi.uio.no/bioinf/Projects/</url>.</p> <p>Conclusion</p> <p>We conclude that platforms based on different technology principles reveal similar aberration patterns, although we observed some unique amplification or deletion peaks at various locations, only detected by one of the platforms. The correct platform choice for a particular study is dependent on whether the appointed research intention is gene, genome, or genotype oriented.</p

Linear and non-linear dependencies between copy number aberrations and mRNA expression reveal distinct molecular pathways in breast cancer

<p>Abstract</p> <p>Background</p> <p>Elucidating the exact relationship between gene copy number and expression would enable identification of regulatory mechanisms of abnormal gene expression and biological pathways of regulation. Most current approaches either depend on linear correlation or on nonparametric tests of association that are insensitive to the exact shape of the relationship. Based on knowledge of enzyme kinetics and gene regulation, we would expect the functional shape of the relationship to be gene dependent and to be related to the gene regulatory mechanisms involved. Here, we propose a statistical approach to investigate and distinguish between linear and nonlinear dependences between DNA copy number alteration and mRNA expression.</p> <p>Results</p> <p>We applied the proposed method to DNA copy numbers derived from Illumina 109 K SNP-CGH arrays (using the log R values) and expression data from Agilent 44 K mRNA arrays, focusing on commonly aberrated genomic loci in a collection of 102 breast tumors. Regression analysis was used to identify the type of relationship (linear or nonlinear), and subsequent pathway analysis revealed that genes displaying a linear relationship were overall associated with substantially different biological processes than genes displaying a nonlinear relationship. In the group of genes with a linear relationship, we found significant association to canonical pathways, including purine and pyrimidine metabolism (for both deletions and amplifications) as well as estrogen metabolism (linear amplification) and BRCA-related response to damage (linear deletion). In the group of genes displaying a nonlinear relationship, the top canonical pathways were specific pathways like PTEN and PI13K/AKT (nonlinear amplification) and Wnt(B) and IL-2 signalling (nonlinear deletion). Both amplifications and deletions pointed to the same affected pathways and identified cancer as the top significant disease and cell cycle, cell signaling and cellular development as significant networks.</p> <p>Conclusions</p> <p>This paper presents a novel approach to assessing the validity of the dependence of expression data on copy number data, and this approach may help in identifying the drivers of carcinogenesis.</p

Allele-specific disparity in breast cancer

Background In a cancer cell the number of copies of a locus may vary due to amplification and deletion and these variations are denoted as copy number alterations (CNAs). We focus on the disparity of CNAs in tumour samples, which were compared to those in blood in order to identify the directional loss of heterozygosity. Methods We propose a numerical algorithm and apply it to data from the Illumina 109K-SNP array on 112 samples from breast cancer patients. B-allele frequency (BAF) and log R ratio (LRR) of Illumina were used to estimate Euclidian distances. For each locus, we compared genotypes in blood and tumour for subset of samples being heterozygous in blood. We identified loci showing preferential disparity from heterozygous toward either the A/B-allele homozygous (allelic disparity). The chi-squared and Cochran-Armitage trend tests were used to examine whether there is an association between high levels of disparity in single nucleotide polymorphisms (SNPs) and molecular, clinical and tumour-related parameters. To identify pathways and network functions over-represented within the resulting gene sets, we used Ingenuity Pathway Analysis (IPA). Results To identify loci with a high level of disparity, we selected SNPs 1) with a substantial degree of disparity and 2) with substantial frequency (at least 50% of the samples heterozygous for the respective locus). We report the overall difference in disparity in high-grade tumours compared to low-grade tumours (p-value < 0.001) and significant associations between disparity in multiple single loci and clinical parameters. The most significantly associated network functions within the genes represented in the loci of disparity were identified, including lipid metabolism, small-molecule biochemistry, and nervous system development and function. No evidence for over-representation of directional disparity in a list of stem cell genes was obtained, however genes appeared to be more often altered by deletion than by amplification. Conclusions Our data suggest that directional loss and amplification exist in breast cancer. These are highly associated with grade, which may indicate that they are enforced with increasing number of cell divisions. Whether there is selective pressure for some loci to be preferentially amplified or deleted remains to be confirmed

arrayMap: A Reference Resource for Genomic Copy Number Imbalances in Human Malignancies

Background: The delineation of genomic copy number abnormalities (CNAs) from cancer samples has been instrumental for identification of tumor suppressor genes and oncogenes and proven useful for clinical marker detection. An increasing number of projects have mapped CNAs using high-resolution microarray based techniques. So far, no single resource does provide a global collection of readily accessible oncoge- nomic array data. Methodology/Principal Findings: We here present arrayMap, a curated reference database and bioinformatics resource targeting copy number profiling data in human cancer. The arrayMap database provides a platform for meta-analysis and systems level data integration of high-resolution oncogenomic CNA data. To date, the resource incorporates more than 40,000 arrays in 224 cancer types extracted from several resources, including the NCBI's Gene Expression Omnibus (GEO), EBIs ArrayExpress (AE), The Cancer Genome Atlas (TCGA), publication supplements and direct submissions. For the majority of the included datasets, probe level and integrated visualization facilitate gene level and genome wide data re- view. Results from multi-case selections can be connected to downstream data analysis and visualization tools. Conclusions/Significance: To our knowledge, currently no data source provides an extensive collection of high resolution oncogenomic CNA data which readily could be used for genomic feature mining, across a representative range of cancer entities. arrayMap represents our effort for providing a long term platform for oncogenomic CNA data independent of specific platform considerations or specific project dependence. The online database can be accessed at http://www.arraymap.org.Comment: 17 pages, 5 inline figures, 3 tables, supplementary figures/tables split into 4 PDF files; manuscript submitted to PLoS ON

Phylogenetics and evolution of Su(var)3-9 SET genes in land plants: rapid diversification in structure and function

<p>Abstract</p> <p>Background</p> <p>Plants contain numerous <it>Su(var)3-9 </it>homologues (<it>SUVH</it>) and related (<it>SUVR</it>) genes, some of which await functional characterization. Although there have been studies on the evolution of plant <it>Su(var)3-9 SET </it>genes, a systematic evolutionary study including major land plant groups has not been reported. Large-scale phylogenetic and evolutionary analyses can help to elucidate the underlying molecular mechanisms and contribute to improve genome annotation.</p> <p>Results</p> <p>Putative orthologs of plant Su(var)3-9 SET protein sequences were retrieved from major representatives of land plants. A novel clustering that included most members analyzed, henceforth referred to as core <it>Su(var)3-9 </it>homologues and related (<it>cSUVHR</it>) gene clade, was identified as well as all orthologous groups previously identified. Our analysis showed that plant Su(var)3-9 SET proteins possessed a variety of domain organizations, and can be classified into five types and ten subtypes. Plant <it>Su(var)3-9 SET </it>genes also exhibit a wide range of gene structures among different paralogs within a family, even in the regions encoding conserved PreSET and SET domains. We also found that the majority of SUVH members were intronless and formed three subclades within the SUVH clade.</p> <p>Conclusions</p> <p>A detailed phylogenetic analysis of the plant <it>Su(var)3-9 SET g</it>enes was performed. A novel deep phylogenetic relationship including most plant <it>Su(var)3-9 SET </it>genes was identified. Additional domains such as SAR, ZnF_C2H2 and WIYLD were early integrated into primordial PreSET/SET/PostSET domain organization. At least three classes of gene structures had been formed before the divergence of <it>Physcomitrella patens </it>(moss) from other land plants. One or multiple retroposition events might have occurred among <it>SUVH </it>genes with the donor genes leading to the V-2 orthologous group. The structural differences among evolutionary groups of plant <it>Su(var)3-9 SET </it>genes with different functions were described, contributing to the design of further experimental studies.</p

Arabidopsis COMPASS-Like Complexes Mediate Histone H3 Lysine-4 Trimethylation to Control Floral Transition and Plant Development

Histone H3 lysine-4 (H3K4) methylation is associated with transcribed genes in eukaryotes. In Drosophila and mammals, both di- and tri-methylation of H3K4 are associated with gene activation. In contrast to animals, in Arabidopsis H3K4 trimethylation, but not mono- or di-methylation of H3K4, has been implicated in transcriptional activation. H3K4 methylation is catalyzed by the H3K4 methyltransferase complexes known as COMPASS or COMPASS-like in yeast and mammals. Here, we report that Arabidopsis homologs of the COMPASS and COMPASS-like complex core components known as Ash2, RbBP5, and WDR5 in humans form a nuclear subcomplex during vegetative and reproductive development, which can associate with multiple putative H3K4 methyltransferases. Loss of function of ARABIDOPSIS Ash2 RELATIVE (ASH2R) causes a great decrease in genome-wide H3K4 trimethylation, but not in di- or mono-methylation. Knockdown of ASH2R or the RbBP5 homolog suppresses the expression of a crucial Arabidopsis floral repressor, FLOWERING LOCUS C (FLC), and FLC homologs resulting in accelerated floral transition. ASH2R binds to the chromatin of FLC and FLC homologs in vivo and is required for H3K4 trimethylation, but not for H3K4 dimethylation in these loci; overexpression of ASH2R causes elevated H3K4 trimethylation, but not H3K4 dimethylation, in its target genes FLC and FLC homologs, resulting in activation of these gene expression and consequent late flowering. These results strongly suggest that H3K4 trimethylation in FLC and its homologs can activate their expression, providing concrete evidence that H3K4 trimethylation accumulation can activate eukaryotic gene expression. Furthermore, our findings suggest that there are multiple COMPASS-like complexes in Arabidopsis and that these complexes deposit trimethyl but not di- or mono-methyl H3K4 in target genes to promote their expression, providing a molecular explanation for the observed coupling of H3K4 trimethylation (but not H3K4 dimethylation) with active gene expression in Arabidopsis

Embryogenic potential and expression of embryogenesis-related genes in conifers are affected by treatment with a histone deacetylase inhibitor

Somatic embryogenesis is used for vegetative propagation of conifers. Embryogenic cultures can be established from zygotic embryos; however, the embryogenic potential decreases during germination. In Arabidopsis, LEAFY COTYLEDON (LEC) genes are expressed during the embryonic stage, and must be repressed to allow germination. Treatment with the histone deacetylase inhibitor trichostatin A (TSA) causes de-repression of LEC genes. ABSCISICACID3 (ABI3) and its Zeamays ortholog VIVIPAROUS1 (VP1) act together with the LEC genes to promote embryo maturation. In this study, we have asked the question whether TSA treatment in a conifer affects the embryogenic potential and the expression of embryogenesis-related genes. We isolated two conifer LEC1-type HAP3 genes, HAP3A and HAP3B, from Picea abies and Pinus sylvestris. A comparative phylogenetic analysis of plant HAP3 genes suggests that HAP3A and HAP3B are paralogous genes originating from a duplication event in the conifer lineage. The expression of HAP3A is high, in both somatic and zygotic embryos, during early embryo development, but decreases during late embryogeny. In contrast, the expression of VP1 is initially low but increases during late embryogeny. After exposure to TSA, germinating somatic embryos of P. abies maintain the competence to differentiate embryogenic tissue, and simultaneously the germination progression is partially inhibited. Furthermore, when embryogenic cultures of P. abies are exposed to TSA during embryo maturation, the maturation process is arrested and the expression levels of PaHAP3A and PaVP1 are maintained, suggesting a possible link between chromatin structure and expression of embryogenesis-related genes in conifers

The degree of segmental aneuploidy measured by total copy number abnormalities predicts survival and recurrence in superficial gastroesophageal adenocarcinoma

Background: Prognostic biomarkers are needed for superficial gastroesophageal adenocarcinoma (EAC) to predict clinical outcomes and select therapy. Although recurrent mutations have been characterized in EAC, little is known about their clinical and prognostic significance. Aneuploidy is predictive of clinical outcome in many malignancies but has not been evaluated in superficial EAC. Methods: We quantified copy number changes in 41 superficial EAC using Affymetrix SNP 6.0 arrays. We identified recurrent chromosomal gains and losses and calculated the total copy number abnormality (CNA) count for each tumor as a measure of aneuploidy. We correlated CNA count with overall survival and time to first recurrence in univariate and multivariate analyses. Results: Recurrent segmental gains and losses involved multiple genes, including: HER2, EGFR, MET, CDK6, KRAS (recurrent gains); and FHIT, WWOX, CDKN2A/B, SMAD4, RUNX1 (recurrent losses). There was a 40-fold variation in CNA count across all cases. Tumors with the lowest and highest quartile CNA count had significantly better overall survival (p = 0.032) and time to first recurrence (p = 0.010) compared to those with intermediate CNA counts. These associations persisted when controlling for other prognostic variables. Significance: SNP arrays facilitate the assessment of recurrent chromosomal gain and loss and allow high resolution, quantitative assessment of segmental aneuploidy (total CNA count). The non-monotonic association of segmental aneuploidy with survival has been described in other tumors. The degree of aneuploidy is a promising prognostic biomarker in a potentially curable form of EAC. © 2014 Davison et al

The ASH1 HOMOLOG 2 (ASHH2) Histone H3 Methyltransferase Is Required for Ovule and Anther Development in Arabidopsis

BACKGROUND:SET-domain proteins are histone lysine (K) methyltransferases (HMTase) implicated in defining transcriptionally permissive or repressive chromatin. The Arabidopsis ASH1 HOMOLOG 2 (ASHH2) protein (also called SDG8, EFS and CCR1) has been suggested to methylate H3K4 and/or H3K36 and is similar to Drosophila ASH1, a positive maintainer of gene expression, and yeast Set2, a H3K36 HMTase. Mutation of the ASHH2 gene has pleiotropic developmental effects. Here we focus on the role of ASHH2 in plant reproduction. METHODOLOGY/PRINCIPAL FINDINGS:A slightly reduced transmission of the ashh2 allele in reciprocal crosses implied involvement in gametogenesis or gamete function. However, the main requirement of ASHH2 is sporophytic. On the female side, close to 80% of mature ovules lack embryo sac. On the male side, anthers frequently develop without pollen sacs or with specific defects in the tapetum layer, resulting in reduction in the number of functional pollen per anther by up to approximately 90%. In consistence with the phenotypic findings, an ASHH2 promoter-reporter gene was expressed at the site of megaspore mother cell formation as well as tapetum layers and pollen. ashh2 mutations also result in homeotic changes in floral organ identity. Transcriptional profiling identified more than 300 up-regulated and 600 down-regulated genes in ashh2 mutant inflorescences, whereof the latter included genes involved in determination of floral organ identity, embryo sac and anther/pollen development. This was confirmed by real-time PCR. In the chromatin of such genes (AP1, AtDMC1 and MYB99) we observed a reduction of H3K36 trimethylation (me3), but not H3K4me3 or H3K36me2. CONCLUSIONS/SIGNIFICANCE:The severe distortion of reproductive organ development in ashh2 mutants, argues that ASHH2 is required for the correct expression of genes essential to reproductive development. The reduction in the ashh2 mutant of H3K36me3 on down-regulated genes relevant to the observed defects, implicates ASHH2 in regulation of gene expression via H3K36 trimethylation in chromatin of Arabidopsis inflorescences

Molecular Characterization of the Region 7q22.1 in Splenic Marginal Zone Lymphomas

Splenic marginal zone lymphomas (SMZL) are an uncommon type of B-cell non-Hodgkin's lymphoma (NHL-B) in which no specific chromosomal translocations have been described. In contrast, the most frequent cytogenetic abnormality is the loss of the long arm of chromosome 7 (7q). Previous reports have located this loss in the 7q32 region. In order to better characterize the genomic imbalances in SMZL, molecular studies were carried out in 73 patients with SMZL. To gain insight into the mapping at 7q a tiling array was also used. The results confirmed the loss of 7q as the most frequent change. In addition, several abnormalities, including 4q22.1, 1q21.3–q22, 6q25.3, 20q13.33, 3q28, 2q23.3–q24.1 and 17p13, were also present. A loss of 7q22.1 at 99925039–101348479 bp was observed in half of the cases. The region of 7q22.1 has not previously been characterised in SMZL. Our results confirmed the presence of a new region of loss on chromosome 7 in these NHL