# Reproductive Biology of Albacore Tuna (Thunnus alalunga) in the Western Indian Ocean

Zahirah Dhurmeea, Iker Zudaire, Emmanuel Chassot, Maria Cedras, Natacha Nikolic, Jérôme Bourjea, Wendy West, Chandani Appadoo, Nathalie Bodin

► To cite this version:

Zahirah Dhurmeea, Iker Zudaire, Emmanuel Chassot, Maria Cedras, Natacha Nikolic, et al.. Reproductive Biology of Albacore Tuna (Thunnus alalunga) in the Western Indian Ocean. PLoS ONE, 2016, 11 (12), pp.e0168605. 10.1371/journal.pone.0168605 . hal-01926616

HAL Id: hal-01926616

https://hal.umontpellier.fr/hal-01926616

Submitted on 26 May 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

# Reproductive Biology of Albacore Tuna
*(Thunnus alalunga)* in the Western Indian Ocean

Zahirah Dhurmeea¹,²*, Iker Zudaire²,³, Emmanuel Chassot²,⁴, Maria Cedras⁴, Natacha Nikolic⁵, Jérôme Bourjea⁶, Wendy West⁷, Chandani Appadoo⁸, Nathalie Bodin²,⁴

¹ Department of Biosciences, Faculty of Science, University of Mauritius, Réduit, Mauritius, ² IRD—research unit MARine Biodiversity, Exploitation & Conservation (MARBEC), Fishing Port, Victoria, Mahé, Seychelles,

³ IKERBASQUE, Basque Foundation for Science, Bilbao, Spain, ⁴ Seychelles Fishing Authority (SFA), Fishing Port, Victoria, Mahé, Seychelles, ⁵ Ifremer DOI La Réunion, Le Port, Reunion Island, France,

⁶ Ifremer—research unit MARine Biodiversity, Exploitation & Conservation (MARBEC), Sète, France,

⁷ Department of Agriculture Forestry and Fisheries (DAFF), Cape Town, South Africa, ⁸ Department of Marine and Ocean Science, Fisheries and Mariculture, Faculty of Ocean Studies, University of Mauritius, Réduit, Mauritius

* dzahirah@hotmail.com

<figure><img src="image_2.png" /></figure>

CrossMark
click for updates

<figure><img src="image_3.png" /></figure>

OPEN ACCESS

**Citation:** Dhurmeea Z, Zudaire I, Chassot E, Cedras M, Nikolic N, Bourjea J, et al. (2016) Reproductive Biology of Albacore Tuna (*Thunnus alalunga*) in the Western Indian Ocean. PLoS ONE 11(12): e0168605. doi:10.1371/journal.pone.0168605

**Editor:** Aldo Corriero, Universita degli Studi di Bari Aldo Moro, ITALY

**Received:** July 25, 2016

**Accepted:** December 2, 2016

**Published:** December 21, 2016

**Copyright:** © 2016 Dhurmeea et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper and its Supporting Information files. All data are also available through the open electronic source Global Biodiversity Information Facility (GBIF; http://www.gbif.org/).

**Funding:** This work was part of the PhD of the first author ZD co-funded by the "Allocations de Recherche pour une Thèse au Sud" (ARTS) programme of the French Institute of Research for Development (IRD), and the Western Indian Ocean Marine Science Association (WIOMSA). Sampling was conducted under the framework of the European Data Collection Framework (EU DCF

## Abstract

The reproductive biology of albacore tuna, *Thunnus alalunga*, in the western Indian Ocean was examined through analysis of the sex ratio, spawning season, length-at-maturity ($L_{50}$), spawning frequency and fecundity. From 2013 to 2015, a total of 923 female and 867 male albacore were sampled. A bias in sex ratio was found in favor of females with fork length ($L_{F}$) < 100 cm. Using histological analyses and gonadosomatic index, spawning was found to occur between 10°S and 30°S, mainly to the east of Madagascar from October to January. Large females contributed more to reproduction through their longer spawning period compared to small individuals. The $L_{50}$ (mean ± standard error) of female albacore was estimated at 85.3 ± 0.7 cm $L_{F}$. Albacore spawn on average every 2.2 days within the spawning region and spawning months, from November to January. Batch fecundity ranged between 0.26 and 2.09 million oocytes and the relative batch fecundity (mean ± standard deviation) was estimated at 53.4 ± 23.2 oocytes g⁻¹ of somatic-gutted weight. The study provides new information on the reproductive development and classification of albacore in the western Indian Ocean. The reproductive parameters will reduce uncertainty in current stock assessment models which will eventually assist the fishery to be sustainable for future generations.

## Introduction

Albacore tuna, *Thunnus alalunga* (Bonnaterre, 1788), is a temperate species widely distributed in the tropical, sub-tropical and temperate zones from 50°N to 40°S [1]. It is one of the most important species of commercial tuna fisheries worldwide representing 6% of the total catch in 2013 [2]. In the Indian Ocean, this species is caught mostly using drifting longlines followed by purse-seines and other gears [3]. Albacore represents a very important food and economic resource for countries in this region. For example, it is the most harvested species by pelagic
project (Reg. 199/2008 and 665/2008) and the
Fisheries European Funds (EU FEP 2007-2013
project GERMON N°759/DMSOI/2013). Laboratory
analysis was co-funded by Seychelles Fishing
Authority (SFA) and University of Mauritius. The
funders had no role in study design, data collection
and analysis, decision to publish, or preparation of
the manuscript.

**Competing Interests:** The authors have declared
that no competing interests exist.

longliners from Reunion Island after swordfish (*Xiphias gladius*) [4]. Baitboats operating along the southern and western coasts of South Africa and Namibia catch around 4,000 tons (t) of albacore annually [5,6]. Albacore is also very important to the economy of small developing island states such as Mauritius where fishing licenses are issued to some hundred foreign longliners annually, most of which target albacore [7]. Here landings of albacore, which are mostly caught in the exclusive economic zone (EEZ), have increased considerably since 2008, reaching 3,580 t in 2011. Moreover, albacore contributes to the revenue of the local artisanal fishermen who target albacore associated to anchored fish aggregating devices which are set along the coast of the island [7]. Catches of albacore in the Indian Ocean were above the maximum sustainable yield (MSY) in 2013 (43,000 t) but were then recorded to be within MSY limits in 2014 (40,981 t) [3]. The last assessment of the albacore tuna stock in the Indian Ocean indicated that the stock is not overfished and not subject to overfishing despite the existence of uncertainties in the methods of assessment [8].

Reproduction plays a major role in a fish population productivity and, as such, in its resilience to fishing and environmental changes [9,10]. Maturity, fecundity, sex ratio and fish condition are the fundamental factors that affect fish population productivity and are thus important parameters for estimating their reproductive potential [10,11], that is their ability to generate viable eggs in relation to the energy available and parental life expectancy [12]. It has recently been recognized that the spawning stock biomass ($B_s$), long used as a biological indicator in stock assessments, is not suitable to estimate a stock's reproductive potential [13] as it fails to account for the plasticity in reproductive traits [14]. As an alternative to $B_s$, reproductive potential indices have been suggested [10,15] that take into account the basic reproductive parameters to improve the stock assessments.

Albacore migrates over vast distances and develops separate groups at particular stages of its life cycle [1]. A boundary has been described at approximately 30°S in the Indian Ocean with immature fish located south of it, while large albacore occurred in warmer waters in the north [16,17]. It was proposed that the distribution of the three life history stages (immature, spawning only, and mature) of albacore coincided with the three main oceanic currents of the Indian Ocean: the monsoon-driven current in the north of 10°S, the subtropical gyre between 10°S and 30°S and the Circumpolar Current south of 30°S [18]. Spawning of albacore in the western Indian Ocean is thought to occur predominantly east of Madagascar between 15°S and 25°S during the first and fourth quarters of each year when sea surface temperature is above 24°C [19]. Both immature and mature individuals exhibit a north-south seasonal migration although immature tend to remain in southern latitudes, while mature ones remain further north [17].

The reproductive biology of albacore tuna is generally poorly known compared with other principal market tuna species [20]. In the Indian Ocean, contrary to the western Atlantic Ocean [21,22], Pacific Ocean [23,24] and Mediterranean Sea [25-27], biological studies are scarce and little new information is available [28]. In the northeast Atlantic, sexual maturity for males and females is reached between 90 and 94 cm fork length ($L_F$) or at age 5, with no difference noted between the sexes [29]. In the South Pacific, the predicted age at which 50% and 100% of the female population reach sexual maturity was found to be 4.5 years (87 cm $L_F$) [30] and 7 years (94 cm $L_F$), respectively [31].

Similar to other tunas, female albacore have asynchronous oocyte development, with oocytes of different development stages present in the ovary at the same time [22]. Albacore has indeterminate annual fecundity [32], that is, the potential annual fecundity is not fixed before the onset of the spawning season and unyolked oocytes continue to be matured and spawned over a protracted spawning period in several batches [33,34]. They were also reported to spawn daily, on average 1.3 days, during peak spawning months in the South Pacific [24].
Currently, reproductive parameters used in stock assessment models for Indian Ocean albacore originate largely from other stocks [8,35]. Sensitivity analyses of related stock assess- ment models have shown that minor modification of these parameter estimates can pro- foundly influence reference points for which management advice is based on [36]. As a result, given that differences may exist in the population biology of albacore stocks between oceans, using reproductive parameters from other stocks may have a large impact on the stock assess- ment of Indian Ocean albacore. It is therefore vital that the knowledge on these parameters is enhanced for albacore in the Indian Ocean so that advanced scientific advice can also be improved for the sustainable management of albacore in the region.

With special focus on females due to the higher importance of maternal role in offspring production [37], the present study aims to fully describe the reproductive biology of albacore tuna from the western Indian Ocean through the assessment of the fundamental traits that characterize the reproductive potential of a fish species: (i) sex ratio, (ii) spawning season and location, (iii) length at 50% maturity (L₅₀), and (iv) spawning frequency and fecundity. Spatial variability in the first two parameters was also investigated.

**Material and Methods**

**Field sampling**

A total of 923 female and 867 male albacore, caught from June 2013 to June 2015 between lati-
tudes 10°N-40°S and longitudes 10°E-70°E, were sampled from different regions of the west-
ern Indian Ocean (Fig 1; Table 1). Based on the observation of immature, spawning only and
mature albacore in the Indian Ocean from previous studies [16,17], the boundaries of latitudes

<figure><img src="image_4.png" /><figcaption>**Fig 1.** Map of the western Indian Ocean showing positions of sampled albacore tunas by latitudinal area and region, and by gear. Seychelles waters located north of 10'S (region A), Mozambique Channel (region B) and Mauritius and Reunion Island waters (region C) both regions located at latitudes 10-30'S, and South Africa waters located south of 30'S (region D). Black crosses: pole-and-liner (PL), grey triangles: purse-seiner (PS), grey crosses: longliner targeting swordfish (ELL), white triangles: deep-freezing longliner (LL), black dots: coastal longliner (LLCO). Regions are separated by grey dotted lines.</figcaption></figure>

doi:10.1371/journal.pone.0168605.g001
<table><caption>Table 1. Number of female and male albacore tuna sampled by region and latitudinal area in the western Indian Ocean.</caption><thead><tr><th rowspan="2">Latitudinal area</th><th rowspan="2">Region</th><th colspan="2">Females</th><th colspan="2">Males</th></tr><tr><th>n</th><th>L<sub>F</sub></th><th>n</th><th>L<sub>F</sub></th></tr></thead><tbody><tr><td rowspan="1">North of 10°S</td><td rowspan="1">A</td><td rowspan="1">355</td><td rowspan="1">94.9 ± 4.4</td><td rowspan="1">275</td><td rowspan="1">96.8 ± 5.4</td></tr><tr><td rowspan="2">10–30°S</td><td rowspan="1">B</td><td rowspan="1">43</td><td rowspan="1">94.1 ± 5.1</td><td rowspan="1">49</td><td rowspan="1">95.8 ± 5.7</td></tr><tr><td rowspan="1">C</td><td rowspan="1">370</td><td rowspan="1">98.7 ± 3.8</td><td rowspan="1">398</td><td rowspan="1">104.0 ± 4.9</td></tr><tr><td rowspan="1">South of 30°S</td><td rowspan="1">D</td><td rowspan="1">155</td><td rowspan="1">84.9 ± 6.9</td><td rowspan="1">145</td><td rowspan="1">85.3 ± 7.8</td></tr><tr><td rowspan="1">Total</td><td></td><td>923</td><td>94.6 ± 6.7</td><td>867</td><td>98.1 ± 8.7</td></tr></tbody></table>

Region A: Seychelles waters, Region B: Mozambique Channel, Region C: Mauritius and Reunion Island waters, Region D: South Africa waters. n: number of fish, L<sub>F</sub>: mean ± standard deviation in cm.

doi:10.1371/journal.pone.0168605.t001

10°S and 30°S were chosen to assess the monthly variation in the ovarian reproductive phases of albacore tuna in the western Indian Ocean by latitudinal area: latitudes north of 10°S which correspond to the waters of Seychelles (region A), latitudes 10–30°S which group fish caught in the Mozambique Channel (region B) and the waters of Mauritius and Reunion Island (region C), and finally latitudes south of 30°S which include fish from South Africa (region D).

Albacore from regions A and B were caught by commercial European purse-seiners (PS) landing in Seychelles. Albacore from Reunion Island (region C) were sampled at sea onboard a commercial longliner chartered from November to December 2013, and additional samples were obtained from the catch of longliners targeting swordfish (ELL) from April to August 2014. All vessels from Reunion Island were licensed under the French Ministry of Fisheries and Agriculture. Albacore sampled in Mauritius (region C) were caught from October 2013 to March 2015 by professional artisanal fishermen along the coast of the island using vertical longlines (LLCO) and commercial longliners (LL), and were licensed by the Ministry of Ocean Economy, Marine Resources, Fisheries, Shipping and Outer Islands in Mauritius. Albacore from region D were caught from November 2013 to January 2014 and from April to May 2014 by commercial pole-and-line fishing boats (PL) licensed by the Department of Agriculture, Forestry and Fisheries in South Africa, and part of the fish were sampled at sea by observers. All commercial fishing vessels were authorized to operate in the Indian Ocean Tuna Commission (IOTC) area under Resolution 15/04 and within the EEZ of coastal countries of the western Indian Ocean through fishing agreements. Apart from fish sampled at sea, all the other albacore were sampled at processing plants in compliance with their respective health and safety practices (regions A and B: Victoria in Seychelles, region C: Saint-Denis in Reunion Island, Port-Louis in Mauritius; region D: Cape Town in South Africa) or at fish landing sites (fish caught by LLCO). Since all fish sampled under the present study were already dead by the time of sampling, no ethical approval was required.

Different parameters were recorded for each fish: fork length (L<sub>F</sub>, cm; projected straight distance measured with a caliper from the tip of the upper jaw (snout) to the fork of the tail), pectoral length (L<sub>D</sub>, cm; projected straight distance measured with a caliper between the cranial insertion of the pectoral fin and the fork of the tail), total fish weight (W<sub>T</sub>, kg), eviscerated fish weight (ie. W<sub>T</sub> excluding the gut, W<sub>E</sub>, kg), whole gonad weight trimmed of fat (W<sub>G</sub>, g), and sex. Somatic-gutted weight (W<sub>S</sub>, g) was calculated as eviscerated weight (g)—W<sub>G</sub>. An estimated W<sub>S</sub> was used where W<sub>T</sub> or W<sub>E</sub> was unavailable but with a known L<sub>F</sub> or L<sub>D</sub> (i.e., regions A and C), by using length-length and length-weight linear regressions estimated by [38]. A 4–5 cm cross-section from the middle third of either ovary was collected randomly for each female fish when available and preserved in neutral 4% buffered formaldehyde for histological analysis [39].
The gonadosomatic index ($I_G$) [40] was calculated as:

$$I_G = \frac{W_G}{W_S} \times 100$$

## Histological classification of females

A subsample (~0.5 cm) of the ovary cross-section was cut, dehydrated in a graded-series of ethyl-alcohol solutions (70% to absolute), cleared with xylene and infiltrated with paraffin. Paraffin blocks were then sectioned at 6 µm and stained with Harris' hematoxylin and eosin prior to being inspected under a Leica light microscope. Based on the oocyte description and classification criteria by [41] for yellowfin tuna (*Thunnus albacares*), the most advanced oocyte stage in each ovary was identified and used to classify the ovaries into different reproductive phases [42]: the immature phase which is characterized by the presence of oocytes in the primary growth stage (chromatin nucleolus and perinucleolar oocytes), the developing phase which includes oocytes in the cortical alveoli, vitellogenenic 1 and vitellogenic 2 stages, and the spawning capable phase which includes vitellogenic 3, germinal vesicle migration, germinal vesicle breakdown and the hydration stage oocytes (Table 2). When ovulation occurs, the granulosa and theca layers surrounding the hydrated oocyte remain in the ovary tissue and are known as postovulatory follicles (POFs). The presence and the estimated age of POFs were recorded from the histological sections, based on the criteria for yellowfin tuna [43]: 0 hours old (new POFs), <12 hours and >12 hours old.

The relative intensity of atresia, i.e. the percentage of yolked oocytes in alpha atresia relative to the total number of yolked oocytes, was estimated for each female from the histological section [24]. The presence of maturity markers was recorded such as: unovulated or residual hydrated oocytes and late stages of atresia (gamma or delta) which are referred to as brown or yellow bodies [24].

Spawning capable females, defined as females capable of spawning at sampling time or the near future, had <50% alpha atresia (atretic stage 1). Females were classified in the regressing phase, at the end of their reproductive cycle, if their ovaries contained ≥50% (atretic stage 2) or 100% (atretic stage 3) alpha atresia of yolked oocytes or if only beta atresia was present

<table><caption>Table 2. Microscopic classification of the ovary of albacore tuna in the western Indian Ocean.</caption><thead><tr><th>Reproductivephase</th><th>MAGO</th><th>POFs</th><th>Description of atresia in yolked oocytes</th><th>Maturity markers</th><th>Atretic stage</th></tr></thead><tbody><tr><td>Immature</td><td>PG</td><td>Absent</td><td>Absent</td><td>Absent</td><td>-</td></tr><tr><td rowspan="3">Developing</td><td>CA</td><td>Absent</td><td>Absent</td><td>Absent</td><td>-</td></tr><tr><td>Vtg1</td><td>Absent</td><td>Absent or some atresia may be present</td><td>Absent</td><td>-</td></tr><tr><td>Vtg2</td><td>Absent</td><td>Absent or some atresia may be present</td><td>Absent</td><td>-</td></tr><tr><td rowspan="4">Spawning capable</td><td>Vtg3</td><td>May be present</td><td>&lt;50% α atresia. β atresia may be present</td><td>May be present</td><td>1</td></tr><tr><td>GVM</td><td>May be present</td><td>&lt;50% α atresia. β atresia may be present</td><td>May be present</td><td>1</td></tr><tr><td>GVBD</td><td>May be present</td><td>&lt;50% α atresia. β atresia may be present</td><td>May be present</td><td>1</td></tr><tr><td>Hyd</td><td>Absent</td><td>&lt;50% α atresia. β atresia may be present</td><td>May be present</td><td>1</td></tr><tr><td rowspan="3">Regressing</td><td>Vtg3</td><td>Absent</td><td>≥50% α atresia. β atresia may be present</td><td>May be present</td><td>2</td></tr><tr><td>PG, CA</td><td>Absent</td><td>100% α atresia. β atresia may be present</td><td>May be present</td><td>3</td></tr><tr><td>PG, CA</td><td>Absent</td><td>No α atresia. β atresia present</td><td>May be present</td><td>4</td></tr><tr><td>Regenerating</td><td>PG, CA</td><td>Absent</td><td>Absent</td><td>Present</td><td>-</td></tr></tbody></table>

MAGO: most advanced group of oocytes, PG: primary growth (chromatin nucleolus and perinucleolar), CA: cortical alveoli, Vtg1: vitellogenic 1, Vtg2: vitellogenic 2, Vtg3: vitellogenic 3, GVM: germinal vesicle migration, GVBD: germinal vesicle breakdown, Hyd: hydration, POFs: postovulatory follicles, α atresia: alpha atresia, β atresia: beta atresia.

doi:10.1371/journal.pone.0168605.t002
(atretic stage 4), while the ovaries of regenerating females contained maturity markers but no atresia.

## Length-at-maturity

In this study, fish were considered mature when their ovaries reached the spawning capable phase containing vitellogenic 3 as the most advanced oocytes in addition to those classified in the spawning capable, regressing and regenerating phases [24,44]. The length at which 50% of the population of female albacore is mature ($L_{50}$) was derived from fitting the proportion of mature female fish ($P$) to $L_F$ using the following logistic equation [45]:

$$P = \frac{\exp(\alpha + \beta \times L_F)}{1 + \exp(\alpha + \beta \times L_F)}$$

where the parameters $\alpha$ and $\beta$ of the regression model were estimated by assuming a binomial distribution for the logit-transformed data of mature female proportion. Fish were grouped into 3-cm classes of $L_F$ to minimize random fluctuations linked to individual variability. $L_{50}$ was computed as $-\alpha\beta^{-1}$. The variance of the estimate of $L_{50}$ was derived from the delta method using a first-order Taylor approximation [46]. A second $L_{50}$ was additionally computed where females having oocytes in the cortical alveoli, vitellogenic 1 or vitellogenic 2 as most advanced stage in their ovaries were also considered mature (second criterion), as also assumed in previous studies [41,42].

## Fecundity and spawning frequency

Batch fecundity ($F_B$; total number of oocytes released per batch) was estimated using a gravimetric method [47] for ovaries in an advanced stage of maturation, i.e. those having migratory nucleus (germinal vesicle migration and germinal vesicle breakdown) or hydration as most advanced oocytes. Oocytes in the hydration stage were large and translucent. Those in the migratory nucleus were observed as being larger than the other yolked oocytes present and had one or several large oil droplets. A hyaline band could also be observed at their border due to the fusion of yolk as stated by [48]. In each preserved portion of ovary, three subsamples of 0.1 ± 0.01 g were retrieved and saturated with glycerin. The number of migratory nucleus or hydrated oocytes in the weighed subsamples was counted under an Olympus stereomicroscope, and the calculated mean number of oocytes per gram was multiplied by $W_G$ to give an estimate of $F_B$. The relative batch fecundity ($F_{RB}$) was estimated by dividing $F_B$ by the $W_S$ of the fish. Assuming that POFs remain visible in the ovary within 24 hours [24,43], the spawning frequency (average interval between successive spawning events) was estimated as the inverse of the spawning fraction. The spawning fraction was calculated as the total number of females with observable POFs divided by the total number of mature females within our sample [49]. In our study, the spawning frequency was estimated during the spawning period excluding frozen samples from deep-freezing longliners in region C since freezing destroys the follicle wall and hence the POFs would be difficult to detect.

## Data analysis

Sex ratio was calculated as the proportion of females to males in the sample and the differences from the expected 1:1 ratio were analyzed by 5-cm $L_F$ class (to obtain representative number of individuals in each $L_F$ class) using binomial tests. A Chi-square test ($\chi^2$) was used to assess the proportion of spawning capable females in small ($\le L_F$ of the largest immature fish) and large ($> L_F$ of the largest immature fish) size classes during the spawning period. Variations in the
$I_G$ with month for each of the three latitudinal areas were tested using analysis of variance (ANOVA) on the log-transformed values to satisfy the assumption of normality. The relationships between the $F_B$ and $F_{RB}$ with other parameters ($L_F$, $W_G$, $W_S$, $I_G$) were analyzed through linear regressions. All statistical analyses were performed using R version 3.2.2 [50]. All maps were also built in R using the package 'rgdal' [51] and data from http://www.diva-gis.org/Data.

# Results

## Length-frequency and sex ratio

The size of the albacore tuna sampled differed between regions and sexes (Fig 2). Most females from regions A (72%) and B (63%) concentrated in the 90–99 cm $L_F$ class (range: 83–108 cm $L_F$) while female albacore caught in region C mainly peaked (82%) in the 95–104 cm $L_F$ class (range: 88–108 cm $L_F$). The majority of male albacore from regions A (85%) and B (82%).

<figure><img src="image_5.png" /><figcaption>**Fig 2.** Length-frequency distributions of female and male albacore tuna by region in the western Indian Ocean. Region A: Seychelles waters, Region B: Mozambique Channel, Region C: Mauritius and Reunion Island waters, Region D: South Africa waters. $L_F$: fork length in cm. Females are represented in light grey bars and males in black dotted lines.</figcaption></figure>

doi:10.1371/journal.pone.0168605.g002
<table><caption>Table 3. Sex ratio of albacore tuna by region in the western Indian Ocean by 5-cm length class.</caption><thead><tr><th rowspan="2">L<sub>F</sub> class<br>(cm)</th><th colspan="2">Region A</th><th colspan="2">Region B</th><th colspan="2">Region C</th><th colspan="2">Region D</th></tr><tr><th>n</th><th>%M</th><th>n</th><th>%M</th><th>N</th><th>%M</th><th>n</th><th>%M</th></tr></thead><tbody><tr><td>65–69</td><td></td><td></td><td></td><td></td><td></td><td></td><td>1</td><td>100.0</td></tr><tr><td>70–74</td><td></td><td></td><td></td><td></td><td></td><td></td><td>5</td><td>20.0</td></tr><tr><td>75–79</td><td></td><td></td><td></td><td></td><td></td><td></td><td>55</td><td>49.1</td></tr><tr><td>80–84</td><td>5</td><td>0.0</td><td></td><td></td><td></td><td></td><td>98</td><td>50.0</td></tr><tr><td>85–89</td><td>64</td><td>34.8*<sup>2</sup></td><td>16</td><td>43.8</td><td>4</td><td>25.0</td><td>83</td><td>45.8</td></tr><tr><td>90–94</td><td>192</td><td>44.0</td><td>31</td><td>51.6</td><td>56</td><td>26.8*<sup>1</sup></td><td>29</td><td>48.3</td></tr><tr><td>95–99</td><td>232</td><td>35.8*<sup>1</sup></td><td>20</td><td>40.0</td><td>226</td><td>28.3*<sup>1</sup></td><td>7</td><td>42.9</td></tr><tr><td>100–104</td><td>115</td><td>57.4</td><td>22</td><td>72.7</td><td>261</td><td>48.7</td><td>12</td><td>41.7</td></tr><tr><td>105–109</td><td>20</td><td>90.0<sup>§1</sup></td><td>3</td><td>66.7</td><td>167</td><td>86.8<sup>§1</sup></td><td>6</td><td>66.7</td></tr><tr><td>110–114</td><td>2</td><td>100.0</td><td></td><td></td><td>39</td><td>100.0<sup>§1</sup></td><td>2</td><td>50.0</td></tr><tr><td>115–119</td><td></td><td></td><td></td><td></td><td>5</td><td>100.0</td><td>1</td><td>100.0</td></tr></tbody></table>

$L_F$: fork length, %M: percentage of observed males, n: total number of fish sampled per $L_F$ class. Sex ratio analysis was performed through binomial tests.

When the *P* value was significant (1: *P* < 0.001, 2: *P* < 0.05), the sex dominating the $L_F$ class is specified

* = females

§ = males.

Region A: Seychelles waters, Region B: Mozambique Channel, Region C: Mauritius and Reunion Island waters, Region D: South Africa waters.

doi:10.1371/journal.pone.0168605.t003

occupied the 90–104 cm $L_F$ class and had a maximum $L_F$ of 112 cm while in region C, most of the males (69%) dominated the 100–109 cm $L_F$ class and had a maximum $L_F$ of 117 cm. In region D, both males and females reached smaller sizes compared to the other regions (range: 67–118 cm for males and 70–110 cm for females) with around 60% of the fish concentrating in the 80–89 cm $L_F$ class.

Analysis of sex ratio showed that female albacore tend to dominate $L_F$ classes below 100 cm while males dominated in $L_F$ classes larger than 105 cm (Table 3). In region C, males significantly dominated the larger $L_F$ classes (105–109 cm and 110–114 cm), while females were dominant in the 90–94 cm and 95–99 cm $L_F$ classes with the sex ratio significantly different from the expected 1:1 (binomial tests, *P* < 0.05; Table 3). Similar sex ratio results were observed in region A with males significantly dominating 105–109 cm $L_F$ class while females significantly dominated $L_F$ classes of 85–89 cm and 95–99 cm. Males and females for regions A and C appeared to occupy a 1:1 ratio in the 100–104 cm $L_F$ range. In regions B and D, no significant difference in sex ratio was identified in any of the length classes (binomial tests *P* > 0.05; Table 3) although Region B had similar percentage of males across $L_F$ classes as in Region A.

## Size-at-maturity

Most (96%) of the females sampled and analyzed histologically in this study were mature (S1 Table). The largest immature fish (immature phase) was 94 cm $L_F$ while the smallest mature fish (regenerating phase) was 83 cm $L_F$. The proportion of mature females fitted the logistic regression model (Fig 3) and allowed for the estimation of $\alpha$ (mean ± standard error: -40.8749 ± 6.1068) and $\beta$ (mean ± standard error: 0.4791 ± 0.0686). The $L_{50}$ was thus estimated at 85.3 ± 0.7 cm $L_F$ (*P* < 0.0001, *n* = 665) for female albacore tuna of the western Indian Ocean at the vitellogenic 3 maturity threshold. For the second criterion (when fish containing cortical alveoli oocytes onwards were mature), $\alpha$ and $\beta$ were estimated at -40.9305 ± 6.3783 and 0.4820 ± 0.0719 respectively, giving an $L_{50}$ of 84.9 ± 0.8 cm $L_F$ (*P* < 0.0001, *n* = 665).
<figure><img src="image_6.png" /><figcaption>**Fig 3.** Proportion of mature female albacore in the western Indian Ocean for two different criteria based on 3-cm fork length (L<sub>F</sub>) interval fitting. The proportion of females mature when their ovaries contained vitellogenous 3 oocytes onwards (first criterion) are represented by circles and the logistic regression curve fitted to the data is indicated by the black solid line. The proportion of females mature when their ovaries contained cortical alveoli oocytes onwards (second criterion) are represented by grey crosses and the logistic regression curve fitted to the data is indicated by the grey solid line. The horizontal dotted line represents the proportion at which 50% of the population of female albacore tuna is mature and the vertical dotted line shows the corresponding L<sub>F</sub> for each criterion (black and grey for the first and second criteria, respectively).</figcaption></figure>

doi:10.1371/journal.pone.0168605.g003

## Histological classification

Examination of the histological sections revealed that the ovary in the spawning capable phase contains a combination of oocytes in all the different stages of development, confirming an asynchronous ovarian development.

The monthly distribution of fish with ovaries at different reproductive phases was grouped by latitudinal area (north of 10°S, 10–30°S, south of 30°S) in the western Indian Ocean (Fig 4). Most (69%) of the fish caught at latitudes north of 10°S were in the regenerating phase while the rest were in the regressing phase. Around 71% of these regressing females had their ovaries in advanced stages of atresia (atretic stages 3–4). The latitudes 10–30°S were characterized by the occurrence (55%) of spawning capable females (atretic stage 1) which were absent from the higher and lower latitudes. The majority (99%) of these fish were sampled in region C, increasing as from October to November-January (95%), although a small proportion of regressing fish (atretic stage 2) was found during the same period in this region. The abundance of spawning capable fish then decreased from January with very few of them observed in April and May. At the same time, the proportion of both regressing and regenerating females increased. Regressing females in atretic stages 3–4 were observed only from April to July and constituted 76% of all regressing females in this period. The proportion of regenerating
<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>
**Fig 4.** Monthly mean gonadosomatic index and proportion of female albacore with different ovarian developmental phases at different latitudinal areas in the western Indian Ocean. (a) north of 10°S (region A), (b) at latitudes 10–30°S (regions B and C) and (c) south of 30°S (region D) in the western Indian Ocean. *I*<sub>G</sub>: gonadosomatic index, *IP*: immature phase, *DP*: developing phase, *SCP*: spawning capable phase, *RGP*: regressing phase and *RGNP*: regenerating phase. Numbers above the bars indicate the numbers of fish analyzed for *I*<sub>G</sub> and for histological development (in bold) when different. Bars represent standard deviation.

doi:10.1371/journal.pone.0168605.g004

females increased progressively from March to dominate completely in August. Finally, south of 30°S, only immature females (with primary growth as most advanced stage) were present except for few regenerating females that occurred mainly during April and May.

Since spawning capable females were found only in region C, the proportion of large (>94 cm *L*<sub>F</sub>, *n* = 223) to small (≤94 cm *L*<sub>F</sub>, *n* = 26) spawning capable females was assessed for this region only. It was found that there was a significantly larger amount of large females in the spawning capable phase (90%) compared to smaller ones (61%) (χ<sup>2</sup> = 13.7, *P* < 0.001).

Log-transformed *I*<sub>G</sub> varied significantly between months for fish caught north of 10°S (*F*<sub>(6,340)</sub> = 36.1, *P* < 0.0001), at latitudes 10–30°S (*F*<sub>(8,306)</sub> = 141.5, *P* < 0.0001) and south of 30°S (*F*<sub>(3,150)</sub> = 3.3, *P* < 0.05). Mean *I*<sub>G</sub> at the latitudes north of 10°S (region A) was found to remain low (*I*<sub>G</sub> < 0.6) throughout all the months sampled with the highest values noted from December (0.52 ± 0.11) to February (0.53 ± 0.08) and the lowest value in April (0.28 ± 0.09) (Fig 4). At latitudes 10–30°S, the highest *I*<sub>G</sub> were observed from October to January and peaking in November (1.71 ± 0.51) and December (1.92 ± 0.78). These *I*<sub>G</sub> values coincided to the high percentage of spawning capable females during the same period (October-January) and which also peaked in November and December. *I*<sub>G</sub> then decreased to <1.0 in March and continued decreasing until July where the lowest mean value was observed (0.38 ± 0.08) corresponding to the presence of regenerating females. Females caught south of 30°S (region D) had the lowest *I*<sub>G</sub> (*I*<sub>G</sub> < 0.3). The mean *I*<sub>G</sub> for those females was lowest in November and December (both 0.14 ± 0.07) but increased in April (0.19 ± 0.10) and May (0.24 ± 0.14) owing to the presence of a higher percentage of large mature regenerating female albacore in the region.

## Spawning frequency and fecundity

POFs were observed in ovaries of 55 spawning capable females only from region C and from November to January except in ovaries of three fish caught in April and May. These POFs were observed almost exclusively in ovaries containing vitellogenic 3 as most advanced group of oocytes (96%) and the remaining in two ovaries containing migratory nucleus oocyte stage. In addition, the POFs observed in these ovaries were either <12 hours or >12 hours old. Most of the POFs were observed (84%) in spawning capable fish that were landed fresh by artisanal longliners (LLCO and ELL) while a minority in frozen fish from the industrial longliners (LL). This confirms our reasoning of excluding frozen ovaries for the spawning frequency calculations. Out of the 159 mature fish caught during the spawning period (January-December) for the whole region (excluding longliners in region C), a total of 46 females were found to contain POFs which gave a spawning fraction of 0.29 and a mean interval of 3.5 days.

In region C, where the main spawning activity occurs, a total of 101 mature fish were found during the spawning period among which 46 females contained POFs in their ovaries. Thus, in region C, the spawning fraction was 0.46 and the mean interval was 2.2 days.

A total of 54 female fish (range: 87–106 cm *L*<sub>F</sub>) were used for estimating *F*<sub>B</sub> of albacore tuna in the western Indian Ocean. Among those females, 11 had ovaries at germinal vesicle migration, 19 had germinal vesicle breakdown and 24 had hydration as most advanced oocyte stage. These fish were caught at latitude 10–30°S and only from region C except for one fish from region B. *F*<sub>B</sub> estimates for the 54 females varied from 0.26 to 2.09 million eggs
(mean ± standard deviation: 1.07 ± 0.49 million eggs). $F_{RB}$ was estimated to be between 11.1 and 100.5 oocytes g⁻¹ of somatic-gutted weight (mean ± standard deviation: 53.4 ± 23.2 oocytes g⁻¹). Linear regressions showed no significant relationship between $F_B$ and $F_{RB}$ with $W_S$ (Table 4). However, a highly significant positive linear relationship was observed between $W_G$ and $I_G$ with both $F_B$ and $F_{RB}$. $F_B$, but not $F_{RB}$, was significantly related to the size of the fish ($L_F$) (Fig 5). It was also found that $W_G$ was significantly dependent on $L_F$ ($F_{(1,52)} = 10.7$, r² = 0.1549, P < 0.01).

## Discussion

### Sex ratio

Sex ratio analysis showed a predominance of males at larger length classes as also observed in the Pacific and Atlantic Oceans [24,29,32,52,53] and for other tuna species such as bigeye tuna (*Thunnus obesus*), yellowfin tuna, Atlantic (*Thunnus thynnus*) and southern bluefin (*Thunnus maccoyii*) tunas [19]. Two main hypotheses have been put forward to explain this: differences in natural mortality [19,54] as a result of higher energetic demand in females when spawning [55], and differential growth rates [56]. Indeed, in South Pacific albacore, length-at-age was found to vary between sexes with males having higher maximum sizes than females [56]. Also similarly to the South Pacific albacore [24], female albacore in the western Indian Ocean tend to dominate the smaller length classes prior to those dominated by males, which is thought to be mostly the effect of the differences in growth rates [57] and could reflect difference in energetic needs for reproduction between sexes [58]. South of 30°S, in South Africa waters, no significant differences in sex ratio were observed at the same length classes where the majority of fish were small and immature, and not yet involved in reproduction.

### Size-at-maturity

The largest fish classified as immature in our study had an $L_F$ of 94 cm which is the same as estimated for albacore in the South and North Pacific [23,30]. The minimum size-at-maturity found was 83 cm $L_F$ which is slightly smaller than that reported for North and South Pacific albacore (78–90 cm $L_F$; [23,24,32,52] and previous estimates for the Indian Ocean albacore (90 cm $L_F$; [59,60]), but larger than in the Mediterranean Sea (56 cm $L_F$; [26]). The present study is the first to calculate the $L_{50}$ for albacore tuna in the Indian Ocean based on histological analyses, with estimates of the $L_{50}$ at 85.3 ± 0.7 cm based on a vitellogenic 3 maturity threshold. Due to the small number of immature fish found in our analysis and the fact that the albacore were caught by different fishing gears, a spatial analysis of size-at-maturity could not be estimated

<table><caption>Table 4. Summary of the regression analyses testing the relationships between batch fecundity ($F_B$ in million oocytes) and relative batch fecundity ($F_{RB}$ in oocytes g⁻¹ of somatic-gutted weight) with fork length ($L_F$, cm), total gonad weight ($W_G$, g), somatic-gutted weight ($W_S$, g) and gonadosomatic index ($I_G$) for albacore tuna in the western Indian Ocean.</caption><thead><tr><th></th><th>Parameters</th><th>r<sup>2</sup></th><th>P value</th><th>F</th></tr></thead><tbody><tr><td rowspan="3">F<sub>B</sub></td><td>L<sub>F</sub></td><td>0.1089</td><td>&lt;0.05</td><td>6.36</td></tr><tr><td>W<sub>G</sub></td><td>0.3517</td><td>&lt;0.0001</td><td>28.2</td></tr><tr><td>W<sub>S</sub></td><td>0.0588</td><td>0.0774</td><td>3.25</td></tr><tr><td rowspan="4">F<sub>RB</sub></td><td>I<sub>G</sub></td><td>0.3196</td><td>&lt;0.0001</td><td>24.4</td></tr><tr><td>L<sub>F</sub></td><td>0.0225</td><td>0.2789</td><td>1.20</td></tr><tr><td>W<sub>G</sub></td><td>0.2667</td><td>&lt;0.0001</td><td>18.9</td></tr><tr><td>W<sub>S</sub></td><td>0.0009</td><td>0.8326</td><td>0.05</td></tr><tr><td></td><td>I<sub>G</sub></td><td>0.2978</td><td>&lt;0.0001</td><td>22.1</td></tr></tbody></table>

doi:10.1371/journal.pone.0168605.t004
<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /><figcaption>Fig 5. Relationships between (a) fork length (L<sub>F</sub>, cm) and batch fecundity (F<sub>B</sub> in million oocytes) and (b) L<sub>F</sub> and total gonad weight (W<sub>G</sub>, g) for albacore tuna in the western Indian Ocean. n = 54.</figcaption></figure>

doi:10.1371/journal.pone.0168605.g005
such as for the South Pacific albacore [30]. However, our *L*₅₀ is slightly lower than theirs which was estimated at around 87 cm also based on a vitellogenic 3 maturity threshold. To obtain more accurate information on the maturation process, further studies could be designed to account for potential variability in maturity spatially.

Since it was previously suggested that the cortical alveoli stage is a sign of maturation [42], we estimated a new *L*₅₀ (84.9 ± 0.8 cm) based on this assumption. However, at present, there is no standard method that allows for the classification of fish as immature or mature [61] and different studies use different oocyte maturity thresholds [24,41,44]. Thus, information on the amount of individual growth that occurs during the time lag between the cortical alveoli and vitellogenic oocyte stage of a species could improve the accuracy of *L*₅₀ estimates by defining the most appropriate level of gonad development to be used as maturity threshold.

## Spawning behavior

Histological analyses revealed that the latitudes south of 30°S are inhabited mainly by immature fish while the regions north of 30°S are mostly occupied by mature fish. Spawning individuals were exclusively located between 10°S and 30°S mainly to the east of Madagascar (region C), although evidence of spawning activity was also observed on the western side in the Mozambique Channel (region B). Spawning concentrated from October to January with peak spawning observed in November and December. Our results agree with previous studies which suggested that spawning of albacore in the Indian Ocean occurs between 10°S and 30°S [16,17,59] and in the waters off eastern Madagascar [62,63]. Spawning of albacore was also found to occur at similar latitudes and season in the South Pacific Ocean [24,52,64]. However, spawning in the northeast Pacific extends from March to September from 10°N to 20°N [23]. Albacore peak spawning period in the present study was further supported by high *I*<sub>G</sub> (>1.5) which corresponded to fish capable of reproducing [65]. The observation of POFs in ovaries of fish from region C was an additional indication of recent spawning activity [66]. The spawning period corresponded to waters with seawater temperatures above 24°C at the latitudes 10–30°S as previously reported for albacore and other tuna species [17,19,67]. The lower *I*<sub>G</sub> (<1.5) observed in October and January could be related to lower spawning activity at the start and at the end of the spawning season as a result of an asynchrony in spawning at the population level. This can be explained by the gradual increase in the percentage of fish spawning (from October) and falling off (from January) at the end of the spawning season, as observed in the South Pacific albacore [24].

The increase in the proportion of females with regressing ovaries is evidence of the end of spawning season [42] when fish species with indeterminate fecundity, such as albacore, are characterized by an increasing amount of oocytes in atresia [58]. During the spawning season, oocytes are continuously being produced which are recruited into the stock of oocytes in vitellogenesis and this creates an excess of oocytes that are reabsorbed through the process of atresia as the fish reaches the end of the reproductive season [68]. This down-regulation enables fish to adjust the number of oocytes produced relative to their energy reserves [69]. As the fish progresses throughout the spawning season, the incidence of the lower atretic stage (1) decreases while the higher atretic stages (2–4) appear as the vitellogenic oocytes are being reabsorbed gradually during the regressing phase. In the present study, regressing fish with atretic stage 2 were seen mostly during or right after the spawning season while those of atretic stages 3–4 occurred mostly after the spawning season, as from April, showing clearly the evolution of atresia with time. Similar observations were observed in the South Pacific Ocean, where alpha and beta atresia occurred in ovaries of albacore mostly during or right after the spawning season, following which they showed a decreasing trend [30]. Oocyte atresia can then be only
<figure><img src="image_12.png" /><figcaption>**Fig 6.** Schematic map showing the distribution and migration of albacore in the western Indian Ocean during austral summer (left) and winter (right). Positions are shown for spawning capable females (red dots), mature inactive (regressing or regenerating phases) females (green crosses) and immature (immature and developing phases) females (blue circles). Diagonal blue dotted lines indicate main location and feeding grounds of immature (thinner lines) and mature (thicker lines) albacore. Arrows indicate movement of mature inactive fish (green arrows) and immature fish (blue arrows).</figcaption></figure>

doi:10.1371/journal.pone.0168605.g006

observed in their advanced stages (gamma/delta) as yellow bodies, used as a maturity marker in regenerating ovaries through histology in both the Indian Ocean and South Pacific albacore [24,30]. This progression of females from the regressing into the regenerating phase increased gradually in favor of the latter after the spawning season in January, dominating from March to August. In the South Pacific albacore [24], a similar trend in proportion of regenerating females could be detected after the spawning season.

Fig 6 shows the spatial distribution of immature, spawning capable and mature (inactive) albacore life-history stages in austral summer (October to March) and winter period (April to August) and their likely migration. The distribution is in agreement with previous studies in the Indian [16,17] and South Pacific [24] Oceans. The length-frequency data of albacore also show the presence of smaller-sized fish south of 30°S while larger ones are situated northwards. However, histology reveals that a small proportion of mature individuals were found south of 30°S and some immature fish were caught at the spawning latitudes (10°S to 30°S) especially during winter, from April to July. [70] also revealed the presence of few immature fish at the spawning latitudes outside the spawning period; the authors suggested a recent migration of immature fish during austral winter to feed or to prepare for their first spawning in the following season. In region D, all mature fish were mature inactive (in the regenerating phase), also observed by [24] south of 30°S in the Pacific Ocean. Despite the small number of samples analyzed histologically in this region, the abundance of these mature inactive females appeared lower during the spawning season, which could suggest that mature females have migrated northward of 30°S to the spawning latitudes to reproduce.

North of 10°S, the percentage of females with the different ovarian phases from August to November remains unknown as no fish from this period and area could be sampled. However, we found that most of the fish had ovaries in the regressing phase, even during the spawning months of December and January, suggesting that large-scale spawning does not occur in this region. The presence of mature inactive (regressing and regenerating) females north of 10°S in
Seychelles waters and south of 30°S in South African waters, suggests that after the completion of the spawning season, some females migrate out of the spawning latitudes towards the north or the south of the Indian Ocean. This suggested migration is consistent with the decreased catch rates of albacore observed after the spawning season in Mauritius [7] and Reunion Island [4], and the increased catch rates around south of 30°S from April to September [71]. It is likely that the regions located north of 10°S, where nutrient-rich waters occur as a result of upwelling [72] could be used as a main feeding ground by mature albacore to replenish their energy reserves after reproduction. The cold and rich waters south to Madagascar and off South Africa could also be a potential feeding zone for albacore as previously thought for swordfish [73]. This could especially be true for the immature which seem to occur mostly in this area. However, the low proportion of inactive but mature female albacore caught throughout the year at latitudes 10-30°S, in addition to fishery statistics from longliners indicating that albacore are caught in this area throughout the year [8], suggest that part of the females also remain at these latitudes even after the spawning season.

The occurrence of regressing or regenerating females during the spawning period in region C suggests that not all fish spawn during the whole period which is consistent with the asynchrony in spawning activity of the population. This asynchrony could be the result of the difference in spawning period between small and large fish. Previous studies have suggested that larger and older female albacore have a longer spawning period [24]. This difference has also been observed in yellowfin tuna [67], skipjack tuna [74], and swordfish [75].

In this study, a higher proportion of smaller-sized females already had their ovaries in the regressing and regenerating phases compared to larger females (which were mostly in the spawning capable phase) during the same period. These results are consistent with the theory that smaller and younger fish that are first maturing will tend to have a lower spawning activity than larger and older fish [34]. This may in turn be due to the energy equilibrium between somatic growth and reproduction which is dependent on fish size [76]. A shorter spawning period for the small and young fish has also been observed in the Pacific mackerel (*Scomber japonicus*) [77] and jack mackerel (*Trachus symmetricus*) [78]. It was also proposed that the adjustment of small and young mature albacore to the tropical waters could result in a lower spawning fraction [24].

## Spawning frequency and fecundity

The mean interval between two successive spawning events for albacore in the entire western Indian Ocean (3.4) and specifically for region C (2.2) during the spawning period (from November to January) was lower than for the South Pacific albacore during its peak spawning months (1.3 days from October to December) from 10-25°S [24]. The values are also lower than that reported for albacore from the North Pacific Ocean during peak spawning activity (1.7 days) [23] although lower values of 3.0 days have been observed from March to April [53]. The Indian Ocean albacore is characterized by a lower spawning frequency compared to skipjack (1.2 days) [49], yellowfin (1.5 days) [55], bigeye (1.1 days) [79] and southern bluefin tunas (1.6 days) [80]. In albacore, POFs have previously been observed in ovaries with oocytes at maturation and hydration stages [23,24]. Usually, the presence of POFs and oocytes in advanced stages of maturation in an ovary implies that the fish is capable of daily spawning [19]. In the present study, the majority of POFs were observed in ovaries of fish with vitellogenic 3 as most advanced oocytes which is different to that for fish from the North and South Pacific Ocean where POFs were observed in ovaries with migratory nucleus and hydrated oocytes [23,24]. The low spawning frequency estimate for Indian Ocean albacore suggests that their oocyte maturation rate could be lower than those of albacore from the South Pacific. This could in turn be
due to lower energy available for reproduction in the region since spawning intervals have been linked to food availability [81]. Further investigations are required to account for the differences in spawning frequencies taking into account possible effects of fishing gear [23].

The mean $F_{RB}$ of 53.4 ± 23.2 oocytes per gram of body weight estimated in this study is in the same range as for the North Pacific albacore (50.5 ± 22.8; [23]), southern bluefin tuna (56.5 ± 16.1; [80]), South Pacific albacore (64.4 ± 24.7; [24]) and the Atlantic bluefin (74.7 ± 45.3; [82]). It is also consistent with $F_{RB}$ estimated for other tunas such as the western Pacific bigeye (16–150; [83]), yellowfin from the western Indian Ocean (9–180; [67]), the eastern Pacific (68.0 ± 20.7; [43]), the western Pacific Ocean (22.2–77.6; [84]) and Hawaii area (31.9–147.1; [84]), but was lower than the Indian Ocean skipjack tuna (140 ± 64; [74]). In the present study, $F_B$ for albacore was 0.26 to 2.09 million oocytes, which is higher than the North Pacific albacore (0.17 to 1.66 million eggs; [23]) but is similar to that from the South Pacific albacore although the latter could reach a higher maximum value (0.26 to 2.83 million oocytes; [24]).

It is known that fecundity varies within individuals, being proportional to fish size [34]. It has been suggested that larger and older fish will tend to produce larger numbers of eggs that will hatch into a larger proportion of viable larvae [85]. In our study, $F_B$ showed a significant relationship with fish size but not with body weight, as is also the case with swordfish [86], and $F_{RB}$ is also not significantly related to $L_F$ as observed in the Atlantic bluefin tuna [82]. This suggests that the increased $F_B$ in larger albacore is mainly due to their larger ovary as gonad growth is proportionally related to fish size [82]. This is consistent with the result found in the present study which shows that the larger fish developed significantly larger gonads. However, the observed relationship between $F_B$ and fish size was variable with $F_B$ of fish of the same size varying to a large extent. This has also been shown in the South Pacific albacore [24] and other tunas [55,79,80] and has been attributed to the period in the spawning season that the fish was caught [24,80]. This is in turn due to the variation in the number of oocytes that are recruited and spawned as the spawning season continues [87]. $F_B$ has been observed to decrease over time for South Pacific albacore, with higher batches occurring at the beginning than at the end of the spawning season [24]. Fish condition is also a factor that can affect fecundity [34] with better condition linked with higher fecundity [88], earlier maturation and higher energy allocation to egg production [89]. However, reduced oocyte growth has been observed in females with poor feeding conditions [90]. Exploration on the performance of albacore offspring in relation to maternal age and condition could provide additional information that may possibly be integrated in the analysis of the albacore population.

## Conclusions

This study describes numerous aspects of the reproductive biology of albacore in the western Indian Ocean and provides a number of parameters for improving future stock assessments of this species in the region. This will, in turn, enhance the confidence of population estimates of Indian Ocean albacore. Most importantly, stock assessments can be improved by taking into account the bias in sex ratio towards males in the larger size classes, and the importance of large females, and their greater contribution to reproduction, and thus recruitment, through their longer spawning period compared to small individuals. In the long-term, the results could be useful to generate advice and development of management strategies to maintain the MSY of the Indian Ocean albacore stock.

## Supporting Information

S1 Table. Number of fish per oocyte developmental stage or ovarian phase by 3-cm interval fork length ($L_F$) used for maturity analysis of female albacore tuna in the western Indian
**Ocean based on histology.** In the first criterion, immature fish were viewed as containing the most advanced oocyte stage the primary growth (PG) stage, cortical alveolar (CA) and vitellogenic 1 (Vtg1). Mature fish were those containing secondary vitellogenic 2 (Vtg2), vitellogenic 3 (Vtg3), germinal vesicle migration (GVM), germinal vesicle breakdown (GVBD) and hydrated (Hyd) oocytes as most advanced stage, as well as those in regressing (RGP) and regenerating phases (RGNP). In the second criterion, the maturity threshold was set at CA oocyte stage.

(PDF)

**Acknowledgments**

This work was conducted under the framework of the EU FEP 2007–2013 project GERMON (N°759/DMSOI/2013) and the EU DCF project (Reg. 199/2008 and 665/2008). It is part of the PhD of the first author, co-funded by the University of Mauritius, the “Allocations de Recherche pour une Thèse au Sud” (ARTS) programme of the French Institute of Research for Development (IRD), and the Western Indian Ocean Marine Science Association (WIOMSA). The authors would like to thank various people for their help during sampling: J. Esparon, A. Tirant, A. Stephen, R. Rose, E. Mathiot, M. Elisabeth, P. Boniface, M. Lesperance, and S. Hollanda from SFA; M. Meyer, J. Cunningham and K. Baloyi from DAFF and C. Heinecken from Capfish; E. Romanov, E. Richard and L. Le Foulgoc from CAP RUN—HYDRORUN, and finally H. Evano, A. Puech and D. Roos from Ifremer DOI La Réunion. We also thank fishermen, skippers and other members of the fishing industry including the workers from the tuna processing factories and canneries in Seychelles (Indian Ocean Tuna Ltd), Mauritius (Thon des Mascareignes Ltd and Pelagic Process Ltd), South Africa (Indian Ocean Tuna Ltd) and Reunion Island (Reunipeche and Manohal team) for their cooperation and assistance during sampling activities. We are grateful for various discussions with J. Farley and H. Pethybridge at CSIRO that helped advance the manuscript. Finally, appreciation goes to K. Constance, C. Belmont, and S. Lawrence from SFA for their help in laboratory preparation of samples, and for the two anonymous reviewers and the editor for their careful readings and precious recommendations that helped improving the manuscript.

**Author Contributions**

**Conceptualization:** ZD NB.

**Data curation:** NB EC ZD IZ.

**Formal analysis:** ZD EC.

**Funding acquisition:** NB.

**Investigation:** ZD NB IZ MC WW NN.

**Methodology:** ZD NB.

**Project administration:** NB.

**Resources:** NB MC CA JB NN WW ZD.

**Supervision:** NB CA.

**Writing – original draft:** ZD.

**Writing – review & editing:** ZD NB IZ EC NN JB.
## References

1. Collette B, Nauen C. FAO species catalogue, Vol 2. Scombrids of the world: an annotated and illustrated catalogue of tunas, mackerels, bonitos and related species known to date. FAO Fish Synop 1983; No. 125(2),137 p.

2. ISSF. ISSF Tuna Stock Status Update, 2015: Status of the world fisheries for tuna. ISSF Technical Report 2015-03A. 2015. Available: http://responsiblequality.riomare.com/download/ISSF-2015-03A-Status-of-the-Stocks-Nov-2015.pdf. Accessed 4 March 2016.

3. IOTC. Report of the Eighteenth Session of the IOTC Scientific Committee. 2015. IOTC-2015-SC18-R [E]: 175. Available: http://www.iotc.org/meetings/18th-scientific-committee-sc18. Accessed 9 July 2016.

4. Coelho R, Nikolic N, Evano H, Miguel N, Bourjea J. Reunion island pelagic longline fishery characterization and standardization of albacore catch rates. Indian Ocean Tuna Commission. 2014. Available: http://archimer.ifremer.fr/doc/00208/31926/. Accessed 5 July 2016.

5. Newcombe HC. Contribution towards the development of a management plan for the baitboat and sport fishery for tuna in South Africa. M.Sc. Thesis, Rhodes University. 2011.

6. Kerwath SE, Winker H, West WM. Standardization of the catch per unit effort for albacore (*Thunnus alalunga*) for the South African tuna-pole (baitboat) fleet for the time series 1999–2010. ICCAT Coll Vol Sci Pap. 2012; 68: 604–614.

7. Dhurmeea Z, Beeharry SP, Sooklall T. Catch/effort and length-frequency data collected on albacore tuna landed in Mauritius. Working Paper IOTC-2012-WPTmT04-12, Indian Ocean Tuna Commission. Working Party on Temperate Tunas. 2012. Available: http://www.iotc.org/documents/catch-effort-and-length-frequency-data-collected-albacore-tuna-landed-mauritius. Accessed 10 October 2012.

8. IOTC. Report of the Fifth Session of the IOTC Working Party on Temperate Tunas. IOTC-2014-WPTmT05-R[E]:53. 2014. Available: http://www.iotc.org/documents/report-fifth-session-iotc-working-party-temperate-tunas. Accessed 9 January 2016.

9. Murua H, Motos L. Reproductive strategy and spawning activity of the European hake *Merluccius merluccius* (L.) in the Bay of Biscay. J Fish Biol. 2006; 69: 1288–1303.

10. Morgan MJ, Murua H, Kraus G, Lambert Y, Marteinsdottir G, Marshall C et al. The evaluation of reference points and stock productivity in the context of alternative indices of stock reproductive potential. Can J Fish Aquat Sci. 2009; 66: 404–414.

11. Lambert Y. Why should we closely monitor fecundity in marine fish populations? J Northwest Atl Fish Sci. 2008; 41: 93–106.

12. Wootton RJ. Introduction: tactics and strategies in fish reproduction. In: Potts GW and Wootton RJ, editors. Fish reproduction: strategies and tactics. London: Academic Press; 1984. pp. 1–12.

13. Marshall TC. Implementing information on stock reproductive potential in fisheries management: the motivation, challenges and opportunities. In: Jakobsen MF, Megrey B and Moksness E, editors. Fish reproductive biology: implications for assessment and management. United Kingdom: Wiley-Blackwell; 2009. pp. 395–420.

14. Morgan M. Integrating reproductive biology into scientific advice for fisheries management. J Northwest Atl Fish Sci. 2008; 41: 37–51.

15. Murua H, Ibaibarriaga L, Álvarez P, Santos M, Korta M, Santurtun M et al. The daily egg production method: A valid tool for application to European hake in the Bay of Biscay? Fish Res. 2010; 104: 100–110.

16. Suda A. Recent status of resources of tuna exploited by longline fishery in the Indian Ocean. Bull Far Seas Fish Res Lab. 1974; 10: 27–62.

17. Chen IC, Lee PF, Tzeng WN. Distribution of albacore (*Thunnus alalunga*) in the Indian Ocean and its relation to environmental factors. Fish Oceanogr. 2005; 14: 71–80.

18. Wyrtki K. Physical oceanography of the Indian Ocean. In: Zeitzschel B, editor. The Biology of the Indian Ocean. Berlin, New York: Spring-Verlag; 1973. pp. 18–36.

19. Schaefer KM. Reproductive biology of tunas. In: Block BA and Stevens ED, editors. Fish Physiology, Vol. 19. San Diego, CA: Academic Press; 2001. pp. 225–270.

20. Juan-Jordá MJ, Mosqueira I, Freire J, Dulvy NK. The conservation and management of tunas and their relatives: setting life history research priorities. PLOS ONE. 2013; 8: e70405. doi: 10.1371/journal.pone.0070405 PMID: 23950930

21. Le Gall JY. Exposé synoptique des données biologiques sur le germon *Thunnus alalunga* (Bonaterre, 1788) de l'Océan Atlantique. FAO Fish Synop. 1974; 109: 1–70.
22. Ortiz de Zárate V, Macías D, Satoh K, Saito H. Information on the reproduction of albacore (*Thunnus alalunga*) in the central and tropical North Atlantic in 2002. ICCAT Coll Vol Sci Pap. 2004; 56: 1450–1462.

23. Chen KS, Crone PR, Hsu CC. Reproductive biology of albacore *Thunnus alalunga*. J Fish Biol. 2010; 77: 119–136. doi: 10.1111/j.1095-8649.2010.02662.x PMID: 20646142

24. Farley JH, Williams AJ, Hoyle SD, Davies CR, Nicol SJ. Reproductive dynamics and potential annual fecundity of South Pacific albacore tuna (*Thunnus alalunga*). PLOS ONE. 2013; 8: e60577. doi: 10.1371/journal.pone.0060577 PMID: 23565258

25. Saber S, Vive MG, Báez JC, Godoy D, Macías D. Estimating batch fecundity of the Mediterranean albacore. A comparison between quantification of postvitellogenic follicles and postovulatory follicles. ICCAT Coll Vol Sci Pap. 2011; 66: 1919–1930.

26. Saber S, de Urbina JO, Gómez-Vives MJ, Macías D. Some aspects of the reproductive biology of alba-core *Thunnus alalunga* from the western Mediterranean Sea. J Mar Biol. Assoc UK. 2015; 95: 1705–1715.

27. Saber S, Macías D, de Urbina JO, Kjesbu OS. Contrasting batch fecundity estimates of albacore (*Thunnus alalunga*), an indeterminate spawner, by different laboratory techniques. Fish Res. 2016; 176: 76–85.

28. IOTC. Report of the Sixteenth Session of the IOTC Scientific Committee. IOTC-2013-SC16-R[E]:312. 2013. Available: http://www.iotc.org/meetings/16th-session-scientific-committee. Accessed 15 March 2016.

29. Bard FX. Le thon germon (Thunnus alalunga Bonnaterre 1788) de l'Océan Atlantique. De la dynamique des populations à la stratégie démographique. Ph.D. Thesis. Université Paris VI. 1981.

30. Farley JH, Hoyle SD, Eveson JP, Williams AJ, Davies CR, Nicol SJ. Maturity ogives for South Pacific albacore tuna (*Thunnus alalunga*) that account for spatial and seasonal variation in the distributions of mature and immature fish. PLOS ONE. 2014; 9: e83017 doi: 10.1371/journal.pone.0083017 PMID: 24416153

31. Farley JH, Williams AJ, Davies CR, Clear NP, Eveson JP, Hoyle SD et al. Population biology of alba-core tuna in the Australian region. Final report for FRDC project 2009/012. Hobart: CSIRO. 124 p; 2012. Available: http://frdc.com.au/research/Final_Reports/2009-012-DLD.pdf. Accessed 9 November 2015.

32. Otsu T, Uchida RN. Sexual maturity and spawning of albacore in the Pacific Ocean. Fish Bull. 1959; 148: 287–305.

33. Hunter JR, Macewicz BJ, Lo NH, Kimbrell CA. Fecundity, spawning and maturity of females Dover Sole, *Microstomus pacificus*, with an evaluation of assumptions and precision. Fish Bull. 1992; 90: 101–128.

34. Murua H, Kraus G, Saborido-Rey F, Witthames PR, Thorsen A, Junquera S. Procedures to estimate fecundity of marine fish species in relation to their reproductive strategy. J Northwest Atl Fish Sci. 2003; 33: 33–54.

35. IOTC. Report of the Fifteenth Session of the IOTC Scientific Committee. IOTC-2012-SC15-R[E]:288. 2012. Available: http://www.iotc.org/meetings/15th-session-scientific-committee. Accessed 20 January 2014.

36. Hoyle S. Adjusted biological parameters and spawning biomass calculations for south Pacific albacore tuna, and their implications for stock assessments. Working paper WCPFC-SC4-2008/ME-WP-2, 4th Scientific Committee meeting of the Western and Central Pacific Fisheries Commission. 2008. Available: www.wcpfc.int/meetings/2008/4th-regular-session-scientificcommittee. Accessed 22 December 2013.

37. Helfman GS, Collette BB, Facey DE. The diversity of fishes. Massachussettes, USA: Blackwell Science; 1997. 529 p.

38. Dhurmeea Z, Chassot E, Augustin E, Assan C, Nikolic N, Bourjea J et al. 2016. Morphometrics of alba-core tuna (Thunnus alalunga) in the Western Indian Ocean. Working Paper IOTC-2016-WPTmT06-28, Indian Ocean Tuna Commission. Working Party on Temperate Tunas. 2016. Available: http://www.iotc.org/documents/morphometrics-albacore-tuna-western-indian-ocean. Accessed 18 July 2016.

39. Hunter JR. Preservation of northern anchovy in formaldehyde solution: An egg production method for estimating spawning biomass of pelagic fish: application to the northern anchovy, *Engraulis mordax*. NOAA Tech Rep NMFS. 1985; 36: 63–65.

40. Schaefer MB, Orange CJ. Studies of the sexual development and spawning of yellowfin tuna (*Neothunnus macroptus*) and skipjack (*Katsuwonus pelamis*) in three areas of the Eastern Pacific Ocean, by examination of gonads. Inter-Am Trop Tuna Comm Bull. 1956; 1: 283–349.

41. Zudaire I, Murua H, Grande M, Korta M, Arrizabalaga H, Areso JJ et al. Fecundity regulation strategy of the yellowfin tuna (*Thunnus albacares*) in the Western Indian Ocean. Fish Res. 2013; 138: 80–88.
42. Brown-Peterson NJ, Wyanski DM, Saborido-Rey F, Macewicz BJ, Lowerre-Barbieri SK. A standardized terminology for describing reproductive development in fishes. Mar Coast Fish: Dyn Manag Ecosyst Sci. 2011; 3: 52–70.

43. Schaefer K. Spawning time, frequency, and batch fecundity of yellowfin tuna, *Thunnus albacares*, near Clipperton Atoll in the eastern Pacific Ocean. Fish Bull. 1996; 94: 98–113.

44. Okochi Y, Abe O, Tanaka S, Ishihara Y, Shimizu A. Reproductive biology of female Pacific bluefin tuna, *Thunnus orientalis*, in the Sea of Japan. Fish Res. 2016; 174: 30–39.

45. Saborido-Rey F, Junquera S. Histological assessment of variations in sexual maturity of cod (*Gadus morhua* L.) at the Flemish Cap (north-west Atlantic). J Mar Sci. 1998; 55: 515–521.

46. Xu J, Long JS. Confidence intervals for predicted outcomes in regression models for categorical outcomes. Stata Journal. 2005; 5(4): 537–559.

47. Hunter JR, Macewicz BJ, Kimbrell CA. Fecundity and other aspects of the reproduction Sablefish, *Anoplopoma fimbria*, in central California water. CalCOFI Rep. 1989; 30: 1–72.

48. Karlou-Riga C, Economidis PS. Spawning frequency and batch fecundity of horse mackerel, *Trachurus trachurus* (L.), in the Saronikos Gulf (Greece). J Appl Ichthyol. 1997; 13: 97–104.

49. Hunter JR, Macewicz BJ, Sibert JR. The spawning frequency of skipjack tuna, *Katsuwonus pelamis*, from the South Pacific. Fish Bull. 1986; 84: 895–903.

50. R Development Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing. Vienna, Austria. 2015 Available: https://www.R-project.org/.

51. Bivanand R, Keitt T, Rowlingson B, Pebesma E, Summer M, Hijmans R et al. Package 'rgdal'. 2016.

52. Ramón D, Bailey K. Spawning seasonality of albacore, *Thunnus alalunga*, in the South Pacific Ocean. Fish Bull. 1996; 94: 724–733.

53. Hsu CC, Chen KS. Reproductive Biology of Albacore (Thunnus alalunga) from Central and Western North Pacific Ocean. (ISC-ALBWG/05/09). Albacore Working Group Meeting, 28 November-2 December 2005. La Jolla, USA: International Scientific Committee.

54. Maunder MN, Watters GM. A-SCALA: an age-structured statistical catch-at-length analysis for assessing tuna stocks in the eastern tropical Pacific Ocean. Inter-Am Trop Tuna Comm Bull. 2003; 22: 435–582.

55. Schaefer KM. Reproductive biology of yellowfin tuna (*Thunnus albacares*) in the eastern Pacific Ocean. Inter-Am Trop Tuna Comm Bull. 1998; 21: 201–272.

56. Williams AJ, Farley JH, Hoyle SD, Davies CR, Nicol SJ. Spatial and sex-specific variation in growth of albacore tuna (*Thunnus alalunga*) across the South Pacific Ocean. PLOS ONE. 2012; 7: e39318. doi: 10.1371/journal.pone.0039318 PMID: 22723993

57. Harley S, Maunder MA. Simple model for age-structured natural mortality based on changes in sex ratios. Inter-American Tropical Tuna Commission. 2003. Available: https://www.iattc.org/PDFFiles/Age%20structured%20natural%20mortality.pdf. Accessed 9 July 2016.

58. Murua H, Saborido-Rey F. Female reproductive strategies of marine fish species of the North Atlantic. J Northwest Atl Fish Sci. 2003; 33: 23–31.

59. Ueyanagi S. Observations on the distribution of tuna larva in the Indo-Pacific Ocean with emphasis on the delineation of spawning areas of albacore, *Thunnus alalunga*. Bull Fish Res Lab. 1969; 2: 177–219.

60. Wu CL, Kuo CL. Maturity and fecundity of albacore, *Thunnus alalunga* (Bonnaterre), from the Indian Ocean. J Fish Soc Taiwan. 1993; 20: 135–151.

61. Lowerre-Barbieri SK, Ganias K, Saborido-Rey F, Murua H, Hunter JR. Reproductive timing in marine fishes: Variability, temporal scales, and methods. Mar Coast Fish. 2011; 3: 71–91.

62. Conand F, Richards WJ. Distribution of tuna larvae between Madagascar and the Equator, Indian Ocean. Biol Oceanogr. 1982; 1: 31–336.

63. Shiohama T. Overall fishing intensity and length composition of albacore caught by long line fishery in the Indian Ocean, 1952–1982. IPTP TWS/85/22. 1985; 91–109.

64. Otsu T, Hansen RJ. Sexual maturity and spawning of albacore in the Central South Pacific Ocean. Fish Bull. 1962; 62: 151–161.

65. Nootmorn P, Yakoh A, Kawises K. Reproductive biology of yellowfin tuna in the Eastern Indian Ocean. Working paper IOTC-2005-WPTT-14. Working party on tropical tunas. Indian Ocean Tuna Commission. 2005. Available: https://iotc.org/sites/default/files/documents/proceedings/2005/wptt/IOTC-2005-WPTT-14.pdf. Accessed 5 November 2015.

66. Fitzhugh GR, Hettler WF. Temperature influence on postovulatory follicle degeneration in Atlantic menhaden, *Brevoortia tyrannus*. Fish Bull. 1995; 93: 568–568.
67. Zudaire I, Murua H, Grande M, Bodin N. Reproductive potential of yellowfin Tuna (*Thunnus albacares*) in the western Indian Ocean. Fish Bull. 2013; 111: 252–263.

68. Guraya SS. Gonadal development and production of gametes in fish. Proc. Indian natn Sci Acad, B60. 1994; 60: 15–32.

69. Kennedy J, Witthames PR, Nash RDM, Fox CJ. Is fecundity in plaice (*Pleuronectes platessa L.*) down-regulated in response to reduced food intake during autumn? J Fish Biol. 2008; 72: 78–92.

70. Farley J, Clear N. Preliminary study of age, growth, and spawning activity of albacore in Australia's Eastern Tuna and Billfish Fishery (WCPFC-SC4-2008/BI-IP-1). Western and Central Pacific Fisheries Commission. 2008. Available: http://www.spc.int/DigitalLibrary/Doc/FAME/Meetings/WCPFC/SC4/BI-IP-1.html. Accessed 20 March 2015.

71. Hoyle SD, Sharma R, Herrera M. Stock assessment of albacore tuna in the Indian Ocean for 2014 using Stock Synthesis. Working paper IOTC-2014-WPTmT05-24, 5th Working Party on temperate tunas. Indian Ocean Tuna Commission. 2014. Available: http://www.iotc.org/documents/indian-ocean-albacore-assessment. Accessed 10 July 2016.

72. Schott FA, McCreary JP Jr. The monsoon circulation of the Indian Ocean. Prog Oceanogr. 2001; 51: 1–123.

73. Poisson F, Fauvel C. Reproductive dynamics of swordfish (*Xiphias gladius*) in the southwestern Indian Ocean (Reunion Island). Part 1: oocyte development, sexual maturity and spawning. Aquat Living Resour. 2009; 22: 45–58.

74. Grande M, Murua H, Zudaire I, Goñi N, Bodin N. Reproductive timing and reproductive capacity of the skipjack tuna (*Katsuwonus pelamis*) in the western Indian Ocean. Fish Res. 2014; 156: 14–22.

75. Poisson F, Fauvel C. Reproductive dynamics of swordfish (*Xiphias gladius*) in the southwestern Indian Ocean (Reunion Island). Part 2: fecundity and spawning pattern. Aquat Living Resour. 2009; 22: 59–68.

76. Claramunt G, Serra R, Castro LR, Cubillos L. Is the spawning frequency dependent on female size? Empirical evidence in *Sardinops sagax* and *Engraulis ringens* off northern Chile. Fish Res. 2007; 85: 248–257.

77. Knaggs EH, Parrish RH. Maturation and growth of Pacific mackerel Scomber japonicus Houttuyn. Long Beach: California Department of Fish and Game, Mar Res Tech Rep 3. 1973. 19 p.

78. MacCall AD, Frey HW, Huppert DD, Knaggs EH, McMillan JA, Stauffer GD. Biology and economics of the fishery for jack mackerel in the northeastern Pacific. NOAA Tech Rep. 1980. NMFS-SWFSC-4.

79. Nikaido H, Miyabe N, Ueyanagi S. Spawning time and frequency of bigeye tuna, *Thunnus obesus*. Bull Nat Res Inst Far Seas Fish. 1991; 28: 47–73.

80. Farley JH, Davis TL. Reproductive dynamics of southern bluefin tuna, *Thunnus maccoyii*. Fish Bull. 1998; 96: 223–236.

81. Tsuruta Y. Reproduction in the Japanese anchovy (*Engraulis japonica*) as related to population fluctua-tion. Bull Nat Res Inst Far Seas Fish. 1992; 13: 129–168.

82. Aranda G, Medina A, Santos A, Abascal FJ, Galaz T. Evaluation of Atlantic bluefin tuna reproductive potential in the western Mediterranean Sea. J Sea Res. 2013; 76: 154–160.

83. Sun CL, Chu SL, Yeh S. The reproductive biology of female bigeye tuna (Thunnus obesus) in the Western Pacific. Working paper WCPFC-SC2-2006/BI WP-1, 2nd Scientific Committee meeting, Western and Central Pacific Fisheries Commission. 2006. Available: www.wcpfc.int/meetings/2006/2nd-regularsession. Accessed 3 November 2015.

84. Itano D. The reproductive biology of yellowfin tuna (Thunnus albacares) in Hawaiian waters and the western tropical Pacific Ocean: Project summary. SOEST 00-01, JIMAR Contribution 00-328. 2000. Available: http://www.soest.hawaii.edu/PFRP/biology/itano/itano_yft.pdf. Accessed 9 July 2016.

85. Solemdal P. Maternal effects—a link between the past and the future. J Sea Res. 1997; 37: 213–227.

86. Arocha F, Lee DW. Maturity at size, reproductive seasonality, spawning frequency, fecundity and sex ratio in the swordfish from the Northwest Atlantic. ICCAT Coll Vol Sci Pap. 1996; 45: 350–357.

87. Arocha F, Lee DW, Grubich JR. Observations on sex ratio, maturity stages, and fecundity estimates of swordfish, *Xiphias gladius*, in the Northwest Atlantic Ocean. ICCAT Coll Vol Sci Pap. 1994; 42: 309–318.

88. Kjesbu OS, Klungsøyr J, Kryvi H, Witthames PR, Walker MG. Fecundity, atresia, and egg size of cap-tive Atlantic cod (*Gadus morhua*) in relation to proximate body composition. Can J Fish Aquat Sci. 1991; 48: 2333–2343.

89. Lambert Y, Yaragina NA, Kraus G, Marteinsdottir G, Wright PJ. Using environmental and biological indi-ces as proxies for egg and larval production of marine fish. J Northwest Atl Fish Sci. 2003; 33: 115–159.

90. Fletcher DA, Wootton RJ. A hierarchical response to differences in ration size in the reproductive perfor-mance of female three-spined sticklebacks. J Fish Biol. 1995; 46: 657–668.