Loggerhead Sea Turtles as Hosts of Diverse Bacterial and Fungal Communities

Research on microbial communities associated with wild animals provides a valuable reservoir of knowledge that could be used for enhancing their rehabilitation and conservation. The loggerhead sea turtle (Caretta caretta) is a globally distributed species with its Mediterranean population categorized as least concern according to the IUCN Red List of Threatened Species as a result of robust conservation efforts. In our study, we aimed to further understand their biology in relation to their associated microorganisms. We investigated epi- and endozoic bacterial and endozoic fungal communities of cloaca, oral mucosa, carapace biofilm. Samples obtained from 18 juvenile, subadult, and adult turtles as well as 8 respective enclosures, over a 3-year period, were analysed by amplicon sequencing of 16S rRNA gene and ITS2 region of nuclear ribosomal gene. Our results reveal a trend of decreasing diversity of distal gut bacterial communities with the age of turtles. Notably, Tenacibaculum species show higher relative abundance in juveniles than in adults. Differential abundances of taxa identified as Tenacibaculum, Moraxellaceae, Cardiobacteriaceae, and Campylobacter were observed in both cloacal and oral samples in addition to having distinct microbial compositions with Halioglobus taxa present only in oral samples. Fungal communities in loggerheads’ cloaca were diverse and varied significantly among individuals, differing from those of tank water. Our findings expand the known microbial diversity repertoire of loggerhead turtles, highlighting interesting taxa specific to individual body sites. This study provides a comprehensive view of the loggerhead sea turtle bacterial microbiota and marks the first report of distal gut fungal communities that contributes to establishing a baseline understanding of loggerhead sea turtle holobiont. Supplementary Information The online version contains supplementary material available at 10.1007/s00248-024-02388-x.


Introduction
Microbial communities associated with animal hosts play crucial roles in various aspects of the host's development, physiology, immune response, metabolism, and reproduction, and may have an impact on the host's evolutionary potential [1,2].While the importance of sea turtle epibiosis with macroepibionts (> 1 mm) such as barnacles has long been recognized [3], research on the microbial epibionts and endobionts of sea turtles has only recently gained attention [4][5][6].Sea turtles hold a unique ecological role as keystone species, connecting terrestrial and coastal habitats, but they are also highly vulnerable to anthropogenic threats, such as climate change, disruption of feeding and breeding habitats, egg poaching, and accidental bycatch [7][8][9].To lessen some of the pressure sea turtles face today, global conservation efforts have focused on safeguarding female turtles and their nesting areas, together with the rehabilitation of injured and sick turtles [8].Motivated by the aim of enhancing the rehabilitation and conservation of wild animals and their associated microbiomes, the studies of sea turtle gut, skin, egg, and nest microbiomes have become a forefront in sea turtle conservation research, building upon cultivation-based and pathogen-oriented research in the veterinary domain [10][11][12].
Loggerhead sea turtles (Caretta caretta) are a widely distributed species and are classified as vulnerable by the IUCN Red List of threatened species [13].However, the Mediterranean subpopulation of loggerheads is categorized as "Least Concern" due to long-term successful conservation efforts [8,13].Currently, loggerhead sea turtles' microbiota is the second most studied, preceded only by green turtles (Chelonia mydas) [6].Previous studies that used culture-dependent approaches have identified the most common pathogens associated with mucosal surfaces, skin lesions (such as bacterial families Aeromonadaceae, Pseudomonadaceae, Enterobacteriaceae), and hatchling failure (Fusarium spp.), with the presence of antibiotic resistance genes indicating loggerheads as sentinels of antibiotic pollution in the Mediterranean [11,[14][15][16][17][18][19][20][21][22].Recent investigations using next-generation sequencing approaches to study the loggerhead microbiota shed light on the bacterial community structure and composition of the gastrointestinal tract [5, 23,24], the impact of rehabilitation on mucosal bacteriomes [24,25], the effects of plastic pollution on the gut bacteria [26], variations in microbial communities driven by localities [27,28] or turtle anatomy [28,29], and host-microbial coevolution patterns [30].On the other hand, the fungal communities associated with marine reptiles, including sea turtles, have received limited attention using cultivation-independent approaches, despite the vulnerability of sea turtles to infections caused by Fusarium spp.during early development [18].Recent work by Guo et al. [31] provided initial insights into the fungal communities found on carapace (healthy and ulcerated), in faeces, and in the seawater of green turtle juveniles undergoing rehabilitation; however, comprehensive surveys of endobiotic fungal communities in loggerhead sea turtles have not yet been conducted.Given the ecological significance of loggerhead sea turtles in the Mediterranean basin ecosystem, their role as sentinels for pollution, and their potential to act as vectors for zoonotic diseases, a comprehensive approach including eukaryotic microorganisms is necessary to understand the loggerhead sea turtle microbiota.This knowledge will contribute to the advancement of current conservation practices and future microbial stewardship efforts [32].
The objective of this study was to investigate the bacterial and fungal communities associated with loggerhead sea turtles found in the Adriatic Sea.More specifically, we aimed to analyse the bacterial communities in the cloacal, oral, and enclosure samples, as well as fungal communities of cloacal and enclosure tank water samples using amplicon sequencing targeting the V3-V4 (V34) region of 16S rRNA gene and the ITS2 region of nuclear ribosomal genes, respectively.Furthermore, when available, we aimed to compare the bacterial communities of carapace biofilm samples corresponding to turtles in this study out of which some were previously analysed as a part of our earlier study on sea turtle epibiosis [28].By combining these datasets, we provide a comprehensive overview of the environmental, surface, and internal microbiota of the loggerhead turtles, establishing a baseline for future holobiont approaches to studying the loggerhead sea turtles.

Loggerhead Sea Turtle Sampling
Loggerhead sea turtles investigated in this study were found at various locations along the Adriatic Sea coast from 2019 to 2021 (Fig. 1a) and transported to two locations where the sampling was conducted: at the Sea Turtle Clinic (STC) of the Department of Veterinary Medicine of University of Bari "Aldo Moro" in Italy and the Sea Turtle Rescue Center Aquarium Pula in Croatia.The turtles were sampled immediately upon their arrival to the rehabilitation centres or during/after rehabilitation, and prior to release (Table S1).A total of 18 loggerhead sea turtles and 8 respective enclosures were included in the sampling.Turtles were classified as juveniles (n = 10), subadult (n = 4), and adults (n = 4) according to their size [33] and sex determination was based on observable physical characteristics when possible (Table 1; Fig. 1b).Additional information about sampling procedures and the loggerhead population surveyed in this study can be found in Supplementary Methods and Table S2.
The endozoic samples were collected from cloacal and oral cavities in triplicate by sterile synthetic swabs (Aptaca Nuova) as described in Filek et al. [25].When available, enclosure tank water was collected in sterile containers, vacuum filtered on 0.2-μm sterile Whatman polycarbonate membrane filters (Sigma-Aldrich), and stored in 2-ml tubes in 96% EtOH.All samples were stored at −20 °C until DNA extraction and further processing.Additionally, corresponding epizoic carapace biofilm samples were obtained by randomly brushing the entire carapace using a toothbrush (Dentalux Classic, hard, Lidl) according to [28,34] and the collected material was resuspended in 96% EtOH, and stored at −20 °C (Table 1).Each endozoic sequencing sample ID that is referred to in this manuscript has a 16S or ITS prefix, sampling event number, and suffix corresponding to sampling site: C -cloaca, O -oral cavity, W -tank water; for example, sample ID ITS0084C represents cloacal sample of fungal ITS2 sequences for sampling event 0084 (turtle ID010;     Epizoic samples have a TB prefix and numbering unrelated to sampling event.

DNA Extraction and Sequencing
Total DNA from swabs and filters was extracted using the DNeasy PowerSoil kit (Qiagen) following the manufacturer's instructions with several modifications: (1) the samples were incubated in C1 solution at 65 °C for 1 h; (2) instead of bead beating, PowerBead Tubes were vortexed horizontally for 10 min at maximum speed; and (3) all downstream incubation times at 2-8 °C were increased to 15 min.The methods used for epizoic carapace biofilm samples that were sequenced for V4 region (ID010-ID074) were described by Kanjer et al. [28].The DNA from biofilm scrapings (ID093-ID122) analysed only in this study was extracted via DNeasy PowerLyzer PowerSoil extraction kit (Qiagen) following the manufacturer's instructions and modified as follows: (1) after addition of C1 solution, the samples were incubated at 70 °C for 10 min; (2) bead beating was performed at 30 Hz for 1 min in TissueLyzer Retsch Qiagen; (3) 50 µl of C6 solution was used for DNA elution and incubated for 5 min at room temperature prior to centrifugation.Nuclease-free water (W4502, Sigma-Aldrich) was used as the negative control for the DNA extraction step and was processed using the DNA extraction kit in parallel to all the samples.The quality and quantity of extracted DNA were evaluated by BioSpec-nano (Shimadzu).
The extracted DNA from both endozoic and epizoic samples was stored at −20 °C and sent for Illumina MiSeq v3 300 × 2 bp paired-end sequencing to Microsynth, Switzerland.Primers used for sequencing the V34 region of 16S rRNA gene were 341F and 805R [35], and primers for fungal ITS2 region of the nuclear ribosomal gene were ITS3 and ITS4 [36].

Bioinformatics and Statistics
The obtained sequences had non-biological sequences trimmed by the sequencing facility and checked for quality with FastQC [37].Sequencing data is available at EMBL ENA at accessions PRJEB62752 and PRJEB68298 for 16S rRNA gene, PRJEB62762 for ITS2 region sequences, and from Kanjer et al. [28] PRJEB51458.The sequences were imported and analysed in QIIME 2 (versions 2021.8 and 2023.2) [38].
Statistical analyses were performed within QIIME 2 environment and with R. Alpha diversity indices, including Shannon's entropy, Pielou's evenness, Faith's phylogenetic diversity, and observed ASVs, were calculated via q2-diversity plugin.The Kruskal-Wallis rank-sum test was employed to determine differences between selected groups, followed by post hoc pairwise comparisons using the Wilcoxon rank-sum exact test.For beta diversity analyses, rarefied data (sequencing depth determined by alpha rarefaction curves) were explored using the q2-diversity plugin with Bray-Curtis, Jaccard, unweighted UniFrac, and weighted UniFrac distances [39].Compositional data analysis on non-rarefied datasets was performed by calculating robust Aitchison distances using the q2-deicode plugin or the R package vegan v.2.6-4 [40][41][42].Principal coordinate analysis (PCoA) was conducted on Bray-Curtis, Jaccard, and all UniFrac distances, while principal component analysis was performed for robust Aitchison (rPCA) using q2-diversity and q2-deicode, respectively.To assess the relative impact of factors (age range, sex, and duration of rehabilitation) on microbial communities, a multi-way permutational multivariate analysis of variance (Adonis2 PERMANOVA) with 9999 permutations was employed (Anderson, 2001) in R by using vegan v.2.6-4 and pairwise Adonis v.0.4.1 packages [43].Resulting p-values from all pairwise tests were adjusted using the Benjamini-Hochberg false-discovery rate (FDR) correction for multiple comparisons (reported as q-values).Differential abundance analysis was used to identify differentially abundant (DA) features in sample site pairs by using the ANCOMBC package "ancombc2" function in R [44,45].Default parameters were used and pairwise testing was enabled (with Holm's method for adjusting p-values, reported as q-values), except for DA testing in V4-trimmed sequences where the "struc_zero" was set to "TRUE" to exclude structural zeros based on sampling sites.Log fold change (LFC) indicates the scale of differential abundance between differentially abundant features.Features with p-value < 0.05 were reported as differentially abundant, and q-value < 0.05 as significantly differentially abundant.

Endozoic and Tank Water Bacterial Communities
Altogether, 50 endozoic and water samples (plus one negative control) were sequenced, and 7,110,067 high-quality sequences were obtained (median frequency per sample was 107,522, min.2, max.911,443).Denoising yielded a total of 11,105 ASVs.After filtering mitochondrial and chloroplast sequences, the number of ASVs decreased to 10,946.Forty-three samples yielded enough high-quality reads (at minimum sequencing depth above 10,000) for downstream analyses (19/21 cloacal, 16/20 oral, and 8/9 tank water samples).

Endozoic and Tank Water Fungal Communities
Overall, 29 samples (plus one negative control) were sequenced, and 2,208,729 high-quality sequences were obtained (median frequency per sample was 69,203, min.11,722, max.150,142).Denoising yielded a total of 9547 ASVs.Based on alpha rarefaction curves, 27 samples yielded enough high-quality reads (at minimum sequencing depth at 25,000 reads) for downstream statistical analyses requiring rarefied data (19/20 cloacal and 8/9 tank water samples).All samples except negative control were used in compositional data analyses.

Epizoic and Endozoic Bacterial Communities
After trimming epizoic, endozoic, and tank water sequences to V4 region of 16S rRNA gene (65 samples in total) and denoising, 10,072,866 high-quality sequences were merged across all sequencing events (median frequency per sample was 126,629, min.0, max.946,183) and yielded 13,263 ASVs.Due to reduced resolution after trimming the longer V34 region to shorter V4 region, the number of ASVs for cloacal, oral, and tank water samples was expectedly lower (7725 in V4 vs. 11,105 in V34 sequences).
The taxonomic composition of endozoic samples and tank water using the V4 sequences resembled the one detected with V34 and it is reported in Table S5.Additional information on taxonomic composition of epizoic bacterial communities in carapace samples is reported in Supplementary Results.Differential abundance analysis on ASVs (structural zeros excluded) detected 17 significantly DA ASVs across all sample sites (Fig. 6).A member of Rhodobacteraceae family (v4ASV8861) and Halioglobus sp.(v4ASV11704) was DA abundant in carapace tested against endozoic samples, but not tank water.In tank water, a member of Flavobacteriaceae (v4ASV2446) and NS3a marine group (v4ASV2817) were consistently DA (Fig. 6).Differential abundance analysis on taxa collapsed to species level is reported in Supplementary Results and Fig. S2.

Discussion
This study provides an overview of the bacterial and fungal communities inhabiting oral and cloacal environments of loggerhead sea turtles.Additionally, we aimed to investigate bacterial communities in carapace samples in relation to the gastrointestinal tract and tank water to explore possible connections between two habitats.The carapace exhibited the highest bacterial diversity, followed by oral samples, influenced by the tank water environment, and then cloacal samples.Each sampling site had distinct microbial communities and cloacal bacterial diversity negatively correlated with turtle size and age.Conversely, fungal communities in the cloaca were distinct from tank water and showed high heterogeneity among individual turtles, with no discernible patterns related to age or sex.

Cloacal Bacterial Diversity and Structure Changes with the Turtle Age
Similarly to previous studies on cloacal microbiota in Adriatic loggerhead sea turtles, we observed changes in bacterial richness, diversity, and structure, negatively correlating with loggerheads' CCL and, consequently, their age [24].According to research on vertebrate gut microbiomes, it is expected that the bacterial diversity increases with body size in animals with complex digestive systems, like ruminants, and decreases in animals with simple guts, often omnivores or carnivores [49].Loggerhead sea turtles are omnivores and exhibit an ontogenetic shift from oceanic-pelagic habitats as juveniles to neritic-benthic feeding grounds as adults.However, in the Mediterranean Sea, they employ an "amphi-habitat strategy", where juveniles, subadults, and adults share feeding grounds and similar feeding behaviours [50,51].Differences in diet, jaw size, bite force, and diving ability between juveniles and adults may result in distinct diets and digestive physiologies, potentially reflected in cloacal bacterial communities as increased richness.Studies on loggerhead diet show no significant differences in stomach contents among adults, subadults, and juveniles, but it is worth noting that softer prey like tunicates or jellyfish may not be as detectable due to easier digestion and lack of hard remains in morphology-based studies [50].
In this study, juveniles exhibited higher abundance and more frequent presence of Flavobacteriaceae and Tenacibaculum spp. in oral and cloacal samples than subadult and adult turtles.Tenacibaculum spp.are marine pathogens possibly carried by cnidarians and ctenophores as vectors [52,53], suggesting a higher proportion of soft prey in juveniles' diet.This study was conducted during non-winter months when jellyfish populations increase due to higher water temperatures potentially influencing prey availability.Adults likely have access to these prey types but also consume more challenging prey inaccessible to juveniles.Notably, in other studies, Tenacibaculum spp.were highly abundant on loggerhead skin as well [28], indicating their propensity to inhabit marine vertebrate surfaces and gastrointestinal tracts, irrespective of animal's diet.While in previous research the Mogibacteraceae family detected in faeces correlated with CCL [24], we did not detect it in our data, possibly due to different taxonomy databases being utilized (Greengenes vs. SILVA in our study).BLAST analysis showed ASVs assigned as Peptostreptococcales-Tissieralles closely related to Mogibacterium kristiansenii but it still did not exhibit a CCL correlation.

Oral Bacterial Communities Harbour Distinct Taxa and Could Reflect Recent Diet
In oral samples, the genera Truepera and Halioglobus showed differential abundance.Trueperaceae, a relatively new family, includes one isolate Truepera radiovictrix, known for radiation resistance and thriving in extreme conditions like other Deinococci members [55].Truepera was also a dominant part of the oral microbiota in splendid japalure lizards [56], but its role and functions in reptile oral microbiomes remain unclear.Tolerance to extreme environments and ability to use diverse carbon sources may allow Truepera spp. to outcompete other taxa on reptilian oral mucosa.The genus Halioglobus is typically found in seawater and marine sediments, and has been associated with dinoflagellate blooms and starving, green-lipped mussels [57,58].We have previously reported Halioglobus in oral samples of loggerhead sea turtles [25].Halioglobus bacteria likely derive from loggerhead prey such as mussels or oysters rather than being an intrinsic property of sea turtle mucosal surfaces, given the limited information on host-associated Halioglobus.
Similarly, the high relative abundance of the genus Exiguobacterium in oral and cloacal samples of one juvenile loggerhead (ID069) sampled upon admission may indicate recent feeding, as some Exiguobacterium species are commonly found on shrimp or algae [59][60][61].

Carapace Microbiota Is Distinct, yet Can Harbour Taxa Specific to Oral or Cloacal Microbiota
The carapace and skin of loggerhead sea turtles harbour diverse microbial biofilms, rich in prokaryotes, microeukaryotes, and macroeukaryotes like barnacles and algae [3,28].These surface microbial communities vary with the turtle's geographical location, carapace condition, and the sampled anatomical site [28,29].They can be considered as microbial reservoirs and "diversity hotspots" in otherwise scarce environments [62,63].In our study, we found expectedly distinct microbial communities in the carapace, oral mucosa, and cloaca, although many microbial taxa were present across all body sites, including potential zoonotic pathogens primarily from the cloaca.Oral and cloacal samples shared several differentially abundant microbial taxa (Tenacibaculum, Cardiobacteraceae, Campylobacter), suggesting co-inhabitation of the gastrointestinal tract.The carapace microbiota differed significantly from tank water, with indication of transfer of certain taxa from carapaces to tank water or vice versa.Taxa like Alteromonadaceae and Colwelliaceae (Thalassotalea), differentially abundant on carapaces in this study, were prominent in enclosure tank water in prior studies that did not examine carapace bacterial communities [24], possibly originating from captive turtles' carapaces.The implications of the interplay between microbes from different body sites on loggerhead sea turtles in their natural habitats and during captivity are not yet clear.However, it is essential to consider surface microbiota and potential opportunistic pathogens, especially when rehabilitating severely injured and possibly immunocompromised individuals.

Fungal Communities of Loggerhead Cloaca Are Highly Heterogeneous and Diverse
Fungal communities in loggerhead sea turtles' cloaca and tank water exhibit high variability among individuals, with greater diversity observed in tank water.In this study, we could not attribute differences in the composition and structure of cloacal mycobiota to the turtles' age, sex, or hospitalization status, possibly due to the limited sample size per sampling site and condition assessed.Conversely, captive juvenile green turtles displayed more consistent mycobiota richness and diversity across various sampling sites and health conditions, unaffected by environmental fungi [31].
The taxonomic composition of cloacal and tank water fungal phyla in our study aligns with previously reported marine fungi groups found in green sea turtle faeces [31], as well as marine algicolous fungi, sediments, and sponges [64].Due to the lower resolution of the ITS2 gene marker and many unassigned fungal ASVs beyond the family level, this study offers just a general overview and serves as a foundation for further exploration of specific groups of interest.We sporadically detected pathogenic fungal genera [18], with nine pathogenic genera found in both cloaca and tank water.Among these, Fusarium species (family Nectriaceae) are recognized sea turtle pathogens linked to reduced hatchling success [65].Our study identified Fusarium oxysporum (a known pathogen), Fusarium neocosmosporellium, and Fusarium waltergamsii, not previously reported as sea turtle pathogens but related to known pathogenic species in the Fusarium solani species complex [18,66].Fusarium ASVs were relatively less abundant, except in one sample where they reached 8%.However, some ASVs assigned to the Nectriaceae genus may belong to Fusarium species as ITS region can be insufficient in detecting Fusarium species [67], potentially affecting our reported abundance of pathogenic fungi in cloacal and tank water samples.
The origin of these fungi -whether they are metabolically active and intrinsic to the sampled turtle population or introduced from the environment as spores or through food -is unclear.The impact of these fungi on the host, the interaction with the turtle's immune system and physiology, and their potential host association remain unknown.Many fungal ASVs detected in this study belong to primarily terrestrial taxa, suggesting possible terrestrial sources during rehabilitation or a lack of marine fungal sequences in taxonomy databases.Additionally, numerous reads were assigned only as "Fungi", indicating either undiscovered fungal taxa, lack of representative sequences, or possible turtle host origin.

Low Biomass of Samples Is a Potential Limitation to Interpreting Results
Microbial composition results should be interpreted with caution due to low biomass collected with swabs and fewer fungal cells compared to bacterial cells in vertebrate guts [68].Negative control (sterile water) sequenced with bacterial primers showed a small number of reads that were also detected in low-read samples, indicating potential contamination during DNA extraction and sequencing specifically for samples with low biomass (e.g.kitome) [69,70].Surprisingly, the negative control sequenced with fungal primers had significantly more reads than the bacterial negative control, containing numerous ASVs also found in cloacal and tank water samples.Although the community structures of cloacal and tank water samples differed from the negative control, indicating genuine fungal communities rather than random contaminants, the high number of reads from seemingly low-biomass samples or the "empty" negative control raises concerns about our sample handling and sequencing approach.To our knowledge, no studies addressed fungal contamination in DNA extraction kits as they do for prokaryotes, making it challenging to pinpoint the exact source of this fungal DNA.For future studies, we suggest including additional negative control samples at various sampling and processing stages when investigating undescribed gut-associated fungal communities.

Conclusion
Loggerhead sea turtle-associated microbial communities are crucial for understanding the biology of these endangered reptiles and supporting conservation.In this study, we examined bacterial and fungal communities in juvenile, subadult, and adult loggerhead sea turtles.Our findings revealed distinct microbial communities in the carapace, cloaca, and oral mucosa, with characteristic taxa for each site and shared taxa between cloacal and oral samples (e.g.Tenacibaculum, Moraxellaceae, Cardiobacteriaceae, and Campylobacter).Cloacal bacterial communities exhibited decreasing diversity and changing composition with turtle age, likely due to shifts in diet as juveniles develop stronger bite force and diving capabilities.Microbial exchange with the environment appears to occur, particularly from turtles to tank water, especially from the carapace.Fungal communities in cloaca and tank water displayed high heterogeneity across individuals, with no age or clinical patterns, possibly due to limited samples and low fungal biomass in cloacal samples.Loggerhead sea turtles host complex microbial communities, including potential bacterial and fungal pathogens, which pose risks to handlers or the general public as changing turtle behaviour leads to increased human-turtle interactions.Despite growing research on loggerhead microbiomes, defining a healthy microbiome beyond bacteria remains challenging.Future research should prioritize establishing a description of healthy loggerhead microbiomes at different developmental stages (from eggs and hatchlings to juveniles and adults) to enhance conservation practices and explore potential probiotics and prebiotics for addressing their current and future needs.

Fig. 1
Fig. 1 Locations and body measurements of loggerhead sea turtles with corresponding IDs.(a) Map of locations where loggerhead sea turtles were found prior to transport to rehabilitation centres and sampling.(b) Relationship of loggerheads' weight in kilogrammes and

Fig. 2
Fig. 2 Bacterial community structure and diversity in loggerhead sea turtles' cloacal, oral, and tank water samples.Alpha diversity boxplots with sample density for Faith's phylogenetic diversity (a) and Shannon's index (b) per sample site.(c) Pearson's correlation

Fig. 4
Fig. 4 Structure and composition of fungal microbial communities in cloaca and tank water samples.Alpha diversity boxplots with sample density for Faith's phylogenetic diversity (a) and Shannon's index (b) per sample site.(c) Differential abundance analysis by ANCOM-BC2 for taxa collapsed to species level in cloacal and tank water samples.Cloacal sample site was used as a denominator for log fold change
of doctoral student K. Filek has been fully supported by the "Young Researchers' Career Development Project -Training of Doctoral Students" of the Croatian Science Foundation funded by the European Union from the European Social Fund.

6.37** -4.05**
[54]ers of the Campilobacterota phyla, particularly Arcobacteriaceae, Campylobacteraceae, and Helicobacteraceae, have recently been found to form unique, cold-adapted communities in ectothermic reptiles[54].As expected, and based on their physiology and lifestyle, the Arcobacter genus, which can thrive at atmospheric oxygen levels and prefers lower temperatures, was present in all sample sites.However, it was more abundant in oral than in cloacal samples and was rarely detected in tank water.The Helicobacter and Campylobacter genera, which are vertebrate-associated, were observed in cloacal and occasionally oral samples.Notably, one adult male (ID057) showed a higher prevalence of Helicobacter, reaching up to 20% relative abundance in the cloaca, represented by a single ASV assigned as Helicobacter sp. and not found in any other sample.This suggests a potential overgrowth or infection despite the relatively good clinical status of the individual observed at the time.Therefore, like other reptiles, loggerhead sea turtles could serve as a valuable source of previously undiscovered cold-adapted Campylobacter and Helicobacter species, as well as mucosal Arcobacter species.