Unraveling the predictive role of temperature in the gut microbiota of the sea urchin Echinometra sp. EZ across spatial and temporal gradients.

Shifts in microbial communities represent a rapid response mechanism for host organisms to respond to changes in environmental conditions. Therefore, they are likely to be important in assisting the acclimatization of hosts to seasonal temperature changes as well as to variation in temperatures across a species' range. The Persian/Arabian Gulf is the world's warmest sea, with large seasonal fluctuations in temperature (20℃ - 37℃) and is connected to the Gulf of Oman which experiences more typical oceanic conditions (<32℃ in the summer). This system is an informative model for understanding how symbiotic microbial assemblages respond to thermal variation across temporal and spatial scales. Here, we elucidate the role of temperature on the microbial gut community of the sea urchin Echinometra sp. EZ and identify microbial taxa that are tightly correlated with the thermal environment. We generated two independent datasets with a high degree of geographic and temporal resolution. The results show that microbial communities vary across thermally variable habitats, display temporal shifts that correlate with temperature, and can become more disperse as temperatures rise. The relative abundances of several ASVs significantly correlate with temperature in both independent datasets despite the >300 km distance between the furthest sites and the extreme seasonal variations. Notably, over 50% of the temperature predictive ASVs identified from the two datasets belonged to the family Vibrionaceae. Together, our results identify temperature as a robust predictor of community-level variation and highlight specific microbial taxa putatively involved in the response to thermal environment.

A Chromosome-level Genome Assembly of the Highly Heterozygous Sea Urchin Echinometra sp. EZ Reveals Adaptation in the Regulatory Regions of Stress Response Genes.

Echinometra is the most widespread genus of sea urchin and has been the focus of a wide range of studies in ecology, speciation, and reproduction. However, available genetic data for this genus are generally limited to a few select loci. Here, we present a chromosome-level genome assembly based on 10x Genomics, PacBio, and Hi-C sequencing for Echinometra sp. EZ from the Persian/Arabian Gulf. The genome is assembled into 210 scaffolds totaling 817.8 Mb with an N50 of 39.5 Mb. From this assembly, we determined that the E. sp. EZ genome consists of 2n = 42 chromosomes. BUSCO analysis showed that 95.3% of BUSCO genes were complete. Ab initio and transcript-informed gene modeling and annotation identified 29,405 genes, including a conserved Hox cluster. E. sp. EZ can be found in high-temperature and high-salinity environments, and we therefore compared E. sp. EZ gene families and transcription factors associated with environmental stress response ("defensome") with other echinoid species with similar high-quality genomic resources. While the number of defensome genes was broadly similar for all species, we identified strong signatures of positive selection in E. sp. EZ noncoding elements near genes involved in environmental response pathways as well as losses of transcription factors important for environmental response. These data provide key insights into the biology of E. sp. EZ as well as the diversification of Echinometra more widely and will serve as a useful tool for the community to explore questions in this taxonomic group and beyond.

Physiological Variation in Response to Vibrio and Hypoxia by Aquacultured Eastern Oysters in the Southeastern United States.

Eastern oysters (Crassostrea virginica) have long been recognized as model organisms of extreme environmental tolerance, showing resilience to variation in temperature, salinity, hypoxia, and microbial pathogens. These phenotypic responses, however, show variability between geographic locations or habitats (e.g., tidal). Physiological, morphological, and genetic differences occur in populations throughout a species' geographical range, which may have been shaped by regional abiotic and biotic variations. Few studies of C. virginica have explored the combined factors of physiological mechanisms of divergent phenotypes between locations and the genetic relationships of individuals between these locations. To characterize genetic relationships of four locations with aquacultured oysters along the North Carolina and Virginia coast, we sequenced a portion of cytochrome oxidase subunit I (COI) that revealed significant variation in haplotype distribution between locations. We then measured mitochondrial physiology and expression of the innate immunity response of hemocytes to lab acclimation and combined stress conditions to compare basal expression and stress response in oysters between these locations. For stress sensing genes, toll-like receptors had the strongest location-specific response to hypoxia and Vibrio, whereas mannose receptor and a stress-receptor were specific to hypoxia and bacteria, respectively. The expression of stress response genes also showed location-specific and stressor-specific changes in expression, particularly for big defensin and the complement gene Cq3. Our results further suggested that genetic similarity of oysters from different locations was not clearly related to physiological and molecular responses. These results are informative for understanding the range of physiological plasticity for stress responses in this commercially important oyster species. They also have implications in the oyster farming industry as well as conservation efforts to restore endangered native oyster beds.

Population Genomic Analyses of the Sea Urchin Echinometra sp. EZ across an Extreme Environmental Gradient.

Extreme environmental gradients represent excellent study systems to better understand the variables that mediate patterns of genomic variation between populations. They also allow for more accurate predictions of how future environmental change might affect marine species. The Persian/Arabian Gulf is extreme in both temperature and salinity, whereas the adjacent Gulf of Oman has conditions more typical of tropical oceans. The sea urchin Echinometra sp. EZ inhabits both of these seas and plays a critical role in coral reef health as a grazer and bioeroder, but, to date, there have been no population genomic studies on this or any urchin species in this unique region. E sp. EZ's life history traits (e.g., large population sizes, large reproductive clutches, and long life spans), in theory, should homogenize populations unless nonneutral processes are occurring. Here, we generated a draft genome and a restriction site-associated DNA sequencing data set from seven populations along an environmental gradient across the Persian/Arabian Gulf and the Gulf of Oman. The estimated genome size of E. sp. EZ was 609 Mb and the heterozygosity was among the highest recorded for an echinoderm at 4.5%. We recovered 918 high-quality SNPs from 85 individuals which we then used in downstream analyses. Population structure analyses revealed a high degree of admixture between all sites, although there was population differentiation and significant pairwise FST values between the two seas. Preliminary results suggest migration is bidirectional between the seas and nine candidate loci were identified as being under putative natural selection, including one collagen gene. This study is the first to investigate the population genomics of a sea urchin from this extreme environmental gradient and is an important contribution to our understanding of the complex spatial patterns that drive genomic divergence.

The complete mitochondrial genome of the sea urchin, Echinometra sp. EZ.

The complete mitogenome of Echinometra sp. EZ has been described and fully annotated in this study. Phylogenetic analysis of cytochrome c oxidase subunit I (COI) from six Echinometra species confirms that our sample is E. sp. EZ. The mitogenome is 15,698 bp in length and contains 13 protein-coding genes, 22 tRNAs, 2 rRNAs, and a non-coding region with an identical organization to other Echinoidea. The E. sp. EZ mitogenome shared ∼99.1% identity to the published Echinometra mathaei mitogenome, differing by 147 SNPs. The E. sp. EZ mitogenome will serve as a resource that can be applied to disentangling the Echinometra species complex and to future population genetic studies of this ecologically important sea urchin species.

Host specificity of Symbiodinium variants revealed by an ITS2 metahaplotype approach.

Analysis of the widely used ITS region is confounded by the presence of intragenomic variants (IGVs). In Symbiodinium, the algal symbionts of reef building corals, deep-sequencing analyses are used to characterise communities within corals, yet these analyses largely overlook IGVs. Here we consider that distinct ITS2 sequences could represent IGVs rather than distinct symbiont types and argue that symbionts can be distinguished by their proportional composition of IGVs, described as their ITS2 metahaplotype. Using our metahaplotype approach on Minimum Entropy Decomposition (MED) analysis of ITS2 sequences from the corals Acropora downingi, Cyphastrea microphthalma and Playgyra daedalea, we show the dominance of a single species-specific Symbiodinium C3 variant within each coral species. We confirm the presence of these species-specific symbionts using the psbA non-coding region. Our findings highlight the importance of accounting for IGVs in ITS2 analyses and demonstrate their capacity to resolve biological patterns that would otherwise be overlooked.

Species-specific trends in the reproductive output of corals across environmental gradients and bleaching histories.

Coral populations in the Persian Gulf have a reputation for being some of the toughest in the world yet little is known about the energetic constraints of living under temperature and salinity extremes. Energy allocation for sexual reproduction in Gulf corals was evaluated relative to conspecifics living under milder environmental conditions in the Oman Sea. Fecundity was depressed at Gulf sites in two Indo-Pacific merulinid species (Cyphastrea microphthalma and Platygyra daedalea) but not in a regionally endemic acroporid (Acropora downingi). Gulf populations of each species experienced high temperature bleaching at the onset of gametogenesis in the study but fecundity was only negatively impacted in P. daedalea and A. downingi. Large population sizes of C. microphthalma and P. daedalea in the Gulf are expected to buffer reductions on colony-level fecundity. However, depleted population sizes of A. downingi at some Gulf sites equate to low reef-wide fecundity and likely impede outcrossing success.

Levels of genetic diversity and taxonomic status of Epinephelus species in United Arab Emirates fish markets.

Understanding the patterns of genetic diversity of fish species is essential for marine conservation and management. This is particularly important in the Arabian Gulf where marine life is subject to extreme environmental conditions that could impact genetic diversity. Here we assess genetic diversity of the most commercially important fish in the United Arab Emirates; groupers (Epinephelus spp.). Sequencing of 973 bp mitochondrial DNA from 140 tissue samples collected in four main fish markets revealed 58 haplotypes clustered within three groups. Data analysis revealed the presence of three distinct Epinephelus species being marketed as one species (hammour): Epinephelus coioides, Epinephelus areolatus and Epinephelus bleekeri. We report species-specific genetic markers and demonstrate that all three species exhibit relatively low levels of genetic variation, reflecting the effect of overfishing and environmental pressures. In light of the genetic evidence presented here, conservation and management of groupers in the UAE warrant the implementation of species-specific measures.