# data[].attributes.body
parrotfish are iteroparous, breeding several times during their lives as different sexes, all year around. Reproduction occurs externally (broadcast spawners) and eggs are left to be spread by the tides where they hatch 25 hours after fertilization. There is no parental care involved

Hays, J. (2023). Parrotfishes: Characteristics, Behavior, Mucus and Sand-Creating Teeth | Sea Life, Islands and Oceania — Facts and Details. Factsanddetails.com. https://ioa.factsanddetails.com/article/entry-126.html

‌

# data[].attributes.body
Sexual Maturation

Some species of parrotfish are hermaphrodite protogeny. Rivulated parrotfish for example go through the early stages of life as females, between the length range of 126-270mm. between the lengths 246-350mm, they turn into males. Female rivulated parrotfish reach sexual maturity at 2 years and undergo a period of sex change to male sex when they are 3 years old

Tuwo, A., Tresnati, J., Huda, N., Yasir, I., Rahmani, P. Y., & Aprianto, R. (2021). Reproductive strategy of rivulated parrotfish Scarus rivulatus Valenciennes, 1840. IOP Conference Series: Earth and Environmental Science, 763(1), 012002. https://doi.org/10.1088/1755-1315/763/1/012002

‌

Mate Selection

Among other parrotfish species such as bumphead parrotfish mate in time with the lunar cycle, possible to coincide their breeding events with tides to ensure their eggs are pushed out to open ocean and away from predators. During breeding events, males position themselves at designated territories and court the females by swimming in bursts and temporarily changing colour. Females, in response, turn their faces white and the pair take their affair near the surface where they release their eggs for fertilization.

Dell’Amore, C. (2015, March 8). Head-Bashing and Other Mating Secrets of Giant Parrotfish. Animals. https://www.nationalgeographic.com/animals/article/150308-parrotfish-science-oceans-animals-mating-moons

‌

# data[].attributes.body
Most species spawn at dawn or dusk to reduce the risk of predation.

Spawning can be pair spawning (one male, one female) or group spawning (one male with multiple females).

Fertilization is external. Females lay their eggs into open waters and males fertilize the eggs by releasing their sperms in their vicinity.

Relationships & Adaptations of a Parrotfish. (2020, November 19). Pets on Mom.com. https://animals.mom.com/relationships-adaptations-parrotfish-8336.html

# data[].attributes.body
most are small in size (can range from 0.7mm - 2.6mm in diameter)

buoyant due to a single oil globule and high water percentage maintaining buoyancy

have higher mortality rates due to fragility of eggs being exposed to higher predation and exposure to harsh environmental conditions

have higher fecundity due to species usually having more spawning seasons

all pelagic eggs are transparent and mostly spherical (exceptions include anchovies which are longer, or others which may be slightly ovoid)

eggs do not consist of adhesive membrane

type of eggs produce do not link between habitat and species of fish (pelagic fish can produce demersal eggs, demersal fish can produce pelagic eggs etc.)

two types of pelagic eggs:

isolated eggs (majority)

agglutinated eggs (lophiidae): eggs form a mass and move together

reference:

IASRI. (n.d.). Fish eggs and their characteristics. eCourses Online. http://ecoursesonline.iasri.res.in/mod/page/view.php?id=38791

# data[].attributes.body
demersal eggs are usually heavy and dense and have filamentous structures to stick to other objects

usually have a thick and more complex chorion and is more resistable to damage

usually has a higher lipid content

develop slowly, can take more than a month to hatch

type of eggs produce do not link between habitat and species of fish (pelagic fish can produce demersal eggs, demersal fish can produce pelagic eggs etc.)

two types of demersal eggs:

adhesive: (exocoetidae, gobiidae) they are adhesive to substratum using adhesive egg membranes / filaments

non-adhesive: (salmonidae)

some examples of species with demersal eggs:

deep sea wolf herring (demersal fish)

references:

Blaxter, J. H. S., & Hempel, G. (1988). Buoyancy, egg size and predation in marine fish eggs. Sarsia, 73(1), 13-21. https://www.tandfonline.com/doi/pdf/10.1080/00364827.1988.10420671#:~:text=Demersal%20eggs%20usually%20have%20a,a%20relatively%20high%20lipid%20content.

IASRI. (n.d.). Fish eggs and their characteristics. eCourses Online. http://ecoursesonline.iasri.res.in/mod/page/view.php?id=38791

# data[].attributes.body
the action of species releasing or depositing reproductive gametes resulting in fertilisation for reproduction by sexually mature male and female species

most species have spawning seasons or exhibit seasonal cycles in the production of gametes during specific seasons of the year. this is usually called breeding season for species

before spawning, female species usually have maturity stages before the release of eggs for fertilisation. stages include resting phase, early maturing phase, advanced maturing phase, mature / prespawning phase, spawning phase and spent phase.

reference:

The Fish Site. (n.d.). Maturation and spawning in fish. Retrieved from https://thefishsite.com/articles/maturation-and-spawning-in-fish

# data[].attributes.body
when an individual starts as a female sexually mature producing female gametes but changes into a male sexually mature producing male gametes

favoured in situations where if large females in a small group changes sex into a large dominant male to control reproductive access to females

reference:

Todd, E. V., Liu, H., Muncaster, S., & Gemmell, N. J. (2016). Bending genders: The biology of natural sex change in fish. BioEssays, 38(11), 1117–1135. https://pmc.ncbi.nlm.nih.gov/articles/PMC5567342/#:~:text=Unique%20among%20vertebrates1%2C%20sex,large%20dominant%20males%20(protogyny).

# data[].attributes.attachment_link.readable_text
Abstract

Sequentially hermaphroditic fish change sex from male to female (protandry) or vice versa (protogyny), increasing their fitness by becoming highly fecund females or large dominant males, respectively. These life-history strategies present different social organizations and reproductive modes, from near-random mating in protandry, to aggregate- and harem-spawning in protogyny. Using a combination of theoretical and molecular approaches, we compared variance in reproductive success ( V k*) and effective population sizes ( N e) in several species of sex-changing fish. We observed that, regardless of the direction of sex change, individuals conform to the same overall strategy, producing more offspring and exhibiting greater V k* in the second sex. However, protogynous species show greater V k*, especially pronounced in haremic species, resulting in an overall reduction of N e compared to protandrous species. Collectively and independently, our results demonstrate that the direction of sex change is a pivotal variable in predicting demographic changes and resilience in sex-changing fish, many of which sustain highly valued and vulnerable fisheries worldwide.

Introduction

Unique among vertebrates1, sex-changing fish develop and reproduce as males first and then grow into highly fecund females (protandry), or reproduce initially as females to later change into large dominant males (protogyny). Sequential hermaphroditism has intrigued evolutionary biologists for decades and a great amount of information has been gathered on sex determination, sex differentiation and the plasticity of sex change in fishes2-5. The main theoretical model proposed to explain its adaptive value, the size advantage model6-8, predicts that sex change should occur when the reproductive success of an individual depends on its size, but more so for one sex than the other. In this scenario, protandry is favoured over fixed separate sexes (gonochorism) when larger females have higher reproductive value than smaller ones (they can produce more eggs), while protogyny is favoured in situations where size allows dominant males to control the reproductive access to females.

At the population level, we can expect that the variance in individual lifetime reproductive success ( V k*) will influence the demographic trajectory of a population (Fig. 1). Strangely, sex-changing populations have seldom been investigated from a population genetic perspective. In those circumstances the focus has been only on one sex-changing mode, mainly protogyny9, 10 or on the comparison between gonochoristic and sex-changing species11, 12. Yet, among sequential hermaphrodites, protandry and protogyny stand as two remarkably different life-history strategies as they are shaped by different social systems and reproductive modes3, 5, 13. In protandrous species, populations are composed of many small males and fewer large, highly fecund females. The reproductive mode is often monogamy or near-random mating2, 6. Protogyny, on the other hand, occurs when there is high potential for polygyny and results in strong social structures dominated by large males6, 8, in some cases antagonised by sneakers. Males are territorial and control harems of females or, in group-spawners, larger males are expected to be the most successful. In either case, in protogyny there is a strong sexual selection on males.

Figure 1.

Interplay among individual fitness, life-history traits and population dynamics. Multiple factors affect effective population size ( N e), including many driven by the mating context of the population. The exemplified protogynous species is drawn in red and the protandrous one in blue; shadows indicate the sex of individuals at any point during their lives (red for females; blue for males). In the equations: N m = number of adult males; N f = number of adult females; N = total number of adults in the population; GL = generation length; V k* = variance in individual lifetime reproductive success.

Mating system and reproductive mode variations are predicted to influence V k*14 (Fig. 1). In particular, V k* should be greater in protogyny than in protandry, as large males can monopolize multiple females and small females will tend to choose the larger males, increasing the reproductive success of a small number of larger males.

As a result of such changes over an individual's lifetime, sequentially hermaphroditic species typically exhibit sex ratios that depart from the balanced ≈ 0.5 observed in gonochoristic species, and are generally skewed towards the 'first sex'. Such a bias is more pronounced in protogynous species than protandrous ones15. Variance in reproductive success and skewed sex ratios are the two most powerful forces that shape a parameter of crucial significance in population genetics, ecology and conservation: the effective population size ( N e), which offers a view of the intensity of genetic drift and the changes in genetic variability in a population and its potential for persistence and resilience16, 17. Wright18 was the first to realize that the "effective population number" should refer only to the "breeding population and not to the total number of individuals of all ages"; indeed, the effective population number in natural populations is generally smaller than the census size of the same population, sometimes by one or more orders of magnitude.

The direct influence of skewed sex ratio applies mainly to gonochoristic species. In sequential hermaphrodites, the sex ratio is dynamic and many individuals contribute to the coming generations as both sexes (even though not all individuals change sex). Biased sex ratios though reinforce reproductive skew, and higher V k* values indicate that female-first sex-changers should have lower N e than male-first sex-changers (Fig. 1). Thus, the direction of sex change is expected to have an influence on N e and, consequently, on population structure.

Calculating N e is extremely difficult because it requires a specific knowledge of the demographic and ecological parameters of the natural population under study. As mentioned, N e can be influenced by a multitude of factors, including fecundity, birth rate, natural mortality, migration, breeding sex ratio, variance in family size, age-structure, spatial and temporal distributions and reproductive mode19. For this reason, direct demographic methods of calculating N e are often replaced by indirect genetic methods20.

Here, we used a combination of theoretical modelling and molecular genetics to compare the effective population size in several species of protandrous and protogynous fish. A life-history model21 was employed to estimate N e, V k*, the effective number of breeders per year ( N b), as well as the annual mean number of offspring () and V k and N b for each sex in two sets of protandrous and protogynous species. At the same time, molecular markers were used to produce empirical estimates for a similar set of protandrous and protogynous species. The combination of these two approaches largely confirms the main expectation of lower N e estimates for protogynous than for protandrous species, and collectively unveils previously unrecognised patterns and trends that are of great relevance to the management of marine living resources.

Results

Life history modelling

A range or realistic life-table scenarios were produced, utilizing data from literature on eight species, four protandrous (PA) and four protogynous (PG; Table 1), to estimate21 multiple key population parameters (Table 2). As predicted, protogynous species overall showed significantly lower N e (mean ± sd: PG = 215.73 ± 28.06; PA = 392.70 ± 83.28; t = 4.03, df = 6, p = 0.006; Fig. 2) and higher lifetime V k* (mean ± sd: PG = 191.61 ± 108.59; PA = 78.90 ± 12.88; W = 0, p = 0.029; Fig. 2) than protandrous species, in the face of a stable total adult N (adult census size) that did not differ significantly between the two groups (mean ± sd: PG = 753.00 ± 213.41; PA = 603.50 ± 120.87; t = −1.22, df = 6, p = 0.280; Fig. 2). Similarly, no significant difference was detected in generation length (mean ± sd: PG = 10.24 ± 5.02; PA = 7.88 ± 1.94; t = −0.87, df = 6, p = 0.416). Indeed, the linear model returned protogyny (t = −7.280; p < 0.001) as the main factor reducing the size effective ratio N e/ N, which allows to control for the variance in population abundance estimates ( N).

Table 1.

Species

Family

Reproductive mode

AgeMax

L inf

t 0

K

t m

α

β

Protandry

Diplodus sargus

Sparidae

Broadcast spawners

1263

45.90

−0.890

0.171

3

0.011

5.1064

Lithognathus mormyrus

Sparidae

Broadcast spawners

1165

38.4466

−1.48366

0.20066

3

0.010†

5.10 †

Sparus aurata

Sparidae

Broadcast spawners

1267

59.7667

−1.71167

0.15367

2

0.010†

5.10 †

Lates calcarifer

Latidae

Broadcast spawners

23[32]68

143.00

−0.860

0.130

4

27.00069

2.8969

Protogyny

Chrysoblephus puniceus

Sparidae

Haremic

1012

44.0070

−0.810

0.180

2

10.253*

2.6071

Pagellus erythrinus

Sparidae

Group spawners

2131

41.7872

−1.210

0.130

3

0.030 ‡

5.49 ‡

Spondyliosoma cantharus

Sparidae

Haremic and nest guarders

1073

47.70

−0.830

0.180

2

0.04074*

4.6074*

Semicossyphus pulcher

Labridae

Haremic and territorial

2947

55.77

−0.710

0.126

4

0.00189

5.4975

Table 2.

Species

Total N

N

GL

V k*

N e

N e/N

N b

N b/N

N b/N e

t c

♀ k ®

♀ V k

♀ N b

♂ k ®

♂ V k

♂ N b

Protandry

Diplodus sargus

1875

560

7.449

85.960

338.7

0.605

251.9

0.450

0.744

6

7.650

32.510

91.7

0.573

3.089

201.3

Protogyny

Chrysoblephus puniceus

1755

458

5.810

111.782

204.2

0.446

155.9

0.340

0.763

5

0.588

1.844

366.9

18.170

18.170

43.6

Figure 2.

Graphical representation of some key population parameter estimates (obtained with AgeNe) by individual species (panel a) and mating system (panel b); protandry in blue vs. protogyny in red. N: estimated number of individuals in the population (census size); V k*: lifetime variance in reproductive success; N e: effective population size; N b: effective number of breeders per year. Key ratios ( N e/ N, N b/ N and N b/ N e) are also represented.

As all the species under study are iteroparous, N b (effective number of breeders per breeding cycle) was also calculated and found to be significantly lower in protogynous than in protandrous species (mean ± sd: PG = 155.98 ± 20.91; PA = 348.20 ± 129.91; W = 16, p = 0.0286; Fig. 2). The linear model returned protogyny (t = −6.179, p = 0.003) as the main factor reducing N b/ N (again, to control for the variance in population abundance estimates, N), with maximum length as a covariate, which significantly increased this size effective ratio (t = 2.954, p = 0.042).

The analysis of the N b/ N e ratio revealed no influence of the mating system (t = −2.074; p = 0.093), but a significant effect of maximum length (t = 4.886; p = 0.004). The male-first sex-changer L. calcarifer showed a ratio larger than 1 ( N b/ N e = 1.059), meaning that in this species the number of breeders per year is actually higher than the overall effective population size (Fig. 2).

It was also possible to analyse , and annual V k and N b for each sex. Graphically (Fig. 3) it is easy to visualize the reproductive strategy for an individual when it reproduces as a different sex during its lifetime, changing from the first sex (male in protandry and female in protogyny) to the second sex (female in protandry and male in protogyny). Individuals belonging to protandrous species have higher annual (mean ± sd: ♀  = 7.17 ± 2.26; ♂ = 0.55 ± 0.06; = 5.06, df = 6, p = 0.002) and V k (♀ = 22.70 ± 14.80; ♂ = 2.24 ± 0.75; t = 5.06, df = 6, p = 0.002) when they are older and larger (second sex, female, Fig. 3a). For this latter stage of life (the female phase), in the population there are fewer breeders per year than for the male phase ( N b♀ = 124.20 ± 52.39; N b ♂ = 299.00 ± 87.71; t = −3.561, df = 6, p = 0.012). A mirrored situation applies to protogynous species, and in this case the magnitude is even stronger (Fig. 3b): larger males second sex) have much higher annual (♂ = 17.82 ± 4.04; ♀ = 0.51 ± 0.09; Welch t-test t = −24.892, df = 5.87, p < 0.0001), higher V k (♂ = 78.03 ± 51.43; ♀ = 1.91 ± 0.36; W = 0, p = 0.029), and lower N b than smaller females, the first sex (♂ = 44.83 ± 7.00; ♀ = 308.78 ± 41.25; W = 16, p = 0.029). The even greater mismatch between the number of breeders for the two sexes reflects nicely the description of haremic species, where few males monopolize the majority of females.

Figure 3.

Graphical representation of some annual key parameter estimates (obtained with AgeNe) in each life history strategy, per sex. k ¯ : annual mean number of offspring; V k: annual variance in reproductive success; N b: annual effective number of breeders. Values are log ₁₀ transformed to improve visualization. Protandry in blue (panel a); protogyny in red (panel b); lighter colour: females - F; darker colour: males - M.

In general, the trend is similar between protandry and protogyny: the second sex has always higher and V k (which confirms the individual advantage of a change of sex in both systems) and lower N b than the first sex (Fig. 3). When comparing the reproductive success of first and second sex across the two systems, it is interesting to see a significant interaction (mating system [PA, PG] crossed with sequential sex [first, second]: t = 26.06; P < 0.0001; Fig. 4): the second sex is comparatively more successful in protogyny, whereas the first sex is slightly more successful in protandry. Protogynous male breeders (second sex in red in Fig. 4) are more successful (but less numerous) than protandrous females (second sex in blue in Fig. 4) but protogynous females (first sex in red in Fig. 4) are not more successful than protandrous males (first sex in blue in Fig. 4), even though at this stage the fitness difference between the first and second sex is not as large as it will become later in life, after sex change.

Figure 4.

Interaction plots of some annual key parameter estimates (obtained with AgeNe) in each life history strategy per sex. k ¯ : annual mean number of offspring; V k: annual variance in reproductive success; N b: annual effective number of breeders. Values are log ₁₀ transformed to improve visualization. Protandry in blue; protogyny in red.

Genetic evidence

We had generated genotypic data from five sparid species: three protandrous, Diplodus sargus, Lithognathus mormyrus 22 and Sparus aurata 23, and two protogynous, Chrysoblephus puniceus 12 and Pagellus erythrinus 24. Furthermore, two more protogynous species, Cephalopholis fulva and Epinephelus guttatus, were included in the analyses, based on datasets made available by Portnoy and colleagues9. No significant signals of large allele drop out, scoring error were detected across each species. Only C. puniceus had one marker (SL35) that had a strong probability for null alleles, and was hence removed from the statistical analysis. Average observed ( ﻿H﻿o)﻿ and expected (﻿ H e)﻿ heterozygosities were respectively 0.756 and 0.823 for D. sargus; 0.844 and 0.836 for L. mormyrus; 0.768 and 0.776 for S. aurata; 0.532 and 0.534 for C. fulva; 0.701 and 0.860 for P. erythrinus; 0.794 and 0.829 for C. puniceus and 0.690 and 0.694 for E. guttatus (Table 3; see Table S1 for genetic diversity parameters by locus and by population). Two protogynous species, P. erythrinus and to a lesser extent C. fulva, showed signs of heterozygotes deficiency, with values of H o smaller than H e. This translates into high and positive F is values recorded especially for P. erythrinus. Overall, H o was on average higher for the protandrous than the protogynous species (W = 99, p = 0.007; Fig. S1), while no significant difference was found for H e (W = 79.5, p = 0.171) or A R (W = 43, p = 0.324). Population structure was also investigated using F st, the Bayesian clustering implemented in STRUCTURE (Fig. S2) and the Discriminant Analysis of Principal Component, DAPC (Fig. S3; Fig. S4). Although some species did show low but significant pairwise F st values, overall only L. mormyrus revealed the presence of substantial population structure, with both analyses. This was expected as the dataset contains Atlantic and Mediterranean populations, which are separated by a strong phylogeographic break22. Cephalopholis fulva and S. aurata also showed some signal of low but significant differentiation with pairwise F st (Table S2), although STRUCTURE failed to identify any sub-structuring. We decided to estimate per single population/location, rather than pooling them together. Estimates of calculated with LDNe, varied between 128 (Cp2) and infinite (Lm1, Lm2, Sa1, Sa3, Pe1, Eg1; Table 3; Fig. 5a). Given the high skew towards large values (i.e., infinite estimates), we calculated also 1/ (Fig. 5b). Overall, we did not detect any significant difference when we compared 1/ across all the populations classified by sex-changing system (protogyny vs. protandry: W = 44, p = 0.345), but, given the high variance (Fig. 5c) recorded for protogynous species, we repeated the analysis using their reproductive mode, and found that haremic species stand out as having significant higher 1/ than the other species (random mating vs. group spawning vs. harem: Kruskall-Wallis χ ²  = 8.0101, df = 2, p = 0.018).

Table 3.

Species

Population

ID

n

loci

H o

H e

F is

A R

N ˆ e

Protandry

Diplodus sargus

France - Sète Italy - Livorno Croatia - Rovinj

Ds1 Ds2 Ds3

49 49 50

10 10 10

0.738 0.754 0.776

0.822 0.826 0.831

0.121 0.089 0.061

15.094 15.222 14.428

239.5 (129.3-1155.7) 1447.1 (293.7-∞) 556.3 (205.9-∞)

Protogyny

Chrysoblephus puniceus

South Africa - Port Edward South Africa - Park Rynie South Africa - Richards Bay

Cp1 Cp2 Cp3

43 43 39

11 11 11

0.782 0.810 0.790

0.827 0.831 0.830

0.048 0.025 0.040

18.005 18.161 18.00

834.9 (286.7-∞) 128.2 (96.3-186.5) 165.8 (113.8-291.6)

Figure 5.

Effective population size () estimates, obtained with LDNe for multiple populations of the seven species under study. Bars in the tones of blue: protandrous species; bars in the tones of orange: protogynous species. ( A) values (up to infinite); ( B) 1/ (note the reverse y axis as an aid to visualize values - higher 1/ values imply lower ); ( C) 1/ (reverse y axis) by reproductive mode.

Discussion

Sequentially hermaphroditic fish are an intriguing group of animals. Their life history includes, for the majority of the individuals in a population (the exception being primary females in protandry and primary males in protogyny), a complete rearrangement of gonads and behaviours which allow them to reproduce as the opposite sex, at a certain stage of their life2.

While fecundity and reproductive success are expected to increase with size for both females and males, under certain circumstances (mainly related to mating context3, 5), the fitness advantage can increase more rapidly for one sex than the other. In this case, it is beneficial to match the sex with the higher reproductive value at a given size, in order to maximise the reproductive output6-8. Since the sexual-transition is mating-system dependent, there is not a single direction to change sex, but two, from male to females or vice-versa; in our study we have not considered bidirectional sex-changers25 given their low occurrence.

To our knowledge, no studies have attempted to calculate and compare the fitness of each sex stage in both sequentially hermaphroditic models (protandry and protogyny). Using a theoretical model based on life history traits, we were able to compare sex-specific reproductive outputs in eight sex-changing species (four protandrous and four protogynous), and their corresponding annual variance. As predicted by theory, in both systems (protandry and protogyny) the second sex is more successful than the first one in terms of the annual average number of offspring, and experiences higher variance in reproductive success, V k. Also, in both cases, the second sex is composed by a smaller number of annual breeders ( N b) compared to the first sex (Fig. 3). The magnitude of these values is different between the two systems and this reflects the underlying reproductive systems that characterize them: in protogyny, sex-changed males (mostly haremic) face a large increase in and V k compared to their previous female stage (and there are less male breeders: they are few and dominant), much larger than the increase occurring during a protandrous transition (Fig. 3).

Recently there has been an increased emphasis on the connections between individual life history traits and population-level responses26, 27 which results in a link between demographic and evolutionary processes28. Lifetime V k* can be translated to the population level, which allows us to address a broadly relevant question, especially in light of conservation and management of fish stocks and biodiversity: are the population trajectories of protandrous and protogynous species different? Our estimates confirmed that the V k* of protogynous species is significantly higher than the V k* of protandrous ones (as hypothesised: many protogynous haremic species have few larger males who can successfully fertilize the majority of the females; Fig. 3). As a consequence, we found that protogynous species have significantly lower N e than protandrous ones (while having similar number of adults in the population), thus supporting the initial hypothesis that the direction of sex change plays an important role in shaping the demographic trajectories of fish population under natural mortality (our life history tables did not include fishing pressure).

For iteroparous species, it has been suggested to estimate the effective number of breeders per reproductive cycle ( N b), which is easier to calculate than the effective population size28, 29. The ability to change sex adds an extra level of complexity to the analysis of N e in organisms which grow indefinitely, reproduce multiple times in their lives and have overlapping generations: in this case the focus on N b is more relevant than the number of reproducing individuals per generation ( N e), as the same individual in different years will reproduce using a different strategy. We thus calculated overall N b estimates and a series of key effective size ratios. By definition, protogynous species have highly skewed sex ratios (unbalanced towards the first sex), which should result in low numbers of male successful breeders ( N b). Also, N e/ N ratios are expected to be low in species with high fecundity but high V k*30 and N e/ N should approximately be similar to N b/ N ratios30. All these theoretical expectations were confirmed in our analyses: N b values, as well as N e/ N and N b/ N ratios were found to be significantly lower in protogynous than in protandrous species. We did not detect differences in N b/ N e in the two groups, but in one case, the protandrous barramundi ( L. calcarifer), N b/ N e exceeded 1. According to theoretical and empirical estimates, a long life span leads to higher values of N b and delayed maturation leads to an increase of N e 28. The barramundi is very long lived but has an early maturation, which explains the somewhat paradoxical (yet not unique28) value obtained.

We were able to strengthen our understanding of the consequences of changing sex in opposite directions by testing the central hypothesis (protogynous species have lower N e than protandrous species) also using independent molecular datasets on multiple populations of five sparid species: higher estimates were obtained for the protandrous species (with the exception of one population of D. sargus), but not all the estimates from the two protogynous species available (which do not fully correspond to the species used for the life-history modelling) conformed to the expectation of low values: indeed, while C. puniceus (with a lek-like system where males defend territories) presented low for the majority of the populations, both populations of P. erythrinus presented high value of , comparable with the ones obtained in protandrous species (Fig. 4; Table 3). Even with estimates obtained from life history tables, the N e estimate for this species was higher than the other species under consideration (Fig. 2; Table 2). Pagellus erythrinus are group-spawners: their V k* is not as high as haremic species, but a protogynous sex change is still advantageous in this case: being a large male can be beneficial in the spawning area31, 32, where the lack of haremic structure removes the typically reduced level of sperm-competition which is the norm in protogynous species33. Furthermore, the pronouncedly lower-than expected V k in P. erythrinus, may to some extent stem from increased levels of female fecundity and/or k ® , or a possible unexplored greater contribution of primary males34, which the standardised approach we employed (see Methods) would not be able to account for.

To test the possible interplay between sex-changing system and reproductive mode, we used two more datasets available for two protogynous groupers, the group-spawner E. guttatus and the haremic C. fulva 9: indeed, the former exhibited values similar to P. erythrinus, while the latter showed estimates comparable to C. puniceus, thereby revealing an additional layer of complexity to the scenario: we detected a decrease of along a gradient of reproductive modes (Fig. 4) and found that harem-spawners protogynous species have lower than group-spawning ones. Thus, behavioural traits contribute to a fuller understanding of the patterns of N e and V k* in sequential hermaphrodites.

Overall, the estimates from life history tables confirmed the theoretical expectation that, in general, regardless of the direction of sex change, individuals use the same overall strategy, producing less offspring () with less variance in reproductive success ( V k) as the first sex and more offspring with higher V k as the second sex, later in life. Concurrently, at the population level, the majority of breeders ( N b) reproduce as the first sex. Although this trend is consistent in the two mating strategies (protandry and protogyny), its magnitude is much higher in protogynous species, resulting in a stronger reduction of N e at the population level. Thus, the direction of sex change has an influence on the overall N e of the populations and combining all sex-changing fish species in one category is ultimately incorrect and can be misleading for conservation and management practices.

A second independent analysis, based on molecular data also indicated that a protogynous mating system is prone to more significant reductions of than a protandrous system. Moreover, this second analysis uncovered the fact that, in protogynous populations, V k* may also be dependent on reproductive modes (from aggregate- to harem-spawning) and can change accordingly. The comparison, within protogynous sex-changing modes, between group-spawners and haremic species confirms the importance of considering and assessing the influence of mating systems and reproductive modes, in order to understand the demography of populations35. The possibility that several natural stocks of haremic protogynous species may owe their persistence to a disproportionately small number of extremely successful old males fits with recent concerns regarding many valuable warm-water fisheries worldwide36, and calls for the practical consideration, not only of sex change, but also of its direction and the behavioural structure underpinning it, in the management arena. Future analyses should focus on the evaluation of the consequences of increased mortality by fishing pressure and the populations' compensatory capacities37. Following this path can significantly expand our ability to predict the viability of many commercially exploited reef-associated species in tropical and warm-temperate areas, and enable us to do something useful for their preservation.

Materials and Methods

Choice of species

Four protandrous and four protogynous species were selected (Table 1), based on existing knowledge of their biology, to produce realistic life history tables. Six of them belong to the Sparidae, a very heterogeneous family in terms of reproductive styles38, which includes both protandrous and protogynous species. We previously generated genetic data for five out of these six sparid species, which we used in our molecular analysis (Table 3), in the attempt to minimize the biases that may arise comparing phylogenetically distant taxa. Two more protogynous species9, characterized by different social structures (aggregate- and harem-spawners) were also included in the analysis.

Theoretical modelling

Key life history parameters (Table 1) were retrieved from the literature and fed into a demographic model21, implemented in the freely available software AgeNe, version 2.0 ( http://conserver.iugo-cafe.org/node/2876). For each species, Leslie matrices were built and age-specific survival ( s x), fecundity and sex ratios were calculated.

Age-specific survival rates

Natural mortalities ( m) were calculated for length classes following Charnov and colleagues39

where L inf is the asymptotic length (in cm), L is the fish standard length (in cm) and K is the growth coefficient of the von Bertalanffy equation:

where L x is the fish total length (in cm) at age x, and x ₀ is the theoretical age at length = 0. Length classes were assigned to age classes (eq. 2) and length values for each class were compared with known data from the literature, to ensure realistic estimates. Natural instantaneous mortalities were then transformed in annual mortality rates per year (d40) using the formula:

To prevent survivorship in the oldest age class, age-specific survival rates were calculated as:

s x =s A −(0.5/ x)−{[e (x−s A )]⋅s A −(0.5/ x)}

4

which is based on a generally accepted type III survivorship curve, where s A is the asymptotic maximum survival ( s = 1 −  d = e −m, from eq. 3) and x is the age. This introduced senescence reduces the adult population to zero at the maximum age (as reported in the literature). Equation 4 adjusts survival for age, following the rationale that natural mortality increases exponentially with age24, 41. Thus, survival increases with size until at a certain age and then decreases in older fish till all fish are dead in the oldest age class. Based on greater availability of bibliographic information, combined data for the two sexes were used. Even though this is a simplification42, the assumption of size-dependent (and not sex-dependent) survival rates should not significantly alter the results. Indeed, for the majority of individuals (excluding the ones which do not change sex, i.e., primary males in protogyny and primary females in protandry) sex is correlated with size. This means that male survival can be assumed to be the same as female survival, but this value changes with the size of the fish39, 43. Average mortality data (across all size classes) were compared to known mortality rates from literature, to ensure realistic estimates (Table S3).

Age-specific fecundity

Female fecundity can be calculated with an allometric (power function) relationship44-46

Female fecundity(N e g g s )=α L β

5

where α and β are constant. Surprisingly, not many data on fecundity are available from the literature. We used the data available and extrapolated the missing data, comparing calculated ranges of fertility with known ranges from literature (when available), to ensure realistic estimates (Table S3). From a practical standpoint (count of gametes), it is more difficult to calculate male fecundity. Following Alonzo and Mangel45, 46 we considered male fecundity as 1000 times female fecundity.

Age-specific sex ratio in the population

The sex ratio in the population changes with the size of the individuals46 and thus it is different among age-classes. Individuals change sex at a specific age ( t c). The proportion of the terminal sex ( S; males in protogynous species and females in protandrous species) was calculated with the formula:

S=S i +(S f −S i )/{1+e [−(t−t c )]}

6

where S i is the proportion of non sex-changing individuals (primary males in protogynous species and primary females in protandrous species) and S f is the final proportion of the terminal sex (overall sex ratios for each species are reported in Table S3).

Given these age-specific parameters (age-specific survival rate s x; birth rate b x and sex ratio), AgeNe can produce complete life history tables for each sex, with estimates of total number of individuals in the population (total N) and estimated number of adults (adult census size N), generation length ( GL i.e., the average age of parents of a newborn cohort, in years), lifetime V k*, N e, and annual number of breeders ( N b). Moreover, the software estimates annual V k, annual mean number of offspring (k ¯ ) and annual N b for each sex. We searched published literature to obtain realistic parameters for our models. We are aware that a great deal of local variation exists in these parameters47, 48 and we acknowledge that geographically targeted predictions should be made for individual populations. Here, we used data from FishBase49 to insure general (and not locality-specific) estimates (Table 1). We checked that our life tables reproduced realistic scenarios and we made changes only when recent articles reported higher maximum length than originally reported in FishBase; Table S3). Our goal was to obtain a general (and not local) estimate of N e, for comparative purposes among species characterized by different mating systems.

Molecular analysis

Samples were genotyped using microsatellite loci (Table S4) run on ABI 3130xl Genetic Analyzer capillary sequencers (Applied Biosystems©). Markers were amplified using fluorescent primers labelled with 6-FAM, NED, PET and VIC dyes (Applied Biosystems©). For S. aurata, we used 12 of the 15 loci used in the original study23, since we removed three possible candidates for directional or balancing selection. Data were analysed with MICRO-CHECKER50 to check for large allele drop out and scoring errors; frequencies of null alleles were estimated using the software FreeNA51. We retained loci with moderate frequency (0.05 < r < 0.20) of null alleles51. Estimates of genetic diversity, including expected ( H e) and observed ( H o) heterozygosities, allelic richness ( A R) and coefficient of inbreeding ( F is), to test for departure from Hardy-Weinberg equilibrium, were calculated using the package hierfstat52 while pairwise F st values and relative 99% confidence intervals were calculated using the package DiveRsity53 both developed in R54. We also checked for the presence of population sub-structuring using two different approaches: STRUCTURE 2.3.455, 56 and the Discriminant Analysis of Principal Component (DAPC) as performed in the R54 package ADEGENET v 2.0.257, 58. Only adult animals were used in the analyses. When data from more than three populations of sparids from our lab were available, we chose truly protandrous and protogynous populations (based on logistic models for size at sex change), i.e., populations with reduced frequencies of non sex-changing individuals (primary males and primary females). Molecular markers were used to calculate N e estimates using the software LDNe59, as implemented in NeEstimator V260 following Gilbert and Whitlock61. Rare alleles with frequencies lower than 0.02 were excluded from our calculations, as suggested by Waples and Do62.

Statistical analyses

Comparisons between protandrous and protogynous species were performed with Student's t-test, Welch t-test, Wilcoxon and Kruskall-Wallis tests depending on the normality and equality of variance of the data. We assessed the effect of mating system (PA and PG) on the effective size ratio N e/ N and N b/ N (to control for abundance), using a series of covariates (after checking for multicollinearity; Fig. S5). For N e/ N the best linear model (Multiple R ² = 0.9396, F 2,5 = 38.92; p = 0.0009; AIC = −18.84) included generation length as a covariate; for N b/ N (Multiple R ² = 0.9446, F 3,4 = 22.73, p = 0.006), the significant covariate was maximum length. A linear model was also run for the ratio N b/ N e (R ² = 0.8929, F 2,5 = 0.85, p = 0.004). Linear models were fitted to annual V k, and N b data considering mating system and sequential sex (first vs second, regardless if male or female) as interaction terms. All analyses and graphical output were performed using the freely available scripts within the R environment54.

Ethical statement

Fish samples were purchased from local fishermen. All methods and protocols were carried out in accordance with relevant guidelines and regulations. Ethical approval AREC-P-10-23 was obtained from UCD, University College Dublin.

Data and materials availability

Microsatellite genotypes can be found on Dryad Digital Repository - doi:10.5061/dryad.n722g4d3 for S. aurata 23; doi:10.5061/dryad.g36ch for C. puniceus 12; doi:10.5061/dryad.sj894 for E. guttatus and C. fulva 9; genotypes for D. sargus, P. erythrinus and L. mormyrus for will be made available upon acceptance.

Electronic supplementary material

Acknowledgements

We would like to thank all researchers that contributed to the empirical studies that generated the genotypes used here for parameter estimation. The manuscript benefited from extensive and constructive comments by Prof Robin Waples. The study was funded by the Irish Research Council and Science Foundation Ireland.

Author Contributions

C.B. and S.M. designed the study, C.B. and J.C. set up the theoretical model, C.B. prepared the life history tables and run the theoretical analysis, C.B. and I.C. analysed the data and produced the figures, S.M. produced the fish drawings for Figure 1, C.B., M.S.-B., I.C. and J.C. collected genetic data, C.B., M.S.-B., I.C. and S.M. commented on study throughout, C.B. wrote the manuscript with the assistance of I.C. and S.M.

Competing Interests

The authors declare that they have no competing interests.

Footnotes

Electronic supplementary material
Supplementary information accompanies this paper at doi:10.1038/s41598-017-09298-8

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

References

10. Buchholz-Sørensen, M. &Vella, A. Population structure, genetic diversity, effective population size, demographic history and regional connectivity patterns of the endangered dusky grouper, Epinephelus marginatus (Teleostei: Serranidae), within Malta's Fisheries Management Zone. PloS one 11, e0159864. doi:0159810.0151371/ journal.pone.0159864 (2016). [DOI] [PMC free article] [PubMed]

14. Shuster, S. M. & Wade, M. J. Mating systems and strategies. Princeton University Press (2003).

24. Chopelet, J. Population Genetic Consequences of Sex Change in Marine Fish. PhD Thesis, University College Dublin (2010).

25. Munday, P. L., Kuwamura, T. & Kroon, F. J. Bidirectional sex change in marine fishes. In: Reproduction and Sexuality in Marine Fishes: Patterns and Processes (ed K.S. Cole) University of California Press, pp. 241-271 (2010).

27. Lowerre-Barbieri, S. et al. Reproductive resilience: a paradigm shift in understanding spawner-recruit systems in exploited marine fish. Fish Fish. doi:10.1111/faf.12180 (2016).

41. Tanasichuk, R. W. Natural mortality rates of adult Pacific herring (Clupea pallasi) from Southern British Columbia. Canadian Stock Assessment Secretariat Research Document 99/168, Nanaimo, B.C., Canada, 34 pp (1999).

49. Froese, R. & Pauly, D. FishBase. World Wide Web electronic publication www.fishbase.org ver. 08/2012 (2011).

Associated Data
This section collects any data citations, data availability statements, or supplementary materials included in this article.
Supplementary Materials

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group

# data[].attributes.attachment_link.title
Ecological and evolutionary consequences of alternative sex-change pathways in fish

# data[].attributes.attachment_link.description
Sequentially hermaphroditic fish change sex from male to female (protandry) or vice versa (protogyny), increasing their fitness by becoming highly fecund females or large dominant males, respectively. These life-history strategies present different ...

# data[].attributes.body
when an individual starts as a male sexually mature producing male gametes but changes into a female sexually mature producing female gametes

this happens in small groups when the largest male would change sex to become the largest female to produce more eggs

reference:

ScienceDirect. (n.d.). Protandry. Retrieved from https://www.sciencedirect.com/topics/biochemistry-genetics-and-molecular-biology/protandry

# data[].attributes.body
Produce offspring by laying eggs

Eggs develop by the fusion of gametes, then the egg hatches

Examples include clownfish and bamboo sharks

# data[].attributes.body
Eggs incubate and hatch inside the mother's body and babies are born live

Egg hatches in the body and develop there too

Examples include sand tiger shark and seahorses

# data[].attributes.body
Young develops inside the mother's placenta are born live

Examples include guppies, mollies and swordtails

# data[].attributes.body
Spawn one and die (single reproductive episode)

Less common

Typically involves a long and exhausting spawning migration from which they cannot recover

Death of parents after single reproductive episode

They overcome this apparent disadvantage by producing more eggs in a single spawning event than closely related species which are iteroparous (high fecundity per reproductive episode)

Survival of young constant/predictable

Examples include sockeye salmon, Pacific lamprey and giant Pacific octopus

# data[].attributes.body
Multiple breeding seasons (annual/monthly/weekly/lunar/daily cycles) (multiple reproductive episodes)

More common

Low fecundity per reproductive episode

Can allow for parental care of eggs and juvenilles (although not all iteroparous species exhibit parental care)

Fish that spawn over multiple annual seasons during their life

Survival of young varied/unpredictable

Examples include trout and cod

# data[].attributes.body
The most common sexual system in fish where each individual is either male or female.

The sex of the fish is determined by genetics (genetic)

Sex can be determined by environmental factors such as water temperature, pH, and oxygen concentration (environmental)

Sex is determined by social factors like the disappearance of a male or female from a group (social)

# data[].attributes.body
Egg Stage

Parrotfish begin as tiny, floating eggs in the water column. These eggs are fertilized externally during spawning events.

They hatch into larvae within a few days, depending on environmental conditions like temperature.

Larval Stage

Newly hatched larvae are planktonic, drifting with ocean currents.

They feed on microscopic plankton and develop rapidly.

This stage lasts a few weeks before they settle on coral reefs.

Juvenile Stage

After settling on a reef, the young parrotfish develop bright colors.

They begin grazing on algae, an essential behavior that helps maintain coral reef health.

Juveniles are often different in color from adults and may resemble different species.

Initial Phase

Most parrotfish start life as females (protogynous hermaphroditism). They will continue growing and reproducing as females.

Some species may remain female their entire lives. (etc. Redband parrotfish)

Terminal Phase

Some dominant females transform into males, displaying dramatic color changes.

These males are larger, more vibrantly colored, and defend territories or harems of females.

This change is influenced by social and environmental factors.

Vedantu. (n.d.). Parrotfish. VEDANTU. https://www.vedantu.com/animal/parrot-fish

Background information. (n.d.). The Parrot Fish. https://parrotfishinformation.weebly.com/background-information.html#:~:text=The%20initial%20phase%20is%20when%20a%20parrot%20fish,They%20also%20will%20experience%20the%20greatest%20reproductive%20success.

# data[].attributes.body
When fish has both male and female reproductive organs and functions

Simultaneous:

Fish have both male and female reproductive organs at the same time

Can alternate between male and female roles during their spawning period

Sequential:

Fish start out as one sex and change to the other later in life

Protandrous - fish develop as male first then change to female

Protogynous - fish develop as female first then change to male

# data[].attributes.body
Structures that transfer sperm during copulation. They can be modified fins, tentacles, or legs.

# data[].attributes.body
A placental analog that helps the embryo absorb nutrients from the mother. There are different types of pseudoplacentas in fish, including trophotaenial placenta, the follicular pseudoplacenta, and the branchial placenta

Allows the developing ambryo to absorb nutrients directly from the mother's ovary in certain species of viviparous fish, essentially providing nourishment to the fetus during gestation

Fishes with pseudoplacenta have a higher likelihood of exhibiting superfetation, which is the ability to gestate multiple broods at different development stages

# data[].attributes.body
A placenta-like structure found in some viviparous fish that allows for the exchange of nutrients and oxygen between the mother and embryo

Allows for the exchange of nutrients and oxygen