Patrick Beaulieu

Genome-wide analysis of transcript isoform variation in humans

We have performed a genome-wide analysis of common genetic variation controlling differential expression of transcript isoforms in the CEU HapMap population using a comprehensive exon tiling microarray covering 17,897 genes. We detected 324 genes with significant associations between flanking SNPs and transcript levels. Of these, 39% reflected changes in whole gene expression and 55% reflected transcript isoform changes such as splicing variants (exon skipping, alternative splice site use, intron retention), differential 5′ UTR (initiation of transcription) use, and differential 3′ UTR (alternative polyadenylation) use. These results demonstrate that the regulatory effects of genetic variation in a normal human population are far more complex than previously observed. This extra layer of molecular diversity may account for natural phenotypic variation and disease susceptibility.

Global patterns of cis variation in human cells revealed by high-density allelic expression analysis.

Cis-acting variants altering gene expression are a source of phenotypic differences. The cis-acting components of expression variation can be identified through the mapping of differences in allelic expression (AE), which is the measure of relative expression between two allelic transcripts. We generated a map of AE associated SNPs using quantitative measurements of AE on Illumina Human1M BeadChips. In 53 lymphoblastoid cell lines derived from donors of European descent, we identified common cis variants affecting 30% (2935/9751) of the measured RefSeq transcripts at 0.001 permutation significance. The pervasive influence of cis-regulatory variants, which explain 50% of population variation in AE, extend to full-length transcripts and their isoforms as well as to unannotated transcripts. These strong effects facilitate fine mapping of cis-regulatory SNPs, as demonstrated by dissection of heritable control of transcripts in the systemic lupus erythematosus–associated C8orf13-BLK region in chromosome 8. The dense collection of associations will facilitate large-scale isolation of cis-regulatory SNPs.

Population genomics in a disease targeted primary cell model

The common genetic variants associated with complex traits typically lie in noncoding DNA and may alter gene regulation in a cell type-specific manner. Consequently, the choice of tissue or cell model in the dissection of disease associations is important. We carried out an expression quantitative trait loci (eQTL) study of primary human osteoblasts (HOb) derived from 95 unrelated donors of Swedish origin, each represented by two independently derived primary lines to provide biological replication. We combined our data with publicly available information from a genome-wide association study (GWAS) of bone mineral density (BMD). The top 2000 BMD-associated SNPs (P < approximately 10(-3)) were tested for cis-association of gene expression in HObs and in lymphoblastoid cell lines (LCLs) using publicly available data and showed that HObs have a significantly greater enrichment (threefold) of converging cis-eQTLs as compared to LCLs. The top 10 BMD loci with SNPs showing strong cis-effects on gene expression in HObs (P = 6 x 10(-10) - 7 x 10(-16)) were selected for further validation using a staged design in two cohorts of Caucasian male subjects. All 10 variants were tested in the Swedish MrOS Cohort (n = 3014), providing evidence for two novel BMD loci (SRR and MSH3). These variants were then tested in the Rotterdam Study (n = 2090), yielding converging evidence for BMD association at the 17p13.3 SRR locus (P(combined) = 5.6 x 10(-5)). The cis-regulatory effect was further fine-mapped to the proximal promoter of the SRR gene (rs3744270, r(2) = 0.5, P = 2.6 x 10(-15)). Our results suggest that primary cells relevant to disease phenotypes complement traditional approaches for prioritization and validation of GWAS hits for follow-up studies.

Characterization of the BclI Polymorphism in the Glucocorticoid Receptor Gene

Glucocorticoids (GCs) have a major antiproliferative effect, which has led to the use of their synthetic homologs for immunosuppression, treatment of inflammation, and induction of cytotoxicity (1)(2). GCs exert their effect by binding to an intracellular GC receptor (GR), forming a complex that translocates to the nucleus, where GCs then regulate the expression of target genes interacting with promoter GC-responsive elements (1). The different GR forms, resulting from GR gene variability, can affect the regulation of many biological functions, such as hypothalamic-pituitary-adrenal axis regulation and GC responsiveness, thereby underlying susceptibility to many diseases. Indeed, GR mutations have been associated with altered cardiovascular function, metabolic disturbances, and hematologic malignancies (3)(4)(5). Likewise, functional GR variability might affect the therapeutic response to corticosteroid drugs (5). Identification of different GR gene variants may thus be helpful in assessing the role of the GR gene in disease susceptibility or in adjudging predisposition to corticosteroid-associated adverse drug reactions. Several polymorphisms of the GR gene, which might have an impact on GC sensitivity, have been reported (6)(7)(8). Among these, the Bcl I polymorphism was identified by Southern blotting using human GR cDNA-specific probes (9) that identified two alleles with fragment lengths of 4.5 and 2.3 kb. Several clinical investigations have subsequently suggested that this GR polymorphism is linked to altered GR function (6)(10)(11)(12)(13)(14)(15). An association between the Bcl I polymorphism and changes in tissue-specific corticosteroid sensitivity, as well as with poor feedback regulation of the hypothalamic-pituitary-adrenal axis, has been reported (6)(10). This was further documented by association of the Bcl I polymorphism with abdominal obesity (11)(12), insulin resistance (6)(13), and development of an atherogenic profile (6)(14). Similarly, …

DNA Variants in Region for Noncoding Interfering Transcript of Dihydrofolate Reductase Gene and Outcome in Childhood Acute Lymphoblastic Leukemia

Purpose: Dihydrofolate reductase (DHFR) is the major target of methotrexate, a key component in childhood acute lymphoblastic leukemia (ALL) treatment. We recently reported an association of DHFR promoter polymorphisms with ALL outcome. Lower event-free survival correlated with haplotype *1, defined by A317 and C1610 alleles. Haplotype *1 was also associated higher DHFR expression. Experimental Design: Here, we analyzed adjacent 400-bp region participating in DHFR regulation as both a major promoter and a noncoding minor transcript. Results: Six polymorphisms were identified, of which five were single nucleotide polymorphisms and one was length polymorphism composed of variable number of 9-bp elements and 9-bp insertion/deletion. Haplotype analysis including all promoter polymorphisms revealed diversification of haplotype *1 into five subtypes (*1a-*1e). DNA variations of major promoter/noncoding transcript region and haplotype *1 subtypes were subsequently analyzed for the association with ALL outcome. Lower event-free survival was associated with an A allele of G308A polymorphism (P = 0.02) and with *1b haplotype (P = 0.01). This association was particularly striking in high-risk patients (P = 0.001) and was subsequently confirmed in independent patient cohort (P = 0.02). Haplotype *1b was the only haplotype *1 subtype associated with higher mRNA levels. Conclusions: The study provides a new insight into DHFR regulatory variations predisposing to an event in ALL patients. (Clin Cancer Res 2009;15(22):69318)

Functional impact of sequence variation in the promoter region of TGFB1.

Pathological deregulation of the transforming growth factor, beta 1 (TGFB1) pathway has been implicated in the development of several major diseases, including cancers. Regulatory variation in the TGFB1 gene may lead to altered TGFB1 expression and activity, and thus, modulate an individuals susceptibility to disease. Here, we performed a study of the functional relevance of cis‐acting regulatory variation in the proximal promoter region of the TGFB1 gene. In a previous study, 9 promoter polymorphisms were identified in the 2kb region upstream of the transcription start site and 9 distinct promoter haplotypes were inferred from a panel of individuals from 5 distinct continental population groups. Following experimental validation, we found that the 2 major haplotypes significantly influenced TGFB1 transcriptional activity in an allele‐specific manner and that 3 of the SNPs (−1886G>A, −509C>T and −1550DEL/AGG) altered DNA‐protein complexe formation. Though the biological relevance of these findings remains to be verified, our study suggests that polymorphisms in the TGFB1 promoter could indeed influence gene expression and potentially contribute to the pathogenesis of TGFB1 related diseases.

Phosphodiesterase Type 4D Gene Polymorphism: Association with the Response to Short-Acting Bronchodilators in Paediatric Asthma Patients

Short-acting b2-adrenergic receptor agonists are commonly used bronchodilators for symptom relief in asthmatics. The aim of this study was to test whether genetic variants in PDE4D gene, a key regulator of b2-adrenoceptor-induced cAMP turnover in airway smooth muscle cells, affect the response to short-acting b2-agonists. Bronchodilator responsiveness was assessed in 133 asthmatic children by % change in baseline forced expiratory volume in one second (FEV1) after administration of albuterol. The analyses were performed in patients with airway obstruction (FEV1/FVC ratio below 90%, n = 93). FEV1  % change adjusted for baseline FEV1 values was significantly different between genotypes of rs1544791 G/A polymorphism (P = 0.006) and −1345 C/T (rs1504982) promoter variation (P = 0.03). The association remained significant with inclusion of age, sex, atopy, and controller medication into multivariate model (P = 0.004 and P = 0.02, resp.). Our work identifies new genetic variants implicated in modulation of asthma treatment, one of them (rs1544791) previously associated with asthma phenotype.

Genomic determinants of long-term cardiometabolic complications in childhood acute lymphoblastic leukemia survivors

BackgroundWhile cure rates for childhood acute lymphoblastic leukemia (cALL) now exceed 80%, over 60% of survivors will face treatment-related long-term sequelae, including cardiometabolic complications such as obesity, insulin resistance, dyslipidemia and hypertension. Although genetic susceptibility contributes to the development of these problems, there are very few studies that have so far addressed this issue in a cALL survivorship context.MethodsIn this study, we aimed at evaluating the associations between common and rare genetic variants and long-term cardiometabolic complications in survivors of cALL. We examined the cardiometabolic profile and performed whole-exome sequencing in 209 cALL survivors from the PETALE cohort. Variants associated with cardiometabolic outcomes were identified using PLINK (common) or SKAT (common and rare) and a logistic regression was used to evaluate their impact in multivariate models.ResultsOur results showed that rare and common variants in the BAD and FCRL3 genes were associated (p<0.05) with an extreme cardiometabolic phenotype (3 or more cardiometabolic risk factors). Common variants in OGFOD3 and APOB as well as rare and common BAD variants were significantly (p<0.05) associated with dyslipidemia. Common BAD and SERPINA6 variants were associated (p<0.05) with obesity and insulin resistance, respectively.ConclusionsIn summary, we identified genetic susceptibility loci as contributing factors to the development of late treatment-related cardiometabolic complications in cALL survivors. These biomarkers could be used as early detection strategies to identify susceptible individuals and implement appropriate measures and follow-up to prevent the development of risk factors in this high-risk population.

Connections Between ETV6-Modulated Genes: Identification of Shared Features

Accumulating genetic and functional evidence point to ETV6 as being the tumour suppressor gene targeted by the deletions at chromosome 12p12–13 found in various cancers, particularly childhood leukemia. ETV6 is a ubiquitously expressed transcription factor (TF) of the ETS family with very few known targeted genes. We recently compiled a list of 87 ETV6-modulated genes that can be classified into a number of subgroups based on their coordinated expression patterns. In the present report, we hypothesized that genes presenting a similar profile of modulation could also share biological features, promoter sequence similarities and/or, common transcription factor binding sites (TFBSs). Using an exploratory approach based on hierarchical clustering of expression data, Gene Ontology (GO) terms, sequence similarity and evolutionary conserved putative TFBSs, we found that many genes presenting a similar expression profile also share biological features and/or conserved predicted TFBSs but rarely show detectable promoter sequence similarities. We also calculated the proportion of ETV6-modulated genes that have any conserved TFBSs of the Jaspar database in their regulatory sequence and compared these proportions to those calculated for two other gene lists, ETV6 non-modulated and ETS-regulated. We found that the NF-kB, c-REL and p65 TFBSs, which all bind TFs of the REL class, were under-represented among the ETV6-modulated genes compared to the ETV6-non-modulated genes, while the Broad-complex 1 TFBS appeared to be over-represented. NF-Y and Chop/cEBP TFBSs were over-represented in the promoters of ETV6-modulated genes compared to ETS-regulated genes. These analyses will help direct further studies intending to understand the role of ETV6 as a transcriptional regulator and aid in constructing the ETV6-regulatory gene network.

KMT2E-ASNS: a novel relapse-specific fusion gene in early T-cell precursor acute lymphoblastic leukemia

To the editor: Early T-cell precursor acute lymphoblastic leukemia (ETP-ALL) is a recently characterized subtype accounting for up to 15% of childhood T-cell acute lymphoblastic leukemia (T-ALL).[1][1][⇓][2]-[3][3] This heterogeneous subgroup shows a distinctive immature immunophenotype