Introduction
Deception Island (62°57'S, 60°38'W), located in the South Shetland Islands archipelago, is one of few active volcanic locations in Antarctica today (Fig. 1). Due to its unique geothermally associated flora, exceptional aesthetic and scientific value and historical importance, compounded by exposure to multiple contemporary sources of human impact from both national research operations and the tourism industry, the entire island has been designated as Antarctic Specially Managed Area (ASMA) 4, and it also contains two Antarctic Specially Protected Areas (ASPAs 140 and 145). The area also includes two Historic Sites or Monuments (HSM 71 and 76), an important reminder that human activity in Antarctic marine exploitation dates back to the early 1820s (Held & Blanchette Reference Held and Blanchette2017).
The sealing industry was very active in the South Shetland Islands during the nineteenth century after the archipelago's discovery in the early 1820s, followed by the whaling industry in the early twentieth century. The first shore-based whaling station on the island was established in 1911 - the Aktieselskabet Hecktor Whaling Station at Whalers Bay - which operated until 1931 when it was no longer commercially viable. In 1944, as part of a military operation towards the end of the Second World War, the UK constructed its Base B. In the post-war years, now operated by the Falkland Islands Dependencies Survey, that base was expanded to include a small aerodrome and associated hangar. Base B included some buildings from the preceding whaling station but also expanded its footprint further with new buildings. Active aircraft operations took place from 1955 to 1957 and again from 1959 to 1969. Base B was decommissioned in 1969, at least in part because many of its structures were heavily damaged during volcanic eruptions in 1967 and 1969 (Smith Reference Smith1984, Morales et al. Reference Morales, Almendros and Carmona2017). The area was partially cleared by the UK in 1990–1992, but many structures dating from both the whaling and research eras remain, with many of these being constructed from or incorporating wood imported to the island (Fig. 2). The Whalers Bay site has been under the protection of the Antarctic Treaty since being designated as HSM 71 in 1995.
Whalers Bay is one of the most visited and popular sites in Antarctica by both national operator staff and the Antarctic tourism industry. Two national operators (Argentina and Spain) currently operate summer-only research stations on the island, with Chile and the UK previously operating year-round stations before the late 1960s' eruptions. Research and logistical support vessels from multiple national operators also routinely visit the island and land personnel in most summers. A total of ~160 000 tourists visited Whalers Bay between 2010 and 2019 (www.iaato.org; Carvalho-Silva et al. Reference Carvalho-Silva, Rosa, Pinto, Da Silva, Henriques, Convey and Câmara2021), the majority of whom visited and explored the ruins of HSM 71. Human impacts on Whalers Bay local terrestrial ecosystems, in the context of the presence of non-native species anthropogenically transferred to the island (Greenslade et al. Reference Greenslade, Potapov, Russell and Convey2012, Hughes et al. Reference Hughes, Pertierra, Molina-Montenegro and Convey2015) or of environmental DNA (eDNA) sequences assigned to exotic species (Rosa et al. Reference Rosa, da Silva, Ogaki, Pinto, Stech and Convey2020, Câmara et al. Reference Câmara, Carvalho-Silva, Pinto, Amorin, Henriques and da Silva2021b, Carvalho-Silva et al. Reference Carvalho-Silva, Rosa, Pinto, Da Silva, Henriques, Convey and Câmara2021), have received some research attention. Many of the wooden structures in Whalers Bay are > 100 years old and are increasingly deteriorating. The effects of decomposer fungal communities on the wood of many Antarctic and Arctic historical sites have been investigated, including on Deception Island (Held et al. Reference Held, Arenz, Blanchette, Barr and Chaplin2011, Held & Blanchette Reference Held and Blanchette2017, Blanchette et al. Reference Blanchette, Held, Jurgens, Stear and Dupont2021). However, little to nothing is known regarding the impacts of other organism groups that may have become established on these exotic wood habitats, such as microalgae and protozoans, as well as regarding the presence of pollen and spores.
Many organisms of these types are difficult to identify conclusively, requiring the development of appropriate culture methodologies, and some may be encysted. Recently developed molecular tools such as DNA metabarcoding using high-throughput sequencing (HTS; Taberlet et al. Reference Taberlet, Coissac, Pompanon, Brochmann and Willerslev2012) are increasingly being applied to infer the presence of active or dormant life forms, propagules, pollen and detritus of plants (Fahner et al. Reference Fahner, Shokralla, Baird and Hajibabaei2016). These methodologies use the total and even sometimes the degraded DNA extracted from environmental samples (e.g. water, soil or air). The methodology has the potential to increase by a factor of ~11 the number of taxa detected when compared with classical morphological approaches (Rippin et al. Reference Rippin, Borchhardt, Williams, Colesie, Jung and Büdel2018), and they can reveal the presence of DNA of taxa not detectable through traditional surveys. In the current study, we used HTS of eDNA to investigate the diversity of non-fungal eukaryotic organisms present in wood samples obtained from various sites in HSM 71 Whalers Bay on Deception Island.
Methods
Study sites and sampling
Four sites were selected for sampling: Biscoe House, an abandoned whaling boat, the Biscoe House Annex and the Hunting Lodge (Fig. 2), all located in close proximity to each other. Three small samples (~1 cm3) of wood from each site were obtained during the Antarctic summer (December 2016) and placed into individual sterilized Whirl-Pak bags (Sigma-Aldrich, USA), which were sealed and kept at -20°C for 3 weeks until being processed at the Microbiology Laboratory at the Federal University of Minas Gerais, Brazil. Data from each sample were obtained separately but are presented together in Table I.
DNA extraction, amplification and sequencing
Total DNA was extracted using a modified sodium dodecyl sulphate (SDS) extraction method (Goldenberger et al. Reference Goldenberger, Perschil, Ritzler and Alwegg1995, Zhou et al. Reference Zhou, Bruns and Tiedje1996, Natarajan et al. Reference Natarajan, Zhang, Morono, Inagaki and Wang2016). Wood fragments of ~1 cm3 were added to plastic tubes each containing 2 ml of SDS extraction buffer (0.1 M ethylenediaminetetraacetic acid (EDTA) at pH 8 and 2% SDS) and ground with sterilized iron beads for 3 min before being incubated at 55°C for 16 h. Then, 330 μl of 5 M NaCl and 330 μl of pre-heated cetrimonium bromide (CTAB) 10% (55°C) were added; the solution was then vortexed, spun down and incubated at 55°C for 10 min. The solution was then transferred to a new tube, and 600 μl of chloroform was added and vortexed at maximum speed for 1 min. After that, the tubes were centrifuged at 13 000 rpm for 10 min and the supernatant transferred to a new tube.
The extracted DNA was cleaned using the Genomic DNA Purification Kit (QIAGEN, USA) following the manufacturer's instructions. Extractions were carried out under strict sterile conditions to avoid contamination (de Menezes et al. Reference de Menezes, Porto, Amorim, Zani, de Almeida Alves and Sales Junior2020). Samples were manipulated inside a sterile flow hood, and all equipment (forceps, spatula, etc.) was previously disinfected. ddH2O samples were used as blanks to ensure the absence of contamination.
DNA quality was analysed using agarose gel electrophoresis (1% agarose in 1 × Trisborate-EDTA) and then quantified using the Quanti- iT™ Pico Green dsDNA Assay (Invitrogen, USA). The internal transcribed spacer 2 region (ITS2) of the nuclear ribosomal DNA was used as a DNA barcode for molecular species identification (Chen et al. Reference Chen, Yao, Han, Liu, Song and Shi2010, Richardson et al. Reference Richardson, Lin, Sponsler, Quijia, Goodell and Johnson2015, Hadi et al. Reference Hadi, Santana, Brunale, Gomes, Oliveira and Matthiensen2016, Ruppert et al. Reference Ruppert, Kline and Rahman2019, Câmara et al. Reference Câmara, de Souza, Pinto, Convey, Amorim, Carvalho-Silva and Rosa2021a,Reference Câmara, Carvalho-Silva, Pinto, Amorin, Henriques and da Silvab, Reference Câmara, de Menezes, Oliveira, Souza, Amorim and Shaefer2022, Carvalho-Silva et al. Reference Carvalho-Silva, Rosa, Pinto, Da Silva, Henriques, Convey and Câmara2021) using the universal primers ITS3 and ITS4 (White et al. Reference White, Bruns, Lee, Taylor, Innis, Gelfand, Sninsky and White1990). Library construction and DNA amplification were performed using the Herculase II Fusion DNA Polymerase Nextera XT Index Kit V2, following the Illumina 16S Metagenomic Sequencing Library Preparation Part #15044223 Rev. B protocol. Paired-end sequencing (2 × 300 bp) was performed commercially on a MiSeq System (Illumina, USA) by Macrogen, Inc. (South Korea), including negative controls.
Data analyses and taxa identification
Quality analysis was carried out using BBDuk v. 38.87 in BBmap software (Bushnell Reference Bushnell2014) with the following parameters: Illumina adapters removing (Illumina artefacts and the PhiX Control v3 Library); ktrim = l; k = 23; mink = 11; hdist = 1; minlen = 50; tpe; tbo; qtrim = rl; trimq = 20; ftm = 5; maq = 20. The remaining sequences were imported into QIIME2 version 2021.4 (https://qiime2.org/) for bioinformatics analyses (Bolyen et al. Reference Bolyen, Rideout, Dillon, Bokulich, Abnet, Al-Ghalith and Alexander2019). The qiime2-dada2 plugin was used for filtering, dereplication, turning paired-end fastq files into merged files, removing chimeras and creating amplicon sequence variants (ASVs) with default parameters (Callahan et al. Reference Callahan, McMurdie, Rosen, Han, Johnson and Holmes2016). Taxonomic assignments of ASVs were determined using the qiime2-feature-classifier (Bokulich et al. Reference Bokulich, Kaehler, Rideout, Dillon, Bolyen and Knight2018) classify-sklearn against different databases; the sequence similarity threshold was set at 97%. Firstly, ASVs were classified against the PLANiTS2 database (Banchi et al. Reference Banchi, Ametrano, Greco, Stanković, Muggia and Pallavicini2020). After this step, ASVs that were not classified were filtered and classify-sklearn classified against the UNITE Eukaryotes ITS database version 8.3 (Abarenkov et al. Reference Abarenkov, Zirk, Piirmann, Pöhönen, Ivanov, Nilsson and Kõljalg2020). Finally, any remaining unclassified ASVs were filtered and aligned against the filtered National Center for Biotechnology Information (NCBI) non-redundant nucleotide sequences (nt) database (October 2021) using BLASTn (Camacho et al. Reference Camacho, Coulouris, Avagyan, Ma, Papadopoulos, Bealer and Madden2009) with default parameters; the nt database was filtered with the following keywords: ‘ITS1’, ‘ITS2’, ‘internal transcribed spacer’ and ‘internal transcribed spacer’. Taxonomic assignments were performed using MEGAN6 (Huson et al. Reference Huson, Beier, Flade, Górska, El-Hadidi and Mitra2016).
For comparative purposes, the number of reads can be used as a proxy for relative abundance (Giner et al. Reference Giner, Forn, Romac, Logares, de Vargas and Massana2016, Deiner et al. Reference Deiner, Bik, Mächler, Seymour, Lacoursière-Roussel and Altermatt2017, Câmara et al. Reference Câmara, de Souza, Pinto, Convey, Amorim, Carvalho-Silva and Rosa2021a,Reference Câmara, Carvalho-Silva, Pinto, Amorin, Henriques and da Silvab, Reference Câmara, de Menezes, Oliveira, Souza, Amorim and Shaefer2022). Classifications and systematic ranks for kingdoms and phyla followed Ruggiero et al. (Reference Ruggiero, Gordon, Orrell, Bailly, Bourgoin and Brusca2015). Venn diagrams were prepared as described by Bardou et al. (Reference Bardou, Mariette, Escudié, Djemiel and Klopp2014).
Results
A total of 2 327 468 reads were obtained, of which 2 218 794 remained after quality control (for detailed data, see Supplemental Table 1). A total of 19 820 reads were DNA sequences representing three kingdoms (Chromista, Protozoa and Viridiplantae) and four phyla (Ciliophora, Perclozoa, Chlorophyta and Magnoliophyta). The remaining reads corresponded to fungi (de Souza et al. Reference De Souza, Teixeira, Coelho, Lopes, Convey and Carvalho-Silva2022). All calculated rarefaction curves (Supplemental Fig. 1) reached a plateau, indicating that the sampling effort was sufficient to represent the taxa analysed in all sampled sites. Sequences representing a total of 43 taxa were detected (Table I). Biscoe House Annex hosted the greatest diversity, with 20 taxa, followed by the whaling boat, Biscoe House and the Hunting Lodge, with 16, 15 and 12 taxa, respectively. Representatives of the ciliate order Sporadotrichida were the most frequently detected sequences, this being the only ciliate taxa present at all sites, followed by the green algal genus Trebouxia, which was present in three of the four sites. The Venn diagram illustrates that only two taxa were detected at all sites (Sporadotrichida and Chlamydomonadales; Fig. 3).
Non-metric multidimensional scaling and cluster analysis did not identify any differences between sampled structures, which were in close proximity, of similar material and subject to the same environmental conditions. Data are available upon request.
Discussion
Among the Chromista detected, although giving the greatest taxonomic resolution of sequence assignments, none were assigned to species level. Fuscheria is a genus of approximately six species and is globally distributed, inhabiting terrestrial, fresh, brackish and marine environments (Petz et al. Reference Petz, Song and Wilbert1995). Two terrestrial species - Fuscheria lacustres and Fuscheria terricola - have been reported from the Antarctic continent near Casey Station (ATCM 2013, Thompson et al. Reference Thompson, Powell and Adams2019), and the Maritime Antarctic species Fuscheria marina has been reported from the Weddell Sea (Petz et al. Reference Petz, Song and Wilbert1995). The genus Homalogastra includes one species record (Homalogastra setosa) from Antarctica (Thompson et al. Reference Thompson, Powell and Adams2019). Trithigmostoma is a cosmopolitan genus with five species inhabiting fresh, marine and brackish water. It has been reported from Antarctica from soils in Wilkes Land (Petz & Foissner Reference Petz and Foissner1997). There are no previous records of these genera from Deception Island.
All other taxa were assigned to higher taxonomic levels. Some of these include taxa with wide geographical distributions, as also reported by Câmara et al. (Reference Câmara, de Souza, Pinto, Convey, Amorim, Carvalho-Silva and Rosa2021a). The very abundant order Sporadotrichida is a highly diverse group present in marine, freshwater and terrestrial habitats, containing four families and 83 genera, some of which are very common, such as Halteria, Oxytricha and Cyrtohymena (Parr et al. Reference Parr, Wilson, Leary, Schulz, Lans and Walley2014), which are all recorded from Antarctica (Thompson et al. Reference Thompson, Powell and Adams2019, Câmara et al. Reference Câmara, de Souza, Pinto, Convey, Amorim, Carvalho-Silva and Rosa2021a,Reference Câmara, Carvalho-Silva, Pinto, Amorin, Henriques and da Silvab,Reference Câmara, Convey, Rangel, Konrath, Barreto and Pintoc). As sequences were assigned only at the order level, and considering that this order contains many very common species, this may be suggestive of the presence of new and as yet undescribed species present on the wood sampled. Ciliates are known to feed on bacteria, which are likely to be present in decaying wood. Protists are also known to play an important role in cellulose degradation (Peterson et al. Reference Peterson, Stewart and Scharf2015). However, functional roles cannot be inferred directly in the current study. Ciliates in the class Spirotrichea and family Acropisthiidae are known to occur in both marine and terrestrial environments in Antarctica (De Broyer et al. Reference De Broyer, Clarke, Koubbi, Pakhomov, Scott, Vanden Berghe and Danis2020). The ciliate groups assigned here are widespread and have been previously recorded in Antarctica, and our data confirm their presence in decaying wood of anthropogenic origin.
Among the Protozoa, the genus Allovahlkampfia includes free-living amoebae that feed primarily on bacteria. It has previously been recorded from Deception and King George islands (Câmara et al. Reference Câmara, Carvalho-Silva, Pinto, Amorin, Henriques and da Silva2021b), and representatives are known to be capable of causing disease in humans (Mohamed et al. Reference Mohamed, Huseein, Farrag, Mostafa and Hassan2016, Tolba et al. Reference Tolba, Huseein, Farrag, Mohamed, Kobayashi and Suzuki2016). In contrast, Neovahlkampfia is a genus that includes only one described marine species that has not been previously recorded from Antarctica. The taxonomy of these genera remains poorly known (Jonckheere et al. Reference Jonckheere, Murase and Opperdoes2011). However, as with the ciliates, bacterivorous organisms are expected to be found on decaying wood.
The assigned Viridiplantae sequences include some very common taxa (genera Trebouxia, Chlamydomonas and Prasiola) previously recorded from multiple locations in Antarctica, including Deception Island (e.g. Câmara et al. Reference Câmara, de Souza, Pinto, Convey, Amorim, Carvalho-Silva and Rosa2021a,Reference Câmara, Carvalho-Silva, Pinto, Amorin, Henriques and da Silvab, Fonseca et al. Reference Fonseca, Câmara, Ogaki, Pinto, Lirio and Coria2022). The close proximity of the Whalers Bay sampling locations to the sea as well as to snow and ice could underlie the taxon composition found. Representatives of Trebouxia are very common and widespread in Antarctica (Câmara et al. Reference Câmara, de Souza, Pinto, Convey, Amorim, Carvalho-Silva and Rosa2021a,Reference Câmara, Carvalho-Silva, Pinto, Amorin, Henriques and da Silvab), and some members are also photobionts in lichens found growing on the wood surface as well as on other natural substrata (de Souza et al. Reference De Souza, Teixeira, Coelho, Lopes, Convey and Carvalho-Silva2022).
In contrast with the algae, the flowering plant (Magnoliophyta) sequence assignments all represented exotic taxa not native to Antarctica, and the presence of their DNA is probably indicative of either association with human activity or the airborne transfer of pollen. Petroselinum crispum (parsley) is a plant originally in the Mediterranean region but now widely cultivated around the planet and used for culinary purposes. Glycine soja (soybean), originally from Asia, is today one of the most important crop plants on the planet and is widely used as food (milk, tofu, oil, etc.). Eucalyptus fulgens (green scent bark) is a tree species endemic to Australia; however, this assignment could be an artefact arising from the quality or completeness of the database consulted, as the genus Eucalyptus includes > 700 species, many widely cultivated in the Americas, Europe, Asia and Africa for timber, cellulose or oil. The ruined wooden structures on Deception Island were primarily constructed of pine (Pinus sp.) and spruce (Picea sp.; Held et al. Reference Held, Arenz, Blanchette, Barr and Chaplin2011), genera that were not detected in the current study. However, even today many wooden materials are imported to Antarctica (Osyczka et al. Reference Osyczka, Mleczko, Karasiński and Chlebicki2012), and ~47% of the garbage found in Antarctica comprises or includes wood (Anfuso et al. Reference Anfuso, Bolívar-Anillo, Asensio-Montesinos, Portantiolo Manzolli, Portz and Villate Daza2020). Lityńska-Zając et al. (Reference Lityńska-Zając, Chwedorzewska, Olech, Korczak-Abshire and Augustyniuk-Kram2012) have also shown that many wood fragments arrive in Antarctica accidentally. Finally, the genus Prunus includes some of the most widely consumed food fruits, including cherries, almonds, plums, apricots and peaches, again suggestive of human influence.
Although we recognize that the possibility exists of pollen or plant propagules reaching Deception Island independently, taking into consideration the high number of visitors to the island (Roura Reference Roura2012, Carvalho-Silva et al. Reference Carvalho-Silva, Rosa, Pinto, Da Silva, Henriques, Convey and Câmara2021) and the ruderal nature of all of the flowering plant taxa assigned, we believe that there is a strong likelihood of the detection of flowering plant sequences here being a consequence of human activity. Carvalho-Silva et al. (Reference Carvalho-Silva, Rosa, Pinto, Da Silva, Henriques, Convey and Câmara2021), in a study of eDNA from Whalers Bay soil samples, similarly reported assignments of many taxa potentially linked with human activities, including foodstuffs.
It is always important to note that assignment of a DNA sequence from eDNA samples does not confirm the presence of the organism itself or viable parts of it. The identification of exotic eDNA highlights the need for protecting these ecosystems from inadvertent and irreversible genetic contamination. The method used here provided valuable and novel data on communities living on exotic decaying wood in Antarctica in addition to the better-studied Fungi and Bacteria.
Conclusions
Historical wooden ruins in Antarctica host a very diverse and largely unknown assembly of organisms, including some not previously reported from the region. DNA metabarcoding using HTS is a useful tool for detecting species that would be otherwise difficult or impossible to identify by morphology alone, including organisms that may be encysted or that cannot at present be grown in culture. The assignment of identities based on DNA sequences alone does not confirm the presence of living organisms or viable propagules, and the assignments themselves rely heavily on the quality of consulted databases, which has been increasing substantially in recent years. In addition, no single DNA marker is effective for all groups of organisms, and further studies using a wider range of markers, backed by traditional direct observation and culturing methods, will ultimately be required to provide a more thorough description of the diversity present. The data presented here will be useful in the development of conservation policies and in studies on colonization processes and human influence in Antarctica.
Author contributions
PEASC, MC-S, LdS and LR performed the fieldwork, PEASC, FAL and FLVB and performed the bioinformatic and data analysis, MC-S analysed the flowering plant data and PEASC analysed the non-flowering plant data. PEASC, PC and LR conceived the paper and its experimental design and provided funds. All of the authors participated in writing the paper.
Acknowledgements
We thank the Brazilian Navy and Air Force for logistical support and the Brazilian Antarctic Program (PROANTAR). We also thank the editor and the reviewers for handling this paper and working to make it better.
Competing interests
The authors declare no conflicts of interest.
Financial support
This study received financial support from the Brazilian Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), Ministério da Ciência, Tecnologia e Inovação (MCTI) and Programa Antártico Brasileiro (PROANTAR). PC is supported by Natural Environment Research Council (NERC) core funding to the British Antarctic Survey (BAS) ‘Biodiversity, Evolution and Adaptation’ Team. Congresswoman Jô Moraes and the Instituto de Ciências Biológicas at Universidade de Brasília provided extra funds.