Identification of ERG Specific Target Genes by Genome-Wide Screening in T-Lymphoblastic Leukemia

Blood ◽  
2008 ◽  
Vol 112 (11) ◽  
pp. 3788-3788
Author(s):  
Liliana H Mochmann ◽  
Konrad Neumann ◽  
Juliane Bock ◽  
Jutta Ortiz Tanchez ◽  
Arend Bohne ◽  
...  

Abstract The Ets related gene, ERG, encodes a transcription factor with a vital role in hematopoiesis. Recent findings have shown that ERG knockout mice require a minimum of one functional allele to ensure embryonic blood development and adult stem cell maintenance. Moreover, it was earlier reported that enforced expression of ERG induced oncogenic transformation in 3T3 cells. Overexpression of ERG, observed in a subset of acute T-lymphoblastic and acute myeloid leukemia patients, was associated with an inferior outcome. However, the impact of ERG contributing to this unfavourable phenotype has yet to be determined, as downstream targets of ERG in leukemia remain unknown. Herein, we conducted a genome-wide analysis of ERG target genes in T-lymphoblastic leukemia. Chromatin immunoprecipitation-on-chip array (ChIP-on-chip) analyses were performed using two ERG specific antibodies for the enrichment of ERG-bound DNA templates in T-lymphoblastic leukemia cells (Jurkat) with input DNA or IgG precipitated DNA as controls. Enriched DNA templates and control DNA were differentially labelled and co-hybridized to high resolution promoter chip arrays with 50–75mer probes (770,000) representing 29,000 annotated human transcripts (NimbleGen). Based on two independent ChIP-on-chip assays, bioinformatic analysis (ACME) yielded statistically significant enriched peaks (using a sliding window of 1000 bp, and a P-value < 0.0001) identifying promoter regions of 365 potential ERG target genes. From these genes, clustering by functional annotation was performed using the DAVID database and subsequently genes related to leukemia were further selected for quantitative PCR validation. The design of promoter primers included the highly conserved ETS GGAA DNA binding site. Genes with greater than two-fold enrichment (ERG ChIP versus control) included WNT2 (17-fold), OLIG2 (14-fold), WNT11 (7-fold), CCND1 (5-fold), WNT9A (4-fold), CD7 (3-fold), EPO (3-fold), ERBB4 (3-fold), RPBJL (3-fold), TRADD (3-fold), PIWIL1 (2-fold), TNFRSF25 (2-fold), TWIST1 (2-fold), and HDAC4 (2-fold). Interestingly, enriched target genes involved in developmental processes (WNT2, WNT9A, WNT11, TWIST1, PIWIL1, ERBB4, and OLIG2) have shown oncogenic potential when mutated or overexpressed. Thus, we hypothesize that overexpression of ERG may contribute to T-cell leukemogenesis by the deregulation of these oncogenic targets. Further disclosure of ERG directed downstream pathways may contribute to the design of specific treatment strategies (such as WNT inhibitors) with particular effectiveness in ERG deregulated leukemia.

Blood ◽  
2021 ◽  
Vol 138 (Supplement 1) ◽  
pp. 3306-3306
Author(s):  
Mabel M. Jung ◽  
Koichi R. Katsumura ◽  
Peng Liu ◽  
Kirby D. Johnson ◽  
Joel P. Mackay ◽  
...  

Abstract Genomic analyses in clinical and experimental contexts have accelerated discoveries of human genetic variants. While elucidating the consequences of conspicuously loss-of-function variants is highly tractable, decoding the impact of missense or non-coding variants is considerably more challenging. Previously, we described a germline variant in GATA2 in a patient with GATA2-deficiency syndrome, which inserts nine amino acids between the two zinc fingers (9aa-Ins), one of which mediates sequence-specific DNA binding (Cavalcante de Andrade Silva M. et al., Leukemia, 2021). Unlike other GATA2 coding region and enhancer variants identified (Bresnick E.H. et al., Blood Adv., 2020), it was unclear whether the altered zinc finger spacing would be inhibitory, stimulatory or of no consequence. The 9aa-Ins variant was defective in activating several target genes (Hdc, Ear2 and Tpsb2) in a genetic complementation assay with Gata2 -77 enhancer-mutant (-77 -/-) primary hematopoietic progenitor cells. As only several target genes were tested, we used RNA-seq to conduct a genome-wide comparison of the capacity of GATA2 and 9aa-Ins to activate and repress transcription. To elucidate mechanisms, we considered the following models: 1) 9aa-Ins fails to regulate all genes normally controlled by GATA2; 2) 9aa-Ins fails to repress all genes normally controlled by GATA2; 3) 9aa-Ins fails to activate genes normally controlled by GATA2; 4) 9aa-Ins ectopically regulates genes not controlled by GATA2. Using a genetic complementation approach with -77 -/- cells that were immortalized by the Hoxb8 transcription factor (hi-77 -/-) (Wang G.G. et al., Nat. Methods, 2006; Johnson K.D. et al., JEM, 2020), we compared GATA2 and 9aa-Ins activities when expressed at a comparable level. This analysis revealed 2,138 GATA2-regulated, 525 GATA2 and 9aa-Ins-regulated, and 414 ectopically-regulated genes (at least two-fold change, adjusted P-value &lt;0.05). A similar number of genes were GATA2-activated (1,061) and repressed (1,077). Only 144 out of the 1,061 (14%) were 9aa-Ins-activated and 381 out of 1,077 (35%) were 9aa-Ins-repressed, illustrating the severe consequences of this mutation and a greater impact on activation versus repression. Statistical analysis with a range of P-values constraints (0.01 to 0.1) verified that activation by 9aa-Ins was more significantly impaired than repression (86% were no longer activated, and 65% were no longer repressed, P = 5.4 x 10 -6). Gene ontology analysis revealed that the 9aa insertion impaired GATA2-mediated activation of genes related to GPCR signaling and GATA2-mediated repression of genes related to innate immune machinery. The ectopically-regulated genes did not conform to a particular mechanism or pathway. Since it was unclear whether the transcriptional defects of the 9aa-Ins mutant reflect a unique inhibitory activity imparted by the 9aa sequence, we systematically varied the inter-zinc finger spacer length to establish whether any alterations can be tolerated. Using the genetic complementation assay, 2, 4, 6, 8, and 9aa spacer variants were compared with GATA2 for their capacity to regulate GATA2-activated genes (Hdc, Il1rl1, Gata1 and S1pr1) and repressed genes (Irf8, Il6st, Il6ra and Tifab). GATA2-mediated activation was compromised by insertions of two amino acids or more, whereas repression tolerated two and four amino acid insertions; 6, 8 and 9 were more inhibitory. Quantitative analyses revealed that a 6aa insertion reduced activation of the GATA2-activated genes by &gt;50% of the wild type value, whereas the GATA2-repressed genes were still repressed by at least 50% (18% retention of activation and 83% retention of repression, P = 0.001). Thus, zinc finger spacing alterations differentially impacted activation versus repression. These results provide a rigorous foundation for interpreting variants that alter zinc finger spacing without disrupting vital finger residues. In vitro and in vivo functional analyses and molecular modeling are ongoing to further dissect the underlying mechanisms and ascertain the importance of genetic networks and circuits that are sensitive or resistant to human disease variants. Disclosures No relevant conflicts of interest to declare.


2021 ◽  
Vol 9 (8) ◽  
pp. 1570
Author(s):  
Chien-Hsun Huang ◽  
Chih-Chieh Chen ◽  
Yu-Chun Lin ◽  
Chia-Hsuan Chen ◽  
Ai-Yun Lee ◽  
...  

The current taxonomy of the Lactiplantibacillus plantarum group comprises of 17 closely related species that are indistinguishable from each other by using commonly used 16S rRNA gene sequencing. In this study, a whole-genome-based analysis was carried out for exploring the highly distinguished target genes whose interspecific sequence identity is significantly less than those of 16S rRNA or conventional housekeeping genes. In silico analyses of 774 core genes by the cano-wgMLST_BacCompare analytics platform indicated that csbB, morA, murI, mutL, ntpJ, rutB, trmK, ydaF, and yhhX genes were the most promising candidates. Subsequently, the mutL gene was selected, and the discrimination power was further evaluated using Sanger sequencing. Among the type strains, mutL exhibited a clearly superior sequence identity (61.6–85.6%; average: 66.6%) to the 16S rRNA gene (96.7–100%; average: 98.4%) and the conventional phylogenetic marker genes (e.g., dnaJ, dnaK, pheS, recA, and rpoA), respectively, which could be used to separat tested strains into various species clusters. Consequently, species-specific primers were developed for fast and accurate identification of L. pentosus, L. argentoratensis, L. plantarum, and L. paraplantarum. During this study, one strain (BCRC 06B0048, L. pentosus) exhibited not only relatively low mutL sequence identities (97.0%) but also a low digital DNA–DNA hybridization value (78.1%) with the type strain DSM 20314T, signifying that it exhibits potential for reclassification as a novel subspecies. Our data demonstrate that mutL can be a genome-wide target for identifying and classifying the L. plantarum group species and for differentiating novel taxa from known species.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Morteza Bitaraf Sani ◽  
Javad Zare Harofte ◽  
Mohammad Hossein Banabazi ◽  
Saeid Esmaeilkhanian ◽  
Ali Shafei Naderi ◽  
...  

AbstractFor thousands of years, camels have produced meat, milk, and fiber in harsh desert conditions. For a sustainable development to provide protein resources from desert areas, it is necessary to pay attention to genetic improvement in camel breeding. By using genotyping-by-sequencing (GBS) method we produced over 14,500 genome wide markers to conduct a genome- wide association study (GWAS) for investigating the birth weight, daily gain, and body weight of 96 dromedaries in the Iranian central desert. A total of 99 SNPs were associated with birth weight, daily gain, and body weight (p-value < 0.002). Genomic breeding values (GEBVs) were estimated with the BGLR package using (i) all 14,522 SNPs and (ii) the 99 SNPs by GWAS. Twenty-eight SNPs were associated with birth weight, daily gain, and body weight (p-value < 0.001). Annotation of the genomic region (s) within ± 100 kb of the associated SNPs facilitated prediction of 36 candidate genes. The accuracy of GEBVs was more than 0.65 based on all 14,522 SNPs, but the regression coefficients for birth weight, daily gain, and body weight were 0.39, 0.20, and 0.23, respectively. Because of low sample size, the GEBVs were predicted using the associated SNPs from GWAS. The accuracy of GEBVs based on the 99 associated SNPs was 0.62, 0.82, and 0.57 for birth weight, daily gain, and body weight. This report is the first GWAS using GBS on dromedary camels and identifies markers associated with growth traits that could help to plan breeding program to genetic improvement. Further researches using larger sample size and collaboration of the camel farmers and more profound understanding will permit verification of the associated SNPs identified in this project. The preliminary results of study show that genomic selection could be the appropriate way to genetic improvement of body weight in dromedary camels, which is challenging due to a long generation interval, seasonal reproduction, and lack of records and pedigrees.


Genetics ◽  
2003 ◽  
Vol 164 (1) ◽  
pp. 247-258 ◽  
Author(s):  
Jinghong Li ◽  
Willis X Li

Abstract Overactivation of receptor tyrosine kinases (RTKs) has been linked to tumorigenesis. To understand how a hyperactivated RTK functions differently from wild-type RTK, we conducted a genome-wide systematic survey for genes that are required for signaling by a gain-of-function mutant Drosophila RTK Torso (Tor). We screened chromosomal deficiencies for suppression of a gain-of-function mutation tor (torGOF), which led to the identification of 26 genomic regions that, when in half dosage, suppressed the defects caused by torGOF. Testing of candidate genes in these regions revealed many genes known to be involved in Tor signaling (such as those encoding the Ras-MAPK cassette, adaptor and structural molecules of RTK signaling, and downstream target genes of Tor), confirming the specificity of this genetic screen. Importantly, this screen also identified components of the TGFβ (Dpp) and JAK/STAT pathways as being required for TorGOF signaling. Specifically, we found that reducing the dosage of thickveins (tkv), Mothers against dpp (Mad), or STAT92E (aka marelle), respectively, suppressed torGOF phenotypes. Furthermore, we demonstrate that in torGOF embryos, dpp is ectopically expressed and thus may contribute to the patterning defects. These results demonstrate an essential requirement of noncanonical signaling pathways for a persistently activated RTK to cause pathological defects in an organism.


2021 ◽  
Vol 11 (1) ◽  
pp. 59
Author(s):  
Kirsten Voorhies ◽  
Joanne E. Sordillo ◽  
Michael McGeachie ◽  
Elizabeth Ampleford ◽  
Alberta L. Wang ◽  
...  

An unaddressed and important issue is the role age plays in modulating response to short acting β2-agonists in individuals with asthma. The objective of this study was to identify whether age modifies genetic associations of single nucleotide polymorphisms (SNPs) with bronchodilator response (BDR) to β2-agonists. Using three cohorts with a total of 892 subjects, we ran a genome wide interaction study (GWIS) for each cohort to examine SNP by age interactions with BDR. A fixed effect meta-analysis was used to combine the results. In order to determine if previously identified BDR SNPs had an age interaction, we also examined 16 polymorphisms in candidate genes from two published genome wide association studies (GWAS) of BDR. There were no significant SNP by age interactions on BDR using the genome wide significance level of 5 × 10−8. Using a suggestive significance level of 5 × 10−6, three interactions, including one for a SNP within PRAG1 (rs4840337), were significant and replicated at the significance level of 0.05. Considering candidate genes from two previous GWAS of BDR, three SNPs (rs10476900 (near ADRB2) [p-value = 0.009], rs10827492 (CREM) [p-value = 0.02], and rs72646209 (NCOA3) [p-value = 0.02]) had a marginally significant interaction with age on BDR (p < 0.05). Our results suggest age may be an important modifier of genetic associations for BDR in asthma.


Endocrinology ◽  
2018 ◽  
Vol 160 (1) ◽  
pp. 38-54 ◽  
Author(s):  
Keiichi Itoi ◽  
Ikuko Motoike ◽  
Ying Liu ◽  
Sam Clokie ◽  
Yasumasa Iwasaki ◽  
...  

Abstract Glucocorticoids (GCs) are essential for stress adaptation, acting centrally and in the periphery. Corticotropin-releasing factor (CRF), a major regulator of adrenal GC synthesis, is produced in the paraventricular nucleus of the hypothalamus (PVH), which contains multiple neuroendocrine and preautonomic neurons. GCs may be involved in diverse regulatory mechanisms in the PVH, but the target genes of GCs are largely unexplored except for the CRF gene (Crh), a well-known target for GC negative feedback. Using a genome-wide RNA-sequencing analysis, we identified transcripts that changed in response to either high-dose corticosterone (Cort) exposure for 12 days (12-day high Cort), corticoid deprivation for 7 days (7-day ADX), or acute Cort administration. Among others, canonical GC target genes were upregulated prominently by 12-day high Cort. Crh was upregulated or downregulated most prominently by either 7-day ADX or 12-day high Cort, emphasizing the recognized feedback effects of GC on the hypothalamic-pituitary-adrenal (HPA) axis. Concomitant changes in vasopressin and apelin receptor gene expression are likely to contribute to HPA repression. In keeping with the pleotropic cellular actions of GCs, 7-day ADX downregulated numerous genes of a broad functional spectrum. The transcriptome response signature differed markedly between acute Cort injection and 12-day high Cort. Remarkably, six immediate early genes were upregulated 1 hour after Cort injection, which was confirmed by quantitative reverse transcription PCR and semiquantitative in situ hybridization. This study may provide a useful database for studying the regulatory mechanisms of GC-dependent gene expression and repression in the PVH.


2017 ◽  
Vol 7 (1) ◽  
Author(s):  
Chunmei He ◽  
Jaime A. Teixeira da Silva ◽  
Jianwen Tan ◽  
Jianxia Zhang ◽  
Xiaoping Pan ◽  
...  

2016 ◽  
Author(s):  
Bethany Signal ◽  
Brian S Gloss ◽  
Marcel E Dinger ◽  
Timothy R Mercer

ABSTRACTBackgroundThe branchpoint element is required for the first lariat-forming reaction in splicing. However due to difficulty in experimentally mapping at a genome-wide scale, current catalogues are incomplete.ResultsWe have developed a machine-learning algorithm trained with empirical human branchpoint annotations to identify branchpoint elements from primary genome sequence alone. Using this approach, we can accurately locate branchpoints elements in 85% of introns in current gene annotations. Consistent with branchpoints as basal genetic elements, we find our annotation is unbiased towards gene type and expression levels. A major fraction of introns was found to encode multiple branchpoints raising the prospect that mutational redundancy is encoded in key genes. We also confirmed all deleterious branchpoint mutations annotated in clinical variant databases, and further identified thousands of clinical and common genetic variants with similar predicted effects.ConclusionsWe propose the broad annotation of branchpoints constitutes a valuable resource for further investigations into the genetic encoding of splicing patterns, and interpreting the impact of common- and disease-causing human genetic variation on gene splicing.


2021 ◽  
Author(s):  
Richard F Oppong ◽  
Pau Navarro ◽  
Chris S Haley ◽  
Sara Knott

We describe a genome-wide analytical approach, SNP and Haplotype Regional Heritability Mapping (SNHap-RHM), that provides regional estimates of the heritability across locally defined regions in the genome. This approach utilises relationship matrices that are based on sharing of SNP and haplotype alleles at local haplotype blocks delimited by recombination boundaries in the genome. We implemented the approach on simulated data and show that the haplotype-based regional GRMs capture variation that is complementary to that captured by SNP-based regional GRMs, and thus justifying the fitting of the two GRMs jointly in a single analysis (SNHap-RHM). SNHap-RHM captures regions in the genome contributing to the phenotypic variation that existing genome-wide analysis methods may fail to capture. We further demonstrate that there are real benefits to be gained from this approach by applying it to real data from about 20,000 individuals from the Generation Scotland: Scottish Family Health Study. We analysed height and major depressive disorder (MDD). We identified seven genomic regions that are genome-wide significant for height, and three regions significant at a suggestive threshold (p-value <1x10^(-5) ) for MDD. These significant regions have genes mapped to within 400kb of them. The genes mapped for height have been reported to be associated with height in humans, whiles those mapped for MDD have been reported to be associated with major depressive disorder and other psychiatry phenotypes. The results show that SNHap-RHM presents an exciting new opportunity to analyse complex traits by allowing the joint mapping of novel genomic regions tagged by either SNPs or haplotypes, potentially leading to the recovery of some of the "missing" heritability.


2017 ◽  
Author(s):  
Tomasz Dzida ◽  
Mudassar Iqbal ◽  
Iryna Charapitsa ◽  
George Reid ◽  
Henk Stunnenberg ◽  
...  

We have developed a machine learning approach to predict context specific enhancer-promoter interactions using evidence from changes in genomic protein occupancy over time. The occupancy of estrogen receptor alpha (ERα), RNA polymerase (Pol II) and histone marks H2AZ and H3K4me3 were measured over time using ChIP-Seq experiments in MCF7 cells stimulated with estrogen. A Bayesian classifier was developed which uses the correlation of temporal binding patterns at enhancers and promoters and genomic proximity as features to predict interactions. This method was trained using experimentally determined interactions from the same system and was shown to achieve much higher precision than predictions based on the genomic proximity of nearest ERα binding. We use the method to identify a genome-wide confident set of ERα target genes and their regulatory enhancers genome-wide. Validation with publicly available GRO-Seq data demonstrates that our predicted targets are much more likely to show early nascent transcription than predictions based on genomic ERα binding proximity alone.


Sign in / Sign up

Export Citation Format

Share Document