Patterns of genetic connectivity among anchialine habitats: a case study of the endemic Hawaiian shrimp Halocaridina rubra on the island of Hawaii

Connectivity modelling identifies sources and sinks of coral recruitment within reef clusters

Connectivity aids the recovery of populations following disturbances, such as coral bleaching and tropical cyclones. Coral larval connectivity is a function of physical connectivity and larval behaviour. In this study, we used OceanParcels, a particle tracking simulator, with 2D and 3D velocity outputs from a high resolution hydrodynamic-biogeochemical marine model (RECOM) to simulate the dispersal and settlement of larvae from broadcast spawning Acropora corals in the Moore Reef cluster, northern Great Barrier Reef, following the annual spawning events in 2015, 2016 and 2017. 3D velocity simulations showed 19.40-68.80% more links and sinks than those of 2D simulations. Although the patterns of connectivity among sites vary over days and years, coral larvae consistently dispersed from east to west in the cluster domain, with some sites consistently acting as sources or sinks for local larval recruitment. Results can inform coral reef intervention plans for climate change, such as the design of marine protected areas and the deployment of proposed interventions within reef clusters. For example, the wider benefits of interventions (e.g., deployment of heat adapted corals) may be optimised when deployed at locations that are a source of larvae to others within comparable habitats across the reef cluster.

Primary microbial succession in the anchialine ecosystem

When new land is created, initial microbial colonization lays the foundation for further ecological succession of plant and animal communities. Primary microbial succession of new aquatic habitats formed during volcanic activity has received little attention. The anchialine ecosystem, which includes coastal ponds in young lava flows, offers an opportunity to examine this process. Here, we characterized microbial communities of anchialine habitats in Hawaii that were created during volcanic eruptions in 2018. Benthic samples from three habitats were collected ∼2 years after their formation and at later time points spanning ∼1 year. Sequence profiling (16S and 18S) of prokaryotic and eukaryotic communities was used to test whether communities were similar to those from older, established anchialine habitats, and if community structure changed over time. Results show that microbial communities from the new habitats were unlike any from established anchialine microbial communities, having higher proportions of Planctomycetota and Chloroflexi but lower proportions of green algae. Each new habitat also harbored its own unique community relative to other habitats. While community composition in each habitat underwent statistically significant changes over time, they remained distinctive from established anchialine habitats. New habitats also had highly elevated temperatures compared to other habitats. These results suggest idiosyncratic microbial consortia form during early succession of Hawaiian anchialine habitats. Future monitoring will reveal whether the early communities described here remain stable after temperatures decline and macro-organisms become more abundant, or if microbial communities will continue to change and eventually resemble those of established habitats. This work is a key first step in examining primary volcanic succession in aquatic habitats and suggests young anchialine habitats may warrant special conservation status.

Extensive gene flow among populations of the cavernicolous shrimp at the northernmost distribution margin in the Ryukyu Islands, Japan

Marine cave habitats in the Ryukyu Islands, Indo-West Pacific, are located at the northern edge of the distribution of many cave-dwelling species. At distribution margins, gene flow is often more restricted than that among core populations due to the smaller effective population size. Here, we used high-throughput sequencing technology to investigate the gene flow pattern among three sampling sites of a marine cave-dwelling species at the margin of its distribution range. We collected individuals of the barbouriid shrimp Parhippolyte misticia from three marine caves in the Ryukyu Islands and performed population genetic analyses by means of multiplexed inter-simple sequence repeat genotyping by sequencing. Based on 62 single-nucleotide polymorphism markers, no clear population structure or directional gene flow pattern was found among the three sites. These results were unexpected because previous studies of other stygobitic shrimps in this region did find significant population genetic structures and northward directional gene flow patterns. Together, these inconsistent findings imply that marine cave-dwelling species in the region have different mechanisms of larval dispersal. Future studies on larval ecology and the biotic and abiotic factors influencing gene flow patterns are needed to clarify the mechanisms underlying the population dynamics of marine cave-dwelling species.

Red Coloration in an Anchialine Shrimp: Carotenoids, Genetic Variation, and Candidate Genes

Red coloration is a widely distributed phenotype among animals, yet the pigmentary and genetic bases for this phenotype have been described in relatively few taxa. Here we show that the Hawaiian endemic anchialine shrimp Halocaridina rubra is red because of the accumulation of astaxanthin. Laboratory colonies of phylogenetically distinct lineages of H. rubra have colony-specific amounts of astaxanthin that are developmentally, and likely genetically, fixed. Carotenoid supplementation and restriction experiments failed to change astaxanthin content from the within-colony baseline levels, suggesting that dietary limitation is not a major factor driving coloration differences. A possible candidate gene product predicted to be responsible for the production of astaxanthin in H. rubra and other crustaceans is closely related to the bifunctional cytochrome P450 family 3 enzyme CrtS found in fungi. However, homologs to the enzyme thought to catalyze ketolation reactions in birds and turtles, CYP2J19, were not found. This work is one of the first steps in linking phenotypic variation in red coloration of H. rubra to genotypic variation. Future work should focus on (1) pinpointing the genes that function in the bioconversion of dietary carotenoids to astaxanthin, (2) examining what genomic variants might drive variation in coloration among discrete lineages, and (3) testing more explicitly for condition-dependent carotenoid coloration in crustaceans.

Dark offshoot: Phylogenomic data sheds light on the evolutionary history of a new species of cave brittle star

Caves are a useful system for testing evolutionary and biogeographic hypotheses, as they are isolated, and their environmental conditions have resulted in adaptive selection across different taxa. Although in recent years many more cave species have been discovered, cave-dwelling members of the class Ophiuroidea (brittle stars) remain scarce. Out of the more than two thousand species of brittle stars described to date, only three are regarded as true cave-dwellers. These occurrences represent rare colonising events, compared to other groups that are known to have successfully diversified in these systems. A third species from an anchihaline cave system in the Yucatan Peninsula, Mexico, has been previously identified from cytochrome oxidase I (COI) barcodes. In this study, we reassess the species boundaries of this putative cave species using a phylogenomic dataset (20 specimens in 13 species, 100 exons, 18.7 kbp). We perform species delimitation analyses using robust full-coalescent methods for discovery and validation of hypotheses on species boundaries, as well as infer its phylogenetic relationships with species distributed in adjacent marine regions, in order to investigate the origin of this cave-adapted species. We assess which hypotheses on the origin of subterranean taxa can be applied to this species by taking into account its placement within the genus Ophionereis and its demographic history. We provide a detailed description of Ophionereis commutabilis n. sp., and evaluate its morphological characters in the light of its successful adaptation to life in caves.

Phylogenetic evidence that both ancient vicariance and dispersal have contributed to the biogeographic patterns of anchialine cave shrimps

Cave shrimps from the genera Typhlatya, Stygiocaris and Typhlopatsa (Atyidae) are restricted to specialised coastal subterranean habitats or nearby freshwaters and have a highly disconnected distribution (Eastern Pacific, Caribbean, Atlantic, Mediterranean, Madagascar, Australia). The combination of a wide distribution and a limited dispersal potential suggests a large-scale process has generated this geographic pattern. Tectonic plates that fragment ancestral ranges (vicariance) has often been assumed to cause this process, with the biota as passive passengers on continental blocks. The ancestors of these cave shrimps are believed to have inhabited the ancient Tethys Sea, with three particular geological events hypothesised to have led to their isolation and divergence; (1) the opening of the Atlantic Ocean, (2) the breakup of Gondwana, and (3) the closure of the Tethys Seaway. We test the relative contribution of vicariance and dispersal in the evolutionary history of this group using mitochondrial genomes to reconstruct phylogenetic and biogeographic scenarios with fossil-based calibrations. Given that the Australia/Madagascar shrimp divergence postdates the Gondwanan breakup, our results suggest both vicariance (the Atlantic opening) and dispersal. The Tethys closure appears not to have been influential, however we hypothesise that changing marine currents had an important early influence on their biogeography.

The Three Domains of Conservation Genetics: Case Histories from Hawaiian Waters

The scientific field of conservation biology is dominated by 3 specialties: phylogenetics, ecology, and evolution. Under this triad, phylogenetics is oriented towards the past history of biodiversity, conserving the divergent branches in the tree of life. The ecological component is rooted in the present, maintaining the contemporary life support systems for biodiversity. Evolutionary conservation (as defined here) is concerned with preserving the raw materials for generating future biodiversity. All 3 domains can be documented with genetic case histories in the waters of the Hawaiian Archipelago, an isolated chain of volcanic islands with 2 types of biodiversity: colonists, and new species that arose from colonists. This review demonstrates that 1) phylogenetic studies have identified previously unknown branches in the tree of life that are endemic to Hawaiian waters; 2) population genetic surveys define isolated marine ecosystems as management units, and 3) phylogeographic analyses illustrate the pathways of colonization that can enhance future biodiversity. Conventional molecular markers have advanced all 3 domains in conservation biology over the last 3 decades, and recent advances in genomics are especially valuable for understanding the foundations of future evolutionary diversity.

Reproduction and development in Halocaridina rubra Holthuis, 1963 (Crustacea: Atyidae) clarifies larval ecology in the Hawaiian anchialine ecosystem

Larvae in aquatic habitats often develop in environments different from those they inhabit as adults. Shrimp in the Atyidae exemplify this trend, as larvae of many species require salt or brackish water for development, while adults are freshwater-adapted. An exception within the Atyidae family is the "anchialine clade," which are euryhaline as adults and endemic to habitats with subterranean fresh and marine water influences. Although the Hawaiian anchialine atyid Halocaridina rubra is a strong osmoregulator, its larvae have never been observed in nature. Moreover, larval development in anchialine species is poorly studied. Here, reproductive trends in laboratory colonies over a 5-y period are presented from seven genetic lineages and one mixed population of H. rubra; larval survivorship under varying salinities is also discussed. The presence and number of larvae differed significantly among lineages, with the mixed population being the most prolific. Statistical differences in reproduction attributable to seasonality also were identified. Larval survivorship was lowest (12% settlement rate) at a salinity approaching fresh water and significantly higher in brackish and seawater (88% and 72%, respectively). Correlated with this finding, identifiable gills capable of ion transport did not develop until metamorphosis into juveniles. Thus, early life stages of H. rubra are apparently excluded from surface waters, which are characterized by lower and fluctuating salinities. Instead, these stages are restricted to the subterranean (where there is higher and more stable salinity) portion of Hawaii's anchialine habitats due to their inability to tolerate low salinities. Taken together, these data contribute to the understudied area of larval ecology in the anchialine ecosystem.

Genetic and Morphological Characterization of Freshwater Shrimps (Caridina africana Kingsley, 1882) Reveals the Presence of Alien Shrimps in the Cape Floristic Region, South Africa

Morphological identification and molecular data (mtDNA COI) were used to resolve the taxonomic identity of a non-native freshwater shrimp in the Cape Floristic Region (CFR) of South Africa and to evaluate levels of genetic diversity and differentiation in the species' core natural distribution. The species was morphologically and genetically identified as Caridina africana Kingsley, 1882, whose main natural distribution is in the KwaZulu-Natal (KZN) Province, more than 1200 km from the point of new discovery. Subsequently, sequence data from natural populations occurring in seven rivers throughout KZN showed the presence of nuclear copies of the mtDNA COI gene (NUMTs) in 46 out of 140 individuals. Upon removal of sequences containing NUMTs, levels of genetic diversity were low in the alien population (possibly as a consequence of a bottleneck event), while varying levels of genetic diversity and differentiation were found in natural populations, indicating habitat heterogeneity, fragmentation and restricted gene flow between rivers. Following the present study, the alien shrimp has survived the Western Cape's winter and dispersed into a nearby tributary of the Eerste River System, hence posing an additional potential threat to endangered endemics. Understanding the biology of this alien species will aid detection and eradication procedures.

Phylogenetic utility, and variability in structure and content, of complete mitochondrial genomes among genetic lineages of the Hawaiian anchialine shrimp Halocaridina rubra Holthuis 1963 (Atyidae:Decapoda)

The Atyidae are caridean shrimp possessing hair-like setae on their claws and are important contributors to ecological services in tropical and temperate fresh and brackish water ecosystems. Complete mitochondrial genomes have only been reported from five of the 449 species in the family, thus limiting understanding of mitochondrial genome evolution and the phylogenetic utility of complete mitochondrial sequences in the Atyidae. Here, comparative analyses of complete mitochondrial genomes from eight genetic lineages of Halocaridina rubra, an atyid endemic to the anchialine ecosystem of the Hawaiian Archipelago, are presented. Although gene number, order, and orientation were syntenic among genomes, three regions were identified and further quantified where conservation was substantially lower: (1) high length and sequence variability in the tRNA-Lys and tRNA-Asp intergenic region; (2) a 317-bp insertion between the NAD6 and CytB genes confined to a single lineage and representing a partial duplication of CytB; and (3) the putative control region. Phylogenetic analyses utilizing complete mitochondrial sequences provided new insights into relationships among the H. rubra genetic lineages, with the topology of one clade correlating to the geologic sequence of the islands. However, deeper nodes in the phylogeny lacked bootstrap support. Overall, our results from H. rubra suggest intra-specific mitochondrial genomic diversity could be underestimated across the Metazoa since the vast majority of complete genomes are from just a single individual of a species.

A complex evolutionary history in a remote archipelago: phylogeography and morphometrics of the Hawaiian endemic Ligia isopods

Compared to the striking diversification and levels of endemism observed in many terrestrial groups within the Hawaiian Archipelago, marine invertebrates exhibit remarkably lower rates of endemism and diversification. Supralittoral invertebrates restricted to specific coastal patchy habitats, however, have the potential for high levels of allopatric diversification. This is the case of Ligia isopods endemic to the Hawaiian Archipelago, which most likely arose from a rocky supralittoral ancestor that colonized the archipelago via rafting, and diversified into rocky supralittoral and inland lineages. A previous study on populations of this isopod from O'ahu and Kaua'i revealed high levels of allopatric differentiation, and suggested inter-island historical dispersal events have been rare. To gain a better understanding on the diversity and evolution of this group, we expanded prior phylogeographic work by incorporating populations from unsampled main Hawaiian Islands (Maui, Moloka'i, Lana'i, and Hawai'i), increasing the number of gene markers (four mitochondrial and two nuclear genes), and conducting Maximum likelihood and Bayesian phylogenetic analyses. Our study revealed new lineages and expanded the distribution range of several lineages. The phylogeographic patterns of Ligia in the study area are complex, with Hawai'i, O'ahu, and the Maui-Nui islands sharing major lineages, implying multiple inter-island historical dispersal events. In contrast, the oldest and most geographically distant of the major islands (Kaua'i) shares no lineages with the other islands. Our results did not support the monophyly of all the supralittoral lineages (currently grouped into L. hawaiensis), or the monophyly of the terrestrial lineages (currently grouped into L. perkinsi), implying more than one evolutionary transition between coastal and inland forms. Geometric-morphometric analyses of three supralittoral clades revealed significant body shape differences among them. A taxonomic revision of Hawaiian Ligia is warranted. Our results are relevant for the protection of biodiversity found in an environment subject to high pressure from disturbances.

Phylogeography of the sponge Suberites diversicolor in Indonesia: insights into the evolution of marine lake populations

The existence of multiple independently derived populations in landlocked marine lakes provides an opportunity for fundamental research into the role of isolation in population divergence and speciation in marine taxa. Marine lakes are landlocked water bodies that maintain a marine character through narrow submarine connections to the sea and could be regarded as the marine equivalents of terrestrial islands. The sponge Suberites diversicolor (Porifera: Demospongiae: Suberitidae) is typical of marine lake habitats in the Indo-Australian Archipelago. Four molecular markers (two mitochondrial and two nuclear) were employed to study genetic structure of populations within and between marine lakes in Indonesia and three coastal locations in Indonesia, Singapore and Australia. Within populations of S. diversicolor two strongly divergent lineages (A & B) (COI: p = 0.4% and ITS: p = 7.3%) were found, that may constitute cryptic species. Lineage A only occurred in Kakaban lake (East Kalimantan), while lineage B was present in all sampled populations. Within lineage B, we found low levels of genetic diversity in lakes, though there was spatial genetic population structuring. The Australian population is genetically differentiated from the Indonesian populations. Within Indonesia we did not record an East-West barrier, which has frequently been reported for other marine invertebrates. Kakaban lake is the largest and most isolated marine lake in Indonesia and contains the highest genetic diversity with genetic variants not observed elsewhere. Kakaban lake may be an area where multiple putative refugia populations have come into secondary contact, resulting in high levels of genetic diversity and a high number of endemic species.

Multiple colonizations lead to cryptic biodiversity in an island ecosystem: comparative phylogeography of anchialine shrimp species in the Ryukyu Archipelago, Japan

Archipelagos of the Indo-West Pacific are considered to be among the richest in the world in biodiversity, and phylogeographic studies generally support either the center of origin or the center of accumulation hypothesis to explain this pattern. To differentiate between these competing hypotheses for organisms from the Indo-West Pacific anchialine ecosystem, defined as coastal bodies of mixohaline water fluctuating with the tides but having no direct oceanic connections, we investigated the genetic variation, population structure, and evolutionary history of three caridean shrimp species (Antecaridina lauensis, Halocaridinides trigonophthalma, and Metabetaeus minutus) in the Ryukyu Archipelago, Japan. We used two mitochondrial genes--cytochrome c oxidase subunit I (COI) and large ribosomal subunit (16S-rDNA)--complemented with genetic examination of available specimens from the same or closely related species from the Indian and Pacific Oceans. In the Ryukyus, each species encompassed 2-3 divergent (9.52%-19.2% COI p-distance) lineages, each having significant population structure and varying geographic distributions. Phylogenetically, the A. lauensis and M. minutus lineages in the Ryukyus were more closely related to ones from outside the archipelago than to one another. These results, when interpreted in the context of Pacific oceanographic currents and geologic history of the Ryukyus, imply multiple colonizations of the archipelago by the three species, consistent with the center of accumulation hypothesis. While this study contributes toward understanding the biodiversity, ecology, and evolution of organisms in the Ryukyus and the Indo-West Pacific, it also has potential utility in establishing conservation strategies for anchialine fauna of the Pacific Basin in general.

On the comparison of population-level estimates of haplotype and nucleotide diversity: a case study using the gene cox1 in animals

Estimates of genetic diversity represent a valuable resource for biodiversity assessments and are increasingly used to guide conservation and management programs. The most commonly reported estimates of DNA sequence diversity in animal populations are haplotype diversity (h) and nucleotide diversity (π) for the mitochondrial gene cytochrome c oxidase subunit I (cox1). However, several issues relevant to the comparison of h and π within and between studies remain to be assessed. We used population-level cox1 data from peer-reviewed publications to quantify the extent to which data sets can be re-assembled, to provide a standardized summary of h and π estimates, to explore the relationship between these metrics and to assess their sensitivity to under-sampling. Only 19 out of 42 selected publications had archived data that could be unambiguously re-assembled; this comprised 127 population-level data sets (n ≥ 15) from 23 animal species. Estimates of h and π were calculated using a 456-base region of cox1 that was common to all the data sets (median h=0.70130, median π=0.00356). Non-linear regression methods and Bayesian information criterion analysis revealed that the most parsimonious model describing the relationship between the estimates of h and π was π=0.0081 h(2). Deviations from this model can be used to detect outliers due to biological processes or methodological issues. Subsampling analyses indicated that samples of n>5 were sufficient to discriminate extremes of high from low population-level cox1 diversity, but samples of n ≥ 25 are recommended for greater accuracy.

Morphological and behavioral evidence for adaptive diversification of sympatric Hawaiian limpets (Cellana spp.)

The endemic Hawaiian limpets (Cellana exarata, Cellana sandwicensis, and Cellana talcosa), reside at different elevations on wave-exposed rocky shores and comprise a monophyletic lineage that diversified within Hawai'i. Here, I report phenotypic differences in shell, soft tissue, and behavioral characters among these limpets and discuss their potential utility in exploiting their respective niches. The high-shore limpet, C. exarata, is characterized by a tall round shell, short mantle tentacles, and long evasion distance when confronted by a predatory gastropod. The mid-shore limpet, C. sandwicensis, is characterized by a shorter oblong shell, long mantle tentacles, and a short evasion distance when confronted by a predatory snail. The low-shore, shallow-subtidal limpet, C. talcosa, is characterized by a flat shell that is thin in juveniles and disproportionately massive in large adults (relative to the other two species), and mantle tentacles of varying lengths (some individuals exhibit short tentacles, some long). These species-specific suites of characters are likely to confer specific fitness advantages on the high shore (C. exarata) where thermal and desiccation stress is severe, on the mid shore (C. sandwicensis) where hydrodynamic forces are severe, and on the low-shallow subtidal shore (C. talcosa) where pelagic predators have free access to the limpets. These data add to the growing body of evidence for adaptive diversification and speciation in the Hawaiian Cellana, and in marine species in general.

Shrimps down under: evolutionary relationships of subterranean crustaceans from Western Australia (Decapoda: Atyidae: Stygiocaris)

BACKGROUND

We investigated the large and small scale evolutionary relationships of the endemic Western Australian subterranean shrimp genus Stygiocaris (Atyidae) using nuclear and mitochondrial genes. Stygiocaris is part of the unique cave biota of the coastal, anchialine, limestones of the Cape Range and Barrow Island, most of whose nearest evolutionary relations are found in coastal caves of the distant North Atlantic. The dominance of atyids in tropical waters and their food resources suggest they are pivotal in understanding these groundwater ecosystems.

METHODOLOGY/PRINCIPLE FINDINGS

Our nuclear and mitochondrial analyses all recovered the Mexican cave genus Typhlatya as the sister taxon of Stygiocaris, rather than any of the numerous surface and cave atyids from Australia or the Indo-Pacific region. The two described Stygiocaris species were recovered as monophyletic, and a third, cryptic, species was discovered at a single site, which has very different physiochemical properties from the sites hosting the two described species.

CONCLUSIONS/SIGNIFICANCE

Our findings suggest that Stygiocaris and Typhlatya may descend from a common ancestor that lived in the coastal marine habitat of the ancient Tethys Sea, and were subsequently separated by plate tectonic movements. This vicariant process is commonly thought to explain the many disjunct anchialine faunas, but has rarely been demonstrated using phylogenetic techniques. The Cape Range's geological dynamism, which is probably responsible for the speciation of the various Stygiocaris species, has also led to geographic population structure within species. In particular, Stygiocaris lancifera is split into northern and southern groups, which correspond to population splits within other sympatric subterranean taxa.

The complete mitochondrial genome of the Hawaiian anchialine shrimp Halocaridina rubra Holthuis, 1963 (Crustacea: Decapoda: Atyidae)

Shrimp of the family Atyidae are important members of nearly all tropical (and most temperate) fresh and brackish water ecosystems in the world. To date, a complete mitochondrial genome from this important crustacean group has not been reported. Here, we present the complete mitochondrial DNA sequence of the Hawaiian atyid Halocaridina rubra [Holthuis, L.B., 1963. On red coloured shrimps (Decapoda, Caridea) from tropical land-locked saltwater pools. Zool. Meded.16, 261-279.] (Crustacea: Decapoda: Atyidae). The genome is a circular molecule of 16,065 bp and encodes the 37 mitochondrial genes (13 protein-coding, 22 tRNAs, and two rRNAs) typically found in the metazoa. Gene order and orientation in the H. rubra mitochondrial genome is syntenic with most malacostracans that have been examined to date. Of special note is the absence of the dihydrouridine (DHU) arm stem from tRNA(Tyr) and the use of CCG as an initiation codon for cytochrome oxidase subunit I (COI); these represent the first reported examples of such phenomena in the Malacostraca. Phylogenetic analyses utilizing complete mitochondrial sequences from other malacostracans place H. rubra as sister to Macrobrachium rosenbergii, which also belongs to the Infraorder Caridea. However, the placement of this infraorder, as well as the Infraorder Dendrobrachiata, in the phylogeny of the Decapoda varied depending on outgroup selection. Data from additional mitochondrial genomes, such as basal decapods like the Stenopodidea, should contribute to a better overall understanding of decapod phylogenetics.