Evidence for an ancient aquatic origin of the RNA viral order Articulavirales.

The emergence of previously unknown disease-causing viruses in mammals is in part the result of a long-term evolutionary process. Reconstructing the deep phylogenetic histories of viruses helps identify major evolutionary transitions and contextualizes the emergence of viruses in new hosts. We used a combination of total RNA sequencing and transcriptome data mining to extend the diversity and evolutionary history of the RNA virus order Articulavirales, which includes the influenza viruses. We identified instances of Articulavirales in the invertebrate phylum Cnidaria (including corals), constituting a novel and divergent family that we provisionally named the "Cnidenomoviridae." We further extended the evolutionary history of the influenza virus lineage by identifying four divergent, fish-associated influenza-like viruses, thereby supporting the hypothesis that fish were among the first hosts of influenza viruses. In addition, we substantially expanded the phylogenetic diversity of quaranjaviruses and proposed that this genus be reclassified as a family-the "Quaranjaviridae." Within this putative family, we identified a novel arachnid-infecting genus, provisionally named "Cheliceravirus." Notably, we observed a close phylogenetic relationship between the Crustacea- and Chelicerata-infecting "Quaranjaviridae" that is inconsistent with virus-host codivergence. Together, these data suggest that the Articulavirales has evolved over at least 600 million years, first emerging in aquatic animals. Importantly, the evolution of the Articulavirales was likely shaped by multiple aquatic-terrestrial transitions and substantial host jumps, some of which are still observable today.

Groundwater is a hidden global keystone ecosystem.

Groundwater is a vital ecosystem of the global water cycle, hosting unique biodiversity and providing essential services to societies. Despite being the largest unfrozen freshwater resource, in a period of depletion by extraction and pollution, groundwater environments have been repeatedly overlooked in global biodiversity conservation agendas. Disregarding the importance of groundwater as an ecosystem ignores its critical role in preserving surface biomes. To foster timely global conservation of groundwater, we propose elevating the concept of keystone species into the realm of ecosystems, claiming groundwater as a keystone ecosystem that influences the integrity of many dependent ecosystems. Our global analysis shows that over half of land surface areas (52.6%) has a medium-to-high interaction with groundwater, reaching up to 74.9% when deserts and high mountains are excluded. We postulate that the intrinsic transboundary features of groundwater are critical for shifting perspectives towards more holistic approaches in aquatic ecology and beyond. Furthermore, we propose eight key themes to develop a science-policy integrated groundwater conservation agenda. Given ecosystems above and below the ground intersect at many levels, considering groundwater as an essential component of planetary health is pivotal to reduce biodiversity loss and buffer against climate change.

Evolutionary responses of a reef-building coral to climate change at the end of the last glacial maximum.

Climate change threatens the survival of coral reefs on a global scale, primarily through mass bleaching and mortality as a result of marine heatwaves. While these short-term effects are clear, predicting the fate of coral reefs over the coming century is a major challenge. One way to understand the longer-term effects of rapid climate change is to examine the response of coral populations to past climate shifts. Coastal and shallow-water marine ecosystems such as coral reefs have been reshaped many times by sea-level changes during the Pleistocene, yet, few studies have directly linked this with its consequences on population demographics, dispersal, and adaptation. Here we use powerful analytical techniques, afforded by haplotype phased whole-genomes, to establish such links for the reef-building coral, Acropora digitifera. We show that three genetically distinct populations are present in northwestern Australia, and that their rapid divergence since the last glacial maximum (LGM) can be explained by a combination of founder-effects and restricted gene flow. Signatures of selective sweeps, too strong to be explained by demographic history, are present in all three populations and overlap with genes that show different patterns of functional enrichment between inshore and offshore habitats. In contrast to rapid divergence in the host, we find that photosymbiont communities are largely undifferentiated between corals from all three locations, spanning almost 1000 km, indicating that selection on host genes and not acquisition of novel symbionts, has been the primary driver of adaptation for this species in northwestern Australia.

A hybrid-capture approach to reconstruct the phylogeny of Scleractinia (Cnidaria: Hexacorallia).

A well-supported evolutionary tree representing most major lineages of scleractinian corals is in sight with the development and application of phylogenomic approaches. Specifically, hybrid-capture techniques are shedding light on the evolution and systematics of corals. Here, we reconstructed a broad phylogeny of Scleractinia to test previous phylogenetic hypotheses inferred from a few molecular markers, in particular, the relationships among major scleractinian families and genera, and to identify clades that require further research. We analysed 449 nuclear loci from 422 corals, comprising 266 species spanning 26 families, combining data across whole genomes, transcriptomes, hybrid capture and low-coverage sequencing to reconstruct the largest phylogenomic tree of scleractinians to date. Due to the large number of loci and data completeness (less than 38% missing data), node supports were high across shallow and deep nodes with incongruences observed in only a few shallow nodes. The "Robust" and "Complex" clades were recovered unequivocally, and our analyses confirmed that Micrabaciidae Vaughan, 1905 is sister to the "Robust" clade, transforming our understanding of the "Basal" clade. Several families remain polyphyletic in our phylogeny, including Deltocyathiidae Kitahara, Cairns, Stolarski & Miller, 2012, Caryophylliidae Dana, 1846, and Coscinaraeidae Benzoni, Arrigoni, Stefani & Stolarski, 2012, and we hereby formally proposed the family name Pachyseridae Benzoni & Hoeksema to accommodate Pachyseris Milne Edwards & Haime, 1849, which is phylogenetically distinct from Agariciidae Gray, 1847. Results also revealed species misidentifications and inconsistencies within morphologically complex clades, such as Acropora Oken, 1815 and Platygyra Ehrenberg, 1834, underscoring the need for reference skeletal material and topotypes, as well as the importance of detailed taxonomic work. The approach and findings here provide much promise for further stabilising the topology of the scleractinian tree of life and advancing our understanding of coral evolution.

Macroalgae exhibit diverse responses to human disturbances on coral reefs.

Scientists and managers rely on indicator taxa such as coral and macroalgal cover to evaluate the effects of human disturbance on coral reefs, often assuming a universally positive relationship between local human disturbance and macroalgae. Despite evidence that macroalgae respond to local stressors in diverse ways, there have been few efforts to evaluate relationships between specific macroalgae taxa and local human-driven disturbance. Using genus-level monitoring data from 1205 sites in the Indian and Pacific Oceans, we assess whether macroalgae percent cover correlates with local human disturbance while accounting for factors that could obscure or confound relationships. Assessing macroalgae at genus level revealed that no genera were positively correlated with all human disturbance metrics. Instead, we found relationships between the division or genera of algae and specific human disturbances that were not detectable when pooling taxa into a single functional category, which is common to many analyses. The convention to use percent cover of macroalgae as an indication of local human disturbance therefore likely obscures signatures of local anthropogenic threats to reefs. Our limited understanding of relationships between human disturbance, macroalgae taxa, and their responses to human disturbances impedes the ability to diagnose and respond appropriately to these threats.

Population connectivity and genetic offset in the spawning coral Acropora digitifera in Western Australia.

Anthropogenic climate change has caused widespread loss of species biodiversity and ecosystem productivity across the globe, particularly on tropical coral reefs. Predicting the future vulnerability of reef-building corals, the foundation species of coral reef ecosystems, is crucial for cost-effective conservation planning in the Anthropocene. In this study, we combine regional population genetic connectivity and seascape analyses to explore patterns of genetic offset (the mismatch of gene-environmental associations under future climate conditions) in Acropora digitifera across 12 degrees of latitude in Western Australia. Our data revealed a pattern of restricted gene flow and limited genetic connectivity among geographically distant reef systems. Environmental association analyses identified a suite of loci strongly associated with the regional temperature variation. These loci helped forecast future genetic offset in gradient forest and generalized dissimilarity models. These analyses predicted pronounced differences in the response of different reef systems in Western Australia to rising temperatures. Under the most optimistic future warming scenario (RCP 2.6), we predicted a general pattern of increasing genetic offset with latitude. Under the extreme climate scenario (RCP 8.5 in 2090-2100), coral populations at the Ningaloo World Heritage Area were predicted to experience a higher mismatch between current allele frequencies and those required to cope with local environmental change, compared to populations in the inshore Kimberley region. The study suggests complex and spatially heterogeneous patterns of climate-change vulnerability in coral populations across Western Australia, reinforcing the notion that regionally tailored conservation efforts will be most effective at managing coral reef resilience into the future.

Phylogeography of recent Plesiastrea (Scleractinia: Plesiastreidae) based on an integrated taxonomic approach.

Scleractinian corals are a diverse group of ecologically important yet highly threatened marine invertebrates, which can be challenging to identify to the species level. An influx of molecular studies has transformed scleractinian systematics, highlighting that cryptic species may be more common than previously understood. In this study, we test the hypothesis that Plesiastrea versipora (Lamarck, 1816), a species currently considered to occur throughout the Indo-Pacific in tropical, sub-tropical and temperate waters, is a single species. Molecular and morphological analyses were conducted on 80 samples collected from 31 sites spanning the majority of the species putative range and twelve mitogenomes were assembled to identify informative regions for phylogenetic reconstruction. Congruent genetic data across three gene regions supports the existence of two monophyletic clades aligning with distinct tropical and temperate provenances. Multivariate macromorphological analyses based on 13 corallite characters provided additional support for the phylogeographic split, with the number of septa and corallite density varying across this biogeographic divide. Furthermore, micromorphological and microstructural analyses identified that the temperate representatives typically develop sub-cerioid corallites with sparse or absent coenosteal features and smooth septal faces. In contrast, tropical representatives typically develop plocoid corallites separated by a porous dissepimental coenosteum and have granulated septal faces. These data suggest that at least two species exist within the genus PlesiastreaMilne Edwards & Haime, 1848. Based on examination of type material, we retain the name Plesiastrea versipora (Lamarck, 1816) for the temperate representatives of the genus and resurrect the name Plesiastrea peroniMilne Edwards & Haime, 1857 for the tropical members. This study highlights how broadly distributed hard coral taxa still need careful re-examination through an integrated systematics approach to better understand their phylogeographic patterns. Furthermore, it demonstrates the utility of integrating micro-, macro-morphological and genetic datasets, and the importance of type specimens when dealing with taxonomic revisions of scleractinian taxa.

eDNA metabarcoding survey reveals fine-scale coral reef community variation across a remote, tropical island ecosystem.

Environmental DNA (eDNA) metabarcoding, a technique for retrieving multispecies DNA from environmental samples, can detect a diverse array of marine species from filtered seawater samples. There is a growing potential to integrate eDNA alongside existing monitoring methods in order to establish or improve the assessment of species diversity. Remote island reefs are increasingly vulnerable to climate-related threats and as such there is a pressing need for cost-effective whole-ecosystem surveying to baseline biodiversity, study assemblage changes and ultimately develop sustainable management plans. We investigated the utility of eDNA metabarcoding as a high-resolution, multitrophic biomonitoring tool at the Cocos (Keeling) Islands, Australia (CKI)-a remote tropical coral reef atoll situated within the eastern Indian Ocean. Metabarcoding assays targeting the mitochondrial 16S rRNA and CO1 genes, as well as the 18S rRNA nuclear gene, were applied to 252 surface seawater samples collected from 42 sites within a 140 km2 area. Our assays successfully detected a wide range of bony fish and elasmobranchs (244 taxa), crustaceans (88), molluscs (37) and echinoderms (7). Assemblage composition varied significantly between sites, reflecting habitat partitioning across the island ecosystem and demonstrating the localisation of eDNA signals, despite extensive tidal and oceanic movements. In addition, we document putative new occurrence records for 46 taxa and compare the efficiency of our eDNA approach to visual survey techniques at CKI. Our study demonstrates the utility of a multimarker metabarcoding approach in capturing multitrophic biodiversity across an entire coral reef atoll and sets an important baseline for ongoing monitoring and management.

Evolutionary biogeography of the reef-building coral genus Galaxea across the Indo-Pacific ocean.

Stony corals (Scleractinia) form the basis for some of the most diverse ecosytems on Earth, but we have much to learn about their evolutionary history and systematic relationships. In order to improve our understanding of species in corals we here investigated phylogenetic relationships between morphologically defined species and genetic lineages in the genus Galaxea (Euphyllidae) using a combined phylogenomic and phylogeographic approach. Previous studies revealed the nominal species G. fascicularis included three genetically well-differentiated lineages (L, S & L+) in the western Pacific, but their distribution and relationship to other species in the genus was unknown. Based on genomic (RAD-seq) and mitochondrial sequence data (non-coding region between cytb and ND2) we investigated whether the morphological taxa represent genetically coherent entities and what is the phylogenetic relationship and spatial distribution of the three lineages of G. fascicularis throughout the observed species range. Using the RAD-seq data, we find that the genus Galaxea is monophyletic and contains three distinct clades: an Indo-Pacific, a Pacific, and a small clade restricted to the Chagos Archipelago. The three lineages of G. fascicularis were associated with different RAD-seq clades, with the 'L' lineage showing some morphological distinction from the other two lineages (larger more asymmetrical polyps). In addition to these, three more genetic lineages in G. fascicularis may be distinguished - a Chagossian, an Ogasawaran, and one from the Indian-Red Sea. Among nominal taxa for which we have multiple samples, G. horrescens was the only monophyletic species. The mitochondrial non-coding region is highly conserved apart of the length polymorphism used to define L, S & L+ lineages and lacks the power to distinguish morphological and genetic groups resolved with genomic RAD-sequencing. The polyphyletic nature of most species warrants a careful examination of the accepted taxonomy of this group with voucher collections and their comparison to type specimens to resolve species boundaries. Further insight to the speciation process in corals will require international cooperation for the sharing of specimens to facilitate scientific discovery.

Genetic signatures through space, time and multiple disturbances in a ubiquitous brooding coral.

The predominance of self-recruitment in many reef-building corals has fundamental and complex consequences for their genetic diversity, population persistence and responses to climate change. Knowledge of genetic structure over local scales needs to be placed within a broad spatial context, and also integrated with genetic monitoring through time to disentangle these consequences. Here, we examined patterns of genetic diversity over multiple spatio-temporal scales across tropical Australia in the ubiquitous brooding coral, Seriatopora hystrix. We also analysed complimentary environmental and demographic data to elucidate the seascape drivers of these patterns. Large genetic differences were detected between the east vs. west coasts of Australia. In northwest Australia, geographic differentiation dominated genetic structure over multiple scales. However, three sympatric lineages were detected at the largest offshore reef system (Scott Reef). Similar to the differences observed among putative species in eastern Australia, these lineages were associated with different levels of wave exposure. Local genetic structure within the Scott Reef system was relatively stable over 10 years, but temporal differences were observed that reflected small but important genetic changes over a few generations during recovery after severe bleaching. These results highlight the importance of self-recruitment together with occasional longer distance connectivity for the persistence of a metapopulation across spatially and temporally variable environments. Our multidimensional research provides a foundation for further long-term genetic monitoring to inform conservation strategies and highlights that sampling scales, ecological effects and cryptic diversity are important considerations to develop realistic understanding of the evolutionary resilience of corals.

Phylogenomics provides new insight into evolutionary relationships and genealogical discordance in the reef-building coral genus Acropora.

Understanding the genetic basis of reproductive isolation is a long-standing goal of speciation research. In recently diverged populations, genealogical discordance may reveal genes and genomic regions that contribute to the speciation process. Previous work has shown that conspecific colonies of Acropora that spawn in different seasons (spring and autumn) are associated with highly diverged lineages of the phylogenetic marker PaxC Here, we used 10 034 single-nucleotide polymorphisms to generate a genome-wide phylogeny and compared it with gene genealogies from the PaxC intron and the mtDNA Control Region in 20 species of Acropora, including three species with spring- and autumn-spawning cohorts. The PaxC phylogeny separated conspecific autumn and spring spawners into different genetic clusters in all three species; however, this pattern was not supported in two of the three species at the genome level, suggesting a selective connection between PaxC and reproductive timing in Acropora corals. This genome-wide phylogeny provides an improved foundation for resolving phylogenetic relationships in Acropora and, combined with PaxC, provides a fascinating platform for future research into regions of the genome that influence reproductive isolation and speciation in corals.

Are Big Food's corporate social responsibility strategies valuable to communities? A qualitative study with parents and children.

OBJECTIVE: Recent studies have identified parents and children as two target groups whom Big Food hopes to positively influence through its corporate social responsibility (CSR) strategies. The current preliminary study aimed to gain an in-depth understanding of parents and children's awareness and interpretation of Big Food's CSR strategies to understand how CSR shapes their beliefs about companies. DESIGN: Community-based qualitative semi-structured interviews. SETTING: New South Wales, Australia. SUBJECTS: Parents (n 15) and children aged 8-12 years (n 15). RESULTS: Parents and children showed unprompted recognition of CSR activities when shown McDonald's and Coca-Cola brand logos, indicating a strong level of association between the brands and activities that target the settings of children. When discussing CSR strategies some parents and most children saw value in the activities, viewing them as acts of merit or worth. For some parents and children, the companies' CSR activities were seen as a reflection of the company's moral attributes, which resonated with their own values of charity and health. For others, CSR strategies were in conflict with companies' core business. Finally, some also viewed the activities as harmful, representing a deceit of the public and a smokescreen for the companies' ultimately unethical behaviour. CONCLUSIONS: A large proportion of participants valued the CSR activities, signalling that denormalising CSR to sever the strong ties between the community and Big Food will be a difficult process for the public health community. Efforts to gain public acceptance for action on CSR may need greater levels of persuasion to gain public support of a comprehensive and restrictive approach.

Assessing implementation fidelity and adaptation in a community-based childhood obesity prevention intervention.

Little research has assessed the fidelity, adaptation or integrity of activities implemented within community-based obesity prevention initiatives. To address this gap, a mixed-method process evaluation was undertaken in the context of the South Australian Obesity Prevention and Lifestyle (OPAL) initiative. An ecological coding procedure assessed fidelity and adaptation of activity settings, targets and strategies implemented in the second year of four communities. Implementation integrity reflected fidelity and adaptation to local context, whereas efforts resulting in significant deviations from the original plan were deemed to lack fidelity and integrity. Staff implemented 284 strategies in 205 projects. Results show that 68.3 and 2.1% of strategies were implemented with fidelity or adapted, respectively. Overall, 70.4% of all strategies were implemented with integrity. Staff experienced barriers with 29.6% of strategies. Chi-square analyses show statistically significant associations between implementation integrity and strategy type, intervention and behavioural targets. These relationships are weak to modest. The strongest relationship was found between implementation integrity and proximal target. Staff experienced implementation barriers at the coalition, policy, organization, interpersonal and community levels. The greatest range of barriers was encountered working with organizations. To overcome these barriers, staff took greater ownership, invested more time, persisted and allocated more financial resources.

The complete mitochondrial genome of the invasive cyanobacteriosponge Terpios hoshinota (Demospongiae, Suberitida, Suberitidae).

The cyanobacteriosponge Terpios hoshinota occurs on tropical reefs throughout the Indo-Pacific. The species encrusts live coral, and other benthos, and is considered a pest species that can threaten the health and productivity of locally native benthic communities on coral reefs. Here we assemble a complete mitochondrial genome to aid further research into the range expansion of this species. The circular genome was 20,504 bp in length and encoded 14 protein-coding genes, two ribosomal RNA (rRNA) genes, and 25 transfer RNA (tRNA) genes. A phylogenetic analysis based on the concatenated sequences of 14 protein-coding genes of 12 members of the subclass Heteroscleromorpha including the newly sequenced T. hoshinota, suggests further taxonomic revisions within the order Suberitida may be warranted.

Aquatic environmental DNA: A review of the macro-organismal biomonitoring revolution.

Environmental DNA (eDNA) is the fastest growing biomonitoring tool fuelled by two key features: time efficiency and sensitivity. Technological advancements allow rapid biodiversity detection at both species and community levels with increasing accuracy. Concurrently, there has been a global demand to standardise eDNA methods, but this is only possible with an in-depth overview of the technological advancements and a discussion of the pros and cons of available methods. We therefore conducted a systematic literature review of 407 peer-reviewed papers on aquatic eDNA published between 2012 and 2021. We observed a gradual increase in the annual number of publications from four (2012) to 28 (2018), followed by a rapid growth to 124 publications in 2021. This was mirrored by a tremendous diversification of methods in all aspects of the eDNA workflow. For example, in 2012 only freezing was applied to preserve filter samples, whereas we recorded 12 different preservation methods in the 2021 literature. Despite an ongoing standardisation debate in the eDNA community, the field is seemingly moving fast in the opposite direction and we discuss the reasons and implications. Moreover, by compiling the largest PCR-primer database to date, we provide information on 522 and 141 published species-specific and metabarcoding primers targeting a wide range of aquatic organisms. This works as a user-friendly 'distillation' of primer information that was hitherto scattered across hundreds of papers, but the list also reflects which taxa are commonly studied with eDNA technology in aquatic environments such as fish and amphibians, and reveals that groups such as corals, plankton and algae are under-studied. Efforts to improve sampling and extraction methods, primer specificity and reference databases are crucial to capture these ecologically important taxa in future eDNA biomonitoring surveys. In a rapidly diversifying field, this review synthetises aquatic eDNA procedures and can guide eDNA users towards best practice.

Spatially varying selection between habitats drives physiological shifts and local adaptation in a broadcast spawning coral on a remote atoll in Western Australia.

At the Rowley Shoals in Western Australia, the prominent reef flat becomes exposed on low tide and the stagnant water in the shallow atoll lagoons heats up, creating a natural laboratory for characterizing the mechanisms of coral resilience to climate change. To explore these mechanisms in the reef coral Acropora tenuis, we collected samples from lagoon and reef slope habitats and combined whole-genome sequencing, ITS2 metabarcoding, experimental heat stress, and transcriptomics. Despite high gene flow across the atoll, we identified clear shifts in allele frequencies between habitats at relatively small linked genomic islands. Common garden heat stress assays showed corals from the lagoon to be more resistant to bleaching, and RNA sequencing revealed marked differences in baseline levels of gene expression between habitats. Our results provide new insight into the complex mechanisms of coral resilience to climate change and highlight the potential for spatially varying selection across complex coral reef seascapes to drive pronounced ecological divergence in climate-related traits.

Extreme seascape drives local recruitment and genetic divergence in brooding and spawning corals in remote north-west Australia.

Management strategies designed to conserve coral reefs threatened by climate change need to incorporate knowledge of the spatial distribution of inter- and intra-specific genetic diversity. We characterized patterns of genetic diversity and connectivity using single nucleotide polymorphisms (SNPs) in two reef-building corals to explore the eco-evolutionary processes that sustain populations in north-west Australia. Our sampling focused on the unique reefs of the Kimberley; we collected the broadcast spawning coral Acropora aspera (n = 534) and the brooding coral Isopora brueggemanni (n = 612) across inter-archipelago (tens to hundreds of kilometres), inter-reef (kilometres to tens of kilometres) and within-reef (tens of metres to a few kilometres) scales. Initial analysis of A. aspera identified four highly divergent lineages that were co-occurring but morphologically similar. Subsequent population analyses focused on the most abundant and widespread lineage, Acropora asp-c. Although the overall level of geographic subdivision was greater in the brooder than in the spawner, fundamental similarities in patterns of genetic structure were evident. Most notably, limits to gene flow were observed at scales <35 kilometres. Further, we observed four discrete clusters and a semi-permeable barrier to dispersal that were geographically consistent between species. Finally, sites experiencing bigger tides were more connected to the metapopulation and had greater gene diversity than those experiencing smaller tides. Our data indicate that the inshore reefs of the Kimberley are genetically isolated from neighbouring oceanic bioregions, but occasional dispersal between inshore archipelagos is important for the redistribution of evolutionarily important genetic diversity. Additionally, these results suggest that networks of marine reserves that effectively protect reefs from local pressures should be spaced within a few tens of kilometres to conserve the existing patterns of demographic and genetic connectivity.

Phylotranscriptomics confirms Alveopora is sister to Montipora within the family Acroporidae.

The genus Alveopora is a scleractinian coral taxon whose phylogenetic classification has recently changed from the family Poritidae to Acroporidae. This change, which was made based on single-locus genetic data, has led to uncertainty about the placement of Alveopora and the ability for deep evolutionary relationships in these groups to be accurately recovered and represented by limited genetic datasets. We sought to characterize the higher-level position of Alveopora using newly available transcriptome data to confirm its placement within Acroporidae and resolve its closest ancestor. Here we present an analysis of a new 2031 gene dataset that confirms the placement of Alveopora within Acroporidae corroborating other single-locus (COI, 16S and ITS) analyses and a mitogenome dataset. We also resolve the position of Alveopora as sister to the genus Montipora. This has allowed the re-interpretation of morphology, and a rediagnosis of the family Acroporidae and the genus Alveopora.

Under the karst: detecting hidden subterranean assemblages using eDNA metabarcoding in the caves of Christmas Island, Australia.

Subterranean ecosystems are understudied and challenging to conventionally survey given the inaccessibility of underground voids and networks. In this study, we conducted a eukaryotic environmental DNA (eDNA) metabarcoding survey across the karst landscape of Christmas Island, (Indian Ocean, Australia) to evaluate the utility of this non-invasive technique to detect subterranean aquatic 'stygofauna' assemblages. Three metabarcoding assays targeting the mitochondrial 16S rRNA and nuclear 18S genes were applied to 159 water and sediment samples collected from 23 caves and springs across the island. Taken together, our assays detected a wide diversity of chordates, cnidarians, porifera, arthropods, molluscs, annelids and bryozoans from 71 families across 60 orders. We report a high level of variation between cave and spring subterranean community compositions which are significantly influenced by varying levels of salinity. Additionally, we show that dissolved oxygen and longitudinal gradients significantly affect biotic assemblages within cave communities. Lastly, we combined eDNA-derived community composition and environmental (water quality) data to predict potential underground interconnectivity across Christmas Island. We identified three cave and spring groups that showed a high degree of biotic and abiotic similarity indicating likely local connectivity. This study demonstrates the applicability of eDNA metabarcoding to detect subterranean eukaryotic communities and explore underground interconnectivity.