Characterization of a genomic signature of pregnancy identified in the breast.

The objective of this study was to comprehensively compare the genomic profiles in the breast of parous and nulliparous postmenopausal women to identify genes that permanently change their expression following pregnancy. The study was designed as a two-phase approach. In the discovery phase, we compared breast genomic profiles of 37 parous with 18 nulliparous postmenopausal women. In the validation phase, confirmation of the genomic patterns observed in the discovery phase was sought in an independent set of 30 parous and 22 nulliparous postmenopausal women. RNA was hybridized to Affymetrix HG_U133 Plus 2.0 oligonucleotide arrays containing probes to 54,675 transcripts, scanned and the images analyzed using Affymetrix GCOS software. Surrogate variable analysis, logistic regression, and significance analysis of microarrays were used to identify statistically significant differences in expression of genes. The false discovery rate (FDR) approach was used to control for multiple comparisons. We found that 208 genes (305 probe sets) were differentially expressed between parous and nulliparous women in both discovery and validation phases of the study at an FDR of 10% and with at least a 1.25-fold change. These genes are involved in regulation of transcription, centrosome organization, RNA splicing, cell-cycle control, adhesion, and differentiation. The results provide initial evidence that full-term pregnancy induces long-term genomic changes in the breast. The genomic signature of pregnancy could be used as an intermediate marker to assess potential chemopreventive interventions with hormones mimicking the effects of pregnancy for prevention of breast cancer.

Evidence for different expression profiles for c-Met, EGFR, PTEN and the mTOR pathway in low and high grade endometrial carcinomas in a cohort of consecutive women. Occurrence of PIK3CA and K-Ras mutations and microsatellite instability.

Molecular and genetic investigations in endometrial carcinogenesis may have prognostic and therapeutic implications. We studied the expression of EGFR, c-Met, PTEN and the mTOR signalling pathway (phospho-AKT/phospho-mTOR/phospho-RPS6) in 69 consecutive tumours and 16 tissue microarrays. We also analysed PIK3CA, K-Ras mutations and microsatellite instability (MSI). We distinguished two groups: group 1 (grade 1 and 2 endometrioid cancers) and group 2 (grade 3 endometrioid and type II clear and serous cell cancers). We hypothesised that these histological groups might have different features. We found that a) survival was higher in group 1 with less aggressive tumours (P⟨0.03); b) EGFR (P=0.01), PTEN and the AKT/mTOR/RPS6 signalling pathway were increased in group 1 versus group 2 (P=0.05 for phospho-mTOR); c) conversely, c-Met was higher (P⟨0.03) in group 2 than in group 1; d) In group 1, EGFR was correlated with c-Met, phospho-mTOR, phospho-RPS6 and the global activity of the phospho-AKT/phospho-mTOR/phospho-RPS6 pathway. In group 2, EGFR was correlated only with the phospho-AKT/phospho-mTOR/phospho-RPS6 pathway, whereas c-Met was correlated with PTEN; e) survival was higher for tumours with more than 50% PTEN-positive cells; f) K-RAS and PIK3CA mutations occurred in 10-12% of the available tumours and MSI in 40.4%, with a loss of MLH1 and PMS2 expression. Our results for endometrial cancers provide the first evidence for a difference in status between groups 1 and 2. The patients may benefit from different targeted treatments, anti-EGFR agents and rapamycin derivatives (anti-mTOR) for group 1 and an anti c-MET/ligand complex for group 2.

Expression of MYC, IgM, As Well As Non-Germinal Centre B-Cell Like Immunophenotype and Positive Immunofish Index Predict a Worse Progression Free Survival and Overall Survival in a Series of 670 De Novo Diffuse Large B-Cell Lymphomas Included in Clinical Trials: A GELA Study of the 2003 Program

Introduction: Diffuse large B-cell lymphomas (DLBCL) represent a heterogeneous disease with variable clinical outcome. Identifying phenotypic biomarkers of tumor cells on paraffin sections that predict different clinical outcome remain an important goal that may also help to better understand the biology of this lymphoma. Differentiating non-germinal centre B-cell-like (non-GCB) from Germinal Centre B-cell-like (GCB) DLBCL according to Hans algorithm has been considered as an important immunohistochemical biomarker with prognostic value among patients treated with R-CHOP although not reproducibly found by all groups. Gene expression studies have also shown that IgM expression might be used as a surrogate for the GCB and ABC subtypes with a strong preferential expression of IgM in ABC DLBCL subtype. ImmunoFISH index based on the differential expression of MUM-1, FOXP1 by immunohistochemistry and on the BCL6 rearrangement by FISH has been previously reported (C Copie-Bergman, J Clin Oncol. 2009;27:5573-9) as prognostic in an homogeneous series of DLBCL treated with R-CHOP. In addition, oncogenic MYC protein overexpression by immunohistochemistry may represent an easy tool to identify the consequences of MYC deregulation in DLBCL. Our aim was to analyse by immunohistochemistry the prognostic relevance of MYC, IgM, GCB/nonGCB subtype and ImmunoFISH index in a large series of de novo DLBCL treated with Rituximab (R)-chemotherapy (anthracyclin based) included in the 2003 program of the Groupe d'Etude des Lymphomes de l'Adulte (GELA) trials. Methods: The 2003 program included patients with de novo CD20+ DLBCL enrolled in 6 different LNH-03 GELA trials (LNH-03-1B, -B, -3B, 39B, -6B, 7B) stratifying patients according to age and age-adjusted IPI. Tumor samples were analyzed by immunohistochemistry using CD10, BCL6, MUM1, FOXP1 (according to Barrans threshold), MYC, IgM antibodies on tissue microarrays and by FISH using BCL6 split signal DNA probes. Considering evaluable Hans score, 670 patients were included in the study with 237 (35.4%) receiving intensive R-ACVBP regimen and 433 (64.6%) R-CHOP/R-mini-CHOP. Results: 304 (45.4%) DLBCL were classified as GCB and 366 (54.6%) as non-GCB according to Hans algorithm. 337/567 cases (59.4%) were positive for the ImmunoFISH index (i.e. two out of the three markers positive: MUM1 protein positive, FOXP1 protein Variable or Strong, BCL6 rearrangement). Immunofish index was preferentially positive in the non-GCB subtype (81.3%) compared to the GCB subtype (31.2%), (p<0.001). IgM was recorded as positive in tumor cells in 351/637 (52.4%) DLBCL cases with a preferential expression in non-GCB 195 (53.3%) vs GCB subtype 100(32.9%), p<0.001). MYC was positive in 170/577 (29.5%) cases with a 40% cut-off and in 44/577 (14.2%) cases with a cut-off of 70%. There was no preferential expression of MYC among GCB or non-GCB subtype (p>0.4) for both cut-offs. Progression-free Survival (PFS) was significantly worse among patients with high IPI score (p<0.0001), IgM positive tumor (p<0.0001), MYC positive tumor with a 40% threshold (p<0.001), ImmunoFISH positive index (p<0.002), non-GCB DLBCL subtype (p<0.0001). Overall Survival (OS) was also significantly worse among patients with high IPI score (p<0.0001), IgM positive tumor (p=0.02), MYC positive tumor with a 40% threshold (p<0.01), ImmunoFISH positive index (p=0.02), non-GCB DLBCL subtype (p<0.0001). All significant parameters were included in a multivariate analysis using Cox Model and in addition to IPI, only the GCB/non-GCB subtype according to Hans algorithm predicted significantly a worse PFS among non-GCB subgroup (HR 1.9 [1.3-2.8] p=0.002) as well as a worse OS (HR 2.0 [1.3-3.2], p=0.003). This strong prognostic value of non-GCB subtyping was confirmed considering only patients treated with R- CHOP for PFS (HR 2.1 [1.4-3.3], p=0.001) and for OS (HR 2.3 [1.3-3.8], p=0.002). Conclusion: Our study on a large series of patients included in trials confirmed the relevance of immunohistochemistry as a useful tool to identify significant prognostic biomarkers for clinical use. We show here that IgM and MYC might be useful prognostic biomarkers. In addition, we confirmed in this series the prognostic value of the ImmunoFISH index. Above all, we fully validated the strong and independent prognostic value of the Hans algorithm, daily used by the pathologists to subtype DLBCL.

Lipoxygenase 2 and the early steps of jasmonic acid biosynthesis in "Arabidopsis thaliana"

Résumé Les oxylipines, telles que l'acide jasmonique (AJ ou jasmonate), jouent un rôle central en réponse à la blessure et à la pathogenèse. De nombreuses études ont montré l'importance de la voie canonique du jasmonate lors de la défense des plantes. De plus, un précurseur cyclopentenone de l'AJ, l'acide oxo-phyto-dienoic (OPDA), a été impliqué comme jouant le rôle d'une molécule signal lors de la défense contre certains pathogènes. En utilisant des mutants bloqués dans la biosynthèse de l'acide jasmonique (aos) ou dans sa perception (coi1-1), nous avons cherché à définir dans quelle mesure l'OPDA joue un rôle de signal induisant l'expression génétique en réponse à la blessure chez Arabidopsis. A l'aide de puces à ADN (microarray), nous avons montré que les transcriptomes d'aos et de coi1-1 sont très semblables après blessure, ce qui suggère que les produits d'AOS sont tous perçus via COI1. Pourtant, lorsqu'on analyse les métabolites présents chez ces mutants, une différence est visible, puisque aos n'accumule pas d'AJ, alors que coi1-1 en accumule encore rapidement après blessure. Nous avons étudié la possibilité qu'un mécanisme de régulation post-traductionnelle sur la voie de biosynthèse du jasmonate explique l'accumulation d'AJ chez coi1-1 après blessure. La lipoxygenase 2 (LOX2) est la première enzyme impliquée dans la biosynthèse de l'AJ et est donc une cible potentielle d'un tel mécanisme. Un indice sur la manière dont l'activité LOX pourrait être régulée vient du mutant fou2 (pour fatty acid oxygenation upregcilated 2) dans lequel l'activité LOX ainsi que le niveau d'AJ sont constitutivement élevés. Cette mutation implique un flux de cation dans la régulation de la production de l'AJ. De plus, il a été montré que plusieurs LOXs, dans des organismes autres que des plantes, peuvent lier le calcium. Nous montrons que l'activité LOX requiert l'addition de cations divalents pour être maximale in vitro, et que non seulement le calcium mais aussi le magnésium joue ce rôle. De plus, nous caractérisons un mutant récessif de LOX2 chez Arabidopsis (lox2-1). Ces plantes sont fertiles, et une analyse quantitative montre qu'elles accumulent toujours un peu d'AJ après blessure. Ceci suggère que LOX2 n'est pas la seule LOX impliquée dans la synthèse d'AJ. Aussi les plantes lox2-1 ne sont pas plus sensibles que les plantes de type sauvage lorsqu'elles sont infectées par la moisissure Botrytis cinerea ou lorsqu'elles sont exposées à un détritivore, néanmoins elles sont plus sensibles lorsqu'elles sont offertes en nourriture à un insecte herbivore. Les insectes et les plantes ont co-évolué conjointement, ainsi une plante ne contenant qu'un niveau réduit d'AJ favorise l'insecte. La disponibilité d'un mutant avec un niveau intermédiaire d'AJ va permettre de mieux comprendre pourquoi les plantes produisent autant de jasmonate. Abstract Oxylipins such as jasmonic acid (JA) play central roles in the wound response and during pathogenesis and many studies have confirmed the important role of the canonical jasmonate pathway in plant defense. Moreover, the cyclopentenone precursor of JA, oxo-phytodienoic acid (OPDA), is also thought to function as a signaling molecule in defense towards some pathogens. Its action was reported to depend on a different signal pathway to JA. By using mutants blocked in the biosynthesis (aos) or perception (coil-1) of JA, we investigated to which extend OPDA works as signaling molecule to trigger gene expression in the wound response of Arabidopsis. Using microarrays, we showed that aos and coil-1 transcriptome are similar in response to wounding, suggesting that products of AOS are all perceived by COI1. However, we found a difference between the two mutants at the metobolomic level, since aos is devoid of JA, but coil-1 can still rapidly accumulate JA upon wounding. We investigated the possibility that the post-translational activation of JA biosynthesis could explain the fast accumulation of JA in coil-1 plants upon wounding. Lipoxygenase (LOX) 2 is the first enzyme implicated in JA synthesis and was thus chosen as a potential target for posttranslational regulation. A clue as to how LOX activity might be regulated came from the fatty acid oxygenation upregulated 2 (foul) mutant in which LOX activity and JA levels are elevated. The foul mutant implicates cations flux in the regulation of JA production, and several LOXs in organisms other than plants have been shown to bind calcium. We showed that Arabidopsis LOX requires divalent cations for full activity in vitro, and that not only calcium but also magnesium can play this role. Moreover, a single recessive mutant of AtLOX2 was characterized. These plants are fully fertile. Quantitative oxylipin analysis showed that lox2-1 can still accumulate some JA after wounding, which suggests that LOX2 is not the only LOX involved in JA biosynthesis. lox2-1 plants do not show altered susceptibility to the fungus Botrytis cinerea or to a detritivore, however, they are more susceptible to an insect herbivore. The insect and plants are closely co-evolved and a reduced ability to synthesize JA favors the insect. The availability of a lox2-1 mutant with intermediate JA levels will further help understanding why plants produce elevated JA levels.

Pseudogenes in the ENCODE regions: consensus annotation, analysis of transcription, and evolution.

Arising from either retrotransposition or genomic duplication of functional genes, pseudogenes are "genomic fossils" valuable for exploring the dynamics and evolution of genes and genomes. Pseudogene identification is an important problem in computational genomics, and is also critical for obtaining an accurate picture of a genome's structure and function. However, no consensus computational scheme for defining and detecting pseudogenes has been developed thus far. As part of the ENCyclopedia Of DNA Elements (ENCODE) project, we have compared several distinct pseudogene annotation strategies and found that different approaches and parameters often resulted in rather distinct sets of pseudogenes. We subsequently developed a consensus approach for annotating pseudogenes (derived from protein coding genes) in the ENCODE regions, resulting in 201 pseudogenes, two-thirds of which originated from retrotransposition. A survey of orthologs for these pseudogenes in 28 vertebrate genomes showed that a significant fraction ( approximately 80%) of the processed pseudogenes are primate-specific sequences, highlighting the increasing retrotransposition activity in primates. Analysis of sequence conservation and variation also demonstrated that most pseudogenes evolve neutrally, and processed pseudogenes appear to have lost their coding potential immediately or soon after their emergence. In order to explore the functional implication of pseudogene prevalence, we have extensively examined the transcriptional activity of the ENCODE pseudogenes. We performed systematic series of pseudogene-specific RACE analyses. These, together with complementary evidence derived from tiling microarrays and high throughput sequencing, demonstrated that at least a fifth of the 201 pseudogenes are transcribed in one or more cell lines or tissues.

The G0/G1 switch gene 2 is a novel PPAR target gene.

PPARs (peroxisome-proliferator-activated receptors) alpha, beta/delta and gamma are a group of transcription factors that are involved in numerous processes, including lipid metabolism and adipogenesis. By comparing liver mRNAs of wild-type and PPARalpha-null mice using microarrays, a novel putative target gene of PPARalpha, G0S2 (G0/G1 switch gene 2), was identified. Hepatic expression of G0S2 was up-regulated by fasting and by the PPARalpha agonist Wy14643 in a PPARalpha-dependent manner. Surprisingly, the G0S2 mRNA level was highest in brown and white adipose tissue and was greatly up-regulated during mouse 3T3-L1 and human SGBS (Simpson-Golabi-Behmel syndrome) adipogenesis. Transactivation, gel shift and chromatin immunoprecipitation assays indicated that G0S2 is a direct PPARgamma and probable PPARalpha target gene with a functional PPRE (PPAR-responsive element) in its promoter. Up-regulation of G0S2 mRNA seemed to be specific for adipogenesis, and was not observed during osteogenesis or myogenesis. In 3T3-L1 fibroblasts, expression of G0S2 was associated with growth arrest, which is required for 3T3-L1 adipogenesis. Together, these data indicate that G0S2 is a novel target gene of PPARs that may be involved in adipocyte differentiation.

Differentiation associated regulation of microRNA expression in vivo in human CD8+ T cell subsets.

BACKGROUND: The differentiation of CD8+ T lymphocytes following priming of naïve cells is central in the establishment of the adaptive immune response. Yet, the molecular events underlying this process are not fully understood. MicroRNAs have been recently shown to play a key role in the regulation of haematopoiesis in mouse, but their implication in peripheral lymphocyte differentiation in humans remains largely unknown. METHODS: In order to explore the potential implication of microRNAs in CD8+ T cell differentiation in humans, microRNA expression profiles were analysed using microarrays and quantitative PCR in several human CD8+ T cell subsets defining the major steps of the T cell differentiation pathway. RESULTS: We found expression of a limited set of microRNAs, including the miR-17~92 cluster. Moreover, we reveal the existence of differentiation-associated regulation of specific microRNAs. When compared to naive cells, miR-21 and miR-155 were indeed found upregulated upon differentiation to effector cells, while expression of the miR-17~92 cluster tended to concomitantly decrease. CONCLUSIONS: This study establishes for the first time in a large panel of individuals the existence of differentiation associated regulation of microRNA expression in human CD8+ T lymphocytes in vivo, which is likely to impact on specific cellular functions.

Human-to-bovine jump of Staphylococcus aureus CC8 is associated with the loss of a β-hemolysin converting prophage and the acquisition of a new staphylococcal cassette chromosome.

Staphylococcus aureus can colonize and infect both humans and animals, but isolates from both hosts tend to belong to different lineages. Our recent finding of bovine-adapted S. aureus showing close genetic relationship to the human S. aureus clonal complex 8 (CC8) allowed us to examine the genetic basis of host adaptation in this particular CC. Using total chromosome microarrays, we compared the genetic makeup of 14 CC8 isolates obtained from cows suffering subclinical mastitis, with nine CC8 isolates from colonized or infected human patients, and nine S. aureus isolates belonging to typical bovine CCs. CC8 isolates were found to segregate in a unique group, different from the typical bovine CCs. Within this CC8 group, human and bovine isolates further segregated into three subgroups, among which two contained a mix of human and bovine isolates, and one contained only bovine isolates. This distribution into specific clusters and subclusters reflected major differences in the S. aureus content of mobile genetic elements (MGEs). Indeed, while the mixed human-bovine clusters carried commonly human-associated β-hemolysin converting prophages, the bovine-only isolates were devoid of such prophages but harbored an additional new non-mec staphylococcal cassette chromosome (SCC) unique to bovine CC8 isolates. This composite cassette carried a gene coding for a new LPXTG-surface protein sharing homologies with a protein found in the environmental bacterium Geobacillus thermoglucosidans. Thus, in contrast to human CC8 isolates, the bovine-only CC8 group was associated with the combined loss of β-hemolysin converting prophages and gain of a new SCC probably acquired in the animal environment. Remaining questions are whether the new LPXTG-protein plays a role in bovine colonization or infection, and whether the new SCC could further acquire antibiotic-resistance genes and carry them back to human.

Immunohistochemistry as a valuable tool to assess CD30 expression in peripheral T-cell lymphomas: high correlation with mRNA levels.

The extended use of brentuximab-vedotin was reported for CD30(+) nonanaplastic peripheral T-cell lymphomas (PTCLs) with promising efficacy. CD30 status assessment is thus a critical factor for therapeutic decision, but the reliability of immunohistochemistry (IHC) in evaluating its expression remains to be defined. This prompted us to investigate the correlation between semiquantitative CD30 protein assessment by IHC and messenger RNA (mRNA) assessment by microarrays in a cohort of 376 noncutaneous PTCLs representative of the main entities. By IHC, CD30 expression was heterogeneous across and within entities and significantly associated with large tumor cell size. In addition to 100% anaplastic large-cell lymphomas, 57% of other PTCL entities were CD30-positive at a 5% threshold. CD30 protein expression was highly correlated to mRNA levels. mRNA levels were bimodal, separating high from low CD30-expressing PTCL cases. We conclude that IHC is a valuable tool in clinical practice to assess CD30 expression in PTCLs.

Networking of differentially expressed genes in human cancer cells resistant to methotrexate

BACKGROUND: The need for an integrated view of data obtained from high-throughput technologies gave rise to network analyses. These are especially useful to rationalize how external perturbations propagate through the expression of genes. To address this issue in the case of drug resistance, we constructed biological association networks of genes differentially expressed in cell lines resistant to methotrexate (MTX). METHODS: Seven cell lines representative of different types of cancer, including colon cancer (HT29 and Caco2), breast cancer (MCF-7 and MDA-MB-468), pancreatic cancer (MIA PaCa-2), erythroblastic leukemia (K562) and osteosarcoma (Saos-2), were used. The differential expression pattern between sensitive and MTX-resistant cells was determined by whole human genome microarrays and analyzed with the GeneSpring GX software package. Genes deregulated in common between the different cancer cell lines served to generate biological association networks using the Pathway Architect software. RESULTS: Dikkopf homolog-1 (DKK1) is a highly interconnected node in the network generated with genes in common between the two colon cancer cell lines, and functional validations of this target using small interfering RNAs (siRNAs) showed a chemosensitization toward MTX. Members of the UDP-glucuronosyltransferase 1A (UGT1A) family formed a network of genes differentially expressed in the two breast cancer cell lines. siRNA treatment against UGT1A also showed an increase in MTX sensitivity. Eukaryotic translation elongation factor 1 alpha 1 (EEF1A1) was overexpressed among the pancreatic cancer, leukemia and osteosarcoma cell lines, and siRNA treatment against EEF1A1 produced a chemosensitization toward MTX. CONCLUSIONS: Biological association networks identified DKK1, UGT1As and EEF1A1 as important gene nodes in MTX-resistance. Treatments using siRNA technology against these three genes showed chemosensitization toward MTX.

Computational and statistical analysis of copy number variants in normal and cancer genomes

AbstractAlthough the genomes from any two human individuals are more than 99.99% identical at the sequence level, some structural variation can be observed. Differences between genomes include single nucleotide polymorphism (SNP), inversion and copy number changes (gain or loss of DNA). The latter can range from submicroscopic events (CNVs, at least 1kb in size) to complete chromosomal aneuploidies. Small copy number variations have often no (lethal) consequences to the cell, but a few were associated to disease susceptibility and phenotypic variations. Larger re-arrangements (i.e. complete chromosome gain) are frequently associated with more severe consequences on health such as genomic disorders and cancer. High-throughput technologies like DNA microarrays enable the detection of CNVs in a genome-wide fashion. Since the initial catalogue of CNVs in the human genome in 2006, there has been tremendous interest in CNVs both in the context of population and medical genetics. Understanding CNV patterns within and between human populations is essential to elucidate their possible contribution to disease. But genome analysis is a challenging task; the technology evolves rapidly creating needs for novel, efficient and robust analytical tools which need to be compared with existing ones. Also, while the link between CNV and disease has been established, the relative CNV contribution is not fully understood and the predisposition to disease from CNVs of the general population has not been yet investigated.During my PhD thesis, I worked on several aspects related to CNVs. As l will report in chapter 3, ! was interested in computational methods to detect CNVs from the general population. I had access to the CoLaus dataset, a population-based study with more than 6,000 participants from the Lausanne area. All these individuals were analysed on SNP arrays and extensive clinical information were available. My work explored existing CNV detection methods and I developed a variety of metrics to compare their performance. Since these methods were not producing entirely satisfactory results, I implemented my own method which outperformed two existing methods. I also devised strategies to combine CNVs from different individuals into CNV regions.I was also interested in the clinical impact of CNVs in common disease (chapter 4). Through an international collaboration led by the Centre Hospitalier Universitaire Vaudois (CHUV) and the Imperial College London I was involved as a main data analyst in the investigation of a rare deletion at chromosome 16p11 detected in obese patients. Specifically, we compared 8,456 obese patients and 11,856 individuals from the general population and we found that the deletion was accounting for 0.7% of the morbid obesity cases and was absent in healthy non- obese controls. This highlights the importance of rare variants with strong impact and provides new insights in the design of clinical studies to identify the missing heritability in common disease.Furthermore, I was interested in the detection of somatic copy number alterations (SCNA) and their consequences in cancer (chapter 5). This project was a collaboration initiated by the Ludwig Institute for Cancer Research and involved other groups from the Swiss Institute of Bioinformatics, the CHUV and Universities of Lausanne and Geneva. The focus of my work was to identify genes with altered expression levels within somatic copy number alterations (SCNA) in seven metastatic melanoma ceil lines, using CGH and SNP arrays, RNA-seq, and karyotyping. Very few SCNA genes were shared by even two melanoma samples making it difficult to draw any conclusions at the individual gene level. To overcome this limitation, I used a network-guided analysis to determine whether any pathways, defined by amplified or deleted genes, were common among the samples. Six of the melanoma samples were potentially altered in four pathways and five samples harboured copy-number and expression changes in components of six pathways. In total, this approach identified 28 pathways. Validation with two external, large melanoma datasets confirmed all but three of the detected pathways and demonstrated the utility of network-guided approaches for both large and small datasets analysis.RésuméBien que le génome de deux individus soit similaire à plus de 99.99%, des différences de structure peuvent être observées. Ces différences incluent les polymorphismes simples de nucléotides, les inversions et les changements en nombre de copies (gain ou perte d'ADN). Ces derniers varient de petits événements dits sous-microscopiques (moins de 1kb en taille), appelés CNVs (copy number variants) jusqu'à des événements plus large pouvant affecter des chromosomes entiers. Les petites variations sont généralement sans conséquence pour la cellule, toutefois certaines ont été impliquées dans la prédisposition à certaines maladies, et à des variations phénotypiques dans la population générale. Les réarrangements plus grands (par exemple, une copie additionnelle d'un chromosome appelée communément trisomie) ont des répercutions plus grave pour la santé, comme par exemple dans certains syndromes génomiques et dans le cancer. Les technologies à haut-débit telle les puces à ADN permettent la détection de CNVs à l'échelle du génome humain. La cartographie en 2006 des CNV du génome humain, a suscité un fort intérêt en génétique des populations et en génétique médicale. La détection de différences au sein et entre plusieurs populations est un élément clef pour élucider la contribution possible des CNVs dans les maladies. Toutefois l'analyse du génome reste une tâche difficile, la technologie évolue très rapidement créant de nouveaux besoins pour le développement d'outils, l'amélioration des précédents, et la comparaison des différentes méthodes. De plus, si le lien entre CNV et maladie a été établit, leur contribution précise n'est pas encore comprise. De même que les études sur la prédisposition aux maladies par des CNVs détectés dans la population générale n'ont pas encore été réalisées.Pendant mon doctorat, je me suis concentré sur trois axes principaux ayant attrait aux CNV. Dans le chapitre 3, je détaille mes travaux sur les méthodes d'analyses des puces à ADN. J'ai eu accès aux données du projet CoLaus, une étude de la population de Lausanne. Dans cette étude, le génome de plus de 6000 individus a été analysé avec des puces SNP et de nombreuses informations cliniques ont été récoltées. Pendant mes travaux, j'ai utilisé et comparé plusieurs méthodes de détection des CNVs. Les résultats n'étant pas complètement satisfaisant, j'ai implémenté ma propre méthode qui donne de meilleures performances que deux des trois autres méthodes utilisées. Je me suis aussi intéressé aux stratégies pour combiner les CNVs de différents individus en régions.Je me suis aussi intéressé à l'impact clinique des CNVs dans le cas des maladies génétiques communes (chapitre 4). Ce projet fut possible grâce à une étroite collaboration avec le Centre Hospitalier Universitaire Vaudois (CHUV) et l'Impérial College à Londres. Dans ce projet, j'ai été l'un des analystes principaux et j'ai travaillé sur l'impact clinique d'une délétion rare du chromosome 16p11 présente chez des patients atteints d'obésité. Dans cette collaboration multidisciplinaire, nous avons comparés 8'456 patients atteint d'obésité et 11 '856 individus de la population générale. Nous avons trouvés que la délétion était impliquée dans 0.7% des cas d'obésité morbide et était absente chez les contrôles sains (non-atteint d'obésité). Notre étude illustre l'importance des CNVs rares qui peuvent avoir un impact clinique très important. De plus, ceci permet d'envisager une alternative aux études d'associations pour améliorer notre compréhension de l'étiologie des maladies génétiques communes.Egalement, j'ai travaillé sur la détection d'altérations somatiques en nombres de copies (SCNA) et de leurs conséquences pour le cancer (chapitre 5). Ce projet fut une collaboration initiée par l'Institut Ludwig de Recherche contre le Cancer et impliquant l'Institut Suisse de Bioinformatique, le CHUV et les Universités de Lausanne et Genève. Je me suis concentré sur l'identification de gènes affectés par des SCNAs et avec une sur- ou sous-expression dans des lignées cellulaires dérivées de mélanomes métastatiques. Les données utilisées ont été générées par des puces ADN (CGH et SNP) et du séquençage à haut débit du transcriptome. Mes recherches ont montrées que peu de gènes sont récurrents entre les mélanomes, ce qui rend difficile l'interprétation des résultats. Pour contourner ces limitations, j'ai utilisé une analyse de réseaux pour définir si des réseaux de signalisations enrichis en gènes amplifiés ou perdus, étaient communs aux différents échantillons. En fait, parmi les 28 réseaux détectés, quatre réseaux sont potentiellement dérégulés chez six mélanomes, et six réseaux supplémentaires sont affectés chez cinq mélanomes. La validation de ces résultats avec deux larges jeux de données publiques, a confirmée tous ces réseaux sauf trois. Ceci démontre l'utilité de cette approche pour l'analyse de petits et de larges jeux de données.Résumé grand publicL'avènement de la biologie moléculaire, en particulier ces dix dernières années, a révolutionné la recherche en génétique médicale. Grâce à la disponibilité du génome humain de référence dès 2001, de nouvelles technologies telles que les puces à ADN sont apparues et ont permis d'étudier le génome dans son ensemble avec une résolution dite sous-microscopique jusque-là impossible par les techniques traditionnelles de cytogénétique. Un des exemples les plus importants est l'étude des variations structurales du génome, en particulier l'étude du nombre de copies des gènes. Il était établi dès 1959 avec l'identification de la trisomie 21 par le professeur Jérôme Lejeune que le gain d'un chromosome supplémentaire était à l'origine de syndrome génétique avec des répercussions graves pour la santé du patient. Ces observations ont également été réalisées en oncologie sur les cellules cancéreuses qui accumulent fréquemment des aberrations en nombre de copies (telles que la perte ou le gain d'un ou plusieurs chromosomes). Dès 2004, plusieurs groupes de recherches ont répertorié des changements en nombre de copies dans des individus provenant de la population générale (c'est-à-dire sans symptômes cliniques visibles). En 2006, le Dr. Richard Redon a établi la première carte de variation en nombre de copies dans la population générale. Ces découvertes ont démontrées que les variations dans le génome était fréquentes et que la plupart d'entre elles étaient bénignes, c'est-à-dire sans conséquence clinique pour la santé de l'individu. Ceci a suscité un très grand intérêt pour comprendre les variations naturelles entre individus mais aussi pour mieux appréhender la prédisposition génétique à certaines maladies.Lors de ma thèse, j'ai développé de nouveaux outils informatiques pour l'analyse de puces à ADN dans le but de cartographier ces variations à l'échelle génomique. J'ai utilisé ces outils pour établir les variations dans la population suisse et je me suis consacré par la suite à l'étude de facteurs pouvant expliquer la prédisposition aux maladies telles que l'obésité. Cette étude en collaboration avec le Centre Hospitalier Universitaire Vaudois a permis l'identification d'une délétion sur le chromosome 16 expliquant 0.7% des cas d'obésité morbide. Cette étude a plusieurs répercussions. Tout d'abord elle permet d'effectuer le diagnostique chez les enfants à naître afin de déterminer leur prédisposition à l'obésité. Ensuite ce locus implique une vingtaine de gènes. Ceci permet de formuler de nouvelles hypothèses de travail et d'orienter la recherche afin d'améliorer notre compréhension de la maladie et l'espoir de découvrir un nouveau traitement Enfin notre étude fournit une alternative aux études d'association génétique qui n'ont eu jusqu'à présent qu'un succès mitigé.Dans la dernière partie de ma thèse, je me suis intéressé à l'analyse des aberrations en nombre de copies dans le cancer. Mon choix s'est porté sur l'étude de mélanomes, impliqués dans le cancer de la peau. Le mélanome est une tumeur très agressive, elle est responsable de 80% des décès des cancers de la peau et est souvent résistante aux traitements utilisés en oncologie (chimiothérapie, radiothérapie). Dans le cadre d'une collaboration entre l'Institut Ludwig de Recherche contre le Cancer, l'Institut Suisse de Bioinformatique, le CHUV et les universités de Lausanne et Genève, nous avons séquencés l'exome (les gènes) et le transcriptome (l'expression des gènes) de sept mélanomes métastatiques, effectués des analyses du nombre de copies par des puces à ADN et des caryotypes. Mes travaux ont permis le développement de nouvelles méthodes d'analyses adaptées au cancer, d'établir la liste des réseaux de signalisation cellulaire affectés de façon récurrente chez le mélanome et d'identifier deux cibles thérapeutiques potentielles jusqu'alors ignorées dans les cancers de la peau.

Mechanistic insights into cytosolic molecular chaperones in protein unfolding and disaggregation

Under optimal non-physiological conditions of low concentrations and low temperatures, proteins may spontaneously fold to the native state, as all the information for folding lies in the amino acid sequence of the polypeptide. However, under conditions of stress or high protein crowding as inside cells, a polypeptide may misfold and enter an aggregation pathway resulting in the formation of misfolded conformers and fibrils, which can be toxic and lead to neurodegenerative illnesses, such as Alzheimer's, Parkinson's or Huntington's diseases and aging in general. To avert and revert protein misfolding and aggregation, cells have evolved a set of proteins called molecular chaperones. Here, I focussed on the human cytosolic chaperones Hsp70 (DnaK) and HspllO, and co-chaperone Hsp40 (DnaJ), and the chaperonin CCT (GroEL). The cytosolic molecular chaperones Hsp70s/Hspll0s and the chaperonins are highly upregulated in bacterial and human cells under different stresses and are involved both in the prevention and the reversion of protein misfolding and aggregation. Hsp70 works in collaboration with Hsp40 to reactivate misfolded or aggregated proteins in a strict ATP dependent manner. Chaperonins (CCT and GroEL) also unfold and reactivate stably misfolded proteins but we found that it needed to use the energy of ATP hydrolysis in order to evict over- sticky misfolded intermediates that inhibited the unfoldase catalytic sites. Ill In this study, we initially characterized a particular type of inactive misfolded monomeric luciferase and rhodanese species that were obtained by repeated cycles of freeze-thawing (FT). These stable misfolded monomeric conformers (FT-luciferase and FT-rhodanese) had exposed hydrophobic residues and were enriched with wrong ß-sheet structures (Chapter 2). Using FT-luciferase as substrate, we found that the Hsp70 orthologs, called HspllO (Sse in yeast), acted similarly to Hsp70 as were bona fide ATP- fuelled polypeptide unfoldases and was much more than a mere nucleotide exchange factor, as generally thought. Moreover, we found that HspllO collaborated with Hsp70 in the disaggregation of stable protein aggregates in which Hsp70 and HspllO acted as equal partners that synergistically combined their individual ATP-consuming polypeptide unfoldase activities to reactivate the misfolded/aggregated proteins (Chapter 3). Using FT-rhodanese as substrate, we found that chaperonins (GroEL and CCT) could catalytically reactivate misfolded rhodanese monomers in the absence of ATP. Also, our results suggested that encaging of an unfolding polypeptide inside the GroEL cavity under a GroES cap was not an obligatory step as generally thought (Chapter 4). Further, we investigated the role of Hsp40, a J-protein co-chaperone of Hsp70, in targeting misfolded polypeptides substrates onto Hsp70 for unfolding. We found that even a large excess of monomeric unfolded a-synuclein did not inhibit DnaJ, whereas, in contrast, stable misfolded a-synuclein oligomers strongly inhibited the DnaK-mediated chaperone reaction by way of sequestering the DnaJ co-chaperone. This work revealed that DnaJ could specifically distinguish, and bind potentially toxic stably aggregated species, such as soluble a-synuclein oligomers involved in Parkinson's disease, and with the help of DnaK and ATP convert them into from harmless natively unfolded a-synuclein monomers (chapter 5). Finally, our meta-analysis of microarray data of plant and animal tissues treated with various chemicals and abiotic stresses, revealed possible co-expressions between core chaperone machineries and their co-chaperone regulators. It clearly showed that protein misfolding in the cytosol elicits a different response, consisting of upregulating the synthesis mainly of cytosolic chaperones, from protein misfolding in the endoplasmic reticulum (ER) that elicited a typical unfolded protein response (UPR), consisting of upregulating the synthesis mainly of ER chaperones. We proposed that drugs that best mimicked heat or UPR stress at increasing the chaperone load in the cytoplasm or ER respectively, may prove effective at combating protein misfolding diseases and aging (Chapter 6).  - Dans les conditions optimales de basse concentration et de basse température, les protéines vont spontanément adopter un repliement natif car toutes les informations nécessaires se trouvent dans la séquence des acides aminés du polypeptide. En revanche, dans des conditions de stress ou de forte concentration des protéines comme à l'intérieur d'une cellule, un polypeptide peu mal se replier et entrer dans un processus d'agrégation conduisant à la formation de conformères et de fibrilles qui peuvent être toxiques et causer des maladies neurodégénératives comme la maladie d'Alzheimer, la maladie de Parkinson ou la chorée de Huntington. Afin d'empêcher ou de rectifier le mauvais repliement des protéines, les cellules ont développé des protéines appelées chaperonnes. Dans ce travail, je me suis intéressé aux chaperonnes cytosoliques Hsp70 (DnaK) et HspllO, la co-chaperones Hsp40 (DnaJ), le complexe CCT/TRiC et GroEL. Chez les bactéries et les humains, les chaperonnes cytosoliques Hsp70s/Hspl 10s et les « chaperonines» sont fortement activées par différentes conditions de stress et sont toutes impliquées dans la prévention et la correction du mauvais repliement des protéines et de leur agrégation. Hsp70 collabore avec Hsp40 pour réactiver les protéines agrégées ou mal repliées et leur action nécessite de 1ATP. Les chaperonines (GroEL) déplient et réactivent aussi les protéines mal repliées de façon stable mais nous avons trouvé qu'elles utilisent l'ATP pour libérer les intermédiaires collant et mal repliés du site catalytique de dépliage. Nous avons initialement caractérisé un type particulier de formes stables de luciférase et de rhodanese monomériques mal repliées obtenues après plusieurs cycles de congélation / décongélation répétés (FT). Ces monomères exposaient des résidus hydrophobiques et étaient plus riches en feuillets ß anormaux. Ils pouvaient cependant être réactivés par les chaperonnes Hsp70+Hsp40 (DnaK+DnaJ) et de l'ATP, ou par Hsp60 (GroEL) sans ATP (Chapitre 2). En utilisant la FT-Luciferase comme substrat nous avons trouvé que HspllO (un orthologue de Hsp70) était une authentique dépliase, dépendante strictement de l'ATP. De plus, nous avons trouvé que HspllO collaborait avec Hsp70 dans la désagrégation d'agrégats stables de protéines en combinant leurs activités dépliase consommatrice d'ATP (Chapitre 3). En utilisant la FT-rhodanese, nous avons trouvé que les chaperonines (GroEL et CCT) pouvaient réactiver catalytiquement des monomères mal repliés en absence d'ATP. Nos résultats suggérèrent également que la capture d'un polypeptide en cours de dépliement dans la cavité de GroEL et sous un couvercle du complexe GroES ne serait pas une étape obligatoire du mécanisme, comme il est communément accepté dans la littérature (Chapitre 4). De plus, nous avons étudié le rôle de Hsp40, une co-chaperones de Hsp70, dans l'adressage de substrats polypeptidiques mal repliés vers Hsp70. Ce travail a révélé que DnaJ pouvait différencier et lier des polypeptide mal repliés (toxiques), comme des oligomères d'a-synucléine dans la maladie de Parkinson, et clairement les différencier des monomères inoffensifs d'a-synucléine (Chapitre 5). Finalement une méta-analyse de données de microarrays de tissus végétaux et animaux traités avec différents stress chimiques et abiotiques a révélé une possible co-expression de la machinerie des chaperonnes et des régulateurs de co- chaperonne. Cette meta-analyse montre aussi clairement que le mauvais repliement des protéines dans le cytosol entraîne la synthèse de chaperonnes principalement cytosoliques alors que le mauvais repliement de protéines dans le réticulum endoplasmique (ER) entraine une réponse typique de dépliement (UPR) qui consiste principalement en la synthèse de chaperonnes localisées dans l'ER. Nous émettons l'hypothèse que les drogues qui reproduisent le mieux les stress de chaleur ou les stress UPR pourraient se montrer efficaces dans la lutte contre le mauvais repliement des protéines et le vieillissement (Chapitre 6).

Molecular and social regulation of worker division of labour in fire ants.

Reproductive and worker division of labour (DOL) is a hallmark of social insect societies. Despite a long-standing interest in worker DOL, the molecular mechanisms regulating this process have only been investigated in detail in honey bees, and little is known about the regulatory mechanisms operating in other social insects. In the fire ant Solenopsis invicta, one of the most studied ant species, workers are permanently sterile and the tasks performed are modulated by the worker's internal state (age and size) and the outside environment (social environment), which potentially includes the effect of the queen presence through chemical communication via pheromones. However, the molecular mechanisms underpinning these processes are unknown. Using a whole-genome microarray platform, we characterized the molecular basis for worker DOL and we explored how a drastic change in the social environment (i.e. the sudden loss of the queen) affects global gene expression patterns of worker ants. We identified numerous genes differentially expressed between foraging and nonforaging workers in queenright colonies. With a few exceptions, these genes appear to be distinct from those involved in DOL in bees and wasps. Interestingly, after the queen was removed, foraging workers were no longer distinct from nonforaging workers at the transcriptomic level. Furthermore, few expression differences were detected between queenright and queenless workers when we did not consider the task performed. Thus, the social condition of the colony (queenless vs. queenright) appears to impact the molecular pathways underlying worker task performance, providing strong evidence for social regulation of DOL in S. invicta.

Revisiting sodium and water reabsorption with functional genomics tools.

PURPOSE OF REVIEW: The kidney plays an essential role in maintaining sodium and water balance, thereby controlling the volume and osmolarity of the extracellular body fluids, the blood volume and the blood pressure. The final adjustment of sodium and water reabsorption in the kidney takes place in cells of the distal part of the nephron in which a set of apical and basolateral transporters participate in vectorial sodium and water transport from the tubular lumen to the interstitium and, finally, to the general circulation. According to a current model, the activity and/or cell-surface expression of these transporters is/are under the control of a gene network composed of the hormonally regulated, as well as constitutively expressed, genes. It is proposed that this gene network may include new candidate genes for salt- and water-losing syndromes and for salt-sensitive hypertension. A new generation of functional genomics techniques have recently been applied to the characterization of this gene network. The purpose of this review is to summarize these studies and to discuss the potential of the different techniques for characterization of the renal transcriptome. RECENT FINDINGS: Recently, DNA microarrays and serial analysis of gene expression have been applied to characterize the kidney transcriptome in different in-vivo and in-vitro models. In these studies, a set of new interesting genes potentially involved in the regulation of sodium and water reabsorption by the kidney have been identified and are currently under detailed investigation. SUMMARY: Characterization of the kidney transcriptome is greatly expanding our knowledge of the gene networks involved in multiple kidney functions, including the maintenance of sodium and water homeostasis.

Proteomic and transcriptomic profiling of Staphylococcus aureus surface LPXTG-proteins: correlation with agr genotypes and adherence phenotypes.

Staphylococcus aureus infections involve numerous adhesins and toxins, which expression depends on complex regulatory networks. Adhesins include a family of surface proteins covalently attached to the peptidoglycan via a conserved LPXTG motif. Here we determined the protein and mRNA expression of LPXTG-proteins of S. aureus Newman in time-course experiments, and their relation to fibrinogen adherence in vitro. Experiments were performed with mutants in the global accessory-gene regulator (agr), surface protein A (Spa), and fibrinogen-binding protein A (ClfA), as well as during growth in iron-rich or iron-poor media. Surface proteins were recovered by trypsin-shaving of live bacteria. Released peptides were analyzed by liquid chromatography coupled to tandem mass-spectrometry. To unambiguously identify peptides unique to LPXTG-proteins, the analytical conditions were refined using a reference library of S. aureus LPXTG-proteins heterogeneously expressed in surrogate Lactococcus lactis. Transcriptomes were determined by microarrays. Sixteen of the 18 LPXTG-proteins present in S. aureus Newman were detected by proteomics. Nine LPXTG-proteins showed a bell-shape agr-like expression that was abrogated in agr-negative mutants including Spa, fibronectin-binding protein A (FnBPA), ClfA, iron-binding IsdA, and IsdB, immunomodulator SasH, functionally uncharacterized SasD, biofilm-related SasG and methicillin resistance-related FmtB. However, only Spa and SasH modified their proteomic and mRNA profiles in parallel in the parent and its agr- mutant, whereas all other LPXTG-proteins modified their proteomic profiles independently of their mRNA. Moreover, ClfA became highly transcribed and active in fibrinogen-adherence tests during late growth (24 h), whereas it remained poorly detected by proteomics. On the other hand, iron-regulated IsdA-B-C increased their protein expression by >10-times in iron-poor conditions. Thus, proteomic, transcriptomic, and adherence-phenotype demonstrated differential profiles in S. aureus. Moreover, trypsin peptide signatures suggested differential protein domain exposures in various environments, which might be relevant for anti-adhesin vaccines. A comprehensive understanding of the S. aureus physiology should integrate all three approaches.