ABSTRACT
The evolution of primates involved increasing body size, brain size and presumably cognitive ability. Cognition is related to neural activity, metabolic rate and rate of blood flow to the cerebral cortex. These parameters are difficult to quantify in living animals. This study shows that it is possible to determine the rate of cortical brain perfusion from the size of the internal carotid artery foramina in skulls of certain mammals, including haplorrhine primates and diprotodont marsupials. We quantify combined blood flow rate in both internal carotid arteries as a proxy of brain metabolism in 34 species of haplorrhine primates (0.116–145 kg body mass) and compare it to the same analysis for 19 species of diprotodont marsupials (0.014–46 kg). Brain volume is related to body mass by essentially the same exponent of 0.70 in both groups. Flow rate increases with haplorrhine brain volume to the 0.95 power, which is significantly higher than the exponent (0.75) expected for most organs according to ‘Kleiber's Law’. By comparison, the exponent is 0.73 in marsupials. Thus, the brain perfusion rate increases with body size and brain size much faster in primates than in marsupials. The trajectory of cerebral perfusion in primates is set by the phylogenetically older groups (New and Old World monkeys, lesser apes) and the phylogenetically younger groups (great apes, including humans) fall near the line, with the highest perfusion. This may be associated with disproportionate increases in cortical surface area and mental capacity in the highly social, larger primates.
INTRODUCTION
The evolution of primates involved an overall trend for increased body size, absolute brain size and relative brain size, although there are reversals in some primate clades (Montgomery et al., 2010). This overall size increase is presumably associated with increased brain metabolic rate and cognitive ability, leading to the great apes and humans. The human brain is an energy-expensive tissue (Aiello and Wheeler, 1995) and much current interest in primate brain evolution concerns how the high cost is compensated in primates compared with other mammals (Navarrete et al., 2011; Pontzer et al., 2014). Recent research shows that the number of neurons is directly proportional to brain size in primates (Herculano-Houzel, 2011c) and that the metabolic rate of the brain appears directly proportional to the number of neurons (Herculano-Houzel, 2011a). Therefore, within primates, brain metabolic rate is thought to be proportional to brain size and the large human brain appears to be simply a scaled-up version of smaller primate brains (Herculano-Houzel, 2012). However, the relationship between brain metabolic rate and brain size in mammals may not follow direct proportionality. When brain metabolic rate in relation to brain mass was analysed with an allometric equation of the form Y=aXb, the exponent b was 0.82 when based on data from rats, dogs and humans (Wang et al., 2001). More recently, a meta-analysis of direct measurements of brain metabolic rate across 10 species of mammal, from mouse to human, showed that it scales with brain volume to the 0.86 power, which is significantly higher than the commonly accepted exponent of 0.75 required by ‘Kleiber's Law’ (Karbowski, 2007). Laboratory and domestic mammals occupy the small end of the brain size distribution (0.35–72 ml) in this analysis and the primates are at the large end, with large brains (100–1389 ml). A larger sample of primate species over a greater range of body and brain size would not only provide a direct test of these important concepts, but would also reflect on the evolution of brain metabolic rate and possibly cognitive ability among primates.
Brain metabolic rate is related to brain size and the intensity of neuronal and synaptic activity (Attwell and Laughlin, 2001). Compared with the rest of the body, the metabolic rate of the brain is high and is similar during mental activity and sleep (Armstrong, 1983), as well as during rest and moderate exercise (Hiura et al., 2014; Wilcox et al., 1970). On a tissue-mass-specific basis, brain metabolic rate in a resting mammal is about two-thirds of the rates of the most active organs (heart, liver, kidney) and more than 10-times higher than that of the rest of the body (Wang et al., 2001). Continuous cranial blood flow is required because of the aerobic nature of the brain and its inability to store glucose or glycogen (Karbowski, 2011; Weisbecker and Goswami, 2010). Blood flow rates in the internal carotid arteries are likely to indicate the levels of activity in the brain, because these vessels mainly supply the cerebrum (with minor supplies to the meninges, eyes and forehead scalp), whereas the vertebral arteries mainly supply the cerebellum in most mammals (Turnquist and Minugh-Purvis, 2012). Moreover, metabolic rate is much greater in a cerebral neuron than a cerebellar neuron (Herculano-Houzel, 2011c). Thus, blood flow rates in the internal carotid arteries of humans are considerably higher than in the vertebral arteries (Schöning et al., 1994). Direct measurements of carotid arterial blood flow rate and brain blood flow rate show a close correlation (van Bel et al., 1994).
Local cerebral blood flow in the internal carotid arteries is considered proportional to the metabolic rate of the cerebral cortex (Changizi, 2001; Hawkins et al., 1983; Lou et al., 1987; Turnquist and Minugh-Purvis, 2012). This is demonstrated quite clearly by similar allometric exponents for brain oxygen consumption rate, glucose utilization rate and cerebral perfusion rate in relation to brain volume, Vbr. The rate of oxygen uptake by the whole brain of seven mammalian species scales with Vbr0.86 and the rate of glucose metabolism in ten species also scales with Vbr0.86 (Karbowski, 2007). The scaling of cortical blood flow rate and cortical brain volume in up to seven species of mammals varies with exponents between 0.81 and 0.87, depending on the particular part of the cortex, but the overall exponent is 0.84 (Karbowski, 2011). This indicates that there is good matching of brain metabolic rate and blood flow rate and confirms the notion that blood flow rate can indicate metabolic rate. The implication of proportionality between perfusion and metabolic rate is that O2 extraction from the cerebral blood is constant. Thus, it is not necessary to know arterial–venous O2 contents, O2 equilibrium curves, haemoglobin levels, transit times, intercapillary distances or other variables associated with O2 delivery.
The metabolic rate of an organ influences the size of the supplying artery. Blood vessels undergo continuous remodelling throughout life in response to local haemodynamic forces, including blood pressure that affects vessel wall thickness, and blood flow rate that induces shear stress and causes changes in vessel lumen diameter (Lehoux et al., 2006). Vessel remodelling ensures that the energy needed for perfusion is minimized by maintaining basal levels of shear stress and producing energy-efficient laminar flow (Ward et al., 2000). Chronic and acute experimental changes in flow rates or blood pressure induce adaptive alteration in vessel geometry and normalize shear stress within narrow limits (Kamiya et al., 1984; Kamiya and Togawa, 1980; Langille, 1999; Smiesko and Johnson, 1993; Tronc et al., 1996). Blood flow rate Q̇ is related to shear stress τ, vessel radius r, and blood viscosity η according to the equation: Q̇ = (τπr3)/(4η) (Lehoux and Tedgui, 2003). Shear stress is the viscous force exerted by the blood flowing along the walls of the vessel, and it is calculated as proportional to the maximum velocity at the axis of the vessel and inversely proportional to vessel radius, assuming a parabolic velocity distribution between the axis and the wall, where velocity is theoretically zero (Jeong and Rosenson, 2013). The vascular epithelium senses blood flow near the wall and induces a molecular mechanism of vascular remodelling that is well understood (Lu and Kassab, 2011). The biophysics involving arterial wall thickness is also well known. Larger arteries have thicker walls because wall cross-sectional tension is normalized according to the principle of Laplace that requires the wall thickness to increase in proportion to blood pressure and the radius of the vessel (Burton, 1965). Mean central arterial blood pressure is essentially invariant across species, except in mammals with long vertical distances between the heart and the head (White and Seymour, 2014), and blood viscosity is also invariant in mammals (Schmid-Schönbein et al., 1969). Therefore, measurements of external arterial dimensions can indicate the rate of blood flow within them. Where an artery passes through a foramen in a bone, the foramen size can be a gauge of arterial size, which is the key innovation that has allowed estimates of bone perfusion in living and extinct vertebrates (Allan et al., 2014; Seymour et al., 2012).
The present study estimates cerebral perfusion from the radius of the internal carotid foramen, measured in skulls of haplorrhine primates and diprotodont marsupials. This comparison was made for a number of reasons. The body sizes of haplorrhine primates and diprotodont marsupials overlap, with humans equivalent to large kangaroos. The blood supply to the forebrain of haplorrhine primates is mainly derived from the internal carotid arteries that service all of the tissues anterior to the hippocampus, whereas the vertebral arteries supply mainly the hindbrain, although there are connections between the two within the cerebral arterial circle (Circle of Willis) (Coceani and Gloor, 1966; Scremin, 2011). The same is true for diprotodonts; the internal carotid arteries supply the forebrain through discrete internal carotid foramina, whereas the vertebral arteries supply the hindbrain through the foramen magnum (Voris, 1928). Haplorrhine primates and marsupials have no additional arterial supplies to the forebrain (Gillilan, 1972). The analysis is not possible in most other groups of placental mammals, because the blood supply is derived from a complex of arteries, often involving retia and entering the skull through various foramina. Strepsirrhine primates have much reduced internal carotid arteries (Schwartz and Tattersall, 1987), nevertheless, we include them for comparison in this study. Examples of other inappropriate groups unfortunately include cetaceans (Geisler and Luo, 1998; Viamonte et al., 1968; Vogl and Fisher, 1981), carnivores (Davis and Story, 1943; Gillilan, 1976) and artiodactyls (Diéguez et al., 1987; Gillilan, 1974). Rodents (Gillilan, 1974) and rabbits (Scremin et al., 1982) appear to have discrete internal carotid arteries, and deserve study, but most are relatively small animals.
The taxonomy of this paper considers that primates comprise two clades, the Suborder Strepsirrhini that includes the lemuriform primates and the Suborder Haplorrhini that includes the tarsiers and simians (=anthropoids). Haplorrhine primate evolution followed a path that gave rise to four clades in succession: the Platyrrhini (New World monkeys) and the Catarrhini, which include the Cercopithecidae (Old World monkeys), the Hylobatidae (lesser apes and gibbons) and the Hominidae (great apes) (Finstermeier et al., 2013; Perelman et al., 2011). We use phylogenetically informed regression analysis to test how the total internal carotid blood flow, which is the sum of flows in both arteries, increases with brain volume, while accounting for the statistical dependence of related species due to their common ancestral history.
RESULTS
Brain volumes (Vbr) for Strepsirrhini, Haplorrhini and Diprotodontia increase with body mass (Mb) with exponents of 0.69, 0.71 and 0.70, respectively (Fig. 1; Table 1). Relative brain volume (body-mass-specific) does not remove the allometric effect of body mass, but is a common measurement in the literature. To obtain it, the equations in Table 1 remain the same, except that the exponent changes to a negative value equal to the given exponent minus 1.0. ANCOVA shows no significant differences in exponent between any of the three groups (P>0.88), but significant differences in elevation between all groups (P<0.001). Thus, for 1 kg animals, Vbr is 22.9 ml in Haplorrhini, 12.9 ml in Strepsirrhini and 5.3 ml in Diprotodontia.
Relationship between brain volume and body mass in primates and marsupials. Data represent 34 species of Haplorrhini, seven species of Strepsirrhini and 19 species of diprotodont marsupials on log–log axes. See Table 1 for equations and statistics for the regression lines.
Relationship between brain volume and body mass in primates and marsupials. Data represent 34 species of Haplorrhini, seven species of Strepsirrhini and 19 species of diprotodont marsupials on log–log axes. See Table 1 for equations and statistics for the regression lines.
Allometric equations of the form Y=aXb, relating body mass, brain volume, total flow through both internal carotid arteries and internal carotid foramen radius in two groups of primates and one group of marsupials

The radius of the internal carotid foramen is our primary measurement for the three groups (Fig. 2; Table 1). ANCOVA shows that exponents of foramen radius rICF on Vbr are all significantly different from each other (P<0.05). Johnson–Neyman (J–N) tests reveal no significant difference between the data for Diprotodontia and Haplorrhini within the range where Vbr overlap, but Strepsirrhini are all significantly lower than the other two groups.
Relationship between internal carotid artery foramen radius and brain volume. See Table 1 for equations and statistics.
Relationship between internal carotid artery foramen radius and brain volume. See Table 1 for equations and statistics.
Total internal carotid blood flow rate (Q̇ICA), which is the sum of flows in both arteries, reflects the relationship of foramen size on Vbr, but is not equivalent to it, because flow rate takes into account the scaling of shear stress with body mass (Fig. 3; Table 1). ANCOVA shows significant differences in exponents between the three groups (P<0.05). There is no significant difference in the data between Haplorrhini and Diprotodontia where the Vbr overlaps, but the smallest three points for Diprotodontia are significantly higher than predicted from Haplorrhini (J–N test). All points for Strepsirrhini are significantly lower than any of the other groups (J–N test).
Relationship between total internal carotid blood flow rate (both arteries) and brain volume. See Table 1 for equations and statistics.
Relationship between total internal carotid blood flow rate (both arteries) and brain volume. See Table 1 for equations and statistics.
Total internal carotid blood flow rate increases with Mb in all three groups (Fig. 4; Table 1). ANCOVA shows that the exponents of Q̇ICA on Mb are not significantly different between Haplorrhini and Diprotodontia (P=0.09), but Haplorrhini have significantly higher elevation (P<0.001). The exponent for Q̇ICA on Mb for Strepsirrhini is significantly higher than both the Haplorrhini and Diprotodontia, but the Q̇ICA data are all significantly lower than those of all Haplorrhini, and all but the lemur Propithecus diadema are lower than Diprotodontia according to the J–N test.
Relationship between total internal carotid blood flow rate (both arteries) and body mass. See Table 1 for equations and statistics.
Relationship between total internal carotid blood flow rate (both arteries) and body mass. See Table 1 for equations and statistics.
Focussing on the Haplorrhini, a tight relationship exists between Q̇ICA and Vbr, with an exponent of 0.95 (Fig. 5). The 95% confidence interval for the exponent is 0.84–1.06. A trend towards larger brains with higher Q̇ICA emerges throughout haplorrhine evolution, from the smaller brained New World monkeys to the large brained apes, including humans, with the highest Q̇ICA. Blood flow to the human brain is not exceptional and falls almost exactly on the regression line for haplorrhine primates in general.
Relationship between total internal carotid blood flow rate (both arteries) and brain volume for Haplorrhini on log–log axes. Symbols represent the Family Hominidae (great apes), Family Hylobatidae (lesser apes), Superfamily Cercopithecidae (Old World monkeys) and Parvorder Platyrrhini (New World monkeys). Regression line and 95% confidence belts for the regression mean are shown. Selected genera are indicated. See supplementary material Table S1 for the full list of species.
Relationship between total internal carotid blood flow rate (both arteries) and brain volume for Haplorrhini on log–log axes. Symbols represent the Family Hominidae (great apes), Family Hylobatidae (lesser apes), Superfamily Cercopithecidae (Old World monkeys) and Parvorder Platyrrhini (New World monkeys). Regression line and 95% confidence belts for the regression mean are shown. Selected genera are indicated. See supplementary material Table S1 for the full list of species.
For haplorrhine primates, the best supported phylogenetic model for explaining variability in log10-transformed Q̇ICA includes Vbr rather than Mb. This model was consistently the top-ranked model across all 1000 phylogenies (Table 2). The extent to which Q̇ICA and Vbr have evolved by Brownian motion (i.e. phylogenetic signal measured by Pagel's Lambda) was very strong among the species of Haplorrhini, indicating a high correlation between brain perfusion derived from the internal carotid arteries and relatedness of species during primate evolution. Across all 1000 phylogenetic trees, the median maximum-likelihood estimate of Pagel's Lambda in this model was 1.00. These results were also confirmed in the Diprodontia (Table 2).
DISCUSSION
The scaling of brain volume Vbr with body mass Mb in the two groups of primates and diprotodonts are parallel, with nearly identical exponents of about 0.70 (Fig. 1; Table 1). However, at any given Mb, haplorrhine primates have brains about 4-times larger, and strepsirrhines 2.5-times larger, than diprotodont marsupials. The implications of these large differences in Vbr are not clear, because the three groups show similarities in diet, locomotion and arboreal habitat, but it may be correlated with a relatively high degree of social interaction and cognitive ability in primates. Within primates, there is a correlation between brain size, innovation and social learning (Reader and Laland, 2002), which is not evident among marsupials.
Links between cognitive ability and brain volume or encephalization quotient have been difficult to sustain (Deaner et al., 2007; Henneberg, 1998; Herculano-Houzel, 2012), but cerebral metabolic rate and blood flow rate may be better, indeed more appropriate, correlates. It is clear that brain metabolism, inferred through internal carotid blood flow increases strongly with increasing brain size in haplorrhine primates with a relatively high scaling exponent, Vbr0.95, which is not significantly different from 1.0, but significantly greater than 0.75 (Fig. 5). The cerebral cortex of primates becomes increasingly convoluted in larger brains, and the surface area increases with an exponent much higher than 0.67 that is expected for geometric expansion with the same appearance (Chenn and Walsh, 2002; Luders et al., 2008; Radinsky, 1975). The anthropoid cerebral cortex is relatively thick and convoluted, with the degree of gyrification increasing significantly as brain volume increases (Zilles et al., 2013). In particular, the surface area and gyrification of the temporal lobes increases with a significantly greater exponent than predicted for geometric growth of the brain (Rilling and Seligman, 2002). Conversely, gyrification in diprotodont species does not vary considerably with increasing brain size (Zilles et al., 2013), and this is reflected in the scaling of blood flow in diprotodont marsupials, which is proportional to a much lower exponent, Vbr0.73 (Fig. 3; Table 1). Although the exponents differ between haplorrhine primates and diprotodonts, the absolute perfusion rates through the internal carotid arteries are similar, where brain volume overlaps (Fig. 3). If perfusion rate is related to cognitive ability, we conclude that there is little difference between the smaller (phylogenetically older) primates and marsupials.
The large difference between the exponent of Q̇ICA on brain volume between haplorrhine primates (0.95) and diprotodont marsupials (0.73) almost exactly straddles the exponent of 0.86 for brain metabolic rate (Karbowski, 2007) and 0.84 for brain perfusion (Karbowski, 2011). ANCOVA shows that neither 0.95 nor 0.73 are significantly different from 0.86, based on the samples and relationships considered here. One explanation for this result may be that the primates represented the high end, and non-primates the low end, of the range in Karbowski's analysis, which could result in an intermediate exponent. Indeed, there is further reason to believe that the scaling of brain metabolism of primates differs from non-primates. The scaling of neuronal density (number of neurons per milligram of brain) on brain mass gives a negative exponent of −0.37 in rodents, but −0.12 in primates (Herculano-Houzel, 2011a). Converting from neuronal density to total neuron number, the exponents in relation to brain mass become 0.63 and 0.88, respectively. Because metabolic rate is related to neuron number, it is clear that primates have a much steeper scaling than rodents.
The significant differences in exponents between primates and diprotodonts may be related to the scaling of neuron density in the two groups. In primates, neuron and non-neuron cell numbers increase linearly with brain volume, maintaining neuron and non-neuron cell densities and cell sizes (Herculano-Houzel, 2011b). Conversely, in diprotodonts, an increase in brain volume results in an increase in neuron size and a decrease in neuron density (Herculano-Houzel et al., 2006; Karbowski, 2007). Glucose and oxygen consumption per neuron is reported to be independent of neuron size (Herculano-Houzel, 2011c; Karbowski, 2011), so the diprotodonts scale lower, with an exponent close to the expected 0.75 value of ‘Kleiber's Law’. It would be informative to determine whether or not other non-primate groups scale similar to diprotodonts.
The data for Strepsirrhini are significantly lower than both of the other groups because the internal carotid arteries are not the main supply to the forebrain. The pattern of brain perfusion in lemurs is different from other primates (Schwartz and Tattersall, 1987). Dwarf lemurs and mouse lemurs (Cheirogaleidae) possess a highly reduced internal carotid artery and rely on the ascending pharyngeal artery. Sportive lemurs (Lepilemuridae) primarily rely on external carotid and vertebral arteries, and have reduced internal carotid, stapedial and ascending pharyngeal arteries. Three other lemur families (Lemuridae, Indriidae, Daubentoniidae) rely on large stapedial arteries and reduced internal carotids. The carotid foramen is therefore a poor indicator of brain perfusion in strepsirrhine primates.
A perplexing feature of the data is that calculated brain perfusion in diprotodont marsupial brains is not significantly different from haplorrhine brains within the range of overlapping brain volume (Fig. 3). This is in stark contrast with measurements of lower cellular and neuronal densities of the South American marsupial, the short-tailed opossum Monodelphis domestica, and estimates that the energy consumption of the neocortex of marsupials is less than 10% of that in similar-sized eutherians (Seelke et al., 2014). At this stage, we cannot provide a satisfactory explanation for why brain perfusion appears relatively high, particularly in the small diprotodonts; however, it could relate to differences in O2 extraction.
The evolution of primate brain perfusion
The observed allometric relationship between brain perfusion and brain volume in primates increased from the small-bodied, small-brained New World monkeys to the larger Old World monkeys, and finally to the great apes and humans (Fig. 5) (Martin, 2001). According to mitochondrial genomes, the Strepsirrhini (lemuriform primates) gave rise to the Haplorrhini (tarsiers and anthropoid Haplorrhini) about 66–69 Mya, the Platyrrhini (New World monkeys) split from the Catarrhini (Old World monkeys and apes) about 46 Mya, the Hominoidea (apes) arose about 32 Mya and the Hylobatidae (lesser apes and gibbons) split from the Hominidae (great apes) about 19 Mya (Finstermeier et al., 2013). Our phylogenetically informed analysis reveals a higher than expected blood flow relative to brain size as these major primate groups evolved. This may reflect an increase in cognitive activity and relative intelligence in primates, and the evolutionary correlates behind this trajectory may involve social behaviour, pair-bonding, dietary changes, among others (Dunbar and Shultz, 2007a,,b). To facilitate the increasing brain size, an increase in energy intake or decrease in energy expenditure on other organ development may have occurred (Aiello and Wheeler, 1995). In more advanced Haplorrhini, including human and chimpanzee, a dietary change from herbivorous to omnivorous resulted in consumption of more energy-dense foods. This shift may have not only increased energy intake without the need to increase foraging time, but also decreased the energy needed for digestive organs (Herculano-Houzel, 2012; Parker, 1996). The invention of cooking, which essentially softens and pre-digests food outside the digestive cavity may have further reduced energy requirements of digestion, allowing for a reduction in the size of the gastrointestinal tract (Herculano-Houzel, 2012). Interestingly, Q̇ICA in humans is exactly as predicted from the other primates for a 1.2 kg brain, indicating that humans are not exceptional in primate evolution other than having a large brain relative to body size.
Critique of methods
It is possible to provide three independent validations of actual blood flow rate and carotid foramen size in humans, rats and mice, because sufficient data are available (Table 3). Here, we arrive at internal carotid blood flow rate in two ways: first, by obtaining actual measured blood flow rates from the literature (L), and second, by calculating flow rate from foramen radius according to the methods of the present (P) study. The L/P ratios of Q̇ICA are 1.28 for humans, 1.00 for rats and 1.35 for mice. The L/P ratios of lumen radius are 1.09, 1.00 and 1.10, respectively. Given that the shear stress equation includes lumen radius to the third power, the agreement between measured and calculated flow rates is fairly good.
Measured and calculated cardiovascular variables relative to blood flow through the internal carotid arteries of humans and rodents

We calculate the rate of blood flow through the carotid foramen according to the shear stress equation, Q̇ = (τπr3)/(4η), in which flow is proportional to the radius cubed. This equation does not require knowledge of blood pressure or vessel length. The equation is less common among circulation physiologists than the Poiseuille–Hagen (P–H) equation, Q̇ = (ΔPπr4)/(8Lη), in which flow rate is proportional to the pressure decrease (ΔP) along the vessel and to r4/L, where L is the length of the vessel. Not only is L difficult to measure for the internal carotid, but ΔP must also be extremely small and most likely never measured. In any case, if r and L scale with body size geometrically (∝Mb0.33), r3 and r4/L are allometrically equivalent.
We assume that brain case volume scales similarly with brain volume and we include the volume of the cerebellum that is perfused largely by vertebral arteries. Although the volume of the cerebellum in relation to the entire brain increases slightly in larger primates, from approximately 9% in small monkeys to 13% in great apes (Rilling and Insel, 1998), the difference is negligible compared with other uncertainties in the study.
Conclusions
Because the metabolic rate of the brain is related to its blood supply, the size of the carotid foramina can be used as an index of brain metabolism and possibly cognitive ability. The evolutionary lineage of Haplorrhini shows a relatively steep progression of increasing brain perfusion, compared with the diprotodont marsupials. The trajectory of primate brain evolution appears to have been determined before the great apes arrived, but they sit near the line and humans on it. Although humans have a large brain in relation to their body mass, brain perfusion in them appears to be simply another step in primate evolution and unexceptional. Whether primates differ from mammals in general is not known, because we are unable to gauge brain perfusion in other groups with the technique of this study.
Our approach to inferring metabolic rates from blood flow and foramen size is not only expedient, but also permits measurements that are not otherwise practical. In the present case, direct determination of the metabolic rate of the brain is technically difficult. However, measuring the carotid foramen is relatively simple and can be used for both extant and extinct species. The size of the carotid foramen in fossilized Sahelanthropus, Australopithecus and extinct Homo may reveal changes in blood flow to the brain that can indicate a trend towards higher-order thinking. A separate analysis for hominids could determine if the advent of cooking did in fact result in an increase in brain volume, and a subsequent increase in carotid foramen area due to increased brain metabolism.
MATERIALS AND METHODS
Morphometric data, including brain case volume and carotid foramen radii, were collected from 34 species of simian primates (Suborder Haplorrhini), 7 species of prosimian primates (Suborder Strepsirrhini) and, for comparison, 19 species of diprotodont marsupials (Order Diprotodontia).
Braincase volume (Vbr) was measured in dried and degreased skulls by filling them with either millet seed (skulls <200 ml) or sorghum seed (>200 ml). After sealing major foramina with cotton wool, skulls were filled to just below the opening of the foramen magnum, and were uniformly agitated to ensure similar packing characteristics. Seed was then poured into an appropriately sized graduated cylinder and agitated to determine volume. Brain case volume is a reliable estimate of brain mass, because the brain has a density of 1.04 g ml−1 (Isler et al., 2008).
We identified the carotid foramina from published anatomical illustrations of primate skulls (Ankel-Simons, 2007) and marsupial skulls (Reig et al., 1987; Travouillon et al., 2014) (Fig. 6). For primates, we measured the external opening of the carotid canal, not the foramen lacerum, which in living humans is sealed by cartilage on the outside of the skull but opens on the inside to pass the internal carotid artery (Tauber et al., 1999).
Diagrams of representative skulls of a diprotodont marsupial and a haplorrhine primate, showing the locations of major foramina. The radius of the external opening of the carotid canal is the focus of this study and is easily measured from a ventral view. The carotid canal of the diprotodont enters obliquely and is best viewed from a posterior–ventral–lateral perspective. Note the foramen lacerum of the primate skull, which is plugged with cartilage on the ventral surface, but transmits the carotid artery from the carotid canal into the brain case in life.
Diagrams of representative skulls of a diprotodont marsupial and a haplorrhine primate, showing the locations of major foramina. The radius of the external opening of the carotid canal is the focus of this study and is easily measured from a ventral view. The carotid canal of the diprotodont enters obliquely and is best viewed from a posterior–ventral–lateral perspective. Note the foramen lacerum of the primate skull, which is plugged with cartilage on the ventral surface, but transmits the carotid artery from the carotid canal into the brain case in life.
The minimum diameter of the carotid foramen, below the flared surface, was measured with callipers directly from large skulls. For small skulls, calibrated digital photographs were taken with a camera (LC20, Olympus Soft Imaging Solutions GmbH, Münster, Germany) mounted on a dissecting microscope (SZ2-ILST, Olympus Corporation, Tokyo, Japan) and foramen area was measured with image analysis software (ImageJ, www.nih.gov) (Fig. 7). Minimum diameter or area was converted to radius, which was the primary measurement that led to calculation of blood flow rate Q̇ (cm3 s−1) according to the shear stress equation: Q̇ = (τπr3)/(4η). The equation is independent of blood pressure, and blood viscosity (η) is taken as a constant (0.04 dyne s cm−1) in the large vessels of mammals (Schmid-Schönbein et al., 1969). The other two variables, vessel lumen radius (r; cm), and shear stress (τ; dynes cm−2) were evaluated as follows.
Examples of internal carotid foramina of Haplorrhini. (A) Homo sapiens (Vbr=1186 ml). (B) Pan troglodytes (Vbr=475 ml). (C) Gorilla gorilla (Vbr=461 ml). (D) Hylobates lar (Vbr=103 ml). (E) Papio hamadryas (Vbr=156.5 ml). (F) Saimiri sciureus (Vbr=20 ml).
Lumen radius (r)
The shear stress equation requires the radius of the lumen of the internal carotid artery, not the radius of the carotid foramen. The internal carotid artery is the only significant blood vessel to pass through the carotid foramen, and it largely fills it. There is a thin layer of fascia, a few small veins, and the sympathetic internal carotid nerve plexus that passes as tiny fibres around and within the vessel wall (Drake et al., 2015). The internal and external radii of the arterial wall almost certainly scale in parallel with one another owing to the effect of the principle of Laplace for cylinders in requiring wall thickness to be directly proportional to lumen radius if blood pressure is constant (Burton, 1965). In fact, arterial blood pressure is practically invariant over the body mass range of primates (White and Seymour, 2014). The ratio of wall thickness (including intima, media and adventitia) to lumen radius is 0.40 in human common carotid arteries (Skilton et al., 2011), 0.40 in dog internal carotid arteries (Orsi et al., 2006) and 0.36 in rat brain posterior communicating arteries (Gules et al., 2002). Thus, we use the conversion: ICA lumen radius=(carotid foramen radius/1.4).
Shear stress (τ)
Although it is widely assumed that the circulatory system is constructed to normalize shear stress at a constant level (∼10–20 dynes cm−2), regardless of the section of the arterial tree or the mass of the animal, it is now clear that shear stress in major arteries decreases in larger mammals. Based on aortic flow rates measured with magnetic resonance angiography in mice, rats and humans, the scaling exponent is −0.38 (Greve et al., 2006). On theoretical grounds, other analyses independently predict an exponent of −0.375 for shear stress scaling in mammalian aortas, but indicate that this value is not necessarily universal for all arteries (Weinberg and Ethier, 2007). We estimated the scaling of shear stress for the internal carotid artery from measured data from 70 kg humans (this study) and 375 g rats (Cai et al., 2012) using the shear stress equation and literature values. Blood flow rate in one internal carotid artery is taken as 3.83 cm3 s−1 in humans (Boysen et al., 1970) and 0.028 cm3 s−1 in rats (mean of four studies: Nakao et al., 2001; Ohno et al., 1979; Sakurada et al., 1978; Takagi et al., 1987). ICA lumen radius is taken as 0.221 cm in humans (mean of three studies: Kane et al., 1996; Krejza et al., 2006; Porsche et al., 2002) and 0.0304 in rats (Cai et al., 2012). Blood viscosity is considered invariant at 0.04 dyne s cm−1 (Schmid-Schönbein et al., 1969). The human and rat data produce a power equation, τ=167Mb−0.20, which we use to calculate ICA flow rate in this study. The exponent is less negative than that for aortas, but is similar to that for the common carotid artery (−0.21) derived from four species of mammal (Weinberg and Ethier, 2007).
Statistics
Internal carotid artery blood flow rate was calculated from the mean foramen size and then doubled to obtain total cerebral perfusion rate, Q̇ICA. A common Q̇ICA datum for each species was obtained as a mean from one to ten adult specimens. Adult body mass of each species was obtained from the literature. All data were analysed allometrically using the power equation Y=aXb, where Y is the variable of interest, X is the body mass or brain volume of the animal, a is the scaling factor (elevation of the curve) and b is the scaling exponent. The exponent indicates the shape of the relationship: Y is proportional to X only if b=1.0; if 1>b>0, the arithmetic curve of Y on X increases with an ever-decreasing slope and the ratio of Y/X decreases; if b>1.0, the arithmetic curve increases with an ever-increasing slope and the ratio of Y/X increases. The power equation is often log10-transformed to a linear one: log Y=log a+(b log X). In this case b is the slope of the line. Ordinary least squares regressions on log-transformed data, 95% confidence intervals of regression means and significant differences between orders were determined using analysis of covariance (ANCOVA) in GraphPad Prism 5 statistical software (GraphPad Software, La Jolla, CA, USA) (Zar, 1998). When ANCOVA revealed differences in the exponent b, a Johnson–Neyman (J–N) test was performed to determine the range of X values over which the individual data for Y are significantly different (White, 2003).
Phylogenetic generalized least-square (PGLS) models were constructed to test the association between log10-transformed values of Q̇ICA, Mb and Vbr in the haplorrhine primates and diprodont marsupials separately. PGLS models were constructed using the ‘pgls’ function in the package ‘caper’ (Orme et al., 2011) in the R software environment for statistical computing and graphics R Core Team (RCoreTeam, 2013). The ‘pgls’ function incorporated the covariance between related taxa into the calculation of estimated coefficients from a generalized least squares model. The covariance matrix of the expected covariance between each pair of tips was calculated using the branch lengths of an estimated phylogeny.
Haplorrhine primate phylogenetic trees were constructed online (http://10ktrees.nunn-lab.org/primates/) using the ‘10kTrees dated phylogeny v3’ with the associated taxonomy from GenBank (Arnold et al., 2010). One thousand phylogenetic trees were constructed, and PGLS models were assessed using Akaike Information Criterion (AICc) corrected for small sample sizes (Burnham and Anderson, 2001). A single Diprodont phylogenetic tree was obtained from a species-level phylogenetic supertree of marsupials (Cardillo et al., 2004). No phylogenetic analysis was done on the Strepsirrhini as the number of species was insufficient. Phylogenetic signal in blood flow index was measured by Pagel's lambda (λ) (Pagel, 1999). Pagel's lambda indicates the strength of the phylogenetic relationship, where values lie between 0 and 1. Lambda values of or near 0 are indicative of phylogenetic independence and values of or near 1 indicate that the variable is related to evolutionary history (Kamilar and Cooper, 2013).
Acknowledgements
We appreciate the efforts of many anonymous referees that improved our paper, especially the primatologists who instructed us to reanalyse our results and communicate them more clearly. We thank Cath Kemper and David Stemmer of the South Australian Museum for access to the primate collection and facilities, Maciej Henneberg for providing advice on the evolution of primate brains and access to his collections, Chris Macgowan for kindly providing original data on internal carotid flow in mice and Qiaohui Hu for help with collecting data and literature and for drawing Fig. 6.
Footnotes
Author contributions
R.S.S., S.E.A. and E.P.S. designed and performed research and wrote the paper; P.C. performed the phylogenetic analysis and helped write the paper.
Funding
The project was funded by the Australian Research Council Discovery Grant (DP-120102081) to R.S.S. P.C. is an ARC Future Fellow (FT0991420).
References
Competing interests
The authors declare no competing or financial interests.