ABSTRACT
Archaeal 16S rRNA was extracted from samples of deep marine subsurface sediments from Peru Margin site 1227, Ocean Drilling Program leg 201. The amounts of archaeal 16S rRNA in each extract were quantified by serial dilution and reverse transcription (RT)-PCR. The results indicated a 1,000-fold variation in rRNA content with depth in the sediment, with the highest concentrations found near the sediment surface and in the sulfate-methane transition zone (SMTZ). The phylogenetic composition of the active archaeal population revealed by cloning and sequencing of RT-PCR products changed with depth. Several phylotypes affiliated with marine benthic group B (MBGB) dominated clone libraries from the upper part of the SMTZ and were detected only in this layer. Members of the miscellaneous crenarchaeotal group (MCG) dominated clone libraries from the other layers. These results demonstrate that archaeal communities change in activity and community composition over short distances in geochemically distinct zones of deep subseafloor sediments and that these changes are traceable in the rRNA pool. It was shown for the first time that members of both the MCG and MBGB Archaea are more active in the SMTZ than in layers above and below. This indicates that they benefit either directly or indirectly from the anaerobic oxidation of methane. They also appear to be ecophysiologically flexible, as they have been retrieved from a wide range of marine sediments of various geochemical properties.
Results obtained through the Deep Sea Drilling Project and Ocean Drilling Program (ODP) have shown that marine subsurface sediment constitutes one of the largest and most widespread reservoirs of biomass on Earth and that the activities of subsurface prokaryotes have profound implications for the global carbon cycle (9, 10, 11, 29, 37, 38, 55). Numerous studies have documented diverse prokaryotic communities in subsurface sediments by extraction and analysis of genomic DNA material (17, 35, 41, 45), but since most of the obtained phylotypes belong to previously unknown and uncultivated groups, the physiological properties of the subsurface organisms remain unknown. DNA may persist in intact but inactive cells as well as in the environment as extracellular material after cell death (8, 34). Since DNA-based gene surveys do not distinguish between living and dead populations, it has remained unclear whether the detected phylotypes reflect metabolically active organisms. This problem has led to the “paleome” concept, wherein deep-subsurface sediments and marine sedimentary rocks are viewed as repositories for inactive cells and ancient DNA (16, 18), and to ensuing debate about this interpretation of microbial community DNA signatures (14). In contrast to genomic DNA, ribosomes are being continuously turned over in cells and the concentration of rRNA is correlated to the growth rate (23, 24, 25, 28, 30). Thus, by analyzing rRNA rather than genomic DNA in sediment samples, it is possible to specifically target the active part of the microbial community and to more directly correlate biogeochemical gradients with the activity of specific phylogenetic groups (32, 33, 52, 53). This could be of particular value in subsurface sediments where a large proportion of the cells may be inactive.
Between one-third and one-tenth of the prokaryotic cells in subsurface sediment samples from ODP site 1227 at the Peru Margin contain sufficient 16S rRNA for staining with catalyzed reporter deposition-fluorescence in situ hybridization (CARD-FISH). Active Archaea were not detected with this approach, but the presence and distribution of archaeal genes for 16S rRNA was documented through quantitative PCR amplification (42). In another study, archaeal genes for 16S rRNA were amplified and characterized from four different depths in the sediment column at site 1227 (19). The detected phylotypes were affiliated with a number of uncultured groups of Archaea, and their physiological status and ecological significance could not be determined. Recently, a combination of intact phospholipid (IPL) extraction, FISH, and rRNA analysis demonstrated that Archaea dominate the active prokaryotic community in subsurface sulfate-methane transition zones (SMTZ) from several different ODP sites on the Peru Margin, including site 1227, where up to 98% of the detected cells hybridized with an Archaea-specific oligonucleotide probe (3). However, depth-related changes of the active archaeal community, for example, between zones of net methanogenesis and methane consumption, were not addressed.
In the present study, the active archaeal communities were characterized throughout the upper 50 m of sediment at Peru Margin site 1227 by extraction and analysis of rRNA. By selecting for metabolically active organisms, changes in archaeal activity and community composition were shown to be related to the strong geochemical gradients that exist at this site.
MATERIALS AND METHODS
Study site and sampling.ODP site 1227 is located on the Peru continental shelf within the upwelling zone in 427 m of water (44). The sediment column is characterized by strong geochemical gradients in concentrations of dissolved sulfate and dissolved methane (Fig. 1) (44). The concentration of dissolved sulfate declines rapidly in the upper 15 m of sediment and more slowly below this level, down to about 41 m below the seafloor (mbsf), where sulfate is depleted. Methane accumulates to several millimolar concentrations below the sulfate-containing zone and is consumed in the layers between 35 and 41 mbsf, presumably with sulfate serving as the electron acceptor. The sediment can thus be divided into three zones, each characterized by different overall methane/sulfate dynamics: an upper zone from 0 to 35 mbsf, where the sulfate concentration decreases gradually with depth and where methane is present only in trace amounts; a transition zone from 35 to 41 mbsf, where sulfate is depleted by microbially mediated reactions with methane; and a deep methanogenic zone below 41 mbsf. Samples were collected during ODP leg 201 in 2001 under strict contamination controls, as have been described elsewhere (15).
Concentration profiles of sulfate (black squares) and methane (open triangles) at site 1227 (redrawn from reference 17).
Extraction of RNA.RNA was extracted from samples of boreholes 1227A (2H-2, 7.35 m; 3H-5, 21.35 m; 5H-3, 37.35 m; 5H-5, 40.35 m; 6H-2, 45.35 m; 6H-5, 49.85 m) and 1227D (1H-5, 6.55 m; 4H-6, 34.25 m; 5H-2, 37.75 m). About 4 ml of sediment was retrieved from the center of each sample core with sterile spatulas and immediately mixed with 10 ml of a 1:1 mixture of extraction buffer (250 mM sodium acetate, 5 mM EDTA, 2% sodium dodecyl sulfate [pH 5]) and phenol (pH 5). The mixture was distributed into 2-ml bead-beating tubes containing 0.5 g of 0.1-mm zirconium beads (BioSpec Products, Inc., Bartlesville, OK). The tubes were homogenized on a FastPrep FP120 homogenizer (Qbiogene, Inc., Carlsbad, CA) for 30 s at level 6.5 and subsequently centrifuged at 10,000 × g for 5 min at 4°C. The aqueous phases were transferred to new 1.5-ml centrifuge tubes and kept on ice. The phenol phases, containing sediment and beads, were subjected to two more rounds of homogenization, centrifugation, and transferal of the water phases, each time with 300 μl of new extraction buffer. The combined aqueous phases from the three rounds of homogenization were extracted three times by vortexing and centrifugation with equal volumes of phenol, phenol-chloroform, and finally chloroform. Extracted RNA was precipitated for 2 h at −20°C with 0.5 volumes of 7.5 M ammonium acetate and 1 volume of isopropanol. The extract was spun down, washed with 80% ethanol, air dried, and finally resuspended 86 μl of water. To digest trace amounts of DNA in the sample, 10 μl of 10× DNase buffer and 4 μl of DNase (Fisher) was added and the sample was digested for 30 min at 37°C. The sample was finally purified with the RNeasy minikit (QIAGEN Inc.) according to the manufacturer's instructions. Parallel blind extractions, in which the entire extraction procedure was performed without addition of sample, were consistently performed, and in case the control yielded a visible reverse transcription (RT)-PCR product, the extraction was redone.
RNA quantification.Each extraction product was serially diluted between 1 and 2,187 times in a 3× dilution series. Template (2 μl) from each dilution was reverse transcribed and amplified in 50-μl reaction mixtures with the archaeal primers ARC8f (5′-TCCGGTTGATCCTGCC-3′) and either ARC915r (5′-GTGCTCCCCCGCCAATTCCT-3′) or ARC518r (5′-GGT(G/A/T)TTACCGCGGC(G/T)GCTG-3′) using the One-Step RT-PCR kit (QIAGEN Inc.). The RT-PCR program consisted of 30 min at 50° and 15 min at 95°C, followed by 35 cycles of 30 s at 92°C, 30 s at 58°C, and 4 min at 72°C. The final elongation step was increased to 10 min. To test for the presence of trace amounts of DNA in the RNA extracts, the extracts were used as templates in PCRs. No PCR products were observed. Upon amplification, the entire volume of each RT-PCR product was loaded onto a 1% low-melting-point agarose gel, run at 70 V for 40 min, and stained with SYBR gold (Bio-Rad). DNA on the gels was visualized with a Dark Reader transilluminator (Clare Chemicals, Dolores, CO), and the smallest amount of template from which a visible amount of DNA could be obtained was recorded. This served as a relative measure of the amount of RNA in the extraction product. Coextraction of unknown PCR inhibitors hindered PCR amplification of undiluted template amounts larger than 2 μl. DNA bands were excised and extracted from the gel with the QIAEXII gel extraction kit, according to the manufacturer's instructions.
Phylogenetic analysis.The excised RT-PCR products were cloned as described previously (45). Sequence data were obtained from the sequencing facility at the Marine Biological Laboratory (Woods Hole, MA). Sequences of more than 98% similarity were regarded as identical phylotypes. The obtained phylotypes were BLAST analyzed against the GenBank 16S rRNA database (1). The sequences were aligned with representative members of the Euryarcheota and Crenarchaeota with the CLUSTAL W program, followed by manual alignment. Bootstrapping of the data set, calculations of Jukes-Cantor distance matrices, and tree construction by neighbor joining or maximum parsimony was performed with the PHYLIP software package (12). The sequence data set was screened for potential chimeric structures by using the CHECK_CHIMERA application at the Ribosomal Database Project website and by constructing phylogenetic trees from three different fragments of the sequence alignment (bp 1 to 250, bp 251 to 500, and bp 501 to the end of the sequence). Sequences found to be of potential chimeric origin either through the screening or due to inconsistent positions in the phylogenetic trees were excluded from this study.
Diversity analysis.The total number of phylotypes in each clone library was estimated by calculating SChao1 estimators (5, 6): $$mathtex$$\[S_{\mathrm{Chaol}}{=}S_{\mathrm{obs}}{+}\frac{F_{1}^{2}}{2(F_{2}{+}1)}{+}\frac{F_{1}\ {\cdot}\ F_{2}}{2(F_{2}{+}1)^{2}}\]$$mathtex$$ where Sobs is the number of different phylotypes detected in a clone library and F1 and F2 are the numbers of phylotypes encountered one and two times, respectively. The estimated number of phylotypes depends on the number of phylotypes detected and is prone to underestimation in insufficiently sampled datasets. To evaluate the stability of the obtained phylotype richness estimates, SChao1 estimators were calculated for randomly sampled subsets, including one-fourth, one-half, or three-fourths of the clone libraries (22).
Nucleotide sequence accession numbers.The phylotypes detected in this study are available from the GenBank nucleotide database, at the National Center for Biotechnology Information (NCBI) website (http://www.ncbi.nlm.nih.gov ), under accession numbers DQ301959 to DQ302007.
RESULTS
RNA quantification.RNA was successfully extracted, reverse transcribed, and amplified from seven of the nine samples (Fig. 2a). The amounts of archaeal 16S rRNA extracted from the sediment samples were highly dependent on depth in the sediment column. The concentration decreased with depth down to 34 mbsf, where no rRNA could be amplified, indicating that the concentration was at least 1,000-fold lower than at the surface. The amount of 16S rRNA then increased through the SMTZ and finally decreased again at greater depths. Also shown in Fig. 2a are relative abundances of Archaea detected in a previous study by quantitative PCR (Q-PCR) for 16S rRNA genes at site 1227 (42). A general trend of decreasing numbers of Archaea with depth and a local maximum at the SMTZ was also found with the Q-PCR approach (42). However, the relative differences between layers were smaller by about an order of magnitude, and the increased rRNA concentrations at 37 mbsf were not registered in the Q-PCR analysis.
Results of 16S rRNA extraction and analysis. (a) Relative amounts of 16S rRNA in extracts from different depths (black squares). A value of 1 corresponds to the detection limit. The data are compared with the amounts of archaeal genes for 16S rRNA extracted from similar depths in the sediment and quantified by Q-PCR (open circles) (redrawn from reference 42). These numbers were normalized such that the amount of DNA found at 37.75 mbsf corresponds to a value of 1. The SMTZ is indicated in the figure by the stippled lines. (b) Phylogenetic composition of RT-PCR clone libraries from each depth. The sample from 37.75 mbsf was extracted twice; the second extract was amplified with two primer sets, yielding three independent clone libraries. Groups shown include MBGB (54), MCG (17), and SAGMEG-1 (51).
Phylogenetic profiling.Clone libraries were constructed from the reverse-transcribed and PCR-amplified archaeal rRNA. The distribution of clones between major phylogenetic groups is summarized in Fig. 2b. Neighbor-joining trees outlining the phylogenetic positions of the detected phylotypes and close relatives are shown in Fig. 3 and Fig. 4. The phylogenetic profile of the active archaeal community changed substantially with depth. In the sediment above the transition layer, the clone libraries were dominated by phylotypes belonging to the miscellaneous crenarchaeotic group (MCG) (17, 50). Smaller numbers of clones were affiliated with South African gold mine euryarchaeotic group 1 (SAGMEG-1) (51) or other groups of uncultured Archaea, including the VAL III (Valkea Kontinen) group (20), the terrestrial miscellaneous euryarchaeotal group (TMEG) (51), marine benthic group D (MBGD) (54), and a cluster of SAGMEG-1 and -2 relatives (51). In the samples from 37.35 and 37.75 mbsf, both located within the SMTZ, the clone libraries were dominated by members of marine benthic group B (MBGB) (54), a deeply rooted sister lineage to the Crenarchaeota that is synonymous with the deep-sea archaeal group (DSAG) (50). Samples from the bottom of the transition zone and below yielded exclusively phylotypes of MCG Archaea. The MCG-MBGB switch in the active archaeal community of the SMTZ was tested and reproduced twice. A second extraction of the sediment from 37.75 mbsf, RNA quantification, and sequencing of clones obtained with the same primers as before (ARC8f and ARC915r) confirmed the original results. To further evaluate the robustness of the results, a second clone library from the material of this second extraction was constructed with an alternative reverse PCR primer, ARC518r. In general, the same phylotypes were detected whether reverse primer ARC915r or ARC518r was used, but the relative frequencies of some phylotypes in the clone libraries varied depending on the primer used (Fig. 3). A new archaeal group, the ancient archaeal group (AAG), was detected with the alternative reverse primer. This group is deeply rooted within the Euryarchaeota and has previously been found at a deep-sea hydrothermal vent (49). Published sequences of AAG Archaea show several mismatches in the target region of primer ARC915r, suggesting that the group was not detected previously due to primer mismatches.
Neighbor-joining tree showing phylogenetic positions of the crenarchaeotal phylotypes retrieved during this study. Phylotype names refer to sample 1H-5 (6.55 mbsf), 2H-2 (7.35 mbsf), 3H-5 (21.35 mbsf), 5H-2 (37.75 mbsf), 5H-3 (37.35 mbsf), 5H-5 (40.35 mbsf), or 6H-2 (45.35 mbsf). Some phylotypes were retrieved from several layers, as indicated. Numbers in parentheses indicates the number of times found in clone libraries. The three numbers given at each phylotype from sample 5H-2 refer to the first extraction, second extraction, and second extraction using an alternative reverse primer for PCR. Bootstrap values equal to or larger than 0.5 are indicated. The values were determined using 1,000 replicates.
Neighbor-joining tree indicating the phylogenetic positions of the euryarchaeotal phylotypes found during this study.
Most of the phylotypes detected in this study were closely related or nearly identical to phylotypes found during 16S rRNA gene surveys from Peru Margin sites, as well as other subsurface marine sediments (19, 38). Several well-supported subclusters (bootstrap values 0.9 to 1.0) were identified within the MCG (Fig. 3). One of these, MCG-1, contained all phylotypes detected in and below the SMTZ. Phylotypes related to known methanogens or putative anaerobic methane oxidizers of the ANME-1 and -2 groups (4, 26, 36) were not detected in any sediment horizon.
Richness analysis.The total numbers of phylotypes were estimated to be 246, 26, and 10 in clone libraries from samples 1H-5 (6.55 mbsf), 2H-2 (7.35 mbsf), and 3H-5 (21.35 mbsf), respectively. To test the stability of these numbers, similar calculations were performed with randomly selected subsets of one-fourth, one-half, and three-fourths of the data. The results are shown in Fig. 5. No plateau was reached with sample 1H-5, suggesting that the total diversity in this clone library may be higher than the estimated 246 phylotypes (22). In contrast, estimated phylotype numbers were stable when 50% or more of the clone libraries from sample 2H-2 or 3H-5 were used, suggesting that phylotype richness had reached saturation in these samples. The number of different phylotypes in the clone libraries thus declined rapidly with depth above the SMTZ.
S Chao1 estimates of total numbers of phylotypes in clone libraries from cores 1H-5, from 6.55 mbsf (squares); 2H-2, from 7.35 mbsf (circles); and 3H-5, from 21.35 mbsf (triangles). Estimators were calculated from randomly sampled subsets as well as from the entire data set from each sample (22). Error bars indicate standard deviations in richness estimates for four resamplings.
DISCUSSION
Archaeal communities in subsurface sediments.The archaeal communities in deep subsurface sediments consist largely of unknown groups of organisms and are only distantly related to cultured strains. With the exception of marine group I (13, 21, 27, 40), organisms belonging to the groups in Table 1 have been detected solely in the form of 16S rRNA gene phylotypes. Some of the groups, including MBGB, MCG, and SAGMEG-1 and -2, appear to be ubiquitous in subsurface environments, while other groups have been detected more sporadically. Since cells may be transported into and within sediments by porewater flow, and since organisms deposited on the sediment surface may become buried in deeper sediment layers during ongoing sedimentation, the distribution of gene phylotypes may or may not be correlated to metabolic activity in the sediments. As a consequence, it has not been clear whether subsurface phylotypes reflect living and actively metabolizing organisms. The finding of RT-PCR-amplifiable rRNA affiliated with MBGB, MBGD, MCG, TMEG, and SAGMEG shows that members of these groups are active in the subsurface sediment at site 1227. As is clear from this and other studies (Table 1), archaeal communities vary considerably from site to site; future research will undoubtedly demonstrate activity of further groups of organisms within subsurface sediments.
Frequencies of uncultured archaeal lineages in clone libraries of genes for 16S rRNA from marine subsurface sediments
MCG.Members of the MCG have been found frequently in 16S rRNA gene surveys of both marine and terrestrial environments, and they appear to be widespread in both surface and subsurface systems (2, 7, 17, 26). The group dominated 16S rRNA gene clone libraries from ODP Peru Margin sites 1227 and 1229, ash layers of Okhotsk Sea sediments, and some sediment layers at the Nankai Trough (Table 1). The majority of the MCG-related rRNA phylotypes found by RT-PCR were closely related or even identical to phylotypes obtained by PCR amplification of genomic DNA extracted during a previously published study of the sediment at site 1227 (19). Others, e.g., the frequently recovered phylotype represented by clone 6H-2D22 within MCG-1, did not have closely related counterparts in rRNA gene clone libraries and are good candidates for highly active Archaea that are not dominant in DNA surveys (Fig. 3).
MBGB.Members of the MBGB were originally found in deep-sea sediments and at hydrothermal vent sites (49, 54), and the group is now detected in a growing number of benthic marine environments, including Atlantic deep-sea sediments, marine mud volcanoes, marine carbonate crusts from the Black Sea, organic-poor deep subsurface sediments of the central oceanic basins, and organic-rich methane- or methane-hydrate-containing sediments near continental margins (3, 26, 38, 45, 54). The group was recently reported to dominate clone libraries obtained from methane-hydrate-containing Peru Margin and Pacific Rim sediments, while being rare in clone libraries from nonhydrate sediments, including site 1227 (19). The RT-PCR results presented here demonstrate that the group may be abundant in certain layers of nonhydrate sediments.
SAGMEG.As the name implies, members of SAGMEG-1 and -2 were originally found in South African gold mines, but they have subsequently been detected numerous times in studies of subsurface and hydrothermal systems. Together with members of the MCG, SAGMEG-related phylotypes dominated clone libraries constructed from 16S rRNA genes from site 1227 (19). While less abundant in the rRNA libraries from the same site, the detection of the group demonstrates an in situ activity.
Activity of archaeal groups with depth in the sediment.The amount of rRNA and the RNA/DNA ratio in actively metabolizing cells may differ between species and strains, reflecting species- or strain-specific responses to environmental stimuli, as well as differences in rRNA gene copy number (31, 46). Furthermore, phylotype frequency in clone libraries of PCR products may be biased by selection during extraction and amplification (39, 47, 48). Consequently, the relative frequencies among different phylogenetic groups in the RT-PCR clone libraries presented here do not necessarily translate into relative levels of metabolic activity. On the other hand, vertical variations in the concentration of rRNA within a specific group or phylotype may provide a measure of its activity at different depths. The rRNA profiles permit evaluation of such vertical variations in activity of two of the detected phylogenetic groups, MCG and MBGB.
A downwardly decreasing richness in 16S rRNA phylotypes affiliated with the MCG group was observed in the clone libraries. The decreasing amount of template rRNA with depth may have diluted rare phylotypes to extinction, resulting in a lower diversity in the deep clone libraries. However, since only one phylotype was detected in sample 5H-5 from 40.35 mbsf, where the rRNA concentration was relatively high, the decreasing phylotype richness probably reflects a lower diversity in the active MCG community rather than a dilution effect. The MCG Archaea constituted more than 80% of all archaeal clone libraries from samples taken from throughout the upper 35 m, and the decreasing amount of rRNA with depth therefore strongly suggests a decreasing activity level of this group. MCG phylotypes were only sporadically detected at the top of the SMTZ, where the large amount of MBGB-related rRNA overshadowed other groups of organisms. Below the MBGB-dominated layer, two related MCG phylotypes accounted for the activity peak in the deepest part of the SMTZ. The phylotype represented by clone 6H-2L10 was detected in clone libraries from throughout the sediment column (Fig. 3), suggesting either that it is not directly involved in methane oxidation or that it is physiologically flexible.
Members of the MBGB had a striking affinity for a narrow sediment horizon at the center of the SMTZ, where they completely dominated the clone libraries. This strongly indicates an increased activity of the group in this specific layer. The increased activity of MBGB at 37.35 and 37.75 mbsf was not reflected in Q-PCR data obtained previously from the same sediment depth at site 1227 (42). This apparent discrepancy between RNA and DNA abundance could be a result of a higher rRNA/DNA ratio in the active MBGB than in the active MCG organisms.
The vertical rRNA profile indicates that the activity of two groups of Crenarchaeota is higher in the SMTZ than in the sediment above and below. High activity of Archaea in the SMTZ is consistent with another study of site 1227, in which Archaea were quantified by FISH counts and extraction of IPL. Of total cells in sediment samples from the SMTZ at site 1227, 98% hybridized with Archaea-specific oligonucleotide probes, and between 31 and 48% of extracted IPL in the SMTZ was of archaeal origin (3). Carbon isotopic data obtained from individual cells and from the IPL extracted from the SMTZ at site 1227 indicated that photosynthetically produced organic carbon, not methane, was the primary carbon source for the dominant archaeal community (3). If MBGB and/or MCG Archaea play any role in methane oxidation in the SMTZ, they do not assimilate methane-derived carbon and must be limited to a strictly dissimilatory methane metabolism (3). Both MCG and MBGB Archaea show a conspicuously cosmopolitan occurrence pattern in a wide spectrum of sediments (Table 1) and vents, indicating considerable ecological flexibility.
The correspondence between the vertical zonation of dissolved chemicals and active archaeal communities indicates a strong coupling between environmental chemistry and community composition in deep subseafloor sediments, with community activities and composition both responding to and mediating the geochemical gradients. Often, the microbial communities of subsurface SMTZs have been evaluated based on sediment from a single depth horizon within the SMTZ (38). The results presented here indicate that considerable stratification of the microbial community may occur within the SMTZ or similar geochemically well-defined sediment horizons. Therefore, detailed habitat studies should include multiple depth horizons within every specific, geochemically defined sediment layer in order to track unforeseen, steep gradients in microbial community composition and activity.
ACKNOWLEDGMENTS
This work was supported by the NASA Astrobiology Institute teams at the University of Rhode Island (“Subsurface Biospheres,” NCC 2-1275) and the Marine Biological Laboratory (“Environmental Genomes,” NCC 2-1054) and by the Danish Natural Science Research Council.
We thank all members of the leg 201 shipboard scientific party for help during the cruise, Barbara McGregor for help with the RNA extraction protocol, the staff of the Josephine Bay Paul Center for Comparative Molecular Biology and Evolution for timely clone library sequencing, and Steve D'Hondt for valuable comments and suggestions regarding the manuscript.
FOOTNOTES
- Received 8 March 2006.
- Accepted 17 April 2006.
- Copyright © 2006 American Society for Microbiology