We define a new balance index for rooted phylogenetic trees based on the symmetry of the evolutive history of every set of 4 leaves. This index makes sense for multifurcating trees and it can be computed in time linear in the number of leaves. We determine its maximum and minimum values for arbitrary and bifurcating trees, and we provide exact formulas for its expected value and variance on bifurcating trees under Ford's α-model and Aldous' β-model and on arbitrary trees under the α-γ-model. of isomorphisms of the restriction of T to them (the rooted quartet they define), and then we add up these values over all 4-tuples of different leaves of T . The idea behind the definition of this balance index is that a highly symmetrical evolutive process should give rise to symmetrical evolutive histories of many small subsets of taxa. In terms of phylogenetic trees, this leads us to expect that, the most symmetrical a phylogenetic tree is, the most symmetrical will be its restrictions to subsets of leaves of a fixed cardinality. Since the smallest number of leaves yielding enough different tree topologies to allow a meaningful comparison of their symmetry is 4, we assess the balance of a tree by measuring the symmetry of all its rooted quartets and adding up the results. And indeed, in Section 4 below we shall find the trees with maximum and minimum values of our rooted quartet index in both the arbitrary and the bifurcating cases, and it will turn out that the minimum value is reached exactly at the combs (see Fig. 1.(a)), which are usually considered the least balanced trees, and the maximum value is reached, in the arbitrary case, exactly at the rooted stars (see Fig. 1.(b)) and, in the bifurcating case, exactly at the maximally balanced trees (cf. Fig. 3 ), which in both cases are considered the most balanced trees.Besides taking its maximum and minimum values at the expected trees, other important features of our index are that it can be easily computed in linear time and that its mean value and variance can be explicitly computed on any probabilistic model of phylogenetic trees satisfying two natural conditions: independence under relabelings and sampling consistency. This allows us to provide these values for two well-known probabilistic models of bifurcating phylogenetic trees, Ford's α-model [13] and Aldous' β-model [2], which include as specific instances the Yule [14,29] and the uniform [6,24,19] models, as well as for Chen-Ford-Winkel's α-γ-model of multifurcating trees [7]. To our knowledge, this is the first shape index for which closed formulas for the expected value and the variance under the α-γ-model have been provided.The rest of this paper is organized as follows. In the next section we introduce the basic notations and facts on phylogenetic trees that will be used in the rest of the paper, and we recall several preliminary results on probabilistic models of phylogenetic trees, proving those results for which we have not been able to find a suitable reference in the literature. Then, in Section 3, we ...