New Avenues for Old Travellers: Phenotypic Evolutionary Trends Meet Morphodynamics, and Both Enter the Global Change Biology Era

Evolutionary trends (ETs) are traditionally defined as substantial changes in the state of traits through time produced by a persistent condition of directional evolution. ETs might also include directional responses to ecological, climatic or biological gradients and represent the primary evolutionary pattern at high taxonomic levels and over long-time scales. The absence of a well-supported operative definition of ETs blurred the definition of conceptual differences between ETs and other key concepts in evolution such as convergence, parallel evolution, and divergence. Also, it prevented the formulation of modern guidelines for studying ETs and evolutionary dynamics related to them. In phenotypic evolution, the theory of morphodynamics states that the interplay between evolutionary factors such as phylogeny, evo-devo constraints, environment, and biological function determines morphological evolution. After introducing a new operative definition, here we provide a morphodynamics-based framework for studying phenotypic ETs, discussing how understanding the impact of these factors on ETs improves the explanation of links between biological patterns and processes underpinning directional evolution. We envisage that adopting a quantitative, pattern-based, and multifactorial approach will pave the way to new potential applications for this field of evolutionary biology. In this framework, by exploiting the catalysing effect of climate change on evolution, research on ETs induced by global change might represent an ideal arena for validating hypotheses about the predictability of evolution.


History and Future of Evolutionary Trends
Evolutionary trends (ETs) are traditionally defined as persistent and directional changes in the state of one or more quantitative traits, resulting in substantial changes through time and representing the primary phenomenon characterising evolution at high taxonomic levels and over long-time scales (Gould, 2002;McNamara, 2006;McShea, 2005). Alroy (2000) described the study of ETs as "one of the oldest and more intriguing topics in evolutionary biology", agreeing with McKinney (1990), who indicated the concept of ET as "arguably the single most important in the study of evolution".
The quest for ETs started early in macroevolution and quickly assumed crucial importance in phenotypic evolution and many other fields of evolutionary biology. The first groundbreaking investigations concerning phenotypic directional evolution led to the description of iconic and pervasive trends throughout the history of a single clade, frequently relying on the evidence available from fossil record (e.g., increasing hypsodonty and gradual acquisition of the monodactyl posture in equids caused by the spread of grasslands - Kowalewsky, 1874). Then, the formulation of the so-called 'biological rules' extended the definition of ETs to include directional responses to ecological, climatic 1 3 or biological gradients, with several clades being simultaneously analysed obtaining the same evolutionary pattern (e.g., latitudinal or elevation gradients, like in Allen's or Bergmann's rules - Allen, 1877;Bergmann, 1847).
During the second half of the twentieth century, an animated debate rose regarding how to identify an ET. In particular, a key point in this debate was whether to find out a directional pattern in evolution suffices to establish the presence of an ET, or if an explanation about the underpinning biological mechanisms is needed (Dayan & Simberloff, 1998;Mayr, 1956;Meiri, 2011). Further theoretical battlegrounds concerned erroneous interpretations of processes involved in the occurrence of directional evolution. In the past, these misinterpretations led to gross oversimplifications of the complexity of evolutionary dynamics and produced "common portrayals of evolution in non-academic settings including outright notions of 'advancement', that held sway within evolutionary biology as well" (Gregory, 2008, p. 259). For example, Edward D. Cope based his concept of kinetogenesis (i.e., theory stating that movements of animals aided in the alteration and development of moving parts of the body) on the idea that key evolutionary novelties originate from generalized phenotypes, in keeping with Lamarck's view of evolution as possessing intrinsic directionality (Bowler, 1977). Whereas kinetogenesis (and its allied concept of orthogenesis) became extinct, the idea that evolution may possess a weak form of directionality (i.e., in contrast to old views of evolution as a movement towards 'perfection of life') in specific circumstances recently found support (Raia & Fortelius, 2013;Raia et al., 2016). This vision relies on the concept that specialisation and adaptation towards optimal phenotypes can only occur after the birth of clades, for example, originating directional trends through time (Ayala, 1988;Gould, 1988a;Raia & Fortelius, 2013;Raia et al., 2016;Rosenzweig & McCord, 1991). Similar weak forms of directionality were suggested to occur in specific fields of evolutionary biology, such as the evolution of organismal diversity and complexity. For example, the zero-force evolutionary law (McShea & Brandon, 2010;McShea et al., 2019) stated that "in any evolutionary system in which there is variation and heredity, there is a tendency for diversity and complexity to increase, one that is always present but may be opposed or augmented by natural selection, other forces or constraints acting on diversity or complexity" (McShea et al., 2019(McShea et al., , p. 1103.
Although the search for new ETs and the validation of the existing ones remain central topics in evolutionary biology (Cardini, 2019;Cardini & Polly, 2013), the formulation of a well-supported and operative definition of ET remains an unsolved issue. The lack of a broadly accepted operative definition prevented evolutionary biologists from unanimously pinpointing the requirements needed for detecting and validating ETs and blurred the definition of conceptual differences between ETs and other independent/ nested key concepts in evolution, whose operative definition had already been formulated (e.g., convergence - Stayton, 2015aStayton, , 2015b. The pivotal importance of ETs in evolutionary biology can be understood by considering many possible contexts in which they represent a central tenet. For instance, ETs constitute an ideal case study to describe and separate the two components of evolution proposed by Simpson (1944), the tempo and mode in evolution, and their mutually independent variations (Felice et al., 2018;Michaud et al., 2018). The strength of an ET can be quantified as the magnitude of a vector, providing a practical way of representing the speed (or tempo) of evolution. The direction of the same vector represents the mode of evolution (i.e., observed pattern of variation). Phenotypic ETs can also be investigated in the framework of evolutionary landscapes, accounting for variations in evolutionary constraints affecting a clade through time or in response to environmental shifts. Phenotypic ETs were frequently found as the resulting evolutionary outcome in studies about the controversial field of evolutionary predictability, that focuses on the occurrence of repetitive and foreseeable patterns in evolution under specific conditions (e.g., recurrent patterns in insular systems- de Visser & Krug, 2014;Whittaker et al., 2017). For instance, directional evolution towards a relative reduction of flight muscles in insular colonies of birds was found to be a predictable and frequent phenotypic specialization produced by an ecological release from mammalian predators (Wright et al., 2016). And it is specialisation that was suggested by recent conservation studies to serve as a fundamental proxy for vulnerability and extinction rate increases of species in future or hypothetical scenarios, supporting the idea that research on phenotypic ETs is conceivable to become additional topics included in the future scientific debate on climate change and conservation biology (Gallagher et al., 2015).
Our main goal is to propose a modern theoretical framework for the study of phenotypic ETs, based on Seilacher's theory of morphodynamics (Briggs, 2017;Seilacher & Gishlick, 2015), that explores the interactions between ETs and the main evolutionary factors involved in their occurrence. The theory of morphodynamics can be visually represented with a tetrahedron whose nodes correspond to the main evolutionary factors influencing morphological evolution: phylogenetic history, developmental constraints, environment, and biological function (Briggs, 2017). Thanks to the inclusion of the environmental component, morphodynamics might be considered an expansion of the concept of 'constructional morphology' proposed by Seilacher (1970Seilacher ( , 1991, which considered the role of phylogenetic history, developmental constraints, and biological function in determining morphological evolution of organisms. Both morphodynamics and constructional morphology "replace the idealistic concept of 'bauplan', in which the concept of organisms as integrated, functional machines was substituted for the sterile notion of organisms as systems of pure morphological form" (Bock, 1991, p. 27).
In order to achieve our goal, after introducing a new operative definition of ETs, our discussion will start by investigating the influence of environment and evo-devo constraints on phenotypic ETs, focusing on the benefits resulting from the adoption of spatially-explicit models and on the role played by the plasticity-first evolution model. Then, we will elaborate on the importance of including phylogenetic information by adopting phylogenetic comparative methods to evaluate evolutionary patterns. Finally, we will describe how recent advances in morphological quantification can be applied to investigate the relationship between (directional) morphological evolution and biological function. In order to suggest new frontiers of research, in the last section we will discuss how, by exploiting the catalysing effect of climate change on evolution, research on directional evolution induced by global change might represent an ideal arena for the validation of hypotheses about the predictability of evolution.

Geometric Definition and Classification of Evolutionary Trends
The description of an ET, traditionally considered as an unambiguous phenomenon, has often been limited to anecdotal lists including self-evident examples, without any further quantitative consideration (Stayton, 2015b). The absence of a widely-supported operative definition led to the formulation of multiple (and partially overlapping or conflicting) theoretical definitions. Furthermore, this absence contributed to generate confusion about the relationship between evolutionary concepts like ETs, convergence, divergent evolution, and parallelism, that were considered nested or independent evolutionary cases according to different authors (Manceau et al., 2010;Scotland, 2011).
The concept of '(non) parallel evolution' (i.e., continuum from convergent through parallel to divergent evolution) was recently proposed by Bolnick et al. (2018) to describe the independent evolution of replicate populations, in reference to bacterial cultures of laboratory experiments. The aim of this term is to adopt a pattern-based (Box 1), geometric, and quantitative definition to include in a single concept all classes of evolutionary patterns existing whenever two or more replicates evolving toward a constant direction in a trait space are taken into account. According to this approach, the distinction between parallelism, convergence, and divergent evolution entirely relies on the amplitude of the angles among these evolutionary trajectories, represented as vectors in the trait space. The transposition of (non) parallel evolution to the macroevolutionary scale ( Fig. 1) represents a perfect geometric definition of ETs. This definition constitutes a key innovation to solve the main theoretical issue related to the study of ETs, that is the description of the existing classes of ETs (and their relationships) occurring whenever multiple clades are simultaneously considered (Stayton, 2015a). Three possible classes (i.e., parallel, convergent, and divergent evolution) are conceivable from this perspective. The criterion of distinction between them relies, even in this case, on the orientation of the evolutionary trajectories shown by the single groups ( Fig. 1B-D). Therefore, parallelism will be the resulting outcome when evolutionary trajectories point in the same direction, otherwise, a condition of convergence or divergent evolution will occur, depending on whether trajectories respectively point 'at' or 'away from' a region of the trait space.
This operative definition of ETs is in line with the idea of 'directional evolution' provided by Hunt and Carrano (2010). The authors described directional evolution, Brownian Motion (BM) evolution, and stasis as the most relevant models in long-term (i.e., paleontological) macroevolutionary studies, even stating that only mean and variance of the distribution of trait values matter over paleontological temporal scales. In particular, Hunt and Carrano (2010) distinguished between directional and BM evolution arguing that whenever deviations from the mean of this distribution occur, the resulting model will be defined as directional evolution; otherwise, we will be in presence of BM evolution. Coherently with this distinction, several methods using BM as a null hypothesis were proposed in order to detect specific ETs, such as convergent evolution (e.g., Wheatsheaf, C1 and θ metrics- Arbuckle et al., 2014;Stayton, 2015a;Castiglione et al., 2019a) and directional trends through time (e.g., Castiglione et al., 2019b;Sherratt et al., 2016). Nevertheless, in the presence of constraints, BM-like diffusions might result in a directional pattern of evolution determined by the position of constraints in the trait space (Fisher, 1986). For example, the presence of a lower boundary in size evolution has been supposed to be the mechanism that accounts for Cope's rule (i.e., tendency for clades to increase their size through time - Stanley, 1973). According to this vision, Cope's rule is simply produced by an increase in size variance through time (Gould, 1988b). Based on these considerations, McShea (1994McShea ( , p. 1747) distinguished between passive and driven ETs, stating that, "in a driven trend, the distribution mean increases on account of a force (which may manifest itself as a bias in the direction of change) that acts on lineages throughout the trait space in which diversification occurs. In a passive system, no pervasive force or bias exists, but the mean increases because change in one direction is blocked by a boundary, or other inhomogeneity, in some limited region of the trait space". McShea also validated the occurrence of driven and passive ETs in the evolution of horses, brachiopods, and rodents, and tried to classify as driven or passive trends several other live hypotheses, such as trends through time in size, complexity, and energy intensiveness (McShea, 1998). Our operative definition of ETs can be applied to both passive and driven ETs.
According to our definition of ETs, parallel evolution, convergence, and divergent evolution can be considered as patterns produced whenever, analysing multiple clades at the same time, a portion of the groups undergoes an episode of directional evolution, whereas other groups are subject to a different evolutionary regime. The adoption of a multiscale approach is therefore recommended to recognise these classes of ETs. Performing evolutionary analyses at different taxonomic levels has also been suggested as a promising avenue to define the evolutionary boundaries acting on a clade more rigorously and the relationships between morphology, evolutionary tempo, and ecological divergence in a broad comparative context (Tamagnini et al., 2017). In this sense, research on ETs would take a major step forward if a pattern-based approach, like the one described above (that can only validate or disprove the presence of a pattern), was systematically accompanied by multifactorial analyses. These analyses would also allow researchers to take into account data relative to underpinning dynamics leading to the occurrence of directional evolution (e.g., ecological variables or life-history traits). When it comes to evolutionary dynamics occurring in phenotypic evolution, the tangled nature of interactions between the factors underpinning morphological evolution has already been explored by the theory of morphodynamics (Briggs, 2017, pp. 203-204), concluding that four components have a predominant role: phylogenetic history, developmental constraints, environment, and biological function. However, the central role of genetic constraints (e.g., scarce genetic variation, low mutation rate, and multivariate genetic correlations) was subsequently recognized as a limiting factor for phenotypic evolution (Futuyma, 2010). Given the role of evo-devo processes as a bridge between genetic and developmental constraints (Müller, 2007), it is in our opinion advisable to adopt them as the fourth factor included in the present theoretical framework (Fig. 2), replacing the factor 'developmental constraints' originally included in the theory of morphodynamics. In the next sections, we will explore the interactions between phenotypic ETs and these evolutionary factors (starting with evo-devo constraints and environment, then shifting onto phylogeny and finally focusing on biological function in the light of new perspectives in morphological quantification), also mentioning applications or methodological approaches that can be used in these contexts.

Process-Based Versus Pattern-Based Approaches and Implications with Evo-Devo Mechanisms
The necessity to identify an underpinning evolutionary process to validate the presence of evolutionary patterns has been criticised and gave rise to a long-standing debate in the history of evolutionary biology (Mayr, 1956). When it comes to ETs, the emerging solutions to this dispute are all stemming from two alternative approaches to individuate the presence of a phenomenon that led to the formulation of process-based and patternbased definitions. Process-based definitions require the a priori demonstration that a common mechanism is acting in response to the same evolutionary conditions in all considered groups. By contrast, pattern-based definitions bind the validation of an evolutionary pattern to the mere fulfilment of mathematical or geometrical criteria, generally based on ratios or angular values: it follows that this type of approach is always coupled with quantitative measures and null models (Rosenblum et al., 2014;Stayton, 2015a).
The improved understanding of the functioning of evodevo mechanisms under different evolutionary conditions initiated an identification process of the linkage between phenotypic variation and specific evo-devo constraints. For instance, Felice and colleagues (2018) hypothesised that phenotypic integration and modularity act as a bridge between developmental, genetic, and functional levels of trait associations, influencing the morphology of organisms, even catalyzing the evolution and the production of biodiversity. Under suitable conditions, phenotypic integration canalises trait variations along limited evolutionary directions, meeting one of the preconditions to the rise of an ET (Goswami et al., 2014;Wagner & Zhang, 2011). Variations in consolidated patterns of integration between morphological modules were also observed to play an essential role in the iterative occurrence of extreme phenotypes, like the skull of saber-toothed cats, classifiable as episodes of convergent evolution (Meloro & Slater, 2012). The recent advances in mechanism investigation might tempt evolutionary biologists to adopt process-based definitions of ETs in every situation, formulating a developmental or physiological explanation for an observed directionality simply assuming its validity on the basis of different and previously observed case studies. However, despite this type of definitions frequently represents the optimal choice in microevolution, the validity of evo-devo processes in broad comparative data still has to be extensively proved, and process-based approaches were frequently demonstrated to produce contradictions or misinterpretations in the explanation of evolutionary patterns when considering large taxonomic scales (Felice et al., 2018;Meiri, 2011;Stayton, 2015a). Furthermore, process-based definitions systematically fail to categorize the occurrence of a specific pattern, like the presence of directionality, whenever an evolutionary mechanism in place cannot be easily found (e.g., Collar et al., 2014).

3
The advantage of adopting pattern-based definitions in macroevolution partially resides in their extreme flexibility, allowing researchers to accommodate additional modifiers easily and clearly define nested cases within an evolutionary phenomenon (Stayton, 2015a). But more importantly, pattern-based validations are suitable for making inferences about different systems in broad comparative data and can serve to steer post hoc studies on the mechanisms responsible for the observed patterns. Inferences based on a process-based approach would require, by contrast, an a priori validation of the putatively involved mechanism for each different system, unless resulting in gross extrapolations beyond the conditions of the system (e.g., Menge, 1992). Obtaining this sort of process corroboration is challenging due to its experimental nature, which limits the number and the scale of analysable samples. The necessity to shift from small scale understanding to large scale predictions is shared in other fields of evolutionary biology, like macroecology, where there is an increasing demand to identify general, predictive, and empirical relationships in natural systems (Currie, 2019).

The Role of Evo-Devo Constraints and Environment in Shaping Directional Macroevolutionary Patterns and Processes: Plasticity-First Evolution and Spatially-Explicit Models
Investigating the interconnection between evo-devo constraints, environment, and phenotypic variations is crucial for better understanding macroevolutionary processes and improving the prediction of patterns produced by evolutionary directionality. The recent formulation of the evolutionary mechanism named plasticity-first evolution (PFE) model (Levis & Pfennig, 2016) has emphasised the influence of environment and its interactions with evo-devo constraints in shaping phenotypes. According to the PFE model, the environmentally initiated phenotypic change, known as phenotypic plasticity, acts as a precursor of evolutionary adaptation. The environmental variability uncovers the cryptic genetic variation shaping the extant phenotypic plasticity. Then, genetic accommodation leads to the adaptive refinement of favoured phenotypes whenever a condition of diffused polyphenism (i.e., condition where different phenotypes can easily coexist) is not favoured by selection. Shifting the focus on the resulting macroevolutionary effects, directional and persistent environmental stimuli might underpin the rise of an ET acting on preexisting phenotypic plasticity. Despite concerns regarding the PFE model being raised due to its (at least apparently) partial inconsistency with the modern synthesis of evolutionary biology (i.e., phenotypic plasticity acting in specific cases as a precursor of evolutionary adaptation), leading to define it as a mere artefact of lab studies (Ho & Zhang, 2018), Levis and colleagues (2018) recently validated the presence of PFE in North American spadefoot toads of the genus Spea demonstrating that a persistent environmental stimulus (i.e., presence of different preys and competitors) caused the evolution of intraspecific dietinduced plasticity into fixed morphs within the genus.
Leaving aside the uncertainties relative to the underlying evolutionary mechanisms, the focus on geographical and ecological information in a theoretical model of directional evolutionary patterns recently became a central topic in macroevolution. For instance, Polly (2018) highlighted that existing statistical models used to study evolutionary patterns rely on a Fisherian view of evolution, whereby species are considered single panmictic populations (i.e., populations in which all individuals are potential partners) whose traits are the same regardless of spatial and environmental influences (Fisher, 1930). According to this vision, evolutionary outcomes disregard drift and habitat-specific selection (and are likely to be affected by a misjudgement of the impact of sexual selection). However, both drift and habitat-specific selection might produce, under neutral evolution, phenotypic paths of evolution similar to those resulting from the occurrence of directionality under Fisherian conditions, like the rise of morphological gradients along species distribution ranges (e.g., cranial shape evolution of European common shrews- Polly, 2018). To this aim, Polly also used computational modelling to produce simulations, supporting the idea that drift and habitat-specific selection represent a potential source of misinterpretations in macroevolution. The inclusion of spatially structured variation in evolutionary models (i.e., Wrightian view of evolution) would allow considering metapopulations composed by interacting demes characterised by different trait values: this would enable the distinction between patterns produced by spatial processes, like drift and habitat-specific selection, and ETs associated with non-neutral and directional evolutionary regimes (Hanski, 1999). A potential solution would be the development of individual-based macroevolutionary models. However, a major obstacle to go in this direction is the discrepancy between the lineagebased focus, typical of evolutionary theories, and the attention for individual organisms and their interactions over timescales of a few generations, which often belongs to ecology (Rosindell et al., 2015). In spite of the temporary absence of available spatially explicit models to distinguish between Fisherian and Wrightian trait evolution, ecological models based on the neutral theory of biodiversity represents a key tool to bridge the gap between individualbased ecology and macroevolution, as well as the inclusion of spatial gradients of selection, metacommunity spatial dynamics, and trait distributions of fossil record (Badgley & Finarelli, 2013;Lyons & Smith, 2013;Maestri et al., 2018;Polly, 2018).

Phylogenetic Comparative Methods and the Necessity to Adopt a Phylogenetically-Informed Approach for the Study of Phenotypic ETs
Macroevolutionary datasets should be collected following a rigorous sampling design. However, evolutionary biologists face issues related to accessibility, incompleteness or sample size, often leading to the implementation of simplified sampling designs that suffer from unknown biases (Albert et al., 2010). Considering taxonomically related groups of species that occur together in space, instead of exclusively analysing monophyletic groups, is often referred to as assemblage approach and is a widespread practice in evolutionary biology, specifically for the study of phenotypic ETs (Meiri, 2011;Stroud et al., 2015). This approach shows, whenever applied to evolutionary biology, the same shortcomings occurring in its original field of application, ecology. For instance, Mittelbach and Schemske (2015) highlighted that assemblagebased analyses can be misleading without a preliminary recognition of the dynamic nature of species pools and the processes underlying species pool formation. The resulting almost inevitable ineffectiveness to describe evolutionary dynamism over large spatial and/or temporal scales originates discrepancies frequently found in assemblage-based studies on ETs (e.g., inconsistent presence of Bergmann's rule in different assemblages of the mammalian order Carnivora- Diniz-Filho et al. 2007) and suggests to always adopt a phylogenetically-informed approach for the study of macroevolutionary patterns and processes, specifically relying on phylogenetic comparative methods (PCMs).
PCMs are statistical models that estimate the evolutionary regime that best approximates the tempo and mode of evolution acting on the considered traits, allowing researchers to correct for biases due to the non-independence of sampled observations in macroevolutionary samples (i.e., phylogenetic relationships - Felsenstein, 1985;Venditti et al., 2011). Although discussing existing PCMs and their assumptions is beyond the aim of the present work, the existing literature about this topic is already impressive (e.g., Adams & Collyer, 2019;Cooper et al., 2016 and references therein). Taking into account paleontological data, the possibility to include the available fossil evidence in phylogenetic reconstructions was recently proved to improve estimations of the best model of evolution, facilitating the detection of changes in tempo and mode of evolution and the reconstruction of the ancestral states, preventing, therefore, erroneous inferences and increasing the statistical power of PCMs (Mitchell et al., 2019;Schnitzler et al., 2017). The inclusion of paleontological datasets in complex modelling of temporal and spatial evolutionary patterns and the application of PCMs to ecological information represent a promising way to proceed for evolutionary biologists in the foreseeable future and would facilitate the detection of ETs and their variations in strength and direction through time (Castiglione et al., 2019a(Castiglione et al., , 2019bMaestri et al., 2018).

New Perspectives in Morphological Quantification: Evolutionary Landscapes, Geometric Morphometrics, and Finite Element Analysis
The tangled interactions between variations in selective advantage deriving from biological function and phenotypic or genotypic responses have been traditionally rendered using evolutionary landscapes since this type of graph is meant to show differences in the level of selective advantage intrinsic to a specific set of trait values (Arnold et al., 2001;Wright, 1932). Nowadays, several versions of evolutionary landscapes are available in phenotypic evolution to investigate these interactions both at clade and single-lineage level, enabling, for instance, to recognize the occurrence of evolution driven by ecological opportunity (Hunt et al., 2015;Caetano and Harmon, 2017;Voje, 2020). Regardless of the adopted version of landscape, these graphs are particularly indicated for the visual detection of ETs, since their conformation is likely to reflect the acting evolutionary pattern. For instance, in the absence of boundaries (i.e., driven ETs sensu McShea, 1994), favoured phenotypes are distributed along a linear line whenever directionality is common to the entire sample. The resulting evolutionary landscape is a rectilinear ridge (Fig. 3A). Parallelism, convergence and divergent evolution originate a series of parallel ridges, a peak or a depression into the landscapes (Fig. 3B-D), respectively. By contrast, non-directional evolution, as in the case of BM evolution, is likely to generate a flat surface (Fig. 3E). In the presence of boundaries (i.e., passive ETs sensu McShea, 1994), their position in the evolutionary landscape contributes to determining the resulting evolutionary pattern.
Recently developed techniques of morphological quantifications (like geometric morphometrics-GMM-and finite element analysis-FEA) allow researchers to quantitatively describe different aspects of phenotypic variation such as size, shape, and functional performance. However, these techniques produce highly multivariate data generally difficult to be implemented for an adequate evolutionary landscape or to be corrected applying PCMs. For instance, Adams and Collyer (2018a) recently demonstrated that the most common methods (i.e., log-likelihood or AIC score-based techniques) used to estimate the best evolutionary model in multivariate GMM analyses become more ill-conditioned as the ratio between considered morphological traits and sampled species increases or as considered models become more complex. Furthermore, paleomorphological studies have to face issues linked to the scarcity and the fragmented nature of preserved material, plus distortions of fossils produced by taphonomic processes (Arbour & Currie, 2012). To overcome these issues, a new generation of PCMs is rising in GMM, allowing operators to ground estimations of the best model and its parameters on methods unbiased by an increase in trait dimensionality (i.e., statistical tests relying on traces of covariance matrices- Adams & Collyer, 2018a), to adopt innovative techniques of randomizing residuals in phylogenetic ANOVAs and regression models (Adams & Collyer, 2018b) or to graphically visualise potential trends in phylogenetic signal (i.e., phylogenetically aligned component analysis- Collyer & Adams, 2020). Furthermore, recent studies demonstrated how FEA analyses can be employed, relying on trait spaces derived from GMM, to obtain quantitative surfaces describing the performances of different morphologies (O'Higgins et al., 2011;Polly et al., 2016). The combination of multiple performance surfaces (e.g., optimisation concepts such as Pareto front) represents an informative tool for the construction of phenotypic evolutionary landscapes that show optimal morphologies able to contemporarily carry out all functions associated with the selected performances (see Jones et al., 2021 for a worked example about the synapsid-mammal locomotor transition). Combining these techniques of morphological quantification was also proved to be useful for retrodeforming altered fossil specimens or reconstructing damaged remains, as well as for hypothesising and analysing non-available intermediate morphologies (Gunz et al., 2009;Schlager et al., 2018;Tseng, 2013). The possibility to consider non-preserved or purely hypothetical evolutionary forms allows researchers to investigate the occurrence of phenotypic ETs filling unexplored regions of morphological trait spaces, obtaining highly detailed landscapes and characterising the response of functional performance produced by different types of continuous variables, ranging from ecological to physiological ones. A further step forward in eco-evolution is the concept of dynamic adaptive landscape proposed by Laughlin and Messier (2015), whose primary focus is detecting changes occurring in an evolutionary landscape along environmental gradients. This technique uses multivariate trait distributions to determine functional trade-offs at different taxonomic scales, improving the understanding of functional trait evolution and its link to the environmental conditions, and therefore represents a powerful tool for studying the interactions between biological function and environment. Adopting a dynamic version of phenotypic evolutionary landscapes (potentially also including paleontological evidence) seems to be a promising target, on the one hand, to visualize selective advantage variations along stratigraphic sequences and, on the other, to graphically show the initial outbreak, the constant transformations, and the final weakening of an ET through time (Fig. 4).

The Search for Phenotypic Evolutionary Trends in Climate Change Scenarios as a Tool for Understanding the Predictability of Evolution
As discussed above, pattern-based methods accompanied by multifactorial approaches deriving from a morphodynamicsbased theoretical framework will pave the way to new potential applications for studies on phenotypic ETs. Research on ETs in the context of global changes can represent a timely and promising endeavour, particularly considering the predictability of evolution.
Clarifying under what conditions evolution proceeds along predictable pathways would enable researchers to increase and refine the use of evolutionary biology in several applied contexts, for instance, providing tools for Fig. 4 Example of dynamic evolutionary landscape including the case study of fossil and living equids. A typical scenario for the occurrence of ETs is represented by an initial condition of null selective advantage differences within a given trait space (that is likely to result in Brownian Motion evolution and a flat evolutionary landscape-t 0 ) that is disrupted by variations in an external factor (e.g., environmental changes). These variations might lead to a temporary phase of directional evolution (potentially reflected by changes in the configuration of the evolutionary landscape-t 1 and t 2 ). The disappearance of the external forcing can terminate this phase (e.g., end of environmental changes and return of a flat evolutionary land-scape configuration-t 3 ). A similar evolutionary scenario might have occurred in limb evolution of equids: the environmental transition from forests to grasslands started during the late Oligocene is likely to have induced an event of directional evolution towards the acquisition of the monodactyl posture within this clade. Before and after this phase, the absence of ETs (e.g., Brownian Motion evolution) can be considered a likely evolutionary outcome, even in the light of the presence of more stable environmental conditions. Yellow dots represent different species of equids that occurred at different (Color figure online) sustainability assessments and reproducibility of laboratory experiments (Rego-Costa et al., 2018). The first pioneering research in this field found strong sources of stochasticity, as erratic fluctuations of the environment or genetic drift, as limiting factors of evolutionary predictability (Saether & Engen, 2015). By contrast, investigating the importance of deterministic chaos under simulated repetitive evolutionary conditions, Rego-Costa and colleagues (2018) recently pointed out that a strong forcing by a changing environment (e.g., rapid shifts in selective optima) can improve the predictability of eco-evolutionary dynamics, otherwise resulting in a chaotic scenario under constant evolutionary conditions. It follows that climate change events, whenever causing directional shifts of selective conditions, might represent an enhancer of evolutionary predictability.
The current phase of anthropogenic climate change, resulting in a pattern of diffused global warming, is threatening the existing biodiversity and is potentially leading the entire planet towards a sixth mass extinction that might require millions of years to restore comparable levels of species richness (Barnosky et al., 2011;Davis et al., 2018;Rosenzweig et al., 2008). These environmental changes force species to shift their geographical distribution, change their phenology or alternatively develop adaptations in their morphology and/or physiology relying on an admixture of genetic change and plasticity (Bellard et al., 2012;Hoffmann & Sgrò, 2011). The rise of a phenotypic ET, frequently regarding variations in organism size or shape and changes in the integration between biological modules, constitutes a frequent and highly repeatable event induced by climate change, probably because reflecting a phenotypic response to a constant and directional change in the environment that fits a linear model (MacLean et al., 2018 and references therein). For this reason, research on the occurrence of phenotypic ETs in currently changing environments has the potential to become, in the near future, an indispensable tool for clarifying the dynamics underpinning evolutionary predictability, allowing researchers to verify hypotheses (by the exploitation of the catalysing effect of climate change on evolution) that would otherwise require geological times to be empirically validated in more stable environments. For this purpose, metrics measuring evolvability (i.e., capacity of a system to produce evolutionary adaptations - Colegrave & Collins, 2008;Kirschner & Gerhart, 1998) might be used to assess the likelihood of a clade to develop phenotypic ETs in changing environments. An investigation of morphological integration in the cranium of the mammalian order Carnivora revealed, for instance, that canids possess a higher amount of phenotypic evolvability if compared to other carnivorans. This outcome is mainly produced by an elevated incidence of pure repetitive sequences, promoting new genetic variants by duplication (Laidlaw et al., 2007;Machado et al., 2018). The conspicuous evolvability of this family was suggested to have underpinned several episodes of directionality in cranial evolution of canids associated with dietary shifts (e.g., increase in relative brachyrostry in hypercarnivore species- Machado et al., 2018). Recent macroevolutionary analyses performed on avian morphology suggest the presence of strong form-function associations that channel phenotypic variation enhancing evolvability towards specific directions of the trait space and produce repeated patterns of morphological convergence that result in highly predictable phenotypes (Felice et al., 2018;Pigot et al., 2020). This evidence seems to confirm the potential that studies on phenotypic ETs in climate change contexts might have for understanding the dynamics underpinning the predictability of evolution. Such studies might open the way for innovative approaches that might improve several inferences about patterns and processes in macroevolution, even producing indirect repercussions in climate change and conservation biology.

Concluding Remarks and Future Directions
ETs channel trait variations and often occur in the presence of directional shifts of selective conditions. Relying on the orientation of evolutionary trajectories in trait spaces, ETs involving multiple clades can be seen as a continuum that goes from convergent through parallel to divergent evolution, making multiscale analyses necessary to distinguish these classes.
Since phylogenetic history, evo-devo constraints, environment, and biological function all play a fundamental role in the occurrence of phenotypic ETs, adopting a morphodynamics-based approach to research on ETs in phenotypic evolution represent a pivotal step forward to fully understand the dynamics that produce this type of evolutionary patterns. Therefore, using methods to assess the influence of each of these factors taken individually, or even performing multifactorial analyses, is of vital importance to shed light on this field (Fig. 5). The possibility to incorporate spatially structured variations, together with the inclusion of fossil record in already existing modelling of evolution, is expected to produce the next major breakthroughs, in terms of analyses on phenotypic ETs, that are likely to be designed in the foreseeable future. An intriguing new frontier of macroevolution is represented by the possibility to perform research on ETs in the presence of extreme environmental shifts, like those resulting from the current phase of climate change. Such research might be crucial for refining predictions of future or hypothetical evolutionary outcomes and might lead to indirect implications in climate change and conservation biology.

Morphodynamics
Conceptual framework used to describe all processes influencing a morphology during its evolutionary and developmental histories, including also the impact of external environment.
Evolutionary developmental biology (Evo-devo) Scientific discipline investigating the mechanistic interactions occurring between phenotypic variation and individual development during evolution. Central questions in evo-devo researches often concern key concepts like integration, modularity, plasticity, and biological innovations.

Fig. 5
Operational framework to investigate ETs using morphodynamics. Any trend observed plotting phenotypic traits in a trait space can be investigated in the light of four factors (i.e., phylogeny, evodevo constraints, environment, and biological function). The impact of each factor is mapped back into the trait space to quantify the amount of selective advantage variation explained over the phenotypic distribution. Optimisation concepts such as Pareto front can be employed to account for multiple constraints and functions within the same surface. Pareto front was applied to phenotypic evolution mainly to investigate the limits of biological function. Typical examples can be found in Shoval et al. (2012), where the authors demonstrated that phenotypic distribution over a morphological trait space generally falls within a triangle whose vertices represent the archetype (i.e., optimal morphological solution for a particular task). Sho-  (2020) demonstrated that functional optimization can also be interpreted applying PCMs, thus covering both biological function and phylogeny within the same framework. The visualisation of horse FEA (light blue box) was modified from Panagiotopoulou et al. (2016), both the spatial distribution plot and the phylogeny of extant equids (pink and yellow boxes, respectively) were modified from Cucchi et al. (2017), and the visualisation of Hox genes in Drosophila, mice, and humans (green box) was modified from Pang and Thompson (2011). Trait space plot (white box) was produced following http:// qingk aikong. blogs pot. com/ 2018/ 05/ (Color figure online)

Modularity
Condition where a biological structure is subdivided into units (i.e., modules) that possess a strong internal covariation but weakly interact with the other units. In morphology, a classic example of modular structure is the mammalian skull, traditionally subdivided into facial and neurocranial modules. Integration Tendency of multiple traits to covary throughout a biological structure. Cryptic genetic variation Alternative developmental variants uncovering genetic variations that have negligible effects on phenotypic variations except under atypical conditions. Genetic accommodation Evolutionary process that refines a phenotype towards an adaptive condition through selection driving quantitative genetic changes. Ecological opportunity Condition where a lineage or a clade can suddenly occupy a new niche space with a wide range of underexploited resources. The sudden disappearance of competitors, the access to new environments or the appearance of key innovations are some of the most common triggers producing ecological opportunity. Geometric Morphometrics (GMM) Technique of morphological quantification able to extract (and separate) size and shape data from a biological structure, relying on the Cartesian coordinates of corresponding anatomical points, known as landmarks, placed on every specimen of the analysed sample. Finite Element Analysis (FEA) Engineering approach frequently employed to quantify the morphological function of a biological structure. The results, expressed in terms of stress and strain, allow estimating the effects produced by applying different loads to the structure and the resulting deformation. Deterministic chaos Dependency of the dynamics of a system on its initial conditions. Phenology Timing of life cycle events of a species, such as flowering, fruiting, and seasonal migrations.
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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence 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 licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.