Introduction
Among all development and phylogeny-related topics, heterochrony has been described consistently as an evolutionary process capable of bringing about macroevolutionary change. Both the term and the concept were coined by Haeckel (Reference Haeckel1866) when he proposed his “biogenetic law,” also called the theory of recapitulation. The initial definition was under the category of cenogenesis, as displacements of developmental time or dislocations of the phylogenetic order of succession, and has oftentimes been associated with the condensation of ancestral adult stages in descendants (Haeckel Reference Haeckel1866; Gould Reference Gould1977).
Heterochrony is here accepted as “change in timing or rate of developmental events, relative to the same events in the ancestor” (McKinney and McNamara Reference McKinney and McNamara1991: p. 387). Different models have been suggested to describe the processes it includes. Besides eight morphogenetic modes raised by De Beer (Reference De Beer1930) and the clock model of Gould (Reference Gould1977), Alberch et al. (Reference Alberch, Gould, Oster and Wake1979) placed heterochrony within a strictly quantitative formalism and introduced an ontogenetic trajectory through a multidimensional space to record shape and size development along with age. Through the 1980s to the present, heterochrony studies have expanded to include all temporal-related morphological changes, including whole body, organs, and morphological traits (McNamara Reference McNamara1986; McKinney and McNamara Reference McKinney and McNamara1991; O'Keefe et al. Reference O'Keefe, Rieppel and Sander1999; Cartolano et al. Reference Cartolano, Pieper, Lempe, Tattersall, Huijser, Tresch, Darrah, Hay and Tsiantis2015; Foth et al. Reference Foth, Hedrick and Ezcurra2016; Dial et al. Reference Dial, Reznick and Brainerd2017; Godoy et al. Reference Godoy, Ferreira, Montefeltro, Vila Nova, Butler and Langer2018), and heterochronic processes have been often simplified and presented in a two-dimensional (2D) shape–size plane.
Time and phylogeny have been two basic components of heterochrony concepts since their original formulation (Haeckel Reference Haeckel1866, Reference Haeckel and McCabe1905; Cope Reference Cope1896; Gould Reference Gould1977; Fink Reference Fink1982; Klingenberg Reference Klingenberg1998). But for fossils, developmental times are often difficult to measure. Thus, body size has frequently been used as a developmental age proxy (e.g., Wei Reference Wei1994; Zollikofer and Leon Reference Zollikofer, De León, Harvati and Harrison2006; Foth et al. Reference Foth, Hedrick and Ezcurra2016; Dial et al. Reference Dial, Reznick and Brainerd2017), despite its inaccuracy as an age indicator (Buddemeier et al. Reference Buddemeier, Maragos and Knutson1974; Shea Reference Shea1983; Emerson Reference Emerson1986; Jones Reference Jones and McKinney1988; Oschmann Reference Oschmann2009). If growth curves are plotted in a 2D shape–size space, distinctions between size and age are necessarily confounded, and misinterpretations can result, as significantly different trajectories in shape–size–age space could appear similar (McKinney Reference McKinney and McKinney1988).
The analysis of heterochronic patterns in a 2D size–shape space is also incompatible with the concept of complex heterochrony. Different from the single heterochronic process that resulted in morphological change, complex heterochrony refers to the combination of heterochronic processes affecting the same parameter, that is, shape or trait (Neige et al. Reference Neige, Marchand and Laurin1997), and has long been referred to as a common observation (Alberch et al. Reference Alberch, Gould, Oster and Wake1979; Dommergues et al. Reference Dommergues, David and Marchand1986; Klingenberg Reference Klingenberg1998). Combinations of two (or more) heterochronic processes have been found to characterize a number of evolutionary changes (Dommergues et al. Reference Dommergues, David and Marchand1986; Klingenberg and Spence Reference Klingenberg and Spence1993; Neige et al. Reference Neige, Marchand and Laurin1997; Klingenberg Reference Klingenberg1998; Godoy et al. Reference Godoy, Ferreira, Montefeltro, Vila Nova, Butler and Langer2018). Unfortunately, when two or more processes are present, the geometry of their effect(s) can be difficult to resolve in this simple space, because size–shape correlated growth patterns may display little apparent conformance to any single process (Godoy et al. Reference Godoy, Ferreira, Montefeltro, Vila Nova, Butler and Langer2018) or may exhibit the mimetic process in others (Neige et al. Reference Neige, Marchand and Laurin1997). Regardless, size and shape are two independent morphological indicators and might suggest the operation of different relative growth processes if they are studied separately in the context of an appropriate age indicator.
In the case of organisms whose shells or bones preserve the aspects of size, shape, and developmental stage, combined heterochronic processes can be explored in their shape–size–age relations thoroughly. The growth rings of mollusks and cephalopods, chamber septa number of ammonites, and chamber and whorl numbers of foraminifera all have been used to estimate organism ages (Doguzhaeva Reference Doguzhaeva1982; Dommergues Reference Dommergues1988; Bucher and Guex Reference Bucher and Guex1990; Neige et al. Reference Neige, Marchand and Laurin1997; Shi and MacLeod Reference Shi and MacLeod2016). Indeed, the chambers or whorls of many fossil taxa contain many valid developmental time indicators such that their shape and size ontogenies can be studied separately to investigate complex patterns of heterochrony.
Fusulinids, the oldest larger benthic foraminifera, lived in shallow-marine late Paleozoic environments (Ross Reference Ross, Scholle, Peryt and Ulmer-Scholle1995; Hohenegger Reference Hohenegger2011) and built multichambered shells, often called tests, with planispirally coiled chambers. Test growth in foraminifera is believed to reflect the growth of the internal cytoplasmic mass (Hohenegger and Briguglio Reference Hohenegger, Briguglio, Kitazato and Bernhard2014), and these chambered tests record the size and shape development at discrete states of ontogeny. The developmental histories of fusulinid tests can be revealed in thin section. Axial thin sections (Fig. 1A,B), produced by grinding the test parallel to the coiling axis until the center of the initial chamber or proloculus is exposed, present the interior test morphology as well as the test's cross-sectional form of each whorl. As a result, the forms of all fusulinid tests can be observed, measured, and studied quantitatively at each whorl-based stage of their ontogenetic development.
Figure 1. The whorl exterior test and axial thin sections of fusulinids. A, Fusulinid individual sample displaying a grind direction of axial section; B, an axial section of sample A; C–M, species used in the analysis: C, Sphaeroschwagerina constans; D, Pseudoschwagerina neotruncata; E, Triticites creekensis; F, Triticites concaviclivis; G, Triticites noinskyi; H, Zellia colaniae; I, Robustoschwagerina yishanensis; J, Pseudoschwagerina subconvexa; K, Triticites stuckenbergi; L, Triticites cf. T. vicrotioensis; M., Triticites cellamagnus. Scale bars, 1 mm.
In the evolutionary history of fusulinids, there is one significant taxonomic radiation event in the early Permian associated with the macroevolutionary occurrence of a new group with large inflated tests, namely Pseudoschwagerininae Zhang, Reference Zhang1963 (Leven Reference Leven1993; Shi and Yang Reference Shi and Yang2005; Yang et al. Reference Yang, Liu, Zhu and Shi2005), and there has long been a suspicion among fusulinid taxonomists that a heterochronic process might have been responsible for this morphological novelty. Pseudoschwagerine and an earlier-evolved genus belonging to the same Schwagerinidae family, namely Triticites, both developed similar fusiform juvenile tests with identical keriotheca-containing test walls and distinct chomata (black parallel ridge pairs deposited to prevent cytoplasmic scattering); therefore, the latter is considered the presumptive pseudoschwagerine ancestor (e.g., Ross Reference Ross1967; Rozovskaya Reference Rozovskaya1975; Sheng et al. Reference Sheng, Zhang and Wang1988; Yang and Hao Reference Yang and Hao1991; Yang Reference Yang1992; Leven Reference Leven1993). While all triticitid species retain a fusiform test shape in their adult life-history stage, pseudoschwagerine species adopted a loosely coiled and hugely inflated fusiform, subspherical or spherical test in the adult stage. The significant shape and size differences exhibited by triticitid and pseudoschwagerine species document a substantial level of morphological divergence in the pseudoschwagerine lineage.
Of the several pseudoschwagerine genera, Robustoschwagerina has been regarded as having evolved through hypermorphosis (Yang and Hao Reference Yang and Hao1991; Yang Reference Yang1992). While the morphological analyses in these studies were incomplete, they raise the questions of whether the pseudoschwagerine lineage resulted from heterochronic modification of ancestral developmental conditions and which heterochronic pattern, or combination of patterns, was responsible for this macroevolutionary event. The present study aims to test several related hypotheses, including: (1) whether geometric morphometric routine could be employed to expose heterochronic pattern during developmental history of fusulinid foraminifera; (2) whether a single or complex heterochronic pattern reflects the mechanism responsible for the pseudoschwagerine evolution in the early Permian; and (3) whether the commonly used shape and size correlation would bias a heterochronic analysis.
Materials
The Pseudoschwagerininae appeared abruptly in the early Permian fusulinid radiation. Among the four genera mainly developed in this radiation, Pseudoschwagerina, Robustoschwagerina, Sphaeroschwagerina, and Zellia, Pseudoschwagerina exhibits inflated fusiform to subspherical adult test forms (Fig. 1D,J), while the other three exhibit subspherical to spherical test forms (Fig. 1C,H,I). Another schwagerinid genus, Triticites, has long been considered the likely genus into which the pseudoschwagerine ancestor would be placed (e.g., Ross Reference Ross1967; Rozovskaya Reference Rozovskaya1975; Yang and Hao Reference Yang and Hao1991; Yang Reference Yang1992; Leven Reference Leven1993).
The present study included 66 specimens belonging to pseudoschwagerine genera (15 Pseudoschwagerina, 24 Robustoschwagerina, 14 Sphaeroschwagerina, and 13 Zellia; see Table 1). Forty Triticites specimens, including 21 from the late Carboniferous as the ancestor group and 19 from the early Permian as descendants, were also recruited. Among all the specimens, 100 were collected or reported from China and 6 were from other countries (Table 1). The Carboniferous triticitid forms are morphologically close and therefore most likely have a phylogenetic connection to the pseudoschwagerine species of this study. Examples include Triticites concaviclivis (Fig. 1F) as the possible ancestor of Sphaeroschwagerina constans (Fig. 1C), both having tightly coiled juvenile whorls with medium-sized proloculi, slim spirotheca, and small but distinct chomata; Triticites creekensis (Fig. 1E) as the possible ancestor of Pseudoschwagerina neotruncata (Fig. 1D) for their loosely coiled juvenile whorls, medium-sized proloculi, and narrow and distinct chomata; and Triticites noinskyi (Fig. 1G) as the possible ancestor of the pseudoschwagerine species exhibiting loosely coiled juvenile whorls, large proloculi, and strong and wide chomata, such as Robustoschwagerina yishanensis (Fig. 1I), Zellia colaniae (Fig. 1H), and Pseudoschwagerina subconvexa (Fig. 1J).
Table 1. Species in the analysis
The Permian triticitid species were selected for analysis to represent the developmental history of descendants with fusiform tests. Most Permian triticitid species are distinct from Carboniferous species for developing strongly folded septa and large fusiform tests with more whorls, but the identical inner whorls between species disclose the phylogenetic relationships. Examples include Permian Triticites stuckenbergi (Fig. 1K) as the possible descendant of Carboniferous Triticites cf. T. vicrotioensis (Fig. 1L) as both developed small proloculus, slim but distinct chomata, and slightly folded septa in the inner whorls; and T. noinskyi as the possible ancestor of Permian Triticites cellamagnus (Fig. 1M) as both developed large proloculi, distinct chomata, and straight septa in the inner whorls.
It is uncertain at present which species are the direct ancestors of the studied pseudoschwagerine and Permian triticitid species, and the morphologies of these taxa might well vary in some small respects from the presumed ancestor forms. Chances would be rare for some South China triticitids as ancestors of the American Permian triticitids. The strategy of mixing multiple species of each genus was therefore conducted to combine patterns of morphological variation into averaged morphologies, to achieve unbiased representation of the ancestor and descendant groups. Altogether, six groups—Carboniferous Triticites, Permian Triticites, Pseudoschwagerina, Robustoschwagerina, Sphaeroschwagerina, and Zellia, were the main focus of comparison and discussion through the entire analysis.
Methods
Thin sections of South China specimens (99 in total) were photographed, and seven previously published specimen images were rescanned to obtain high-resolution representations of the axial-sectional morphology. Ontogenetic test-form growth patterns were quantified using the tpsDig program via representations of changes in the axial-sectional outline of successive growth stages (see Shi and MacLeod Reference Shi and MacLeod2016). To facilitate comparison, the growth interval was set to coincide with chamber whorl formation, irrespective of the actual time interval between successive chamber secretions. The final chamber of each half whorl was exposed in the axial section and therefore numbered to represent the growth interval as a whorl stage. Each axial section image was edited to expose a series of outlines illustrating the overall form of the test at each whorl stage (Fig. 2). Test outline at each whorl stage was digitized by specification of four type 2 landmarks, located at the maxima of curvatures (Bookstein Reference Bookstein1986), and 36 semilandmarks representing the interlandmark curve segments (Fig. 2). The number of semilandmarks required to represent these test outlines was determined by shape complexity analysis, with the outline fidelity index set to ≥95% (MacLeod Reference MacLeod1999; Shi and MacLeod Reference Shi and MacLeod2016). In total, 1396 test outlines representing ontogenetic form changes in the 106 fusulinid specimens were acquired (Supplementary Material).
Figure 2. Digitization protocol for a specimen image. Outline digitized by four type 2 landmarks (black dots) and 36 other equally spaced semilandmarks (gray dots) representing the four divided segments. A, Test outlines at the successive 13th and 14th whorl stages. B,C, Outline of the odd-numbered stages (i.e., 13th) were inverted to ensure the ultimately added chamber was in an orientation consistent with those of the even-numbered stages (i.e., 14th). C,D, Standardized outlines after centroid size value has been removed.
To study the size and shape variation aspects of fusulinid test form, generalized least-squares Procrustes superposition (Rohlf and Slice Reference Rohlf and Slice1990) was employed to align all the outlines and normalize them for the extraneous factors of position, orientation, and size. Test size variation was quantified using the centroid size index (square root of the sum of squared distances of the landmarks in configuration to their mean location; Bookstein Reference Bookstein1986; MacLeod Reference MacLeod2008) and extracted from the outlines before they were aligned. Following outline alignment, principal component analysis (PCA) based on the shape covariance structure of the pooled outline dataset was used to facilitate geometric interpretation (MacLeod Reference MacLeod2009; Shi and MacLeod Reference Shi and MacLeod2016). When the projected scores of the outline data along the principal component axes captured various aspects of their shape variation, these scores could be used to represent the outline shapes. As a result, the specimens’ ontogenetic shape trajectory could be constructed with suitable principal component scores representing the majority shape variation of the outline dataset at comparable developmental stages. Test size estimates were also added to each trajectory as a third, independent variable. The shape and size trajectories for each studied group were composed with the means of shape scores or size values of the group specimens.
Whorl stage based on half-whorl counting was regarded as the developmental age indicator in the current plots of ontogenetic variation. Chamber number is believed to be an accurate age proxy for some benthic foraminifera (Hemleben et al. Reference Hemleben, Spindler and Erson1989; Speijer et al. Reference Speijer, Van Loo, Masschaele, Vlassenbroeck, Cnudde and Jacobs2008; Hohenegger and Briguglio Reference Hohenegger, Briguglio, Kitazato and Bernhard2014), but whorl number is regarded a valid alternative for fusulinids and could be traced in the axial sections (Yang and Hao Reference Yang and Hao1991; Yang Reference Yang1992; Shi and MacLeod Reference Shi and MacLeod2016). Although the true developmental time of individual whorls cannot be assumed to be equal, for the ancestor–descendant comparisons (e.g., Pseudoschwagerininae and Triticites), patterns of their whorl development and the indicated time or age interval would be similar.
Ontogenetic trajectories of the studied generic group were summarized by patterns of variation in shape (PC 1 scores), size (centroid size), and age (whorl stage), with their mean values, in a three-dimensional (3D) Cartesian space. Trajectories were projected in the 2D shape–size, shape–whorl stage, and size–whorl stage spaces to compare episodes of morphological development. Ontogenetic trajectories of all 106 specimens and six groups were plotted in the shape–whorl stage and size–whorl stage planes to isolate those aspects of the morphological development processes, including potential heterochronic pattern, involved between triticitids and pseudoschwagerines.
Results
A Fusulinid Shape Space
A pooled shape space for the studied specimens was established through a PCA of the Procrustes-aligned outline coordinates for all specimens at each whorl stage. In total, 80 shape dimensions account for 100% of the observed shape variance, and the first principal component summarized the overwhelming majority (75.6%) of shape structure (Fig. 3A). The first two dimensions of the PC shape space, as well as the shape models along PC 1, are shown in Figure 3.
Figure 3. Patterns of fusulinid test shape variation represented by the first two principal components. A, Principal component analysis based on the pooled sample shape covariance matrix for all developmental whorl stages of all specimens belonging to Triticites and the four genera of Pseudoschwagerinae. Test shapes of different whorl stages have been marked with circles with different radii. The positions of shape models in B are marked with black dots. B, Shape models. The models were calculated with coordinates along the PC 1 axis in the principal component shape spaces of A and are displayed with numbers representing their PC 1 coordinates.
Shape models reconstructed from the PC 1 scores identified the major shape variation among all the test outlines (Fig. 3B). Outline shapes that project to positions near the negative end along the PC 1 axis represent highly elongate fusiform shapes (e.g., close to the characteristic adult test shapes of triticitids), whereas those projecting to positions near the positive end demonstrate the inflated fusiform to spherical shapes similar to adult pseudoschwagerine tests. While juvenile tests (at 1–2 whorl stages) of most specimens have high positive PC 1 scores, the tests that developed 2–4 whorls (at 3–9 whorl stages) mostly projected in the middle area of the PC 1 axis, and mature tests of triticitids and pseudoschwagerines finally diverged to the two ends. Therefore, the PC 1 score of the outline shape can be used to quantify the primary shape change during developmental history.
Ontogenetic Trajectories and Their Projections
Ontogenetic trajectories in the 3D shape–size–whorl stage space demonstrate the nonlinear growth pattern of the studied generic groups (Fig. 4A). Roughly two bundles of trajectories could be recognized in this space, one is composed of four pseudoschwagerine genera and another of the triticitid groups. The four pseudoschwagerine genera exhibit similar patterns, and three subintervals describing significant morphological change could be recognized in each of them. The triticitid groups present patterns formed by only two distinct subintervals, similar to the first two of the pseudoschwagerines.
Figure 4. Ontogenetic trajectories of four pseudoschwagerine genera and two triticitid groups in the 3D shape–size–whorl stage space (A) and in 2D shape–size, shape–whorl stage, and size–whorl stage planes (B). Shape value is represented by PC 1 scores (no unit) and size value is represented by centroid size index (no unit). Shape and size value means of the studied species in each group are used. Error bars are not displayed to allow trajectories in 3D space and 2D planes to be distinguished.
The first subinterval (S1) encompasses the developmental period when all six trajectories expanded along the direction of test size, whorl stage counting increasing and shape value decreasing (Fig. 4A). It covers the first several whorl stages for pseudoschwagerines and most stages through triticitid development. S1 is followed by the second subinterval (S2), when whorl stage, test size, and shape value increased concurrently, and the curves therefore turned into a different direction. For pseudoschwagerines, S2 subintervals consist of 8 to 10 whorl stages while those of triticitids only contain two stages. The last two or three stages of pseudoschwagerines constitute the third subinterval (S3) in their development, during which whorl stage number kept increasing, test sizes slowly increased, and shape values became stable or even dropped slightly. Because of these changes, the curves’ extension direction turned again and formed the distinct “tails” of the pseudoschwagerine trajectories.
These subintervals could be easily recognized in the trajectory projections in the shape–size plane (Fig. 4B). In the beginning of their development history, the six groups clustered together in the S1 subinterval, and all experienced a significant shape value drop from ca. 0.2 to ca. −0.1 within a minor size rise. After that, the divergence became evident when pseudoschwagerine shape score started to increase and their S2 subinterval started. Triticitid groups stayed in S1, and their shape value kept dropping. In S2, the sphaeroschwagerinid, robustoschwagerinid, and zellid curves clustered as their tests gradually became spherical, while test size increased significantly. The pseudoschwagerinid shape value only increased to ca. 0.05, much lower than those of the other three, describing the inflated fusiform shape of their mature tests. S3 subintervals of pseudoschwagerines are composed of the zigzag “tails” that are actually false signals resulting from variation of the mean size value.
The trajectories projections in the shape–whorl stage and size–whorl stage planes separated the shape and size variations and exhibited more regular patterns (Fig. 4C). The three subintervals could also be recognized in the shape–whorl stage but are less significant in the size–whorl stage plane. The trajectories of all species revealed more details during their developmental histories.
Ontogenetic Trajectories in the Shape–Whorl Stage Plane
Because whorl stages were common to all fusulinids, ontogenetic trajectories can be constructed such that the test shape–phase intervals can be recognized for all the studied specimens (Fig. 5). As mentioned earlier, the ontogenetic trajectories of pseudoschwagerine species contain three recognizable form change subintervals, while triticitid species only include two and exhibit a decreasing shape value trend over most of their ontogeny.
Figure 5. Ontogenetic shape and size trajectories of pseudoschwagerines and triticitids. A–F, Species-based shape trajectories in shape–whorl stage space. Trajectories of all studied species belonging to each genus are plotted in colored lines with pseudoschwagerine generic and triticitid group trajectories represented by mean values (black lines) of the species with 95% confidence level (black whiskers). Species list refers to Table 1. G, Assembled shape-average trajectories for four generic and two triticitid groups in the shape–whorl stage space illustrated by colored lines and whiskers. H, Size-average trajectories for pseudoschwagerine generic and two triticitid groups represented by mean values (colored lines) with 95% confidence level (colored whiskers) in size–whorl stage plane.
The S1 intervals of pseudoschwagerines display a decreasing trend, with shape values dropping from 0.2 to −0.1. The proloculus is spherical in almost all foraminiferal species (Alve and Goldstein Reference Alve and Goldstein2003; Murray Reference Murray2006), and in this analysis, its shape is represented by a shape value of ca. 0.2. From this point, the developmental curves of all pseudoschwagerine species shape values decrease for four to eight whorl stages, indicating these tests altered to a fusiform shape as new chambers were added. In triticitids, S1 shape values generally continue decreasing beyond the states characteristic of pseudoschwagerines, which indicates their test shapes became more fusiform.
S2 and S3 constitute two phases of an overall inflation trend. At about the fourth (for robustoschwagerinids and zellids) or eighth (for pseudoschwagerinids and sphaeroschwagerinids) whorl stage, shape values of pseudoschwagerine tests reached their lowest, indicating the most fusiform shape the species have developed. Then pseudoschwagerines began to exhibit test inflation. During S2, their shape scores turned back toward 0.1 at about the 11th or 14th whorl stage, reflecting tests growing into spherical forms. In triticitid species, S2 does not occur until the very end of their development stages, that is, at the 12th or 13th whorl stage, where the shape score increased. This suggests that their typical fusiform shape was maintained throughout the adult period, followed by an incipient interval of test inflation. S3 appears at the very end of the developmental sequence for all pseudoschwagerine species and is marked by a visible turn in the developmental trajectory over the last one or two whorl stages. During this interval, the shape score of most species remained steady or exhibited a slight drop in value, suggesting test inflation ceased when individuals entered the final (possibly the mature) stage of development.
The average shape trajectories were constructed by shape value means of the species in each group along with developmental stages. Six curves exhibit divergences in the later developmental intervals, and three assemblages can be recognized (Fig. 5G). The triticitid assemblage consists of the two triticitid groups, with a long S1 interval including about 13 whorl stages and a short S2 duration consisting of the last one or two whorl stages. The “middle” assemblage includes pseudoschwagerinid and sphaeroschwagerinid species, both with long S1 and S2 intervals that occupy nearly half of their life histories. The third assemblage consists of robustoschwagerinid and zellid species, both of which exhibit a relatively short S1 subinterval with a longer S2.
The Unique Inflation Growth of Pseudoschwagerininae Suggested by Size Change
Size trajectories for the four pseudoschwagerine genera, as well as triticitids, all displayed exponential growth trends (Fig. 5H). Two groups are recognized based on their test size trajectories, one composed of robustoschwagerinid, pseudoschwagerinid, and zellid species (= the RPZ group) and the other composed of sphaeroschwagerinid and triticitid morphologies (= the ST group). During early ontogeny, all species increased in size slowly relative to whorl stages, and the RPZ group exhibited a markedly higher growth rate. At the fourth or fifth whorl stage, the rate of whorl-marked growth in this group increased rapidly, with zellids outpacing the other two genera. During this interval, the ST group was characterized by a notably slower whorl growth rate. As a result, the size trajectories of these groups diverged substantially.
During the 4th to 10th whorl stage, growth rates in the RPZ consistently remained several times higher than in the ST group, and test size of the zellids turned out to be the largest at their mature period. Sphaeroschwagerinid whorl growth rates are as slow as those of triticitids in their early ontogeny, until around the 10th whorl stage, but increased in apparent lockstep with test inflation after that point. In late ontogeny, that is, the 14th whorl stage, sphaeroschwagerinid test size and shape change rates both became very similar to those of both the robustoschwagerinids and zellids, while triticitids retained a slow growth.
Discussion
In the evolution history of fusulinids, the early Permian pseudoschwagerines comprise a special branch whose taxa exhibited significant shape shifts during their ontogeny (Ross Reference Ross1967; Rozovskaya Reference Rozovskaya1975; Leven Reference Leven1993). These spherical–fusiform–spherical test shape shifts have been quantitatively described through a geometric morphometric routine so that further comparison regarding the shift timing and process between different taxa or groups could be conducted. The advantage of the PC 1 score representing the majority of fusulinid test outline shape variation might not always be the case, but through PCA the major shape variation can be captured and investigated with sensitive procedures. Ontogenetic trajectories of fusulinid foraminifera established here provided sufficient information about their shape and size changes across their life spans and are therefore well suited to various research aims, including heterochronic analysis.
Growth curves reflecting the shape–size correspondence have been frequently used in morphological studies, including heterochronic analyses (Mitteroecker et al. Reference Mitteroecker, Gunz and Bookstein2005; Zollikofer and Leon Reference Zollikofer, De León, Harvati and Harrison2006; Gerber et al. Reference Gerber, Neige and Eble2007; Foth et al. Reference Foth, Hedrick and Ezcurra2016; Dial et al. Reference Dial, Reznick and Brainerd2017). They have been established here (Fig. 4B), and the patterns are close to the ontogenetic trajectories in the 3D shape–size–whorl stage space. However, the shape–size plane failed to provide either an accurate or an adequate record of trends in size variation. In the plane, the growth curve of zellids overlapped with those of robustoschwagerinids and sphaeroschwagerinids when the size value of their tests was from two to nine, but their trajectories clearly separated in the 3D space, due to the different size growth rates (Fig. 4). To avoid this type of bias, an independent age indicator is a necessity. Creatures containing an independent, countable time record, such as chambered shells or growth lines, represent the only candidate fossil groups suitable for the detailed study of heterochrony.
The well-established theoretical relationships among size, shape, and age in the context of heterochrony analyses (Gould Reference Gould1977; Alberch et al. Reference Alberch, Gould, Oster and Wake1979; Raff Reference Raff1996), suggest that the 3D developmental space plays an important role in the analysis and comparison of heterochronic patterns, as well as in the testing of various pattern-level heterochronic hypotheses. Retention and use of all three of these critical heterochronic variables are necessary to achieve a complete understanding of any organism's developmental pathway, irrespective of taxonomic/phylogenetic level. When shape development of pseudoschwagerines, along with size increases, was quantified and compared with those of their presumptive triticitid ancestors during their life-history stages, the changes documented suggested the operation of heterochronic processes. In particular, information on developmental time in the ontogenetic trajectories could further determine whether condensation (acceleration) or retardation is displayed in pseudoschwagerine evolutionary history.
Peramorphic Heterochronic Growth of Pseudoschwagerines
The consistent shape changes of Carboniferous triticitids and pseudoschwagerines during their early developmental stages could be observed in the current results in the form of the overlapped shape trajectories during the S1 subinterval (Fig. 5G). The divergence in their post-S1 developmental shape change can now be seen as a pseudoschwagerine apomorphy and results in the morphological distinctions between the two lineages, triticitid-represented fusiform schwagerines and inflated pseudoschwagerines.
Carboniferous triticitids developed a rather long S1 and a very short S2 (Fig. 5A). Pseudoschwagerine species repeated these two intervals, but S1 became relatively shortened and S2 extended substantially (Fig. 5B,D,F). Their shape trajectories overlapped with those of Carboniferous triticitids during the S1 interval, and S2 started at the “middle” of the Carboniferous triticitids S1, though with variation in timing. The shape trajectory of sphaeroschwagerinids diverted a bit from the others during the middle of its S1 interval but was soon back on track. S2 intervals of pseudoschwagerines developed earlier than those of their triticitid ancestors, fitting the criterion of peramorphic predisplacement heterochrony. Their S2 intervals are much longer than the two whorl stages of Carboniferous triticitids, suggesting retardation from peramorphic hypermorphosis. Moreover, the shape changes of pseudoschwagerines during the S2 interval are faster than those of their triticitid ancestors, as indicated by the higher slopes of their S2 curves. This phenomenon indicates the peramorphic acceleration of shape change. The recapitulation of S1 and S2 ontogenetic stages of the (presumed) ancestors by descendant pseudoschwagerines, along with the developmental truncation indicated by the shortened S1, retardation indicated by the prolonged S2, and acceleration indicated by the S2 higher rate of shape change, can be attributed to the extreme shape novelties displayed by the pseudoschwagerines. In particular, the shortened S1 developmental time interval indicates peramorphic predisplacement, the developmentally prolonged S2 suggests peramorphic hypermorphosis, and higher shape change rate reveals peramorphic acceleration.
The maturity of foraminifera, both benthic and planktonic, is oftentimes indicated by the slow growth of last two whorls or chambers (Hemleben et al. Reference Hemleben, Spindler and Erson1989; Speijer et al. Reference Speijer, Van Loo, Masschaele, Vlassenbroeck, Cnudde and Jacobs2008; Hohenegger Reference Hohenegger2011). In fusulinids, growth termination is coincident with the S3 subinterval. Along with the shape value plateau or drop, size growth rate of pseudoschwagerines slowed down in their end stages, that is, the 12th to 14th whorl stages (Figs. 4B, 5H). In the size–time space, the growth rates of triticitid species apparently slowed down after about the 11th whorl stage. By analogy, the decreasing or plateauing growth rate trend should be regarded as a signal of fusulinid reproductive maturity or onset of the reproduction phase of their life history. Therefore, compared with the triticitid ancestor, maturity in pseudoschwagerine genera occurred in a later stage, indicating a developmental change consistent with hypermorphosis.
Predisplacement shortened the early developmental periods of the characteristic triticitid fusiform shapes, while acceleration, expressed as an increased rate of shape change (Gould Reference Gould1977; Alberch et al. Reference Alberch, Gould, Oster and Wake1979; Raff Reference Raff1996), sped up the inflated shape change trend. Hypermorphosis is expressed by prolongation of ontogeny, which causes delayed maturation (De Beer Reference De Beer1930; Gould Reference Gould1977; Raff Reference Raff1996). This developmental change, along with the predisplacement and acceleration, led to the expression of a radically new morphology in these triticitid descendants, such as larger body size, larger test size, and new shapes (Gould Reference Gould1977; McNamara Reference McNamara1986).
Compared with their Carboniferous counterparts, Permian triticitids developed a similar S1 and an S2 with identical developmental time but higher shape change rate (Fig. 5A,C,H). This slight acceleration resulted in visible changes to these descendants’ test shapes.
Size Rate Modification of Pseudoschwagerines
Size change was not defined as a heterochronic process (Gould Reference Gould1977; Alberch et al. Reference Alberch, Gould, Oster and Wake1979), but gigantism and dwarfism have long been regarded as resulting from the action of heritable changes in developmental timing (Gould Reference Gould1977). Later, peramorphic acceleration was explained as increased rate of morphological development that could affect the whole organism, organs, or structures. Thus, overall size change has long been thought to be an important, though largely unacknowledged, element of heterochrony (McNamara Reference McNamara1986; Raff Reference Raff1996; Webster and Zelditch Reference Webster and Zelditch2005) and often represents an important factor responsible for complex morphological changes (Neige et al. Reference Neige, Marchand and Laurin1997; McNamara and Long Reference McNamara, Long, Brett-Surman, Thomas, Holtz, Farlow and Walters2012; Godoy et al. Reference Godoy, Ferreira, Montefeltro, Vila Nova, Butler and Langer2018).
Ontogenetic trajectories of Carboniferous triticitid and pseudoschwagerines in the size–whorl stage space displayed exponential growth patterns (Fig. 5H). This type of size trajectory is common to most (possibly all) larger benthic foraminifera (Briguglio et al. Reference Briguglio, Metscher and Hohenegger2011; Hohenegger and Briguglio Reference Hohenegger, Briguglio, Kitazato and Bernhard2014). However, the growth rates recorded by the slopes of these curves are variable. Robustoschwagerinids, sphaeroschwagerinids, and zellids exhibited higher values than their Carboniferous triticitid ancestor(s) across the whole developmental history, while sphaeroschwagerinids did not get beyond their ancestors’ size growth until the 12th whorl stage. Higher size growth suggests descendants grow faster than ancestors through certain developmental stages, a pattern that can be assigned to peramorphic developmental acceleration (McNamara Reference McNamara1986; Raff Reference Raff1996) or recast as rate modification (Webster and Zelditch Reference Webster and Zelditch2005). Based on these preliminary results, accelerated size growth may have existed from the beginning of development for robustoschwagerinids, pseudoschwagerinids, and zellids. However, in sphaeroschwagerinids, this trend appeared much later in the developmental sequence, although the rate increased quickly and soon exceeded that characteristic of triticitids (Fig. 5H). As a result, size acceleration appears responsible for the pronounced test-size enlargement characteristic of pseudoschwagerine fusulinids as a group.
Complex Heterochrony in Pseudoschwagerine Development
Based on the fusulinid shape–size–whorl stage trajectories documented in this investigation, three different peramorphic heterochronies—predisplacement, hypermorphosis, and acceleration—appear to characterize pseudoschwagerine developmental history. An overall peramorphic pattern has been proposed for the development of Robustoschwagerina but without further detail (Yang and Hao Reference Yang and Hao1991), and effects of pairwise combinations of heterochronic processes have been illustrated by Klingenberg (Reference Klingenberg1998). The current morphometric routine exposed a threefold complex heterochronic pattern involving all three single peramorphic processes being responsible for the unique morphology of the four pseudoschwagerine genera. Regarding the development of a certain trait (i.e., test shape), there are altogether eight possible combinations for a threefold complex heterochronic pattern (Fig. 6), and in fact, development rate (acceleration or neoteny) and time (progenesis or hypermorphosis) are the two factors determining the trait's morphology, while the timing (predisplacement or postdisplacement) decides its relationship with other traits. In the current example, hypermorphosis indicated by the prolonged test-inflation interval (S2) and acceleration indicated by the high growth rate of both test inflation and size can explain the large, loosely coiled spherical tests of pseudoschwagerines. Predisplacement, indicated by an early start of S2, resulted in the early termination of fusiform shape development. In the commonly used shape–size space, because another important aspect of organism developmental history, size growth, is totally neglected, the size-related heterochronic process would be either buried or difficult to recognize. Complex heterochrony might have been overlooked in the previous research on the development history of other fossil organisms, if body size was regarded as an age indicator.
Figure 6. Ontogenetic trajectories effected by the threefold complex heterochronic processes and the S2 shape trajectories in the current investigation. A, Example of ancestor and descendant ontogenetic trajectories. B, Effects of the eight modes of threefold complex heterochronic processes. The simple heterochronic process involved is represented by the number. C, Simplified S2 subinterval shape trajectories of Carboniferous triticitids (C-t) and pseudoschwagerines (p).
Pseudoschwagerine Thrived as K-Strategists
The early Permian Asselian to Sakmarian stages bore witness to the fusulinid radiation in which the Pseudoschwagerinae first appeared (Leven Reference Leven1993; Shi and Yang Reference Shi and Yang2005; Yang et al. Reference Yang, Liu, Zhu and Shi2005). Although the fusiform schwagerines remained major components of fusulinid shape diversity throughout this interval, pseudoschwagerines thrived and accounted for 14.2% of all fusulinid species occurring throughout this interval in well-studied South China stratigraphic successions (Shi et al. Reference Shi, Yang and Huang2009b). As the schwagerines that flourished in this interval remained morphologically close to those of late Carboniferous species with fusiform adult tests (e.g., Triticites), pseudoschwagerines constituted a considerable morphological departure from this standard. The successful radiation of pseudoschwagerines in the early Permian might be related to this radical morphological change.
Gould (Reference Gould1977) connected heterochronic variations in life-history strategies with r- and K-selection. Pseudoschwagerines have been regarded as representing fusulinid K-strategists (Leven Reference Leven1993; Shi and Yang Reference Shi and Yang2005; Yang et al. Reference Yang, Liu, Zhu and Shi2005), a designation they share with numerous other coeval species. The Asselian to Sakmarian interval is believed to have been a favorable time for fusulinids, with its high dissolved oxygen levels and the widespread occurrence of suitable, shallow-marine habitats (Berner Reference Berner2006; Groves and Wang Reference Groves and Wang2009; Payne et al. Reference Payne, Groves, Jost, Nguyen, Moffitt, Hill and Skotheim2012). In such a stable and crowded environment, heterochronic hypermorphosis likely resulted in the retarded maturation typical of K-strategists (Gould Reference Gould1977), within the pseudoschwagerine lineage, resulting ultimately in the appearance of forms characterized by large, spherical tests. The spherical test could also facilitate transportation to and colonization of new sites via selective transportation by turbulent bottom currents (Shi and MacLeod Reference Shi and MacLeod2016).
Fusulinids are considered to have been symbiont-bearing foraminifera because of their similarities in size and internal structure to modern symbiont-bearing benthic foraminifera (Ross Reference Ross, Scholle, Peryt and Ulmer-Scholle1995; Groves et al. Reference Groves, Pike and Westley2012). For modern symbiont-bearing foraminifera, flat and fusiform tests are usually interpreted as strategies to enlarge their test surface so more symbionts could be held within the test and exposed to sunlight (Hohenegger Reference Hohenegger2011). The spherical tests that pseudoschwagerine fusulinids developed could enhance their ability to become widely distributed, and the objective of benefiting symbionts was accomplished through large tests with large surface areas. These are probably some of the reasons for their rapid diversification, irrespective of the fact that large size would also pose biomechanical and metabolic challenges. Regardless, this diversification constituted a significant episode in the increase of diversity in the entire late Paleozoic marine invertebrate community (Fan et al. Reference Fan, Shen, Erwin, Sadler, MacLeod, Cheng, Hou, Yang, Wang, Wang, Zhang, Chen, Li, Zhang, Shi, Yuan, Chen, Zhang, Li and Zhao2020).
Conclusions
Inspection of ontogenetic trajectories of organisms in the 3D shape–size–age (time) space is required to adequately represent, assess, and analyze the morphological signatures of heterochronic processes. Organisms with shells or bones that record information pertaining to specimen age are ideal subjects for such investigation. If a complete representation of this space is combined with a morphometric approach to shape characterization, morphological analyses can be used to describe the developmental histories of organisms quantitatively and so facilitate the comparisons between taxa.
With whorl stage as an independent age indicator for fusulinids, heterochronic patterns were confirmed in the development of early Permian pseudoschwagerines. Keeping in mind that both size and shape change are necessary elements for documenting complete morphological development history, ontogenetic trajectories of pseudoschwagerine fusulinids suggest that complex patterns of heterochrony were part of and served as the mechanistic basis for these morphological transitions. Peramorphic predisplacement, hypermorphosis, and acceleration of morphological development dominated the developmental history of these groups and resulted in the novel appearance of large, inflated fusiform and spherical tests in larger benthic foraminifera. Compared with their presumptive triticitid ancestor(s), pseudoschwagerine species compressed the developmental time interval in which their test morphology was characterized by a fusiform shape and added a much longer terminal stage in which the test quickly expanded to adopt a characteristic inflated or spherical shape. Accelerated test size growth was also among the trends responsible for the appearance of large pseudoschwagerine tests.
The morphometric approach presented here can successfully summarize morphological variations among organisms such as foraminifera, thus providing a common routine for various developmental investigations. The 3D shape, size, and age (time) relationship established through this procedure can quantitatively illustrate the differences and similarities among the developmental histories of organisms and facilitate further morphological studies.
Acknowledgments
This work is supported by the National Scientific Foundation of China (grant nos. 41772017 and 91955201) and the National Key R&D Program of China (no. 2018YFE0204201). Valuable suggestions from N. MacLeod of Nanjing University, Associate Editor M. Hopkins and two anonymous reviewers are greatly appreciated. Morphometric analysis programs from Norman MacLeod have been used in the current study.
Open access
