Latent physiological factors of complex human diseases revealed by independent component analysis of clinarrays
© Butte et al; licensee BioMed Central Ltd. 2010
Published: 28 October 2010
Diagnosis and treatment of patients in the clinical setting is often driven by known symptomatic factors that distinguish one particular condition from another. Treatment based on noticeable symptoms, however, is limited to the types of clinical biomarkers collected, and is prone to overlooking dysfunctions in physiological factors not easily evident to medical practitioners. We used a vector-based representation of patient clinical biomarkers, or clinarrays, to search for latent physiological factors that underlie human diseases directly from clinical laboratory data. Knowledge of these factors could be used to improve assessment of disease severity and help to refine strategies for diagnosis and monitoring disease progression.
Applying Independent Component Analysis on clinarrays built from patient laboratory measurements revealed both known and novel concomitant physiological factors for asthma, types 1 and 2 diabetes, cystic fibrosis, and Duchenne muscular dystrophy. Serum sodium was found to be the most significant factor for both type 1 and type 2 diabetes, and was also significant in asthma. TSH3, a measure of thyroid function, and blood urea nitrogen, indicative of kidney function, were factors unique to type 1 diabetes respective to type 2 diabetes. Platelet count was significant across all the diseases analyzed.
The results demonstrate that large-scale analyses of clinical biomarkers using unsupervised methods can offer novel insights into the pathophysiological basis of human disease, and suggest novel clinical utility of established laboratory measurements.
Factorial methods have been successfully applied to diverse forms of multidimensional biological and clinical data to uncover both biological and clinical phenomena. Principal Component Analysis (PCA) and Independent Component Analysis (ICA) have been applied to gene expression microarrays to uncover latent factors in the expression data related to cellular processes, transcriptional control programs, and disease subtypes [1–4]. ICA has also been successfully applied to fMRI brain imaging data to uncover nondeterministic signals of interest .
ICA is a powerful extension of PCA originally developed as a solution to signal separation problems. PCA can only impose statistical independence of components up to the second order; therefore it can only identify directions that are uncorrelated and orthogonal to each other. ICA is capable of exploiting higher-order statistics to relax the orthogonality assumption and identify components that are mutually statistically independent from each other, which is a stronger condition than lack of correlation . Given that ICA aims to find components that are non-Gaussian, it has the desirable side effect of ignoring variance resulting from noise in the data. It is likely that this higher-order model is more reflective of biological phenomenon, and thus offers an explanation for the many successful applications of ICA in the biomedical domain.
In previous studies, we demonstrated the utility for patient clinical biomarker vectors, which we termed clinarrays, in distinguishing between different severities of disease ; in integrative analysis with gene expression measurements to elucidate genes related to maturation and aging ; and to build a model of pediatric aging . Clinarrays are a vector representation of a person’s physiological biomarker state across all visits at a hospital. Each value in this vector represents a summary statistic of a particular biomarker across a period of time. Aggregations of clinarrays enable the use of methods developed for microarray analysis to gain new insights into the clinical characteristics of disease-affected patients.
We demonstrate that ICA can be applied to clinarrays of patient lab measurements to uncover physiological factors known to be concomitant in common chronic diseases. Furthermore, we identify a novel latent physiological factor for cystic fibrosis that serves as a promising lead for further investigation into diagnostic or prognostic biomarkers for CF.
Creation of disease-specific clinarrays
Number of patients and biomarkers remaining after pruning
Type 1 diabetes
Type 2 diabetes
Duchenne muscular dystrophy
Independent component analysis of clinarrays
Significant biomarkers after ICA analysis
Type 1 Diabetes
Blood Urea Nitrogen
Type 2 Diabetes
Duchenne Muscular Dystrophy
IgE antibody response
We applied ICA analysis to biomarker clinarrays in an effort to elucidate latent physiological factors concomitant in five chronic disorders. Overall we find associations that are consistent with known physiological factors of the diseases studied, showing that the factors extracted by the ICA algorithm are reflective of actual physiological processes underlying each disease. It is interesting to note that platelet count, which we use as a marker for thrombogenesis, appears to be the only factor shared across all of the diseases studied. Recent studies have revealed physiological roles for platelets that extend well beyond thrombogenesis, including roles in innate immunity, microbial defense, and angiogenesis. Our findings suggest that platelet activation may be a salient pathophysiological factor for a broad range of chronic disorders.
We investigated biomarkers of patients with Cystic Fibrosis (CF), a hereditary multi-organ inflammatory disorder caused by mutations in the CFTR gene. The results identify three significant physiological factors latent in the disease pathology: IgE antibody response (Total IgE), dephosphorylation (alkaline phosphatase) and thrombogenesis (platelet count). Platelet activation has been implicated as a major factor in the pathogenesis and progression of CF. Falco and colleagues suggest that platelets participate in the pathogenesis of CF by increasing levels of soluble CD40 ligand, whose secretion has been linked to coagulation activity . These findings elucidate a putative role for platelet activation in disease progression, and inform the course of treatment for patients exhibiting high degrees of platelet activation.
A clear pathophysiological explanation for alkaline phosphatase (ALKP) is not immediately evident. Differential levels of ALKP in CF may be attributed to bone disease or biliary obstruction. The interpretation is confounded by the fact that the data was derived from pediatric patients who may have greater variance of ALKP, due to differential rates of bone growth. However, no other disease resulted in significance for ALKP, even though the range of ages was similar for the other diseases. ALKP remains a potentially novel finding.
One of the more interesting findings from a clinical perspective is the identification of an Immunoglobulin E (IgE) component in CF. IgE plays an important role in allergic sensitivity and response . IgE levels are known to be elevated in a subset of CF patients due to allergic bronchopulmonary aspergillosis (ABPA), an allergic reaction to a secondary fungal infection of an aspergillius species. The prevalence ABPA in pediatric CF patients is estimated to be anywhere from 2% to 15%. Despite the fact that detection and treatment of fungus has been known to improve prognosis, it is not common practice to screen for fungus. Testing for ABPA in CF patients is not usually conducted until the patient fails to respond to antibiotic therapy for an extended period. Our findings suggest that ABPA may be more of a significant factor among patients seeking treatment for CF symptom which suggests a more routine screening for IgE antibodies to mitigate prolonged lung damage and improve patient prognosis.
Type 1 and type 2 diabetes, though similar in regards to causing an abnormal increase in blood sugar, have etiological roots in distinctly different pathophysiological mechanisms. Both types of diabetes are known to cause hyperglycemia, which decreases serum sodium concentrations due to water exiting the intracellular to the extracellular space . Hyperglycemia is known to be associated with decreased measured serum sodium . Therefore the identification of serum sodium as a significant factor for both forms of diabetes by our analysis is concordant with expectations from established pathophysiology. Sterner and colleagues showed that an increased platelet count is associated with type 1 diabetes both with and before renal impairment . Jesri and colleagues showed that platelet counts significantly increased as patients accrued more risk factors for metabolic syndrome , a precursor for type 2 diabetes. Therefore while platelet count was significant in all diseases, there is an established basis for its association with both types of diabetes analyzed.
Thyroid stimulating hormone (TSH) was found to be significant in distinguishing patients with type 1 diabetes, while it was conspicuously absent from our analysis of type 2 diabetes. Thyroid autoimmunity affects approximately 20% of patients diagnosed with type 1 diabetes . For patients with thyroid autoimmunity TSH levels were found to be significantly higher. Moreover, patients with two thyroid antibodies had an even higher TSH level . As expected, we did not find this biomarker in type 2 diabetes. Similarly, we recovered a serum alanine aminotransferase (ALT) association with type 2 diabetes. It has been previously shown that elevated ALT was associated with insulin resistance, the main component of the pathogenesis of type 2 diabetes . We also note that certain biomarkers with known associations to type 1 diabetes and type 2 diabetes, such as hemoglobin A1c, are not returned in our analysis. This may be due to the inability of these biomarkers to differentiate between patient subsets within these particular diseases, or a lack of statistical independence in the associated data.
We also examined Duchenne muscular dystrophy (DMD), a disease whose pathogenesis stems from a mutation in the dystrophin gene. It was previously shown that patients with DMD may have elevated levels of AST, which is corroborated by our analysis . These high levels of AST have been attributed to muscle breakdown rather than liver pathology, which is traditionally associated with AST levels. As more muscle begins to break down in the progression of DMD, we would expect to see an increase in levels markers for tissue damage, such as AST. Therefore we suggest that AST may serve as a marker for disease progression in DMD. It also has been shown that lactate dehydrogenase (LDH), also identified by our analysis, is increased in patients affected with DMD .
A relationship whose basis is less evident is the implication of triglyceride levels as a factor in DMD. Young and colleagues have shown that mice with hereditary muscular dystrophy have significantly increased amounts of triglyceride in skeletal muscle biopsies . There has been a conspicuous absence of examining triglyceride levels in the blood for DMD patients, which is captured in clinical lipid profiles. A recent study by Wren and colleagues showed that muscle adiposity values are accurate in determining disease severity of DMD patients . Our analysis suggests that triglyceride levels in blood could also serve as a novel biomarker that indicates the severity of DMD.
We acknowledge some limitations in our approach. Foremost, the range of lab measurements found in the patient data limited the identification of significant biomarkers for each disease. These biomarkers were also manually mapped to general physiological processes and thus are limited to our current understanding of physiology. Biomarkers may be involved in multiple physiological processes and may be part of processes we do not capture. We also acknowledge that clinarrays do not take into consideration the temporal aspect of disease progression. However, as with microarrays, the temporal aspect of gene expression within a given tissue is also often ignored unless specifically part of the protocol. ICA-based clinarray analysis could be extended to incorporate a temporal component in future work. We also note that while we currently do not use the longitudinal data available at a patient-specific level in electronic medical records, this type of data could be incorporated in future work to potentially increase the power of our analysis. In several cases, lab measurements could not be included in the clinarray for a disease as the underlying patient data was too sparse. Thus, the set of laboratory measurements was not the same for each clinarray. In addition, the use of ICA is only appropriate when the assumption of combinatorial linearity of components is satisfied. Further testing is needed to determine if this assumption holds generally for clinical biomarker data.
In future work we hope to expand our data set to include biomarkers for additional diseases and formulate clinical validation of the results as well as integrate longitudinal and temporal information to produce more extensive models of disease pathophysiology.
We have developed a novel approach that incorporates independent component analysis of patient laboratory biomarkers represented as vectors in a clinarray. This approach seeks to separate out statistically independent signals of disease pathophysiology as measured by clinical laboratory tests on patient populations. We applied this approach to a data set characterizing five distinct disease conditions in a pediatric patient population. The results of this analysis yielded both unique and common physiological factors associated with the diseases analyzed. Several of these associations can be explained by known disease pathophysiology, whereas some of the associations are novel, and therefore suggest novel disease biology as well as novel uses for established laboratory biomarkers in clinical management of these diseases. This work provides a novel basis for unsupervised, data-driven analysis of disease pathophysiology from aggregate clinical laboratory measurements of patient populations.
Building clinarrays from patient lab tests
Quantitative clinical biomarker data, consisting of 893,956 measurements across 972 distinct biomarkers, obtained at the Lucile Packard Children’s Hospital (LPCH), were collected in a de-identified manner from the Stanford Translational Research Database Environment (STRIDE; http://stride.stanford.edu). This data represented 4085 patients diagnosed with one or more of five chronic diseases: asthma, type 1 diabetes, type 2 diabetes, cystic fibrosis, and Duchenne muscular dystrophy. The Institutional Review Board of the Stanford University School of Medicine approved the use of de-identified clinical data in this manner.
Clinarrays were generated from patients diagnosed with one or more of the diseases previously mentioned. The median value for each biomarker across all measurements of that biomarker for an individual patient was calculated. The median values for all biomarkers were aggregated to create the clinarray. Due to the paucity of data, thresholds were implemented to exclude patients and biomarkers that were poorly represented. Patients with less than 10 distinct biomarkers were excluded. Biomarkers without measurements for at least fifty percent of the remaining patients were also excluded.
Independent component analysis
For our analysis we consider a clinarray data matrix X where the columns correspond to individual patient clinarrays, one clinarray per patient, and whose rows correspond to the variable biomarker measurements for each patient. Missing values were imputed using the K nearest neighbor imputation (KNNImpute) algorithm with K = 10. The biomarker values were centered to a mean of zero and scaled to unit variance. The ICA model considers X as the matrix product of statistically independent components or signals S, and a mixing matrix A.
X = AS
Thus the observed physiological measures in the clinarray X are considered to be a linear combination of statistically independent latent physiological modes in various combinations as defined by the unknown mixing matrix A. Thus the goal of ICA is to estimate an optimal un-mixing matrix W that attempts to satisfy the equation:
Y = WX = WAS
Where Y is an approximation of S and if W=A -1 then Y is a perfect reconstruction of S.
In this analysis we employed the FastICA algorithm proposed by Hyvärinen. FastICA employs a contrast function that uses approximations of neg-entropy to identify non-Gaussian components. We determined the number of components to extract using the screen test method, which is used in exploratory factor analysis to estimate the putative number of latent principal factors.
Extracting physiological factors
The results produced by the FastICA algorithm are influenced by random initializations from which the objective function of minimizing statistical independence of components is optimized. Thus the algorithm is known to find local minima depending on the initial conditions, and successive runs of the algorithm will produce slightly different results. In addition, ICA algorithms do not extract the components in order. Therefore we designed an iterative ICA analysis in which the FastICA algorithm was applied to the clinarray data over 500 iterations.
We investigated two different techniques for extracting significant laboratory biomarkers from the ICA iterations. First we examined the distribution of the absolute signal loading scores across all components and iterations. We then identified significant lab tests by identifying, for each disease, lab tests whose mean absolute component loading score was in the top 1% of the distribution. In a second approach, we counted number of times each lab was found to have the highest absolute loading score within a component across all iterations, and selected lab tests found to load highest on a component in all iterations. Both approaches yielded exactly the same results.
The work was supported by the Lucile Packard Foundation for Children's Health, the William R. Hewlett Stanford Graduate Fellowship, NIBIB and NLM Grant (T15 LM007033), National Library of Medicine (R01 LM009719), National Institute of General Medical Sciences (R01 GM079719), and the Howard Hughes Medical Institute. We thank Alex Skrenchuk and Boris Oskotsky from Stanford University for computer support.
This article has been published as part of BMC Bioinformatics Volume 11 Supplement 9, 2010: Selected Proceedings of the 2010 AMIA Summit on Translational Bioinformatics. The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2105/11?issue=S9.
- Raychaudhuri S, Stuart JM, Altman RB: Principal components analysis to summarize microarray experiments: application to sporulation time series. Pacific Symposium on Biocomputing 2000, 455–466.Google Scholar
- Liebermeister W: Linear modes of gene expression determined by independent component analysis. Bioinformatics 2002, 18(1):51–60. 10.1093/bioinformatics/18.1.51View ArticlePubMedGoogle Scholar
- Frigyesi A, Veerla S, Lindgren D, Höglund M: Independent component analysis reveals new and biologically significant structures in micro array data. BMC Bioinformatics 2006, 7: 290. 10.1186/1471-2105-7-290PubMed CentralView ArticlePubMedGoogle Scholar
- Saidi SA, Holland CM, Kreil DP, MacKay DJ, Charnock-Jones DS, Print CG, Smith SK: Independent component analysis of microarray data in the study of endometrial cancer. Oncogene 2004, 23(39):6677–6683. 10.1038/sj.onc.1207562View ArticlePubMedGoogle Scholar
- Kiviniemi V, Kantola JH, Jauhiainen J, Hyvarinen A, Tervonen O: Independent component analysis of nondeterministic fMRI signal sources. Neuroimage 2003, 19(2 Pt 1):253–260. 10.1016/S1053-8119(03)00097-1View ArticlePubMedGoogle Scholar
- Comon P: Independent Component Analysis, a New Concept. Signal Processing 1994, 36(3):287–314. 10.1016/0165-1684(94)90029-9View ArticleGoogle Scholar
- Chen DP, Weber SC, Constantinou PS, Ferris TA, Lowe HJ, Butte AJ: Clinical Arrays of Laboratory Measures, or “Clinarrays”, Built from an Electronic Health Record Enable Disease Subtyping by Severity. 2007: American Medical Informatics Association 2007, 115.Google Scholar
- Chen DP, WS, Constantinou PS, Ferris TA, Lowe HJ, Butte AJ: Novel Integration of hospital electronic medical records and gene expression measurements to identify genetic markers of maturation. Pacific Symposium on Biocomputing 2008, 243–254.Google Scholar
- Chen DP, Morgan AA, Butte AJ: Validating pathophysiological models of aging using clinical electronic medical records. Journal of biomedical informatics 2009.Google Scholar
- Fischbach FT, Dunning MB: A manual of laboratory and diagnostic tests. 7th edition. Philadelphia: Williams & Wilkins; 2004.Google Scholar
- Leslie M: Cell biology. Beyond clotting: the powers of platelets. Science 2010, 328(5978):562–564. 10.1126/science.328.5978.562View ArticlePubMedGoogle Scholar
- Falco A, Romano M, Iapichino L, Collura M, Davi G: Increased soluble CD40 ligand levels in cystic fibrosis. Journal of Thrombosis and Haemostasis 2004, 2(4):557–560. 10.1111/j.1538-7836.2004.00683.xView ArticlePubMedGoogle Scholar
- O'Sullivan BP, Linden MD, Frelinger AL, Barnard MR, Spencer-Manzon M, Morris JE, Salem RO, Laposata M, Michelson AD: Platelet activation in cystic fibrosis. Blood 2005, 105(12):4635–4641. 10.1182/blood-2004-06-2098View ArticlePubMedGoogle Scholar
- Haworth CS, Selby PL, Webb AK, Dodd ME, Musson H, Niven RM, Economou G, Horrocks AW, Freemont AJ, Mawer EB, et al.: Low bone mineral density in adults with cystic fibrosis. Thorax 1999, 54(11):961–967. 10.1136/thx.54.11.961PubMed CentralView ArticlePubMedGoogle Scholar
- Aris RM, Merkel PA, Bachrach LK, Borowitz DS, Boyle MP, Elkin SL, Guise TA, Hardin DS, Haworth CS, Holick MF, et al.: Consensus statement: Guide to bone health and disease in cystic fibrosis. Journal of Clinical Endocrinology and Metabolism 2005, 90(3):1888–1896. 10.1210/jc.2004-1629View ArticlePubMedGoogle Scholar
- Gould HJ, Sutton BJ, Beavil AJ, Beavil RL, McCloskey N, Coker HA, Fear D, Smurthwaite L: The biology of IgE and the basis of allergic disease. Annual Review of Immunology 2003, 21: 579–628. 10.1146/annurev.immunol.21.120601.141103View ArticlePubMedGoogle Scholar
- Tillie-Leblond I, Tonnel AB: Allergic bronchopulmonary aspergillosis. Allergy 2005, 60(8):1004–1013. 10.1111/j.1398-9995.2005.00887.xView ArticlePubMedGoogle Scholar
- Mastella G, Rainisio M, Harms HK, Hodson ME, Koch C, Navarro J, Strandvik B, McKenzie SG: Allergic bronchopulmonary aspergillosis in cystic fibrosis. European Respiratory Journal 2001, 17(5):1052–1053. 10.1183/09031936.01.17510520View ArticlePubMedGoogle Scholar
- Seldin DW, Tarail R: Effect of Hypertonic Solutions on Metabolism and Excretion of Electrolytes. American Journal of Physiology 1949, 159(1):160–174.PubMedGoogle Scholar
- Hillier TA, Abbott RD, Barrett EJ: Hyponatremia: Evaluating the correction factor for hyperglycemia. American Journal of Medicine 1999, 106(4):399–403. 10.1016/S0002-9343(99)00055-8View ArticlePubMedGoogle Scholar
- Sterner G, Carlson J, Ekberg G: Raised platelet levels in diabetes mellitus complicated with nephropathy. Journal of Internal Medicine 1998, 244(6):437–441. 10.1111/j.1365-2796.1998.00349.xView ArticlePubMedGoogle Scholar
- Jesri A, Okonofua EC, Egan BM: Platelet and white blood cell counts are elevated in patients with the metabolic syndrome. In Journal of clinical hypertension. Volume 7. Greenwich, Conn; 2005:705–711. quiz 712–703 quiz 712-703 10.1111/j.1524-6175.2005.04809.xGoogle Scholar
- Leshin M: Polyglandular Autoimmune Syndromes. American Journal of the Medical Sciences 1985, 290(2):77–88. 10.1097/00000441-198508000-00007View ArticlePubMedGoogle Scholar
- Kordonouri O, Klinghammer A, Lang EB, Grabert M, Holl RW: Thyroid Autoimmunity in children and adolescents with type 1 diabetes mellitus: A multi-center survey. Diabetes 2002, 51: A52-A52.Google Scholar
- Hanley AJ, Wagenknecht LE, Festa A, D'Agostino RB Jr., Haffner SM: Alanine aminotransferase and directly measured insulin sensitivity in a multiethnic cohort: the Insulin Resistance Atherosclerosis Study. Diabetes care 2007, 30(7):1819–1827. 10.2337/dc07-0086View ArticlePubMedGoogle Scholar
- EH R: Muscular Dystrophy. Pediatr Rev 2000, 21(7):233–237. 10.1542/pir.21-7-233View ArticleGoogle Scholar
- Perloff JK, Deleon AC, Odoherty D: Cardiomyopathy of Progressive Muscular Dystrophy. Circulation 1966, 33(4):625-&.View ArticlePubMedGoogle Scholar
- Young HL, Young W, Edelman IS: Electrolyte and lipid composition of skeletal and cardiac muscle in mice with hereditary muscular dystrophy. The American journal of physiology 1959, 197: 487–490.PubMedGoogle Scholar
- Wren TA, Bluml S, Tseng-Ong L, Gilsanz V: Three-point technique of fat quantification of muscle tissue as a marker of disease progression in Duchenne muscular dystrophy: preliminary study. Ajr 2008, 190(1):W8–12. 10.2214/AJR.07.2732View ArticlePubMedGoogle Scholar
- Troyanskaya O, Cantor M, Sherlock G, Brown P, Hastie T, Tibshirani R, Botstein D, Altman RB: Missing value estimation methods for DNA microarrays. Bioinformatics 2001, 17(6):520–525. 10.1093/bioinformatics/17.6.520View ArticlePubMedGoogle Scholar
- Hyvarinen A: Fast and robust fixed-point algorithms for independent component analysis. Ieee Transactions on Neural Networks 1999, 10(3):626–634. 10.1109/72.761722View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.