De Lussanet claims that our model that accounts for the degree of folding of the cerebral cortex based on the product of cortical surface area and the square root of cortical thickness is better reduced to the product of gray-matter proportion and folding index. Lewitus et al., in turn, claim that the assumptions of our model are in conflict with experimental data; that the model does not accurately fit the data; and that the ancestral mammalian brain was gyrencephalic. Here, we show that both claims are inappropriate.
In our original Report (1), we showed for the first time that the relationship between total surface area (A T ) and exposed surface area (A E ) of the mammalian cerebral cortex can be universally described, across lissencephalic and gyrencephalic species alike, by the power, where T is the average thickness of the cortical gray matter (Fig. 1A). Our model predicts the degree of folding of the mammalian cerebral cortex (that is, the relationship between A T and A E ) with an r 2 of 0.998; that is, 99.8% of the variance in A T T 1/2 across species is explained by the variation in kA E 1.305 . The empirical function is very close to our theoretical function A T T 1/2 = kA E 1.25 that expresses the relationship between A T and A E , given T, that minimizes the effective free energy of a growing cortical volume-that is, the function that describes the most stable conformation of a growing, elastic, deformable cerebral cortex given its surface area and thickness. We thus propose that folding occurs as the expanding cortex, during development, settles at each moment in time into the most energetically favorable conformation, one that depends simply on its current combination of A T and T, not number of neurons.Based on the calculation from our data of values of k for the theoretical exponent of 1.25, de Lussanet notes that the residuals of the function, that is, k, vary systematically with increasing A E across species (2). This indicates a (presumably hidden) significant discrepancy between the exponent for A E predicted by our effective free energy-minimization theory and the empirical exponent 1.305 (excluding cetaceans) that is found in the data. Our original Report makes it clear that, with a value of 1.305 ± 0.010, the exponent for A E that best fits the data is close to, but statistically distinguishable (given the standard error) from, the theoretical exponent of 1.250, which is not surprising for a mean-field model based on a number of simplifying assumptions.The origin of the discrepancy of 0.055 in the scaling exponent is an important question. In purely theoretical terms, the use of A E as an order parameter was a necessary simplification, based on the expected relationship V T~AE 3/2 [made explicitly in the supplementary information in (1), and a good empirical approximation, given that we obtain V T~AE 1.52 ]. Whether modifications to correct for Gaussian curvature of the surface and the presence of deep nuclei would suffice to close the gap of 0.055 between theoretical and empirical e...