Genetic influences on peripheral blood cell counts: a study in baboons
http://www.100md.com
《血液学杂志》
the Department of Genetics, Southwest Foundation for Biomedical Research
Southwest National Primate Research Center, San Antonio, TX
the Department of Laboratory Medicine, Children's Hospital, Boston, MA.
Abstract
Interperson differences in peripheral blood cell counts in healthy individuals result from genetic and environmental influences. We used multivariate genetic analyses to assess the relative impact of genes and environment on baseline blood cell counts and indices using a pedigreed colony of baboons, an animal with well-documented analogies to human blood physiology. After accounting for age, sex, and weight, we found that genetic influences explain a significant proportion of the remaining variability, ranging from a low of 13.7% for mean corpuscular hemoglobin concentration (MCHC) to a high of 72.4% for red blood cell (RBC) number. Genes influence 38.5% of the variation in baseline white blood cell (WBC) count, a characteristic that correlates with mortality in both the general human population and clinically defined subgroups such as individuals with sickle-cell disease. We examined the interaction between pairs of traits and identified those that share common genetic influences (pleiotropy). We unexpectedly observed that the same gene or group of genes influences both WBC count and mean platelet volume (MPV). We anticipate that this approach will ultimately lead to discovery of novel insights into the biology of related traits, and ultimately identify important genes that affect hematopoiesis.
Introduction
Peripheral blood cell counts represent a complex clinical phenotype—the end result of interactions between a multitude of genetic and environmental factors. Increasingly, baseline blood cell count phenotypes are emerging as important markers of risk for disease and disease severity. examples include the correlation of steady-state white blood cell (WBC) count with mortality in the general population,1-6 and with disease severity in patients with sickle-cell anemia (SS).7-10 efforts to dissect the genetic underpinnings of these complex traits have the potential to greatly enhance our understanding of hematopoiesis, pathophysiology, and disease modifier genes.
The primary genetic questions about quantitative traits such as peripheral blood cell counts (ie, traits measured on continuous rather than dichotomous scales) are formulated in terms of variation in the traits.11 The amount of variation is expressed as the variance, or in statistical terms, the mean of the squared deviations from the population mean for the trait. The fundamental objective of statistical genetic analyses of any quantitative trait is to partition the total variance for that trait (the phenotype) into 2 components: those due to the effects of genes (the genetic component), and those influenced by environmental factors (the environmental component). In this study, we examine the patterns of normal quantitative variation in blood cell counts in a pedigreed colony of baboons, a nonhuman primate with well-documented analogies to human blood physiology,12-15 in order to partition the blood cell phenotypes into their genetic and environmental components. Our specific aims are to quantify the genetic and environmental components, and determine which genetic effects on one cell count may also account for variation in other cell counts (pleiotropy).
Materials and methods
Subjects
We analyzed the blood counts and indices of 600 individual healthy pedigreed baboons (Papio hamadryas)16,17 to quantify the range of values observed in this normal population. Housed at the Southwest National Primate Research Center, located at the Southwest Foundation for Biomedical Research (San Antonio, TX), this sample consisted mainly of olive baboons (Papio hamadryas anubis), yellow baboons (Papio hamadryas cynocephalus), and their hybrids, and initially included 415 females and 185 males, ranging in age from 5.5 to 30 years. While strict genetic management is used to control inbreeding, all nonfounder animals in this study were the results of matings that were random with respect to phenotypes analyzed here. The subjects were distributed among 11 pedigrees with 67 to 171 members each. Within these pedigrees there were 4714 relative pairs (Table 1).
Since birth or arrival at the colony, all animals had been housed out of doors in social group cages and maintained on monkey chow diets to which they had ad libitum access. Animal care personnel and staff veterinarians provided daily maintenance and regular, urgent, and emergent health care to all animals in accordance with the Guide for the Care and Use of Laboratory Animals.18 All procedures related to their treatment during the conduct of this study were approved by the Institutional Animal Care and Use Committee in accordance with the established guidelines.18
Hematologic studies
After sedation with intramuscular ketamine and intravenous xylazine hydrochloride, atropine, acepromazine, and ketamine to assure relaxation, blood samples were obtained by femoral venipuncture. Blood was collected in eDTA (ethylenediaminetetraacetic acid), stored on ice, and analyzed within 4 hours using a Coulter JT (Beckman Coulter, Fullerton, CA), conditions that maintained the stability of the measured parameters. We measured hematocrit (HCT), hemoglobin (Hb), red blood cell (RBC) count, mean corpuscular volume (MCV), red cell distribution width (RDW), mean corpuscular hemoglobin (MCH), mean corpuscular hemoglobin concentration (MCHC), white blood cell (WBC) count, platelet (Plt) count, and mean platelet volume (MPV). Details regarding protocols for these assays, and measurements are reported elsewhere.19
Statistical analysis
Statistical genetic analyses were conducted by means of a maximum likelihood–based variance decomposition approach implemented in the computer package SOLAR (Sequential Oligogenic Linkage Analysis Routines; Southwest Foundation for Biomedical Research, San Antonio, TX).20 This approach, which was developed following methodology originally proposed by Hopper and Mathews21 and Lange and Boehnke,22 has been described in detail elsewhere.23
In short, we used this approach to partition the phenotypic variance (2P) into components corresponding to the additive genetic (2G) and nongenetic (ie, environmental) (2e) effects. Because these components are additive, such that 2P =2G +2e, we estimated the heritability, or proportion of the phenotypic variance attributable to additive genetic effects, as h2 =2G /2P. We estimated the proportion of the phenotypic variance attributable to nongenetic factors as e2 = 1 – h2. When covariate effects are included in the maximum likelihood model, the h2 and e2 estimates are obtained as proportions of the residual phenotypic variance in the trait (ie, that remaining after the accounting for the effects of the covariates) rather than as proportions the total phenotypic variance. Multiplication of h2 by the proportion of the total variance not explained by covariates (ie, the quantity 1 – c2, where "c2" refers to the proportion of the total variance attributable to covariate effects) yields a heritability value that is rescaled to the total phenotypic variance. We also simultaneously estimated the mean effects of sex, sex-specific age, and sex-specific age,2 and body weight on each of the hematological phenotypes. We model the phenotype of an individual (yi) as linear function of the population mean of the trait (μ), the covariates (x) and their regression coefficients (), plus the additive genetic values (g) and random environmental deviations (e), such that yi = μ + 1 x1+ 2 x2 + nxn + g + e.
Using extensions to univariate genetic analysis that encompass the multivariate state,21,22,24 we followed an approach presented in detail previously25 to similarly model multivariate phenotypes in these baboons (ie, to simultaneously analyze measurements on pairs of traits in individuals within each pedigree). By maximizing these multivariate models, we obtained estimates of 2 additional parameters: the additive genetic correlation, G, and the environmental correlation, e, between trait pairs. Respectively, these correlations are estimates of the shared additive effects of genes (ie, pleiotropy) and shared environmental (ie, unmeasured, nongenetic) factors on the variance in a trait. Genetic correlations with values equal to (or not significantly different from) 1 or –1 are interpreted as evidence for complete pleiotropy (ie, where the additive genetic contributions to the variance in the 2 correlated traits are due entirely to the effects of the same gene or genes). Genetic correlations that are significantly different from both 0 and 1 or – 1 are evidence of incomplete pleiotropy (ie, where some but not all of the additive genetic effects on a pair of traits are due to the same gene or genes). The genetic and environmental components of the phenotypic correlation matrix are additive, like those of the corresponding variance-covariance matrix, so we could use the maximum likelihood estimates of the additive genetic and environmental correlations, to obtain the total phenotypic correlation between 2 traits, P, as .
Using this approach as implemented in SOLAR,20 we conducted a series of bivariate quantitative genetic analysis of all hematologic trait pairs.
Significance of the maximum likelihood estimates for heritability and other parameters was assessed by means of likelihood ratio tests.26 The maximum likelihood for the general model in which all parameters were estimated was compared to that for restricted models in which the value of the parameter to be tested was held constant at some value (usually 0). For tests of parameters like h2, for which a fixed value of 0 in a restricted model is at a boundary of the parameter space, twice the difference in the natural logarithm (ln) of the likelihoods of the 2 models compared is distributed asymptotically approximately as a 0.5:0.5 mixture of 2 with a point mass at 0. For tests of covariates, for which 0 is not a boundary value, this statistic is distributed asymptotically approximately as a 2 variate.21 Degrees of freedom can be obtained as the difference in the number of estimated parameters in the 2 models.21 However, in tests of parameters like h2, whose values may be fixed at a boundary of their parameter space in the null model, the appropriate significance level is obtained by halving the P value.24
Results
Inheritance of the individual peripheral blood cell counts
In these healthy pedigreed baboons, who live in a highly controlled environment designed to maximize the detection of genetic signals, a significant fraction of the phenotypic variance of all of the peripheral blood cell counts and indices is under genetic control (Table 2). After accounting for the individual and combined effects of age, sex, and weight, the additive effects of genes (h2) account for a low of 13.7% of the observed variability in MCHC, to a high of 72.4% of the variability of RBC counts.
Phenotypic correlations between the individual cell traits
Many of the cell counts and indices are correlated with others. As shown in Table 3, the erythrocyte traits are highly correlated, some positively (eg, Hct and Hb), and some negatively (eg, RBC counts and MCV). Platelet count and MPV are negatively correlated. WBC count and MPV are positively correlated; WBC count and Plt are negatively correlated.
evidence that shared genes influence more than one trait
We had not anticipated finding significant genetic influences shared between WBC count and MPV (P < .05). In fact, similar to what we found in the case of Hct and Hb, the genetic effects shared between WBC count and MPV is virtually 100% (ie, complete pleiotropy (G = 1).
Discussion
Most of what we know about the genetics of the peripheral blood counts and indices comes from the study of patients and experimental animals whose counts lie outside the normal range. Humans and animals with too many or too few cells that are too big or too small have long been the focus of attention of hematologists, and have established the pathophysiologic importance of a variety of genes that control phenotype in a Mendelian fashion. Some of these genes include those that encode or control globins, membrane proteins, growth factors, receptors, ion channels, and adhesion molecules. In much the same way, clinical investigators have discovered critical environmental (ie, noninherited) modulators of blood counts such as normal aging, maturation and growth, pregnancy, diet (iron, folate, cobalamin), hypoxia, bleeding, smoking, infection, inflammation, and exogenous growth factors and cytotoxic agents.
In this study, we shift focus from the classical clinical study of hematologically abnormal individuals to a statistical genetic methodology that analyzes the normal variation of cell counts in a healthy pedigreed population (ie, a population of individuals that were not ascertained on the basis of a disease state). Conceptually, this method is very similar to studies of healthy human twins. In twin studies, the genetic influence is estimated by comparing degrees of similarity between 2 types of relative pairs: monozygotic and dizygotic twins. Studies of more complex, multigenerational pedigrees such as ours make use of information from a much larger range of relative pairs comprising parents and offspring, grandparents and grandchildren, siblings, cousins, aunts/uncles, and nieces/nephews (the "avuncular" relationship), etc. Although not unique to nonhuman primates, the baboon pedigrees from which data for this study were obtained also include large numbers of pairs exhibiting half relationships (eg, half-sibling and half-avuncular pairs). Because it uses information from more relative classes, this strategy offers demonstrably greater statistical power to both detect genetic effects and, when combined with genotype data, efficiently localize quantitative trait loci to specific chromosomal regions.27,28 This approach considers the peripheral blood count a complex phenotype that is controlled by multiple underlying genetic and environmental factors, and is designed to quantify the impact and relatedness of those factors. We studied the baboon because of its close analogies to the human, and the ability to minimize the number and impact of environmental factors on interindividual variation.
We found that in the baboon, the steady-state set points for the number and size of peripheral blood cells are under significant genetic control (Table 1). These results are extremely similar to findings in humans by investigators working with major twin registries in the United Kingdom and Australia,29-31 and underscore the advantages of using such populations to discover genetic influences.
Our analyses further allowed us to examine the patterns of interrelatedness among these traits and, more specifically, to dissect out the genetic and random environmental components underlying these patterns. The positive phenotypic correlations between the various Hb, RBC count, and Hct pairings, and negative correlation of Hb, Hct, and RBC count with MCV are generally similar to the patterns of correlations observed between these traits in healthy human30-32 and other mammalian32 populations, although they have not been examined in detail. Many of these pairs demonstrated significant genetic contribution to the relationships—consistent with the findings that emerged from the Australian twin registry study.30 A similar inverse phenotypic relationship between Plt and MPV reported by Whitfield and Martin30 and others33-40 was also found in this baboon study, but other studies did not decompose the phenotypic covariance between platelet count and MPV into genetic and environmental components.
The positive phenotypic correlation that we found between WBC count and MPV has not been previously reported; and the observation of complete pleiotropy (ie, that all or most of the genetic effects on WBC count [responsible for 37% of the total phenotypic variance] and all of the genetic effects on MPV [49% of the total phenotypic variance] are due to the effects of the same gene or group of genes) is remarkable. Conceptually, the MPV resembles the MCV. Like the reticulocyte, the newly formed platelet is large, contains mRNA,41 and as an analogy—the increased number of reticulocytes in immune hemolytic anemia causes an elevated MCV, just as the high number of young platelets in idiopathic thrombocytopenic purpura elevates the MPV.42
Genetic abnormalities may make red cells small (thalassemia) or large (stomatocytosis), but we do not know what genes are responsible for that portion of variation in MCV (or MPV) that is due to the effects of genes in the healthy population. It may well be that these genes have little or no overlap with the globin genes that so dramatically influence cell size in abnormal populations.
It is interesting to note that high MPV43-45 and high WBC count1-6 both correlate with heart disease and stroke risk in the healthy population. While these may represent 2 genetically independent risk markers, it is tempting to speculate that genes exerting pleiotropic effects on both traits underlie the risk. Perhaps the genetic correlation between WBC count and platelet volume in relatively healthy individuals is attributable to the effects of genes that influence steps in hematopoiesis common to megakaryocytes and white cells, or genes that allow more mature platelets and white cells to respond to or enhance inflammation. It is important to point out that we do not know how much of the heart disease/stroke risk is influenced by the genetic component of the variation in WBC count and MPV, and how much by environmental factors. Similarly, we know that patients with sickle-cell anemia and high severity have high WBC counts, but do not know whether the high WBC count is an "environmental" (ie, nongenetic) reaction to disease severity, or whether the patients with a genetic predisposition to having a high WBC count are more likely to have severe disease as an expression of the role white blood cells play in disease pathophysiology.
Multivariate statistical genetic analyses of the type reported here offer a means to identify pairs and, in some cases, networks of traits that are bound together by the shared effects of the same gene or group of genes. In many instances, the detected pleiotropic relationships will be entirely consonant with our current expectations regarding the phenotypes under study (eg, RBC count and hemoglobin-related measures). In other cases however, previously unappreciated relationships, such as that between WBC count and MPV in this study, may be revealed or suggested by these sorts of analyses. The pleiotropic relationships detected in this study suggest a rationale for subsequent multivariate statistical analyses in both this pedigreed animal population, and in human families. Such analyses not only offer potentially novel insights into the biology of related traits but, in comparison to analyses of single traits, they also provide substantially greater amounts of information and improved statistical power to detect, localize, and identify quantitative trait loci contributing to the pleiotropic effects. Our plans to combine the data presented in this paper with genotypic data we are developing, will allow us to map the quantitative trait loci underlying patterns of covariation among these phenotypes and identify positional candidate genes in this nonhuman primate model for human hematological variation.
Footnotes
Prepublished online as Blood First edition Paper, May 3, 2005; DOI 10.1182/blood-2004-12-4863.
Supported by grants from the National Institutes of Health (NIH) grants R01 HL068922 (O.S.P.), R01 HL054141 (M.C.M), P01 HL028972 (M.C.M.), and P51 RR013086 (M.C.M.).
The publication costs of this article were defrayed in part by page charge payment. Therefore, and solely to indicate this fact, this article is hereby marked "advertisement" in accordance with 18 U.S.C. section 1734.
References
de Labry LO, Campion eW, Glynn RJ, Vokonas PS. White blood cell count as a predictor of mortality: results over 18 years from the Normative Aging Study. J Clin epidemiol. 1990;43: 153-157.
James AL, Knuiman MW, Divitini ML, Musk AW, Ryan G, Bartholomew HC. Associations between white blood cell count, lung function, respiratory illness and mortality: the Busselton Health Study. eur Respir J. 1999;13: 1115-1119.
Burr ML, Holliday RM, Fehily AM, Whitehead PJ. Haematological prognostic indices after myocardial infarction: evidence from the diet and reinfarction trial (DART). eur Heart J. 1992;13: 166-170.
Gillum RF, Ingram DD, Makuc DM. White blood cell count, coronary heart disease, and death: the NHANeS I epidemiologic Follow-up Study. Am Heart J. 1993;125: 855-863.
Grimm RH, Jr., Neaton JD, Ludwig W. Prognostic importance of the white blood cell count for coronary, cancer, and all-cause mortality. JAMA. 1985;254: 1932-1937.
Weijenberg MP, Feskens eJ, Kromhout D. White blood cell count and the risk of coronary heart disease and all-cause mortality in elderly men. Arterioscler Thromb Vasc Biol. 1996;16: 499-503.
Castro O, Brambilla DJ, Thorington B, et al. The acute chest syndrome in sickle cell disease: incidence and risk factors. The Cooperative Study of Sickle Cell Disease. Blood. 1994;84: 643-649.
Kinney TR, Sleeper LA, Wang WC, et al. Silent cerebral infarcts in sickle cell anemia: a risk factor analysis. The Cooperative Study of Sickle Cell Disease. Pediatrics. 1999;103: 640-645.
Platt OS, Brambilla DJ, Rosse WF, et al. Mortality in sickle cell disease: life expectancy and risk factors for early death. N engl J Med. 1994;330: 1639-1644.
Miller ST, Sleeper LA, Pegelow CH, et al. Prediction of adverse outcomes in children with sickle cell disease. N engl J Med. 2000;342: 83-89.
Falconer D. Introduction to Quantitative Genetics. 3rd ed. essex, United Kingdom: Longman; 1989.
Hampton JW, Matthews C. Similarities between baboon and human blood clotting. J Appl Physiol. 1966;21: 1713-1716.
Todd Me, McDevitt e, Goldsmith eI. Blood-clotting mechanisms of nonhuman primates: choice of the baboon model to simulate man. J Med Primatol. 1972;1: 132-141.
Harker LA, Hanson SR. experimental arterial thromboembolism in baboons: mechanism, quantitation, and pharmacologic prevention. J Clin Invest. 1979;64: 559-560.
Hanson SR, Harker LA. Baboon models of acute arterial thrombosis. Thromb Haemost. 1987;58: 801-805.
Jolly CJ. Species, subspecies, and baboon systematics. In: Kimbel WH, Martin LB, eds. Species, Species Concepts, and Primate evolution. New York and London: Plenum Press; 1993: 67-107.
Williams-Blangero S, Vandeb, Vandeberg J, Blangero L, Konigsbewrg L, Dyke B. Genetic differentiation between baboon subspecies: relevance for biomedical research. Am J Primatol. 1990;20: 67-81.
Council NR. Guide for care and use of laboratory animals. Washington, D.C.: National academy of sciences; 1996.
Hainsey BM, Hubbard GB, Leland MM, Brasky KM. Clinical parameters of the normal baboons (Papio species) and chimpanzees (Pan troglodytes). Lab Anim Sci. 1993;43: 236-243.
Almasy L, Blangero J. Multipoint quantitative-trait linkage analysis in general pedigrees. Am J Hum Genet. 1998;62: 1198-1211.
Hopper JL, Mathews JD. extensions to multivariate normal models for pedigree analysis. Ann Hum Genet. 1982;46: 373-383.
Lange K, Boehnke M. extensions to pedigree analysis, IV: covariance components models for multivariate traits. Am J Med Genet. 1983;14: 513-524.
Wang XL, Mahaney MC, Sim AS, et al. Genetic contribution of the endothelial constitutive nitric oxide synthase gene to plasma nitric oxide levels. Arterioscler Thromb Vasc Biol. 1997;17: 3147-3153.
Boehnke M, Moll PP, Kottke BA, Weidman WH. Partitioning the variability of fasting plasma glucose levels in pedigrees: genetic and environmental factors. Am J epidemiol. 1987;125: 679-689.
Mahaney MC, Blangero J, Comuzzie AG, Vande-Berg JL, Stern MP, MacCluer JW. Plasma HDL cholesterol, triglycerides, and adiposity: a quantitative genetic test of the conjoint trait hypothesis in the San Antonio Family Heart Study. Circulation. 1995;92: 3240-3248.
edwards A. Likelihood. Baltimore: Johns Hopkins Press; 1992.
Williams JT, Van eerdewegh P, Almasy L, Blangero J. Joint multipoint linkage analysis of multivariate qualitative and quantitative traits, I: likelihood formulation and simulation results. Am J Hum Genet. 1999;65: 1134-1147.
Williams JT, Begleiter H, Porjesz B, et al. Joint multipoint linkage analysis of multivariate qualitative and quantitative traits, II: alcoholism and event-related potentials. Am J Hum Genet. 1999; 65: 1148-1160.
Garner C, Tatu T, Reittie Je, et al. Genetic influences on F cells and other hematologic variables: a twin heritability study. Blood. 2000;95: 342-346.
Whitfield JB, Martin NG. Genetic and environmental influences on the size and number of cells in the blood. Genet epidemiol. 1985;2: 133-144.
evans DM, Frazer IH, Martin NG. Genetic and environmental causes of variation in basal levels of blood cells. Twin Res. 1999;2: 250-257.
Norton JM, Rand PW. A relationship between erythrocyte volume and concentration in humans and other mammals. Blut. 1980;41: 421-426.
Graham SS, Traub B, Mink IB. Automated platelet-sizing parameters on a normal population. Am J Clin Pathol. 1987;87: 365-369.
Levin J, Bessman JD. The inverse relation between platelet volume and platelet number: abnormalities in hematologic disease and evidence that platelet size does not correlate with platelet age. J Lab Clin Med. 1983;101: 295-307.
Bancroft AJ, Abel eW, McLaren M, Belch JJ. Mean platelet volume is a useful parameter: a reproducible routine method using a modified Coulter thrombocytometer. Platelets. 2000;11: 379-387.
Ruan C, Wu Y, Okada T, et al. Studies on reference intervals for platelet counts in three cities in China and one in Japan. Thromb Haemost. 2002; 88: 111-114.
Jansson L. Relationship of mean platelet volume to platelet count in morphologic evaluation of thrombopoiesis. Scand J Haematol. 1985;35: 363-366.
Giles C. The platelet count and mean platelet volume. Br J Haematol. 1981;48: 31-37.
Bessman JD, Williams LJ, Gilmer PR, Jr. Mean platelet volume: the inverse relation of platelet size and count in normal subjects, and an artifact of other particles. Am J Clin Pathol. 1981;76: 289-293.
Bain BJ. Platelet count and platelet size in males and females. Scand J Haematol. 1985;35: 77-79.
Ingram M, Coopersmith A. Reticulated platelets following acute blood loss. Br J Haematol. 1969; 17: 225-229.
Niethammer AG, Forman eN. Use of the platelet histogram maximum in evaluating thrombocytopenia. Am J Hematol. 1999;60: 19-23.
Bath P, Algert C, Chapman N, Neal B. Association of mean platelet volume with risk of stroke among 3134 individuals with history of cerebrovascular disease. Stroke. 2004;35: 622-626.
endler G, Klimesch A, Sunder-Plassmann H, et al. Mean platelet volume is an independent risk factor for myocardial infarction but not for coronary artery disease. Br J Haematol. 2002;117: 399-404.
Martin JF, Bath PM, Burr ML. Influence of platelet size on outcome after myocardial infarction. Lancet. 1991;338: 1409-1411.(Michael C. Mahaney, Carlo)
Southwest National Primate Research Center, San Antonio, TX
the Department of Laboratory Medicine, Children's Hospital, Boston, MA.
Abstract
Interperson differences in peripheral blood cell counts in healthy individuals result from genetic and environmental influences. We used multivariate genetic analyses to assess the relative impact of genes and environment on baseline blood cell counts and indices using a pedigreed colony of baboons, an animal with well-documented analogies to human blood physiology. After accounting for age, sex, and weight, we found that genetic influences explain a significant proportion of the remaining variability, ranging from a low of 13.7% for mean corpuscular hemoglobin concentration (MCHC) to a high of 72.4% for red blood cell (RBC) number. Genes influence 38.5% of the variation in baseline white blood cell (WBC) count, a characteristic that correlates with mortality in both the general human population and clinically defined subgroups such as individuals with sickle-cell disease. We examined the interaction between pairs of traits and identified those that share common genetic influences (pleiotropy). We unexpectedly observed that the same gene or group of genes influences both WBC count and mean platelet volume (MPV). We anticipate that this approach will ultimately lead to discovery of novel insights into the biology of related traits, and ultimately identify important genes that affect hematopoiesis.
Introduction
Peripheral blood cell counts represent a complex clinical phenotype—the end result of interactions between a multitude of genetic and environmental factors. Increasingly, baseline blood cell count phenotypes are emerging as important markers of risk for disease and disease severity. examples include the correlation of steady-state white blood cell (WBC) count with mortality in the general population,1-6 and with disease severity in patients with sickle-cell anemia (SS).7-10 efforts to dissect the genetic underpinnings of these complex traits have the potential to greatly enhance our understanding of hematopoiesis, pathophysiology, and disease modifier genes.
The primary genetic questions about quantitative traits such as peripheral blood cell counts (ie, traits measured on continuous rather than dichotomous scales) are formulated in terms of variation in the traits.11 The amount of variation is expressed as the variance, or in statistical terms, the mean of the squared deviations from the population mean for the trait. The fundamental objective of statistical genetic analyses of any quantitative trait is to partition the total variance for that trait (the phenotype) into 2 components: those due to the effects of genes (the genetic component), and those influenced by environmental factors (the environmental component). In this study, we examine the patterns of normal quantitative variation in blood cell counts in a pedigreed colony of baboons, a nonhuman primate with well-documented analogies to human blood physiology,12-15 in order to partition the blood cell phenotypes into their genetic and environmental components. Our specific aims are to quantify the genetic and environmental components, and determine which genetic effects on one cell count may also account for variation in other cell counts (pleiotropy).
Materials and methods
Subjects
We analyzed the blood counts and indices of 600 individual healthy pedigreed baboons (Papio hamadryas)16,17 to quantify the range of values observed in this normal population. Housed at the Southwest National Primate Research Center, located at the Southwest Foundation for Biomedical Research (San Antonio, TX), this sample consisted mainly of olive baboons (Papio hamadryas anubis), yellow baboons (Papio hamadryas cynocephalus), and their hybrids, and initially included 415 females and 185 males, ranging in age from 5.5 to 30 years. While strict genetic management is used to control inbreeding, all nonfounder animals in this study were the results of matings that were random with respect to phenotypes analyzed here. The subjects were distributed among 11 pedigrees with 67 to 171 members each. Within these pedigrees there were 4714 relative pairs (Table 1).
Since birth or arrival at the colony, all animals had been housed out of doors in social group cages and maintained on monkey chow diets to which they had ad libitum access. Animal care personnel and staff veterinarians provided daily maintenance and regular, urgent, and emergent health care to all animals in accordance with the Guide for the Care and Use of Laboratory Animals.18 All procedures related to their treatment during the conduct of this study were approved by the Institutional Animal Care and Use Committee in accordance with the established guidelines.18
Hematologic studies
After sedation with intramuscular ketamine and intravenous xylazine hydrochloride, atropine, acepromazine, and ketamine to assure relaxation, blood samples were obtained by femoral venipuncture. Blood was collected in eDTA (ethylenediaminetetraacetic acid), stored on ice, and analyzed within 4 hours using a Coulter JT (Beckman Coulter, Fullerton, CA), conditions that maintained the stability of the measured parameters. We measured hematocrit (HCT), hemoglobin (Hb), red blood cell (RBC) count, mean corpuscular volume (MCV), red cell distribution width (RDW), mean corpuscular hemoglobin (MCH), mean corpuscular hemoglobin concentration (MCHC), white blood cell (WBC) count, platelet (Plt) count, and mean platelet volume (MPV). Details regarding protocols for these assays, and measurements are reported elsewhere.19
Statistical analysis
Statistical genetic analyses were conducted by means of a maximum likelihood–based variance decomposition approach implemented in the computer package SOLAR (Sequential Oligogenic Linkage Analysis Routines; Southwest Foundation for Biomedical Research, San Antonio, TX).20 This approach, which was developed following methodology originally proposed by Hopper and Mathews21 and Lange and Boehnke,22 has been described in detail elsewhere.23
In short, we used this approach to partition the phenotypic variance (2P) into components corresponding to the additive genetic (2G) and nongenetic (ie, environmental) (2e) effects. Because these components are additive, such that 2P =2G +2e, we estimated the heritability, or proportion of the phenotypic variance attributable to additive genetic effects, as h2 =2G /2P. We estimated the proportion of the phenotypic variance attributable to nongenetic factors as e2 = 1 – h2. When covariate effects are included in the maximum likelihood model, the h2 and e2 estimates are obtained as proportions of the residual phenotypic variance in the trait (ie, that remaining after the accounting for the effects of the covariates) rather than as proportions the total phenotypic variance. Multiplication of h2 by the proportion of the total variance not explained by covariates (ie, the quantity 1 – c2, where "c2" refers to the proportion of the total variance attributable to covariate effects) yields a heritability value that is rescaled to the total phenotypic variance. We also simultaneously estimated the mean effects of sex, sex-specific age, and sex-specific age,2 and body weight on each of the hematological phenotypes. We model the phenotype of an individual (yi) as linear function of the population mean of the trait (μ), the covariates (x) and their regression coefficients (), plus the additive genetic values (g) and random environmental deviations (e), such that yi = μ + 1 x1+ 2 x2 + nxn + g + e.
Using extensions to univariate genetic analysis that encompass the multivariate state,21,22,24 we followed an approach presented in detail previously25 to similarly model multivariate phenotypes in these baboons (ie, to simultaneously analyze measurements on pairs of traits in individuals within each pedigree). By maximizing these multivariate models, we obtained estimates of 2 additional parameters: the additive genetic correlation, G, and the environmental correlation, e, between trait pairs. Respectively, these correlations are estimates of the shared additive effects of genes (ie, pleiotropy) and shared environmental (ie, unmeasured, nongenetic) factors on the variance in a trait. Genetic correlations with values equal to (or not significantly different from) 1 or –1 are interpreted as evidence for complete pleiotropy (ie, where the additive genetic contributions to the variance in the 2 correlated traits are due entirely to the effects of the same gene or genes). Genetic correlations that are significantly different from both 0 and 1 or – 1 are evidence of incomplete pleiotropy (ie, where some but not all of the additive genetic effects on a pair of traits are due to the same gene or genes). The genetic and environmental components of the phenotypic correlation matrix are additive, like those of the corresponding variance-covariance matrix, so we could use the maximum likelihood estimates of the additive genetic and environmental correlations, to obtain the total phenotypic correlation between 2 traits, P, as .
Using this approach as implemented in SOLAR,20 we conducted a series of bivariate quantitative genetic analysis of all hematologic trait pairs.
Significance of the maximum likelihood estimates for heritability and other parameters was assessed by means of likelihood ratio tests.26 The maximum likelihood for the general model in which all parameters were estimated was compared to that for restricted models in which the value of the parameter to be tested was held constant at some value (usually 0). For tests of parameters like h2, for which a fixed value of 0 in a restricted model is at a boundary of the parameter space, twice the difference in the natural logarithm (ln) of the likelihoods of the 2 models compared is distributed asymptotically approximately as a 0.5:0.5 mixture of 2 with a point mass at 0. For tests of covariates, for which 0 is not a boundary value, this statistic is distributed asymptotically approximately as a 2 variate.21 Degrees of freedom can be obtained as the difference in the number of estimated parameters in the 2 models.21 However, in tests of parameters like h2, whose values may be fixed at a boundary of their parameter space in the null model, the appropriate significance level is obtained by halving the P value.24
Results
Inheritance of the individual peripheral blood cell counts
In these healthy pedigreed baboons, who live in a highly controlled environment designed to maximize the detection of genetic signals, a significant fraction of the phenotypic variance of all of the peripheral blood cell counts and indices is under genetic control (Table 2). After accounting for the individual and combined effects of age, sex, and weight, the additive effects of genes (h2) account for a low of 13.7% of the observed variability in MCHC, to a high of 72.4% of the variability of RBC counts.
Phenotypic correlations between the individual cell traits
Many of the cell counts and indices are correlated with others. As shown in Table 3, the erythrocyte traits are highly correlated, some positively (eg, Hct and Hb), and some negatively (eg, RBC counts and MCV). Platelet count and MPV are negatively correlated. WBC count and MPV are positively correlated; WBC count and Plt are negatively correlated.
evidence that shared genes influence more than one trait
We had not anticipated finding significant genetic influences shared between WBC count and MPV (P < .05). In fact, similar to what we found in the case of Hct and Hb, the genetic effects shared between WBC count and MPV is virtually 100% (ie, complete pleiotropy (G = 1).
Discussion
Most of what we know about the genetics of the peripheral blood counts and indices comes from the study of patients and experimental animals whose counts lie outside the normal range. Humans and animals with too many or too few cells that are too big or too small have long been the focus of attention of hematologists, and have established the pathophysiologic importance of a variety of genes that control phenotype in a Mendelian fashion. Some of these genes include those that encode or control globins, membrane proteins, growth factors, receptors, ion channels, and adhesion molecules. In much the same way, clinical investigators have discovered critical environmental (ie, noninherited) modulators of blood counts such as normal aging, maturation and growth, pregnancy, diet (iron, folate, cobalamin), hypoxia, bleeding, smoking, infection, inflammation, and exogenous growth factors and cytotoxic agents.
In this study, we shift focus from the classical clinical study of hematologically abnormal individuals to a statistical genetic methodology that analyzes the normal variation of cell counts in a healthy pedigreed population (ie, a population of individuals that were not ascertained on the basis of a disease state). Conceptually, this method is very similar to studies of healthy human twins. In twin studies, the genetic influence is estimated by comparing degrees of similarity between 2 types of relative pairs: monozygotic and dizygotic twins. Studies of more complex, multigenerational pedigrees such as ours make use of information from a much larger range of relative pairs comprising parents and offspring, grandparents and grandchildren, siblings, cousins, aunts/uncles, and nieces/nephews (the "avuncular" relationship), etc. Although not unique to nonhuman primates, the baboon pedigrees from which data for this study were obtained also include large numbers of pairs exhibiting half relationships (eg, half-sibling and half-avuncular pairs). Because it uses information from more relative classes, this strategy offers demonstrably greater statistical power to both detect genetic effects and, when combined with genotype data, efficiently localize quantitative trait loci to specific chromosomal regions.27,28 This approach considers the peripheral blood count a complex phenotype that is controlled by multiple underlying genetic and environmental factors, and is designed to quantify the impact and relatedness of those factors. We studied the baboon because of its close analogies to the human, and the ability to minimize the number and impact of environmental factors on interindividual variation.
We found that in the baboon, the steady-state set points for the number and size of peripheral blood cells are under significant genetic control (Table 1). These results are extremely similar to findings in humans by investigators working with major twin registries in the United Kingdom and Australia,29-31 and underscore the advantages of using such populations to discover genetic influences.
Our analyses further allowed us to examine the patterns of interrelatedness among these traits and, more specifically, to dissect out the genetic and random environmental components underlying these patterns. The positive phenotypic correlations between the various Hb, RBC count, and Hct pairings, and negative correlation of Hb, Hct, and RBC count with MCV are generally similar to the patterns of correlations observed between these traits in healthy human30-32 and other mammalian32 populations, although they have not been examined in detail. Many of these pairs demonstrated significant genetic contribution to the relationships—consistent with the findings that emerged from the Australian twin registry study.30 A similar inverse phenotypic relationship between Plt and MPV reported by Whitfield and Martin30 and others33-40 was also found in this baboon study, but other studies did not decompose the phenotypic covariance between platelet count and MPV into genetic and environmental components.
The positive phenotypic correlation that we found between WBC count and MPV has not been previously reported; and the observation of complete pleiotropy (ie, that all or most of the genetic effects on WBC count [responsible for 37% of the total phenotypic variance] and all of the genetic effects on MPV [49% of the total phenotypic variance] are due to the effects of the same gene or group of genes) is remarkable. Conceptually, the MPV resembles the MCV. Like the reticulocyte, the newly formed platelet is large, contains mRNA,41 and as an analogy—the increased number of reticulocytes in immune hemolytic anemia causes an elevated MCV, just as the high number of young platelets in idiopathic thrombocytopenic purpura elevates the MPV.42
Genetic abnormalities may make red cells small (thalassemia) or large (stomatocytosis), but we do not know what genes are responsible for that portion of variation in MCV (or MPV) that is due to the effects of genes in the healthy population. It may well be that these genes have little or no overlap with the globin genes that so dramatically influence cell size in abnormal populations.
It is interesting to note that high MPV43-45 and high WBC count1-6 both correlate with heart disease and stroke risk in the healthy population. While these may represent 2 genetically independent risk markers, it is tempting to speculate that genes exerting pleiotropic effects on both traits underlie the risk. Perhaps the genetic correlation between WBC count and platelet volume in relatively healthy individuals is attributable to the effects of genes that influence steps in hematopoiesis common to megakaryocytes and white cells, or genes that allow more mature platelets and white cells to respond to or enhance inflammation. It is important to point out that we do not know how much of the heart disease/stroke risk is influenced by the genetic component of the variation in WBC count and MPV, and how much by environmental factors. Similarly, we know that patients with sickle-cell anemia and high severity have high WBC counts, but do not know whether the high WBC count is an "environmental" (ie, nongenetic) reaction to disease severity, or whether the patients with a genetic predisposition to having a high WBC count are more likely to have severe disease as an expression of the role white blood cells play in disease pathophysiology.
Multivariate statistical genetic analyses of the type reported here offer a means to identify pairs and, in some cases, networks of traits that are bound together by the shared effects of the same gene or group of genes. In many instances, the detected pleiotropic relationships will be entirely consonant with our current expectations regarding the phenotypes under study (eg, RBC count and hemoglobin-related measures). In other cases however, previously unappreciated relationships, such as that between WBC count and MPV in this study, may be revealed or suggested by these sorts of analyses. The pleiotropic relationships detected in this study suggest a rationale for subsequent multivariate statistical analyses in both this pedigreed animal population, and in human families. Such analyses not only offer potentially novel insights into the biology of related traits but, in comparison to analyses of single traits, they also provide substantially greater amounts of information and improved statistical power to detect, localize, and identify quantitative trait loci contributing to the pleiotropic effects. Our plans to combine the data presented in this paper with genotypic data we are developing, will allow us to map the quantitative trait loci underlying patterns of covariation among these phenotypes and identify positional candidate genes in this nonhuman primate model for human hematological variation.
Footnotes
Prepublished online as Blood First edition Paper, May 3, 2005; DOI 10.1182/blood-2004-12-4863.
Supported by grants from the National Institutes of Health (NIH) grants R01 HL068922 (O.S.P.), R01 HL054141 (M.C.M), P01 HL028972 (M.C.M.), and P51 RR013086 (M.C.M.).
The publication costs of this article were defrayed in part by page charge payment. Therefore, and solely to indicate this fact, this article is hereby marked "advertisement" in accordance with 18 U.S.C. section 1734.
References
de Labry LO, Campion eW, Glynn RJ, Vokonas PS. White blood cell count as a predictor of mortality: results over 18 years from the Normative Aging Study. J Clin epidemiol. 1990;43: 153-157.
James AL, Knuiman MW, Divitini ML, Musk AW, Ryan G, Bartholomew HC. Associations between white blood cell count, lung function, respiratory illness and mortality: the Busselton Health Study. eur Respir J. 1999;13: 1115-1119.
Burr ML, Holliday RM, Fehily AM, Whitehead PJ. Haematological prognostic indices after myocardial infarction: evidence from the diet and reinfarction trial (DART). eur Heart J. 1992;13: 166-170.
Gillum RF, Ingram DD, Makuc DM. White blood cell count, coronary heart disease, and death: the NHANeS I epidemiologic Follow-up Study. Am Heart J. 1993;125: 855-863.
Grimm RH, Jr., Neaton JD, Ludwig W. Prognostic importance of the white blood cell count for coronary, cancer, and all-cause mortality. JAMA. 1985;254: 1932-1937.
Weijenberg MP, Feskens eJ, Kromhout D. White blood cell count and the risk of coronary heart disease and all-cause mortality in elderly men. Arterioscler Thromb Vasc Biol. 1996;16: 499-503.
Castro O, Brambilla DJ, Thorington B, et al. The acute chest syndrome in sickle cell disease: incidence and risk factors. The Cooperative Study of Sickle Cell Disease. Blood. 1994;84: 643-649.
Kinney TR, Sleeper LA, Wang WC, et al. Silent cerebral infarcts in sickle cell anemia: a risk factor analysis. The Cooperative Study of Sickle Cell Disease. Pediatrics. 1999;103: 640-645.
Platt OS, Brambilla DJ, Rosse WF, et al. Mortality in sickle cell disease: life expectancy and risk factors for early death. N engl J Med. 1994;330: 1639-1644.
Miller ST, Sleeper LA, Pegelow CH, et al. Prediction of adverse outcomes in children with sickle cell disease. N engl J Med. 2000;342: 83-89.
Falconer D. Introduction to Quantitative Genetics. 3rd ed. essex, United Kingdom: Longman; 1989.
Hampton JW, Matthews C. Similarities between baboon and human blood clotting. J Appl Physiol. 1966;21: 1713-1716.
Todd Me, McDevitt e, Goldsmith eI. Blood-clotting mechanisms of nonhuman primates: choice of the baboon model to simulate man. J Med Primatol. 1972;1: 132-141.
Harker LA, Hanson SR. experimental arterial thromboembolism in baboons: mechanism, quantitation, and pharmacologic prevention. J Clin Invest. 1979;64: 559-560.
Hanson SR, Harker LA. Baboon models of acute arterial thrombosis. Thromb Haemost. 1987;58: 801-805.
Jolly CJ. Species, subspecies, and baboon systematics. In: Kimbel WH, Martin LB, eds. Species, Species Concepts, and Primate evolution. New York and London: Plenum Press; 1993: 67-107.
Williams-Blangero S, Vandeb, Vandeberg J, Blangero L, Konigsbewrg L, Dyke B. Genetic differentiation between baboon subspecies: relevance for biomedical research. Am J Primatol. 1990;20: 67-81.
Council NR. Guide for care and use of laboratory animals. Washington, D.C.: National academy of sciences; 1996.
Hainsey BM, Hubbard GB, Leland MM, Brasky KM. Clinical parameters of the normal baboons (Papio species) and chimpanzees (Pan troglodytes). Lab Anim Sci. 1993;43: 236-243.
Almasy L, Blangero J. Multipoint quantitative-trait linkage analysis in general pedigrees. Am J Hum Genet. 1998;62: 1198-1211.
Hopper JL, Mathews JD. extensions to multivariate normal models for pedigree analysis. Ann Hum Genet. 1982;46: 373-383.
Lange K, Boehnke M. extensions to pedigree analysis, IV: covariance components models for multivariate traits. Am J Med Genet. 1983;14: 513-524.
Wang XL, Mahaney MC, Sim AS, et al. Genetic contribution of the endothelial constitutive nitric oxide synthase gene to plasma nitric oxide levels. Arterioscler Thromb Vasc Biol. 1997;17: 3147-3153.
Boehnke M, Moll PP, Kottke BA, Weidman WH. Partitioning the variability of fasting plasma glucose levels in pedigrees: genetic and environmental factors. Am J epidemiol. 1987;125: 679-689.
Mahaney MC, Blangero J, Comuzzie AG, Vande-Berg JL, Stern MP, MacCluer JW. Plasma HDL cholesterol, triglycerides, and adiposity: a quantitative genetic test of the conjoint trait hypothesis in the San Antonio Family Heart Study. Circulation. 1995;92: 3240-3248.
edwards A. Likelihood. Baltimore: Johns Hopkins Press; 1992.
Williams JT, Van eerdewegh P, Almasy L, Blangero J. Joint multipoint linkage analysis of multivariate qualitative and quantitative traits, I: likelihood formulation and simulation results. Am J Hum Genet. 1999;65: 1134-1147.
Williams JT, Begleiter H, Porjesz B, et al. Joint multipoint linkage analysis of multivariate qualitative and quantitative traits, II: alcoholism and event-related potentials. Am J Hum Genet. 1999; 65: 1148-1160.
Garner C, Tatu T, Reittie Je, et al. Genetic influences on F cells and other hematologic variables: a twin heritability study. Blood. 2000;95: 342-346.
Whitfield JB, Martin NG. Genetic and environmental influences on the size and number of cells in the blood. Genet epidemiol. 1985;2: 133-144.
evans DM, Frazer IH, Martin NG. Genetic and environmental causes of variation in basal levels of blood cells. Twin Res. 1999;2: 250-257.
Norton JM, Rand PW. A relationship between erythrocyte volume and concentration in humans and other mammals. Blut. 1980;41: 421-426.
Graham SS, Traub B, Mink IB. Automated platelet-sizing parameters on a normal population. Am J Clin Pathol. 1987;87: 365-369.
Levin J, Bessman JD. The inverse relation between platelet volume and platelet number: abnormalities in hematologic disease and evidence that platelet size does not correlate with platelet age. J Lab Clin Med. 1983;101: 295-307.
Bancroft AJ, Abel eW, McLaren M, Belch JJ. Mean platelet volume is a useful parameter: a reproducible routine method using a modified Coulter thrombocytometer. Platelets. 2000;11: 379-387.
Ruan C, Wu Y, Okada T, et al. Studies on reference intervals for platelet counts in three cities in China and one in Japan. Thromb Haemost. 2002; 88: 111-114.
Jansson L. Relationship of mean platelet volume to platelet count in morphologic evaluation of thrombopoiesis. Scand J Haematol. 1985;35: 363-366.
Giles C. The platelet count and mean platelet volume. Br J Haematol. 1981;48: 31-37.
Bessman JD, Williams LJ, Gilmer PR, Jr. Mean platelet volume: the inverse relation of platelet size and count in normal subjects, and an artifact of other particles. Am J Clin Pathol. 1981;76: 289-293.
Bain BJ. Platelet count and platelet size in males and females. Scand J Haematol. 1985;35: 77-79.
Ingram M, Coopersmith A. Reticulated platelets following acute blood loss. Br J Haematol. 1969; 17: 225-229.
Niethammer AG, Forman eN. Use of the platelet histogram maximum in evaluating thrombocytopenia. Am J Hematol. 1999;60: 19-23.
Bath P, Algert C, Chapman N, Neal B. Association of mean platelet volume with risk of stroke among 3134 individuals with history of cerebrovascular disease. Stroke. 2004;35: 622-626.
endler G, Klimesch A, Sunder-Plassmann H, et al. Mean platelet volume is an independent risk factor for myocardial infarction but not for coronary artery disease. Br J Haematol. 2002;117: 399-404.
Martin JF, Bath PM, Burr ML. Influence of platelet size on outcome after myocardial infarction. Lancet. 1991;338: 1409-1411.(Michael C. Mahaney, Carlo)