Taxonomic distinctness and diversity patterns of a polychaete (Annelida) community on the continental shelf of the Southern Gulf of Mexico

The spatial patterns of taxonomic diversity of annelid polychaete species from the continental shelf in the Southern Gulf of Mexico were examined in this study. We used taxonomic distinctness and its spatial variations to explore the diversity patterns and how they change between Southern Gulf of Mexico regions. In addition, using taxonomic distinctness as a dissimilarity measure and Ward’s Clustering, we characterized three distinct faunal assemblages. We also investigated patterns of richness, taxonomic distinctness, and distance decay of similarity between sampling stations as a ß-diversity measure. Finally, we examined the spatial relationships between polychaete assemblages and environmental variables to test the relative importance of spatial and environmental components in annelid polychaete community structure from the Southern Gulf of Mexico. We used a combination of eigenvector-based multivariate analyses (dbMEMs) and distance-based redundancy analysis (dbRDA) to quantify the relative importance of these explanatory variables on the spatial variations of taxonomic distinctness. The significance level of spatial and environmental components to the distribution of polychaete species showed that the combined effect of spatial processes and sediment characteristics explained a higher percentage of the variance than those parameters could alone.

This statement is required for submission and will appear in the published article if the submission is accepted.Please make sure it is accurate.

Unfunded studies
Enter: The author(s) received no specific funding for this work.General guidance is provided below.
Consult the submission guidelines for detailed instructions.Make sure that all information entered here is included in the Methods section of the manuscript.
Samples were collected during six oceanographic expeditions, on board the R/V "Justo Sierra" (Universidad Nacional Autónoma de México), from 1988 to 1996, as part of the institutional projects "IMCA" and "DINAMO" and contact information or URL).This text is appropriate if the data are owned by a third party and authors do not have permission to share the data.

Introduction
Explaining the causes of biodiversity variability at various scales and correlating these variations to changes in environmental conditions are two of ecology's primary goals (Benedetti-Cecchi et al., 2010;Economo & Keitt, 2008).To estimate how biodiversity may respond to climate change, it is essential to understand the patterns of species distribution and the ecological and evolutionary processes that lead to these patterns.Abiotic factors influence the borders of the species ranges based on their physiological preferences (Chase & Leibold, 2003).Numerous complex mechanisms have determined the distribution of species at diverse geographical and temporal scales (Leibold et al., 2004).
Over the past few years, there has been a remarkable growth in research on marine biodiversity patterns, and understanding the spatial distribution of species in soft-bottom habitats is essential to comprehend how species interact with one another and with the environment (Dauvin & Bellan-Santini, 2004) and how benthic communities and ecosystems function (Checon & Amaral, 2017;Legendre et al., 2005).Community ecology studies frequently address the connection between species composition and environmental factors (Alves et al., 2020;Giangrande & Gravina, 2015;Mikac et al., 2020;Quiroz-Martínez et al., 2021;Semprucci et al., 2019;Thrush et al., 2010).Spatially structured environmental variation plays a vital role in species distribution (Dray et al., 2012); consequently, discriminating spatial ecological structures in species composition, distribution, and diversity is a pivotal question in community ecology (Legendre, 1993).However, the spatial scale can also incorporate puzzling factors; thus, statistical models should include them as predictors or covariable parameters (Dray et al., 2006).The local scale is generally considered the spatial autocorrelation created by community dynamics, while the regional scale corresponds to the scale of the environmental drivers (Dray et al., 2012).
The spatial component in community analysis has been incorporated using a variety of methodologies (Dormann et al., 2007;Keitt et al., 2002).For example, ordination analysis (Borcard et al., 1992;Legendre, 2018) can partition the variance across different variables.This approach entails assessing environmental variables, creating a spatial coordinates matrix, and characterizing these elements' individual and combined contributions to the community structure.
Due to its ability to explain the regional and local dynamics that shape communities, variation partitioning has become a crucial exploratory tool (Checon & Amaral, 2017;Gilbert & Bennett, 2010;Heino et al., 2012).Processes associated with sediment properties, temperature, salinity, dissolved oxygen, and nutrient content are among the drivers that structure the geographical distribution of biodiversity (Levin et al., 2001;Quiroz-Martínez et al., 2021;Schnurr et al., 2018).
Annelid polychaetes are among the most diverse and characteristic groups of benthic macroinvertebrates (Hutchings, 1998;Knox, 1977) and one of the richest concerning species numbers (Arvanitidis et al., 2002;Capa & Hutchings, 2021); it includes nearly 11,500 accepted nominal species and is often the dominant group in benthic macrofauna (Pamungkas et al., 2019), representing up to 60% of the species and over 70% of the abundance, mainly in soft-bottom communities (Dogan et al., 2005;Knox, 1977).They are ubiquitous in soft sediments and show significant spatial distribution patterns under the influence of environmental gradients (Gopal et al., 2020).Polychaetes can be used as "indicators" of different ecological conditions (Gambi & Giangrande, 1986;Quiroz-Martínez et al., 2012), and they can be particularly informative in assessing the health of benthic environments.Therefore, the analysis of the diversity patterns of this fauna can be crucial in understanding the whole ecosystem functioning since the biological processes that determine these patterns can reflect those of the entire ecosystem (Belan, 2003;Ergen et al., 2006;Giangrande et al., 2005;Glasby & Read, 1998;Mackie et al., 1997;Mutlu et al., 2010;Olsgard & Somerfield, 2000).In the continental shelf of the Gulf of Mexico, the polychaetes have been the focus of comprehensive research for decades as they represent a crucial group in abundance and diversity (Fauchald et al., 2009); however, most of these studies only present taxonomic lists and are often presented in a fragmented or isolated way; thus, information about their distribution patterns related to environmental factors is still scarce.A comprehensive polychaete database was therefore built using the results of several research projects examining the occurrence of this group in the continental shelf of the Southern Gulf of Mexico.
Using distance-based redundancy analysis (dbRDA) and variation partitioning, this study aims to analyze the environmental factors that influence polychaete diversity at small and intermediate scales along the continental shelf of the Southern Gulf of Mexico: 1) assessing the relative impact of spatial and environmental components on polychaete community structure; and 2) examining the spatial relationships between polychaete assemblages and environmental variables.A common approach to examining spatial ß-diversity is considering a directional turnover in community composition with distance (Distance Decay).In this framework, the similarity amongst pairs of biological communities typically decreases (decays) with increasing spatial and environmental distance.

Study area
The Gulf of Mexico is a marginal sea shared by three countries, Mexico, Cuba, and the United States of America.It is recognized as a large marine ecosystem (Dorr et al., 2019), where activities with significant economic importance, such as fisheries, transportation, and gas and oil extraction, constantly increase (Dorr et al., 2019).The Southern Gulf of Mexico (18° 30′22° 20′ N; 89° 41′-97° 4′ W), delimited by an imaginary line that goes from the port of Tampico, in the State of Tamaulipas to Progreso in the State of Yucatan covers approximately 46 000 km 2 .Its southeastern portion (Fig. 1, upper panel) is composed of two distinct regions: The Campeche Bank, an expansive shelf characterized by a gentle slope and uneven bottoms and sandbanks, coral reefs, and native biogenic and autogenic sediments along most of the Yucatan Peninsula's coast, is to the east of the Bay of Campeche, which is a narrow shelf with its upper limit at 130 m depth and between 45 and 65 km wide (Fig. 1) (Antoine, 1972;Martin & Bouma, 1978).Three seasons characterize the southern Gulf of Mexico, the dry season from March to May, the rainy season from June to October, and the winter storms season, locally called the "Nortes" season, from November to February (López-Cabello et al., 2022).The southern Gulf of Mexico has combined influences from river run-off to the west and winter storms ("Nortes") (Fuentes-Yaco et al., 2001;Hernández-Arana et al., 2005).The study area has been classified according to the sediments' carbonate content: the "terrigenous" region to the west with less than 25% and the "carbonate" region to the east and north, with more than 75% with a transitional carbonate-terrigenous region, between the first two, whose limits change throughout the year, with intermediate values (Carranza-Edwards et al., 1993;Hernández-Arana et al., 2003;Ortiz-Zamora et al., 2002).The terrigenous sediments are predominantly silty with variable amounts of gravel, sand, and clay, while fragments of shells, corals, and algae mainly constitute the carbonate ones.The Southern Gulf of Mexico is a complex tropical domain with intense river discharges to the west from the Grijalva-Usumacinta River system, which modifies the temperature and salinity, inducing coastal frontogenetic processes, in addition to being the primary source of terrigenous sediments (Carranza-Edwards et al., 1993;Monreal-Gómez et al., 1992).Discharges from the Grijalva-Usumacinta and San Pedro-San Pablo rivers, "Nortes", summer rains, and intrusion of oceanic water from the Caribbean Sea are the physical driving processes on the continental shelf of the southern Gulf of Mexico (Hernández-Arana et al., 2005;Monreal-Gómez et al., 1992).The Caribbean current drives the water circulation pattern in the spring and summer, with a southsouthwest flow; however, in the autumn and winter, the flow shifts to an east-northeast direction (Boicourt et al., 1998;Hernandez Arana et al., 2005;Salas-Monreal et al., 2018)

Data source
The polychaete database, constructed from extensive studies carried out on the continental shelf (< 200 m depth) of the Southern Gulf of Mexico (18° 30′-22° 20′ N; 89° 41′-97° 4′ W) was the primary source of information: the samples were collected during six oceanographic expeditions, on board the R/V "Justo Sierra" (Universidad Nacional Autónoma de México), from 1988 to 1996, as part of two institutional projects "IMCA" and "DINAMO" (Solís-Weiss et al., 1995).The collected specimens from 61 stations were identified to species level and deposited in the "Colección Nacional de Anélidos Poliquetos" from the Instituto de Ciencias del Mar y Limnología, Universidad Nacional Autónoma de México (CNAP-ICML, UNAM; DFE.IN.061.0598).The database includes information on all identified polychaete species, the geographical position of the sampling stations (latitude and longitude), and environmental variables such as sediment composition, depth, temperature, and salinity where they were collected.In addition, the taxonomic information was carefully checked.We verified the validity of names and synonymies and omitted species whose taxonomic identification was doubtful.We used the more recent systematic reviews and the World Polychaeta Database (Read & Fauchald, 2021) accessed through the World Register of Marine Species (WORMS) database.

Statistical analyses
The information in the database was linked to the corresponding supra-generic levels for every species; their relationships to genus, family, and class were included.We then calculated two indices to compare polychaete diversity between sampling stations: 1) the Average taxonomic distinctness (Δ + ), and 2) the variation in taxonomic distinctness (Λ + ), as defined by Warwick & Clarke (2001) and Clarke & Warwick (1998).
The Average taxonomic distinctness (Δ + ) was calculated using the function: where ωij is the phylogenetic/taxonomic path length between species i and j, and S is the number of species).
While the variation in taxonomic distinctness (Λ+) was calculated as follows: The average taxonomic distinctness (Δ + ) measures the average taxonomic distance between different species in an assemblage; the higher the values of Δ + , the more significant the average taxonomic difference between species in the assemblage.Conversely, Λ + is the variation in branch lengths amongst all pairs of species and measures the distribution of branch lengths within a sample; the greater the variation in branch lengths, the greater the value of Λ + .
To assess the differences in taxonomic distinctness (Δ + ) from expected Δ + values determined from the species list of polychaetes from the Gulf of Mexico, we used a randomization approach (as suggested by Clarke & Warwick (1998, 2001) and by Warwick & Clarke (2001)).A simulated distribution was created, leading to a theoretical mean (a horizontal line displaying the taxonomic distinctness for all polychaete species shown in the graph of Δ + or Λ + for all the sampling stations against richness in each station) and to a confidence funnel for each, Δ + , and Λ + , from random subsamples of the polychaete species as suggested by Bhat and Magurran (2006).
The simulation creates numerous random subsets of species, each of size m, from the entire species list, computes the corresponding Δ + and Λ + values, and establishes a range in which 95% of these values are included.Values of Δ + and Λ + located within the 95% probability funnel indicate that species diversity in the corresponding areas falls within the expected range, thus allowing both for sample size and effort-free diversity comparisons.
To evaluate the effects of geographical distance on taxonomic composition (Distance Decay), we calculated the distance in km between all pairs of sampling stations.Next, we used the geographical coordinates to calculate the distance between all pairs of sampling stations.Next, the relationship between pairwise faunal dissimilarity and spatial distance was assessed by fitting an exponential function describing the increase in faunal dissimilarity with spatial distance (Nekola & McGill, 2014;Nekola & White, 1999); the R package "betapart" (functions "decay.model"and "boot.coefs")was used; as it adjusts a General Linear Model (GLM) with dissimilarity as the response variable, spatial distance as the predictor, log link, and Gaussian error.Finally, the intercept and slope parameters were bootstrapped (1000 replicates) (Baselga et al., 2013;Gómez-Rodríguez & Baselga, 2018).
The indices were computed using the polychaete species' taxonomic hierarchy based on the Linnaean classification using the Plymouth Routines in Multivariate Ecological Research PRIMER-E (Clarke et al., 2014;Clarke & Gorley, 2015).The resulting matrices were examined to derive dissimilarity patterns through Ward's hierarchical cluster analysis on the taxonomic dissimilarity matrix.The optimal number of clusters was determined using the average silhouette width and the diagnostic group criteria suggested by Borcard et al. (2018).The cluster analysis was run on the θ + taxonomic dissimilarity index, a presence/absence "beta diversity" coefficient (Clarke et al., 2014;Clarke et al., 2006).The significant differences in the assemblage structure between the two studied regions were tested with a one-way analysis of similarity (ANOSIM).
In addition, we calculated the geographic distance between each pair of stations using the

Taxonomic Distinctness and Distance Decay
One hundred and seventy-three species from 40 families were recorded from the continental shelf of the Southern Gulf of Mexico.The eastern continental shelf region (Campeche Bank) had the highest number of species (164).In contrast, 125 species were collected in the western region of the study area (Bay of Campeche).The estimated Δ + values resulting from the polychaete species list in each of the 61 sampling stations from the Southern Gulf of Mexico showed that most of the calculated Δ + values were close to the expected simulated mode of the funnel (Fig. 2).According to these Δ + values, most of the stations were as diverse as expected since they were located inside the funnel.
The variation in taxonomic distinctness (Λ + ) showed that most stations were placed inside the funnel, close to the expected values.However, a few stations fell above the simulated funnel's upper limit, indicating higher-than-expected Λ + values (Fig. 2).The exponential fit between the calculated value of θ + taxonomic dissimilarity index as a function of the distance between all pairs of sampling stations showed an exponential relationship between compositional dissimilarity and distance.The distance among stations exhibited high variability, but most differed at less than 200 km.The dissimilarity between the polychaete taxonomic composition from all stations exponentially increased with increasing geographical distance (Figure 3).

Cluster Analysis
Ward's Clustering, using the θ + taxonomic dissimilarity index relating to the polychaete composition of the 61 sampling stations, showed three main faunal assemblages with high taxonomic dissimilarity (>0.75) between them (Figure 4).The ANOSIM analysis confirmed that spatial differences between the three polychaete assemblages were significant (RANOSIM(C1-C2)= 0.374, p= 0.001; RANOSIM(C1-C3)= 0.45, p= 0.001; RANOSIM(C2-C3)= 0.508, p= 0.001).The  The distance-based redundancy analysis (dbRDA) showed that the first axis explained 43.56% of the total variation, and the second explained 24.08% (Figure 5).The first axis represented a large-scale (with a high correlation with dbMEM1) sedimentary and thermal gradient characterized by increasing sand content and temperature.The "Carbonate" assemblage (blue) in the eastern region was associated with warm water and sandy sediments.Stations belonging to the "Terrigenous" assemblage (red) in the western study area are associated with a medium-scale (dbMEM7) salinity gradient.The "Transitional" assemblage (green) in front of Laguna de Términos was related to salinity variability and increasing mud content in the sediment (Figure 5).
All the environmental variables seem to be spatially structured.The unexplained variation may be linked to describers not used in the present study or simply due to random factors.The set of environmental (X1) and spatial variables (dbMEM1, X3) explains 5% of the variation in taxonomic dissimilarity across the study area.The environmental variables alone (X1 in the partitioning results) explain almost 10% of the variation, of which 5% is not spatially structured.This fraction represents species-environment relationships associated with local environmental conditions.Medium-scale spatial processes (dbMEM7, X4) explained 3.4% of the taxonomic variability.Finally, the linear spatial structure (x and y coordinates; X2) contributions were less significant (1.5%) (Figure 6).

Discussion
In the Southern Gulf of Mexico, intensive oil extraction and other hydrocarbons have been carried out since the 70s; thus, the ability to assess this industry's naturally polluting and associated anthropic impacts is important (Hernández-Arana et al., 2003, 2005).Therefore, understanding the distributional patterns of the polychaetes species present can help us understand the ecological processes occurring in the benthic ecosystems and how these organisms could be used as indicators of different ecological conditions at a larger scale (Gambi & Giangrande, 1986;Quiroz-Martínez et al., 2021).
In the Southern Gulf of Mexico, polychaete taxonomic distinctness changed across the longitudinal gradient, increasing from west to east, linked to a decline in species richness towards the terrigenous shelf in the west.The low taxonomic distinctness observed in terrigenous environments could be related to its more homogenous conditions, allowing the coexistence of taxonomically close species.In contrast, in the carbonated environment, the habitats are more varied, with highly heterogenous bottoms having different grain sizes, enabling the occurrence of taxonomically distinct species with different ecological preferences.Sediment composition is one of the main factors influencing the settling and distribution of benthic species.Sediment type and grain size variations determine organic matter content, dissolved gases, and nutrient concentrations (Granados Barba, 2001;Gray, 1974Gray, , 1981;;Magni et al., 2004Magni et al., , 2008Magni et al., , 2009)).The mixed sandy-mud bottoms dominate the Campeche Bank, favoring higher densities and species richness.This contrasts with the terrigenous sediments from the western Gulf, which negatively affect the benthic communities' diversity.Our analyses allowed us to distinguish the occurrence of three distinct polychaete species assemblages with the sedimentary distribution in the southern Gulf: one located in the western region on the terrigenous shelf, another associated with the carbonate shelf to the east, and an intermediate zone, in front of Laguna de Términos.
In the Bay of Campeche, the circulation is dominated by the mesoscale cyclonic eddy, which develops in the eastern bay and moves westward, enhancing small-scale anticyclonic and cyclonic eddies (Salas- de-León et al., 2008).In the study area, a coastal current moves, southward or south-westward, along the east coast mainly; approximately at 93°50W longitude, it converges with an eastward current along the coast to form a northward current along the shelf break (Salas- de-León et al., 2008).The sediment distribution in the southern Gulf of Mexico (Bay of Campeche) is mainly driven by wind-induced currents and mesoscale cyclonic eddies (Salas- de-León et al., 2008).Coarse sediments are located near the coast in the convergence region and the core of the anticyclonic eddy.In turn, fine sediments are situated off the shelf coast and in areas where the weakest currents are recorded (Salas- de-León et al., 2008).Carbonates dominate sediments on the eastern side of the Bay of Campeche since no river discharges exist or any other systems that could contribute terrigenous sediments, whereas, in the southernmost part of the bay, where major rivers are located, the bottom sediments are mainly of terrigenous origin (Salas- de-León et al., 2008).
The composition and structure of benthic assemblages can be influenced by numerous physical and biological factors (Hernández-Alcántara et al., 2011;Warwick, 1984) that define their zonation along environmental gradients (Bergen et al., 2001).In the study area, sediment variations were essential in structuring the polychaete assemblages.In addition, bathymetric fluctuations usually promote changes in other depth-dependent physical variables (Carvalho et al., 2012;Gray, 1981;Probert et al., 2001;Van Hoey et al., 2004).Particularly in the Southern Gulf, the distance between the sampling stations and their geographical proximity directly affected the faunal differences since the dissimilarity in polychaete species composition was a function of the distance separating them.Distance decay curves at regional or local scales are primarily shaped by the organisms' biological characteristics (niche requirements, dispersal ability) and environmental interactions.Nearby locations should typically have highly different species compositions when a species has a narrow geographic range, and power-law or negative exponential functions can be used to approximate distance-decay curves.However, close-by biological communities may be very similar when species have large spatial ranges relative to the study extent.This leads to a low or zero rate of change in community composition at the closest distances, while the actual decrease in community similarity would start at intermediate distances (Martín-Devasa et al., 2022).Several studies have shown that sediment composition can be crucial in controlling benthic communities in marine environments (Snelgrove & Butman, 1994).Notably, in the Gulf of Mexico, it has been demonstrated that sediment type and depth were leading factors in determining polychaete species' distribution and abundance patterns.Granados- Barba (2001) and Quiroz-Martínez et al. (2021, 2022) also showed that, in the southern Gulf, the polychaete species richness increased with depth.
Several studies have shown that depth is the major factor in structuring macrofauna (Hyland et al., 1991;Oug et al., 2017).
In contrast, in other studies, it has been found that sediment composition is the main factor affecting the benthic fauna (Hernández-Alcántara et al., 2011;Quiroz-Martínez et al., 2021;Snelgrove & Butman, 1994).In this sense, our results, founded on distance-based Redundancy analysis, showed that sediment type (and the dbMEMs) were the primary factor, rather than depth, to determine the polychaete faunal structure in the Southern Gulf of Mexico.For example, the syllids, usually small and highly motile, constituted the most diverse polychaete family and were mainly found in carbonated bottoms, where coarse sands mixed with coral and shell rubble dominate.However, the syllids become uncommon in mud and sandy-mud terrigenous sediments.
In the soft bottoms of the eastern Campeche Bay, syllids species richness generally followed a decreasing gradient from carbonate to terrigenous sediments (east to west) (Granados Barba et al., 2003).In addition, dbMEMs spatial processes explained a higher proportion of the variance in polychaete community structure than the environmental variables.Although polychaete larvae can disperse long distances, adult dispersal is often restricted and occurs over local/small spatial scales, mainly in intertidal areas and during local disturbance events (Checon & Amaral, 2017).Our results suggest that species composition was spatially correlated at shorter distances and that dispersal limitation and species response to spatially structured environmental gradients might be involved in determining polychaete distribution patterns.
The observed changes in sediment type affecting the establishment and development of benthic organisms can be associated with the water-sediment interphase dynamics that modify the organisms' excavation processes and limit the number of species that can settle (Díaz-Castañeda & Harris, 2004).According to Santibañez-Aguascalientes et al. ( 2020), the carbonate bottoms to the east comprise medium and fine calcarenites and oxides that provide suitable habitats for benthic communities (Carranza-Edwards et al., 1993;Cuevas Madrid et al., 2018;Granados Barba et al., 2003).In turn, the terrigenous bottoms are influenced by fluvial discharges dragging mud and sand sediments and harbor a high diversity of polychaete species adapted to turbid bottoms and terrigenous material (Hernández-Alcántara et al., 2011;Santibañez-Aguascalientes et al., 2018).
Characterizing distinct groups of invertebrates can be a complex process since, along biological gradients, a gradual transition between different assemblages frequently occurs (Gray, 1981;Van Hoey et al., 2004).In the present study, the separation of species groups based on multivariate analysis and taxonomic distinctness allowed us to distinguish three large polychaete assemblages; this indicates a gradual change in the composition of the polychaete fauna along the longitudinal gradient.Furthermore, polychaete assemblages distributed in the transitional region were different from the other sedimentary environments; this could be related to the high environmental variability of the transitional area, characterized by the transition from terrigenous to carbonated sediments, modulated by the Grijalva-Usumacinta plume and the oil extraction activities that affect the benthic fauna by reducing species diversity and richness (Salas- de-León et al., 2008;Santibañez-Aguascalientes et al., 2021, 2023).Ecological research in the Southern Gulf of Mexico shows that "Nortes" and summer rains provoking high river run-off are the most important physical processes influencing large-scale benthic community structure across the carbonate and terrigenous shelves (Cornell & Karlson, 2000;Granados Barba et al., 2003;Hernández-Arana et al., 2003, 2005;Soto & Escobar-Briones, 1995;Yáñez-Arancibia & Sánchez-Gil, 1983).
Finally, connectivity is another process that could explain the results presented in this study.Connectivity is an essential process in any marine ecosystem because it influences, among other things, the colonization of habitats, the maintenance of environmental services, and the persistence and survival of marine benthic populations (Gerber et al., 2014).According to Treml et al. (2015), there are 4-stage processes that explain the connectivity between populations; the second one (their Figure 1) includes the dispersion, transport, and movement of the organisms, which is strongly determined by the advection potential, which in turn is related to the magnitude of the currents.This could explain the three distinct faunal assemblages observed in this study, considering that the Southern Gulf of Mexico is highly dynamic, with the confluence of mesoscale physical processes and intense coastal currents, as mentioned above, which could have influenced the polychaete patterns of distribution presented in this study.For example, the connectivity between coral reef systems from the Gulf of Mexico has already been demonstrated using biological data (Carricart-Ganivet, 2004;Johnston & Akins, 2016;Mumby & Hastings, 2007;Murphy et al., 1999;Ortiz-Lozano et al., 2013) to the traditional coastal alongshore path, which is constrained to the continental shelf, with the first linking the Campeche Reef System with the Veracruz System.The water circulation pattern is driven by the Caribbean current during the spring and summer, with south to southwestern direction; however, during the autumn-winter seasons, the flow reverts to the east to a northeastern direction (Boicourt et al., 1998;Salas-Monreal et al., 2018).Various studies on connectivity in the southern Gulf of Mexico indicate a low east-west flow from the Campeche Bank to the Veracruz reefs; connectivity is stronger going west to east (Sanvicente-Añorve et al., 2014).However, ocean currents in this region are complex, with eddies and seasonal shifts in direction within the inner shelf, showing fluxes from east to west during summer and west to east during winter (Salas-Monreal et al., 2018).The pattern of taxonomic distribution observed in this study suggests connectivity between the reef systems In this study, the environment was characterized by spatial processes (dbMEMs), temperature, depth, salinity, and sediment variability; however, other variables, such as concentration of organic matter or nutrients, influence the distribution of polychaetes.Therefore, further studies should explore whether these relationships may produce higher percentages of variance explained by environmental structuring.
1. Thank you for submitting the above manuscript to PLOS ONE.During our internal evaluation of the manuscript, we found that this manuscript is very closely related to your previous published article: https://www.mdpi.com/1424-2818/13/9/425

Response
Our work attempts to describe the species richness patterns for this group, using multivariate analysis such as Cluster, Distance-based Redundancy Analysis, and Variation partitioning.We also analyzed the distance decay in similarity between sampling sites to measure beta diversity.Also, we define the key parameters that control the diversity in the study area.Although this manuscript relates to our previously published article (https://www.mdpi.com/1424-2818/13/9/425), it is substantially different.For one, we used the average taxonomic distinctness (Δ + , a measure of the average degree to which species are related to each other) as a diversity measure in the present manuscript; this diversity index differs from the more classical Sorensen index because it includes taxonomic relationships between the species in an assemblage and it has been shown that indices based on the taxonomic relatedness of species may give additional information to other univariate measures of biodiversity.In addition, the clustering method used in this study was also different, as we used Ward's Clustering instead of the UPGMA clustering method.The new diversity index and clustering method combined allowed us to characterize the polychaete fauna from the Southern Gulf of Mexico in more detail, as shown by the three faunal assemblages described in our study instead of the two assemblages described in our previous study.Finally, the spatial variables were simulated by the analysis of the Distance-based Moran's Eigenvector maps (dbMEMs) to reveal the influence of spatial variables on the community structure of annelid polychaetes from the Gulf of Mexico.An analysis which was not performed on the previous article.All the analysis therefore provided different, complementary information to the one obtained on our previous manuscript and allowed us to refine the patterns of distribution of the polychaetes.
or funders play any role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript?• NO -Include this sentence at the end of your statement: The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscriptbelow to enter a competing interest statement for this submission.On behalf of all authors, disclose any competing interests that could be perceived to bias this work-acknowledging all financial support and any other relevant financial or nonfinancial competing interests.This statement is required for submission and will appear in the published article if the submission is accepted.Please make sure it is accurate and that any funding sources listed in your Funding Information later in the submission form are also declared in your Financial Disclosure statement.View published research articles from PLOS ONE for specific examples.The authors have declared that no competing interests exist.NO authors have competing interestsEnter: The authors have declared that no competing interests exist.Authors with competing interestsEnter competing interest details beginning with this statement: I have read the journal's policy and the authors of this manuscript have the following competing interests/A" if the submission does not require an ethics statement.

Format
or tissues) Provide the name of the Institutional Animal Care and Use Committee (IACUC) or other relevant ethics board that reviewed the study protocol, and indicate whether they approved this research or granted a formal waiver of ethical approval • Include an approval number if one was obtained • If the study involved non-human primates, add additional details about animal welfare and steps taken to ameliorate suffering • If anesthesia, euthanasia, or any kind of animal sacrifice is part of the study, include briefly which substances and/or methods were applied • Field Research Include the following details if this study involves the collection of plant, animal, or other materials from a natural setting: to make all data underlying the findings described fully available, without restriction, and from the time of publication.PLOS allows rare exceptions to address legal and ethical concerns.See the PLOS Data Policy and FAQ for detailed information.No -some restrictions will apply Powered by Editorial Manager® and ProduXion Manager® from Aries Systems Corporation A Data Availability Statement describing where the data can be found is required at submission.Your answers to this question constitute the Data Availability Statement and will be published in the article, if accepted.Important: Stating 'data available on request from the author' is not sufficient.If your data are only available upon request, select 'No' for the first question and explain your exceptional situation in the text box.Do the authors confirm that all data underlying the findings described in their manuscript are fully available without restriction?Describe where the data may be found in full sentences.If you are copying our sample text, replace any instances of XXX with the appropriate details.If the data are held or will be held in a public repository, include URLs, accession numbers or DOIs.If this information will only be available after acceptance, indicate this by ticking the box below.For example: All XXX files are available from the XXX database (accession number(s) XXX, XXX.).• If the data are all contained within the manuscript and/or Supporting Information files, enter the following: All relevant data are within the manuscript and its Supporting Information files.• If neither of these applies but you are able to provide details of access elsewhere, with or without limitations, please do so.For example: Data cannot be shared publicly because of [XXX].Data are available from the XXX Institutional Data Access / Ethics Committee (contact via XXX) for researchers who meet the criteria for access to confidential data.The data underlying the results presented in the study are available from (include the name of the third party • Data will be posted in Obis in the future but is available from request from the author Powered by Editorial Manager® and ProduXion Manager® from Aries Systems Corporation . As the Caribbean Current flows northward through the Yucatan Channel, it is renamed the Yucatan Current.The Caribbean Tropical Surface Water (CTSW: T>28 °C and S>36.4) and the Caribbean Subtropical Underwater (CSUW: 22.2<T<26 °C and 36.4<S<36.7)influence the Campeche Bank since the Yucatan current flows into the Gulf, not only contributing significantly to the Loop Current, but another branch flows westward on the Yucatan continental shelf and Campeche Bank (Aldeco-Ramírez et al., 2009).

Figure 1 .
Figure 1.Sampling stations in the Southern Gulf of Mexico.
longitude and latitude coordinates of each station and the environmental distance between each pair of stations with Euclidean distances based on all environmental variables measured.We used each sampling station's sediment composition, depth, temperature, and salinity as environmental variables.Distance-based Moran's eigenvector maps (dbMEMs), redundancy analysis (dbRDA), and variation partitioning were used to quantify the effect of spatial structure and biotic and abiotic variables on community dissimilarity (taxonomic distinctness).Before the dbRDA analysis, forward selection procedures were conducted to select the most significant spatial and environmental variables; only the dbMEMs corresponding to positive autocorrelation as spatial variables were used in the analyses.Distance-based Moran eigenvector maps (dbMEMs) for the matrix of the sites' geographical coordinates were used to describe the dataset's spatial structure.The initial dbMEMs model the large-scale spatial correlation, while the last dbMEMs correspond to fine-scale spatial correlations, which may capture variation at the sampling site scale.All nonsignificant (P > 0.05) variables were eliminated from further analyses.Afterward, we performed a variation partitioning analysis to assess the relative contribution of the significant spatial and environmental variables in explaining the community dissimilarity(Peres-Neto et al., 2006).The variation was partitioned into four fractions: (a) influence of environmental variables alone (sediment composition, depth, temperature, and salinity); (b) influence of spatial variables; (c) influence of spatial variables alone (distance-based Moran eigenvector maps; dbMEMs); and (d) unexplained variation.The variation partitioning analysis included only the variables selected by this procedure.The significance of each fraction was tested by Monte Carlo permutation tests.In addition, the relative contribution of each significant spatial and environmental variable in determining community dissimilarity was assessed by partial dbRDA to partial out the effect of each significant environmental/spatial variable.A Monte Carlo permutation test was run to test the significance of each significant spatial and environmental variable.All analyses were run in R-4.0.3 with the "vegan" and other built-in packages(Oksanen et al., 2020;R Development Core Team, 2020) Results

Figure 2 .
Figure 2. Simulated distribution of average taxonomic distinctness (Δ + ) and the variation in taxonomic distinctness Λ + (theoretical mean, horizontal, dashed line) for random subsets of species from the complete species list of 173 polychaete species from the Southern Gulf of Mexico and the 95% confidence limits (the funnel).

Figure 3 .
Figure 3. Taxonomic dissimilarity index θ + plotted against distance for all pairwise comparisons between sampling stations.The blue line represents the exponential function describing the increase in dissimilarity (a.intercept = 26.56;b.slope = 0.001).

"
Terrigenous assemblage" (Cluster 1 in red) was integrated by stations located in the southernmost region of the terrigenous shelf of the Gulf of Mexico.Interestingly, sites 114, 126, and 140, associated with Cluster 1, are geographically located on the Carbonate Shelf.The "Transitional" assemblage (Cluster 2 in green) grouped the stations from the central part of the study area in front of Laguna de Términos, considered a sedimentary transitional region.Finally, the "Carbonated" assemblage (Cluster 3 in blue) included stations located on the carbonate shelf on the eastern Gulf of Mexico.

Figure 4 .
Figure 4. Cluster analysis based on taxonomic dissimilarity and Ward's clustering method

Figure
Figure 5. Distance-based redundancy analysis (dbRDA) ordination plot showing the relationships between sampling stations based on polychaete taxonomic distinctness and spatial and environmental constraints.
. In addition, Salas-Monreal et al. (2018) elucidated the different pathways that connect the Western Gulf of Mexico reef systems using monthly geostrophic velocities, sea surface temperature, and chlorophyll-a values.They found high connectivity between the reef systems that bring substrates, suspended matter, and organisms from the Campeche Reef System (CRS) in the Campeche Bank to the Veracruz Reef System (VRS) in the Southwestern Gulf of Mexico.According to them, the geostrophic current velocities suggested various paths connecting the various coral reef habitats, with two open ocean pathways in addition