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.

The functions of DNA methylation by CcrM in Caulobacter crescentus: a global approach.

DNA methylation is involved in a diversity of processes in bacteria, including maintenance of genome integrity and regulation of gene expression. Here, using Caulobacter crescentus as a model, we exploit genome-wide experimental methods to uncover the functions of CcrM, a DNA methyltransferase conserved in most Alphaproteobacteria. Using single molecule sequencing, we provide evidence that most CcrM target motifs (GANTC) switch from a fully methylated to a hemi-methylated state when they are replicated, and back to a fully methylated state at the onset of cell division. We show that DNA methylation by CcrM is not required for the control of the initiation of chromosome replication or for DNA mismatch repair. By contrast, our transcriptome analysis shows that >10% of the genes are misexpressed in cells lacking or constitutively over-expressing CcrM. Strikingly, GANTC methylation is needed for the efficient transcription of dozens of genes that are essential for cell cycle progression, in particular for DNA metabolism and cell division. Many of them are controlled by promoters methylated by CcrM and co-regulated by other global cell cycle regulators, demonstrating an extensive cross talk between DNA methylation and the complex regulatory network that controls the cell cycle of C. crescentus and, presumably, of many other Alphaproteobacteria.

Polarity and specific sequence requirements of peroxisome proliferator-activated receptor (PPAR)/retinoid X receptor heterodimer binding to DNA. A functional analysis of the malic enzyme gene PPAR response element.

The malic enzyme (ME) gene is a target for both thyroid hormone receptors and peroxisome proliferator-activated receptors (PPAR). Within the ME promoter, two direct repeat (DR)-1-like elements, MEp and MEd, have been identified as putative PPAR response elements (PPRE). We demonstrate that only MEp and not MEd is able to bind PPAR/retinoid X receptor (RXR) heterodimers and mediate peroxisome proliferator signaling. Taking advantage of the close sequence resemblance of MEp and MEd, we have identified crucial determinants of a PPRE. Using reciprocal mutation analyses of these two elements, we show the preference for adenine as the spacing nucleotide between the two half-sites of the PPRE and demonstrate the importance of the two first bases flanking the core DR1 in 5'. This latter feature of the PPRE lead us to consider the polarity of the PPAR/RXR heterodimer bound to its cognate element. We demonstrate that, in contrast to the polarity of RXR/TR and RXR/RAR bound to DR4 and DR5 elements respectively, PPAR binds to the 5' extended half-site of the response element, while RXR occupies the 3' half-site. Consistent with this polarity is our finding that formation and binding of the PPAR/RXR heterodimer requires an intact hinge T region in RXR while its integrity is not required for binding of the RXR/TR heterodimer to a DR4.

Acetylation of cell wall is required for structural integrity of the leaf surface and exerts a global impact on plant stress responses.

The epidermis on leaves protects plants from pathogen invasion and provides a waterproof barrier. It consists of a layer of cells that is surrounded by thick cell walls, which are partially impregnated by highly hydrophobic cuticular components. We show that the Arabidopsis T-DNA insertion mutants of REDUCED WALL ACETYLATION 2 (rwa2), previously identified as having reduced O-acetylation of both pectins and hemicelluloses, exhibit pleiotrophic phenotype on the leaf surface. The cuticle layer appeared diffused and was significantly thicker and underneath cell wall layer was interspersed with electron-dense deposits. A large number of trichomes were collapsed and surface permeability of the leaves was enhanced in rwa2 as compared to the wild type. A massive reprogramming of the transcriptome was observed in rwa2 as compared to the wild type, including a coordinated up-regulation of genes involved in responses to abiotic stress, particularly detoxification of reactive oxygen species and defense against microbial pathogens (e.g., lipid transfer proteins, peroxidases). In accordance, peroxidase activities were found to be elevated in rwa2 as compared to the wild type. These results indicate that cell wall acetylation is essential for maintaining the structural integrity of leaf epidermis, and that reduction of cell wall acetylation leads to global stress responses in Arabidopsis.

Restriction site-associated DNA sequencing, genotyping error estimation and de novo assembly optimization for population genetic inference.

Restriction site-associated DNA sequencing (RADseq) provides researchers with the ability to record genetic polymorphism across thousands of loci for nonmodel organisms, potentially revolutionizing the field of molecular ecology. However, as with other genotyping methods, RADseq is prone to a number of sources of error that may have consequential effects for population genetic inferences, and these have received only limited attention in terms of the estimation and reporting of genotyping error rates. Here we use individual sample replicates, under the expectation of identical genotypes, to quantify genotyping error in the absence of a reference genome. We then use sample replicates to (i) optimize de novo assembly parameters within the program Stacks, by minimizing error and maximizing the retrieval of informative loci; and (ii) quantify error rates for loci, alleles and single-nucleotide polymorphisms. As an empirical example, we use a double-digest RAD data set of a nonmodel plant species, Berberis alpina, collected from high-altitude mountains in Mexico.

A glycine-arginine domain in control of the human MRE11 DNA repair protein.

Human MRE11 is a key enzyme in DNA double-strand break repair and genome stability. Human MRE11 bears a glycine-arginine-rich (GAR) motif that is conserved among multicellular eukaryotic species. We investigated how this motif influences MRE11 function. Human MRE11 alone or a complex of MRE11, RAD50, and NBS1 (MRN) was methylated in insect cells, suggesting that this modification is conserved during evolution. We demonstrate that PRMT1 interacts with MRE11 but not with the MRN complex, suggesting that MRE11 arginine methylation occurs prior to the binding of NBS1 and RAD50. Moreover, the first six methylated arginines are essential for the regulation of MRE11 DNA binding and nuclease activity. The inhibition of arginine methylation leads to a reduction in MRE11 and RAD51 focus formation on a unique double-strand break in vivo. Furthermore, the MRE11-methylated GAR domain is sufficient for its targeting to DNA damage foci and colocalization with gamma-H2AX. These studies highlight an important role for the GAR domain in regulating MRE11 function at the biochemical and cellular levels during DNA double-strand break repair.

Results of a collaborative study of the EDNAP group regarding mitochondrial DNA heteroplasmy and segregation in hair shafts.

A collaborative exercise was carried out by the European DNA Profiling Group (EDNAP) in order to evaluate the distribution of mitochondrial DNA (mtDNA) heteroplasmy amongst the hairs of an individual who displays point heteroplasmy in blood and buccal cells. A second aim of the exercise was to study reproducibility of mtDNA sequencing of hairs between laboratories using differing chemistries, further to the first mtDNA reproducibility study carried out by the EDNAP group. Laboratories were asked to type 2 sections from each of 10 hairs, such that each hair was typed by at least two laboratories. Ten laboratories participated in the study, and a total of 55 hairs were typed. The results showed that the C/T point heteroplasmy observed in blood and buccal cells at position 16234 segregated differentially between hairs, such that some hairs showed only C, others only T and the remainder, C/T heteroplasmy at varying ratios. Additionally, differential segregation of heteroplasmic variants was confirmed in independent extracts at positions 16093 and the poly(C) tract at 302-309, whilst a complete A-G transition was confirmed at position 16129 in one hair. Heteroplasmy was observed at position 16195 on both strands of a single extract from one hair segment, but was not observed in the extracts from any other segment of the same hair. Similarly, heteroplasmy at position 16304 was observed on both strands of a single extract from one hair. Additional variants at positions 73, 249 and the HVII poly(C) region were reported by one laboratory; as these were not confirmed in independent extracts, the possibility of contamination cannot be excluded. Additionally, the electrophoresis and detection equipment used by this laboratory was different to those of the other laboratories, and the discrepancies at position 249 and the HVII poly(C) region appear to be due to reading errors that may be associated with this technology. The results, and their implications for forensic mtDNA typing, are discussed in the light of the biology of hair formation.

Liprin (beta)1 is highly expressed in lymphatic vasculature and is important for lymphatic vessel integrity.

The lymphatic vasculature is important for the regulation of tissue fluid homeostasis, immune response, and lipid absorption, and the development of in vitro models should allow for a better understanding of the mechanisms regulating lymphatic vascular growth, repair, and function. Here we report isolation and characterization of lymphatic endothelial cells from human intestine and show that intestinal lymphatic endothelial cells have a related but distinct gene expression profile from human dermal lymphatic endothelial cells. Furthermore, we identify liprin beta1, a member of the family of LAR transmembrane tyrosine phosphatase-interacting proteins, as highly expressed in intestinal lymphatic endothelial cells in vitro and lymphatic vasculature in vivo, and show that it plays an important role in the maintenance of lymphatic vessel integrity in Xenopus tadpoles.

Use of DNA profiles for investigation using a simulated national DNA database: Part I. Partial SGM Plus profiles.

In traditional criminal investigation, uncertainties are often dealt with using a combination of common sense, practical considerations and experience, but rarely with tailored statistical models. For example, in some countries, in order to search for a given profile in the national DNA database, it must have allelic information for six or more of the ten SGM Plus loci for a simple trace. If the profile does not have this amount of information then it cannot be searched in the national DNA database (NDNAD). This requirement (of a result at six or more loci) is not based on a statistical approach, but rather on the feeling that six or more would be sufficient. A statistical approach, however, could be more rigorous and objective and would take into consideration factors such as the probability of adventitious matches relative to the actual database size and/or investigator's requirements in a sensible way. Therefore, this research was undertaken to establish scientific foundations pertaining to the use of partial SGM Plus loci profiles (or similar) for investigation.

A New Large-DNA-Fragment Delivery System Based on Integrase Activity from an Integrative and Conjugative Element.

During the past few decades, numerous plasmid vectors have been developed for cloning, gene expression analysis, and genetic engineering. Cloning procedures typically rely on PCR amplification, DNA fragment restriction digestion, recovery, and ligation, but increasingly, procedures are being developed to assemble large synthetic DNAs. In this study, we developed a new gene delivery system using the integrase activity of an integrative and conjugative element (ICE). The advantage of the integrase-based delivery is that it can stably introduce a large DNA fragment (at least 75 kb) into one or more specific sites (the gene for glycine-accepting tRNA) on a target chromosome. Integrase recombination activity in Escherichia coli is kept low by using a synthetic hybrid promoter, which, however, is unleashed in the final target host, forcing the integration of the construct. Upon integration, the system is again silenced. Two variants with different genetic features were produced, one in the form of a cloning vector in E. coli and the other as a mini-transposable element by which large DNA constructs assembled in E. coli can be tagged with the integrase gene. We confirmed that the system could successfully introduce cosmid and bacterial artificial chromosome (BAC) DNAs from E. coli into the chromosome of Pseudomonas putida in a site-specific manner. The integrase delivery system works in concert with existing vector systems and could thus be a powerful tool for synthetic constructions of new metabolic pathways in a variety of host bacteria.

Are there two species of Pygmy Shrews (Sorex)? Revisiting the question using DNA sequence data

Pygmy Shrews in North America have variously been considered to be one species (Sorex hoyi) or two species (S. hoyi and S. thompsoni). Currently, only S. hoyi is recognized. In this study, we examine mitochondrial DNA sequence data for the cytochrome b gene to evaluate the level of differentiation and phylogeographic relationships among eleven samples of Pygmy Shrews from across Canada. Pygmy Shrews from eastern Canada (i.e., Ontario, Quebec, New Brunswick, Nova Scotia, and Prince Edward Island) are distinct from Pygmy Shrews from western Canada (Alberta, Yukon) and Alaska. The average level of sequence divergence between these clades (3.3%) falls within the range of values for other recognized pairs of sister species of shrews. A molecular clock based on third position transversion substitutions suggests that these two lineages diverged between 0.44 and 1.67 million years ago. These molecular phylogenetic data. combined with a reinterpretation of previously published morphological data, are suggestive of separate species status for S. hoyi and S. thompsoni as has been previously argued by others. Further analysis of specimens from geographically intermediate areas (e.g., Manitoba. northern Ontario) is required to determine if there is secondary contact and/or introgression between these two putative species.

Tightening of DNA knots by supercoiling facilitates their unknotting by type II DNA topoisomerases.

Using numerical simulations, we compare properties of knotted DNA molecules that are either torsionally relaxed or supercoiled. We observe that DNA supercoiling tightens knotted portions of DNA molecules and accentuates the difference in curvature between knotted and unknotted regions. The increased curvature of knotted regions is expected to make them preferential substrates of type IIA topoisomerases because various earlier experiments have concluded that type IIA DNA topoisomerases preferentially interact with highly curved DNA regions. The supercoiling-induced tightening of DNA knots observed here shows that torsional tension in DNA may serve to expose DNA knots to the unknotting action of type IIA topoisomerases, and thus explains how these topoisomerases could maintain a low knotting equilibrium in vivo, even for long DNA molecules.

Visualisation of human rad52 protein and its complexes with hRad51 and DNA.

The human Rad52 protein stimulates joint molecule formation by hRad51, a homologue of Escherichia coli RecA protein. Electron microscopic analysis of hRad52 shows that it self-associates to form ring structures with a diameter of approximately 10 nm. Each ring contains a hole at its centre. hRad52 binds to single and double-stranded DNA. In the ssDNA-hRad52 complexes, hRad52 was distributed along the length of the DNA, which exhibited a characteristic "beads on a string" appearance. At higher concentrations of hRad52, "super-rings" (approximately 30 nm) were observed and the ssDNA was collapsed upon itself. In contrast, in dsDNA-hRad52 complexes, some regions of the DNA remained protein-free while others, containing hRad52, interacted to form large protein-DNA networks. Saturating concentrations of hRad51 displaced hRad52 from ssDNA, whereas dsDNA-Rad52 complexes (networks) were more resistant to hRad51 invasion and nucleoprotein filament formation. When Rad52-Rad51-DNA complexes were probed with gold-conjugated hRad52 antibodies, the presence of globular hRad52 structures within the Rad51 nucleoprotein filament was observed. These data provide the first direct visualisation of protein-DNA complexes formed by the human Rad51 and Rad52 recombination/repair proteins.

The twist, writhe and overall shape of supercoiled DNA change during counterion-induced transition from a loosely to a tightly interwound superhelix. Possible implications for DNA structure in vivo.

A cryo-electron microscopy study of supercoiled DNA molecules freely suspended in cryo-vitrified buffer was combined with Monte Carlo simulations and gel electrophoretic analysis to investigate the role of intersegmental electrostatic repulsion in determining the shape of supercoiled DNA molecules. It is demonstrated here that a decrease of DNA-DNA repulsion by increasing concentrations of counterions causes a higher fraction of the linking number deficit to be partitioned into writhe. When counterions reach concentrations likely to be present under in vivo conditions, naturally supercoiled plasmids adopt a tightly interwound conformation. In these tightly supercoiled DNA molecules the opposing segments of interwound superhelix seem to directly contact each other. This form of supercoiling, where two DNA helices interact laterally, may represent an important functional state of DNA. In the particular case of supercoiled minicircles (178 bp) the delta Lk = -2 topoisomers undergo a sharp structural transition from almost planar circles in low salt buffers to strongly writhed "figure-eight" conformations in buffers containing neutralizing concentrations of counterions. Possible implications of this observed structural transition in DNA are discussed.

Structure of a multisubunit complex that promotes DNA branch migration.

The RuvA and RuvB proteins of Escherichia coli, which are induced in response to DNA damage, are important in the formation of heteroduplex DNA during genetic recombination and related recombinational repair processes. In vitro studies show that RuvA binds Holiday junctions and acts as a specificity factor that targets the RuvB ATPase, a hexameric ring protein, to the junction. Together, RuvA and RuvB promote branch migration, an ATP-dependent reaction that increases the length of the heteroduplex DNA. Electron microscopic visualization of RuvAB now provides a new insight into the mechanism of this process. We observe the formation of a tripartite protein complex in which RuvA binds the crossover and is sandwiched between two hexameric rings of RuvB. The Holliday junction within this complex adopts a square-planar structure. We propose a molecular model for branch migration, a unique feature of which is the role played by the two oppositely oriented RuvB ring motors.