Next Article in Journal
Germination Stimulant Activity of Isothiocyanates on Phelipanche spp.
Next Article in Special Issue
A Critical Gap in Seagrass Protection: Impact of Anthropogenic Off-Shore Nutrient Discharges on Deep Posidonia oceanica Meadows
Previous Article in Journal
Development and Application of a Cultivar-Specific Sequence-Characterized Amplified Region (SCAR) Marker for the Detection of Chrysanthemum morifolium Ramat. ‘Daboju’
Previous Article in Special Issue
Taxonomic Revision of Hook-Forming Acrosorium (Delesseriaceae, Rhodophyta) from the Northwestern Pacific Based on Morphology and Molecular Data
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Plants
Volume 11
Issue 5
10.3390/plants11050605
plants-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

More than What Meets the Eye: Differential Spatiotemporal Distribution of Cryptic Intertidal Bangiales

by
Fernanda P. Cid Alda
1,*,
Nelson Valdivia
2,3 and
Marie-Laure Guillemin
3,4,5,*
1
Scientific and Technological Bioresource Nucleus (BIOREN), Universidad de La Frontera, Avenida Francisco Salazar 1145, Temuco 4780000, Chile
2
Instituto de Ciencias Marinas y Limnológicas, Facultad de Ciencias, Campus Isla Teja, Universidad Austral de Chile, Valdivia 5090000, Chile
3
Centro FONDAP de Investigación de Ecosistemas Marinos de Altas Latitudes (IDEAL), Valdivia 5090000, Chile
4
Instituto de Ciencias Ambientales y Evolutivas, Universidad Austral de Chile, Casilla 567, Valdivia 5090000, Chile
5
CNRS, Sorbonne Université, IRL 3614, Evolutionary Biology and Ecology of Algae, Station Biologique de Roscoff, CS 90074, F-29688 Roscoff, France
*
Authors to whom correspondence should be addressed.
Plants 2022, 11(5), 605; https://doi.org/10.3390/plants11050605
Submission received: 6 December 2021 / Revised: 31 January 2022 / Accepted: 10 February 2022 / Published: 24 February 2022
(This article belongs to the Special Issue Biodiversity in Marine Plants)
Browse Figures
Versions Notes

Abstract

:
Morphologically similar but genetically distinct species have been termed cryptic and most have been assumed to be ecologically similar. However, if these species co-occur at a certain spatial scale, some niche differences at finer scales should be expected to allow for coexistence. Here, we demonstrate the existence of a disjointed distribution of cryptic bladed Bangiales along spatial (intertidal elevations) and temporal (seasons) environmental gradients. Bladed Bangiales were identified and quantified across four intertidal elevations and four seasons for one year, at five rocky intertidal sites (between 39° S and 43° S) in southern Chile. Species determination was based on partial sequences of the mitochondrial cytochrome c oxidase 1 (COI) gene amplification. To assess species gross morphology, thallus shape, color, and maximum length and width were recorded. Hundreds of organisms were classified into nine Bangiales species belonging to three genera (i.e., Fuscifolium, Porphyra, and Pyropia), including five frequent (>97% of specimens) and four infrequent species. All species, except for Pyropia saldanhae, had been previously reported along the coasts of Chile. The thallus shape and color were very variable, and a large overlap of the maximum width and length supported the cryptic status of these species. Multivariate analyses showed that the main variable affecting species composition was intertidal elevation. Species such as Py. orbicularis were more abundant in low and mid intertidal zones, while others, such as Po. mumfordii and Po. sp. FIH, were principally observed in high and spray elevations. Despite all numerically dominant species being present all year long, a slight effect of seasonal variation on species composition was also detected. These results strongly support the existence of spatial niche partitioning in cryptic Bangiales along the Chilean rocky intertidal zone.

1. Introduction

The study of species co-occurring along environmental gradients allows us to improve our mechanistic understanding of biodiversity. For instance, cryptic species are commonly found co-occurring in organisms such as corals [1,2], nematodes [3], rotifers [4], and red algae [5,6] living in aquatic environments subjected to strong physical or chemical gradients (e.g., light intensity, desiccation stress, or salinity). Cryptic species have been classically defined as genetically distinct taxa that have been erroneously classified under a single nominal species name because they are, at least superficially, morphologically indistinguishable [7,8]. These cryptic taxa have been shown to be especially common in the marine realm [9]. Cryptic species that are morphologically similar have been assumed to be ecologically similar, which implies that they require very similar resources [7,8]. However, if these species are actually coexisting, some niche differences would be expected [1,2,9,10,11].
Contemporary coexistence theory predicts that niche partitioning allows for stable long-term coexistence of co-occurring species if their differences in competitive abilities (i.e., “fitness differences”) are not extremely large [12,13]. Therefore, cryptic species being truly ecologically similar and competing for the same resources should be observed in the same habitat only transiently [14,15]. Alternatively, Hubbell’s [16] unified neutral theory of biodiversity and biogeography posits that coexistence could be due to stochastic events of reproduction, death, and dispersal along with speciation but without any effect of niche partitioning or selection. Yet, such an extreme ecological equivalence between coexisting species has proven to be elusive in nature, and minor niche differences have been shown to underpin the slow competitive exclusion of inferior competitors (e.g., [17]). Cryptic species, despite largely similar morphologies, have been shown to exhibit subtle differential physiological tolerances, leading to niche partitioning and stable coexistence [4].
The rocky intertidal zone provides an ideal habitat to assess the association between the composition of cryptic species and environmental gradients over space and time. In the intertidal zone, stress increases with elevation, due to desiccation and sun exposure, favoring rapid vertical changes in community composition [18,19]. The strength of desiccation, jointly with other factors, such as the distribution of consumers, can affect the temporal and spatial distribution and composition of resource species (e.g., macroalgae [20,21]). Moreover, intertidal habitats in temperate biogeographic zones can be exposed to important seasonal variations in abiotic environmental conditions, affecting sessile invertebrate and macroalga performance and community structure [22,23,24,25,26]. These sharp environmental intertidal gradients offer a variety of environmental niches for locally co-occurring (i.e., a few meters) species. If intertidal cryptic species are not fully ecologically similar, therefore, we should expect a significant association between the occurrence and abundance (i.e., community composition) of cryptic species and environmental gradients over space and time [27].
Multiple macroalgae cryptic species have been identified co-occurring along rocky intertidal habitats. For example, cryptic species of the genus Ectocarpus (brown algae) show differences in attachment substrate and tidal zonation at a very small scale (i.e., a few meters [28,29]). Regarding other groups of cryptic macroalgae, such as bladed Bangiales (red algae), only a few studies have been performed along the intertidal gradient. In New Zealand, nine species were recorded along a three-year study in a single location, but no differences in seasonal growth or small-scale distribution were detected between the two most common species [30]. In central Chile, a small-scale study conducted at Maitencillo beach reported the presence of the two most common cryptic Bangiales (i.e., Porphyra luchea and Pyropia variabilis) across several types of habitats, such as mid and high intertidal habitats and rocky platforms [31]. However, Porphyra longissima was restricted to rock boulders surrounded by sand in Maitencillo, a habitat where no other species were recorded, hinting at a possible differentiation of ecological niches among bladed Bangiales [31]. This idea is supported by the slightly distinct gene expression and ecophysiological responses of Pyropia orbicularis and Py. variabilis to desiccation [32]. Similarly, variations in ecophysiological responses have been used to explain differences in species distribution among cryptic red algae of the genus Bostrychia [33]. In addition, strong differences in the temporal patterns of appearance and growth were detected between the seven species of bladed Bangiales present in New Hampshire (USA), with some aseasonal species present all year round (as Po. ‘umbilicalis’), while others were highly seasonal (as Po. ‘yezoensis’) [34]. Taken jointly, these previous studies suggest that ecophysiological differences and/or differences in competitive abilities could explain the long-term coexistence of groups of closely related cryptic bladed Bangiales along temperate coasts.
In Chile, as in other part of the word, bladed Bangiales have attracted considerable interest because of their ecological and economical importance. Indeed, this group includes species commonly known as “nori”, which represent the most valuable macroalgal crop, with a market value estimated to be more than USD 1.3 billion per year in China, Japan, and Korea [35]. In Chile, the bladed Bangiales are commonly known as “Luche”. They are harvested and sold fresh or dry for local consumption and represent an important economic resource for the Chilean artisanal fishermen [31,36,37]. A first genetic analysis carried out along the southeastern Pacific coast has uncovered a high number of genetically distinct taxa, most of them still unnamed [6]. Chilean Bangiales lack conspicuous diagnostic morphological differences that could allow unambiguous species determination. Indeed, even when taking into account microscopic and macroscopic traits and the developmental features of reproductive tissues, morphology did not allow the four more common species sampled in Maitencillo to be distinguished effectively [31]. Accurate species identification based on morphology is extremely difficult, which explains why Chilean Bangiales have been misidentified for decades. For example, many specimens sampled along the Chilean coast were, for many years, incorrectly assumed to be Porphyra columbina, a species that we now know to be endemic of the sub-Antarctic waters of New Zealand [38,39]. Moreover, critical basic information, such as the ecological requirement for growth and reproduction, is still lacking for these species, limiting the development of aquaculture programs for these important resources for both human food and new bioproducts.
Here, we test the hypothesis that niche differences will lead to distinct spatiotemporal patterns of several cryptic bladed Bangiales along spatial (intertidal elevations) and temporal (seasons) environmental gradients. As an alternative hypothesis, we predict that ecological similarity and, thus, competitive exclusion due to fitness differences generates transient patterns in the abundance of these species. Under a neutral scenario, with species lacking major niche and fitness differences, a random distribution of cryptic bladed Bangiales along intertidal elevations and seasons of the year is expected. To test these hypotheses, bladed Bangiales were collected from five rocky intertidal sites (platforms and/or boulders located between the 39° S and 43° S) across four intertidal elevations and four seasons for one year. Cryptic species determination was based on partial sequences of the mitochondrial cytochrome c oxidase 1 (COI) gene and we assembled a large dataset of morphological attributes and molecular data for these key primary producers of the Chilean intertidal zone.

2. Results

A total of 1990 individuals were sampled. No Bangiales could be observed in the spray-elevation zone during summer in Niebla (two unsampled spray-elevation plots; five samples each). Of these 1990 individuals, 1727 (i.e., 87%) were identified at the species level using COI sequences. For the remaining 263 individuals, low-quantity and/or low-quality DNA extraction did not allow for COI amplification and sequencing. Nine Bangiales species were identified (Figure 1A). The most abundant species were Pyropia sp. CHJ (579 individuals, 34% of the sequenced Bangiales), Pyropia orbicularis (523 individuals, 30% of the sequenced Bangiales), Porphyra mumfordii (351 individuals, 20% of the sequenced Bangiales), Porphyra sp. FIH (170 individuals, 10% of the sequenced Bangiales), and Pyropia sp. CHH (66 individuals, 4% of the sequenced Bangiales). Four rare species of Bangiales were detected in our study area: Porphyra longissima (24 individuals identified), Porphyra luchea (7 individuals identified), Fuscifolium sp. CHA (4 individuals identified), and Pyropia saldanhae (3 individuals identified) (Table 1).
These 4 species corresponded to less than 2.5% of the sequenced Bangiales. No clear difference in gross morphology was detected between the nine Bangiales species for which measures of maximum width and length largely overlapped (Figure 1B). In the most abundant species, Po. mumfordii, Po. sp. FIH, and Py. orbicularis, thallus shape was very variable, with specimens showing lancelolate, elongated rosette, and rosette-like blades. Some species, such as Py. sp. CHH and Py. sp. CHJ, were characterized by more rosette-like or elongated rosette blades. The highest homogeneity in terms of thallus shape was observed in Po. longissima, for which 19 of the 21 measured specimens were characterized by lancelolate blades. Most sampled specimens were of brownish color. However, thallus color was highly variable in Po. longissima and Py. sp. CHJ, where thalli varied from red to brown, green, and yellowish, depending on the sampled specimen (Supplementary Table S1).
The first axis of the principal coordinates analysis (PCoA) separated the samples according to intertidal elevation, with low- and mid-elevation samples located on the left and high- and spray-elevation samples located on the right of the axis (Figure 2). Nevertheless, a relatively high overlap between both groups was observed, in particular for the mid-elevation samples obtained during autumn (light-blue squares in Figure 2). Accordingly, the global PERMANOVA showed an elevation by season interactive effect on species composition (Table 2, R2 = 0.31). Pairwise comparisons indicated that species composition significantly varied along the vertical intertidal stress gradient across the year (Table 2). The exception was summer, when species composition at the mid and low elevations did not vary statistically (Table 2, p = 0.339).
The abundance of Py. sp. CHJ increased from low to high intertidal elevations in spring and summer but decreased along the same environmental gradient during autumn and winter (Figure 3; see also post-hoc comparisons in Table 3). Elevation and season accounted for 27% of the spatial variation in the abundance of Py. sp. CHJ (R2c = 0.27). The entire model, including sampling sites, accounted for 40% of the variation in the abundance of this species (R2m = 0.40). The model significantly fit the observed data of Py. sp. CHJ abundance (LRT: χ216 = 179.67, p < 0.001). Pyropia orbicularis decreased in abundance from low to spray elevations (Figure 3). The magnitude of this variation, however, differed among seasons, with broader variations observed during spring and summer than autumn and winter (Figure 3; Table 3). The fixed factors accounted for 89% of the variation in Py. orbicularis abundance and the entire model explained an 89.5% of the variation, indicating a small effect of the differences among sites on the abundance of this species (R2c and R2m = 0.890 and 0.895, respectively; LRT: χ216 = 402.9, p < 0.001). The abundance of Po. mumfordii tended to increase with increasing intertidal vertical environmental stress across the year (Figure 3). These differences, however, were statistically significant only between the high and low elevations during summer and winter (Table 3). Accordingly, R2c and R2m were 0.88 and 0.96, respectively (LRT: χ216 = 167, p < 0.001). Porphyra sp. FIH increased in abundance from low to spray elevations (Figure 3). The species was absent at low elevations during spring and summer and occurred with minimal abundances at low elevations during autumn and winter (Figure 3). The statistical model significantly fit the observed data of Porphyra sp. FIH abundance (R2c and R2m = 0.96 and 0.99, respectively; LRT: χ216 = 108.4, p < 0.001). The comparatively low abundances of Porphyra sp. FIH resulted in statistically significant comparisons only during autumn (Table 3). Pyropia sp. CHH occurred mainly at low and mid elevations in spring and summer (Figure 3; R2c and R2m = 0.57 and 0.78, respectively; LRT: χ216 = 145.73, p < 0.001). As above, the low abundances of this species led to a statistically non-significant post-hoc test (Table 3).

3. Discussion

Using the molecular and morphological data of hundreds of specimens that were collected seasonally in the southern part of the Chilean coast (39–43° S) for one year, we were able to detect patterns of niche partitioning among cryptic bladed Bangiales along the intertidal gradient. Indeed, differences in species composition according to intertidal elevation were detected, with Pyropia orbicularis being more abundant in the low and mid intertidal zones and Porphyra mumfordii and Porphyra sp. FIH being more abundant in the high and spray elevations. In total, we identified nine Bangiales species belonging to three genera (i.e., Fuscifolium, Porphyra, and Pyropia) in the study area, including five dominant (i.e., those species accounted for >97% of the sequenced specimens) and four infrequent species. Measurements of maximum thallus width and length showed overlapping patterns of gross morphology, supporting the cryptic status of these nine Bangiales species. Excepting Porphyra longissima, which was characterized by more lanceolate thalli, thallus shape was highly variable between the specimens of the same species, confirming a high intraspecific morphological variability. Even if we detected variation in community composition between seasons, all dominant species were a-seasonal and present all year long in our study area. As hypothesized for other cryptic species, we propose that niche partitioning could allow for the stable long-term coexistence of morphologically similar co-occurring species of bladed Bangiales.

3.1. Responses of Bladed Bangiales to Vertical Stress Gradients in Intertidal Habitats

Our results support the hypothesis that niche differences lead to distinct spatiotemporal patterns of cryptic bladed Bangiales in the southern part of the Chilean coast. Intertidal elevation accounted for most of the variation in species composition. Thus, and despite the observed overlaps in species abundances, our results support the existence of spatial niche partitioning in cryptic Bangiales along the intertidal zone.
The constant alternation of low and high tides determines that high intertidal elevations are characterized by strong physiological constraints for organisms living in intertidal habitats in terms of desiccation [40], high temperatures [41], nutrient shortages [42], and high UV radiation [43,44]. Thus, the upper distribution boundary of intertidal species is related to species physiological limits. Contrastingly, the lower limit of intertidal species distribution seems to be mainly regulated by biological interactions (e.g., consumption and competition) [45,46,47,48]. In southern Chile, for instance, small-sized grazers have been shown to strongly control the growth of intertidal macroalgae and sessile invertebrates over succession [49,50]. Thus, the interplay of biotic and abiotic factors generates a sharp gradient of environmental conditions and a broad niche space in intertidal habitats [51,52,53].
Even though bladed Bangiales were morphologically very similar, some species occurred at different intertidal elevations. For example, while Py. orbicularis was mainly present at low/mid-intertidal elevations, both Po. mumfordii and Po. sp. FIH were distributed higher on the shore. These patterns were not detected in the previous studies focusing on the center part of the Chilean coast [31,32,39,54]. The species Py. orbicularis was first described as colonizing the upper and mid intertidal area of the Maitencillo beach (central Chile, [39]), but the species’ abundance across intertidal elevations was not estimated in that study. In the same site, various studies also reported Py. orbicularis as commonly encountered in the high intertidal elevation [31,32]. It is possible that ecotypes with distinct preferences in intertidal elevations exist within Py. orbicularis, a species for which clear genetic differences have been detected between the samples from extreme south, south, central and north Chile [6]. In the same way, Po. mumfordii was reported from the low to the upper intertidal elevations in Montemar (central Chile, [54]) but without any estimation of the species abundance. However, and confirming our results, higher abundances of Po. mumfordii at high intertidal elevations have also been observed in the northern hemisphere [55,56].
As for most bladed Bangiales, Chilean species have been classically described as extremely resistant to the stressing conditions characteristic of the high intertidal zone (e.g., desiccation, temperature, and UV; [40,44,57,58,59]) and the existence of species-specific tolerance limits acting as key determinants of the upper vertical limits of the distribution in these species is still unclear. For example, experimental studies have shown that Chilean bladed Bangiales are highly tolerant to desiccation, being able to lose more than 90% of their water contents and recover cellular activities after only five minutes of rehydration [59]. Moreover, these species present a low sensitivity to UV radiation and a rapid recovery from solar stress [44]. In other coastal regions, differences in physiological response to stress have been detected between Bangiales species occupying distinct tidal levels. Indeed, strong differences have been detected between the high- and low-intertidal species in terms of nutrient uptake and growth rate (northwest Atlantic; [60,61]), membrane damage (North East Pacific; [41]), and photosynthetic activities (North West Pacific, [62]). Collectively, these works suggest that fine-scale niche partitioning between bladed Bangiales species, mostly due to dissimilarities in physiological tolerance to stress, can generate distinct distribution patterns. These results have also been confirmed in filamentous Bangiales species for which distinct adaptation to temperature and salinity could explain, at least in part, species distribution [63]. Whether the cryptic bladed Bangiales identified in our study respond differently to abiotic environmental factors is still an open question, and further manipulative experiments are now needed to understand the physiological responses to environmental stress in these species assemblages.

3.2. Seasonal Variation in Intertidal Species Assemblages

Even if the effect of seasonal variation on species composition was much weaker than that imposed by intertidal elevation, our results show that seasonal changes affect bladed Bangiales communities in south Chile. For example, the distribution patterns observed for Py. sp. CHJ varied between seasons, with an increase in abundance from low to high intertidal elevations in spring and summer but a decrease along the same environmental gradient during autumn and winter.
Nonetheless, all the relatively abundant species observed in the study area were present across the four seasons sampled. While we might need a multiyear sampling of species occurrences and abundances to confirm these results, the present study suggests that the assemblage of species is very persistent over time (i.e., there are no transient species). The relatively strong temporal persistence of the species assemblage observed in our study contrasts with previous studies demonstrating the existence of both persistent and transient cryptic bladed Bangiales in intertidal habitats (e.g., New Hampshire, USA [34]; New Zealand [30]. For rare species, it is much harder to conclude on a possible persistent or transient status. However, even for species for which only a few specimens were sampled, no clear seasonal patterns were detected. For example, Po. longissima (24 specimens collected in total) was observed during the four seasons, Po. luchea (6 specimens collected in total) was observed during both winter and summer, F. sp. CHA (4 specimens collected in total) was observed during spring and summer, and Py. saldanhae (3 specimens collected in total) was observed during spring and winter. In northern Atlantic populations of Po. umbilicalis, the production of spores (i.e., asexual spores or neutral spores) has been shown to be highly seasonal, with a higher spore release from fall to early spring [57,64,65]. These results are concordant with observations of the bladed Bangiales recruitment peak during autumn in New Zealand [30]. It is possible that most species of Chilean bladed Bangiales also reproduce during the colder periods, affecting the cover and biomass of Bangiales within each site but not necessarily the structure of the entire species assemblage [65].

3.3. Diversity and Biogeography of Bladed Bangiales

In the present study, measures of maximum blade width and length exhibited an overlap between the nine identified species. These results support the cryptic status of these species, at least at the level of gross morphology. Similar results were obtained for the bladed Bangiales species community located in central Chile, where only Po. longissima could be clearly differentiated morphologically from the rest of the bladed Bangiales, due to its long and thin blades [31]. At our study sites, Po. longissima was also characterized by long, thin, lanceolate thalli. However, Po. longissima shared these characteristics with some specimens from other species. A high level of variability in terms of thallus color, shape, and size was observed within species, probably due to phenotypic plasticity. A previous study conducted in central Chile detected clear phenotypic plasticity in Po. mumfordii, with specimens from the high and mid intertidal habitats showing long, thin, lanceolate thalli, while the ones from the low intertidal habitats were characterized by much wider elongated rosettes [54]. Taken together, these results are in agreement with the previous studies that showed that morphological characters alone could lead to inaccurate species determination in bladed Bangiales [66,67,68,69]. In many taxa, most detected cryptic species represent recently diverged entities that still share the same gross morphology [70]. However, molecular studies have also revealed the existence of non-monophyletic complexes of cryptic species, and morphological similarities in these cases have been associated with evolutionary convergence, morphological stasis, or developmental constraints [70]. In our species assemblage, Py. sp. CHJ and Py. orbicularis have been shown to be close genetic groups [6] and could represent recently diverged species that have not accumulated any morphological differences yet. However, the nine species sampled in the present study are part of three very divergent genera (i.e., Fuscifolium, Porphyra, and Pyropia; [68]), and the evolutionary convergence or morphological stasis linked to life in the highly stressful intertidal zone (e.g., [71]) could be hypothesized in bladed Bangiales.
Eight of the nine species sampled in our study sites had previously been reported in Chile, with six of them sampled along the southern part of the coast [6]. Two rare species, F. sp. CHA and Po. luchea, were previously known only in the northern part of Chile (F. sp. CHA was reported in Puerto Oscuro, 31° S and Po. luchea was reported in Maitencillo, 32° S and Chañaral de Aceituno, 29° S; [6,31]). Thus, our results demonstrate that the distribution areas of both F. sp. CHA and Po. luchea extend well into the southern part of the country. Another dominant and commonly encountered species, Po. mumfordii, also presents a much wider distribution than the one proposed by Guillemin and collaborators (i.e., reported from 39° S up to 41° S, 2016), extending between 8° S Salaverry, Peru [72], and 43° S (Melinka, Chile, present study).
The present study is the first report of Py. saldanhae in South America, a species considered to be endemic to the west coast of South Africa and Namibia [73,74,75]. Genetic data show that Py. saldanhae forms a well-supported clade with various species of the Falkland Islands and New Zealand, and the authors suggested that Py. saldanhae has speciated along the South African coasts after a Pleistocene-Pliocene west-to-east colonization event [75,76]. Actually, rafting-mediated transport along the Antarctic Circumpolar Current has been shown to favor macroalgae dispersal between isolated landmasses in the Southern Ocean [77,78,79]. The presence of Py. saldanhae in south Chile, though, opens questions about the possible origin of this population. Even if we cannot completely rule out a recent introduction by human marine transport, the geographical position of our study region does not link to any major South Africa to south east Pacific shipping routes and is not located near any major Chilean harbor. It is also possible that the species is locally rare but presents an ample distribution, including the Chile Patagonia, Falkland Islands, and South Africa. Renewed efforts, including extensive sampling in the cold waters of the southern hemisphere and the acquisition of a multigenic data set, are now needed to better resolve the evolutionary history of the bladed Bangiales in the region.

4. Materials and Methods

4.1. Study Area

Samples were collected during each season in five sampling sites: San Carlos (39°51′46.52″ S/73°26′28.24″ W), Los Liles (39°53′31.19″ S/73°29′6.42″ W), Niebla (39°52′6.97″ S/73°24′7.05″ W), Pilolcura (39°40′20.54″ S/73°21′10.72″ W), and Melinka (43°53′ 56.13″ S/73°44′13.41″ W). The first four sites are located in the Los Ríos region while Melinka is located in the Aysén region (Ascension Island; Figure 4). Both regions are separated by approximately 500 km (Figure 4).
Abiotic environmental data were downloaded from Bio-ORACLE [80,81] (www.bio-oracle.org, accessed on 4 January 2022; Supplementary Table S1 in the Supplementary Materials). Present-day climate layers were obtained at a 5 armin spatial resolution. The maximum surface current velocity ranged from 0.398 m s−1 (San Carlos and Los Liles) to 0.095 m s−1 (Melinka, Supplementary Table S1). The maximal salinity was relatively higher in Los Ríos (between 34.12 and 34.14 PSU) than in Aysén (33. 03 PSU). The maximal surface seawater temperature was highest at Pilolcura (15.23 °C, Los Ríos) and lowest at Melinka (13.88 °C, Aysén). The maximal photosynthetically active radiation (PAR) peaked at Niebla (64.02 μE m−2 d−1) and was lowest at Melinka (47.59 μE m−2 d−1; Supplementary Table S1).
The sites in the Los Ríos region correspond to wave semi-exposed rocky shores. However, water circulation has been reported as highly complex in the area, due to the effect of local upwelling, freshwater plumes, and offshore eddies [82]. Two of the four sites, San Carlos and Niebla, are located at the mouth of the Valdivia River. Los Liles and Pilolcura, on the other hand, correspond to more open costal habitats [83].
In the Aysén region, Ascension Island is part of the large archipelago of Las Guaitecas, an area characterized by an intricate coastline with many small islands separated by channels. Melinka is located along the island’s coast facing the Corcovado gulf and represents a fairly protected area from the prevailing winds, flowing from west to east (see also Supplementary Table S1).
Rocky intertidal habitats in both study regions are dominated in terms of biomass by sessile species, such as the acorn barnacles Jehlius cirratus and Notochthamalus scabrosus (high intertidal), the purple mussel Perumytilus purpuratus, the red corticated alga Mazzaella laminarioides (mid intertidal), and large brown kelps, such as Lessonia spicata, Durvillaea incurvata, and Macrocystis pyrifera (low intertidal) [44,84,85,86]. The assemblage of consumers comprises macrograzers, such as Fissurella picta, Chiton granosus, and Scurria zebrina, and mesograzers, such as the pulmonate gastropod Siphonaria lessoni, littorinids, and amphipods [84,86]. Mesograzers are particularly abundant in the Los Ríos region, likely stemming from the overexploitation and functional extinction of predators and larger grazers; consequently, mesograzers have been demonstrated to exert a strong top-down control on the assembly of sessile communities [50,87].
The samplings were conducted between late May and early June 2020 for the autumn season, between late July and early August 2020 for the winter season, during December 2020 for the spring season, and between late February and early March 2021 for the summer season. A slight difference in sampling timing was due to the limited access to Ascension Island, especially during bad weather conditions.

4.2. Bladed Bangiales Sample Collection

Samples were taken from rocky platforms and boulders within four intertidal elevations (Figure 5). The upper and lower limits of the Mazzaella laminarioides belt were used as the major local biological markers to define the distinct intertidal levels sampled (Figure 5). The first intertidal level, named “low”, corresponded to samples encountered mixed with species such as Ahnfeltiopsis spp. or Sarcothalia crispata and completely embedded within the M. laminarioides belt. These samples were located at the limit between the midlittoral and infralittoral zones and generally corresponded to the bladed Bangiales specimens observed lowest at the shore. The second level, named “mid”, corresponded to samples encountered at the upper limit of the M. laminarioides belt or slightly higher. The third level, named “high”, corresponded to samples located at the limit between the supralittoral and midlittoral zones, where no M. laminarioides were observed. The fourth level, named “splash zone”, corresponded to the highest elevation of the supralittoral zone where bladed Bangiales could still be found. Samples were collected by group of five nearby specimens, all located within a circle of 25 cm and bagged together. When possible, neighboring specimens presenting different phenotypes (i.e., distinct color, form, or size) were collected. For the first three levels, six samples were collected and noted A to F, corresponding to “low”, G to L to “mid” and M to R to the “high” intertidal levels (Figure 5). For the “splash zone” two samples were collected, noted S and T. The samples collected within the same level were separated by at least 2 m. When available, at each season and in each locality, 100 algae specimens were collected in total (i.e., 20 samples of five specimens each).
Algal blades were spread flat to the best of our ability and photographed jointly with a size scale before cutting a small portion of tissue that was then stored in silica gel for further molecular work. Only complete blades, still wet, were photographed.

4.3. Morphological Traits

Using the photographs taken of the blades spread flat, the maximum lengths and maximum widths of the thalli were measured using the software ImageJ (version Java 1.8.0) (http://imagej.nih.gov/ij/, accessed on 24 May 2021). Some other characteristics such as color (brown, green, yellow, and red), and form (rosette, elongated rosette, and lanceolate blade) were also noted.

4.4. Species Determination Based on Molecular Data

Algal material was subjected to cell disruption for 30 s in a Mini-BeadBeater (BioSpec Inc., Bartlesville, OK, USA) and DNA was extracted using the Plant DNA kit (Omega BIO-TEK, Norcross, GA, USA) following the provider’s protocol. Quantity and quality of the extracted DNA was determined using a NanoDrop Lite spectrophotometer (Thermo Scientific, Waltham, MA, USA). A fragment of the mitochondrial cytochrome c oxidase 1 (COI) gene was amplified using the primer pair GazF1 (5′-TCA ACA AAT CAT AAA GAT ATT GG -3′) and GazR1 (5′-ACT TCT GGA TGT CCA AAA AAY CA -3′) following the PCR protocol described by Saunders [88], with a melting temperature of 56 °C. The PCR products were checked using 1.5% agarose electrophoresis gels. Sanger sequencing was performed using the GazF1 primer in the AUSTRAL-omics Core-Facility (Universidad Austral de Chile). The sequences were checked by hand using Chromas v.2.6.6 (http://technelysium.com.au/wp/chromas/, accessed on 10 May 2021), and species identification was performed using the GenBank Basic Local Alignment Search Tool (BLAST) to compare the newly obtained COI sequences with the ones already deposited in the public databases in NCBI GenBank (https://www.ncbi.nlm.nih.gov/, accessed on 20 December 2021). Species determination was based on the percent identity BLAST metric and a minimum identity threshold of 3% was used (i.e., reflecting the congeneric interspecific divergences observed for Chilean bladed Bangiales for the COI) [6]. The nucleotide sequences obtained in this study were deposited in the GenBank database; all accession numbers are available in the Supplementary Table S2.

4.5. Statistical Analyses

Composition of cryptic species (i.e., the combination of species occurrences and abundances) was analyzed with permutational multivariate analyses of variance (PERMANOVA; [89]. The model included intertidal elevation (low, mid, high, and spray level) and season (spring, summer, autumn, and winter) as fixed and crossed factors. The species abundance data (number of individuals) were transformed to Bray–Curtis dissimilarities before the analysis, and 999 permutations of raw data were constrained within sites (Pilolcura, San Carlos, Niebla, Los Liles, and Melinka). After fitting the global model, we conducted between-elevation PERMANOVA post-hoc pairwise comparisons within each season. To control the type I error rate, we used a “treatment” contrast in which each elevation was compared against the low-tide elevation (reference group). Moreover, the spatial patterns of the community structure were described in a constrained principal coordinates analysis (PCoA) ordination [90]. This method allowed us to show the multivariate axes that accounted for most of the variation in the species abundance dataset according to our predictive model (i.e., separate and joint elevation by season effects on species composition).
Individual counts of numerically dominant species were separately analyzed with zero-inflated hierarchical generalized linear models. The error structure was modeled as a negative binomial distribution and model parameters were estimated by means of maximum likelihood. The models included elevation and season as fixed and crossed factors and site as a hierarchical factor. Model goodness-of-fit was assessed by means of conditional and marginal pseudo-coefficients of determination (R2c and R2m, respectively; [83,84]). Conditional R2c represents the fit of the fixed factors (elevation and season); marginal R2m represents that of the entire model (including the random effects). Plots of fitted vs. observed values were used as model diagnostics. Statistical fit was tested with a likelihood ratio test (LRT) in which a chi-squared test (χ2) was used to contrast the explained deviance of the model against that of a null model that included only the intercept. After fitting the global models, we conducted Tukey-corrected post-hoc comparisons. As above, we used a “treatment” contrast with low elevation as a reference group.
The vegan R package was used for running PERMANOVA and PCoA, glmmTMB was used for the zero-inflated models, MuMIn was used for pseudo-R2, and the tidyverse, cowplot, RColorBrewer, and wesanderson packages were used for data preparation and plotting [90,91,92,93,94,95,96,97,98,99,100]. All statistical analyses were conducted in R version 3.0.3 [93].

5. Conclusions

Our study provides useful insights into the mechanistic understanding of the growing diversity of cryptic species. The observed differences in the occurrence and abundance of cryptic bladed Bangiales along the rocky intertidal zone in south Chile hint at niche partitioning within this assemblage. This, in turn, has probably favored the co-existence in time of these Bangiales species and indicates that they are not fully ecologically similar. Future research on physiological responses to biotic and abiotic environmental stress would help us to understand the mechanisms underpinning the observed spatial patterns. In bladed Bangiales, as in other species characterized by a complex life cycle, both the microscopic sporophyte (also known as the conchocelis phase) and the macroscopic gametophytic phases are important to study. Indeed, laboratory-based experimental studies have shown that the sporophyte and gametophyte phases respond differently to abiotic stresses (light and temperature [101,102]) and grazing pressure by intertidal mollusks [103]). We suggest that similar morphologies do not equal similar environmental tolerances, which challenges the widely accepted use of comparative anatomy to define functional types, e.g., [104] and, thus, to understand how individual species influence the community structure and ecosystem functioning.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/plants11050605/s1, Table S1: Abiotic environmental variables at each sampling site. Data were downloaded from Bio-ORACLE (www.bio-oracle.org), Table S2: Bladed Bangiales specimens sampled for this study. Season, locality, intertidal elevation, morphological traits, GenBank accession numbers, and percentage of identity between each sequence and the best BLAST in GenBank are given.

Author Contributions

Conceptualization, F.P.C.A. and M.-L.G.; data curation, F.P.C.A., N.V. and M.-L.G.; formal analysis, N.V.; writing—original draft, F.P.C.A., N.V. and M.-L.G.; writing—review and editing, F.P.C.A., N.V. and M.-L.G. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the FONDECYT Postdoctoral Grant no. 3200728 (ANID, ex CONICYT, Chile) awarded to F.P.C.A. N.V. and M.-L.G. were supported by the FONDAP grant no. 15150003 (FONDAP IDEAL). While writing, N.V. was supported by FONDECYT grant no. 1190529.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Reported results can be found at ncbi (https://www.ncbi.nlm.nih.gov/, accessed on 20 December 2021), accession numbers OM141383-OM141394 and OM215205-OM216838; see Table S2 for more details.

Acknowledgments

We thank Enrique Higuera, Alvaro Aguilar, and Oscar Alvarado for their support in sampling the Melinka site; Suany Quesada, María Paz Melo Gaymer, Francisca Heise, and Lorena Gallardo for their help in the field; and Sarai Morales, Suany Quezada, and María Paz Melo Gaymer for helping with DNA extraction and sequencing. F.P.C.A. would like to give special thanks to Michel Abanto Marin, postdoc advisor, and María de La Luz Mora for guidance and support. Three anonymous reviewers provided useful criticism that improved a previous version of this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Johnston, E.C.; Wyatt, A.S.J.; Leichter, J.J.; Burgess, S.C. Niche differences in co-occurring cryptic coral species (Pocillopora spp.). Coral Reefs 2021, 1–12. [Google Scholar] [CrossRef]
  2. Bongaerts, P.; Cooke, I.R.; Ying, H.; Wels, D.; Haan, S.D.; Hernandez-Agreda, A.; Brunner, C.A.; Dove, S.; Englebert, N.; Eyal, G.; et al. Morphological stasis masks ecologically divergent coral species on tropical reefs. Curr. Biol. 2021, 31, 2286–2298.e8. [Google Scholar] [CrossRef]
  3. de Meester, N.; Derycke, S.; Bonte, D.; Moens, T. Salinity effects on the coexistence of cryptic species: A case study on marine nematodes. Mar Biol. 2011, 158, 2717–2726. [Google Scholar] [CrossRef]
  4. Montero-Pau, J.; Ramos-Rodríguez, E.; Serra, M.; Gómez, A. Long-Term Coexistence of Rotifer Cryptic Species. PLoS ONE 2011, 6, e21530. [Google Scholar] [CrossRef]
  5. Pardo, C.; Lopez, L.; Peña, V.; Hernandez, J.; Le Gall, L.; Bárbara, I.; Barreiro, R. A Multilocus Species Delimitation Reveals a Striking Number of Species of Coralline Algae Forming Maerl in the OSPAR Maritime Area. PLoS ONE 2014, 9, e104073. [Google Scholar] [CrossRef] [Green Version]
  6. Guillemin, M.-L.; Contreras-Porcia, L.; Ramírez, M.E.; Macaya, E.; Contador, C.B.; Woods, H.; Wyatt, C.; Brodie, J. The bladed Bangiales (Rhodophyta) of the South Eastern Pacific: Molecular species delimitation reveals extensive diversity. Mol. Phylogenetics Evol. 2016, 94, 814–826. [Google Scholar] [CrossRef] [PubMed]
  7. Chenuil, A.; Cahill, A.E.; Délémontey, N.; Du Salliant du Luc, E.; Fanton, H. Problems and Questions Posed by Cryptic Species. A Framework to Guide Future Studies. In From Assessing to Conserving Biodiversity: Conceptual and Practical Challenges; Casetta, E., da Silva, J.M., Vecchi, D., Eds.; Springer Nature: Berlin, Germany, 2019; p. 452. [Google Scholar]
  8. Bickford, D.; Lohman, D.J.; Sodhi, N.S.; Ng, P.K.; Meier, R.; Winker, K.; Ingram, K.K.; Das, I. Cryptic species as a window on diversity and conservation. Trends Ecol. Evol. 2007, 22, 148–155. [Google Scholar] [CrossRef]
  9. Knowlton, N. Molecular genetic analyses of species boundaries in the sea. Hydrobiologia 2000, 420, 73–90. [Google Scholar] [CrossRef]
  10. Trontelj, P.; Fišer, C. Perspectives: Cryptic species diversity should not be trivialised. Syst. Biodivers. 2009, 7, 1–3. [Google Scholar] [CrossRef]
  11. Mayfield, M.M.; Levine, J.M. Opposing effects of competitive exclusion on the phylogenetic structure of communities. Ecol. Lett. 2010, 13, 1085–1093. [Google Scholar] [CrossRef]
  12. Van Campenhout, J.; Derycke, S.; Moens, T.; Vanreusel, A. Differences in life-histories refute ecological equivalence of cryptic species and provide clues to the origin of bathyal Halomonhystera (nematoda). PLoS One. 2014, 9, e111889. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Chesson, P. Mechanisms of Maintenance of Species Diversity. Annu. Rev. Ecol. Syst. 2000, 31, 343–366. [Google Scholar] [CrossRef] [Green Version]
  14. Spaak, J.W.; Godoy, O.; De Laender, F. Mapping species niche and fitness differences for communities with multiple interaction types. Oikos 2021, 130, 2065–2077. [Google Scholar] [CrossRef]
  15. Zhang, D.-Y.; Lin, K.; Hanski, I. Coexistence of cryptic species. Ecol. Lett. 2004, 7, 165–169. [Google Scholar] [CrossRef]
  16. Hubbell, S.P. The Unified Neutral Theory of Biodiversity and Biogeography. In Monographs in Population Biology Vol. 83; Levin, S.A., Horn, H.S., Eds.; Princeton University Press: Princeton, NJ, USA, 2001. [Google Scholar]
  17. Shinen, J.L.; Navarrete, S.A. Lottery coexistence on rocky shores: Weak niche differentiation or equal competitors engaged in neutral dynamics? Am. Nat. 2014, 183, 342–362. [Google Scholar] [CrossRef] [PubMed]
  18. Stephenson, T.A.; Stephenson, A. The Universal Features of Zonation Between Tide-Marks on Rocky Coasts. J. Ecol. 1949, 37, 289. [Google Scholar] [CrossRef]
  19. Menge, B.A.; Branch, G.M. Rocky Intertidal Communities. In Marine Community Ecology; Bertness, M.D., Gaines, S.D., Hay, M.E., Eds.; Sinauer Associates: Sunderland, MA, USA, 2001; pp. 221–251. [Google Scholar]
  20. Branch, G.M.; Moreno, C.A. Intertidal and subtidal grazers. In Rocky Shores: Exploitation in Chile and South Africa Ecological Studies; Siegfried, W.R., Ed.; Springer Verlag: Berlin, Germany, 1994; pp. 75–100. [Google Scholar]
  21. Aguilera, M.A. The functional roles of herbivores in the rocky intertidal systems in Chile: A review of food preferences and consumptive effects. Rev. Chil. de Hist. Nat. 2011, 84, 241–261. [Google Scholar] [CrossRef] [Green Version]
  22. Dethier, M.N.; Williams, S.L. Seasonal stresses shift optimal intertidal algal habitats. Mar. Biol. 2009, 156, 555–567. [Google Scholar] [CrossRef]
  23. Benincà, E.; Ballantine, B.; Ellner, S.P.; Huisman, J. Species fluctuations sustained by a cyclic succession at the edge of chaos. Proc. Natl. Acad. Sci. USA 2015, 112, 6389–6394. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Valdivia, N.; Aguilera, M.; Navarrete, S.; Broitman, B. Disentangling the effects of propagule supply and environmental filtering on the spatial structure of a rocky shore metacommunity. Mar. Ecol. Prog. Ser. 2015, 538, 67–79. [Google Scholar] [CrossRef] [Green Version]
  25. Kanamori, Y.; Fukaya, K.; Noda, T. Seasonal changes in community structure along a vertical gradient: Patterns and processes in rocky intertidal sessile assemblages. Popul. Ecol. 2017, 59, 301–313. [Google Scholar] [CrossRef]
  26. Suárez, J.L.; Hansen, M.; Urtubia, U.; Lenz, M.; Valdivia, N.; Thiel, M. Season-dependent effects of ocean warming on the physiological performance of a native and a non-native sea anemone. J. Exp. Mar. Biol. Ecol. 2019, 522, 151229. [Google Scholar] [CrossRef]
  27. Vellend, M. The Theory of Ecological Communities (MPB-57); Princeton University Press: Princeton, NJ, USA, 2016. [Google Scholar] [CrossRef]
  28. Couceiro, L.; Le Gac, M.; Hunsperger, H.M.; Mauger, S.; Destombe, C.; Cock, J.M.; Ahmed, S.; Coelho, S.M.; Valero, M.; Peters, A.F. Data from: Evolution and maintenance of haploid-diploid life cycles in natural populations: The case of the marine brown alga Ectocarpus. Evolution 2015, 69, 1808–1822. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Montecinos, A.E.; Couceiro, L.; Peters, A.F.; Desrut, A.; Valero, M.; Guillemin, M.L. Species delimitation and phylogeographic analyses in the Ectocarpus subgroup siliculosi (Ectocarpales, Phaeophyceae). J Phycol. 2017, 53, 17–31. [Google Scholar] [CrossRef]
  30. Schweikert, K.; Sutherland, J.E.; Burritt, D.J.; Hurd, C.L. Analysis of spatial and temporal diversity and distribution of porphyra (rhodophyta) in southeastern new zealand supported by the use of molecular tools1. J. Phycol. 2012, 48, 530–538. [Google Scholar] [CrossRef] [PubMed]
  31. Meynard, A.; Zapata, J.; Salas, N.; Betancourtt, C.; Pérez-Lara, G.; Castañeda, F.; Ramírez, M.E.; Contador, C.B.; Guillemin, M.L.; Porcia, L.C. Genetic and morphological differentiation of Porphyra and Pyropia species (Bangiales, Rhodophyta) coexisting in a rocky intertidal in Central Chile. J. Phycol. 2019, 55, 297–313. [Google Scholar] [CrossRef] [PubMed]
  32. Zapata, J.; Meynard, A.; Anguita, C.; Espinoza, C.; Alvear, P.; Kumar, M.; Porcia, L.C. Non-random distribution and ecophysiological differentiation of Pyropia species (Bangiales, Rhodophyta) through environmental gradients. J. Phycol. 2019, 55, 1140–1153. [Google Scholar] [CrossRef] [PubMed]
  33. Muangmai, N.; Preuss, M.; Zuccarello, G.C. Comparative physiological studies on the growth of cryptic species of Bostrychia intricata (Rhodomelaceae, Rhodophyta) in various salinity and temperature conditions. Phycol. Res. 2015, 63, 300–306. [Google Scholar] [CrossRef]
  34. West, A.L.; Mathieson, A.C.; Klein, A.S.; Neefus, C.D.; Bray, T.L. Molecular ecological studies of New England species of Porphyra (Rhodophyta, Bangiales). Nova Hedwig. 2005, 80, 1–24. [Google Scholar] [CrossRef]
  35. Yang, L.; Lu, Q.; Brodie, J. A review of the bladed Bangiales (Rhodophyta) in China: History, culture and taxonomy. Eur. J. Phycol. 2017, 52, 251–263. [Google Scholar] [CrossRef]
  36. Avila, M.; Seguel, M. An overview of seaweed resources in Chile. J. Appl. Phycol. 1993, 5, 133–139. [Google Scholar] [CrossRef]
  37. Yáñez, E.; Lagos, N.A.; Norambuena, R.; Silva, C.; Letelier, J.; Muck, K.-P.; Martin, G.S.; Benítez, S.; Broitman, B.R.; Contreras, H.; et al. Impacts of Climate Change on Marine Fisheries and Aquaculture in Chile. Clim. Change Impacts Fish. Aquac. A Glob. Anal. 2017, 1, 239–332. [Google Scholar] [CrossRef]
  38. Nelson, W.A.; Broom, J.E.S. The identity of Porphyra columbina (Bangiales, Rhodophyta) originally described from the New Zealand subantarctic islands. Aust. Syst. Bot. 2010, 23, 16–26. [Google Scholar] [CrossRef]
  39. Ramírez, M.E.; Contreras-Porcia, L.; Guillemin, M.-L.; Brodie, J.; Valdivia, C.; Flores-Molina, M.R.; Núñez, A.; Contador, C.B.; Lovazzano, C. Pyropia orbicularis sp. nov. (Rhodophyta, Bangiaceae) based on a population previously known as Porphyra columbina from the central coast of Chile. Phytotaxa 2014, 158, 133. [Google Scholar] [CrossRef]
  40. Contreras-Porcia, L.; Callejas, S.; Thomas, D.; Sordet, C.; Pohnert, G.; Contreras, A.; Lafuente, A.; Flores-Molina, M.R.; Correa, J.A. Seaweeds early development: Detrimental effects of desiccation and attenuation by algal extracts. Planta 2011, 235, 337–348. [Google Scholar] [CrossRef] [PubMed]
  41. Smith, C.M.; A Berry, J. Recovery of photosynthesis after exposure of intertidal algae to osmotic and temperature stresses: Comparative studies of species with differing distributional limits. Oecologia 1986, 70, 6–12. [Google Scholar] [CrossRef] [PubMed]
  42. Schonbeck, M.; Norton, T.A. The effects of brief periodic submergence on intertidal fucoid algae. Estuar. Coast. Mar. Sci. 1979, 8, 205–211. [Google Scholar] [CrossRef]
  43. Mansilla, A.; Werlinger, C.; Palacios, M.; Navarro, N.P.; Cuadra, P. Effects of UVB Radiation on the Initial Stages of Growth of Gigartina Skottsbergii, Sarcothalia Crispata and Mazzaella Laminarioides (Gigartinales, Rhodophyta). J. Appl. Phycol. 2006, 18, 451–459. [Google Scholar] [CrossRef]
  44. Gómez, I.; Huovinen, P. Morpho-functional patterns and zonation of South Chilean seaweeds: The importance of photosynthetic and bio-optical traits. Mar. Ecol. Prog. Ser. 2011, 422, 77–91. [Google Scholar] [CrossRef] [Green Version]
  45. Underwood, A.; Jernakoff, P. The effects of tidal height, wave-exposure, seasonality and rock-pools on grazing and the distribution of intertidal macroalgae in New South Wales. J. Exp. Mar. Biol. Ecol. 1984, 75, 71–96. [Google Scholar] [CrossRef]
  46. Schonbeck, M.W.; Norton, T.A. Factors controlling the lower limits of fucoid algae on the shore. J. Exp. Mar. Biol. Ecol. 1980, 43, 131–150. [Google Scholar] [CrossRef]
  47. Mangialajo, L.; Chiantore, M.; Susini, M.-L.; Meinesz, A.; Cattaneo-Vietti, R.; Thibaut, T. Zonation patterns and interspecific relationships of fucoids in microtidal environments. J. Exp. Mar. Biol. Ecol. 2012, 412, 72–80. [Google Scholar] [CrossRef]
  48. Kéfi, S.; Berlow, E.L.; Wieters, E.A.; Joppa, L.N.; Wood, S.A.; Brose, U.; Navarrete, S.A. Network structure beyond food webs: Mapping non-trophic and trophic interactions on Chilean rocky shores. Ecology 2015, 96, 291–303. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Moreno, C.A.; Jaramillo, E. The Role of Grazers in the Zonation of Intertidal Macroalgae of the Chilean Coast. Oikos 1983, 41, 73. [Google Scholar] [CrossRef]
  50. Tejada-Martinez, D.; López, D.N.; Bonta, C.C.; Sepúlveda, R.D.; Valdivia, N. Positive and negative effects of mesograzers on early-colonizing species in an intertidal rocky-shore community. Ecol. Evol. 2016, 6, 5761–5770. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Paine, R.T. Food Webs: Linkage, Interaction Strength and Community Infrastructure. J. Anim. Ecol. 1980, 49, 666. [Google Scholar] [CrossRef]
  52. Scrosati, R.A.; Heaven, C. Spatial trends in community richness, diversity, and evenness across rocky intertidal environmental stress gradients in eastern Canada. Mar. Ecol. Prog. Ser. 2007, 342, 1–14. [Google Scholar] [CrossRef] [Green Version]
  53. Valdivia, N.; Segovia-Rivera, V.; Fica, E.; Bonta, C.C.; Aguilera, M.A.; Broitman, B. Context-dependent functional dispersion across similar ranges of trait space covered by intertidal rocky shore communities. Ecol. Evol. 2017, 7, 1882–1891. [Google Scholar] [CrossRef] [PubMed]
  54. Muñoz-Muga, P.; Romo, H.; Calderón, C.; Evrard, O.; Díaz, H. Amplificación del gen rbcL revela primer registro de Porphyra mumfordii (Bangiales, Rhodophyta) en la bahía de Valparaíso, Chile central. Rev. Biol. Mar. Oceanogr. 2018, 53, 131–140. [Google Scholar] [CrossRef] [Green Version]
  55. Lindstrom, S.C.; Cole, K.M. The Porphyra lanceolata–P. pseudolanceolata (Bangiales, Rhodophyta) complex unmasked: Recognition of new species based on isozymes, morphology, chromosomes and distributions. Phycologia 1992, 31, 431–448. [Google Scholar] [CrossRef]
  56. Lindstrom, S.C.; Hughey, J.R.; Aguilar Rosas, L.E. Four new species of Pyropia (Bangiales, Rhodophyta) from the west coast of North America: The Pyropia lanceolata species complex updated. PhytoKeys. 2015, 52, 1–22. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Blouin, N.A.; Brodie, J.A.; Grossman, A.C.; Xu, P.; Brawley, S.H. Porphyra: A marine crop shaped by stress. Trends Plant Sci. 2011, 16, 29–37. [Google Scholar] [CrossRef] [PubMed]
  58. Contreras-Porcia, L.; Thomas, D.; Flores, V.; Correa, J.A. Tolerance to oxidative stress induced by desiccation in Porphyra columbina (Bangiales, Rhodophyta). J. Exp. Bot. 2010, 62, 1815–1829. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Contreras-Porcia, L.; López-Cristoffanini, C.; Meynard, A.; Kumar, M. Tolerance Pathways to Desiccation Stress in Seaweeds. In Systems Biology of Marine Ecosystems; Springer: Cham, Switzerland, 2017; pp. 13–33. [Google Scholar] [CrossRef]
  60. Kim, J.K.; Kraemer, G.P.; Yarish, C. Physiological activity of Porphyra in relation to eulittoral zonation. J. Exp. Mar. Biol. Ecol. 2008, 365, 75–85. [Google Scholar] [CrossRef]
  61. Kim, J.K.; Kraemer, G.P.; Yarish, C. Comparison of growth and nitrate uptake by New England Porphyra species from different tidal elevations in relation to desiccation. Phycol. Res. 2009, 57, 152–157. [Google Scholar] [CrossRef]
  62. Abe, S.; Kurashima, A.; Yokohama, Y.; Tanaka, J. The Cellular Ability of Desiccation Tolerance in Japanese Intertidal Seaweeds. Bot. Mar. 2001, 44, 125–131. [Google Scholar] [CrossRef]
  63. Boedeker, C.; Farr, T.J.; Nelson, W.A. Comparative culture experiments with filamentous members of the Bangiales (Rhodophyta) from New Zealand: Insights into ecologic adaptation and biogeography. Phycol. Res. 2008, 56, 183–192. [Google Scholar] [CrossRef]
  64. Royer, C.J.; Redmond, S.; Lai, C.S.; Brawley, S.H. Porphyra umbilicalis in applied and basic research: Reproductive phenology, development, seed stock culture, and a field trial for aquaculture. J. Appl. Phycol. 2018, 31, 547–560. [Google Scholar] [CrossRef]
  65. Holmes, M.J.; Brodie, J. Morphology, seasonal phenology and observations on some aspects of the life history in culture of Porphyra dioica (Bangiales, Rhodophyta) from Devon, UK. Phycologia 2004, 43, 176–188. [Google Scholar] [CrossRef]
  66. Broom, J.E.; Nelson, W.A.; Yarish, C.; Jones, W.A.; Rosas, R.A.; Rosas, L.E.A. A reassessment of the taxonomic status of Porphyra suborbiculata, Porphyra carolinensis and Porphyra lilliputiana (Bangiales, Rhodophyta) based on molecular and morphological data. Eur. J. Phycol. 2002, 37, 227–235. [Google Scholar] [CrossRef]
  67. Lindstrom, S.C. Cryptic diversity, biogeography and genetic variation in Northeast Pacific species of Porphyra sensu lato (Bangiales, Rhodophyta). J. Appl. Phycol. 2008, 20, 951–962. [Google Scholar] [CrossRef]
  68. Sutherland, J.E.; Lindstrom, S.C.; Nelson, W.A.; Brodie, J.; Lynch, M.D.J.; Hwang, M.S.; Choi, H.-G.; Miyata, M.; Kikuchi, N.; Oliveira, M.C.; et al. A new look at an ancient order: Generic revision of the bangiales (rhodophyta)1. J. Phycol. 2011, 47, 1131–1151. [Google Scholar] [CrossRef] [PubMed]
  69. Kucera, H.; Saunders, G.W. A a survey of bangiales (rhodophyta) based on multiple molecular markers reveals cryptic diversity1. J. Phycol. 2012, 48, 869–882. [Google Scholar] [CrossRef] [PubMed]
  70. Egea, E.; David, B.; Choné, T.; Laurin, B.; Féral, J.P.; Chenuil, A. Morphological and genetic analyses reveal a cryptic species complex in the echinoid Echinocardium cordatum and rule out a stabilizing selection explanation. Mol. Phylogenetics Evol. 2016, 94, 207–220. [Google Scholar] [CrossRef] [PubMed]
  71. Pfingstl, T.; Lienhard, A.; Baumann, J.; Koblmüller, S. A taxonomist‘s nightmare – Cryptic diversity in Caribbean intertidal arthropods (Arachnida, Acari, Oribatida). Mol. Phylogenet. Evol. 2021, 163, 107240. [Google Scholar] [CrossRef] [PubMed]
  72. Calderon, M.S.; Bustamante, D.E.; Boo, S.M. Red algal diversity (Rhodophyta) from Peru based on molecular analysis. Phytotaxa 2020, 454, 1–23. [Google Scholar] [CrossRef]
  73. Stegenga, H.; Bolton, J.J.; Anderson, R.J. Seaweeds of the South African West Coast. In Contributions from the Bolus Herbarium 18; Hall, A.V., Ed.; Bolus Herbarium, University of Cape Town: Cape Town, South Africa, 1997; pp. 1–655. [Google Scholar]
  74. Jones, W.A.; Griffin, N.J.; Jones, D.T.; Nelson, W.A.; Farr, T.J.; Broom, J.E. Phylogenetic diversity in South African Porphyra (Bangiales, Rhodophyta) determined by nuclear SSU sequence analyses. Eur. J. Phycol. 2004, 39, 197–211. [Google Scholar] [CrossRef]
  75. Reddy, M.M.; De Clerck, O.; Leliaert, F.; Anderson, R.J.; Bolton, J.J. An appraisal of the genus Pyropia (Bangiales, Rhodophyta) from southern Africa based on a multi-gene phylogeny, morphology and ecology, including the description of Pyropia meridionalis sp. nov. South Afr. J. Bot. 2020, 131, 18–32. [Google Scholar] [CrossRef]
  76. Reddy, M.M. Taxonomy and Systematics of the Bangiales (Rhodophyta) in South Africa using an Integrative Approach. Ph.D. Thesis, University of Cape Town, Cape Town, South Africa, 2018. [Google Scholar]
  77. Hommersand, M.H.; Fredericq, S. Biogeography of the marine red algae of the South African West Coast: A molecular approach. In Proceedings of the XVIIth International Seaweed Symposium, Cape Town, South Africa, 28 January–2 February 2001; pp. 325–336. [Google Scholar]
  78. Fraser, C.I.; Nikula, R.; Spencer, H.G.; Waters, J.M. Kelp genes reveal effects of subantarctic sea ice during the Last Glacial Maximum. Proc. Natl. Acad. Sci. USA 2009, 106, 3249–3253. [Google Scholar] [CrossRef] [Green Version]
  79. Broom, J.E.S.; Nelson, W.A.; Farr, T.J.; Phillips, L.E.; Clayton, M. Relationships of the Porphyra (Bangiales, Rhodophyta) flora of the Falkland Islands: A molecular survey using rbcL and nSSU sequence data. Aust. Syst. Bot. 2010, 23, 27–37. [Google Scholar] [CrossRef]
  80. Tyberghein, L.; Verbruggen, H.; Pauly, K.; Troupin, C.; Mineur, F.; De Clerck, O. Bio-ORACLE: A global environmental dataset for marine species distribution modelling. Glob. Ecol. Biogeogr. 2011, 21, 272–281. [Google Scholar] [CrossRef]
  81. Assis, J.; Tyberghein, L.; Bosch, S.; Verbruggen, H.; Serrão, E.A.; De Clerck, O.; Tittensor, D. Bio-ORACLE v2.0: Extending marine data layers for bioclimatic modelling. Glob. Ecol. Biogeogr. 2017, 27, 277–284. [Google Scholar] [CrossRef]
  82. Atkinson, L.P.; Valle-Levinson, A.; Figueroa, D.; De Pol-Holz, R.; Gallardo, V.A.; Schneider, W.; Blanco, J.L.; Schmidt, M. Oceanographic observations in Chilean coastal waters between Valdivia and Concepción. J. Geophys. Res. Earth Surf. 2002, 107, 18-1–18-13. [Google Scholar] [CrossRef]
  83. Pérez, C.A.; Lagos, N.A.; Saldías, G.S.; Waldbusser, G.; Vargas, C.A. Riverine discharges impact physiological traits and carbon sources for shell carbonate in the marine intertidal mussel Perumytilus purpuratus. Limnol. Oceanogr. 2016, 61, 969–983. [Google Scholar] [CrossRef] [Green Version]
  84. Moreno, C.A. Community patterns generated by human harvesting on Chilean shores: A review. Aquat. Conserv. Mar. Freshw. Ecosyst. 2001, 11, 19–30. [Google Scholar] [CrossRef]
  85. Häussermann, V.; Försterra, G. Marine Benthic Fauna of Chilean Patagonia; Nature in Focus: Puerto Montt, Chile, 2009; p. 1000. [Google Scholar]
  86. Velásquez, C.; Jaramillo, E.; A Camus, P.; Manzano, M.; Sánchez, R. Biota del intermareal rocoso expuesto de la Isla Grande de Chiloé, Archipiélago de Chiloé, Chile: Patrones de diversidad e implicancias ecológicas y biogeográficas. Rev. De Biol. Mar. Oceanogr. 2016, 51, 33–50. [Google Scholar] [CrossRef] [Green Version]
  87. Jara, F.H.; Moreno, C.A. Herbivory and Structure in a Midlittorial Rocky Community: A Case in Southern Chile. Ecol. Soc. Am. 1984, 65, 28–38. [Google Scholar] [CrossRef]
  88. Saunders, G.W. Applying DNA barcoding to red macroalgae: A preliminary appraisal holds promise for future applications. Philos. Trans. R. Soc. B Biol. Sci. 2005, 360, 1879–1888. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  89. Anderson, M.J. A new method for non-parametric multivariate analysis of variance. Austral. Ecol. 2001, 26, 32–46. [Google Scholar] [CrossRef]
  90. Borcard, D.; Gillet, F.; Legendre, P. Numerical Ecology with R; Springer: New York, NY, USA, 2011. [Google Scholar]
  91. Johnson, P.C.D. Extension of Nakagawa & Schielzeth’s R2GLMM to random slopes models. Methods Ecol. Evol. 2014, 5, 944–946. [Google Scholar] [CrossRef] [Green Version]
  92. Nakagawa, S.; Johnson, P.C.D.; Schielzeth, H. The coefficient of determination R 2 and intra-class correlation coefficient from generalized linear mixed-effects models revisited and expanded. J. R. Soc. Interface 2017, 14, 20170213. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  93. R Core Team. R: A Language and Environment for Statistical Computing. Available online: https://www.r-project.org/ (accessed on 16 August 2021).
  94. Neuwirth, E. RColorBrewer: ColorBrewer Palettes. Package version 1. Available online: https://CRAN.R-project.org/package=RColorBrewer (accessed on 16 August 2021).
  95. Wilke, C.O. Cowplot: Streamlined Plot Theme and Plot Annotations for ggplot2. R Package Version 100. Available online: https://cran.r-project.org/package=cowplot (accessed on 16 August 2021).
  96. Brooks, M.E.; Kristensen, K.; van Benthem, K.J.; Magnusson, A.; Berg, C.W.; Nielsen, A.; Skaug, H.J.; Mächler, M.; Bolker, B.M. glmmTMB balances speed and flexibility among packages for zero-inflated generalized linear mixed modeling. R. J. 2017, 9, 378–400. [Google Scholar] [CrossRef] [Green Version]
  97. Ram, K.; Wickham, H. Wesanderson: A Wes Anderson Palette Generator. R Package Version 036. Available online: https://rdrr.io/cran/wesanderson (accessed on 16 August 2021).
  98. Oksanen, J.; Blanchet, F.G.; Friendly, M.; Kindt, R.; Legendre, P.; McGlinn, D.; Minchin, P.R.; O’Hara, R.B.; Simpson, G.L.; Solymos, P.; et al. Vegan: Community Ecology Package. R Packag Version 2. Available online: https://CRAN.R-project.org/package=vegan (accessed on 16 August 2021).
  99. Wickham, H.; Averick, M.; Bryan, J.; Chang, W.; McGowan, L.D.A.; François, R.; Grolemund, G.; Hayes, A.; Henry, L.; Hester, J.; et al. Welcome to the Tidyverse. J. Open Source Softw. 2019, 4, 1686. [Google Scholar] [CrossRef]
  100. Bartoń, K. MuMIn: Multi-Model Inference. R package version 1.43.17. Available online: https://CRAN.R-project.org/package=MuMIn (accessed on 16 August 2021).
  101. Luo, Q.; Zhu, Z.; Zhu, Z.; Yang, R.; Qian, F.; Chen, H.; Yan, X. Different Responses to Heat Shock Stress Revealed Heteromorphic Adaptation Strategy of Pyropia haitanensis (Bangiales, Rhodophyta). PLoS ONE 2014, 9, e94354. [Google Scholar] [CrossRef]
  102. Watanabe, Y.; Yamada, H.; Mine, T.; Kawamura, Y.; Nishihara, G.N.; Terada, R. Photosynthetic responses of Pyropia yezoensis f. narawaensis (Bangiales, Rhodophyta) to a thermal and PAR gradient vary with the life-history stage. Phycologia 2016, 55, 665–672. [Google Scholar] [CrossRef]
  103. Santelices, B.; Ugarte, R. Algal life-history strategies and resistance to digestion. Mar. Ecol. Prog. Ser. 1987, 35, 267–275. [Google Scholar] [CrossRef]
  104. Steneck, R.S.; Dethier, M.N. A Functional Group Approach to the Structure of Algal-Dominated Communities. Oikos 1994, 69, 476. [Google Scholar] [CrossRef] [Green Version]
Plants 11 00605 g001 550
Figure 1. (A) Total number of individuals identified through COI sequencing to one of the nine Bangiales species detected in our study area. (B) Max length (mm ± SD) and max width (mm ± SD) for each of the nine Bangiales species studied. Number of individuals measured per species: Pyropia sp. CHJ (n = 356), Pyropia orbicularis (n = 224), Porphyra mumfordii (n = 204), Porphyra sp. FIH (n = 63), Pyropia sp. CHH (n = 60), Porphyra longissima (n = 21), Porphyra luchea (n = 7), Fuscifolium sp. CHA (n = 4), and Pyropia saldanhae (n = 3). Only complete wet thalli were measured. See the text for more details.
Figure 1. (A) Total number of individuals identified through COI sequencing to one of the nine Bangiales species detected in our study area. (B) Max length (mm ± SD) and max width (mm ± SD) for each of the nine Bangiales species studied. Number of individuals measured per species: Pyropia sp. CHJ (n = 356), Pyropia orbicularis (n = 224), Porphyra mumfordii (n = 204), Porphyra sp. FIH (n = 63), Pyropia sp. CHH (n = 60), Porphyra longissima (n = 21), Porphyra luchea (n = 7), Fuscifolium sp. CHA (n = 4), and Pyropia saldanhae (n = 3). Only complete wet thalli were measured. See the text for more details.
Plants 11 00605 g001
Plants 11 00605 g002 550
Figure 2. Principal coordinates analysis (PCoA) showing changes in species composition depending on the elevation and season. The first two axes explain more than 89% of the total variance.
Figure 2. Principal coordinates analysis (PCoA) showing changes in species composition depending on the elevation and season. The first two axes explain more than 89% of the total variance.
Plants 11 00605 g002
Plants 11 00605 g003 550
Figure 3. Proportion of Bangiales species identified at each intertidal elevation level (low, mid, high, and spray) during the four seasons. Bar width is proportional to number of samples at each intertidal elevation.
Figure 3. Proportion of Bangiales species identified at each intertidal elevation level (low, mid, high, and spray) during the four seasons. Bar width is proportional to number of samples at each intertidal elevation.
Plants 11 00605 g003
Plants 11 00605 g004 550
Figure 4. Location of the five sampling sites with Pilolcura (pink), Niebla (orange), San Carlos (light blue), and Los Liles (green) situated at the mouth of the Valdivia River estuary and Melinka (dark blue) in Ascension Island.
Figure 4. Location of the five sampling sites with Pilolcura (pink), Niebla (orange), San Carlos (light blue), and Los Liles (green) situated at the mouth of the Valdivia River estuary and Melinka (dark blue) in Ascension Island.
Plants 11 00605 g004
Plants 11 00605 g005 550
Figure 5. Intertidal zonation of bladed Bangiales sampling. The four levels that were sampled (“splash zone”, “high”, “mid”, and “low”) are indicated. Capital letters correspond to samples encompassing 5 nearby specimens collected in a 25 cm radius. See the text for more details.
Figure 5. Intertidal zonation of bladed Bangiales sampling. The four levels that were sampled (“splash zone”, “high”, “mid”, and “low”) are indicated. Capital letters correspond to samples encompassing 5 nearby specimens collected in a 25 cm radius. See the text for more details.
Plants 11 00605 g005
Table
Table 1. Raw counts of the species identified in the five sites sampled for each season.
Table 1. Raw counts of the species identified in the five sites sampled for each season.
Pyropia sp. CHJPyropia orbicularisPorphyra mumfordiiPorphyra sp. FIHPyropia sp. CHHPorphyra longissimaPorphyra lucheaFuscifolium sp. CHAPyropia saldanhaeƩ
AutumnPilolcura 28232274000084
Niebla4811802300072
San Carlos15342180000078
Los Liles44323190000089
Melinka16132130000044
WinterPilolcura 4220071300073
Niebla39183101300092
San Carlos5351540700185
Los Liles5523080000095
Melinka10304530060094
SpringPilolcura 29480111000089
Niebla55110021500092
San Carlos33282700201293
Los Liles40825170001091
Melinka0583720000097
SummerPilolcura 1276011011092
Niebla23251035000084
San Carlos4647200000095
Los Liles312114200101088
Melinka25642000000100
Ʃ57952335117066247431727
Table
Table 2. Results of the PERMANOVA on the effects of intertidal elevation and season on Bangiales species composition. Pairwise tests were performed within each season with the low-tide elevation used as the reference group for the comparisons between the intertidal elevations.
Table 2. Results of the PERMANOVA on the effects of intertidal elevation and season on Bangiales species composition. Pairwise tests were performed within each season with the low-tide elevation used as the reference group for the comparisons between the intertidal elevations.
EffectdfFR2
Elevation342.91 **0.19
Season315.82 **0.07
Elevation: season93.59 **0.05
Residuals 464 0.68
Total479 1
Pairwise test
SeasonIntertidal elevation
SpringMid **
SpringHigh **
SpringSpray **
SummerMid
SummerHigh **
SummerSpray **
Autumn Mid **
AutumnHigh **
AutumnSpray **
WinterMid **
WinterHigh **
WinterSpray **
df = degrees of freedom; F = F ratio; R2 = coefficient of determination; Asterisks denote statistical significance: ** p < 0.01.
Table
Table 3. Tukey-corrected post-hoc comparisons after the zero-inflated hierarchical generalized linear models. All comparisons were conducted with a “treatment” contrast with low elevation as the reference group and 462 degrees of freedom. The t-ratio of each comparison is shown.
Table 3. Tukey-corrected post-hoc comparisons after the zero-inflated hierarchical generalized linear models. All comparisons were conducted with a “treatment” contrast with low elevation as the reference group and 462 degrees of freedom. The t-ratio of each comparison is shown.
SeasonIntertidal Elevation Pyropia sp. CHJPyropia orbicularisPorphyra mumfordiiPorphyra sp. FIHPyropia sp. CHH
SpringMid2.40−1.5<−0.01<−0.01−1.5
SpringHigh3.16 **−6.2 ***<−0.01<−0.01−2.5
SpringSpray−1.31−4.3 ***<−0.01<−0.01<−0.01
SummerMid2.3−0.31<−0.01<−0.01−0.9
SummerHigh3.49 **−4.9 ***3.9 **<−0.01<−0.01
SummerSpray0.04<−0.011.5<−0.01<−0.01
AutumnMid1.25−4.1 ***1.61.7<−0.01
AutumnHigh−0.35<−0.010.93.9 **<−0.01
AutumnSpray−4.41 ***<−0.01−0.92.7 *<−0.01
Winter Mid1.24−3.8 **2.41.2<−0.01
WinterHigh−0.33−4.5 ***5.0 ***1.6<−0.01
WinterSpray−4.31 ***<−0.012.7 *1.9<−0.01
Asterisks denote statistical significance: * p < 0.05; ** p < 0.01; *** p < 0.001.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cid Alda, F.P.; Valdivia, N.; Guillemin, M.-L. More than What Meets the Eye: Differential Spatiotemporal Distribution of Cryptic Intertidal Bangiales. Plants 2022, 11, 605. https://doi.org/10.3390/plants11050605

AMA Style

Cid Alda FP, Valdivia N, Guillemin M-L. More than What Meets the Eye: Differential Spatiotemporal Distribution of Cryptic Intertidal Bangiales. Plants. 2022; 11(5):605. https://doi.org/10.3390/plants11050605

Chicago/Turabian Style

Cid Alda, Fernanda P., Nelson Valdivia, and Marie-Laure Guillemin. 2022. "More than What Meets the Eye: Differential Spatiotemporal Distribution of Cryptic Intertidal Bangiales" Plants 11, no. 5: 605. https://doi.org/10.3390/plants11050605

APA Style

Cid Alda, F. P., Valdivia, N., & Guillemin, M. -L. (2022). More than What Meets the Eye: Differential Spatiotemporal Distribution of Cryptic Intertidal Bangiales. Plants, 11(5), 605. https://doi.org/10.3390/plants11050605

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop