Integrative analysis of RUNX1 downstream pathways and target genes

BACKGROUND: The RUNX1 transcription factor gene is frequently mutated in sporadic myeloid and lymphoid leukemia through translocation, point mutation or amplification. It is also responsible for a familial platelet disorder with predisposition to acute myeloid leukemia (FPD-AML). The disruption of the largely unknown biological pathways controlled by RUNX1 is likely to be responsible for the development of leukemia. We have used multiple microarray platforms and bioinformatic techniques to help identify these biological pathways to aid in the understanding of why RUNX1 mutations lead to leukemia. RESULTS: Here we report genes regulated either directly or indirectly by RUNX1 based on the study of gene expression profiles generated from 3 different human and mouse platforms. The platforms used were global gene expression profiling of: 1) cell lines with RUNX1 mutations from FPD-AML patients, 2) over-expression of RUNX1 and CBFbeta, and 3) Runx1 knockout mouse embryos using either cDNA or Affymetrix microarrays. We observe that our datasets (lists of differentially expressed genes) significantly correlate with published microarray data from sporadic AML patients with mutations in either RUNX1 or its cofactor, CBFbeta. A number of biological processes were identified among the differentially expressed genes and functional assays suggest that heterozygous RUNX1 point mutations in patients with FPD-AML impair cell proliferation, microtubule dynamics and possibly genetic stability. In addition, analysis of the regulatory regions of the differentially expressed genes has for the first time systematically identified numerous potential novel RUNX1 target genes. CONCLUSION: This work is the first large-scale study attempting to identify the genetic networks regulated by RUNX1, a master regulator in the development of the hematopoietic system and leukemia. The biological pathways and target genes controlled by RUNX1 will have considerable importance in disease progression in both familial and sporadic leukemia as well as therapeutic implications.

Fish possess multiple copies of fgfrl1, the gene for a novel FGF receptor

FGFRL1 is a novel FGF receptor that lacks the intracellular tyrosine kinase domain. While mammals, including man and mouse, possess a single copy of the FGFRL1 gene, fish have at least two copies, fgfrl1a and fgfrl1b. In zebrafish, both genes are located on chromosome 14, separated by about 10 cM. The two genes show a similar expression pattern in several zebrafish tissues, although the expression of fgfrl1b appears to be weaker than that of fgfrl1a. A clear difference is observed in the ovary of Fugu rubripes, which expresses fgfrl1a but not fgfrl1b. It is therefore possible that subfunctionalization has played a role in maintaining the two fgfrl1 genes during the evolution of fish. In human beings, the FGFRL1 gene is located on chromosome 4, adjacent to the SPON2, CTBP1 and MEAEA genes. These genes are also found adjacent to the fgfrl1a gene of Fugu, suggesting that FGFRL1, SPON2, CTBP1 and MEAEA were preserved as a coherent block during the evolution of Fugu and man.

Selecting control genes for RT-QPCR using public microarray data

BACKGROUND: Gene expression analysis has emerged as a major biological research area, with real-time quantitative reverse transcription PCR (RT-QPCR) being one of the most accurate and widely used techniques for expression profiling of selected genes. In order to obtain results that are comparable across assays, a stable normalization strategy is required. In general, the normalization of PCR measurements between different samples uses one to several control genes (e.g. housekeeping genes), from which a baseline reference level is constructed. Thus, the choice of the control genes is of utmost importance, yet there is not a generally accepted standard technique for screening a large number of candidates and identifying the best ones. RESULTS: We propose a novel approach for scoring and ranking candidate genes for their suitability as control genes. Our approach relies on publicly available microarray data and allows the combination of multiple data sets originating from different platforms and/or representing different pathologies. The use of microarray data allows the screening of tens of thousands of genes, producing very comprehensive lists of candidates. We also provide two lists of candidate control genes: one which is breast cancer-specific and one with more general applicability. Two genes from the breast cancer list which had not been previously used as control genes are identified and validated by RT-QPCR. Open source R functions are available at http://www.isrec.isb-sib.ch/~vpopovic/research/ CONCLUSION: We proposed a new method for identifying candidate control genes for RT-QPCR which was able to rank thousands of genes according to some predefined suitability criteria and we applied it to the case of breast cancer. We also empirically showed that translating the results from microarray to PCR platform was achievable.

The evolutionary diversification of parrots supports a taxon pulse model with multiple trans-oceanic dispersal events and local radiations

Vicariance is thought to have played a major role in the evolution of modern parrots. However, as the relationships especially of the African taxa remained mostly unresolved, it has been difficult to draw firm conclusions about the roles of dispersal and vicariance. Our analyses using the broadest taxon sampling of old world parrots ever based on 3219 bp of three nuclear genes revealed well-resolved and congruent phylogenetic hypotheses. Agapornis of Africa and Madagascar was found to be the sister group to Loriculus of Australasia and Indo-Malayasia and together they clustered with the Australasian Loriinae, Cyclopsittacini and Melopsittacus. Poicephalus and Psittacus from mainland Africa formed the sister group Of the Neotropical Arini and Coracopsis from Madagascar and adjacent islands may be the closest relative of Psittrichas from New Guinea. These biogeographic relationships are best explained by independent colonization of the African continent via trans-oceanic dispersal from Australasia and Antarctica in the Paleogene following what may have been vicariance events in the late Cretaceous and/or early Paleogene. Our data support a taxon pulse model for the diversification of parrots whereby trans-oceanic dispersal played a more important role than previously thought and was the prerequisite for range expansion into new continents. (C) 2009 Elsevier Inc. All rights reserved

Unravelling the multiple roles of FUS in RNA processing and on ALS pathogenesis

Two genes with related functions in RNA biogenesis were recently reported in patients with familial ALS: the FUS/TLS gene at the ALS6 locus and the TARDBP/TDP-43 gene at the ALS10 locus [1, 2]. FUS has been implicated to function in several steps of gene expression, including transcription regulation [3], RNA splicing [4, 5], mRNA transport in neurons [6] and, interestingly, in microRNA (miRNA) processing [7]. The goal of this project is to identify the molecular mechanisms leading to the development of FUS mutations-associated ALS. Specifically, we want to test the hypothesis that these FUS mutations misregulate miRNA levels that in turn affect the expression of genes critical for motor neuron survival. In addition we want to test whether misregulation of the miRNA profile is a common feature in ALS. We have performed immunoprecipitations from total extracts of 293T cells expressing FLAG-tagged FUS to characterize its interactome by mass spectrometry. This proteomic study not only revealed a strong interaction of FUS with splicing factors, but shows that FUS might be involved in many, quite different pathways. To map which parts of the FUS protein contribute to the interaction with splicing factors, we have performed a set of experiments with a series of missense and deletion mutants. With this approach, we will not only gain information on the binding partners of FUS along with a map of the required domains for the interactions, but it will also help to unravel whether certain ALS-associated FUS mutations lead to a loss or gain of function due to gain or loss of interactors. Additionally, we have performed quantitative interactomics using SILAC to identify interactome differences of ALS-associated FUS mutants. To this end we have performed immunoprecipitations of total extract from 293T cells, stably transduced with constructs expressing wild-type FUS-FLAG as well as three different ALS-associated mutants (G156E, R244C, P525L). First results indicate striking differences in the interactome with certain RNA binding proteins. We are now validating these candidates in order to reveal the importance of these differential interactions in the context of ALS.

Cycling in the nucleus: regulation of RNA 3' processing and nuclear organization of replication-dependent histone genes.

The histones which pack new DNA during the S phase of animal cells are made from mRNAs that are cleaved at their 3' end but not polyadenylated. Some of the factors used in this reaction are unique to it while others are shared with the polyadenylation process that generates all other mRNAs. Recent work has begun to shed light on how the cell manages the assignment of these common components to the two 3' processing systems, and how it achieves their cell cycle-regulation and recruitment to the histone pre-mRNA. Moreover, recent and older findings reveal multiple connections between the nuclear organization of histone genes, their transcription and 3' end processing as well as the control of cell proliferation.

Deletion of exon 8 from the EXT1 gene causes multiple osteochondromas (MO) in a family with three affected members.

Multiple osteochondromas (also called hereditary multiple exostoses) is an autosomal dominant disorder characterized by multiple cartilaginous tumors, which are caused by mutations in the genes for exostosin-1 (EXT1) and exostosin-2 (EXT2). The goal of this study was to elucidate the genetic alterations in a family with three affected members. Isolation of RNA from the patients' blood followed by reverse transcription and PCR amplification of selected fragments showed that the three patients lack a specific region of 90 bp from their EXT1 mRNA. This region corresponds to the sequence of exon 8 from the EXT1 gene. No splice site mutation was found around exon 8. However, long-range PCR amplification of the region from intron 7 to intron 8 indicated that the three patients contain a deletion of 4318 bp, which includes exon 8 and part of the flanking introns. There is evidence that the deletion was caused by non-homologous end joining because the breakpoints are not located within a repetitive element, but contain multiple copies of the deletion hotspot sequence TGRRKM. Exon 8 encodes part of the active site of the EXT1 enzyme, including the DXD signature of all UDP-sugar glycosyltransferases. It is conceivable that the mutant protein exerts a dominant negative effect on the activity of the EXT glycosyltransferase since it might interact with normal copies of the enzyme to form an inactive hetero-oligomeric complex. We suggest that sequencing of RNA might be superior to exome sequencing to detect short deletions of a single exon.

Multiple regulatory variants located in cell type-specific enhancers within the PKP2 locus form major risk and protective haplotypes for canine atopic dermatitis in German shepherd dogs.

BACKGROUND Canine atopic dermatitis (CAD) is a chronic inflammatory skin disease triggered by allergic reactions involving IgE antibodies directed towards environmental allergens. We previously identified a ~1.5 Mb locus on canine chromosome 27 associated with CAD in German shepherd dogs (GSDs). Fine-mapping indicated association closest to the PKP2 gene encoding plakophilin 2. RESULTS Additional genotyping and association analyses in GSDs combined with control dogs from five breeds with low-risk for CAD revealed the top SNP 27:19,086,778 (p = 1.4 × 10(-7)) and a rare ~48 kb risk haplotype overlapping the PKP2 gene and shared only with other high-risk CAD breeds. We selected altogether nine SNPs (four top-associated in GSDs and five within the ~48 kb risk haplotype) that spanned ~280 kb forming one risk haplotype carried by 35 % of the GSD cases and 10 % of the GSD controls (OR = 5.1, p = 5.9 × 10(-5)), and another haplotype present in 85 % of the GSD cases and 98 % of the GSD controls and conferring a protective effect against CAD in GSDs (OR = 0.14, p = 0.0032). Eight of these SNPs were analyzed for transcriptional regulation using reporter assays where all tested regions exerted regulatory effects on transcription in epithelial and/or immune cell lines, and seven SNPs showed allelic differences. The DNA fragment with the top-associated SNP 27:19,086,778 displayed the highest activity in keratinocytes with 11-fold induction of transcription by the risk allele versus 8-fold by the control allele (pdifference = 0.003), and also mapped close (~3 kb) to an ENCODE skin-specific enhancer region. CONCLUSIONS Our experiments indicate that multiple CAD-associated genetic variants located in cell type-specific enhancers are involved in gene regulation in different cells and tissues. No single causative variant alone, but rather multiple variants combined in a risk haplotype likely contribute to an altered expression of the PKP2 gene, and possibly nearby genes, in immune and epithelial cells, and predispose GSDs to CAD.

Novel multiple sclerosis susceptibility loci implicated in epigenetic regulation.

We conducted a genome-wide association study (GWAS) on multiple sclerosis (MS) susceptibility in German cohorts with 4888 cases and 10,395 controls. In addition to associations within the major histocompatibility complex (MHC) region, 15 non-MHC loci reached genome-wide significance. Four of these loci are novel MS susceptibility loci. They map to the genes L3MBTL3, MAZ, ERG, and SHMT1. The lead variant at SHMT1 was replicated in an independent Sardinian cohort. Products of the genes L3MBTL3, MAZ, and ERG play important roles in immune cell regulation. SHMT1 encodes a serine hydroxymethyltransferase catalyzing the transfer of a carbon unit to the folate cycle. This reaction is required for regulation of methylation homeostasis, which is important for establishment and maintenance of epigenetic signatures. Our GWAS approach in a defined population with limited genetic substructure detected associations not found in larger, more heterogeneous cohorts, thus providing new clues regarding MS pathogenesis.

Analysis of Plasminogen Genetic Variants in Multiple Sclerosis Patients.

Multiple sclerosis (MS) is a prevalent neurological disease of complex etiology. Here, we describe the characterization of a multi-incident MS family that nominated a rare missense variant (p.G420D) in plasminogen (PLG) as a putative genetic risk factor for MS. Genotyping of PLG p.G420D (rs139071351) in 2160 MS patients, and 886 controls from Canada, identified 10 additional probands, two sporadic patients and one control with the variant. Segregation in families harboring the rs139071351 variant, identified p.G420D in 26 out of 30 family members diagnosed with MS, 14 unaffected parents, and 12 out of 30 family members not diagnosed with disease. Despite considerably reduced penetrance, linkage analysis supports cosegregation of PLG p.G420D and disease. Genotyping of PLG p.G420D in 14446 patients, and 8797 controls from Canada, France, Spain, Germany, Belgium, and Austria failed to identify significant association with disease (P = 0.117), despite an overall higher prevalence in patients (OR = 1.32; 95% CI = 0.93-1.87). To assess whether additional rare variants have an effect on MS risk, we sequenced PLG in 293 probands, and genotyped all rare variants in cases and controls. This analysis identified nine rare missense variants, and although three of them were exclusively observed in MS patients, segregation does not support pathogenicity. PLG is a plausible biological candidate for MS owing to its involvement in immune system response, blood-brain barrier permeability, and myelin degradation. Moreover, components of its activation cascade have been shown to present increased activity or expression in MS patients compared to controls; further studies are needed to clarify whether PLG is involved in MS susceptibility.