Spreading of diseases through comorbidity networks across life and gender

We use a database of the Main Association of Austrian Social Security Institutions that contains pseudonymized claims data of all persons receiving outpatient and inpatient care in Austria between 1 January 2006 and 31 December 2007. The data provides a comprehensive, nationwide collection of the medical condition of the vast majority of Austriaʼs population of 8.3 million people. Information on diagnoses is available for all persons receiving inpatient care, so-called inpatients. The total sample of inpatients consists of 1,862,258 patients (1,064,952 females and 797,306 males) and includes their year of birth, gender, date of death, their drug prescriptions, and main- and side-diagnoses. Patients are grouped according to their age. The age groups, labeled by a, contain all patients whose age is between a and $a+8$ years. Diagnosis are provided in the 10th revision of the International Statistical Classification of Diseases and Related Health Problems (ICD 10) [24], a medical classification list by the World Health Organization (WHO). Not all ICD codes represent disorders—they may also indicate general examinations, injuries, pregnancies, or collections of symptoms that could be indicative of a large number of diseases. We exclude these categories and work with the D = 1,055 diagnoses on the three-digit ICD levels in chapters A–N. We will use the words disease, disorder, and diagnosis interchangeably whenever referring to a specific ICD entry. It was found that such data can generally be considered to be complete with only a small amount of gross diagnostic miscoding [25]. Prognostic models for hospitalizations built on comorbidity scores derived from medical claims data perform similarly well as survey-derived data [26].

We use the following notation. In our data we have ${{N}^{m/f}}(a)$ men/women in age group a. Let us consider index pairs, where index i indicates patients suffering disease i and $\neg i$ indicates patients not having the diagnosis i. Each of the ${{N}^{m/f}}(a)$ patients belongs to one of four index pairs, $\{(i,j),(\neg i,j),(i,\neg j),(\neg i,\neg j)\}$, meaning that the patient has diseases i and j, has j but not i, has i but not j, or has neither i nor j, respectively. The numbers $N_{i,j}^{m/f}(a),N_{i,\neg j}^{m/f}(a),N_{\neg i,j}^{m/f}(a),N_{\neg i,\neg j}^{m/f}(a)$ denote the number of male/female patients in age group a in these respective groups. The prevalence for each disease i is defined as $p_{i}^{m/f}(a)=\frac{N_{i}^{m/f}(a)}{{{N}^{m/f}}(a)}$. In the following, we sometimes suppress the indices for age and gender for clarity.

2.2.1. Layer 1: statistical significance of comorbidity relations

The first layer of the PHDN multiplex provides the statistical significance of relations between any two diseases i and j. With the above notation the contingency coefficient ${{\phi }_{ij}}$ for a specific age group is defined by

$$\begin{eqnarray}&&{{\phi }_{ij}}=\frac{{{N}_{i,j}}{{N}_{\neg i,\neg j}}-{{N}_{i,\neg j}}{{N}_{\neg i,j}}}{\sqrt{{{N}_{i}}{{N}_{\neg i}}{{N}_{j}}{{N}_{\neg j}}}}.\end{eqnarray} \tag{ 1 }$$

A value of ${{\phi }_{ij}}(a)$ is computed for the patient numbers in each age group a. These values are directly related to the chi-square-coefficient ${{({{\chi }_{ij}})}^{2}}(a)=N(a){{({{\phi }_{ij(a)}})}^{2}}$, where N(a) is the number of patients in the considered age group. A chi-square test can be used to reject the null hypothesis for a given confidence level that the occurrence of disease i is independent of the occurrence of disease j in the same patient. $\phi =({{\phi }_{ij}})$ is a symmetric matrix containing the test statistics for diseases i and j. We discard negative entries by setting ${{\phi }_{ij}}\to 0$, whenever ${{\phi }_{ij}}\lt 0.$ For later use we define an unweighted adjacency matrix J derived from ϕ by defining ${{J}_{ij}}=1$ if ${{\phi }_{ij}}\gt 0,$ and ${{J}_{ij}}=0$ otherwise.

There are known biases in the test statistic ${{\phi }_{ij}}$ that depend on the prevalences $p_{i}^{m/f}(a)$ and $p_{j}^{m/f}(a)$ of the diseases i and j. While the contingency coefficient is a reliable measure for correlations whenever diseases i and j have similar prevalences $p_{i}^{m/f}(a)$ and $p_{j}^{m/f}(a)$, ϕ underestimates correlations between very rare and highly frequent diseases [11, 19]. It is therefore unclear if low values for ${{\phi }_{ij}}$ for two diseases i and j imply weak correlations or large differences in their respective prevalences.

Assume that disease i has a very low prevalence $p_{i}^{m/f}(a)$ with respect to the prevalence values of the majority of other diseases. Then each of its ϕ-values can be assumed to underestimate the true correlations, with respect to similar ϕ-values obtained from more frequent diseases. A simple way to correct these biases is to compare values of ${{\phi }_{ij}}$ for disease i with the typical correlation strength for disease i [20]. To this end, consider the matrix

$$\begin{eqnarray}&&\tilde{{{\phi }_{ij}}}=\frac{{{\phi }_{ij}}}{{{\langle \phi \rangle }_{i}}},\end{eqnarray} \tag{ 2 }$$

where ${{\langle \phi \rangle }_{i}}=\frac{\sum _{j=1}^{D}{{\phi }_{ij}}}{\sum _{j=1}^{D}{{J}_{ij}}}$. While ϕ is an undirected network, $\tilde{\phi }$ is typically directed. To obtain a unique and bias-corrected measure for the strength of the correlation between diseases i and j, an undirected version ${{\phi }^{*}}$ of the network $\tilde{{{\phi }_{ij}}}$ is obtained by $\phi _{ij}^{*}={\rm max} \left\{ \tilde{{{\phi }_{ij}}},\tilde{{{\phi }_{ji}}} \right\}$.

2.2.2. Layer 2: conditional disease risk

The second layer of the PHDN multiplex is given by the risk for a patient of a certain age group and gender to have disease i, given that she has disease j. The conditional probability $P(i|j)$ for a patient who suffers disease j to also suffer disease i is given by $P(i|j)=\frac{{{N}_{i,j}}}{{{N}_{j}}}$. A thresholded, undirected, and unweighted version P^* of $P(i|j)$ is

$$\begin{eqnarray}P_{ij}^{*}=\left\{ \begin{array}{ccccccccccccccc} 1 & {} & {\rm max} \left\{ {{P}_{ij}},{{P}_{ji}} \right\}\gt {{p}_{c}} \\ 0 & {} & {\rm otherwise} \\ \end{array} \right.,\end{eqnarray} \tag{ 3 }$$

for a given threshold ${{p}_{c}}\gt 0$.

2.2.3. Overlap of the multiplex network layers

The overlap (intersection) of the two multiplex network layers, $\phi *$ and P^*, the network O, is given by

$$\begin{eqnarray}{{O}_{i,j}}=\left\{ \begin{array}{ccccccccccccccc} 1 & {} & {\rm if}\ \phi _{ij}^{*}\gt {{\phi }_{c}}\;{\rm and}\;P_{ij}^{*}\gt {{p}_{c}} \\ 0 & {} & {\rm otherwise} \\ \end{array} \right..\end{eqnarray} \tag{ 4 }$$

The network O_ij links two diseases i and j if the link weights in both network layers exceed given thresholds—namely, ${{\phi }_{c}}$ for the statistical significance layer ${{\Phi }^{*}}$, and p_c for the conditional risk layer P^*. Note that up to first order of $\frac{1}{\sum _{j=1}^{D}{{J}_{ij}}}$ the thresholding of the statistical significance layer ${{\phi }^{*}}$ by the value ${{\phi }_{c}}$, as done in equation (4), yields the same results as the application of the multiscale filter proposed in [20]. From now on we refer to O_ij simply as the 'disease network.' It is a variant of a PHDN. In contrast to the PHDN studied in [11], the disease network O_ij is built around a combination of measures (conditional probabilities and contingency coefficients) that have been corrected for biases that result from the comparison of very rare and frequent diseases.

Let d_a be the number of diseases diagnosed in patients with a hospital stay of at least one day in a given age group during the two years 2006 and 2007. Figure 1 shows the cumulative distribution functions (cdf) for the numbers of diseases d_a for three different age groups, a = 10–20 years, a = 50–60 years, and $a\gt 80$ years, for males and females. The cdf for all age groups are presented in supplementary information (SI) figures 1 and 2. In the first age group the distribution suggests an approximate power-law behavior up to about 20 diseases. With increasing age the disease distributions takes on a more exponential character. The average number of diseases per patient in each age group, $\langle {{d}_{a}}\rangle$, is shown in the inset of figure 1. Children in the first age group have a higher average number of diagnoses than teenagers. After adolescence there is a continuous increase in the number of diseases per patient in the population, from two (at the age of around 20) to more than five for patients older than 80. Men between 30 and 70 have generally more diseases then women, which reverses at older ages. The number of diagnoses for females around 20 is higher than those for males.

Zoom In Zoom Out Reset image size

Zoom In Zoom Out Reset image size

The large increase in diagnoses per patient over age is related to substantial topological changes of the ${{\phi }^{*}}$ and P^* networks. Network properties for the two layers are shown in figure 2. The number of nodes N (figure 2(A)), the average degrees $\langle k\rangle$ (figure 2(B)), and the average nearest neighbour degrees ${{\langle k\rangle }_{nn}}$ (figure 2(D)), increase with age in both layers. For patients aged above 65 the degrees are larger for women than for men. The clustering coefficients C for each network are given in units of the clustering as expected from a random graph, C_r, in figure 2(C). $C/{{C}_{r}}$ is larger for ϕ than for the conditional risk layer for all age groups except for the youngest children. This ratio shows a clear maximum around age 25 in the ${{\phi }^{*}}$ layer, but it decreases monotonically in the P^* layer. Correlation coefficients between the degree sequences of ${{\phi }^{*}}$ and P^* are given for two choices of thresholds in the SI; see table S1. In one case we set the thresholds in the ${{\phi }^{*}}$ and P^* layers to p_c = 0.01 and ${{\phi }_{c}}=2$, respectively, and in the other case the thresholds are both zero, ${{p}_{c}}={{\phi }_{c}}=0$. We find higher degree correlations across the two layers for the non-zero threshold values, compared to the case where the thresholds are both zero. This suggests that the thresholding tends to filter out links that appear in only one of the layers. With increasing patient age, the difference between the degree correlations also increases.

Network properties for the disease network O for the threshold values ${{\phi }_{c}}=2$ and p_c = 0.01 are shown in figure 3. The numbers of nodes N (figure 3(A)), average degrees $\langle k\rangle$ (figure 3(B)), and the ratios of clustering coefficients $C/{{C}_{r}}$ (figure 3(C)), show a similar behavior with increasing age. They decrease from the first to second age group, and increase until about age 50. For higher ages the increase levels off and the network measures decrease again. The average nearest neighbour degrees do not level off at high ages but continue to increase; see figure 3(D). The disease network is most dense for both males and females for the age range 48–56. Results for the modularities of ${{\phi }^{*}}$, P^*, and O as a function of patient age and gender are shown in the SI; see figure S6. The disease network O has typically higher modularities than the layers ${{\phi }^{*}}$ and P^*, which suggests that disease clusters can be easier identified in O than in its individual layers. We confirmed that the network properties shown in figure 3 are to a large extend independent from the actual choices of the threshold values p_c and ${{\phi }_{c}}$; see SI figures S7–10.

Zoom In Zoom Out Reset image size

Visualizations of the disease networks O across lifetime are presented in figure 4 for males and females. Massive structural reorganization in the disease networks are clearly visible across age. Nodes represent diseases i; their size is proportional to the disease prevalence, $N_{i}^{m/f}(a)$. The diseases type, as represented by the first letters from the ICD 10 code, is indicated by node colors. Link colors are identical to the node colors if both diseases share the same type; otherwise, it is a mix of the two colors. Disease clusters belonging to the same ICD type (first letter of ICD code) are highlighted by colored patches in figure 4. Alternative visualizations of the disease networks O are shown in the SI, figures S21–28, for each age group and gender. In figures S21–28 we set p_c = 0.02 and adjust ${{\phi }_{c}}$ such that each network has the same average degree $\langle k\rangle$.

Zoom In Zoom Out Reset image size

We find three distinct phases in the evolution of the diseases networks across lifetime.

Phase I. The disease network for children, phase I (age 0–16), shows one cluster containing diseases of many different types, such as diseases of the respiratory system (letter J in the ICD code), infectious diseases (A), diseases of the eye and ear (H), and endocrine, nutritional, and metabolic diseases (E). Highly prevalent are viral and other specified intestinal infections (A08), diarrhea and gastroenteritis of presumed infectious origin (A09), chronic diseases of tonsils and adenoids (J35), and non-suppurative otitis media (H65). These diseases correspond to nodes that are characterized by large degrees and tend to be connected to each other. There is a small cluster formed by diseases of the nervous system (G), with a local hub being epilepsy (G40), and a cluster for mental and behavioral disorders (F). For male children there is a cluster containing disorders of psychological development that is not visible for females. For females in phase I there are small clusters containing diseases of the digestive system (K) and of the genitourinary system (N), whereas for males phimosis and paraphimosis (N47) belong to the main cluster. The network for young males is in general more dense than for females in phase I, but the formation of disease clusters and the prevalences are similar. The network for patients aged 8–16 is less dense than the network of the former age group, but the same clusters are still visible (see SI figure S12). There are obvious changes, however, for the next age group, 16–24, (see figure 4, second row), where the cluster for mental diseases (F) dominates and contains more different nodes than in the previous networks. This can be seen as the onset of the structural phase II of the disease networks.

Phase II. In phase II (age 16–64), a cluster of mental and behavioral disorders (F) appears more clearly and contains addiction to nicotine, alcohol, and other substance abuses. It has a higher density for women than for men. For both genders in the age group 16–24, the prevalences of diseases are smaller than in phase I, but infections (A08 and A09), tonsil problems (J35), and disorders of the prepuce for men (N47) are still highly prevalent. Children and young patients suffer similarly from obesity (E66), which becomes increasingly prevalent starting with patients aged 16–24y. Hypertension (I10) becomes also visible at this age, but it is not yet a dominant disease.

At ages 32–40 the cluster of diseases of the digestive system (K) emerges for men (see figure 4, third row). Hemorrhoids (I84) are found at the periphery of this cluster with links exclusively to digestive diseases; there are no links from hemorrhoids to any other circulatory disease. The cluster of diseases of the genitourinary system (N) exists still only for women and grows with age. It is connected to benign neoplasms and blood diseases (D). For patients over 40, leiomyoma of the uterus (D25) are very common. This disease has many links to both genitourinary disorders (N) and to malignant neoplasms (C), especially related to female genital organs. The most prevalent neoplasms are breast cancer (C50) for women and prostate cancer (C61) for men. Neoplasms are located at the periphery of the network and show many inter-cluster connections, such as links to secondary neoplasms.

The cluster of circulatory diseases (I) appears for men older than 32 and women older than 40. It densifies for older age groups. This cluster is also much denser for men than for women; see figure 4, fourth row. Another gender difference for this age group is the absence of genitourinary diseases (N) for men, while for women these diseases become increasingly prevalent until the age of 40. The most common circulatory disease is hypertension (I10). Already at age 40–48 it has a large number of links inside and outside the cluster of circulatory diseases.

The origin of the cluster related to diseases of the musculoskeletal system (M) is visible for both genders at age 24–32 (see SI figure S14); for older groups, this cluster becomes denser more quickly for females than for males. This may be related to osteoporosis (M81), which is a hub only for women over 48 and serves as a gateway to many musculoskeletal diseases. An important structural change begins at age 64, where some clusters start to shrink, such as mental disorders (F) and genitourinary diseases (N) for women. This marks the beginning of phase III (age 64–80) of the evolution of the disease network.

Phase III. A contraction of the network toward its center can be observed in figure 4, bottom row, where it becomes increasingly dominated by circulatory (I) and metabolic (E) diseases. Global hubs appear—in particular, hypertension (I10), depression (F32), ischemic heart diseases (I25), chronic obstructive pulmonary disease (COPD, J44), the metabolic syndrome including disorders of lipoprotein metabolism (E78), type 2 diabetes (E11), and obesity (E66). There are also some local hubs—for example, cataracts (H25). Clearly discernible clusters disappear in phase III, with the exception of musculoskeletal disorders (M) for females. Neoplasms (C) become increasingly interconnected. The effect of the network contraction progresses in the oldest groups (age 72–80), where it is even more difficult to distinguish specific clusters (see SI figure S19).

The three phases are also apparent in the evolution of the network properties of the disease network; see figure 3. In phase I (childhood) we find large clustering and degrees at a comparably small network size N. This suggests a tightly interconnected cluster of diseases spanning only a small part of the entire possible network. Phase II presents itself as a steady increase in clustering and network density. This trend levels off once substantial parts of the network show high levels of disease prevalences. In phase III clustering and density decrease, whereas the average nearest neighbor degree increases. This suggests the emergence of hubs that become connected to an increasing number of other diseases, whereas the local clusters that characterize phase II disappear. On a more quantitative level, we fitted the degree distributions with q-exponentials [27, 28] and report the values of q for males and females across age groups in figures S2–4 in the SI. In general, q increases over lifetime, but there is a prominent peak at ages 16 and 24, indicating again a massive change of regime in network topology.