Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Comparative Assessment of the Reproductive Status of Female Atlantic Bluefin Tuna from the Gulf of Mexico and the Mediterranean Sea

  • Jessica M. Knapp ,

    jessieknapp@gmail.com

    Affiliations Department of Biological Sciences, University of New Hampshire, Durham, New Hampshire, United States of America, Large Pelagics Research Center, University of Massachusetts Amherst, Gloucester, Massachusetts, United States of America

  • Guillermo Aranda,

    Affiliation Departamento de Biología, Facultad de Ciencias del Mar y Ambientales, Campus de Excelencia Internacional del Mar (CEI-MAR), Universidad de Cádiz, Puerto Real, Cádiz, Spain

  • Antonio Medina,

    Affiliation Departamento de Biología, Facultad de Ciencias del Mar y Ambientales, Campus de Excelencia Internacional del Mar (CEI-MAR), Universidad de Cádiz, Puerto Real, Cádiz, Spain

  • Molly Lutcavage

    Affiliation Large Pelagics Research Center, University of Massachusetts Amherst, Gloucester, Massachusetts, United States of America

Abstract

Despite attention focused on the population status and rebuilding trajectory of Atlantic bluefin tuna (Thunnus thynnus), the reproduction and spawning biology remains poorly understood, especially in the NW Atlantic. At present, the eastern and western spawning populations are believed to exhibit different reproductive characteristics and, consequently, stock productivity. However, our study suggests that the two spawning populations, the Gulf of Mexico and the Mediterranean Sea, could show similar reproductive features and spawning strategies. Between 2007 and 2009, gonad samples from female Atlantic bluefin tuna were collected in the northern Gulf of Mexico (n = 147) and in the western Mediterranean Sea (n = 40). The histological and stereological analysis confirmed that sampled eastern and western bluefin tuna exhibit the same spawning duration (three months) but the spawning in the Gulf of Mexico begins one month earlier than in the Mediterranean Sea. Western bluefin tuna caught in the peak of the spawning season (May) showed a similar spawning frequency (60%) to the spawning peak observed in the Mediterranean Sea (June). Fecundity for the Gulf of Mexico fish () was lower but not significantly different than for fish sampled in the Mediterranean Sea (). Our study represents the first comparative histological analysis of the eastern and western spawning stocks whose findings, combined with new determinations of size/age at maturity and possible alternative spawning areas, might suggest basic life history attributes warrant further scientific and management attention.

Introduction

The reproductive biology of Atlantic bluefin tuna (ABFT; Thunnus thynnus, L. 1758) remains poorly understood despite the high economic value of this fishery and its exploitation throughout the Atlantic Ocean and the Mediterranean Sea. These uncertainties directly affect our understanding of the recruitment and productivity of the stock and could result in inefficient management of the fishery. The International Commission for the Conservation of Atlantic Tunas (ICCAT) manages the population as two stocks (eastern and western), separated at W based on two known spawning areas, the Gulf of Mexico (GMX) and the Mediterranean Sea (MED). This division assumes ABFT exhibit natal homing [1], [2] and have different maturation schedules [3], [4]. Nevertheless, electronic tagging and genetic/chemical marker studies have shown stock mixing on the foraging grounds indicating more complex population dynamics [1], [2], [5][7]. The eastern stock is estimated to be tenfold larger than the western stock [8]; consequently, the mixing rates are unbalanced with the eastern stock having greater influence on the western population, and any management action aimed at the eastern stock may indirectly affect the western stock [2]. Understanding the reproductive potential of both stocks is essential as it influences recruitment and hence the sustainability of the stocks and their capacity for supporting commercial fisheries [3],[9].

Electronic tagging and macroscopic examination of gonads are useful tools but, alone, are not sufficient for assessing population reproductive dynamics. Studies based on the histological analysis of gonads have allowed the characterization of the reproductive cycle of eastern spawners by the identification of maturation stages and estimations of several reproductive parameters [4], [10][13]. In the western and central Atlantic, these studies are scarce or impaired by large uncertainties [3], [9], [14]. Past studies examined gonad histology in fish from the Bahamas and the mid-Atlantic Bight and found a significant proportion of post-spawning females, but no ovaries with ovulated oocytes were observed [9], [15], [16]. However, these studies provided limited information since the samples were taken far from the known spawning areas and could not provide information about the reproductive condition of fish classified as active spawners [17], such as spawning frequency and/or fecundity [18].

While age at maturity in the eastern stock has been estimated as 3–4 y [4], [10], this parameter for the western spawning stock is the object of intense debate. ICCAT assumes an age at 50% maturity for the western stock of 8 y, but other studies provide estimations ranging from 4–16 y [9], [14], [19][21]. The younger maturation age observed in the eastern population [4], [9][11] leads to higher lifetime reproductive output and, consequently, larger productivity in the eastern than in the western stock.

Fecundity estimates allow the quantification of the reproductive capacity of individual fish and are essential for accurate assessment of the spawning stock biomass [22]. Assuming environmental characteristics between the MED and GMX spawning grounds are different, and ABFT exhibit natal homing, fecundity must be calculated separately for each stock.

An intensive evaluation of the maturity status of ABFT in the GMX compared with the MED would improve our comprehension of the reproductive connections between both stocks. In the present study, a histological and stereological analysis of gonads from ABFT caught in the GMX and the MED was undertaken to determine the reproductive status of female ABFT in the GMX through qualitative and quantitative histological studies, and compare the reproductive parameters estimated between both spawning grounds.

Materials and Methods

Sample collection and biometry

Female ABFT were sampled from commercial fisheries in the GMX and the MED between 2007 and 2009. The National Marine Fisheries Service Pelagic Observer Program sampled fish from longline vessels in the northern region of the Gulf of Mexico from February–July (). Eight samples obtained in February () and March () were not included in the stereological analysis because of low monthly sample size. In the MED, ABFT were sampled from longline fishing vessels on the western MED spawning grounds from mid June to mid July, 2008 (). For all samples, curved fork length () of each individual was measured to the nearest cm and converted to straight fork length () using the formula: [23]. Body mass () was calculated from by location and timing of catch according to ICCAT conversions [23]. Ovaries were immediately removed and weighed, and the volume of the pair of ovaries () was estimated from their mass () according to the equation [12]. The gonadosomatic index () was calculated from .

Histology

A subsample (0.5–1 ) was removed from the central portion of one ovary and fixed for at least 24 h in 4% formaldehyde (10% formalin) in phosphate buffer, 0.1 M, pH 7.2. Tissue samples were dehydrated through increasing concentrations of ethyl alcohol, cleared, and embedded in paraffin wax. Samples were cut into 5–6 m sections using a microtome, and stained with either haematoxylin-eosin (GMX) or haematoxylin-VOF (MED) [23]. All images were taken on a Leica DMI 6000b using Image-Pro Plus.

Based on the oocyte development, eight distinct types of developing follicles were distinguished: perinucleolar (PNF), lipid-stage (LSF), vitellogenic (VF), and oocyte maturation follicles (OMF), which consisted of migratory-nucleus (MNF) and hydrated (HF) follicles. Additionally, - and -stage atretic follicles (AF and AF, respectively) and postovulatory follicles (POF) were counted (Figure 1). Based on the most advanced follicle type in the ovary and the extent of atresia, fish were classified as inactive (IN), active non-spawning (ANS), or active spawning (AS) [17], [25].

thumbnail
Figure 1. Stages of oocyte development observed in Atlantic bluefin tuna sampled from the Gulf of Mexico and Mediterranean Sea.

A: lipid stage follicles (LSF); B: LSF, and early/late vitellogenic stage follicles (VF); C, D: migratory-nucleus follicles (MNF); E: hydrated follicle (HF); F, G: post-ovulatory follicles (POF); H: alpha (AF) and beta(AF) atretic follicles. All scale bars are 100 m.

https://doi.org/10.1371/journal.pone.0098233.g001

Spawning frequency

Spawning frequency was estimated by the postovulatory follicle method [26] as adapted by [27]. This method calculates the mean spawning fraction as the total number of spawning females whose ovaries show postovulatory follicles (POF) divided by the total number of mature females sampled.

Stereology

A stereological model-based method was applied to estimate numbers of the different categories of follicles according to the formula, where is the numerical density (number per unit volume) of the considered follicle type, is a shape coefficient, K is a size distribution coefficient, is the number of follicle transections per unit area, and is the partial area (volume density). was calculated by analyzing 10 digital micrographs from each ovary using ImageJ [18], [28]. Values of both and K previously calculated for ABFT were used for this study [11], [12], [29]. The total number of follicles (N) was calculated by extrapolating to the total ovarian volume, . Ovarian volume loss through processing was measured, and corrections were applied, 34.8% for MED and 43.3% for GMX. Finally, the relative number of follicles (number of follicles per gram) was estimated as .

Statistical analysis

Comparisons of means of the stereological and biometrical parameters among years were performed using the Kruskal-Wallis test, and parameters with no significant difference were regrouped by month (). Monthly variation was also analyzed using the Kruskal-Wallis test (). The Mann-Whitney U-test with Bonferroni correction was used to assess significant differences between pairs of months () [30].

Results

The for fish sampled in the GMX (n = 147) was 172–326 cm, and 120–240 cm for fish sampled in the MED (n = 40). The was 0.32–6.9 in the GMX and 0.3–5.8 in the MED. Mean of fish sampled in the GMX ( cm) was significantly larger (Mann-Whitney, ) than those sampled in the MED ( cm; Table 1). No significant differences in were observed within the GMX throughout the sampling period (Kruskal-Wallis, ). The mean was higher in the MED than in the GMX (Mann-Whitney, ), but this difference was only significant in June. Within the GMX, the for females sampled in June was significantly lower than for those sampled in the other months (Kruskal-Wallis, ; Table 2).

thumbnail
Table 1. Biometric data from Atlantic bluefin tuna, Thunnus thynnus, caught on longline fishing vessels operating in the Gulf of Mexico and Mediterranean Sea.

https://doi.org/10.1371/journal.pone.0098233.t001

thumbnail
Table 2. Stereological data for Atlantic bluefin tuna, Thunnus thynnus, sampled in the Gulf of Mexico and Mediterranean Sea.

https://doi.org/10.1371/journal.pone.0098233.t002

Histology

Histological analysis of ovarian tissue from the GMX sampled from April to June showed no significant differences in the gonad development between years. The number of inactive (IN) females was less than 20% of our sample except in 2007 (28.0%). The proportion of active non-spawning (ANS) females was 20.0% (2007), 50.0% (2008) and 30.0% (2009). The proportion of active spawning (AS) females was consistent at about 50% for all years (Figure 2). Given the lack of annual variation, GMX samples were pooled for all years for subsequent analyses.

thumbnail
Figure 2. Maturity stages for female Atlantic bluefin tuna sampled in the Gulf of Mexico from 2007–2009.

No significant difference was found between years for maturity stages, and thus, all years were pooled for future analyses. IN = inactive, ANS = active non-spawning, AS = active spawning.

https://doi.org/10.1371/journal.pone.0098233.g002

When ovary samples were arranged by month, consistent differences in the reproductive condition were observed with increasing maturation throughout the sampling period (Figure 3). In the GMX, the samples collected in February and March did not include any AS females. The proportion of AS females increased from April to June until a small decline was observed near the end of June. Near the end of the presumed spawning season in the MED (mid-June–mid-July), proportions of AS, ANS, and IN fish were similar to those observed in the GMX in May (Figure 3).

thumbnail
Figure 3. Maturity stages for female Atlantic bluefin tuna sampled in the Gulf of Mexico and Mediterranean Sea separated by month.

No significant difference was found between years for maturity stages, and thus, all years were pooled for future analyses. IN = inactive, ANS = active non-spawning, AS = active spawning.

https://doi.org/10.1371/journal.pone.0098233.g003

Spawning frequency

Spawning fraction in the GMX was estimated to be 0.45 throughout the sampling period as 49 out of 108 mature females contained POFs in their ovaries. When this parameter was calculated by month, the proportion of females with POFs caught in April was lower than that in May and June. When the spawning frequency was calculated considering only AS females, its value increased significantly and remained similar among months. The spawning fraction in the MED samples was higher (0.60) than in the GMX (Table 3).

thumbnail
Table 3. Spawning frequency and spawning interval estimates for Atlantic bluefin tuna, Thunnus thynnus, sampled from the Gulf of Mexico and the Mediterranean Sea.

https://doi.org/10.1371/journal.pone.0098233.t003

Stereology

No significant differences were found between months for the stereological parameters of atretic follicles (AF and AF) in fish sampled in the GMX (Table 2). No significant differences were found for the numerical density of LSF (LSF) in the GMX despite the decrease observed throughout the reproductive season. Nevertheless, the total number of LSF (NLSF) and the relative number of LSF (NgLSF) estimated in GMX fish in April was significantly higher than in June. In general, stereological counts of LSF quantified in the MED showed significantly lower values than those from the GMX sampled in April and May. Although the mean number of VF per (VF) and the mean number of VF (NVF) remained unchanged in the GMX, the relative number of VF (NgVF) was significantly lower in June than in May. VF and NVF estimated in the MED were similar to the GMX values in June; however, because eastern fish were smaller than western fish on average, the relative number of VF (NgVF) was much higher for MED fish. The low number of females with MNF and HF is likely the cause for finding no significant differences for these stages throughout the sampling period in the GMX. The highest number of POF (POF, NPOF, and NgPOF) occurred at the beginning of the sampling period (April) and the lowest values were observed in June; however, no significant differences were found either among months or between the spawning areas. The realized fecundity estimated as the mean absolute number of POFs (SD) was for fish from the GMX. This corresponds to relative batch fecundity () of . In the MED, the estimated mean number of POF was and the relative fecundity (Table 2).

Discussion

Although not spatially and temporally exhaustive, this study represents the first attempt since the 1970s to accurately assess the spawning condition of ABFT sampled on the GMX spawning grounds and is the first histological and stereological comparison of the eastern and western spawning stocks.

Histological analysis of gonad samples from the spawning grounds throughout the spawning season is essential for evaluating the reproductive condition and performance of ABFT. Additionally, systematic sampling across the extent of the spawning grounds allows the study of temporal variations in key reproductive parameters, such as sex ratio, proportion of mature fish, spawning frequency, and spawning periodicity and fecundity. This information was lacking for the western stock resulting in large uncertainties for stock assessment and evaluation of productivity [8]. Since the implementation of the moratorium on directed fishing for ABFT in the GMX, federal fisheries observers have sampled ABFT caught as bycatch in the yellowfin tuna and swordfish longline fisheries (Thunnus albacares and Xiphias gladius, respectively). Such restrictions prevent comprehensive size sampling of spawning ABFT in the GMX and hinder the determination of the spatial and temporal extent of western spawning [31].

As a result of the bycatch sampling, previous studies have lacked small/medium fish ( cm) leading to larger/older size and age at maturity estimates [20], [21]. Despite electronic tagging data that show presumed mature fish outside of known spawning areas [1], [32][34], these studies assumed western bluefin only spawn in the GMX. A New England study suggested previous ABFT sampling did not accurately represent the spawning size range of the western population because it only included fish sampled by longliners on known spawning grounds rather than all size classes sampled throughout their range [14]. Gear type, size selectivity, and vertical distribution of tuna by size also influence the size of spawners sampled by commercial fishing fleets [12], [35].

In this study, the smallest tuna sampled from the GMX, 172 cm and estimated age 7–8 y [36], had ripe ovaries with numerous recent POFs. This is consistent with an earlier study finding mature fish at 8 y [14] but not the recently proposed 12–16 y [20], [21]. In order to fully understand the reproductive dynamics of the western spawning stock, the maturity ogive should be revised using comprehensive size sampling over larger temporal scales including histological examination of the ovaries and endocrine profiling. New studies utilizing endocrine hormones developed and calibrated in captive eastern ABFT [37] indicate that western ABFT become sexually mature at less than 8 y [38].

Prior to the US moratorium, catch records indicated the presence of only giant ABFT ( cm) in the GMX [3]. As opposed to the current management paradigm of western ABFT maturing at an older age than the eastern stock, fish may exhibit size and temporal segregation on the spawning grounds as is seen with the eastern spawning stock [3], [32], [39] and in Pacific bluefin tuna (Thunnus orientalis) [40]. There is indirect evidence that smaller fish may utilize alternative spawning locations, such as the Caribbean Sea, the Bahamas, or the Gulf Stream margins because ripe fish have been sampled there [9], [15], [41]. Given this evidence, it is possible that smaller bluefin spawn in alternative locations within the GMX as previous sampling has been concentrated in the north/central Gulf where US longline vessels operate. An ABFT life history model predicts that smaller/younger maturing fish should have shorter migration routes and spawn in areas closer to feeding areas than larger, older fish with higher energy reserves [42].

Electronic tagging results have consistently shown annual migration patterns of giant ABFT not entering either known spawning ground before returning to northern feeding grounds [1], [32][34]. It is possible that western ABFT spawn over a broader area of regions with oceanographic conditions appropriate for larval development than previously assumed [3], [14], [32], [34], [43]. Recent larval cruises found bluefin larvae outside the GMX [44], but spawning areas beyond the northern GMX await histological validation. Pop-up satellite tagged 2–5 year old ABFT did not enter the GMX or MED during presumed spawning periods (April–June) [45]. Nevertheless, some fish lingered in subtropical seas north of the Bahamas and in the southern mid-Atlantic Bight, areas visited by tagged adults [1], [5], [33], [46].

Our histological examination of the GMX bluefin ovaries revealed differences in follicle maturation between months throughout the spawning period. While samples collected in February and March contained no active spawning (AS) individuals, 31% of samples collected in April were AS individuals. As the spawning season progressed, the number of AS females increased and peaked in May (60%). The observed in June was significantly lower than other months indicating a decrease in ovarian size, and thus, the impending end of the spawning season for this region. The statistical results of the stereological analysis are consistent with previous findings and with the progression of the spawning season[11]. As the spawning season progresses, LSFs become less frequent indicating high levels of recruitment to VFs thereby compensating for losses caused by atresia or spawning. Similarly, the relative number of VFs (VF) was significantly higher in April and May than in June indicating a decrease in the recruitment of VFs as the end of the spawning season approaches.

In spawning fish, atresia is a natural mechanism for regulating the number of eggs spawned. Alternatively, massive atresia can indicate a cessation of oocyte maturation and/or spawning activity[47]. The GMX samples showed relatively low and stable levels of AF throughout the spawning season indicating ABFT found in the GMX are in favorable condition for oocyte maturation and spawning. While these fish appear to be actively spawning, tagged ABFT of presumed reproductively mature size observed outside the known spawning areas during the spawning period could be skipping spawning due to unfavorable body condition [14], [34]. ABFT sampled on the New England and Canadian foraging grounds have had periods of reduced somatic condition [48], [49] possibly accounting for an increased incidence of skipped spawning [50]. The incidence of skipped spawning in ABFT, however, is unknown, and modeling results show it is less likely to occur in larger, older fish in positive energy balance [42]. Giant ABFT sampled on western foraging grounds in the fall [14], [48] and in the MED in the spring and early summer [51] have extensive perigonadal fat and somatic lipid stores, and thus seem unlikely candidates for skipped spawning.

While ABFT have been observed on the western spawning grounds for months during the spawning season [34], individuals are believed to be actively spawning for only a few weeks [9], [52], [53]. These findings, albeit fishery-dependent, define the temporal borders of the reproductive events occurring in the north/central part of the GMX indicating the spawning season runs from April to June with maximum spawning activity in May. However, ABFT begin entering the GMX in late November, and those arriving in winter experience warm water masses of C in the lower GMX [34]. Reproductive sampling has been primarily conducted in the northern GMX and US territorial seas in late spring [9], [52]. Spawning activity occurring earlier in other areas of the GMX awaits confirmation by broader sampling, especially in Mexican territorial seas.

The relatively low proportion of GMX females with POFs in their ovaries () contrasts with the high spawning frequency (60%) observed in the western MED [12], [13]. The lower spawning frequency observed in the GMX could be the result of bias associated with utilizing the yellowfin and swordfish fisheries as the only sampling method. As long as bluefin reproductive studies rely on bycatch in commercial fisheries, it is not possible to obtain an unbiased, accurate assessment of ABFT reproduction. Given these constraints, it is important to note the temporal and spatial aspects of the sampling as well as the fishing gear used for any bluefin maturity study.

Stereological methods have often been used as an accurate tool for estimating fecundity in fishes, including eastern ABFT [11], [22], [29], [54], [55]. Realized fecundity can be estimated through stereological counts of POFs, whereas the number of MNF is an estimation of the potential fecundity [29], [56]. In this study, the mean relative batch fecundity was calculated directly from stereological counts of POFs and showed a decrease as the season progressed. This is atypical for indeterminate spawners [57], [58] but is likely due to the selective nature of sampling bluefin as bycatch in a longline fishery. One study suggested that monthly variation of fecundity may be masked by a decrease in the condition factor of fish appearing later on the spawning ground [59]. Although significant differences were not found between months, the highest value in the relative fecundity occurred early in the season (April), even though the number of AS females was still quite low. ABFT entering the GMX early might exhibit higher reproductive potential than those arriving later due to the good condition acquired on the foraging grounds. Otherwise, the lower spawning frequency observed on the western spawning grounds could be a consequence of migration distance [42] or decreased body condition observed on the western foraging grounds [48]. While being the first to arrive on the spawning grounds might provide increased resource availability for offspring, arriving in poor condition could decrease larval survival rates [51], [60].

The fecundity of eastern spawning bluefin was estimated at [12] and [13] for potential and realized fecundity, respectively. Our results show a lower realized fecundity for western spawners () than eastern spawners; however, this difference was not statistically significant. This is in agreement with previous work indicating realized fecundity is not proportional to body size [13]. Although the differences in fecundity and maturity schedules of both ABFT stocks can be masked by sampling bias and complex population dynamics, fecundity has been shown to vary within a given species as a result of different adaptations to environmental habitats [61]. The unexpected variation of fecundity shown among large ABFT could be related to the environmental adaptations and the balance between body condition and energetic expenses during migration [42], [51].

Eastern and western ABFT spawning sites seem to exhibit the same periodicity (three months), but spawning in the northern GMX occurs one month earlier than in the western Mediterranean spawning ground. This is likely due to specific oceanographic conditions and the early warmer temperatures observed in the GMX (data from CCA-UNAM). In this study, we have observed similar values in bluefin reproductive parameters showing that the spawning condition of Mediterranean spawners from mid-June to mid-July is comparable with the reproductive peak observed in the GMX in May. Depth and temperature associations of electronically tagged ABFT entering the GMX in winter also suggest that it also serves as a foraging ground [34], as it does for other tunas and billfish species.

A more holistic view of the population dynamics of ABFT requires life history characteristics, reproductive profiles, and spawning areas and periodicity to be well defined, especially since they undoubtedly will change with shifts in climate and ocean productivity [62]. The extent and quality of lipids acquired by tunas before they arrive in spawning areas will affect eggs and larvae, and therefore, overall stock biomass [42], [51]. Future work should address energetic relationships between reproduction, migration, and early life history through modelling, biological sampling, and the development of smart tags to detect actual spawning events.

Acknowledgments

We thank Gilad Heinisch, John Logan, Jim Franks, Ed Chesney, Larry Beerkircher, Steve Turner, Sascha Cushner, and the NMFS Pelagic Observer Program for helping with sampling in the Gulf of Mexico. We also thank Captain Carl Roby (FV Cynthia Renee) and Captain Tim (FV Orion) for sampling help in the Gulf of Mexico. Barbara Brothwell and the UNH Veterinary Diagnostic Lab provided significant histology support. We also thank David Macias (IEO) for providing gonad samples from the Balearic Sea and Agustín Santos (UCA) for technical support in stereological and histological procedures.

Author Contributions

Conceived and designed the experiments: JK GA AM ML. Performed the experiments: JK GA. Analyzed the data: JK GA. Contributed reagents/materials/analysis tools: JK GA AM ML. Wrote the paper: JK GA AM ML.

References

  1. 1. Block BA, Teo SLH, Walli A, Boustany A, Stokesbury MJW, et al. (2005) Electronic tagging and population structure of Atlantic bluefin tuna. Nature 434: 1121–1127.
  2. 2. Rooker JR, Secor DH, De Metrio G, Kaufman AJ, Belmonte RA, et al. (2008) Evidence of trans-Atlantic movement and natal homing of bluefin tuna from stable isotopes in otoliths. Marine Ecology Progress Series 368.
  3. 3. Mather FJ, Mason Jr JM, Jones AC (1995) Historical document: Life history and fisheries of Atlantic bluefin tuna. NOAA Technical Memorandum NMFS-SEFSC-370.
  4. 4. Corriero A, Karakulak S, Santamaria N, Deorio M, Spedicato D, et al. (2005) Size and age at sexual maturity of female bluefin tuna (Thunnus thynnus, L. 1758) from the Mediterranean Sea. Journal of Applied Ichthyology 21: 483–486.
  5. 5. Block BA, Dewar H, Blackwell SB, Williams TD, Prince ED, et al. (2001) Migratory movements, depth preferences, and thermal biology of Atlantic bluefin tuna. Science 293: 1310–1314.
  6. 6. Pujolar J, Roldán M, Pla C (2003) Genetic analysis of tuna populations, Thunnus thynnus thynnus and T-alalunga. Marine Biology 143: 613–621.
  7. 7. Dickhut RM, Deshpande AD, Cincinelli A, Cochran MA, Corsolini S, et al. (2009) Atlantic bluefin tuna (Thunnus thynnus) population dynamics delineated by organochlorine tracers. Environmental Science & Technology 43: 8522–8527.
  8. 8. Fromentin JM, Powers JE (2005) Atlantic bluefin tuna: population dynamics, ecology, fisheries and management. Fish and Fisheries 6: 281–306.
  9. 9. Baglin RE (1982) Reproductive biology of western Atlantic bluefin tuna. Fishery Bulletin 80: 121–134.
  10. 10. Rodríguez-Roda J (1967) Fecundidad del atún, Thunnus thynnus (L.), de la costa sudatlántica de España. Investigación Pesquera 31: 33–52.
  11. 11. Medina A, Abascal FJ, Megina C, Garcia A (2002) Stereological assessment of the reproductive status of female Atlantic northern bluefin tuna during migration to Mediterranean spawning grounds through the Strait of Gibraltar. Journal of Fish Biology 60: 203–217.
  12. 12. Medina A, Abascal FJ, Aragón L, Mourente G, Aranda G, et al. (2007) Influence of sampling gear in assessment of reproductve parameters for bluefin tuna in the western Mediterranean. Marine Ecology Progress Series 337: 221–230.
  13. 13. Aranda G, Medina A, Santos A, Abascal FJ, Galaz T (2013) Evaluation of Atlantic bluefin tuna reproductive potential from ovarian histology and its use in stock assessments. Journal of Sea Research 76: 154–160.
  14. 14. Goldstein JL, Heppell SA, Cooper A, Brault S, Lutcavage ME (2007) Reproductive status and nutritional condition of Atlantic bluefin tuna in the Gulf of Maine. Marine Biology 151: 2063–2075.
  15. 15. Rivas LR (1954) A preliminary report on the spawning of the western north Atlantic bluefin tuna (Thunnus thynnus) in the Straits of Florida. Bulletin of Marine Science of the Gulf of Mexico and Caribbean 4: 302–322.
  16. 16. Baglin RE (1976) A preliminary study of the gonadal development and fecundity of the western Atlantic bluefin tuna. International Commission for the Conservation of Atlantic Tunas Collective Volume of Scientific Papers 5: 279–289.
  17. 17. Schaefer KM (1998) Reproductive biology of yellowfin tuna (Thunnus albacares) in the eastern Pacific Ocean. Bulletin of the Inter-American Tropical Tuna Commission 21: 201–249.
  18. 18. Clay D (1991) Atlantic bluefin tuna (Thunnus thynnus thynnus (L.)): a review. In: World meeting on stock assessment of bluefin tunas: strengths and weaknesses, Inter-American Tropical Tuna Commission Special Report No. 7. pp. 91–179.
  19. 19. Westman JR, Neville WC (1942) The tuna fishery of Long Island, New York. A survey conducted by the United States Department of the Interior, Fish and Wildlife Service, in cooperation with the County Executive and the Board of Supervisors, Nassau County, Long Island, New York.
  20. 20. Diaz GA, Turner SC (2007) Size frequency distribution analysis, age composition, and maturity of western bluefin tuna in the Gulf of Mexico from the U.S. (1981–2005) and Japanese (1975–1981) longline fleets. International Commission for the Conservation of Atlantic Tunas Collective Volume of Scientific Papers 60: 1160–1170.
  21. 21. Diaz GA (2011) A revision of western Atlantic bluefin tuna age of maturity derived from size samples collected by the Japanese longline fleet in the Gulf of Mexico (1975–1980). International Commission for the Conservation of Atlantic Tunas Collective Volume of Scientific Papers 66.
  22. 22. Murua H, Kraus G, Saborido-Rey F, Witthames P, Thorsen A, et al. (2003) Procedures to estimate fecundity of marine fish species in relation to their reproductive strategy. Journal of Northwest Atlantic Fishery Science 33: 33–54.
  23. 23. ICCAT Website. Available: http://www.iccat.int/documents/stats/convers.pdf. Accessed 2014 May 6.
  24. 24. Gutiérrez M (1967) Coloración histológica para ovarios de peces, crustáceos y moluscos. Investigación Pesquera 31: 265–271.
  25. 25. Schaefer KM (2001) Reproductive biology of tunas. In: Block BA, Stevens ED, editors, Tuna: physiology, ecology, and evolution, London: Academic Press, volume 19 of Fish Physiology. pp. 225–270.
  26. 26. Stauffer G, Picquelle S (1981) The 1980 and 1981 egg production estimates of anchovy spawning biomass. Southwest Fisheries Center La Jolla Laboratory, National Marine Fisheries Service, NOAA.
  27. 27. Hunter JR, Macewicz BJ (1985) Measurement of spawning frequency in multiple spawning fishes. In: Lasker R, editor, An egg production method for estimating spawning biomass of pelagic fish: application to the northern anchovy, Engraulis mordax, NOAA Technical Report NMFS. pp. 79–94.
  28. 28. Rasband WS (1997–2011) ImageJ, U.S. National Institutes of Health, Bethesda, Maryland, USA. http://imagejnihgov/ij/.
  29. 29. Aragón L, Aranda G, Santos A, Medina A (2010) Quantification of ovarian follicles in bluefin tuna Thunnus thynnus by two stereological methods. Journal of Fish Biology 77: 719–730.
  30. 30. Sokal RR, Rohlf FJ (1998) Biometry: the principles and practice of statistics in biological research. New York: W.H. Freeman and Company.
  31. 31. Rosenfeld H, Mylonas CC, Bridges CR, Heinisch G, Corriero A, et al. (2012) GnRHa-mediated stimulation of the reproductive endocrine axis in captive Atlantic bluefin tuna, Thunnus thynnus. General and Comparative Endocrinology 175: 55–64.
  32. 32. Lutcavage ME, Brill RW, Skomal GB, Chase BC, Howey PW (1999) Results of pop-up satellite tagging of spawning size class fish in the Gulf of Maine: do North Atlantic bluefin tuna spawn in the mid-Atlantic? Canadian Journal of Fisheries and Aquatic Sciences 56: 173–177.
  33. 33. Sibert JR, Lutcavage ME, Nielsen A, Brill RW, Wilson SG (2006) Interannual variation in largescale movement of Atlantic bluefin tuna (Thunnus thynnus) determined from pop-up satellite archival tags. Canadian Journal of Fisheries and Aquatic Sciences 63: 2154–2166.
  34. 34. Galuardi B, Royer F, Golet W, Logan J, Lutcavage ME (2010) Complex migration routes of Atlantic bluefin tuna question current population structure paradigm. Canadian Journal of Fisheries and Aquatic Sciences 67: 966–976.
  35. 35. Davis TL, Farley JH (2001) Size distribution of southern bluefin tuna (Thunnus maccoyii) by depth on their spawning ground. Fishery Bulletin 99: 381–386.
  36. 36. Restrepo VR, Diaz GA, Walter JF, Neilson JD, Campana SE, et al. (2010) Updated estimate of the growth curve of western Atlantic bluefin tuna. Aquatic Living Resources 23: 335–342.
  37. 37. Berkovich N, Corriero A, Santamaria N, Mylonas CC, Vassallo-Aguis R, et al. (2013) The intrapituitary relationship of follicle stimulating hormone and luteinizing hormone during pubertal development in Atlantic bluefin tuna (Thunnus thynnus). General and Comparative Endocrinology 194: 10–23.
  38. 38. Heinisch G, Gordin H, Rosenfeld H, Knapp JM, Lutcavage ME (2013) Sexual maturity in western Atlantic bluefin tuna. In prep.
  39. 39. Heinisch G, Corriero A, Medina A, Abascal FJ, de la Serna JM, et al. (2008) Spatial-temporal pattern of bluefin tuna (Thunnus thynnus L. 1758) gonad maturation across the Mediterranean Sea. Marine Biology 154: 623–630.
  40. 40. Itoh T (2006) Sizes of adult bluefin tuna, Thunnus orientalis, in different areas of the western Pacific Ocean. Fisheries Science 72: 53–62.
  41. 41. Wilson PC, Bartlett MR (1967) Inventory of U.S. exploratory longline fishing effort and catch rates for tunas and swordfish in the western Atlantic.
  42. 42. Chapman EW, Jørgensen C, Lutcavage ME (2011) Atlantic bluefin tuna (Thunnus thynnus): a state-dependent energy allocation model for growth, maturation, and reproductive investment. Canadian Journal of Fisheries and Aquatic Sciences 68: 1934–1951.
  43. 43. Lutcavage M, Galuardi B, Lam CH (2012) Predicting potential alternative spawning grounds of western Atlantic bluefin tuna based on electronic tagging results. In prep.
  44. 44. Muhling BA, Lamkin JT, Quattro JM, Smith RH, Roberts MA, et al. (2011) Collection of larval bluefin tuna (Thunnus thynnus) outside documented western Atlantic spawning grounds. Bulletin of Marine Science
  45. 45. Galuardi B, Lutcavage M (2012) Dispersal routes and habitat utilization of juvenile Atlantic bluefin tuna, Thunnus thynnus, tracked with mini PSAT and archival tags. PLoS One 7: e37829
  46. 46. Wilson SG, Lutcavage ME, Brill RW, Genovese MP, Cooper AB, et al. (2005) Movements of bluefin tuna (Thunnus thynnus) in the northwestern Atlantic Ocean recorded by pop-up satellite archival tags. Marine Biology 146: 409–423.
  47. 47. Tyler CR, Sumpter JP (1996) Oocyte growth and development in teleosts. Reviews in Fish Biology and Fisheries 6: 287–318.
  48. 48. Golet WJ, Cooper AB, Campbell R, Lutcavage ME (2007) Decline in condition of northern bluefin tuna (Thunnus thynnus) in the Gulf of Maine. Fishery Bulletin 105: 390–395.
  49. 49. Paul SD, Hanke A, Vanderlaan ASM, Busawon D, Neilson JD (2011) Indices of stock status from the 2009 Canadian bluefin tuna fishery. International Commission for the Conservation of Atlantic Tunas Collective Volume of Scientific Papers 66: 1170–1203.
  50. 50. Rideout RM, Rose GA, Burton MPM (2005) Skipped spawning in female iteroparous fishes. Fish and Fisheries 6: 50–72.
  51. 51. Mourente G, Megina C, Díaz-Salvago E (2002) Lipids in female northern bluefin tuna (Thunnus thynnus thynnus L.) during sexual maturation. Fish Physiology and Biochemistry 24: 351–363.
  52. 52. Richards WJ (1976) Spawning of bluefin tuna (Thunnus thynnus) in the Atlantic Ocean and adjacent seas. International Commission for the Conservation of Atlantic Tunas Collective Volume of Scientific Papers 5: 267–278.
  53. 53. Teo SLH, Boustany A, Dewar H, Stokesbury MJW, Weng KC, et al. (2007) Annual migrations, diving behavior, and thermal biology of Atlantic bluefin tuna, Thunnus thynnus, on their Gulf of Mexico breeding grounds. Marine Biology 151: 1–18.
  54. 54. Coward K, Bromage NR (2002) Quantification of ovarian condition in fish: a safer, more precise alternative to established methodology. Aquatic Living Resources 15: 259–261.
  55. 55. Kjesbu OS, Fonn M, Gonzáles BD, Nilsen T (2010) Stereological calibration of the profile method to quickly estimate atresia levels in fish. Fisheries Research 104: 8–18.
  56. 56. Aranda G, Aragón L, Corriero A, Mylonas CC, de la Gándara F, et al. (2011) GnRHa induced spawning in cage-reared Atlantic bluefin tuna: An evaluation using stereological quantification of ovarian post-ovulatory follicles. Aquaculture 317: 255–259.
  57. 57. Murua H, Motos L, Lucio P (1998) Reproductive modality and batch fecundity of the European hake (Merluccius merluccius L.) in the Bay of Biscay. California Cooperative Oceanic Fisheries Investigations Report 39: 196–203.
  58. 58. Murua H, Motos L (2006) Reproductive strategy and spawning activity of the European hake Merluccius merluccius (L.) in the Bay of Biscay. Journal of Fish Biology 69: 1288–1303.
  59. 59. Kjesbu OS, Witthames PR, Solemdal P, Greer Walker M (1998) Temporal variations in the fecundity of Arcto-Norwegian cod (Gadus morhua) in response to natural changes in food and temperature. Journal of Sea Research 40: 303–321.
  60. 60. Donelson JM, Munday PL, McCormick MI (2009) Parental effects on offspring life histories: when are they important? Biology Letters 5: 262–265.
  61. 61. Witthames PR, Greer-Walker M, Dinis MT, Whiting CL (1995) The geographical variation in the potential fecundity of Dover sole Solea solea L. from European Shelf edge waters during 1991. Journal of Sea Research 34: 45–58.
  62. 62. Muhling BA, Lee SK, Lamkin JT, Liu Y (2011) Predicting the effects of climate change on bluefin tuna (Thunnus thynnus) spawning habitat in the Gulf of Mexico. ICES Journal of Marine Science 68: 1051–1062.