Ecology

Suborder Glossiphoniiformes Tessler & de Carle, 2018Family Glossiphoniidae Vaillant, 1890Comments. Our two-locus phylogeny reveals the presence of two large clades, corresponding to the subfamilies Glossiphoniinae and Haementeriinae (Fig. 1). The subfamily Theromyzinae Sawyer, 1986, delineated by some authors2,10,22, was not supported as a distant phylogenetic clade and their representatives are clustered within the monophyletic Glossiphoniinae. The same pattern was recovered by earlier phylogenetic reconstructions3,30,33,34. These data indicate that Theromyzinae may represent a synonym of the latter subfamily. However, a subfamily-level revision of the Glossiphoniidae is beyond the framework of the present study.Subfamily Glossiphoniinae Vaillant, 1890Genus Alboglossiphonia Lukin, 1976Type species: Alboglossiphonia heteroclita (Linnaeus, 1761) (= Hirudo heteroclita Linnaeus, 1761; by original designation).Arctic occurrences. Our results reveal that members of this genus are not common inhabitants of the Arctic but two species, A. heteroclita (Linnaeus, 1761) and A. sibirica sp. nov., cross the Arctic Circle on the Yamal Peninsula through the Ob and Taz rivers (Table 1). Previously, it was shown that A. heteroclita occurs in the lower Ob Basin, northern edge of Western Siberia23.Comments. This genus contains inconspicuous minute leeches and is characterized by a nearly global distribution1. It definitely requires an integrative taxonomic revision. Available genetic evidence (Fig. 1 and Supplementary Fig. S1) reveals that the North American populations of what was traditionally assigned to A. heteroclita should be considered a separate species, A. pallida (Verrill, 1872) (type locality: West River near New Haven, Connecticut, USA)35,36. Other species, which occurs in Siberia and the Far East, was tentatively assigned to Alboglossiphonia cf. papillosa (Braun, 1805) based on a darker pigmentation of its dorsum37,38 but it represents a separate North Asian species, which is described here.
Alboglossiphonia sibirica Bolotov, Eliseeva, Klass & Kondakov sp. nov = Alboglossiphonia heteroclita Lukin (1957): 27339 (identification error). = Alboglossiphonia heteroclita papillosa Kaygorodova et al. (2014): 337; Kaygorodova (2015): 4140 (identification error). = Alboglossiphonia cf. papillosa Klass et al. (2018): 2638 (identification error).Figures 4a, 5a, 7a, Supplementary Figs. S2a, S3a, S4, Supplementary Table S2.LSID: https://zoobank.org/urn:lsid:zoobank.org:act:19B581C3-E912-487C-B9EC-8E50DDEFD380.Holotype. RMBH Hir_0542_2-H (non-sequenced), RUSSIA: Lake Torfyanka, 43.0761° N, 131.9620° E, Vladivostok, Primorye, August 12, 2021, Y. E. Chapurina leg.Paratypes (N = 13). RUSSIA: 1 specimen RMBH Hir_0542_2 (sequenced: COI sequence acc. No. ON873332), the type locality, the same date, and collector; 1 specimen RMBH Hir_0396 (non-sequenced), an oxbow lake of Taz River, near Tazovsky settlement, 67.5063° N, 78.6751° E, Yamal-Nenets Region, August 22, 2019, E. S. Babushkin leg.; 1 specimen RMBH Hir_0394 (DNA voucher; sequenced: COI sequence acc. No. ON548508), Vitim River, 57.2010° N, 116.4300° E, Lena River basin, Vitimsky Nature Reserve, Irkutsk Region, July 12, 2019, E. S. Babushkin leg.; 4 specimens RMBH Hir_0013 (3 sequenced with DNA vouchers and one placed on 36 permanent slides as a series of slices; COI sequence acc. No. MH286267, MH286268, and MH286269; 18S rRNA sequence acc. No. MH286273), between zooids of a bryozoan colony, small floodplain lake of the Lena River near Yakutsk, 62.3076° N, 129.8999° E, Yakutia Republic, August 20, 2017, I. N. Bolotov leg.; 1 specimen RMBH Hir_0417_2 (DNA voucher; sequenced: COI sequence acc. No. ON548511), Oron Lake, Gnilaya Kurya Bay, 57.1750° N, 116.4031° E, Lena River basin, Vitimsky Nature Reserve, Irkutsk Region, July 1, 2019, E. S. Babushkin leg.; 1 specimen RMBH Hir_0409_1 (sequenced: COI sequence acc. No. ON548509), a roadside ditch in Knevichi settlement, 43.3886° N, 132.1880° E, Primorye, September 10, 2020, O. V. Aksenova et al. leg.; 1 specimen RMBH Hir_0413 (sequenced: COI sequence acc. No. ON548510), a puddle near railway at Artem city, 43.3794° N, 132.2188° E, Primorye, September 10, 2020, O. V. Aksenova et al. leg.; 1 specimen RMBH Hir_0003_3 (DNA voucher; sequenced: COI sequence acc. No. MN393256), Tumnin River, 49.9451° N, 139.9181° E, Khabarovsk Region, July 14, 2014, I. N. Bolotov & I. V. Vikhrev leg.; 1 specimen RMBH Hir_0509_1 (sequenced: COI sequence acc. No. ON548516), a reservoir on the Bolshoy Alim River, near Tolstovka settlement, 50.1981° N, 127.9431° E, Amur Region, July 3, 2021, O. V. Aksenova et al. leg.; 1 specimen RMBH Hir_0510_1 (DNA voucher; sequenced: COI sequence acc. No. ON548517), an oxbow lake of Bureya River, near Novospassk, 49.6756° N, 129.7343° E, Amur Region, July 3, 2021, O. V. Aksenova et al. leg.Etymology. The name of this species reflects its broad distribution in Siberia.Differential diagnosis. Small leech, which could be distinguished from other congeners by a combination of the following characters: dorsum covered by numerous small, shallow, and indistinct papillae, light yellow, with multiple dark spots and short dashes arranged to 18–24 longitudinal rows; these spots and dashes merged into longitudinal lines in the anterior half of the body (the dark markings pattern often lost in ethanol-preserved animals). Externally, the new species is similar to A. heteroclita, A. hyalina (O. F. Müller, 1773), and A. striata (Apáthy, 1888). However, all these species do not have numerous dark spots and short dashes arranged to multiple longitudinal rows. Additionally, A. heteroclita differs from the new species by having a median row of segmentally arranged dark spots and a smooth dorsum without papillae. A. hyalina differs from A. sibirica sp. nov. by the general lack of dark pigmentation. A. striata differs from the new species by having a median row of segmentally arranged dark transverse stripes and a smooth dorsum without papillae.Molecular diagnosis. The new species represents a separate genetic lineage but is more closely related to A. heteroclita (mean pairwise COI p-distance = 5.1%; range 4.9–5.4%). The intraspecific pairwise COI p-distance ranges from 0.0 to 2.1% (mean ± s.e.m. = 1.31 ± 0.10%; N = 14 sequences and 91 pairwise distance values). The GenBank acc. numbers of reference DNA sequences (COI and 18S rRNA) are given in Supplementary Table S2 and Supplementary Datasets S1–S2.Description. Small leech (body length up to 11.9 mm). Measurements of the type series are given in Supplementary Table S2. Body broad, leaf-like, ovate. Dorsum with numerous small, shallow, and indistinct papillae. Posterior sucker small, circular (maximum diameter of 2.25 mm), ventrally directed. Proboscis pore in the center of anterior sucker. Coloration of living animals: body dirty yellow with multiple brown spots and dashes arranged to longitudinal rows; in the anterior half of the body, these spots and dashes merged into longitudinal lines. Coloration of ethanol-preserved animals: body light yellow; dorsum with multiple dark spots and short dashes arranged to 18–24 longitudinal rows; these spots and dashes merged into longitudinal lines in the anterior half of the body but the dark markings pattern often lost due to preservation. Three pairs of eyespots; the eyespots of the distal pair joined into a single spot; the eyespots of the next two pairs are spaced apart and fused together. Venter light yellow or whitish. Total number of annuli: 70. Somites I–IV joined to form a head region, somites V–XXIV triannulate, somites XXV–XXVII uniannulate. Gonopores joined and open in the furrow XIIa1/a2. Reproductive system: 6 pairs of large, bag-like testisacs inter-segmentally from XIII/XIV to XIX/XX; atrium small, spherical, the atrial cornua twisted anteriorly; paired ejaculatory ducts twisted, short; paired ovisacs narrow, very short. Digestive system: proboscis sheath massive, long, thick; salivary glands diffuse; crop with 6 pairs of crop caeca: 1st-5th uniform, bag-like, 6th pair (posterior caeca) with 3 blind processes; intestine with 4 pairs of rather short processes and an ovate extention after the last pair of processes.Distribution. North Asia: Western Siberia, Eastern Siberia, the Russian Far East, and Mongolia39.Habitats and ecology. This species is known to occur in a broad range of freshwater environments such as rivers, oxbow lakes, large to small natural lakes, reservoirs, road ditches, and even puddles (Supplementary Dataset S2). An unusual example of its association with a bryozoan species was described from Eastern Siberia38. Probably, the record of an Alboglossiphonia leech in the mantle cavity of an unidentified lymnaeid snail from the Altai Mountains, Russia41 could also be attributed to this species. The life cycle of the new species is unknown.Genus Glossiphonia Johnson, 1816Type species: Glossiphonia complanata (Linnaeus, 1758) (= Hirudo complanata Linnaeus, 1758; by subsequent designation).Arctic occurrences. Representatives of this genus are the most remarkable component of the Arctic Glossiphoniidae fauna. Altogether seven species were recorded north of the Arctic Circle, two of which are new to science and described here (Table 1).Comments. In general, sequenced representatives of the genus Glossiphonia could phylogenetically be delineated to three species groups (or subgenera): (1) the complanata-group (= subgenus Glossiphonia s. str.); (2) the verrucata-group (= subgenus Boreobdella Johansson, 1929; type species: Clepsine verrucata Müller, 1844); and (3) the concolor-group (= subgenus Paratorix Lukin & Epstein, 1960; type species: Torix baicalensis Stschegolew, 1922) (Table 1, Fig. 1 and Supplementary Fig. S1).
Glossiphonia arctica Bolotov, Eliseeva, Klass & Kondakov sp. novFigures 4B, 5b,c, 7c, Supplementary Figs. S2b, S3b, S5, Supplementary Table S2.LSID: https://zoobank.org/urn:lsid:zoobank.org:act:FADF0993-A946-413A-9680-25BA0F9BE90D.Holotype. RMBH Hir_0457_2_1-H (sequenced: COI sequence acc. No. ON810735; 18S rRNA sequence acc. No. ON819028), RUSSIA: a large lake near Sob’ railway station, 67.0480°N, 65.6316°E, Polar Urals, June 23, 2021, A. V. Kondakov et al. leg.Paratypes (N = 18). 18 specimens RMBH Hir_0457 (two specimens sequenced: COI sequence acc. No. ON810736 and ON810737; 18S rRNA sequence acc. No. ON819029; one specimen placed on 18 permanent slides as a series of slices), the type locality, the same date, and collectors.Etymology. The name of the new species indicates that its type locality is situated in the Arctic Region.Differential diagnosis. Medium-sized leech, which could be distinguished from other congeners by a combination of the following characters: dorsum with four rows of ovate, broad but very shallow and indistinct papillae on annulus a2 (outer paramedian and inner paramarginal series); each papilla bears ovate light yellow or white spot; dorsal black markings pattern absent. Externally, the new species is similar to G. mollissima. However, the latter species differs from G. arctica sp. nov. by having larger papillae and a well-developed black markings pattern dorsally.Molecular diagnosis. The new species represents a separate genetic lineage belonging to the verrucata-group (Fig. 1). The pairwise COI p-distance of the new species from other congeners varies from 7.0 to 12.4%. The intraspecific pairwise COI p-distance ranges from 0.0 to 0.2% (mean ± s.e.m. = 0.10 ± 0.05%; N = 3 sequences and 3 pairwise distance values). The GenBank acc. numbers of reference DNA sequences (COI and 18S rRNA) are given in Supplementary Table S2 and Supplementary Datasets S1–S2.Description. Medium-sized leech (body length up to 13.3 mm). Measurements of the type series are given in Supplementary Table S2. Body broad, leaf-like, ovate. Dorsum with four rows of ovate, broad but very shallow and indistinct papillae on annulus a2 (outer paramedian and inner paramarginal series). Posterior sucker small, circular (maximum diameter of 1.9 mm), ventrally directed. Proboscis pore in the center of anterior sucker. Coloration of living animals: body almost transparent, light brown, with multiple yellowish pigment cells. Coloration of ethanol-preserved animals: dorsum beige to light brown, with darker broad inner paramedian lines and light yellowish areas laterally and anteriorly; ovate light yellow or white spots at each papillae on annulus a2 arranged into four longitudinal rows (outer paramedian and inner paramarginal), sometimes with a few white spots between them. Three pairs of ovate eyespots arranged to two parallel rows; in some specimens eyes on each side are joined to a single large spot. Venter whitish to light brown, sometimes with irregular brownish shading. Total number of annuli: 70. Somites I–III uniannulate, IV biannulate, V–XXIV triannulate, XXV biannulate, XXVII uniannulate. The male and female genital pores are separated by two annuli and are located in the furrows XIa3/XIIa1 and XIIa2/a3, respectively. Reproductive system: 6 pairs of spherical testisacs inter-segmentally from XIII/XIV to XVIII/XIX; atrium spherical, the atrial cornua large, twisted anteriorly; paired ejaculatory ducts very long, extending to XVIII; paired ovisacs massive, long, with multiple lobes, arranged as loops, extending to XVIII (pregnant specimen with eggs). Digestive system: proboscis sheath massive, thick, elongated; esophagus narrow; salivary glands diffuse; crop with 7 pairs of crop caeca: 1st-6th uniform, bag-like, 7th pair (posterior caeca) with 4 blind processes and several smaller lobes; intestine enlarged, with 4 pairs of large, long, bag-like processes, expanding distally, each with several short lobes; a large circular extension after the last pair of processes.Distribution. Polar Urals (not known beyond the type locality).Habitats and ecology. The type series of this species was collected from a natural mountain lake with stony bottom. The leeches were recorded beneath flat stones (Fig. 3b); their feeding behavior and life cycle remain unknown.
Glossiphonia taymyrensis Bolotov, Eliseeva, Klass & Kondakov sp. novFigures 4E, 5d, 7b, Supplementary Figs. S2h, S3c, S6, Supplementary Table S2.LSID: https://zoobank.org/urn:lsid:zoobank.org:act:40269BF4-FE1C-4269-A7CC-41020789DC44.Holotype. RMBH Hir_0258_1-H (sequenced: COI sequence acc. No. ON810695), RUSSIA: small lake near Dudinka on Taymyr Peninsula, 69.4008°N, 86.3384°E, July 16, 2018, O. V. Aksenova et al. leg.Paratypes (N = 8). RUSSIA: 2 specimens RMBH Hir_0263_1 and RMBH Hir_0264_3 (sequenced: COI sequence acc. No. ON810701 and ON810705; 18S rRNA sequence acc. No. ON819017), the type locality, the same date, and collectors; 2 specimens RMBH Hir_0256_1 (one sequenced and one placed on 20 permanent slides as a series of slices; COI sequence acc. No. ON810693), small lake near Dudinka on Taymyr Peninsula, 69.3987° N, 86.3505° E, July 16, 2018, O. V. Aksenova et al. leg.; 1 specimen RMBH Hir_0261_2 (sequenced: COI sequence acc. No. ON810699; 18S rRNA sequence acc. No. ON819016), small lake near Dudinka on Taymyr Peninsula, 69.4014° N, 86.3250° E, July 16, 2018, O. V. Aksenova et al. leg.; 1 specimen RMBH Hir_0265_2 (sequenced: COI sequence acc. No. ON810706; 18S rRNA sequence acc. No. ON819021), Bolgokhtokh River near Dudinka, Taymyr Peninsula, 69.3780° N, 87.2215° E, July 21, 2018, O. V. Aksenova et al. leg.; 1 specimen RMBH Hir_0488 (sequenced: COI sequence acc. No. ON810755), a lake on Putorana Plateau, 68.7607° N, 91.9014° E, July, 2021, E. S. Chertoprud leg.; 1 specimen RMBH Hir_0449 (sequenced: COI sequence acc. No. ON810731), Pyzas River near Ust-Kabyrza settlement, 52.8277° N, 88.3973° E, Tashtagolsky District, Kemerovo Region, July 23, 2020, E. S. Babushkin & M. V. Vinarski leg.Etymology. The new species is named after the Taymyr Peninsula, where the majority of the type specimens were collected.Differential diagnosis. Small leech with broad, leaf-like, ovate body; three pairs of eyespots (distal pair joined; next two pairs separate); dorsal papillae absent; dorsal coloration with two inner paramedian rows of black spots, sometimes joining into unclear dashed lines; two annuli between the male (XIa3/XIIa1) and female (XIIa2/a3) genital pores. The new species largely resembles G. complanata but could be distinguished from it by having a smooth dorsum, without clear papillae. These taxa seem to have non-overlapping, allopatric ranges and, hence, could be separated on the basis of geographic criteria. However, the DNA approach seems to be the most appropriate way to distinguish these two species.Molecular diagnosis. The new species represents a separate genetic lineage belonging to the complanata-group (Fig. 1). The pairwise COI p-distance of the new species from other congeners varies from 6.0 to 12.2%. The intraspecific pairwise COI p-distance ranges from 0.0 to 1.1% (mean ± s.e.m. = 0.52 ± 0.07%; N = 8 sequences and 28 pairwise distance values). The GenBank acc. numbers of reference DNA sequences (COI and 18S rRNA) are given in Supplementary Table S2 and Supplementary Datasets S1–S2.Description. Small leech (body length up to 11.3 mm). Measurements of the type series are given in Supplementary Table S2. Body broad, leaf-like, ovate. Dorsum smooth, without clear papillae. Posterior sucker ovate (maximum diameter of 3.0 mm), ventrally directed. Proboscis pore in the center of anterior sucker. Coloration of living animals: not examined. Coloration of ethanol-preserved animals: (1) typical form having beige to light brown ground color without light spots but with darker brown coloration between inner paramedian lines; (2) melanic forms having dark brown ground color with four rows of large yellow spots (outer paramedian and marginal series) and yellow median stripe anteriorly (f. ‘maculosa’) or with strongly reduced yellow markings pattern. In all forms, there are two inner paramedian rows of black spots, sometimes joining into unclear dashed lines. Three pairs of ovate eyespots; the eyespots of the distal pair joined into a single spot; the eyespots of the next two pairs separate and are spaced apart. In the typical form, venter light yellow, with paired brown median and outer paramedian lines, which may be reduced to series of narrow brown longitudinal stripes. In melanic forms, ventral markings is more developed, with a series of brown longitudinal lines from median to inner paramarginal position and outer paramarginal brown spots. Posterior sucker with dense brown spots in melanic forms and with scarce brown spots in typical form. Total number of annuli: 68. Somites I–IV uniannulate, V–XXIV triannulate, XXV biannulate, XXVI–XXVII uniannulate. The male and female genital pores are separated by two annuli and are located in the furrows XIa3/XIIa1 and XIIa2/a3, respectively. Reproductive system: 6 pairs of spherical testisacs inter-segmentally from XII/XIII to XVIII/XIX; atrium ovate, the atrial cornua directed laterally; paired ejaculatory ducts twisted, short; paired ovisacs short, thick (undeveloped). Digestive system: salivary glands diffuse; proboscis sheath moderately thick; esophagus ovate; crop with 6 pairs of massive, bag-like, uniform crop caeca; intestine with 4 pairs of processes.Distribution. Western and Eastern Siberia.Habitats and ecology. The new species was recorded from natural lakes and rivers (Supplementary Dataset S2); its feeding behavior and life cycle are unknown.Genus Hyperboreomyzon Bolotov, Eliseeva, Klass & Kondakov gen. novLSID: https://zoobank.org/urn:lsid:zoobank.org:act:298FF41E-AF0D-4442-9F82-3022B8094A67.Type species: Hyperboreomyzon polaris gen. & sp. nov.Etymology. This name is compiled using two Greek words: ‘Hyperborea’ (meaning a mythical far northern land) and ‘myzon’ (meaning sucking).Diagnosis. Medium-sized, elongate, sub-fusiform glossiphoniid leeches; body and posterior sucker densely covered by shallow, ‘fish-scale’-like papillae; somite V biannulate; somites XII–XXIII secondarily sexannulate dorsally and ventrally due to the presence of very deep, prominent furrows separating each annulus to two semi-annuli; six rows of prominent dorsal tubercle-like papillae at a2 (inner paramedian, inner paramarginal, and marginal series) from V to XXVI; two pairs of circular eyespots on II and Va1 at inner paramedian position; gonopores at the furrows XIa3/XIIa1 (male) and XIIa2/a3 (female) and separated by two annuli; male atrium spherical; proboscis pore opens in a thick velar fold in the anterior half of oral sucker; one pair of compact, massive, elongated, incurved salivary glands, each gland with a bunch of a few short processes apically; 9 crop caeca pairs. Comparison of the new genus with other genera in the family based on morphological and anatomical features is presented in
Supplementary Table S3. Sexannulate condition was also recorded in the genus Actinobdella Moore, 1901 from North America36, but it differs from Hyperboreomyzon gen. nov. by having one pair of eyespots, diffuse salivary glands, and an apical position of proboscis pore (Supplementary Table S3).Comments. This genus is established for a single species, which is described below.
Hyperboreomyzon polaris Bolotov, Eliseeva, Klass & Kondakov gen. & sp. nov.Figures 4J, 5e, 6a-j, 7c, Supplementary Figs. S21, S8, S9, S10, S11, Supplementary Table S2.LSID: https://zoobank.org/urn:lsid:zoobank.org:act:503A9A26-CEDE-4747-952D-8416AE4EF4EB.Holotype. RMBH Hir_0486-H (sequenced: COI sequence acc. No. ON810753; 18S rRNA sequence acc. No. ON819030), RUSSIA: small alpine lake on Putorana Plateau, 68.9008°N, 94.1599°E, July, 2021, E. S. Chertoprud leg.Paratypes (N = 2). RUSSIA: 1 specimen RMBH Hir_0689 (dissected and placed on 60 permanent slides as a series of slices), small alpine lake on Putorana Plateau, 68.6659° N, 93.1365° E, August 11, 2021, E. S. Chertoprud leg.; 1 specimen RMBH Hir_0216 (sequenced and dissected; COI sequence acc. No. ON810677; 18S rRNA sequence acc. No. ON819005), water puddle on Kolguev Island, 68.9300° N, 49.0303° E, August 12, 2018, O. V. Travina & V. M. Spitsyn leg.Etymology. The name of the new species reflects its occurrences in polar (Arctic) areas of Eurasia.Differential diagnosis. As for the genus.Molecular diagnosis. None of congeneric species is known. Based on uncorrected pairwise COI p-distances between a haplotype of the new taxon and selected species-level haplotypes in each genus (Supplementary Table S1), Hyperboreomyzon seems to be more closely related to members of Hemiclepsis (mean distance ± s.e.m. = 11.62 ± 0.15%, range = 9.75–14.08%, N = 9) and Theromyzon (mean distance ± s.e.m. = 11.37 ± 0.07%, range = 10.47–12.64%, N = 9) without significant differences between distances from these two genera (Mann–Whitney test: p = 0.72). Other Glossiphoniidae genera are more distantly related, with a mean pairwise uncorrected COI p-distance of  > 13.0% (Mann–Whitney test: p  More

Study siteThe Cariaco Basin, located on the continental shelf off Venezuela, is a large (about 160 km long and about 65 km wide) depression, composed of two approximately 1,400-m-deep sub-basins. It is partially isolated from the Caribbean Sea by a series of sills with depths of less than 150 m (ref. 47). This limits renewal of deep water in the basin and, paired with the high oxygen demand resulting from intense surface primary productivity, leads to anoxic waters below a depth of about 275 m at present47,48.The marked seasonality in the Cariaco Basin, combined with anoxic bottom waters that effectively prevent bioturbation, results in the accumulation of annually laminated (varved) sediments. As sediments are varved for the last deglaciation and the Holocene, and because of the sensitivity of the area to climate change, they are considered to be one of the most valuable high-resolution marine climate archives and have been successfully used to study climate variability in the tropics3,11,16,17,18. Varve thickness is about 1 mm or more during the YD–Holocene transition18.Core and age modelCore MD03-2621 was retrieved during IMAGES cruise XI (PICASSO) aboard R/V Marion Dufresne in 2003 (Laj and Shipboard Party 2004). Cariaco cores have been collected under the regulations of the Ocean Drilling Program and the IMAGES coring programme. In this study, data from depths between 480 and 540 cm below the seafloor are presented, encompassing the YD–Holocene transition. A detailed age model for core MD03-2621 was established by Deplazes et al.11 and is based on the cross-correlation of total reflectance to dated colour records from the Cariaco Basin49,50. For the studied interval, the original age model is based on a floating varve chronology anchored to tree ring data by matching 14C data49. The age model for core MD03-2621 was further fine-tuned by correlation of reflectance data to the NGRIP ice core δ18O record on the GICC05 age scale11. The transition from the YD to the Holocene is characterized by a decrease in the sedimentation rate from 1.4 to 0.5 mm year−1.To account for possible depth offsets during storage and subsampling, we matched sediment colour data expressed as greyscale (GS) to the reflectance data from Deplazes et al.11 with the software QAnalySeries51. To enable comparison with our record, ages in Lea et al.3 were corrected for the age difference between the sediment-colour-based midpoint of the YD–Holocene transition in their record (11.56 kyr b2k) and in data from Deplazes et al.11 (11.673 kyr b2k). The start and end of the change in reflectance were determined by the RAMPFIT software52.Sample preparationSamples for MSI of molecular proxies were prepared as described in Alfken et al.53: the original core was subsampled by LL channels, from which X-ray pictures (Hewlett-Packard Faxitron 43855A X-ray cabinet) and high-resolution digital images (smart-CIS 1600 Line Scanner) were obtained. The LL channels were then cut into 5-cm pieces, which were subsequently freeze-dried, embedded in a gelatin:carboxymethyl cellulose (4%:1%) mixture and thin-sectioned on a Microm HM 505 E cryomicrotome. From each piece, one 60-µm-thick and one 100-µm-thick, longitudinal slice (spanning the whole 5 cm piece) were prepared and affixed to indium-tin-oxide-coated glass slides (Bruker Daltonik, Bremen, Germany) for MSI and elemental mapping, respectively. Slices for MSI were further amended with a fullerite matrix54.For all slices, a high-resolution picture was taken on a M4 Tornado micro-X-ray fluorescence spectroscopy system (Bruker Nano Analytics). This picture was used as a reference to set up elemental mapping and MSI analysis, and also for the 2D comparison of elemental and proxy data to sediment colour. Sediment colour is expressed as GS value. To account for differences between single slices, ΔGS was calculated as the difference between a value and the median GS of each individual slice. Very low GS values corresponding to areas devoid of sediment, identified by a black background, were excluded from analysis.Elemental mappingElemental mapping of 100-µm-thick slices was performed on a M4 Tornado micro-X-ray fluorescence spectroscopy system (Bruker Nano Analytics) equipped with a micro-focused Rh source (50 kV, 600 µA) with a polycapillary optic. Measurements were conducted under vacuum, with a resolution of 50 µm, two scans per spot and a scan time of 5 ms per scan. Data were initially processed and visualized with M4 Tornado Software version 1.3. XY matrices of relevant elements and sediment colour were imported into MATLAB (R2016b) for further processing. To assess the correspondence between sediment colour and elemental composition, for each 5-cm piece, signal intensities of Ca, Fe, Ti and Si in single spots were binned according to ΔGS and average intensities were calculated for each bin (Extended Data Fig. 5). The bin size was 5 GS units.Molecular proxy analysis by MSIMSI was carried out on a 7T solariX XR Fourier transform ion cyclotron resonance mass spectrometer coupled to a matrix-assisted laser desorption/ionization source equipped with a Smartbeam II laser (Bruker Daltonik, Bremen, Germany). Analyses were performed in positive ionization mode selecting for a continuous accumulation of selective ions window of m/z 554 ± 12. Spectra were acquired with 25% data reduction to limit data size. Spatial resolution was obtained by rastering the ionizing laser across the sample in a defined rectangular area at a 100-µm spot distance. Considering laminae thickness in the millimetre range18, such raster resolution is suited for seasonally resolved SST reconstruction. Settings for laser power, frequency and number of shots were adjusted for optimal signal intensities before each measurement; typical values were 250 shots with 200 Hz frequency and 60% laser power. External mass calibration was performed in electrospray ionization mode with sodium trifluoroacetate (Sigma-Aldrich). Each spectrum was also calibrated after data acquisition by an internal lock mass calibration using the Na+ adduct of pyropheophorbide a (m/z 557.2523), a chlorophyll a derivative generally present in relatively young marine sediments. Around 20,000 individual spots were thereby obtained for every 5-cm slice, each spot containing information on the abundance of diunsaturated and triunsaturated C37 alkenones needed to calculate the ({{rm{U}}}_{37}^{{{rm{K}}}^{{prime} }}) SST proxy.We provide a two-pronged approach to decode SST proxy information: (1) a downcore ({{rm{U}}}_{37}^{{{rm{K}}}^{{prime} }}) profile is obtained by pooling alkenone data from coeval horizons, and results in SST reconstructions with annual resolution, and (2) 2D images of alkenone distribution are examined in conjunction with maps of sediment colour and elemental distribution to filter single-spot alkenone data for season of deposition.SST reconstruction with yearly resolutionFor the downcore profile, MSI data were referenced to the X-ray image by the identification of three teaching points per 5-cm piece. Afterwards, the X-ray image was corrected for tilting of laminae in the LL channels. This was done by identification of single laminae in the X-ray image and selection of a minimum of four tie points per lamina. A detailed description can be found in Alfken et al.9. After applying the corresponding age model, downcore profiles were established with 1-year resolution: the intensity of the two alkenone species relevant to the ({{rm{U}}}_{37}^{{{rm{K}}}^{{prime} }}) proxy (C37:2 and C37:3) were recorded for each individual laser spot and filtered for a signal-to-noise threshold of 3. Only spots in which both compounds were detected were further considered. Intensity values were then summed over the depth corresponding to 1 year. By pooling proxy data into 1-year horizons, the effect of changing sedimentation rate and, thereby, changing downcore resolution is minimized. If at least ten spots presenting both compounds were available for a single horizon, data quality criteria were satisfied54 and a ({{rm{U}}}_{37}^{{{rm{K}}}^{{prime} }}) value was calculated as defined by Prahl and Wakeham22:$${{rm{U}}}_{37}^{{{rm{K}}}^{{prime} }}=frac{{text{C}}_{37:2}}{{text{C}}_{37:2}+{text{C}}_{37:3}}$$
(1)
To apply the gas chromatography (GC)-based calibrations for the ({{rm{U}}}_{37}^{{{rm{K}}}^{{prime} }}) proxy, MSI-based data were converted to GC equivalents. Therefore, after MSI, sediment slices were extracted for conventional proxy analysis. Sediment was scraped off the slide and extracted following a modified Bligh and Dyer procedure55,56. Extracts were evaporated under a stream of nitrogen, re-dissolved in n-hexane and analysed on a Thermo Finnigan Trace GC-FID equipped with a Restek Rxi-5ms capillary column (30 m × 0.25 mm ID). For each 5-cm piece, a ratio between the ({{rm{U}}}_{37}^{{{rm{K}}}^{{prime} }}) values obtained by GC flame ionization detector analysis and MSI of the whole piece was calculated. The average ratio of all pieces for which GC-based values could be obtained was 1.194, with a standard deviation of 0.021.$${{rm{U}}}_{37,{rm{G}}{rm{C}}-{rm{F}}{rm{I}}{rm{D}}}^{{{rm{K}}}^{{prime} }}=1.194times {{rm{U}}}_{37,{rm{M}}{rm{S}}{rm{I}}}^{{{rm{K}}}^{{prime} }}$$
(2)
This ratio was used to calculate GC-equivalent ({{rm{U}}}_{37}^{{{rm{K}}}^{{prime} }}) values, which were then translated into SST using the BAYSPLINE calibration57. The average standard error of the BAYSPLINE model is 0.049 ({{rm{U}}}_{37}^{{{rm{K}}}^{{prime} }}) units (corresponding to 1.4 °C) for samples with SST below 23.4 °C, but increases at higher values (to up to 4.4 °C)57. This is explained by the fact that sensitivity of the ({{rm{U}}}_{37}^{{{rm{K}}}^{{prime} }}) to SST (that is, the slope of the regression) declines at higher values. In the current dataset, the 95% confidence interval is, on average, ±3.6 °C. The analytical precision of MSI-based SST reconstructions for the ({{rm{U}}}_{37}^{{{rm{K}}}^{{prime} }}), using at least ten data points, according to Alfken et al.9, is about 0.3 °C. Sources of uncertainty are summarized in Extended Data Fig. 10a.For frequency analysis, a continuous, annually spaced record was constructed by linearly interpolating 49 missing values. The record was subsequently detrended. Spectral analysis was performed with the REDFIT module58 using a Hanning window (oversample 2, segments 2). Continuous wavelet transforms were applied to investigate changes in cyclicity over time, using the Morlet wavelet with code provided by Torrence and Compo59 for MATLAB. All steps, except for the wavelet analysis, were performed with the PAST software60.For the assessment of the interannual variability, the SST record was band-pass-filtered for periods between 2 and 8 years. The record is based on 1-year binned data; seasonality is thereby nullified and the highest frequency to be evaluated (Nyquist frequency) corresponds to a period of 2 years. Variability of this time series was quantified by calculating the standard deviation of the band-pass-filtered ({{rm{U}}}_{37}^{{{rm{K}}}^{{prime} }}) signal in 25-year intervals. To account for the potential impact of analytical precision on the observed signal (Methods, section titled ‘The effect of changing sedimentation rate on reconstructed interannual SST variability during the YD–Holocene transition’), the variability experiment from Alfken et al.9 was revisited. A sediment extract had been sprayed on an ITO slide and analysed by MSI. We then randomly selected n spots and obtained a ({{rm{U}}}_{37}^{{{rm{K}}}^{{prime} }}) value for the summed intensities of these spots. Precision was calculated as the standard deviation of five replicate experiments for n = 1, 5, 10, 20, 30, 40, 50 and 60. Decreasing analytical variability with increasing number of observations was fitted to a curve (R2 = 0.838) described by the equation$${rm{Analytical}},{rm{variability}}=0.0741times {rm{number}},{rm{of}},{{rm{spots}}}^{-0.558}$$
(3)
On the basis of this equation, analytical variability for each horizon could be calculated on the basis of the number of values included (Extended Data Fig. 10b). The mean variability for each 25-year window was then subtracted from the observed variability in the band-pass-filtered signal and the resulting proxy values were translated to SST following the equation by Müller et al.61. Statistical significance of the change in corrected SST variability after 11.66 kyr b2k was assessed with a t-test.Assessment of SST seasonalityFor the assessment of SST seasonality, alkenone intensities from individual spots were binned according to ΔGS, with a bin size of 1 unit. Spots were then separated into the categories upwelling season and non-upwelling season by identifying the threshold ΔGS value that maximized the difference between average SST in the bins above and below it. Furthermore, this value had to fulfill three conditions: (1) be higher (lighter) than the bins with the highest relative abundance of Ca, Ti and Fe, which is indicative of the dark sediments associated to non-upwelling season, (2) be lower (darker) than the bin with highest relative abundance for Si indicative of light sediment associated to the upwelling season and (3) the number of spots categorized as upwelling and non-upwelling had to account for at least 25% of total spots. If criteria 1 and 2 prevented criteria 3 from being fulfilled, a limit of 15% was set. After separating data into these two categories, data were processed separately as described above for the unfiltered dataset and a downcore temporal resolution of 5 years was applied. Seasonality was calculated as the difference between both records and thus represents the difference between 5-year average SST in the non-upwelling and upwelling seasons.Shift in seasonality was fitted to two different ramps with the RAMPFIT software52. An unconstrained approach and a constrained approach (in which the start and end points of the ramp were restricted to the intervals 11.725–11.8 kyr b2k and 11.6–11.675 kyr b2k) were applied. Negative values were excluded from this fitting. The resulting groups of data were compared by a Mann–Whitney rank test.SST seasonality in the modern Cariaco Basin was calculated for the years 1980 to 2020 based on the HadISST dataset62 by dividing monthly data from each year into two groups and searching for the largest difference between the average temperatures of both groups. Each group had to include at least three consecutive months. In 36 out of 41 years, the warm season was defined from May to November or from July to November.Decadal-scale to centennial-scale SST changes during the YD–Holocene transition and in the early HoloceneAnnually reconstructed SST (average SST = 24.3 °C) remains relatively stable during the YD–Holocene transition. At around 11.4 kyr b2k, a warming trend is observed. Averaging all data before 11.39 kyr and after 11.37 kyr results in a warming from 23.9 ± 1.6 °C to 25.5 ± 1.4 °C. Trends identified by MSI are consistent with conventional ({{rm{U}}}_{37}^{{{rm{K}}}^{{prime} }}) analyses performed in the present study and those previously reported by Herbert and Schuffert23 on Ocean Drilling Program core 165-1002C (Extended Data Fig. 1). These authors observed a slight warming several hundred years after the transition into the Holocene, between about 11.53 and 11.32 kyr b2k.Three prominent SST maxima are observed between about 11.50 and 11.45 kyr b2k. The average SST in these 50 years is 1.3 °C higher than in the 50 years before and after. These maxima are synchronous with the 11.4-ka cold event or PBO characterized by a negative excursion in δ18O and reduced snow accumulation rates in Greenland ice cores63 (Extended Data Fig. 2). The PBO coincides with the oldest of the Bond events, that is, pulses of ice rafting in the Northern Atlantic indicative of climatic deterioration64.A warm tropical response to the PBO would be supported by the lower-resolution foraminiferal SST record of Lea et al.3, which shows two data points of increased SST shortly after the end of the YD–Holocene transition. To enable direct comparison, ages in Lea et al.3 were corrected for the age difference between the sediment-colour-based YD termination midpoint in their record and in data from Deplazes et al.11. After this correction, these maxima correspond to 11.43 and 11.50 kyr b2k (Extended Data Fig. 2). Further, the SST maxima coincide with a short-lived change to lighter-coloured sediments. Hughen et al.19 described a correlation between brief North Atlantic cold events, such as the PBO, and changes in tropical primary productivity mediated by stronger upwelling that result in lighter sediments in the Cariaco Basin. Far-reaching effects of the PBO have previously been described in West Asia, with increased dust plumes being related to a southward shift of the westerlies65.The identification of the mechanisms behind a potential TNA response to the PBO is beyond the scope of this study. However, we wish to point out that high-resolution records are crucial to identify such events and to differentiate between underlying changes coinciding in time and, as in the present case, sharp signals that act on the same multidecadal timescales and can potentially be triggered by the same processes66.The effect of changing sedimentation rate on reconstructed interannual SST variability during the YD–Holocene transitionPooling proxy data into 1-year horizons establishes a constant sampling rate and thereby prevents potential effects of changing sedimentation rates. The onset of the Holocene in the Cariaco Basin sediments is characterized by a sharp decrease in sedimentation rates from 1.4 to 0.5 mm year−1 (refs. 11,19). Consequently, in the yearly pooled data, we observe a reduction in the number of values summed for each horizon (Extended Data Fig. 10b), as fewer laser spots fit into the thinner Holocene annual layers. At the same time, the mean intensity in each of these spots slightly increases, consistent with a relative increase of the contribution of haptophytes to primary production20.We have previously shown that the precision of MSI-based molecular proxy analysis is dependent on both the number of spots pooled per data point and the signal intensity in these spots54. All horizons used in the downcore record are above the established threshold of ten spots and proxy variability was shown to stabilize above this threshold9,54. However, as a decrease in the number of values per horizon might still result in lower analytical precision and contribute to higher signal variability, we corrected variability in the 2–8-year window with the estimated analytical variability (see equation (3)). With this correction, the magnitude of the described variability decreases across the record, but the trend towards higher interannual variability in the Holocene persists (Fig. 2c).Varve formation and alkenone deposition in the sediments of the Cariaco Basin during the YD–Holocene transitionComparison of elemental maps and sediment colour (Extended Data Fig. 5) shows a consistent pattern of lamination across the YD–Holocene transition that results from the seasonal interplay of precipitation, upwelling and dominant phytoplankton community composition. Darker laminae represent the rainy, non-upwelling (summer/fall) season and are enriched in Fe and Ti from terrigenous material and Ca sourced from biogenic CaCO3 produced by foraminifera or coccolithophores. Lighter laminae are characterized by high abundance of Si and correspond to the increased production of biogenic opal by diatoms during the upwelling (winter/spring) season67. This is in agreement with observations by Hughen et al.18, who described the laminae couplets in the Cariaco Basin as representing annual cycles, whereby light laminae are an indicator of high productivity associated with the winter/spring upwelling season and dark laminae are an indicator of summer/fall runoff and accumulation of terrigenous material. Deplazes et al.68 described a divergent origin of lamination for a deeper section of the YD, with light laminae being rich in calcareous and terrigenous elements characteristic for the summer season, whereas dark layers were enriched in Si and Br, indicative of diatoms and organic-walled primary producers characteristic for the more productive winter season. Such an alteration of the characteristic pattern of lamination is not observed in the late YD investigated here.This blueprint of seasonality was used to assess the seasonal behaviour of alkenones. Alkenones were deposited throughout the year, as evidenced by the fact that the number of spots containing detectable amounts of both alkenone species are not restricted to the upwelling or non-upwelling seasons but distributed across a relatively wide range of GS values to both sides of the median (Extended Data Fig. 6). Average alkenone signal intensity is higher in the non-upwelling season, pointing to a preference of alkenone producers for this season and/or to a stronger dilution of the signal in the upwelling season. In regards to the ({{rm{U}}}_{37}^{{{rm{K}}}^{{prime} }}) SST proxy distribution in light versus dark layers, our observations are in agreement with the ability to capture the seasonal SST cycle with alkenones in sinking particles in the modern Cariaco Basin69.Effect of changing seasonality on YD and early Holocene SST records from the western TNAChanging seasonality can contribute to explaining contrasting lower-resolution SST records in the western TNA during the YD and the early Holocene. The strong warming during the YD–Holocene transition recorded in the foraminiferal Mg/Ca record of the Cariaco Basin (Lea et al.3; Extended Data Fig. 1) might be reflecting the more robust thermohaline stratification and increasingly warmer non-upwelling seasons, given the preference of Globigerinoides ruber for this season.Globigerinoides ruber (white), as used by Lea et al.3, is considered to be a dominant species in the tropics, with a relatively uniform annual distribution. However, in the modern Cariaco Basin, upwelling leads to a distinct foraminiferal community composition and seasonal turnover70, consistent with the notion of warm-water foraminifera narrowing their occurrence to the warmest season71. The relative abundance of G. ruber increases in the non-upwelling (warm) season but rarely exceeds 15%, whereas the upwelling season is clearly dominated by Globigerina bulloides72,73. Globigerinoides ruber fluxes are consistently lowest when upwelling is most vigorous, as expressed in annual minima in SST (Extended Data Fig. 9b). As upwelling during the YD and early Holocene was more intense than in the present70, the preference of G. ruber for the summer (non-upwelling) season might have been even more pronounced.The development of a stronger seasonality in the early Holocene would thus have led to a narrower temporal occurrence of G. ruber in the non-upwelling season, during which it would also be exposed to higher SST. The average SST difference between seasons obtained in our analysis can be converted into annual SST amplitude by assuming a sinusoidal curve. By doing so, we observe an increase in the seasonal amplitude of 1.5 to 1.9 °C (depending on the ramp fitted), which is similar to the warming described by Lea et al.3.This interpretation is in agreement with Bova et al.46, who observed that most Holocene climate reconstructions are biased towards the boreal summer/fall and reflect the evolution of seasonal rather than annual temperatures. As discussed above, this is probably not true for the ({{rm{U}}}_{37}^{{{rm{K}}}^{{prime} }}) index in the Cariaco Basin, as alkenones are deposited throughout the year. The suggested weakening of summer stratification during the YD (as compared with the Holocene) might, however, explain why the lower-resolution ({{rm{U}}}_{37}^{{{rm{K}}}^{{prime} }}) records from the semi-enclosed Cariaco Basin show no or weaker warming23 than other, open-ocean, tropical YD records4, where the interplay of upwelling, freshwater input and stratification are less relevant to the SST signal. More

Hawley, W. A. The biology of Aedes albopictus. J. Am. Mosq. Control Assoc. 1, 1–39 (1988).CAS 

Google Scholar 
Benedict, M. Q., Levine, R. S., Hawley, W. A. & Lounibos, L. P. Spread of the tiger: global risk of invasion by the mosquito Aedes albopictus. Vector-Borne Zoonotic Dis. 7, 76–85 (2007).Article 
PubMed 

Google Scholar 
Paupy, C., Delatte, H., Bagny, L., Corbel, V. & Fontenille, D. Aedes albopictus, an arbovirus vector: From the darkness to the light. Microb. Infect. 11, 1177–1185 (2009).Article 
CAS 

Google Scholar 
Delatte, H. et al. Blood-feeding behavior of Aedes albopictus, a vector of Chikungunya on La Réunion. Vector-Borne Zoonotic Dis. 10, 249–258 (2010).Article 
PubMed 

Google Scholar 
Pereira-dos-Santos, T., Roiz, D., Lourenço-de-Oliveira, R. & Paupy, C. A systematic review: Is Aedes albopictus an efficient bridge vector for zoonotic arboviruses? Pathogens 9, 266 (2020).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Gratz, N. Critical review of the vector status of Aedes albopictus. Med. Vet. Entomol. 18, 215–227 (2004).Article 
CAS 
PubMed 

Google Scholar 
Grard, G. et al. Zika virus in Gabon (Central Africa)—2007: A new threat from Aedes albopictus? PLoS Negl. Trop. Dis. 8, e2681 (2014).Article 
PubMed 
PubMed Central 

Google Scholar 
Lambrechts, L., Scott, T. W. & Gubler, D. J. Consequences of the expanding global distribution of Aedes albopictus for dengue virus transmission. PLoS Negl. Trop. Dis. 4, e646 (2010).Article 
PubMed 
PubMed Central 

Google Scholar 
Lounibos, L. P. & Kramer, L. D. Invasiveness of Aedes aegypti and Aedes albopictus and vectorial capacity for chikungunya virus. J. Infect. Dis. 214, S453–S458 (2016).Article 
PubMed 
PubMed Central 

Google Scholar 
European Centre for Disease Prevention and Control (ECDC). Vector Control with a Focus on Aedes aegypti and Aedes albopictus Mosquitoes: Literature Review and Analysis of Information (ECDC, Stockholm, Sweden, 2017).Tatem, A. J., Hay, S. I. & Rogers, D. J. Global traffic and disease vector dispersal. PNAS 103, 6242–6247 (2006).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Lowe, S., Browne, M., Boudjelas, S. & De Poorter, M. 100 of the World’s Worst Invasive Alien Species: A Selection From the Global Invasive Species Database, Vol. 12 (Invasive Species Specialist Group, 2000).Diagne, C. et al. High and rising economic costs of biological invasions worldwide. Nature 592, 571–576 (2021).Article 
CAS 
PubMed 

Google Scholar 
Hulme, P. E. Trade, transport and trouble: Managing invasive species pathways in an era of globalization. J. Appl. Ecol. 46, 10–18 (2009).Article 

Google Scholar 
Marini, F., Caputo, B., Pombi, M., Tarsitani, G. & Della-Torre, A. Study of Aedes albopictus dispersal in Rome, Italy, using sticky traps in mark–release–recapture experiments. Med. Vet. Entomol. 24, 361–368 (2010).Article 
CAS 
PubMed 

Google Scholar 
Bonizzoni, M., Gasperi, G., Chen, X. & James, A. A. The invasive mosquito species Aedes albopictus: current knowledge and future perspectives. Trends Parasitol. 29, 460–468 (2013).Article 
PubMed 
PubMed Central 

Google Scholar 
Collantes, F. et al. Review of ten-years presence of Aedes albopictus in Spain 2004–2014: Known distribution and public health concerns. Parasit Vectors 8, 1–11 (2015).Article 

Google Scholar 
Aranda, C., Eritja, R. & Roiz, D. First record and establishment of the mosquito Aedes albopictus in Spain. Med. Vet. Entomol. 20, 150–152 (2006).Article 
CAS 
PubMed 

Google Scholar 
Giménez, N. et al. Introduction of Aedes albopictus in Spain: A new challenge for public health. Gac. Sanit. 21, 25–28 (2007).Article 
PubMed 

Google Scholar 
European Centre for Disease Prevention and Control and European Food Safety Authority. Mosquito maps [internet]. Stockholm: ECDC. https://ecdc.europa.eu/en/disease-vectors/surveillance-and-disease-data/mosquito-maps (2022).Shigesada, N. & Kawasaki, K. Biological Invasions: Theory and Practice (Oxford University Press, 1997).
Google Scholar 
Puth, L. M. & Post, D. M. Studying invasion: Have we missed the boat? Ecol. Lett. 8, 715–721 (2005).Article 

Google Scholar 
Leung, B. et al. An ounce of prevention or a pound of cure: Bioeconomic risk analysis of invasive species. Proc. R Soc. Lond. Ser. B Biol. Sci. 269, 2407–2413 (2002).Article 

Google Scholar 
Lounibos, L. P. Invasions by insect vectors of human disease. Annu. Rev. Entomol. 47, 233–266 (2002).Article 
CAS 
PubMed 

Google Scholar 
Manni, M. et al. Genetic evidence for a worldwide chaotic dispersion pattern of the arbovirus vector, Aedes albopictus. PLoS Negl. Trop. Dis. 11, e0005332 (2017).Article 
PubMed 
PubMed Central 

Google Scholar 
Roiz, D. et al. Integrated Aedes management for the control of Aedes-borne diseases. PLoS Negl. Trop. Dis. 12, e0006845 (2018).Article 
PubMed 
PubMed Central 

Google Scholar 
Lühken, R. et al. Microsatellite typing of Aedes albopictus (Diptera: Culicidae) populations from Germany suggests regular introductions. Infect. Genet. Evol. 81, 104237 (2020).Article 
PubMed 

Google Scholar 
Battaglia, V. et al. The worldwide spread of the tiger mosquito as revealed by mitogenome haplogroup diversity. Front. Genet. 7, 208 (2016).Article 
PubMed 
PubMed Central 

Google Scholar 
Medley, K. A., Jenkins, D. G. & Hoffman, E. A. Human-aided and natural dispersal drive gene flow across the range of an invasive mosquito. Mol. Ecol. 24, 284–295 (2015).Article 
PubMed 

Google Scholar 
Eritja, R., Palmer, J. R., Roiz, D., Sanpera-Calbet, I. & Bartumeus, F. Direct evidence of adult Aedes albopictus dispersal by car. Sci. Rep. 7, 1–15 (2017).Article 
CAS 

Google Scholar 
Sherpa, S. et al. Unravelling the invasion history of the Asian tiger mosquito in Europe. Mol. Ecol. 28, 2360–2377 (2019).Article 
PubMed 

Google Scholar 
Swan, T. et al. A literature review of dispersal pathways of Aedes albopictus across different spatial scales: Implications for vector surveillance. Parasit Vectors 15, 1–13 (2022).Article 

Google Scholar 
Ballard, J. W. O. & Whitlock, M. C. The incomplete natural history of mitochondria. Mol. Ecol. 13, 729–744. https://doi.org/10.1046/j.1365-294X.2003.02063.x (2004).Article 
PubMed 

Google Scholar 
Toews, D. P. L. & Brelsford, A. The biogeography of mitochondrial and nuclear discordance in animals. Mol. Ecol. 21, 3907–3930. https://doi.org/10.1111/j.1365-294X.2012.05664.x (2012).Article 
CAS 
PubMed 

Google Scholar 
Hurst, G. D. & Jiggins, F. M. Problems with mitochondrial DNA as a marker in population, phylogeographic and phylogenetic studies: The effects of inherited symbionts. Proc. R. Soc. B: Biol. Sci. 272, 1525–1534 (2005).Article 
CAS 

Google Scholar 
Cariou, M., Duret, L. & Charlat, S. The global impact of Wolbachia on mitochondrial diversity and evolution. J. Evol. Biol. 30, 2204–2210 (2017).Article 
CAS 
PubMed 

Google Scholar 
Zug, R. & Hammerstein, P. Still a host of hosts for Wolbachia: analysis of recent data suggests that 40% of terrestrial arthropod species are infected. PLoS ONE 7, e38544 (2012).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Weinert, L. A., Araujo-Jnr, E. V., Ahmed, M. Z. & Welch, J. J. The incidence of bacterial endosymbionts in terrestrial arthropods. Proc. R. Soc. B: Biol. Sci. 282, 20150249 (2015).Article 

Google Scholar 
Goubert, C., Minard, G., Vieira, C. & Boulesteix, M. Population genetics of the Asian tiger mosquito Aedes albopictus, an invasive vector of human diseases. Heredity 117, 125–134 (2016).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Western, D. Human-modified ecosystems and future evolution. PNAS 98, 5458–5465 (2001).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Pech-May, A. et al. Population genetics and ecological niche of invasive Aedes albopictus in Mexico. Acta Trop. 157, 30–41 (2016).Article 
PubMed 

Google Scholar 
Vargo, E. L. et al. Hierarchical genetic analysis of German cockroach (Blattella germanica) populations from within buildings to across continents. PLoS ONE 9, e102321 (2014).Article 
PubMed 
PubMed Central 

Google Scholar 
von Beeren, C., Stoeckle, M. Y., Xia, J., Burke, G. & Kronauer, D. J. Interbreeding among deeply divergent mitochondrial lineages in the American cockroach (Periplaneta americana). Sci. Rep. 5, 1–7 (2015).
Google Scholar 
Tseng, S.-P. et al. Genetic diversity and Wolbachia infection patterns in a globally distributed invasive ant. Front. Genet. 10, 838 (2019).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Wesson, D. M., Porter, C. H. & Collins, F. H. Sequence and secondary structure comparisons of ITS rDNA in mosquitoes (Diptera: Culicidae). Mol. Phylogen. Evol. 1, 253–269 (1992).Article 
CAS 

Google Scholar 
Mishra, S., Sharma, G., Das, M. K., Pande, V. & Singh, O. P. Intragenomic sequence variations in the second internal transcribed spacer (ITS2) ribosomal DNA of the malaria vector Anopheles stephensi. PLoS ONE 16, e0253173 (2021).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Artigas, P. et al. Aedes albopictus diversity and relationships in south-western Europe and Brazil by rDNA/mtDNA and phenotypic analyses: ITS-2, a useful marker for spread studies. Parasit Vectors 14, 1–23 (2021).Article 

Google Scholar 
Armbruster, P. et al. Infection of New-and Old-World Aedes albopictus (Diptera: Culicidae) by the intracellular parasite Wolbachia: implications for host mitochondrial DNA evolution. J. Med. Entomol. 40, 356–360 (2003).Article 
PubMed 

Google Scholar 
Maia, R., Scarpassa, V. M., Maciel-Litaiff, L. & Tadei, W. P. Reduced levels of genetic variation in Aedes albopictus (Diptera: Culicidae) from Manaus, Amazonas State, Brazil, based on analysis of the mitochondrial DNA ND5 gene. Gen. Mol. Res. 2000, 998–1007 (2009).Article 

Google Scholar 
Birungi, J. & Munstermann, L. E. Genetic structure of Aedes albopictus (Diptera: Culicidae) populations based on mitochondrial ND5 sequences: Evidence for an independent invasion into Brazil and United States. Ann. Entomol. Soc. Am. 95, 125–132 (2002).Article 
CAS 

Google Scholar 
Kambhampati, S. & Rai, K. S. Mitochondrial DNA variation within and among populations of the mosquito Aedes albopictus. Genome 34, 288–292 (1991).Article 
CAS 
PubMed 

Google Scholar 
Werren, J. H., Baldo, L. & Clark, M. E. Wolbachia: master manipulators of invertebrate biology. Nat. Rev. Microbiol. 6, 741–751 (2008).Article 
CAS 
PubMed 

Google Scholar 
Wiwatanaratanabutr, I. Geographic distribution of wolbachial infections in mosquitoes from Thailand. J. Invertebr. Pathol. 114, 337–340 (2013).Article 
PubMed 

Google Scholar 
Carvajal, T. M., Hashimoto, K., Harnandika, R. K., Amalin, D. M. & Watanabe, K. Detection of Wolbachia in field-collected Aedes aegypti mosquitoes in metropolitan Manila, Philippines. Parasit. Vectors 12, 1–9 (2019).Article 

Google Scholar 
Atyame, C. M., Delsuc, F., Pasteur, N., Weill, M. & Duron, O. Diversification of Wolbachia endosymbiont in the Culex pipiens mosquito. Mol. Biol. Evol. 28, 2761–2772 (2011).Article 
CAS 
PubMed 

Google Scholar 
Damiani, C. et al. Wolbachia in Aedes koreicus: Rare detections and possible implications. Insects 13, 216 (2022).Article 
PubMed 
PubMed Central 

Google Scholar 
Jiggins, F. M. Male-killing Wolbachia and mitochondrial DNA: Selective sweeps, hybrid introgression and parasite population dynamics. Genetics 164, 5–12 (2003).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Schuler, H. et al. The hitchhiker’s guide to Europe: The infection dynamics of an ongoing Wolbachia invasion and mitochondrial selective sweep in Rhagoletis cerasi. Mol. Ecol. 25, 1595–1609 (2016).Article 
PubMed 
PubMed Central 

Google Scholar 
Ross, P. A., Ritchie, S. A., Axford, J. K. & Hoffmann, A. A. Loss of cytoplasmic incompatibility in Wolbachia-infected Aedes aegypti under field conditions. PLoS Negl. Trop. Dis. 13, e0007357 (2019).Article 
PubMed 
PubMed Central 

Google Scholar 
Avise, J. C. Phylogeography: The history and formation of species (Harvard University Press, 2000).Book 

Google Scholar 
Rokas, A., Atkinson, R. J., Brown, G. S., West, S. A. & Stone, G. N. Understanding patterns of genetic diversity in the oak gallwasp Biorhiza pallida: Demographic history or a Wolbachia selective sweep? Heredity 87, 294–304 (2001).Article 
CAS 
PubMed 

Google Scholar 
Porretta, D., Mastrantonio, V., Bellini, R., Somboon, P. & Urbanelli, S. Glacial history of a modern invader: Phylogeography and species distribution modelling of the Asian tiger mosquito Aedes albopictus. PLoS ONE 7, e44515. https://doi.org/10.1371/journal.pone.0044515 (2012).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Motoki, M. T. et al. Population genetics of Aedes albopictus (Diptera: Culicidae) in its native range in Lao People’s Democratic Republic. Parasit. Vectors 12, 1–12 (2019).Article 
CAS 

Google Scholar 
Zhong, D. et al. Genetic analysis of invasive Aedes albopictus populations in Los Angeles County, California and its potential public health impact. PLoS ONE 8, e68586 (2013).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Usmani-Brown, S., Cohnstaedt, L. & Munstermann, L. E. Population genetics of Aedes albopictus (Diptera: Culicidae) invading populations, using mitochondrial nicotinamide adenine dinucleotide dehydrogenase subunit 5 sequences. Ann. Entomol. Soc. Am. 102, 144–150 (2009).Article 
CAS 
PubMed 

Google Scholar 
Mousson, L. et al. Phylogeography of Aedes (Stegomyia) aegypti (L.) and Aedes (Stegomyia) albopictus (Skuse) (Diptera: Culicidae) based on mitochondrial DNA variations. Genet. Res. 86, 1–11 (2005).Article 
CAS 
PubMed 

Google Scholar 
Bazin, E., Glémin, S. & Galtier, N. Population size does not influence mitochondrial genetic diversity in animals. Science 312, 570–572. https://doi.org/10.1126/science.1122033 (2006).Article 
CAS 
PubMed 

Google Scholar 
Dowling, D. K., Friberg, U. & Lindell, J. Evolutionary implications of non-neutral mitochondrial genetic variation. Ecol. Evol. 23, 546–554 (2008).Article 

Google Scholar 
Montero-Pau, J., Gómez, A. & Muñoz, J. Application of an inexpensive and high-throughput genomic DNA extraction method for the molecular ecology of zooplanktonic diapausing eggs. Limnol. Oceanogr. Methods 6, 218–222 (2008).Article 
CAS 

Google Scholar 
Porter, C. H. & Collins, F. H. Species-diagnostic differences in a ribosomal DNA internal transcribed spacer from the sibling species Anopheles freeborni and Anopheles hermsi (Diptera: Culicidae). Am. J. Trop. Med. 45, 271–279 (1991).Article 
CAS 

Google Scholar 
Folmer, O., Black, M., Hoeh, W., Lutz, R. & Vrijenhoek, R. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol. Mar. Biol. Biotechnol. 3, 294–299 (1994).CAS 
PubMed 

Google Scholar 
Prosser, S., Martínez-Arce, A. & Elías-Gutiérrez, M. A new set of primers for COI amplification from freshwater microcrustaceans. Mol. Ecol. Resour. 13, 1151–1155 (2013).CAS 
PubMed 

Google Scholar 
Ivanova, N. V., Zemlak, T. S., Hanner, R. H. & Hebert, P. D. Universal primer cocktails for fish DNA barcoding. Mol. Ecol. Notes 7, 544–548 (2007).Article 
CAS 

Google Scholar 
Kumar, S., Stecher, G. & Tamura, K. MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 33, 1870–1874 (2016).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Zhou, W., Rousset, F. & O’Neill, S. Phylogeny and PCR–based classification of Wolbachia strains using wsp gene sequences. Proc. R Soc. Lond. Ser. B Biol. Sci. 265, 509–515 (1998).Article 
CAS 

Google Scholar 
Braig, H. R., Zhou, W., Dobson, S. L. & O’Neill, S. L. Cloning and characterization of a gene encoding the major surface protein of the bacterial endosymbiont Wolbachia pipientis. J. Bacteriol. 180, 2373–2378 (1998).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Hu, Y. et al. Identification and molecular characterization of Wolbachia strains in natural populations of Aedes albopictus in China. Parasit. Vectors 13, 1–14 (2020).
Google Scholar 
Heddi, A., Grenier, A.-M., Khatchadourian, C., Charles, H. & Nardon, P. Four intracellular genomes direct weevil biology: Nuclear, mitochondrial, principal endosymbiont, and Wolbachia. PNAS 96, 6814–6819 (1999).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Rozas, J. et al. DnaSP 6: DNA sequence polymorphism analysis of large data sets. Mol. Biol. Evol. 34, 3299–3302 (2017).Article 
CAS 
PubMed 

Google Scholar 
Salzburger, W., Ewing, G. B. & Von Haeseler, A. The performance of phylogenetic algorithms in estimating haplotype genealogies with migration. Mol. Ecol. 20, 1952–1963 (2011).Article 
PubMed 

Google Scholar 
Stamatakis, A. RAxML-VI-HPC: Maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics 22, 2688–2690 (2006).Article 
CAS 
PubMed 

Google Scholar 
Darriba, D., Taboada, G. L., Doallo, R. & Posada, D. jModelTest 2: More models, new heuristics and parallel computing. Nat. Methods 9, 772. https://doi.org/10.1038/nmeth.2109 (2012).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Gower, J. C. Some distance properties of latent root and vector methods used in multivariate analysis. Biometrika 53, 325–338 (1966).Article 
MathSciNet 
MATH 

Google Scholar 
R Core Team. R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing. http://www.R-project.org/. (2021).Hijmans, R. J., Williams, E., Vennes, C. & Hijmans, M. R. J. Package ‘geosphere’. Spher. Trigon. 1, 5 (2017).
Google Scholar 
Palmer, J. R. et al. Citizen science provides a reliable and scalable tool to track disease-carrying mosquitoes. Nat. Commun. 8, 1–13 (2017).Article 

Google Scholar 
Mantel, N. & Valand, R. S. A technique of nonparametric multivariate analysis. Biometrics 1970, 547–558 (1970).Article 

Google Scholar 
Goslee, S. C. & Urban, D. L. The ecodist package for dissimilarity-based analysis of ecological data. J. Stat. Softw. 22, 1–19 (2007).Article 

Google Scholar 
Stewart, C. Zero-inflated beta distribution for modeling the proportions in quantitative fatty acid signature analysis. J. Appl. Stat. 40, 985–992 (2013).Article 
MathSciNet 
MATH 

Google Scholar 
Figueroa-Zúñiga, J. I., Arellano-Valle, R. B. & Ferrari, S. L. Mixed beta regression: A Bayesian perspective. Comput. Stat. Data Anal. 61, 137–147 (2013).Article 
MathSciNet 
MATH 

Google Scholar 
Branscum, A. J., Johnson, W. O. & Thurmond, M. C. Bayesian beta regression: Applications to household expenditure data and genetic distance between foot-and-mouth disease viruses. Aust. N. Z. J. Stat. 49, 287–301 (2007).Article 
MathSciNet 
MATH 

Google Scholar 
Ospina, R. & Ferrari, S. L. Inflated beta distributions. Stat. Pap. 51, 111–126 (2010).Article 
MathSciNet 
MATH 

Google Scholar 
Chung, H. & Beretvas, S. N. The impact of ignoring multiple membership data structures in multilevel models. Br. J. Math. Stat. Psychol. 65, 185–200 (2012).Article 
MathSciNet 
PubMed 
MATH 

Google Scholar  More

Author notesMikayla M. ClarkPresent address: University of Tennessee, Knoxville, TN, USAMichael W. SneddonPresent address: Predicine, Inc., Hayward, CA, USARoman SutorminPresent address: Google, Inc., San Francisco, CA, USAAuthors and AffiliationsLawrence Berkeley National Laboratory, Berkeley, CA, USADylan Chivian, Sean P. Jungbluth, Paramvir S. Dehal, Elisha M. Wood-Charlson, Richard S. Canon, Gavin A. Price, William J. Riehl, Michael W. Sneddon, Roman Sutormin & Adam P. ArkinOak Ridge National Laboratory, Oak Ridge, TN, USABenjamin H. Allen, Mikayla M. Clark, Miriam L. Land & Robert W. CottinghamArgonne National Laboratory, Lemont, IL, USATianhao Gu, Qizhi Zhang & Chris S. HenryAuthorsDylan ChivianSean P. JungbluthParamvir S. DehalElisha M. Wood-CharlsonRichard S. CanonBenjamin H. AllenMikayla M. ClarkTianhao GuMiriam L. LandGavin A. PriceWilliam J. RiehlMichael W. SneddonRoman SutorminQizhi ZhangRobert W. CottinghamChris S. HenryAdam P. ArkinCorresponding authorsCorrespondence to
Dylan Chivian or Adam P. Arkin. More

McFall-Ngai, M. et al. Animals in a bacterial world, a new imperative for the life sciences. Proc. Natl. Acad. Sci. 110, 3229–3236 (2013).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Archibald, J. M. Endosymbiosis and eukaryotic cell evolution. Curr. Biol. 25, R911–R921 (2015).Article 
CAS 
PubMed 

Google Scholar 
Moran, N. A. Symbiosis as an adaptive process and source of phenotypic complexity. Proc. Natl. Acad. Sci. U.S.A. 104(Suppl 1), 8627–8633 (2007).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Hurst, G. D. D. Extended genomes: Symbiosis and evolution. Interface Focus. https://doi.org/10.1098/rsfs.2017.0001 (2017).Article 
PubMed 
PubMed Central 

Google Scholar 
Moran, N. A., McCutcheon, J. P. & Nakabachi, A. Genomics and evolution of heritable bacterial symbionts. Annu. Rev. Genet. 42, 165–190 (2008).Article 
CAS 
PubMed 

Google Scholar 
Kikuchi, Y., Hosokawa, T. & Fukatsu, T. Insect-microbe mutualism without vertical transmission: A stinkbug acquires a beneficial gut symbiont from the environment every generation. Appl. Environ. Microbiol. 73, 4308–4316 (2007).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Kikuchi, Y., Hosokawa, T. & Fukatsu, T. An ancient but promiscuous host-symbiont association between Burkholderia gut symbionts and their heteropteran hosts. ISME J. 5, 446–460 (2011).Article 
PubMed 

Google Scholar 
Hu, Y. et al. Herbivorous turtle ants obtain essential nutrients from a conserved nitrogen-recycling gut microbiome. Nat. Commun. 9, 2440. https://doi.org/10.1038/s41467-018-03357-y (2018).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Salem, H. et al. Drastic genome reduction in an herbivore’s pectinolytic symbiont. Cell 171, 1520–1531 (2017).Article 
CAS 
PubMed 

Google Scholar 
Bennett, G. M. & Moran, N. A. Heritable symbiosis: The advantages and perils of an evolutionary rabbit hole. Proc. Natl. Acad. Sci. 112, 10169–10176 (2015).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Campbell, M. A. et al. Changes in endosymbiont complexity drive host-level compensatory adaptations in cicadas. MBio 9, e02104-18 (2018).Article 
PubMed 
PubMed Central 

Google Scholar 
Buchner, P. Symbiosis in animals which suck plant juices. In Endosymbiosis of Animals with Plant Microorganisms 210–432 (Interscience, 1965).
Google Scholar 
McCutcheon, J. P., McDonald, B. R. & Moran, N. A. Convergent evolution of metabolic roles in bacterial co-symbionts of insects. Proc. Natl. Acad. Sci. 106, 15394–15399 (2009).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Christensen, H. & Fogel, M. L. Feeding ecology and evidence for amino acid synthesis in the periodical cicada (Magicicada). J. Insect Physiol. 57, 211–219 (2011).Article 
CAS 
PubMed 

Google Scholar 
McCutcheon, J. P., McDonald, B. R. & Moran, N. A. Origin of an alternative genetic code in the extremely small and GC-rich genome of a bacterial symbiont. PLoS Genet. 5, e1000565 (2009).Article 
PubMed 
PubMed Central 

Google Scholar 
Campbell, M. A. et al. Genome expansion via lineage splitting and genome reduction in the cicada endosymbiont Hodgkinia. Proc. Natl. Acad. Sci. 112, 10192–10199 (2015).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Müller, H. J. Neuere vorstellungen über verbreitung und phylogenie der endosymbiosen der zikaden. Z. Morphol. Oekol. Tiere 61, 190–210 (1962).Article 

Google Scholar 
Müller, H. J. Zur systematik und phylogenie der zikaden-endosymbiosen. Biol. Zent. 68, 343–368 (1949).
Google Scholar 
Matsuura, Y. et al. Recurrent symbiont recruitment from fungal parasites in cicadas. Proc. Natl. Acad. Sci. 115, E5970–E5979 (2018).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Zhou, W. et al. Analysis of inter-individual bacterial variation in gut of cicada Meimuna mongolica (Hemiptera: Cicadidae). J. Insect Sci. 15, 1–6 (2015).Article 

Google Scholar 
Zheng, Z., Wang, D., He, H. & Wei, C. Bacterial diversity of bacteriomes and organs of reproductive, digestive and excretory systems in two cicada species (Hemiptera: Cicadidae). PLoS One 12, 1–21 (2017).
Google Scholar 
Wang, D., Huang, Z., He, H. & Wei, C. Comparative analysis of microbial communities associated with bacteriomes, reproductive organs and eggs of the cicada Subpsaltria yangi. Arch. Microbiol. 200, 227–235 (2018).Article 
CAS 
PubMed 

Google Scholar 
Dillon, R. J. & Dillon, V. M. The gut bacteria of insects: Nonpathogenic interactions. Annu. Rev. Entomol. 49, 71–92 (2004).Article 
CAS 
PubMed 

Google Scholar 
Ng, S. H., Stat, M., Bunce, M. & Simmons, L. W. The influence of diet and environment on the gut microbial community of field crickets. Ecol. Evol. 8, 4704–4720 (2018).Article 
PubMed 
PubMed Central 

Google Scholar 
Nishida, A. H. & Ochman, H. Rates of gut microbiome divergence in mammals. Mol. Ecol. 27, 1884–1897 (2018).Article 
PubMed 
PubMed Central 

Google Scholar 
Douglas, A. E. & Werren, J. H. Holes in the hologenome: Why host–microbe symbioses are not holobionts. MBio 7, e02099 (2016).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Grueneberg, J., Engelen, A. H., Costa, R. & Wichard, T. Macroalgal morphogenesis induced by waterborne compounds and bacteria in coastal seawater. PLoS One 11, e0146307 (2016).Article 
PubMed 
PubMed Central 

Google Scholar 
Lin, J. D., Lemay, M. A. & Parfrey, L. W. Diverse bacteria utilize alginate within the microbiome of the giant kelp Macrocystis pyrifera. Front. Microbiol. 9, 1914 (2018).Article 
PubMed 
PubMed Central 

Google Scholar 
Coon, K. L., Vogel, K. J., Brown, M. R. & Strand, M. R. Mosquitoes rely on their gut microbiota for development. Mol. Ecol. 23, 2727–2739 (2014).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Coon, K. L., Brown, M. R. & Strand, M. R. Mosquitoes host communities of bacteria that are essential for development but vary greatly between local habitats. Mol. Ecol. 25, 5806–5826 (2016).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Kwong, W. K. et al. Dynamic microbiome evolution in social bees. Sci. Adv. 3, 1–17 (2017).Article 

Google Scholar 
Brooks, A. W., Kohl, K. D., Brucker, R. M., van Opstal, E. J. & Bordenstein, S. R. Phylosymbiosis: Relationships and functional effects of microbial communities across host evolutionary history. PLoS Biol. 14, 1–29 (2016).Article 

Google Scholar 
Kropáčková, L. et al. Codiversification of gastrointestinal microbiota and phylogeny in passerines is not explained by ecological divergence. Mol. Ecol. 26, 5292–5304 (2017).Article 
PubMed 

Google Scholar 
Hird, S. M., Sánchez, C., Carstens, B. C. & Brumfield, R. Comparative gut microbiota of 59 neotropical bird species. Front. Microbiol. 6, 1403. https://doi.org/10.3389/fmicb.2015.01403 (2015).Article 
PubMed 
PubMed Central 

Google Scholar 
Hu, Y., Lukasik, P., Moreau, C. S. & Russell, J. A. Correlates of gut community composition across an ant species (Cephalotes varians) elucidate causes and consequences of symbiotic variability. Mol. Ecol. 23, 1284–1300 (2014).Article 
PubMed 

Google Scholar 
Hammer, T. J., Sanders, J. G. & Fierer, N. Not all animals need a microbiome. FEMS Microbiol. Lett. https://doi.org/10.1093/femsle/fnz117 (2019).Article 
PubMed 

Google Scholar 
Hammer, T. J., Janzen, D. H., Hallwachs, W., Jaffe, S. P. & Fierer, N. Caterpillars lack a resident gut microbiome. Proc. Natl. Acad. Sci. 114, 9641–9646 (2017).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Shapira, M. Gut microbiotas and host evolution: Scaling up symbiosis. Trends Ecol. Evol. 31, 539–549 (2016).Article 
PubMed 

Google Scholar 
Marshall, D. C. et al. Inflation of molecular clock rates and dates: Molecular phylogenetics, biogeography, and diversification of a global cicada radiation from Australasia (Hemiptera: Cicadidae: Cicadettini). Syst. Biol. 65, 16–34 (2016).Article 
PubMed 

Google Scholar 
Lane, D. H. The recognition concept of speciation applied in an analysis of putative hybridization in New Zealand cicadas of the genus Kikihia (Insects: Hemiptera: Tibicinidae). Speciation and the Recognition Concept: Theory and Application (The Johns Hopkins Univ Press, 1995).
Google Scholar 
Cooley, J. R. & Marshall, D. C. Sexual signaling in periodical cicadas, Magicicada spp. (Hemiptera: Cicadidae). Behaviour 138, 827–855 (2001).Article 

Google Scholar 
Fleming, C. A. Adaptive Radiation in New Zealand Cicadas (American Philosophical Society, 1975).
Google Scholar 
Dugdale, J. S. & Fleming, C. A. New Zealand cicadas of the genus Maoricicada (Homoptera: Tibicinidae). N. Z. J. Zool. 5, 295–340 (1978).Article 

Google Scholar 
Marshall, D. C., Hill, K. B. R., Cooley, J. R. & Simon, C. Hybridization, mitochondrial DNA phylogeography, and prediction of the early stages of reproductive isolation: Lessons from New Zealand cicadas (genus Kikihia). Syst. Biol. 60, 482–502 (2011).Article 
PubMed 

Google Scholar 
Bolyen, E. et al. QIIME 2: Reproducible, interactive, scalable, and extensible microbiome data science. Report No.: e27295v2. PeerJ https://doi.org/10.7287/peerj.preprints.27295v2 (2018).Article 

Google Scholar 
Callahan, B. J. et al. DADA2: High-resolution sample inference from Illumina amplicon data. Nat. Methods 13, 581–583 (2016).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Katoh, K. & Standley, D. M. MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Mol. Biol. Evol. 30, 772–780 (2013).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
McMurdie, P. J. & Holmes, S. phyloseq: An R package for reproducible interactive analysis and graphics of microbiome census data. PLoS One 8, e61217 (2013).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Davis, N. M., Proctor, D., Holmes, S. P., Relman, D. A. & Callahan, B. J. Simple statistical identification and removal of contaminant sequences in marker-gene and metagenomics data. Microbiome 6, 226. https://doi.org/10.1101/221499 (2018).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Lemmon, A. R., Emme, S. A. & Lemmon, E. M. Anchored hybrid enrichment for massively high-throughput phylogenomics. Syst. Biol. 61, 727–744 (2012).Article 
CAS 
PubMed 

Google Scholar 
Simon, C. et al. Off-target capture data, endosymbiont genes and morphology reveal a relict lineage that is sister to all other singing cicadas. Biol. J. Linn. Soc. Lond. https://doi.org/10.1093/biolinnean/blz120 (2019).Article 

Google Scholar 
Owen, C. L. et al. Detecting and removing sample contamination in phylogenomic data: An example and its implications for Cicadidae phylogeny (Insecta: Hemiptera). Syst. Biol. 71, 1504–1523 (2022).Article 
PubMed 

Google Scholar 
Bushnell, B. BBMap: A Fast, Accurate, Splice-Aware Aligner. Report No.: LBNL-7065E. https://www.osti.gov/biblio/1241166-bbmap-fast-accurate-splice-aware-aligner (Lawrence Berkeley National Lab. (LBNL), 2014).Bolger, A. M., Lohse, M. & Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120 (2014).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Bankevich, A. et al. SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 19, 455–477 (2012).Article 
MathSciNet 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Nurk, S. et al. Assembling genomes and mini-metagenomes from highly chimeric reads. In Research in Computational Molecular Biology 158–170 (Springer, 2013).Chapter 

Google Scholar 
Łukasik, P. et al. One hundred mitochondrial genomes of cicadas. J. Hered. 110, 247–256 (2019).Article 
PubMed 

Google Scholar 
Li, H. & Durbin, R. Fast and accurate short read alignment with Burrows–Wheeler transform. Bioinformatics 25, 1754–1760 (2009).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Katoh, K., Rozewicki, J. & Yamada, K. D. MAFFT online service: Multiple sequence alignment, interactive sequence choice and visualization. Brief. Bioinform. https://doi.org/10.1093/bib/bbx108 (2017).Article 
PubMed Central 

Google Scholar 
Kearse, M. et al. Geneious Basic: An integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics 28, 1647–1649 (2012).Article 
PubMed 
PubMed Central 

Google Scholar 
Hahn, C., Bachmann, L. & Chevreux, B. Reconstructing mitochondrial genomes directly from genomic next-generation sequencing reads—A baiting and iterative mapping approach. Nucleic Acids Res. 41, e129 (2013).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Stamatakis, A. RAxML version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 30, 1312–1313 (2014).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Miller, M. A., Pfeiffer, W., Schwartz, T. Creating the CIPRES Science Gateway for inference of large phylogenetic trees. 2010 Gateway Computing Environments Workshop (GCE) 1–8 (2010).Buckley, T. R., Cordeiro, M., Marshall, D. C. & Simon, C. Differentiating between hypotheses of lineage sorting and introgression in New Zealand alpine cicadas (Maoricicada Dugdale). Syst. Biol. 55, 411–425 (2006).Article 
PubMed 

Google Scholar 
Marshall, D. C., Slon, K., Cooley, J. R., Hill, K. B. R. & Simon, C. Steady Plio-Pleistocene diversification and a 2-million-year sympatry threshold in a New Zealand cicada radiation. Mol. Phylogenet. Evol. 48, 1054–1066 (2008).Article 
PubMed 

Google Scholar 
Bator, J., Marshall, D. C., Leston, A., Cooley, J. & Simon, C. Phylogeography of the endemic red-tailed cicadas of New Zealand (Hemiptera: Cicadidae: Rhodopsalta): Molecular, morphological and bioacoustical confirmation of the existence of Hudson’s Rhodopsalta microdora. Zool. J. Linn. Soc. 195, 1219–1244 (2022).Article 

Google Scholar 
Brumfield, K. D. et al. Gut microbiome insights from 16S rRNA analysis of 17-year periodical cicadas (Hemiptera: Magicicada spp.) Broods II, VI, and X. Sci. Rep. 12, 16967. https://doi.org/10.1038/s41598-022-20527-7 (2022).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Rakitov, R. A. Structure and function of the Malpighian tubules, and related behaviors in juvenile cicadas: Evidence of homology with spittlebugs (Hemiptera: Cicadoidea & Cercopoidea). Zool. Anz. 241, 117–130 (2002).Article 

Google Scholar 
Andersen, P. C., Brodbeck, B. V. & Mizell, R. F. Feeding by the leafhopper, Homalodisca coagulata, in relation to xylem fluid chemistry and tension. J. Insect Physiol. 38, 611–622 (1992).Article 
CAS 

Google Scholar 
Cheung, W. W. K. & Marshall, A. T. Water and ion regulation in cicadas in relation to xylem feeding. J. Insect Physiol. 19, 1801–1816 (1973).Article 
CAS 

Google Scholar 
Williams, K. S. & Simon, C. The ecology, behavior, and evolution of periodical cicadas. Annu. Rev. Entomol. 40, 269–295 (1995).Article 
CAS 

Google Scholar 
Logan, D. P., Rowe, C. A. & Maher, B. J. Life history of chorus cicada, an endemic pest of kiwifruit (Cicadidae: Homoptera). N. Z. Entomol. 37, 96–106 (2014).Article 

Google Scholar 
Buckley, T. R. & Simon, C. Evolutionary radiation of the cicada genus Maoricicada Dugdale (Hemiptera: Cicadoidea) and the origins of the New Zealand alpine biota. Biol. J. Linn. Soc. Lond. 91, 419–435 (2007).Article 

Google Scholar 
Banker, S. E., Wade, E. J. & Simon, C. The confounding effects of hybridization on phylogenetic estimation in the New Zealand cicada genus Kikihia. Mol. Phylogenet. Evol. 116, 172–181 (2017).Article 
PubMed 

Google Scholar 
Groussin, M. et al. Unraveling the processes shaping mammalian gut microbiomes over evolutionary time. Nat. Commun. 8, 14319 (2017).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Sanders, J. G. et al. Stability and phylogenetic correlation in gut microbiota: Lessons from ants and apes. Mol. Ecol. 23, 1268–1283 (2014).Article 
PubMed 

Google Scholar 
Wang, J. et al. Analysis of intestinal microbiota in hybrid house mice reveals evolutionary divergence in a vertebrate hologenome. Nat. Commun. 6, 6440 (2015).Article 
CAS 
PubMed 

Google Scholar 
Brucker, R. M. & Bordenstein, S. R. The hologenomic basis of speciation. Science 466, 667–669 (2013).Article 

Google Scholar 
Chandler, J. A. & Turelli, M. Comment on “The hologenomic basis of speciation: Gut bacteria cause hybrid lethality in the genus Nasonia”. Science 345, 1011 (2014).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Li, Z. et al. Changes in the rumen microbiome and metabolites reveal the effect of host genetics on hybrid crosses. Environ. Microbiol. Rep. 8, 1016–1023 (2016).Article 
CAS 
PubMed 

Google Scholar 
Weintraub, P. G. & Beanland, L. Insect vectors of phytoplasmas. Annu. Rev. Entomol. 51, 91–111 (2006).Article 
CAS 
PubMed 

Google Scholar 
Hopkins, D. L. Xylella fastidiosa: Xylem-limited bacterial pathogen of plants. Annu. Rev. Phytopathol. 27, 271–290 (1989).Article 

Google Scholar 
Karban, R. Why cicadas (Hemiptera: Cicadidae) develop so slowly. Biol. J. Linn. Soc. Lond. 135, 291–298 (2021).Article 

Google Scholar 
Krell, R. K., Boyd, E. A., Nay, J. E., Park, Y.-L. & Perring, T. M. Mechanical and insect transmission of Xylella fastidiosa to Vitis vinifera. Am. J. Enol. Vitic. 58, 211–216 (2007).Article 
CAS 

Google Scholar 
Paião, F., Meneguim, A. M., Casagrande, E. C., Lovato, L. & Leite, R. P. Levantamento de espécies de cigarras e transmissão de Xylella fastidiosa em cafeeiro. http://www.sbicafe.ufv.br/handle/123456789/1457 (2003).Elbeaino, T. et al. Identification of three potential insect vectors of Xylella fastidiosa in southern Italy. Phytopathol. Mediterr. 53, 328–332 (2014).
Google Scholar  More

Thornton, I. Island Colonization: The Origin and Development of Island Communities (Cambridge University Press, 2007).Book 

Google Scholar 
Weigelt, P., Jetz, W. & Kreft, H. Bioclimatic and physical characterization of the world’s islands. Proc. Natl Acad. Sci. 110, 15307–15312 (2013).Article 
PubMed 
PubMed Central 

Google Scholar 
Vitousek, P., Adsersen, H. & Loope, L. Introduction. In Islands: Biological Diversity and Ecosystem Function (eds Vitousek, P. et al.) 1–6 (Berlin, 1995).Chapter 

Google Scholar 
Whittaker, R. J. & Fernández-Palacios, J. M. Island Biogeography: Ecology, Evolution, and Conservation (Oxford University Press, 2007).
Google Scholar 
Lomolino, M., Brown, J. & Sax, D. Island biogeography theory. In The Theory of Island Biogeography Revisited (eds Losos, J. & Ricklefs, R.) 13–51 (Princeton University Press, 2010).
Google Scholar 
Colom, P., Carreras, D. & Stefanescu, C. Long-term monitoring of Menorcan butterfly populations reveals widespread insular biogeographical patterns and negative trends. Biodivers. Conserv. 28, 1837–1851 (2019).Article 

Google Scholar 
Preston, F. W. The canonical distribution of commonness and rarity, part II. Ecology 43, 410–432 (1962).Article 

Google Scholar 
Rosenzweig, M. L. Species Diversity in Space and Time (Cambridge University Press, 1995).Book 

Google Scholar 
Drakare, S., Lennon, J. J. & Hillebrand, H. The imprint of the geographical, evolutionary and ecological context on species–area relationships. Ecol. Lett. 9(2), 215–227 (2006).Article 
PubMed 

Google Scholar 
Field, R. et al. Spatial species-richness gradients across scales: A meta-analysis. J. Biogeogr. 36, 132–147 (2009).Article 

Google Scholar 
Brehm, G., Süssenbach, D. & Fiedler, K. Unique elevational diversity patterns of geometrid moths in an Andean montane rainforest. Ecography 26, 456–466 (2003).Article 

Google Scholar 
Brehm, G., Colwell, R. K. & Kluge, J. The role of environment and mid-domain effect on moth species richness along a tropical elevational gradient. Glob. Ecol. Biogeogr. 16, 205–219 (2007).Article 

Google Scholar 
Beck, J. & Kitching, I. J. Drivers of moth species richness on tropical altitudinal gradients: A cross-regional comparison. Glob. Ecol. Biogeogr. 18, 361–371 (2009).Article 

Google Scholar 
Ashton, L. A. et al. Altitudinal patterns of moth diversity in tropical and subtropical A ustralian rainforests. Aust. Ecol. 41, 197–208 (2016).Article 

Google Scholar 
Maunsell, S. C. et al. Elevational turnover in the composition of leaf miners and their interactions with host plants in Australian subtropical rainforest. Aust. Ecol. 41, 238–247 (2016).Article 

Google Scholar 
McCain, C. M. Global analysis of reptile elevational diversity. Glob. Ecol. Biogeogr. 19, 541–553 (2010).
Google Scholar 
Yu, X. D., Lü, L., Luo, T. H. & Zhou, H. Z. Elevational gradient in species richness pattern of epigaeic beetles and underlying mechanisms at east slope of Balang Mountain in Southwestern China. PLoS ONE 8, e69177 (2013).Article 
PubMed 
PubMed Central 

Google Scholar 
Beck, J. et al. Elevational species richness gradients in a hyperdiverse insect taxon: A global meta-study on geometrid moths. Glob. Ecol. Biogeogr. 26, 412–424 (2017).Article 

Google Scholar 
Szewczyk, T. & McCain, C. M. A systematic review of global drivers of ant elevational diversity. PLoS ONE 11, e0155404 (2016).Article 
PubMed 
PubMed Central 

Google Scholar 
Rahbek, C. The elevational gradient of species richness: A uniform pattern?. Ecography 18, 200–205 (1995).Article 

Google Scholar 
Vitousek, P. M. Oceanic islands as model systems for ecological studies. J. Biogeogr. 29, 573–582 (2002).Article 

Google Scholar 
Kidane, Y. O., Steinbauer, M. J. & Beierkuhnlein, C. Dead end for endemic plant species? A biodiversity hotspot under pressure. Global Ecol. Conserv. 19, e00670 (2019).Article 

Google Scholar 
Meyer, W. M. III. et al. Ground-dwelling arthropod communities of a sky island mountain range in Southeastern Arizona, USA: Obtaining a baseline for assessing the effects of climate change. PLoS ONE 10, e0135210 (2015).Article 
PubMed 
PubMed Central 

Google Scholar 
Kong, W. S. Biogeography of Korean plants 335 (Geobook, 2007) (in Korean).
Google Scholar 
Kitching, R. L. et al. Moth assemblages as indicators of environmental quality in remnants of upland Australian rain forest. J. Appl. Ecol. 37, 284–297 (2000).Article 

Google Scholar 
Froidevaux, J. S., Broyles, M. & Jones, G. Moth responses to sympathetic hedgerow management in temperate farmland. Agric. Ecosyst. Environ. 270, 55–64 (2019).Article 
PubMed 

Google Scholar 
Fox, R. The decline of moths in Great Britain: A review of possible causes. Insect Conserv. Div. 6, 5–19 (2013).Article 

Google Scholar 
Keret, N. M., Mutanen, M. J., Orell, M. I., Itämies, J. H. & Välimäki, P. M. Climate change-driven elevational changes among boreal nocturnal moths. Oecologia 192, 1085–1098 (2020).Article 
PubMed 
PubMed Central 

Google Scholar 
Wenzel, M., Schmitt, T., Weitzel, M. & Seitz, A. The severe decline of butterflies on western German calcareous grasslands during the last 30 years: A conservation problem. Biol. Conserv. 128, 542–552 (2006).Article 

Google Scholar 
Dirzo, R. et al. Defaunation in the anthropocene. Science 345, 401–406 (2014).Article 
PubMed 

Google Scholar 
Hallmann, C. A. et al. More than 75 percent decline over 27 years in total flying insect biomass in protected areas. PLoS ONE 12, e0185809 (2017).Article 
PubMed 
PubMed Central 

Google Scholar 
Sánchez-Bayo, F. & Wyckhuys, K. A. G. Worldwide decline of the entomofauna: A review of its drivers. Biol. Conserv. 232, 8–27 (2019).Article 

Google Scholar 
Zenker, M. M. et al. Diversity and composition of Arctiinae moth assemblages along elevational and spatial dimensions in Brazilian Atlantic Forest. J. Insect Conserv. 19, 129–140 (2015).Article 

Google Scholar 
Brehm, G. & Fiedler, K. Faunal composition of geometrid moths changes with altitude in an Andean montane rain forest. J. Biogeogr. 30, 431–440 (2003).Article 

Google Scholar 
McCain, C. M. & Grytnes, J. A. Elevational gradients in species richness. In Encyclopedia of Life Sciences (Wiley, Chichester, 2010).
Google Scholar 
Heinrich, B. The Hot-Blooded Insects: Strategies and Mechanisms of Thermoregulation 601 (Harvard University Press, 1993).Book 

Google Scholar 
Heinrich, B. Thermoregulation in Endothermic Insects: Body temperature is closely attuned to activity and energy supplies. Science 185, 747–756 (1974).Article 
PubMed 

Google Scholar 
May, M. L. Insect thermoregulation. Annu. Rev. Entomol. 24, 313–349 (1979).Article 

Google Scholar 
Heidrich, L. et al. Noctuid and geometrid moth assemblages show divergent elevational gradients in body size and color lightness. Ecography 44, 1169–1179 (2021).Article 
MathSciNet 

Google Scholar 
Holloway, J. D. Macrolepidoptera diversity in the Indo-Australian tropics, geographic, biotopic and taxonomic variations. Biol. J. Linn. Soc. 30, 325–341 (1987).Article 

Google Scholar 
Axmacher, J. C. et al. Diversity of geometrid moths (Lepidoptera: Geometridae) along an Afrotropical elevational rainforest transect. Divers. Distrib. 10, 293–302 (2004).Article 

Google Scholar 
Heinrich, B. & Mommsen, T. P. Flight of winter moths near 0°C. Science 228, 177–179 (1985).Article 
PubMed 

Google Scholar 
Rydell, J. & Lancaster, W. C. Flight and thermoregulation in moths were shaped by predation from bats. Oikos 88, 13–18 (2000).Article 

Google Scholar 
Skou, P. The geometroid moths of North Europe. Entomonograph, Vol. 6. Brill, Leiden. (1986).Zahiri, R. et al. Molecular phylogenetics of Erebidae (Lepidoptera, noctuoidea). Syst. Entomol. 37, 102–124 (2012).Article 

Google Scholar 
Fiedler, K., Brehm, G., Hilt, N., Sussenbach, D. & Hauser, C. L. Variation of diversity patterns across moth families along a tropical altitudinal gradient. Ecol. Stud. 198, 167–179 (2008).Article 

Google Scholar 
Longino, J. T. & Colwell, R. K. Density compensation, species composition, and richness of ants on a neotropical elevational gradient. Ecosphere 2, 1–20 (2011).Article 

Google Scholar 
Beck, J. & Chey, V. K. Explaining the elevational diversity pattern of geometrid moths from Borneo: A test of five hypotheses. J. Biogeogr. 35, 1452–1464 (2008).Article 

Google Scholar 
Nogués-Bravo, D., Araújo, M. B., Romdal, T. & Rahbek, C. Scale effects and human impact on the elevational species richness gradients. Nature 453, 216–219 (2008).Article 
PubMed 

Google Scholar 
Kwon, T. S. Ants foraging on grasses in South Korea: High diversity in Jeju Island and negative correlation with aphids. J. Asia-Pac. Biodivers. 10, 465–471 (2017).Article 

Google Scholar 
Han, E. K. et al. A disjunctive marginal edge of evergreen broad-leaved oak (Quercus gilva) in East Asia: The high genetic distinctiveness and unusual diversity of Jeju island populations and insight into a massive, independent postglacial colonization. Genes 11, 1114 (2020).Article 
PubMed 
PubMed Central 

Google Scholar 
Chi, Y., Shi, H., Wang, Y., Guo, Z. & Wang, E. Evaluation on island ecological vulnerability and its spatial heterogeneity. Mar. Pollut. Bull. 125, 216–241 (2017).Article 
PubMed 

Google Scholar 
Vehviläinen, H., Koricheva, J. & Ruohomäki, K. Tree species diversity influences herbivore abundance and damage: Meta-analysis of long-term forest experiments. Oecologia 152, 287–298 (2007).Article 
PubMed 

Google Scholar 
Root, R. B. Organization of plant–arthropod association in simple and diverse habitats: The fauna of collards (I. Brassica oleracea). Ecol. Monogr. 43, 95–124 (1973).Article 

Google Scholar 
Otway, S. J., Hector, A. & Lawton, J. H. Resource dilution effects on specialist insect herbivores in a grassland biodiversity experiment. J. Anim. Ecol. 74, 234–240 (2005).Article 

Google Scholar 
Hawkins, B. A. et al. Energy, water, and broad-scale geographic patterns of species richness. Ecology 84, 3105–3117 (2003).Article 

Google Scholar 
Qian, H. Environment–richness relationships for mammals, birds, reptiles, and amphibians at global and regional scales. Ecol. Res. 25, 629–637 (2010).Article 

Google Scholar 
Major, J. A climatic index to vascular plant activity. Ecology 44, 485–498 (1963).Article 

Google Scholar 
Latham, R. E. & Ricklefs, R. E. Global patterns of tree species richness in moist forests: Energy-diversity theory does not account for variation in species richness. Oikos 67, 325–333 (1993).Article 

Google Scholar 
Francis, A. P. & Currie, D. J. A globally consistent richness-climate relationship for angiosperms. Am. Nat. 161, 523–536 (2003).Article 
PubMed 

Google Scholar 
Storch, D. et al. Energy, range dynamics and global species richness patterns: Reconciling mid-domain effects and environmental determinants of avian diversity. Ecol. Lett. 9, 1308–1320 (2006).Article 
PubMed 

Google Scholar 
Intachat, J., Holloway, J. D. & Staines, H. Effects of weather and phenology on the abundance and diversity of geometroid moths in a natural Malaysian tropical rain forest. J. Trop. Ecol. 17, 411–429 (2001).Article 

Google Scholar 
Choi, S. W. Effects of weather factors on the abundance and diversity of moths in a temperate deciduous mixed forest of Korea. Zool. Sci. 25, 53–58 (2008).Article 

Google Scholar 
Kreft, H. & Jetz, W. Global patterns and determinants of vascular plant diversity. Proc. Natl Acad. Sci. 104, 5925–5930 (2007).Article 
PubMed 
PubMed Central 

Google Scholar 
Lennon, J. J., Koleff, P., Greenwood, J. J. D. & Gaston, K. J. The geographical structure of British bird distributions: Diversity, spatial turnover and scale. J. Anim. Ecol. 70, 966–979 (2001).Article 

Google Scholar 
Choi, S. W. A high mountain moth assemblage quickly recovers after fire. Ann. Entomol. Soc. Am. 111, 304–311 (2018).
Google Scholar 
van Swaay, C., Warren, M. & Loïs, G. Biotope use and trends of European butterflies. J. Insect Conserv. 10, 189–209 (2006).Article 

Google Scholar 
De Frenne, P. et al. Microclimate moderates plant responses to macroclimate warming. Proc. Natl Acad. Sci. 110, 18561–18565 (2013).Article 
PubMed 
PubMed Central 

Google Scholar 
Seibold, S. et al. Arthropod decline in grasslands and forests is associated with landscape-level drivers. Nature 574, 671–674 (2019).Article 
PubMed 

Google Scholar 
Conrad, K. F., Warren, M. S., Fox, R., Parsons, M. S. & Woiwod, I. P. Rapid declines of common, widespread British moths provide evidence of an insect biodiversity crisis. Biol. Conserv. 132, 279–291 (2006).Article 

Google Scholar 
White, E. R. Minimum time required to detect population trends: The need for long-term monitoring programs. Bioscience 69, 40–46 (2019).Article 

Google Scholar 
Forister, M. L., Pelton, E. M. & Black, S. H. Declines in insect abundance and diversity: We know enough to act now. Conserv. Sci. Pract. 1, e80 (2019).
Google Scholar 
Didham, R. K. et al. Interpreting insect declines: Seven challenges and a way forward. Insect Conserv. Div. 13, 103–114 (2020).Article 

Google Scholar 
Kim, J. W., Boo, K. O., Choi, J. T. & Byun, Y. H. Climate Change of 100 Years on the Korean Peninsula (National Institute of Meteorological Science, 2018).
Google Scholar 
Kim, S. S., Beljaev, E. A. & Oh, S. H. Illustrated Catalogue of Geometridae in Korea (Lepidoptera: Geometrinae, Ennominae) (Korea Research Institute of Bioscience and Biotechnology & Center for Insect Systematics, 2001).
Google Scholar 
Kononenko, V.S., Ahn, S.B. & Ronkay, L. Illustrated catalogue of Noctuidae in Korea (Lepidoptera). Insects of Korea 3. KRIBB & CIS, Junghaengsa (1998).Shin, Y.H. Coloured illustrations of the moths of Korea. Academybook (2001).Kim, S.S., Choi, S.W., Sohn, J.C., Kim, T. & Lee, B.W. The Geometrid moths of Korea (Lepidoptera: Geometridae). Junghaengsa (2016).Kim, C. G. & Kim, N. W. Altitudinal pattern of evapotranspiration and water need for upland crops in Jeju Island. J. Korea Water Resour. Assoc. 48, 915–923 (2015).Article 

Google Scholar 
Zuur, A. F., Ieno, E. N. & Elphick, C. S. A protocol for data exploration to avoid common statistical problems. Methods Ecol. Evol. 1, 3–14 (2010).Article 

Google Scholar 
Magurran, A. E. Ecological Diversity and its Measurement (Princeton University Press, 1988).Book 

Google Scholar 
Hammer, Ø., Harper, D. A. & Ryan, P. D. PAST: Paleontological statistics software package for education and data analysis. Palaeontol. Electron. 4, 9 (2001).
Google Scholar 
Burnham, K. P. & Anderson, D. R. Model Selection and Multimodel Inference: A Practical Information-Theoretic Approach 2nd edn. (Springer, 2002).MATH 

Google Scholar 
R Development Core Team. R 4.0.3. R: A language and environment for statistical computing. R Foundation for statistical computing Vienna. Austria. URL http://www.R-project.org. (2020).Pohlert, T. Non-parametric trend tests and change-point detection. R-package version 0.0.1. (2020).Hipel, K. W. & McLeod, A. I. Time Series Modelling of Water Resources and Environmental Systems (Elsevier, 1994).
Google Scholar 
Chao, A., Chazdon, R. L., Colwell, R. K. & Shen, T. J. Abundance-based similarity indices and their estimation when there are unseen species in samples. Biometrics 62, 361–371 (2006).Article 
MathSciNet 
PubMed 
MATH 

Google Scholar 
Colwell, R. K. Estiamtes, Version 91: Statistical Estimation of Species Richness and Shared Species from Samples (University of Connecticut, 2013).
Google Scholar 
Baselga, A. Partitioning the turnover and nestedness components of beta diversity. Global Ecol. Biogeogr. 19, 134–143 (2010).Article 

Google Scholar 
Baselga, A. The relationship between species replacement, dissimilarity derived from nestedness, and nestedness. Glob. Ecol. Biogeogr. 21, 1223–1232 (2012).Article 

Google Scholar  More

Time using the rock as refugeTemperature had an effect in the refuge usage of both species when analysed together (lme.zig: F3,192 = 7.97, p = 0.0001; Fig. 1A). However, species behaved differently (lme.zig: F1,192 = 14.79, p = 0.0004; Fig. 1A). As hypothesised, there was an interaction between temperature and species (lme.zig: F3,192 = 11.90, p  0.14, Fig. 1B).Size had an effect in the time exploring the rock (lme: F1,192 = 6.91, p = 0.012, Fig. 3) when species were analysed together, but there was no interaction with temperatures (lme: F3,192 = 0.42, p = 0.74, Fig. 3). We found that the interaction between species and size was close to be significant (lme: F1,192 = 3.62, p = 0.064, Fig. 3), implying that possibly smaller fish spent more time exploring the rock than bigger fish. However, when analysed separately, we did not find an effect of size in the exploring behaviour neither for twoline skiffias (lme: F1,96 = 2.99, p = 0.099, Fig. 3) nor for guppies (lme: F1,96 = 0.33, p = 0.569, Fig. 3).Figure 3Proportion of the total time observed (600 s) fish of different sizes spent exploring the rock. Lines represent the areas where the density of data is higher.Full size imageTime spent swimmingTemperature had an effect in the time spent swimming for both species when analysed together (lme: F3,192 = 23.48, p  More

MacPherson, E. Ontogenetic shifts in habitat use and aggregation in juvenile sparid fishes. J. Exp. Mar. Bio. Ecol. 220, 127–150 (1998).Article 

Google Scholar 
Freitas, C., Olsen, E. M., Knutsen, H., Albretsen, J. & Moland, E. Temperature-associated habitat selection in a cold-water marine fish. J. Anim. Ecol. 85, 628–637 (2016).Article 
PubMed 

Google Scholar 
Michelot, C. et al. Seasonal variation in coastal marine habitat use by the European shag: Insights from fine scale habitat selection modeling and diet. Deep. Res. Part II Top. Stud. Oceanogr. 141, 224–236 (2017).Article 

Google Scholar 
Davoren, G. K., Montevecchi, W. A. & Anderson, J. T. Distributional patterns of a marine bird and its prey: Habitat selection based on prey and conspecific behaviour. Mar. Ecol. Prog. Ser. 256, 229–242 (2003).Article 

Google Scholar 
Chiarello, A. G. et al. A translocation experiment for the conservation of maned sloths, Bradypus torquatus (Xenarthra, Bradypodidae). Biol. Conserv. 118, 421–430 (2004).Article 

Google Scholar 
Fukuda, Y. et al. Environmental resistance and habitat quality influence dispersal of the saltwater crocodile. Mol. Ecol. 31, 1076–1092 (2022).Article 
PubMed 

Google Scholar 
O’Leary, S. J., Dunton, K. J., King, T. L., Frisk, M. G. & Chapman, D. D. Genetic diversity and effective size of Atlantic sturgeon, Acipenser oxyrhinchus oxyrhinchus river spawning populations estimated from the microsatellite genotypes of marine-captured juveniles. Conserv. Genet. 15, 1173–1181 (2014).Article 

Google Scholar 
Brüniche-Olsen, A. et al. Genetic data reveal mixed-stock aggregations of gray whales in the North Pacific Ocean. Biol. Lett. 14, 1–4 (2018).Article 

Google Scholar 
Carroll, E. L. et al. Genetic diversity and connectivity of southern right whales (Eubalaena australis) found in the Brazil and Chile-Peru wintering grounds and the South Georgia (Islas Georgias Del Sur) feeding ground. J. Hered. 111, 263–276 (2020).Article 
PubMed 
PubMed Central 

Google Scholar 
Bowen, A. B. W. et al. Origin of hawksbill turtles in a Caribbean feeding area as indicated by genetic markers. Ecol. Appl. 6, 566–572 (1996).Article 

Google Scholar 
Paxton, K. L., Yau, M., Moore, F. R. & Irwin, D. E. Differential migratory timing of western populations of Wilson’s Warbler (Cardellina pusilla) revealed by mitochondrial DNA and stable isotopes. Auk 130, 689–698 (2013).Article 

Google Scholar 
Anderson, E. C., Waples, R. S. & Kalinowski, S. T. An improved method for predicting the accuracy of genetic stock identification. Can. J. Fish. Aquat. Sci. 65, 1475–1486 (2008).Article 

Google Scholar 
Debevec, E. M. SPAM (version 3.2): Statistics program for analyzing mixtures. J. Hered. 91, 509–511 (2000).Article 
PubMed 

Google Scholar 
Bolker, B. M., Okuyama, T., Bjorndal, K. A. & Bolten, A. B. Incorporating multiple mixed stocks in mixed stock analysis: ‘Many-to-many’ analyses. Mol. Ecol. 16, 685–695 (2007).Article 
PubMed 

Google Scholar 
Neaves, P. I., Wallace, C. G., Candy, J. R. & Beacham, T. D. CBayes: Computer Program for Mixed Stock Analysis of Allelic Data. Free Program Distributed by the Authors Over the Internet. at (2005).Pella, J. & Masuda, M. Bayesian methods for analysis of stock mixtures from genetic characters. Fish. Bull. 99, 151–167 (2001).
Google Scholar 
Bolker, B., Okuyama, T., Bjorndal, K. A. & Bolten, A. B. Sea turtle stock estimation using genetic markers: Accounting for sampling error of rare genotypes. Ecol. Appl. 13, 763–775 (2003).Article 

Google Scholar 
Okuyama, T. & Bolker, B. M. Combining genetic and ecological data to estimate sea turtle origins. Ecol. Appl. 15, 315–325 (2005).Article 

Google Scholar 
Nishizawa, H. et al. Composition of green turtle feeding aggregations along the Japanese archipelago: Implications for changes in composition with current flow. Mar. Biol. 160, 2671–2685 (2013).Article 

Google Scholar 
Naro-Maciel, E. et al. Predicting connectivity of green turtles at Palmyra Atoll, central Pacific: A focus on mtDNA and dispersal modelling. J. R. Soc. Interface 11, 20130888 (2014).Proietti, M. C. et al. Green turtle Chelonia mydas mixed stocks in the western South Atlantic, as revealed by mtDNA haplotypes and drifter trajectories. Mar. Ecol. Prog. Ser. 447, 195–209 (2012).Article 

Google Scholar 
van der Zee, J. P. et al. Population recovery changes population composition at a major southern Caribbean juvenile developmental habitat for the green turtle, Chelonia mydas. Sci. Rep. 9, 1–11 (2019).
Google Scholar 
Shamblin, B. M. et al. Mexican origins for the Texas green turtle foraging aggregation: A cautionary tale of incomplete baselines and poor marker resolution. J. Exp. Mar. Bio. Ecol. 488, 111–120 (2017).Article 

Google Scholar 
Seminoff, J. A. et al. Status Review of the Green Turtle (Chelonia mydas) Under the Endangered Species Act. (NOAA Technical Memorandum, NOAA-NMFS-SWFSC, 2015).Chaloupka, M. et al. Encouraging outlook for recovery of a once severely exploited marine megaherbivore. Glob. Ecol. Biogeogr. 17, 297–304 (2008).Article 

Google Scholar 
Bjorndal, K. A. & Bolten, A. B. Annual variation in source contributions to a mixed stock: Implications for quantifying connectivity. Mol. Ecol. 17, 2185–2193 (2008).Article 
PubMed 

Google Scholar 
Roland, J., Keyghobadi, N. & Fownes, S. Alpine Parnassius butterfly dispersal: Effects of landscape and population size. Ecology 81, 1642–1653 (2000).Article 

Google Scholar 
Vanschoenwinkel, B., De Vries, C., Seaman, M. & Brendonck, L. The role of metacommunity processes in shaping invertebrate rock pool communities along a dispersal gradient. Oikos 116, 1255–1266 (2007).Article 

Google Scholar 
Shamblin, B. M. et al. Mitogenomic sequences better resolve stock structure of southern Greater Caribbean green turtle rookeries. Mol. Ecol. 21, 2330–2340 (2012).Article 
PubMed 

Google Scholar 
Witherington, B., Hirama, S. & Hardy, R. Young sea turtles of the pelagic Sargassum-dominated drift community: Habitat use, population density, and threats. Mar. Ecol. Prog. Ser. 463, 1–22 (2012).Article 

Google Scholar 
Putman, N. F. & Mansfield, K. L. Direct evidence of swimming demonstrates active dispersal in the sea turtle ‘lost years’. Curr. Biol. 25, 1221–1227 (2015).Article 
PubMed 

Google Scholar 
Mansfield, K. L., Wyneken, J. & Luo, J. First Atlantic satellite tracks of ‘lost years’ green turtles support the importance of the Sargasso Sea as a sea turtle nursery. Proc. R. Soc. B Biol. Sci. 288, 20210057 (2021).Putman, N. F. et al. Predicted distributions and abundances of the sea turtle ‘lost years’ in the western North Atlantic Ocean. Ecography (Cop.) 43, 506–517 (2020).Article 

Google Scholar 
Putman, N. F. & Naro-Maciel, E. Finding the ‘lost years’ in green turtles: Insights from ocean circulation models and genetic analysis. Proc. R. Soc. B Biol. Sci. 280, 20131468 (2013).Naro-Maciel, E., Hart, K. M., Cruciata, R. & Putman, N. F. DNA and dispersal models highlight constrained connectivity in a migratory marine megavertebrate. Ecography (Cop.) 40, 586–597 (2017).Article 

Google Scholar 
Ehrhart, L. M., Redfoot, W. E. & Bagley, D. A. Marine turtles of the central region of the Indian River Lagoon system, Florida. Florida Sci. 70, 415–434 (2007).
Google Scholar 
Redfoot, W. & Ehrhart, L. Trends in size class distribution, recaptures, and abundance of juvenile green turtles (Chelonia mydas) utilizing a rock riprap lined embayment at Port Canaveral, Florida, USA, as developmental habitat. Chelonian Conserv. Biol. 12, 252–261 (2013).Article 

Google Scholar 
Ehrhart, L., Redfoot, W., Bagley, D. & Mansfield, K. Long-term trends in loggerhead (Caretta caretta) nesting and reproductive success at an important western Atlantic rookery. Chelonian Conserv. Biol. 13, 173–181 (2014).Article 

Google Scholar 
Bolten, A. B. Techniques for measuring sea turtles. in Research and Management Techniques for the Conservation of Sea Turtles. (eds. Eckert, K. L., Bjorndal, K. A., Abreu-Grobois, F. A. & Donnelly, M.). 1–5 (1999).Bagley, D. A. Characterizing Juvenile Green Turtles, (Chelonia mydas), from Three East Central Florida Developmental Habitats. (University of Central Florida, 2003).Rohland, N. & Reich, D. Cost-effective, high-throughput DNA sequencing libraries for multiplexed target capture. Genome Res. 22, 939–946 (2012).Article 
PubMed 
PubMed Central 

Google Scholar 
Faircloth, B. & Glenn, T. Preparation of an AMPure XP Substitute. AKA Serapure https://doi.org/10.6079/J9MW2F26 (2016).Article 

Google Scholar 
Abreu-Grobois, F. A. et al. New mtDNA Dloop primers which work for a variety of marine turtle species may increase the resolution of mixed stock analyses. in Proceedings of the 26th Annual Symposium on Sea Turtle Biology. 179 (International Sea Turtle Society, 2006).Kearse, M. et al. Geneious basic: An integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics 28, 1647–1649 (2012).Article 
PubMed 
PubMed Central 

Google Scholar 
Leigh, J. W. & Bryant, D. PopART: Full-feature software for haplotype network construction. Methods Ecol. Evol. 6, 1110–1116 (2015).Article 

Google Scholar 
Excoffier, L. & Lischer, H. E. L. Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Windows. Mol. Ecol. Resour. 10, 564–567 (2010).Article 
PubMed 

Google Scholar 
Wright, S. Evolution and the Genetics of Populations. Vol. 4. Variability Within and Among Natural Populations. (University of Chicago Press, 1978).Hays, G. C. Ocean currents and marine life. Curr. Biol. 27, R470–R473 (2017).Article 
PubMed 

Google Scholar 
Engstrom, T. N., Meylan, P. A. & Meylan, A. B. Origin of juvenile loggerhead turtles (Caretta caretta) in a tropical developmental habitat in Caribbean Panamá. Anim. Conserv. 5, 125–133 (2002).Article 

Google Scholar 
Florida Fish and Wildlife Conservation Commission-Fish and Wildlife Research Institute, F. W. C. F. W. R. I. Index Nesting Beach Survey (INBS). (2021).Cuevas Flores, E. A., Guzmán Hernández, V., Guerra Santos, J. J. & Rivas Hernández, G. A. El uso del Conocimiento de las Tortugas Marinas Como Herramienta para la Restauración de sus Poblaciones y Hábitats Asociados. (Universidad Autónoma del Carmen, 2019).Pineda, O. G. & Rocha, A. R. B. Las Tortugas Marinas en México: Logros y Perspectivas para su Conservación. (CONANP, 2016).Varela, R. G., Quílez, G. Z. & Harrison, E. Report on the 2014 Green Turtle Program at Tortuguero, Costa Rica. (2015).Azanza Ricardo, J. et al. Nesting ecology of Chelonia mydas (Testudines: Cheloniidae) on the Guanahacabibes Peninsula. Cuba. Rev. Biol. Trop. 61, 1935–1945 (2013).PubMed 

Google Scholar 
Nalovic, M. A. et al. Sea Turtles in the North Atlantic & Wider Caribbean Region. (2020).Makowski, D., Ben-Shachar, M. & Lüdecke, D. bayestestR: Describing effects and their uncertainty, existence and significance within the Bayesian framework. J. Open Source Softw. 4, 1541 (2019).Article 

Google Scholar 
Kruschke, J. K. Doing Bayesian Data Analysis: A Tutorial with R, JAGS, and Stan. https://doi.org/10.1016/B978-0-12-405888-0.09999-2 (Academic Press, 2015).Ruiz-Urquiola, A. et al. Population genetic structure of greater Caribbean green turtles (Chelonia mydas) based on mitochondrial DNA sequences, with an emphasis on rookeries from southwestern Cuba. Rev. Investig. Mar. 31, 33–52 (2010).
Google Scholar 
Long, C. A. et al. Incongruent long-term trends of a marine consumer and primary producers in a habitat affected by nutrient pollution. Ecosphere 12, e03553 (2021).Article 

Google Scholar 
Phillips, K. F., Stahelin, G. D., Chabot, R. M. & Mansfield, K. L. Long-term trends in marine turtle size at maturity at an important Atlantic rookery. Ecosphere 12, 7 (2021).Article 

Google Scholar 
Bjorndal, K. A., Bolten, A. B. & Chaloupka, M. Y. Evaluating trends in abundance of immature green turtles, Chelonia mydas, in the Greater Caribbean. Ecol. Appl. 15, 304–314 (2005).Article 

Google Scholar 
Naro-Maciel, E. et al. The interplay of homing and dispersal in green turtles: A focus on the southwestern atlantic. J. Hered. 103, 792–805 (2012).Article 
PubMed 

Google Scholar 
Monzón-Argüello, C. et al. Evidence from genetic and Lagrangian drifter data for transatlantic transport of small juvenile green turtles. J. Biogeogr. 37, 1752–1766 (2010).Article 

Google Scholar 
Luke, K., Horrocks, J. A., LeRoux, R. A. & Dutton, P. H. Origins of green turtle (Chelonia mydas) feeding aggregations around Barbados, West Indies. Mar. Biol. 144, 799–805 (2004).Article 

Google Scholar 
Bass, A. L., Epperly, S. P. & Braun-McNeill, J. Green turtle (Chelonia mydas) foraging and nesting aggregations in the Caribbean and Atlantic: Impact of currents and behavior on dispersal. J. Hered. 97, 346–354 (2006).Article 
PubMed 

Google Scholar 
Lahanas, P. N. et al. Genetic composition of a green turtle (Chelonia mydas) feeding ground population: Evidence for multiple origins. Mar. Biol. 130, 345–352 (1998).Article 

Google Scholar 
Foley, A. M. et al. Characteristics of a green turtle (Chelonia mydas) assemblage in northwestern Florida determined during a hypothermic stunning event. Gulf Mex. Sci. 25, 131–143 (2007).
Google Scholar 
Bass, A. L., Lagueux, C. J. & Bowen, B. W. Origin of green turtles, Chelonia mydas, at ‘Sleeping Rocks’ off the Northeast coast of Nicaragua. Copeia 1998, 1064 (1998).Article 

Google Scholar 
Bass, A. L. & Witzell, W. N. Demographic composition of immature green turtles (Chelonia mydas) from the East Central Florida Coast: Evidence from mtDNA markers. Herpetologica 56, 357–367 (2000).
Google Scholar 
Bjorndal, K. A., Parsons, J., Mustin, W. & Bolten, A. B. Threshold to maturity in a long-lived reptile: Interactions of age, size, and growth. Mar. Biol. 160, 607–616 (2013).Article 

Google Scholar 
Perrault, J. R. et al. Maternal health status correlates with nest success of leatherback sea turtles (Dermochelys coriacea) from Florida. PLoS ONE 7, e31841 (2012).Article 
PubMed 
PubMed Central 

Google Scholar 
Montero, N. et al. Warmer and wetter conditions will reduce offspring production of hawksbill turtles in Brazil under climate change. PLoS ONE 13, 1–16 (2018).Article 

Google Scholar 
Shamblin, B. M. et al. Geographic patterns of genetic variation in a broadly distributed marine vertebrate: New insights into loggerhead turtle stock structure from expanded mitochondrial DNA sequences. PLoS ONE 9, 85956 (2014).Article 

Google Scholar 
Anderson, J. D., Shaver, D. J. & Karel, W. J. Genetic Diversity and Natal Origins of Green Turtles (Chelonia mydas) in the Western Gulf of Mexico. J. Herpetol. 47, 251–257 (2013).Article 

Google Scholar  More