Sporadic nesting reveals long distance colonisation in the philopatric loggerhead sea turtle (Caretta caretta)

The colonisation of new suitable habitats is crucial for species survival at evolutionary scale under changing environmental conditions. However, colonisation potential may be limited by philopatry that facilitates exploiting successful habitats across generations. We examine the mechanisms of long distance dispersal of the philopatric loggerhead sea turtle (Caretta caretta) by analysing 40 sporadic nesting events in the western Mediterranean. The analysis of a fragment of the mitochondrial DNA and 7 microsatellites of 121 samples from 18 of these nesting events revealed that these nests were colonising events associated with juveniles from distant populations feeding in nearby foraging grounds. Considering the temperature-dependent sex determination of the species, we simulated the effect of the incubation temperature and propagule pressure on a potential colonisation scenario. Our results indicated that colonisation will succeed if warm temperature conditions, already existing in some of the beaches in the area, extend to the whole western Mediterranean. We hypothesize that the sporadic nesting events in developmental foraging grounds may be a mechanism to overcome philopatry limitations thus increasing the dispersal capabilities of the species and the adaptability to changing environments. Sporadic nesting in the western Mediterranean can be viewed as potential new populations in a scenario of rising temperatures.

Philopatry has been found in very different taxa [11][12][13][14] but marine turtles are one of the best examples of this strategy 15 . Early tag studies demonstrated that female turtles return to the same beaches to lay their eggs, sometimes with differences of only a few metres among nesting seasons 16,17 . Posterior genetic studies 18,19 demonstrated that this site fidelity was in fact a true philopatry, as the reused nesting beaches corresponded to the beaches where the nesting females hatched, thus generating a strong female-mediated genetic structuring 20 . Recent studies have shown that males would also show high degree of philopatry 21,22 . The combination of natal imprinting and accurate navigation mechanism has been proposed as the key elements maintaining the philopatry in marine turtles and in other taxa 2 . Newborn hatchlings would memorize different chemical and magnetic cues of the nesting beaches where they hatch and use this information in adulthood to find the natal nesting beaches to reproduce 15,23 .
The potential limitations of philopatry are especially accentuated in marine turtle species, most of them of conservation concern 24 . Furthermore, many researchers have suggested the possible impact of the predicted climate change on marine turtles 25 due to their Temperature Sex Determination (TSD) 26 . Global warming could increase feminisation of the populations and philopatry might prevent dispersal to colder nesting beaches to counteract this effect [27][28][29] . Nonetheless, the circumtropical distribution of most marine turtle species suggests the existence of mechanisms for long distance dispersal 30 . Consequently it has been proposed for marine turtles that 'non-philopatric exploratory behaviours are needed to colonize new nesting environments on evolutionary time scales' 31 , 'strays and wandering must occur, and are no doubt adaptively advantageous aberrations, necessary for colony proliferation' 16 as 'absolute natal homing, over the 100-million-year history of this group, would be a strategy for extinction' 20 .
Tagging has revealed that the typical distance between successive nesting sites of loggerhead turtle (Caretta caretta) individual females is 5 km or less 32 , although interchange of nesting females among more distant localities has also been reported [33][34][35] . In fact, the females nesting in the north-western Atlantic have a remigration rate close to 70% 35 , meaning that a significant portion of the nesting females are not strictly philopatric and lay their clutches in other nesting beaches. These deviations are of tens to hundreds of kilometres from the original nesting beach and could easily explain the lateral spread of the nesting areas along continuous or semi-continuous nesting habitats and thus the existence of Regional Management Units (RMU) 36 . However, these strays are not enough to explain the presence of very distant nesting beaches separated by vast marine areas. Thus, colonisation through long distance dispersal across oceans or seas is the only likely explanation for the current circumtropical distribution of most marine turtle species but these transoceanic dispersal strays have not yet been described in detail.
The Mediterranean Sea offers a unique scenario to ascertain how this long distance dispersal might operate, since the major nesting aggregations of the loggerhead sea turtle in the central and eastern part of this sea are the result of at least two independent colonization events during the late Pleistocene and the Holocene 37 . Furthermore, some sporadic nesting events of this species have been recently detected in the western Mediterranean [38][39][40][41][42] , defined as rare nesting events in an area where low or no nesting activity has been recorded to date. These clutches have always been found in the vicinity of developmental habitats for juveniles of Atlantic and Mediterranean nesting populations [43][44][45][46] and most of them thousands of kilometres away from any known regular nesting area. One possible explanation for these sporadic nesting events is that they are relicts of ancient nesting populations in these locations, as some sporadic nesting had been reported in the past 47,48 . However, another possibility could be that these nests are examples of contemporary long distance dispersal events from distant nesting populations 38 .
The aim of the present study is to contextualize these sporadic nesting events in Caretta caretta under a philopatric scenario and examine the role of sporadic nesting in the long distance dispersal of this species. We thus aim to reconcile a philopatric strategy with a circumtropical distribution and evaluate the possible adaptability of this species under different global warming scenarios through long distance colonisation.

Results
Genetic analysis. Evidence of a total of 40 sporadic nesting events were reported in the western and central Mediterranean in 1870,1990 and from 2001 to 2015 (Table 1, Fig. 1). Almost all nesting events with biometric data (86%, Table 1) produced viable hatchlings, although almost half of the nests with information about the incubation duration (54%, Table 1) suggested a null or low production of females. Furthermore, it was not possible to obtain samples from all the clutches, considering that some were laid before starting the present study, while in other cases no samples for genetic studies were collected by the discoverers or were preserved in formalin. For the genetic analyses we obtained 121 samples (120 hatchlings plus one nesting female) from 18 different clutches in the western Mediterranean (Table 2, Supplementary Data S1). A total of six mitochondrial DNA haplotypes were found among the samples, all of them described in previous studies and found in the GenBank (Table 2). From these six haplotypes, two of them are reported as exclusive from the Atlantic nesting beaches (CC-A1.1 and CC-A9.1) and three of them are common in both Atlantic and Mediterranean nesting beaches (CC-A2.1, CC-A3.1 and CC-A20.1) 49 . The remaining haplotype (CC-A10.4) has only been reported from the nesting population of Melbourne beach (Florida, USA) 50 and from an adult individual foraging in Drini bay (Albania) 51 . However, the short (~380 bp) sequence of this haplotype (CC-A10) had been also found in the nesting population of Zakynthos island (Greece) 52 where no long version of this haplotype has yet been described. Hence, this haplotype may be present as well in the eastern Mediterranean nesting populations. Individual assignments using microsatellites revealed different origins for the samples from different clutches ( Table 2, Figs 1 and 2). Hatchlings from the nests with Atlantic mtDNA haplotypes were assigned to the Atlantic while those from the nest with the CC-A10.4 haplotype were assigned to the Mediterranean. Hatchlings from nests with common mtDNA haplotypes were assigned either to the Atlantic or Mediterranean nesting beaches ( origin. However, we cannot discard that the mixed probabilities found are due to the lack of resolution of the markers as observed in previous studies 46,53 . We collected samples from more than two hatchlings for 8 of the 18 clutches. Multiple paternities were found by GERUD in all but two of these clutches, resulting in a minimum number of 2 fathers per nest ( Table 2). The inferred genotypes of all possible mothers and fathers presented high assignment probabilities to the same nesting area as the offspring. The assignation of the parents were generally better than the offspring (e.g. nest N18, Fig. 2a) probably because all the possible combinations of parents integrated the allele information of all the sampled hatchlings. For instance, rare alleles with high discriminating power may appear only in some hatchlings, and thus only the hatchlings with these alleles will have good assignations, but they will always be present in any parent combination. The nest N13 (Fig. 2a) was an exception as both the hatchlings and the possible parents yielded mixed assignment values resulting in inconclusive mean values. In this case, we separated the putative mothers that were assigned to the Atlantic from those assigned to the Mediterranean and we then reassigned the corresponding pair of fathers inferred by GERUD. If the mother was assigned to the Mediterranean with a good probability, the fathers were assigned to the Atlantic and vice versa (Fig. 2b), thus indicating that the parents had different origin. The relatedness analysis showed that the values obtained among samples within the same clutch were higher than the values obtained among samples of different clutches (Fig. 3a). However, the values obtained within samples of some specific pairs of sporadic nesting events were of the same magnitude than the values obtained within a clutch thus suggesting some level of relatedness. Specifically, the samples of the nests N32 and N33 showed high levels of the relatedness and shared the same maternal haplotype (Fig. 3a, Table 2). Furthermore, the genotype of the mother sampled in Pulpi (nest N33) was compatible with the offspring from the nest in Torrevieja (nest N32) and the two possible fathers inferred for Pulpi nest were also found among the possible father combinations of Torrevieja nest for the genotype of the Pulpi mother. This female, assigned to be of Atlantic origin ( Table 2, Fig. 3a) measured 74 cm of curved carapace length (CCL), slightly below of the mean maturation size of Atlantic turtles visiting the Mediterranean of 80 cm CCL 54 . Considering our results and that the period between the two nests (14 days, Table 1) is within the interesting interval of the species, we hypothesize that these two nests were laid by the same female. The relatedness values among hatchlings of some of the nests of Conigli Beach (Fig. 3b) were generally high (specifically between the nests N8-N9 and N8-N10) suggesting that the same or related females have been laying different clutches in that area.
Effect of incubation temperature and propagule pressure. A model was performed under different scenarios to test the impact of the incubation temperature and propagule pressure rate (Pr) on colonisation success. Colonisation was possible only at temperatures high enough to produce some females in all scenarios (>28 °C) and was generally faster at warmer temperatures, as the percentage of hatchling females produced increases allowing for descendent females to identify the new area as a suitable nesting site and return as adults to lay their eggs (Fig. 4). However, colonisation decreased again at very high temperatures (33 °C), mainly associated to a decrease of the emergence success due to excess of heating 55 . As expected, Pr values impacted the colonisation  Table 1. Stars indicate nesting events with genetic data and reliable assignation to the Atlantic (grey), to the Mediterranean (white) or mixed (black). White dots indicate nesting events with genetic data but with no reliable assignment due to resulting low (<0.80) probabilities or due to not having microsatellite data when published 64  process, by accelerating the time when the first returning females start to contribute to the growing population (Fig. 4). Even at very low Pr values (1 nest per 100 years) colonisation was possible if temperature conditions were appropriate but hundreds of years would be needed to establish a new nesting population (Fig. 4a).

Discussion
Marine turtles have existed for at least 110 million years 20,56 . Considering that all existing species exhibit at least some degree of philopatry, it is reasonable to conclude that it originated at least several tens of millions years ago, before the current marine turtle families diverged 57 . The Earth temperature has experienced drastic oscillations during this period 58,59 , and thus the temporal heterogeneity of the nesting beaches has been necessarily very high. Considering the constraints imposed by TSD and philopatry, it is surprising how this taxon has survived to these oscillations while other lineages went extinct 60 . Furthermore, in almost all marine turtle species with an oceanic stage, philopatry coexists with a circumtropical distribution 61 . As a consequence, long dispersal mechanisms must be present to escape from philopatry constraints and allow the colonisation of distant areas. These mechanisms would explain the apparent adaptability of these species to historical temperature changes, their wide distribution and the existence of widespread common haplotypes 30 . By examining the rare loggerhead sea turtle sporadic nesting events in the western Mediterranean we argue that any delay of the migration from developmental foraging habitats to the adult nesting habitats may result in the establishment of a new population in the vicinity of juvenile foraging grounds, if the conditions of beach temperature are appropriate.

Sporadic nesting events: relicts or colonisers.
Our results suggest that the sporadic nests found in the central and western Mediterranean are not remnants of a past population but the result of an ongoing process of colonisation from distant nesting beaches. Many hatchlings were assigned to Mediterranean or Atlantic populations with high probability values while previous studies have shown a deep genetic differentiation between Atlantic and eastern Mediterranean nesting beaches suggesting a profound isolation of these two distant nesting areas 46 . Furthermore, nests with very different origins were found at close distances, suggesting independent colonisation and not consistent with the hypothesis of an historical population in the western Mediterranean.
The second line of evidence comes from the habitats where these clutches were laid. Beach temperature data derived from temperature recorders deployed at 50 cm below beach surface 62 , our data collected in some of the reported nesting events (Table 1) and simulations from air temperature data 63 indicate that almost all available nesting habitats along the European shore of the western Mediterranean remain below the pivotal temperature and hence are too cold to support a viable population. There are some exceptions, as some of the beaches within   the western Mediterranean could potentially host a nesting population (Table 2) 63 , including southern Spain, southern Italy and the African coast. However, although they may be suitable nesting habitats, they are not of the highest quality 63 as egg viability is predicted to be much lower than at regular nesting sites, as our in situ data collection confirmed (Table 1). Thus the expected production of females would be generally low. Furthermore, temperature conditions can be highly variable among years 62 , and all clutches were, to some extent, manipulated as a consequence of their discovery leading in some cases to artificially increase the proportion of females (Table 1) 64 . Thus, the expected production of females in some areas could be greater in warmer years as estimated at Conigli Beach in Lampedusa, where nests from two different years were obtained: 2006 was warm enough to produce a significant amount of females, while almost only males were produced in 2002. The disparity of nest temperatures among years was also found in another set of records from the Italian coast 65 , with an expected higher production of females on warmer years. As a final line of evidence, beach patrolling was done in all the areas where these sporadic nests were found in the same period of the year and subsequent years, with no female return detected 38 Table 1. the data collected in situ and by the habitat suitability ( Table 1). The repeated nesting found in different years, and the evidence of female returns 66,67 suggests that perhaps these locations hosts a very small and new growing population, being in Lampedusa first reported in 1975 68 and in the Tyrrhenian in 2002 64 . These new populations may be maintained by the females produced in warmer years along with independent recurrent sporadic nesting, probably similar to the nesting events in Calabria 69 or in Sicily island 70 . Our tests of relatedness performed in Conigli nests point in that direction, as some of the nests seem to be related. Finally the clutches from nests N32 and N33 were laid by the same nesting female, within the same nesting season and with a renesting interval normal for the species. Considering all these lines of evidence, we can thus conclude that these sporadic nets are the result of accidental nesting events from individuals originated in distant populations and that mating is possible in foraging areas far away from regular nesting areas. Long distance colonisation mechanism. Considering our results, we propose a long distance colonisation mechanism for the philopatric loggerhead sea turtle. If the individuals undertaking developmental migrations mature before returning to their nesting areas of origin, they would produce sporadic nesting events near the developmental foraging area, as their nesting beaches of origin would be too far to lay the eggs. Under this hypothesis, it is not surprising that the only female that we could sample and measure as a sporadic nester had an Atlantic origin and a size of 74 cm CCL. This size is slightly below the maturation size of individuals of Atlantic origin developed within the Mediterranean 54 , despite being large enough to return to her natal area through the Straits of Gibraltar 71 . Once a sporadic nest is laid, all the females born would be imprinted by the same mechanisms that maintain the philopatry in the species, causing them to return to the new nesting area if they survive to adulthood. This would fix the new nesting location as the place to return in only one generation and without losing philopatry as previously suggested 2 . This is not unprecedented, as this imprinting mechanism has been the basis of reintroduction projects, like the establishment of a new Kemp's Ridley nesting area in Texas (USA) by imprinting hatchlings from Rancho Nuevo (Mexico) nesting area 72 .
Furthermore, the process of colonisation is subject to the production of females and thus it is dependent on temperature. Mean incubation temperatures below 28.3 °C during the middle third of incubation period 73 , would produce no females to return to the new nesting beaches and the colonisation would fail independently of the number of sporadic nests laid in the new area. However, if the temperature increases, the sporadic nesting events, as those observed in the western Mediterranean would produce more females and potentially allow the establishment of new nesting areas under optimal conditions. Consequently, sporadic nesting nearby developmental feeding areas could act as a dormant mechanism of colonisation that activates when the environmental conditions change. Propagule pressure was long identified as the major single parameter predicting the success of biotic introductions 74 , and hence isolated nesting events will hardly result into the establishment of new populations of loggerhead turtles. Our results showed that low colonisation rates could potentially establish a new population but hundreds of years would be needed, while faster rates would produce the same effect in only a few decades. Thus, all factors affecting the propagule pressure would have an impact on the probability of colonisation success. For instance, it has been hypothesized that warming of the regular nesting areas would imply a significant population growth through the production of more females, thus increasing the number of females and thus the rate of potential colonisers 75 . A reduced number of sporadic nests is likely to produce a bottleneck in the potential new population leading to high levels of inbreeding. However, the multiple paternity found in all but two of the nests, in one case even involving parents from very distant nesting areas (nest N13, and among Conigli beach nests), could help to increase the overall genetic variability, reducing the bottleneck effect and thus increasing the potential success of the colonisation process.

Colonisation of the Mediterranean. The Mediterranean loggerhead nesting populations were originated from
Atlantic individuals during the Pleistocene before the last glacial maximum (20,000-200,000 years ago) and some of the nesting populations have been suggested to be older than others 37,76 . In some cases, the genetic structure among populations have been explained by recolonisations within the basin as the species distribution range retracted during the colder phases of the Pleistocene to warmer refugia (such as Libya, Greece and Turkey) and expanded again within the eastern Mediterranean, when the thermal conditions changed. However, this hypothesis hardly explains the genetic structuring of some of the nesting areas in the central Mediterranean 69 . Previous studies concluded that the high levels of genetic diversity and the presence of one haplotype (CC-A20.1) shared with Atlantic nesting populations in the nesting population of Calabria (Italy) results from an independent colonisation event from Atlantic individuals that migrated into the Mediterranean during the Holocene (<10,000 years) 37,69 . This has been explained by the fact that central Mediterranean has mean temperatures colder than the eastern basin but warmer than the western basin 63 . Considering our results, we suggest that the general warming produced during the Holocene could have activated the colonisation process of the central Mediterranean area from sporadic nesting events close to developmental feeding grounds. A similar process could be happening now in Conigli Beach in Lampedusa island, having higher temperatures than all the other reported sporadic nesting sites and probably some of the nests reported already are females born in these areas that are returning after reaching maturity. If Conigli Beach is really an ongoing colonisation, the Atlantic haplotype CC-A9.1 would be the most recent acquired haplotype of the Mediterranean nesting populations.  Table 3. Summary of the demographic parameters used in the present study. Six different temperatures during the middle third of the incubation were considered. Some of the parameters were assumed to be constant in the range of temperatures considered. *First value for turtles of Mediterranean origin and second value for turtles of Atlantic origin. We thus suggest the hypothesis of a sequential colonisation of the Mediterranean; the eastern Mediterranean was colonized before the last glacial maximum 37 , the Calabria (Southern Italy) was colonized during the Holocene 37 and the central and western Mediterranean may be an ongoing colonisation process (present study), depending on the environmental conditions including global warming 28,63,77 . Future monitoring of the western Mediterranean potential nesting beaches along the next decades would clarify if we are really facing this ongoing colonisation process and how the origin of the new colonizers impacts on the whole basin.
Currents in the Straits of Gibraltar and the southwestern Mediterranean could play an important role in this potential new colonisation as they trap loggerhead turtles of Atlantic origin within the Mediterranean for long periods, thus probably increasing the chances that they lay a sporadic nest before they return to their beaches of origin 71 . These currents are reinforced during cold periods 78 , but if its strength lowers under warmer temperatures, the Atlantic individuals would not be trapped anymore within the western Mediterranan, thus lowering the colonisation pressure in this area. Furthermore, stochastic factors are expected to play an important role in the colonisation process, not in vain we have only detected a minimum of two independent colonisation processes in the Mediterranean in the last 65,000 years 37,69 . For instance, sea temperature may have an impact in hatchling survival, as it is predicted that the temperature at sea is not going to rise in the same proportion that at the beaches 64 . Thus, low winter sea surface temperatures may limit the survival and self recruitment of hatchlings.
Conservation implications. The mechanism for long distance expansion of nesting habitat proposed here would be an alternative method for the species to expand its distribution to new suitable distant habitats, as the old nesting habitats become suboptimal. For instance, a warming of Mediterranean sea surface temperatures from east to west has been shown, starting in 1970 and continuing in this century 28 . Considering our results, a general warming would favour the colonisation of the western Mediterranean during the following decades 28 , while the eastern Mediterranean nesting areas would become too hot thus decreasing hatchlings survival as predicted in previous studies 77 . The combination of contiguous range expansion and long distance dispersal mechanisms has proved to be effective in an evolutionary scale, as these animals have survived the drastic climate changes of the last 110 million years. However, whether or not these mechanisms would be fast enough to counterbalance some of the effects of the current climate change is something that remains to be tested. In any case, several potential conservation and research actions may be considered in the light of our results. The first one would be the monitoring of the effect of warming temperatures on sex ratio, inbreeding and nestling survival on current nesting beaches. The detection of sporadic nesting events through extensive monitoring of potential suitable habitats, coupled with protection and conservation of newly colonised sites, might facilitate the possible expansion and long term survival of the species. Furthermore, assisted migration through egg translocation might be an effective action to promote the creation of new populations in more suitable habitats 72 . The detection of these rare events through extensive monitoring on potential suitable habitats, coupled with its protection and conservation may be crucial to facilitate the possible expansion and long term survival of the species.

Methods
Genetic analysis. Loggerhead turtle sporadic nesting events are very rare and thus their sampling is necessarily opportunistic (Table 1, Fig. 1). When a nesting event occurred, we collected the basic biometric data of the nest and at least one sample per clutch from a dead hatchling and/or embryo found after emergence. When possible, several hatchlings were sampled in order to test for multiple paternities. The collection of samples was conducted in strict accordance with Spanish and European regulations. The Ethics Comitttee of Animal Experimentation of the University of Barcelona stated that the analysed procedures fit the essential ethical rules and the legislation in force according to the article RD2013 of B.O.E from 8th of February 2013. This sort of studies is excluded from the area of application of the legislation, and therefore, the corresponding authorization of this Ethics Comitee is not needed. Furthermore, although the study species is listed in CITES, transportation of samples within the European Union does not require CITES permits. Muscle or skin samples of 120 hatchlings and one female were collected from 18 clutches ( Table 2) and stored in 95% ethanol. DNA was extracted using the QIAamp extraction kit (QIAGEN ® ) following the manufacturer's instructions.
Individual assignments of all samples to either Atlantic or Mediterranean nesting beaches were done using a combination of a fragment of the mtDNA control region and seven microsatellites 46,71 . As a first step, we amplified a long (~800 bp) fragment of the mitochondrial DNA (mtDNA) control region of one sample per clutch, using the primers LCM15382 and H950 79 as it has been proven to be much more informative 49 than the shorter fragment (~380 bp) of the same region used in previous studies 44,80 . We used the same PCR conditions as in previous studies (e.g. 37 ). Sequences were aligned using BioEdit v7.1.11 81 and compared to known loggerhead haplotypes found in the database maintained by the Archie Carr Centre for Sea Turtle Research (http://accstr.ufl.edu/). Published haplotype frequencies on nesting populations 37,49,50,[82][83][84] were compared with our samples and the origin of the nesting female was directly determined if it had an exclusive haplotype. Additionally, we genotyped seven microsatellites of all the samples using primers previously used for this species: Cm84, Cc117, Cm72 and Ei8 85 ; Cc141 and Cc7 86 ; and a modified version 80 of Ccar176 87 using the same protocols described in the literature 80 . The combination of the seven microsatellites was checked against the baseline of Atlantic and Mediterranean individuals used in a previous study to perform individual assignments 46 using the program STRUCTURE v 2.3.4 88 . This baseline comprised a total of 112 individuals of Mediterranean origin and 56 individuals of Atlantic origin (Supplementary Dataset S1). A previous study 46 indicated that the best number of clusters of this baseline was K = 2 and that Atlantic and Mediterranean samples were highly differentiated (F ST = 0.029, P < 0.001). Five samples of this baseline were genotyped again for the present study in order to check for changes in allele sizing and thus all microsatellite data was compatible after a correction for allele size changes (data not shown). An input file was prepared including both the baseline and the sporadic nests genotypes. All the samples of the baseline Scientific RePoRTS | (2018) 8:1435 | DOI:10.1038/s41598-018-19887-w were labelled as belonging either to the Mediterranean (1) or to the Atlantic (2) while all sporadic nests samples were labelled as belonging to an additional group (3). Only samples from the baseline were used as locprior and thus no prior was assumed for the sporadic nests samples. STRUCTURE was run under the assumption of no admixture, considering the differentiation between the Atlantic and Mediterranean samples of the baseline 46 , performing 1,000,000 repetitions after a 100,000 burn-in period and assuming K = 2. Each assignation was replicated 10 times and the results were averaged using CLUMPP v1.2 69 . Hence, a probability of being either Atlantic or Mediterranean was ascribed to each sample. Samples with assignment values higher than 0.8 to one of the two groups were assumed to be reliable as previous studies recommended 46 .
Multiple paternity was tested for all clutches sampled for more than two hatchlings, using the microsatellite data and the program GERUD v 2.0. This software allows the reconstruction of parental genotypes from half-sib progeny with unknown parents and infers multiple paternity by identifying more than four different parental alleles from a clucth 90 . Additionally, we generated the genotypes of all possible mothers and fathers of our offspring and we made individual assignments of up to the 100 most probable parents using the same methodology described above. Furthermore, we calculated the Lynch & Ritland relatedness index 91 between all hatchling pairs using GenAlex v 6.5 92 . Values were multiplied by two as indicated in the program to allow the index to vary between −1 and 1. Relatedness values among all sample pairs were used to create a heatmap using 'ggplots2' 93 . Effect of incubation temperature and propagule pressure. The population dynamics in a colonisation process under ideal conditions was modelled using a simplification of the population models developed in the literature [94][95][96] and implemented on an Excel TM spreadsheet in order to assess the impact of the temperature at incubation and propagule pressure. We defined the propagule pressure rate (Pr) as the number of sporadic clutches laid in a new area per year by individuals from distant populations (colonisers). Considering the number and frequency of the sporadic nests found in the western Mediterranean, we assumed that accidental colonisers lay a single nest before returning to their original nesting area and that the colonisation rate is constant in the whole area. These sporadic colonisers are expected to produce in a year (t) a number of hatchling females (Fc (t) ) equivalent to the mean annual number of eggs laid per nest (En), the emergence success (Es) and the expected percentage of females produced per nest (Of) as in equation 1.

Fc
EnEsOf The females born in the new nesting area will potentially return in the future as adults to lay their own nests as a result of philopatry through natal imprinting as in equation 2. Thus the annual number of hatchling females produced by the females originally hatched in the new nesting area (Fp (t) ) would be: Which is basically the same as equation 1 but considering N (t) as the number of reproductive females in a given year (t) originated in the new population and Cy the mean number of clutches per female laid each year. As a consequence, the total number of female hatchlings produced at any year (t) in the new nesting area (F (t) ) is the sum of the hatchling females produced by the colonisers and those produced by the reproductive females originally hatched in the new population (Equation 3).

= +
At the start of the colonisation process, N (t) value is zero (and so is Fp (t) ), as no resident females are found in the new area. Thus, the production of hatchling females in the area would be caused only by the sporadic colonisers until some of the females born in the area settle back there to reproduce. This will happen when the newborn hatchling females reach maturity after Ma years. We thus consider the number of recruited females in any time (R (t) ) as the result of the female hatchlings produced Ma years before in the area and that survived to reproduce as adults (Equation 4).
where Sm is the survival rate from hatch to maturity. These recruits will continue to reproduce in the area as a consequence of the philopatry during the following years. Thus, the nesting female population size (N (t) ) at any moment, measured as the number of reproductively active females in the population, is the sum of all the females recruited up to Tr years ago that survived to this moment to reproduce (Equation 5).
where T r is the mean duration of the reproductive period (years) and Sa is the annual survival for adults. As a consequence, the temporal changes in population size from the start of the colonisation can be obtained. In order to test the colonisation process under different scenarios, we used a different combination of parameters that are dependent on the temperature of the middle third of incubation (Table 3). A recent study showed that these reproductive parameters are likely to be independent of population size 97 so we assumed that they would be similar in a colonising scenario. Furthermore, four different propagule pressure rates (Pr) were modelled: 1 sporadic nest per a hundred years (Pr = 0.01), 1 sporadic nest per ten years (Pr = 0.1) and 1 sporadic nest per year (Pr = 1) and 10 sporadic nests per year (Pr = 10). Finally, the upper threshold for population size is supposed to be produced by limiting factors related to the nesting area (probably density dependent). As population sizes range from 20 to 687 females in Mediterranean nesting rookeries 98 , we considered these values as a rough guide for what would be normal population sizes in the Mediterranean after a colonisation.