Hidden diversity and phylogeographic history provide conservation insights for the edible seaweed Sargassum fusiforme in the Northwest Pacific

Multigene phylogenetics of Sargassum (Phaeophyceae) revealed low molecular diversity in contrast to high morphological variability in the NE Atlantic Ocean

Sargassum species play a key role in habitat formation in tropical and subtropical regions; however, species identification has been hampered by the phenological plasticity exhibited in response to environmental conditions and life history. Molecular phylogenetics has challenged taxa circumscriptions and proven critical in delimiting species in this genus. Yet, the Atlantic species of Sargassum remain poorly understood, and recent studies have shown low molecular diversity between the species in the NW Atlantic. Here, we expand the taxon sampling to the NE Atlantic to assess the diversity of Sargassum in the Atlantic basin, based on a comprehensive morphological and multigene approach. We selected genes commonly used in delineating species in this genus (ITS2, rbcLS, cox3, mtsp) and explored additional markers (cox2, nad6, psbC, clpC, atpB) to infer the phylogenetic relationships between the morphospecies observed in the NE Atlantic. Phylogenetic analyses using single-gene and multigene alignments including 185 new sequences confirmed the low molecular diversity and supported the distinction of a single clade in Sargassum section Sargassum of N Atlantic benthic species. In contrast, morphological analysis resulted in the identification of 10 species and three new morphospecies that we described here but opt not to equate to species until further molecular evidence is available. Our results were congruent with previous findings from the NW Atlantic and highlight the morphological and ecological diversity of Sargassum in the Atlantic. These results suggest a recent colonization and incipient speciation of Sargassum in the Atlantic basin and showcase the need of further high-throughput analyses.

Will climate change cause Sargassum beds in temperate waters to expand or contract? Evidence from the range shift pattern of Sargassum

Understanding the range shift patterns of foundation species (e.g., macroalgae) under future climatic conditions is critical for biodiversity conservation in coastal ecosystems. These predictions are typically made using species distribution models (SDMs), and severe habitat loss has been predicted for most brown algal forests. Nevertheless, some models showed that local adaptation within species can reduce range loss projections. In this study, we used the brown algae Sargassum fusiforme and Sargassum thunbergii, which are distributed in the Northwest Pacific, to determine whether climate change will cause the Sargassum beds in Northwest Pacific temperate waters to expand or contract. We divided S. fusiforme and S. thunbergii into northern and southern lineages, considering the temperature gradients and phylogeographic structures. We quantified the realized niches of the two lineages using an n-dimensional hypervolume. Significant niche differentiation was detected between lineages for both species, suggesting the existence of local adaptation. Based on these results, lineage-level SDMs were constructed for both species. The prediction results showed the different responses of different lineages to climate change. The suitable distribution area for both species was predicted to move northward, retaining part of the suitable habitat at low latitudes (along the East China Sea). Unfortunately, this expansion could not compensate for losing middle-low latitude areas. Our results have important implications for the future management and protection of macroalgae and emphasize the importance of incorporating intraspecific variation into species distribution predictions.

Contrasting population genetic structure of three semi-terrestrial brachyuran crabs on the coast of the Japanese archipelago

Anthropogenic activities have reduced ecotones between the ocean and land, which is likely to threaten the population of brackish-water brachyuran crabs. To assess the current status of these crabs, we examine the population genetic structures of three semi-terrestrial brachyuran crabs widely distributed along the coast of the Japan and to clarify factors determining their genetic structures. We collected 184 Orisarma dehaani, 252 Chiromantes haematocheir, and 151 Helice tridens crabs from 36 localities of the Japanese archipelago. Genome-wide SNP data from these crabs were analyzed using MIG-seq. Bayesian clustering of STRUCTURE and DAPC analysis were used to identify genetically disturbed populations and to visualize genetic differentiation between local populations. Genetic population structure showed clear differentiation between populations on the Pacific coast of the Tohoku region and on other Japanese coasts in O. dehaani, but not in C. haematocheir or H. tridens. The inbreeding coefficient of O. dehaani was significantly higher on the Pacific coast of the Tohoku region compared to other Japanese coasts. C. haematocheir and H. tridens had homogeneous genetic structures along the Japanese coast, but showed genetic differentiation of a local population at their range limits. Thus, O. dehaani showed little gene flow and clear genetic differentiation between populations in the Tohoku Pacific region and those on other Japanese coasts due to ocean currents. Although such a regional differentiation was not found in C. haematocheir and H. tridens, one population of C. haematocheir was genetically isolated at the edge of its distribution range and likely vulnerable to environmental changes.

Multiple genetic sources facilitate the northward range expansion of an intertidal oyster along China's coast

Coastal artificial structures on the former mudflats provide available habitats for the rocky intertidal species which can establish new populations in these emerging habitats over their former distribution range limits. As a former southern species, the oyster Crassostrea sikamea has become a pioneer and rapidly invaded the artificial shorelines in northern China. We used a seascape genomics approach to investigate the population structure and genetic sources of C. sikamea on the coastal artificial structures, which is crucial for understanding the genetic mechanisms driving species distribution range expansion and invasion pathway of intertidal species. Five C. sikamea populations, including two artificial substrate populations (WGZ and ZAP), one oyster reef population (LS), and two natural rocky shore populations (ZS and XM), were measured using single nucleotide polymorphism (SNPs) obtained from double digest restriction-site associated DNA sequencing (ddRAD-Seq). Redundancy analyses (RDA) were implemented for investigating the relationship between local temperature variables and the temperature adaptability of C. sikamea. Genetic diversity, direction and strength of gene flow, and population structure all revealed that the LS and ZS populations were the genetic sources for the oyster populations on the emerging northern coastal artificial structures. Results of RDA showed that there were different adaptive potentials for northern and southern populations to local temperature variables and the oyster reef population which frequently suffers from heat stress owned high heat adaptability. The ZS population as a genetic source nearby the Yangtze River estuary provided mass larvae for the northern populations, and the other genetic source, the heat-tolerant LS population, in the oyster reef played an important role in the post-settlement success by providing preadapted genotypes. These results highlight the importance of multiple sources with divergent adaptative capabilities for biological invasion, and also emphasize the importance of the oyster reef in coastal biodiversity and conservation.

ddRAD Sequencing-Based Scanning of Genetic Variants in Sargassum fusiforme

Sargassum fusiforme is a commercially important brown seaweed that has experienced significant population reduction both from heavy exploitation and degradation of the environment. Cultivated breed strains are also in a state of population mixing. These population stressors make it necessary to investigate the population genetics to discover best practices to conserve and breed this seaweed. In this study, the genetic diversity and population structure of S. fusiforme were investigated by the genome-wide SNP data acquired from double digest restriction site-associated DNA sequencing (ddRAD-seq). We found a low genetic diversity and a slight population differentiation within and between wild and cultivated populations, and the effective population size of S. fusiforme had experienced a continuous decline. Tajima’s D analysis showed the population contraction in wild populations may be related to copper pollution which showed a consistent trend with the increase of the sea surface temperature. The potential selection signatures may change the timing or level of gene expression, and further experiments are needed to investigate the effect of the mutation on relevant pathways. These results suggest an urgent need to manage and conserve S. fusiforme resources and biodiversity considering the accelerating change of the environment.

Phylogeographic pattern of a cryptoviviparous mangrove, Aegiceras corniculatum, in the Indo-West Pacific, provides insights for conservation actions

This study identified the historical geoclimatic factors which caused low genetic diversity and strong phylogeographic structure in a cryptoviviparous mangrove. The phylogeographic pattern was used to suggest conservation actions. Phylogeographic studies are used to understand the spatial distribution and evolution of genetic diversity, and have major conservation implications, especially for threatened taxa like the mangroves. This study aimed to assess the phylogeographic pattern of Aegiceras corniculatum, a cryptoviviparous mangrove, across its distribution range in the Indo-West Pacific (IWP) region. We genotyped 398 samples, collected from 37 populations, at four chloroplast DNA (cpDNA) loci, and identified the influence of historical processes on the contemporary population structure of the species. Low genetic diversity at the population level was observed. The evolutionary relationship between 12 cpDNA haplotypes suggested a strong phylogeographic structure, which was further validated by the clustering algorithms and proportioning of maximum variation among hierarchical population groups. The magnitude and direction of historical gene flow indicated that the species attained its wide distribution from its likely ancestral area of the Malay Archipelago. The divergence time estimates of the haplotypes indicated that the geoclimatic changes during the Pleistocene, especially the glacial sea-level changes and emergence of landmasses, hindered genetic exchange and created genetic differentiation between the phylogenetic groups. The species overwintered the last glacial maxima in multiple refugia in the IWP, as identified by the environmental niche modelling. Overall, our findings indicated that ancient glacial vicariance had influenced the present genetic composition of A. corniculatum, which was maintained by the current demographic features of this region. We discussed how these findings can be used to prioritize areas for conservation actions, restore disturbed habitats and prevent further genetic erosion.

Vegetative Reproduction Is More Advantageous Than Sexual Reproduction in a Canopy-Forming Clonal Macroalga under Ocean Warming Accompanied by Oligotrophication and Intensive Herbivory

Ocean warming and the associated changes in fish herbivory have caused polarward distributional shifts in the majority of canopy-forming macroalgae that are dominant in temperate Japan, but have little effect on the alga Sargassum fusiforme. The regeneration ability of new shoots from holdfasts in this species may be advantageous in highly grazed environments. However, little is known about the factors regulating this in Sargassum species. Moreover, holdfast tolerance to high-temperature and nutrient-poor conditions during summer has rarely been evaluated. In the present study, S. fusiforme holdfast responses to the combined effects of temperature and nutrient availability were compared to those of sexually reproduced propagules. The combined effects of holdfast fragmentation and irradiance on regeneration were also evaluated. Propagule growth rate values changed from positive to negative under the combination of elevated temperature (20 °C–30 °C) and reduced nutrient availability, whereas holdfasts exhibited a positive growth rate even at 32 °C in nutrient-poor conditions. The regeneration rate increased with holdfast fragmentation (1 mm segments), but was unaffected by decreased irradiance. These results suggest that S. fusiforme holdfasts have a higher tolerance to high-temperature and nutrient-poor conditions during summer than propagules, and regenerate new shoots even if 1-mm segments remain in shaded refuges for fish herbivory avoidance.

Climate-induced range shifts shaped the present and threaten the future genetic variability of a marine brown alga in the Northwest Pacific

Glaciation-induced environmental changes during the last glacial maximum (LGM) have strongly influenced species' distributions and genetic diversity patterns in the northern high latitudes. However, these effects have seldom been assessed on sessile species in the Northwest Pacific. Herein, we chose the brown alga Sargassum thunbergii to test this hypothesis, by comparing present population genetic variability with inferred geographical range shifts from the LGM to the present, estimated with species distribution modelling (SDM). Projections for contrasting scenarios of future climate change were also developed to anticipate genetic diversity losses at regional scales. Results showed that S. thunbergii harbours strikingly rich genetic diversity and multiple divergent lineages in the centre-northern range of its distribution, in contrast with a poorer genetically distinct lineage in the southern range. SDM hindcasted refugial persistence in the southern range during the LGM as well as post-LGM expansion of 18 degrees of latitude northward. Approximate Bayesian computation (ABC) analysis further suggested that the multiple divergent lineages in the centre-northern range limit stem from post-LGM colonization from the southern survived lineage. This suggests divergence due to demographic bottlenecks during range expansion and massive genetic diversity loss during post-LGM contraction in the south. The projected future range of S. thunbergii highlights the threat to unique gene pools that might be lost under global changes.

Contrasting Phylogeographic Patterns in Lumnitzera Mangroves Across the Indo-West Pacific

Mangroves are ecologically important forest communities in tropical and subtropical coasts, the effective management of which requires understanding of their phylogeographic patterns. However, these patterns often vary among different species, even among ecologically similar taxa or congeneric species. Here, we investigated the levels and patterns of genetic variation within Lumnitzera consisting of two species (L. racemosa and L. littorea) with nearly sympatric ranges across the Indo-West Pacific (IWP) region by sequencing three chloroplast DNA regions (for both species) and genotyping 11 nuclear microsatellite loci (for L. littorea). Consistent with findings in studies on other mangrove species, we found that both L. racemosa and L. littorea showed relatively high genetic variation among populations but low genetic variation within populations. Haplotype network and genetic clustering analyses indicated two well-differentiated clades in both L. racemosa and L. littorea. The relationship between geographic and genetic distances and divergence time estimates of the haplotypes indicated that limited dispersal ability of the propagules, emergence of land barriers during ancient sea-level changes, and contemporary oceanic circulation pattern in the IWP influenced the current population structure of the two species. However, the position of genetic break was found to vary between the two species: in L. racemosa, strong divergence was observed between populations from the Indian Ocean and the Pacific Ocean possibly due to land barrier effect of the Malay Peninsula; in L. littorea, the phylogeographic pattern was created by a more eastward genetic break along the biogeographic barrier identified as the Huxley's line. Overall, our findings strongly supported previous hypothesis of mangrove species divergence and revealed that the two Lumnitzera species have different phylogeographic patterns despite their close genetic relationship and similar current geographic distribution. The findings also provided references for the management of Lumnitzera mangroves, especially for the threatened L. littorea.

The Complete Plastid Genomes of Seven Sargassaceae Species and Their Phylogenetic Analysis

Sargassum is one of the most important genera of the family Sargassaceae in brown algae and is used to produce carrageenan, mannitol, iodine, and other economic substances. Here, seven complete plastid genomes of Sargassum ilicifolium var. conduplicatum, S. graminifolium, S. phyllocystum, S. muticum, S. feldmannii, S. mcclurei, and S. henslowianum were assembled using next-generation sequencing. The sizes of the seven circular genomes ranged from 124,258 to 124,563 bp, with two inverted regions and the same set of plastid genes, including 139 protein-coding genes (PCGs), 28 transfer (t)RNAs, and 6 ribosomal (r)RNAs. Compared with the other five available plastid genomes of Fucales, 136 PCGs were conserved, with two common ones shared with Coccophora langsdorfii, and one with S. fusiforme and S. horneri. The co-linear analysis identified two inversions of trnC(gca) and trnN(gtt) in ten Sargassum species, against S. horneri and C. langsdorfii. The phylogenetic analysis based on the plastid genomes of 55 brown algae (Phaeophyceae) showed four clades, whose ancient ancestor lived around 201.42 million years ago (Mya), and the internal evolutionary branches in Fucales started to be formed 92.52 Mya, while Sargassum species were divided into two subclades 14.33 Mya. Our novel plastid genomes provided evidence for the speciation of brown algae and plastid genomic evolution events.

Limited Genetic Connectivity Among Sargassum horneri (Phaeophyceae) Populations in the Chinese Marginal Seas Despite Their high Dispersal Capacity

Sargassum horneri is a habitat-forming species in the Northwest Pacific and an important contributor to seaweed rafts. In this study, 131 benthic samples and 156 floating samples were collected in the Yellow Sea and East China Sea (ECS) to test the effects of seaweed rafts on population structure and connectivity. Our results revealed high levels of genetic diversity in both benthic and floating samples based on concatenated mitochondrial markers (rpl5-rps3, rnl-atp9, and cob-cox2). Phylogenetic analyses consistently supported the existence of two lineages (lineages I and II), with divergence dating to c. 0.692 Mya (95% HPD: 0.255-1.841 Mya), indicating that long-term isolation may have occurred during the mid-Pleistocene (0.126-0.781 Mya). Extended Bayesian skyline plots demonstrated a constant population size over time in lineage I and slight demographic expansion in lineage II. Both lineages were found in each marginal sea (including both benthic and floating samples), but PCoA, FST , and AMOVA analyses consistently revealed deep genetic variation between regions. Highly structured phylogeographic pattern supports limited genetic connectivity between regions. IMA analyses demonstrated that asymmetric gene flow between benthic populations in the North Yellow Sea (NYS) and ECS was extremely low (ECS→NYS, 2Nm = 0.6), implying that high dispersal capacity cannot be assumed to lead to widespread population connectivity, even without dispersal barriers. In addition, there were only a few shared haplotypes between benthic and floating samples, suggesting the existence of hidden donors for the floating masses in the Chinese marginal seas.

Phylogeographic diversification and postglacial range dynamics shed light on the conservation of the kelp Saccharina japonica

Studies of postglacial range shifts could enhance our understanding of seaweed species' responses to climate change and hence facilitate the conservation of natural resources. However, the distribution dynamics and phylogeographic diversification of the commercially and ecologically important kelp Saccharina japonica in the Northwest Pacific (NWP) are still poorly surveyed. In this study, we analyzed the evolutionary history of S. japonica using two mitochondrial markers and 24 nuclear microsatellites. A STRUCTURE analysis revealed two partially isolated lineages: lineage H, which is scattered along the coast of Japan; and lineage P, which occurs along the west coast of the Japan Sea. Ecological niche modeling projections to the Last Glacial Maximum (LGM) revealed that the southern coasts of the Japan Sea and the Pacific side of the Oshima and Honshu Peninsulas provided the most suitable habitats for S. japonica, implying that these regions served as ancient refugia during the LGM. Ancient isolation in different refugia may explain the observed divergence between lineages P and H. An approximate Bayesian computation analysis indicated that the two lineages experienced post-LGM range expansion and that postglacial secondary contact occurred in Sakhalin. Model projections into the year 2,100 predicted that S. japonica will shift northwards and lose its genetic diversity center on the Oshima Peninsula in Hokkaido and Shimokita Peninsula in Honshu. The range shifts and evolutionary history of S. japonica improve our understanding of how climate change impacted the distribution range and diversity of this species and provide useful information for the conservation of natural resources under ongoing environmental change in the NWP.

Phylogeography and genetic connectivity of the marine macro-alga Sargassum ilicifolium (Phaeophyceae, Ochrophyta) in the northwestern Pacific1

The evolutionary influences of historical and contemporary factors on the population connectivity and phylogeographic structure of a brown seaweed, Sargassum ilicifolium, were elucidated using the nuclear ITS2 and mitochondrial COI markers for the collections newly sampled within its distribution range in the northwestern Pacific (NWP). Significant genetic structure at variable levels was identified between populations (pairwise FST ) and among populations grouped by geographical proximity (ΦCT among regions). The adjacent groups of populations with moderate structure revealed from AMOVA appeared to have high genetic connectivity. However, a lack of genealogical concordance with the geographic distribution was uncovered for S. ilicifolium from the NWP. Such genetic homogeneity is interpreted as a result of the interaction between postglacial recolonization and dynamic oceanic current regimes in the region. Two separated glacial refugia, the South China Sea and the Okinawa Trough, in the marginal seas of east China were recognized based on the presence of endemic haplotypes and high haplotype diversity in the populations at southern China and northeast of Taiwan. Populations persisting in these refugia may have served as the source for recolonization in the NWP with the rise of sea level during the warmer interglacial periods. The role of oceanic currents in maintaining genetic connectivity of S. ilicifolium in the region was further corroborated by the coherence between the direction of oceanic currents and that of gene flow, especially along the eastern coast of Taiwan. This study underlines the interaction between historical postglacial recolonization and contemporary coastal hydrodynamics in contributing to population connectivity and distribution for this tropical seaweed in the NWP.

Historical isolation and contemporary gene flow drive population diversity of the brown alga Sargassum thunbergii along the coast of China

BACKGROUND

Long-term survival in isolated marginal seas of the China coast during the late Pleistocene ice ages is widely believed to be an important historical factor contributing to population genetic structure in coastal marine species. Whether or not contemporary factors (e.g. long-distance dispersal via coastal currents) continue to shape diversity gradients in marine organisms with high dispersal capability remains poorly understood. Our aim was to explore how historical and contemporary factors influenced the genetic diversity and distribution of the brown alga Sargassum thunbergii, which can drift on surface water, leading to long-distance dispersal.

RESULTS

We used 11 microsatellites and the plastid RuBisCo spacer to evaluate the genetic diversity of 22 Sargassum thunbergii populations sampled along the China coast. Population structure and differentiation was inferred based on genotype clustering and pairwise F ST and allele-frequency analyses. Integrated genetic analyses revealed two genetic clusters in S. thunbergii that dominated in the Yellow-Bohai Sea (YBS) and East China Sea (ECS) respectively. Higher levels of genetic diversity and variation were detected among populations in the YBS than in the ECS. Bayesian coalescent theory was used to estimate contemporary and historical gene flow. High levels of contemporary gene flow were detected from the YBS (north) to the ECS (south), whereas low levels of historical gene flow occurred between the two regions.

CONCLUSIONS

Our results suggest that the deep genetic divergence in S. thunbergii along the China coast may result from long-term geographic isolation during glacial periods. The dispersal of S. thunbergii driven by coastal currents may facilitate the admixture between southern and northern regimes. Our findings exemplify how both historical and contemporary forces are needed to understand phylogeographical patterns in coastal marine species with long-distance dispersal.

Hidden diversity and phylogeographic history provide conservation insights for the edible seaweed Sargassum fusiforme in the Northwest Pacific

Understanding the evolutionary processes that have created diversity and the genetic potential of species to adapt to environmental change is an important premise for biodiversity conservation. Herein, we used mitochondrial trnW‐L and cox3 and plastid rbcL‐S data sets to analyze population genetic variation and phylogeographic history of the brown alga Sargassum fusiforme, whose natural resource has been largely exterminated in the Asia–Northwest Pacific in the past decades. Phylogenetic trees and network analysis consistently revealed three major haplotype groups (A, B, and C) in S. fusiforme, with A and B distributed in the Japan‐Pacific coast. Group C consisted of three subgroups (C1, C2, and C3) which were distributed in the Sea of Japan, the Yellow–Bohai Sea, and East China Sea, respectively. Isolation‐with‐migration (IMa) analysis revealed that the three groups diverged approximately during the mid‐Pleistocene (c. 756–1,224 ka). Extended Bayesian skyline plots (EBSP) showed that groups A and B underwent relatively long‐term stable population size despite a subsequent rapid demographic expansion, while subgroups C2 and C3 underwent a sudden expansion at c. 260 ka. F_ST and AMOVA detected low population‐level genetic variation and high degrees of divergence between groups. The cryptic diversity and phylogeographic patterns found in S. fusiforme not only are essential to understand how environmental shifts and evolutionary processes shaped diversity and distribution of coastal seaweeds but also provide additional insights for conserving and managing seaweed resources and facilitate predictions of their responses to future climate change and habitat loss.