# Brain evolution in Proboscidea (Mammalia, Afrotheria) across the Cenozoic

Julien Benoit, Lucas Legendre, Rodolphe Tabuce, Théodor Obada, Vladislav Mararescul, Paul R. Manger

► To cite this version:

Julien Benoit, Lucas Legendre, Rodolphe Tabuce, Théodor Obada, Vladislav Mararescul, et al.. Brain evolution in Proboscidea (Mammalia, Afrotheria) across the Cenozoic. Scientific Reports, 2019, 9 (1), 10.1038/s41598-019-45888-4. http://02408988

HAL Id: hal-02408988

https://hal.science/hal-02408988

Submitted on 27 May 2021

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

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

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

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC REPORTS

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

OPEN

# Brain evolution in Proboscidea (Mammalia, Afrotheria) across the Cenozoic

Received: 31 January 2019

Accepted: 7 June 2019

Published online: 27 June 2019

Julien Benoit<fnref n="1" />, Lucas J. Legendre<fnref n="1" />,², Rodolphe Tabuce³, Theodor Obada⁴, Vladislav Mararescul⁴ & Paul Manger⁵

As the largest and among the most behaviourally complex extant terrestrial mammals, proboscideans (elephants and their extinct relatives) are iconic representatives of the modern megafauna. The timing of the evolution of large brain size and above average encephalization quotient remains poorly understood due to the paucity of described endocranial casts. Here we created the most complete dataset on probosccidean endocranial capacity and analysed it using phylogenetic comparative methods and ancestral character states reconstruction using maximum likelihood. Our analyses support that, in general, brain size and body mass co-evolved in proboscideans across the Cenozoic; however, this pattern appears disrupted by two instances of specific increases in relative brain size in the late Oligocene and early Miocene. These increases in encephalization quotients seem to correspond to intervals of important climatic, environmental and faunal changes in Africa that may have positively selected for larger brain size or body mass.

In mammals, brains larger than 700 g have evolved infrequently, being found primarily in extant humans, cete-
ceans and elephants<fnref n="1" />. The temporal/evolutionary history of human brain size and shape has been studied in
detail²⁻⁵, with the 700 g barrier on brain size evolution being surpassed with the appearance of *Homo ergaster/erec-
tus*, approximately 1.8 million years ago<fnref n="1" />. The suborder Cetacea has many species that also exhibit large brains,
both absolute and relative to body mass⁶,⁷. It has been established that large cetacean brains emerged with the
evolution of the Neoceti approximately 32 million years ago⁸,⁹. In contrast to humans and cetaceans, the temporal
and evolutionary history of brain size in Proboscidea is poorly understood<fnref n="1" />,¹⁰.

Elephants are the largest and among the most behaviourally complex terrestrial mammals. They are known for their extensive long-term memory, problem-solving abilities, behavioural adaptability, their ability to recognize themselves in a mirror, to manipulate their environment, and to manufacture tools with their trunk<sup>11,12</sup>. On the one hand, given that the metabolic cost of brain tissue is high<sup>13</sup>, and following the assumption that natural selection would not maintain a costly organ that brings no benefit, the presence of a large brain in elephants suggests there is some adaptive value, potentially related to their cognitive capacities. It has been proposed that evolutionary selective pressures requiring increased intelligence (e.g. cognitive and sensory processing, memory, behavioural flexibility) to cope with a variety of environmental variables (e.g. sociality, gregariousness, diet, habitat) has driven the evolution of large brains in both Primates and 'ungulates'<sup>14-17</sup>. This is based on the retrospective assignation of implied cognitive needs centred on the apparent behavioural flexibility and intelligence of extant species. Similarly, this retrospectively applied cognitive explanatory paradigm has been argued to apply to the evolution of a large brain in cetaceans<sup>18-20</sup>, and elephants<sup>10,12</sup>; however, alternative explanations of brain size evolution have been forwarded, focussing on features such as the structural laws of form linking body and brain mass over evolutionary time<sup>1,7,21,22</sup>. Under this assumption, the elephantine brain is a typical mammalian brain that appears "to scale consistently with the scaling laws governing brain and body mass relationships across mammals"<fnref n="1" />. In parallel, the assumption that complex cognition requires an increase in brain size has been heavily

<footnote n="1">¹Evolutionary Studies Institute (ESI), University of the Witwatersrand, Braamfontein, 2050, Johannesburg, South Africa. ²Jackson School of Geosciences, The University of Texas at Austin, 2275 Speedway Stop C9000, Austin, TX, United States. ³Institut des Sciences de L'Evolution de Montpellier, Université Montpellier 2, Place Eugène Batillon, F-34095 Montpellier, cedex 05, Montpellier, France. ⁴Academy of Sciences of Moldova, Institute of Zoology, Chișinău, Moldova. ⁵School of Anatomical Sciences, University of the Witwatersrand, 7 York Road, Parktown, 2193, Johannesburg, South Africa. Correspondence and requests for materials should be addressed to J.B. (email: **Julien.benoit@wits.ac.za**)</footnote>
<table><caption>**Table 1.** Endocast, calculated brain size and Manger's EQ⁷ for the proboscideans used in the analyses. The sirenian *Prorastomus* and the hyracoid *Seggeurius* were added for comparative purposes but were not included in the analyses. Sources of data: a¹⁰, b³⁶, c³⁴, d³⁷.</caption><thead><tr><th></th><th>Endocast volume (cm³)</th><th>Brain size (g)</th><th>Body mass (g)</th><th>EQ Manger</th></tr></thead><tbody><tr><td>Elephas maximus</td><td>4550<sup>b</sup></td><td>4181</td><td>2267430<sup>b</sup></td><td>1,83</td></tr><tr><td>Elephas maximus</td><td>5220<sup>b</sup></td><td>4798</td><td>3216000<sup>b</sup></td><td>1,62</td></tr><tr><td>Elephas maximus</td><td>5000<sup>b</sup></td><td>4595</td><td>3450400<sup>b</sup></td><td>1,48</td></tr><tr><td>Elephas maximus</td><td>6075<sup>b</sup></td><td>5584</td><td>3190098<sup>b</sup></td><td>1,90</td></tr><tr><td>Loxodonta africana</td><td>5712<sup>b</sup></td><td>5250</td><td>6654000<sup>b</sup></td><td>1,05</td></tr><tr><td>Loxodonta africana</td><td>9000<sup>b</sup></td><td>8274</td><td>4380000<sup>b</sup></td><td>2,23</td></tr><tr><td>Loxodonta africana</td><td>4000<sup>b</sup></td><td>3676</td><td>5174400<sup>b</sup></td><td>0,88</td></tr><tr><td>Loxodonta africana</td><td>4050<sup>b</sup></td><td>3722</td><td>1793300<sup>b</sup></td><td>1,93</td></tr><tr><td>Loxodonta africana</td><td>4420<sup>b</sup></td><td>4062</td><td>3505000<sup>b</sup></td><td>1,29</td></tr><tr><td>Loxodonta africana</td><td>5300<sup>b</sup></td><td>4871</td><td>5550000<sup>b</sup></td><td>1,11</td></tr><tr><td>Loxodonta africana</td><td>4480<sup>b</sup></td><td>4117</td><td>2750000<sup>b</sup></td><td>1,56</td></tr><tr><td>Loxodonta africana</td><td>4210<sup>b</sup></td><td>3869</td><td>4000000<sup>b</sup></td><td>1,12</td></tr><tr><td>Loxodonta africana</td><td>4100<sup>b</sup></td><td>3768</td><td>2160000<sup>b</sup></td><td>1,70</td></tr><tr><td>Loxodonta africana</td><td>4000<sup>b</sup></td><td>3676</td><td>2537000<sup>b</sup></td><td>1,48</td></tr><tr><td>Palaeoloxodon falconeri</td><td>1800<sup>a</sup></td><td>1652</td><td>168000<sup>a</sup></td><td>4,81</td></tr><tr><td>Palaeoloxodon antiquus</td><td>5446<sup>a</sup></td><td>5005</td><td>3649880<sup>a</sup></td><td>1,54</td></tr><tr><td>Palaeoloxodon antiquus</td><td>9000<sup>a</sup></td><td>8274</td><td>11000000<sup>c</sup></td><td>1,14</td></tr><tr><td>Mammuthus primigenius</td><td>4687<sup>a</sup></td><td>4307</td><td>6000000<sup>d</sup></td><td>0,92</td></tr><tr><td>Mammuthus meridionalis</td><td>5828<sup>a</sup></td><td>5357</td><td>11000000<sup>d</sup></td><td>0,74</td></tr><tr><td>Mammuthus columbi</td><td>6232<sup>a</sup></td><td>5728</td><td>9800000<sup>d</sup></td><td>0,86</td></tr><tr><td>Mammut americanum</td><td>3862<sup>a</sup></td><td>3549</td><td>6384056<sup>a</sup></td><td>0,73</td></tr><tr><td>Mammut americanum</td><td>4630<sup>a</sup></td><td>4255</td><td>8000000<sup>d</sup></td><td>0,74</td></tr><tr><td>Stegodon insignis</td><td>3838</td><td>3527</td><td>2000000<sup>d</sup></td><td>1,69</td></tr><tr><td>Zygolophodon borsoni</td><td>5133</td><td>4718</td><td>16000000<sup>d</sup></td><td>0,50</td></tr><tr><td>Palaeomastodon beadnelli</td><td>771</td><td>706</td><td>2500000<sup>d</sup></td><td>0,29</td></tr><tr><td>Moeritherium lyonsi</td><td>240<sup>a</sup></td><td>218</td><td>810000<sup>d</sup></td><td>0,20</td></tr><tr><td>Prorastomus sirenoides</td><td>87<sup>b</sup></td><td>90</td><td>98156<sup>b</sup></td><td>0,39</td></tr><tr><td>Seggeurius amourensis</td><td>5<sup>b</sup></td><td>5</td><td>2932<sup>b</sup></td><td>0,29</td></tr></tbody></table>

criticized¹,²,¹,²,³,²⁴. Therefore whether the evolution of a large brain is an adaptive trait (e.g.)¹⁰, or simply accompanied the increase in body mass in proboscideans (e.g.)¹, remains an open question.

While proboscidean evolution is well documented and understood²⁵⁻³⁰, studies of the brains and endocranial casts of extinct proboscidean genera are limited¹⁰,¹²,³¹⁻³⁵. These studies have indicated that brain size, both in absolute and relative terms, appears to have increased dramatically in Elephantiformes, sometime between the Oligocene and Pliocene, a poorly constrained interval of more than 30 million years¹⁰. Thus, the timing and rate of changes in brain size, and the taxa involved in the evolution of a large brain during proboscidean evolution remains unclear. Building on the database of Benoit¹⁰, by adding data from three key taxa, *Zygolophodon (Mammut) borsoni*, *Palaeomastodon beadnelli* and *Stegodon insignis*, we apply statistical and phylogenetic comparative methodology to changes in brain size in proboscideans throughout their evolutionary history. To address whether a large brain in proboscideans evolved as an adaptation or as a consequence of increased body mass, we examined the temporal scaling relationships of brain with body size and with associated environmental changes that have occurred during proboscidean evolutionary history.

## Methods

**Data collection.** Data analysed (endocranial volume and body mass) are summarized in Table 1. Data for extant species of elephants are from Benoit *et al.*³⁶ (note that this dataset contained two duplicated lines that are removed here) and those for extinct species are from Benoit¹⁰. The body masses of some extinct taxa were updated using the recent estimates of Larramendi and Palombo³² and Larramendi³⁷ (see Table 1).

This dataset was also updated with the addition of endocranial volume and body mass of three taxa: *Zygolophodon (Mammut) borsoni* (MCFFM-CLB-1), from the Pliocene (Mammal Neogene zone 15) of Otman Hill (Moldova)³⁸, *Stegodon insignis* (MNHN-A952) from an unknown Miocene locality in the Himalayas³⁹, and *Palaeomastodon beadnelli* (NHMUK PV M 8464) from the early Oligocene of the Fayum deposits (Egypt)⁴⁰. As one of the largest proboscideans that existed³⁷, *Z. borsoni* is important for discussing the effect of body mass on the evolution of brain mass in proboscideans. Stegodontids may be the possible outgroup of modern elephants, and *Palaeomastodon beadnelli* is the basal-most known elephantiform and the closest relative of the clade Elephantimorpha in the dataset²⁶,²⁸⁻³⁰,⁴¹. For this study, the braincase of *P. beadnelli* (NHMUK PV M 8464), which is partially exposed⁴⁰, was carefully examined and drawn by a palaeo-artist (Sophie Vrard) in dorsal and ventral
views under the supervision of J.B. (Supplementary Fig. 1). The endocranial volume was estimated to 771 cm³ using double-graphic integration on the drawings³¹. Radinsky (ref. ⁴², p. 84) performed comparisons with the water displacement technique that validated the accuracy of double graphic integration and suggested that this method should be preferred in case of incomplete preservation or when direct measurement is impossible, as is the case here with NHMUK PV M 8464. As originally described, specimen MNHN-A952 is the cast of the right side of the braincase of a *Stegodon insignis* skull assembled with the left endocast of a modern elephant to enable Gervais³⁹ to reconstruct a complete endocranial cast (Supplementary Fig. 2). The endocranial size of *Zygolophodon (Mammut) borsoni* is based on a latex endocranial cast of specimen MCFFM-CLB-1 (Fossil Fauna Complex Museum of Moldova, Chisinau) (Supplementary Fig. 3). To measure its volume, the endocast was digitized using photogrammetry AgiSoft PhotoScan Professional Version 1.3.2), a method that has proven as accurate as CT scanning⁴³ and sometimes more accurate than laser-scanning⁴⁴,⁴⁵. The volume of the resulting mesh was then measured using Avizo (FEI VSG, Hillsboro OR, USA). The resulting volume is 5133 cm³ (the *Zygolophodon* endocast was also measured at 4910 cm³ using water displacement, but the digital measure was preferred as it less prone to errors)². To measure the volume of the endocast of *Stegodon insignis*, the complete endocast was digitized using photogrammetry, its right half isolated digitally, and the volume of the resulting mesh was measured using the same software programs. The resulting volume, 1919 cm³, was then multiplied by two to obtain the volume of a complete endocast. Body masses for these three taxa are from Larramendi³⁷. The encephalization quotient (EQ) used is that of Manger⁷, which is very similar to Eisenberg's EQ⁴⁶, although it excludes outliers such as cetaceans and primates. Brain volumes used to calculate EQ are based on endocranial volumes corrected for meningeal thickness using Benoit's¹⁰ method. This is a necessary step as the thick meninges in proboscideans occupy some 14% of the endocranial cavity with non-neural tissue on average¹⁰. The resulting brain volumes were then multiplied by 1.036, which represents the average specific gravity of brain tissue³²,⁴⁷.

**Ancestral character state reconstruction.** In order to analyse temporal patterns of endocranial volume, body mass and EQ changes across proboscidean phylogenetic history, an ancestral character state reconstruction using maximum likelihood (ACSRML) was performed using the function fastAnc in the R package 'phytools'⁴⁸ (excluding the sirenian *Prorastomus* and the hyracoid *Seggeurius*). The phylogeny used for the ACSRML was assembled following the current consensus²⁸⁻³⁰,⁴¹ (Supplementary Fig. 4). The resulting dates are consistent with the recently published palaeogenomic study by Palkopoulou *et al.*⁴⁹. There have been many discrepancies concerning the phylogenetic placement of the genus *Palaeoloxodon*. Palaeogenomic studies have consistently placed *P. antiquus* within the African elephant genus *Loxodonta*⁴⁹,⁵⁰, whereas studies of the mitochondrial DNA have suggested that *P. falconeri* belongs to the Asian elephant genus *Elephas*⁵¹,⁵². Given the morphological support for the validity of the genus *Palaeoloxodon*²⁵,⁵³,⁵⁴ and the burden of evidence in favour of a direct ancestor to descendant relationship between *P. antiquus* and *P. falconeri* respectively²⁹,⁵⁵,⁵⁶, we preferred to remain conservative and preserved the monophyly of the genus *Palaeoloxodon*. Nevertheless, the branch leading to *Palaeoloxodon* is in an unresolved position within Elephantida in order to reflect these phylogenetic inconsistencies (Supplementary Fig. 4). Divergence ages and branch lengths were determined based on published literature²⁹,⁵⁷.

**Phylogenetic regressions.** In order to test for the presence of a functional relationship between brain mass and body mass, linear regressions were performed on our datasets (excluding the sirenian *Prorastomus* and the hyracoid *Seggeurius*). In such statistical analyses, there are several assumptions made about the data, such as normality, homoscedasticity and independence of observations. When comparing traits of different species, independence is considered potentially violated because of the evolutionary relationships of species to each other, and thus phylogeny needs to be considered as part of the analysis⁵⁸. All analyses were performed in R v. 3.5.3 (R Core Team, 2019).

Phylogenetic generalized least squares (PGLS)⁵⁹ are a very well-known and frequently used phylogenetic comparative method (see Garamszegi⁶⁰, for an extensive review), designed to consider the influence of evolutionary relationships between species (i.e., phylogenetic signal *sensu* Blomberg and Garland)⁶¹ in a generalized linear model, using a regression of one or more continuous predictors over a response variable. After a log-conversion to better approximate the allometric relationship between variables of interest⁶²,⁶³, endocranial volume was regressed on body mass using a calibrated phylogeny of Proboscidea in the dataset (Supplementary Fig. 4; Supplementary Material 5). Since the evolutionary model followed by this regression is not known, the analysis was performed several times using different evolutionary models, and the goodness of fit of each model determined from their AICc score, i.e., the Akaike Information Criterion corrected for finite sample size⁶⁴. Using the R package 'AICcmodavg'⁶⁵, AICc scores were compared for five different evolutionary models: Brownian motion, i.e., a purely neutral model of evolution⁶⁶; Ornstein-Uhlenbeck, i.e., a model with selective optima along the branches⁶⁷; lambda model, i.e. modified Brownian motion model with all branch lengths weighted using Pagel's lambda, a common estimator of phylogenetic signal⁶⁸, here estimated from the dataset through maximum likelihood; Early Burst, i.e., model with high initial evolutionary rates that decelerate through time after a cladogenesis event⁶⁹; and Ordinary Least Squares (i.e. simple, non-phylogenetic least squares regression). The model with the lowest AICc score was found to be the lambda model (AICc = -18.2, vs a range of -6.8-10.11 for other evolutionary models), with λ = 0.90, indicating a strong phylogenetic signal in the regression. Thus, the lambda model was selected to perform further analyzes on the dataset. PGLS were performed using the function gls in the R package 'nlme'⁷⁰, and the phylogenetic correlation structure corPagel (for the lambda model) specified in the package 'ape'⁷¹.

Normality and homoscedasticity of the residuals in the model were assessed using a Lilliefors test in the package 'nortest'⁷², as well as graphically checked through plotting against corresponding fitted values, to conform to the parametric assumptions of PGLS⁷³. Residuals were found to be non-normal, and seven outliers were graphically detected and removed from the dataset (Supplementary Material 6). These outliers consist of *Moeritherium*
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Regressions of log brain mass over log body mass corrected for phylogeny excluding outliers.</figcaption></figure>

lyonsi, *Palaeomastodon beadnelli*, *Palaeoloxodon falconeri*, one specimen of *Palaeoloxodon antiquus* as well as one specimen of *Elephas maximus* and two of *Loxondonta africana*. After pruning the corresponding branches in the tree, regression analyses were performed again, and the new residuals were found to be normal and homoscedastic (Fig. 1), allowing for a straightforward interpretation of the new results.

## Results

According to the ACSRML, which was performed on the whole dataset (Fig. 2), the EQ in the last common ancestor (LCA) of *Meretherium* and *Elephantiformes* is 0.24. It increases slightly in the next node to reach 0.3 in the LCA of *Elephantiformes* (*Palaeomastodon* and more derived proboscideans). The reconstructed brain and body masses at this node are about 1.8 times larger than in the previous node, while the EQ remains similar, which indicates that during the Paleogene, increases in brain mass were proportional to increases in body mass in basal proboscideans (Fig. 2). The reconstructed EQ doubles in the node corresponding to the common ancestor of *Elephantimorpha* that lived some 26 Ma ago in the late Oligocene, where it reaches a value of 0.73. At this node, the ACSRML indicates that body mass has increased by 2.3 times while brain mass has increased by 3.6 times compared to the previous node (Fig. 2). The reconstructed EQ then decreases slightly in the lineage leading to *Mammutidae* and *Mammut* (0.64 and 0.68 respectively) in which the reconstructed body mass exceeds 8 tons while brain mass gains only 1 kg (Fig. 2). In contrast, the EQ is nearly doubled again in the LCA of the *Elephantoidea* that lived some 20 Ma ago in the Early Miocene, to reach 1.09. At this node, the reconstructed body mass decreases slightly while brain mass is reconstructed as being about 1.2 times larger than in the previous node (Fig. 2). Reconstructed ancestral EQs, and brain and body masses then remain relatively stable across *Elephantoidea* and *Elephantidae* phylogeny (Fig. 2). Reconstructed ancestral body masses vary between about 4 and 7 tons, while brain masses vary between 4.5 and 6.4 kg in *Elephantidae* (Fig. 2). Reconstructed variations in the EQ nevertheless appear mostly due to variations in body mass. In the *Mammuthus* clade, the elevated reconstructed body mass (9.1 tons) coincides with a decrease in the reconstructed EQ compared to the ancestral value in the LCA of *Elephantoidea* to 0.89 (Fig. 2). In contrast, the increase in EQ to 1.68 in the *Elephas* clade is accompanied by a decrease in ancestral body mass to 3.2 tons (Fig. 2). The ancestral value for the *Palaeoloxodon* clade is surprisingly close to that for the LCA of *Elephantidae* (1.51).

The EQ varies around 1.47 ± 0.76 among all *Elephantimorpha* except in two taxa that display extreme EQ values. On one end of the spectrum is the small bodied and highly encephalized dwarf elephant of Sicily, *Palaeoloxodon falconeri*, which displays an EQ of 4.81 (Table 1). This species is well studied and its large EQ has been extensively discussed in the literature<sup>10,32,34</sup>. Its body size is between one tenth and one hundredth that of its close relative, *Palaeoloxodon antiquus*, whereas its endocast volume is about one half to one quarter that of *P. antiquus*<sup>10,32,34</sup> (Table 1). On the opposite end of the spectrum is *Zygolophodon (Mammut) borsoni* which, with an estimated body mass of 16 tons<sup>37</sup>, is the heaviest species in the dataset. The *Zygolophodon* EQ of 0.50 also makes it the least encephalized of all *Elephantimorpha*, although the clade to which it belongs, the *Mammutidae*, ancestrally display a low EQ value (0.64) similar to that reconstructed for the ancestral *Elephantimorpha* (Fig. 2).

The regression of brain mass over body mass corrected for phylogeny and excluding outliers (Fig. 1) shows that the relationship between the two variables is significant (p-value = 0.0004), but their correlation is relatively weak (R² = 0.50) and the slope of the regression line is low (0.20). This indicates that brain and body mass are significantly correlated, but this correlation and the effect of the brain and body mass on each other are moderate.

When outliers are included (Supplementary Material 6), the correlation between the two variables is slightly higher (R² = 0.61) and remains significant (P-value < 0.0001); however, since the assumptions of normality and homoscedasticity of the sample are not met for this analysis, this result is of dubious reliability<sup>73</sup>.
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Results of the Ancestral Character State Reconstruction using Maximum Likelihood. Geological time scale and δ18O curve after Zachos *et al.*<sup>84</sup>. Abbreviations: BoM, body mass; BrM, brain mass; EQ, encephalization quotient.</figcaption></figure>

Discussion

The oldest known endocranial cast of a proboscidean belongs to *Moeritherium lyonsi* from the late Eocene Qasr el Sagha Formation of the Fayum, Egypt<sup>40</sup>. The EQ of *Moeritherium* equals 0.2 (Table 1), which is an order of magnitude smaller than in any Elephantimorpha (Table 1). A comparably low EQ value is also found in *Palaeomastodon* (EQ = 0.29), which was contemporaneous with *Moeritherium*, and in the paenungulates *Seggeurius* (EQ = 0.29) and *Prorastomus* (EQ = 0.39), the basal-most hyracoid and sirenian respectively, and the closest relatives of proboscideans that lived during the early Eocene (~50 Ma)<sup>36</sup> (Table 1). As a consequence, it is not surprising that a rather small EQ is reconstructed as primitive for Proboscidea and Elephantiformes by the ancestral character state analysis (Fig. 2). The relatively small brain cavity of *Phosphatherium*, a stem proboscidean from the early Eocene (−55 Ma)<sup>74</sup> supports this assumption.

The analysis also suggests that brain mass evolved proportionally to body size in non-elephantimorph proboscideans ('Plesielephantiformes') and among more derived elephantids (Fig. 2), which is supported by the significant correlation between brain and body mass in the dataset (p-value = 0.0004). This would support the idea that brain mass evolved as a 'simple passenger' of body mass in these two groups, as hypothesized by Manger¹. Given the relationship between brain and body mass variation, it can be predicted that non-elephantimorph taxa with body masses already approaching or exceeding two tons, such as *Barytherium* and the Deinotheriidae<sup>27,34,75,76</sup> would have had an absolutely larger brain than the basal most proboscideans. This assumption might be addressed in the future, when more data on the endocasts of *Barytherium* and the Deinotheriidae becomes available.

In contrast, the ACSRML suggests that brain mass increases much faster than body mass in Elephantimorpha, resulting in two pulses of increase in the EQ: (i) the EQ doubles in the LCA of Elephantimorpha and (ii) it almost doubles again in the LCA of Elephantoidea (Fig. 2). Subsequently, the EQ essentially remains the same among Elephantoidea, including modern proboscideans (Fig. 2). According to the ACSRML (Fig. 2), brain mass began increasing in the basal-most elephantimorphs, the Mammutida (*M. americanum* and *Z. borsoni*), where the EQ is about twice that of *Moeritherium* and *Palaeomastodon* (Table 1). A brain mass equivalent to that of most modern elephants is achieved in the LCA of the Elephantoidea (Fig. 2). These two pulses suggest that, at least in the LCA of Elephantimorpha and Elephantoidea, increases in brain mass were decoupled from increases in body mass, and that an increase in EQ was positively selected for.

What could have driven such increases in EQ? According to both the fossil record and molecular studies, the LCA of Elephantimorpha lived some 26 million years ago in the late Oligocene<sup>28,29,57,77</sup>, and the LCA of the
Elephantoidea lived around 20 million years ago in the early Miocene²⁸,²⁹. These time points are contemporaneous with two periods of increased aridity on the African continent, accompanied by glaciations, arrival of invasive faunal elements from Eurasia, dispersal of proboscideans out of Africa and the spread of C4 plants²⁸,⁷⁸⁻⁸⁴ (Fig. 2). Gregariousness¹⁷,⁸⁵, greater behavioural flexibility⁸⁶ and the ability to produce complex long distance, infrasonic calls¹²,²⁵,⁸⁷,⁸⁸ have been documented in Miocene elephantiforms and may have played a role in increasing the fitness of relatively larger brained individuals in the context of ecological changes¹⁰,¹⁴,¹⁵. An enhanced long term memory may have facilitated the recalling of the location of distant water holes during increasingly frequent droughts¹¹, though evidence from the fossil record is lacking¹⁰. A shift in proboscidean diets (as a result of the advent of the conveyor-belt type of dental replacement and increasing hypsodonty)⁷⁷,⁸⁶ seems an unlikely driver for the evolution of a higher EQ because (i) proboscideans remained primarily browsers until elephantids became grazers some 7 million years ago⁸⁶,⁸⁹,⁹⁰, and (ii) no undisputable correlation is evident between diet and brain size in 'ungulates'¹⁴,¹⁵.

Alternatively, given that brain and body size variations appear significantly correlated in proboscideans (Fig. 1), it is possible that body mass, not brain mass, could have been the main trait under selective pressure during these two time intervals. Positive selection for larger body size could have provided early elephantimorphs with protection against the novel invasive predators and competitors from Eurasia, and large herbivores are known to be less sensitive to droughts and increased seasonality as they can store more fat and water, and can subsist on lower quality food (e.g. grasslands) due to their larger digestive tract and lower metabolic rate⁷⁸,⁹¹,⁹². In this respect, it is noteworthy that other 'ungulate' taxa contemporaneously show an overall increase of their brain and body mass in the early Neogene¹,⁹³,⁹⁴. This would make the increase in EQ merely a correlate of the increase of body mass, as already hypothesized by Manger *et al.*¹. However, brain mass increases faster than body mass in this portion of the tree, and the latter even decreases slightly in the LCA of Elephantoidea compared to the LCA of Elephantimorpha (Fig. 2), which does not match what would be expected under the assumption that body mass was the main driver of EQ increase.

## Conclusion

Overall, there does not seem to be one single hypothesis accounting for the complete evolution of the brain in proboscideans. Part of the evolution of a larger absolute brain mass may be explained by the co-evolution of brain and body mass (in 'plesielephantiforms' and elephantids) and part of this process seems to better explained by potential selective pressures (in the LCAs of Elephantoidea and Elephantimorpha), perhaps as a result of climatic and environmental changes in Africa, as well as with dispersal events. Despite this, brain size evolution remains significantly correlated to body size in proboscideans, which cautions against an explanatory scenario involving intelligence alone as the main focus of natural selection. The precise timing of the two identified pulsed increases in relative brain size remains poorly temporally constrained due to the small sample size that makes it difficult to obtain reasonable phylogenetic resolution. Considerable effort will have to be put into collecting data for key taxa that will fill the temporal voids in the current data set. The inclusion of Barytheriidae and Deinotheriidae - as some of the closest relatives to Elephantiformes and the earliest large proboscideans- and investigation of whether they evolved a larger brain independently of other 'Plesielephantiformes' will prove crucial in this respect¹.

## Data Availability

All the data necessary to reproduce this research are provided with the manuscript.

## References

1. Manger, P. R., Spocter, M. A. & Patzke, N. The Evolutions of Large Brain Size in Mammals: The 'Over-700-Gram Club Quartet'. *Brain Behav. Evol.* **82**, 68–78 (2013).

2. Holloway, R. L. On the Making of Endocasts: The New and the Old in Paleoneurology in Digital Endocasts (eds Bruner, E., Ogihara, N. & Tanabe, H.), 1–8. (Springer Japan KK, 2018).

3. Carlson, K. J. *et al.* The Endocast of MH1, *Australopithecus sediba*. *Science* **333**, 1402–1407 (2011).

4. Zollikofer, C. P. E. & De León, M. S. P. Pandora's growing box: Inferring the evolution and development of hominin brains from endocasts. *Evol. Anthropol.* **22**, 20–33 (2013).

5. Beaudet, A. *et al.* The endocast of StW 573 ('Little Foot') and hominin brain evolution. *J Hum Evol* **126**, 112–123 (2019).

6. Marino, L. A Comparison of Encephalization between Odontocete Cetaceans and Anthropoid Primates. *Brain Behav. Evol.* **51**, 230–238 (1998).

7. Manger, P. R. An examination of cetacean brain structure with a novel hypothesis correlating thermogenesis to the evolution of a big brain. *Biol. Rev.* **81**, 293 (2006).

8. Marino, L., McShea, D. W. & Uhen, M. D. Origin and evolution of large brains in toothed whales. *Anat. Rec.* **281A**, 1247–1255 (2004).

9. Montgomery, S. H. *et al.* The evolutionary history of Cetacean brain and body size. *Evol.* **67**, 3339–3353 (2013).

10. Benoit, J. A new method of estimating brain mass through cranial capacity in extinct proboscideans to account for the non-neural tissues surrounding their brain. *J. Vertebr. Paleontol.* **35**, e991021 (2015).

11. Hart, B. L., Hart, L. A. & Pinter-Wollman, N. Large brains and cognition: Where do elephants fit in? *Neurosci. Biobehav. Rev.* **32**, 86–98 (2007).

12. Shoshani, J., Kupsky, W. J. & Marchant, G. H. Elephant brain. *Brain Res. Bull.* **70**, 124–157 (2006).

13. Aiello, L. C. & Wheeler, P. The expensive-tissue hypothesis: The brain and the digestive system in human and primate evolution. *Curr. Anthropol.* **36**, 199–221 (1995).

14. Shultz, S. & Dunbar, R. I. M. Both social and ecological factors predict ungulate brain size. *Proc. Roy. Soc. B* **273**, 207–215 (2006).

15. Pérez-Barbería, F. J. & Gordon, I. J. Gregariousness increases brain size in ungulates. *Oecologia* **145**, 41–52 (2005).

16. DeCasien, A. R., Williams, S. A. & Higham, J. P. Primate brain size is predicted by diet but not sociality. *Nat. Ecol. Evol.* **1**(5), 112 (2017).

17. Powell, L. E., Isler, K. & Barton, R. A. Re-evaluating the link between brain size and behavioural ecology in primates. *Proc. Roy. Soc. B* **284**, 20171765 (2017).

18. Marino, L. *et al.* A claim in search of evidence: reply to Manger's thermogenesis hypothesis of cetacean brain structure. *Biol. Rev. Camb. Philos. Soc.* **83**, 417–440 (2008).

19. Fox, K. C. R., Muthukrishna, M. & Shultz, S. The social and cultural roots of whale and dolphin brains. *Nat. Ecol. Evol.* **1**, 1699–1705 (2017).
20. Xu, S. et al. Genetic basis of brain size evolution in cetaceans: insights from adaptive evolution of seven primary microcephaly (MCPh) genes. BMC Evol. Biol. **17**(1), 206 (2017).

21. Manger, P. R. Questioning the interpretations of behavioral observations of cetaceans: Is there really support for a special intellectual status for this mammalian order? *Neurosci.* **250**, 664–696 (2013).

22. Manger, P. R., Hemingway, J., Spocter, M. A. & Gallagher, A. The mass of the human brain: is it a spandrel? In African Genesis (eds Reynolds, S. C., Gallagher, A.), 181–204. (Cambridge University Press, 2012).

23. Deacon, T. W. Fallacies of progression in theories of brain-size evolution. *Int. J. Primatol.* **11**, 193–236 (1990).

24. Healy, S. D. & Rowe, C. A critique of comparative studies of brain size. *Proc. Roy. Soc. B* **274**(1609), 453–464 (2012).

25. Shoshani, J. & Tassy, P. Advances in proboscidean taxonomy & classification, anatomy & physiology, and ecology & behavior. *Quat. Int.* **126–128**, 5–20 (2005).

26. Gheerbrant, E. & Tassy, P. L'origine et l'évolution des éléphants. *C. R. Palevol* **8**, 281–294 (2009).

27. Sanders, W. J., Gheerbrant, E., Harris, J. M., Saegusa, H. & Delmer, C. Proboscidea. In *Cenozoic Mammals of Africa* (eds Werdelin, L., Sanders, W. J.), 161–251 (University of California Press, 2010).

28. Van Der Made, J. The evolution of elephants and their relatives in the context of a changing climate and geography in Elefantenreich – Eine Fossilwelt in Europa (ed. Meller, H.), 341–360. (Landesamt für Denkmalpflege und Archäologie Sachsen-Anhalt & Landesmuseum Vorgeschichte, 2010).

29. Shauer, K. Commentary and sources for the evolutionary diagram of the Proboscidea in Africa and Eurasia in Elefantenreich – Eine Fossilwelt in Europa (ed. Meller, H.), 631–650 (Landesamt für Denkmalpflege und Archäologie Sachsen-Anhalt & Landesmuseum Vorgeschichte, 2010).

30. Fisher, D. C. Paleobiology of Pleistocene Proboscideans. *Annu. Rev. Earth. Planet. Sci.* **46**, 229–260 (2018).

31. Jerison, H. J. *Evolution of the Brain and Intelligence* (Academic Press, 1973).

32. Palombo, M. R. & Giovinazzo, C. *Elephas falconeri* from Spinagallo Cave (South-Eastern Sicily, Hyblean Plateau, Siracusa): brain to body weight comparison. *Monogr. Soc. Hist. Nat. Balears* **12**, 255–264A (2005).

33. Fisher, D. C. et al. X-ray computed tomography of two mammoth calf mummies. *J. Paleontol.* **88**, 664–675 (2014).

34. Larramendi, A. & Palombo, M. R. Body size, biology and encephalization quotient of *Palaeoloxodon ex gr. P. falconeri* from Spinagallo Cave (Hyblean plateau, Sicily). *Hystrix Ital. J. Mammal.* **26**, 102–109 (2015).

35. Kharlamova, A. S. et al. The mummified brain of a Pleistocene woolly mammoth (*Mammuthus primigenius*) compared with the brain of the extant African elephant (*Loxodonta africana*). *J. Comp. Neurol.* **523**, 2326–2343 (2015).

36. Benoit, J., Crumpton, N., Mérigeaud, S. & Tabuce, R. A Memory Already like an Elephant's? The Advanced Brain Morphology of the Last Common Ancestor of Afrotheria (Mammalia). *Brain Behav. Evol.* **81**, 154–169 (2013).

37. Larramendi, A. Shoulder height, body mass, and shape of proboscideans. *Acta Palaeontol. Pol.* **61**, 537–574 (2016).

38. Obada, T. Preliminary data on the *Mammut borsoni* (Hays, 1834) from Otman Hill (Colibaş, Republic of Moldova). *Abstract Book of the VIth International Conference on Mammoths and their Relatives. S.A.S.G.* **102**, 143–144 (2014).

39. Gervais, P. Mémoire sur les formes cérébrales propres à différents groupes de mammifères. *J. Zool.* **1**, 425–469 (1872).

40. Andrews, C. W. A Descriptive Catalogue of the Tertiary Vertebrata of Fayum, Egypt. (British Museum Natural History press, 1906).

41. Tassy, P. The classification of Proboscidea: How many cladistic classifications? *Cladistics* **4**, 43–57 (1988).

42. Radinsky, L. Early primate brains: Facts and fiction. *J. Hum. Evol.* **6**, 79–86 (1977).

43. Fahlke, J. M. & Autenrieth, M. Photogrammetry vs. Micro-CT scanning for 3D surface generation of a typical vertebrate fossil – A case study. *J. Paleontol. Tech.* **14**, 1–18 (2016).

44. Sutton, M., Rahman, I. & Garwood, R. Virtual paleontology – an overview. *Paleontol. Soc. Pap.* **22**, 1–20 (2016).

45. Evin, A. et al. The use of close-range photogrammetry in zooarchaeology: Creating accurate 3D models of wolf crania to study dog domestication. *J. Archaeol. Sci. Rep.* **9**, 87–93 (2016).

46. Eisenberg, J. F. *The mammalian radiations* (University of Chicago Press, 1981).

47. Stephan, H., Frahm, H. & Baron, G. New and revised data on volumes of brain structures in insectivores and primates. *Folia Primatol.* **35**, 1–29 (1981).

48. Revell, L. J. phytools: an R package for phylogenetic comparative biology (and other things). *Meth. Ecol. Evol.* **3**, 217–223 (2011).

49. Palkopoulou, E. et al. A comprehensive genomic history of extinct and living elephants. *Proc. Nat. Acad. Sci.* **115**, E2566–E2574 (2018).

50. Meyer, M. et al. Palaeogenomes of Eurasian straight-tusked elephants challenge the current view of elephant evolution. *eLife* **6**, e25413 (2017).

51. Poulakakis, N., Theodorou, G. E., Zouros, E. & Mylonas, M. Molecular Phylogeny of the Extinct Pleistocene Dwarf Elephant *Palaeoloxodon antiquus* falconeri from Tilos Island, Dodekanisa, Greece. *J. Mol. Evol.* **55**, 364–374 (2002).

52. Poulakakis, N. et al. Ancient DNA forces reconsideration of evolutionary history of Mediterranean pygmy elephantids. *Biol. Lett.* **2**, 451–454 (2006).

53. Osborn, H. F. *Palaeoloxodon antiquus italicus* sp. nov., final stage in the 'Elephas antiquus' phylum. *Am. Mus. Novit.* **460**, 1–24 (1931).

54. Osborn, H. F. Proboscidea: a monograph of the discovery, evolution, migration and extinction of the mastodonts and elephants of the world. Vol. II: Stegodontoidea, Elephantoidea (The American Museum Press, 1942).

55. Palombo, M. R. Paedomorphic features and allometric growth in the skull of Elephas falconeri from Spinagallo (Middle Pleistocene, Sicily) in Proceedings of the 1st International Congress of "La Terra degli Elefanti", The World of Elephants (eds Cavarretta, G., Gioia, P., Mussi, M., Palombo, M. R.), 492–496 (Consiglio Nazionale delle Ricerche, 2001).

56. Palombo, M. R. & Giovinazzo, C. *Elephas falconeri* from Spinagallo Cave (south-eastern Sicily, Hyblean Plateau, Siracusa): A preliminary report on brain to body weight comparison. *Monogr. Soc. Hist. Nat. Balears* **12**, 255–264 (2005).

57. Rohland, N. et al. Proboscidean Mitogenomics: Chronology and Mode of Elephant Evolution Using *Mastodon* as Outgroup. *PLoS Biol.* **5**, e207 (2007).

58. Felsenstein, J. Phylogenies and the Comparative Method. *Am. Nat.* **125**, 1–15 (1985).

59. Grafen, A. The phylogenetic regression. *Philos. Trans. Royal Soc. B* **326**, 119–157 (1989).

60. Garamszegi, L. Z. & Mundry, R. Multimodel-inference in comparative analyses in Modern Phylogenetic Comparative Methods and Their Application in Evolutionary Biology (ed. Garamszegi, L. Z.), 305–331 (Springer Berlin Heidelberg, 2014).

61. Blomberg, S. P. & Garland, T. Jr. Tempo and mode in evolution: phylogenetic inertia, adaptation and comparative methods. *J. Evol. Biol.* **15**, 899–910 (2002).

62. Mascaro, J., Litton, C. M., Hughes, R. F., Uowolo, A. & Schnitzer, S. A. Is logarithmic transformation necessary in allometry? Ten, one-hundred, one-thousand-times yes. *Biol. J. Linn. Soc.* **111**, 230–233 (2014).

63. Montgomery, S. H., Mundy, N. I. & Barton, R. A. Brain evolution and development: adaptation, allometry and constraint. *Philos. Trans. Royal Soc. B* **283**, 20160433 (2016).

64. Burnham, K. P., Anderson, D. R. & Huyvaert, K. P. AIC model selection and multimodel inference in behavioral ecology: some background, observations, and comparisons. *Behav. Ecol. Sociobiol.* **65**, 23–35 (2010).

65. Mazerolle, M. J. AICcmodav: Model selection and multimodel inference based on (Q)AIC(c). R package version 2.1-1 (2017).

66. Felsenstein, J. Maximum likelihood estimation of evolutionary trees from continuous characters. *Am. J. Hum. Genet.* **25**, 471–492 (1973).

67. Butler, M. A. & King, A. A. Phylogenetic Comparative Analysis: A Modeling Approach for Adaptive Evolution. *Am. Nat.* **164**, 683–695 (2004).

68. Pagel, M. Inferring the historical patterns of biological evolution. *Nature* **401**, 877–884 (1999).
69. Harmon, L. J. et al. Early bursts of body size and shape evolution are rare in comparative data. *Evol.* **64**(8), 2385–2396 (2010).

70. Pinheiro, J., Bates, D., DebRoy, S., Sarkar, D. & R. Core Team. nlme: Linear and Nonlinear Mixed Effects Models. R package version 3.1-137 (2018).

71. Paradis, E., Claude, J. & Strimmer, K. APE: Analyses of Phylogenetics and Evolution in R language. *Bioinform.* **20**, 289–290 (2004).

72. Gross, J. & Ligges, U. nortest: Tests for Normality. R package version 1.0-3 (2015).

73. Mundry, R. Statistical Issues and Assumptions of Phylogenetic Generalized Least Squares in Modern Phylogenetic Comparative Methods and Their Application in Evolutionary Biology (eds Garamszegi, L.), 131–153 (Springer Berlin Heidelberg, 2014).

74. Gheerbrant, E. et al. Nouvelles données sur *Phosphatherium escuilliei* (Mammalia, Proboscidea) de l'Eocène inférieur du Maroc, apports à la phylogénie des Proboscidea et des ongulés lophodontes. *Geodiversitas* **27**, 239–333 (2005).

75. Christiansen, P. Body size in proboscideans, with notes on elephant metabolism. *Zool. J. Linn. Soc.* **140**, 523–549 (2004).

76. Evans, A. R. et al. The maximum rate of mammal evolution. *Proc. Nat. Acad. Sci.* **109**, 4187–4190 (2012).

77. Shoshani, J. et al. A proboscidean from the late Oligocene of Eritrea, a 'missing link' between early Elephantiformes and Elephantimorpha, and biogeographic implications. *Proc. Nat. Acad. Sci.* **103**, 17296–17301 (2006).

78. Kappelman, J. et al. Oligocene mammals from Ethiopia and faunal exchange between Afro-Arabia and Eurasia. *Nature* **426**, 549–552 (2003).

79. Antoine, P.-O. et al. First record of Paleogene Elephantoidea (Mammalia, Proboscidea) from the Bugti Hills of Pakistan. *J. Vertebr. Paleontol.* **23**, 977–980 (2003).

80. Feakins, S. J. & Demenocal, P. B. Global and African Regional Climate during the Cenozoic in Cenozoic Mammals of Africa (eds Werdelin, L., Sanders, J. S.), 45–56. (University of California Press, 2010).

81. Jacobs, B. F., Pan, A. D. & Scotese, C. R. A Review of the Cenozoic Vegetation History of Africa in Cenozoic Mammals of Africa (eds Werdelin, L. & Sanders, J. S.), 57–72. (University of California Press, 2010).

82. Sen, S. Dispersal of African mammals in Eurasia during the Cenozoic: Ways and whys. *Geobios* **46**, 159–172 (2013).

83. Tabuce, R. & Marivaux, L. Mammalian interchanges between Africa and Eurasia: an analysis of temporal constraints on plausible anthropoid dispersals during the Paleogene. *Anthropol. Sci.* **113**, 27–32 (2005).

84. Zachos, J. Trends, Rhythms, and Aberrations in Global Climate 65 Ma to Present. *Science* **292**, 686–693 (2001).

85. Bibi, F. et al. Early evidence for complex social structure in Proboscidea from a late Miocene trackway site in the United Arab Emirates. *Biol. Let.* **8**, 670–673 (2012).

86. Lister, A. M. The role of behaviour in adaptive morphological evolution of African proboscideans. *Nature* **500**, 331–334 (2013).

87. O'Connell-Rodwell, C. E. Keeping an 'ear' to the ground: seismic communication in elephants. *Physiology* **22**, 287–294 (2007).

88. Benoit, J., Merigeaud, S. & Tabuce, R. Homoplasy in the ear region of Tethytheria and the systematic position of Embrithopoda (Mammalia, Afrotheria). *Geobios* **46**, 357–370 (2013).

89. Cerling, T. E., Harris, J. M. & Leaky, M. G. Browsing and Grazing in Elephants: The Isotope Record of Modern and Fossil Proboscideans. *Oecologia* **120**, 364–374 (1999).

90. Cerling, T. E., Harris, J. M. & Leaky, M. G. Environmentally Driven Dietary Adaptations in African Mammals in A History of Atmospheric CO₂ and its Effects on Plants, Animals, and Ecosystems. Ecological Studies (eds Ehleringer, J. R. et al.), 258–272 (Springer, New York, 2005).

91. Peters, R. H. *The Ecological Implications of Body Size* (Cambridge University Press, 1983).

92. Sander, P. M. et al. Biology of the saupod dinosaurs: the evolution of gigantism. *Biol. Rev.* **86**, 117–155 (2010).

93. Radinsky, L. Evolution of Brain Size in Carnivores and Ungulates. *Am. Nat.* **112**(987), 815–831 (1978).

94. Huang, S. et al. Mammal body size evolution in North America and Europe over 20 Myr: similar trends generated by different processes. *Proc. Roy. Soc. B* **284**, 20162361 (2017).

## Acknowledgements

The authors would like to thank S. Vrard, C. Bibal, M.R. Palombo, P. Tassy, P. Brewer, C. Argot, M. Cazenave, C. Hendrickx; and to the financial support of the DST-NRF Centre of Excellence in Palaeosciences, the Claude Leon Foundation and the des Treilles Foundation (Paris): “La Fondation des Treilles, créée par Anne Gruner Schlumberger, a notamment pour vocation d'ouvrir et de nourrir le dialogue entre les sciences et les arts afin de faire progresser la création et la recherche contemporaines. Elle accueille également des chercheurs et des écrivains dans le domaine des Treilles (Var, France). www.lestreilles.com”.

## Author Contributions

All authors had full access to all the data in the study and gave final approval for publication. Study concept and design: J.B. and P.M. Acquisition of data: J.B., R.T., T.O. and V.M. Analysis and interpretation of data: J.B. and L.J.L. Drafting of the manuscript: J.B., L.J.L., R.T. and P.M. Critical revision of the manuscript for important intellectual content: T.O., R.T. and P.M. Obtained funding: J.B.

## Additional Information

Supplementary information accompanies this paper at https://doi.org/10.1038/s41598-019-45888-4.

**Competing Interests:** The authors declare no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2019