Contrasting Early Ordovician assembly patterns
highlight the complex initial stages of the Ordovician
Radiation

Farid Saleh, Pauline Guenser, Corentin Gibert, Diego Balseiro, Fernanda

Serra, Beatriz G. Waisfeld, Jonathan B. Antcliffe, Allison C. Daley, M.

Gabriela Mángano, Luis A. Buatois, et al.

► To cite this version:

Farid Saleh, Pauline Guenser, Corentin Gibert, Diego Balseiro, Fernanda Serra, et al. . Contrasting Early Ordovician assembly patterns highlight the complex initial stages of the Ordovician Radiation. Scientific Reports, 2022, 12, 10.1038/s41598-022-07822-z . insu-03710171v2

HAL Id: insu-03710171

https://insu.hal.science/insu-03710171v2

Submitted on 30 Jun 2022

**HAL** is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for updates

OPEN

# Contrasting Early Ordovician assembly patterns highlight the complex initial stages of the Ordovician Radiation

Farid Saleh<fnref n="1" />,²✉, Pauline Guenser³, Corentin Gibert⁴✉, Diego Balseiro⁵,⁶, Fernanda Serra⁵,⁶, Beatriz G. Waisfeld⁵,⁶, Jonathan B. Antcliffe⁷, Allison C. Daley⁷, M. Gabriela Mángano⁸, Luis A. Buatois⁸, Xiaoya Ma<fnref n="1" />,²,⁹✉, Daniel Vizcaïno¹⁰ & Bertrand Lefebvre¹¹

The Early Ordovician is a key interval for our understanding of the evolution of life on Earth as it lays at the transition between the Cambrian Explosion and the Ordovician Radiation and because the fossil record of the late Cambrian is scarce. In this study, assembly processes of Early Ordovician trilobite and echinoderm communities from the Central Anti-Atlas (Morocco), the Montagne Noire (France), and the Cordillera Oriental (Argentina) are explored. The results show that dispersal increased diachronically in trilobite communities during the Early Ordovician. Dispersal did not increase for echinoderms. Dispersal was most probably proximally triggered by the planktic revolution, the fall in seawater temperatures, changes in oceanic circulation, with an overall control by tectonic frameworks and phylogenetic constraints. The diachronous increase in dispersal within trilobite communities in the Early Ordovician highlights the complexity of ecosystem structuring during the early stages of the Ordovician Radiation. As Early Ordovician regional dispersal was followed by well-documented continental dispersal in the Middle/Late Ordovician, it is possible to consider that alongside a global increase in taxonomic richness, the Ordovician Radiation is also characterized by a gradual increase in dispersal.

The fossil record of the Cambrian Period documents the earliest evolution and diversification of eumetazoans on Earth<fnref n="1" />, during which animals developed a wide range of morphological features and characteristics²,³. The ancient ecology of these animals can be reconstructed based on the exceptional fidelity of anatomy preserved in Burgess Shale (BS)-type fossil sites⁴, such as Burgess Shale, Chengjiang, and Qingjiang biotas⁴⁻⁶. The Cambrian Explosion was followed by the Ordovician Radiation². During the Ordovician, diversity significantly increased and was accompanied by structural changes in ecosystems⁷⁻⁹. However, most palaeoecological studies have focused on the Middle Ordovician and onwards. As a result, the Early Ordovician remains an understudied interval, despite having the potential to help elucidate the transition between the Cambrian Explosion and the Ordovician Radiation, particularly given that the fossil record of the late Cambrian is scarce (i.e., the Furongian gap)¹⁰.

This work investigates the ecological processes that structure community composition during the Early Ordovician of three regions: The Central Anti-Atlas (CAA; Morocco), Montagne Noire (MN; France), and Cordillera

<footnote n="1">¹Yunnan Key Laboratory for Palaeobiology, Institute of Palaeontology, Yunnan University, Kunming, China. ²MEC International Joint Laboratory for Palaeobiology and Palaeoenvironment, Institute of Palaeontology, Yunnan University, Kunming, China. ³Université Claude Bernard Lyon 1, CNRS, UMR5023, LEHNA, Université de Lyon, 69622 Villeurbanne, France. ⁴Laboratoire de la Préhistoire à l'Actuel: Culture, Environnement et Anthropologie (PACEA, UMR 5199 CNRS, INEE), University of Bordeaux, Bordeaux, France. ⁵Facultad de Ciencias Exactas, Físicas y Naturales, Universidad Nacional de Córdoba, Córdoba, Argentina. ⁶Centro de Investigaciones en Ciencias de la Tierra (CICTERRA), Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Av. Vélez Sarsfield 1611, CP X5016GCA Córdoba, Argentina. ⁷Institute of Earth Sciences, University of Lausanne, Géopolis, 1015 Lausanne, Switzerland. ⁸Department of Geological Sciences, University of Saskatchewan, Saskatoon, SK S7N 5E2, Canada. ⁹Centre for Ecology and Conservation, University of Exeter, Penryn, UK. ¹⁰Carcassonne, France. ¹¹Université Claude Bernard Lyon 1, École Normale Supérieure de Lyon, CNRS, UMR5276, LGL-TPE, Université de Lyon, Villeurbanne, France. ¹²email: farid.nassim.saleh@gmail.com; corentingibert@gmail.com; X.Ma2@exeter.ac.uk</footnote>
<figure><img src="image_3.png" /><figcaption>**Figure 1.** (a) Stratigraphic correlations between the Central Anti-Atlas (CAA), the Montagne Noire (MN), and the Cordillera Oriental (CO). The studied intervals within this work are highlighted in orange. Note that for the Fezouata Shale orange intervals were not selectively chosen, as those are the only ones yielding echinoderm and trilobite fossils. Se. = Series, St. = Stage, TS = Time Slices²⁷; Morocco²⁸; Cordillera Oriental¹⁹; Montagne Noire¹⁴. (b) Paleogeography of the three regions; modified from Saleh et al.²³.</figcaption></figure>

Oriental (CO; Argentina) (Fig. 1a,b). Because not all regions preserve all types of taxa owing to taphonomic biases¹¹,¹², this study is focusing on echinoderms and trilobites as these are readily preserved and are well studied in the three regions (Fig. 2)¹³⁻²³. Moreover, trilobites are the most abundant and diverse group in the three studied regions followed by echinoderms (in the CAA and the MN)¹³⁻²³. Communities can be niche-assembled (Fig. 3a) when they are at equilibrium, closed, and balanced, with a stable taxonomic composition maintained by deterministic processes like positive (e.g., mutualism) and negative (e.g., competition) interactions between individuals²⁴⁻²⁶. Communities are dispersal-assembled (Fig. 3b) when they are open, in non-equilibrium state with their compositions depending on continuous dispersion between communities, and subjected to random originations, and extinctions²⁴⁻²⁶. Once seen as incompatible hypotheses, niche- and dispersal-perspectives are now seen as opposite ends of a continuum²⁴⁻²⁶. Within this continuum, the taxonomic composition of an assemblage may result to a greater or lesser extent of each process, or from an equal combination of the two²⁴⁻²⁶.

Identifying niche and dispersal processes is a fundamental question in modern ecology and has become a central theme for conservation studies²⁴. Novel approaches have been developed that can identify assembly-processes within fossil datasets, and describe changes in processes over evolutionary timescales, and their potential links with abiotic and biotic events²⁴. In order to quantify niche and dispersal processes within trilobite and echinoderm assemblages in the three regions, this work applies the novel Dispersal and Niche Continuum index (DNCI) analyses to understand changes in community assembly between the Tremadocian and the Floian²⁴, i.e., the two Early Ordovician stratigraphic stages. The DNCI quantifies the relative contribution of niche and dispersal processes to the taxonomic dissimilarity between assemblages, which are defined as statistical clusters of fossiliferous levels in each region²⁴. The DNCI ranges from significantly positive (i.e., niche-assembled community) to significantly negative (i.e., dispersal-assembled community) values. The lower the value of the DNCI, the more important dispersal was in a certain assemblage²⁴. A community was equally niche- and dispersal-assembled when the mean DNCI and its confidence interval overlap with zero²⁴. Once niche and dispersal contributions were quantified, they are then placed in the larger evolutionary framework of trilobites and echinoderms in addition to the geological context of the three localities.

From a terminological perspective, the generic term “assemblage” is used to refer to all taxa inhabiting a particular area, regardless of the spatial scale considered (e.g., community, meta-community, bioregions). “Assembly-processes” (i.e., niche and dispersal) are the ecological forces bringing together and maintaining these taxa within assemblages. Niche- and dispersal-assembly should be differentiated as terms from cosmopolitanism and endemism. For instance, within a community model, if dispersal is the only ecological force considered, this can lead to an endemic and/or a cosmopolitan signal.

## Results and discussion

Trilobites in the CAA are assembled equally by both niche- and dispersal-processes during the Early Ordovician (i.e., the range of the DNCI and its confidence interval overlap with zero in the Tremadocian and the Floian; Fig. 4b). Trilobites in the MN are dispersal-assembled in the Tremadocian and the Floian (i.e., DNCI is negative and does not overlap with zero; Fig. 4b). Trilobite assemblages in the CO are equally driven by both niche and dispersal processes in the Tremadocian and mainly driven by dispersal processes in the Floian (Fig. 4b). Trilobites in the MN and CO show little overlap between the Tremadocian and the Floian (Fig. 4b). This indicates a clear difference between the Tremadocian and the Floian records of trilobites in the CO and the MN. This difference is characterized by a transition towards more negative DNCI values (Fig. 4b) associated with a decrease in niche processes, and an increase in dispersal²⁴. The contrasting results for trilobites between the three regions (Fig. 4b) do not reflect a differential sampling effort bias because trilobites in the CAA, the CO, and the MN are equally sampled (i.e., S<sub>obs</sub>/S<sub>ex</sub> ~ 0.9; Fig. 4a).

The increase in trilobite regional dispersal during the Early Ordovician (at least in the MN and the CO) correlates with environmental changes during this time interval. During the Furongian, the planktic revolution started, and increased the availability of a globally distributed food source for animals²⁹. Furthermore,
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Early Ordovician trilobites from the Central Anti-Atlas (a–d), the Montagne Noire (e,f), and the Cordillera Oritental (g,h) in addition to Early Ordovician echinoderms from the Central Anti-Atlas (i–l), and the Montagne Noire (m–p). (a) *Euloma* sp., late Tremadocian; AA.TGR1c.OI.232, (b) *Platypeltoides magrebiensis*, late Tremadocian; AA.TGR0a.OI.132, (c) *Bavarilla zemmourensis*, late Tremadocian; AA.BIZ15.OI.16, (d) *Anacheirurus adserai*, late Tremadocian; AA.BIZ15.OI.291, (e) *Euloma filacovi*, late Tremadocian; UCBL-FSL 740109, (f) *Paramegalaspis immarginata*, late Tremadocian; UCBL-FSL 740027, (g) *Pliomeridius sulcatus*, middle Floian; CEGH-UNC-16311, (h) *Leptoplastides granulosa*; middle Tremadocian, CEGH-UNC-21324, (i) *Villebrunaster* sp., late Tremadocian; FSL 424 961, (j) *Macrocystella bohemica*, late Tremadocian, MHNM.15690.196, (k) Eocrinoidea indet., late Tremadocian; MHNM.15690.141, (l) Cornuta n. gen. n. sp. cf. "*Phyllocystis* jingxiensis; FSL.712971, (m), *Villebrunaster thorali*, late Tremadocian; UCBL-FSL 168673, (n) *Macrocystella azaisi*, late Tremadocian; MBB-GG50, (o) *Aethocrinus moorei*, late Tremadocian; MNHN.A49684, (p) *Phyllocystis blayaci*, late Tremadocian; UCBL-FSL 168703. 1 cm scale bars (a–f; i–p), magnification = 0.8× (g) and 1.8× (h). Abbreviations (repositories): AA Cadi-Ayyad University, Marrakesh (Morocco), MBB Musée du Biterrois, Béziers (France), MNHN Muséum National d'Histoire Naturelle, Paris (France), UCBL-FSL palaeontological collections, Lyon 1 University, Villeurbanne (France), CEGH-UNC Universidad Nacional de Córdoba, Córdoba (Argentina).</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Niche- versus Dispersal-assembly perspectives. (a) In a niche-assembly scenario, communities are closed and balanced, built by a limited number of species at equilibrium (similar assemblages at t + 1 and t + 2). Niche taxonomic composition is controlled by deterministic (i.e., non-random) biotic interactions between taxa and their environment. (b) Dispersal-assembly scenario in which communities are constantly changing (i.e., non-equilibrium; different assemblages between t + 1 and t + 2) owing to random processes such as migrations.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Figure 4.** (a) Estimation of the sampling effort [observed sampling ($S_{obs}$) over extrapolated sampling ($S_{ex}$)] of echinoderms (echino) and trilobites (trilo) in the Central Anti-Atlas (CAA), the Montagne Noire (MN), and the Cordillera Oriental (CO). $S_{obs}$ and $S_{ex}$ were also calculated while separating between the Floian and the Tremadocian of the three regions (Supplementary Material 1). (b) DNCI analyses for trilobites and echinoderms during the Tremadocian (Trema) and Floian of the Central Anti-Atlas, the Montagne Noire, and the Cordillera Oriental.</figcaption></figure>

temperatures achieved a global hothouse at the onset of the Furongian (also known as the Steptoean Positive Carbon Isotope Excursion or SPICE event) and started to decrease right afterward in the early/middle Furongian up until the Late Ordovician³⁰. The global cooling trend established a more prominent global oceanic circulation in the Middle/Late Ordovician than in the late Cambrian³¹,³². This means that at least transitional conditions of oceanic ventilation and circulation existed in the Tremadocian and the Floian (periods of high sea level)³³. Thus, it is most likely that oceanic circulation, cooling, and the planktic revolution facilitated dispersal because these environmental changes are linked to animals' specific evolutionary patterns and ecological behavior²⁹. For instance, an increased colonization of deep-sea floors is indicated by ichnologic evidence, and took place during the Early Ordovician³⁴. This increased animal activity in deep seas has been attributed to augmented oceanic ventilation³⁴. Moreover, in the Furongian, for the first time, trilobites started to disperse into new habitats by exploring brackish water environments, a trend that continued well into the Middle Ordovician²⁰,³⁵. Additionally, most of the known early to middle Cambrian trilobite larvae were not spherical but flat in shape (adult-like protaspides) and are interpreted as having lived on or near the seafloor³⁶⁻³⁹. In the Furongian, some trilobite groups developed bulbous configurations³⁶. Bulbous larvae (non-adult-like protaspides) were adapted to a planktic mode of life⁴⁰. Although bulbous larvae have so far been found in few Ordovician lineages³⁶,⁴⁰, Bignon et al.⁴¹ argued that this type of larvae must have been more abundant. Thus, some of the lineages yielding either bulbous larvae or juveniles of this type might explain the observed increase in trilobite dispersal. However, this evolutionary trait only explains the increase in dispersal but it does not explain why trilobite dispersal was non-uniform between regions. The non-uniform dispersal pattern might be resulting from the limited temporal sampling in the CAA (Fig. 1a), which represents a significantly more restricted window in comparison to the CO, and the MN (Fig. 1a). It can also result from differences in regional settings (e.g., geodynamic framework). During the Early Ordovician, the CAA was a restricted area⁴²,⁴³. As such, the CAA was isolated from the global oceanic circulation (i.e., the newly established global currents did not increase inter-communication between CAA sites), which may explain why trilobites did not disperse within this area as early as trilobites in other regions (Fig. 4b). Moreover, the CO was partially protected by a volcanic arc to the west which might explain why the composition of the
assemblages of this area remained equally niche- and dispersal-assembled in the Tremadocian in comparison to the dispersal-assembled MN⁴⁴.

Another aspect of these data is that echinoderms in the CAA and the MN are assembled equally by both niche- and dispersal-processes during the Early Ordovician, without showing a change in community assembly processes (Fig. 4b). Due to sampling biases, trilobite and echinoderm data cannot be directly compared [i.e., S<sub>obs</sub>/S<sub>ex</sub>(echinoderm) ~ 0.6 < S<sub>obs</sub>/S<sub>ex</sub>(trilobite) ~ 0.9; Fig. 4a]. However, the similarity between echinoderms in the CAA and the MN is not biased by sampling (i.e., S<sub>obs</sub>/S<sub>ex</sub> ~ 0.6 for both regions; Fig. 4a). This essentially means that echinoderm assembly patterns were uniform between regions and were not differentially affected by the environmental changes happening at the global scale. Unfortunately, it is impossible to further discuss this finding owing to the absence of fossilized echinoderm larvae. It is assumed that the possession of planktotrophic larvae is a plesiomorphic condition in the five extant classes⁴⁵, which appeared and diversified either in the Early Ordovician (asteroids, crinoids, ophiuroids) or in the Darriwilian (echinoids, holothuroids)⁴⁶. However, confirming that echinoderm planktotrophic larvae were acquired by the Early Ordovician, is not possible⁴⁷

The relevance of dispersal on assemblages is also dependent on the possibility for new immigrants to enter local communities (i.e., openness of communities). Therefore, an increase in dispersal in community assembly, also suggests that local communities became more susceptible to invasion by immigrants. Such a change can be driven by an increase in the amount of unused resources⁴⁸ related to the diversification of phytoplankton during the plankton revolution²⁹. The increase in dispersal could also result from transient taxa⁴⁹ as new dispersive morphologies (e.g., bulbous larvae) enable the development of new demographics (i.e., source/sink)⁵⁰ and evolutionary (increased taxa longevity) dynamics⁵¹. For instance, non-specialized taxa could survive in remote and hostile environments if they possessed an adequate dispersive morphology that allowed for a continuous influx of new individuals from their source community (i.e., competition-colonization trade-off).

With these new data, it is now possible to place the transition between the Cambrian Explosion and the Ordovician Radiation in an ecological context. The Cambrian Explosion is primarily linked with increased niche packing and partitioning⁵². Previous studies made at a global scale suggested that niche-processes remained high during the first two Ordovician stages⁵², indicating that the remarkable decrease in niche-assembly was only acquired later, starting in the Middle Ordovician⁵³. The overlap of the DNCI with zero in the Tremadocian, and to some extent in the Floian, indicates that both niche and dispersive regimes are equally detected in the analyses, highlighting that the Early Ordovician was a time of ecological transition between the Cambrian Explosion and the Ordovician Radiation. This study provides a new perspective, complementing previously described dynamics by showing that changes in community assembly (i.e., increase in dispersal in comparison to the Cambrian) can already be observed in the Early Ordovician, especially at a regional scale.

## Conclusion

The Early Ordovician is a complex and poorly understood time in the evolution of life on Earth. Dispersal increased during this time and this was likely triggered by evolutionary dynamics (e.g., the planktonic revolution) with underlying environmental forcing (e.g., oceanic cooling and circulation changes). However, dispersal increase was not homogenous for all animal groups and regions (Fig. 4b) as it was controlled by a combination of clade inherent features and the tectonic framework of each region. The importance of the results presented herein is that they (1) were obtained using a novel method that is independent of classical taxonomic richness metrics, and (2) highlight the complexity of the early stages of the Ordovician Radiation. Considering that diachronous regional dispersal in the Early Ordovician was followed by episodes of continental dispersal during the Middle Ordovician⁵⁴⁻⁵⁷, it is possible to suggest that the Ordovician Radiation can be envisaged as facilitated by a set of diachronous (i.e., in one region prior to another) increases (i.e., from regional to continental) in dispersal. This increase in dispersal could also represent a transition from fragmented habitat patches separated by non-livable oxygen-depleted zones in the Early Ordovician to a more homogenous ocean that is dominantly cooler and oxygen-rich in the Middle Ordovician⁵⁸. In the future, it is essential to apply DNCI analyses to a wider range of regions in order to investigate whether either latitudinal or longitudinal dispersal gradients exist during this time interval. These future studies will bring insights on the spreading of regional dispersal during the Early Ordovician, and on how to calibrate the diachronicity of dispersal during this critical time interval.

## Materials and methods

**Studied fossils and regions.** The echinoderm and trilobite data in this study are primarily field-based and result from intensive sampling in the three investigated areas over several decades. Data from the Montagne Noire are based on the public collections of Lyon 1 University, France. The Central Anti-Atlas data was extracted from the public collections of Cadi Ayyad University, Morocco, while data from the Cordillera Oriental is from both the literature and from the public collections at the National University of Córdoba, Argentina. This work comprises more than 300 well-dated levels (Fig. 1a). “Levels” represent spatial-temporal limited intervals. For instance, two stratigraphically distinct fossiliferous horizons occurring at the same site were considered as two different levels. Moreover, two coeval fossiliferous layers were considered as distinct “levels”, when they occur in two different sites of a same region. Additional detailed information on the studied fossils is provided in Supplementary Material 1 (SM1).

In order to avoid local depositional biases on the analyses, occurrence data within the same range of depositional conditions between the fair-weather wave base and the storm-weather wave base (i.e., offshore) were included in this study. Moreover, there is no dominant oxygen or salinity stress during deposition of most investigated levels and the three regions have a very comparable substrate¹⁹,²²,⁵⁹. Therefore, levels in which an oxygen stress is evident in the CAA²² were excluded from the analyses. During the Early Ordovician, the CAA was a relatively restricted sea⁴²,⁴³. The CO was partially protected by a volcanic arc to the west⁴⁴. The MN was
an open sea in direct contact with the ocean⁴²,⁴³. More information on the geological context is provided in the
expanded materials and methods section in SM1.

**Statistical and ecological approaches.** Incidence matrices reflecting occurrences (absence/presence) in all layers for each region were built (i.e., five datasets; Supplementary Material 2). The sampling effort observed as a ratio between observed and extrapolated sampling (S<sub>obs</sub>/S<sub>ex</sub>) was calculated using the iNEXT package on R software⁶⁰ and its confidence interval was obtained using the function "Ratio of counts CI" of the software PAST⁶¹ and is reported on Fig. 4a. Each dataset (e.g., echinoderms in the CAA) was then split into two subsets (e.g., Tremadocian and Floian echinoderms in the CAA; Supplementary Material 3). Cluster analyses (UPGMA) were performed on each subset, and results are present in SM3. The significance of these clusters was then tested and confirmed as statistically significant with an analysis of similarity (i.e., ANOSIM; SM3). These analyses were performed in PAST using the Jaccard index, which is the most commonly used index in paleoecological investigations⁶². DNCI identifies the main assembly process driving species composition between clusters by comparing empirical and null SIMPER profiles⁶³. The SIMPER method has been designed to compute and depict the contribution of taxa to the overall taxonomic dissimilarity (OAD) between two clusters/groups of assemblages/communities⁶³. The contribution of taxa to OAD depends on the relationship of β-diversity within clusters of assemblages and β-diversity between clusters of assemblages (Supplementary Material 1). In other words, SIMPER considers: (1) taxa that are largely distributed within one cluster, (2) taxa that are largely and evenly distributed between two clusters; (3), rare taxa that are evenly distributed and (4) rare taxa that are only distributed within one cluster. Then an empirical SIMPER profile is produced by ordering the contribution of each taxa to the OAD. This profile is then compared to three null permuted SIMPER (PERSIMPER) profiles (only niche, only dispersal or the equal contribution of niche and dispersal processes). The closer the null profile, the more important is the respective assembly model (Supplementary Material 1). In this sense, the taxonomic dissimilarity and the β-diversity are at the core of PER-SIMPER/DNCI. However, there is no correlation between dissimilarity indices or β-diversity variations and DNCI variations. To perform DNCI analyses, sheets in SM3 (i.e., 10 in total) were exported as CSV files and are present in a compressed folder (Supplementary Material 4). Three levels that do not belong to any cluster (i.e., highlighted in yellow in SM3) were removed from SM4 because a cluster has to have by definition more than one level to be used in the following DNCI analyses²⁴. Clusters also need to be composed of a similar number of levels in order to avoid bias on niche and dispersal processes identification⁶³. If the compared clusters are uneven, DNCI will shift toward zero, leading to a false identification of an equal contribution of niche and dispersal processes²⁴. A subsampling procedure can be used to make the clusters of equal size and correct this bias. In this case, when clusters differ in size (e.g., cluster 1 has 40% more levels than cluster 2), a subsampling randomly reducing the larger cluster(s) to the size of the smallest cluster was implemented 100 times to minimize the sampling bias effect and produce robust mean DNCI values. DNCI analyses were done using DNCImper⁶⁴ R package available on Github (https://github.com/Coren tin-Gibert-Paleontology/DNCImper). A customized R Markdown script of DNCI analyses is provided in Supplementary Material 5. More details regarding the methods used, the choices made, in addition to explanations on how the DNCI operates are provided in SM1.

**Data availability**

All data needed to evaluate the conclusions in the paper are present in the Supplementary Materials. Additional
data related to this paper may be requested from the authors.

Received: 11 November 2021; Accepted: 25 February 2022

Published online: 09 March 2022

**References**

1. Marshall, C. R. Explaining the Cambrian "explosion" of animals. *Annu. Rev. Earth Planet. Sci.* **34**, 355–384 (2016).

2. Bush, A. M. & Bambach, R. K. Paleoecologic megatrends in marine metazoa. *Annu. Rev. Earth Planet. Sci.* **39**, 241–269 (2011).

3. Smith, M. P. & Harper, D. A. Causes of the Cambrian explosion. *Science* **341**(6152), 1355–1356 (2013).

4. Daley, A. C., Antcliffe, J. B., Drage, H. B. & Pates, S. Early fossil record of Euarthropoda and the Cambrian explosion. *Proc. Natl. Acad. Sci.* **115**(21), 5323–5331 (2018).

5. Fu, D. et al. The Qingjiang biota—A Burgess Shale–type fossil Lagerstätte from the early Cambrian of South China. *Science* **363**(6433), 1338–1342 (2019).

6. Nanglu, K., Caron, J. B. & Gaines, R. R. The Burgess Shale paleocommunity with new insights from Marble Canyon, British Columbia. *Paleobiology* **46**(1), 58–81 (2020).

7. Sepkoski, J. J. Jr. The Ordovician radiations: Diversification and extinction shown by global genus-level taxonomic data. In *Ordovician Odyssey: Short Papers, 7th International Symposium on the Ordovician System* (eds Cooper, J. D. *et al*. ) 393–396 (Pacific Section Society for Sedimentary Geology (SEPM), 1995).

8. Servais, T., Cascales-Miñana, B. & Harper, D. A. The Great Ordovician Biodiversification event (GOBE) is not a single event. *Paleontol. Res.* **25**(4), 315–328 (2021).

9. Harper, D. A., Cascales-Miñana, B., Kroeck, D. M. & Servais, T. The palaeogeographical impact on the biodiversity of marine faunas during the Ordovician radiations. *Glob. Planet. Change* **207**, 103665 (2021).

10. Harper, D. A. et al. The Furongian (late Cambrian) biodiversity gap: Real or apparent? *Palaeoworld* **28**(1–2), 4–12 (2019).

11. Saleh, F. et al. Taphonomic bias in exceptionally preserved biotas. *Earth Planet. Sci. Lett.* **529**, 115873 (2020).

12. Saleh, F. et al. A novel tool to untangle the ecology and fossil preservation knot in exceptionally preserved biotas. *Earth Planet. Sci. Lett.* **569**, 117061 (2021).

13. Vizcaíno, D. & Lefebvre, B. Les échinodermes du Paléozoïqueinférieur de Montagne Noire: Biostratigraphie et paléodiversité. *Geobios* **32**(2), 353–364 (1995).

14. Vizcaíno, D. & Álvaro, J. J. Adequacy of the Early Ordovician trilobite record in the southern Montagne Noire (France): Biases for biodiversity documentation. *Earth Environ. Sci. Trans. R. Soc. Edinb.* **93**(4), 393–401 (2002).
15. Lefebvre, B. et al. Palaeoecological aspects of the diversification of echinoderms in the Lower Ordovician of central Anti-Atlas, Morocco. *Palaeogeogr. Palaeoclimatol. Palaeoecol.* **460**, 97–121 (2016).

16. Lefebvre, B. et al. Exceptionally preserved soft parts in fossils from the Lower Ordovician of Morocco clarify stylophoran affinities within basal deuterostomes. *Geobios* **52**, 27–36 (2019).

17. Martin, E. L. O. et al. Biostratigraphic and palaeoenvironmental controls on the trilobite associations from the Lower Ordovician Fezouata Shale of the central Anti-Atlas, Morocco. *Palaeogeogr. Palaeoclimatol. Palaeoecol.* **460**, 142–154 (2016).

18. Waisfeld, B. G. & Balseiro, D. Decoupling of local and regional dominance in trilobite assemblages from northwestern Argentina: New insights into Cambro-Ordovician ecological changes. *Lethaia* **49**(3), 379–392 (2016).

19. Serra, F., Balseiro, D. & Waisfeld, B. G. Diversity patterns in upper Cambrian to Lower Ordovician trilobite communities of north-western Argentina. *Palaeontology* **62**(4), 677–695 (2019).

20. Serra, F., Balseiro, D., Vaucher, R. & Waisfeld, B. G. Structure of trilobite communities along a delta-marine gradient (lower Ordovician; Northwestern Argentina). *Palaios* **36**(2), 39–52 (2021).

21. Saleh, F., Lefebvre, B., Hunter, A. W. & Nohejlová, M. Fossil weathering and preparation mimic soft tissues in eocrinoid and somasteroid echinoderms from the Lower Ordovician of Morocco. *Microsc. Today* **28**(1), 24–28 (2020).

22. Saleh, F. et al. Insights into soft-part preservation from the Early Ordovician Fezouata Biota. *Earth Sci. Rev.* **213**, 103464 (2021).

23. Saleh, F. et al. Large trilobites in a stress-free Early Ordovician environment. *Geol. Mag.* **158**(2), 261–270 (2021).

24. Vilmi, A. et al. Dispersal–niche continuum index: A new quantitative metric for assessing the relative importance of dispersal versus niche processes in community assembly. *Ecography* **44**(3), 370–379 (2021).

25. Hubbell, S. P. *A Unified Theory of Biodiversity and Biogeography* (Princeton University Press, 2001).

26. Gravel, D., Canham, C. D., Beaudet, M. & Messier, C. Reconciling niche and neutrality: The continuum hypothesis. *Ecol. Lett.* **9**(4), 399–409 (2006).

27. Bergström, S. M., Chen, X., Gutiérrez-Marco, J. C. & Dronov, A. The new chronostratigraphic classification of the Ordovician system and its relations to major regional series and stages and to δ13C chemostratigraphy. *Lethaia* **42**, 97–107 (2009).

28. Lefebvre, B. et al. Age calibration of the Lower Ordovician Fezouata Lagerstätte, Morocco. *Lethaia* **51**(2), 296–311 (2018).

29. Servais, T. et al. The onset of the ‘Ordovician Plankton Revolution’in the late Cambrian. *Palaeogeogr. Palaeoclimatol. Palaeoecol.* **458**, 12–28 (2016).

30. Lee, J. H. & Riding, R. Marine oxygenation, lithistid sponges, and the early history of Paleozoic skeletal reefs. *Earth Sci. Rev.* **181**, 98–121 (2018).

31. Servais, T., Danelian, T., Harper, D. A. T. & Munnecke, A. Possible oceanic circulation patterns, surface water currents and upwelling zones in the Early Palaeozoic. *GFF* **136**(1), 229–233 (2014).

32. Rasmussen, C. M. et al. Onset of main Phanerozoic marine radiation sparked by emerging Mid Ordovician icehouse. *Sci. Rep.* **6**(1), 1–9 (2016).

33. Edwards, C. T. Links between early Paleozoic oxygenation and the Great Ordovician Biodiversification Event (GOBE): A review. *Palaeoworld* **28**(1–2), 37–50 (2019).

34. Buatois, L. A. et al. Quantifying ecospace utilization and ecosystem engineering during the early Phanerozoic—The role of bio-turbation and bioerosion. *Sci. Adv.* **6**(33), eabb0618 (2020).

35. Mängano, M. G. et al. Were all trilobites fully marine? Trilobite expansion into brackish water during the early Palaeozoic. *Proc. R. Soc. B* **288**(1944), 20202263 (2021).

36. Park, T. Y. S. et al. Ontogeny of the Furongian (late Cambrian) trilobite *Proceratopyge* cf. *P. Lata* Whitehouse from northern Victoria Land, Antarctica, and the evolution of metamorphosis in trilobites. *Palaeontology* **59**(5), 657–670 (2016).

37. Laibl, L. & Fatka, O. Review of early developmental stages of trilobites and agnostids from the Barrandian area (Czech Republic). *J. Natl. Mus. (Prague) Nat. Hist. Ser.* **186**(1), 103–112 (2017).

38. Laibl, L., Cederström, P. & Ahlberg, P. Early post-embryonic development in *Ellipostrenua* (Trilobita, Cambrian, Sweden) and the developmental patterns in Ellipsocephaloidea. *J. Paleontol.* **92**(6), 1018–1027 (2018).

39. Laibl, L., Maletz, J. & Olschewski, P. Post-embryonic development of *Fritzolenellus* suggests the ancestral morphology of the early developmental stages in Trilobita. *Pap. Palaeontol.* **7**(2), 839–859 (2021).

40. Chatterton, B. D. E. & Speyer, S. E. Ontogeny in *Treatise on Invertebrate Paleontology*. Part O, Arthropoda 1, Trilobita 1, revised. 7–11 (Geological Society of America and University of Kansas Press, Lawrence, 1997).

41. Bignon, A., Waisfeld, B. G., Vaccari, N. E. & Chatterton, B. D. Reassessment of the order Trinucleida (Trilobita). *J. Syst. Palaeontol.* **18**(13), 1061–1077 (2020).

42. Torsvik, T. H. & Cocks, L. R. M. The Palaeozoic palaeogeography of central Gondwana. *Geol. Soc. Lond. Spec. Publ.* **357**(1), 137–166 (2011).

43. Torsvik, T. H. & Cocks, L. R. M. New global palaeogeographical reconstructions for the early Palaeozoic and their generation. *Geol. Soc. Lond. Mem.* **38**(1), 5–24 (2013).

44. Bahlburg, H., Moya, M. C. & Zeil, W. Geodynamic evolution of the early Palaeozoic continental margin of Gondwana in the Southern Central Andes of Northwestern Argentina and Northern Chile. In *Tectonics of the Southern Central Andes*. 293–302 (Springer, 1994).

45. McEdward, L. R. & Miner, B. G. Larval and life-cycle patterns in echinoderms. *Can. J. Zool.* **79**(7), 1125–1170 (2001).

46. Lefebvre, B. et al. Palaeobiogeography of Ordovician echinoderms. *Geol. Soc. Lond. Mem.* **38**(1), 173–198 (2013).

47. Signor, P. W. & Vermeij, G. J. The plankton and the benthos: Origins and early history of an evolving relationship. *Paleobiology* **20**, 297–319 (1994).

48. Davis, M. A., Grime, J. P. & Thompson, K. Fluctuating resources in plant communities: A general theory of invasibility. *J. Ecol.* **88**(3), 528–534 (2000).

49. Franeck, F. Perspectives on the Great Ordovician Biodiversification Event-local to global patterns (2020).

50. Pulliam, H. R. Sources, sinks, and population regulation. *Am. Nat.* **132**(5), 652–661 (1988).

51. Kröger, B., Franeck, F. & Rasmussen, C. M. The evolutionary dynamics of the early Palaeozoic marine biodiversity accumulation. *Proc. R. Soc. B* **286**(1909), 20191634 (2019).

52. Penny, A. & Kröger, B. Impacts of spatial and environmental differentiation on early Palaeozoic marine biodiversity. *Nat. Ecol. Evol.* **3**(12), 1655–1660 (2019).

53. Rasmussen, C. M., Kröger, B., Nielsen, M. L. & Colmenar, J. Cascading trend of Early Paleozoic marine radiations paused by Late Ordovician extinctions. *Proc. Natl. Acad. Sci.* **116**(15), 7207–7213 (2019).

54. Stigall, A. L. The invasion hierarchy: Ecological and evolutionary consequences of invasions in the fossil record. *Annu. Rev. Ecol. Evol. Syst.* **50**, 355–380 (2019).

55. Stigall, A. L., Edwards, C. T., Freeman, R. L. & Rasmussen, C. M. Coordinated biotic and abiotic change during the Great Ordovician Biodiversification Event: Darriwilian assembly of early Paleozoic building blocks. *Palaeogeogr. Palaeoclimatol. Palaeoecol.* **530**, 249–270 (2019).

56. Stigall, A. L. How is biodiversity produced? Examining speciation processes during the GOBE. *Lethaia* **51**(2), 165–172 (2018).

57. Servais, T. & Harper, D. A. T. The great Ordovician biodiversification event (GOBE): Definition, concept and duration. *Lethaia* **51**(2), 151–164 (2018).

58. Trotter, J. A. et al. Did cooling oceans trigger Ordovician biodiversification? Evidence from conodont thermometry. *Science* **321**(5888), 550–554 (2008).
59. Vizcaíno, D., Álvaro, J. J. & Lefebvre, B. The lower Ordovician of the southern Montagne Noire. *Ann. Soc. Géol. Nord* 8(4), 213–220 (2001).

60. Hsieh, T. C., Ma, K. H. & Chao, A. iNEXT: An R package for rarefaction and extrapolation of species diversity (Hill numbers). *Methods Ecol. Evol.* 7(12), 1451–1456 (2016).

61. Suchéras-Marx, B., Escarguel, G., Ferreira, J. & Hammer, Ø. Statistical confidence intervals for relative abundances and abundance-based ratios: Simple practical solutions for an old overlooked question. *Mar. Micropaleontol.* **151**, 101751 (2019).

62. Hammer, Ø., Harper, D. A. T. & Ryan, P. D. PAST-palaeontological statistics, ver. 1.89. *Palaeontol. Electron* **4**(1), 1–9 (2001).

63. Gibert, C. & Escarguel, G. PER-SIMPER—A new tool for inferring community assembly processes from taxon occurrences. *Glob. Ecol. Biogeogr.* **28**(3), 374–385 (2019).

64. Gibert, C. DNCImper: *Assembly Process Identification Based on SIMPER Analysis*. R package ver. 0.0.1.0000. https://github.com/ Corentin-Gibert-Paleontology/DNCImper (2019).

## Acknowledgements

The authors thank Alycia Stigall and Christian Klug for their reviews and Lukáš Laibl and Francesc Pérez-Peris for fruitful discussions. This paper is supported by grant no. 2020M683388 from the Chinese Postdoctoral Science Foundation awarded to F. Saleh and by grant no. 205321_179084 from the Swiss National Science Foundation, awarded to A.C.D. as principal investigator. D.B., B.G.W. and F. Serra acknowledge supporting facilities from the Consejo Nacional de Investigaciones Científicas y Técnicas (PUE 2016 CONICET -CICTERRA) and grant PICT 2016 0588 provided by the Agencia Nacional de Promoción Científica y Tecnológica to B.G.W. Research by M.G.M. and L.A.B. was supported by Natural Sciences and Engineering Research Council (NSERC) Discovery Grants 311727-20 and 422931-20, respectively. M.G.M. thanks funding by the George J. McLeod Enhancement Chair in Geology. This paper is a contribution of the IGCP653 and IGCP 735 programs.

## Author contributions

F Saleh and BL designed the project. F Saleh, BL, DB, BGW, F Serra, DV provided raw data. F Saleh, PG and CG ran the analyses. All authors interpreted and discussed the results. F Saleh wrote the manuscript and made the figures with the help of all coauthors.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-022-07822-z.

**Correspondence** and requests for materials should be addressed to F.S., C.G. or X.M.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_7.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2022