A family of short retroposons (Squam1) from squamate reptiles (Reptilia: Squamata): Structure,evolution, and correlation with phylogeny

Sequences of the SINE family specific to squamate reptiles have been isolated from the genomes of lacertid lizards and sequenced. These retroposons, which we called Squam1, are 360–390 bp long and contain a region similar to the tRNA gene sequence at the 5’ end. This family has also been detected in representatives of other reptile families (varanids, iguanids (Anolis), gekkonids, and snakes), being absent from the genomes of crocodiles as well as amphibians, birds, and mammals. The primary structures of Squam1 copies have been comprehensively analyzed and compared with GenBank sequences. The genomes of most taxa contain two to three SINE subfamilies with specific diagnostic features in their primary structures. Individual similarity between the copies within each taxon is about 85%, with intrageneric similarity being only slightly higher. A comparison of consensus sequences between different lizard families has shown that Squam1 may be a convenient phylogenetic marker for this group of reptiles, having a number of both apomorphic and more or less pronounced synapomorphic features. By this criterion, snakes slightly differ from lizards but obviously belong to the same clade. However, they show no special affinity to varanids as the putative closest relatives of snakes, compared to other lizards.     

Development of a Simple and Rapid Method Based on Polymerase Chain Reaction–Based Restriction Fragment Length Polymorphism Analysis to Differentiate Helicobacter, Campylobacter, and Arcobacter Species

Polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) analysis of amplified DNA fragment of the 16S and 23S rRNA genes was performed on 35 Helicobacter, 24 Campylobacter, and 15 Arcobacter strains. PCR amplification generated a 1004-bp fragment of 16S rDNA and a 2.6-Kbp fragment of 23S rDNA from each strain. The amplicons were digested with DdeI and HpaII, respectively. For both assays, distinctive profiles were obtained for each genus. 23S rDNA PCR-RFLP analysis with HpaII enzyme identified Campylobacter and Helicobacter strains at the species level. Analysis of 16S rRNA gene with DdeI enzyme was not useful for the specific identification of Campylobacter and Arcobacter, although it discriminated among Helicobacter species. The PCR-RFLP technique allowed for the discrimination among these three related genus with only one restriction enzyme; therefore it can be a simple, rapid, and useful method for routine identification.     

Detection of Arcobacter spp. in the Coastal Environment of the Mediterranean Sea

The occurrence of Arcobacter spp. was studied in seawater and plankton samples collected from the Straits of Messina, Italy, during an annual period of observation by using cultural and molecular techniques. A PCR assay with three pairs of primers targeting the 16S and 23S rRNA genes was used for detection and identification of Arcobacter butzleri, Arcobacter cryaerophilus, and Arcobacter skirrowii in cultures and environmental samples. Only one of the Arcobacter species, A. butzleri, was isolated from seawater and plankton samples. With some samples the A. butzleri PCR assay gave amplified products when cultures were negative. A. cryaerophilus and A. skirrowii were never detected by culture on selective agar plates; they were detected only by PCR performed directly with environmental samples. Collectively, our data suggest that culturable and nonculturable forms of Arcobacter are present in marine environments. The assay was useful for detecting Arcobacter spp. both as free forms and intimately associated with plankton. This is the first report showing both direct isolation of A. butzleri and the presence of nonculturable Arcobacter spp. in the coastal environment of the Mediterranean Sea.     

Intestinal Microbiota and Species Diversity of Campylobacter and Helicobacter spp. in Migrating Shorebirds in Delaware Bay

Using 16S rRNA gene sequencing analysis, we examined the bacterial diversity and the presence of opportunistic bacterial pathogens (i.e., Campylobacter and Helicobacter) in red knot (Calidris canutus; n = 40), ruddy turnstone (Arenaria interpres; n = 35), and semipalmated sandpiper (Calidris pusilla; n = 22) fecal samples collected during a migratory stopover in Delaware Bay. Additionally, we studied the occurrence of Campylobacter spp., enterococci, and waterfowl fecal source markers using quantitative PCR (qPCR) assays. Of 3,889 16S rRNA clone sequences analyzed, the bacterial community was mostly composed of Bacilli (63.5%), Fusobacteria (12.7%), Epsilonproteobacteria (6.5%), and Clostridia (5.8%). When epsilonproteobacterium-specific 23S rRNA gene clone libraries (i.e., 1,414 sequences) were analyzed, the sequences were identified as Campylobacter (82.3%) or Helicobacter (17.7%) spp. Specifically, 38.4%, 10.1%, and 26.0% of clone sequences were identified as C. lari (>99% sequence identity) in ruddy turnstone, red knot, and semipalmated sandpiper clone libraries, respectively. Other pathogenic species of Campylobacter, such as C. jejuni and C. coli, were not detected in excreta of any of the three bird species. Most Helicobacter-like sequences identified were closely related to H. pametensis (>99% sequence identity) and H. anseris (92% sequence identity). qPCR results showed that the occurrence and abundance of Campylobacter spp. was relatively high compared to those of fecal indicator bacteria, such as Enterococcus spp., E. faecalis, and Catellicoccus marimammalium. Overall, the results provide insights into the complexity of the shorebird gut microbial community and suggest that these migratory birds are important reservoirs of pathogenic Campylobacter species.     

Captive Reptile Mortality Rates in the Home and Implications for the Wildlife Trade

The trade in wildlife and keeping of exotic pets is subject to varying levels of national and international regulation and is a topic often attracting controversy. Reptiles are popular exotic pets and comprise a substantial component of the live animal trade. High mortality of traded animals raises welfare concerns, and also has implications for conservation if collection from the wild is required to meet demand. Mortality of reptiles can occur at any stage of the trade chain from collector to consumer. However, there is limited information on mortality rates of reptiles across trade chains, particularly amongst final consumers in the home. We investigated mortality rates of reptiles amongst consumers using a specialised technique for asking sensitive questions, additive Randomised Response Technique (aRRT), as well as direct questioning (DQ). Overall, 3.6% of snakes, chelonians and lizards died within one year of acquisition. Boas and pythons had the lowest reported mortality rates of 1.9% and chameleons had the highest at 28.2%. More than 97% of snakes, 87% of lizards and 69% of chelonians acquired by respondents over five years were reported to be captive bred and results suggest that mortality rates may be lowest for captive bred individuals. Estimates of mortality from aRRT and DQ did not differ significantly which is in line with our findings that respondents did not find questions about reptile mortality to be sensitive. This research suggests that captive reptile mortality in the home is rather low, and identifies those taxa where further effort could be made to reduce mortality rates.     

Prevalence and molecular identification of Cryptosporidium isolates from pet lizards and snakes in Italy

In order to acquire prevalence and genetic data on Cryptosporidium infections in captive lizards and snakes kept as pets, a survey was conducted on 150 individual reptiles from southern Italy. Fecal samples were preserved in 5% formalin and analyzed using a commercial immunofluorescence assay (IFA) for the detection of Cryptosporidium oocysts and Giardia cysts. IFA revealed the presence of Cryptosporidium oocysts in nine of the 150 samples examined (6.0%), precisely in 6/125 snakes (4.8%) and in 3/25 lizards (12.0%); all fecal samples tested negative for the presence of Giardia cysts. Molecular characterization based on nested PCR amplification and sequencing of the SSU-rRNA gene, revealed the presence of Cryptosporidium serpentis in three samples from snakes (Boa constrictor constrictor, Elapheguttata guttata guttata and Python molurus).     

A Molecular Enrichment Strategy Based on <Emphasis Type="Italic">cpn</Emphasis>60 for Detection of Epsilon-Proteobacteria in the Dog Fecal Microbiome

Members of the rare microbiome can be important components of complex microbial communities. For example, pet dog ownership is a known risk factor for human campylobacteriosis, and Campylobacter is commonly detected in dog feces by targeted assays. However, these organisms have not been detected by metagenomic methods. The goal of this study was to characterize fecal microbiota from healthy and diarrheic pet dogs using two different levels of molecular detection. PCR amplification and pyrosequencing of the universal cpn60 gene target was used to obtain microbial profiles from each dog. To investigate the relatively rare epsilon-proteobacteria component of the microbiome, a molecular enrichment was carried out using a PCR that first amplified the cpn10–cpn60 region from epsilon-proteobacteria, followed by universal cpn60 target amplification and pyrosequencing. From the non-enriched survey, the major finding was a significantly higher proportion of Bacteroidetes, notably Bacteroides vulgatus, in healthy dogs compared to diarrheic dogs. Epsilon-proteobacteria from the genera Helicobacter and Campylobacter were also detected at a low level in the non-enriched profiles of some dogs. Molecular enrichment increased the proportion of epsilon-proteobacteria sequences detected from each dog, as well as identified novel, presumably rare sequences not seen in the non-enriched profiles. Enriched profiles contained known species of Arcobacter, Campylobacter, Flexispira, and Helicobacter and identified two possibly novel species. These findings add to our understanding of the canine fecal microbiome in general, the epsilon-proteobacteria component specifically, and present a novel modification to traditional metagenomic approaches for study of the rare microbiome.     

Arcobacter molluscorum sp. nov., a new species isolated from shellfish

Nineteen bacteria isolates recovered from shellfish samples (mussels and oysters) showed a new and specific 16S rDNA-RFLP pattern with an Arcobacter identification method designed to recognize all species described up to 2008. These results suggested that they could belong to a new species. ERIC-PCR revealed that the 19 isolates belonged to 3 different strains. The sequence of the 16S rRNA gene of a representative strain (F98-3^T) showed 97.6% similarity with the closest species Arcobacter marinus followed by Arcobacter halophilus (95.6%) and Arcobacter mytili (94.7%). The phylogenetic analysis with the16S rRNA, rpoB, gyrB and hsp60 genes placed the shellfish strains within the same cluster as the three species mentioned (also isolated from saline habitats) but they formed an independent phylogenetic line. The DDH results between strain F98-3^T and A. marinus (54.8% ± 1.05), confirmed that it represents a new species. Several biochemical tests differentiated the shellfish isolates from all other Arcobacter species. Although the new species was different from A. mytili, they shared not only the same habitat (mussels) but also the characteristic of being so far the only Arcobacter species that are simultaneously negative for urea and indoxyl acetate hydrolysis. All results supported the classification of the shellfish strains as a new species, for which the name Arcobacter molluscorum sp. nov. with the type strain F98-3^T is proposed (=CECT 7696^T = LMG 25693^T).     

The Microaerophilic Microbiota of De-Novo Paediatric Inflammatory Bowel Disease: The BISCUIT Study

Introduction

Children presenting for the first time with inflammatory bowel disease (IBD) offer a unique opportunity to study aetiological agents before the confounders of treatment. Microaerophilic bacteria can exploit the ecological niche of the intestinal epithelium; Helicobacter and Campylobacter are previously implicated in IBD pathogenesis. We set out to study these and other microaerophilic bacteria in de-novo paediatric IBD.

Patients and Methods

100 children undergoing colonoscopy were recruited including 44 treatment naïve de-novo IBD patients and 42 with normal colons. Colonic biopsies were subjected to microaerophilic culture with Gram-negative isolates then identified by sequencing. Biopsies were also PCR screened for the specific microaerophilic bacterial groups: Helicobacteraceae, Campylobacteraceae and Sutterella wadsworthensis.

Results

129 Gram-negative microaerophilic bacterial isolates were identified from 10 genera. The most frequently cultured was S. wadsworthensis (32 distinct isolates). Unusual Campylobacter were isolated from 8 subjects (including 3 C. concisus, 1 C. curvus, 1 C. lari, 1 C. rectus, 3 C. showae). No Helicobacter were cultured. When comparing IBD vs. normal colon control by PCR the prevalence figures were not significantly different (Helicobacter 11% vs. 12%, p = 1.00; Campylobacter 75% vs. 76%, p = 1.00; S. wadsworthensis 82% vs. 71%, p = 0.312).

Conclusions

This study offers a comprehensive overview of the microaerophilic microbiota of the paediatric colon including at IBD onset. Campylobacter appear to be surprisingly common, are not more strongly associated with IBD and can be isolated from around 8% of paediatric colonic biopsies. S. wadsworthensis appears to be a common commensal. Helicobacter species are relatively rare in the paediatric colon.

Trial Registration

This study is publically registered on the United Kingdom Clinical Research Network Portfolio (9633).     

Presence and Analysis of Plasmids in Human and Animal Associated Arcobacter Species

In this study, we report the screening of four Arcobacter species for the presence of small and large plasmids. Plasmids were present in 9.9% of the 273 examined strains. One Arcobacter cryaerophilus and four Arcobacter butzleri plasmids were selected for further sequencing. The size of three small plasmids isolated from A. butzleri and the one from A. cryaerophilus strains ranged between 4.8 and 5.1 kb, and the size of the large plasmid, isolated from A. butzleri, was 27.4 kbp. The G+C content of all plasmids ranged between 25.4% and 26.2%. A total of 95% of the large plasmid sequence represents coding information, which contrasts to the 20 to 30% for the small plasmids. Some of the open reading frames showed a high homology to putative conserved domains found in other related organisms, such as replication, mobilization and genes involved in type IV secretion system. The large plasmid carried 35 coding sequences, including seven genes in a contiguous region of 11.6 kbp that encodes an orthologous type IV secretion system found in the Wolinella succinogenes genome, Helicobacter pylori and Campylobacter jejuni plasmids, which makes this plasmid interesting for further exploration.     

Facultative parthenogenesis discovered in wild vertebrates

Facultative parthenogenesis (FP)—asexual reproduction by bisexual species—has been documented in a variety of multi-cellular organisms but only recently in snakes, varanid lizards, birds and sharks. Unlike the approximately 80 taxa of unisexual reptiles, amphibians and fishes that exist in nature, FP has yet to be documented in the wild. Based on captive documentation, it appears that FP is widespread in squamate reptiles (snakes, lizards and amphisbaenians), and its occurrence in nature seems inevitable, yet the task of detecting FP in wild individuals has been deemed formidable. Here we show, using microsatellite DNA genotyping and litter characteristics, the first cases of FP in wild-collected pregnant females and their offspring of two closely related species of North American pitviper snakes—the copperhead (Agkistrodon contortrix) and cottonmouth (Agkistrodon piscivorus). Our findings support the view that non-hybrid origins of parthenogenesis, such as FP, are more common in squamates than previously thought. With this confirmation, FP can no longer be viewed as a rare curiosity outside the mainstream of vertebrate evolution. Future research on FP in squamate reptiles related to proximate control of induction, reproductive competence of parthenogens and population genetics modelling is warranted.     

Morphology,Reproduction and Diet in Australian and Papuan Death Adders (Acanthophis,Elapidae)

Death adders (genus Acanthophis) differ from most other elapid snakes, and resemble many viperid snakes, in their thickset morphology and ambush foraging mode. Although these snakes are widely distributed through Australia and Papua New Guinea, their basic biology remains poorly known. We report morphological and ecological data based upon dissection of >750 museum specimens drawn from most of the range of the genus. Female death adders grow larger than conspecific males, to about the same extent in all taxa (20% in mean adult snout-vent length,  =  SVL). Most museum specimens were adult rather than juvenile animals, and adult males outnumbered females in all taxa except A. pyrrhus. Females have shorter tails (relative to SVL) than males, and longer narrower heads (relative to head length) in some but not all species. The southern A. antarcticus is wider-bodied (relative to SVL) than the other Australian species. Fecundity of these viviparous snakes was similar among taxa (mean litter sizes 8 to 14). Death adders encompass a broad range of ecological attributes, taking a wide variety of vertebrate prey, mostly lizards (55%), frogs and mammals (each 21%; based on 217 records). Dietary composition differed among species (e.g. frogs were more common in tropical than temperate-zone species), and shifted with snake body size (endotherms were taken by larger snakes) and sex (male death adders took more lizards than did females). Overall, death adders take a broader array of prey types, including active fast-moving taxa such as endotherms and large diurnal skinks, than do most other Australian elapids of similar body sizes. Ambush foraging is the key to capturing such elusive prey.     

Antimicrobial resistance and virulence genes profiles of Arcobacter butzleri strains isolated from back yard chickens and retail poultry meat in Chile

This research aims to investigate the presence and pathogenic potential of Arcobacter in poultry meat samples purchased in the retail market of Valdivia (South of Chile) as well as in faecal samples from backyard chickens from rural areas around this city. The isolates obtained were identified by molecular methods. Furthermore, putative virulence genes were assessed by PCR and the antimicrobial resistance was tested by phenotypic methods. Arcobacter was present in 41·6% of the samples, with the highest value in retail poultry meat (55·7%) followed by backyard production (28·0%). Arcobacter butzleri was the most prevalent species (75·6%) followed by Arcobacter skirrowii (14·8%) and Arcobacter cryaerophilus (9·6%). An 8·5% of A. butzleri strains from meat were resistant to both ciprofloxacin and tetracycline and 6·1% were resistant to erythromycin, while none was resistant to gentamycin, unlike strains from domestic chickens, which showed no resistance. Furthermore, A. butzleri strains from chicken meat presented a higher prevalence of virulence genes than strains from domestic chickens. In fact, in this last group, some genes (hecA, hecB and irgA) were completely absent. Therefore, this study provides insight on the epidemiology of Arcobacter in Chilean poultry and suggests that under traditional breeding conditions strains are, apparently, less pathogenic and drug resistant.     

Adherence to and Invasion of Human Intestinal Cells by Arcobacter Species and Their Virulence Genotypes

The genus Arcobacter is composed of 17 species which have been isolated from various sources. Of particular interest are A. butzleri, A. cryaerophilus, and A. skirrowii, as these have been associated with human cases of diarrhea, the probable transmission routes being through the ingestion of contaminated drinking water and food. To date, only limited studies of virulence traits in this genus have been undertaken. The present study used 60 Arcobacter strains isolated from different sources, representing 16 of the 17 species of the genus, to investigate their ability to adhere to and invade the human intestinal cell line Caco-2. In addition, the presence of five putative virulence genes (ciaB, cadF, cj1349, hecA, and irgA) was screened for in these strains by PCR. All Arcobacter species except A. bivalviorum and Arcobacter sp. strain W63 adhered to Caco-2 cells, and most species (10/16) were invasive. The most invasive species were A. skirrowii, A. cryaerophilus, A. butzleri, and A. defluvii. All invasive strains were positive for ciaB (encoding a putative invasion protein). Other putative virulence genes were present in other species, i.e., A. butzleri (cadF, cj1349, irgA, and hecA), A. trophiarum (cj1349), A. ellisii (cj1349), and A. defluvii (irgA). No virulence genes were detected in strains which showed little or no invasion of Caco-2 cells. These results indicate that many Arcobacter species are potential pathogens of humans and animals.     

Genetic Diversity of Cryptosporidium spp. in Captive Reptiles

The genetic diversity of Cryptosporidium in reptiles was analyzed by PCR-restriction fragment length polymorphism and sequence analysis of the small subunit rRNA gene. A total of 123 samples were analyzed, of which 48 snake samples, 24 lizard samples, and 3 tortoise samples were positive for Cryptosporidium. Nine different types of Cryptosporidium were found, including Cryptosporidium serpentis, Cryptosporidium desert monitor genotype, Cryptosporidium muris, Cryptosporidium parvum bovine and mouse genotypes, one C. serpentis-like parasite in a lizard, two new Cryptosporidium spp. in snakes, and one new Cryptosporidium sp. in tortoises. C. serpentis and the desert monitor genotype were the most common parasites and were found in both snakes and lizards, whereas the C. muris and C. parvum parasites detected were probably the result of ingestion of infected rodents. Sequence and biologic characterizations indicated that the desert monitor genotype was Cryptosporidium saurophilum. Two host-adapted C. serpentis genotypes were found in snakes and lizards.     

Partitioning of Bacterial Communities between Seawater and Healthy, Black Band Diseased, and Dead Coral Surfaces

Distinct partitioning has been observed in the composition and diversity of bacterial communities inhabiting the surface and overlying seawater of three coral species infected with black band disease (BBD) on the southern Caribbean island of Curaçao, Netherlands Antilles. PCR amplification and sequencing of bacterial 16S rRNA genes (rDNA) with universally conserved primers have identified over 524 unique bacterial sequences affiliated with 12 bacterial divisions. The molecular sequences exhibited less than 5% similarity in bacterial community composition between seawater and the healthy, black band diseased, and dead coral surfaces. The BBD bacterial mat rapidly migrates across and kills the coral tissue. Clone libraries constructed from the BBD mat were comprised of eight bacterial divisions and 13% unknowns. Several sequences representing bacteria previously found in other marine and terrestrial organisms (including humans) were isolated from the infected coral surfaces, including Clostridium spp., Arcobacter spp., Campylobacter spp., Cytophaga fermentans, Cytophaga columnaris, and Trichodesmium tenue.     

Arcobacter cloacae sp. nov. and Arcobacter suis sp. nov., two new species isolated from food and sewage

Three strains recovered from mussels (F26), sewage (SW28-13^T) and pork meat (F41^T) were characterized as Arcobacter. They did not appear to resemble any known species on the basis of their 16S rDNA-RFLP patterns and the rpoB gene analyses. However, strains F26 and SW28-13^T appeared to be the same species. The 16S rRNA gene sequence similarity of strains SW28-13^T and F41^T to the type strains of all other Arcobacter species ranged from 94.1% to 99.6% and 93.4% to 98.8%, respectively. Phenotypic characteristics and the DNA–DNA hybridization (DDH) results showed that they belonged to 2 new Arcobacter species. A multilocus phylogenetic analysis (MLPA) with the concatenated sequences of 5 housekeeping genes (gyrA, atpA, rpoB, gyrB and hsp60) was used for the first time in the genus, showing concordance with the 16S rRNA gene phylogenetic analysis and DDH results. The MALDI-TOF mass spectra also discriminated these strains as two new species. The names proposed for them are Arcobacter cloacae with the type strain SW28-13^T (=CECT 7834^T = LMG 26153^T) and Arcobacter suis with the type strain F41^T (=CECT 7833^T = LMG 26152^T).     

<Emphasis Type="Italic">Arcobacter acticola</Emphasis> sp. nov., isolated from seawater on the East Sea in South Korea

A Gram-stain-negative, facultative aerobic, non-flagellated, and rod-shaped bacterium, designated AR-13^T, was isolated from a seawater on the East Sea in South Korea, and subjected to a polyphasic taxonomic study. Strain AR-13^T grew optimally at 30°C, at pH 7.0–8.0 and in the presence of 0–0.5% (w/v) NaCl. The phylogenetic trees based on 16S rRNA gene sequences showed that strain AR-13^T fell within the clade comprising the type strains of Arcobacter species, clustering coherently with the type strain of Arcobacter venerupis. Strain AR-13^T exhibited 16S rRNA gene sequence similarity values of 98.1% to the type strain of A. venerupis and of 93.2–96.9% to the type strains of the other Arcobacter species. Strain AR-13^T contained MK-6 as the only menaquinone and summed feature 3 (C_16:1ω7c and/or C_16:1ω6c), C_16:0, C_18:1ω7c, and summed feature 2 (iso-C_16:1 I and/or C_14:0 3-OH) as the major fatty acids. The major polar lipids detected in strain AR-13^T were phosphatidylethanolamine, phosphatidylglycerol, and one unidentified aminophospholipid. The DNA G+C content was 28.3 mol% and its mean DNA-DNA relatedness value with the type strain of A. venerupis was 21%. Differential phenotypic properties, together with its phylogenetic and genetic distinctiveness, revealed that strain AR-13^T is separated from recognized Arcobacter species. On the basis of the data presented, strain AR-13^T is considered to represent a novel species of the genus Arcobacter, for which the name Arcobacter acticola sp. nov. is proposed. The type strain is AR-13^T (=KCTC 52212^T =NBRC 112272^T).     

Snakes on the Balearic Islands: An Invasion Tale with Implications for Native Biodiversity Conservation

Biological invasions are a major conservation threat for biodiversity worldwide. Islands are particularly vulnerable to invasive species, especially Mediterranean islands which have suffered human pressure since ancient times. In the Balearic archipelago, reptiles represent an outstanding case with more alien than native species. Moreover, in the last decade a new wave of alien snakes landed in the main islands of the archipelago, some of which were originally snake-free. The identification of the origin and colonization pathways of alien species, as well as the prediction of their expansion, is crucial to develop effective conservation strategies. In this study, we used molecular markers to assess the allochthonous status and the putative origin of the four introduced snake species (Hemorrhois hippocrepis, Malpolon monspessulanus, Macroprotodon mauritanicus and Rhinechis scalaris) as well as ecological niche models to infer their patterns of invasion and expansion based on current and future habitat suitability. For most species, DNA sequence data suggested the Iberian Peninsula as the potential origin of the allochthonous populations, although the shallow phylogeographic structure of these species prevented the identification of a restricted source-area. For all of them, the ecological niche models showed a current low habitat suitability in the Balearic, which is however predicted to increase significantly in the next few decades under climate change scenarios. Evidence from direct observations and spatial distribution of the first-occurrence records of alien snakes (but also lizards and worm lizards) suggest the nursery trade, and in particular olive tree importation from Iberian Peninsula, as the main pathway of introduction of alien reptiles in the Balearic islands. This trend has been reported also for recent invasions in NE Spain, thus showing that olive trees transplantation may be an effective vector for bioinvasion across the Mediterranean. The combination of molecular and ecological tools used in this study reveals a promising approach for the understanding of the complex invasion process, hence guiding conservation management actions.