# Effect of Genome and Environment on Metabolic and Inflammatory Profiles

Marina Sirota, Gonneke Willemsen, Purnima Sundar, Steven Pitts, Shobha Potluri, Edi Prifti, Sean Kennedy, S. Dusko Ehrlich, Jacoline Neuteboom, Cornelis Kluft, et al.

► To cite this version:

Marina Sirota, Gonneke Willemsen, Purnima Sundar, Steven Pitts, Shobha Potluri, et al. Effect of Genome and Environment on Metabolic and Inflammatory Profiles. PLoS ONE, 2015, 10 (4), pp.e0120898. 10.1371/journal.pone.0120898. pasteur-02552107

HAL Id: pasteur-02552107

https://pasteur.hal.science/pasteur-02552107v1

Submitted on 23 Apr 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
RESEARCH ARTICLE

# Effect of Genome and Environment on Metabolic and Inflammatory Profiles

Marina Sirota¹*, Gonneke Willemsen²*, Purnima Sundar¹, Steven J. Pitts¹, Shobha Potluri¹, Edi Prifti³, Sean Kennedy³, S. Dusko Ehrlich³, Jacoline Neuteboom⁴, Cornelis Kluft⁴, Karen E. Malone⁴, David R. Cox¹, Eco J. C. de Geus², Dorret I. Boomsma²

1. Rinat-Pfizer, South San Francisco, California, United States of America,
2. Department of Biological Psychology, VU University Amsterdam, Amsterdam, Netherlands,
3. INRA, Institut National de la Recherche Agronomique, Metagenopolis, Jouy-en-Josas, France,
4. Good Biomarker Sciences, Leiden, Netherlands

© These authors contributed equally to this work.

* a.h.m.willemsen@vu.nl

<figure><img src="image_2.png" /></figure>

CrossMark

click for updates

<figure><img src="image_3.png" /></figure>

## OPEN ACCESS

**Citation:** Sirota M, Willemsen G, Sundar P, Pitts SJ, Potluri S, Prifti E, et al. (2015) Effect of Genome and Environment on Metabolic and Inflammatory Profiles. PLoS ONE 10(4): e0120898. doi:10.1371/journal.pone.0120898

**Academic Editor:** Suryaprakash Sambhara, Centers for Disease Control and Prevention, UNITED STATES

**Received:** August 7, 2014

**Accepted:** January 27, 2015

**Published:** April 8, 2015

**Copyright:** © 2015 Sirota et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** The data were collected as part of the NTR Biobank, they can be released by application from NTR. Data are available from the VU Institutional Data Access / Ethics Committee for researchers who meet the criteria for access to confidential data. The data can be requested from ntr@vu.nl.

**Funding:** Funding was obtained from INRA-Pfizer, the Netherlands Organisation for Scientific Research (NWO 911-09-032; NWO 480-04-004), EMGO+ Institute for Health and Care Research, Neuroscience Campus Amsterdam, BBMRI-NL (Biobanking and

## Abstract

Twin and family studies have established the contribution of genetic factors to variation in metabolic, hematologic and immunological parameters. The majority of these studies analyzed single or combined traits into pre-defined syndromes. In the present study, we explore an alternative multivariate approach in which a broad range of metabolic, hematologic, and immunological traits are analyzed simultaneously to determine the resemblance of monozygotic (MZ) twin pairs, twin-spouse pairs and unrelated, non-cohabiting individuals. A total of 517 participants from the Netherlands Twin Register, including 210 MZ twin pairs and 64 twin-spouse pairs, took part in the study. Data were collected on body composition, blood pressure, heart rate, and multiple biomarkers assessed in fasting blood samples, including lipid levels, glucose, insulin, liver enzymes, hematological measurements and cytokine levels. For all 51 measured traits, pair-wise Pearson correlations, correcting for family relatedness, were calculated across all the individuals in the cohort. Hierarchical clustering techniques were applied to group the measured traits into sub-clusters based on similarity. Sub-clusters were observed among metabolic traits and among inflammatory markers. We defined a phenotypic profile as the collection of all the traits measured for a given individual. Average within-pair similarity of phenotypic profiles was determined for the groups of MZ twin pairs, spouse pairs and pairs of unrelated individuals. The average similarity across the full phenotypic profile was higher for MZ twin pairs than for spouse pairs, and lowest for pairs of unrelated individuals. Cohabiting MZ twins were more similar in their phenotypic profile compared to MZ twins who no longer lived together. The correspondence in the phenotypic profile is therefore determined to a large degree by familial, mostly genetic, factors, while household factors contribute to a lesser degree to profile similarity.
Biomolecular Resources Research Infrastructure), and the European Research Council (Genetics of Mental Illness 230374). Pfizer provided support in the form of salaries for authors PS, SP and SJP, but did not have any additional role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript. The specific roles of these authors are articulated in the author contributions section.

**Competing Interests:** PS, SP and SJP are employees of Pfizer, whose company provide funding towards this study. There are no patents, products in development or marketed products to declare. This does not alter the authors' adherence to all the PLOS ONE policies on sharing data and materials.

# INTRODUCTION

Twin and family studies have been historically used to examine the effects of genetics and environment on a wide range of complex human phenotypic traits [1-4]. There has been an abundance of studies examining the influence of genetic factors on a large spectrum of cardiovascular, immunological and metabolic traits. For most physiological parameters significant heritability is demonstrated. For instance, heritability of blood pressure ranges from 30 to 60% [5,6], while the heritability estimate for BMI can be even higher, above 70% [7-9], though it may differ as a function of age and social situation [10]. Previous studies have also shown genetic influences on indices of pro-inflammatory state, with heritability estimates ranging between 20 and 45% [11,12]. These measurements, however, are not independent from one-another. For instance, the metabolic syndrome has been characterized as consisting of several well established risk factors for cardiovascular disease, whose co-variation is largely due to genetic factors [13]. A syndrome score reflecting all risk factors is often used, which can result in a loss of information. Moreover, the metabolic syndrome (obesity, hyperlipidemia, hyperglycemia and hypertension) is associated with many other traits that may also explain cardiovascular disease, but are not formally part of the syndrome. In general, many studies that assess a large range of correlated biomarkers tend to analyze them in a one-by-one fashion, which fails to take into account their often-substantial inter-dependence.

Over the last decade, increasing attention is given to the methodology for the joint analyses of multiple traits, including principal components, factor analyses and other structural equation modeling approaches and recently the application of item response theory [14-17]. However, some of these techniques are difficult to implement in studies of familial resemblance, when for example the number of traits is large relative to the number of families. Here, we propose an approach, in which we define a phenotypic profile of an individual as the collection of **all** quantitatively measured traits for that individual and apply correlation metrics to capture the pairwise similarity across a large number of quantitative traits. By comparing the pairwise similarity of the phenotypic profiles within specific groups, which share genes and environment to a different extent, the underlying cause of variation in the phenotypic profile can be explored. To illustrate this, we applied our approach (see Methods) to a range of metabolic, hematological, and immunological traits (Table 1), which were collected in a large cohort of monozygotic (MZ) twins and their spouses from the Netherlands Twin Register. We explored positive and negative linear relationships between the traits correcting for age, sex and the family structure (see Methods) and used a hierarchical clustering approach to further test and visualize the relationships. These relationships can be used in the future to inform study design and account for the necessary covariates. They also provide insight into which traits are most informative to measure.

In addition to studying relationships between the traits, we examine the similarity of specific phenotypic traits as well as phenotypic profiles across individuals. The inclusion of spouses, who share a household, in addition to MZ twin pairs, who share their genome, allow for an exploration of genetic and environmental influences on the phenotypic profile by determining the similarity in profile for MZ pairs, twin-spouse pairs and unrelated individuals (who share neither genes nor environment). Family members may resemble one-another due to shared genes and shared environments. In the classical twin design the similarity in MZ and DZ twins is compared to disentangle the effect of these genetic and the shared environmental influences. Here, we use an alternative method by contrasting MZ twin similarity to the similarity in spouse pairs, who share a household but can be considered not to share a genetic background [18]. A higher profile similarity in MZ twins compared to unrelated individuals indicates familial, genetic and/or shared—childhood- environment influences (Fig. 1A). A staircase pattern
<table><caption>Table 1. List of measured traits with their median, mean and standard deviation.</caption><thead><tr><td>Variable</td><td>Description</td><td>Median</td><td>Mean</td><td>Standard Deviation</td></tr></thead><tbody><tr><td>Age</td><td>Age</td><td>34.0</td><td>33.88</td><td>11.58</td></tr><tr><td>ALT</td><td>Alanine transaminase (U/L)</td><td>18.0</td><td>19.14</td><td>7.48</td></tr><tr><td>AST</td><td>Aspartate aminotransferase (U/L)</td><td>23.0</td><td>24.90</td><td>8.30</td></tr><tr><td>basophils</td><td>Basophils (10*9/L)</td><td>0.0</td><td>0.02</td><td>0.01</td></tr><tr><td>BMI</td><td>Body mass index (kg/m2)</td><td>23.6</td><td>24.25</td><td>3.79</td></tr><tr><td>cholesterol</td><td>Total cholesterol (mmol/L)</td><td>4.8</td><td>4.90</td><td>0.91</td></tr><tr><td>cotinine</td><td>Cotinine</td><td>0.0</td><td>83.51</td><td>254.53</td></tr><tr><td>creatinine</td><td>Creatinine (mol/L)</td><td>80.0</td><td>80.13</td><td>12.17</td></tr><tr><td>CRP</td><td>C-reactive protein (mg/L)</td><td>1.2</td><td>2.25</td><td>2.81</td></tr><tr><td>DBP</td><td>Diastolic blood pressure (mmHg)</td><td>80.0</td><td>80.39</td><td>10.79</td></tr><tr><td>eosinophils</td><td>Eosinophils (10*9/L)</td><td>0.2</td><td>0.19</td><td>0.13</td></tr><tr><td>fibrinogen</td><td>Fibrinogen (g/L)</td><td>3.1</td><td>3.20</td><td>0.63</td></tr><tr><td>GGT</td><td>Gamma glutamyl transpeptidase (U/L)</td><td>19.0</td><td>23.96</td><td>15.98</td></tr><tr><td>glucose</td><td>Glucose (mmol/L)</td><td>4.7</td><td>4.76</td><td>0.47</td></tr><tr><td>GMCSF</td><td>Granulocyte-macrophage colony-stimulating factor (pg/ml)</td><td>0.0</td><td>0.34</td><td>0.82</td></tr><tr><td>HCT</td><td>Hematocrit (ratio)</td><td>0.4</td><td>0.42</td><td>0.03</td></tr><tr><td>HDL</td><td>HDL cholesterol (mmol/L)</td><td>1.4</td><td>1.44</td><td>0.37</td></tr><tr><td>height</td><td>Height (cm)</td><td>173.0</td><td>173.40</td><td>9.37</td></tr><tr><td>HGB</td><td>Hemoglobin (mmol/L)</td><td>8.7</td><td>8.80</td><td>0.77</td></tr><tr><td>hip</td><td>Hip circumference (cm)</td><td>102.0</td><td>102.72</td><td>7.74</td></tr><tr><td>HR</td><td>Heart rate (bpm)</td><td>71.0</td><td>71.82</td><td>11.81</td></tr><tr><td>IFNg</td><td>Interferon gamma (pg/ml)</td><td>1.9</td><td>2.37</td><td>2.22</td></tr><tr><td>IL10</td><td>Interleukin-10 (pg/ml)</td><td>1.7</td><td>4.53</td><td>10.89</td></tr><tr><td>IL12p70</td><td>Interleukin-12p70 (pg/ml)</td><td>2.3</td><td>10.14</td><td>38.60</td></tr><tr><td>IL1b</td><td>Interleukin-1 beta (pg/ml)</td><td>0.0</td><td>0.44</td><td>0.74</td></tr><tr><td>IL2</td><td>Interleukin-2 (pg/ml)</td><td>0.0</td><td>0.22</td><td>0.54</td></tr><tr><td>IL6</td><td>Interleukin-6 (pg/ml)</td><td>0.8</td><td>0.97</td><td>1.34</td></tr><tr><td>IL6R</td><td>Interleukin-6 receptor</td><td>380.5</td><td>392.57</td><td>117.44</td></tr><tr><td>IL8</td><td>Interleukin-8 (pg/ml)</td><td>2.7</td><td>2.64</td><td>1.21</td></tr><tr><td>insulin</td><td>Insulin (microlU/ml)</td><td>5.0</td><td>6.10</td><td>5.32</td></tr><tr><td>LDL</td><td>LDL cholesterol (mmol/L)</td><td>2.8</td><td>2.85</td><td>0.82</td></tr><tr><td>lymphocytes</td><td>Lymphocytes (10*9/L)</td><td>2.2</td><td>2.25</td><td>0.69</td></tr><tr><td>MCH</td><td>Mean corpuscular hemoglobin (fmol)</td><td>1.8</td><td>1.84</td><td>0.10</td></tr><tr><td>MCHC</td><td>Mean corpuscular hemoglobin concentration (mmol/L)</td><td>20.8</td><td>20.82</td><td>0.50</td></tr><tr><td>MCV</td><td>Mean corpuscular volume (fl)</td><td>88.3</td><td>88.32</td><td>4.15</td></tr><tr><td>monocytes</td><td>Monocytes (10*9/L)</td><td>0.5</td><td>0.55</td><td>0.17</td></tr><tr><td>MPV</td><td>Mean platelet volume (fl)</td><td>11.1</td><td>11.12</td><td>0.91</td></tr><tr><td>neutrophils</td><td>Neutrophils (10*9/L)</td><td>3.0</td><td>3.19</td><td>1.13</td></tr><tr><td>PCT</td><td>Platelet Count</td><td>0.3</td><td>0.28</td><td>0.06</td></tr><tr><td>PDW</td><td>Platelet distribution width</td><td>13.3</td><td>13.59</td><td>2.06</td></tr><tr><td>PLT</td><td>Platelets (10*9/L)</td><td>244.0</td><td>251.53</td><td>55.72</td></tr><tr><td>RBC</td><td>Red blood cell count (10*12/L)</td><td>4.8</td><td>4.79</td><td>0.42</td></tr><tr><td>RDWCV</td><td>Red cell distribution width, coefficient variation</td><td>12.9</td><td>13.00</td><td>0.83</td></tr><tr><td>RDWSD</td><td>Red cell distribution width, standard variation</td><td>40.9</td><td>41.06</td><td>2.38</td></tr><tr><td>SBP</td><td>Systolic blood pressure (mmHg)</td><td>127.0</td><td>128.87</td><td>15.15</td></tr><tr><td>telomere length</td><td>Telomere length</td><td>2.6</td><td>2.66</td><td>0.40</td></tr><tr><td>TNFA</td><td>Tumor necrosis factor alpha (pg/ml)</td><td>6.3</td><td>6.66</td><td>1.96</td></tr></tbody></table>

(Continued)
<table><caption>**Table 1. (Continued)**</caption><thead><tr><td>Variable</td><td>Description</td><td>Median</td><td>Mean</td><td>Standard Deviation</td></tr></thead><tbody><tr><td>triglycerides</td><td>Triglycerides (mmol/L)</td><td>1.1</td><td>1.34</td><td>0.78</td></tr><tr><td>waist</td><td>Waist circumference (cm)</td><td>78.0</td><td>80.02</td><td>11.49</td></tr><tr><td>WBC</td><td>White blood cell count (10*9/L)</td><td>6.0</td><td>6.18</td><td>1.54</td></tr><tr><td>weight</td><td>Weight (kg)</td><td>71.4</td><td>73.17</td><td>14.42</td></tr><tr><td>WHR</td><td>Waist to hip ratio (cm/cm)</td><td>0.8</td><td>0.78</td><td>0.08</td></tr></tbody></table>

doi:10.1371/journal.pone.0120898.t001

with similarity being highest in MZ twins, lower in spouses, and much lower among unrelated individuals, would indicate influence of shared environment, where spouse resemblance provides an upper limit for household effects if phenotypic assortative mating is also present. In

<figure><img src="image_4.png" /><figcaption>**Fig 1.** Study overview and participant characterization. A) Study Design. Pairs of twins were recruited for this study. A proportion of the cohort consists of twins who are co-habiting. Spouses were also recruited to index resemblance through shared household factors; B) Age distribution in twin study participants stratified by co-habitation status; C) Age distribution in study participants stratified by sex.</figcaption></figure>

doi:10.1371/journal.pone.0120898.g001
contrast, a pattern in which the similarity for spouses is comparable to unrelated individuals, with MZ twin similarity being higher, would indicate that household effects do not play a role. The importance of sharing a household is further explored by determining profile similarity in MZ twins who cohabit and in MZ twins who live separate from each other. If non-cohabiting twins are less alike in their profile than cohabiting twins, this points to an influence of shared household factors, although we recognize that cohabiting twins tend to be younger.

# MATERIALS AND METHODS

## Participants

The participants in this study were recruited from the Adult Netherlands Twin Register (ANTR), a large twin-family based register for the study of individual differences in health, lifestyle and personality across the lifespan [19,20]. Adult participants were invited in the present study when they did not report any diseases related to the immune system, endocrine system or gastrointestinal tract in earlier survey studies. The study mainly targeted MZ twin pairs and their spouses, but in the initial pilot phase of the study other healthy family members were also invited. Invitation letters to NTR participants were followed by a phone call. Of those participants reached by phone, 71% agreed to take part in the study. The final cohort consisted of 517 individuals (67% female): 447 twins (210 pairs and 27 unpaired twins), 64 spouses, 1 sibling, 2 fathers, 2 mothers and 1 singleton. The age of the participants ranges from 19 to 68 (mean 33.9, standard deviation 11.9). The age bimodal distribution stratified on sex (more females than males for the younger individuals) and cohabitation status are shown in Fig. 1C and 1B. All data are used for the analysis of relationships among traits. The heritability analysis focuses on the 210 MZ twin pairs, and 64 twin-spouse pairs included within this cohort. Within the 210 pairs of twins, 57 were cohabiting at the time of the sampling and 153 pairs were not. The zygosity was confirmed by DNA markers. Ethical approval was obtained from the Central Ethics Committee on Research Involving Human Subjects of the VU University Medical Center, Amsterdam. Written informed consent was obtained from all participants. The study was approved by the Central Ethics Committee on Research Involving Human Subjects of the VU University Medical Centre, Amsterdam, an Institutional Review Board certified by the U.S. Office of Human Research Protections (IRB number IRB-2991 under Federal-wide Assurance-3703; IRB/institute codes, NTR 03-180).

## Data collection

Participants were visited between 7 and 10 AM at home or, in a few cases, at work. They were instructed to be fasting and avoid exercise as of the night before and to refrain from smoking and, when possible, medication prior to the visit. Fertile women were bled on a fixed day (day 3) of their menstrual cycle, or during the placebo week of their birth control hormone regimen, if feasible. During the visit, participants were asked to provide information on fasting status, general health, medication use and lifestyle (such as a smoking status, diet).

**Body composition indices.** Height, weight, waist and hip circumference were determined, which were used to calculate body mass index (BMI) and waist-to-hip ratio (WHR).

**Cardiovascular measurements.** Using an OMRON HEM907 blood pressure device, systolic blood pressure (SBP), diastolic blood pressure (DBP) and heart rate (HR) were measured twice, before and after the blood sampling. In the present study we used the average SBP, DBP and HR across the two measurements.

**Blood Sample Collection.** Using the Vacutainer system, 8 blood tubes were collected: 2 x 9 ml EDTA, 1 x 9 ml Li heparine, 1 x 9 ml Na heparine, 1 x 4.5 ml CTAD, 1 x 2.5 ml PAX, 1 x 5 ml
serum and 1 x 2 ml EDTA. To prevent clotting, all tubes were inverted gently 8–10 times immediately after collection [21].

**Lipids.** Total cholesterol, HDL-cholesterol and triglycerides were determined in plasma from a lithium heparin tube, which was stored in melting ice during transport and processed at the laboratory within 6 hours of transport. After centrifugation of the tube for 15 minutes at 1000x g at 4° C, heparin plasma was obtained and divided into 8 subsamples of 0.5 ml, snap-frozen and stored at—30° C. One these subsamples provided lipid levels, using Vitros 250 total cholesterol, Vitros 250 direct HDL cholesterol and Vitros 250 Triglycerides assays (Johnson & Johnson, Rochester, USA) respectively. LDL-cholesterol was calculated using the Friedewald Equation: $LDL = Total cholesterol—HDL—(Triglycerides/5)$ [22].

**Glucose metabolism.** Insulin was measured using the Immulite 1000 Insulin Method (Diagnostic Product Corporation, Los Angeles, USA) and glucose was measured using the Vitros 250 Glucose assay (Johnson & Johnson, Rochester, USA).

*Liver enzymes.* GGT, ALT and AST levels in units per liter (U/L) were determined for one of the lithium heparin plasma aliquots, using Vitros assays (Vitros 250, Ortho-Phenotypic Diagnostics; Johnson & Johnson, Rochester). Please see Table 1 for a complete list.

**Hematology.** The 2 ml EDTA tube was stored at room temperature during transport and used to obtain a hematological profile via the Coulter system (Coulter Corporation, Miami, USA). Parameters provided were white blood cell count (WBC), numbers of neutrophils, lymphocytes, monocytes, eosinophils, basophils as well as red blood cell count, hemoglobin, hematocrit, mean corpuscular volume (MCV), mean corpuscular hemoglobin (MCH), mean corpuscular hemoglobin concentration (MCHC), red cell distribution width (RDW), platelet count (PCT) and mean platelet volume (MPV). Please see Table 1 for a complete list.

**Cytokines.** Nine cytokines were determined in blood from one of the EDTA tubes, which were stored in melting ice during transport. Upon arrival at the laboratory, the EDTA tubes were centrifuged for 20 minutes at 2000x g at 4° C. EDTA plasma was harvested and aliquoted (0.5 ml), snapfrozen in dry ice, and stored at—30° C. Using the Meso Scale Discovery platform and one of the EDTA plasma aliquots, the levels of IL6, TNF-α, GMCSF, IFN-γ, IL10, IL12p70, IL1b, IL2, and IL8 were determined. Please see Table 1 for a complete list.

**C-reactive protein (CRP).** CRP was obtained from one of the lithium heparin aliquots using the Immulite 1000 CRP assay (Diagnostic Product Corporation, USA).

*Fibrinogen.* Fibrinogen levels were obtained from CTAD plasma. The 4.5 ml CTAD tube was stored during transport in melting ice and upon arrival at the laboratory, centrifuged for 20 minutes at 2000x g at 4° C, after which citrated plasma was harvested from the buffy coat and red blood cells, aliquoted (0.5 ml), snap frozen in dry ice, and stored at—30° C. One of these aliquots was used to obtain fibrinogen levels on a STA Compact Analyzer (Diagnostica Stago, France), using STA Fibrinogen (Diagnostica Stago, France).

**Data Analysis.** Data were inspected and extreme outliers (defined below) were removed. For individuals who were taking antihypertensive medication (N = 24), systolic blood pressure was increased by 14 mmHg and diastolic blood pressure by 10 mmHg [23]. Triglyceride levels were set to missing when non-fasting (N = 8), and 0.74 mmol/l was added to LDL cholesterol and total cholesterol when LDL-lowering medication (statins, N = 8) were used [24]. Glucose and insulin levels were set to missing when non-fasting (N = 8) and additionally when on diabetes medication (N = 4). In the case of glucocorticoid or anti-inflammatory medication, cytokine levels were set to missing (N = 8).

Based on inspection of the distribution of the hematology data we set data to missing when white blood cell count was above 12 (N = 7), neutrophil count above 8 (N = 5), lymphocyte count above 5 (N = 2), monocyte count above 1.1 (N = 6) and basophil count above 10 (N = 1). We further set to missing MCHC when above 24 (N = 1), HCT above. 6 (N = 1), HGB above
14 (N = 1) and red blood cell count above 7 (N = 1). Values were set to missing when individuals were using anti-inflammatory or endocrine medication or had values above 5 for CRP (N = 21), IL6 (N = 9), TNFalpha (N = 22), IFNgamma (N = 17), or above 6 for fibrinogen (N = 7), above 110 for IL10 (N = 24), above 500 for IL12p70 (N = 16), and above 5 for IL1b (N = 11), IL2 (N = 10), and IL8 (N = 14).

Table 1 provides an overview of the 51 quantitatively measured traits included in the analyses. It's important to note that several cytokines have skewed distributions. The phenomenon of such tailed distributions might be due to transient or chronic immunological or inflammatory processes. Many of the observed values are below detection limit (29% for Il-6; 77% for GMCSF, 19% for IFN-gamma, 65% for Il-1-beta and 83% for Il-2) resulting in a difference between means and median measurements, which are also reported in Table 1. We recognize that this is a limitation of the measurement modalities that are applied here. In order to correct for age and sex variability in the cohort, we computed linear regression residuals for each trait. The distributions of these residual values per trait are shown in S1 Fig.

First, we tested the correlation between all traits in the full sample of 517 individuals by computing partial pairwise Pearson correlations across all the individuals in the cohort (ppcor Package in R), in which we controlled for the relatedness between individuals by taking into account individuals from the same family and their relationship to each other (twin or spouse). In total, 1275 pairwise trait correlations were generated. Based on Bonferroni correction for multiple hypothesis testing, the threshold for significance was < 3.92e-05. We applied a hierarchical clustering technique using the complete linkage method [25] to group the measured traits into sub-clusters based on 1—correlation as a distance metric. Complete linkage method is a type of agglomerative clustering where each phenotypic trait starts in its own cluster, and consequently pairs of clusters are merged as one moves up the hierarchy based on the 1—partial Pearson correlation value as a distance metric.

To establish the relative contribution of genetic and household factors to the variance in all individual traits, we computed the pair-wise correlation between MZ twin pairs (living together and living apart), twin-spouse pairs as well as randomly sampled unrelated individuals as a measure of heritability for each trait. The unrelated individual pairs were constructed by randomly sampling from the cohort for each person excluding their twin and spouse from the cohort. For each of the 51 measured traits we computed: **twin_correlation_cohab** = correlation(t[twin1], t[twin2]) where twin1 and twin2 are cohabiting, **twin_correlation_not_cohab** = correlation(t[twin1], t[twin2]), where twin1 and twin2 are not cohabiting, **twin_spouse_correlation** = correlation(t[twin], t[spouse]) and **unrelated_individuale** = correlation(t[indiv1], t[indiv2]) where indiv1 and indiv2 are chosen at random.

The phenotypic profile p of an individual was defined as a vector of all the phenotypic measurements (n = 51) for that individual. Pearson correlation, capturing the resemblance between pairs of vectors, is used a measure of similarity of two phenotypic profiles. Pearson correlation coefficient is a measure of the linear correlation between two variables X and Y, giving a value between +1 and -1 inclusive, where 1 is total positive correlation, 0 is no correlation, and -1 is total negative correlation. It is furthermore defined as the covariance of the two variables divided by the product of their standard deviations. Pairwise profile similarity between all pairs of individuals was determined for four subgroups: MZ twin pairs sharing (57) or not-sharing a household (153), twin-spouse pairs (64) and pairs of unrelated individuals (132,596). Here for each within pair comparison in our four groups we compute a similarity score s = correlation (p_indiv1, p_indiv2).

All statistical analyses were performed in R.
# RESULTS

## Characterization and Clustering of the Phenotypic Traits

Table 1 includes all the measured traits used for the analysis with their median, mean and standard deviation in the dataset. Using the data from the 517 individuals in our sample, we computed partial Pearson correlations, correcting for the effect of family relatedness, for the 51 measured traits. Out of the 1,275 relationships tested, 169 were significant at $p < 0.05/1275 \sim 4e05$ (correlations ranged between -0.558 to +0.996, see S1 Table).

We applied hierarchical clustering to elucidate and visualize relationships between the measured traits (Fig. 2). Many of the body composition measurements such as weight, waist, hip, waist to hip ratio were strongly correlated with blood pressure, cholesterol levels, insulin and glucose measurements and inversely correlated with HDL levels, defining the metabolic syndrome. As to be expected, hematocrit and hemoglobin measurements were highly correlated (correlation coefficient 0.89–0.93, $p < 4e05$) with other types of red blood cell count measurements and inversely related to measures of red blood cell volume and size. Inflammatory markers such as white blood cell count (WBC), neutrophil, lymphocyte, eosinophil, monocyte and basophil counts were correlated with IL6 cytokine levels (correlation coefficients 0.24–0.28, $p < 4e05$). Those traits also correlated positively with platelet counts (PCT, PLT, correlation coefficient 0.27–0.35, $p < 4e05$). A positive correlation was observed between IFN-γ and TNF-α (correlation coefficient = 0.25, $p < 4e05$), which are known to act synergistically during inflammation [26]. IFN-γ and TNF-α also cluster together with IL2, IL1b and GMCSF, which are known to have common mechanisms in response to co-stimulatory signals in T-cells [27]. Levels of IL10 and IL12, two cytokines which are both involved in Th1 T-cell differentiation [28], were highly correlated (correlation coefficient = 0.95, $p < 4E-05$).

Interestingly, we observed a weak, but significant positive correlation between cotinin, an alkaloid found in tobacco and a metabolite of nicotine, and various inflammatory markers (correlation coefficient 0.23–0.32, $p < 4e05$). We also observed a weak positive correlation between heart rate and glucose as well as insulin levels (correlation coefficient = 0.19, $p < 4e05$). Other studies have demonstrated that an increased heart rate at rest is associated with glucose levels [29,30].

A strong correlation was seen between fibrinogen, a coagulation marker, and CRP, a general measure of inflammation (correlation coefficient = 0.46, $p < 4e05$). The liver synthesizes both CRP and fibrinogen and in this study they were both positively correlated with body composition (weight, waist and BMI, correlation coefficient 0.24–0.3, $p < 4e05$). Fibrinogen was also positively correlated with other inflammatory markers (WBC and neutrophil counts, IL6 levels, correlation coefficients 0.19–0.26, $p < 4e05$).

Weak positive correlations were also seen for white blood cell count, an inflammatory marker, with insulin and triglyceride levels (correlation coefficients 0.27–0.29, $p < 4E-05$) and negatively correlated with HDL (correlation coefficient -0.22, $p < 4e05$).

## Trait similarity across twins, spouses and unrelated individuals

In order to establish the relative importance of genetic and household factors for each trait, we examined the pairwise similarity for each of the separate traits in four groups—unrelated individuals, spouse pairs, non-cohabiting MZ twins and cohabiting MZ twins respectively. Fig. 3 shows the distribution of the correlation coefficients for the 51 traits in MZ twin pairs living together (mean = 0.56), MZ twin pairs living apart (mean = 0.48), spouse pairs (mean = 0.08) and pairs of unrelated individuals (mean = 0.00). Although we do see a higher average correlation across the traits in MZ twins who are living together than for those living apart, the
<figure><img src="image_5.png" /><figcaption>**Fig 2.** Hierarchical clustering of phenotypic measurements. This figure shows the clustering of the phenotypic measurements based on their correlations with each other. Correlations between all pairs of phenotypes are computed, while partiaLLing out the effect related to family membership. Positive correlations are shown in red, negative correlations in purple. The clusters are shown in green.</figcaption></figure>

doi:10.1371/journal.pone.0120898.g002

difference is not statistically significant. The traits are much more correlated in both twin
groups than in twin-spouse pairs or unrelated individuals (all comparisons were highly signifi-
cant, p < 2.2e16). In comparison to unrelated individuals (mean = 0.00) the distribution of cor-
relations between twin-spouse (mean = 0.08) was significantly higher (p-value = 0.0001)
demonstrating the presence of a small environmental effect.
Phenotype Correlations (51 Traits)

<figure><img src="image_6.png" /><figcaption>**Fig 3.** Trait similarity across twins, spouses and unrelated individuals. This figure shows the boxplot with distribution across the 51 measured traits showing correlation in twins (cohabiting and not co-habiting), twin-spouse pairs and unrelated individuals. Average within pair correlation for the MZ twins living together is 0.56 (living together), for MZ twins living apart 0.48 (living apart), for twin-spouse pairs 0.08 and for non-related individuals 0.</figcaption></figure>

doi:10.1371/journal.pone.0120898.g003

MZ pair resemblance for certain traits is more striking than for others (Fig. 4A, Table 2). As expected [31], highest resemblance was found for height (correlation coefficient = 0.9, p-value < 2.2e16). In addition, telomere length, which has been previously associated with normal ageing[32,33], was highly correlated (correlation coefficient = 0.875, p-value < 2.2e16) in MZ twins (Fig. 4B). This is in line with a previous study by Broer et al. [34], who also demonstrated high heritability for telomere length. We further observed a high MZ twin correlation for the plasma level of the soluble IL6 receptor (correlation coefficient = 0.8, p-value < 2.2e16),
<figure><img src="image_7.png" /><figcaption>**Fig 4.** Examples of traits and their concordance in twins. MZ twin pair concordance for A) height, B) telomere length, C) IL6 receptor, D) IL6. In black the correlations for twins living apart, in red the correlations for twins living together. Values shown are sex and age corrected residuals (see Methods).</figcaption></figure>

doi:10.1371/journal.pone.0120898.g004

whereas the resemblance for plasma levels of the cytokine itself was very low (correlation coefficient = 0.09, p-value = ns) Fig. 4C and 4D.

Evidence for a moderate influence of shared genetic background is seen for other inflammatory markers, platelet measurements, HDL and LDL cholesterol, as well as total cholesterol and creatinine levels. We observed that blood pressure was more correlated in the MZ twins (correlation coefficient > 0.5, p-value < 2.2e16) than in spouses (correlation coefficient = 0.21,
<table><caption>Table 2. Pair-wise correlations (confidence intervals are provided in S2 Table) for each trait in each of the four comparison groups (MZ twin pairs living together, MZ twin pairs living apart, twin-spouse pairs, unrelated individuals) ordered by concordance in twins living together.</caption><thead><tr><th rowspan="2">Variable</th><th rowspan="2">Description</th><th colspan="2">Twin-Twin Living Together</th><th colspan="2">Twin-Twin Living Apart</th><th colspan="2">Twin-Spouse</th><th colspan="2">Unrelated</th></tr><tr><th>Correlation</th><th>P-Value</th><th>Correlation</th><th>P-Value</th><th>Correlation</th><th>P-Value</th><th>Correlation</th><th>P-Value</th></tr></thead><tbody><tr><td>telomere length</td><td>Telomere length</td><td>0.9133</td><td>0.0000*</td><td>0.8561</td><td>0.0000*</td><td>-0.0279</td><td>0.8270</td><td>0.0003</td><td>0.9963</td></tr><tr><td>cotinine</td><td>Cotinine</td><td>0.9001</td><td>0.0000*</td><td>0.3469</td><td>0.0000*</td><td>0.2950</td><td>0.0189</td><td>0.0739</td><td>0.2897</td></tr><tr><td>height</td><td>Height (cm)</td><td>0.8841</td><td>0.0000*</td><td>0.9126</td><td>0.0000*</td><td>0.1451</td><td>0.2526</td><td>0.0440</td><td>0.5265</td></tr><tr><td>triglycerides</td><td>Triglycerides (mmol/L)</td><td>0.8838</td><td>0.0000*</td><td>0.3777</td><td>0.0000*</td><td>0.0783</td><td>0.5454</td><td>-0.0491</td><td>0.4885</td></tr><tr><td>hip</td><td>Hip circumference (cm)</td><td>0.8035</td><td>0.0000*</td><td>0.6716</td><td>0.0000*</td><td>0.3027</td><td>0.0159</td><td>-0.0613</td><td>0.3812</td></tr><tr><td>BMI</td><td>Body mass index (kg/m2)</td><td>0.8007</td><td>0.0000*</td><td>0.6239</td><td>0.0000*</td><td>0.2560</td><td>0.0464</td><td>-0.0252</td><td>0.7193</td></tr><tr><td>weight</td><td>Weight (kg)</td><td>0.7901</td><td>0.0000*</td><td>0.6832</td><td>0.0000*</td><td>0.2053</td><td>0.1125</td><td>-0.0968</td><td>0.1661</td></tr><tr><td>IL6R</td><td>Interleukin-6 receptor</td><td>0.7794</td><td>0.0000*</td><td>0.8038</td><td>0.0000*</td><td>0.2900</td><td>0.0259</td><td>-0.1104</td><td>0.1189</td></tr><tr><td>waist</td><td>Waist circumference (cm)</td><td>0.7671</td><td>0.0000*</td><td>0.5477</td><td>0.0000*</td><td>0.1174</td><td>0.3595</td><td>-0.0626</td><td>0.3705</td></tr><tr><td>PCT</td><td>Platelet Count</td><td>0.7504</td><td>0.0000*</td><td>0.6870</td><td>0.0000*</td><td>0.1156</td><td>0.3628</td><td>-0.0487</td><td>0.4870</td></tr><tr><td>MCV</td><td>Mean corpuscular volume (fl)</td><td>0.7396</td><td>0.0000*</td><td>0.7561</td><td>0.0000*</td><td>0.0733</td><td>0.5646</td><td>0.1344</td><td>0.0542</td></tr><tr><td>basophils</td><td>Basophils (10*9/L)</td><td>0.7245</td><td>0.0000*</td><td>0.6079</td><td>0.0000*</td><td>0.0254</td><td>0.8423</td><td>0.0044</td><td>0.9500</td></tr><tr><td>LDL</td><td>LDL cholesterol (mmol/L)</td><td>0.7096</td><td>0.0000*</td><td>0.7039</td><td>0.0000*</td><td>0.0651</td><td>0.6093</td><td>-0.0026</td><td>0.9699</td></tr><tr><td>MPV</td><td>Mean platelet volume (fl)</td><td>0.6957</td><td>0.0000*</td><td>0.7703</td><td>0.0000*</td><td>-0.0880</td><td>0.4893</td><td>-0.0399</td><td>0.5691</td></tr><tr><td>RDWSD</td><td>Red cell distribution width, standard variation</td><td>0.6936</td><td>0.0000*</td><td>0.6657</td><td>0.0000*</td><td>0.0309</td><td>0.8086</td><td>0.0558</td><td>0.4257</td></tr><tr><td>PDW</td><td>Platelet distribution width</td><td>0.6918</td><td>0.0000*</td><td>0.7456</td><td>0.0000*</td><td>-0.0893</td><td>0.4828</td><td>-0.0138</td><td>0.8437</td></tr><tr><td>PLT</td><td>Platelets (10*9/L)</td><td>0.6847</td><td>0.0000*</td><td>0.6708</td><td>0.0000*</td><td>0.1011</td><td>0.4265</td><td>-0.0413</td><td>0.5556</td></tr><tr><td>cholesterol</td><td>Total cholesterol (mmol/L)</td><td>0.6684</td><td>0.0000*</td><td>0.6464</td><td>0.0000*</td><td>-0.0132</td><td>0.9176</td><td>0.0593</td><td>0.3961</td></tr><tr><td>MCH</td><td>Mean corpuscular hemoglobin (fmol)</td><td>0.6662</td><td>0.0000*</td><td>0.7019</td><td>0.0000*</td><td>0.0012</td><td>0.9925</td><td>-0.0252</td><td>0.7189</td></tr><tr><td>creatinine</td><td>Creatinine (mol/L)</td><td>0.6657</td><td>0.0000*</td><td>0.6441</td><td>0.0000*</td><td>0.3282</td><td>0.0081</td><td>0.0972</td><td>0.1635</td></tr><tr><td>HR</td><td>Heart rate (bpm)</td><td>0.6395</td><td>0.0000*</td><td>0.3131</td><td>0.0001*</td><td>0.0378</td><td>0.7686</td><td>0.0119</td><td>0.8643</td></tr><tr><td>DBP</td><td>Diastolic blood pressure (mmHg)</td><td>0.6323</td><td>0.0000*</td><td>0.5213</td><td>0.0000*</td><td>0.0518</td><td>0.6868</td><td>-0.0558</td><td>0.4214</td></tr><tr><td>eosinophils</td><td>Eosinophils (10*9/L)</td><td>0.6233</td><td>0.0000*</td><td>0.4949</td><td>0.0000*</td><td>0.0134</td><td>0.9163</td><td>-0.0014</td><td>0.9838</td></tr><tr><td>WHR</td><td>Waist to hip ratio (cm/cm)</td><td>0.6026</td><td>0.0000*</td><td>0.4095</td><td>0.0000*</td><td>-0.1052</td><td>0.4121</td><td>-0.0575</td><td>0.4118</td></tr><tr><td>WBC</td><td>White blood cell count (10*9/L)</td><td>0.6001</td><td>0.0000*</td><td>0.5301</td><td>0.0000*</td><td>0.2986</td><td>0.0165</td><td>-0.0797</td><td>0.2634</td></tr><tr><td>lymphocytes</td><td>Lymphocytes (10*9/L)</td><td>0.5975</td><td>0.0000*</td><td>0.5353</td><td>0.0000*</td><td>0.2148</td><td>0.0909</td><td>-0.0489</td><td>0.4875</td></tr><tr><td>SBP</td><td>Systolic blood pressure (mmHg)</td><td>0.5948</td><td>0.0000*</td><td>0.5810</td><td>0.0000*</td><td>0.2139</td><td>0.0923</td><td>-0.0775</td><td>0.2636</td></tr><tr><td>insulin</td><td>Insulin (microlU/ml)</td><td>0.5925</td><td>0.0000*</td><td>0.3254</td><td>0.0001*</td><td>0.1678</td><td>0.2041</td><td>-0.0608</td><td>0.3940</td></tr><tr><td>fibrinogen</td><td>Fibrinogen (g/L)</td><td>0.5886</td><td>0.0000*</td><td>0.4423</td><td>0.0000*</td><td>-0.0080</td><td>0.9522</td><td>-0.0248</td><td>0.7350</td></tr><tr><td>monocytes</td><td>Monocytes (10*9/L)</td><td>0.5873</td><td>0.0000*</td><td>0.6474</td><td>0.0000*</td><td>0.0928</td><td>0.4693</td><td>0.1068</td><td>0.1314</td></tr><tr><td>neutrophils</td><td>Neutrophils (10*9/L)</td><td>0.5800</td><td>0.0000*</td><td>0.3687</td><td>0.0000*</td><td>0.2512</td><td>0.0453</td><td>-0.0529</td><td>0.4556</td></tr><tr><td>CRP</td><td>C-reactive protein (mg/L)</td><td>0.5730</td><td>0.0000*</td><td>0.3200</td><td>0.0002*</td><td>-0.1983</td><td>0.1287</td><td>0.0936</td><td>0.2012</td></tr><tr><td>RDWCV</td><td>Red cell distribution width, coefficient variation</td><td>0.5708</td><td>0.0000*</td><td>0.5665</td><td>0.0000*</td><td>-0.0895</td><td>0.4821</td><td>0.0752</td><td>0.2826</td></tr><tr><td>HDL</td><td>HDL cholesterol (mmol/L)</td><td>0.5697</td><td>0.0000*</td><td>0.7372</td><td>0.0000*</td><td>0.0498</td><td>0.6960</td><td>-0.0147</td><td>0.8331</td></tr><tr><td>GGT</td><td>Gamma glutamyl transpeptidase (U/L)</td><td>0.5630</td><td>0.0000*</td><td>0.5487</td><td>0.0000*</td><td>-0.0277</td><td>0.8277</td><td>0.0488</td><td>0.4853</td></tr><tr><td>MCHC</td><td>Mean corpuscular hemoglobin concentration (mmol/L)</td><td>0.5510</td><td>0.0000*</td><td>0.5076</td><td>0.0000*</td><td>-0.0445</td><td>0.7289</td><td>0.0966</td><td>0.1670</td></tr><tr><td>RBC</td><td>Red blood cell count (10*12/L)</td><td>0.5446</td><td>0.0000*</td><td>0.6880</td><td>0.0000*</td><td>0.2382</td><td>0.0580</td><td>0.0722</td><td>0.3035</td></tr><tr><td>glucose</td><td>Glucose (mmol/L)</td><td>0.5342</td><td>0.0000*</td><td>0.4509</td><td>0.0000*</td><td>0.2283</td><td>0.0794</td><td>-0.1023</td><td>0.1496</td></tr><tr><td>HGB</td><td>Hemoglobin (mmol/L)</td><td>0.5285</td><td>0.0000*</td><td>0.4730</td><td>0.0000*</td><td>0.1138</td><td>0.3708</td><td>0.0395</td><td>0.5735</td></tr><tr><td>HCT</td><td>Hematocrit (ratio)</td><td>0.4834</td><td>0.0002*</td><td>0.5137</td><td>0.0000*</td><td>0.0950</td><td>0.4553</td><td>0.0799</td><td>0.2548</td></tr></tbody></table>

(Continued)
<table><caption>**Table 2. (Continued)**</caption><thead><tr><th rowspan="2">Variable</th><th rowspan="2">Description</th><th colspan="2">Twin-Twin Living Together</th><th colspan="2">Twin-Twn Living Apart</th><th colspan="2">Twin-Spouse</th><th colspan="2">Unrelated</th></tr><tr><th>Correlation</th><th>P-Value</th><th>Correlation</th><th>P-Value</th><th>Correlation</th><th>P-Value</th><th>Correlation</th><th>P-Value</th></tr></thead><tbody><tr><td>ALT</td><td>Alanine transaminase (U/L)</td><td>0.4445</td><td>0.0006*</td><td>0.3756</td><td>0.0000*</td><td>0.2899</td><td>0.0201</td><td>0.0149</td><td>0.8313</td></tr><tr><td>AST</td><td>Aspartate aminotransferase (U/L)</td><td>0.3495</td><td>0.0077</td><td>0.3411</td><td>0.0000*</td><td>0.1013</td><td>0.4297</td><td>0.0137</td><td>0.8452</td></tr><tr><td>IFNg</td><td>Interferon gamma (pg/ml)</td><td>0.3270</td><td>0.0205</td><td>0.2404</td><td>0.0051</td><td>0.1498</td><td>0.2704</td><td>-0.0675</td><td>0.3611</td></tr><tr><td>IL8</td><td>Interleukin-8 (pg/ml)</td><td>0.3237</td><td>0.0193</td><td>0.4329</td><td>0.0000*</td><td>0.4530</td><td>0.0005*</td><td>-0.0256</td><td>0.7236</td></tr><tr><td>GMCSF</td><td>Granulocyte-macrophage colony-stimulating factor (pg/ml)</td><td>0.2370</td><td>0.0940</td><td>-0.0499</td><td>0.5600</td><td>-0.0292</td><td>0.8294</td><td>0.1685</td><td>0.0202</td></tr><tr><td>IL2</td><td>Interleukin-2 (pg/ml)</td><td>0.1461</td><td>0.2967</td><td>0.0960</td><td>0.2662</td><td>-0.1095</td><td>0.4263</td><td>0.1014</td><td>0.1651</td></tr><tr><td>IL6</td><td>Interleukin-6 (pg/ml)</td><td>0.1244</td><td>0.3846</td><td>0.0931</td><td>0.2703</td><td>-0.0508</td><td>0.7024</td><td>-0.0115</td><td>0.8736</td></tr><tr><td>TNFA</td><td>Tumor necrosis factor alpha (pg/ml)</td><td>0.0628</td><td>0.6651</td><td>0.1573</td><td>0.0655</td><td>0.1206</td><td>0.3851</td><td>0.0296</td><td>0.6872</td></tr><tr><td>IL1b</td><td>Interleukin-1 beta (pg/ml)</td><td>-0.0058</td><td>0.9682</td><td>0.1149</td><td>0.1894</td><td>0.1727</td><td>0.2117</td><td>-0.1439</td><td>0.0540</td></tr><tr><td>IL10</td><td>Interleukin-10 (pg/ml)</td><td>-0.0603</td><td>0.6773</td><td>-0.0322</td><td>0.7129</td><td>-0.0331</td><td>0.8107</td><td>-0.0425</td><td>0.5672</td></tr><tr><td>IL12p70</td><td>Interleukin-12p70 (pg/ml)</td><td>-0.1040</td><td>0.4769</td><td>-0.0482</td><td>0.5888</td><td>-0.1040</td><td>0.4632</td><td>-0.0093</td><td>0.9018</td></tr></tbody></table>

Significant p-values (p < 0.001) are denoted with a *.

doi:10.1371/journal.pone.0120898.t002

p-value = ns) and unrelated individuals (correlation coefficient = -0.08, p-value = ns), indicating genetic influences. These findings are consistent with previous work by Xu et al who studied the genetics and environmental effects on BP variability [35]. MZ correlations in a similar range were observed for liver enzymes, heart rate and GGT, but correlations were lower for the other liver enzymes and cytokines (Table 2).

Interestingly, several variables were more correlated in twins living together than those living apart, signifying the role that the environment plays for those traits. The strongest such finding is for cotinine, a biomarker indicating tobacco exposure, for which the correlation coefficient is 0.9 (p-value < 2.2e16) in twins living together and 0.34 (p-value < 2.2e16) in twins living apart. Several other traits including triglyceride levels, waist circumference, heart rate and insulin levels showed a higher concordance in twins living together than those living apart and might have a larger environmental impact.

When examining how correlated these individual traits are in twin-spouse pairs and unrelated individuals (Table 2), only IL8 levels were significantly correlated in twin-spouse pairs (correlation coefficient = 0.45 and p = 0.0005). A few other traits were nominally significantly correlated in twin-spouse pairs including creatinine levels (correlation coefficient = 0.29, p-value = 0.02), hip circumference (correlation coefficient = 0.3, p-value = 0.02), white blood cell count (correlation coefficient = 0.3, p-value = 0.02) and cotinine levels (correlation coefficient = 0.3, p-value = 0.02).

## Phenotype similarity across twins, spouses and unrelated individuals

Fig. 5 summarizes the phenotypic profile similarity for unrelated individuals, spouse pairs, non-cohabiting MZ twins and cohabiting MZ twins respectively. The 57 cohabiting MZ pairs are younger and more often male than the MZ twin pairs who live separately (Fig. 1C, mean 23.3 vs. 35.4 years old, 63% vs. 73% female). However, all the measurements in the analysis have been corrected for sex and age.

We observe a significantly higher similarity in the profiles of twins than in twin-spouse pairs (average profile correlation = 0.54 versus 0.06, p-value < 2.2e16), which suggests a strong
<figure><img src="image_8.png" /><figcaption>**Fig 5. Correlation of phenotypic profiles across twins, spouses and unrelated individuals.** This boxplot shows the pairwise correlations of phenotypic profiles between individuals: cohabiting MZ twins (57 pairs), MZ twins living apart (153 pairs), twin-spouse pairs (64 pairs) and non-related individuals (132596 pairs).</figcaption></figure>

doi:10.1371/journal.pone.0120898.g005

genetic influence on the observed phenotypic profiles. In order to validate this finding, we examined how closely twins would cluster together based on their phenotypic profiles and compared results on the actual data to those computed on a randomly permuted dataset (S2 Fig.). For nearly half of the twin pairs, the two members have most correlated profiles in the cohort relative to other individuals. In the majority of the pairs, starting with twin 1, twin 2 is in the top five most similar with regard to the phenotypic profile. Comparing the twin pairs who live together with the twin pairs who live apart, we observe a higher phenotypic profile correlation
in the cohabiting group (average correlation 0.58 and 0.52 respectively, p-value = 0.01) suggest-
ing an additional role for the environment in the phenotypic profile. The higher phenotypic
profile correlations for the spouse pairs versus the unrelated pairs (average correlation = 0.06
versus -0.0023, p-value = 0.0007) provide further support for this small role for household fac-
tors in influencing one's phenotypic profile.

DISCUSSION

In this study we use correlation approaches to explore the similarity of genetically related indi-
viduals or individuals sharing a household across a large set of metabolic, hematological and
immunological traits. In addition we identify known and novel relationships between the mea-
sured phenotypes. The unique design of the study, including twins as well as their spouses, al-
lows us to study the genetic and household effects of the measured phenotypic traits. We
define the concept of a phenotypic profile of an individual allowing us to look at multivariate
relationships simultaneously to determine the resemblance of monozygotic twins in compari-
son to their spouses and unrelated individuals.

Looking at the overall relationships among all traits several interesting clusters could be
identified. Many of the inflammatory markers were highly correlated across this healthy popu-
lation based cohort. The overall white blood cell count as well as levels of various types of white
blood cells (neutrophils, lymphocytes, monocytes) correlated with levels of IL6, a cytokine se-
creted by T-cells and macrophages to stimulate immune response. As expected, many of the
metabolic phenotypes including weight, BMI, and waist circumference, were also correlated
with each other and other traits such as LDL, total cholesterol, glucose and insulin levels.

Some traits such as heart rate, liver enzyme levels, fibrinogen and CRP were correlated with both metabolic (BMI, waist, weight, triglyceride and glucose levels) and inflammatory markers (IL6, WBC, neutrophils), which have also been seen by Il'yasova et al in 2008 [36]. The biological links between metabolism and inflammation, mostly through lipid-activated nuclear receptors (LXR) and peroxisome proliferator-activated receptors (PPARs) have been explored extensively in previous studies [37,38]. The relationship that we observe between fibrinogen and CRP has also been previously noted[39]. There are several novel relationships that are observed such as cotinine and inflammation (WBC, neutrophils, lymphocytes, monocytes). This study thus captures associations between metabolism and inflammation through phenotypic measurements of healthy individuals. This type of analysis can also be extended to individuals with disease to see how the network of phenotypic measurements and the relationships between them differ as a function of disease status or by therapeutic interventions.

The unique design of this study allows us to determine the relative contribution of genetic and household factors to familial resemblance in the metabolic, hematological and immuno-logical traits. In general, when comparing the intra-pair resemblances of the traits of cohabiting MZ twins, resemblance can be a function of both their genes and shared environment, whereas the resemblance of non-cohabiting MZ twins provides an upper limit of the heritability of the trait. In twin-spouse pairs, who share environment only, the resemblance is caused by house-hold factors. Finally, we compare twin-twin and twin-spouse correlations to those obtained in unrelated individuals who share neither genes nor environment.

Results of the comparison of similarity across MZ twin pairs, spouse pairs and pairs of unre-
lated individuals clearly demonstrated the importance of genetic factors on phenotypic profiles
as well as many of the individual phenotypic traits in healthy individuals. The profiles among
MZ twins, either cohabiting or non-cohabiting, were significantly higher correlated than those
of twin-spouse pairs or unrelated individuals, pointing to the importance of genetic factors al-
though a lingering role for childhood environmental factors in non-cohabiting twins cannot be
excluded. Some of the most concordant traits included height, telomere length, IL6 receptor levels, and platelet measures. One of the few traits that showed high correspondence in spouse pairs was cotinine level. This may not be surprising, since cotinine levels are reflective of smoking behavior. Concordance rates for smoking are generally high in spouses, both before and after marriage [40] and smoking is one of the cardiovascular risk factors with highest concordance among spouses [41].

The results demonstrate that substantial differences may exist in the relative contribution of genetic factors to the various components of a single biological system. The interaction between circulating IL6 and its soluble receptor is a case in point. IL6 is a potent pleiotropic cytokine secreted by T-cells and macrophages that regulates cell growth and differentiation and plays an important role in immune response. Interleukin 6 receptor (IL6R) forms a ligand-receptor complex with IL-6 that is capable of stimulating a variety of cellular responses including proliferation, differentiation and activation of inflammatory processes [42,43]. Raggi et al previously showed similar high heritability for IL6R [44], while other studies reported low to moderate heritability for IL6. IL6R level seems to be a stable individual characteristic with a genetic component [45], while IL6 level may be more responsive to environmental conditions.

A few extensions and limitations of our approach should be mentioned. The current paper provides a proof of principle in the proposed methodology combined with the unique study design, but includes only a snapshot of the variables available. However, the number of phenotypic traits that are examined here can be extended in number as well as towards categorical variables. The development of integrated phenotypic databases and increased use of electronic medical record systems will allow for richer and more comprehensive analyses to be carried out. In addition, longitudinal data can be incorporated. Following phenotypic markers in twins from birth, through adolescence, adulthood and old age and studying the change in their phenotypic landscape though time will be a very valuable extension of this work. Furthermore, validating some of the relationships and correlations presented here in a larger independent cohort would strengthen the findings. The current study compared MZ twin pairs, who share their genes and childhood environment, with spouse pairs, who share their adult environment. The impact of the shared environment in spouse pairs may differ depending on the time the spouses have been together and the extent to which they take part in similar activities. It would be interesting to investigate the similarity of the phenotypic profile as function of shared interests and cohabitation duration in a larger sample. Similarly we assume that sharing a household in adulthood implies that shared environment can affect phenotype. The effect of the shared as well as non-shared environment of the individuals at a young age or in utero is not captured by our study and is a limitation of the approach. Because there is such limited evidence from studies using the classical twin design (comparing MZ to DZ resemblance) that early common environment has an influence on the traits we study, we explored if the current shared environment has an influence. Both approaches seem to lead to similar conclusions: a main effect of shared environment is very hard to detect.

In conclusion, we present here a comprehensive analysis of the phenotypic profile and relationships between anthropometric, hematological, metabolic, and inflammatory traits of interest in twin pairs and spouse pairs. We show complex relationships between many of the phenotypic variables, forming metabolic and inflammatory clusters. The increased twin similarity demonstrates the importance of genetic factors while the larger concordance in spouse pairs compared to unrelated pairs and the comparison of cohabiting and non-cohabiting twins indicates that specific household factors may also play a role, albeit small, in the metabolic and immunological profile.
**Supporting Information**

**S1 Fig. The distribution of phenotypic measurements (shown as sex-and age corrected residuals) in the total sample (N = 51).**
(PDF)

**S2 Fig. How closely do twins cluster together in comparison to unrelated individuals?** For each twin pair the rank of how similar the other twin's profile is in comparison to all the other phenotypic profiles in the cohort is displayed. Panel A shows the results on the actual data, panel B shows the results on a representative permutation of the data. For the majority of the twins, their twin pair falls within the top five most similar individuals based on their phenotypic profile.
(PDF)

**S1 Table. Partial Pairwise Correlations between Traits.** This table contains all the significant partial pairwise correlation values between the measured phenotypic traits. The non-significant values were set to 0.
(XLSX)

**S2 Table. Pair-wise correlations and confidence intervals for each trait in each of the four comparison groups (MZ twin pairs living together, MZ twin pairs living apart, twin-spouse pairs, unrelated individuals) ordered by concordance in twins living together.**
(XLSX)

**Acknowledgments**

We would like to dedicate this manuscript to the memory of David R. Cox, our dear friend and colleague, a visionary and passionate scientist who has had a profound influence on the fields of genetics and biomedicine. We will always be inspired by David's lifelong dedication to biomedical research and the use of human genetics to improve human health. We would like to thank all twins and family members for their participation in the Netherlands Twin Register Biobank project. We thank Gareth Davies, Erik Ehli and Patrizia Huizenga from Avera Institute for Human Genetics, Sioux Falls, USA for their assistance with the zygosity determination.

**Author Contributions**

Conceived and designed the experiments: MS GW DRC EJCG DIB. Performed the experiments: MS GW JN CK KEM. Analyzed the data: MS GW. Contributed reagents/materials/ analysis tools: PS SP SJP EP. Wrote the paper: MS GW PS SP SJP EP SK SDE JN CK KEM DRC EJCG DIB.

**REFERENCES**

1. McCartney K, Harris MJ, Bernieri F. Growing up and growing apart: a developmental meta-analysis of twin studies. Psychological bulletin. 1990; 107(2):226-37. PMID: 2138795

2. Pilia G, Chen WM, Scuteri A, Orru M, Albai G, Dei M, et al. Heritability of cardiovascular and personality traits in 6,148 Sardinians. PLoS genetics. 2006; 2(8):e132. PMID: 16934002

3. van Dongen J, Slagboom PE, Draisma HH, Martin NG, Boomsma DI. The continuing value of twin studies in the omics era. Nature reviews Genetics. 2012; 13(9):640-53. doi: 10.1038/nrg3243 PMID: 22847273

4. Boomsma D, Busjahn A, Peltonen L. Classical twin studies and beyond. Nature reviews Genetics. 2002; 3(11):872-82. PMID: 12415317

5. van Rijn MJ, Schut AF, Aulchenko YS, Deinum J, Sayed-Tabatabaei FA, Yazdanpanah M, et al. Heritability of blood pressure traits and the genetic contribution to blood pressure variance explained by four blood-pressure-related genes. Journal of hypertension. 2007; 25(3):565-70. PMID: 17278972
6. Wu T, Treiber FA, Snieder H. Genetic influence on blood pressure and underlying hemodynamics measured at rest and during stress. Psychosomatic medicine. 2013; 75(4):404–12. doi: 10.1097/PSY. 0b013e31828d3cb6 PMID: 23576770

7. Schousboe K, Willemsen G, Kyvik KO, Mortensen J, Boomsma DI, Cornes BK, et al. Sex differences in heritability of BMI: a comparative study of results from twin studies in eight countries. Twin research: the official journal of the International Society for Twin Studies. 2003; 6(5):409–21. PMID: 14624725

8. Korkeila M, Kaprio J, Rissanen A, Koskenvuo M. Consistency and change of body mass index and weight. A study on 5967 adult Finnish twin pairs. International journal of obesity and related metabolic disorders: journal of the International Association for the Study of Obesity. 1995; 19(5):310–7.

9. Turula M, Kaprio J, Rissanen A, Koskenvuo M. Body weight in the Finnish Twin Cohort. Diabetes research and clinical practice. 1990; 10 Suppl 1:S33–6. PMID: 2286148

10. Min J, Chiu DT, Wang Y. Variation in the heritability of body mass index based on diverse twin studies: a systematic review. Obesity reviews: an official journal of the International Association for the Study of Obesity. 2013; 14(11):871–82.

11. Neijs M, van Dongen J, Kluft C, Boomsma DI, Willemsen G, de Geus EJ. Genetic architecture of the pro-inflammatory state in an extended twin-family design. Twin research and human genetics: the official journal of the International Society for Twin Studies. 2013; 16(5):931–40.

12. Su S, Snieder H, Miller AH, Ritchie J, Bremner JD, Goldberg J, et al. Genetic and environmental influences on systemic markers of inflammation in middle-aged male twins. Atherosclerosis. 2008; 200 (1):213–20. doi: 10.1016/j.atherosclerosis.2007.12.009 PMID: 18243214

13. van Dongen J, Willemsen G, Chen WM, de Geus EJ, Boomsma DI. Heritability of metabolic syndrome traits in a large population-based sample. Journal of lipid research. 2013; 54(10):2914–23. doi: 10. 1194/jlr.P041673 PMID: 23918046

14. Fragoso TM, Giolo SR, Pereira AC, de Andrade M, Soler JM. Using item response theory to model multiple phenotypes and their joint heritability in family data. Genetic epidemiology. 2014; 38(2):152–61. doi: 10.1002/gepi.21784 PMID: 24415554

15. Burgel PR, Paillasseur JL, Caillaud D, Tillie-Leblond I, Chanez P, Escamilla R, et al. Clinical COPD phenotypes: a novel approach using principal component and cluster analyses. The European respiratory journal. 2010; 36(3):531–9. doi: 10.1183/09031936.00175109 PMID: 20075045

16. Moore WC, Meyers DA, Wenzel SE, Teague WG, Li H, Li X, et al. Identification of asthma phenotypes using cluster analysis in the Severe Asthma Research Program. American journal of respiratory and critical care medicine. 2010; 181(4):315–23. doi: 10.1164/rccm.200906-0896OC PMID: 19892860

17. Patrawalla P, Kazeros A, Rogers L, Shao Y, Liu M, Fernandez-Beros ME, et al. Application of the asthma phenotype algorithm from the Severe Asthma Research Program to an urban population. PloS one. 2012; 7(9):e44540. doi: 10.1371/journal.pone.0044540 PMID: 23028556

18. Draisma HH, Beekman M, Pool R, van Ommen GJ, Adamski J, Prehn C, et al. Familial resemblance for serum metabolite concentrations. Twin research and human genetics: the official journal of the International Society for Twin Studies. 2013; 16(5):948–61.

19. Boomsma DI, de Geus EJ, Vink JM, Stubbe JH, Distel MA, Hottenga JJ, et al. Netherlands Twin Register: from twins to twin families. Twin research and human genetics: the official journal of the International Society for Twin Studies. 2006; 9(6):849–57.

20. Willemsen G, Vink JM, Abdellaoui A, den Braber A, van Beek JH, Draisma HH, et al. The Adult Netherlands Twin Register: twenty-five years of survey and biological data collection. Twin research and human genetics: the official journal of the International Society for Twin Studies. 2013; 16(1):271–81.

21. Willemsen G, de Geus EJ, Bartels M, van Beijsterveldt CE, Brooks AI, Estourgie-van Burk GF, et al. The Netherlands Twin Register biobank: a resource for genetic epidemiological studies. Twin research and human genetics: the official journal of the International Society for Twin Studies. 2010; 13(3):231–45.

22. Friedewald WT, Levy RI, Fredrickson DS. Estimation of the concentration of low-density lipoprotein cholesterol in plasma, without use of the preparative ultracentrifuge. Clinical chemistry. 1972; 18 (6):499–502. PMID: 4337382

23. Hottenga JJ, Whitfield JB, Posthuma D, Willemsen G, de Geus EJ, Martin NG, et al. Genome-wide scan for blood pressure in Australian and Dutch subjects suggests linkage at 5P, 14Q, and 17P. Hypertension. 2007; 49(4):832–8. PMID: 17325240

24. Ray KK, Cannon CP, Braunwald E. Recent trials of lipid lowering. International journal of clinical practice. 2007; 61(7):1145–59. PMID: 17577297

25. Defays D. An efficient algorithm for a complete link method. The Computer Journal. 1977; 20(4):364–6.

26. Ohmori Y, Schreiber RD, Hamilton TA. Synergy between interferon-gamma and tumor necrosis factor-alpha in transcriptional activation is mediated by cooperation between signal transducer and activator
of transcription 1 and nuclear factor kappaB. The Journal of biological chemistry. 1997; 272 (23):14899-907. PMID: 9169460

27. Shannon MF, Himes SR, Coles LS. GM-CSF and IL-2 share common control mechanisms in response to costimulatory signals in T cells. Journal of leukocyte biology. 1995; 57(5):767-73. PMID: 7759956

28. O'Garra A, Murphy KM. From IL-10 to IL-12: how pathogens and their products stimulate APCs to in- duce T(H)1 development. Nature immunology. 2009; 10(9):929-32. doi: 10.1038/ni0909-929 PMID: 19692989

29. Brito Diaz B, Aleman Sanchez JJ, Cabrera de Leon A. [Resting heart rate and cardiovascular disease.]. Medicina clinica. 2013.

30. Valensi P, Extramiana F, Lange C, Cailleau M, Haggui A, Maison Blanche P, et al. Influence of blood glucose on heart rate and cardiac autonomic function. The DESIR study. Diabetic medicine: a journal of the British Diabetic Association. 2011; 28(4):440-9. doi: 10.1111/j.1464-5491.2010.03222.x PMID: 21204961

31. Silventoinen K. Determinants of variation in adult body height. Journal of biosocial science. 2003; 35 (2):263-85. PMID: 12664962

32. Benetos A, Okuda K, Lajemi M, Kimura M, Thomas F, Skurnick J, et al. Telomere length as an indicator of biological aging: the gender effect and relation with pulse pressure and pulse wave velocity. Hyper- tension. 2001; 37(2 Pt 2):381-5. PMID: 11230304

33. Insel KC, Merkle CJ, Hsiao CP, Vidrine AN, Montgomery DW. Biomarkers for cognitive aging part I: telomere length, blood pressure and cognition among individuals with hypertension. Biological research for nursing. 2012; 14(2):124-32. doi: 10.1177/1099800411406433 PMID: 21586494

34. Broer L, Codd V, Nyholt DR, Deelen J, Mangino M, Willemsen G, et al. Meta-analysis of telomere length in 19,713 subjects reveals high heritability, stronger maternal inheritance and a paternal age effect. Eu- ropean journal of human genetics: EJHG. 2013; 21(10):1163-8. doi: 10.1038/ejhg.2012.303 PMID: 23321625

35. Xu X, Ding X, Zhang X, Su S, Treiber FA, Vlietinck R, et al. Genetic and environmental influences on blood pressure variability: a study in twins. Journal of hypertension. 2013; 31(4):690-7. doi: 10.1097/ HJH.0b013e32835e2a4a PMID: 23470779

36. Il'yasova D, Ivanova A, Morrow JD, Cesari M, Pahor M. Correlation between two markers of inflamma- tion, serum C-reactive protein and interleukin 6, and indices of oxidative stress in patients with high risk of cardiovascular disease. Biomarkers: biochemical indicators of exposure, response, and susceptibili- ty to chemicals. 2008; 13(1):41-51.

37. Bensinger SJ, Tontonoz P. Integration of metabolism and inflammation by lipid-activated nuclear recep- tors. Nature. 2008; 454(7203):470-7. doi: 10.1038/nature07202 PMID: 18650918

38. Hong C, Tontonoz P. Coordination of inflammation and metabolism by PPAR and LXR nuclear recep- tors. Current opinion in genetics & development. 2008; 18(5):461-7.

39. Schnabel RB, Lunetta KL, Larson MG, Dupuis J, Lipinska I, Rong J, et al. The relation of genetic and environmental factors to systemic inflammatory biomarker concentrations. Circulation Cardiovascular genetics. 2009; 2(3):229-37. doi: 10.1161/CIRCGENETICS.108.804245 PMID: 20031590

40. Merline AC, Schulenberg JE, O'Malley PM, Bachman JG, Johnston LD. Substance use in marital dyads: premarital assortment and change over time. Journal of studies on alcohol and drugs. 2008; 69 (3):352-61. PMID: 18432377

41. Di Castelnuovo A, Quacquaruccio G, Donati MB, de Gaetano G, Iacoviello L. Spousal concordance for major coronary risk factors: a systematic review and meta-analysis. American journal of epidemiology. 2009; 169(1):1-8. doi: 10.1093/aje/kwn234 PMID: 18845552

42. Jones SA, Horiuchi S, Topley N, Yamamoto N, Fuller GM. The soluble interleukin 6 receptor: mecha- nisms of production and implications in disease. FASEB journal: official publication of the Federation of American Societies for Experimental Biology. 2001; 15(1):43-58.

43. Mihara M, Hashizume M, Yoshida H, Suzuki M, Shiina M. IL-6/IL-6 receptor system and its role in physi- ological and pathological conditions. Clinical science. 2012; 122(4):143-59. doi: 10.1042/CS20110340 PMID: 22029668

44. Raggi P, Su S, Karohl C, Veledar E, Rojas-Campos E, Vaccarino V. Heritability of renal function and in- flammatory markers in adult male twins. American journal of nephrology. 2010; 32(4):317-23. doi: 10. 1159/000319449 PMID: 20720405

45. van Dongen J, Jansen R, Smit D, Hottenga JJ, Mbarek H, Willemsen G, et al. The Contribution of the Functional IL6R Polymorphism rs2228145, eQTLs and Other Genome-Wide SNPs to the Heritability of Plasma sIL-6R Levels. Behavior genetics. 2014.