Tutorials › High-Dimensional Genomic Data Analysis › Nodes, Edges, And Co-Expression Network Estimands

Co-Expression Networks · Tutorial 484 of 795

Nodes, Edges, And Co-Expression Network Estimands

Learn to state what a co-expression network’s nodes and edges represent, and why an observed expression association is not automatically a regulatory or causal relationship.

Advanced 11 min read

What You'll Learn

  • Define the gene universe represented by a network’s nodes.
  • Distinguish a population edge estimand from an edge estimated from observed expression data.
  • Explain how signed, unsigned, threshold-based, and significance-based edge definitions differ.
  • Identify how the target population and sampling mixture affect co-expression.
  • Separate expression association from regulatory, causal, and physical-interaction claims.

What a Co-Expression Network Represents

A co-expression network is a representation of relationships among genes based on their measured expression across biological units, such as samples, individuals, or experimental specimens. Its apparent simplicity can hide several different statistical questions: which genes are eligible to appear, what relationship between two genes counts as an edge, and in what population or experimental context that relationship is defined. A network is interpretable only when those choices are made explicit.

The central distinction is between the estimand and the network calculated from data. The estimand is the population quantity or relationship of scientific interest. An estimator uses observed expression values to approximate it. A plotted edge is therefore not, by itself, a complete description of what has been estimated.

Definition: A node is an element of the specified gene universe represented in a network. An edge is a defined relationship between a pair of nodes. In a co-expression network, that relationship describes an expression association in a stated target population or experimental context; it does not, on its own, assert regulation or causation.

Let \(X_g\) denote the expression measurement for gene \(g\) on a randomly selected biological unit from a specified target population. The population distribution of the expression vector depends on what units are sampled, the biological conditions represented, and how measurements are defined. For a pair of genes \(g\) and \(h\), one possible association estimand is their population correlation:

$$ \rho_{gh} = \frac{\operatorname{Cov}(X_g,X_h)} {\sqrt{\operatorname{Var}(X_g)\operatorname{Var}(X_h)}} $$

This quantity is defined when both genes have finite, nonzero population variances. It is one possible basis for an edge, not the universal meaning of an edge. Other definitions may use covariance, a rank-based association, or another specified pairwise relationship. Choosing among correlation measures is a separate methodological decision; here the important point is that the chosen quantity defines what the edge is intended to mean.

A target population must be specified to make an estimand meaningful. For example, association across untreated samples from one tissue is not necessarily the same estimand as association across samples pooled from several tissues and treatment groups. A gene universe also needs to be defined: it may include all genes passing stated measurement and filtering rules, or a narrower set chosen for a scientific purpose. The node list is thus a substantive part of the network definition, not just a visual convenience.

Edge Definitions and Their Interpretation

A useful formal definition separates the pairwise population quantity from the rule that turns it into an edge. For instance, if the target is correlation, an edge rule could declare a pair connected when the magnitude of its population correlation exceeds a specified cutoff \(c\). This defines a population network. In practice, the population correlation is unknown, so a network is commonly constructed by applying a stated rule to an estimate or test result from the data. That data-based network is an estimated object, and its edges may differ from the population edges.

Definition: An edge estimand specifies (1) the target population and expression context, (2) the pairwise relationship of interest, and (3) the rule determining whether that relationship qualifies as an edge. An estimated edge is the result of applying an estimation and network-building procedure to observed data.

Common edge meanings are not interchangeable. A signed edge retains whether the measured association is positive or negative. An unsigned edge uses association magnitude, so strong positive and strong negative associations can both qualify. A threshold-based rule asks whether an estimated association meets a specified magnitude criterion. A significance-based rule instead uses a statistical test and a stated decision rule; its result depends on the test, error-rate procedure, and chosen level. These rules answer different questions, even if they produce networks with similar-looking diagrams.

Some network displays encode the magnitude of an association as an edge weight; others retain only whether an edge is present. A weight is meaningful only with its definition and scale stated. For instance, a weight equal to correlation preserves sign, whereas a weight equal to absolute correlation does not. In either case, it represents the specified association summary, not the strength of a causal mechanism.

Worked Example: Applying a Signed Correlation Edge Rule

Suppose five hypothetical biological units have centered expression measurements for three genes, \(A\), \(B\), and \(C\). Let the values for \(A\) be \((-2,-1,0,1,2)\), for \(B\) be \((-4,-2,0,2,4)\), and for \(C\) be \((2,-1,-2,-1,2)\). Consider an estimated network that connects a pair when the absolute sample correlation \( |r| \) is at least \(0.70\).

For \(A\) and \(B\), both sample means are zero. The sums of squares are \(\sum A_i^2=10\) and \(\sum B_i^2=40\), and the cross-product sum is \(\sum A_iB_i=20\). Thus:

$$ r_{AB} = \frac{\sum_i A_iB_i} {\sqrt{\left(\sum_i A_i^2\right)\left(\sum_i B_i^2\right)}} = \frac{20}{\sqrt{10(40)}} = 1 $$

Since \( |r_{AB}|=1\geq 0.70\), the estimated network contains an \(A\)-\(B\) edge. For \(A\) and \(C\), \(\sum_i A_iC_i=0\), while their centered sums of squares are \(10\) and \(14\). Therefore \(r_{AC}=0/\sqrt{10(14)}=0\), so this rule does not include an \(A\)-\(C\) edge. This calculation establishes what the specified sample rule does with these data. It does not establish that the corresponding population correlations are exactly 1 and 0.

Why Population and Context Matter

A co-expression association is relative to the distribution of units being analyzed. If the samples represent several conditions, a pooled association can reflect differences between conditions even when the genes show no association within any one condition. Conversely, analyzing only one condition targets a different distribution and can yield a different edge. The estimand is not fully defined by naming two genes and a statistic; the population and sampling context matter as well.

Worked Example: A Pooled Edge Created by Group Differences

Consider eight hypothetical units divided equally between two conditions. In condition 1, the values for genes \(A\) and \(B\) are \((-1,-1,1,1)\) and \((-1,1,-1,1)\). In condition 2, add 10 to every value for both genes, giving \(A=(9,9,11,11)\) and \(B=(9,11,9,11)\). Within each condition, the centered cross-product sum is zero: in condition 1 it is \(1-1-1+1=0\), and adding a constant to both variables does not change centered cross-products. Thus the within-condition sample correlation is zero in each group.

For the pooled data, each gene has mean 5. The four condition-1 observations have centered values with a group offset of \(-5\), and the four condition-2 observations have a group offset of \(+5\). Within each group, the residual cross-products sum to zero. The pooled centered cross-product sum is therefore \(4(25)+4(25)=200\). For each gene, the pooled centered sum of squares is \(200+8=208\): 200 comes from the group offsets and 8 from the within-group values. Hence:

$$ r_{AB,\text{pooled}} = \frac{200}{\sqrt{208(208)}} = \frac{200}{208} \approx 0.9615 $$

A pooled threshold rule with cutoff \(0.70\) would include an estimated edge, despite zero within-condition correlations. The pooled edge summarizes the combined sample distribution, where condition membership shifts both genes together. It does not show that the genes are associated within either condition. Whether the pooled or within-condition relationship is scientifically relevant depends on the question and target population.

This example is a reason to report the sampling context and, where appropriate, examine whether associations differ across conditions. It is not a reason to assume that either a pooled or stratified network is always preferable. Each targets a different relationship. Expression normalization, gene filtering, batch structure, and the definition of a biological replicate can also affect the data used to estimate edges and therefore belong in a reproducible network description.

Association Is Not Regulation or Causation

An undirected co-expression edge states that two expression measurements satisfy a specified association rule in the target population. It does not identify which gene influences the other, whether either gene directly regulates the other, or whether the association is produced by a third variable. Shared cell composition, a common upstream response, technical effects, or condition differences could all contribute to an observed association.

A regulatory network makes a stronger claim: its directed links represent regulatory relationships, with direction carrying meaning. A causal network likewise concerns effects under interventions or a causal model with defensible assumptions. Establishing those interpretations generally requires evidence and design beyond pairwise co-expression. A physical-interaction network is different again: its edges represent molecular interactions, not necessarily correlated expression. The same genes can therefore appear in several network types with different edge meanings.

A common mistake is to look at a positive co-expression edge and conclude, “Gene \(A\) activates gene \(B\).” The correction is to say, “Their measured expression values meet the stated positive-association edge rule in this sample and target context.” The latter claim describes the analysis. Directional regulation or causation requires additional evidence and a suitable design.

A Short Justification: Why the Target Mixture Changes Correlation

The pooled example follows from the covariance decomposition by group. Let \(G\) identify the condition, and let \(X\) and \(Y\) be the expression values of two genes. Writing \(\mu_{X,G}=E(X\mid G)\) and \(\mu_{Y,G}=E(Y\mid G)\), decompose each centered value into a within-group deviation and a group-mean deviation:

$$ X-E(X) = \{X-\mu_{X,G}\}+\{\mu_{X,G}-E(X)\} $$

The analogous decomposition holds for \(Y-E(Y)\). On multiplying these expressions and taking expectations, the cross terms have expectation zero because \(E\{X-\mu_{X,G}\mid G\}=0\) and \(E\{Y-\mu_{Y,G}\mid G\}=0\). Consequently:

$$ \operatorname{Cov}(X,Y) = E\{\operatorname{Cov}(X,Y\mid G)\} + \operatorname{Cov}\{\mu_{X,G},\mu_{Y,G}\} $$

The pooled covariance contains both average within-group covariance and covariance between group means. In the worked example, the first term is zero and the second is positive. This explains why the pooled association is strong even though neither group has an internal association. Correlation additionally divides covariance by the pooled standard deviations, which also reflect within- and between-group variation.

Defining and Checking a Network

Before interpreting a network, make its estimand and construction auditable. A concise description should identify:

  • the node universe and any expression or quality filters;
  • the biological units, conditions, and target population represented;
  • the pairwise association quantity and whether edges retain sign;
  • the rule used to include or weight edges, including any cutoff or testing decision;
  • the expression preprocessing and the distinction between the population target and its data-based estimate.

A diagnostic check is whether the interpretation matches the edge rule. If the rule uses absolute association, the resulting edge does not encode positive versus negative direction. If it uses a pooled sample, it describes the pooled distribution and may reflect group structure. If the rule is based on statistical significance, an edge represents a testing decision under that procedure; it is not a direct measure of biological importance. In small or unrepresentative samples, estimated relationships may also be unstable, so a displayed edge should not be treated as certain merely because it appears in a graph.

There is no edge definition that is best for every scientific aim. A signed association can be useful when direction of co-variation matters; an unsigned one may be appropriate when the aim is to group genes with strong expression relationships regardless of sign. A significance-based decision addresses evidence against a specified null under a testing framework, while a magnitude rule directly encodes a chosen association scale. The appropriate choice depends on the scientific question, and different choices may produce different networks from the same measurements.

Key takeaway: A co-expression network has meaning only after its nodes, target population, pairwise relationship, and edge rule are specified. Its edges describe expression association as defined by those choices; they do not by themselves establish regulation, physical interaction, or causation.

Check Your Understanding

Answer each question by stating the relevant network definition or calculation. Brief solutions follow.

  1. A network includes pairs when the absolute sample correlation is at least \(0.8\). What does an edge say, and what does it not say?
  2. In a signed network, a pair has estimated correlation \(-0.9\). What is the association direction, and would an unsigned rule based on absolute correlation retain it at a cutoff of \(0.8\)?
  3. Why can a pooled network contain an edge even when within-condition associations are zero?
  4. Name two pieces of information, besides the list of genes, needed to interpret a co-expression network.
  5. A researcher interprets an undirected co-expression edge as proof that one gene regulates the other. Identify the inferential error.

Solutions. (1) The observed pair meets the stated sample-correlation magnitude rule; that alone does not establish its exact population association or a regulatory effect. (2) The estimated association is negative. An unsigned rule retains it because \( |-0.9|=0.9\geq0.8\). (3) Covariance between group means can contribute to pooled covariance even when average within-group covariance is zero. (4) Examples include the target population and conditions, the association measure, the edge rule, and whether sign is retained. (5) The researcher has treated an undirected observational association as directional causal evidence, which the edge definition does not provide.

Clinical Biostats