First record of sexual size dimorphism in fossil Strombidae (Mollusca, Gastropoda) from the Miocene of Kutch, western India and its evolutionary implications

Persististrombus deperditus (Sowerby) from the Lower Miocene of Kutch, Gujarat, western India is represented by two size classes in our collection. Statistical analyses discriminate the size morphs. Large size variations generally result from either (1) sexual differences or (2) ecophenotypic causes. All the living species of the family Strombidae, wherever examined, are characterized by sexual size dimorphism (SSD). Persististrombus deperditus shares all the characters of SSD in these recent species. Size variations due to difference in ecological factors generally occur in allopatric populations. Similar variations are known to characterize sympatric sub-populations of molluscs living only in the intertidal zone, where upper and lower shorefaces differ significantly in physico-chemical and biological properties. Persististrombus deperditus comes from a stable shelf setting that received less siliciclastic input in response to transgression. Hence, its size dimorphism is considered to have sexual origin. This is the first report of SSD in a fossil strombid gastropod. It is argued that fecundity selection was the primary driving force behind the evolution of SSD in this gonochoristic gastropod species. Hence, the larger morph is the female.


Introduction
The family Strombidae is present in the Holocene exclusively in shallow tropical and sub-tropical seas. Wherever strombid species are found they are generally quite abundant [1]. The family is characterized by determinate growth. Once the shell 2. Dimorphic Persististrombus deperditus (Sowerby, 1840) Strombus deperditus Sowerby, 1840 was known from the Lower Miocene of Kutch and Pakistan [10,11]. Harzhauser et al. [12] included it (misspelled as depertitus) in the genus Persististrombus Kronenberg & Lee, 2007 and placed Strombus nodosus Sowerby, 1840 into its synonymy. The latter comes from the same stratigraphic level and geographical locality as it. The species is abundant in Kutch. Harzhauser et al. [12] provided a formal systematic description of the species.
The species is fusiform, anomphalous, dextrally coiled with elevated and step-like spire. The spire is shorter than the last whorl. The last whorl overlaps about 70% of the preceding whorl. The base is conical and extended into a short and slightly curved anterior siphonal canal. An angular shoulder separates a moderately wide and slightly sloping shelf from the rather flat whorl side. Its surface bears characteristic ornament having prominent spiral ribbing and shoulder tubercles. The species is easily identifiable because of these prominent morphological features. Interestingly, it occurs in a wide range of adult sizes in Kutch. On visual inspection, they can be broadly grouped into two size classes, which do not differ significantly in other features. The number of whorls in all the adult specimens is about eight. The larger and the smaller size groups are referred to here as the macro-and the microconch, respectively (table 1 and figure 1).
The height of the largest macroconch is twice that of the smallest microconch although overlapping is present in size ranges of the two groups. Therefore, it is difficult to assign some specimens to either of the groups. Hence, we have performed statistical analyses to see the viability of these groups.

Material studied
The Lower Miocene is represented in Kutch, Gujarat by the Khari Nadi Formation and the Chhasra Formation [13]. The Khari Nadi Formation was deposited in the Aquitanian and the Chhasra Formation represents the Burdigalian [14]. The succession comprises a thick pile of sandstone, shale and siltstone with intervening thin layers of marl and limestone. These marl and limestone bands are richly fossiliferous and mainly yield varieties of benthic molluscs-bivalves and gastropods. Persististrombus deperditus is quite common at certain levels of this succession mainly from the uppermost part of the Khari Nadi Formation and lower part of the Chhasra Formation. In this study, out of several hundred specimens in our collection we have used only 72 sub-adult to adult specimens. These are chosen because of preservation of major parts of their apical and basal sides so that faithful measurements can be taken. Several of these specimens are preserved with shell. The specimens were collected randomly from three localities. The specimens come from the Khari Nadi Formation exposed near village Aida (Loc 1: 23824 0 48.5 00 N, 68848 0 58 00 E), and Chhasra Formation exposed about 2.5 km north of village royalsocietypublishing.org/journal/rsos R. Soc. open sci. 6: 181320  2). All these specimens approach adulthood, which has been confirmed using one or more of the features discussed below.

Methodology
The dimension of a gastropod is traditionally represented by the shell height (figure 3), especially in the recent forms. In fossil strombids, either the apex or the anterior extremity of the siphonal canal or both are commonly broken. For this analysis, we have selected specimens with relatively less truncated ends to minimize error in the measurement of height. Still, measurement of the shell height (H) could have incorporated some error due to reconstruction. Hence, we have also measured the shell diameter (D) of the last whorl near the final position of aperture where the suture turns up adapically (figure 3). Measurements were taken using a digital slide caliper and used up to two decimal digits. Histograms were constructed with D and H data. The data distribution of designated macro-and microconch specimens were tested for normality by Shapiro -Wilk test ( p , 0.05). Equality of mean, variance and distribution were tested using t-, F-and Kolmogorov-Smirnov tests, respectively ( p , 0.05). Normal probability plots were also obtained for the data. We have performed discriminant analysis to get a quantitative evaluation of success of our visual classification into two size classes. All the statistical analyses and plots were obtained using the PAST 3.18 software [15].
10 mm 10 mm    Representative specimens were coated with magnesium oxide before they were photographed by Nikon D7000 DSLR camera (figure 1).

Determining adulthood
Adulthood in the family Strombidae is generally marked by an apical extension of the posterior canal and an upturned suture. Another feature commonly associated with the adult aperture in this family is a flared and thickened outer lip having plicate interior. In Persististrombus deperditus, apical extension of the suture is relatively short and rarely crosses the preceding whorl ( figure 1a,b,f,i,k,l; however, 1c for an exception). Flared outer lip with plicate interior is preserved in only a few specimens of our collection (figure 1b,c). However, a few other features facilitated recognition of the final size here. Strength of the shoulder tubercles and surface ornamentation reduce considerably about one whorl before the final aperture (figure 1g,h). This area, which lies at the same level with the adult aperture, was in contact with the substrate in the living individual. Apparently, this reduction of ornamentation is an adaptation for smooth movement of the adult individual. This is also economical because anti-predatory ornaments were not required in this part. The strength of the shoulder tubercles abruptly increases after this stage (figure 1g,h,k,n). This increase is accompanied by sudden appearance of spirally arranged tubercles at the middle of the whorl side (figure 1d,k,m). We have restricted our analysis to specimens where one or more of these features could be studied.

Results
Visual inspection of our collection of Persististrombus deperditus from the Lower Miocene of Kutch indicated large variability in the adult shell size and the presence of broadly two size classes. The larger size class, referred as the macroconch, has shell diameter of about 20 mm and above, and shell height of 50 mm and above. The smaller size group, referred as the microconch, comprises of individuals with D mostly less than 15 mm, and H less than 40 mm (table 1). The histograms that have been prepared from the data of H and D broadly show two-peak distributions (figure 4). This supports the division of the assemblage into the macro-and the microconch. The Shapiro -Wilk test reveals that H and D for both the macro-and the microconch specimens are distributed normally ( p , 0.05) (table 2). The normal probability plots arranged the data points for both the morphs largely along the reduced major axis regression lines, which were drawn for reference and comparison ( figure 5). This also exhibits normal distribution of shell height and shell diameter for the two morphs. The probability of their means, variances and distributions being equal appears to be negligible ( p , 0.05) (table 3). Because the variances of D and H for the two morphs can be significantly different, unequal variance t-test scores are also given in the table 3. However, p-values for both equal and unequal variance t-tests show significant difference between the macro-and the microconch. The discriminant analysis reveals that 95.83% of the specimens were correctly classified along the discriminant axis ( figure 6). Only three of the macroconchs are classified as microconchs in the discriminant analysis. H has a higher loading (4.4091) on the  royalsocietypublishing.org/journal/rsos R. Soc. open sci. 6: 181320 discriminant function than D (1.7803). The normal distribution, and significantly different means, variances and distributions indicate that the macro-and the microconch specimens represent natural populations, which differ significantly in their diameter and height.
The longitudinal data gathered from measurement of successive ontogenetic stages are given in table 4. PCA of the data and the plot of PC2 versus PC1 illustrate moderate degree of separation between the macro-and the microconch specimens along PC2 (figure 7a). PC1 accounts for more than 98% variability, whereas PC2 describes more than 1% (table 5). Distribution of loadings of the measured parameters demonstrates that PC1 is essentially the size axis, whereas PC2 is defined mainly by the ratio between Sp H and Wh D (table 5 and figure 7b,c). This ratio is plotted against Wh D ( figure 8). This scatter plot indicates that the microconch is more slender than the macroconch.

Discussion
The analyses reveal significant difference in the adult size between the macro-and the microconch, which were separated by visual inspection. The microconch is significantly smaller and also more slender than  Figure 5. Normal probability plots with correlation co-efficient for shell diameter of macroconchs (a), shell height of macroconchs (b), shell diameter of microconchs (c), and shell height of microconchs (d ) of P. deperditus (Sowerby). All the plots show strong correlation indicating normal distribution. the macroconch. Size dimorphism in a gonochoristic species is generally attributed to differences due to sexes. However, intraspecific size dimorphism due to non-sexual causes are known. Such dimorphism is generally linked to ecophenotypic differences due to difference in external physico-chemical and biological conditions of the habitat such as temperature, latitude, moisture content, food availability, predation pressure, wave-current energy, etc. [16][17][18][19][20][21]. These factors commonly vary with geography. However, within-habitat variation in these factors is possible depending on local differences. Here, we examine whether sexual or non-sexual causes were responsible for the evolution of size dimorphism in Persististrombus deperditus by trying to refute them one by one.

Sexual dimorphism?
SSD in fossil mega-invertebrates is best exemplified and almost omnipresent in the Jurassic-Cretaceous ammonites (Ammonitida, Cephalopoda). There, the antidimorphs differ in size and other morphologies in adulthood, whereas they share juvenile features. They are accepted as sexual consorts if found from the same stratigraphic level and overlapping geographical areas [22,23]. The study of SSD in ammonites has a long history [24][25][26][27]. The phenomenon is now almost unequivocally accepted by ammonoid researchers in spite of the fact that nothing is known about anatomical differences in the soft part of the antidimorphs [23]. The macroconch is often twice to four times the microconch in size. However, overlapping size ranges of the macro-and the microconch are known [22]. They also generally differ in adult apertural modifications, ornamentation and some shape parameters. However, in contradiction to the above statement, dimorphic pair with identical coiling and ornamentation is also known [28]. The relative abundance of the macro-and the microconch also commonly differs and may be in the range of 100 to one in either way. The difference in numerical ratios between the  Figure 6. Histogram of discriminant projection values for the two morphs of P. deperditus (Sowerby). There is slight overlapping between the two morphs defined by visual examination.   antidimorphs is generally attributed to facies and taphonomy [29]. The evolution and continuation of SSD in ammonoid lineages of the Jurassic and the Cretaceous appear to be phylogenetically linked [22,30]. Overwhelming and global presence of size dimorphism of similar nature in almost all the Jurassic-Cretaceous ammonite lineages and its evolutionary sustenance led to the large-scale acceptance of its sexual origin among ammonite researchers. SSD is present in all the living species of the family Strombidae wherever examined from this point of view. This fact prompted Cob et al. [9] to consider SSD the general condition of the lineage. SSD in living strombids has also been known for a long time [6][7][8]. In these living species, it is generally characterized primarily by moderate difference in size between the macro-and the microconch with overlapping in  their size ranges [1,2,31]. They also differ in shape, with males having relatively slimmer shells than females [5]. Males are also known to have thinner lips than females [1]. Females are also more abundant than males [4,5]. Strombid SSD is essentially similar to that of ammonoids except that it is characterized by relatively smaller difference in size between the sexes. Significant difference in ornamentation has also not been reported in any of the dimorphic strombids. SSD has never been reported in fossil strombids. Persististrombus deperditus resembles recent strombids in exhibiting moderate size dimorphism with some overlapping. Microconchs have slimmer shells than macroconchs. The lip thickness in most of the fossil specimens could not be properly studied and compared between the two varieties. The two morphs come from the same stratigraphic levels and from the same localities. Macroconchs are much more abundant than microconchs in all localities of collection. High abundance of the species wherever it was encountered, random orientation of its shells, low degree of fragmentation and abrasion, and presence of the specimens only in plane-bedded marl and limestone indicate that the fossils are autochthonous to parautochthonous. They apparently had substrate preference and flourished in abundance whenever substrate conditions were suitable. Independent determination of sex is not possible in P. deperditus as is possible in the living representatives of Strombidae. However, given the identical nature of dimorphism in it with those in living strombids, sexual origin of the size dimorphism in P. deperditus is difficult to refute. This also supports evolutionary sustenance of the identical SSD in living strombids.

Ecophenotypic difference?
Closely related species sometimes show distinct geographical pattern of morphological variation, which is commonly guided by latitude [18]. Intra-and interspecific size variations are known to occur in many animals, including molluscs, due to difference in ecological conditions of the habitats. However, significant size difference has also been reported between closely related species and between subpopulations of a species, which are separated by only about 10 m distance [18,19,21]. Johannesson et al. [19] reported wide variation in the size of Littorina saxatilis (Olivi, 1792) (Littorinidae, Gastropoda) population from Galicia. It is a direct developer species having indeterminate growth. It lives in the intertidal zone. They found large difference in size between individuals living in the lower shore and those in the upper shore [19]. There is a steep environmental gradient between the lower and the upper shores in terms of temperature, salinity, subaerial exposure, substrate, topography, wave energy, predation pressure and several other factors ( [21] and references therein). These differences are reflected in large differences in the shape and the size of species and sub-populations that thrive in these spatially close areas [18,21]. Faunal diversity is generally low in the intertidal areas because of strong fluctuation in ecological factors. But, there are a number of examples that record differences-including in size-between sub-populations living in sympatry in the intertidal habitat [18,19,21]. However, prominent non-sexual size difference between sub-populations living in sympatry are not known from other marine habitats.
Molluscs that lack a planktotrophic dispersive larval ontogeny are known to display more prominent differences between their populations in terms of size than those with planktotrophic larval ontogeny royalsocietypublishing.org/journal/rsos R. Soc. open sci. 6: 181320 ([19] and references therein). However, strombids with planktotrophic larval development are also known to vary significantly in size. Savazzi [2] observed that conspecific strombids living in the intertidal zone are much smaller than those living in the subtidal part a few tens of metres away. The intertidal and the subtidal areas differ markedly with respect to hydrodynamic condition, substrate type, temperature variation, food availability and predation pressure. Within the same strombid population, however, size varies only between 20 and 30% [2]. The specimens of Persististrombus deperditus from Kutch were mostly collected from limestone beds near the top of the Khari Nadi Formation and the basal part of the Chhasra Formation. These are biostromal composite shell concentrations ( [32] and references therein). The shells belong to a diverse benthic community that colonized repeatedly. The community flourished when siliciclastic supply was poor [32,33]. Lower part of the Khari Nadi Formation is dominated by fine sandstone and siltstone. These were deposited in the intertidal zone and are largely barren. With transgression and deepening, the basin became starved of clastic sediment supply. This allowed colonization of this relatively quiet and stable shelf by varieties of benthic animals. The shells suffered only minor within-habitat transportation [32]. Persististrombus deperditus is one of the common members of this community. They are not encountered in the lower part of the Khari Nadi Formation, which represents intertidal zone.
In summary, distinct intraspecific non-sexual size difference generally occurs (1) between sub-populations that live in upper and lower shores of the intertidal zone or (2) between populations that thrive in distinct habitats separated by some barrier-small or large. Persististrombus deperditus flourished in the Kutch basin only when transgression curtailed silisiclastic supply to the deepened shelf. The population belonged to a community that lived in a stable habitat, which was not segmented into distinct zones based on difference in ecological factors. Hence, the size difference between the micro-and the macroconch is difficult to explain by ecophenotypic causes.
It appears that large size variation in Persististrombus deperditus as ecophenotypic difference is more easily refuted than as sexual difference. Hence, difference in sexes is considered the most plausible reason for the development of size dimorphism in P. deperditus. Savazzi's [2] observation that fossil strombids of the same assemblage may display a size ratio of 4 : 1 can also include SSD.

Evolution of SSD in Persististrombus deperditus
The family Strombidae perhaps appeared in the Upper Cretaceous but became a common component of benthic mollusc faunas in shallow tropical and sub-tropical seas only from the later part of the Palaeogene [2,3]. The present report indicates that SSD, comparable to that of the living strombids, had appeared quite early in the history of this lineage. Another stromboid species belonging to the genus Tibia Rö ding, 1798 (Rostellariidae, Stromboidea, Gastropoda) co-occurs with Persististrombus deperditus in the Lower Miocene of Kutch. It also shows large variability in adult size in spite of a similar number of whorls in the adulthood, indicating the presence of possible SSD (K Halder & S Paira 2018, unpublished data).
In all examined living species of Strombidae, SSD is characterized by larger adult females [1,4,5,9,34]. Cob et al. [9] observed that while size difference between the sexes in Laevistrombus canarium is small, most of the sampled large individuals are females and the small individuals are males. In the four recent strombid species from the USA, live-caught specimens of which were used for captive breeding experiment, the female shells were observed to be higher and heavier than the male shells although their ranges overlap to a great extent (see table 2 in [31]). In P. deperditus we have found significant difference between the size of the macro-and the microconch although there is overlapping in their size ranges. Wherever sexual consorts of a living mollusc species are observed to have prominent difference in their adult size, the females are usually larger [35]. Fossil molluscs including ammonoids [22,23] and nautiloids [36][37][38] followed the same pattern of SSD. Hence, the macroconch of P. deperditus has greater probability of being the female.
It has been observed that the macroconch of P. deperditus is considerably more abundant than the microconch. Cob et al. [4] also observed abundance of females over males throughout the year of study in Laevistrombus canarium from Malaysia (see also [5]). Disparate sex ratios are a commonplace observation in sexually dimorphic species of ammonites, where it has been seen to vary from nearly one to 100 in either way [29]. Disparate sex ratio has also been observed in nautiloids [37,38]. In these cephalopods, the disparity has been attributed to seasonal segregation of the sexes or taphonomic bias [22,23,29,39]. Seasonal niche partitioning between the sexes is also known in living cephalopods [40]. It has been documented that in living strombids, sexes sometimes segregate during burial behaviour [41]. This may be a part of their reproductive behaviour pattern ( [42] and references therein). At the present royalsocietypublishing.org/journal/rsos R. Soc. open sci. 6: 181320 state of knowledge, it is uncertain whether the relative abundance of the morphs of P. deperditus is real or due to niche separation.

Why dimorphism: general remarks
SSD can evolve in an organism in response to one or more of the following causes.

Sexual selection
A male can father more offspring than a female can bear by mating with more than one female. This increases the reproductive success of the male but a female's reproductive success does not increase by mating with multiple males. This ensues in selection pressure on males. In response to this selection pressure males either fight among themselves to win over a female or mate a female coercively. As a consequence larger males are generally selected because strength is the key factor in this selection [43,44]. Larger males are also commonly selected in response to this selection pressure where epigamic choice by females is involved. There may be other responses to this selection pressure that result in smaller males. A smaller male with higher agility and manoeuvrability is selected because of greater success in search of mate, especially in fluid media [45,46]. Smaller males that mature early can spend maximum energy in reproduction [47]. Such progenetic growth of the male due to reproductive haste is common in molluscs [40,48].

Fecundity selection
The female grows to a larger size than the male and often differs in body shape to physically accommodate a large clutch of eggs. SSD in some bivalves and gastropods has been attributed to this factor [48][49][50].

Niche partitioning
If males and females occupy different ecological niches-either geographically or based on resource utilization-occupation of the eco-space increases and intersexual competition decreases [43,51]. As a consequence, sexes adapt and evolve independently. Such dimorphism is reported in molluscs. The dimorphism produced in such cases is generally seen to be characterized by larger females [52,53]. Niche partitioning itself does not explain the larger size of females. This may be in response to some other factor, such as dietary habits [54].

Non-adaptive causes
SSD can result as a by-product of change in some other character, e.g. feeding rates in juveniles or paedomorphic growth in one sex relative to the other [55]. It is worth mentioning that slower feeding rate and progenetic growth with early maturity in the male had been reported and explained in adaptive terms also [56].

Possible causes of SSD in Persististrombus deperditus
Sexual selection that selects larger males due to male-male fighting or involvement of coercion in copulation is not known in any mollusc or any other invertebrate organism. Epigamic choice by females is also primarily a phenomenon known in more complex groups of vertebrates. In the recent Strombidae, often several males including those that belong to different species are seen to attempt copulation with a female individual simultaneously [57,58]. A male attempting copulation with another male has also been observed [57,58]. Epigamic choice in such a group is almost impossible. Small size would not provide much agility and manoeuvrability to a benthic crawler, like a strombid gastropod. Although early maturity in males is common in molluscs, progenetic growth has not been reported in males of living Strombidea. Instead, neotenic growth of the male with respect to the female has been observed in Laevistrombus canarium [4]. Possibility of reproductive haste as the driving mechanism for evolution of SSD in P. deperditus seems low.
The larger and wider shell of the female gastropod to accommodate more eggs definitely provides higher success in producing large number of progeny. Hence, fecundity selection can be a probable cause of SSD in strombids. Cob et al. [9] also suggested that females of Laevistrombus canarium royalsocietypublishing.org/journal/rsos R. Soc. open sci. 6: 181320 allocated more energy to gonad production. Similar observations were made in other strombids and also other gastropods ( [5,35] and references therein).
Ecological niche separation between males and females leading to independent adaptation and evolution of the sexes can be a possible cause behind SSD in P. deperditus. This can also explain the significantly discrepant relative abundance of the morphs. However, information on living strombids is inadequate to support this hypothesis unequivocally. Living strombids are gregarious, and herbivorous or detritivorous in habit [2]. They commonly show a burial behaviour, during which sometimes sexes segregate [41]. But, this separation does not adequately explain the size difference.
Evolution of SSD as a fabricational noise (sensu [59]), such as by paedomorphic growth of one sex compared to the other, can only be established if such modification does not involve any adaptive advantage. In the case of living strombids with SSD the role of growth process could not be dissociated from that of adaptation.
It appears from the above discussion that fecundity selection is the most plausible cause for the evolution of SSD in Persististrombus deperditus. The only other factor that might have played some role was niche separation between the sexes.

Conclusion
SSD is present in all the living species of Strombidae wherever examined but was unknown in fossil strombids. The oldest SSD in the gastropod family Strombidae is reported here in Persististrombus deperditus from the Lower Miocene of Kutch, western India. The two morphs are distinguished mainly by the adult size, which is represented in this analysis by the shell diameter and the shell height. The SSD seems to have served the purpose of accommodating a larger clutch of eggs in the larger female. Discrepant relative abundance of the sexes, if not actual, might have resulted from niche separation between the two sexes. Appearance of SSD in such an early stage of evolution of the family suggests that further investigation may reveal more dimorphic species in this lineage.
Data accessibility. The specimens used in this study are deposited in the fossil collection of the Department of Geology, Presidency University, Kolkata, India. Measurement data used here are provided in the tables within the text. All the statistical analyses were done using the PAST 3.18 software [15].