Diversity of feeding strategies in loggerhead sea turtles from the Cape Verde archipelago

As biodiversity worldwide is decreasing, to preserve adaptive potential, the importance of maintaining species’ genetic and trait diversities is increasing. An efficient foraging strategy is a critical trait for an organism’s fitness, as it affects its physiology and reproduction. Understanding such strategies is especially relevant for species with long feeding migrations such as sea turtles. Using carbon and nitrogen stable isotopes combined with mitochondrial sequencing, we explored the diversity of feeding strategies in genetically differentiated nesting groups of loggerhead sea turtles (Caretta caretta) within the Cape Verde Archipelago. Here, we reveal a pattern where turtles from most islands use two distinct oceanic feeding strategies, including one putatively linked to a 15N-enriched zone of the West African upwelling area. On the Eastern island of Boavista, an additional third strategy exists used by turtles feeding mostly neritically. Contrary to previous paradoxical assumptions, oceanic turtles, that represent the vast majority of the population, are not smaller than neritic turtles and therefore do not seem to feed in a suboptimal environment. Our results also suggest that the number of feeding strategies may correlate with demography, whereby a greater feeding strategy diversity matches demographic signs of recent expansion after a population bottleneck for turtles nesting on the island of Boavista. Overall, the feeding ecology of Cape Verde loggerhead turtles is complex and likely shaped by an interaction between environmental and population parameters. Our results stress the importance of conservation efforts to prevent loss of critical diversity in endangered species.


Introduction
As the number of extant species continues to fall at an unprecedented rate (Pimm et al. 2014), the importance of preserving diversity worldwide has never been higher, whether to protect species' genetic diversity, traits, or 1 3 130 Page 2 of 12 behaviours (Isaac and Cowlishaw 2004;Myers et al. 2000). An efficient foraging strategy is one of the most important components for the survival and Darwinian fitness of an organism (Le Galliard et al. 2004;Stephens and Krebs 1986). The evolution of foraging strategies is the result of various selective pressures to maximize resource acquisition for growth, maintenance, or reproduction as well as to reduce encounters with predators (Börger et al. 2008;Sims et al. 2008;Stephens and Krebs 1986). Optimal feeding strategies are therefore likely to be habitat or location specific.
Several species have evolved different feeding strategies which may be linked to their jaw or beak morphology (e.g., Abzhanov et al. 2004;Albertson et al. 2005) but also to their size, allowing them to access specific food items (e.g., Hawkes et al. 2006). Understanding the origin of traits underlying feeding strategies, as well as their variation, can offer unique insights into the ecology of organisms and help maintain their adaptive potential for improved conservation measures (Eizaguirre and Baltazar-Soares 2014).
Stable isotopes have become important tools for conservation biologists. Particularly, carbon (C, δ 13 C) and nitrogen (N, δ 15 N) provide information on the specific habitats used by the target species as well as on their trophic position and the source of primary production (Hobson 1999;Post 2002;Reich et al. 2007). Relatively slow turnover of isotopes in certain tissues also means that migrating consumers continue to exhibit stable isotope values representative of their feeding grounds even after leaving them. This characteristic makes δ 13 C and δ 15 N excellent markers to identify the presence of various feeding strategies in a population of migrating and difficult to access organisms.
Sea turtles are prime examples of marine species whose migrations to and from feeding grounds can span thousands of kilometres (Godley et al. 2008;Hawkes et al. 2006) or entire ocean basins (Bolten and Witherington 2003). Adult female sea turtles are capital breeders and, as such, spend most of their time gaining energy for reproduction at their feeding grounds, only returning to their natal area every so often to breed and lay eggs (Plotkin 2003;Stearns 1992;Zbinden et al. 2011). During the nesting period, female turtles are accessible for sampling on land with minimal stress imposed on them.
The loggerhead sea turtle (Caretta caretta) represents the most abundant sea turtle species in warm temperate regions (Bolten and Witherington 2003;Dodd 1988;Eder et al. 2012). At a global level, its population size is low despite a recent upgrade on the IUCN Red List from 'Endangered' to 'Vulnerable'. The species is, however, still classified as "Endangered" (Cape Verde) and even "Critically Endangered" (North Indian Ocean & South Pacific) in certain other rookeries (IUCN 2016).
Recent stable isotope studies, bolstered by improved satellite tracking technology, have allowed monitoring of loggerhead turtles throughout their feeding migrations. Numerous studies report a dichotomy of foraging strategies with individuals of the same population following either an oceanic or neritic feeding strategy (Cape Verde Rookery-Eder et al. 2012;Hawkes et al. 2006, East USA Rookery-McClellan et al. 2010McClellan and Read 2007;Japanese Rookery-Hatase et al. 2002. Oceanic individuals feed epipelagically and opportunistically on neustonic organisms, e.g., jellyfish or crustaceans (Frick et al. 2009). Neritic individuals feed upon benthic prey items such as gastropods and other slow-moving and/or sessile organisms (Frick et al. 2009;Hopkins-Murphy et al. 2003;Plotkin 2003), along with pelagic prey inhabiting coastal habitats such as jellyfish (McClellan et al. 2010;Wallace et al. 2009).
Due to the abundance of nutrient-rich benthic fauna in neritic feeding sites compared to those in often low productivity oceanic habitats, neritic turtles are thought to be larger than their oceanic counterparts (Cape Verde Rookery-Eder et al. 2012;Hawkes et al. 2006;Japanese Rookery-Hatase et al. 2002Tsukamoto 2008, Western Australian Rookery-Thomson et al. 2012). Whilst this prediction seems to hold in the Northwest Atlantic, Japanese and central Mediterranean populations (Hatase et al. 2010;Hawkes et al. 2011;Zbinden et al. 2011), most individuals sampled from the Cape Verde archipelago and Oman are oceanic foragers (Cape Verde = e.g. 87.5%, Eder et al. 2012; Oman = 80%, Rees et al. 2010). This reveals an evolutionary paradox whereby most of the turtles from those rookeries may be using an a priori suboptimal feeding strategy that is assumed to be heritable, which raises the question about how such a behaviour has evolved and is maintained (Smith and Price. 1973).
The Cape Verde archipelago consists of ten volcanic islands ( Fig. 1), located approximately 570 km off the West coast of Africa in the North Atlantic Ocean. The archipelago sustains the third largest aggregation of loggerhead turtles on the planet (after Florida and Oman, Marco et al. 2012). There, turtles are classified as "Endangered" (IUCN 2016), due to traditional poaching, fisheries bycatch and coastal development (López-Jurado et al. 2003;Taylor and Cozens 2010). While the majority of nesting occurs on the most Eastern islands (Boavista, Sal, Maio), each island supports independent nesting groups with unique genetic diversity of conservation value (Stiebens et al. 2013a). Previous research undertaken on the island of Boavista revealed the existence of at least two distinct feeding strategies (Cardona et al. 2017;Eder et al. 2012), oceanic and neritic, with turtles smaller than 90 cm being oceanic, while turtles > 90 cm were considered neritic feeders (Eder et al. 2012). Those strategies are hence characterized by different feeding locations and different prey items. More recently, δ 13 C and δ 15 N readings from scutes of 14 individuals showed that turtles can shift between foraging habitats (Cardona et al. 2017). Here, we determined diversity of feeding ecology in genetically distinct nesting groups using stable isotope analyses from skin tissue across the entire archipelago. We also endeavour to address the paradox of the use of suboptimal feeding strategies by the majority of turtles. We then combined this information with mitochondrial DNA sequencing to test for correlations between population demography and diversity of feeding strategies.

Materials and methods
Overall, sampling took place during the nesting seasons 2012-2014 on six different islands of the Cape Verde archipelago: Boavista, Fogo, Maio, São Vicente, Santa Luzia and Santo Antão (Fig. 1). Sampling consisted of collecting one or two pieces of 3 mm non-keratinized tissue from the front flippers of nesting females using a single-use sterile scalpel, immediately following egg deposition (Stiebens et al. 2013b). One was used for genetic analyses, while the other was used for stable isotopes. At that time, females were tagged with coded inconel and/or PIT (AVID) tags to avoid sampling duplication. Curved carapace length (CCL) was also measured (± 0.1 cm).

Stable isotopes: preparation and readings
We used samples from our first encounter with an individual in the nesting season to avoid possible confounding effects linked to local feeding, if it exists. We collected 419 samples, 287 of which were also sampled for mtDNA sequencing. Samples were initially washed in distilled water for approximately 1 min to remove sand. They were dried at 60 °C for 48 h. Between 0.7 and 1.3 µg of sample was cut and weighed into tin capsules (4 mm), before being combusted using a continuous flow isotope ratio mass spectrometer (Integra2, Sercon) and analysed for nitrogen and carbon elements simultaneously. Repeated run of an internal standard (casein) showed analytical error to be ca. ± 0.1‰ for both C and N.

DNA sampling, and mitochondrial sequencing of the control region
Skin samples were extracted using the DNeasy ® 96 Blood & Tissue Kit (QIAGEN, Hilden, Germany) following the manufacturer's protocol except for the final elution step which was conducted in twice 75 μl of AE Buffer. The long fragment (here, 750 bp) of control region of the mitochondrial DNA was amplified using the primers LCM15382 (5′-GCT TAA CCC TAA AGC ATT GG-3′) and H950 (5′-GTC  Stiebens et al. (2013a). PCR products were cleaned with ExoSAP-IT ® following the manufacturer's protocol. Cycle sequencing reactions were performed with Big Dye ® Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, Darmstadt, Germany). Sequences were obtained from the forward direction (primer LCM15382). Sequencing was performed with an ABI 3730 Genetic Analyzer (Applied Biosystems, Darmstadt, Germany).

Demographic history of nesting groups
Sequences were assembled in Codon Code Aligner v5.0 (CodonCode Corporation, Dedham, Massachusetts) and ambiguities were corrected by hand. All the amplified mitochondrial sequences were classified accordingly to the standardized nomenclature of the Archie Carr Centre for Sea Turtle Research (http://accst r.ufl.edu). The entire data set was aligned in Muscle v8.3.1 (Edgar 2004). Population structure of the archipelago was also analysed using mtDNA with Wright's fixation index (F ST ) implemented within the software Arlequin v3.5 (Excoffier and Lischer 2010). We estimated the relationship among mtDNA haplotypes throughout the archipelago using a median-joining network (Bandelt et al. 1999), with the haplotype frequencies, using Network, version 5.0.1.1 and post-processing maximum parsimony (Polzin and Daneschmand 2003). The demographic history of each nesting group was investigated using the software DNAsp (Rozas and Librado 2009), first through moment estimates of Tajima's D (computed with 1000 coalescent simulations), measurement of goodness of fit, the raggedness index r (Harpending et al. 1993, Harpending 1994, neutrality tests Fu's Fs (Fu 1997) and Ramos-Onsins and Rozas' R 2 (Ramos-Onsins and Rozas 2002). All were assessed for significance with 1000 coalescent simulations.

Determining feeding strategies
Statistical analyses of stable isotopes (SIA) were performed in R Studio, version 3.2.5 (The R Foundation for Statistical Computing, 2016). We classified individual turtles into putative feeding groups based on δ 13 C and δ 15 N using affinity propagation clustering (apcluster, Bodenhofer et al. 2011). Affinity propagation relies on machine learning and, unlike k-means, can cluster data without pre-assuming the number of clusters. To allow different numbers of clusters to exist across nesting islands, analyses were performed within each nesting group. SIBER (Stable Isotope Bayesian Ellipses in R) was used to visualise the extent of the isotopic niche of feeding groups, identified using affinity propagation clusters, for turtles on each island. Additionally, the standard ellipse areas representing the core isotopic niche for each group, corrected for sample size (SEAc), were calculated.
Turtles were also categorised into size-related clusters where those smaller than 90 cm (curved carapace length, CCL) were considered a priori oceanic and those larger than 90 cm (N = 26, CCL) were considered neritic (Eder et al. 2012;Hawkes et al. 2006;Cardona et al. 2017), independently of their SIA values. To evaluate the reliability of this grouping and test whether size really correlates with feeding ecology, linear discriminant analyses (LDA) with jackknifed classification were used to assign individuals from each island into the size-related (CCL) oceanic or neritic groups based on their δ 15 N and δ 13 C. We then tested whether size-based feeding strategies actually matched SIA-based feeding strategies-hence testing for the reliability of size as a proxy for feeding ecology.
The relationship between δ 15 N and δ 13 C overall in all the datasets was assessed with Pearson's product moment correlation coefficient test. Independent linear mixed-effects models (lme4 and lmertest) were used to explore the relationship between CCL and isotopic cluster, as well as their interaction on δ 15 N and δ 13 C for turtles from each nesting island. This was to test for separate vs. common slopes between the isotope values and CCL across isotope clusters. Possible variation across years was accounted for by including "year" as a random effect. Analyses had to be conducted within island, as not all detected strategies were present on all islands and therefore would bias the evaluation of the interaction between CCL and feeding strategy. Post hoc analyses were performed using multiple comparisons of means of Tukey contrast (multcomp). Furthermore, differences in CCL among isotopic clusters were assessed using independent linear mixed-effects models, with "year" as a random effect. Lastly, Pearson correlations were calculated and tested for significance to investigate possible links between the number of isotopic clusters and demographic indices (Tajima's D, raggedness index r, Fu's F and R 2 ) across islands.

Results
For both SIA and mtDNA sequencing, the Eastern island of Boavista contributed a large proportion of collected samples [Electronic supplementary material (Table S1)]. This reflects the fact that most (~ 60 to 70%) nesting activity of the archipelago occurs on that island  and decreases in a westerly direction. Our sample sizes are therefore mostly representative of nesting densities. Overall, δ 15 N and δ 13 C were positively correlated with one another (cor = 0.140, p = 0.004).

Matching isotopic clusters with feeding strategies
The affinity propagation clustering approach revealed different numbers of clusters across different islands. On Santa Luzia, turtles were grouped into a single feeding strategy (mean ± SD δ 15 N = 14.4 ± 0.6‰ and δ 13 C = − 14.1 ± 0.3‰), corresponding to an oceanic niche. On the islands of Fogo, Maio, Santo Antão and São Vicente, we identified two isotopic clusters where all turtles had depleted 13 C values, also indicative of oceanic feeding strategies (Fig. 1, Table 1): δ 13 C varied between − 16.9 (± 0.8)‰ in the southern island of Fogo and − 14.5 (± 1.5)‰ for turtles nesting on São Vicente. In the oceanic isotopic cluster 1, δ 15 N varied between 10.5 (± 1.0)‰ for turtles on the island of São Vicente and 11.3 (± 2.4)‰ for those on the island of Maio. This pattern contrasts with the second group (isotopic cluster 2) with increased δ 15 N values varying between 13.2 (± 0.9)‰ and 16.8 (± 2.3)‰ (Table 1). Overall, assuming turtles do not feed across a very large number of trophic positions, the 15 N-enriched values of cluster 2 turtles suggest that feeding took place in a 15 N-enriched area.
Turtles nesting on Boavista island were classified into three putative feeding clusters (Fig. 1, Table 1). The two oceanic clusters present on most of the other islands were detected along with a third group. This third group (cluster 3) was relatively 15 N-depleted (mean ± SD δ 15 N = 11.0 ± 1.6‰) as well as significantly 13 C-enriched compared to the other clusters (δ 13 C = -12.0 ± 2.0‰) and therefore formed by turtles using a neritic feeding strategy.

Clusters based on turtle size
We evaluated whether a size-based definition of feeding strategy does match the true SIA-based feeding ecology used by turtles. To this end, we reassigned individuals into our affinity propagated isotopic clusters based on their body size (CCL, Electronic Supplementary Material, Table S2). Based on size, we could reassign correctly a maximum of 46% of the turtles on the island of Fogo and as few as 3.9% for turtles nesting on Santa Luzia to their SIA-based feeding strategy. Even in Boavista, where previous studies have been conducted, the reassignment rate was only 20% overall, notably with only 4 (among 13 larger than 90 cm on that island) individuals reassigned to the neritic group in total. This ultimately demonstrates that size does not represent a sufficiently strong proxy for feeding strategy in this system.

Testing for size by feeding strategy interactions
Having determined the most likely number of feeding clusters representing distinct feeding strategies, we tested if isotope signatures change with CCL, and whether such Table 1 Affinity propagation clusters obtained for nesting female loggerhead turtles by island based on SIA values as well as determined feeding strategy Standard ellipse area corrected (SEAc) for sample size for each cluster is also shown and expressed as ‰. NA appears for the island supporting only one feeding group and hence, SEAc could not be calculated changes were the same in each isotope cluster. Noteworthy, because turtles nesting on Santa Luzia island belong to only one feeding strategy, the analysis could not be conducted.

Linking stable isotopes and demographic parameters
We detected 17 distinct haplotypes in the mtDNA control region (Electronic Supplementary Table S3). The haplotype network mostly reflects that of the species across its entire distribution area in the Atlantic (Shamblin et al. 2014), where two haplotypes are observed (CC-A.1 and CC-A.2), with increased diversity around the CC-A.1 and CC-A17.1 haplotypes (Electronic Supplementary Figure S1, Table S3). In line with Stiebens et al. 2013a, we observed a similar genetic structure in the archipelago with the strongest genetic differentiation observed with turtles nesting on Sao Vicente and Santa Luzia islands (Electronic Supplementary  Table S4). This result reinforces the view of relative independent populations being present across the archipelago (Stiebens et al. 2013a).
Analysing the control region of the mitochondrial DNA, we found that turtles nesting on the island of Boavista showed a significant negative Tajima's D (D = − 1.897, p = 0.003) suggesting a rapid expansion after a bottleneck event. On Boavista, R 2 also showed to be significant (R 2 = 0.031, p = 0.033). Santa Luzia turtles had a significant Fu's F value (Fu's F = − 1.889, p = 0.029), yet this result must be taken with caution given the small sample size for this island. All demographic indices for the other islands were not significant (Electronic supplementary material, Table S1).

Discussion
Maintaining diversity, whether genetic, phenotypic or behavioural is a major goal of modern conservation biology that increasingly focuses on species' adaptive potential (Eizaguirre and Baltazar-Soares 2014). Using stable isotope analyses, we show that several feeding strategies are used by turtles nesting at the Cape Verde archipelago-the only significant nesting aggregation for loggerhead turtles in the Eastern Atlantic ). Contrary to expectations, CCL is not revealed as a strong proxy to define a turtle's feeding strategy (illustrated in Fig. 3). With our data, we also address one of the major previous paradoxes that postulated that the vast majority of sea turtles from Cape Verde used a suboptimal feeding strategy, because oceanic turtles were thought to be smaller than their neritic counterpart. When defining feeding strategies based on stable isotopes, we did not detect this pattern, suggesting that turtles use adaptive ecological niches. Lastly, our results suggest there may be a correlation between diversity of feeding strategy and demography.
Our findings showed that the four islands of Fogo, Maio, Santo Antão and São Vicente support two groups of turtles that diverge in their δ 15 N values. Despite some variation, these two groups had δ 13 C lower than − 14‰, which shows that both groups are formed of oceanic feeders when using previously established criteria for this population (Eder et al. 2012;Pinela et al. 2010, Cardona et al. 2017. The large variation in δ 15 N seen among individuals nesting in the archipelago (~ 15‰) is unlikely to represent variation in trophic position, as such a variation would cover over four different trophic levels using a typical trophic enrichment factor of 3.4‰ (Post 2002)-an improbable scenario for turtles. Instead, this reflects baseline variation in δ 15 N across the vast oceanic area utilised by Cape Verde turtles who, during their foraging time, remain in the open water between Cape Verde and the African continent, covering ~ 177,325 km 2 (Pikesley et al. 2015). Localised oceanic upwelling in part of this area (Hood et al. 2004;Prange and Schulz 2004) results in consumers incorporating a 15 N-enriched label in their tissues (Bergman et al. 2013). Conversely, most pelagic habitats in the region are largely unaffected by upwelling, and turtles feeding in these areas are relatively 15 N depleted compared to their counterparts feeding in upwelling zones. Differences in δ 15 N across feeding grounds have been observed in other loggerhead turtle rookeries (the Pacific and Northwest Atlantic), other sea turtle species such as the leatherback (Seminoff et al. 2012;Tucker et al. 2014;Wallace et al. 2006) and even other marine species such as Atlantic bluefin tuna (Thunnus thynnus, Estrada et al. 2005) and humpback whales (Megaptera novaeangliae, Witteveen et al. 2009). Like in our study, δ 15 N fluctuations in plankton have been linked to upwelling activity in the Eastern Atlantic (Montoya et al. 2002). Average nitrogen fixation rates are up to seven times higher in the surface waters of the upwelling region of the Eastern Equatorial Atlantic because of blooms of non-diazotrophic phytoplankton (Christian and Murtugudde 2003;Subramaniam et al. 2013). Subsequently, this can cause enrichment of 15 N in organisms feeding at higher trophic positions such as sea turtles (Bergman et al. 2013;Goebel et al. 2010;Montoya et al. 2002).
Our data indicate that turtles nesting on Boavista island belong to three isotopic clusters. There, these isotopic clusters differ from one another for both their δ 15 N and δ 13 C values. The distinct locations in δ 15 N-δ 13 C isotope space reveal the presence of the two separate oceanic clusters, seen in the other islands, in addition to a third cluster where turtles show more 13 C-enriched values. This third cluster likely represents individuals utilising a more neritic habitat and more benthic feeding strategy (Cardona et al. 2017;Eder et al. 2012;Hawkes et al. 2006). Turtles from this group are most likely those described, using satellite telemetry, to feed around the coasts of Sierra Leone (Hawkes et al. 2006;Pikesley et al. 2015). Far fewer turtles were reported to make the journey to this neritic ground than the oceanic (12.5% neritic vs. 87.5% oceanic-Eder et al. 2012; Table 2) and neritic feeders utilise a small range (~ 500 km 2 - Pikesley et al. 2015), consistent with the reduced observed δ 15 N isotopic variation seen in this group compared to the oceanic groups. We interpret individuals within these different clusters to differ in their feeding strategies, either differing in both their feeding locations and their prey items (i.e., neritic-oceanic) or only their feeding locations (oceanic-oceanic upwelling).
Variation in the number of feeding clusters among nesting groups raises questions about the factors driving such apparent trophic diversity. In Cape Verde, the highest nesting density occurs on Boavista , followed by the island of Sal (not sampled here), Maio and then the western set of northern and southern islands. While the most distant islands away from Boavista support smaller nesting groups, Fig. 3 Scatterplots illustrating the relationship between stable isotopes δ 15 N and δ 13 C (‰) and the size range of turtles for each affinity propagated isotopic cluster on each island. The black line on each plot represents the 90 cm cutoff between neritic (> 90 cm) and oceanic (< 90 cm) size categories it was demonstrated that those groups harbour unique genetic diversity, which are identified as a conservation priority (Stiebens et al. 2013a). Here, we found that the turtles from Boavista have a significant negative Tajima's D value, suggestive of population expansion after a bottleneck event, supported by a significant R 2 . Our results further suggest a match between population expansion and increased diversity of feeding strategies, linking large population size with increased trait diversity. Noteworthy, our hypothesis that links demography and trait diversity must be treated with the necessary caution, because the occurrence of rare feeding strategies in smaller populations might be overlooked despite exhaustive sampling. As a precautionary conservation measure, however, we encourage increased protection in islands with small nesting groups. With this study, we do not suggest that the control region of the mtDNA mediates variation of feeding strategies and we cannot exclude this proxy to be correlated with any other genetically coded trait. Here, we imply that mtDNA of the sampled females describes the population history of the nesting groups. Future studies may focus on association mapping using more complete genomic information and may provide further insights into the evolution of feeding strategies whether linked to adaptive diversity or better resolution of demography.
Alternatively, the feeding strategy followed by individual turtles from the Cape Verde archipelago likely reflects some environmental determinism. A new paradigm was proposed where upon natal dispersion linked to ocean currents, hatchlings imprint on feeding grounds, allowing them to return to these grounds as adults (Hays et al. 2010). Hence, under this alternative hypothesis, we suggest that the majority of Cape Verde islands are exposed to ocean currents that propel hatchlings ultimately towards oceanic habitats. Because the nesting density on Boavista is high and consequently the nesting season long, there are likely climatic/oceanic conditions (e.g., Monzón-Argüello et al. 2012) that will propel hatchlings towards the neritic habitats. The smaller nesting groups in the West may be less likely to be exposed to those conditions because of their geographic location and shorter nesting period as a result of reduced population size. This would explain why turtles from Boavista island showed an additional neritic feeding strategy compared to the rest of the nesting groups we sampled. Because turtles are mostly faithful to their natal island and evolve local adaptation under philopatry (Stiebens et al. 2013a), this would explain why variation in feeding strategies also appears to be island specific. Overall, independently of the relative contribution of the genetic and environmental determinisms of feeding ecology, conservation efforts should not neglect turtles from any of the nesting groups to guarantee long-term maintenance of behavioural and trait diversities.
Adult female loggerhead turtles from Cape Verde were thought to utilise both productive neritic habitats and putatively suboptimal oceanic habitats due to variable sizes of observed individuals (Eder et al. 2012;Hawkes et al. 2006). Evidence for the multiple oceanic feeding groups shown in this study existed in a previous study focusing solely on Boavista. Eder et al. (2012) actually referred to group 2a and 2b within their oceanic groups, however, sample size likely prevented the statistical detection of two distinct groups. By drastically increasing sample size in our study, and using a non-biased clustering approach, we show that oceanic turtles that feed in a 15 N-enriched environment are larger than their counterparts from both other clusters. Hence, oceanic turtles may not actually be utilising suboptimal habitats, but the presence of variable oceanic conditions and the existence  Rees et al. (2010) of upwelling may have driven the evolution of this feeding strategy, explaining the increased abundance of oceanic turtles in the Cape Verde rookery compared to other locations (Table 2). Consequently, this suggests that turtle size cannot be considered a universal indicator of habitat use in Cape Verde. This is contrary to results reported from other rookeries such as the Japanese Rookery (Hatase et al. 2007(Hatase et al. , 2013, the Northeast Atlantic Rookery (Eder et al. 2012) or the Northwest Atlantic Rookery (McClellan et al. 2010;Reich et al. 2010, Table 2). Noteworthy, there is no correlation between CCL and δ 15 N of loggerhead turtles from the Mediterranean Sea (Zbinden et al. 2011) or Florida (Mansfield et al. 2009), where there is substantial overlap in the sizes of two groups distinguished by δ 13 C value (Reich et al. 2010), further suggesting that the correlation between turtle size and feeding ecology is rookery specific. From a conservation perspective, studies on feeding ecology of endangered species are important, as they allow the indirect identification of habitats that need particular protection. Here, the protection of feeding grounds off the coast of Sierra Leone necessitates local and international efforts similar to those habitats exposed to the upwelling events where fishing rates and turtle bycatch are important. At the local scale of Cape Verde, the different size-feeding ecology relationships are likely associated with the evolution of nesting group-specific traits, matching the rather independent functioning of the nesting groups in Cape Verde (Stiebens et al. 2013a). This strongly supports the perspective that functional trait diversity and demography should not be interpreted as independent of one another.