A compendium of myeloma-associated chromosomal copy number abnormalities and their prognostic value.

To obtain a comprehensive genomic profile of presenting multiple myeloma cases we performed high-resolution single nucleotide polymorphism mapping array analysis in 114 samples alongside 258 samples analyzed by U133 Plus 2.0 expression array (Affymetrix). We examined DNA copy number alterations and loss of heterozygosity (LOH) to define the spectrum of minimally deleted regions in which relevant genes of interest can be found. The most frequent deletions are located at 1p (30%), 6q (33%), 8p (25%), 12p (15%), 13q (59%), 14q (39%), 16q (35%), 17p (7%), 20 (12%), and 22 (18%). In addition, copy number-neutral LOH, or uniparental disomy, was also prevalent on 1q (8%), 16q (9%), and X (20%), and was associated with regions of gain and loss. Based on fluorescence in situ hybridization and expression quartile analysis, genes of prognostic importance were found to be located at 1p (FAF1, CDKN2C), 1q (ANP32E), and 17p (TP53). In addition, we identified common homozygously deleted genes that have functions relevant to myeloma biology. Taken together, these analyses indicate that the crucial pathways in myeloma pathogenesis include the nuclear factor-κB pathway, apoptosis, cell-cycle regulation, Wnt signaling, and histone modifications. This study was registered at http://isrctn.org as ISRCTN68454111.

Homozygous deletion mapping in myeloma samples identifies genes and an expression signature relevant to pathogenesis and outcome.

PURPOSE: Myeloma is a clonal malignancy of plasma cells. Poor-prognosis risk is currently identified by clinical and cytogenetic features. However, these indicators do not capture all prognostic information. Gene expression analysis can be used to identify poor-prognosis patients and this can be improved by combination with information about DNA-level changes. EXPERIMENTAL DESIGN: Using single nucleotide polymorphism-based gene mapping in combination with global gene expression analysis, we have identified homozygous deletions in genes and networks that are relevant to myeloma pathogenesis and outcome. RESULTS: We identified 170 genes with homozygous deletions and corresponding loss of expression. Deletion within the "cell death" network was overrepresented and cases with these deletions had impaired overall survival. From further analysis of these events, we have generated an expression-based signature associated with shorter survival in 258 patients and confirmed this signature in data from two independent groups totaling 800 patients. We defined a gene expression signature of 97 cell death genes that reflects prognosis and confirmed this in two independent data sets. CONCLUSIONS: We developed a simple 6-gene expression signature from the 97-gene signature that can be used to identify poor-prognosis myeloma in the clinical environment. This signature could form the basis of future trials aimed at improving the outcome of poor-prognosis myeloma.

Deletions of CDKN2C in Multiple Myeloma: Biological and Clinical Implications

Purpose: Deletions of chromosome 1 have been described in 7% to 40% of cases of myeloma with inconsistent clinical consequences. CDKN2C at 1p32.3 has been identified in myeloma cell lines as the potential target of the deletion. We tested the clinical impact of 1p deletion and used high-resolution techniques to define the role of CDKN2C in primary patient material.Experimental Design: We analyzed 515 cases of monoclonal gammopathy of undetermined significance (MGUS), smoldering multiple myeloma (SMM), and newly diagnosed multiple myeloma using fluorescence in situ hybridization (FISH) for deletions of CDKN2C. In 78 myeloma cases, we carried out Affymetrix single nucleotide polymorphism mapping and U133 Plus 2.0 expression arrays. In addition, we did mutation, methylation, and Western blotting analysis.Results: By FISH we identified deletion of 1p32.3 (CDKN2C) in 3 of 66 MGUS (4.5%), 4 of 39 SMM (10.3%), and 55 of 369 multiple myeloma cases (15%). We examined the impact of copy number change at CDKN2C on overall survival (OS), and found that the cases with either hemizygous or homozygous deletion of CDKN2C had a worse OS compared with cases that were intact at this region (22 months versus 38 months; P = 0.003). Using gene mapping we identified three homozygous deletions at 1p32.3, containing CDKN2C, all of which lacked expression of CDKN2C. Cases with homozygous deletions of CDKN2C were the most proliferative myelomas, defined by an expression-based proliferation index, consistent with its biological function as a cyclin-dependent kinase inhibitor.Conclusions: Our results suggest that deletions of CDKN2C are important in the progression and clinical outcome of myeloma.

Gene mapping and expression analysis of 16q loss of heterozygosity identifies WWOX and CYLD as being important in determining clinical outcome in multiple myeloma.

We performed fluorescent in situ hybridization (FISH) for 16q23 abnormalities in 861 patients with newly diagnosed multiple myeloma and identified deletion of 16q [del(16q)] in 19.5%. In 467 cases in which demographic and survival data were available, del(16q) was associated with a worse overall survival (OS). It was an independent prognostic marker and conferred additional adverse survival impact in cases with the known poor-risk cytogenetic factors t(4;14) and del(17p). Gene expression profiling and gene mapping using 500K single-nucleotide polymorphism (SNP) mapping arrays revealed loss of heterozygosity (LOH) involving 3 regions: the whole of 16q, a region centered on 16q12 (the location of CYLD), and a region centered on 16q23 (the location of the WW domain-containing oxidoreductase gene WWOX). CYLD is a negative regulator of the NF-kappaB pathway, and cases with low expression of CYLD were used to define a "low-CYLD signature." Cases with 16q LOH or t(14;16) had significantly reduced WWOX expression. WWOX, the site of the translocation breakpoint in t(14;16) cases, is a known tumor suppressor gene involved in apoptosis, and we were able to generate a "low-WWOX signature" defined by WWOX expression. These 2 genes and their corresponding pathways provide an important insight into the potential mechanisms by which 16q LOH confers poor prognosis.

Integration of global SNP-based mapping and expression arrays reveals key regions, mechanisms, and genes important in the pathogenesis of multiple myeloma.

Multiple myeloma is characterized by genomic alterations frequently involving gains and losses of chromosomes. Single nucleotide polymorphism (SNP)-based mapping arrays allow the identification of copy number changes at the sub-megabase level and the identification of loss of heterozygosity (LOH) due to monosomy and uniparental disomy (UPD). We have found that SNP-based mapping array data and fluorescence in situ hybridization (FISH) copy number data correlated well, making the technique robust as a tool to investigate myeloma genomics. The most frequently identified alterations are located at 1p, 1q, 6q, 8p, 13, and 16q. LOH is found in these large regions and also in smaller regions throughout the genome with a median size of 1 Mb. We have identified that UPD is prevalent in myeloma and occurs through a number of mechanisms including mitotic nondisjunction and mitotic recombination. For the first time in myeloma, integration of mapping and expression data has allowed us to reduce the complexity of standard gene expression data and identify candidate genes important in both the transition from normal to monoclonal gammopathy of unknown significance (MGUS) to myeloma and in different subgroups within myeloma. We have documented these genes, providing a focus for further studies to identify and characterize those that are key in the pathogenesis of myeloma.

Genome-wide association study reveals greater polygenic loading for schizophrenia in cases with a family history of illness

Genome-wide association studies (GWAS) of schizophrenia have yielded more than 100 common susceptibility variants, and strongly support a substantial polygenic contribution of a large number of small allelic effects. It has been hypothesized that familial schizophrenia is largely a consequence of inherited rather than environmental factors. We investigated the extent to which familiality of schizophrenia is associated with enrichment for common risk variants detectable in a large GWAS. We analyzed single nucleotide polymorphism (SNP) data for cases reporting a family history of psychotic illness (N = 978), cases reporting no such family history (N = 4,503), and unscreened controls (N = 8,285) from the Psychiatric Genomics Consortium (PGC1) study of schizophrenia. We used a multinomial logistic regression approach with model-fitting to detect allelic effects specific to either family history subgroup. We also considered a polygenic model, in which we tested whether family history positive subjects carried more schizophrenia risk alleles than family history negative subjects, on average. Several individual SNPs attained suggestive but not genome-wide significant association with either family history subgroup. Comparison of genome-wide polygenic risk scores based on GWAS summary statistics indicated a significant enrichment for SNP effects among family history positive compared to family history negative cases (Nagelkerke's R(2 ) = 0.0021; P = 0.00331; P-value threshold <0.4). Estimates of variability in disease liability attributable to the aggregate effect of genome-wide SNPs were significantly greater for family history positive compared to family history negative cases (0.32 and 0.22, respectively; P = 0.031). We found suggestive evidence of allelic effects detectable in large GWAS of schizophrenia that might be specific to particular family history subgroups. However, consideration of a polygenic risk score indicated a significant enrichment among family history positive cases for common allelic effects. Familial illness might, therefore, represent a more heritable form of schizophrenia, as suggested by previous epidemiological studies.

Association of VEGF-1499C -> T polymorphism with diabetic nephropathy in type 1 diabetes mellitus

Vascular endothelial growth factor (VEGF) is reported to be implicated in the development of diabetic nephropathy. We performed a case-control study to determine if VEGF-2578C -> A, VEGF-1499C -> T, and VEGF-635G -> C single-nucleotide polymorphisms (SNPs) in the VEGF gene are associated with predisposition to diabetic nephropathy in type I diabetes. Genomic DNA was obtained from Irish type I diabetic individuals with nephropathy (cases, n=242) and those without nephropathy (controls, n=301), in addition to 400 healthy control samples. These samples were genotyped for the three SNPs using TaqMan or Pyrosequencing technology. Chi-squared analyses revealed no significant differences in genotype or allele frequencies in cases versus controls for VEGF-2578C -> A (genotype, P=.58; allele, P=.52) and VEGF-635G -> C (genotype, P=.58; allele, P=.33). However, a positive association with diabetic nephropathy was observed for the VEGF-1499T allele in the Northern Ireland population (P

Trait directed de novo population transcriptome dissects genetic regulation of a balanced polymorphism in phosphorus nutrition/arsenate tolerance in a wild grass Holcus lanatus

The aim of this study was to characterize the transcriptome of a balanced polymorphism, under the regulation of a single gene, for phosphate fertilizer responsiveness/arsenate toler- ance in wild grass Holcus lanatus genotypes screened from the same habitat.

De novo transcriptome sequencing, RNAseq (RNA sequencing) and single nucleotide poly- morphism (SNP) calling were conducted on RNA extracted from H.lanatus. Roche 454 sequencing data were assembled into c. 22 000 isotigs, and paired-end Illumina reads for phosphorus-starved (P) and phosphorus-treated (P+) genovars of tolerant (T) and nontoler- ant (N) phenotypes were mapped to this reference transcriptome.

Heatmaps of the gene expression data showed strong clustering of each P+/P treated genovar, as well as clustering by N/T phenotype. Statistical analysis identified 87 isotigs to be significantly differentially expressed between N and T phenotypes and 258 between P+ and P treated plants. SNPs and transcript expression that systematically differed between N and T phenotypes had regulatory function, namely proteases, kinases and ribonuclear RNA- binding protein and transposable elements.

A single gene for arsenate tolerance led to distinct phenotype transcriptomes and SNP pro- files, with large differences in upstream post-translational and post-transcriptional regulatory genes rather than in genes directly involved in P nutrition transport and metabolism per se.

The Cyclin D1 (CCND1) A870G polymorphism predicts clinical outcome to lapatinib and capecitabine in HER2-positive metastatic breast cancer

BACKGROUND: Lapatinib plus capecitabine emerged as an efficacious therapy in metastatic breast cancer (mBC). We aimed to identify germline single-nucleotide polymorphisms (SNPs) in genes involved in capecitabine catabolism and human epidermal receptor signaling that were associated with clinical outcome to assist in selecting patients likely to benefit from this combination.

PATIENTS AND METHODS: DNA was extracted from 240 of 399 patients enrolled in EGF100151 clinical trial (NCT00078572; clinicaltrials.gov) and SNPs were successfully evaluated in 234 patients. The associations between SNPs and clinical outcome were analyzed using Fisher's exact test, Kaplan-Meier curves, log-rank tests, likelihood ratio test within logistic or Cox regression model, as appropriate.

RESULTS: There were significant interactions between CCND1 A870G and clinical outcome. Patients carrying the A-allele were more likely to benefit from lapatinib plus capecitabine versus capecitabine when compared with patients harboring G/G (P = 0.022, 0.024 and 0.04, respectively). In patients with the A-allele, the response rate (RR) was significantly higher with lapatinib plus capecitabine (35%) compared with capecitabine (11%; P = 0.001) but not between treatments in patients with G/G (RR = 24% and 32%, respectively; P = 0.85). Time to tumor progression (TTP) was longer in patients with the A-allele treated with lapatinib plus capecitabine compared with capecitabine (median TTP = 7.9 and 3.4 months; P < 0.001), but not in patients with G/G (median TTP = 6.1 and 6.6 months; P = 0.92).

CONCLUSION: Our findings suggest that CCND1A870G may be useful in predicting clinical outcome in HER2-positive mBC patients treated with lapatinib plus capecitabine.

The cyclin D1 (CCND1) rs9344 G>A polymorphism predicts clinical outcome in colon cancer patients treated with adjuvant 5-FU-based chemotherapy

Recent evidence indicates a potential prognostic and predictive value for germline polymorphisms in genes involved in cell cycle control. We investigated the effect of cyclin D1 (CCND1) rs9344 G>A in stage II/III colon cancer patients and validated the findings in an independent study cohort. For evaluation and validation set, a total of 264 and 234 patients were included. Patients treated with 5-fluorouracil-based chemotherapy, carrying the CCND1 rs9344 A/A genotype had significantly decreased time-to-tumor recurrence (TTR) in univariate analysis and multivariate analysis (hazard ratio (HR) 2.47; 95% confidence interval (CI) 1.16-5.29; P=0.019). There was no significant association between CCND1 rs9344 G>A and TTR in patients with curative surgery alone. In the validation set, the A allele of CCND1 rs9344 G>A remained significantly associated with decreased TTR in univariate and multivariate analyses (HR 1.94; 95% CI 1.05-3.58; P=0.035). CCND1 rs9344 G>A may be a predictive and/or prognostic biomarker in stage II/III colon cancer patients, however, prospective trials are warranted to confirm our findings.

Exome sequencing of gastric adenocarcinoma identifies recurrent somatic mutations in cell adhesion and chromatin remodelling genes

Gastric cancer is a major cause of global cancer mortality. We surveyed the spectrum of somatic alterations in gastric cancer by sequencing the exomes of 15 gastric adenocarcinomas and their matched normal DNAs. Frequently mutated genes in the adenocarcinomas included TP53 (11/15 tumors), PIK3CA (3/15) and ARID1A (3/15). Cell adhesion was the most enriched biological pathway among the frequently mutated genes. A prevalence screening confirmed mutations in FAT4, a cadherin family gene, in 5% of gastric cancers (6/110) and FAT4 genomic deletions in 4% (3/83) of gastric tumors. Frequent mutations in chromatin remodeling genes (ARID1A, MLL3 and MLL) also occurred in 47% of the gastric cancers. We detected ARID1A mutations in 8% of tumors (9/110), which were associated with concurrent PIK3CA mutations and microsatellite instability. In functional assays, we observed both FAT4 and ARID1A to exert tumor-suppressor activity. Somatic inactivation of FAT4 and ARID1A may thus be key tumorigenic events in a subset of gastric cancers.

Conservation genetics of Ireland's sole population of the River water crowfoot (Ranunculus fluitans Lam.)

Populations of many freshwater species are becoming increasingly threatened as a result of a wide range of anthropogenically mediated factors. In the present study, we wanted to assess levels and patterns of genetic diversity in Ireland's sole population of the River water crowfoot (Ranunculus fluitans), which is restricted to a 12 km stretch of a single river, to assist the formation of conservation strategies. Analysis using amplified fragment length polymorphism (AFLP) indicated comparable levels of genetic diversity to those exhibited by a more extensive population of the species in England, and revealed no evidence of clonal reproduction. Allele-specific PCR analysis of five nuclear single nucleotide polymorphisms (SNPs) indicated no evidence of hybridization with its more abundant congener Ranunculus penicillatus, despite previous anecdotal reports of the occurrence of hybrids. Although the population currently exhibits healthy levels of genetic diversity and is not at risk of genetic assimilation via hybridization with R. penicillatus, it still remains vulnerable to other factors such as stochastic events and invasive species. © 2013 Elsevier B.V. All rights reserved.

Association Study of 167 Candidate Genes for Schizophrenia Selected by a Multi-Domain Evidence-Based Prioritization Algorithm and Neurodevelopmental Hypothesis

Integrating evidence from multiple domains is useful in prioritizing disease candidate genes for subsequent testing. We ranked all known human genes (n = 3819) under linkage peaks in the Irish Study of High-Density Schizophrenia Families using three different evidence domains: 1) a meta-analysis of microarray gene expression results using the Stanley Brain collection, 2) a schizophrenia protein-protein interaction network, and 3) a systematic literature search. Each gene was assigned a domain-specific p-value and ranked after evaluating the evidence within each domain. For comparison to this
ranking process, a large-scale candidate gene hypothesis was also tested by including genes with Gene Ontology terms related to neurodevelopment. Subsequently, genotypes of 3725 SNPs in 167 genes from a custom Illumina iSelect array were used to evaluate the top ranked vs. hypothesis selected genes. Seventy-three genes were both highly ranked and involved in neurodevelopment (category 1) while 42 and 52 genes were exclusive to neurodevelopment (category 2) or highly ranked (category 3), respectively. The most significant associations were observed in genes PRKG1, PRKCE, and CNTN4 but no individual SNPs were significant after correction for multiple testing. Comparison of the approaches showed an excess of significant tests using the hypothesis-driven neurodevelopment category. Random selection of similar sized genes from two independent genome-wide association studies (GWAS) of schizophrenia showed the excess was unlikely by chance. In a further meta-analysis of three GWAS datasets, four candidate SNPs reached nominal significance. Although gene ranking using integrated sources of prior information did not enrich for significant results in the current experiment, gene selection using an a priori hypothesis (neurodevelopment) was superior to random selection. As such, further development of gene ranking strategies using more carefully selected sources of information is warranted.

Genetic variation in the 6p22.3 gene DTNBP1, the human ortholog of the mouse dysbindin gene, is associated with schizophrenia

Prior evidence has supported the existence of multiple susceptibility genes for schizophrenia. Multipoint linkage analysis of the 270 Irish high-density pedigrees that we have studied, as well as results from several other samples, suggest that at least one such gene is located in region 6p24-21. In the present study, family-based association analysis of 36 simple sequence-length-polymorphism markers and of 17 SNP markers implicated two regions, separated by approximately 7 Mb. The first region, and the focus of this report, is 6p22.3. In this region, single-nucleotide polymorphisms within the 140-kb gene DTNBP1 (dystrobrevin-binding protein 1, or dysbindin) are strongly associated with schizophrenia. Uncorrected, empirical P values produced by the program TRANSMIT were significant (P

Identification of a high-risk haplotype for the dystrobrevin binding protein 1 (DTNBP1) gene in the Irish study of high-density schizophrenia families

A recent report showed significant associations between several SNPs in a previously unknown EST cluster with schizophrenia. (1). The cluster was identified as the human dystrobrevin binding protein 1 gene (DTNBP1) by sequence database comparisons and homology with mouse DTNBP1. (2). However, the linkage disequilibrium (LD) among the SNPs in DTNBP1 as well as the pattern of significant SNP-schizophrenia association was complex. This raised several questions such as the number of susceptibility alleles that may be involved and the size of the region where the actual disease mutation(s) could be located. To address these questions, we performed different single-marker tests on the 12 previously studied and 2 new SNPs in DTNBP1 that were re-scored using an improved procedure, and performed a variety of haplotype analyses. The sample consisted of 268 Irish multiplex families selected for high density of schizophrenia. Results suggested a simple structure where the LD in the target region could be explained by 6 haplotypes that together accounted for 96% of haplotype diversity in the whole sample. From these six, a single high-risk haplotype was identified that showed a significant association with schizophrenia and explained the pattern of significant findings in the analyses with individual markers. This haplotype was 30 kb long, had a large effect, could be measured with two tag SNPs only, had a frequency of 6% in our sample, seemed to be of relatively recent origin in evolutionary terms, and was equally distributed over Ireland. Implications of these findings for follow-up and replication studies are discussed.