Human bladder uroepithelial cells synergize with monocytes to promote IL-10 synthesis and other cytokine responses to Uropathogenic Escherichia coli

Urinary tract infections are a major source of morbidity for women and the elderly, with Uropathogenic Escherichia coli (UPEC) being the most prevalent causative pathogen. Studies in recent years have defined a key anti-inflammatory role for Interleukin-10 (IL-10) in urinary tract infection mediated by UPEC and other uropathogens. We investigated the nature of the IL-10-producing interactions between UPEC and host cells by utilising a novel co-culture model that incorporated lymphocytes, mononuclear and uroepithelial cells in histotypic proportions. This co-culture model demonstrated synergistic IL-10 production effects between monocytes and uroepithelial cells following infection with UPEC. Membrane inserts were used to separate the monocyte and uroepithelial cell types during infection and revealed two synergistic IL-10 production effects based on contact-dependent and soluble interactions. Analysis of a comprehensive set of immunologically relevant biomarkers in monocyte-uroepithelial cell co-cultures highlighted that multiple cytokine, chemokine and signalling factors were also produced in a synergistic or antagonistic fashion. These results demonstrate that IL-10 responses to UPEC occur via multiple interactions between several cells types, implying a complex role for infection-related IL-10 during UTI. Development and application of the co-culture model described in this study is thus useful to define the degree of contact dependency of biomarker production to UPEC, and highlights the relevance of histotypic co-cultures in studying complex host-pathogen interactions.

Mary Robinson's declaration of climate justice: Climate change, human rights and fossil fuel divestment

In her biography, Everybody Matters: My Life Giving Voice, Mary Robinson explained how she became interested in the topic of human rights and climate change, after hearing testimony from African farmers, with Archbishop Desmond Tutu.

p53 and DNA damage checkpoint responses in the human prostate

Prostate cancer is the most common noncutaneous malignancy and the second leading cause of cancer mortality in men. In 2004, 5237 new cases were diagnosed and altogether 25 664 men suffered from prostate cancer in Finland (Suomen Syöpärekisteri). Although extensively investigated, we still have a very rudimentary understanding of the molecular mechanisms leading to the frequent transformation of the prostate epithelium. Prostate cancer is characterized by several unique features including the multifocal origin of tumors and extreme resistance to chemotherapy, and new treatment options are therefore urgently needed. The integrity of genomic DNA is constantly challenged by genotoxic insults. Cellular responses to DNA damage involve elegant checkpoint cascades enforcing cell cycle arrest, thus facilitating damage repair, apoptosis or cellular senescence. Cellular DNA damage triggers the activation of tumor suppressor protein p53 and Wee1 kinase which act as executors of the cellular checkpoint responses. These are essential for genomic integrity, and are activated in early stages of tumorigenesis in order to function as barriers against tumor formation. Our work establishes that the primary human prostatic epithelial cells and prostatic epithelium have unexpectedly indulgent checkpoint surveillance. This is evidenced by the absence of inhibitory Tyr15 phosphorylation on Cdk2, lack of p53 response, radioresistant DNA synthesis, lack of G1/S and G2/M phase arrest, and presence of persistent gammaH2AX damage foci. We ascribe the absence of inhibitory Tyr15 phosphorylation to low levels of Wee1A, a tyrosine kinase and negative regulator of cell cycle progression. Ectopic Wee1A kinase restored Cdk2-Tyr15 phosphorylation and efficiently rescued the ionizing radiation-induced checkpoints in the human prostatic epithelial cells. As variability in the DNA damage responses has been shown to underlie susceptibility to cancer, our results imply that a suboptimal checkpoint arrest may greatly increase the accumulation of genetic lesions in the prostate epithelia. We also show that small molecules can restore p53 function in prostatic epithelial cells and may serve as a paradigm for the development of future therapeutic agents for the treatment of prostate cancer We hypothesize that the prostate has evolved to activate the damage surveillance pathways and molecules involved in these pathways only to certain stresses in extreme circumstances. In doing so, this organ inadvertently made itself vulnerable to genotoxic stress, which may have implications in malignant transformation. Recognition of the limited activity of p53 and Wee1 in the prostate could drive mechanism-based discovery of preventative and therapeutic agents.

Role of Basement Membrane and Extracellular Matrix Proteins in the Adhesion and Spreading of Immortalized Human Corneal Epithelial Cells

The repair of corneal wounds requires both epithelial cell adhesion and migration. Basement membrane (BM) and extracellular matrix (ECM) proteins function in these processes via integrin and non-integrin receptors. We have studied the adhesion, spreading and migration of immortalized human corneal epithelial (HCE) cells and their interactions with the laminins (Lms), fibronectins and tenascins produced. Human corneal BM expresses Lms-332 and -511, while Lm-111 was not found in these experiments. HCE cells produced both processed and unprocessed Lm-332, whereas neither Lm-111 nor Lm-511 was produced. Because HCE cells did not produce Lm-511, although it was present in corneal BM, we suggest that Lm-511 is produced by stromal keratocytes. The adhesion of HCE cells to Lms-111, -332 and -511 was studied first by determining the receptor composition of HCE cells and then by using quantitative cell adhesion assays. Immunofluorescence studies revealed the presence of integrin α2, α3, α6, β1 and β4 subunits. Among the non-integrin receptors, Lutheran (Lu) was found on adhering HCE cells. The cells adhered via integrin α3β1 to both purified human Lms-332 and -511 as well as to endogenous Lm-332. However, only integrin β1 subunit functioned in HCE cell adhesion to mouse Lm-111. The adhesion of HCE cells to Lm-511 was also mediated by Lu. Since Lm-511 did not induce Lu into focal adhesions in HCE cells, we suggest that Lm-511 serves as an ECM ligand enabling cell motility. HCE cells produced extradomain-A fibronectin, oncofetal fibronectin and tenascin-C (Tn-C), which are also found during corneal wound healing. Monoclonal antibodies (MAbs) against integrins α5β1 and αvβ6 as well as the arginine-glycine-aspartic acid (RGD) peptide inhibited the adhesion of HCE cells to fibronectin. Although the cells did not adhere to Tn-C, they adhered to the fibronectin/Tn-C coat and were then more efficiently inhibited by the function-blocking MAbs and RGD peptide. During the early adhesion, HCE cells codeposited Lm-332 and the large subunit of tenascin-C (Tn-CL) beneath the cells via the Golgi apparatus and microtubules. Integrin β4 subunit, which is a hemidesmosomal component, did not mediate the early adhesion of HCE cells to Lm-332 or Lm-332/Tn-C. Based on these results, we suggest that the adhesion of HCE cells is initiated by Lm-332 and modulated by Tn-CL, as it has been reported to prevent the assembly of hemidesmosomes. Thereby, Tn-CL functions in the motility of HCE cells during wound healing. The different distribution of processed and unprocessed Lm-332 in adhering, spreading and migrating HCE cells suggests a distinct role for these isoforms. We conclude that the processed Lm-332 functions in cell adhesion, whereas the unprocessed Lm-332 participates in cell spreading and migration.

Associations among Emdogain®, human oral carcinoma cells and oral proteolytic enzymes

The Enamel matrix derivative Emdogain® (EMD) is a commercially available tissue extract preparation of porcine enamel origin. Studies have shown EMD to be clinically useful in promoting periodontal regeneration. EMD has been widely used in periodontal therapy for over ten years, but the mechanism of its action and the exact composition are not completely clear. EMD is predominantly amelogenin (>90%). However, unlike amelogenin, EMD has a number of growth factor-like effects and it has been shown to enhance the proliferation, migration and other cellular functions of periodontal ligament fibroblasts and osteoblasts. In contrast, the effects of EMD on epithelial cell lines and in particular on oral malignant cells have not been adequately studied. In addition, EMD has effects on the production of cytokines by several oral cell lines and the product is in constant interaction with different oral enzymes. Regardless of the various unknown properties of EMD, it is said to be clinically safe in regenerative procedures, also in medically compromised patients. The aim of the study was to examine whether gingival crevicular fluid (GCF), which contains several different proteolysis enzymes, could degrade EMD and alter its biological functions. In addition, the objective was to study the effects of EMD on carcinogenesis-related factors, in particular MMPs, using in vitro and in vivo models. This study also aimed to contribute to the understanding of the composition of EMD. GCF was capable of degrading EMD, depending on the periodontal status, with markedly more degradation in all states of periodontal disease compared to healthy controls. EMD was observed to stimulate the migration of periodontal ligament fibroblasts (PLF), whereas EMD together with GCF could not stimulate this proliferation. In addition, recombinant amelogenin, the main component of EMD, decreased the migration of PLFs. A comparison of changes induced by EMD and TGF-β1 in the gene profiles of carcinoma cells showed TGF-β1 to regulate a greater number of genes than EMD. However, both of the study reagents enhanced the expression of MMP-10 and MMP-9. Furthermore, EMD was found to induce several factors closely related to carcinogenesis on gene, protein, cell and in vivo levels. EMD enhanced the production of MMP-2, MMP-9 and MMP-10 proteins by cultured carcinoma cells. In addition, EMD stimulated the migration and in vitro wound closure of carcinoma cells. EMD was also capable of promoting metastasis formation in mice. In conclusion, the diseased GCF, containing various proteases, causes degradation of EMD and decreased proliferation of PLFs. Thus, this in vitro study suggests that the regenerative effect of EMD may decrease due to proteases present in periodontal tissues during the inflammation and healing of the tissues in vivo. Furthermore, EMD was observed to enhance several carcinoma-related factors and in particular the production of MMPs by benign and malignant cell lines. These findings suggest that the clinical safety of EMD with regard to dysplastic mucosal lesions should be further investigated.

Characterization of cytokines, matrix metalloproteinases and toll-like receptors in human periodontal tissue destruction

Periodontal Disease affects the supporting structures of the teeth and is initiated by a microbial biofilm called dental plaque. Severity ranges from superficial inflammation of the gingiva (gingivitis) to extensive destruction of connective tissue and bone leading to tooth loss (periodontitis). In periodontitis the destruction of tissue is caused by a cascade of microbial and host factors together with proteolytic enzymes. Matrix metalloproteinases (MMPs) are known to be central mediators of the pathologic destruction in periodontitis. Initially plaque bacteria provide pathogen-associated molecular patterns (PAMPs) which are sensed by Toll-like receptors (TLRs), and initiate intracellular signaling cascades leading to host inflammation. Our aim was to characterize TNF-α (tumor necrosis factor-alpha) and its type I and II receptors in periodontal tissues, as well as, the effects of TNF-α, IL-1β (interleukin-1beta) and IL-17 on the production and/or activation of MMP-3, MMP-8 and MMP-9. Furthermore we mapped the TLRs in periodontal tissues and assessed how some of the PAMPs binding to the key TLRs found in periodontal tissues affect production of TNF-α and IL-1β by gingival epithelial cells with or without combination of IL-17. TNF-α and its receptors were detected in pericoronitis. Furthermore, increased expression of interleukin-1β and vascular cell adhesion molecule-1 was found as a biological indicator of TNF-α ligand-receptor interaction. MMP-3, -8, and 9 were investigated in periodontitis affected human gingival crevicular fluid and gingival fibroblasts produced pro-MMP-3. Following that, the effect of IL-17 was studied on MMP and pro-inflammatory cytokine production. IL-17 was increased in periodontitis and up-regulated IL-1β, TNF-α, MMP-1 and MMP-3. We continued by demonstrating TLRs in gingival tissues, in which significant differences between patients with periodontitis and healthy controls were found. Finally, enzyme-linked immunosorbent assays were performed to show that the gingival cells response to inflammatory responses in a TLR-dependent manner. Briefly, this thesis demonstrates that TLRs are present in periodontal tissues and present differences in periodontitis compared to healthy controls. The cells of gingival tissues respond to inflammatory process in a TLR-dependent manner by producing pro-inflammatory cytokines. During the destruction of periodontal tissues, the release (IL-1β and TNF-α) and co-operation with other pro-inflammatory cytokines (IL-17), which in turn increase the inflammation and thus be more harmful to the host with the increased presence of MMPs (MMP-1, MMP-3, MMP-8, MMP-9) in diseased over healthy sites.

Periodontitis and peri-implantitis biomarkers in human oral fluids and the null-allele mouse model

Tissue destruction associated with the periodontal disease progression is caused by a cascade of host and microbial factors and proteolytic enzymes. Aberrant laminin-332 (Ln-332), human beta defensin (hBD), and matrix metalloproteinase (MMP) functions have been found in oral inflammatory diseases. The null-allele mouse model appears as the next step in oral disease research. The MMP-8 knock-out mouse model allowed us to clarify the involvement of MMP-8 in vivo in oral and related inflammatory diseases where MMP-8 is suggested to play a key role in tissue destruction. The cleaved Ln-332 γ2-chain species has been implicated in the apical migration of sulcular epithelial cells during the formation of periodontal pockets. We demonstrated that increased Ln-332 fragment levels in gingival crevicular fluid (GCF) are strongly associated with the severity of inflammation in periodontitis. Porphyromonas gingivalis trypsin-like proteinase can cleave an intact Ln-332 γ2-chain into smaller fragments and eventually promote the formation of periodontal pockets. hBDs are components of an innate mucosal defense against pathogenic microbes. Our results suggest that P. gingivalis trypsin-like proteinase can degrade hBD and thus reduce the innate immune response. Elevated levels and the increased activity of MMPs have been detected in several pathological tissue-destructive conditions where MMPs are shown to cleave extracellular matrix (ECM) and basement membrane (BM) molecules and to facilitate tissue destruction. Elevated levels of MMP-8 have been reported in many inflammatory diseases. In periodontitis, MMP-8 levels in gingival crevicular fluid (GCF) and in peri-implant sulcular fluid (PISF) are elevated at sites of active inflammation, and the increased levels of MMP-8 are mainly responsible for collagenase activity, which leads to tissue destruction. MMP-25, expressed by neutrophils, is involved in inflammatory diseases and in ECM turnover. MMP-26 can degrade ECM components and serve as an activator of other MMP enzymes. We further confirmed that increased levels and activation of MMP-8, -25, and -26 in GCF, PISF, and inflamed gingival tissue are associated with the severity of periodontal/peri-implant inflammation. We evaluated the role of MMP-8 in P. gingivalis-induced periodontitis by comparing MMP-8 knock-out (MMP8-/-) and wild-type mice. Surprisingly, MMP-8 significantly attenuated P. gingivalis-induced site-specific alveolar bone loss. We also evaluated systemic changes in serum immunoglobulin and lipoprotein profiles among these mouse groups. P. gingivalis infection increased HDL/VLDL particle size in the MMP-8-/- mice, which is an indicator of lipoprotein responses during systemic inflammation. Serum total LPS and IgG antibody levels were enhanced in both mice groups. P. gingivalis-induced periodontitis, especially in MMP-8-/- mice, is associated with severe alveolar bone loss and with systemic inflammatory and lipoprotein changes that are likely to be involved in early atherosclerosis.

Defining the role of common variation in the genomic and biological architecture of adult human height

Using genome-wide data from 253,288 individuals, we identified 697 variants at genome-wide significance that together explained one-fifth of the heritability for adult height. By testing different numbers of variants in independent studies, we show that the most strongly associated approximately 2,000, approximately 3,700 and approximately 9,500 SNPs explained approximately 21%, approximately 24% and approximately 29% of phenotypic variance. Furthermore, all common variants together captured 60% of heritability. The 697 variants clustered in 423 loci were enriched for genes, pathways and tissue types known to be involved in growth and together implicated genes and pathways not highlighted in earlier efforts, such as signaling by fibroblast growth factors, WNT/beta-catenin and chondroitin sulfate-related genes. We identified several genes and pathways not previously connected with human skeletal growth, including mTOR, osteoglycin and binding of hyaluronic acid. Our results indicate a genetic architecture for human height that is characterized by a very large but finite number (thousands) of causal variants.

Mast cells and endothelial erosion: Implications for human atherothrombotic disease

More than 40% of all deaths in Finland are caused by atherosclerosis. The complications of atherosclerosis are due to either detachment of the luminal endothelium (erosion) or rupture of the fibrous cap of an atherosclerotic plaque (rupture). As a result, a thrombus is formed at the site of the intimal lesion. Indeed, erosions cause roughly 40% of sudden atherothrombotic deaths and 25% of all atherothrombotic deaths. Erosions are overrepresented in young subjects, diabetics, smokers and women. This dissertation focuses on endothelial erosion. Endothelial erosions were studied in the context of arterial grafting and vascular inflammation. Special attention was given to the role of intimal mast cells and the methodological viewpoints of reliable identification of endothelial erosions. Mast cells are inflammatory cells mostly known for their ability to cause allergic symptoms. In addition to occurring in skin and mucosal surfaces, mast cells are abundant in arterial intima and adventitia. In this study, mast cells were found to associate with endothelial erosions in non-lesional and atherosclerotic human coronary arteries. Thus, mast cells may participate in atherogenesis at the initial phases of the disease process already. We also showed that the mast cell proteases tryptase, chymase, and cathepsin G are all capable of cleaving molecules essential for endothelial cell-to-cell and cell-to-extracellular matrix interactions, such as VE-cadherin and fibronectin. Symptom-causing carotid plaques were found to contain more inflammatory cells, especially mast cells, than non-symptom-causing plaques. Furthermore, the atherogenic serum lipid profile and the degree of carotid stenosis turned out to correlate with the density of carotid plaque mast cells. Apoptotic and proliferating cells were more abundant in non-symptom causing plaques (active renewal of endothelial cells), but erosions were larger in symptom-causing plaques (capacity of endothelial regeneration exceeded). The process of identifying endothelial erosions with immunostainings has been ambiguous, since both endothelial cells and platelets express largely the same antigens. This may have caused inaccurate interpretations of the presence of endothelial erosion. In the last substudy of this thesis we developed a double immunostaining method for simultaneous identification of endothelial cells and platelets. This method enables more reliable identification of endothelial erosions.

Identification of two novel human neuronal ceroid lipofuscinosis genes

The neuronal ceroid lipofuscinoses (NCLs) are a group of mostly autosomal recessively inherited neurodegenerative disorders. The aim of this thesis was to characterize the molecular genetic bases of these, previously genetically undetermined, NCL forms. Congenital NCL is the most aggressive form of NCLs. Previously, a mutation in the cathepsin D (CTSD) gene was shown to cause congenital NCL in sheep. Based on the close resemblance of the phenotypes between congenital NCLs in sheep and human, CTSD was considered as a potential candidate gene in humans as well. When screened for mutations by sequencing, a homozygous nucleotide duplication creating a premature stop codon was identified in CTSD in one family with congenital NCL. While in vitro the overexpressed truncated mutant protein was stable although inactive, the absence of CTSD staining in brain tissue samples of patients indicated degradation of the mutant CTSD in vivo. A lack of CTSD staining was detected also in another, unrelated family with congenital NCL. These results imply that CTSD deficiency underlies congenital NCL. While initially Turkish vLINCL was considered a distinct genetic entity (CLN7), mutations in the CLN8 gene were later reported to account for the disease in a subset of Turkish patients with vLINCL. To further dissect the genetic basis of the disease, all known NCL genes were screened for homozygosity by haplotype analysis of microsatellite markers and/or sequenced in 13 mainly consanguineous, Turkish vLINCL families. Two novel, family-specific homozygous mutations were identified in the CLN6 gene. In the remaining families, all known NCL loci were excluded. To identify novel gene(s) underlying vLINCL, a genomewide single nucleotide polymorphism scan, homozygosity mapping, and positional candidate gene sequencing were performed in ten of these families. On chromosome 4q28.1-q28.2, a novel major facilitator superfamily domain containing 8 (MFSD8) gene with six family-specific homozygous mutations in vLINCL patients was identified. MFSD8 transcript was shown to be ubiquitously expressed with a complex pattern of alternative splicing. Our results suggest that MFSD8 is a novel lysosomal integral membrane protein which, as a member of the major facilitator superfamily, is predicted to function as a transporter. Identification of MFSD8 emphasizes the genetic heterogeneity of Turkish vLINCL. In families where no MFSD8 mutations were detected, additional NCL-causing genes remain to be identified. The identification of CTSD and MFSD8 increases the number of known human NCL-causing genes to eight, and is an important step towards the complete understanding of the genetic spectrum underlying NCLs. In addition, it is a starting point for dissecting the molecular mechanisms behind the associated NCLs and contributes to the challenging task of understanding the molecular pathology underlying the group of NCL disorders.

Development of regulatory T cells in the human thymus

The role of the immune system is to protect an organism against pathogens while maintaining tolerance against self. T cells are an essential component of the immune system and they develop in the thymus. The AIRE (autoimmune regulator) gene product plays an important role in T cell development, as it promotes expression of peripheral tissue antigens in the thymus. Developing T cells, thymocytes, which recognize self-antigens with high affinity are deleted. However, this deletion process is not perfect and not all autoreactive T cells are destroyed. When the distinction between self and non-self fails, tolerance breaks and the immune system attacks the host s own tissues. This results in autoimmunity. Regulatory T cells contribute to the maintenance of self-tolerance. They can actively suppress the function of autoreactive cells. Several populations of cells with regulatory properties have been described, but the best characterized population is the natural regulatory T cells (Treg cells), which develop in the thymus and express the transcription factor FOXP3. The thymic development of Treg cells in humans is the subject of this thesis. Thymocytes at different developmental stages were analyzed using flow cytometry. The CD4-CD8- double-negative (DN) thymocytes are the earliest T cell precursors in the T cell lineage. My results show that the Treg cell marker FOXP3 is up-regulated already in a subset of these DN thymocytes. FOXP3+ cells were also found among the more mature CD4+CD8+ double-positive (DP) cells and among the CD4+ and CD8+ single-positive (SP) thymocytes. The different developmental stages of the FOXP3+ thymocytes were isolated and their gene expression examined by quantitative PCR. T cell receptor (TCR) repertoire analysis was used to compare these different thymocyte populations. My data show that in humans commitment to the Treg cell lineage is an early event and suggest that the development of Treg cells follows a linear developmental pathway, FOXP3+ DN precursors evolving through the DP stage to become mature CD4+ Treg cells. Most T cells have only one kind of TCR on their cell surface, but a small fraction of cells expresses two different TCRs. My results show that the expression of two different TCRs is enriched among Treg cells. Furthermore, both receptors were capable of transmitting signals when bound by a ligand. By extrapolating flow cytometric data, it was estimated that the majority of peripheral blood Treg cells are indeed dual-specific. The high frequency of dual-specific cells among human Treg cells suggests that dual-specificity has a role in directing these cells to the Treg cell lineage. It is known that both genetic predisposition and environmental factors influence the development of autoimmunity. It is also known that the dysfunction or absence of Treg cells leads to the development of autoimmune manifestations. APECED (autoimmune polyendocrinopathy-candidiasis-ectodermal dystrophy) is a rare monogenic autoimmune disease, caused by mutations in the AIRE gene. In the absence of AIRE gene product, deletion of self-specific T cells is presumably disturbed and autoreactive T cells escape to the periphery. I examined whether Treg cells are also affected in APECED. I found that the frequency of FOXP3+ Treg cells and the level of FOXP3 expression were significantly lower in APECED patients than in controls. Additionally, when studied in cell cultures, the suppressive capacity of the patients' Treg cells was impaired. Additionally, repertoire analysis showed that the TCR repertoire of Treg cells was altered. These results suggest that AIRE contributes to the development of Treg cells in humans and the selection of Treg cells is impaired in APECED patients. In conclusion, my thesis elucidates the developmental pathway of Treg cells in humans. The differentiation of Tregs begins early during thymic development and both the cells dual-specificity and AIRE probably affect the final commitment of Treg cells.

Examination of Genetic Components Affecting Human Obesity-Related Quantitative Traits

Obesity increases the risk for several conditions, including type 2 diabetes mellitus, cardiovascular disease, hypertension, osteoarthirits and certain types of cancer. Twin- and family studies have shown that there is a major genetic component in the determination of body mass. In recent years several technological and scientific advance have been made in obesity research. For instance, novel replicated loci have been revealed by a number of genome wide association studies. This thesis aimed to investigate the association of genetic factors and obesity-related quantitative traits. The first study investigated the role of the lactase gene in anthropometric traits. We genetically defined lactose persistence by genotyping 31 720 individuals of European descent. We found that lactase persistence was significantly correlated with weight and body mass index but not with height. In the second study we performed the largest whole genome linkage scan for body mass index to date. The sample consisted of 4401 twin families and 10 535 individuals from six European countries. We found supporting evidence for two loci (3q29 and 7q36). We observed that the heritability estimate increased substantially when additional family members were removed from the analyses, which suggests reduced environmental variance in the twin sample. In the third study we assessed metabonomic, transcriptomic and genomic variation in a Finnish population cohort of 518 individuals. We formed gene expression networks to portray pathways and showed that a set of highly correlated genes of an inflammatory pathway associated with 80 serum metabolites (of 134 quantified measures). Strong association was found, for example, with several lipoprotein subclasses. We inferred causality by using genetic variation as anchors. The expression of the network genes was found to be dependent on the circulatory metabolite concentrations.