skip to main content

Explore Research Products in the NSF-PAR

Title: Islands Within Islands: Bacterial Phylogenetic Structure and Consortia in Hawaiian Lava Caves and Fumaroles
Lava caves, tubes, and fumaroles in Hawai‘i present a range of volcanic, oligotrophic environments from different lava flows and host unexpectedly high levels of bacterial diversity. These features provide an opportunity to study the ecological drivers that structure bacterial community diversity and assemblies in volcanic ecosystems and compare the older, more stable environments of lava tubes, to the more variable and extreme conditions of younger, geothermally active caves and fumaroles. Using 16S rRNA amplicon-based sequencing methods, we investigated the phylogenetic distinctness and diversity and identified microbial interactions and consortia through co-occurrence networks in 70 samples from lava tubes, geothermal lava caves, and fumaroles on the island of Hawai‘i. Our data illustrate that lava caves and geothermal sites harbor unique microbial communities, with very little overlap between caves or sites. We also found that older lava tubes (500–800 yrs old) hosted greater phylogenetic diversity (Faith's PD) than sites that were either geothermally active or younger (<400 yrs old). Geothermally active sites had a greater number of interactions and complexity than lava tubes. Average phylogenetic distinctness, a measure of the phylogenetic relatedness of a community, was higher than would be expected if communities were structured at random. This suggests that bacterial communities of Hawaiian volcanic environments are phylogenetically over-dispersed and that competitive exclusion is the main driver in structuring these communities. This was supported by network analyses that found that taxa (Class level) co-occurred with more distantly related organisms than close relatives, particularly in geothermal sites. Network “hubs” (taxa of potentially higher ecological importance) were not the most abundant taxa in either geothermal sites or lava tubes and were identified as unknown families or genera of the phyla, Chloroflexi and Acidobacteria. These results highlight the need for further study on the ecological role of microbes in caves through targeted culturing methods, metagenomics, and long-read sequence technologies.

  more » « less
Award ID(s):
1711856
NSF-PAR ID:
10485778
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ;
Publisher / Repository:
Frontiers in Microbiology
Date Published:
Journal Name:
Frontiers in Microbiology
Volume:
13
ISSN:
1664-302X
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ; ; ( , ARPHA Conference Abstracts)
    Chi Fru, Ernest ; Chik, Alex ; Colwell, Fredrick ; Dittrich, Maria ; Engel, Annette ; Keenan, Sarah ; Meckenstock, Rainer ; Omelon, Christopher ; Purkamo, Lotta ; Weisener, Chris (Ed.)

    Roots are common features in basaltic lava tube caves on the island of Hawai‘i. For the past 50 years, new species of cave-adapted invertebrates, including cixiid planthoppers, crickets, thread-legged bugs, and spiders, have been discovered from root patches in lava tubes on different volcanoes and across variable climatic conditions. Assessing vegetation on the surface above lava tube passages, as well as genetic characterization of roots from within lava tubes, suggest that most roots belong to the native pioneer tree, ‘ōhi‘a lehua (Metrosideros polymorpha). Planthoppers are the primary consumers of sap at the base of the subsurface food web. However, root physicochemistry and rhizobiome microbial diversity and functional potential have received little attention. This study focuses on characterizing the ‘ōhi‘a rhizobiome, accessed from free-hanging roots inside lava tubes. Using these results, we can begin to evaluate the development and evolution of plant-microbe-invertebrate relationships.

    We explored lava tubes formed in flows of differing elevations and ages, from about 140 to 3000 years old, on Mauna Loa, Kīlauea, and Hualālai volcanoes on Hawai‘i Island. Invertebrate diversity was evaluated from root galleries and non-root galleries, in situ fluid physicochemistry was measured, and root and bare rock fluids (e.g., water, sap) were collected to determine major ion concentrations, as well as non-purgeable organic carbon (NPOC) and total nitrogen (TN) content. To verify root identity, DNA was extracted, and three sets of primers were used. After screening for onlyMetrosiderosspp., the V4 region of the 16S rRNA gene was sequenced and taxonomy was assigned.

    Root fluids were viscous and ranged in color from clear to yellow to reddish orange. Root fluids had 2X to 10X higher major ion concentrations compared to rock water. The average root NPOC and TN concentrations were 192 mg/L and 5.2 mg/L, respectively, compared to rock water that had concentrations of 6.8 mg/L and 1.8 mg/L, respectively. Fluids from almost 300 root samples had pH values that ranged from 2.2 to 5.6 (average pH 4.63) and were lower than rock water (average pH 6.39). Root fluid pH was comparable to soil pH from montane wet forests dominated by ‘ōhi‘a (Selmants et al. 2016), which can grow in infertile soil with pH values as low as 3.6. On Hawai‘i, rain water pH averages 5.2 at sea level and systematically decreases with elevation to pH 4.3 at 2500 m (Miller and Yoshinaga 2012), but root fluid pH did not correlate with elevation, temperature, relative humidity, inorganic and organic constituents, or age of flow. Root fluid acidity is likely due to concentrated organic compounds, sourced as root exudates, and this habitat is acidic for the associated invertebrates.

    From 62 root samples, over 66% were identified to the genusMetrosideros. A few other identifications of roots from lava tube systems where there had been extensive clear-cutting and ranching included monkey pod tree, coconut palm,Ficusspp., and silky oak.

    The 16S rRNA gene sequence surveys revealed that root bacterial communities were dominated by few groups, including Burkholderiaceae, as well as Acetobacteraceae, Sphingomonadaceae, Acidobacteriaceae, Gemmataceae, Xanthobacteraceae, and Chitinophagaceae. However, most of the reads could not be classified to a specific genus, which suggested that the rhizobiome harbor novel diversity. Diversity was higher from wetter climates. The root communities were distinct from those described previously from ‘ōhi‘a flowers and leaves (Junker and Keller 2015) and lava tube rocky surfaces (Hathaway et al. 2014) where microbial groups were specifically presumed capable of heterotrophy, methanotrophy, diazotrophy, and nitrification. Less can be inferred for the rhizobiome metabolism, although most taxa are likely aerobic heterotrophs. Within the Burkholderiaceae, there were high relative abundances of sequences affiliated with the genusParaburkholderia, which includes known plant symbionts, as well as the acidophilic generaAcidocellaandAcidisomafrom the Acetobacteraceae, which were retrieved predominately from caves in the oldest lava flows that also had the lowest root pH values. It is likely that the bacterial groups are capable of degrading exudates and providing nutritional substrates for invertebrate consumers that are not provided by root fluids (i.e., phloem) alone.

    As details about the biochemistry of ‘ōhi‘a have been missing, characterizing the rhizobiome from lava tubes will help to better understand potential plant-microbe-invertebrate interactions and ecological and evolutionary relationships through time. In particular, the microbial rhizobiome may produce compounds used by invertebrates nutritionally or that affect their behavior, and changes to the rhizobiome in response to environmental conditions may influence invertebrate interactions with the roots, which could be important to combat climate change effects or invasive species introductions.

     
    more » « less
  2. ; ; ; ; ; ; ; ; ; ; et al ( , Proceedings of the International Ocean Discovery Program Expedition reports)
    Hotspot tracks (chains of seamounts, ridges, and other volcanic structures) provide important records of plate motions, as well as mantle geodynamics, magma flux, and mantle source compositions. The Tristan-Gough-Walvis Ridge (TGW) hotspot track, extending from the active volcanic islands of Tristan da Cunha and Gough through a province of guyots and then along Walvis Ridge to the Etendeka flood basalt province, forms one of the most prominent and complex global hotspot tracks. The TGW hotspot track displays a tight linear age progression in which ages increase from the islands to the flood basalts (covering ~135 My). Unlike Pacific tracks, which are often simple, nearly linear chains of seamounts, the TGW track is alternately a steep-sided narrow ridge, an oceanic plateau, subparallel linear ridges and chains of seamounts (most are flat-topped guyots). The track displays isotopic zonation over the last ~70 My. The zonation appears near the middle of the track just before it splits into two to three chains of ridge- and guyot-type seamounts. Walvis Ridge, forming the older part of the track, is also overprinted with age-progressive late-stage volcanism, which was emplaced ~30–40 My after the initial eruptions and has a distinct isotopic composition. The plan for Expedition 391 was to drill at six sites, three along Walvis Ridge and three in the seamounts of the Guyot Province, to collect igneous rocks to better understand the formation of volcanic edifices, the temporal and geochemical evolution of the hotspot, and the variation in paleolatitudes at which the volcanic edifices formed. After a delay of 18 days to address a shipboard Coronavirus (COVID-19) outbreak, Expedition 391 proceeded to drill at four of the proposed sites: three sites on Walvis Ridge around Valdivia Bank, an ocean plateau within the ridge, and one site on the lower flank of a guyot in the Center track of the Guyot Province, a ridge located between the Tristan subtrack (which extends from the end of Walvis Ridge to the islands of Tristan da Cunha) and the Gough subtrack (which extends from Walvis Ridge to Gough Island). The first hole was drilled at Site U1575, located on a low portion of the northeastern Walvis Ridge just north of Valdivia Bank. At this location, 209.9 m of sediments and 122.4 m of igneous basement were cored. The sediments ranged in age from Late Pleistocene (~0.43–1.24 Ma) to Late Cretaceous (Campanian; 72–78 Ma). The igneous basement comprised 10 submarine lava units consisting of pillow, lobate, sheet, and massive lava flows, the thickest of which was ~21 m. Most lavas are tholeiitic, but some alkalic basalts were recovered. A portion of the igneous succession consists of low-Ti basalts, which are unusual because they appear in the Etendeka flood basalts but have not been previously found on Walvis Ridge. Two holes were drilled at Site U1576 on the west flank of Valdivia Bank. The first of these holes was terminated because a bit jammed shortly after entering the igneous basement. Hole U1576A recovered a remarkable ~380 m thick sedimentary section consisting mostly of chalk covering a nearly complete sequence from Late Pleistocene (~0.43–1.24 Ma) to Late Cretaceous (Campanian; ~79–81.38 Ma). These sediments display short and long cyclic color changes that imply astronomically forced and longer term paleoenvironmental changes. The igneous basement recovered in Hole U1576B yielded 11 submarine lava units (total thickness = ~65 m). The flows range from pillows to massive flows with compositions varying from tholeiitic basalt to basaltic andesite, only the second occurrence of the latter composition recovered from the TGW track thus far. These units are separated by seven sedimentary chalk units that range 0.1–11.6 m in thickness, implying a long-term interplay of sedimentation and lava eruptions. These intercalated sediments revealed Upper Cretaceous (Campanian) ages of ~77–79 Ma for the upper two interbeds and ~79–81.38 Ma for the lower beds. Coring at Site U1577, on the extreme eastern flank of Valdivia Bank, penetrated a 154.8 m thick sedimentary section ranging from the Paleocene (Thanetian; ~58.8 Ma) to Upper Cretaceous (Campanian; ~81.43–83.20 Ma). Igneous basement coring progressed only 39.1 m below the sediment/basalt contact, recovering three massive submarine tholeiitic basalt lava flows that are 4.1, 15.5, and >19.1 m thick, respectively. Paleomagnetic data from Sites U1577 and U1576 indicate that the former volcanic basement formed just before the end of the Cretaceous Normal Superchron and the latter during Chron 33r, shortly afterward. Biostratigraphic and paleomagnetic data suggest that Valdivia Bank becomes younger from east to west. Site U1578, located on a Center track guyot, provided a long and varied igneous section. After coring through 184.3 m of pelagic carbonate sediments mainly consisting of Eocene and Paleocene chalk (~55.64–63.5 Ma), Hole U1578A cored 302.1 m of igneous basement. Basement lavas are largely pillows but are interspersed with sheet and massive flows. Lava compositions are mostly alkalic basalts with some hawaiite. Several intervals contain abundant olivine (some fresh), and some of the pillow stacks consist of basalt with remarkably high Ti content. The igneous sequence is interrupted by 10 sedimentary interbeds consisting of chalk and volcaniclastics and ranging 0.46–10.19 m in thickness. Investigations of toothpick samples from the intercalated sediments were examined, each revealing the same age range of ~63.5–64.81 Ma (lower Paleocene; Danian). Paleomagnetic data display a change in basement magnetic polarity ~100 m above the base of the hole. Combining magnetic stratigraphy with biostratigraphic data, the igneous section is inferred to span >1 My. Nearly 7 months after Expedition 391, JOIDES Resolution transited from Cape Town to the north Atlantic. During this transit (Expedition 397T), 7.9 days of ship time were used to drill two holes (U1584A and U1585A) at sites on the Gough and Tristan tracks that had been omitted because of COVID-19–related time loss on the earlier cruise. For both, coring was begun only a short distance above the igneous basement to save time. The 75.2 m thick section drilled in Hole U1584A contains two sedimentary units: clay-rich carbonate sediments overlie a pumice-dominated volcaniclastic deposit containing basalt fragments. Because the goal was to core basalt and the base of the volcaniclastic deposit was not imaged in the seismic profile, the hole was terminated early to save operation time for the next site. In Hole U1585A, coring penetrated a 273.5 m thick sediment section overlying an 81.2 m thick pile of massive basalt flows. The sediment section is divided into four units: The uppermost unit consists of nannofossil chalk; The two intermediate units contain alternating chalk and volcaniclastic sediments containing several breccia units; and The lowermost unit consists of volcanic breccia containing juvenile blocks, bombs, and accretionary lapilli. This thick sedimentary section documents a transition from shallow-water volcanism to open-ocean sedimentation as the seamount subsided. The thick underlying basalt section is made up of four sparsely to highly phyric massive flows, the thickest of which is >43 m thick. Samples of these units are mostly basalt with a few trachybasalts and one trachyandesite. Although the igneous penetration was less than planned, coring during Expeditions 391 and 397T obtained samples that clearly will lead to an improved understanding of the evolution of the TGW hotspot and its track. Reasonable recovery of fresh basalt in some holes provides ample samples for geochemical, geochronologic, and paleomagnetic studies. Good recovery of Late Cretaceous and early Cenozoic chalk successions provides samples for paleoenvironmental study. 
    more » « less
  3. ; ; ( , Preliminary report)
    Hotspot tracks (quasilinear chains of seamounts, ridges, and other volcanic structures) provide important records of plate motions, as well as mantle geodynamics, magma flux, and mantle source compositions. The Tristan-Gough-Walvis Ridge (TGW) hotspot track, extending from the active volcanic islands of Tristan da Cunha and Gough through a province of guyots and then along Walvis Ridge to the Etendeka flood basalt province, forms one of the most prominent and complex global hotspot tracks. The TGW hotspot track displays a tight linear age progression in which ages increase from the islands to the flood basalts (covering ~135 My). Unlike Pacific tracks, which are simple chains of seamounts that are often compared to chains of pearls, the TGW track is alternately a steep-sided narrow ridge, an oceanic plateau, subparallel linear ridges and chains of seamounts, and areas of what appear to be randomly dispersed seamounts. The track displays isotopic zonation over the last ~70 My. The zonation appears near the middle of the track just before it splits into two to three chains of ridge- and guyot-type seamounts. The older ridge is also overprinted with age-progressive late-stage volcanism, which was emplaced ~30–40 My after the initial eruptions and has a distinct isotopic composition. The plan for Expedition 391 was to drill at six sites, three along Walvis Ridge and three in the seamount (guyot) province, to gather igneous rocks to better understand the formation of track edifices, the temporal and geochemical evolution of the hotspot, and the variation in paleolatitudes at which the volcanic edifices formed. After a delay of 18 days to address a shipboard outbreak of the coronavirus disease 2019 (COVID-19) virus, Expedition 391 proceeded to drill at four of the proposed sites: three sites on the eastern Walvis Ridge around Valdivia Bank, an ocean plateau within the ridge, and one site on the lower flank of a guyot in the Center track, a ridge located between the Tristan subtrack (which extends from the end of Walvis Ridge to the island of Tristan da Cunha) and the Gough subtrack (which extends from Walvis Ridge to the island of Gough). One hole was drilled at Site U1575, located on a low portion of the northeastern Walvis Ridge north of Valdivia Bank. At this location, 209.9 m of sediments and 122.4 m of igneous basement were cored. The latter comprised 10 submarine lava units consisting of pillow, lobate, sheet, and massive lava flows, the thickest of which was ~21 m. Most lavas are tholeiitic, but some alkalic basalts were recovered. A portion of the igneous succession consists of low-Ti basalts, which are unusual because they appear in the Etendeka flood basalts but have not been previously found on Walvis Ridge. Two holes were drilled at Site U1576 on the west flank of Valdivia Bank. The first hole was terminated because a bit jammed shortly after penetrating igneous basement. Hole U1576A recovered a remarkable ~380 m thick sedimentary section consisting mostly of chalk covering a nearly complete sequence from Paleocene to Late Cretaceous (Campanian). These sediments display short and long cyclic color changes that imply astronomically forced and longer term paleoenvironmental changes. The igneous basement yielded 11 submarine lava units ranging from pillows to massive flows, which have compositions varying from tholeiitic basalt to basaltic andesite, the first occurrence of this composition recovered from the TGW track. These units are separated by seven sedimentary chalk units that range in thickness from 0.1 to 11.6 m, implying a long-term interplay of sedimentation and lava eruptions. Coring at Site U1577, on the extreme eastern flank of Valdivia Bank, penetrated a 154 m thick sedimentary section, the bottom ~108 m of which is Maastrichtian–Campanian (possibly Santonian) chalk with vitric tephra layers. Igneous basement coring progressed only 39.1 m below the sediment-basalt contact, recovering three massive submarine tholeiite basalt lava flows that are 4.1, 15.5, and >19.1 m thick, respectively. Paleomagnetic data from Sites U1577 and U1576 indicate that their volcanic basements formed just before the end of the Cretaceous Normal Superchron and during Chron 33r, shortly afterward, respectively. Biostratigraphic and paleomagnetic data suggest an east–west age progression across Valdivia Bank, becoming younger westward. Site U1578, located on a Center track guyot, provided a long and varied igneous section. After coring through 184.3 m of pelagic carbonate sediments mainly consisting of Eocene and Paleocene chalk, Hole U1578A cored 302.1 m of igneous basement. Basement lavas are largely pillows but are interspersed with sheet and massive flows. Lava compositions are mostly alkalic basalts with some hawaiite. Several intervals contain abundant olivine, and some of the pillow stacks consist of basalt with remarkably high Ti content. The igneous sequence is interrupted by 10 sedimentary interbeds consisting of chalk and volcaniclastics and ranging in thickness from 0.46 to 10.19 m. Paleomagnetic data display a change in basement magnetic polarity ~100 m above the base of the hole. Combining magnetic stratigraphy with biostratigraphic data, the igneous section is inferred to span >1 My. Abundant glass from pillow lava margins was recovered at Sites U1575, U1576, and U1578. Although the igneous penetration was only two-thirds of the planned amount, drilling during Expedition 391 obtained samples that clearly will lead to a deeper understanding of the evolution of the Tristan-Gough hotspot and its track. Relatively fresh basalts with good recovery will provide ample samples for geochemical, geochronologic, and paleomagnetic studies. Good recovery of Late Cretaceous and early Cenozoic chalk successions provides samples for paleoenvironmental study. 
    more » « less
  4. ; ; ; ; ; ; ; ; ; ; et al ( , Journal of Geophysical Research: Planets)
    Abstract

    Craters of the Moon National Monument and Preserve (CotM) is a strong terrestrial analog for lava tube formations on Mars. The commonality of its basalt composition to Martian lava tubes makes it especially useful for probing how interactions between water, rock, and life have developed over time, and what traces of these microbial communities may be detectable by current flight‐capable instrumentation. Our investigations found that secondary mineral deposits within these caves contain a range of underlying compositions that support diverse and active microbial communities. Examining the taxonomy, activity, and metabolic potential of these communities revealed largely heterotrophic life strategies supported by contributions from chemolithoautotrophs that facilitate key elemental cycles. Finally, traces of these microbial communities were detectable by flight‐capable pyrolysis and wet chemistry gas chromatography‐mass spectrometry methods comparable to those employed by the Sample Analysis at Mars (SAM) instrument aboard the Curiosity rover and the Mars Organic Molecule Analyzer (MOMA) on the upcoming Rosalind Franklin rover. Using a suite of methods for chemical derivatization of organic compounds is beneficial for resolving the greatest variety of biosignatures. Tetramethylammonium hydroxide (TMAH), for example, allowed for optimal resolution of long chain fatty acids. Taken together, these results have implications for the direction of mass spectrometry as a tool for biosignature detection on Mars, as well as informing the selection of sampling sites that could potentially host biosignatures.

     
    more » « less
  5. ; ; ; ; ( , Ecosphere)
    Abstract

    Habitat fragmentation resulting in habitat loss and increased isolation is a dominant driver of global species declines. Habitat isolation and connectivity vary across scales, and understanding how connectivity affects biodiversity can be challenging because the relevant scale depends on the taxa involved. A multiscale analysis can provide insight in biodiversity patterns across spatial scale when information on dispersal ability is not available, in particular for community‐level studies focusing on multiple taxa. In this study, we examine the relationship between arthropod diversity, patch area, and connectivity using a multiscale approach. We make use of a natural experiment on Hawai‘i Island, where historic volcanic activity has transformed contiguous native forests to lava matrix and discrete forest patches. This landscape of patches has persisted for 150 yr, and we selected 10,000 ha consisting of 863 patches to analyze landscape connectivity using a graph theory approach. We collected arthropod samples fromMetrosideros polymorpha tree canopies in 34 forest patches during multiple years. We analyzed the relationship of arthropod diversity with area, as well as with connectivity across increasing scales, or dispersal threshold distances. In contrast to well‐established ecological theory as well as prior work on birds and fungi in this system, we did not find support for a canonical species–area relationship. Next, we calculated connectivity across spatial scales and found lower Shannon diversity with higher connectivity at small scales, but no effect at increased dispersal threshold distances. We examined the landscape structure and found all habitat patches connected into three subnetworks at a 350 m threshold distance. All patches were connected at 700 m threshold distance, indicating structural dispersal limitation only at small scales. Our findings suggest that canopy arthropods are not dispersal limited at scales shown to impact both soil fungi and birds in this system. Instead, Hawaiian canopy arthropods may perceive the landscape as a connected area where discrete forest patches and the early‐successional matrix contribute resources that vary spatially with regard to habitat quality. We argue for the utility of multiscale approaches, and the importance of examining maintenance of biodiversity in fragmented landscapes that persist for hundreds of years.

     
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email