Previous Article | Next Article ![]()
Applied and Environmental Microbiology, October 2004, p. 5963-5972, Vol. 70, No. 10
0099-2240/04/$08.00+0 DOI: 10.1128/AEM.70.10.5963-5972.2004
Copyright © 2004, American Society for Microbiology. All Rights Reserved.
British Antarctic Survey, Natural Environment Research Council, Cambridge CB3 0ET, United Kingdom
Received 21 January 2004/ Accepted 3 June 2004
|
|
|---|
|
|
|---|
The majority of previous studies of Antarctic eukaryotic microbiota in the terrestrial environment have been based on traditional morphological and/or culture approaches (6, 44, 61). However, the environmental range of some studies is relatively limited, and it is often difficult to compare results because of methodological variations (6). In recent years a selection of molecular methods which enable the identification of diverse community members in any given environment have become available. These methods do not require extraction of individuals or use of culture protocols and can be designed to target a very broad or very narrow range of organisms. To date, the majority of this molecular work has been directed towards prokaryote communities, often in extreme environments (38, 60, 66). Molecular studies have also been utilized to look at population differences in various environments or along stress gradients (10, 19, 20, 46). Recently, these methods have been applied to eukaryotic populations in a range of environments, although the majority have focused on aquatic microscopic eukaryotes (10, 18, 19, 31, 35, 37, 54, 65). These studies have identified unexpectedly high levels of diversity and suggest the possibility that a similarly hidden diversity of microscopic eukaryotes may exist in terrestrial environments.
The Antarctic forms a natural laboratory providing wide environmental gradients over which the effects of increasing stress on diversity and community structure can be observed. Conventionally, the region is separated into three terrestrial biogeographical zones (sub-Antarctic, maritime Antarctic, and continental Antarctic), based on consistent biological and climatic differences (14, 30, 48). The study described here utilizes molecular techniques to study diversity patterns across latitude, including material obtained across almost the entire latitudinal range of the maritime Antarctic, from the South Orkney Islands (ca. 60°S) to southern Alexander Island (ca. 72°S). Further data were obtained from the LaGorce Mountains (ca. 87°S), among the most southern terrestrial habitats existing in the continental Antarctic, and from Sky Hi Nunataks (ca. 75°S) in the currently undefined and largely unstudied "boundary region" between maritime and continental Antarctica. In addition to the obvious biogeographical significance of understanding processes across the boundary region between maritime and continental Antarctica, this study represents the first opportunity that biologists have had to work south of Alexander Island in the Antarctic Peninsula region.
The aim of this study was to utilize molecular methods to identify a wide range of eukaryotes from a selection of soils across the latitudinal and environmental gradient between the northern maritime and deep continental Antarctic, in order to test the hypothesis that diversity would decrease as latitude (and stress) increased. Prokaryote diversity in the soil community was not addressed in this study, due to the prohibitive number of clones that would have been required to perform an extensive survey. By using standard methods across all samples and study sites and targeting a specific type of soil habitat, we were able to assess the degree of similarity between the sites.
|
|
|---|
![]() View larger version (47K): [in a new window] |
FIG. 1. Map showing the study region and indicating locations of collection sites.
|
|
View this table: [in a new window] |
TABLE 1. Summary of the physical and chemical properties of soils obtained from each study site
|
Five separate soil samples (2 to 10 g) were collected from each site by one of two methods. If enough unfrozen soil was available, sterile monovette corers (Sarstedt Ltd., Leicester, United Kingdom) were plunged into soil fines to a depth of 10 to 15 mm. Alternatively, where little soil was available, larger particles precluded the use of monovettes, or the soil was frozen, samples were scraped into sterile whirlpack plastic bags with a sterile spatula. Soils were stored under ambient conditions in the field following collection before they could be returned to a research station, where they were then held at 20°C until use.
Physical and chemical analyses.
Physical analyses were undertaken to ascertain substratum similarities and any major differences. Soil samples (5 to 10 g) from each site were air dried at 80°C for 48 h, and the water content was estimated gravimetrically. A subsample of the air-dried soil was then incinerated at 400°C for 24 h to obtain organic matter content. The remainder of the air-dried sample was analyzed commercially (Centre for Ecology and Hydrology, Merlewood, United Kingdom) for pH and contents of C, N, P, K, Fe, As, Cd, Cr, Ni, Pb, V, and Zn. In order to provide sufficient material, soil from the five samples obtained at each site was combined, precluding statistical comparisons between sites.
Estimation of molecular diversity. (i) DNA extraction, amplification and clone library construction.
In order to provide comparable data across a wide range of organisms from all sites, a clone library-sequence approach was taken. Clone libraries were constructed by extracting total genomic DNA, amplifying the near-full-length eukaryotic small-subunit ribosomal DNA (SSU rRNA gene), and then cloning the products. Total genomic DNA was extracted separately from the five individual soil samples obtained from each site, and each of the five isolates was then used as a template in at least three amplification reactions. DNA was extracted from soil samples by using a combined physical and chemical treatment. Approximately 0.5 g of soil was placed in an Eppendorf tube with 500 µl of sterile 0.1-mm-diameter glass beads (Stratech Scientific, Soham, United Kingdom) and 700 µl of lysis buffer (0.1 M Tris [pH 7.4], 0.25 M EDTA, 2% [wt/vol] glucose, 0.5% [wt/vol] sodium dodecyl sulfate, and 2 mg of proteinase K [Sigma, Gillingham, United Kingdom] ml1). Samples were vortexed at full speed for 5 min and then subjected to five cycles of freezing and thawing (80°C for 2 min and +80°C for 2 min). Cell debris was then removed by centrifugation (12,000 x g, 2 min). The supernatant was phenol extracted, and the DNA was ethanol precipitated and resuspended in 400 µl of water. Aliquots from each sample (80 µl) were cleaned for PCR amplification with GFX DNA purification columns according to the instructions of the manufacturer (Amersham Pharmacia Biotech, Little Chalfont, United Kingdom).
PCR amplifications were performed with 43 µl of ABgene 1.1x ReddyMix PCR master mix (ABgene, Epsom, United Kingdom), 20 pmol of each primer, and approximately 10 to 100 ng of template DNA in a final volume of 50 µl. The amplification program consisted of 1 cycle of 95°C for 5 min; 30 cycles of 94°C for 60 s, 55°C for 60 s, and 72°C for 70 s; and a final extension step of 72°C for 10 min. To amplify the near-full-length eukaryotic SSU rRNA gene, the primers E4F (5'-CTGGTTGATTCTGCCAGT-3') and E1628R (5'-CGACGGGCGGTGTGTA-3') (53) were used.
PCR products from the same sample were pooled and purified by using GFX DNA purification columns. Cleaned products originating from each sample from a single site were pooled and cloned into the pGEM-T Easy vector (Promega, Southampton, United Kingdom). Ligation products were transformed into JM109 (Promega) according to the manufacturer's instructions. Transformants were screened by color selection on Luria agar containing 100 µg of ampicillin ml1, IPTG (isopropyl-ß-D-thiogalactopyranoside), and either S-Gal or X-Gal (5-bromo-4-chloro-3-indolyl-ß-D-galactopyranoside) (Sigma) as an indicator. White colonies were picked from the plates and resuspended in 50 of µl water. The cells were subjected to two freeze-thaw cycles (80°C and +80°C), and an aliquot of 2 µl was used as the template in a screening PCR with M13 forward (TGTAAAACGACGGCCAGT) and reverse (CAGGAAACAGCTATGACC) universal primers. The amplification program consisted of 1 cycle of 95°C for 5 min; 30 cycles of 94°C for 45 s, 58.5°C for 45 s, and 72°C for 60 s; and a final extension step of 72°C for 10 min. Aliquots of the amplified products were electrophoresed in 1% (wt/vol) agarose gels, and products showing bands of the correct size were cleaned with GFX DNA purification columns.
(ii) Controls.
The methods described above are potentially open to contamination at a number of stages, and thus a series of controls were put in place. First, all samples were taken aseptically, and aliquots for DNA extraction were removed from the sampling containers in flow hoods. The hoods and balances were washed down with 70% ethanol and Virkon (main active ingredient, potassium monopersulfate; Antek International, Sudbury, United Kingdom) between samples. A control DNA extraction (containing no soil) was carried out in parallel with the soil DNA extractions and was used as the template in subsequent PCR amplification. Negative control reactions with mixtures containing water in place of genomic or plasmid DNA were also carried out at each step of the protocol. None of these control reactions generated positive results.
(iii) RFLP screening and sequencing.
Cleaned PCR products of the correct size were separately digested with HaeIII and RsaI for 2 to 3 h at 37°C. Digests were separated on 3% (wt/vol) high-resolution agarose (Appligene, Middlesex, United Kingdom) and stained with ethidium bromide, and images were captured and exported to the Gelcompar II software package (Applied Maths, Saint-Martens-Latem, Belgium). Patterns from all gels were normalized against a size standard, and all clones showing unique patterns were sequenced. Where groups of clones displayed the same restriction fragment length polymorphism (RFLP) patterns, two to four isolates were chosen, at random, for sequencing. Sequence reactions were carried out with the ET-Terminator sequence kit (Amersham Biosciences, Little Chalfont, United Kingdom) and M13 forward and reverse primers and were run on a Megabace 500 sequencer (Amersham Biosciences). For the purpose of this study, sequences corresponding to the first 700 bp of the eukaryotic SSU rRNA gene were used.
Operational taxonomic units (OTUs) were defined as clones sharing
97.5% sequence identity, regardless of whether or not they shared RFLP patterns. This value has been suggested for bacterial studies (3, 36). For eukaryotes no comparable cutoff has been proposed, and likely values appear to vary considerably depending on the DNA region used and the organisms studied. Hebert et al. (26) found that intraspecific divergences were rarely greater than 2% for the mitochondrial cytochrome c oxidase subunit I gene in animals, and a value of 0 to 3% divergence between species of basidiomycete fungi has been reported for internal transcribed spacer and 5.8S regions (64). Generally, the small-subunit regions of eukaryotes are very conserved, but given the potential for errors in both sequencing and amplification protocols from environmental samples (9), this value should ensure that differences are representative of biological differences and are not artifactual. Samples were taken until the sampling coverage (see below) for each clone library was greater than 85%.
(iv) Sequence identification.
All sequences were analyzed with the rRNA Database Project CHECK_CHIMERA program (32); six potential chimeric sequences were detected and removed from the analysis. All remaining clone sequences were compared with the GenBank database of genetic sequences by using gapped BLAST searches (1) to determine their approximate phylogenetic affiliations. Several sequences were most closely related to prokaryotes, and some sequences had no affiliations to any sequences in the database. These did not contain any of the SSU rRNA gene conserved regions and, along with the prokaryote sequences, were not included in further analyses. Sequences obtained from each site were aligned through the JalView utility of ClustalW run at the European Bioinformatics Institute (http://www.ebi.ac.uk), and a consensus phylogenetic tree was obtained from DNA distances and neighbor joining in PHYLIP. FASTA searches were then used to compare sequences against the appropriate EMBL subset database (e.g., EMBL-FUNGI) indicated by the position at which the sequence was recovered in the tree and the BLAST results. Sequences were identified to at least the subdivision level.
(v) Statistical analyses.
Library coverage, an estimate of the proportion of a clone library of infinite size that has been sampled, was calculated by using the relative distribution of OTUs and the equation described by Good (24). Collectors' curves, or species abundance curves, were constructed to compare the diversity of OTUs at each site. Curves were calculated by using the freeware program Analytic Rarefaction 1.3 (available at http://www.uga.edu/
strata/software/). The program uses the rarefaction equations described by Hurlbert (28) and Heck et al. (27). The probability of drawing a new OTU (
n), on the next draw was estimated by using the method described by Clayton and Frees (13).
Between-site relationships based on numbers and types of taxa present were investigated by principal-coordinate analyses. Data in each group were transformed into presence-absence data and converted to a percentage of sequences recovered at each site. A distance matrix was constructed with the average taxonomic distance coefficient, and a principal-coordinate plot was made based on the first three axes. All computations used the MVSP package (Kovach Computing, Anglesea, United Kingdom).
Pairwise comparisons of sites were made by using Z scores, which indicate whether an estimated number of OTUs (E) at a given number of clones (n) is significantly different between pairs of sites. E was obtained by normalizing libraries at n = 77 by using rarefaction (as described above). The rarefaction calculations also gave the upper and lower 95% confidence levels, which were used to calculate the standard error for each site.
|
|
|---|
Clone libraries.
Clone library coverage ranged from 87 to 99%, and all libraries showed very low
n values. To reach this level of coverage, between 78 (Sky Hi Nunataks) and 169 (Signy Island) clones were required (Table 2). Across all sites, 712 clones produced 204 RFLP patterns that were divided into 111 OTUs. Several OTUs were comprised of more than one RFLP pattern, and in two cases the same RFLP pattern gave rise to different OTUs. Alexander Island sites (Mars Oasis and Coal Nunatak) had numbers of OTUs similar to or higher than those of northern maritime Antarctic Signy Island, while Rothera Point (central maritime Antarctic) generated fewer. Sky Hi Nunataks and LaGorce Mountains generated the lowest numbers of OTUs (Tables 2 and 3). The latter libraries also contained the lowest number of clones and the highest coverage values, supporting the conclusion that they are the least diverse of the sites studied.
|
View this table: [in a new window] |
TABLE 2. Estimates of sampling coverage and n (the unconditional probability of a new OTU discovery if another clone was taken) for each study sitea
|
|
View this table: [in a new window] |
TABLE 3. Identification and grouping of sequences. The table lists the major higher taxonomic groupings (phylum or division) into which sequences obtained were separated, the number (N) of separate operational taxonomic units (OTUs) obtained in each group from each study site, and the overall percentage contribution made by each kingdom or sub-kingdom to the biota of each study site.
|
![]() View larger version (14K): [in a new window] |
FIG. 2. Example of consensus tree of aligned complete sequences, showing major kingdom and subkingdom groups present in Signy Island samples. Bootsrap replications are shown for major groupings.
|
|
View this table: [in a new window] |
TABLE 4. Taxonomic breakdown by class of OTUs identified at each study site by comparison with the GenBank sequence databasea
|
Site comparisons.
Rarefaction curves created for each of the sampling sites are shown in Fig. 3. A plateau, indicating close to full biodiversity coverage, was approached only for the LaGorce and Sky Hi Nunataks soils. While clone library coverage values were also high for the other sites (Table 2), the curves did not reach a plateau. The slopes of the Signy Island and Rothera Point curves were shallower than those of the curves obtained for Mars Oasis and Coal Nunatak, suggesting that a small number of OTUs were dominating these libraries. Indeed, the Signy Island clone library contained a single OTU corresponding to a nematode worm sequence that accounted for 55% of clones. No other library showed such a level of dominance by a single OTU; however, all had between one and three OTUs which each contributed 10% or more of the library. Removal of the dominating nematode OTU from the Signy Island data generated a rarefaction curve that was more similar to those of Mars Oasis and Coal Nunatak, although the number of OTUs found remained low in comparison (Fig. 3). Rarefaction was used to normalize each of the clone libraries, and comparisons were made between the sites based on the number of OTUs observed (Table 5). This comparison showed the two Alexander Island sites (Mars Oasis and Coal Nunatak) to be more diverse than the other maritime Antarctic sites, while the boundary region site at Sky Hi Nunataks was less diverse and the continental site at LaGorce Mountains had by far the lowest diversity.
![]() View larger version (17K): [in a new window] |
FIG. 3. Rarefaction curves constructed for eukaryotic clone libraries from each of the six collection sites. Clones were grouped into OTUs at a level of sequence similarity of >97.5%. The "Signy without nematodes" curve has had a dominating clone removed (see text).
|
|
View this table: [in a new window] |
TABLE 5. Diversity comparisons between the six study sites
|
![]() View larger version (9K): [in a new window] |
FIG. 4. Principal-coordinate analysis of sites based on the percentage of sequences obtained for each taxonomic group at each site (original data are in Table 3).
|
![]() View larger version (12K): [in a new window] |
FIG. 5. Numbers of sites at which different OTUs were found. Fungal, algal, and nonalgal protist groups and a combination of all groups (including metazoans, plants, and uncultured clones) are shown. Numbers above the bars indicate the total number of OTUs in each group found at a given number of sites.
|
|
|
|---|
Identification of sequences obtained is dependent on the quality of the sample data, the available reference data, and the discrimination accorded by the sequence used. In this study sample, sequence quality was good (sequences of ca. 1,000 bp were routinely recovered) and was clearly sufficient to separate OTUs into the major biological kingdoms and subkingdoms. Below this taxonomic level, discrimination was more variable, with some species and genera being very distinct, while other assignments could be made only at a family or order level. This illustrates a limitation in using a single sequence to screen and identify a wide diversity of organisms. Reference data quality remains problematic. Mislabeled and incorrect sequences appear in databases (5), while the proportion of organisms for which sequences are available is still relatively low, although coverage varies between taxonomic groups (see, e.g., reference 5). Our cutoffs of 97.5% homology for OTUs and identification to order level only are justified by the limitations of the reference information.
Latitudinal patterns.
It is generally accepted that Antarctic terrestrial diversity decreases as latitude (and stresses) increase (6, 29, 47). In contrast, the present data provide support for this pattern only when a large (regional)-scale comparison is made between the maritime and continental Antarctic. Instead, the two southern maritime Antarctic sites (Mars Oasis and Coal Nunatak) generated diversity measures comparable to or higher than those of the central and northern sites (Rothera Point and Signy Island). Greatest diversity was recorded for Coal Nunatak. While Coal Nunatak is only ca. 30 km from Mars Oasis and they are part of the same sedimentary geological formation, the two sites contrast markedly in their degree of colonization by macroscopic biota. Mars Oasis, which is ca. 30 m above sea level, has a relatively diverse fauna and flora (16), while those of Coal Nunatak, which is ca. 750 m above sea level, are extremely limited. The finding that the eukaryotic microbial diversity of Coal Nunatak soils is approximately 50% greater than that of Mars Oasis (Table 3) is, at least, surprising and may suggest a longer duration of exposure through recent glacial cycles.
The macroscopic biota of Alexander Island are closely related to those of the maritime Antarctic (16), as are the soil meiofauna as indicated by the very few published records (34). However, there are also endemic components, including a springtail and several currently undescribed nematode and tardigrade species (16, 34; N. R. Maslen and P. Convey, unpublished data; S. J. McInnes, personal communication). Although it is generally thought that contemporary terrestrial habitats of this island have been exposed for, at most, 5,000 to 6,000 years (16, 51), when combined with the current finding of high microeukaryotic diversity, these various observations suggest the existence of terrestrial habitats beyond the time scale allowed by the last glacial maximum.
Lower diversity was present at the two more southern sites examined, at Sky Hi Nunataks and the LaGorce Mountains. Given that the present study could include only these two sampling points, it is not possible to conclude whether their lower diversity is a function of increasing latitude or is simply representative of a different Antarctic biogeographical zone. The range of organisms isolated from the LaGorce soils is similar to that found by Broady and Weinstein (8), who used microscopy and culture-based methods to identify cyanobacteria, chlorophytes, and fungi in the same soils.
The organisms recovered are also likely to be influenced by the environmental conditions at each site. For instance, at the sites with the lowest availability of liquid water (LaGorce Mountains and Sky Hi Nunataks), largely water-associated stramenopilan taxa were not present or were poorly represented. Unlike in temperate and tropical solids, where basidiomycete and zygomycete fungi generally dominate, in this study ascomycetes predominated, a group that includes many lichenized taxa (Table 4). As the primarily lichen-forming chlorophyte Trebouxiaceae also contributed a high proportion of the sequences obtained, this suggests that the soils sampled were strongly lichen influenced (e.g., by propagules, early colonizing stages). Samples from four sites (Rothera Point, Signy Island, Coal Nunatak, and La Gorce Mountains) contained sequences that were conclusively identified as being plasmodiophoran. This group has not previously been reported from Antarctica.
Endemism or cosmopolitanism?
Comparison of molecular and traditional surveys is difficult. In studies of prokaryotes, they often show very little overlap, although greater similarity is seen in more extreme and less diverse environments (4, 40, 49). Some eukaryote groups (e.g., nematodes and arthropods) are relatively well documented in the Antarctic (2, 14). However, although sequences showed affiliation to nematodes or arthropods in the GenBank database, none of the latter species or their close relatives are present in Antarctica. Classical taxonomic data suggest that many metazoans are endemic in the continent (2, 25, 41, 42). Thus, while this study has most likely detected the presence of known Antarctic species, it is currently impossible to match the sequence data with known taxa. These difficulties are compounded for some groups (particularly the protists) by a lack of consistency within classical taxonomy.
Endemism, especially applied to microorganisms, is a controversial topic. The term generally implies that an organism is native or confined to a certain region, or, as expressed by Vincent (55) "genotypes of bacteria, protists or other microorganisms specific to a geographical region." Vincent further postulates that if microbial endemism is possible, the Antarctic should be among the most likely places in which such organisms may be found. However, there is currently limited evidence for microbial endemism in Antarctic terrestrial, or indeed any, polar environments (6, 22, 50, 56).
The alternative to endemism is cosmopolitanism. It has long been suspected that microbial dispersal is not restricted by geographical boundaries (21, 63): "everything is everywhere, but the environment selects" (Baas-Becking [1934], cited in reference 50). Darling et al. (17) report evidence for bipolar genetic flow but do not confirm that cosmopolitan microorganisms exist in both polar regions. However, bipolar fungi, lichens, and mosses are known (48). Microorganisms may be transported into and around the Antarctic via various routes, including air currents (33), ocean currents (39), and animals (including birds [45] and humans [7]). Propagule banks in soil or snow are reported to contain cosmopolitan species of algae, fungi, protists, bacteria, and archaea (39, 52, 56, 61). However, while many of the organisms detected in our study show reasonable homology with isolates from other geographic locations, it remains impossible to categorically confirm that they are conspecific.
The clone library data from this study suggest that a minority of the organisms may be ubiquitous (Fig. 5), as a small number of OTUs were identified at multiple, widely separated sites and/or are closely related to sequences from cosmopolitan organisms (e.g., the alga Stichococcus bacillarus and the protists Heteromita globosa and Oxytricha granulifera). Many other isolates also show a high degree of similarity (>98% sequence identity) with non-Antarctic organisms, although 18S ribosomal gene sequence similarity does not necessarily imply functional similarity and species identity. Conversely, some sequences are considerably different from any others present in the GenBank database, an observation which may be indicative either of endemism or simply that these organisms are yet to be detected or sequenced in other environments. Although proof of endemism is not possible, our results suggest that there is very little effective transfer of biota between the study sites. Further insight will require isolation and genetic analysis of identified Antarctic organisms and comparison with known Antarctic and non-Antarctic relatives.
This study forms part of the BAS "Biological Responses to Environmental Stress in Antarctica" core project and also contributes to the SCAR RiSCC (Regional Sensitivity to Climate Change in Antarctica) Program.
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»