Ru-Fang Yeh

miR-126 Regulates Angiogenic Signaling and Vascular Integrity

Precise regulation of the formation, maintenance, and remodeling of the vasculature is required for normal development, tissue response to injury, and tumor progression. How specific microRNAs intersect with and modulate angiogenic signaling cascades is unknown. Here, we identified microRNAs that were enriched in endothelial cells derived from mouse embryonic stem (ES) cells and in developing mouse embryos. We found that miR-126 regulated the response of endothelial cells to VEGF. Additionally, knockdown of miR-126 in zebrafish resulted in loss of vascular integrity and hemorrhage during embryonic development. miR-126 functioned in part by directly repressing negative regulators of the VEGF pathway, including the Sprouty-related protein SPRED1 and phosphoinositol-3 kinase regulatory subunit 2 (PIK3R2/p85-beta). Increased expression of Spred1 or inhibition of VEGF signaling in zebrafish resulted in defects similar to miR-126 knockdown. These findings illustrate that a single miRNA can regulate vascular integrity and angiogenesis, providing a new target for modulating vascular formation and function.

MicroRNA Regulation of Cell Lineages in Mouse and Human Embryonic Stem Cells

Cell fate decisions of pluripotent embryonic stem (ES) cells are dictated by activation and repression of lineage-specific genes. Numerous signaling and transcriptional networks progressively narrow and specify the potential of ES cells. Whether specific microRNAs help refine and limit gene expression and, thereby, could be used to manipulate ES cell differentiation has largely been unexplored. Here, we show that two serum response factor (SRF)-dependent muscle-specific microRNAs, miR-1 and miR-133, promote mesoderm formation from ES cells but have opposing functions during further differentiation into cardiac muscle progenitors. Furthermore, miR-1 and miR-133 were potent repressors of nonmuscle gene expression and cell fate during mouse and human ES cell differentiation. miR-1s effects were in part mediated by translational repression of the Notch ligand Delta-like 1 (Dll-1). Our findings indicate that muscle-specific miRNAs reinforce the silencing of nonmuscle genes during cell lineage commitment and suggest that miRNAs may have general utility in regulating cell-fate decisions from pluripotent ES cells.

Model-based clustering of DNA methylation array data: a recursive-partitioning algorithm for high-dimensional data arising as a mixture of beta distributions.

BackgroundEpigenetics is the study of heritable changes in gene function that cannot be explained by changes in DNA sequence. One of the most commonly studied epigenetic alterations is cytosine methylation, which is a well recognized mechanism of epigenetic gene silencing and often occurs at tumor suppressor gene loci in human cancer. Arrays are now being used to study DNA methylation at a large number of loci; for example, the Illumina GoldenGate platform assesses DNA methylation at 1505 loci associated with over 800 cancer-related genes. Model-based cluster analysis is often used to identify DNA methylation subgroups in data, but it is unclear how to cluster DNA methylation data from arrays in a scalable and reliable manner.ResultsWe propose a novel model-based recursive-partitioning algorithm to navigate clusters in a beta mixture model. We present simulations that show that the method is more reliable than competing nonparametric clustering approaches, and is at least as reliable as conventional mixture model methods. We also show that our proposed method is more computationally efficient than conventional mixture model approaches. We demonstrate our method on the normal tissue samples and show that the clusters are associated with tissue type as well as age.ConclusionOur proposed recursively-partitioned mixture model is an effective and computationally efficient method for clustering DNA methylation data.

Chromatin remodelling complex dosage modulates transcription factor function in heart development

Dominant mutations in cardiac transcription factor genes cause human inherited congenital heart defects (CHDs); however, their molecular basis is not understood. Interactions between transcription factors and the Brg1/Brm-associated factor (BAF) chromatin remodelling complex suggest potential mechanisms; however, the role of BAF complexes in cardiogenesis is not known. In this study, we show that dosage of Brg1 is critical for mouse and zebrafish cardiogenesis. Disrupting the balance between Brg1 and disease-causing cardiac transcription factors, including Tbx5, Tbx20 and Nkx2–5, causes severe cardiac anomalies, revealing an essential allelic balance between Brg1 and these cardiac transcription factor genes. This suggests that the relative levels of transcription factors and BAF complexes are important for heart development, which is supported by reduced occupancy of Brg1 at cardiac gene promoters in Tbx5 haploinsufficient hearts. Our results reveal complex dosage-sensitive interdependence between transcription factors and BAF complexes, providing a potential mechanism underlying transcription factor haploinsufficiency, with implications for multigenic inheritance of CHDs.

Epigenetic profiles distinguish pleural mesothelioma from normal pleura and predict lung asbestos burden and clinical outcome

Mechanisms of action of nonmutagenic carcinogens such as asbestos remain poorly characterized. As pleural mesothelioma is known to have limited numbers of genetic mutations, we aimed to characterize the relationships among gene-locus-specific methylation alterations, disease status, asbestos burden, and survival in this rapidly fatal asbestos-associated tumor. Methylation of 1505 CpG loci associated with 803 cancer-related genes were studied in 158 pleural mesotheliomas and 18 normal pleura. After false-discovery rate correction, 969 CpG loci were independently associated with disease status (Q < 0.05). Classifying samples based on CpG methylation profile with a mixture model approach, methylation classes discriminated tumor from normal pleura (permutation P < 0.0001). In a random forests classification, the overall misclassification error rate was 3.4%, with <1% (n = 1) of tumors misclassified as normal (P < 0.0001). Among tumors, methylation class membership was significantly associated with lung tissue asbestos body burden (P < 0.03), and significantly predicted survival (likelihood ratio P < 0.01). Consistent with prior work, asbestos burden was associated with an increased risk of death (hazard ratio, 1.4; 95% confidence interval, 1.1-1.8). Our results have shown that methylation profiles powerfully differentiate diseased pleura from nontumor pleura and that asbestos burden and methylation profiles are independent predictors of mesothelioma patient survival. We have added to the growing body of evidence that cellular epigenetic dysregulation is a critical mode of action for asbestos in the induction of pleural mesothelioma. Importantly, these findings hold great promise for using epigenetic profiling in the diagnosis and prognosis of human cancers.

Differentially expressed genes are marked by histone 3 lysine 9 trimethylation in human cancer cells

Histone H3 lysine 9 trimethylation (H3K9Me3) has been associated with transcriptional repression, but recent findings implicate this chromatin modification in transcriptional activation and mRNA elongation by RNA polymerase II. Here, we applied immunoprecipitation (IP) with a custom DNA tiling microarray containing many transcription factors important in development and cancer (for example homeotic genes; N=683 total genes) to explore the relationship between H3K9Me3 and other histone modifications with the differential expression of genes. Cancer cell lines derived from different tissues (2 leukemia, 2 medulloblastoma) were characterized with IP antibodies to H3K9Me3, H3K4 dimethylation (H3K4Me2) and H3K9 acetylation (H3K9Ac). MV4-11 is known to overexpress the HOXA9 and MEIS1 genes, whereas D283 overexpresses the OTX2 homeobox gene. Gene expression was assessed by Affymetrix U133 array. Mapping the number and size of histone markings demonstrated significant colocalization of H3K9Ac and H3K4Me2 with H3K9Me3, indicating a pattern of putative ‘activating’ and ‘repressive’ markings. The median site size was 600–821 bp and 72–95% or 53–80% of chromatin signal sites were located within 1 kb or 500 bp of transcription start sites (TSS), respectively. A relatively small number of genes displayed additional H3K9Me3 sites in the 5′-region distant from the TSS. Comparing genes with modification sites to those without sites in their promoters confirmed the positive associations of H3K9Ac and H3K4Me2 with gene expression and revealed that H3K9Me3 is associated with active genes rather than being a repressive marking as previously thought. The positive regulatory effect of all three types of modifications were quantitatively correlated with site size, and applied to absolute gene expression within a single cell line as well as relative expression among pairs of cell lines. Extended patterns of H3K9Me3 upstream of some genes (for example HOXA9 and OTX2) may result from the action of multiple promoter elements. We found an inverse relationship between promoter DNA hypermethylation and H3K9Me3 in three studied genes (HOXA9, TMS1, RASSF1A). The localization of H3K9Me3 downstream of the TSSs of expressed genes and not within promoter regions of hypermethylated and silenced genes is consistent with the proposed coupling of H3K9Me3 with RNA polymerase II. Our results indicate a need for revising aspects of the histone code involving H3 lysine methylation. Awareness of H3K9Me3 as a mark of gene activity, not repression, is especially important for the classification of human cancer using chromatin and histone profiles.

Transcriptional Analysis of Pluripotency Reveals the Hippo Pathway as a Barrier to Reprogramming

Pluripotent stem cells are derived from culture of early embryos or the germline and can be induced by reprogramming of somatic cells. Barriers to reprogramming that stabilize the differentiated state and have tumor suppression functions are expected to exist. However, we have a limited understanding of what such barriers might be. To find novel barriers to reprogramming to pluripotency, we compared the transcriptional profiles of the mouse germline with pluripotent and somatic cells, in vivo and in vitro. There is a remarkable global expression of the transcriptional program for pluripotency in primordial germ cells (PGCs). We identify parallels between PGC reprogramming to pluripotency and human germ cell tumorigenesis, including the loss of LATS2, a tumor suppressor kinase of the Hippo pathway. We show that knockdown of LATS2 increases the efficiency of induction of pluripotency in human cells. LATS2 RNAi, unlike p53 RNAi, specifically enhances the generation of fully reprogrammed iPS cells without accelerating cell proliferation. We further show that LATS2 represses reprogramming in human cells by post-transcriptionally antagonizing TAZ but not YAP, two downstream effectors of the Hippo pathway. These results reveal transcriptional parallels between germ cell transformation and the generation of iPS cells and indicate that the Hippo pathway constitutes a barrier to cellular reprogramming.

Pathway Analysis of Single-Nucleotide Polymorphisms Potentially Associated with Glioblastoma Multiforme Susceptibility Using Random Forests

Glioma is a complex disease that is unlikely to result from the effect of a single gene. Genetic analysis at the pathway level involving multiple genes may be more likely to capture gene-disease associations than analyzing genes one at a time. The current pilot study included 112 Caucasians with glioblastoma multiforme and 112 Caucasian healthy controls frequency matched to cases by age and gender. Subjects were genotyped using a commercially available (ParAllele/Affymetrix) assay panel of 10,177 nonsynonymous coding single-nucleotide polymorphisms (SNP) spanning the genome known at the time the panel was constructed. For this analysis, we selected 10 pathways potentially involved in gliomagenesis that had SNPs represented on the panel. We performed random forests (RF) analyses of SNPs within each pathway group and logistic regression to assess interaction among genes in the one pathway for which the RF prediction error was better than chance and the permutation P < 0.10. Only the DNA repair pathway had a better than chance classification of case-control status with a prediction error of 45.5% and P = 0.09. Three SNPs (rs1047840 of EXO1, rs12450550 of EME1, and rs799917 of BRCA1) of the DNA repair pathway were identified as promising candidates for further replication. In addition, statistically significant interactions (P < 0.05) between rs1047840 of EXO1 and rs799917 or rs1799966 of BRCA1 were observed. Despite less than complete inclusion of genes and SNPs relevant to glioma and a small sample size, RF analysis identified one important biological pathway and several SNPs potentially associated with the development of glioblastoma. (Cancer Epidemiol Biomarkers Prev 2008;17(6):1368–73)

TEL-AML1 regulation of survivin and apoptosis via miRNA-494 and miRNA-320a

There is increasing evidence that miRNA and transcription factors interact in an instructive fashion in normal and malignant hematopoiesis. We explored the impact of TEL-AML1 (ETV6-RUNX1), the most common fusion protein in childhood leukemia, on miRNA expression and the leukemic phenotype. Using RNA interference, miRNA expression arrays, and quantitative polymerase chain reaction, we identified miRNA-494 and miRNA-320a to be up-regulated upon TEL-AML1 silencing independently of TEL expression. Chromatin immunoprecipitation analysis identified miRNA-494 as a direct miRNA target of the fusion protein TEL-AML1. Using bioinformatic analysis as well as functional luciferase experiments, we demonstrate that survivin is a target of the 2 miRNAs. miRNA-494 and miRNA-320a were introduced to the cells by transfection and survivin expression determined by Western blot analysis. These miRNAs blocked survivin expression and resulted in apoptosis in a similar manner as TEL-AML1 silencing by itself; this silencing was also shown to be Dicer-dependent. miRNAs-494 and -320a are expressed at lower levels in TEL-AML1+ leukemias compared with immunophenotype-matched nonTEL-AML1 acute lymphoblastic leukemia subtypes, and within TEL-AML1+ leukemias their expression is correlated to survivin levels. In summary our data suggest that TEL-AML1 might exert its antiapoptotic action at least in part by suppressing miRNA-494 and miRNA-320a, lowering their expression causing enhanced survivin expression.

SNPLogic: an interactive single nucleotide polymorphism selection, annotation, and prioritization system

SNPLogic (http://www.snplogic.org) brings together single nucleotide polymorphism (SNP) information from numerous sources to provide a comprehensive SNP selection, annotation and prioritization system for design and analysis of genotyping projects. SNPLogic integrates information about the genetic context of SNPs (gene, chromosomal region, functional location, haplotypes tags and overlap with transcription factor binding sites, splicing sites, miRNAs and evolutionarily conserved regions), genotypic data (allele frequencies per population and validation method), coverage of commercial arrays (ParAllele, Affymetrix and Illumina), functional predictions (modeled on structure and sequence) and connections or established associations (biological pathways, gene ontology terms and OMIM disease terms). The SNPLogic web interface facilitates construction and annotation of user-defined SNP lists that can be saved, shared and exported. Thus, SNPLogic can be used to identify and prioritize candidate SNPs, assess custom and commercial arrays panels and annotate new SNP data with publicly available information. We have found integration of SNP annotation in the context of pathway information and functional prediction scores to be a powerful approach to the analysis and interpretation of SNP-disease association data.