Integrated Analysis of High-Resolution Gene Expression and Copy Number Profiling Identified Biallelic Deletion of CDKN2A/2B Tumor Suppressor Locus As the Most Frequent and Unique Genomic Abnormality in Diffuse Large B-Cell Lymphoma (DLBCL) with Strong Prognostic Value in Both GCB and ABC Subtypes and Not Overcome by a Dose-Intensive Immunochemotherapy Regimen Plus Rituximab. Results of a Prospective GELA Clinical Trial Program

Background and aim of the study: Genomic gains and losses play a crucial role in the development and progression of DLBCL and are closely related to gene expression profiles (GEP), including the germinal center B-cell like (GCB) and activated B-cell like (ABC) cell of origin (COO) molecular signatures. To identify new oncogenes or tumor suppressor genes (TSG) involved in DLBCL pathogenesis and to determine their prognostic values, an integrated analysis of high-resolution gene expression and copy number profiling was performed. Patients and methods: Two hundred and eight adult patients with de novo CD20+ DLBCL enrolled in the prospective multicentric randomized LNH-03 GELA trials (LNH03-1B, -2B, -3B, 39B, -5B, -6B, -7B) with available frozen tumour samples, centralized reviewing and adequate DNA/RNA quality were selected. 116 patients were treated by Rituximab(R)-CHOP/R-miniCHOP and 92 patients were treated by the high dose (R)-ACVBP regimen dedicated to patients younger than 60 years (y) in frontline. Tumour samples were simultaneously analysed by high resolution comparative genomic hybridization (CGH, Agilent, 144K) and gene expression arrays (Affymetrix, U133+2). Minimal common regions (MCR), as defined by segments that affect the same chromosomal region in different cases, were delineated. Gene expression and MCR data sets were merged using Gene expression and dosage integrator algorithm (GEDI, Lenz et al. PNAS 2008) to identify new potential driver genes. Results: A total of 1363 recurrent (defined by a penetrance > 5%) MCRs within the DLBCL data set, ranging in size from 386 bp, affecting a single gene, to more than 24 Mb were identified by CGH. Of these MCRs, 756 (55%) showed a significant association with gene expression: 396 (59%) gains, 354 (52%) single-copy deletions, and 6 (67%) homozygous deletions. By this integrated approach, in addition to previously reported genes (CDKN2A/2B, PTEN, DLEU2, TNFAIP3, B2M, CD58, TNFRSF14, FOXP1, REL...), several genes targeted by gene copy abnormalities with a dosage effect and potential physiopathological impact were identified, including genes with TSG activity involved in cell cycle (HACE1, CDKN2C) immune response (CD68, CD177, CD70, TNFSF9, IRAK2), DNA integrity (XRCC2, BRCA1, NCOR1, NF1, FHIT) or oncogenic functions (CD79b, PTPRT, MALT1, AUTS2, MCL1, PTTG1...) with distinct distribution according to COO signature. The CDKN2A/2B tumor suppressor locus (9p21) was deleted homozygously in 27% of cases and hemizygously in 9% of cases. Biallelic loss was observed in 49% of ABC DLBCL and in 10% of GCB DLBCL. This deletion was strongly correlated to age and associated to a limited number of additional genetic abnormalities including trisomy 3, 18 and short gains/losses of Chr. 1, 2, 19 regions (FDR < 0.01), allowing to identify genes that may have synergistic effects with CDKN2A/2B inactivation. With a median follow-up of 42.9 months, only CDKN2A/2B biallelic deletion strongly correlates (FDR p.value < 0.01) to a poor outcome in the entire cohort (4y PFS = 44% [32-61] respectively vs. 74% [66-82] for patients in germline configuration; 4y OS = 53% [39-72] vs 83% [76-90]). In a Cox proportional hazard prediction of the PFS, CDKN2A/2B deletion remains predictive (HR = 1.9 [1.1-3.2], p = 0.02) when combined with IPI (HR = 2.4 [1.4-4.1], p = 0.001) and GCB status (HR = 1.3 [0.8-2.3], p = 0.31). This difference remains predictive in the subgroup of patients treated by R-CHOP (4y PFS = 43% [29-63] vs. 66% [55-78], p=0.02), in patients treated by R-ACVBP (4y PFS = 49% [28-84] vs. 83% [74-92], p=0.003), and in GCB (4y PFS = 50% [27-93] vs. 81% [73-90], p=0.02), or ABC/unclassified (5y PFS = 42% [28-61] vs. 67% [55-82] p = 0.009) molecular subtypes (Figure 1). Conclusion: We report for the first time an integrated genetic analysis of a large cohort of DLBCL patients included in a prospective multicentric clinical trial program allowing identifying new potential driver genes with pathogenic impact. However CDKN2A/2B deletion constitutes the strongest and unique prognostic factor of chemoresistance to R-CHOP, regardless the COO signature, which is not overcome by a more intensified immunochemotherapy. Patients displaying this frequent genomic abnormality warrant new and dedicated therapeutic approaches.

Whose DNA is this? How relevant a question? (a note for forensic scientists).

This communication seeks to draw the attention of researchers and practitioners dealing with forensic DNA profiling analyses to the following question: is a scientist's report, offering support to a hypothesis according to which a particular individual is the source of DNA detected during the analysis of a stain, relevant from the point of view of a Court of Justice? This question relates to skeptical views previously voiced by commentators mainly in the judicial area, but is avoided by a large majority of forensic scientists. Notwithstanding, the pivotal role of this question has recently been evoked during the international conference "The hidden side of DNA profiles. Artifacts, errors and uncertain evidence" held in Rome (April 27th to 28th, 2012). Indeed, despite the fact that this conference brought together some of the world's leading forensic DNA specialists, it appeared clearly that a huge gap still exists between questions lawyers are actually interested in, and the answers that scientists deliver to Courts in written reports or during oral testimony. Participants in the justice system, namely lawyers and jurors on the one hand and forensic geneticists on the other, unfortunately talk considerably different languages. It thus is fundamental to address this issue of communication about results of forensic DNA analyses, and open a dialogue with practicing non-scientists at large who need to make meaningful use of scientific results to approach and help solve judicial cases. This paper intends to emphasize the actuality of this topic and suggest beneficial ways ahead towards a more reasoned use of forensic DNA in criminal proceedings.

Fission yeast rad51 and dmc1, two efficient DNA recombinases forming helical nucleoprotein filaments.

Homologous recombination is important for the repair of double-strand breaks during meiosis. Eukaryotic cells require two homologs of Escherichia coli RecA protein, Rad51 and Dmc1, for meiotic recombination. To date, it is not clear, at the biochemical level, why two homologs of RecA are necessary during meiosis. To gain insight into this, we purified Schizosaccharomyces pombe Rad51 and Dmc1 to homogeneity. Purified Rad51 and Dmc1 form homo-oligomers, bind single-stranded DNA preferentially, and exhibit DNA-stimulated ATPase activity. Both Rad51 and Dmc1 promote the renaturation of complementary single-stranded DNA. Importantly, Rad51 and Dmc1 proteins catalyze ATP-dependent strand exchange reactions with homologous duplex DNA. Electron microscopy reveals that both S. pombe Rad51 and Dmc1 form nucleoprotein filaments. Rad51 formed helical nucleoprotein filaments on single-stranded DNA, whereas Dmc1 was found in two forms, as helical filaments and also as stacked rings. These results demonstrate that Rad51 and Dmc1 are both efficient recombinases in lower eukaryotes and reveal closer functional and structural similarities between the meiotic recombinase Dmc1 and Rad51. The DNA strand exchange activity of both Rad51 and Dmc1 is most likely critical for proper meiotic DNA double-strand break repair in lower eukaryotes.

Meta-analysis of estrogen response in MCF-7 distinguishes early target genes involved in signaling and cell proliferation from later target genes involved in cell cycle and DNA repair.

ABSTRACT: BACKGROUND: Many studies have been published outlining the global effects of 17 beta-estradiol (E2) on gene expression in human epithelial breast cancer derived MCF-7 cells. These studies show large variation in results, reporting between ~100 and ~1500 genes regulated by E2, with poor overlap. RESULTS: We performed a meta-analysis of these expression studies, using the Rank product method to obtain a more accurate and stable list of the differentially expressed genes, and of pathways regulated by E2. We analyzed 9 time-series data sets, concentrating on response at 3-4 hrs (early) and at 24 hrs (late). We found >1000 statistically significant probe sets after correction for multiple testing at 3-4 hrs, and >2000 significant probe sets at 24 hrs. Differentially expressed genes were examined by pathway analysis. This revealed 15 early response pathways, mostly related to cell signaling and proliferation, and 20 late response pathways, mostly related to breast cancer, cell division, DNA repair and recombination. CONCLUSIONS: Our results show that meta-analysis identified more differentially expressed genes than the individual studies, and that these genes act together in networks. These results provide new insight into E2 regulated mechanisms, especially in the context of breast cancer.

Using graphical probability analysis (Bayes Nets) to evaluate a conditional DNA inclusion

This paper discusses the analysis of cases in which the inclusion or exclusion of a particular suspect, as a possible contributor to a DNA mixture, depends on the value of a variable (the number of contributors) that cannot be determined with certainty. It offers alternative ways to deal with such cases, including sensitivity analysis and object-oriented Bayesian networks, that separate uncertainty about the inclusion of the suspect from uncertainty about other variables. The paper presents a case study in which the value of DNA evidence varies radically depending on the number of contributors to a DNA mixture: if there are two contributors, the suspect is excluded; if there are three or more, the suspect is included; but the number of contributors cannot be determined with certainty. It shows how an object-oriented Bayesian network can accommodate and integrate varying perspectives on the unknown variable and how it can reduce the potential for bias by directing attention to relevant considerations and distinguishing different sources of uncertainty. It also discusses the challenge of presenting such evidence to lay audiences.

Genes encoding tumor necrosis factor alpha and granzyme A are expressed during development of autoimmune diabetes.

Progressive destruction of the insulin-producing beta cells in nonobese diabetic mice is observed after infiltration of the pancreas with lymphocytes [Makino, S., Kunimoto, K., Muraoka, Y., Mizushima, Y., Katagiri, K. &amp; Tochino, Y. (1980) Exp. Anim. (Tokyo) 29, 1-13]. We show that the genes for tumor necrosis factor alpha and granzyme A, a serine protease associated with cytoplasmic granules of cytotoxic cells, are expressed during the development of spontaneous diabetes mellitus in the nonobese diabetic mouse. Granzyme A-positive cells are found both in and surrounding the islets, implying induction prior to islet infiltration. Tumor necrosis factor alpha expression is exclusively observed in the intra-islet infiltrate, predominantly in lymphocytes adjacent to insulin-producing beta cells, the targets of the autoimmune destruction, implying that tumor necrosis factor alpha expression is induced locally--i.e., in the islet. A considerable portion of cells expressing tumor necrosis factor alpha appear to be CD4+ T cells. This T-cell subset was previously shown to be necessary for development of the disease. Thus, these findings may be important for understanding the pathogenesis of autoimmune diabetes mellitus and potentially also for that of other T-cell-mediated autoimmune diseases.

Construction and electrophoretic migration of single-stranded DNA knots and catenanes.

In recent years there has been growing interest in the question of how the particular topology of polymeric chains affects their overall dimensions and physical behavior. The majority of relevant studies are based on numerical simulation methods or analytical treatment; however, both these approaches depend on various assumptions and simplifications. Experimental verification is clearly needed but was hampered by practical difficulties in obtaining preparative amounts of knotted or catenated polymers with predefined topology and precisely set chain length. We introduce here an efficient method of production of various single-stranded DNA knots and catenanes that have the same global chain length. We also characterize electrophoretic migration of the produced single-stranded DNA knots and catenanes with increasing complexity.

Phosphorylation modulates FOXC2 transcriptional program by controlling the high-order interaction with DNA.

Phosphorylation of transcription factors is a rapid and reversible process linking cell signaling and control of gene expression, therefore understanding how it controls the transcription factor functions is one of the challenges of functional genomics. We performed such analysis for the forkhead transcription factor FOXC2 mutated in human hereditary disease lymphedemadistichiasis and important for the development of venous and lymphatic valves and lymphatic collecting vessels. We found that FOXC2 is phosphorylated in a cell-cycle dependent manner on eight evolutionary conserved serine/threonine residues, seven of which are clustered within a 70 amino acid domain. Surprisingly, the mutation of phosphorylation sites or a complete deletion of the domain did not affect the transcriptional activity of FOXC2 in a synthetic reporter assay. However, overexpression of the wild type or phosphorylation-deficient mutant resulted in overlapping but distinct gene expression profiles suggesting that binding of FOXC2 to individual sites under physiological conditions is affected by phosphorylation. To gain a direct insight into the role of FOXC2 phosphorylation, we performed comparative genome-wide location analysis (ChIP-chip) of wild type and phosphorylation-deficient FOXC2 in primary lymphatic endothelial cells. The effect of loss of phosphorylation on FOXC2 binding to genomic sites ranged from no effect to nearly complete inhibition of binding, suggesting a mechanism for how FOXC2 transcriptional program can be differentially regulated depending on FOXC2 phosphorylation status. Based on these results, we propose an extension to the enhanceosome model, where a network of genomic context-dependent DNA-protein and protein-protein interactions not only distinguishes a functional site from a nonphysiological site, but also determines whether binding to the functional site can be regulated by phosphorylation. Moreover, our results indicate that FOXC2 may have different roles in quiescent versus proliferating lymphatic endothelial cells in vivo.

Evaluation of heterogeneity of DNA ploidy in early gastric cancers.

DNA ploidy has been shown to be a predictive parameter for prognosis in various solid tumours. The prognostic value of DNA-ploidy in gastric cancers is still a matter of controversy. A possible explanation for the discrepant results reported in the literature could be sampling error in tumours with multiple stemlines differing in DNA-ploidy. In order to determine whether or not such heterogeneity exists in early gastric carcinoma, we have performed DNA cytophotometry on multiple samples of a group of 17 early gastric carcinomas, of which 8 were pure intramucosal and 9 were infiltrating into the submucosa. We found an aneuploid DNA-stemline in 8 (47%) early gastric cancers, more often in tumours invading into the submucosa (5/9) than in purely mucosal tumours (3/8). Multiple DNA-stemlines were found more frequently in submucosally infiltrating tumours (4/5). These results confirm the presence of DNA-aneuploid early gastric carcinoma which are frequently heterogeneous and suggest that heterogeneity occurs more frequently in tumours invading the submucosa. This heterogeneity is best detected by analysing multiple samples of tumours for DNA-ploidy.

Typing Candida albicans oral isolates from human immunodeficiency virus-infected patients by multilocus enzyme electrophoresis and DNA fingerprinting.

A total of 189 Candida albicans isolates have been typed by multilocus enzyme electrophoresis. The results obtained confirm the clonal mode of reproduction of C. albicans. The C. albicans populations found in the oropharynx of human immunodeficiency virus (HIV)-infected patients, in the oropharynx of healthy carriers, or in association with invasive candidiasis could not be distinguished. No clone or group of clones could be associated with the appearance of clinical disorders or with a reduced in vitro susceptibility to the antifungal agent fluconazole. Multiple and sequential oral isolates from 24 HIV-infected patients were also typed by restriction enzyme analysis with the enzymes EcoRI and HinfI and by use of the Ca3 repetitive probe. The results obtained by the combination of all three typing methods show that all but one patient each carried a unique major C. albicans clone in their oropharynx. The 21 patients with sequential isolates had the same C. albicans clones in their throats during recurrent oropharyngeal candidiasis episodes, independently of clinical status or of changes of in vitro susceptibility to fluconazole. Finally, several isolates of the same C. albicans clone found simultaneously in the oropharynx of a patient may present different levels of susceptibility to fluconazole.

Transcription factor CTF1 acts as a chromatin domain boundary that shields human telomeric genes from silencing.

Telomeres are associated with chromatin-mediated silencing of genes in their vicinity. However, how epigenetic markers mediate mammalian telomeric silencing and whether specific proteins may counteract this effect are not known. We evaluated the ability of CTF1, a DNA- and histone-binding transcription factor, to prevent transgene silencing at human telomeres. CTF1 was found to protect a gene from silencing when its DNA-binding sites were interposed between the gene and the telomeric extremity, while it did not affect a gene adjacent to the telomere. Protein fusions containing the CTF1 histone-binding domain displayed similar activities, while mutants impaired in their ability to interact with the histone did not. Chromatin immunoprecipitation indicated the propagation of a hypoacetylated histone structure to various extents depending on the telomere. The CTF1 fusion protein was found to recruit the H2A.Z histone variant at the telomeric locus and to restore high histone acetylation levels to the insulated telomeric transgene. Histone lysine trimethylations were also increased on the insulated transgene, indicating that these modifications may mediate expression rather than silencing at human telomeres. Overall, these results indicate that transcription factors can act to delimit chromatin domain boundaries at mammalian telomeres, thereby blocking the propagation of a silent chromatin structure.

Identification of HIV integration sites in infected host genomic DNA.

The integration of the Human Immunodeficiency Virus (HIV) genetic information into the host genome is fundamental for its replication and long-term persistence in the host. Isolating and characterizing the integration sites can be useful for obtaining data such as identifying the specific genomic location of integration or understanding the forces dictating HIV integration site selection. The methods outlined in this article describe a highly efficient and precise technique for identifying HIV integration sites in the host genome on a small scale using molecular cloning techniques and standard sequencing or on a massive scale using 454 pyrosequencing.

DNA methylation profiles of human active and inactive X chromosomes.

X-chromosome inactivation (XCI) is a dosage compensation mechanism that silences the majority of genes on one X chromosome in each female cell. To characterize epigenetic changes that accompany this process, we measured DNA methylation levels in 45,X patients carrying a single active X chromosome (X(a)), and in normal females, who carry one X(a) and one inactive X (X(i)). Methylated DNA was immunoprecipitated and hybridized to high-density oligonucleotide arrays covering the X chromosome, generating epigenetic profiles of active and inactive X chromosomes. We observed that XCI is accompanied by changes in DNA methylation specifically at CpG islands (CGIs). While the majority of CGIs show increased methylation levels on the X(i), XCI actually results in significant reductions in methylation at 7% of CGIs. Both intra- and inter-genic CGIs undergo epigenetic modification, with the biggest increase in methylation occurring at the promoters of genes silenced by XCI. In contrast, genes escaping XCI generally have low levels of promoter methylation, while genes that show inter-individual variation in silencing show intermediate increases in methylation. Thus, promoter methylation and susceptibility to XCI are correlated. We also observed a global correlation between CGI methylation and the evolutionary age of X-chromosome strata, and that genes escaping XCI show increased methylation within gene bodies. We used our epigenetic map to predict 26 novel genes escaping XCI, and searched for parent-of-origin-specific methylation differences, but found no evidence to support imprinting on the human X chromosome. Our study provides a detailed analysis of the epigenetic profile of active and inactive X chromosomes.

DNA mismatch repair gene MSH6 implicated in determining age at natural menopause.

The length of female reproductive lifespan is associated with multiple adverse outcomes, including breast cancer, cardiovascular disease and infertility. The biological processes that govern the timing of the beginning and end of reproductive life are not well understood. Genetic variants are known to contribute to ∼50% of the variation in both age at menarche and menopause, but to date the known genes explain <15% of the genetic component. We have used genome-wide association in a bivariate meta-analysis of both traits to identify genes involved in determining reproductive lifespan. We observed significant genetic correlation between the two traits using genome-wide complex trait analysis. However, we found no robust statistical evidence for individual variants with an effect on both traits. A novel association with age at menopause was detected for a variant rs1800932 in the mismatch repair gene MSH6 (P = 1.9 × 10(-9)), which was also associated with altered expression levels of MSH6 mRNA in multiple tissues. This study contributes to the growing evidence that DNA repair processes play a key role in ovarian ageing and could be an important therapeutic target for infertility.