Introduction
Neurons, the primary components of central nervous system are electrically excitable cells that receive, process, and transmit information through electrical and chemical signals between each other. Digital reconstructions of neurons provide information for quantitative measurements of neuronal morphology. With digital reconstructions, many studies have highlighted the importance of neuronal morphology in its biological function. Mainen and Sejnowski [1] have shown a causal relationship between dendritic structures and intrinsic firing patterns observed from in vitro electrical recordings for a wide variety of cell types. Vetter et al. have shown that branching pattern strongly affected the propagation of action potentials which links information processing at different regions of the dendritic tree [2]. Ferrante et al. have shown that even subtle membrane readjustments at branch point could drastically alter the ability of synaptic input to generate and propagate the action potentials [3]. It has also been shown that the inputoutput response function of neuron’s dendritic arbor grows with the tree size [4, 5]. More recently, Yi et al. have demonstrated a crucial role of neuronal morphology in determining fieldinduced neural response [6]. These results, among others, have contributed to a now widespread acceptance that neuronal morphology plays a critical role in its activity and function.
The formation of neurons through branching is driven by complex interactions of intracellular and extracellular signaling cascades which are proving difficult to be completely understood by molecular biology alone. Mathematical or computational modelling approaches instead provides an alternative and complementary approach to uncover mechanism underlying neuronal morphology. Simulators LNeuron [7, 8] and NeuGen [9] were developed to create virtual neuronal structures in silico by means of iteratively sampling experimental statistical distributions of shape parameters (including e.g. branch diameter, length, ect). In contrast to LNeuron or NeuGen simulators which require a large number of experimental neuron samples to obtain reliable distributions of shape parameters, Van Ooyen et al. used models in which the morphology of a single dendritic tree was represented in a highly abstract manner [10, 11, 12, 13] where growing neurons were modeled as expanding, circular neuritic fields. Based on such phenomenological dendritic growth models, NETMORPH was developed to simulate 3D neuronal networks from the perspective of individual growth cone [14]
, using simple rules for neurite branching probability at each terminal segment. Other modelling approach with simple rules to account for the spatial embedding of tree structures, includes optimal rewiring and particlebased diffusion limited aggregation approach: the optimal rewiring approach generates branching geometry by minimizing the wiring length and the pathlength to root between branch points and synapses in dendrite trees
[15, 16]; the particlebased diffusion limited aggregation approach provides modelbased measures to estimate “diffusive” shape of neuronal treelike structures [17, 18].However, analysis and modelling of neuron morphology to the best of our knowledge so far are largely based on neuron digital reconstructions from certain regions or layers of a brain (e.g. [7, 14]). Moreover, a majority of these neuronal morphology analysis and modelling concentrate on neuron dendrites, with only a few studies of axonal branching structure in certain layers recently (e.g. [19, 20]). Knowledge about how dendrites and axons branch in an entire brain is still limited. Recently, breakthroughs in imaging [21, 22] and molecular labeling [23, 24] techniques have provided tools to trace and digital reconstruct the almost complete morphology of neuronal populations at a singleaxon resolution through a whole brain [25, 26]; this offers opportunity to quantify brainwide dendritic and axonal branching morphology. In this manuscript, we take advantage of current availability of brainwide neuron digital reconstructions of pyramidal cells from a mouse brain [25] and analyze the complete axonal and dendritic branching morphology. In particular we show that axons are topologically selfaffine whereas dendrites are topologically selfsimilar, and topological structures of both axon and dendrites are far away from symmetry and appear to be random. We also develop a selforganized probabilistic model for the entire axonal and dendritic branching structures. In contrast to growth models [13] which require a simulation/growth time to create virtual trees of finite size, we use the generation of rooted 3Cayley trees (where nonterminal nodes are linked to 3 neighbors; in neuron trees, 3 linkages of a node represent one mother branch and two daughter branches) via a stochastic branching process. Using estimated branching probability from brainwide axonal and dendritic trees, we show that this simple probabilistic model is sufficient to quantitative recapture several statistical properties of neuronal morphology including distributions of neuron size and topological width/length.
Results
Brainwide axonal and dendritic morphology quantifications
In this section, we analyse neuron morphology from 35 brainwide digital reconstructions of pyramidal neurons in a mouse brain [25]. A typical neuron consists of a cell body (soma), dendrites, and an axon. Here we analyze axon and dendrite tree morphology in aspects of size, asymmetry, and shape of trees as well as correlation among axon and dendrites in single neurons.
Topological structure of axonal and dendritic arbors
We first focus on topological structure of dendritic and axonal branching from soma to neuron terminals, despite the spatial position of each branch. Fig 1 illustrates an example of a single brainwide neuronal morphology and its topological branching structure. We seperate individual neurons into axon and dendrites for analysis; see Fig 1 as an example. From 35 actual neuron reconstructions, we find that on average there are number of dendrites besides one axon. Neurons have been widely viewed as binary trees (e.g.[27]), here we first examine this by looking into the degree (i.e. the number of branches each node links to) of each node. Clearly, each terminal is connected to one branching node (i.e. degree 1) and branching nodes are mainly linked with three branches in both axons and dendrites  with only of branching nodes in axons and in dendrites linking to more than 3 branches; those could be due to limited spatial resolution in imaging which leads to possible errors in digital reconstructions. Overall, this suggests that brainwide neurons are well approximated as binary trees.
We quantify the number of branching nodes (which is equivalent to the number of terminals minus one) for both axonal and dendritic trees, and refer as the topological size. On average, we find that axons are of size , whereas dendrites are significantly smaller, of size , i.e. about
of the axonal tree size. We also classify the branching nodes into three types
[28]: Btype if both child branches themselves bifurcate; Mtype if only one child branch bifurcate; Stype if branching with two terminals. By such classification, we find that axons have smaller proportion of Stype nodes whereas dendrite trees have larger proportion of Stype branching nodes; see Fig 2 (a).We characterize the topological shape of neuron trees by measuring its topological width and length as illustrated in Fig 1 (c). To do this, we first give an order (denoted as ) for each node as its topological distance to the soma, i.e. the number of branches in its path to the soma. The topological length (denoted as ) of a tree is then defined as the maximal order among nodes. For the topological width, we calculate the number of branching nodes (denoted as ) for each order and define the topological width (denoted as ) of a tree as the maximum among all [29]. This is similar to Sholl analysis [30] when topological distance is applied. In the population of our 35 actual brainwide neuron reconstructions, we observe that on average and the scaling exponents are different between axon () and dendrite (); see Fig 2·(b). Note that in axons, exponents and are significantly different (F test; ) where in dendrites there is no significant difference between the exponents (Ftest ); this suggests that in aspects of topological shape, dendrites are approximately selfsimilar (i.e. ) whereas axons are selfaffine (i.e. ).
Next, we quantify morphological asymmetry of neuron trees. Following the definition in [31], tree asymmetry describes the average of local partition asymmetry for all branching nodes (here are the number of terminals in two subtrees from a branching node) except for those with partition
(i.e. branching to two terminals). This tree asymmetry reflects the relative balance of branching within a tree based on the distribution of terminals between two subtrees birfucated from a node. This asymmetry index ranges from 0 (symmetry) to 1 (asymmetry). The axon and dendrite from brainwide neuron digital reconstructions show no significant difference (student ttest
) on tree asymmetry with for axonal trees and for dendritic trees and both are far away to symmetry. These measured partition asymmetry values are close to reported values in http://neuromorpho.org/. Another asymmetry measurement  excess partition asymmetry introduced by Samsonovich and Ascoli [32] considers the difference between partition asymmetry actually measured at one branch and the average of partition asymmetry computed for the same branch after all possible shuffling of the granddaughter branches. This excess partition asymmetry is useful to test the randomness in branching or the existence of control/regulation process within neurons. Interestingly, in contrast to ’regulation mechanism’ suggested in [32] based on systematically positive excess partition asymmetry measured from pyramidal cells on some slice of a brain [33], we show that our brainwide axonal and dendritic trees appear to have almost no excess partition asymmetry with (ttest to zero mean with ) for axon and (ttest ). This suggests branching in both axonal and dendritic trees appears to be random.Geometric size of axon and dendrites
Regarding geometric size of neuron trees, we measure the total branch length of a tree (i.e. the sum of all branch length in a tree). On average, axons have a total branch length of (median ) and dendrites have a total branch length (median ). Moreover, we observe that the total length strongly correlates with tree size (Pearson for axons and for dendrites) and increase approximately linear as the tree size for both axons and dendrites; see Fig 3 top panels. The mean total length is expected to be the mean branch length (say ) multiplied by the number of branches for a given tree size . The data show a best linear fit to the function with a mean branch length for axons and for dendrites. Such strong correlation and similar linear slope between axon length and its tree size were also observed in vitro axons of human neurons across all layers of medial temporal cortex, though total length of axons in human brains is much larger than that in mouse [20]; this suggests a potential universal mechanism controlling branch length cross species. We also note that in our brainwide neurons, though the total branch length approximately grows linearly with the number of branches, the branch length is not uniform along orders; see Fig 3 bottom panels for both axons and dendrites. In particular, the branch length of axons is larger at intermediate orders than those at smaller or larger branch orders.
Correlation among axon and dendrites in single neurons
It has been suggested that in one given portion of a neuron, dendritic morphology may be under intrinsic homeostatic control which regulates tree size fluctuations systematically by counterbalancing the remaining dendrites in the same cell [32]. We test such a control mechanism in axons and dendrites from our actual brainwide neurons. If trees of a neuron were mainly up(down)regulated by common factors, then neurons with larger axons in terms of tree size would have larger dendrites; if axons and dendrites were mutually regulated by competition, then neurons with larger axons would have smaller dendrites. Based on calculated tree size of axons, average tree size among dendrites and the number of dendrites for each of 36 neurons, we perform Pearson correlation analysis among these quantities. Interestingly, we find no significant correlation between axon tree size and average dendrite tree size (), between axon tree size and number of dendrites (), or between average dendrite tree size and the number of trees (
). These suggest that axon tree size, dendrite tree size and dendrite number in single neurons are likely to be independent. Moreover, if tree size is regulated within single neurons, then tree shuffling would lead to an increase in the variance of average tree size in the neuron population; if tree size is random among populations, then shuffling would not lead to significant change on the variance
[32]. To further test such randomness, we consider neurons with the same number of dendrites. Then we randomly shuffle the dendrites in these neurons and calculate the standard deviation (std) of the average dendrite tree size among single neurons. Fig
4 (c) show that 35.5% of 1000 random shuffling gives higher std of the average dendrite tree size than the std from actual neurons, suggesting tree size of dendrites within single neurons is random.Modelling topology of axonal and dendritic trees
As indicated by above statistical analysis, the tree size of dendrites and axons appear to be random independent of the number of dendrites, indicating that the formation of brainwide neuron trees is likely to be a stochastic process. In this section, we use stochastic generation of rooted zCayley trees [34] to model the topological structure of axonal and dendritic arbors. A zCayley tree is a tree where nonterminals have zneighbours. Here to model neuron trees, we use as nonterminal nodes in neuron trees are shown to mainly have 3 branches (3 neighboring nodes): one in previous order and two in the next order. In the generation of a rooted 3Cayley tree, we start with an arbitrary node as the root (which represents the soma in a neuron tree), and assign an link from this root connecting to another node (viewed as the 1th order node); the generation of tree proceeds from this 1th order node, by 1) bifurcating to two daughter branches (linking to 2 neighbouring nodes in the next order) and 2) randomly selecting its neighboring nodes in the next order with certain probability; if one node in the th order is selected, then it bifurcate into two daughter branches to next order and random selection on its two neighbours applies independently; if a neighbouring node is not selected, then this neighboring node is considered as a terminal (which could not bifurcate into the next order). The procedure continues until no more nodes are selected to bifurcate and eventually generates a rooted 3Cayley tree.
We first investigate homogeneous branching model where all nodes are selected with the equal probability . In this case it is well known that beyond a critical probability , a positive probability exists that an infinite tree is generated, whereas below this threshold, the probability for such infinite spanning is zero and a finite tree is generated [35]. Moreover, the mean tree size of generated finite trees increases with branching probability as
(1) 
In fact if we denote by the mean number of branches in one branch from the 1th order node, then by statistical equivalence between subbranches , which gives , and the mean tree size where corresponds to the 1th order branching node.
To model the axon and dendrite trees, we estimate the constant branching probability from Eq. (1) using measured average tree size, which gives for axon and for dendrite. Indeed, the model with estimated branching probability produce similar mean tree size as the actual data; see Fig 5 top two panels. We then test this simple homogeneous model for other quantities described in the above section, such as proportions of 3 types of branching nodes (we take the Btype branching node as an example for the comparison), topological shape indexes ( and ), and tree asymmetry . Note that the model takes independent probability for each node and thus no excess asymmetry is expected. Fig 5 top panels show the relative difference of the modelling results to actual data on these quantities; in particular, we see that proportion of Btype branching node and topological width () from homogeneous model from both axons and dendrites are far different ( relative difference) to measurements from actual data. Moreover, distributions on tree size (), tree topological length () and width () as shown in Fig 5 bottom panels for both axon and dendrite are significantly different to corresponding distributions from actual axonal and dendritic trees ( for both axons and dendrites). We also notice that distributions from the homogeneous model illustrate higher frequency at lower values and larger variation on these quantities for both axon and dendrite, indicating that it has high probability to generate smaller tree than the expected mean. We thus conclude that homogeneous model is not sufficient to model the topology of neuron trees.
We next consider inhomogeneous model where branching probability differs between nodes. In particular we investigate order dependent branching probability ; such dependence was considered in the growth model[12] and the computational model [36]. Fig 6 indeed shows that the branching probability does depend on the order in particular at low orders. In the growth model [12] for dendritic topology, branching probability is assumed to be decay exponentially to as the order . However, the observed branching frequencies for both axon and dendrites in our brainwide neurons do not fit well to exponential decay with zero plateau, instead, the data exhibit steady branching frequencies at high orders and statistically the data fit better to exponential decay with nonzero plateau in the form of (Ftest; for both axons and dendrites). In the best fitting curve, as , the branching probability (the critical probability); this ensures a finite tree generation in our model [35]. The plateau also reflects the almost constant branching probability at large orders where subtrees branch in an approximately “homogeneous” manner.
With the best fitting branching probability in the form of exponential decay to a nonzero plateau, the inhomogeneous model leads to mean tree sizes and for axonal and dendritic trees respectively, which are close to measured values from actual data. We also test this inhomogeneous model on other topological quantities and show in Fig 5 that the model not only recaptures the mean tree topological measurements (with relative difference to measurement from actual neurons), but also agrees well on their distributions in particular the tree size and topological shapes and (twosample Kolmogorov–Smirnov test with , except for test on of dendrite where ). We next examine the scaling of topological shapes in relation to tree size. Fig 7 shows that the conditional averages of tree length and width scale with the tree size as and with exponents (significant difference between and ; Ftest; ) for axon and (no significant difference between and ; Ftest; ) for dendrite. These exponents from simulations are close to those estimated from actual neurons, showing that as similar as in actual neurons, simulated virutal axons are selfaffine while simulated virtual dendritic trees are selfsimilar. These suggest that inhomogeneous model with a simple order dependent rule on branching probability is sufficient to statistically capture a number of topological features of both axonal and dendritic morphology in a quantitative manner.
Statistical properties of the probability model
In this section, we explore how topological characteristics of artificial trees generated from the inhomogeneous model with branching probability vary with parameters. Note that we take a operation in here, considering a probability cannot be larger than . Also note that, if we denote as the branching probability for order , then , as the first branching is initially selected in the model.
Simply, the mean tree size which can be calculated as the total number of branching nodes from all orders, and read as
(2) 
For homogeneous model , it reduces to Eq (1). For inhomogeneous model, this estimation agrees well with simulated data; see Fig 8 top panels. As expected, tree size , tree length and width increase with decreasing decay rate or increasing . Moreover, we note that for a similar variation of parameters, the decay rate allows to give a board range of values on topological quantities in particular the tree size . Furthermore, from the simulated data, we observe that on average tree length is larger than tree width when tree size is small, whereas when tree size is big (e.g. see simulations with small on the top left panel in Fig 8) tree length is smaller than tree width on average. This is further confirmed in Fig 8 bottom panel, which illustrates changes of conditional average of tree length and width for various tree size . Note that the model has three parameters and and different combinations of the parameters could give a same tree size depending on the branching probability express on . However, from Fig 8 bottom panel, we see that for the same tree size , the corresponding mean tree length and width from different combinations of parameters are similar. This suggests that the tree shape largely depends on the tree size regardless of the branching probability pattern along the orders in the model. Moreover, Fig 8 bottom panel shows that for small magnitude of tree size , the mean tree length and width grow with the tree size in a similar scale, i.e. with ; the dendrite trees are particular of this case. Whereas for a large magnitude of tree size , the mean tree width grow faster than tree length when increasing tree size , i.e. ; axons are particular of this case. An intersection between tree length and tree width in relation with tree size occurs at around ; for tree size much larger than this, the mean tree width is larger than mean tree length (i.e. ). These could be used as predictions for topological tree length and width when knowing tree size.
Conclusion and Discussion
In this manuscript we examine several neuronal morphological features and develop an inhomogeneous model to generate finite realistic virtual axonal and dendritic tree structures. By analyzing topological shape (length and width) in relation with tree size, we find that axons exhibit an selfaffine patten while dendrites are selfsimilar. Moreover, analyses on excess asymmetry and dendrite shuffling, suggest that tree size appear to be random; this supports the idea of a probability model for neutron tree structures. Using inhomogeneous branching probability, our probability model generates finite virtual trees of statistically similar to brainwide Pyramidal neutrons in a number of features (including axonal and dendritic tree size and shape). In contrast to published simulation tools [9, 7, 8, 13], this modelling approach provides a tool with simple rules in a selforganized manner to generate virtual topological neuron structures of finite size. We remark here that models using different order dependent branching probability functions (e.g. polynomial decay, modelling data not shown) which fit the measured branching frequency as shown in Fig 6, are also able to recapture measured data. Mechanisms underlying the decay pattern of branching probability remains unclear; however order dependent branching frequency alone is thus sufficient to capture the tree topology.
Our modelling of brainwide neuron branching topology can be extended to take geometric information of branches into account, e.g. branch length and branch angles as studied in [14, 36]. Note that a branch in a real neuron is not a straight line between two branching points, but exhibit complicated curvature structures; see Fig 1 (a). Including curvature structures would thus be important for a complete and better understanding of mechanism underlying brainwide neuron morphology and its relation to specific functions at a single neuron scale. One possibility for such extending could be the inclusion of a stochastic process for the outgrowth direction of each branch as considered in the NETMORPH [14]. Single neurons can expand towards different regions in the brain and different areas in the brain are specialized for different functions; e.g. left hemisphere of the brain is dedicated to language while the right hemisphere is involved in more creative activities such as drawing. Thus morphological difference between different regions in a single neuron scale is likely to link at with different regionspecified functions [37, 38, 39]. Indeed, neuron morphology and function in different layers of brains have been studied [40, 41]. However, we are still far away to understand neuron morphologyfunction relationship at a single neuron scale.
As a theoretical approach in studying morphological role of neurons in their function, our brainwide modelling of neuron branching could be useful to study the enhancement of dynamic range [42], by creating virtual realistic neuron topologies. Traditional Cayley trees (where terminals are of the same order) considered as an excitable media for input signal propagation has been recently used to investigate dynamic range [42, 4]. With our modelling approach for realistic neuron tree topology, it would be interesting to explore statistical properties of neuron trees in relation with dynamic ranges of neurons. In particular, it has been suggested that larger tree size could give larger dynamic ranges [42]; this can be tested with realistic virtual neuron trees using our modelling approach instead of traditional Cayley trees. Moreover, using more realistic virtual neuron trees also allow investigating the impact of tree asymmetry as well as tree shape in dynamic range.
Methods
In this manuscript, we use 35 neuron digital reconstructions of pyramidal cells for analysis. Pyramidal neurons digital reconstructions are provided by authors in [25]. Briefly, neuron digita reconstructions were obtained as follows. The brain is from C57BL/6J mouse line, and its pyramid neurons are sparsely labelled with Adenoassociated virus (AAV). Fluorescence microoptical sectioning tomography microscopy (fMOST) was used to image these labelled neurons which span different brain regions or even the whole brain. From the imaging dataset, GTree method was used to reconstruct these labelled neurons, all of which contained axonal and dendrite tree morphology. The animal experiments were approved by the Institutional Animal Ethics Committee of Huazhong University of Science and Technology, and all experiments were performed in accordance with relevant guidelines and regulations.
The digital reconstruction of each brainwide neuron is stored in a morphological file (in “SWC” format) which includes position information of each traced point and its linkage points. For analysis, each empirically neuron is decomposed into one axonal tree and several dendritic trees; each tree stems from the same soma. The axon part is identified as the largest component (in terms of the number of branching points) among all subcomponents and the rest are considered as dendrites.
A subcomponent of digital reconstruction stemmed from the soma is excluded for analysis if there is no branching points in this subcomponent; this in fact corresponds to a tree of size . Note that there is a small proportion of branching points which bifurcate to more than 2 branches; this could possibly due to errors in digital reconstruction. We manually adjust this by splitting such branching points so that all branching points bifurcate to 2 branches.
References
 [1] Mainen, Z. F. & Sejnowski, T. J. Influence of dendritic structure on firing pattern in model neocortical neurons. Nature 382, 363–366 (1996).
 [2] Vetter, P., Roth, A. & Häusser, M. Propagation of action potentials in dendrites depends on dendritic morphology. Journal of neurophysiology 85, 926–937 (2001).
 [3] Ferrante, M., Migliore, M. & Ascoli, G. A. Functional Impact of Dendritic BranchPoint Morphology. Journal of Neuroscience 33, 2156–2165 (2013).
 [4] Gollo, L. L., Kinouchi, O. & Copelli, M. Singleneuron criticality optimizes analog dendritic computation. Scientific Reports 3, 3222 (2013).
 [5] Publio, R., Ceballos, C. C. & Roque, A. C. Dynamic Range of Vertebrate Retina Ganglion Cells: Importance of Active Dendrites and Coupling by Electrical Synapses. PLoS ONE 7, e48517 (2012).
 [6] Yi, G. S., Wang, J., Deng, B. & Wei, X. L. Morphology controls how hippocampal CA1 pyramidal neuron responds to uniform electric fields: A biophysical modeling study. Scientific Reports 7, 3210 (2017).
 [7] Ascoli, G. A. & Krichmar, J. L. Lneuron: A modeling tool for the efficient generation and parsimonious description of dendritic morphology. Neurocomputing 3233, 1003–1011 (2000).
 [8] Ascoli, G. A., Krichmar, J. L., Nasuto, S. J. & Senft, S. L. Generation, description and storage of dendritic morphology data. Philosophical Transactions of the Royal Society B: Biological Sciences 356, 1131–1145 (2001).

[9]
Eberhard, J. P., Wanner, A. &
Wittum, G.
NeuGen: A tool for the generation of realistic morphology of cortical neurons and neural networks in 3D.
Neurocomputing 70, 327–342 (2006).  [10] Van Pelt, J. & Verwer, R. W. Growth models (including terminal and segmental branching) for topological binary trees. Bulletin of Mathematical Biology 47, 323–336 (1985).
 [11] Van Pelt, J. & Verwer, R. W. Topological properties of binary trees grown with orderdependent branching probabilities. Bulletin of Mathematical Biology 48, 197–211 (1986).
 [12] Van Pelt, J. & Schierwagen, A. Morphological analysis and modeling of neuronal dendrites. Mathematical Biosciences 188, 147–155 (2004).
 [13] Van Pelt, J. & Uylings, H. B. M. Modeling Neuronal Growth and Shape. In Laublicher, M. D. & Müller, G. B. (eds.) Modeling Biology – Structures, Behaviors, Evolution, 195–215 (The MIT Press, Cambridge, Massachusetts, 2007).
 [14] Koene, R. A. et al. NETMORPH: A framework for the stochastic generation of large scale neuronal networks with realistic neuron morphologies. Neuroinformatics 7, 195–210 (2009).
 [15] Cuntz, H., Forstner, F., Borst, A. & Häusser, M. One rule to grow them all: A general theory of neuronal branching and its practical application. PLoS Computational Biology 6, e1000877 (2010).
 [16] Cuntz, H., Mathy, A. & Hausser, M. A scaling law derived from optimal dendritic wiring. Proceedings of the National Academy of Sciences 109, 11014–11018 (2012).
 [17] Luczak, A. Measuring neuronal branching patterns using modelbased approach. Frontiers in Computational Neuroscience 4, 135 (2010).
 [18] Luczak, A. Spatial embedding of neuronal trees modeled by diffusive growth. Journal of Neuroscience Methods 157, 132–141 (2006).
 [19] Gillette, T. A. & Ascoli, G. A. Topological characterization of neuronal arbor morphology via sequence representation: II  global alignment. BMC Bioinformatics 16, 209 (2015).
 [20] Mohan, H. et al. Dendritic and Axonal Architecture of Individual Pyramidal Neurons across Layers of Adult Human Neocortex. Cerebral Cortex 25, 4839–4853 (2015).
 [21] Li, A. et al. Microoptical sectioning tomography to obtain a highresolution atlas of the mouse brain. Science 330, 1404–1408 (2010).
 [22] Ragan, T. et al. Serial twophoton tomography for automated ex vivo mouse brain imaging. Nature Methods 9, 255–258 (2012).
 [23] Jefferis, G. S. & Livet, J. Sparse and combinatorial neuron labelling. Current Opinion in Neurobiology 22, 101–110 (2012).
 [24] Chung, K. & Deisseroth, K. CLARITY for mapping the nervous system. Nature Methods 10, 508–513 (2013).
 [25] Zhou, H. et al. Dense reconstruction of brainwide neuronal population close to the ground truth. bioRxiv 223834 (2018).
 [26] Quan, T. et al. NeuroGPSTree: automatic reconstruction of largescale neuronal populations with dense neurites. Nature Methods 13, 51–54 (2016).
 [27] Vormberg, A., Effenberger, F., Muellerleile, J. & Cuntz, H. Universal features of dendrites through centripetal branch ordering. PLoS Computational Biology 13, e1005615 (2017).
 [28] Gillette, T. A. & Ascoli, G. A. Topological characterization of neuronal arbor morphology via sequence representation: I  motif analysis. BMC Bioinformatics 16, 216 (2015).
 [29] Molnar, P. On geometrical scaling of Cayley trees and river networks. Journal of Hydrology 322, 199–210 (2006).
 [30] Sholl, D. A. Dendritic organization in the neurons of the visual and motor cortices of the cat. Journal of Anatomy 87, 387–406 (1953).
 [31] Van Pelt, J., Uylings, H. B., Verwer, R. W., Pentney, R. J. & Woldenberg, M. J. Tree asymmetryA sensitive and practical measure for binary topological trees. Bulletin of Mathematical Biology 54, 759–784 (1992).
 [32] Samsonovich, A. V. & Ascoli, G. A. Morphological homeostasis in cortical dendrites. Proceedings of the National Academy of Sciences of the United States of America 103, 1569–1574 (2006).

[33]
Samsonovich, A. V. & Ascoli, G. A.
Statistical determinants of dendritic morphology in hippocampal pyramidal neurons: A hidden Markov model.
Hippocampus 15, 166–183 (2005).  [34] Ostilli, M. Cayley Trees and Bethe Lattices: A concise analysis for mathematicians and physicists. Physica A: Statistical Mechanics and its Applications 391, 3417–3423 (2012).
 [35] Albert, R. & Barabási, A.L. Statistical mechanics of complex networks. Reviews of Modern Physics 74, 47–97 (2002).
 [36] Donohue, D. E. & Ascoli, G. A. A comparative computer simulation of dendritic morphology. PLoS Computational Biology 4, e1000089 (2008).
 [37] Zeng, H. & Sanes, J. R. Neuronal celltype classification: Challenges, opportunities and the path forward (2017).
 [38] Joo, W. J., Sweeney, L. B., Liang, L. & Luo, L. Linking cell fate, trajectory choice, and target selection: Genetic analysis of sema2b in olfactory axon targeting. Neuron 78, 673–686 (2013).
 [39] Han, Y. et al. The logic of singlecell projections from visual cortex. Nature 556, 51–56 (2018).
 [40] Schubert, D., Kötter, R., Luhmann, H. J. & Staiger, J. F. Morphology, electrophysiology and functional input connectivity of pyramidal neurons characterizes a genuine layer Va in the primary somatosensory cortex. Cerebral Cortex 16, 223–236 (2006).
 [41] Marx, M. & Feldmeyer, D. Morphology and physiology of excitatory neurons in layer 6b of the somatosensory rat barrel cortex. Cerebral cortex 23, 2803–2817 (2013).
 [42] Gollo, L. L., Kinouchi, O. & Copelli, M. Active dendrites enhance neuronal dynamic range. PLoS Computational Biology 5, e1000402 (2009).
Acknowledgements
CL acknowledges financial support from National foundation of Science in China (NFSC, grant No. 11701201 and No. 11871061), YZ acknowledges financial support from by NFSC (grant No. 11701200 and 11871262) and TQ acknowledges financial support from Science Fund for Young and Middleaged Creative Research Group of the Universities in Hubei Province (Grant No. T201520) and NSFC (grant No. 81771913). We also thank the Optical Bioimaging Core Facility of WNLOHUST for the support in data acquisition, and the Analytical and Testing Center of HUST for spectral measurements. We thank Prof. Shaoqun Zeng for helpful discussions on actual neuron data.
Author contributions statement
YH, CL and TQ performed the analysis, CL and YH performed modelling, and YZ and CL conceived and designed the project. All authors reviewed the manuscript.
Additional information
Competing interests: The authors declare no competing interests.
Data availability
The data that support the findings of this study are available from the corresponding author on request.
Comments
There are no comments yet.