Profiling of signal transduction in human memory T cells

Résumé pour un large public: La vaccination a eu un impact énorme sur la santé mondiale. Mais, quel est le principe d'un vaccin? Il est basé sur la 'mémoire immunologique', qui est une particularité exclusive des systèmes immunitaires des organismes évolués. Suite à une infection par un pathogène, des cellules spécialisées de notre système immunitaire (les lymphocytes) le reconnaissent et initient une réaction immunitaire qui a pour but son élimination. Pendant cette réaction se développent aussi des cellules, appelées cellules lymphocytaires mémoire, qui persistent pour longue durée et qui ont la capacité de stimuler une réaction immunitaire très efficace immédiatement après une seconde exposition à ce même pathogène. Ce sont ces cellules mémoires (lymphocytes B et T) qui sont à la base de la 'mémoire immunologique' et qui sont stimulées lors de la vaccination. Chez l'homme, deux populations distinctes des lymphocytes T mémoires ont été identifiées: les cellules centrales (CM) et effectrices (EM) mémoires. Ces populations sont fonctionnellement hétérogènes et exercent des rôles distincts et essentiels dans l'immunité protectrice. Typiquement, les cellules effectrices mémoires sont capables de tuer immédiatement le pathogène tandis que les cellules centrales mémoires sont responsables d'initier une réponse immunitaire complète. Pourtant, les mécanismes biochimiques qui contrôlent les fonctions de ces cellules ont été jusqu'à présent peu étudiés à cause de la faible fréquence de ces cellules et de la quantité limitée de tissus humains disponibles pour les analyses. La compréhension de ces mécanismes est cruciale pour la réalisation de vaccins efficaces et pour le développement de nouveaux médicaments capables de moduler la réponse immunitaire lymphocytaire. Dans cette thèse, nous avons d'abord développé et amélioré une technologie appelée 'protéine array en phase inverse' qui possède un niveau de sensibilité beaucoup plus élevé par rapport aux technologies classiquement utilisées dans l'étude des protéines. Grâce à cette technique, nous avons pu comparer la composition protéique du système de transmission des signaux d'activation des cellules CM et EM humaines. L'analyse de 8 à 13 sujets sains a montré que ces populations des cellules mémoires possèdent un système de signalisation protéique différent. En effet, les cellules EM possèdent, par rapport aux cellules CM, des niveaux réduits d'une protéine régulatrice (appelée c-Cbl) que nous avons démontré comme étant responsable des fonctions spécifiques de ces cellules. En effet, en augmentant artificiellement l'expression de cette protéine régulatrice dans les cellules EM jusqu'au niveau de celui des cellules CM, nous avons induit dans les cellules EM des capacités fonctionnelles caractéristiques des cellules CM. En conclusion, notre étude a identifié, pour la première fois chez l'homme, un mécanisme biochimique qui contrôle les fonctions des populations des cellules mémoires. Résumé en Français: Les cellules mémoires persistent inertes dans l'organisme et produisent des réactions immunitaires rapides et robustes contre les pathogènes précédemment rencontrés. Deux populations distinctes des cellules mémoires ont été identifiées chez l'homme: les cellules centrales (CM) et effectrices (EM) mémoires. Ces populations sont fonctionnellement hétérogènes et exercent des rôles distincts et critiques dans l'immunité protectrice. Les mécanismes biochimiques qui contrôlent leurs fonctions ont été jusqu'à présent peu étudiés, bien que leur compréhension soit cruciale pour le développement des vaccins et des nouveaux traitements/médicaments. Les limites majeures à ces études sont la faible fréquence de ces populations et la quantité limitée de tissus humains disponibles. Dans cette thèse nous avons d'abord développé et amélioré la technologie de 'protéine array en phase inverse' afin d'analyser les molécules de signalisation des cellules mémoires CD4 et CD8 humaines isolées ex vivo. L'excellente sensibilité, la reproductibilité et la linéarité de la détection, ont permis de quantifier des variations d'expression protéiques supérieures à 20% dans un lysat équivalent à 20 cellules. Ensuite, grâce à l'analyse de 8 à 13 sujets sains, nous avons prouvé que les cellules mémoires CD8 ont une composition homogène de leur système de signalisation tandis que les cellules CD4 EM expriment significativement de plus grandes quantités de SLP-76 et des niveaux réduits de c-Cbl, Syk, Fyn et LAT par rapport aux cellules CM. En outre, l'expression réduite du régulateur négatif c-Cbl est corrélée avec l'expression des SLP-76, PI3K et LAT uniquement dans les cellules EM. L'évaluation des propriétés fonctionnelles des cellules mémoires a permis de démontrer que l'expression réduite du c-Cbl dans les cellules EM est associé à une diminution de leur seuil d'activation. En effet, grâce a la technique de transduction cytosolique, nous avons augmenté la quantité de c-Cbl des cellules EM à un niveau comparable à celui des cellules CM et constaté une réduction de la capacité des cellules EM à proliférer et sécréter des cytokines. Ce mécanisme de régulation dépend principalement de l'activité d'ubiquitine ligase de c-Cbl comme démontré par l'impact réduit du mutant enzymatiquement déficient de c-Cbl sur les fonctions de cellules EM. En conclusion, cette thèse identifie c-Cbl comme un régulateur critique des réponses fonctionnelles des populations de cellules T mémoires et fournit, pour la première fois chez l'homme, un mécanisme contrôlant l'hétérogénéité fonctionnelle des ces cellules. De plus, elle valide l'utilisation combinée des 'RPP arrays' et de la transduction cytosolique comme outil puissant d'analyse quantitative et fonctionnel des protéines de signalisation. Summary : Memory cells persist in a quiescent state in the body and mediate rapid and vigorous immune responses toward pathogens previously encountered. Two subsets of memory cells, namely central (CM) and effector (EM) memory cells, have been identified in humans. These subsets display high functional heterogeneity and assert critical and distinct roles in the control of protective immunity. The biochemical mechanisms controlling their functional properties remain so far poorly investigated, although their clarification is crucial for design of effective T-cell vaccine and drug development. Major limitations to these studies lie in the low frequency of memory T cell subsets and the limited amount of human specimen available. In this thesis we first implemented the innovative reverse phase protein array approach to profile 15 signalling components in human CD8 and CD4 memory T cells isolated ex vivo. The high degree of sensitivity, reproducibility and linearity achieved, allowed an excellent quantification of variations in protein expression higher than 20% in as few as 20-cell equivalent per spot. Based on the analysis of 8 to 13 healthy subjects, we showed that CD8 memory cells have a homogeneous composition of their signaling machinery while CD4 EM cells express statistically significant increased amounts of SLP-76 and reduced levels of c- Cbl, Syk, Fyn and LAT as compared to CM cells. Moreover, in EM but not CM cells, reduced expression of negative regulator c-Cbl correlated with the expression of SLP-76, PI3K and LAT. Subsequently, we demonstrated that the higher functional properties and the lower functional threshold of EM cells is associated with reduced expression of c-Cbl. Indeed, by increasing c-Cbl content of EM cells to the same level of CM cells using cytosolic transduction, we impaired their proliferation and cytokine production. This regulatory mechanism was primarily dependent on c-Cbl E3 ubiquitin ligase activity as evidenced by the weaker impact of enzymatically deficient c-Cbl C381A mutant on EM cell functions. Together, these results identify c-Cbl as a critical regulator of the functional responses of memory T cell subsets and provides, for the first time in humans, a mechanism controlling the functional heterogeneity of memory CD4 cells. Moreover it validates the combined use of RPP arrays and cytosolic transduction approaches as a powerful tool to quantitatively analyze signalling proteins and functionally assess their roles.

Development of a harmonised method for the profiling of amphetamine

Résumé : Cette thèse de doctorat est le fruit d'un projet de recherche européen financé par le quatrième programme cadre de la Commission Européenne (DG XII, Standards, Measurement and Testing). Ce projet, dénommé SMT-CT98-2277, a été financé pour la partie suisse par l'Office Fédéral de l'Education et de la Science (OFES, Berne, Suisse). Le but de ce projet était de développer une méthode harmonisée et collaborativement testée pour le profilage des impuretés de l'amphétamine illicite par chromatographie capillaire en phase gazeuse. Le travail a été divisé en sept phases majeures qui concernaient la synthèse de l'amphétamine, l'identification d'impuretés, l'optimisation de la préparation de l'échantillon et du système chromatographique, la variabilité des résultats, l'investigation de méthodes mathématiques pour la classification et la comparaison de profils et finalement l'application de la méthode à des réels échantillons illicites. La méthode résultant de ce travail n'a pas seulement montré que les données étaient interchangeables entre laboratoires mais aussi qu'elle était supérieure en de nombreux points aux méthodes préalablement publiées dans la littérature scientifique. Abstract : This Ph.D. thesis was carried out in parallel to an European project funded by the fourth framework program of the European Commission (DG XII, Standards, Measurement and Testing). This project, named SMT-CT98-2277 was funded, for the Swiss part, by the Federal Office of Education and Science (OFES, Bern, Switzerland). The aim of the project was to develop a harmonised, collaboratively tested method for the impurity profiling of illicit amphetamine by capillary gas chromatography. The work was divided into seven main tasks which deal with the synthesis of amphetamine, identification of impurities, optimization of sample preparation and of the chromatographic system, variability of the results, investigation of numerical methods for the classification and comparison of profiles and finally application of the methodology to real illicit samples. The resulting method has not only shown to produce interchangeable data between different laboratories but was also found to be superior in many aspects to previously published methods.

Implication of DNA methylation, CTCF and BORIS factors in the transcriptional regulation of the human telomerase catalytic subunit hTERT

RESUME La télomérase confère une durée de vie illimitée et est réactivée dans la plupart des cellules tumorales. Sa sous-unité catalytique hTERT est définie comme le facteur limitant pour son activation. De l'identification de facteurs liant la région régulatrice d'hTERT, au rôle de la méthylation de l'ADN et de la modification des histones, de nombreux modèles de régulation ont été suggérés. Cependant, aucun de ces modèles n'a pu expliquer l'inactivation de la télomérase dans la plupart des cellules somatiques et sa réactivation dans la majorité des cellules tumorales. De plus, les observations contradictoires entre le faible niveau d'expression d'ARN messager d'hTERT dans les cellules télomérase-positives et la très forte activité transcriptionnelle du promoteur d'hTERT en transfection restent incomprises. Dans cette étude, nous avons montré que la région proximale du gène hTERT (exon 1 et 2) était impliquée dans la répression de l'activité de son promoteur. Nous avons identifié le facteur CTCF comme étant un inhibiteur du promoteur d'hTERT, en se liant au niveau de son premier exon. La méthylation de l'exon 1 du gène hTERT, couramment observée dans les tumeurs mais pas dans les cellules normales, empêcherait la liaison de CTCF. L'étude du profil de méthylation du promoteur d'hTERT indique qu'une partie du promoteur reste déméthylée et qu'elle semble suffisante pour permettre une faible activité transcriptionnelle du gène hTERT. Ainsi, la méthylation particulière des régions régulatrices d'hTERT inhibe la liaison de CTCF tout en permettant une faible transcription du gène. Cependant, dans certaines cellules tumorales, le promoteur et la région proximale du gène hTERT ne sont pas méthylés. Dans les lignées cellulaires tumorales de tesitcules et d'ovaires, l'inhibition de CTCF est contrée par son paralogue BORIS, qui se lie aussi au niveau de l'exon 1 d'hTERT, mais permet ainsi l'activation du promoteur. L'étude de l'expression du gène BORIS montre qu'il est exclusivement exprimé dans les tissus normaux de testicules et d'ovaires jeunes, ainsi qu'à différents niveaux dans la plupart des tumeurs. Sa transcription est sous le contrôle de deux promoteurs. Le promoteur proximal est régulé par méthylation et un transcrit alternatif majoritaire, délété de l'exon 6, est trouvé lorsque ce promoteur est actif. Tous ces résultats conduisent à un modèle de régulation du gène hTERT qui tient compte du profil épigénétique du gène et qui permet d'expliquer le faible taux de transcription observé in vivo. De plus, l'expression de BORIS dans les cancers et son implication dans l'activation du gène hTERT pourrait permettre de comprendre les phénomènes de dérégulation épigénétique et d'immortalisation qui ont lieu durant la tumorigenèse. SUMMARY Telomerase confers an unlimited lifespan, and is reactivated in most tumor cells. The catalytic subunit of telomerase, hTERT, is defined as the limiting factor for telomerase activity. Between activators and repressors that bind to the hTERT 5' regulatory region, and the role of CpG methylation and histone acetylation, an abundance of regulatory models have been suggested. None of these models can explain the silence of telomerase in most somatic cells and its reactivation in tumor cells. Moreover, the contradictory observations of the low level of hTERT mRNA in telomerase-positive cells and the high transcriptional activity of the hTERT promoter in transfection experiments remain unresolved. In this study, we demonstrated that the proximal exonic region of the hTERT gene (exon 1 and 2) is involved in the inhibition of its promoter. We identified the protein CTCF as the inhibitor of the hTERT promoter, through its binding to the first exon. The methylation of the first exon region, which is often observed in cancer cells but not in noimal cells, represses CTCF binding. Study of hTERT promoter methylation shows a partial demethylation sufficient to activate the transcription of the hTERT gene. Therefore, we demonstrated that the particular methylation profile of the hTERT regulatory sequences inhibits the binding of CTCF, while it allows a low transcription of the gene. Nevertheless, in some tumor cells, the promoter and the proximal exonic region of hTERT are unmethylated. In testicular and ovarian cancer cell lines, CTCF inhibition is counteracted by its BORIS paralogue that also binds the hTERT first exon but allows the promoter activation. The study of BORIS gene regulation showed that this factor is exclusively expressed in normal tissue of testis and ovary of young woman, as well as in almost all tumors with different levels. Two promoters were found to induce its transcription. The proximal promoter was regulated by methylation. Moreover, a major alternative transcript, deleted of the exon 6, is detected when this promoter is active. All these results lead to a model for hTERT regulation that takes into account the epigenetic profile of the gene and provides an explanation for the low transcriptional level observed in vivo. BORIS expression in cancers and its implication in hTERT activation might also permit the understanding of epigenetic deregulation and immortalization phenomena that occur during tumorigenesis.

Can eukaryotic cells monitor the presence of unreplicated DNA?

Completion of DNA replication before mitosis is essential for genome stability and cell viability. Cellular controls called checkpoints act as surveillance mechanisms capable of detecting errors and blocking cell cycle progression to allow time for those errors to be corrected. An important question in the cell cycle field is whether eukaryotic cells possess mechanisms that monitor ongoing DNA replication and make sure that all chromosomes are fully replicated before entering mitosis, that is whether a replication-completion checkpoint exists. From recent studies with smc5–smc6 mutants it appears that yeast cells can enter anaphase without noticing that replication in the ribosomal DNA array was unfinished. smc5–smc6 mutants are proficient in all known cellular checkpoints, namely the S phase checkpoint, DNA-damage checkpoint, and spindle checkpoint, thus suggesting that none of these checkpoints can monitor the presence of unreplicated segments or the unhindered progression of forks in rDNA. Therefore, these results strongly suggest that normal yeast cells do not contain a DNA replication-completion checkpoint.

First systematic chemical profiling of cocaine police seizures in Finland in the framework of an intelligence-led approach

For the first time in Finland, the chemical profiling of cocaine specimens was performed at the National Bureau of Investigation (NBI). The main goals were to determine the chemical composition of cocaine specimens sold in the Finnish market and to study the distribution networks of cocaine in order to provide intelligence related to its trafficking. An analytical methodology enabling through one single GC-MS injection the determination of the added cutting agents (adulterants and diluents), the cocaine purity and the chemical profile (based on the major and minor alkaloids) for each specimen was thus implemented and validated. The methodology was found to be efficient for the discrimination between specimens coming from the same source and specimens coming from different sources. The results highlighted the practical utility of the chemical profiling, especially for supporting the investigation through operational intelligence and improving the knowledge related to the cocaine trafficking through strategic intelligence.

EndoG links Bnip3-induced mitochondrial damage and caspase-independent DNA fragmentation in ischemic cardiomyocytes

Mitochondrial dysfunction, caspase activation and caspase-dependent DNA fragmentation are involved in cell damage in many tissues. However, differentiated cardiomyocytes repress the expression of the canonical apoptotic pathway and their death during ischemia is caspase-independent. The atypical BH3-only protein Bnip3 is involved in the process leading to caspase-independent DNA fragmentation in cardiomyocytes. However, the pathway by which DNA degradation ensues following Bnip3 activation is not resolved. To identify the mechanism involved, we analyzed the interdependence of Bnip3, Nix and EndoG in mitochondrial damage and DNA fragmentation during experimental ischemia in neonatal rat ventricular cardiomyocytes. Our results show that the expression of EndoG and Bnip3 increases in the heart throughout development, while the caspase-dependent machinery is silenced. TUNEL-positive DNA damage, which depends on caspase activity in other cells, is caspase-independent in ischemic cardiomyocytes and ischemia-induced DNA high and low molecular weight fragmentation is blocked by repressing EndoG expression. Ischemia-induced EndoG translocation and DNA degradation are prevented by silencing the expression of Bnip3, but not Nix, or by overexpressing Bcl-xL. These data establish a link between Bnip3 and EndoG-dependent, TUNEL-positive, DNA fragmentation in ischemic cardiomyocytes in the absence of caspases, defining an alternative cell death pathway in postmitotic cells.

Casein Kinase 1γ Ensures Monopolar Growth Polarity under Incomplete DNA Replication Downstream of Cds1 and Calcineurin in Fission Yeast.

Cell polarity is essential for various cellular functions during both proliferative and developmental stages, and it displays dynamic alterations in response to intracellular and extracellular cues. However, the molecular mechanisms underlying spatiotemporal control of polarity transition are poorly understood. Here, we show that fission yeast Cki3 (a casein kinase 1γ homolog) is a critical regulator to ensure persistent monopolar growth during S phase. Unlike the wild type, cki3 mutant cells undergo bipolar growth when S phase is blocked, a condition known to delay transition from monopolar to bipolar growth (termed NETO [new end takeoff]). Consistent with this role, Cki3 kinase activity is substantially increased, and cells lose their viability in the absence of Cki3 upon an S-phase block. Cki3 acts downstream of the checkpoint kinase Cds1/Chk2 and calcineurin, and the latter physically interacts with Cki3. Autophosphorylation in the C terminus is inhibitory toward Cki3 kinase activity, and calcineurin is responsible for its dephosphorylation. Cki3 localizes to the plasma membrane, and this localization requires the palmitoyltransferase complex Erf2-Erf4. Membrane localization is needed not only for proper NETO timing but also for Cki3 kinase activity. We propose that Cki3 acts as a critical inhibitor of cell polarity transition under S-phase arrest.

The Smc5/6 complex is required for dissolution of DNA-mediated sister chromatid linkages

Mitotic chromosome segregation requires the removal of physical connections between sister chromatids. In addition to cohesin and topological entrapments, sister chromatid separation can be prevented by the presence of chromosome junctions or ongoing DNA replication. We will collectively refer to them as DNA-mediated linkages. Although this type of structures has been documented in different DNA replication and repair mutants, there is no known essential mechanism ensuring their timely removal before mitosis. Here, we show that the dissolution of these connections is an active process that requires the Smc5/6 complex, together with Mms21, its associated SUMO-ligase. Failure to remove DNA-mediated linkages causes gross chromosome missegregation in anaphase. Moreover, we show that Smc5/6 is capable to dissolve them in metaphase-arrested cells, thus restoring chromosome resolution and segregation. We propose that Smc5/6 has an essential role in the removal of DNA-mediated linkages to prevent chromosome missegregation and aneuploidy.

Mitochondrial DNA damage and animal longevity: insights from comparative studies

Chemical reactions in living cells are under strict enzyme control and conform to a tightly regulated metabolic program. However, uncontrolled and potentially deleterious endogenous reactions occur, even under physiological conditions. Aging, in this chemical context, could be viewed as an entropic process, the result of chemical side reactions that chronically and cumulatively degrade the function of biological systems. Mitochondria are a main source of reactive oxygen species (ROS) and chemical sidereactions in healthy aerobic tissues and are the only known extranuclear cellular organelles in animal cells that contain their own DNA (mtDNA). ROS can modify mtDNA directly at the sugar-phosphate backbone or at the bases, producing many different oxidatively modified purines and pyrimidines, as well as single and double strand breaks and DNA mutations. In this scenario, natural selection tends to decrease the mitochondrial ROS generation, the oxidative damage to mtDNA, and the mitochondrial mutation rate in long-lived species, in agreement with the mitochondrial oxidative stress theory of aging.

Genome-wide profiling of the cardiac transcriptome after myocardial infarction identifies novel heart-specific long non-coding RNAs.

AIM: Heart disease is recognized as a consequence of dysregulation of cardiac gene regulatory networks. Previously, unappreciated components of such networks are the long non-coding RNAs (lncRNAs). Their roles in the heart remain to be elucidated. Thus, this study aimed to systematically characterize the cardiac long non-coding transcriptome post-myocardial infarction and to elucidate their potential roles in cardiac homoeostasis. METHODS AND RESULTS: We annotated the mouse transcriptome after myocardial infarction via RNA sequencing and ab initio transcript reconstruction, and integrated genome-wide approaches to associate specific lncRNAs with developmental processes and physiological parameters. Expression of specific lncRNAs strongly correlated with defined parameters of cardiac dimensions and function. Using chromatin maps to infer lncRNA function, we identified many with potential roles in cardiogenesis and pathological remodelling. The vast majority was associated with active cardiac-specific enhancers. Importantly, oligonucleotide-mediated knockdown implicated novel lncRNAs in controlling expression of key regulatory proteins involved in cardiogenesis. Finally, we identified hundreds of human orthologues and demonstrate that particular candidates were differentially modulated in human heart disease. CONCLUSION: These findings reveal hundreds of novel heart-specific lncRNAs with unique regulatory and functional characteristics relevant to maladaptive remodelling, cardiac function and possibly cardiac regeneration. This new class of molecules represents potential therapeutic targets for cardiac disease. Furthermore, their exquisite correlation with cardiac physiology renders them attractive candidate biomarkers to be used in the clinic.

Identificação de microssatélites para o mamoeiro por meio da exploração do banco de dados de DNA

Os marcadores microssatélites são ferramentas úteis em diversas análises genéticas em plantas. No caso do mamoeiro (Carica papaya L.), poucos locos de microssatélites foram descritos até o momento. Assim, o objetivo deste trabalho foi explorar a base de dados do GenBank / NCBI (National Center of Biotechnoloy Information) à procura de microssatélites de mamoeiro, visando a seu futuro uso em estudos genéticos e moleculares aplicados ao melhoramento genético. As seqüências foram obtidas no GenBank / NCBI, no formato FASTA, e analisadas para a presença de microssatélites com um mínimo de 20; 7 e 5 repetições dos motivos de mono-, di- e trinucleotídeos, respectivamente, e acima de 4 repetições para tetra- e pentanucleotídeos. Seqüências com mais de 90% de similaridade foram consideradas redundantes e, portanto, eliminadas das análises. Foram analisadas 44.591 seqüências, das quais 3.180 foram não-redundantes e apresentaram 3.947 microssatélites. Desse total, 3.587 foram classificados como microssatélites perfeitos, 8 imperfeitos, 65 interrompidos, 239 compostos-perfeitos, 8 compostos-imperfeitos e 40 compostos-interrompidos. As repetições de di- e trinucleotídeos representaram 65,7 e 14,4% do total de seqüências analisadas, respectivamente. Somente os motivos do tipo AT/TA representaram 44,1% dos microssatélites encontrados. Os motivos mais comuns de tri-, tetra- e pentanucleotídeos foram AAT, AATT e TTTAA, respectivamente. Observou-se que, nas seqüências disponíveis, o genoma do mamoeiro apresenta, em média, um microssatélite a cada 5,65 kb.

Chromosome 7 gain and DNA hypermethylation at the HOXA10 locus are associated with expression of a stem cell related HOX-signature in glioblastoma.

BACKGROUND: HOX genes are a family of developmental genes that are expressed neither in the developing forebrain nor in the normal brain. Aberrant expression of a HOX-gene dominated stem-cell signature in glioblastoma has been linked with increased resistance to chemo-radiotherapy and sustained proliferation of glioma initiating cells. Here we describe the epigenetic and genetic alterations and their interactions associated with the expression of this signature in glioblastoma. RESULTS: We observe prominent hypermethylation of the HOXA locus 7p15.2 in glioblastoma in contrast to non-tumoral brain. Hypermethylation is associated with a gain of chromosome 7, a hallmark of glioblastoma, and may compensate for tumor-driven enhanced gene dosage as a rescue mechanism by preventing undue gene expression. We identify the CpG island of the HOXA10 alternative promoter that appears to escape hypermethylation in the HOX-high glioblastoma. An additive effect of gene copy gain at 7p15.2 and DNA methylation at key regulatory CpGs in HOXA10 is significantly associated with HOX-signature expression. Additionally, we show concordance between methylation status and presence of active or inactive chromatin marks in glioblastoma-derived spheres that are HOX-high or HOX-low, respectively. CONCLUSIONS: Based on these findings, we propose co-evolution and interaction between gene copy gain, associated with a gain of chromosome 7, and additional epigenetic alterations as key mechanisms triggering a coordinated, but inappropriate, HOX transcriptional program in glioblastoma.

DnaSP Version 5. DNA Sequence Polymorphism

DnaSP, DNA Sequence Polymorphism, is a software package for the analysis of nucleotide polymorphism from aligned DNA sequence data. DnaSP can estimate several measures of DNA sequence variation within and between populations (in noncoding, synonymous or nonsynonymous sites, or in various sorts of codon positions), as well as linkage disequilibrium, recombination, gene flow and gene conversion parameters. DnaSP can also carry out several tests of neutrality: Hudson, Kreitman and Aguadé (1987), Tajima (1989), McDonald and Kreitman (1991), Fu and Li (1993), and Fu (1997) tests. Additionally, DnaSP can estimate the confidence intervals of some test-statistics by the coalescent. The results of the analyses are displayed on tabular and graphic form.