<title>Introduction</title> <p>Many biological networks share topological properties. Common global properties include modular organization ##REF##12202830##[1]##,##REF##15044247##[2]##, the presence of highly connected hub nodes, and approximate ‘scale free topology’ ##REF##14735121##[3]##,##REF##16254242##[4]##. Common local topological properties include the presence of recurring patterns of interconnections (‘network motifs’) in regulation networks ##REF##12399590##[5]##–##REF##15908506##[7]##.</p> <p>One goal of this article is to describe existing and novel network concepts (also known as network statistics or indices ##REF##17547772##[8]##) that can be used to describe local and global network properties. For example, the clustering coefficient ##REF##9623998##[9]## is a network concept, which measures the cohesiveness of the neighborhood of a node. We are particularly interested in network concepts that are defined with regard to a ‘gene significance measure’. Gene significance measures are of great practical importance since they allow one to incorporate external gene information into the network analysis. In functional enrichment analysis, a gene significance measure could indicate pathway membership. In gene knock-out experiments, gene significance could indicate knock-out essentiality. We study gene significance measures since a microarray sample trait (e.g., case control status) gives rise to a statistical measure of gene significance. For example, the Student <italic>t</italic>-test of differential expression leads to a gene significance measure. Many traditional microarray data analysis methods focus on the relationship between the microarray sample trait and the gene expression data. For example, gene filtering methods aim to find a list of (differentially expressed) genes that are significantly associated with the microarray sample trait; another example are microarray-based prediction methods that aim to accurately predict the sample trait on the basis of the gene expression data.</p> <p>Gene expression profiles across microarray samples can be highly correlated and it is natural to describe their pairwise relations using network language. Genes with similar expression patterns may form complexes, pathways, or participate in regulatory and signaling circuits ##REF##9843981##[10]##–##REF##17646300##[12]##. Gene coexpression networks have been used to describe the transcriptome in many organisms, e.g., yeast, flies, worms, plants, mice, and humans ##REF##11027309##[13]##–##UREF##1##[23]##. Gene coexpression network methods have also been used for typical microarray data analysis tasks such as gene filtering ##UREF##0##[19]##, ##REF##17090670##[24]##–##REF##17668265##[26]## and outcome prediction ##REF##16926220##[27]##,##REF##17940530##[28]##.</p> <p>While the utility of network methods for analyzing microarray data has been demonstrated in numerous publications, the utility of microarray data analysis techniques for solving network theoretic problems has not yet been fully appreciated. One goal of this article is to show that simple geometric arguments can be used to derive network theoretic results if the networks are defined on the basis of a correlation matrix.</p> <title>Definition of Gene Coexpression Networks</title> <p>Although many of our network concepts will be useful for general networks, we are particularly interested in gene coexpression networks (also known as association-, influence-, relevance-, or correlation networks). Gene coexpression networks are built on the basis of a gene coexpression measure. The network nodes correspond to genes—or more precisely to gene expression profiles. The <italic>i</italic>th gene expression profile <italic>x<sub>i</sub></italic> is a vector whose components report the gene expression values across <italic>m</italic> microarrays. We define the coexpression similarity <italic>s<sub>ij</sub></italic> between genes <italic>i</italic> and <italic>j</italic> as the absolute value of the correlation coefficient between their expression profiles:
</p> <p>Using a thresholding procedure, this coexpression similarity is transformed into a measure of connection strength (adjacency). An unweighted network adjacency <italic>a<sub>ij</sub></italic> between gene expression profiles <italic>x<sub>i</sub></italic> and <italic>x<sub>j</sub></italic> can be defined by hard thresholding the coexpression similarity <italic>s<sub>ij</sub></italic> as followswhere τ is the “hard” threshold parameter. Thus, two genes are linked (<italic>a<sub>ij</sub></italic> = 1) if the absolute correlation between their expression profiles exceeds the (hard) threshold τ. Hard thresholding of the correlation leads to simple network concepts (e.g., the gene connectivity equals the number of direct neighbors) but it may lead to a loss of information: if τ has been set to 0.8, there will be no link between two genes if their correlation equals 0.799. To preserve the continuous nature of the coexpression information, one could simply define a weighted adjacency matrix as the absolute value of the gene expression correlation matrix, i.e., [<italic>a<sub>ij</sub></italic>] = [<italic>s<sub>ij</sub></italic>]. However, since microarray data can be noisy and the number of samples is often small, we and others have found it useful to emphasize strong correlations and to punish weak correlations. It is natural to define the adjacency between two genes as a power of the absolute value of the correlation coefficient ##UREF##0##[19]##,##REF##17090670##[24]##:with <italic>β</italic>≥1. This soft thresholding approach leads to a weighted gene coexpression network. We present empirical results for weighted and unweighted networks in the main text, ##SUPPL##0##Text S1##, ##SUPPL##1##Text S2##, and ##SUPPL##2##Text S3##.</p> <title>Social Network Analogy: Affection Network</title> <p>Since humans are organized into social networks, social network analogies should be intuitive to many readers. Therefore, we will refer to the following ‘affection network’ throughout this article. Assume that <italic>n</italic> individuals filled out an interest questionnaire, which was used to define a pairwise similarity score <italic>s<sub>ij</sub></italic>. For convenience, we assume that the similarity measure takes on values between 0 and 1. Our definition of the affection network is based on the following assumption: the more similar the interests between two individuals, the more affection they feel for each other. More specifically, we assume that the affection (adjacency) <italic>a<sub>ij</sub></italic> between two individuals is proportional to their similarity on a logarithmic scale, i.e.,This is equivalent to our soft thresholding approach <italic>a<sub>ij</sub></italic> = <italic>s<sub>ij</sub><sup>β</sup></italic> (Equation 2). A soft threshold <italic>β</italic> = 2 implies that the affection <italic>a<sub>ij</sub></italic> equals 0.25 if the similarity <italic>s<sub>ij</sub></italic> equals 0.5.</p>
<title>Materials and Methods</title> <title>Network Concept Functions and Fundamental Network Concepts</title> <p>Analogous to ##REF##17547772##[8]##, we define a <italic>network concept function</italic> to be function of a square matrix <italic>M</italic> = [<italic>M<sub>ij</sub></italic>] (1≤<italic>i</italic>,<italic>j</italic>≤<italic>n</italic>) and/or a corresponding vector <italic>G</italic> = (<italic>G</italic>
<sub>1</sub>,…,<italic>G<sub>n</sub></italic>). For example, <italic>M</italic> could be the adjacency matrix (with diagonal set to 0) and <italic>G</italic> could be a corresponding gene significance measure.</p> <p>We make use of the following network concept functions:where the components of matrix <italic>B<sub>M</sub></italic> in the denominator of the clustering coefficient function are given by <italic>b<sub>ij</sub></italic> = 1 if <italic>i</italic>≠<italic>j</italic> and <italic>b<sub>ii</sub></italic> = Ind(<italic>m<sub>ii</sub></italic>&gt;0). Here the indicator function Ind(·) takes on the value 1 if the condition is satisfied and 0 otherwise.</p> <p>According to our convention, the diagonal elements of the adjacency matrix are set to 1. Therefore, the diagonal elements of <italic>A–I</italic> (where <italic>I</italic> denotes the identity matrix) equal 0. Now we are ready to define the (fundamental) network concepts that are studied in this article.</p> <p>
<bold>Definition of Fundamental Network Concepts:</bold>
<italic>The fundamental network concepts of a network A are defined by evaluating the network functions (</italic>
<italic>Equation 42</italic>
<italic>) on A–I and the gene significance measure GS, i.e.,</italic>

</p> <p>For example, the connectivity is given by
</p> <p>We define an <bold>intramodular network concept</bold>
<italic>NCF</italic>(<italic>A</italic>
<sup>(<italic>q</italic>)</sup>−<italic>I</italic>,<italic>GS</italic>
<sup>(<italic>q</italic>)</sup>) by evaluating the network concept function on the restricted adjacency matrix <italic>A</italic>
<sup>(<italic>q</italic>)</sup> and the restricted gene significance measure <italic>GS</italic>
<sup>(<italic>q</italic>)</sup>.</p> <p>We will now define eigengene-based network concepts. Using the eigengene-based adjacency matrix <italic>A<sub>E</sub></italic>
<sup>(<italic>q</italic>)</sup> = <italic>a<sub>e</sub></italic>
<sup>(<italic>q</italic>)</sup>(<italic>a<sub>e</sub></italic>
<sup>(<italic>q</italic>)</sup>)<italic><sup>T</sup></italic> (Equation 28) and the eigengene-based gene significance measure <italic>GS<sub>E,i</sub></italic>
<sup>(<italic>q</italic>)</sup> = <italic>a<sub>e,i</sub></italic>
<sup>(<italic>q</italic>)</sup>
<italic>a<sub>e,t</sub></italic>
<sup>(<italic>q</italic>)</sup> (Equation 29), we define an <bold>eigengene-based network concept</bold> as <italic>NCF</italic>(<italic>A<sub>E</sub></italic>
<sup>(<italic>q</italic>)</sup>,<italic>GS<sub>E</sub></italic>
<sup>(<italic>q</italic>)</sup>).</p> <p>As example, consider the eigengene-based connectivity given by
</p> <title>Deriving Observation 1: Expression Data with High Eigengene Factorizability Lead to Approximately Factorizable Networks</title> <p>Here we derive Observation 1, which characterizes approximately factorizable gene coexpression module networks. To simplify the presentation, we omit the superscripts (<italic>q</italic>) in the following, e.g., we will write <italic>EF</italic>(<italic>X</italic>) instead of <italic>EF</italic>(<italic>X</italic>
<sup>(<italic>q</italic>)</sup>). We will argue that if the eigengene factorizability <italic>EF</italic>(<italic>X</italic>) is close to 1, the adjacencies of the weighted coexpression module network <italic>A</italic> = |cor(<italic>X</italic>)|<italic><sup>β</sup></italic> and the trait-based gene significance measure <italic>GS<sub>i</sub></italic> = |cor(<italic>x<sub>i</sub></italic>,<italic>T</italic>)|<italic><sup>β</sup></italic> can be factored as followswhere

</p> <p>Since our gene coexpression networks are defined with respect to the correlation matrix [cor(<italic>x<sub>i</sub></italic>,<italic>x<sub>j</sub></italic>)], which is scale-invariant, we can assume that the gene expression profiles have been scaled as follows:  where <italic>m</italic> is the number of microarray samples. Then one can derive the following relationshipsNote that <italic>u</italic>
<sub>1,<italic>i</italic></sub>|<italic>d</italic>
<sub>1</sub>|<sup>2</sup>
<italic>u</italic>
<sub>1,<italic>j</italic></sub>/<italic>m</italic> = cor(<italic>x<sub>i</sub></italic>,<italic>E</italic>)cor(<italic>x<sub>j</sub></italic>,<italic>E</italic>). Using the fact that <italic>U</italic> is an orthogonal matrix, it is straightforward to show that
</p> <p>This equation motivates us to propose the following measure of <bold>eigengene factorizability</bold>:Note that 0≤<italic>EF</italic>(<italic>E</italic>)≤1. By definition <italic>EF</italic>(<italic>E</italic>)≈1 implies thatBy raising both sides of this equation to a power <italic>β</italic>, we findThe last step highlights an important theoretical advantage of the soft thresholding method: it preserves the approximate factorizability of the underlying correlation matrix.</p> <p>An alternative, possibly more direct way of motivating the observation is based on the insight that the squared singular values <italic>|d<sub>l</sub>|<sup>2</sup></italic> correspond to the eigenvalues of the correlation matrix <italic>COR</italic> = [cor(<italic>x<sub>i</sub></italic>,<italic>x<sub>j</sub></italic>)]. For high values of <italic>EF</italic>(<italic>E</italic>), the correlation matrix can be factored as followswhere <underline>u</underline>
<sub>1</sub> denotes an eigenvector of length 1.</p> <title>Relationships among Network Concepts When the Maximum Conformity Assumption Does Not Hold</title> <p>Here we describe relationships among eigengene-based network concepts if the maximum conformity assumption does not hold (i.e., <italic>a<sub>e</sub></italic>
<sub>,max</sub>
<sup>(<italic>q</italic>)</sup>&lt;&lt;1). It is straightforward to derive the following relationships among eigengene-based network concepts:Observation 2 can be used to derive the following</p> <title>Observation 4</title> <p>
<italic>If A<sup>(q)</sup> = |cor(X<sup>(q)</sup>)|<sup>β</sup> and the eigengene factorizability is close to 1 (EF(X<sup>(q)</sup>)≈1), the relationships among eigengene-based concepts approximately apply to their network analogs as well.</italic>
</p> <p>For example, we find that
</p> <title>Deriving the Geometric Interpretation of Factorizability</title> <p>In the following we provide details on our geometric interpretation of the factorizability. To simplify the notation, we sometimes drop the superscript (<italic>q</italic>) in the following expressions. We denote by <italic>θ<sub>l,i</sub></italic> the angle between the right singular vector <italic>v<sub>l</sub></italic> (Equation 20) and the <italic>i</italic>th gene expression profile <italic>x<sub>i</sub></italic>. The smaller the angle <italic>θ<sub>l,i</sub></italic>, the bigger the correlation cor(<italic>v<sub>l</sub></italic>,<italic>x<sub>i</sub></italic>) = cos(<italic>θ<sub>l</sub></italic>
<sub>,<italic>i</italic></sub>). Using , one can reexpress the eigengene factorizability (Equation 24) as follows
</p> <p>Thus, <italic>EF</italic>(<italic>X</italic>
<sup>(<italic>q</italic>)</sup>)≈1 if the module gene expressions <italic>x<sub>i</sub></italic> are approximately orthogonal (cos(<italic>θ<sub>l</sub></italic>
<sub>,<italic>i</italic></sub>)≈0) to the right singular vectors <bold><italic>v</italic></bold>
<italic><sub>l</sub></italic> for <italic>l</italic>≥2, i.e., if on average the gene expression profiles point in the direction of the module eigengene <italic>v</italic>
<sub>1</sub> = <italic>E</italic>.</p> <p>Under this assumption, we provide a rough geometric intuition of <italic>a<sub>ij</sub></italic>≈<italic>a<sub>e,i</sub>a<sub>e,j</sub></italic> (Equation 25) depicted in ##FIG##4##Figure 5A##. We denote by <italic>θ<sub>i</sub></italic> = <italic>θ</italic>
<sub>1,<italic>i</italic></sub> the angle between the module eigengene <italic>E</italic> and the <italic>i</italic>th gene expression profile and by <italic>θ<sub>ij</sub></italic> the angle between gene expression profiles <italic>i</italic> and <italic>j</italic>. Using the assumptions described in ##FIG##4##Figure 5A##, <italic>θ<sub>ij</sub></italic>≈|<italic>θ<sub>i</sub></italic>±<italic>θ<sub>j</sub></italic>| and sin(<italic>θ<sub>i</sub></italic>) sin(<italic>θ<sub>j</sub></italic>)≈0, we find thati.e., the correlation matrix is approximately factorizable.</p> <title>Heterogeneity Increases with the Soft Threshold <italic>β</italic>
</title> <p>Here we prove that the eigengene-based heterogeneity increases with the soft threshold <italic>β</italic> (Equation 2). Recall that  (Equation 30) which implies that it is a decreasing function ofNote that <italic>a<sub>i</sub></italic> = |cor(<italic>x<sub>i</sub></italic>,<italic>E</italic>)| is a nonnegative number.</p> <p>To prove that the heterogeneity increases with <italic>β</italic>, it suffices to prove the following</p> <p>
<bold>Proposition:</bold>
<italic> Let {a<sub>i</sub>, i = 1,…,n} be a group of nonnegative number and β&gt;1 then the following inequality holds</italic>:
</p> <p>To prove the Proposition, we will make use of the following</p> <p>
<bold>Lemma:</bold>
<italic>Let {u<sub>i</sub>, i = 1,…,n} and {v<sub>i</sub>, i = 1,…,n} be groups of nonnegative numbers, and θ be a number 0≤θ&lt;</italic>1. <italic>Then the following inequality holds:</italic>
The Lemma can be proved with Hölder's inequality, which is given byWe use the Lemma with <italic>θ</italic>
<sub>1</sub> = <italic>β</italic>/(<italic>2β</italic>−1), <italic>u<sub>i</sub></italic> = <italic>a<sub>i</sub></italic>, and <italic>v<sub>i</sub></italic> = <italic>a<sub>i</sub><sup>2β</sup></italic> to deriveFurther, we use the Lemma with <italic>θ</italic>
<sub>2</sub> = (<italic>2β</italic>−2)/(2<italic>β</italic>−1), <italic>u<sub>i</sub></italic> = <italic>a<sub>i</sub></italic>, and <italic>v<sub>i</sub></italic> = <italic>a<sub>i</sub><sup>2β</sup></italic> to deriveBy squaring the first inequality and multiplying it with the second inequality, we arrive atsince 2<italic>θ</italic>
<sub>1</sub>+<italic>θ</italic>
<sub>2</sub> = 2 and 3−(2<italic>θ</italic>
<sub>1</sub>+<italic>θ</italic>
<sub>2</sub>) = 1. The last inequality completes the proof since it is equivalent to the inequality in Equation 53.</p>
<title>Results</title> <title>Gene Significance Based on a Microarray Sample Trait</title> <p>Many network applications use at least one gene significance measure. Abstractly speaking, we define a gene significance measure as a function <italic>GS</italic> that assigns a nonnegative number to each gene; the higher <italic>GS<sub>i</sub></italic> the more <italic>biologically</italic> significant is gene <italic>i</italic>. We assume that the minimum gene significance is 0. For example, if a statistical significance level (<italic>p</italic>-value) is available for each gene, the gene significance of the <italic>i</italic>th gene can be defined as minus log of the <italic>p</italic>-value, i.e., <italic>GS<sub>i</sub></italic> = −log(<italic>p<sub>i</sub></italic>). In this article, we are particularly interested in gene significance measures that are based on a microarray sample trait, e.g., a clinical outcome. The microarray sample trait <italic>T</italic> = (<italic>T</italic>
<sub>1</sub>,…,<italic>T<sub>m</sub></italic>) may be quantitative (e.g., body weight) or binary (e.g., case control status). Since our goal is to provide a simple geometric interpretation of coexpression network analysis, we define the <bold>trait-based gene significance measure</bold> by raising the correlation between the <italic>i</italic>th gene expression profile <italic>x<sub>i</sub></italic> and the clinical trait <italic>T</italic> to a power <italic>β</italic>
Although any power <italic>β</italic> could be used in Equation 4, we use the same power as in Equation 2 to facilitate a simple geometric interpretation.</p> <title>Geometric Interpretation Using a Hypersphere</title> <p>We find it convenient to express network quantities in terms of correlation coefficients since the correlation between two vectors can be interpreted as the cosine of the angle between them (measured in radians) if the vectors are scaled to have a mean of 0. Since the correlation is scale-invariant, i.e., cor(<italic>ax<sub>i</sub></italic>+<italic>b</italic>, <italic>cx<sub>j</sub></italic>+<italic>d</italic>) = cor(<italic>x<sub>i</sub></italic>,<italic>x<sub>j</sub></italic>), we can assume without loss of generality that the vectors <italic>x<sub>i</sub></italic> have a mean 0 and are of the same length. In other words, they correspond to points on a hypersphere.</p> <p>The network adjacency <italic>a<sub>ij</sub></italic> is a monotonically decreasing function of the angle <italic>θ<sub>ij</sub></italic> between the two scaled expression profiles if 0≤<italic>θ<sub>ij</sub></italic>≤<italic>π</italic>/2. When the angle <italic>θ<sub>ij</sub></italic> equals 0 or <italic>π</italic>/2, the adjacency equals 1 or 0, respectively. The network adjacency is a monotonically decreasing function of the length of the shortest path (geodesic) between the two points on the hypersphere. Soft thresholding methods (Equation 2) preserve the continuous nature of these distances. The higher the soft threshold <italic>β</italic>, the more weight is assigned to short geodesic distances compared to large distances.</p> <p>Since the trait-based gene significance measure <italic>GS<sub>i</sub></italic> = |cor(<italic>x<sub>i</sub></italic>,<italic>T</italic>)|<italic><sup>β</sup></italic>, (Equation 4) is scale-invariant, the sample trait <italic>T</italic> can also be considered a point on the hypersphere. Analogous to the network adjacency, the smaller the geodesic distance between the <italic>i</italic>th gene expression profile and the trait <italic>T</italic>, the higher the gene significance of the <italic>i</italic>th gene. In other words, the smaller the angle between the sample trait and the expression profile, the more significant is the gene.</p> <title>A Motivational Example</title> <p>As a motivational example, we study the pairwise correlations among 498 genes that had previously been found to form a sub-network related to mouse body weight. The microarray data measure the expression levels in multiple tissue samples (liver, adipose, brain, muscle) from male and female mice of an F2 intercross. Approximately 100 tissue samples are available for each gender/tissue combination. The biological significance of this subnetwork is described in ##UREF##1##[23]##,##REF##17668265##[26]##. Here we focus on the mathematical and topological properties of the pairwise absolute correlations <italic>a<sub>ij</sub></italic> = |cor(<italic>x<sub>i</sub></italic>,<italic>x<sub>j</sub></italic>)| between the genes. For each gender and tissue type ##FIG##0##Figure 1A## depicts a hierarchical cluster tree of the genes. ##FIG##0##Figure 1B## shows the corresponding heat maps, which color-code the absolute pairwise correlations <italic>a<sub>ij</sub></italic>. As can be seen from the color bar underneath the heat maps, red and green in the heat map indicate high and low absolute correlation, respectively. The genes in the rows and columns of each heat map are sorted by the corresponding cluster tree.</p> <p>It is visually obvious that the heat maps and the cluster trees of different gender/tissue combinations can look quite different. Network theory offers a wealth of intuitive concepts for describing the pairwise relationships among genes that are depicted in cluster trees and heat maps. To illustrate this point, we describe several such concepts in the following. By visual inspection of ##FIG##0##Figure 1B##, genes appear to be more highly correlated in liver than in adipose (a lot of red versus green color in the corresponding heat maps). This property can be captured by the concept of network density (defined below). The density of the female liver network is 0.39 while it is only 0.23 for the female adipose network. Another example for the use of network concepts is to quantify the extent of cluster (module) structure. In this example, branches of a cluster tree (##FIG##0##Figure 1A##) correspond to modules in the corresponding network. The cluster structure is also reflected in the corresponding heat maps: modules correspond to large red squares along the diagonal. Network theory provides a concept for quantifying the extent of module structure in a network: the mean clustering coefficient (defined below). The female liver, male liver and female brain networks have high mean clustering coefficients (mean <italic>ClusterCoef</italic> = 0.42, 0.43, 0.41, respectively). In contrast, the female adipose, male adipose, and male brain networks have lower mean clustering coefficients (mean <italic>ClusterCoef</italic> = 0.27, 0.27, 0.25, respectively). Difference in module structure may reflect true biological differences or they may reflect noise (e.g. technical artifacts or tissue contaminations).</p> <p>As another example for the use of network concepts, compare the cluster tree of the female brain network with that of the male brain network. The cluster tree of the female network appears to be comprised of a single large branch, i.e., a highly connected hub gene at the tip of the branch forms the center in this network. In contrast, the cluster tree corresponding to the male brain network appears to split into multiple smaller branches, i.e., no single gene forms the center. To measure whether a highly connected hub gene forms the center in a network, one can use the concept of centralization (defined below). The female brain and male brain networks have centralization 0.34 and 0.21, respectively.</p> <p>These examples illustrate that graph theory contains a wealth of network concepts that can be used to describe microarray data. But we will argue that microarray data analysis techniques can also be used to derive network theoretic results. For example, network theorists have long studied the relationship between gene significance and connectivity. Several network articles have pointed out that highly connected hub nodes are central to the network architecture ##REF##15130938##[17]##, ##REF##10935628##[29]##–##REF##15190252##[32]## but hub genes may not always be biologically significant ##REF##16515682##[33]##. To define a sample trait based gene significance measure (Equation 4), we define the gene significance of gene <italic>i</italic> as the absolute correlation between the gene expression profile <italic>x<sub>i</sub></italic> and body weight <italic>T</italic>, i.e., <italic>GS<sub>i</sub></italic> = |cor(<italic>x<sub>i</sub></italic>,<italic>T</italic>)|. ##FIG##0##Figure 1C## shows the relationship between this gene significance measure and connectivity in the different gender/tissue type networks. We find a strong positive relationship between gene significance and connectivity in the female and the male mouse liver networks. The positive relationship between gene significance and connectivity suggests that both variables could be used to implicate genes related to body weight. For example, we used connectivity as a variable in a systems biologic gene screening method ##REF##17668265##[26]##. While most network theorists would agree that connectivity is an important variable for finding important genes in a network ##REF##15130938##[17]##,##UREF##0##[19]##, the statistical advantages of combining gene significance and connectivity are not clear. Below, we use the geometric interpretation of coexpression network analysis to argue that intramodular connectivity can be interpreted as a fuzzy measure of module membership. Thus, a systems biologic gene screening method that combines a gene significance measure with intramodular connectivity amounts to a pathway based gene screening method. Empirical evidence shows that the resulting systems biologic gene screening methods can lead to important biological insights ##UREF##1##[23]##–##REF##17668265##[26]##. Before combining gene significance and connectivity in a systems biologic gene screening approach, it is important to study their relationship. Toward this end, we propose a measure of hub gene significance <italic>HGS</italic> as slope of a regression line (through the origin) between gene significance and scaled connectivity. As can be seen from ##FIG##0##Figure 1C##, the hub gene significance is high in liver and adipose tissues but it is low in brain and muscle tissues. Below, we use the geometric interpretation of coexpression networks to characterize coexpression networks that have high hub gene significance if the gene significance measure is based on a microarray sample trait <italic>T</italic>.</p> <title>Network Concepts</title> <title>Abstract definition of network concepts</title> <p>We define network concepts for (weighted) undirected networks that can be represented by a symmetric adjacency matrix <italic>A</italic> = [<italic>a<sub>ij</sub></italic>], where 1≤<italic>i</italic>,<italic>j</italic>≤<italic>n</italic>. We assume that the pairwise adjacency (connection strength) <italic>a<sub>ij</sub></italic> takes on values in the unit interval, i.e., 0≤<italic>a<sub>ij</sub></italic>≤1. For notational convenience, we set the diagonal elements to 1. In the <xref ref-type="sec" rid="s4">Methods</xref> section, we define a network concept <italic>NCF</italic>(<italic>A</italic>,<italic>GS</italic>) by evaluating a network concept function <italic>NCF</italic>(·,·) on the adjacency matrix <italic>A</italic> and/or a corresponding gene significance measure <italic>GS</italic>. This abstract definition will be useful in defining intramodular network concepts (e.g., Equation 17) and eigengene-based analogs of network concepts (e.g., Equation 30). In the following, we describe several network concepts including the connectivity, the maximum adjacency ratio, the density, and the centralization.</p> <title>Connectivity and related concepts</title> <p>The <italic>connectivity</italic> (also known as degree) of the <italic>i</italic>th gene is defined byIn unweighted networks, the connectivity <italic>k<sub>i</sub></italic> equals the number of genes that are directly linked to gene <italic>i</italic>. In weighted networks, the connectivity equals the sum of connection weights between gene <italic>i</italic> and the other genes.</p> <p>The <italic>maximum connectivity</italic> is defined asThe <italic>scaled connectivity K<sub>i</sub></italic> of the <italic>i</italic>-th gene is defined byBy definition, 0≤<italic>K<sub>i</sub></italic>≤1. Note that we distinguish the scaled from the unscaled connectivity by using an upper case “<italic>K</italic>” and a lower case “<italic>k</italic>”, respectively.</p> <p>
<italic>Social Network Interpretation of the Connectivity:</italic> For the aforementioned affection network (Equation 3), assume that the affection (adjacency) <italic>a<sub>ij</sub></italic> equals 1 if two individuals strongly like each other; it equals 0.5 if they are neutral towards each other, and it equals 0 if they strongly dislike each other. Then the scaled connectivity <italic>K<sub>i</sub></italic> is a measure of relative popularity: high values of <italic>K<sub>i</sub></italic> indicate that the <italic>i</italic>th person is well liked by many others.</p> <p>
<italic>Potential Uses of the Connectivity:</italic> The connectivity is the most widely used concept for distinguishing the nodes of a network. As described in the motivational example and detailed below, intramodular connectivity can be used to define a systems biologic gene screening strategy that keeps track of module membership information ##REF##17090670##[24]##.</p> <title>Maximum adjacency ratio</title> <p>For weighted networks, we define the <italic>maximum adjacency ratio</italic> of gene <italic>i</italic> as followswhich is defined if <italic>k<sub>i</sub></italic> = <italic>Σ<sub><italic>j≠i</italic></sub><italic>a<sub>ij</sub></italic>&gt;0. One can easily verify that 0≤<italic>a<sub>ij</sub></italic>≤1 implies 0≤<italic>MAR<sub>i</sub></italic>≤1. Note that <italic>MAR<sub>i</sub></italic> = 1 if all nonzero adjacencies take on their maximum value of 1, which justifies the name “maximum adjacency ratio.” By contrast, if all nonzero adjacencies take on a small (but constant) value <italic>a<sub>ij</sub></italic> = <italic>ε</italic>, then <italic>MAR<sub>i</sub></italic> = <italic>ε</italic> will be small.</italic>
</p> <p>
<italic>Social Network Interpretation of the Maximum Adjacency Ratio: MAR<sub>i</sub></italic> = 1 suggests that the <italic>i</italic>th individual does not form neutral relationships; this individual either strongly likes or dislikes others. In contrast, <italic>MAR<sub>i</sub></italic> = 0.5 suggests the <italic>i</italic>th individual forms less intense relationships with others.</p> <p>
<italic>Potential Uses of the Maximum Adjacency Ratio:</italic> Since <italic>MAR<sub>i</sub></italic> = 1 for all genes in an unweighted network, the maximum adjacency ratio is only useful for weighted networks. The <italic>MAR</italic> can be used to determine whether a hub gene forms moderate relationships with a lot of genes or very strong relationships with relatively few genes. To illustrate this point, we show in the following simple example that the <italic>MAR</italic> can be used to distinguish nodes that have the same connectivity. Assume a network (labeled by <italic>I</italic>) for which the adjacency between node 1 and every other node equals <italic>a</italic>
<sub>1,<italic>j</italic></sub>
<sup>(I)</sup> = 1/(<italic>n</italic>−1). Then <italic>k</italic>
<sub>1</sub>
<sup>(I)</sup> = (<italic>n</italic>−1)/(<italic>n</italic>−1) = 1 and <italic>MAR</italic>
<sub>1</sub>
<sup>(I)</sup> = 1/(<italic>n</italic>−1). For a different network (labeled by II) where <italic>a</italic>
<sub>1,2</sub>
<sup>(II)</sup> = 1 and <italic>a</italic>
<sub>1,<italic>j</italic></sub>
<sup>(II)</sup> = 0 for <italic>j</italic>≥3, the connectivity <italic>k</italic>
<sub>1</sub>
<sup>(II)</sup> still equals 1 but <italic>MAR</italic>
<sub>1</sub>
<sup>(II)</sup> = 1.</p> <p>In weighted coexpression networks, we find empirically that <italic>MAR<sub>i</sub></italic> is often highly correlated with the connectivity <italic>K<sub>i</sub></italic> (see also Equation 36). As we demonstrate in ##FIG##1##Figure 2##, the <italic>MAR<sub>i</sub></italic> is sometimes (but not always) superior to <italic>K<sub>i</sub></italic> when it comes to identifying biologically important intramodular hub genes. As aside, we mention that a directed network analog of <italic>MAR<sub>i</sub></italic> has been used in the analysis of metabolic fluxes ##REF##14985762##[34]##.</p> <title>Network density</title> <p>The <italic>network density</italic> (also known as line density ##UREF##3##[35]##) is defined as the mean off-diagonal adjacency and is closely related to the mean connectivity.where <italic>k</italic> = (<italic>k</italic>
<sub>1</sub>,…,<italic>k<sub>n</sub></italic>) denotes the vector of connectivities and the function vector <italic>v</italic> is defined by <italic>S<sub>p</sub></italic>(<italic>v</italic>) = <italic>Σ<sub><italic>i</italic></sub><italic>v<sub>i</sub><sup>p</sup></italic>.</italic>
</p> <p>
<italic>Social Network Interpretation of the Density:</italic> The density measures the overall affection among individuals. A density close to 1 indicates that all individuals strongly like each other while a density of 0.5 suggests the presence of more ambiguous relationships.</p> <p>
<italic>Potential Uses of the Density:</italic> The density of genes in a subnetwork (e.g., a pathway) can be used to measure whether this sub-network is tight or cohesive. In our motivational mouse tissue example, we find that a network of genes has high density in liver tissue but low density in adipose tissue. The goal of many module detection methods is to find clusters of genes with high density.</p> <title>Network centralization</title> <p>The <italic>network centralization</italic> (also known as degree centralization ##UREF##4##[36]##) is given byThe centralization is 1 for a network with star topology; by contrast, it is 0 for a network where each node has the same connectivity. A regular grid network such as a square has centralization 0.</p> <p>
<italic>Social Network Interpretation of the Centralization:</italic> The centralization of the affection network is close to 1, if one individual has loving relationships with all others who in turn strongly dislike each other. In contrast, a centralization of 0 indicates that all individuals are equally popular.</p> <p>
<italic>Potential Uses of the Centralization:</italic> While the centralization is a widely used measure in social network studies, it has only rarely been used to describe structural differences of metabolic networks ##REF##15603590##[37]##. As described in our motivational example, the centralization can be used to describe properties of cluster trees, see also ##REF##17547772##[8]##.</p> <title>Network heterogeneity</title> <p>The <italic>network heterogeneity</italic> measure is based on the variance of the connectivity. Authors differ on how to scale the variance ##UREF##3##[35]##. We define it as the coefficient of variation of the connectivity distribution, i.e.This heterogeneity measure is invariant with respect to multiplying the connectivity by a scalar.</p> <p>
<italic>Social Network Interpretation of the Heterogeneity:</italic> The heterogeneity can be used to measure the variation of popularity (connectivity) across the individuals.</p> <p>
<italic>Potential Uses of the Heterogeneity:</italic> Describing the reasons for and the meaning of the heterogeneity of complex networks has been the focus of considerable research in recent years ##REF##10935628##[29]##,##REF##16578874##[38]##. Many complex networks have been found to exhibit an approximate scale-free topology, which implies that these networks are very heterogeneous ##REF##14735121##[3]##.</p> <title>Clustering coefficient</title> <p>The <italic>clustering coefficient</italic> of gene <italic>i</italic> is a density measure of local connections, or “cliquishness” ##REF##9623998##[9]##. Specifically,In unweighted networks, <italic>ClusterCoef<sub>i</sub></italic> equals 1 if and only if all neighbors of <italic>i</italic> are also linked to each other. For weighted networks, 0≤<italic>a<sub>ij</sub></italic>≤1 implies that 0≤<italic>ClusterCoef<sub>ij</sub></italic>≤1 ##UREF##0##[19]##.</p> <p>
<italic>Social Network Interpretation of the Clustering Coefficient:</italic> The higher the clustering coefficient of an individual, the higher is the affection among his friends. The clustering coefficient is zero if all of his friends strongly dislike each other.</p> <p>
<italic>Potential Uses of the Clustering Coefficient:</italic> As described in our motivational example, the mean clustering coefficient has been used to measure the extent of module structure present in a network. The relationship between the clustering coefficient and connectivity has been used to describe structural (hierarchical) properties of networks ##REF##12202830##[1]##.</p> <title>Hub gene significance</title> <p>To measure the association between connectivity and gene significance, we propose the following measure of <italic>hub gene significance</italic>:When <italic>GS<sub>i</sub></italic> is proportional to the scaled connectivity (<italic>GS<sub>i</sub></italic> = <italic>cK<sub>i</sub></italic>), the hub gene significance equals the constant of proportionality: <italic>HubGeneSignif</italic> = <italic>c</italic>. The hub gene significance equals the slope of the regression line between <italic>GS<sub>i</sub></italic> and <italic>K<sub>i</sub></italic> if the intercept term is set to 0 (##FIG##2##Figure 3D and 3E##).</p> <p>
<italic>Social Network Interpretation of the Hub Gene Significance:</italic> Assume that the node significance measures the grade point average of the <italic>i</italic>th individual. Then the hub node significance can be used to assess whether there is a relationship between popularity (connectivity) and grade point average.</p> <p>
<italic>Potential Uses of the Hub Gene Significance:</italic> Several studies have shown that the relationship between connectivity and gene significance (i.e., the hub gene significance) carries important biological information. For example, in the analysis of yeast networks, highly connected hub genes were found to be essential for yeast survival and there is evidence that hub genes are preserved across species ##REF##15130938##[17]##, ##REF##17101986##[25]##, ##REF##10935628##[29]##–##REF##15190252##[32]##. A detailed analysis shows that the positive relationship between connectivity and knockout essentiality cannot always be observed ##REF##16515682##[33]##, i.e., the hub gene significance can be close to 0.</p> <title>Network significance measure</title> <p>We define the <italic>network significance measure</italic> as the average gene significance of the genes:
</p> <p>
<italic>Social Network Interpretation of the Network Significance:</italic> The network significance simply measures the average grade point average among the individuals.</p> <p>
<italic>Potential Uses of the Network Significance:</italic> We refer to the network significance of a module network as “module significance.” The module significance measure can be used to address a major goal of gene network analysis: the identification of biologically significant subnetworks or pathways.</p> <title>Centroid significance and centroid conformity</title> <p>We define the <italic>centroid significance</italic> as the gene significance of a suitably chosen representative node (centroid) in the network.where <italic>i.centroid</italic> denotes the index associated with the centroid. A centroid can be defined in many different ways, e.g., based on connectivity or other centrality measures. In our applications, we define the centroid as the most highly connected gene in the network. If multiple genes attain the maximum connectivity, we define the centroid significance by their average gene significance.</p> <p>We define the <italic>centroid conformity</italic> of the <italic>i</italic>th gene as the adjacency between the centroid and the <italic>i</italic>th geneIf multiple genes attain the maximum connectivity, we define the centroid conformity as their average adjacency with the <italic>i</italic>th gene.</p> <p>
<italic>Social Network Interpretation of the Centroid Conformity:</italic> In our affection network, we choose the most popular individual as centroid; then his or her grade point average is the centroid significance. The centroid conformity of the <italic>i</italic>th individual equals his or her affection (connection strength) with the most popular individual.</p> <p>
<italic>Potential Uses of the Centroid Conformity:</italic> Below, we will characterize coexpression networks for which the adjacency <italic>a<sub>ij</sub></italic> can be approximated by a product of the centroid conformities: <italic>a<sub>ij</sub></italic>≈<italic>CentroidConformity<sub>i</sub> CentroidConformity<sub>j</sub></italic>. We will use this insight to derive relationships among seemingly disparate network concepts. For example, the mean clustering coefficient (Equation 12), the density (Equation 9), and the heterogeneity (Equation 11) measure different network properties but we show that they satisfy a simple relationship (Equation 31) in coexpression modules. Further, we will use the centroid significance to derive a simple relationship (Equation 37) between module significance (Equation 14) and hub gene significance (Equation 13).</p> <title>Overview of Weighted Gene Coexpression Network Analysis</title> <p>One of the many biological applications of gene coexpression networks is the identification of pathways (modules) and centrally located genes (referred to as module centroids). In our applications, we define highly connected intramodular hub genes as module centroids. Weighted gene coexpression network analysis (WGCNA, ##UREF##0##[19]##,##REF##17090670##[24]##) can be considered a step-wise microarray data reduction technique, which starts from the level of thousands of genes, identifies clinically interesting gene modules, and finally represents the modules by their centroids. The module centric analysis alleviates the multiple testing problem inherent in microarray data analysis. Instead of relating thousands of genes to a sample trait, it focuses on the relationship between a few (usually less than 10) modules and the sample trait.</p> <p>An outline of WGCNA is presented in ##FIG##2##Figure 3A##. The module definition does not make use of <italic>a priori</italic> defined gene sets. Instead, modules are constructed from the expression data by using a tight clustering procedure. Although it is advisable to relate the resulting modules to gene ontology information to assess their biological plausibility, it is not required. Because the modules may correspond to biological pathways, focusing the analysis on modules (and corresponding centroids) amounts to a biologically motivated data reduction method. Intramodular hub genes are centrally located in the module and thus lend themselves as candidates for biomarkers. Examples of biological studies that show the importance of intramodular hub genes can be found reported in ##UREF##1##[23]##–##REF##17101986##[25]##,##REF##16515682##[33]##,##REF##16912112##[39]##. Because the expression profiles of intramodular hub genes are highly correlated (in our data, <italic>r</italic>&gt;0.90), typically dozens of candidates result. Although these candidates are statistically equivalent, they may differ in terms of biological plausibility or clinical utility.</p> <title>Network Modules</title> <p>Roughly speaking, we define network modules as groups of highly interconnected genes. As detailed in ##SUPPL##0##Text S1##, ##SUPPL##1##Text S2##, ##SUPPL##2##Text S3##, and in our online R tutorials, we use a hierarchical clustering procedure to identify modules (clusters) as branches of the resulting cluster tree. A common but inflexible branch cutting method uses a constant height cutoff value. Alternatively, dynamic branch cutting adaptively chooses cutting values depending on the shape of the branch ##REF##18024473##[40]##. Each module is assigned a unique color label (##FIG##2##Figure 3B##). Our branch cutting algorithm only assigns module colors to branches whose size exceeds a user-specified threshold parameter. In practice, it is advisable to vary the minimum module size and other branch cutting parameters to determine how the results are affected by different parameter choices. An iterative approach for choosing the parameters could be defined by optimizing the module significance. This module detection approach has led to biologically meaningful modules in several applications ##REF##12202830##[1]##, ##REF##17547772##[8]##, ##UREF##1##[23]##–##REF##17101986##[25]##, ##REF##16515682##[33]##, ##REF##16912112##[39]##–##REF##17110366##[43]## but our theoretical results transcend this particular module detection method. Any module detection method that results in clusters of highly correlated gene expressions could be used.</p> <title>Intramodular Network Concepts</title> <p>In the following, we assume that a module detection method (e.g., a clustering procedure) has found <italic>Q</italic> modules. We denote the adjacency matrix of the genes inside the <italic>q</italic>th module by <italic>A</italic>
<sup>(<italic>q</italic>)</sup>. Thus, <italic>A</italic>
<sup>(<italic>q</italic>)</sup> represents a subnetwork comprised of the genes in the <italic>q</italic>th module. Analogously, we define <italic>GS</italic>
<sup>(<italic>q</italic>)</sup> as the gene significance measure restricted to the module genes. Denote by <italic>n</italic>
<sup>(<italic>q</italic>)</sup> the number of genes inside the <italic>q</italic>th module. Throughout the manuscript, we use the superscript (<italic>q</italic>) to denote quantities associated with the <italic>q</italic>th module. But for notational convenience, we sometimes omit (<italic>q</italic>) when the context is clear.</p> <p>We define an intramodular network concept <italic>NCF</italic>(<italic>A</italic>
<sup>(<italic>q</italic>)</sup>,<italic>GS</italic>
<sup>(<italic>q</italic>)</sup>) by evaluating a network concept function <italic>NCF</italic>(·,·) on the adjacency matrix <italic>A</italic>
<sup>(<italic>q</italic>)</sup> and/or a corresponding gene significance measure <italic>GS</italic>
<sup>(<italic>q</italic>)</sup>.</p> <p>For example, the intramodular connectivity is defined bywhere the <italic>j</italic> indexes the genes in the <italic>q</italic>th module. Intramodular connectivity has been found to be an important complementary gene screening variable for finding biologically important genes ##REF##17090670##[24]##,##REF##17101986##[25]##,##REF##16912112##[39]##.</p> <p>We refer to the network significance (Equation 14) of a module network simply as the <bold>module significance measure</bold>, i.e., the module significance is the average gene significance of the module genes:
</p> <title>Data Reduction Methods for Microarray Data</title> <p>The high dimensionality of gene expression data has inspired two broad categories of data reduction techniques. The first category, often used by network theorists, is to reduce the gene coexpression networks into modules. Each module can be represented by a centroid, e.g., an intramodular hub gene. The second category, often used by microarray data analysts, reduces the gene expression data to a small number of components that capture the essential behavior of the expression profiles ##REF##16926220##[27]##, ##REF##10963673##[44]##–##REF##17389406##[51]##. One of our goals is to understand how the two categories of data reduction methods relate to each other. Here we use the singular value decomposition ##REF##10963673##[44]##,##REF##10890920##[45]##,##REF##11983907##[48]## since this will allow us to define a simple measure of factorizability (Equation 24).</p> <title>Singular value decomposition</title> <p>For the <italic>q</italic>th module, denote by <italic>X</italic>
<sup>(<italic>q</italic>)</sup> the <italic>n</italic>
<sup>(<italic>q</italic>)</sup>×<italic>m</italic> matrix of <italic>n</italic>
<sup>(<italic>q</italic>)</sup> gene expression profiles across <italic>m</italic> microarrays:where <italic>x<sub>i</sub></italic> denotes the gene expression vector of the <italic>i</italic>th gene.</p> <p>The singular value decomposition (SVD) of <italic>X</italic>
<sup>(<italic>q</italic>)</sup> is given by <italic>X</italic>
<sup>(<italic>q</italic>)</sup> = <italic>U</italic>
<sup>(<italic>q</italic>)</sup>
<italic>D</italic>
<sup>(<italic>q</italic>)</sup>(<italic>V</italic>
<sup>(<italic>q</italic>)</sup>)<italic><sup>T</sup></italic>, where <italic>U</italic>
<sup>(<italic>q</italic>)</sup> is an <italic>n</italic>
<sup>(<italic>q</italic>)</sup>×<italic>m</italic> matrix with orthonormal columns, <italic>V</italic>
<sup>(<italic>q</italic>)</sup> is an <italic>m</italic>×<italic>m</italic> orthogonal matrix, and <italic>D</italic>
<sup>(<italic>q</italic>)</sup> is an <italic>m</italic>×<italic>m</italic> diagonal matrix of the singular values {|<italic>d<sub>l</sub></italic>
<sup>(<italic>q</italic>)</sup>|}. Specifically, <italic>V</italic>
<sup>(<italic>q</italic>)</sup> and <italic>D</italic>
<sup>(<italic>q</italic>)</sup> are given by
</p> <p>The singular value decomposition of <italic>X</italic>
<sup>(<italic>q</italic>)</sup> is closely related to the principal component analysis of the correlation matrix <italic>COR</italic> = [cor(<italic>x<sub>i</sub></italic>
<sup>(<italic>q</italic>)</sup>,<italic>x<sub>j</sub></italic>
<sup>(<italic>q</italic>)</sup>)] whose entries correspond to the pairwise correlations between the rows (genes) of <italic>X</italic>
<sup>(<italic>q</italic>)</sup>. For example, the eigenvalues of the correlation matrix <italic>COR</italic> are squares of corresponding singular values |<italic>d<sub>l</sub></italic>
<sup>(<italic>q</italic>)</sup>|.</p> <p>We assume that the singular values |<italic>d<sub>l</sub></italic>
<sup>(<italic>q</italic>)</sup>| are arranged in decreasing order. Adapting terminology from ##REF##10963673##[44]##, we refer to the first column of <italic>V</italic>
<sup>(<italic>q</italic>)</sup> as the <bold>Module Eigengene</bold>:
</p> <p>For brevity, we sometimes drop the superscript (<italic>q</italic>) and simply refer to <italic>E</italic> as the eigengene. The module eigengene can be used to summarize and represent the expression profiles of the module genes, see ##FIG##3##Figure 4B##. The proportion of variance explained by the module eigengene <italic>E</italic>
<sup>(<italic>q</italic>)</sup> is defined as
</p> <title>High Level View of Gene Coexpression Networks and Eigengene Networks</title> <p>The module eigengenes of different modules can be highly correlated (##FIG##3##Figure 4A##). Detecting a high correlation between module eigengenes may either be of biological interest (suggesting interactions between pathways) or it may be a methodological artifact (suggesting poorly defined modules that should be merged). The correlations between two eigengenes can be used to define eigengene coexpression networks ##REF##18031580##[52]##, e.g., a weighted eigengene coexpression network can be defined as followswhere <italic>E</italic>
<sup>(<italic>q</italic>)</sup> and <italic>E</italic>
<sup>(<italic>p</italic>)</sup> represent the eigengenes of two distinct modules. Apart from correlating the module eigengenes of different modules to each other, one can relate the module eigengenes to an external microarray sample trait <italic>T</italic> to identify trait related modules. Thus, eigengene network analysis can be viewed as a network reduction scheme that reduces a gene coexpression network involving thousands of genes to an orders of magnitude smaller metanetwork involving module representatives (one eigengene per module).</p> <p>Unlike traditional microarray data reduction methods that impose orthogonality (e.g., principal component analysis) or independence (e.g., independent component analysis), gene coexpression network analysis can be considered a pathway-based data reduction method that allows dependencies between the modules. When focusing on the use of module eigengenes, network analysis can be considered a variant of oblique factor analysis.</p> <title>Low Level View of a Single Module and Factorizable Networks</title> <p>While a high level view of modular gene coexpression networks can be viewed as a data reduction technique, many network analyses focus on the pairwise relationships of relatively few (hundreds) of correlated genes, i.e., genes that form a single module in a larger network. For example, the 498 genes of our motivational example were part of a body weight related module, which was found in a large gene coexpression network based on the female mouse liver samples ##UREF##1##[23]##.</p> <p>The low-level analysis of a single network module may help identify key genes that may be used as therapeutic targets or candidate biomarkers. An important question of low level analysis is to efficiently describe the connection strengths between interacting module genes. We have provided <italic>empirical</italic> evidence that many module adjacency matrices, i.e., networks comprised of genes of a single module, are approximately factorizable ##REF##17547772##[8]##. In such networks, the adjacency between module genes <italic>i</italic> and <italic>j</italic> can approximately be factored into gene specific contributions, i.e., <italic>a<sub>ij</sub></italic>
<sup>(<italic>q</italic>)</sup>≈<italic>CF<sub>i</sub></italic>
<sup>(<italic>q</italic>)</sup>
<italic>CF<sub>j</sub></italic>
<sup>(<italic>q</italic>)</sup> with <italic>CF<sub>i</sub></italic>
<sup>(<italic>q</italic>)</sup> defined as the conformity of gene <italic>i</italic>. Thus, the adjacency matrix of an approximately factorizable network can be approximated using the rank 1 matrix [<italic>CF<sub>i</sub></italic>
<sup>(<italic>q</italic>)</sup>
<italic>CF<sub>j</sub></italic>
<sup>(<italic>q</italic>)</sup>]. The conformity vector <italic>CF</italic>
<sup>(<italic>q</italic>)</sup> can be estimated in several ways ##REF##17547772##[8]##; it is highly related to a single factor nonnegative matrix decomposition of <italic>A</italic>
<sup>(<italic>q</italic>)</sup>
##REF##17389406##[51]## and it is highly related to the connectivity .</p> <title>Characterizing approximately factorizable coexpression modules</title> <p>An open theoretical research question is to characterize microarray data that lead to factorizable coexpression networks. Here we solve this problem for the case of modules in a gene coexpression network. Toward this end, we propose the following measure of <bold>eigengene factorizability</bold>:Note that 0≤<italic>EF</italic>(<italic>X</italic>
<sup>(<italic>q</italic>)</sup>)≤1 and the close resemblance to the proportion of variance explained by the module eigengene (Equation 22). In the <xref ref-type="sec" rid="s4">Methods</xref> section, we argue that <italic>EF</italic>(<italic>X</italic>
<sup>(<italic>q</italic>)</sup>)≈1 implies that the correlation matrix factors as followsFurther, we derive the following</p> <title>Observation 1</title> <p>
<italic>If the eigengene factorizability EF(X<sup>(q)</sup>) is close to 1, the adjacencies of the weighted coexpression module network A<sup>(q)</sup> = |cor(X<sup>(q)</sup>)|<sup>β</sup> and the trait-based gene significance measure GS<sub>i</sub><sup>(q)</sup> = |cor(x<sub>i</sub><sup>(q)</sup>,T)|<sup>β</sup> can be factored as follows</italic>

<italic>where</italic>

<italic>is referred to as the </italic>
<bold><italic>eigengene conformity</italic></bold>
<italic> of the ith gene, and</italic>

<italic>is referred to as the qth module </italic>
<bold><italic>eigengene significance</italic></bold>
<italic> with respect to T, also denoted as EigengeneSignif<sup>(q)</sup>.</italic>
</p> <p>As described in ##TAB##0##Table 1##, the eigengene significance and the eigengene conformity are the eigengene-based counterparts of the centroid significance (Equation 15) and centroid conformity (Equation 16), respectively.</p> <p>The eigengene-based approximations on the right hand side of Equation 25 motivate us to define the eigengene-based adjacency matrix <italic>A<sub>E</sub></italic>
<sup>(<italic>q</italic>)</sup> and gene significance measure <italic>GS<sub>E</sub></italic>
<sup>(<italic>q</italic>)</sup> as follows:
For our coexpression modules, we find empirically that the eigengene factorizability is close to 1 (see ##TAB##1##Table 2##, ##SUPPL##0##Text S1##, ##SUPPL##1##Text S2##, and ##SUPPL##2##Text S3##).</p> <p>Abstractly speaking, Observation 1 allows us to characterize coexpression networks for which the adjacency <italic>a<sub>ij</sub></italic> can be approximated by a product of the centroid conformities (Equation 16): <italic>a<sub>ij</sub></italic>≈<italic>CentroidConformity<sub>i</sub> CentroidConformity<sub>j</sub></italic>.</p> <title>Geometric interpretation of factorizability</title> <p>In the <xref ref-type="sec" rid="s4">Methods</xref> section, we argue that <italic>EF</italic>(<italic>X</italic>
<sup>(<italic>q</italic>)</sup>)≈1 if the module gene expressions <italic>x<sub>i</sub></italic>
<sup>(<italic>q</italic>)</sup> are approximately orthogonal to the right singular vectors <italic>v<sub>l</sub></italic>
<sup>(<italic>q</italic>)</sup> for <italic>l</italic>≥2, i.e., if on average the gene expression profiles point in the direction of the module eigengene <italic>v</italic>
<sub>1</sub>
<sup>(<italic>q</italic>)</sup> = <italic>E</italic>
<sup>(<italic>q</italic>)</sup>. A rough geometric intuition of <italic>a<sub>ij</sub></italic>
<sup>(<italic>q</italic>)</sup>≈<italic>a<sub>e,j</sub></italic>
<sup>(<italic>q</italic>)</sup>
<italic>a<sub>e,j</sub></italic>
<sup>(<italic>q</italic>)</sup> (Equation 25) is presented in ##FIG##4##Figure 5A##. The angle between the module eigengene <italic>E</italic>
<sup>(<italic>q</italic>)</sup> and the <italic>i</italic>th gene expression profile is denoted by <italic>θ<sub>i</sub></italic>. The angle between gene expression profiles <italic>i</italic> and <italic>j</italic> is denoted by <italic>θ<sub>ij</sub></italic>. In the <xref ref-type="sec" rid="s4">Methods</xref> section, we show that <italic>θ<sub>ij</sub></italic>≈|<italic>θ<sub>i</sub></italic>±<italic>θ<sub>j</sub></italic>| and sin(<italic>θ<sub>i</sub></italic>) sin(<italic>θ<sub>j</sub></italic>)≈0 imply approximate factorizability of the correlation matrix.</p> <title>Eigengene-Based Analogs of Network Concepts</title> <p>Here we define <bold>eigengene-based network concepts</bold> as a step towards a geometric interpretation of network concepts. Analogous to the case of intramodular network concepts, we define eigengene-based network concepts by evaluating the network concept function <italic>NCF</italic>(<italic>A<sub>E</sub></italic>
<sup>(<italic>q</italic>)</sup>,<italic>GS<sub>E</sub></italic>
<sup>(<italic>q</italic>)</sup>) on the eigengene-based adjacency matrix <italic>A<sub>E</sub></italic>
<sup>(<italic>q</italic>)</sup> (Equation 28) and the eigengene-based gene significance measure <italic>GS<sub>E</sub></italic>
<sup>(<italic>q</italic>)</sup> (Equation 29). One can easily derive the following formulas for eigengene-based network concepts:where . Under the assumptions of Observation 1, we find that <italic>A</italic>
<sup>(<italic>q</italic>)</sup>≈<italic>A<sub>E</sub></italic>
<sup>(<italic>q</italic>)</sup> and <italic>GS<sub>i</sub></italic>≈<italic>GS<sub>E,i</sub></italic>. For a continuous network concept function <italic>NCF</italic>(·,·) this implies <italic>NCF</italic>(<italic>A</italic>
<sup>(<italic>q</italic>)</sup>,<italic>GS</italic>)≈<italic>NCF</italic>(<italic>A<sub>E</sub></italic>
<sup>(<italic>q</italic>)</sup>,<italic>GS<sub>E</sub></italic>). We summarize this observation as follows</p> <title>Observation 2</title> <p>
<italic>If A<sup>(q)</sup> = |cor(X<sup>(q)</sup>)|<sup>β</sup> and the eigengene factorizability EF(X<sup>(q)</sup>) is close to 1, the network concepts can be approximated by their eigengene-based analogs.</italic>
</p> <p>This observation is illustrated in ##FIG##5##Figure 6##.</p> <title>Using the eigengene-based heterogeneity to study the effect of soft thresholding</title> <p>It can be advantageous to replace network concepts by their eigengene-based analogs when studying theoretical properties. To illustrate this point, we briefly describe the effect of soft thresholding <italic>a<sub>ij</sub></italic> = <italic>s<sub>ij</sub><sup>β</sup></italic> (Equation 2) on the network heterogeneity. Using extensive simulation studies reported on our webpage, we found that for the vast majority of networks, the heterogeneity increases with the soft threshold <italic>β</italic>. Thus, for most coexpression networks, increasing <italic>β</italic> makes it easier to discern highly connected genes from less connected genes. However, one can construct networks for which increasing <italic>β</italic> leads to a lower heterogeneity. The situation is much simpler for the eigengene-based heterogeneity <italic>Heterogeneity<sub>E</sub></italic>
<sup>(<italic>q</italic>)</sup> (Equation 30). In the <xref ref-type="sec" rid="s4">Methods</xref> section, we prove that the eigengene-based heterogeneity is a monotonically increasing function of the soft threshold <italic>β</italic>. Thus, the heterogeneity will be an increasing function of <italic>β</italic> if it can be approximated by its eigengene based analog (Observation 2).</p> <title>Relationships among eigengene-based network concepts</title> <p>A major theoretical advantage of eigengene-based network concepts is that they reveal simple relationships amongst each other. For example, it is straightforward to derive
</p> <p>To arrive at particular simple relationships among network concepts, we make use of the following terminology. We denote the <bold>maximum eigengene conformity</bold> as <italic>a<sub>e</sub></italic>
<sub>,max</sub>
<sup>(<italic>q</italic>)</sup> = max<italic><sub>j</sub></italic>(<italic>a<sub>e,j</sub></italic>
<sup>(<italic>q</italic>)</sup>), where <italic>a<sub>e,j</sub></italic>
<sup>(<italic>q</italic>)</sup> = |cor(<italic>x<sub>j</sub></italic>
<sup>(<italic>q</italic>)</sup>,<italic>E</italic>
<sup>(<italic>q</italic>)</sup>)|<italic><sup>β</sup></italic> (Equation 26). In most modules, we find genes that have very high correlations (<italic>r</italic>≈<italic>0.99</italic>) with the module eigengene. For a low power <italic>β</italic>, this implies that the maximum eigengene conformity is approximately equal to 1:
</p> <p>We refer to Equation 32 as the <bold>maximum conformity assumption</bold>. With the results in the <xref ref-type="sec" rid="s4">Methods</xref> section, one can show that the maximum conformity assumption implies the following</p> <title>Observation 3</title> <p>
<italic>If A<sup>(q)</sup> = |cor(X<sup>(q)</sup>)|<sup>β</sup>, EF(X<sup>(q)</sup>)≈1 and the maximum conformity assumption applies, intramodular network concepts satisfy the following relationships</italic>








<italic>where mean(ClusterCoef<sup>(q)</sup>) denotes the mean clustering coefficient, ClusterCoef<sub>max</sub><sup>(q)</sup> = max<sub>j</sub>(ClusterCoef<sub>j</sub><sup>(q)</sup>) and MAR<sub>max</sub><sup>(q)</sup> = max<sub>j</sub>(MAR<sub>j</sub><sup>(q)</sup>).</italic>
</p> <p>In practice, we find that the maximum conformity assumption holds well for low values of <italic>β</italic>. Below, we study the robustness of our results with respect to higher powers and alternative network construction methods.</p> <title>Geometric interpretation of network concepts</title> <p>Observations 2 and 3 allow us to provide a geometric interpretation of intramodular network concepts.</p> <p>The relationship between the scaled intramodular connectivity <italic>K<sub>i</sub></italic>
<sup>(<italic>q</italic>)</sup> and its eigengene based analog <italic>a<sub>e,i</sub></italic>
<sup>(<italic>q</italic>)</sup> = |cor(<italic>x<sub>i</sub></italic>
<sup>(<italic>q</italic>)</sup>,<italic>E</italic>
<sup>(<italic>q</italic>)</sup>)|<italic><sup>β</sup></italic> (Equation 33) facilitates a geometric interpretation of the intramodular connectivity: the smaller the angle <italic>θ<sub>i</sub></italic> between the <italic>i</italic>th gene expression profile and the module eigengene, the larger is |cos(<italic>θ<sub>i</sub></italic>)|<italic><sup>β</sup></italic> = <italic>a<sub>e,i</sub></italic>
<sup>(<italic>q</italic>)</sup>, i.e., the larger is the scaled intramodular connectivity. Since the module eigengene summarizes the overall behavior of the module, <italic>a<sub>e,i</sub></italic>
<sup>(<italic>q</italic>)</sup> measures how well gene <italic>i</italic> conforms to the overall module. Thus, a tongue-in-cheek social network interpretation of Equation 33 is that group-conforming behavior leads to high popularity.</p> <p>We provide two geometric interpretations of the density. The first makes use of the relationship <italic>a<sub>ij</sub></italic>
<sup>(<italic>q</italic>)</sup> = |cos(<italic>θ<sub>ij</sub></italic>)|<italic><sup>β</sup></italic> where <italic>θ<sub>ij</sub></italic> denotes the angle between gene expression profiles <italic>i</italic> and <italic>j</italic>. By definition (Equation 9), the smaller the pairwise angles <italic>θ<sub>ij</sub></italic> between the gene expression profiles, the higher is the module density. Equation 39 provides another interpretation: the smaller the angles <italic>θ<sub>i</sub></italic> between the module gene expression profiles and the module eigengene, the higher is the density. Thus, the density can be interpreted as a measure of average closeness between the gene expression profiles and the module eigengene. By definition, coexpression module networks have a relatively high density (see ##TAB##1##Table 2##, ##SUPPL##0##Text S1##, ##SUPPL##1##Text S2##, and ##SUPPL##2##Text S3##).</p> <p>The eigengene-based heterogeneity equals the coefficient of variation of the <italic>a<sub>E</sub></italic>
<sup>(<italic>q</italic>)</sup>, i.e., it is a measure of variability of the angles <italic>θ<sub>i</sub></italic> between the gene expression profiles and the module eigengene. The heterogeneity equals 0 if the angles <italic>θ<sub>i</sub></italic> are all equal.</p> <p>The <italic>i</italic>th gene has high eigengene-based significance <italic>GS<sub>E,i</sub></italic>
<sup>(<italic>q</italic>)</sup> (Equation 29) if the eigengene has a small angle with the sample trait and <italic>θ<sub>i</sub></italic> is small. Similarly, the geometric interpretation of the hub gene significance (Equation 13) is straightforward: the smaller the angle between the module eigengene and the sample trait, the higher is the hub gene significance (Equation 34).</p> <p>We provide two geometric interpretations of the module significance (Equation 14). The first interpretation is based on the definition of the module significance as average gene significance; a module has high module significance if on average the angles between the module expression profiles and the sample trait tend to be small. The second interpretation of the module significance is based on Equation 37: a module has high significance if the module density is high and the angle between the module eigengene and the sample trait is small.</p> <title>What Can Microarray Data Analysts Learn from the Geometric Interpretation?</title> <p>Here we illustrate how the geometric interpretation of gene coexpression networks can be used to derive results, which may be interesting to microarray data analysts.</p> <title>Summarizing the expression profiles of a module</title> <p>Multiple approaches are conceivable for summarizing the expression profiles of the genes inside a single module. One approach (popular with statisticians) applies a singular value decomposition to the expression data and summarizes the module with the module eigengene. Another approach (popular with network theorists) is to construct a module network and to use the most highly connected hub gene as centroid. Since Equation 33 implies that hub genes are highly correlated with the module eigengene, we find that the two seemingly different approaches will lead to very similar results in practice (##FIG##3##Figure 4C##).</p> <title>Intramodular connectivity is a measure of module membership</title> <p>Since module construction is computationally intensive, one often restricts the module detection analysis to a subset of the original genes on the microarray, e.g., the most varying and/or the most connected genes. To counter this loss of information, generalizing the intramodular connectivity to extramodular genes, i.e. genes outside the module, is an important problem. Our solution is motivated by the relationship between the intramodular connectivity and its eigengene based analog (Equation 33). Specifically, the <italic>q</italic>th module eigengene gives rise to an eigengene-based scaled intramodular connectivity measure
</p> <p>Under the assumptions of Observation 3, Equation 33 implies that <italic>K<sub>i</sub></italic>
<sup>(<italic>q</italic>)</sup>≈|<italic>K</italic>
<sub>cor,<italic>i</italic></sub>
<sup>(<italic>q</italic>)</sup>|<italic><sup>β</sup></italic> for the subset of genes that are in the <italic>q</italic>th module. The larger <italic>K</italic>
<sub>cor,<italic>i</italic></sub>
<sup>(<italic>q</italic>)</sup>, the more similar is gene <italic>i</italic> is to the summary profile of the <italic>q</italic>th module. Thus, <italic>K</italic>
<sub>cor,<italic>i</italic></sub>
<sup>(<italic>q</italic>)</sup> can be used to measure module membership. A theoretical advantage of <italic>K</italic>
<sub>cor,<italic>i</italic></sub>
<sup>(<italic>q</italic>)</sup> over <italic>K<sub>i</sub></italic>
<sup>(<italic>q</italic>)</sup> is that its definition can be easily extended to expression profiles <italic>x<sub>i</sub></italic> outside the <italic>q</italic>th module. Another advantage of <italic>K</italic>
<sub>cor,<italic>i</italic></sub>
<sup>(<italic>q</italic>)</sup> is that a simple correlation test <italic>p</italic>-value can be used to assess the statistical significance of the correlation between <italic>x<sub>i</sub></italic> and <italic>E</italic>
<sup>(<italic>q</italic>)</sup>.</p> <title>Fuzzy module annotation of genes</title> <p>Module detection usually involves certain parameter choices. For some genes, it may be difficult to decide whether they belong to a particular module or whether they belong to more than one module. Instead of reporting a binary indicator of module membership, it can be advantageous to report a fuzzy measure of module membership, which takes on values in the unit interval [0,1]. A natural choice for a fuzzy measure of module membership is the eigengene-based scaled intramodular connectivity measure <italic>K</italic>
<sub>cor,<italic>i</italic></sub>
<sup>(<italic>q</italic>)</sup> (Equation 41). The fuzzy module membership measures <italic>K</italic>
<sub>cor,<italic>i</italic></sub>
<sup>(<italic>q</italic>)</sup> specify how close gene <italic>i</italic> is to modules <italic>q</italic> = 1,…,<italic>Q</italic>. It is straightforward to use these measures for finding genes that are close to two modules, i.e., intermediate genes. In ##FIG##6##Figure 7##, we show the pairwise relationships among different <italic>K</italic>
<sub>cor,<italic>i</italic></sub>
<sup>(<italic>q</italic>)</sup> measures where the genes are colored by their original module assignment. Note that many of the nonmodule (grey) genes lie intermediate between the proper module genes.</p> <title>What Can Network Theorists Learn from the Geometric Interpretation?</title> <p>In the following, we provide several examples that illustrate potential uses of the geometric interpretation.</p> <title>Statistical significance of network concepts</title> <p>While fundamental network concepts are defined as functions of the network adjacency matrix, their eigengene-based analogs are often simple monotonic functions of correlation coefficients. This insight can be used to attach significance levels (<italic>p</italic>-values) to several eigengene-based network concepts. For example, the eigengene-based hub gene significance is a monotonic function of the correlation between the eigengene and the sample trait (Equation 34). Thus, one can use a correlation test <italic>p</italic>-value ##UREF##7##[53]## or a regression-based <italic>p</italic>-value for assessing the statistical significance between <italic>E</italic>
<sup>(<italic>q</italic>)</sup> and the sample trait <italic>T</italic>. Analogously, one can attach a significance level to the fuzzy module membership measures <italic>K</italic>
<sub>cor,<italic>i</italic></sub>
<sup>(<italic>q</italic>)</sup> (Equation 41).</p> <p>Since the gene coexpression network concepts are based on correlations between quantitative variables, one can use permutation test procedures to attach significance levels to network concepts. By randomly permuting the gene expression values of each gene, it is possible to noise up the correlation structure inherent in the original data. We find that the resulting permuted data lead to networks with low density and low mean clustering coefficients (reflecting the lack of large modules).</p> <title>Relationship between centralization and density</title> <p>The relationship between centralization and density (Equation 40) is surprisingly simple for coexpression networks but it does not hold in general networks. For a general network, one can only derive an upper bound for the centralization in terms of the density ##UREF##3##[35]##. As a caveat, we mention that our empirical studies (described below) show that Equation 40 is not very robust with regard to deviations from our theoretical assumptions.</p> <title>Intramodular hub genes cannot be intermediate genes in coexpression networks</title> <p>The geometric interpretation of gene coexpression network analysis can be used to argue that a gene that lies “intermediate” between two distinct modules cannot be a highly connected intramodular hub gene in either module (see ##FIG##4##Figure 5B##). More precisely, we refer to gene <italic>i</italic> as hub gene in module 1 if its scaled connectivity <italic>K<sub>i</sub></italic>
<sup>(1)</sup> is very high (say larger than 0.9). Further, we refer to two modules as distinct if their respective eigengenes have a low correlation, say |cor(<italic>E</italic>
<sup>(1)</sup>,<italic>E</italic>
<sup>(2)</sup>)|&lt;0.3. We refer to gene <italic>i</italic> as intermediate between modules 1 and 2 if it has a moderately high connectivity with both modules, say <italic>K<sub>i</sub></italic>
<sup>(1)</sup>&gt;0.5 and <italic>K<sub>i</sub></italic>
<sup>(2)</sup>&gt;0.5.</p> <p>Equation 33 allows us to translate statements about the scaled intramodular connectivity into statements about the angles between genes and module eigengenes. A gene is an intermediate gene if it has a moderately small angle with both module eigengenes. If the eigengenes are distinct (i.e., the angle between them is large), the intermediate gene cannot have a very small angle with either module eigengene, i.e., it cannot be an intramodular hub gene in either module. A geometric interpretation of this example can be found in ##FIG##4##Figure 5B##.</p> <p>As an important caveat, we mention that intermediate network genes may well be highly connected “hub” genes if the factorizability property does not hold such as in the entire network comprised of multiple distinct modules.</p> <title>Characterizing module networks where hub genes are significant</title> <p>For a trait-based gene significance measure, the striking relationship between module significance and hub gene significance (Equation 37) suggests a positive relationship between connectivity and gene significance (high hub gene significance) in modules that are enriched with significant genes (high module significance).</p> <p>Further, Equation 34 shows that the hub gene significance of a module network is determined by the angle between the module eigengene and the sample trait. This allows us to describe situations when a module has high hub gene significance, i.e., when there is a strong positive relationship between gene significance and intramodular connectivity. In the example provided in ##FIG##4##Figure 5C and 5D##, the angle between <italic>E</italic> and <italic>T2</italic> is small which implies that the hub gene significance with regard to <italic>GS2<sub>i</sub></italic> = |cor(<italic>x<sub>i</sub></italic>,<italic>T2</italic>)| is high. By contrast, the angle between <italic>E</italic> and <italic>T1</italic> is large, which implies that the hub gene significance with regard to <italic>GS1<sub>i</sub></italic> = |cor(<italic>x<sub>i</sub></italic>,<italic>T1</italic>)| is low.</p> <title>Dictionary for Translating between Network Concepts and Their Eigengene-Based Analogs</title> <p>To facilitate the communication between microarray data analysts and network theorists, we provide a short dictionary for translating between microarray data analysis and network theory terminology. More specifically, for a subset (module) of genes that have high expression factorizability, ##TAB##0##Table 1## describes the correspondence between general network terms and their eigengene-based counterparts. While our theoretical derivations assume a weighted gene coexpression network, our robustness studies show empirically that many of the findings apply to unweighted networks as well. The summary of empirical robustness studies is described below.</p> <p>In general, eigengene-based concepts are no substitute for network concepts. It is natural to use network concepts when describing the pairwise relationships between genes and to use eigengene-based network concepts when relating the gene expression profiles to a module eigengene. Since eigengene-based network concepts tend to be relatively simple, they often simplify theoretical derivations. Further, many of them allow one to calculate a statistical significance level (<italic>p</italic>-value) using a correlation or regression based test statistic.</p> <title>Real Data Applications</title> <p>To illustrate the theoretical results we report 4 different microarray data applications. The underlying data sets and R software code can be found on our webpage <ext-link ext-link-type="uri" xlink:href="http://www.genetics.ucla.edu/labs/horvath/ModuleConformity/GeometricInterpretation/">http://www.genetics.ucla.edu/labs/horvath/ModuleConformity/GeometricInterpretation/</ext-link>.</p> <title>Brain cancer network application</title> <p>Here we describe a weighted gene coexpression network that was constructed on the basis of 55 microarray samples of glioblastoma (brain cancer) patients. A detailed description of the data, modules, and biological implications can be found in ##REF##17090670##[24]##. We defined 6 modules as branches of an average linkage hierarchical cluster tree (##FIG##2##Figure 3B##). Module membership in the 6 “proper” modules is color-coded by turquoise, blue, brown, yellow, green and red. Grey denotes the color of genes that were not grouped into any of the 6 proper modules. To allow for a comparison, we also report results for the “improper” module comprised of grey genes.</p> <p>We used the patient survival time as microarray sample trait <italic>T</italic>. We defined a gene significance measure as the absolute value of the correlation between <italic>T</italic> and the gene expression profiles (Equation 4). The module significance was defined as average gene significance (Equation 14). ##FIG##2##Figure 3C## shows that the brown module had the highest module significance. This module was previously found to be enriched with genes that are prognostic of patient survival ##REF##17090670##[24]##.</p> <p>By relating the gene significance measure <italic>GS<sub>i</sub></italic> to the scaled connectivity <italic>K<sub>i</sub></italic>, we arrive at a hub gene significance measure (Equation 13). As illustrated in ##FIG##2##Figure 3D and 3E##, the hub gene significance is defined as the slope of a regression model without intercept term. The brown module had the highest hub gene significance, see ##TAB##1##Table 2##.</p> <p>We defined the module eigengene significance (Equation 27) as the absolute value of the correlation between the module eigengene and patient survival time. The brown module eigengene also had the highest eigengene significance: <italic>a<sub>e,t</sub></italic>
<sup>brown</sup> = |cor(<italic>E</italic>
<sup>brown</sup>,<italic>T</italic>)| = 0.202. An advantage of the eigengene-based hub gene significance (the eigengene significance) is that it allows one to compute a corresponding p-value. Using a correlation test, we find that the value of the eigengene significance <italic>a<sub>e,t</sub></italic>
<sup>brown</sup> is statistically insignificant (<italic>p</italic> = 0.30) in this dataset. However, when we combined these data with an additional data set, we found that the brown module eigengene is significantly related to survival time ##REF##17090670##[24]##.</p> <p>We visualize the gene expression profiles of module genes with a heat map plot (##FIG##3##Figure 4B##) where rows correspond to the genes, the columns to the samples, and the gene expression profiles have been standardized to a mean of 0 and a variance of 1. The heat map colors high and low expression values by red and green, respectively. For a given module, the heat map exhibits characteristic vertical bands that reflect the high correlation among module gene expression profiles. For the 6 proper modules of our brain cancer application, the proportion of variance explained by the first eigengene ranges from 0.59 to 0.71 (##TAB##1##Table 2##). For the improper grey module genes (defined as genes outside of all proper modules) the proportion of variance explained by the first eigengene is only 0.28. Similarly, when all network genes are used to define an improper module, the proportion of variance explained by the first eigengene is only 0.32. As expected by module construction, we find that the gene expression data of proper modules have high eigengene factorizabilities <italic>EF</italic>(<italic>X</italic>)≥0.97 (##TAB##1##Table 2##). By contrast, the factorizability of the grey genes (i.e., the genes outside of proper modules) is relatively low (<italic>EF</italic>(<italic>X</italic>) = 0.66).</p> <p>For each module, ##TAB##1##Table 2## reports network properties including network size, density, centralization, heterogeneity, mean clustering coefficient, module significance, hub gene significance, and eigengene significance. For the proper (nongrey) modules, we find that the numerical values of the intramodular network concepts and their eigengene-based analogs support our theoretical derivations.</p> <p>Our empirical results illustrate Observation 2 regarding the relationship between intramodular network concepts and their eigengene-based analogs. ##FIG##5##Figure 6A–E## depict the relationships among centralization, heterogeneity, clustering coefficient, module significance, hub gene significance and their respective eigengene-based analogs when a soft threshold of <italic>β</italic> = 1 is used for the weighted network construction (Equation 2). The analogous results for <italic>β</italic> = 6 are depicted in ##FIG##5##Figure 6G–K##. ##FIG##5##Figure 6F and 6L## depicts the relationship between hub gene significance (Equation 13) and module eigengene significance (Equation 27) for <italic>β</italic> = 1 and <italic>β</italic> = 6, respectively. For completeness, we also report the results for the grey, nonmodule genes in the figures. But since our theoretical results assume proper modules, we exclude the grey genes from the calculation of the squared correlation coefficient <italic>R</italic>
<sup>2</sup>. The summary of a robustness analysis with regard to different soft thresholds <italic>β</italic> and hard thresholds <italic>τ</italic> is reported in ##TAB##2##Table 3## and ##SUPPL##0##Text S1##. Overall, we find very high squared correlations (<italic>R</italic>
<sup>2</sup>&gt;0.85), which confirm our theoretical results. Only the <italic>R</italic>
<sup>2</sup> values for the relationship between clustering coefficient and its eigengene-based analog is decreased if <italic>β</italic>&gt;3.</p> <p>
##FIG##7##Figure 8## illustrate the implications of Observation 3 regarding the relationships among network concepts in the cancer coexpression module networks. ##FIG##7##Figure 8A## shows that the scaled connectivity <italic>K<sub>i</sub></italic>
<sup>(<italic>q</italic>)</sup> is highly correlated (<italic>R</italic>
<sup>2</sup>&gt;0.99) with <italic>a<sub>e,i</sub></italic>
<sup>(<italic>q</italic>)</sup>, which illustrates Equation 33. This relationship is highly robust with regard to high soft thresholds <italic>β</italic> as can be seen from ##TAB##2##Table 3##.</p> <p>
##FIG##7##Figure 8B## illustrates the relationship between the clustering coefficient (the mean corresponds to the short horizontal line) and (1+<italic>Heterogeneity</italic>
<sup>2</sup>)<sup>2</sup>×<italic>Density</italic> (Equation 31). This relationship is diminished for soft thresholds <italic>β</italic>&gt;3 as can be seen from ##TAB##2##Table 3##.</p> <p>
##FIG##7##Figure 8C## illustrates the relation  (Equation 37), which is highly robust with regard to different choices of <italic>β</italic> (##TAB##2##Table 3##). ##FIG##7##Figure 8D## illustrates  (Equation 40). This relationship is <italic>not</italic> robust with regard to <italic>β</italic>: the <italic>R<sup>2</sup></italic> value is only 0.058 for <italic>β</italic> = 3. ##FIG##7##Figure 8E## illustrates  (Equation 38), which is highly robust with regard to <italic>β</italic> (##TAB##2##Table 3##).</p> <p>Although our theoretical results were derived using relatively restrictive assumptions, we find that most results are robust in the weighted networks, see ##FIG##8##Figure 9##, ##TAB##2##Table 3##, and ##SUPPL##0##Text S1##. However, in unweighted networks, several relationships have lower <italic>R</italic>
<sup>2</sup> values and show a strong dependence on the hard threshold <italic>τ</italic> (##TAB##2##Table 3##).</p> <title>Motivational example: Mouse tissues of an F2 intercross</title> <p>The mouse tissues came from an F2 intercross between two mouse strains C3H/HeJ and C57BL/6J. The data were already described above and in ##FIG##0##Figure 1##. The 498 genes were part of a body weight related module in liver tissue (the Blue module described in reference ##UREF##1##[23]##). ##TAB##3##Table 4## presents network concepts and their eigengene-based analogs in the different tissue networks. As predicted by Observation 2, we find a close relationship between the two types of network concepts if the eigengene factorizability of the corresponding network is close to 1. This example also illustrates that our results apply to coexpression networks comprised of relatively few genes (here 498 genes).</p> <title>Mouse gene coexpression network application</title> <p>Here we focus on the female mouse liver tissues of the above-mentioned F2 mouse cross. Specifically, 135 female mice were used to construct a weighted network comprised of 3,400 highly connected genes. The biological significance and gene ontology enrichment analysis of the 12 modules in this large network is described in ##UREF##1##[23]##. In ##SUPPL##1##Text S2##, ##TAB##4##Table 5##, and ##FIG##8##Figure 9##, we focus on the relationships among the network concepts. We find that many of our theoretical results hold approximately even if the expression factorizability is low. ##TAB##4##Table 5## shows how the relationship (<italic>R</italic>
<sup>2</sup> values) between network concepts and their eigengene-based analogs depend on the soft threshold <italic>β</italic>. Overall, we find that our theoretical results are highly robust in weighted networks. The relationship between the clustering coefficient and its eigengene-based analog is diminished (down to 0.44) for <italic>β</italic>&gt;3. The relationship between heterogeneity and its eigengene-based analog is diminished (down to 0.54 when <italic>β</italic>&lt;3).</p> <p>The relation  (Equation 40) has a relatively low <italic>R</italic>
<sup>2</sup> value (down to 0.21) for low values of <italic>β</italic>≤3 but the other relationships among network concepts are highly robust with respect to <italic>β</italic>. For unweighted networks, the <italic>R</italic>
<sup>2</sup> values tend to be lower and several relationships show a marked dependency on the hard threshold <italic>τ</italic> (##TAB##4##Table 5##).</p> <title>Yeast gene coexpression network application</title> <p>In ##SUPPL##2##Text S3##, we illustrate our theoretical derivations using a yeast gene coexpression network. The yeast microarray data were derived from experiments designed to study the cell cycle ##REF##9843569##[54]##. A detailed biological description of the modules and the importance of intramodular connectivity can be found in previous work ##REF##16515682##[33]##. In ##SUPPL##2##Text S3## and in ##FIG##8##Figure 9##, we use a gene significance measure that encodes knock-out essentiality, i.e., <italic>GS<sub>i</sub></italic> = 1 if the <italic>i</italic>th gene is known to be essential and 0 otherwise. In contrast to the other applications, <italic>this gene significance measure is not based on a sample trait</italic>. Our theoretical derivations for relating module significance to hub gene significance (Equation 37) assumed a sample-trait based gene significance measure. Although this important assumption is violated for knock-out essentiality, it is striking that the relationship between hub gene significance and module significance can still be observed empirically (##FIG##8##Figure 9##).</p> <p>
##TAB##5##Table 6## shows how the relationship (squared correlation <italic>R</italic>
<sup>2</sup>) between network concepts and their eigengene-based analogs depend on the soft threshold <italic>β</italic>. Overall, we find that our theoretical results are highly robust in weighted networks. The relation  (Equation 40) breaks down for <italic>β</italic> = 3 or 4 but the other relationships among network concepts are highly robust with respect to <italic>β</italic>. For unweighted networks, the <italic>R</italic>
<sup>2</sup> values tend to be lower and several relationships show a marked dependency on the hard threshold <italic>τ</italic> (##TAB##5##Table 6##).</p>
<title>Discussion</title> <p>Network theoretic methods and concepts are increasingly used for the systems biologic analysis of microarray data. We illustrate how network concepts can be used for describing large correlation matrices and for arriving at biologically plausible data reduction techniques. Many alternative approaches for defining gene coexpression networks are possible, e.g., ##REF##11027309##[13]##, ##REF##11099257##[55]##–##REF##17784779##[61]##. Here we define the network adjacency and the gene significance measure in terms of correlations since this allows us to interpret pairwise relations in terms of angles between scaled versions of the variables. For example, the sample trait based gene significance measure of the <italic>i</italic>th gene is determined by the angle between the <italic>i</italic>th gene expression profile and the sample trait <italic>T</italic> (Equation 4); the scaled intramodular connectivity of the <italic>i</italic>th gene (Equation 33) is determined by the angle between the <italic>i</italic>th gene expression profile and the module eigengene; the hub gene significance (Equation 34) is determined by the angle between module eigengene and the sample trait.</p> <p>The geometric interpretation of gene coexpression network analysis reveals a deep connection to other statistical methods. Since it projects the gene expressions profiles onto the hypersphere in an <italic>m</italic>-dimensional Euclidean space, network analysis can be considered a special case of directional statistics. When focusing on the use of module eigengenes, network analysis can be considered a variant of oblique factor analysis.</p> <p>A high level view of modules and their centroids (eigengenes) can be used to define eigengene networks ##REF##18031580##[52]##. High correlations (small angles) between module eigengenes may suggest close relationships between the corresponding pathways. A low level view of a single module allows us to provide a geometric interpretation of intramodular network concepts. We use the singular value decomposition of module expression data to characterize approximately factorizable gene coexpression networks, i.e., adjacency matrices that satisfy <italic>a<sub>ij</sub></italic>
<sup>(<italic>q</italic>)</sup>≈<italic>CF<sub>i</sub></italic>
<sup>(<italic>q</italic>)</sup>
<italic>CF<sub>j</sub></italic>
<sup>(<italic>q</italic>)</sup>. We provide an intuitive formula of the conformity <italic>CF<sub>i</sub></italic>
<sup>(<italic>q</italic>)</sup>≈|cor(<italic>x<sub>i</sub></italic>
<sup>(<italic>q</italic>)</sup>,<italic>E</italic>
<sup>(<italic>q</italic>)</sup>)|<italic><sup>β</sup></italic>. Since the module eigengene <italic>E</italic>
<sup>(<italic>q</italic>)</sup> summarizes the overall behavior of the module, the eigengene conformity |cor(<italic>x<sub>i</sub></italic>
<sup>(<italic>q</italic>)</sup>,<italic>E</italic>
<sup>(<italic>q</italic>)</sup>)|<italic><sup>β</sup></italic> measures how well gene <italic>i</italic> conforms to the overall module. This insight led us to coin the term “conformity”. Using the singular values, we propose a measure of eigengene factorizability (Equation 24) that is analogous to the proportion of variance explained by the module eigengene (Equation 22). We provide a geometric interpretation of network factorizability in ##FIG##4##Figure 5A##.</p> <p>The derivation of Observation 1 in the <xref ref-type="sec" rid="s4">Methods</xref> section highlights a theoretical advantage of the soft-thresholding approach (Equation 2); the resulting weighted network maintains the approximate factorizability of the underlying correlation matrix: <italic>a<sub>ij</sub></italic>
<sup>(<italic>q</italic>)</sup> = |cor(<italic>x<sub>i</sub></italic>
<sup>(<italic>q</italic>)</sup>,<italic>x<sub>j</sub></italic>
<sup>(q)</sup>)|<italic><sup>β</sup></italic>≈|cor(<italic>x<sub>i</sub></italic>
<sup>(<italic>q</italic>)</sup>,<italic>E</italic>
<sup>(<italic>q</italic>)</sup>)cor(<italic>x<sub>j</sub></italic>
<sup>(<italic>q</italic>)</sup>,<italic>E</italic>
<sup>(<italic>q</italic>)</sup>)|<italic><sup>β</sup></italic> = |cor(<italic>x<sub>i</sub></italic>
<sup>(<italic>q</italic>)</sup>,<italic>E</italic>
<sup>(<italic>q</italic>)</sup>)|<italic><sup>β</sup></italic>|cor(<italic>x<sub>j</sub></italic>
<sup>(<italic>q</italic>)</sup>,<italic>E</italic>
<sup>(<italic>q</italic>)</sup>)|<italic><sup>β</sup></italic>.</p> <p>Using multiple different gene coexpression networks from mouse tissues, brain cancer, and yeast, we provide empirical evidence that coexpression modules tend to have high eigengene factorizability and that the maximum conformity assumption (Equation 32) is satisfied for low powers of <italic>β</italic>.</p> <p>We propose eigengene-based analogs of network concepts (Equation 30). While network concepts are functions of the adjacency matrix, eigengene-based network concepts are analogous functions of the eigengene conformities |cor(<italic>x<sub>i</sub></italic>
<sup>(<italic>q</italic>)</sup>,<italic>E</italic>
<sup>(<italic>q</italic>)</sup>)|<italic><sup>β</sup></italic>. Algebraically, eigengene-based network concepts are closely related to “approximate conformity based” network concepts ##REF##17547772##[8]## but they allow for a geometric interpretation.</p> <p>We use the correspondence between intramodular network concepts and their eigengene-based analogs to provide a geometric interpretation of network concepts. Observation 2 states that network concepts in weighted gene coexpression module networks are approximately equal to their eigengene-based analogs. A major theoretical advantage of eigengene-based network concepts is that they reveal simple relationships. To arrive at particularly simple relationships, we make the maximum conformity assumption (Equation 32) for the results presented in the main text. ##TAB##0##Table 1## provides a rough dictionary for translating between gene coexpression network analysis and the singular value decomposition if the underlying expression data have high eigengene factorizability (say <italic>EF</italic>(<italic>X</italic>
<sup>(<italic>q</italic>)</sup>)&gt;0.95) and if the maximum conformity assumption (Equation 32) is satisfied. However, even if the maximum conformity assumption does not hold, one can still find simple relationships among the network concepts (Equation 49).</p> <p>The geometric interpretation of gene coexpression networks facilitates the derivation of several results that should be interesting to network theorists. For example, we argue that highly connected intramodular hub genes cannot be intermediate between two distinct coexpression modules (##FIG##4##Figure 5B##). The geometric interpretation is particularly useful when studying gene significance and module significance measures that are based on a microarray sample trait (Equation 4). To study the relationship between connectivity and gene significance, we propose a novel measure of hub gene significance (Equation 13). We find that the hub gene significance of a module network is determined by the angle between the module eigengene and the microarray sample trait (Equation 34). Our geometric interpretation of coexpression networks allows us to describe situations when a module has low hub gene significance (##FIG##4##Figure 5C and 5D##). Our theoretical derivations for relating module significance to hub gene significance (Equation 37) assumes a gene significance measure based on a sample trait. Although this important assumption is violated for the gene significance measure (knock-out essentiality) in the yeast network, it is striking that the relationship between hub gene significance and module significance can still be observed in this application (##FIG##8##Figure 9##).</p> <p>We provide a robustness analysis that shows that many of our theoretical results apply even if our underlying assumptions are not satisfied (##FIG##5##Figures 6## and ##FIG##8##9##, ##TAB##2##Tables 3##, ##TAB##4##5##, and ##TAB##5##6##, ##SUPPL##0##Text S1##, ##SUPPL##1##Text S2##, and ##SUPPL##2##Text S3##). We find that the correspondence between network concepts and their eigengene-based analogs is often better in weighted networks than in unweighted networks. Further, we find that the results in weighted networks tend to be more robust than those in unweighted networks with regard to changing the network construction thresholds <italic>β</italic> and <italic>τ</italic>, respectively. Thus, weighted coexpression networks are preferable over unweighted networks when a geometric interpretation of network concepts is desirable.</p> <p>The correspondence between coexpression module networks and the singular value decomposition (##TAB##0##Table 1##) can break down when a high soft threshold is used for constructing a weighted network or when dealing with an unweighted network. Thus, eigengene-based concepts do not replace network concepts when describing interaction patterns among genes.</p> <p>While this article has a theoretical bent, we illustrate the results on three different microarray data sets (human, mouse, and yeast) that are described in our online R software tutorials, in ##SUPPL##0##Text S1##, ##SUPPL##1##Text S2##, and ##SUPPL##2##Text S3##. Our theoretical results also apply to networks comprised of genes that are highly correlated with a sample trait. The key assumption underlying our results is high eigengene factorizability <italic>EF</italic>(<italic>X</italic>
<sup>(<italic>q</italic>)</sup>). To illustrate this point, ##SUPPL##3##Text S4## describes a brain cancer network comprised of the 500 genes with highest absolute correlation with brain cancer survival time. Our results illustrate that the geometric interpretation of gene coexpression networks has important theoretical and practical implications that may guide the development and application of network methods.</p>
<p>Conceived and designed the experiments: SH. Analyzed the data: JD. Wrote the paper: SH JD.</p> <p>The merging of network theory and microarray data analysis techniques has spawned a new field: gene coexpression network analysis. While network methods are increasingly used in biology, the network vocabulary of computational biologists tends to be far more limited than that of, say, social network theorists. Here we review and propose several potentially useful network concepts. We take advantage of the relationship between network theory and the field of microarray data analysis to clarify the meaning of and the relationship among network concepts in gene coexpression networks. Network theory offers a wealth of intuitive concepts for describing the pairwise relationships among genes, which are depicted in cluster trees and heat maps. Conversely, microarray data analysis techniques (singular value decomposition, tests of differential expression) can also be used to address difficult problems in network theory. We describe conditions when a close relationship exists between network analysis and microarray data analysis techniques, and provide a rough dictionary for translating between the two fields. Using the angular interpretation of correlations, we provide a geometric interpretation of network theoretic concepts and derive unexpected relationships among them. We use the singular value decomposition of module expression data to characterize approximately factorizable gene coexpression networks, i.e., adjacency matrices that factor into node specific contributions. High and low level views of coexpression networks allow us to study the relationships among modules and among module genes, respectively. We characterize coexpression networks where hub genes are significant with respect to a microarray sample trait and show that the network concept of intramodular connectivity can be interpreted as a fuzzy measure of module membership. We illustrate our results using human, mouse, and yeast microarray gene expression data. The unification of coexpression network methods with traditional data mining methods can inform the application and development of systems biologic methods.</p> <title>Author Summary</title> <p>Similar to natural languages, network language is ever evolving. While some network terms (concepts) are widely used in gene coexpression network analysis, others still need to be developed to meet the ever increasing demand for describing the system of gene transcripts. There is a need to provide an intuitive geometric explanation of network concepts and to study their relationships. For example, we show that certain seemingly disparate network concepts turn out to be synonyms in the context of coexpression modules. We show how coexpression network language affects our understanding of biology. For example, there are geometric reasons why highly connected hub genes in important coexpression modules tend to be important, and why hub genes in one module cannot be hubs in another distinct module. We provide a short dictionary for translating between microarray data analysis language and network theory language to facilitate communication between the two fields. We describe several examples that illustrate how the two data analysis fields can inform each other.</p>
<title>Supporting Information</title>
<p>We are grateful for discussions with Andy Yip, Lora Bagryanova, Dan Geschwind, Peter Langfelder, Tova Fuller, Jake Lusis, Tom Drake, Paul Mischel, Stan Nelson, Mike Oldham, Anja Presson, Lin Wang, and Nan Zhang.</p>
<fig id="pcbi-1000117-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000117.g001</object-id><label>Figure 1</label><caption><title>This motivational example explores the pairwise absolute correlations <italic>a<sub>ij</sub></italic> = |cor(<italic>x<sub>i</sub></italic>,<italic>x<sub>j</sub></italic>)| among 498 genes in different mouse tissues.</title><p>The biological significance of this network is described in ##UREF##1##[23]##,##REF##17668265##[26]##. Each figure panel contains 8 subfigures for different genders and tissue types (liver, adipose, brain, muscle). (A) An average linkage hierarchical cluster tree of the genes. (B) The corresponding heat maps, which color-code the absolute pairwise correlations <italic>a<sub>ij</sub></italic>: red and green in the heat map indicate high and low absolute correlation, respectively. The genes in the rows and columns of each heat map are sorted by the corresponding cluster tree. (C) The relationship between gene significance <italic>GS</italic> (<italic>y</italic>-axis) and connectivity (<italic>x</italic>-axis). The gene significance of the <italic>i</italic>th gene was defined as the absolute correlation between the <italic>i</italic>th gene expression profile and mouse body weight. The hub gene significance <italic>HGS</italic> (Equation 13) is defined as the slope of the red line, which results from a regression model without an intercept term.</p></caption></fig> <fig id="pcbi-1000117-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000117.g002</object-id><label>Figure 2</label><caption><title>Relationships among maximum adjacency ratio, scaled connectivity, and gene significance.</title><p>(A) The relationship between <italic>MAR<sub>i</sub></italic> (y-axis) and scaled connectivity <italic>K<sub>i</sub></italic> using the female mouse muscle tissue network described in the motivational example. The genes are colored red or black depending on whether they are significantly (<italic>p</italic>-value&lt;0.05) related to mouse body weight. (B) Boxplots and a Kruskal-Wallis test <italic>p</italic>-value (<italic>p</italic> = 0.00072) for studying whether <italic>MAR<sub>i</sub></italic> differs between significant (red) and non-significant (black) genes. (C) The analogous boxplots and <italic>p</italic>-value for the scaled connectivity <italic>K<sub>i</sub></italic>. In this female muscle tissue application, <italic>MAR<sub>i</sub></italic> is more significantly (<italic>p</italic> = 0.00072) related to <italic>GS<sub>i</sub></italic> than is <italic>K<sub>i</sub></italic> (<italic>p</italic> = 0.051). (D,E,F) The analogous relationships for male muscle. Here, the <italic>MAR<sub>i</sub></italic> is more significantly (<italic>p</italic> = 0.00014) related to <italic>GS<sub>i</sub></italic> than is <italic>K<sub>i</sub></italic> (<italic>p</italic> = 0.0034). (G,H,I) The analogous relationships for the brown module of the brain cancer application. Here, the <italic>MAR<sub>i</sub></italic> is slightly more significantly (<italic>p</italic> = 1.6E-8) related to <italic>GS<sub>i</sub></italic> than is <italic>K<sub>i</sub></italic> (<italic>p</italic> = 2.6E-7). As a caveat, we mention that in other applications (e.g., the yeast network), we have found that <italic>K<sub>i</sub></italic> is more significantly related to <italic>GS<sub>i</sub></italic> than <italic>MAR<sub>i</sub></italic>.</p></caption></fig> <fig id="pcbi-1000117-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000117.g003</object-id><label>Figure 3</label><caption><title>Overview and an example application of gene coexpression network analysis.</title><p>(A) Outline of an analysis flow chart. Gene coexpression network analysis aims to identify pathways (modules) and their key drivers (e.g., intramodular hub genes). (B) The hierarchical cluster tree of genes in the brain cancer network. Modules correspond to branches of the tree. The branches and module genes are assigned a color as can be seen from the color-bands underneath the tree. Grey denotes genes outside of proper modules. A functional enrichment analysis of these modules can be found in Horvath et al. (2006). (C) The module significance (average gene significance) of the modules. The underlying gene significance is defined with respect to the patient survival time (Equation 4). (D,E) Scatter plots of gene significance <italic>GS</italic> (<italic>y</italic>-axis) versus scaled connectivity <italic>K</italic> (<italic>x</italic>-axis) in the brown and blue module, respectively. The hub gene significance (Equation 13) is defined as the slope of the red line, which results from a regression model without an intercept term.</p></caption></fig> <fig id="pcbi-1000117-g004" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000117.g004</object-id><label>Figure 4</label><caption><title>Module eigengenes in the brain cancer gene coexpression network.</title><p>(A) The pairwise scatter plots among the module eigengenes <italic>E</italic>
<sup>(<italic>q</italic>)</sup> of different modules and cancer survival time <italic>T</italic>. Each dot represents a microarray sample. ME.blue denotes the module eigengene <italic>E</italic>
<sup>(blue)</sup> of the blue module. Numbers below the diagonal are the absolute values of the corresponding correlations. Note that the module eigengenes of different modules can be highly correlated. The brown module eigengene has the highest absolute correlation (<italic>r</italic> = 0.20) with survival time. Frequency plots (histograms) of the variables are plotted along the diagonal. (B) Upper panel: heat map plot of the brown module gene expression profiles (rows) across the microarray samples (columns). Red corresponds to high- and green to low- expression values. Since the genes of a module are highly correlated, one observes vertical bands. (B) Lower panel: the values of the components of the module eigengene (<italic>y</italic>-axis) versus microarray sample number (<italic>x</italic>-axis). Note that vertical bands of red (green) in the upper panel correspond to high (low) values of the eigengene in the lower panel. (C) The expression profile of the module eigengene (<italic>y</italic>-axis) is highly correlated with that of the most highly connected hub gene (<italic>x</italic>-axis). A linear regression line has been added.</p></caption></fig> <fig id="pcbi-1000117-g005" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000117.g005</object-id><label>Figure 5</label><caption><title>Using vectors to illustrate results in gene coexpression network analysis.</title><p>(A) A geometric interpretation of factorizability if the gene expression profiles and the module eigengene lie in a Euclidean plane. Then the angle <italic>θ</italic>
<sub>12</sub> between gene expressions profiles 1 and 2 can be expressed in terms of angles with the module eigengene, i.e., <italic>θ</italic>
<sub>12</sub> = <italic>θ</italic>
<sub>1</sub>−<italic>θ</italic>
<sub>2</sub>. Similarly, <italic>θ</italic>
<sub>23</sub> = <italic>θ</italic>
<sub>2</sub>+<italic>θ</italic>
<sub>3</sub>. Under the assumptions stated in the text, we find <italic>θ<sub>ij</sub></italic>≈|<italic>θ<sub>i</sub></italic>±<italic>θ<sub>j</sub></italic>|. Using a trigonometric formula (Equation 51) this implies that the correlation matrix is approximately factorizable. (B) Illustrating why intramodular hub genes cannot be “intermediate” genes between two distinct coexpression modules. The large angle between module eigengenes <italic>E1</italic> and <italic>E2</italic> reflects that the corresponding modules are distinct. Since intermediate gene 1 does not have a small angle with either eigengene, it is not an intramodular hub gene. By contrast, intramodular hub gene 2 has a small angle with eigengene <italic>E1</italic> but is not close to module eigengene <italic>E2</italic>. (C,D) Illustrating that the hub gene significance of a module depends on the relationship between the module eigengene and the underlying microarray sample trait (Equation 34). For sample traits <italic>T2</italic> and <italic>T1</italic> the hub gene significance (and corresponding eigengene significance cor(<italic>E</italic>,<italic>T</italic>)) are high and low, respectively. The geometry of (C) implies relationships between the connectivity <italic>k</italic> of a gene (determined by its angle with the eigengene E) and gene significance measure <italic>GS1</italic> (its angle with trait <italic>T1</italic>) and <italic>GS2</italic> (its angle with trait <italic>T2</italic>). As shown in (D), the gene significance measure <italic>GS2</italic> increases with <italic>k</italic> since the small angle between <italic>E</italic> and <italic>T2</italic> implies that genes with high <italic>k</italic> (small angle with <italic>E</italic>) also have a small angle with <italic>T2</italic>. In contrast, high connectivity <italic>k</italic> implies a large angle with <italic>T1</italic> and thus <italic>GS1</italic> decreases as a function of <italic>k</italic>.</p></caption></fig> <fig id="pcbi-1000117-g006" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000117.g006</object-id><label>Figure 6</label><caption><title>Illustrating Observation 2 regarding the relationship between a network concept (<italic>y</italic>-axis) and its eigengene-based analog (<italic>x</italic>-axis) in the brain cancer data.</title><p>Each point corresponds to a module. (A–F) Corresponding to a weighted network constructed with a soft threshold (Equation 2) of <italic>β</italic> = 1. (G–L) Analogous plots for <italic>β</italic> = 6. (A,G) Centralization (<italic>y</italic>-axis) versus eigengene-based Centralization<sub>E</sub> (<italic>x</italic>-axis). The following are analogous plots for (B,H): heterogeneity; (C,I) clustering coefficient; (D,J) module significance; and (E,K) hub gene significance. (F,L) Illustrating Equation 13 regarding the relationship between eigengene significance and hub gene significance. The blue line is the regression line through the points representing proper modules (i.e., the grey, nonmodule genes are left out). While the red reference line (slope 1, intercept 0) does not always fit well, we observe high squared correlations <italic>R</italic>
<sup>2</sup> between network concepts and their analogs. Since the grey point corresponds to the genes outside properly defined modules, we did not include it in calculations.</p></caption></fig> <fig id="pcbi-1000117-g007" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000117.g007</object-id><label>Figure 7</label><caption><title>Fuzzy module annotation of genes in the brain cancer network.</title><p>A natural choice for a fuzzy measure of module membership is the generalized scaled connectivity measure <italic>K</italic>
<sub>cor,<italic>i</italic></sub>
<sup>(<italic>q</italic>)</sup> = |cor(<italic>x<sub>i</sub></italic>,<italic>E</italic>
<sup>(<italic>q</italic>)</sup>)| (Equation 41). (A) Scatterplot of the brown module membership measure (<italic>y</italic>-axis) versus that of the blue module (<italic>x</italic>-axis). Note that grey dots corresponding to genes outside of properly defined modules can be intermediate between module genes. (B) The corresponding plot for blue versus turquoise module membership. (C) Brown versus turquoise module membership. (D) The relationship between gene significance based on survival time (<italic>y</italic>-axis) and brown module membership (<italic>x</italic>-axis).</p></caption></fig> <fig id="pcbi-1000117-g008" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000117.g008</object-id><label>Figure 8</label><caption><title>Using the brain cancer data to illustrate Observation 3 regarding the relationships among network concepts.</title><p>(A) Illustrating Equation 33 regarding the relationship between scaled intramodular connectivity <italic>K<sub>i</sub></italic>
<sup>(<italic>q</italic>)</sup> (<italic>y</italic>-axis) and eigengene conformity <italic>a<sub>e</sub></italic>
<sub>,<italic>i</italic></sub> (<italic>x</italic>-axis). Each dot corresponds to a gene colored by its module membership. We find a high squared correlation <italic>R</italic>
<sup>2</sup> even for the grey genes outside properly defined modules. (B) Illustrating Equation 31 regarding the relationship between the clustering coefficient and (1+<italic>Heterogeneity</italic>
<sup>2</sup>)<sup>2</sup>×<italic>Density</italic>. Again each dot represents a gene. The clustering coefficients of grey genes vary more than those of genes in proper modules. The short horizontal lines correspond to the mean clustering coefficient of each module. (C) Illustrating  (Equation 37); here each dot corresponds to a module. Since the grey dot corresponds to genes outside of properly defined modules, we have excluded it from the calculation of the squared correlation <italic>R</italic>
<sup>2</sup>. (D) Illustrating  (Equation 40); (E) Illustrating  (Equation 38). A reference line (red) with intercept 0 and slope 1 has been added to each plot. The blue line is the regression line through the points representing proper modules (i.e., the grey, non-module genes are left out). A robustness analysis with regard to different network construction methods, e.g., <italic>β</italic>&gt;1, can be found in ##SUPPL##0##Text S1##.</p></caption></fig> <fig id="pcbi-1000117-g009" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000117.g009</object-id><label>Figure 9</label><caption><title>Using three different data (brain cancer, mouse liver, and yeast cell cycle) and three different network construction methods to illustrate Equation 37 regarding the relationship between module significance (<italic>y</italic>-axis) and  (<italic>x</italic>-axis).</title><p>Points correspond to modules. The square of the correlation coefficient <italic>R</italic>
<sup>2</sup> was computed without the grey, improper module. (A,D,G) Corresponding to the brain cancer gene coexpression networks. (B,E,H) Corresponding to mouse liver networks. (C,F,I) Corresponding to yeast networks. (A–C) Corresponding to a weighted network (Equation 2) constructed with soft thresholds <italic>β</italic> = 1. (D–F) Corresponding to <italic>β</italic> = 6. (G–I) Corresponding to an unweighted network (Equation 1) that results from thresholding the correlation matrix at <italic>τ</italic> = 0.5. Overall, we find that the reported relationship is quite robust with respect to our theoretical assumptions (e.g., factorizability). The blue line is the regression line through the points representing proper modules (i.e., the grey, nonmodule genes are left out). A reference line with slope 1 and intercept 0 is shown in red. Additional details can be found in ##SUPPL##0##Text S1##, ##SUPPL##1##Text S2##, and ##SUPPL##2##Text S3##.</p></caption></fig>
<table-wrap id="pcbi-1000117-t001" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000117.t001</object-id><label>Table 1</label><caption><title>Dictionary for translating between general network terms and their eigengene-based counterparts.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Term</td><td align="left" rowspan="1" colspan="1">General network</td><td align="left" rowspan="1" colspan="1">Gene expression</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Adjacency matrix <italic>A</italic>
<sup>(<italic>q</italic>)</sup> = [<italic>a<sub>ij</sub></italic>]</td><td align="left" rowspan="1" colspan="1">Microarray data <italic>X</italic>
<sup>(<italic>q</italic>)</sup>
</td></tr><tr><td align="left" rowspan="1" colspan="1">Decomposition</td><td align="left" rowspan="1" colspan="1">Factor analysis of <italic>A</italic>
</td><td align="left" rowspan="1" colspan="1">Singular value decomposition of <italic>X</italic> = <italic>UDV<sup>T</sup></italic>
</td></tr><tr><td align="left" rowspan="1" colspan="1">Centroid</td><td align="left" rowspan="1" colspan="1">Intramodular hub gene</td><td align="left" rowspan="1" colspan="1">Module eigengene <italic>E</italic>
</td></tr><tr><td align="left" rowspan="1" colspan="1">Conformity(i)</td><td align="left" rowspan="1" colspan="1">
<italic>CF<sub>i</sub></italic> defined as 1st factor of <italic>A</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>a<sub>e</sub></italic>
<sub>,<italic>i</italic></sub> = |cor(<italic>x<sub>i</sub></italic>,<italic>E</italic>)|<italic><sup>β</sup></italic>
</td></tr><tr><td align="left" rowspan="1" colspan="1">Approximately factorizable means</td><td align="left" rowspan="1" colspan="1">
<italic>a<sub>ij</sub></italic>≈<italic>CF<sub>i</sub>CF<sub>j</sub></italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>x<sub>i</sub></italic>≈<italic>u</italic>
<sub>1</sub>(<italic>i</italic>)|<italic>d</italic>
<sub>1</sub>|<italic>E</italic>
</td></tr><tr><td align="left" rowspan="1" colspan="1">Factorizability measure</td><td align="left" rowspan="1" colspan="1">

</td><td align="left" rowspan="1" colspan="1">

</td></tr><tr><td align="left" rowspan="1" colspan="1">CentroidSignif(i)</td><td align="left" rowspan="1" colspan="1">
<italic>GS<sub>i,centroid</sub></italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>a<sub>e</sub></italic>
<sub>,<italic>t</italic></sub> = |cor(<italic>E</italic>,<italic>T</italic>)|<italic><sup>β</sup></italic>
</td></tr><tr><td align="left" rowspan="1" colspan="1">CentroidConformity(i)</td><td align="left" rowspan="1" colspan="1">
<italic>a<sub>i.centroid,i</sub></italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>a<sub>e</sub></italic>
<sub>,<italic>i</italic></sub> = |cor(<italic>E</italic>,<italic>x<sub>i</sub></italic>)|<italic><sup>β</sup></italic>
</td></tr><tr><td colspan="3" align="left" rowspan="1">Weighted gene coexpression network and its eigengene-based approximation if <italic>EF</italic>(<italic>X</italic>
<sup>(<italic>q</italic>)</sup>)≈1</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Coexpression network</td><td align="left" rowspan="1" colspan="1">Eigengene-based counterpart</td></tr><tr><td align="left" rowspan="1" colspan="1">Network</td><td align="left" rowspan="1" colspan="1">
<italic>A</italic> = |cor(<italic>X</italic>)|<italic><sup>β</sup></italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>A<sub>E</sub></italic> = <italic>a<sub>e</sub>a<sub>e</sub><sup>T</sup></italic>
</td></tr><tr><td align="left" rowspan="1" colspan="1">Gene significance(i)</td><td align="left" rowspan="1" colspan="1">
<italic>GS<sub>i</sub></italic> = |cor(<italic>x<sub>i</sub></italic>,<italic>T</italic>)|<italic><sup>β</sup></italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>GS<sub>E</sub></italic>
<sub>,<italic>i</italic></sub> = <italic>a<sub>e,i</sub>a<sub>e,t</sub></italic>
</td></tr><tr><td align="left" rowspan="1" colspan="1">Connectivity(<italic>i</italic>)</td><td align="left" rowspan="1" colspan="1">
<italic>k<sub>i</sub></italic> = <italic>Σ<sub>j≠i</sub>a<sub>ij</sub></italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>k<sub>E,i</sub></italic> = <italic>a<sub>e</sub></italic>
<sub>,<italic>i</italic></sub>
<italic>Σ<sub>j</sub>a<sub>e</sub></italic>
<sub>,<italic>j</italic></sub>
</td></tr><tr><td colspan="3" align="left" rowspan="1">Network concepts based on a network concept function <italic>NCF</italic>(·,·) if <italic>EF</italic>(<italic>X</italic>
<sup>(<italic>q</italic>)</sup>)≈1 and max<italic><sub>j</sub></italic>(<italic>a<sub>e</sub></italic>
<sub>,<italic>j</italic></sub>)≈1</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Intramodular</td><td align="left" rowspan="1" colspan="1">Eigengene-based</td></tr><tr><td align="left" rowspan="1" colspan="1">Concept</td><td align="left" rowspan="1" colspan="1">
<italic>NCF(A,GS)</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>NCF(A<sub>E</sub>,GS<sub>E</sub>)</italic>
</td></tr><tr><td align="left" rowspan="1" colspan="1">Scaled connectivity(<italic>i</italic>)</td><td align="left" rowspan="1" colspan="1">

</td><td align="left" rowspan="1" colspan="1">
<italic>K<sub>E,i</sub></italic>≈<italic>a<sub>e</sub></italic>
<sub>,<italic>i</italic></sub>
</td></tr><tr><td align="left" rowspan="1" colspan="1">Density</td><td align="left" rowspan="1" colspan="1">

</td><td align="left" rowspan="1" colspan="1">

</td></tr><tr><td align="left" rowspan="1" colspan="1">Centralization</td><td align="left" rowspan="1" colspan="1">

</td><td align="left" rowspan="1" colspan="1">

</td></tr><tr><td align="left" rowspan="1" colspan="1">Heterogeneity</td><td align="left" rowspan="1" colspan="1">

</td><td align="left" rowspan="1" colspan="1">

</td></tr><tr><td align="left" rowspan="1" colspan="1">Clustering Coefficient(<italic>i</italic>)</td><td align="left" rowspan="1" colspan="1">

</td><td align="left" rowspan="1" colspan="1">

</td></tr><tr><td align="left" rowspan="1" colspan="1">Max. adjacency ratio(<italic>i</italic>)</td><td align="left" rowspan="1" colspan="1">

</td><td align="left" rowspan="1" colspan="1">

</td></tr><tr><td align="left" rowspan="1" colspan="1">Hub gene significance</td><td align="left" rowspan="1" colspan="1">

</td><td align="left" rowspan="1" colspan="1">
<italic>a<sub>e</sub></italic>
<sub>,<italic>t</italic></sub>
</td></tr><tr><td align="left" rowspan="1" colspan="1">Module significance</td><td align="left" rowspan="1" colspan="1">

</td><td align="left" rowspan="1" colspan="1">

</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pcbi-1000117-t002" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000117.t002</object-id><label>Table 2</label><caption><title>Values of network concepts in weighted gene coexpression module networks (brain cancer data).</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Module</td><td align="left" rowspan="1" colspan="1">Blue</td><td align="left" rowspan="1" colspan="1">Brown</td><td align="left" rowspan="1" colspan="1">Green</td><td align="left" rowspan="1" colspan="1">Grey</td><td align="left" rowspan="1" colspan="1">Red</td><td align="left" rowspan="1" colspan="1">Turquoise</td><td align="left" rowspan="1" colspan="1">Yellow</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">Module size (<italic>n</italic>
<sup>(<italic>q</italic>)</sup>)</td><td align="left" rowspan="1" colspan="1">606</td><td align="left" rowspan="1" colspan="1">185</td><td align="left" rowspan="1" colspan="1">136</td><td align="left" rowspan="1" colspan="1">1313</td><td align="left" rowspan="1" colspan="1">105</td><td align="left" rowspan="1" colspan="1">1112</td><td align="left" rowspan="1" colspan="1">143</td></tr><tr><td align="left" rowspan="1" colspan="1">Eigengene factorizability (<italic>EF</italic>(<italic>X</italic>
<sup>(<italic>q</italic>)</sup>))</td><td align="left" rowspan="1" colspan="1">0.97</td><td align="left" rowspan="1" colspan="1">0.99</td><td align="left" rowspan="1" colspan="1">0.99</td><td align="left" rowspan="1" colspan="1">0.66</td><td align="left" rowspan="1" colspan="1">0.98</td><td align="left" rowspan="1" colspan="1">0.98</td><td align="left" rowspan="1" colspan="1">0.99</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>VarExplained(E<sup>(q)</sup>)</italic>
</td><td align="left" rowspan="1" colspan="1">0.59</td><td align="left" rowspan="1" colspan="1">0.66</td><td align="left" rowspan="1" colspan="1">0.70</td><td align="left" rowspan="1" colspan="1">0.28</td><td align="left" rowspan="1" colspan="1">0.68</td><td align="left" rowspan="1" colspan="1">0.57</td><td align="left" rowspan="1" colspan="1">0.71</td></tr><tr><td align="left" rowspan="1" colspan="1">Max. conformity max(<italic>a<sub>e</sub></italic>
<sub>,<italic>i</italic></sub>)</td><td align="left" rowspan="1" colspan="1">0.97</td><td align="left" rowspan="1" colspan="1">0.97</td><td align="left" rowspan="1" colspan="1">0.98</td><td align="left" rowspan="1" colspan="1">0.91</td><td align="left" rowspan="1" colspan="1">0.95</td><td align="left" rowspan="1" colspan="1">0.98</td><td align="left" rowspan="1" colspan="1">0.98</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>Density</italic>
</td><td align="left" rowspan="1" colspan="1">0.58</td><td align="left" rowspan="1" colspan="1">0.65</td><td align="left" rowspan="1" colspan="1">0.69</td><td align="left" rowspan="1" colspan="1">0.29</td><td align="left" rowspan="1" colspan="1">0.67</td><td align="left" rowspan="1" colspan="1">0.55</td><td align="left" rowspan="1" colspan="1">0.70</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>Density<sub>E</sub></italic>
</td><td align="left" rowspan="1" colspan="1">0.58</td><td align="left" rowspan="1" colspan="1">0.65</td><td align="left" rowspan="1" colspan="1">0.70</td><td align="left" rowspan="1" colspan="1">0.23</td><td align="left" rowspan="1" colspan="1">0.68</td><td align="left" rowspan="1" colspan="1">0.55</td><td align="left" rowspan="1" colspan="1">0.71</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>Centralization</italic>
</td><td align="left" rowspan="1" colspan="1">0.16</td><td align="left" rowspan="1" colspan="1">0.13</td><td align="left" rowspan="1" colspan="1">0.12</td><td align="left" rowspan="1" colspan="1">0.15</td><td align="left" rowspan="1" colspan="1">0.11</td><td align="left" rowspan="1" colspan="1">0.17</td><td align="left" rowspan="1" colspan="1">0.12</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>Centralization<sub>E</sub></italic>
</td><td align="left" rowspan="1" colspan="1">0.16</td><td align="left" rowspan="1" colspan="1">0.13</td><td align="left" rowspan="1" colspan="1">0.12</td><td align="left" rowspan="1" colspan="1">0.21</td><td align="left" rowspan="1" colspan="1">0.11</td><td align="left" rowspan="1" colspan="1">0.18</td><td align="left" rowspan="1" colspan="1">0.12</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>Heterogeneity</italic>
</td><td align="left" rowspan="1" colspan="1">0.14</td><td align="left" rowspan="1" colspan="1">0.10</td><td align="left" rowspan="1" colspan="1">0.11</td><td align="left" rowspan="1" colspan="1">0.17</td><td align="left" rowspan="1" colspan="1">0.091</td><td align="left" rowspan="1" colspan="1">0.17</td><td align="left" rowspan="1" colspan="1">0.11</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>Heterogeneity<sub>E</sub></italic>
</td><td align="left" rowspan="1" colspan="1">0.14</td><td align="left" rowspan="1" colspan="1">0.10</td><td align="left" rowspan="1" colspan="1">0.11</td><td align="left" rowspan="1" colspan="1">0.44</td><td align="left" rowspan="1" colspan="1">0.091</td><td align="left" rowspan="1" colspan="1">0.17</td><td align="left" rowspan="1" colspan="1">0.11</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>Mean(ClusterCoef)</italic>
</td><td align="left" rowspan="1" colspan="1">0.60</td><td align="left" rowspan="1" colspan="1">0.66</td><td align="left" rowspan="1" colspan="1">0.71</td><td align="left" rowspan="1" colspan="1">0.32</td><td align="left" rowspan="1" colspan="1">0.68</td><td align="left" rowspan="1" colspan="1">0.59</td><td align="left" rowspan="1" colspan="1">0.72</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>ClusterCoef<sub>E</sub></italic>
</td><td align="left" rowspan="1" colspan="1">0.60</td><td align="left" rowspan="1" colspan="1">0.66</td><td align="left" rowspan="1" colspan="1">0.71</td><td align="left" rowspan="1" colspan="1">0.33</td><td align="left" rowspan="1" colspan="1">0.68</td><td align="left" rowspan="1" colspan="1">0.59</td><td align="left" rowspan="1" colspan="1">0.72</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>ModuleSignif</italic>
</td><td align="left" rowspan="1" colspan="1">0.088</td><td align="left" rowspan="1" colspan="1">0.12</td><td align="left" rowspan="1" colspan="1">0.21</td><td align="left" rowspan="1" colspan="1">0.11</td><td align="left" rowspan="1" colspan="1">0.16</td><td align="left" rowspan="1" colspan="1">0.14</td><td align="left" rowspan="1" colspan="1">0.065</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>ModuleSignif<sub>E</sub></italic>
</td><td align="left" rowspan="1" colspan="1">0.018</td><td align="left" rowspan="1" colspan="1">0.093</td><td align="left" rowspan="1" colspan="1">0.21</td><td align="left" rowspan="1" colspan="1">0.008</td><td align="left" rowspan="1" colspan="1">0.16</td><td align="left" rowspan="1" colspan="1">0.13</td><td align="left" rowspan="1" colspan="1">0.039</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>HubGeneSignif</italic>
</td><td align="left" rowspan="1" colspan="1">0.11</td><td align="left" rowspan="1" colspan="1">0.14</td><td align="left" rowspan="1" colspan="1">0.25</td><td align="left" rowspan="1" colspan="1">0.15</td><td align="left" rowspan="1" colspan="1">0.18</td><td align="left" rowspan="1" colspan="1">0.19</td><td align="left" rowspan="1" colspan="1">0.074</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>HubGeneSignif<sub>E</sub></italic>
</td><td align="left" rowspan="1" colspan="1">0.023</td><td align="left" rowspan="1" colspan="1">0.11</td><td align="left" rowspan="1" colspan="1">0.25</td><td align="left" rowspan="1" colspan="1">0.014</td><td align="left" rowspan="1" colspan="1">0.18</td><td align="left" rowspan="1" colspan="1">0.17</td><td align="left" rowspan="1" colspan="1">0.045</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>EigengeneSignif = a<sub>e,i</sub><sup>(q)</sup></italic>
</td><td align="left" rowspan="1" colspan="1">0.024</td><td align="left" rowspan="1" colspan="1">0.12</td><td align="left" rowspan="1" colspan="1">0.25</td><td align="left" rowspan="1" colspan="1">0.016</td><td align="left" rowspan="1" colspan="1">0.19</td><td align="left" rowspan="1" colspan="1">0.18</td><td align="left" rowspan="1" colspan="1">0.046</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pcbi-1000117-t003" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000117.t003</object-id><label>Table 3</label><caption><title>Robustness analysis of the brain cancer gene coexpression network results.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col width="187*" align="left" span="1"/><col width="88*" align="center" span="1"/><col width="89*" align="center" span="1"/><col width="88*" align="center" span="1"/><col width="89*" align="center" span="1"/><col width="88*" align="center" span="1"/><col width="89*" align="center" span="1"/><col width="88*" align="center" span="1"/><col width="89*" align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">
<italic/>
</td><td colspan="6" align="left" rowspan="1">Weighted networks</td><td colspan="2" align="left" rowspan="1">Unweighted networks</td></tr><tr><td align="left" rowspan="1" colspan="1">Squared correlation R 2 across modules</td><td colspan="6" align="left" rowspan="1">Soft threshold <italic>β</italic>
</td><td colspan="2" align="left" rowspan="1">Hard threshold <italic>τ</italic>
</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">Relation</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">4</td><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1">6</td><td align="left" rowspan="1" colspan="1">0.7</td><td align="left" rowspan="1" colspan="1">0.5</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>Centralization</italic>≈<italic>Centralization<sub>E</sub></italic>
</td><td align="left" rowspan="1" colspan="1">1.0</td><td align="left" rowspan="1" colspan="1">1.0</td><td align="left" rowspan="1" colspan="1">0.97</td><td align="left" rowspan="1" colspan="1">0.90</td><td align="left" rowspan="1" colspan="1">0.87</td><td align="left" rowspan="1" colspan="1">0.88</td><td align="left" rowspan="1" colspan="1">0.07</td><td align="left" rowspan="1" colspan="1">0.93</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>Heterogeneity</italic>≈<italic>Heterogeneity<sub>E</sub></italic>
</td><td align="left" rowspan="1" colspan="1">1.0</td><td align="left" rowspan="1" colspan="1">1.0</td><td align="left" rowspan="1" colspan="1">0.99</td><td align="left" rowspan="1" colspan="1">0.98</td><td align="left" rowspan="1" colspan="1">0.97</td><td align="left" rowspan="1" colspan="1">0.96</td><td align="left" rowspan="1" colspan="1">0.89</td><td align="left" rowspan="1" colspan="1">0.87</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>ClusterCoef<sub>i</sub></italic>≈<italic>ClusterCoef<sub>E</sub></italic>
</td><td align="left" rowspan="1" colspan="1">0.99</td><td align="left" rowspan="1" colspan="1">0.96</td><td align="left" rowspan="1" colspan="1">0.88</td><td align="left" rowspan="1" colspan="1">0.74</td><td align="left" rowspan="1" colspan="1">0.58</td><td align="left" rowspan="1" colspan="1">0.45</td><td align="left" rowspan="1" colspan="1">0.04</td><td align="left" rowspan="1" colspan="1">0.32</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>ModuleSignif</italic>≈<italic>ModuleSignif<sub>E</sub></italic>
</td><td align="left" rowspan="1" colspan="1">0.98</td><td align="left" rowspan="1" colspan="1">0.91</td><td align="left" rowspan="1" colspan="1">0.87</td><td align="left" rowspan="1" colspan="1">0.85</td><td align="left" rowspan="1" colspan="1">0.85</td><td align="left" rowspan="1" colspan="1">0.86</td><td align="left" rowspan="1" colspan="1">0.98</td><td align="left" rowspan="1" colspan="1">0.98</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>HubGeneSignif</italic>≈<italic>HubGeneSignif<sub>E</sub></italic>
</td><td align="left" rowspan="1" colspan="1">0.96</td><td align="left" rowspan="1" colspan="1">0.91</td><td align="left" rowspan="1" colspan="1">0.89</td><td align="left" rowspan="1" colspan="1">0.90</td><td align="left" rowspan="1" colspan="1">0.92</td><td align="left" rowspan="1" colspan="1">0.94</td><td align="left" rowspan="1" colspan="1">0.93</td><td align="left" rowspan="1" colspan="1">0.87</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>EigengeneSignif</italic>≈<italic>HubGeneSignif</italic>
</td><td align="left" rowspan="1" colspan="1">0.96</td><td align="left" rowspan="1" colspan="1">0.89</td><td align="left" rowspan="1" colspan="1">0.87</td><td align="left" rowspan="1" colspan="1">0.88</td><td align="left" rowspan="1" colspan="1">0.90</td><td align="left" rowspan="1" colspan="1">0.92</td><td align="left" rowspan="1" colspan="1">0.93</td><td align="left" rowspan="1" colspan="1">0.87</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>ClusterCoef</italic> = (1+<italic>Heterogeneity</italic>
<sup>2</sup>)<sup>2</sup>×<italic>Density</italic>
</td><td align="left" rowspan="1" colspan="1">0.99</td><td align="left" rowspan="1" colspan="1">0.96</td><td align="left" rowspan="1" colspan="1">0.89</td><td align="left" rowspan="1" colspan="1">0.76</td><td align="left" rowspan="1" colspan="1">0.61</td><td align="left" rowspan="1" colspan="1">0.49</td><td align="left" rowspan="1" colspan="1">0.006</td><td align="left" rowspan="1" colspan="1">0.32</td></tr><tr><td align="left" rowspan="1" colspan="1">

</td><td align="left" rowspan="1" colspan="1">1.0</td><td align="left" rowspan="1" colspan="1">0.99</td><td align="left" rowspan="1" colspan="1">0.99</td><td align="left" rowspan="1" colspan="1">0.98</td><td align="left" rowspan="1" colspan="1">0.97</td><td align="left" rowspan="1" colspan="1">0.95</td><td align="left" rowspan="1" colspan="1">0.85</td><td align="left" rowspan="1" colspan="1">0.99</td></tr><tr><td align="left" rowspan="1" colspan="1">

</td><td align="left" rowspan="1" colspan="1">0.90</td><td align="left" rowspan="1" colspan="1">0.68</td><td align="left" rowspan="1" colspan="1">0.058</td><td align="left" rowspan="1" colspan="1">0.016</td><td align="left" rowspan="1" colspan="1">0.16</td><td align="left" rowspan="1" colspan="1">0.35</td><td align="left" rowspan="1" colspan="1">0.20</td><td align="left" rowspan="1" colspan="1">1.0</td></tr><tr><td align="left" rowspan="1" colspan="1">

</td><td align="left" rowspan="1" colspan="1">0.94</td><td align="left" rowspan="1" colspan="1">0.94</td><td align="left" rowspan="1" colspan="1">0.94</td><td align="left" rowspan="1" colspan="1">0.94</td><td align="left" rowspan="1" colspan="1">0.93</td><td align="left" rowspan="1" colspan="1">0.92</td><td align="left" rowspan="1" colspan="1">0.95</td><td align="left" rowspan="1" colspan="1">0.98</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>K<sub>i</sub></italic>≈<italic>a<sub>e</sub></italic>
<sub>,<italic>i</italic></sub> (median <italic>R</italic>
<sup>2</sup>)</td><td align="left" rowspan="1" colspan="1">1.0</td><td align="left" rowspan="1" colspan="1">1.0</td><td align="left" rowspan="1" colspan="1">1.0</td><td align="left" rowspan="1" colspan="1">1.0</td><td align="left" rowspan="1" colspan="1">1.0</td><td align="left" rowspan="1" colspan="1">0.99</td><td align="left" rowspan="1" colspan="1">0.95</td><td align="left" rowspan="1" colspan="1">0.83</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pcbi-1000117-t004" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000117.t004</object-id><label>Table 4</label><caption><title>Values of network concepts in the different mouse gender/tissue networks reported in ##FIG##0##Figure 1##.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1"/><td colspan="2" align="left" rowspan="1">Female liver</td><td colspan="2" align="left" rowspan="1">Female adipose</td><td colspan="2" align="left" rowspan="1">Female brain</td><td colspan="2" align="left" rowspan="1">Female muscle</td></tr><tr><td align="left" rowspan="1" colspan="1">Network concept</td><td align="left" rowspan="1" colspan="1">Network</td><td align="left" rowspan="1" colspan="1">Eigengene</td><td align="left" rowspan="1" colspan="1">Network</td><td align="left" rowspan="1" colspan="1">Eigengene</td><td align="left" rowspan="1" colspan="1">Network</td><td align="left" rowspan="1" colspan="1">Eigengene</td><td align="left" rowspan="1" colspan="1">Network</td><td align="left" rowspan="1" colspan="1">Eigengene</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">Factorizability</td><td align="left" rowspan="1" colspan="1">0.92</td><td align="left" rowspan="1" colspan="1">0.91</td><td align="left" rowspan="1" colspan="1">0.72</td><td align="left" rowspan="1" colspan="1">0.46</td><td align="left" rowspan="1" colspan="1">0.89</td><td align="left" rowspan="1" colspan="1">0.82</td><td align="left" rowspan="1" colspan="1">0.79</td><td align="left" rowspan="1" colspan="1">0.68</td></tr><tr><td align="left" rowspan="1" colspan="1">Density</td><td align="left" rowspan="1" colspan="1">0.39</td><td align="left" rowspan="1" colspan="1">0.39</td><td align="left" rowspan="1" colspan="1">0.23</td><td align="left" rowspan="1" colspan="1">0.14</td><td align="left" rowspan="1" colspan="1">0.32</td><td align="left" rowspan="1" colspan="1">0.27</td><td align="left" rowspan="1" colspan="1">0.24</td><td align="left" rowspan="1" colspan="1">0.19</td></tr><tr><td align="left" rowspan="1" colspan="1">Centralization</td><td align="left" rowspan="1" colspan="1">0.19</td><td align="left" rowspan="1" colspan="1">0.19</td><td align="left" rowspan="1" colspan="1">0.11</td><td align="left" rowspan="1" colspan="1">0.19</td><td align="left" rowspan="1" colspan="1">0.34</td><td align="left" rowspan="1" colspan="1">0.23</td><td align="left" rowspan="1" colspan="1">0.17</td><td align="left" rowspan="1" colspan="1">0.22</td></tr><tr><td align="left" rowspan="1" colspan="1">Heterogeneity</td><td align="left" rowspan="1" colspan="1">0.18</td><td align="left" rowspan="1" colspan="1">0.19</td><td align="left" rowspan="1" colspan="1">0.22</td><td align="left" rowspan="1" colspan="1">0.59</td><td align="left" rowspan="1" colspan="1">0.36</td><td align="left" rowspan="1" colspan="1">0.54</td><td align="left" rowspan="1" colspan="1">0.32</td><td align="left" rowspan="1" colspan="1">0.57</td></tr><tr><td align="left" rowspan="1" colspan="1">Mean cluster coef</td><td align="left" rowspan="1" colspan="1">0.42</td><td align="left" rowspan="1" colspan="1">0.42</td><td align="left" rowspan="1" colspan="1">0.27</td><td align="left" rowspan="1" colspan="1">0.26</td><td align="left" rowspan="1" colspan="1">0.41</td><td align="left" rowspan="1" colspan="1">0.46</td><td align="left" rowspan="1" colspan="1">0.30</td><td align="left" rowspan="1" colspan="1">0.33</td></tr></tbody></table><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1"/><td colspan="2" align="left" rowspan="1">Male liver</td><td colspan="2" align="left" rowspan="1">Male adipose</td><td colspan="2" align="left" rowspan="1">Male brain</td><td colspan="2" align="left" rowspan="1">Male muscle</td></tr><tr><td align="left" rowspan="1" colspan="1">Network concept</td><td align="left" rowspan="1" colspan="1">Network</td><td align="left" rowspan="1" colspan="1">Eigengene</td><td align="left" rowspan="1" colspan="1">Network</td><td align="left" rowspan="1" colspan="1">Eigengene</td><td align="left" rowspan="1" colspan="1">Network</td><td align="left" rowspan="1" colspan="1">Eigengene</td><td align="left" rowspan="1" colspan="1">Network</td><td align="left" rowspan="1" colspan="1">Eigengene</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">Factorizability</td><td align="left" rowspan="1" colspan="1">0.93</td><td align="left" rowspan="1" colspan="1">0.92</td><td align="left" rowspan="1" colspan="1">0.76</td><td align="left" rowspan="1" colspan="1">0.59</td><td align="left" rowspan="1" colspan="1">0.73</td><td align="left" rowspan="1" colspan="1">0.46</td><td align="left" rowspan="1" colspan="1">0.76</td><td align="left" rowspan="1" colspan="1">0.48</td></tr><tr><td align="left" rowspan="1" colspan="1">Density</td><td align="left" rowspan="1" colspan="1">0.37</td><td align="left" rowspan="1" colspan="1">0.36</td><td align="left" rowspan="1" colspan="1">0.23</td><td align="left" rowspan="1" colspan="1">0.16</td><td align="left" rowspan="1" colspan="1">0.21</td><td align="left" rowspan="1" colspan="1">0.13</td><td align="left" rowspan="1" colspan="1">0.25</td><td align="left" rowspan="1" colspan="1">0.16</td></tr><tr><td align="left" rowspan="1" colspan="1">Centralization</td><td align="left" rowspan="1" colspan="1">0.19</td><td align="left" rowspan="1" colspan="1">0.21</td><td align="left" rowspan="1" colspan="1">0.15</td><td align="left" rowspan="1" colspan="1">0.21</td><td align="left" rowspan="1" colspan="1">0.21</td><td align="left" rowspan="1" colspan="1">0.20</td><td align="left" rowspan="1" colspan="1">0.12</td><td align="left" rowspan="1" colspan="1">0.21</td></tr><tr><td align="left" rowspan="1" colspan="1">Heterogeneity</td><td align="left" rowspan="1" colspan="1">0.28</td><td align="left" rowspan="1" colspan="1">0.32</td><td align="left" rowspan="1" colspan="1">0.28</td><td align="left" rowspan="1" colspan="1">0.65</td><td align="left" rowspan="1" colspan="1">0.27</td><td align="left" rowspan="1" colspan="1">0.67</td><td align="left" rowspan="1" colspan="1">0.30</td><td align="left" rowspan="1" colspan="1">0.64</td></tr><tr><td align="left" rowspan="1" colspan="1">Mean cluster coef</td><td align="left" rowspan="1" colspan="1">0.43</td><td align="left" rowspan="1" colspan="1">0.44</td><td align="left" rowspan="1" colspan="1">0.27</td><td align="left" rowspan="1" colspan="1">0.31</td><td align="left" rowspan="1" colspan="1">0.25</td><td align="left" rowspan="1" colspan="1">0.26</td><td align="left" rowspan="1" colspan="1">0.31</td><td align="left" rowspan="1" colspan="1">0.31</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pcbi-1000117-t005" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000117.t005</object-id><label>Table 5</label><caption><title>Robustness analysis of the mouse coexpression network.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col width="208*" align="left" span="1"/><col width="85*" align="center" span="1"/><col width="86*" align="center" span="1"/><col width="85*" align="center" span="1"/><col width="86*" align="center" span="1"/><col width="85*" align="center" span="1"/><col width="86*" align="center" span="1"/><col width="85*" align="center" span="1"/><col width="86*" align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">
<italic/>
</td><td colspan="6" align="left" rowspan="1">Weighted networks</td><td colspan="2" align="left" rowspan="1">Unweighted networks</td></tr><tr><td align="left" rowspan="1" colspan="1">Squared correlation <italic>R<sup>2</sup></italic> across modules</td><td colspan="6" align="left" rowspan="1">Soft threshold <italic>β</italic>
</td><td colspan="2" align="left" rowspan="1">Hard threshold <italic>τ</italic>
</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">Relation</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">4</td><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1">6</td><td align="left" rowspan="1" colspan="1">0.65</td><td align="left" rowspan="1" colspan="1">0.5</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>Centralization</italic>≈<italic>Centralization<sub>E</sub></italic>
</td><td align="left" rowspan="1" colspan="1">0.69</td><td align="left" rowspan="1" colspan="1">0.74</td><td align="left" rowspan="1" colspan="1">0.90</td><td align="left" rowspan="1" colspan="1">0.95</td><td align="left" rowspan="1" colspan="1">0.94</td><td align="left" rowspan="1" colspan="1">0.92</td><td align="left" rowspan="1" colspan="1">0.007</td><td align="left" rowspan="1" colspan="1">0.66</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>Heterogeneity</italic>≈<italic>Heterogeneity<sub>E</sub></italic>
</td><td align="left" rowspan="1" colspan="1">0.54</td><td align="left" rowspan="1" colspan="1">0.59</td><td align="left" rowspan="1" colspan="1">0.71</td><td align="left" rowspan="1" colspan="1">0.82</td><td align="left" rowspan="1" colspan="1">0.88</td><td align="left" rowspan="1" colspan="1">0.86</td><td align="left" rowspan="1" colspan="1">0.30</td><td align="left" rowspan="1" colspan="1">0.33</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>ClusterCoef<sub>i</sub></italic>≈<italic>ClusterCoef<sub>E</sub></italic>
</td><td align="left" rowspan="1" colspan="1">0.94</td><td align="left" rowspan="1" colspan="1">0.84</td><td align="left" rowspan="1" colspan="1">0.70</td><td align="left" rowspan="1" colspan="1">0.59</td><td align="left" rowspan="1" colspan="1">0.50</td><td align="left" rowspan="1" colspan="1">0.44</td><td align="left" rowspan="1" colspan="1">0.09</td><td align="left" rowspan="1" colspan="1">0.33</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>ModuleSignif</italic>≈<italic>ModuleSignif<sub>E</sub></italic>
</td><td align="left" rowspan="1" colspan="1">0.96</td><td align="left" rowspan="1" colspan="1">0.96</td><td align="left" rowspan="1" colspan="1">0.96</td><td align="left" rowspan="1" colspan="1">0.97</td><td align="left" rowspan="1" colspan="1">0.98</td><td align="left" rowspan="1" colspan="1">0.99</td><td align="left" rowspan="1" colspan="1">0.96</td><td align="left" rowspan="1" colspan="1">0.96</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>HubGeneSignif</italic>≈<italic>HubGeneSignif<sub>E</sub></italic>
</td><td align="left" rowspan="1" colspan="1">0.98</td><td align="left" rowspan="1" colspan="1">0.98</td><td align="left" rowspan="1" colspan="1">0.98</td><td align="left" rowspan="1" colspan="1">0.99</td><td align="left" rowspan="1" colspan="1">1.0</td><td align="left" rowspan="1" colspan="1">1.0.</td><td align="left" rowspan="1" colspan="1">0.88</td><td align="left" rowspan="1" colspan="1">0.91</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>EigengeneSignif</italic>≈<italic>HubGeneSignif</italic>
</td><td align="left" rowspan="1" colspan="1">0.98</td><td align="left" rowspan="1" colspan="1">0.98</td><td align="left" rowspan="1" colspan="1">0.98</td><td align="left" rowspan="1" colspan="1">0.99</td><td align="left" rowspan="1" colspan="1">1.0</td><td align="left" rowspan="1" colspan="1">1.0</td><td align="left" rowspan="1" colspan="1">0.89</td><td align="left" rowspan="1" colspan="1">0.92</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>ClusterCoef</italic> = (1+<italic>Heterogeneity</italic>
<sup>2</sup>)<sup>2</sup>×<italic>Density</italic>
</td><td align="left" rowspan="1" colspan="1">0.89</td><td align="left" rowspan="1" colspan="1">0.78</td><td align="left" rowspan="1" colspan="1">0.70</td><td align="left" rowspan="1" colspan="1">0.62</td><td align="left" rowspan="1" colspan="1">0.54</td><td align="left" rowspan="1" colspan="1">0.48</td><td align="left" rowspan="1" colspan="1">0.08</td><td align="left" rowspan="1" colspan="1">0.31</td></tr><tr><td align="left" rowspan="1" colspan="1">

</td><td align="left" rowspan="1" colspan="1">0.99</td><td align="left" rowspan="1" colspan="1">0.99</td><td align="left" rowspan="1" colspan="1">0.99</td><td align="left" rowspan="1" colspan="1">0.99</td><td align="left" rowspan="1" colspan="1">0.99</td><td align="left" rowspan="1" colspan="1">0.99</td><td align="left" rowspan="1" colspan="1">0.90</td><td align="left" rowspan="1" colspan="1">0.96</td></tr><tr><td align="left" rowspan="1" colspan="1">

</td><td align="left" rowspan="1" colspan="1">0.52</td><td align="left" rowspan="1" colspan="1">0.21</td><td align="left" rowspan="1" colspan="1">0.43</td><td align="left" rowspan="1" colspan="1">0.73</td><td align="left" rowspan="1" colspan="1">0.82</td><td align="left" rowspan="1" colspan="1">0.84</td><td align="left" rowspan="1" colspan="1">0.60</td><td align="left" rowspan="1" colspan="1">0.82</td></tr><tr><td align="left" rowspan="1" colspan="1">

</td><td align="left" rowspan="1" colspan="1">0.95</td><td align="left" rowspan="1" colspan="1">0.97</td><td align="left" rowspan="1" colspan="1">0.97</td><td align="left" rowspan="1" colspan="1">0.98</td><td align="left" rowspan="1" colspan="1">0.98</td><td align="left" rowspan="1" colspan="1">0.98</td><td align="left" rowspan="1" colspan="1">0.93</td><td align="left" rowspan="1" colspan="1">0.80</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>K<sub>i</sub></italic>≈<italic>a<sub>e</sub></italic>
<sub>,<italic>i</italic></sub> (median <italic>R</italic>
<sup>2</sup>)</td><td align="left" rowspan="1" colspan="1">1.0</td><td align="left" rowspan="1" colspan="1">0.99</td><td align="left" rowspan="1" colspan="1">0.98</td><td align="left" rowspan="1" colspan="1">0.96</td><td align="left" rowspan="1" colspan="1">0.95</td><td align="left" rowspan="1" colspan="1">0.94</td><td align="left" rowspan="1" colspan="1">0.74</td><td align="left" rowspan="1" colspan="1">0.86</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pcbi-1000117-t006" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000117.t006</object-id><label>Table 6</label><caption><title>Robustness analysis of the yeast coexpression network.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col width="194*" align="left" span="1"/><col width="88*" align="center" span="1"/><col width="90*" align="center" span="1"/><col width="88*" align="center" span="1"/><col width="90*" align="center" span="1"/><col width="88*" align="center" span="1"/><col width="90*" align="center" span="1"/><col width="88*" align="center" span="1"/><col width="90*" align="center" span="1"/><col width="88*" align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1"/><td colspan="7" align="left" rowspan="1">Weighted networks</td><td colspan="2" align="left" rowspan="1">Unweighted networks</td></tr><tr><td align="left" rowspan="1" colspan="1">Squared correlation <italic>R<sup>2</sup></italic> across modules</td><td colspan="7" align="left" rowspan="1">Soft threshold <italic>β</italic>
</td><td colspan="2" align="left" rowspan="1">Hard threshold <italic>τ</italic>
</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">Relation</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">4</td><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1">6</td><td align="left" rowspan="1" colspan="1">7</td><td align="left" rowspan="1" colspan="1">0.65</td><td align="left" rowspan="1" colspan="1">0.5</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>Centralization</italic>≈<italic>Centralization<sub>E</sub></italic>
</td><td align="left" rowspan="1" colspan="1">0.99</td><td align="left" rowspan="1" colspan="1">0.97</td><td align="left" rowspan="1" colspan="1">0.97</td><td align="left" rowspan="1" colspan="1">0.98</td><td align="left" rowspan="1" colspan="1">0.98</td><td align="left" rowspan="1" colspan="1">0.98</td><td align="left" rowspan="1" colspan="1">0.98</td><td align="left" rowspan="1" colspan="1">0.53</td><td align="left" rowspan="1" colspan="1">0.60</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>Heterogeneity</italic>≈<italic>Heterogeneity<sub>E</sub></italic>
</td><td align="left" rowspan="1" colspan="1">0.86</td><td align="left" rowspan="1" colspan="1">0.92</td><td align="left" rowspan="1" colspan="1">0.94</td><td align="left" rowspan="1" colspan="1">0.94</td><td align="left" rowspan="1" colspan="1">0.93</td><td align="left" rowspan="1" colspan="1">0.92</td><td align="left" rowspan="1" colspan="1">0.91</td><td align="left" rowspan="1" colspan="1">0.13</td><td align="left" rowspan="1" colspan="1">0.006</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>ClusterCoef<sub>i</sub></italic>≈<italic>ClusterCoef<sub>E</sub></italic>
</td><td align="left" rowspan="1" colspan="1">0.98</td><td align="left" rowspan="1" colspan="1">0.97</td><td align="left" rowspan="1" colspan="1">0.94</td><td align="left" rowspan="1" colspan="1">0.92</td><td align="left" rowspan="1" colspan="1">0.89</td><td align="left" rowspan="1" colspan="1">0.86</td><td align="left" rowspan="1" colspan="1">0.82</td><td align="left" rowspan="1" colspan="1">0.18</td><td align="left" rowspan="1" colspan="1">0.25</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>ClusterCoef</italic> = (1+<italic>Heterogeneity</italic>
<sup>2</sup>)<sup>2</sup>×<italic>Density</italic>
</td><td align="left" rowspan="1" colspan="1">0.99</td><td align="left" rowspan="1" colspan="1">0.97</td><td align="left" rowspan="1" colspan="1">0.95</td><td align="left" rowspan="1" colspan="1">0.92</td><td align="left" rowspan="1" colspan="1">0.89</td><td align="left" rowspan="1" colspan="1">0.86</td><td align="left" rowspan="1" colspan="1">0.83</td><td align="left" rowspan="1" colspan="1">0.21</td><td align="left" rowspan="1" colspan="1">0.27</td></tr><tr><td align="left" rowspan="1" colspan="1">

</td><td align="left" rowspan="1" colspan="1">1.0</td><td align="left" rowspan="1" colspan="1">0.99</td><td align="left" rowspan="1" colspan="1">0.98</td><td align="left" rowspan="1" colspan="1">0.98</td><td align="left" rowspan="1" colspan="1">0.97</td><td align="left" rowspan="1" colspan="1">0.95</td><td align="left" rowspan="1" colspan="1">0.94</td><td align="left" rowspan="1" colspan="1">0.99</td><td align="left" rowspan="1" colspan="1">0.99</td></tr><tr><td align="left" rowspan="1" colspan="1">

</td><td align="left" rowspan="1" colspan="1">0.51</td><td align="left" rowspan="1" colspan="1">0.24</td><td align="left" rowspan="1" colspan="1">0.04</td><td align="left" rowspan="1" colspan="1">0.06</td><td align="left" rowspan="1" colspan="1">0.33</td><td align="left" rowspan="1" colspan="1">0.53</td><td align="left" rowspan="1" colspan="1">0.68</td><td align="left" rowspan="1" colspan="1">0.76</td><td align="left" rowspan="1" colspan="1">0.98</td></tr><tr><td align="left" rowspan="1" colspan="1">

</td><td align="left" rowspan="1" colspan="1">0.89</td><td align="left" rowspan="1" colspan="1">0.91</td><td align="left" rowspan="1" colspan="1">0.93</td><td align="left" rowspan="1" colspan="1">0.94</td><td align="left" rowspan="1" colspan="1">0.95</td><td align="left" rowspan="1" colspan="1">0.96</td><td align="left" rowspan="1" colspan="1">0.96</td><td align="left" rowspan="1" colspan="1">0.51</td><td align="left" rowspan="1" colspan="1">0.20</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>K<sub>i</sub></italic>≈<italic>a<sub>e</sub></italic>
<sub>,<italic>i</italic></sub> (median <italic>R</italic>
<sup>2</sup>)</td><td align="left" rowspan="1" colspan="1">1.0</td><td align="left" rowspan="1" colspan="1">0.99</td><td align="left" rowspan="1" colspan="1">0.99</td><td align="left" rowspan="1" colspan="1">0.98</td><td align="left" rowspan="1" colspan="1">0.97</td><td align="left" rowspan="1" colspan="1">0.97</td><td align="left" rowspan="1" colspan="1">0.96</td><td align="left" rowspan="1" colspan="1">0.93</td><td align="left" rowspan="1" colspan="1">0.92</td></tr></tbody></table></alternatives></table-wrap>
<disp-formula></disp-formula> <disp-formula><label>(1)</label></disp-formula> <disp-formula><label>(2)</label></disp-formula> <disp-formula><label>(3)</label></disp-formula> <disp-formula><label>(4)</label></disp-formula> <disp-formula><label>(5)</label></disp-formula> <disp-formula><label>(6)</label></disp-formula> <disp-formula><label>(7)</label></disp-formula> <disp-formula><label>(8)</label></disp-formula> <disp-formula><label>(9)</label></disp-formula> <disp-formula><label>(10)</label></disp-formula> <disp-formula><label>(11)</label></disp-formula> <disp-formula><label>(12)</label></disp-formula> <disp-formula><label>(13)</label></disp-formula> <disp-formula><label>(14)</label></disp-formula> <disp-formula><label>(15)</label></disp-formula> <disp-formula><label>(16)</label></disp-formula> <disp-formula><label>(17)</label></disp-formula> <disp-formula><label>(18)</label></disp-formula> <disp-formula><label>(19)</label></disp-formula> <disp-formula><label>(20)</label></disp-formula> <disp-formula><label>(21)</label></disp-formula> <disp-formula><label>(22)</label></disp-formula> <disp-formula><label>(23)</label></disp-formula> <inline-formula></inline-formula> <disp-formula><label>(24)</label></disp-formula> <disp-formula></disp-formula> <disp-formula><label>(25)</label></disp-formula> <disp-formula><label>(26)</label></disp-formula> <disp-formula><label>(27)</label></disp-formula> <disp-formula></disp-formula> <disp-formula></disp-formula> <disp-formula></disp-formula> <disp-formula></disp-formula> <disp-formula></disp-formula> <disp-formula></disp-formula> <disp-formula></disp-formula> <disp-formula></disp-formula> <disp-formula></disp-formula> <disp-formula></disp-formula> <disp-formula></disp-formula> <disp-formula></disp-formula> <disp-formula></disp-formula> <disp-formula></disp-formula> <disp-formula></disp-formula> <disp-formula></disp-formula> <disp-formula><label>(28)</label></disp-formula> <disp-formula><label>(29)</label></disp-formula> <disp-formula><label>(30)</label></disp-formula> <inline-formula></inline-formula> <disp-formula><label>(31)</label></disp-formula> <disp-formula><label>(32)</label></disp-formula> <disp-formula><label>(33)</label></disp-formula> <disp-formula><label>(34)</label></disp-formula> <disp-formula><label>(35)</label></disp-formula> <disp-formula><label>(36)</label></disp-formula> <disp-formula><label>(37)</label></disp-formula> <disp-formula><label>(38)</label></disp-formula> <disp-formula><label>(39)</label></disp-formula> <disp-formula><label>(40)</label></disp-formula> <disp-formula><label>(41)</label></disp-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <disp-formula><label>(42)</label></disp-formula> <disp-formula></disp-formula> <disp-formula><label>(43)</label></disp-formula> <disp-formula><label>(44)</label></disp-formula> <disp-formula><label>(45)</label></disp-formula> <disp-formula><label>(46)</label></disp-formula> <disp-formula><label>(47)</label></disp-formula> <inline-formula></inline-formula> <disp-formula></disp-formula> <disp-formula></disp-formula> <disp-formula><label>(48)</label></disp-formula> <disp-formula></disp-formula> <disp-formula></disp-formula> <disp-formula></disp-formula> <disp-formula><label>(49)</label></disp-formula> <disp-formula></disp-formula> <inline-formula></inline-formula> <disp-formula><label>(50)</label></disp-formula> <disp-formula><label>(51)</label></disp-formula> <inline-formula></inline-formula> <disp-formula><label>(52)</label></disp-formula> <disp-formula><label>(53)</label></disp-formula> <disp-formula><label>(54)</label></disp-formula> <disp-formula><label>(55)</label></disp-formula> <disp-formula></disp-formula> <disp-formula></disp-formula> <disp-formula></disp-formula>
<supplementary-material content-type="local-data" id="pcbi.1000117.s001"><label>Text S1</label><caption><p>Robustness Analysis of the Brain Cancer Gene Coexpression Network. This supporting text provides a detailed analysis of the brain cancer gene coexpression network. The robustness analysis illustrates how the results change with regard to different network construction methods.</p><p>(3.83 MB PDF)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000117.s002"><label>Text S2</label><caption><p>Robustness Analysis of the Mouse Gene Coexpression Network. This supporting text provides a detailed analysis of the mouse tissue gene coexpression network. The robustness analysis illustrates how the results change with regard to different network construction methods.</p><p>(3.76 MB PDF)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000117.s003"><label>Text S3</label><caption><p>Robustness Analysis of the Yeast Gene Coexpression Network. This supporting text provides a detailed analysis of the yeast cell cycle gene coexpression network. The robustness analysis illustrates how the results change with regard to different network construction methods.</p><p>(2.62 MB PDF)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000117.s004"><label>Text S4</label><caption><p>Brain Cancer Network Comprised of 500 Prognostic Genes. Here we analyze a brain cancer network comprised of the 500 genes with highest absolute correlation with brain cancer survival time. The results illustrate that our theoretical results also apply to small networks comprised of sample trait related genes. The robustness analysis illustrates how the results change with regard to different network construction methods.</p><p>(0.38 MB PDF)</p></caption></supplementary-material>
<table-wrap-foot><fn id="nt101"><p>Here we omit the reference to the <italic>q</italic>th module.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt102"><p>Here we report the results for soft thresholding with <italic>β</italic> = 1 (Equation 2). The results for higher powers <italic>β</italic> and for unweighted networks can be found in ##SUPPL##0##Text S1##. Grey colors genes outside the 6 properly defined modules. The table shows that network concepts in the proper modules are close to their eigengene based analogs.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt103"><p>The table reports how the relationships among network concepts change as function of different soft threshold parameters <italic>β</italic> (Equation 2) or hard thresholds (Equation 1) used in the network construction. For each relationship and each network construction method, the table entry reports the squared correlation <italic>R</italic>
<sup>2</sup> across the proper modules. For within module comparisons the table reports median <italic>R</italic>
<sup>2</sup> values. Additional details can be found in ##SUPPL##0##Text S1##.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt104"><p>For each network, the table reports the network factorizability <italic>F</italic>(<italic>A</italic>), the eigengene factorizability <italic>EF</italic>(<italic>X</italic>), network concepts, and their eigengene-based analogs. Here we use a soft threshold <italic>β</italic> = 1 (Equation 2).</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt105"><p>The table reports how the relationships among network concepts change as function of different soft threshold parameters <italic>β</italic> (Equation 2) or hard thresholds <italic>τ</italic> (Equation 1) used in the network construction. For each relationship and each network construction method, the table entry reports the squared correlation <italic>R</italic>
<sup>2</sup> across the proper modules. For within module comparisons the table reports median <italic>R</italic>
<sup>2</sup> values. Additional details can be found in ##SUPPL##1##Text S2##.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt106"><p>The table reports how the relationships among network concepts change as function of different soft threshold parameters <italic>β</italic> (Equation 2) or hard thresholds <italic>τ</italic> (Equation 1) used in the network construction. For each relationship and each network construction method, the table entry reports the squared correlation <italic>R</italic>
<sup>2</sup> across the proper modules. For within module comparisons the table reports median <italic>R</italic>
<sup>2</sup> values. Additional details can be found in ##SUPPL##2##Text S3##.</p></fn></table-wrap-foot> <fn-group><fn fn-type="COI-statement"><p>The authors have declared that no competing interests exist.</p></fn><fn fn-type="financial-disclosure"><p>We acknowledge grant support from 1U19AI063603-01, P50CA092131, 1U24NS043562-01, 5P30CA016042-28, and HL28481.</p></fn></fn-group>
<graphic xlink:href="pcbi.1000117.e001.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e002.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e003.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e004.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e005.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.g001"/> <graphic xlink:href="pcbi.1000117.e006.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e007.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e008.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e009.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.g002"/> <graphic xlink:href="pcbi.1000117.e010.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e011.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e012.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e013.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e014.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.g003"/> <graphic xlink:href="pcbi.1000117.e015.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e016.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e017.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e018.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e019.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e020.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e021.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e022.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e023.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.g004"/> <graphic xlink:href="pcbi.1000117.e024.jpg" mimetype="image" position="float"/> <inline-graphic xlink:href="pcbi.1000117.e025.jpg" mimetype="image"/> <graphic xlink:href="pcbi.1000117.e026.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e027.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e028.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e029.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e030.jpg" mimetype="image" position="float"/> <graphic id="pcbi-1000117-t001-1" xlink:href="pcbi.1000117.t001"/> <graphic xlink:href="pcbi.1000117.e031"/> <graphic xlink:href="pcbi.1000117.e032"/> <graphic xlink:href="pcbi.1000117.e033"/> <graphic xlink:href="pcbi.1000117.e034"/> <graphic xlink:href="pcbi.1000117.e035"/> <graphic xlink:href="pcbi.1000117.e036"/> <graphic xlink:href="pcbi.1000117.e037"/> <graphic xlink:href="pcbi.1000117.e038"/> <graphic xlink:href="pcbi.1000117.e039"/> <graphic xlink:href="pcbi.1000117.e040"/> <graphic xlink:href="pcbi.1000117.e041"/> <graphic xlink:href="pcbi.1000117.e042"/> <graphic xlink:href="pcbi.1000117.e043"/> <graphic xlink:href="pcbi.1000117.e044"/> <graphic xlink:href="pcbi.1000117.e045"/> <graphic xlink:href="pcbi.1000117.e046"/> <graphic xlink:href="pcbi.1000117.e047.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e048.jpg" mimetype="image" position="float"/> <graphic id="pcbi-1000117-t002-2" xlink:href="pcbi.1000117.t002"/> <graphic xlink:href="pcbi.1000117.g005"/> <graphic xlink:href="pcbi.1000117.e049.jpg" mimetype="image" position="float"/> <inline-graphic xlink:href="pcbi.1000117.e050.jpg" mimetype="image"/> <graphic xlink:href="pcbi.1000117.g006"/> <graphic xlink:href="pcbi.1000117.e051.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e052.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e053.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e054.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e055.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e056.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e057.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e058.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e059.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e060.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e061.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.g007"/> <graphic id="pcbi-1000117-t003-3" xlink:href="pcbi.1000117.t003"/> <inline-graphic xlink:href="pcbi.1000117.e062.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000117.e063.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000117.e064.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000117.e065.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000117.e066.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000117.e067.jpg" mimetype="image"/> <graphic xlink:href="pcbi.1000117.g008"/> <inline-graphic xlink:href="pcbi.1000117.e068.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000117.e069.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000117.e070.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000117.e071.jpg" mimetype="image"/> <graphic xlink:href="pcbi.1000117.g009"/> <graphic id="pcbi-1000117-t004-4" xlink:href="pcbi.1000117.t004"/> <graphic id="pcbi-1000117-t005-5" xlink:href="pcbi.1000117.t005"/> <inline-graphic xlink:href="pcbi.1000117.e072.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000117.e073.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000117.e074.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000117.e075.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000117.e076.jpg" mimetype="image"/> <graphic id="pcbi-1000117-t006-6" xlink:href="pcbi.1000117.t006"/> <inline-graphic xlink:href="pcbi.1000117.e077.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000117.e078.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000117.e079.jpg" mimetype="image"/> <graphic xlink:href="pcbi.1000117.e080.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e081.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e082.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e083.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e084.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e085.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e086.jpg" mimetype="image" position="float"/> <inline-graphic xlink:href="pcbi.1000117.e087.jpg" mimetype="image"/> <graphic xlink:href="pcbi.1000117.e088.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e089.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e090.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e091.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e092.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e093.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e094.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e095.jpg" mimetype="image" position="float"/> <inline-graphic xlink:href="pcbi.1000117.e096.jpg" mimetype="image"/> <graphic xlink:href="pcbi.1000117.e097.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e098.jpg" mimetype="image" position="float"/> <inline-graphic xlink:href="pcbi.1000117.e099.jpg" mimetype="image"/> <graphic xlink:href="pcbi.1000117.e100.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e101.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e102.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e103.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e104.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e105.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000117.e106.jpg" mimetype="image" position="float"/>
<media xlink:href="pcbi.1000117.s001.pdf"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000117.s002.pdf"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000117.s003.pdf"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000117.s004.pdf"><caption><p>Click here for additional data file.</p></caption></media>
[{"label": ["19"], "element-citation": ["\n"], "surname": ["Zhang", "Horvath"], "given-names": ["B", "S"], "year": ["2005"], "article-title": ["A general framework for weighted gene co-expression network analysis."], "source": ["Stat Appl Genet Mol Biol"], "volume": ["4"], "fpage": ["17"]}, {"label": ["23"], "element-citation": ["\n"], "surname": ["Ghazalpour", "Doss", "Zhang", "Plaisier", "Wang"], "given-names": ["A", "S", "B", "C", "S"], "year": ["2006"], "article-title": ["Integrating genetics and network analysis to characterize genes related to mouse weight."], "source": ["PloS Genet"], "volume": ["2"], "fpage": ["8"], "comment": ["doi:10.1371/journal.pgen.0020130"]}, {"label": ["31"], "element-citation": ["\n"], "surname": ["Albert", "Barabasi"], "given-names": ["R", "A"], "year": ["2002"], "article-title": ["Statistical mechanics of complex networks."], "source": ["Rev Mod Phys"], "volume": ["74"], "fpage": ["47"], "lpage": ["97"]}, {"label": ["35"], "element-citation": ["\n"], "surname": ["Snijders"], "given-names": ["T"], "year": ["1981"], "article-title": ["The degree variance: an index of graph heterogeneity."], "source": ["Soc Networks"], "volume": ["3"], "fpage": ["163"], "lpage": ["174"]}, {"label": ["36"], "element-citation": ["\n"], "surname": ["Freeman"], "given-names": ["L"], "year": ["1978"], "article-title": ["Centrality in social networks: conceptual clarification."], "source": ["Soc Networks"], "volume": ["1"], "fpage": ["215"], "lpage": ["239"]}, {"label": ["41"], "element-citation": ["\n"], "surname": ["Ye", "Godzik"], "given-names": ["Y", "A"], "year": ["2004"], "article-title": ["Comparative analysis of protein domain organization."], "source": ["Genome Biol"], "volume": ["14"], "fpage": ["343"], "lpage": ["353"]}, {"label": ["50"], "element-citation": ["\n"], "surname": ["Adrian", "Chris", "Beatrix", "Joseph", "Guang"], "given-names": ["D", "H", "J", "R", "Y"], "year": ["2004"], "article-title": ["Sparse graphical models for exploring gene expression data."], "source": ["J Multivar Anal"], "volume": ["90"], "fpage": ["196"], "lpage": ["212"]}, {"label": ["53"], "element-citation": ["\n"], "surname": ["Fisher"], "given-names": ["RA"], "year": ["1915"], "article-title": ["On the \u2018probable error\u2019 of a coefficient of correlation deduced from a small sample."], "source": ["Metron"], "volume": ["1"], "fpage": ["1"], "lpage": ["32"]}, {"label": ["56"], "element-citation": ["\n"], "surname": ["Perkins", "Jaeger", "Reinitz", "Glass"], "given-names": ["TJ", "J", "J", "L"], "year": ["2005"], "article-title": ["Reverse engineering the gap gene network of Drosophila melanogaster."], "source": ["PLoS Comput Biol"], "volume": ["2"], "fpage": ["e51"], "comment": ["doi:10.1371/journal.pcbi.0020051"]}, {"label": ["60"], "element-citation": ["\n"], "surname": ["Price", "Dehal", "Arkin"], "given-names": ["MN", "PS", "AP"], "year": ["2007"], "article-title": ["Orthologous transcription factors in bacteria have different functions and regulate different genes."], "source": ["PLoS Comput Biol"], "volume": ["3"], "fpage": ["e175"], "comment": ["doi:10.1371/journal.pcbi.0030175"]}]
PLoS Comput Biol. 2008 Aug 15; 4(8):e1000117