Integrative Meta-Analysis of Differential Gene Expression in Acute Myeloid Leukemia

BACKGROUND: Acute myeloid leukemia (AML) is a heterogeneous disease with an overall poor prognosis. Gene expression profiling studies of patients with AML has provided key insights into disease pathogenesis while exposing potential diagnostic and prognostic markers and therapeutic targets. A systematic comparison of the large body of gene expression profiling studies in AML has the potential to test the extensibility of conclusions based on single studies and provide further insights into AML. METHODOLOGY/PRINCIPAL FINDINGS: In this study, we systematically compared 25 published reports of gene expression profiling in AML. There were a total of 4,918 reported genes of which one third were reported in more than one study. We found that only a minority of reported prognostically-associated genes (9.6%) were replicated in at least one other study. In a combined analysis, we comprehensively identified both gene sets and functional gene categories and pathways that exhibited significant differential regulation in distinct prognostic categories, including many previously unreported associations. CONCLUSIONS/SIGNIFICANCE: We developed a novel approach for granular, cross-study analysis of gene-by-gene data and their relationships with established prognostic features and patient outcome. We identified many robust novel prognostic molecular features in AML that were undetected in prior studies, and which provide insights into AML pathogenesis with potential diagnostic, prognostic, and therapeutic implications. Our database and integrative analysis are available online (http://gat.stamlab.org)

The role of mutation rate variation and genetic diversity in the architecture of human disease

Background We have investigated the role that the mutation rate and the structure of genetic variation at a locus play in determining whether a gene is involved in disease. We predict that the mutation rate and its genetic diversity should be higher in genes associated with disease, unless all genes that could cause disease have already been identified. Results Consistent with our predictions we find that genes associated with Mendelian and complex disease are substantially longer than non-disease genes. However, we find that both Mendelian and complex disease genes are found in regions of the genome with relatively low mutation rates, as inferred from intron divergence between humans and chimpanzees, and they are predicted to have similar rates of non-synonymous mutation as other genes. Finally, we find that disease genes are in regions of significantly elevated genetic diversity, even when variation in the rate of mutation is controlled for. The effect is small nevertheless. Conclusions Our results suggest that gene length contributes to whether a gene is associated with disease. However, the mutation rate and the genetic architecture of the locus appear to play only a minor role in determining whether a gene is associated with disease

Integrative analysis of 3604 GWAS reveals multiple novel cell type-specific regulatory associations

BACKGROUND: Genome-wide association study (GWAS) single nucleotide polymorphisms (SNPs) are known to preferentially co-locate to active regulatory elements in tissues and cell types relevant to disease aetiology. Further characterisation of associated cell type-specific regulation can broaden our understanding of how GWAS signals may contribute to disease risk. RESULTS: To gain insight into potential functional mechanisms underlying GWAS associations, we developed FORGE2 ( https://forge2.altiusinstitute.org/ ), which is an updated version of the FORGE web tool. FORGE2 uses an expanded atlas of cell type-specific regulatory element annotations, including DNase I hotspots, five histone mark categories and 15 hidden Markov model (HMM) chromatin states, to identify tissue- and cell type-specific signals. An analysis of 3,604 GWAS from the NHGRI-EBI GWAS catalogue yielded at least one significant disease/trait-tissue association for 2,057 GWAS, including > 400 associations specific to epigenomic marks in immune tissues and cell types, > 30 associations specific to heart tissue, and > 60 associations specific to brain tissue, highlighting the key potential of tissue- and cell type-specific regulatory elements. Importantly, we demonstrate that FORGE2 analysis can separate previously observed accessible chromatin enrichments into different chromatin states, such as enhancers or active transcription start sites, providing a greater understanding of underlying regulatory mechanisms. Interestingly, tissue-specific enrichments for repressive chromatin states and histone marks were also detected, suggesting a role for tissue-specific repressed regions in GWAS-mediated disease aetiology. CONCLUSION: In summary, we demonstrate that FORGE2 has the potential to uncover previously unreported disease-tissue associations and identify new candidate mechanisms. FORGE2 is a transparent, user-friendly web tool for the integrative analysis of loci discovered from GWAS

Induction of endogenous γ-globin gene expression with decoy oligonucleotide targeting Oct-1 transcription factor consensus sequence

Human β-globin disorders are relatively common genetic diseases cause by mutations in the β-globin gene. Increasing the expression of the γ-globin gene has great benefits in reducing complications associated with these diseases. The Oct-1 transcription factor is involved in the transcriptional regulation of the γ-globin gene. The human γ-globin genes (both Aγ and Gγ-globin genes) carry three Oct-1 transcription factor consensus sequences within their promoter regions. We have studied the possibility of inducing γ-globin gene expression using decoy oligonucleotides that target the Oct-1 transcription factor consensus sequence. A double-stranded 22 bp decoy oligonucleotide containing the Oct-1 consensus sequence was synthesized. The results obtained from our in vitro binding assay revealed a strong competitive binding of the decoy oligonucleotide for the Oct-1 transcription factor. When K562 human erythroleukemia cells were treated with the Oct-1 decoy oligonucleotide, significant increases in the level of the γ-globin mRNA were observed. The results of our western blots further demonstrated significant increases of the fetal hemoglobin (HbF, α2γ2) in the Oct-1 decoy oligonucleotide-treated K562 cells. The results of our immunoprecipitation (IP) studies revealed that the treatment of K562 cells with the Oct-1 decoy oligonucleotide significantly reduced the level of the endogenous γ-globin gene promoter region DNA co-precipitated with the Oct-1 transcription factor. These results suggest that the decoy oligonucleotide designed for the Oct-1 transcription factor consensus sequence could induce expression of the endogenous γ-globin gene through competitive binding of the Oct-1 transcription factor, resulting in activation of the γ-globin genes. Therefore, disrupting the bindings of the Oct-1 transcriptional factors with the decoy oligonucleotide provides a novel approach for inducing expression of the γ-globin genes. It also provides an innovative strategy for the treatment of many disease conditions, including sickle cell anemia and β-thalassemia

Patterns of regulatory activity across diverse human cell types predict tissue identity, transcription factor binding, and long-range interactions

Regulatory elements recruit transcription factors that modulate gene expression distinctly across cell types, but the relationships among these remains elusive. To address this, we analyzed matched DNase-seq and gene expression data for 112 human samples representing 72 cell types. We first defined more than 1800 clusters of DNase I hypersensitive sites (DHSs) with similar tissue specificity of DNase-seq signal patterns. We then used these to uncover distinct associations between DHSs and promoters, CpG islands, conserved elements, and transcription factor motif enrichment. Motif analysis within clusters identified known and novel motifs in cell-type-specific and ubiquitous regulatory elements and supports a role for AP-1 regulating open chromatin. We developed a classifier that accurately predicts cell-type lineage based on only 43 DHSs and evaluated the tissue of origin for cancer cell types. A similar classifier identified three sex-specific loci on the X chromosome, including the XIST lincRNA locus. By correlating DNase I signal and gene expression, we predicted regulated genes for more than 500K DHSs. Finally, we introduce a web resource to enable researchers to use these results to explore these regulatory patterns and better understand how expression is modulated within and across human cell types

Genetic, environmental and stochastic factors in monozygotic twin discordance with a focus on epigenetic differences

PMCID: PMC3566971This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited

The role of transcriptional activator GATA-1 at human β-globin HS2

GATA-1 is an erythroid activator that binds β-globin gene promoters and DNase I hypersensitive sites (HSs) of the β-globin locus control region (LCR). We investigated the direct role of GATA-1 interaction at the LCR HS2 enhancer by mutating its binding sites within minichromosomes in erythroid cells. Loss of GATA-1 in HS2 did not compromise interaction of NF-E2, a second activator that binds to HS2, nor was DNase I hypersensitivity at HS2 or the promoter of a linked ε-globin gene altered. Reduction of NF-E2 using RNAi confirmed the overall importance of this activator in establishing LCR HSs. However, recruitment of the histone acetyltransferase CBP and RNA pol II to HS2 was diminished by GATA-1 loss. Transcription of ε-globin was severely compromised with loss of RNA pol II from the transcription start site and reduction of H3 acetylation and H3K4 di- and tri-methylation in coding sequences. In contrast, widespread detection of H3K4 mono-methylation was unaffected by loss of GATA-1 in HS2. These results support the idea that GATA-1 interaction in HS2 has a prominent and direct role in co-activator and pol II recruitment conferring active histone tail modifications and transcription activation to a target gene but that it does not, by itself, play a major role in establishing DNase I hypersensitivity

eFORGE v2.0: updated analysis of cell type-specific signal in epigenomic data

SUMMARY: The Illumina Infinium EPIC BeadChip is a new high-throughput array for DNA methylation analysis, extending the earlier 450k array by over 400,000 new sites. Previously, a method named eFORGE was developed to provide insights into cell type-specific and cell composition effects for 450k data. Here, we present a significantly updated and improved version of eFORGE that can analyse both EPIC and 450k array data. New features include analysis of chromatin states, TF motifs and DNase I footprints, providing tools for EWAS interpretation and epigenome editing. AVAILABILITY: eFORGE v2.0 is implemented as a web tool available from https://eforge.altiusinstitute.org and https://eforge-tf.altiusinstitute.org/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online

Defining functional DNA elements in the human genome

With the completion of the human genome sequence, attention turned to identifying and annotating its functional DNA elements. As a complement to genetic and comparative genomics approaches, the Encyclopedia of DNA Elements Project was launched to contribute maps of RNA transcripts, transcriptional regulator binding sites, and chromatin states in many cell types. The resulting genome-wide data reveal sites of biochemical activity with high positional resolution and cell type specificity that facilitate studies of gene regulation and interpretation of noncoding variants associated with human disease. However, the biochemically active regions cover a much larger fraction of the genome than do evolutionarily conserved regions, raising the question of whether nonconserved but biochemically active regions are truly functional. Here, we review the strengths and limitations of biochemical, evolutionary, and genetic approaches for defining functional DNA segments, potential sources for the observed differences in estimated genomic coverage, and the biological implications of these discrepancies. We also analyze the relationship between signal intensity, genomic coverage, and evolutionary conservation. Our results reinforce the principle that each approach provides complementary information and that we need to use combinations of all three to elucidate genome function in human biology and disease