Cited by 255Pluripotent Stem Cells ResearchBioinformatics and Genomic NetworksGene Regulatory Network Analysis

The DOI matches its Crossref record; 76 of the 77 references it cites were found in a registry or an index.

Signed weighted gene co-expression network analysis of transcriptional regulation in murine embryonic stem cells

Open access
Mason, Mike J, Fan, Guoping, Plath, Kathrin, Zhou, Qing, Horvath, Steve·2009·BMC Genomics

Abstract

Background: Recent work has revealed that a core group of transcription factors (TFs) regulates the key characteristics of embryonic stem (ES) cells: pluripotency and self-renewal. Current efforts focus on identifying genes that play important roles in maintaining pluripotency and self-renewal in ES cells and aim to understand the interactions among these genes. To that end, we investigated the use of unsigned and signed network analysis to identify pluripotency and differentiation related genes. Results: We show that signed networks provide a better systems level understanding of the regulatory mechanisms of ES cells than unsigned networks, using two independent murine ES cell expression data sets. Specifically, using signed weighted gene co-expression network analysis (WGCNA), we found a pluripotency module and a differentiation module, which are not identified in unsigned networks. We confirmed the importance of these modules by incorporating genome-wide TF binding data for key ES cell regulators. Interestingly, we find that the pluripotency module is enriched with genes related to DNA damage repair and mitochondrial function in addition to transcriptional regulation. Using a connectivity measure of module membership, we not only identify known regulators of ES cells but also show that Mrpl15, Msh6, Nrf1, Nup133, Ppif, Rbpj, Sh3gl2, and Zfp39, among other genes, have important roles in maintaining ES cell pluripotency and self-renewal. We also report highly significant relationships between module membership and epigenetic modifications (histone modifications and promoter CpG methylation status), which are known to play a role in controlling gene expression during ES cell self-renewal and differentiation. Conclusion: Our systems biologic re-analysis of gene expression, transcription factor binding, epigenetic and gene ontology data provides a novel integrative view of ES cell biology.

Overview

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Integrity signals

Registry

Crossref · 10.1186/1471-2164-10-327

Cites

76 of 76 references verified0 citation errors · 0 not checkable

Statistics

Figures

Figure 1 Network Connection Strength Versus Expression Correlation. Network adjacency (y-axis) versus correlation (x-axis) for an unweighted network (black step function with τ = 0.8) and weighted networks (dashed lines corresponding to different powers, β) in an unsigned network (a) and a signed network (b). Note that cor(xi, xj) = -1 leads to adjacency = 0 in the signed network. The weighted network preserves the continuous nature of the co-expression information while an unweighted network dichotomizes the correlation.
Figure 1 Network Connection Strength Versus Expression Correlation. Network adjacency (y-axis) versus correlation (x-axis) for an unweighted network (black step function with τ = 0.8) and weighted networks (dashed lines corresponding to different powers, β) in an unsigned network (a) and a signed network (b). Note that cor(xi, xj) = -1 leads to adjacency = 0 in the signed network. The weighted network preserves the continuous nature of the co-expression information while an unweighted network dichotomizes the correlation.CC BY · PubMed Central
Figure 2 Unsigned and Signed Mouse ES cell Networks in Ivanova et al. a, left, Dendrogram of the unsigned network of the Ivanova et al (2006) data set with color bands below indicating module membership for the unsigned network (U) and the signed network (S). a, right, heat map for visualizing standardized gene expressions (rows) across samples (columns) for genes in the turquoise module in the unsigned network. b, left, Dendrogram of the signed network of the Ivanova et al (2006) data set with color bands below indicating module membership for the signed network (S) and the unsigned network (U). b, right, heat map of expression profiles across samples for genes in the turquoise, black, and blue modules in the signed network. Note, modules are not scaled to reflect the number of genes in each module. c, scatter plot of module membership, kME, (x-axis) plotted against gene significance, GS, (y-axis) for the black and blue modules in the signed network with known ES cell regulators and differentiation genes labelled.
Figure 2 Unsigned and Signed Mouse ES cell Networks in Ivanova et al. a, left, Dendrogram of the unsigned network of the Ivanova et al (2006) data set with color bands below indicating module membership for the unsigned network (U) and the signed network (S). a, right, heat map for visualizing standardized gene expressions (rows) across samples (columns) for genes in the turquoise module in the unsigned network. b, left, Dendrogram of the signed network of the Ivanova et al (2006) data set with color bands below indicating module membership for the signed network (S) and the unsigned network (U). b, right, heat map of expression profiles across samples for genes in the turquoise, black, and blue modules in the signed network. Note, modules are not scaled to reflect the number of genes in each module. c, scatter plot of module membership, kME, (x-axis) plotted against gene significance, GS, (y-axis) for the black and blue modules in the signed network with known ES cell regulators and differentiation genes labelled.CC BY · PubMed Central
Figure 3 Relating Module Membership to Epigenetic Regulation. The top 1000 genes with highest module membership in the black module (top row) and blue module (bottom row) are related to 3 epigenetic variables (corresponding to the 3 columns). The y-axis reports the proportion of top 1000 genes that are known to belong to the group of genes defined on the x-axis. Histone H3K4me3 trimethylation status is abbreviated K4, H3K27me3 trimethylation statys is abbreviated by K27. Promoters that are both H3K4 and H3K27 trimethilated in ES cells (denoted K4&K27) are thought to poise key developmental genes for activation upon differentiation [50,51]. Note that genes with promoter CpG methylation are significantly (p = 2.0 × 10-14) under-enriched with respect to the top 1000 black module genes.
Figure 3 Relating Module Membership to Epigenetic Regulation. The top 1000 genes with highest module membership in the black module (top row) and blue module (bottom row) are related to 3 epigenetic variables (corresponding to the 3 columns). The y-axis reports the proportion of top 1000 genes that are known to belong to the group of genes defined on the x-axis. Histone H3K4me3 trimethylation status is abbreviated K4, H3K27me3 trimethylation statys is abbreviated by K27. Promoters that are both H3K4 and H3K27 trimethilated in ES cells (denoted K4&K27) are thought to poise key developmental genes for activation upon differentiation [50,51]. Note that genes with promoter CpG methylation are significantly (p = 2.0 × 10-14) under-enriched with respect to the top 1000 black module genes.CC BY · PubMed Central
Figure 4 Unsigned and Signed Networks of the Zhou et al ES Expression Data. a, left, Dendrogram of the unsigned network of the Zhou et al (2007) data with color bands below indicating module membership for the unsigned network (U) and the signed network (S). a, right, A heat map shows microarray expression profiles accross samples for genes in the blue and black modules in the unsigned network. b, left, Dendrogram of the signed network of the Zhou et al (2007) data with color bands below indicating module membership for the signed network (S) and the unsigned network (U). b, right, heatmap of expression profiles (rows) across samples (columns) for genes in the blue and black modules in the signed network.
Figure 4 Unsigned and Signed Networks of the Zhou et al ES Expression Data. a, left, Dendrogram of the unsigned network of the Zhou et al (2007) data with color bands below indicating module membership for the unsigned network (U) and the signed network (S). a, right, A heat map shows microarray expression profiles accross samples for genes in the blue and black modules in the unsigned network. b, left, Dendrogram of the signed network of the Zhou et al (2007) data with color bands below indicating module membership for the signed network (S) and the unsigned network (U). b, right, heatmap of expression profiles (rows) across samples (columns) for genes in the blue and black modules in the signed network.CC BY · PubMed Central
Figure 5 Expression Changes Versus Module Membership in the Black and Blue Modules (Zhou et al). Module membership, kME, is plotted against log2 expression fold change (FC) for the black (a) and blue (b) modules in the unsigned network of the Zhou et al (2007) data. FC is the ratio between the average expression in Oct4 positive samples and Oct4 negative microarray samples. Known ES cell regulators are labeled. Genes are colored by module membership in the signed network. c and d are analogous to a and b but module membership is with regard to the signed black and blue modules.
Figure 5 Expression Changes Versus Module Membership in the Black and Blue Modules (Zhou et al). Module membership, kME, is plotted against log2 expression fold change (FC) for the black (a) and blue (b) modules in the unsigned network of the Zhou et al (2007) data. FC is the ratio between the average expression in Oct4 positive samples and Oct4 negative microarray samples. Known ES cell regulators are labeled. Genes are colored by module membership in the signed network. c and d are analogous to a and b but module membership is with regard to the signed black and blue modules.CC BY · PubMed Central
Figure 6 Comparison of Genes Ranked by Network Connectivity and Differential Expression in the Ivanova et al data set. Ingenuity Pathway Analysis of functional enrichments in the set of genes ranked within the top 1000 by Student's t-test and kME and yet do not overlap with each other. Venn diagrams show the amount of gene overlap between the top 1000 black (pluripotency) module genes and the top 1000 genes most significantly down-regulated upon Oct4 RNAi (left); gene overlap between the top 1000 blue (differentiation) module genes and the 1000 genes most significantly up-regulated with Oct4 RNAi (right). Significance of differential expression was determined using Student's t-statistic. p-values have been corrected for multiple hypothesis tests (Benjamini-Hochberg). Only significantly enriched functional groups are shown.
Figure 6 Comparison of Genes Ranked by Network Connectivity and Differential Expression in the Ivanova et al data set. Ingenuity Pathway Analysis of functional enrichments in the set of genes ranked within the top 1000 by Student's t-test and kME and yet do not overlap with each other. Venn diagrams show the amount of gene overlap between the top 1000 black (pluripotency) module genes and the top 1000 genes most significantly down-regulated upon Oct4 RNAi (left); gene overlap between the top 1000 blue (differentiation) module genes and the 1000 genes most significantly up-regulated with Oct4 RNAi (right). Significance of differential expression was determined using Student's t-statistic. p-values have been corrected for multiple hypothesis tests (Benjamini-Hochberg). Only significantly enriched functional groups are shown.CC BY · PubMed Central
Figure 7 Transcriptional Regulators related to Pluripotency and Differentiation in the Zhou Network. TF and Suz12 binding in the promoter regions of highly connected genes related to ES cell pluripotency and self-renewal with GO terms of transcriptional regulation or chromatin structure. Genes are listed by black kME (positive, left and negative, right) along with their corresponding significance level (log10 of the Bonferroni corrected p-value generated by a correlation test). Binding data from Chen et al (2008), Boyer et al (2006), and Loh et al (2006), are marked in blue (bound) and beige (unbound). For Oct4, Sox2, and Suz12, where binding is given by two studies, binding will be blue if it is found in both studies and light blue if found in only one.
Figure 7 Transcriptional Regulators related to Pluripotency and Differentiation in the Zhou Network. TF and Suz12 binding in the promoter regions of highly connected genes related to ES cell pluripotency and self-renewal with GO terms of transcriptional regulation or chromatin structure. Genes are listed by black kME (positive, left and negative, right) along with their corresponding significance level (log10 of the Bonferroni corrected p-value generated by a correlation test). Binding data from Chen et al (2008), Boyer et al (2006), and Loh et al (2006), are marked in blue (bound) and beige (unbound). For Oct4, Sox2, and Suz12, where binding is given by two studies, binding will be blue if it is found in both studies and light blue if found in only one.CC BY · PubMed Central
Figure 8 Non-Transcriptional Regulators Related to Pluripotency and Differentiation. TFs and Suz12 binding of highly connected genes related to ES cell pluripotency and self-renewal lacking GO terms for transcriptional regulation or chromatin structure. Genes are tablulated in the same format as Figure 7.
Figure 8 Non-Transcriptional Regulators Related to Pluripotency and Differentiation. TFs and Suz12 binding of highly connected genes related to ES cell pluripotency and self-renewal lacking GO terms for transcriptional regulation or chromatin structure. Genes are tablulated in the same format as Figure 7.CC BY · PubMed Central
Figure 9 Pluripotency Transcriptional Regulators that are not Bound by the Core TF Machinery in ES Cells. Genes related to ES cell pluripotency and self-renewal with GO terms of transcriptional regulation or chromatin structure and little pluripotency TF binding. Genes are listed by black kME along with their corresponding significance level (log10 of the Bonferroni corrected p-value generated by a correlation test). Binding data from Chen et al, Boyer et al and Loh et al, are marked in blue (bound) and beige (unbound). For Suz12, light blue indicates that a genes is called bound in Chen et al or Boyer et al but not in both.
Figure 9 Pluripotency Transcriptional Regulators that are not Bound by the Core TF Machinery in ES Cells. Genes related to ES cell pluripotency and self-renewal with GO terms of transcriptional regulation or chromatin structure and little pluripotency TF binding. Genes are listed by black kME along with their corresponding significance level (log10 of the Bonferroni corrected p-value generated by a correlation test). Binding data from Chen et al, Boyer et al and Loh et al, are marked in blue (bound) and beige (unbound). For Suz12, light blue indicates that a genes is called bound in Chen et al or Boyer et al but not in both.CC BY · PubMed Central

Verification

Claims the paper makes

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

The self-check needs the paper’s full text.

Citations

Cites

Show all 77 references

Cited by

Sorted by citations

Metadata

Authors
Mason, M. J., Fan, G., Plath, K., Zhou, Q., & Horvath, S.
Affiliations
University of California, Los Angeles (US)
Corresponding author
  • Zhou, Qing
  • Horvath, Steve
Venue
BMC Genomics
Published
2009-07-20
Volume
10
Issue
1
Pages
327-327
ISSN
1471-2164 (ISSN-L)
OpenAlex
W2138720219
Cited by
243
References
77
Citations by year
2026: 8 · 2025: 10 · 2024: 12 · 2023: 12 · 2022: 18 · 2021: 34 · 2020: 18 · 2019: 18 · 2018: 18 · 2017: 12 · 2016: 10 · 2015: 8 · 2014: 14 · 2013: 22 · 2012: 11
FWCI
3.89
Type
journal article
Language
en
Open access
Licence
cc-by
Locations
  • BMC Genomics (publishedVersion, cc-by)
  • PubMed Central (publishedVersion, cc-by)
  • PubMed (publishedVersion)
  • DOAJ (DOAJ: Directory of Open Access Journals) (submittedVersion)
  • eScholarship (California Digital Library) (submittedVersion)
  • The HKU Scholars Hub (University of Hong Kong) (submittedVersion)
Funders
  • National Science Foundation (0805491)
  • National Institutes of Health (1dp2od001686-01)
  • National Cancer Institute (P50CA092131, P30CA016042)
  • National Institute of General Medical Sciences (P01GM081621)
  • National Institute of Allergy and Infectious Diseases (U19AI063603)
  • National Institute of Diabetes and Digestive and Kidney Diseases (R01DK072206, DK-072206)
Topics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Pluripotent Stem Cells Research
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Bioinformatics and Genomic Networks
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Gene Regulatory Network Analysis
Keywords
mouse embryonic stem cells, pluripotency, self-renewal, gene co-expression network, transcriptional regulation, transcription factor binding, DNA repair, mitochondrial function, histone modifications, promoter CpG methylation, epigenetic regulation, module membership
MeSH terms
Animals, Chromatin, DNA, DNA Damage, DNA Repair, Gene Expression Regulation, Protein Binding, Transcription Factors, Transcription, Genetic, CpG Islands, DNA Methylation, Gene Expression Profiling, Pluripotent Stem Cells, Mice, Gene Regulatory Networks, Embryonic Stem Cells

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review