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.

Transgenic expression of Ly49A on T cells impairs a specific antitumor response.

Inhibitory MHC receptors determine the reactivity and specificity of NK cells. These receptors can also regulate T cells by modulating TCR-induced effector functions such as cytotoxicity, cytokine production, and proliferation. Here we have assessed the capacity of mouse T cells expressing the inhibitory MHC class I receptor Ly49A to respond to a well-defined tumor Ag in vivo using Ly49A transgenic mice. We find that the presence of Ly49A on the vast majority of lymphocytes prevents the development of a significant Ag-specific CD8+ T cell response and, consequently, the rejection of the tumor. Despite minor alterations in the TCR repertoire of CD8+ T cells in the transgenic lines, precursors of functional tumor-specific CD8+ T cells exist but could not be activated most likely due to a lack of appropriate CD4+ T cell help. Surprisingly, all of these effects are observed in the absence of a known ligand for the Ly49A receptor as defined by its ability to regulate NK cell function. Indeed, we found that the above effects on T cells may be based on a weak interaction of Ly49A with Kb or Db class I molecules. Thus, our data demonstrate that enforced expression of a Ly49A receptor on conventional T cells prevents a specific immune response in vivo and suggest that the functions of T and NK cells are differentially sensitive to the presence of inhibitory MHC class I receptors.

Emergence of secretion-defective sublines of Pseudomonas aeruginosa PAO1 resulting from spontaneous mutations in the vfr global regulatory gene.

Pseudomonas aeruginosa undergoes spontaneous mutation that impairs secretion of several extracellular enzymes during extended cultivation in vitro in rich media, as well as during long-term colonization of the cystic fibrosis lung. A frequent type of strong secretion deficiency is caused by inactivation of the quorum-sensing regulatory gene lasR. Here we analyzed a spontaneously emerging subline of strain PAO1 that exhibited moderate secretion deficiency and partial loss of quorum-sensing control. Using generalized transduction, we mapped the secretion defect to the vfr gene, which is known to control positively the expression of the lasR gene and type II secretion of several proteases. We confirmed this secretion defect by sequencing and complementation of the vfr mutation. In a reconstruction experiment conducted with a 1:1 mixture of wild-type strain PAO1 and a vfr mutant of PAO1, we observed that the vfr mutant had a selective advantage over the wild type after growth in static culture for 4 days. Under these conditions, spontaneous vfr emerged in a strain PAO1 population after four growth cycles, and these mutants accounted for more than 40% of the population after seven cycles. These results suggest that partial or complete loss of quorum sensing and secretion can be beneficial to P. aeruginosa under certain environmental conditions.

Shedding light on proteolytic cleavage of CD44: the responsible sheddase and functional significance of shedding.

CD44 is the major cell-surface receptor for hyaluronan, which is implicated in cell-cell and cell-matrix adhesion, cell migration, and signaling. Studies have shown that CD44-dependent migration requires CD44 to be shed from the cell surface and that matrix metalloproteinase-mediated cleavage may provide an underlying mechanism. However, the full spectrum of proteases that may participate in CD44 shedding has yet to be defined. In this issue, Anderegg et al. demonstrate that ADAM10, but not ADAM17 or MMP14, mediates constitutive shedding of CD44 in human melanoma cells and that knockdown of ADAM10 blocks the antiproliferative activity of the soluble proteolytic cleavage product of CD44.

Validated prediction of clinical outcome in sarcomas and multiple types of cancer on the basis of a gene expression signature related to genome complexity.

Sarcomas are heterogeneous and aggressive mesenchymal tumors. Histological grading has so far been the best predictor for metastasis-free survival, but it has several limitations, such as moderate reproducibility and poor prognostic value for some histological types. To improve patient grading, we performed genomic and expression profiling in a training set of 183 sarcomas and established a prognostic gene expression signature, complexity index in sarcomas (CINSARC), composed of 67 genes related to mitosis and chromosome management. In a multivariate analysis, CINSARC predicts metastasis outcome in the training set and in an independent 127 sarcomas validation set. It is superior to the Fédération Francaise des Centres de Lutte Contre le Cancer grading system in determining metastatic outcome for sarcoma patients. Furthermore, it also predicts outcome for gastrointestinal stromal tumors (GISTs), breast carcinomas and lymphomas. Application of the signature will permit more selective use of adjuvant therapies for people with sarcomas, leading to decreased iatrogenic morbidity and improved outcomes for such individuals.

T cell receptor beta chain gene rearrangements and V beta gene usage in horse cytochrome c-specific T cell hybridomas.

Comparison of T cell receptor alpha and beta-chain genes in murine major histocompatibility complex (MHC) class I and class II-restricted T cell clones and hybridomas recognizing different antigens indicates that no simple correlation exists between the observed antigen/MHC specificity and the expression of certain alpha and beta-chain heterodimers. We have attempted to establish a possible correlation by analyzing T cell receptor beta chain gene rearrangements and V beta gene usage in five T cell hybridomas with identical antigen/MHC specificity and another hybridoma recognizing a different antigenic determinant in association with the same restriction molecule. We report here that in each of the five clones a uniquely rearranged beta chain gene is expressed in combination with at least two different V beta gene segments. The presence of the differently rearranged T cell receptor beta chain genes correlated with the finding of distinct fine specificity pattern of antigen recognition in each of the hybridomas. Interestingly, two hybridomas specific for different epitopes showed identical beta chain D-J rearrangements indicating that the differences might be encoded by the alpha chain gene or/and the V beta gene element.

Human natural Treg microRNA signature: role of microRNA-31 and microRNA-21 in FOXP3 expression.

Treg are the main mediators of dominant tolerance. Their mechanisms of action and applications are subjects of considerable debate currently. However, a human microRNA (miR) Treg signature has not been described yet. We investigated human natural Treg and identified a signature composed of five miR (21, 31, 125a, 181c and 374). Among those, two were considerably under-expressed (miR-31 and miR-125a). We identified a functional target sequence for miR-31 in the 3' untranslated region (3' UTR) of FOXP3 mRNA. Using lentiviral transduction of fresh cord blood T cells, we demonstrated that miR-31 and miR-21 had an effect on FOXP3 expression levels. We showed that miR-31 negatively regulates FOXP3 expression by binding directly to its potential target site in the 3' UTR of FOXP3 mRNA. We next demonstrated that miR-21 acted as a positive, though indirect, regulator of FOXP3 expression. Transduction of the remaining three miR had no direct effect on FOXP3 expression or on the phenotype and will remain the subject of future investigations. In conclusion, not only have we identified and validated a miR signature for human natural Treg, but also unveiled some of the mechanisms by which this signature was related to the control of FOXP3 expression in these cells.

Inhibition of estrogen-responsive gene activation by the retinoid X receptor beta: evidence for multiple inhibitory pathways.

The retinoid X receptor beta (RXR beta; H-2RIIBP) forms heterodimers with various nuclear hormone receptors and binds multiple hormone response elements, including the estrogen response element (ERE). In this report, we show that endogenous RXR beta contributes to ERE binding activity in nuclear extracts of the human breast cancer cell line MCF-7. To define a possible regulatory role of RXR beta regarding estrogen-responsive transcription in breast cancer cells, RXR beta and a reporter gene driven by the vitellogenin A2 ERE were transfected into estrogen-treated MCF-7 cells. RXR beta inhibited ERE-driven reporter activity in a dose-dependent and element-specific fashion. This inhibition occurred in the absence of the RXR ligand 9-cis retinoic acid. The RXR beta-induced inhibition was specific for estrogen receptor (ER)-mediated ERE activation because inhibition was observed in ER-negative MDA-MB-231 cells only following transfection of the estrogen-activated ER. No inhibition of the basal reporter activity was observed. The inhibition was not caused by simple competition of RXR beta with the ER for ERE binding, since deletion mutants retaining DNA binding activity but lacking the N-terminal or C-terminal domain failed to inhibit reporter activity. In addition, cross-linking studies indicated the presence of an auxiliary nuclear factor present in MCF-7 cells that contributed to RXR beta binding of the ERE. Studies using known heterodimerization partners of RXR beta confirmed that RXR beta/triiodothyronine receptor alpha heterodimers avidly bind the ERE but revealed the existence of another triiodothyronine-independent pathway of ERE inhibition. These results indicate that estrogen-responsive genes may be negatively regulated by RXR beta through two distinct pathways.

Comparative modular analysis of gene expression in vertebrate organs.

ABSTRACT: BACKGROUND: The degree of conservation of gene expression between homologous organs largely remains an open question. Several recent studies reported some evidence in favor of such conservation. Most studies compute organs' similarity across all orthologous genes, whereas the expression level of many genes are not informative about organ specificity. RESULTS: Here, we use a modularization algorithm to overcome this limitation through the identification of inter-species co-modules of organs and genes. We identify such co-modules using mouse and human microarray expression data. They are functionally coherent both in terms of genes and of organs from both organisms. We show that a large proportion of genes belonging to the same co-module are orthologous between mouse and human. Moreover, their zebrafish orthologs also tend to be expressed in the corresponding homologous organs. Notable exceptions to the general pattern of conservation are the testis and the olfactory bulb. Interestingly, some co-modules consist of single organs, while others combine several functionally related organs. For instance, amygdala, cerebral cortex, hypothalamus and spinal cord form a clearly discernible unit of expression, both in mouse and human. CONCLUSIONS: Our study provides a new framework for comparative analysis which will be applicable also to other sets of large-scale phenotypic data collected across different species.

Characterization of the PHO1 gene family and the responses to phosphate deficiency of Physcomitrella patens.

PHO1 was previously identified in Arabidopsis (Arabidopsis thaliana) as a protein involved in loading inorganic phosphate (Pi) into the xylem of roots and its expression was associated with the vascular cylinder. Seven genes homologous to AtPHO1 (PpPHO1;1-PpPHO1;7) have been identified in the moss Physcomitrella patens. The corresponding proteins harbor an SPX tripartite domain in the N-terminal hydrophilic portion and an EXS domain in the conserved C-terminal hydrophobic portion, both common features of the plant PHO1 family. Northern-blot analysis showed distinct expression patterns for the PpPHO1 genes, both at the tissue level and in response to phosphate deficiency. Transgenic P. patens expressing the beta-glucuronidase reporter gene under three different PpPHO1 promoters revealed distinct expression profiles in various tissues. Expression of PpPHO1;1 and PpPHO1;7 was specifically induced by Pi starvation. P. patens homologs to the Arabidopsis PHT1, DGD2, SQD1, and APS1 genes also responded to Pi deficiency by increased mRNA levels. Morphological changes associated with Pi deficiency included elongation of caulonemata with inhibition of the formation of side branches, resulting in colonies with greater diameter, but reduced mass compared to Pi-sufficient plants. Under Pi-deficient conditions, P. patens also increased the synthesis of ribonucleases and of an acid phosphatase, and increased the ratio of sulfolipids over phospholipids. These results indicate that P. patens and higher plants share some common strategies to adapt to Pi deficiency, although morphological changes are distinct, and that the PHO1 proteins are well conserved in bryophyte despite the lack of a developed vascular system.

Voltage-gated sodium channel expression in mouse DRG after SNI leads to re-evaluation of projections of injured fibers.

BACKGROUND: Dysregulation of voltage-gated sodium channels (Na(v)s) is believed to play a major role in nerve fiber hyperexcitability associated with neuropathic pain. A complete transcriptional characterization of the different isoforms of Na(v)s under normal and pathological conditions had never been performed on mice, despite their widespread use in pain research. Na(v)s mRNA levels in mouse dorsal root ganglia (DRG) were studied in the spared nerve injury (SNI) and spinal nerve ligation (SNL) models of neuropathic pain. In the SNI model, injured and non-injured neurons were intermingled in lumbar DRG, which were pooled to increase the tissue available for experiments. RESULTS: A strong downregulation was observed for every Na(v)s isoform expressed except for Na(v)1.2; even Na(v)1.3, known to be upregulated in rat neuropathic pain models, was lower in the SNI mouse model. This suggests differences between these two species. In the SNL model, where the cell bodies of injured and non-injured fibers are anatomically separated between different DRG, most Na(v)s were observed to be downregulated in the L5 DRG receiving axotomized fibers. Transcription was then investigated independently in the L3, L4 and L5 DRG in the SNI model, and an important downregulation of many Na(v)s isoforms was observed in the L3 DRG, suggesting the presence of numerous injured neurons there after SNI. Consequently, the proportion of axotomized neurons in the L3, L4 and L5 DRG after SNI was characterized by studying the expression of activating transcription factor 3 (ATF3). Using this marker of nerve injury confirmed that most injured fibers find their cell bodies in the L3 and L4 DRG after SNI in C57BL/6 J mice; this contrasts with their L4 and L5 DRG localization in rats. The spared sural nerve, through which pain hypersensitivity is measured in behavioral studies, mostly projects into the L4 and L5 DRG. CONCLUSIONS: The complex regulation of Na(v)s, together with the anatomical rostral shift of the DRG harboring injured fibers in C57BL/6 J mice, emphasize that caution is necessary and preliminary anatomical experiments should be carried out for gene and protein expression studies after SNI in mouse strains.

Uncoupling phosphate deficiency from its major effects on growth and transcriptome via PHO1 expression in Arabidopsis.

Inorganic phosphate (Pi) is one of the most limiting nutrients for plant growth in both natural and agricultural contexts. Pi-deficiency leads to a strong decrease in shoot growth, and triggers extensive changes at the developmental, biochemical and gene expression levels that are presumably aimed at improving the acquisition of this nutrient and sustaining growth. The Arabidopsis thaliana PHO1 gene has previously been shown to participate in the transport of Pi from roots to shoots, and the null pho1 mutant has all the hallmarks associated with shoot Pi deficiency. We show here that A. thaliana plants with a reduced expression of PHO1 in roots have shoot growth similar to Pi-sufficient plants, despite leaves being strongly Pi deficient. Furthermore, the gene expression profile normally triggered by Pi deficiency is suppressed in plants with low PHO1 expression. At comparable levels of shoot Pi supply, the wild type reduces shoot growth but maintains adequate shoot vacuolar Pi content, whereas the PHO1 underexpressor maintains maximal growth with strongly depleted Pi reserves. Expression of the Oryza sativa (rice) PHO1 ortholog in the pho1 null mutant also leads to plants that maintain normal growth and suppression of the Pi-deficiency response, despite the low shoot Pi. These data show that it is possible to unlink low shoot Pi content with the responses normally associated with Pi deficiency through the modulation of PHO1 expression or activity. These data also show that reduced shoot growth is not a direct consequence of Pi deficiency, but is more likely to be a result of extensive gene expression reprogramming triggered by Pi deficiency.

Leishmania major: differential regulation of the surface metalloprotease in amastigote and promastigote stages.

During its life cycle, the protozoan parasite Leishmania major alternates from an intracellular amastigote form in the mammalian host to a flagellated promastigote form in the insect vector. The expression of the surface metalloprotease (PSP) during differentiation in vitro was investigated by Western and Northern blots, by immunoprecipitation of cells metabolically labeled with [35S]methionine or labeled at the surface with radioactive iodine, and by quantification of the proteolytic activity in substrate-containing polyacrylamide gels. We report that the surface metalloprotease is down-regulated at both the mRNA and the protein level in amastigotes, where it represents less than 1% of the equivalent proteolytic activity detected in promastigotes. A significant amount of mRNA is detected 4 hr after the onset of differentiation. The expression of the protease begins at that time and reaches steady state 8 hr later. The synthesis of PSP precedes the complete morphological differentiation to the promastigote stage and the appearance of the lipophosphoglycan, another major promastigote surface component. In contrast to PSP, a family of mercaptoethanol-activated proteases present in the amastigote exists only at a reduced level in the promastigote. The confinement of the surface metalloprotease to the insect stage of the parasite suggests that it has no physiological function in the parasitism maintenance of mammalian host macrophages.

Expression analyses of three members of the AtPHO1 family reveal differential interactions between signaling pathways involved in phosphate deficiency and the responses to auxin, cytokinin, and abscisic acid.

The PHO1 protein is involved in loading inorganic phosphate (Pi) to the root xylem. Ten genes homologous to AtPHO1 are present in the Arabidopsis thaliana (L.) Heyn genome. From this gene family, transcript levels of only AtPHO1, AtPHO1;H1 and AtPHO1;H10 were increased by Pi-deficiency. While the up-regulation of AtPHO1;H1 and AtPHO1;H10 by Pi deficiency followed the same rapid kinetics and was dependent on the PHR1 transcription factor, phosphite only strongly suppressed the expression of AtPHO1;H1 and had a minor effect on AtPHO1;H10. Addition of sucrose was found to increase transcript levels of both AtPHO1 and AtPHO1;H1 in Pi-sufficient or Pi-deficient plants, but to suppress AtPHO1:H10 under the same conditions. Treatments of plants with auxin or cytokinin had contrasting effect depending on the gene and on the Pi status of the plants. Thus, while both hormones down-regulated expression of AtPHO1 independently of the plant Pi status, auxin and cytokinin up-regulated AtPHO1;H1 and AtPHO1;H10 expression in Pi-sufficient plants and down-regulated expression in Pi-deficient plants. Treatments with abscisic acid inhibited AtPHO1 and AtPHO1;H1 expression in both Pi-sufficient and Pi-deficient plants, but increased AtPHO1;H10 expression under the same conditions. The inhibition of expression by abscisic acid of AtPHO1 and AtPHO1;H1, and of the Pi-starvation responsive genes AtPHT1;1 and AtIPS1, was dependant on the ABI1 type 2C protein phosphatase. These results reveal that various levels of cross talk between the signal transduction pathways to Pi, sucrose and phytohormones are involved in the regulation of expression of the three AtPHO1 homologues.

ExpressionView--an interactive viewer for modules identified in gene expression data.

SUMMARY: ExpressionView is an R package that provides an interactive graphical environment to explore transcription modules identified in gene expression data. A sophisticated ordering algorithm is used to present the modules with the expression in a visually appealing layout that provides an intuitive summary of the results. From this overview, the user can select individual modules and access biologically relevant metadata associated with them. AVAILABILITY: http://www.unil.ch/cbg/ExpressionView. Screenshots, tutorials and sample data sets can be found on the ExpressionView web site.