Doulaye Dembélé

Fuzzy C-means method for clustering microarray data

MOTIVATION Clustering analysis of data from DNA microarray hybridization studies is essential for identifying biologically relevant groups of genes. Partitional clustering methods such as K-means or self-organizing maps assign each gene to a single cluster. However, these methods do not provide information about the influence of a given gene for the overall shape of clusters. Here we apply a fuzzy partitioning method, Fuzzy C-means (FCM), to attribute cluster membership values to genes. RESULTS A major problem in applying the FCM method for clustering microarray data is the choice of the fuzziness parameter m. We show that the commonly used value m = 2 is not appropriate for some data sets, and that optimal values for m vary widely from one data set to another. We propose an empirical method, based on the distribution of distances between genes in a given data set, to determine an adequate value for m. By setting threshold levels for the membership values, genes which are tigthly associated to a given cluster can be selected. Using a yeast cell cycle data set as an example, we show that this selection increases the overall biological significance of the genes within the cluster. AVAILABILITY Supplementary text and Matlab functions are available at http://www-igbmc.u-strasbg.fr/fcm/

Novel insights into the relationships between dendritic cell subsets in human and mouse revealed by genome-wide expression profiling

BackgroundDendritic cells (DCs) are a complex group of cells that play a critical role in vertebrate immunity. Lymph-node resident DCs (LN-DCs) are subdivided into conventional DC (cDC) subsets (CD11b and CD8α in mouse; BDCA1 and BDCA3 in human) and plasmacytoid DCs (pDCs). It is currently unclear if these various DC populations belong to a unique hematopoietic lineage and if the subsets identified in the mouse and human systems are evolutionary homologs. To gain novel insights into these questions, we sought conserved genetic signatures for LN-DCs and in vitro derived granulocyte-macrophage colony stimulating factor (GM-CSF) DCs through the analysis of a compendium of genome-wide expression profiles of mouse or human leukocytes.ResultsWe show through clustering analysis that all LN-DC subsets form a distinct branch within the leukocyte family tree, and reveal a transcriptomal signature evolutionarily conserved in all LN-DC subsets. Moreover, we identify a large gene expression program shared between mouse and human pDCs, and smaller conserved profiles shared between mouse and human LN-cDC subsets. Importantly, most of these genes have not been previously associated with DC function and many have unknown functions. Finally, we use compendium analysis to re-evaluate the classification of interferon-producing killer DCs, lin-CD16+HLA-DR+ cells and in vitro derived GM-CSF DCs, and show that these cells are more closely linked to natural killer and myeloid cells, respectively.ConclusionOur study provides a unique database resource for future investigation of the evolutionarily conserved molecular pathways governing the ontogeny and functions of leukocyte subsets, especially DCs.

Identification of genes associated with tumorigenesis and metastatic potential of hypopharyngeal cancer by microarray analysis

Head and neck squamous cell carcinoma (HNSCC) is the sixth most common cancer among men in the developed world. There is a need, for both clinical and scientific reasons, to find markers to identify patients with aggressive disease as early as possible, and to understand the events leading to malignant transformation and susceptibility to metastasis. We report the first large-scale gene expression analysis of a unique HNSCC location, the hypopharynx. Four normal and 34 tumour samples were analysed with 12 600 gene microarrays. Clusters of differentially expressed genes were identified in the chromosomal regions 3q27.3, 17q21.2–q21.31, 7q11.22–q22.1 and 11q13.1–q13.3, which, interestingly, have already been identified by comparative genomic hybridization (CGH) as major regions of gene amplification. We showed that six overexpressed genes (EIF4G1, DVL3, EPHB4, MCM7, BRMS1 and SART1) located in these regions are indeed amplified. We report 119 genes that are highly differentially expressed between ‘early’ tumours and normal samples. Of these, we validated by quantitative PCR six novel poorly characterized genes. These genes are potential new markers of HNSCC. Comparing patients with relatively nonaggressive and aggressive tumours (without or with clinical evidence of metastasis 3 years after surgery), we identified 164 differentially expressed genes potentially involved in the acquisition of metastatic potential. This study contributes to the understanding of HNSCC, staging patients into prognostic groups and identifying high-risk patients who may benefit from more aggressive treatment.

SOX2 is an oncogene activated by recurrent 3q26.3 amplifications in human lung squamous cell carcinomas.

Squamous cell carcinoma (SCC) of the lung is a frequent and aggressive cancer type. Gene amplifications, a known activating mechanism of oncogenes, target the 3q26-qter region as one of the most frequently gained/amplified genomic sites in SCC of various types. Here, we used array comparative genomic hybridization to delineate the consensus region of 3q26.3 amplifications in lung SCC. Recurrent amplifications occur in 20% of lung SCC (136 tumors in total) and map to a core region of 2 Mb (Megabases) that encompasses SOX2, a transcription factor gene. Intense SOX2 immunostaining is frequent in nuclei of lung SCC, indicating potential active transcriptional regulation by SOX2. Analyses of the transcriptome of lung SCC, SOX2-overexpressing lung epithelial cells and embryonic stem cells (ESCs) reveal that SOX2 contributes to activate ESC-like phenotypes and provide clues pertaining to the deregulated genes involved in the malignant phenotype. In cell culture experiments, overexpression of SOX2 stimulates cellular migration and anchorage-independent growth while SOX2 knockdown impairs cell growth. Finally, SOX2 over-expression in non-tumorigenic human lung bronchial epithelial cells is tumorigenic in immunocompromised mice. These results indicate that the SOX2 transcription factor, a major regulator of stem cell function, is also an oncogene and a driver gene for the recurrent 3q26.33 amplifications in lung SCC.

Misregulation of miR-1 processing is associated with heart defects in myotonic dystrophy.

Myotonic dystrophy is an RNA gain-of-function disease caused by expanded CUG or CCUG repeats, which sequester the RNA binding protein MBNL1. Here we describe a newly discovered function for MBNL1 as a regulator of pre-miR-1 biogenesis and find that miR-1 processing is altered in heart samples from people with myotonic dystrophy. MBNL1 binds to a UGC motif located within the loop of pre-miR-1 and competes for the binding of LIN28, which promotes pre-miR-1 uridylation by ZCCHC11 (TUT4) and blocks Dicer processing. As a consequence of miR-1 loss, expression of GJA1 (connexin 43) and CACNA1C (Cav1.2), which are targets of miR-1, is increased in both DM1- and DM2-affected hearts. CACNA1C and GJA1 encode the main calcium- and gap-junction channels in heart, respectively, and we propose that their misregulation may contribute to the cardiac dysfunctions observed in affected persons.

Misregulated alternative splicing of BIN1 is associated with T tubule alterations and muscle weakness in myotonic dystrophy

Myotonic dystrophy is the most common muscular dystrophy in adults and the first recognized example of an RNA-mediated disease. Congenital myotonic dystrophy (CDM1) and myotonic dystrophy of type 1 (DM1) or of type 2 (DM2) are caused by the expression of mutant RNAs containing expanded CUG or CCUG repeats, respectively. These mutant RNAs sequester the splicing regulator Muscleblind-like-1 (MBNL1), resulting in specific misregulation of the alternative splicing of other pre-mRNAs. We found that alternative splicing of the bridging integrator-1 (BIN1) pre-mRNA is altered in skeletal muscle samples of people with CDM1, DM1 and DM2. BIN1 is involved in tubular invaginations of membranes and is required for the biogenesis of muscle T tubules, which are specialized skeletal muscle membrane structures essential for excitation-contraction coupling. Mutations in the BIN1 gene cause centronuclear myopathy, which shares some histopathological features with myotonic dystrophy. We found that MBNL1 binds the BIN1 pre-mRNA and regulates its alternative splicing. BIN1 missplicing results in expression of an inactive form of BIN1 lacking phosphatidylinositol 5-phosphate–binding and membrane-tubulating activities. Consistent with a defect of BIN1, muscle T tubules are altered in people with myotonic dystrophy, and membrane structures are restored upon expression of the normal splicing form of BIN1 in muscle cells of such individuals. Finally, reproducing BIN1 splicing alteration in mice is sufficient to promote T tubule alterations and muscle weakness, a predominant feature of myotonic dystrophy.

Systematic Gene Expression Mapping Clusters Nuclear Receptors According to Their Function in the Brain

Nuclear receptors (NRs) compose a large family of transcription factors that operate at the interface between genes and environment, acting as sensors and effectors that translate endocrine and metabolic cues into well-defined gene expression programs. We report here on a systematic quantitative and anatomical expression atlas of the 49 NR genes in 104 regions of the adult mouse brain, organized in the interactive MousePat database. MousePat defines NR expression patterns to cellular resolution, a requirement for functional genomic strategies to understand the function of a highly heterogeneous and complex organ such as the brain. Using MousePat data, NR expression patterns can be clustered into anatomical and regulatory networks that delineate the role of NRs in brain functions, like the control of feeding and learning/memory. Mining the MousePat resource will improve the understanding of NR function in the brain and elucidate hierarchical networks that control behavior and whole body homeostasis.

Head and neck squamous cell carcinoma transcriptome analysis by comprehensive validated differential display.

Head and neck squamous cell carcinoma (HNSCC) is common worldwide and is associated with a poor rate of survival. Identification of new markers and therapeutic targets, and understanding the complex transformation process, will require a comprehensive description of genome expression, that can only be achieved by combining different methodologies. We report here the HNSCC transcriptome that was determined by exhaustive differential display (DD) analysis coupled with validation by different methods on the same patient samples. The resulting 820 nonredundant sequences were analysed by high throughput bioinformatics analysis. Human proteins were identified for 73% (596) of the DD sequences. A large proportion (>50%) of the remaining unassigned sequences match ESTs (expressed sequence tags) from human tumours. For the functionally annotated proteins, there is significant enrichment for relevant biological processes, including cell motility, protein biosynthesis, stress and immune responses, cell death, cell cycle, cell proliferation and/or maintenance and transport. Three of the novel proteins (TMEM16A, PHLDB2 and ARHGAP21) were analysed further to show that they have the potential to be developed as therapeutic targets.

Apoptosis and differentiation commitment: novel insights revealed by gene profiling studies in mouse embryonic stem cells.

Mouse embryonic stem (ES) cells remain pluripotent in vitro when grown in the presence of leukemia inhibitory factor (LIF). LIF starvation leads to apoptosis of some of the ES-derived differentiated cells, together with p38α mitogen-activated protein kinase (MAPK) activation. Apoptosis, but not morphological cell differentiation, is blocked by a p38 inhibitor, PD169316. To further understand the mechanism of action of this compound, we have identified its specific targets by microarray studies. We report on the global expression profiles of genes expressed at 3 days upon LIF withdrawal (d3) compared to pluripotent cells and of genes whose expression is modulated at d3 under anti-apoptotic conditions. We showed that at d3 without LIF cells express, earlier than anticipated, specialized cell markers and that when the apoptotic process was impaired, expression of differentiation markers was altered. In addition, functional tests revealed properties of anti-apoptotic proteins not to alter cell pluripotency and a novel role for metallothionein 1 gene, which prevents apoptosis of early differentiated cells.

Mu-opioid receptor activation induces transcriptional plasticity in the central extended amygdala

Addiction develops from the gradual adaptation of the brain to chronic drug exposure, and involves genetic reprogramming of neuronal function. The central extended amygdala (EAc) is a network formed by the central amygdala and the bed nucleus of the stria terminalis. This key site controls drug craving and seeking behaviors, and has not been investigated at the gene regulation level. We used Affymetrix microarrays to analyze transcriptional activity in the murine EAc, with a focus on mu‐opioid receptor‐associated events because these receptors mediate drug reward and dependence. We identified 132 genes whose expression is regulated by a chronic escalating morphine regimen in the EAc from wild‐type but not mu‐opioid receptor knockout mice. These modifications are mostly EAc‐specific. Gene ontology analysis reveals an overrepresentation of neurogenesis, cell growth and signaling protein categories. A separate quantitative PCR analysis of genes in the last of these groups confirms the dysregulation of both orphan (Gpr88) and known (DrD1A, Adora2A, Cnr1, Grm5, Gpr6) G protein‐coupled receptors, scaffolding (PSD95, Homer1) and signaling (Sgk, Cap1) proteins, and neuropeptides (CCK, galanin). These transcriptional modifications do not occur following a single morphine injection, and hence result from long‐term adaptation to excessive mu receptor activation. Proteins encoded by these genes are classically associated with spine modules function in other brain areas, and therefore our data suggest a remodeling of EAc circuits at sites where glutamatergic and monoaminergic afferences interact. Together, mu receptor‐dependent genes identified in this study potentially contribute to drug‐induced neural plasticity, and provide a unique molecular repertoire towards understanding drug craving and relapse.