Ecology

Hydrological changes in high-latitude peatlandsWe observed three hydrological pathways, i.e., drying, wetting, and fluctuating, for both peatland clusters, non-permafrost and permafrost peatlands (Fig. 2). Approximately 54% of the studied peatlands have shifted towards drier surface conditions since 1800 CE and more intensively since 1900 CE (Fig. 2a, d), which is in line with the post-LIA warming. The overall change point to drier conditions was dated to ca. 1950 CE for non-permafrost sites and ca. 1890 CE for the permafrost compiled group. Approximately 32% of the studied peatlands have shifted towards wetter conditions (Fig. 2b, e). The overall shifting point to wetter conditions was dated to ca. 1995 CE for non-permafrost peatlands and to ca. 1990 CE for permafrost peatlands. Wetting has been especially intensive since 1900 CE for non-permafrost peatlands and since 1950 CE for permafrost peatlands. Interestingly, the data showed that in permafrost peatlands a significant dry shift always preceded a wet shift (Fig. 2e). Approximately 14% of the studied peatlands indicated no clear trend, with fluctuating hydrological conditions (Fig. 2c, f).Non-permafrost peatlands generally showed spatially extensive drying across the northern high latitudes, apart from northeastern Canada, where a wetting trend was more frequently observed. Permafrost peatlands, however, were more variable, with some drying, some wetting, and no overall coherent regional pattern was visible (Fig.1a, b). It should be noted that peatlands synthesized here have undergone little or no direct human impact, i.e., their surface hydrology was not significantly affected by human disturbances such as drainage, when compared to, for example, central European peatlands discussed in Swindles et al. (2019)5. This implies that climate and/or local topographical forcing are the predominant hydrological drivers in this study. The dataset is to some extent biased as there are more non-permafrost records from northeastern Canada but more permafrost records from northern Sweden and this might result in regional overestimation to either wetting or drying trends. Nevertheless, the pattern of diverse timing of the hydrological shifts between the individual coring points (Fig. 2) indicates the variability in sensitivity of different regions/peatlands to climate changes.Potential links to climate change and permafrost dynamicsThe comparison of the reconstructed water table and climate data suggests that climate, especially summer temperature, has played an important role in shaping the peatland water table (Fig. 1c–f). The pattern detected here for non-permafrost peatlands, an extensive drying, is comparable to that observed for central European peatlands5. In addition to direct climate forcing, a recent acceleration of peat accumulation might partly explain the drying trend by disconnecting the peatland surface from the water table17. However, for northeastern Canada a wetting trend has been observed more often, possibly regulated by the regional climate that shows clearly less warming in the focused period compared to other regions (Fig. 1c, d).Permafrost initiation in the past caused a peat surface uplift and is probably detected as a dry phase (Fig. 2d, e). Post-LIA warming-induced increase in evapotranspiration may have strengthened the surface drying which originally resulted from surface-uplift and probably mitigated the gradual wetting related to permafrost thawing11. The level of warming has varied among the regions. In some areas such as northeastern Canada temperature has increased less and, when combined with higher precipitation or higher effective moisture level, may have caused surface wetting in permafrost peatlands. This is a direct climate forcing rather than permafrost thawing, which is a consequence of climate warming, i.e., more indirect climate forcing. To date, it is yet challenging to estimate any one tipping point of warming that might trigger permafrost thawing, as the local conditions vary from bottom ground soil conditions to hydrology and vegetation. The consequent wetting or drying depends on evapotranspiration and ice richness etc., which further challenges the prediction of hydrological conditions of permafrost peatlands.The divergent three moisture patterns may occur in the same region and even in the same peatland, especially if the permafrost is present. This complex response pattern is well supported by the records from the Abisko region, Sweden, where replicated sampling was carried out, and captured different successional stages of local permafrost peatlands7,18. In contrast to permafrost peatlands, non-permafrost peatlands are more likely to experience a more consistent ecosystem response pattern19 as supported by the replicated records suggesting the same pattern happening simultaneously in several regions (Fig. 1a). The fluctuating pattern of many records reported here suggests that the past and recent climate has not yet caused a state change in hydrological conditions.Insights into carbon dynamics and future perspectiveGenerally, our results suggest that the recent climate warming has caused hydrological shifts in most high-latitude peatlands, highlighting its pronounced effect on shaping peatland moisture balance, and further on driving peatland C balance. It has been reported that a 1-cm water-table drawdown would increase 3.3–5.0 mg CO2 m−2 h−1 and decrease 2.2–3.6 mg CO2-eq m−2 h−1 (CH4) to the atmosphere, and the average sensitivity of CO2 and CH4 combined was 0.8–2.3 mg CO2-eq m−2 h−1 cm−1 according to a global scale analysis, including sites from high-latitudes3. However, it should be noted that the sensitivity of greenhouse gas fluxes to the magnitude of hydrological changes might vary among different regions and peatland types. It appears that most pan-Arctic peatlands are undergoing a drying trend, that may lead to a decreased C sink capacity3,19, if not compensated by increased C uptake from the atmosphere20.It is very likely that over the 21st century warming in high latitudes will continue to be more pronounced than the global average21. Precipitation is projected to increase, albeit with large regional variability. Also, extreme events with heavy rainfall and drought are becoming more frequent and intense22. It is estimated that about 20% of permafrost zone is experiencing accelerated and abrupt permafrost thaw that is likely causing wetting conditions4, while gradual permafrost thaw has been observed across the circumpolar regions23. Both an increase in precipitation and permafrost thaw might mitigate the drying pressure caused by warming and increased evapotranspiration. However, abrupt permafrost thaw in peatlands can result in a rapid (over years to decades) loss of C from the formerly frozen permafrost peat, causing these peatlands to be a net source of C to the atmosphere before post-thaw accumulation returns them to a net sink (centuries to millennia)12,13. The future C sink and source function of peatlands is a key element in contributing to climate change, but the observed divergent pathways of peatland hydrological successions further challenge the projections of high-latitude peatland C sink and source dynamics. Conversely, it clearly highlights the importance of climate forcing in peatland succession scenarios. Our study reveals that the response of high-latitude peatlands to changing climate conditions is complex. We detect variable ecohydrological trajectories, and in the future, these will determine the C sink capacity of northern peatlands. The observed patterns inevitably create challenges for the climate change modelling community. How to capture the highly heterogenic successional pathways of northern peatlands needs to be a key research focus. More

To assess coral holobiont health, we used the commercially available Urinalysis strips (Accutest URS-14 100 strips Urinalysis Reagent Test Paper; ca. $11–35 USD/ bottle) that are designed to detect disease markers in humans, potentially indicating diabetes, metabolic abnormalities, liver diseases, kidney function, and urinary tract infections. There are ten tests on each strip that provide an initial assessment of health using the following standardized markers: leukocytes, nitrite, urobilinogen, protein, pH, blood, specific gravity, ketone, bilirubin, and glucose. To apply this test to coral extracts, we used ambient and stressed nubbins in a comparative approach to identify trends in the data. We presumed that some of the tests are likely not applicable to corals, whereas others such as ketones and leukocytes, which target conserved animal pathways (see below) could prove useful. Samples from three different coral experiments (two bleaching experiments and one environmental survey) were used to assess the response and utility of the different tests on each strip. The first bleaching experiments was done in 2019 on two Hawaiian species: Montipora capitata (stress resistant) and Pocillopora acuta (stress sensitive); this research has been previously described in detail11,12. Briefly, three M. capitata colonies (different genotypes) were selected, with samples collected from each in triplicate at three timepoints (referred to as T1, T3, and T5) that span a period of prolonged thermal stress (16 days). The second bleaching experiment was done in 2021 and included nubbins from M. capitata, P. acuta, and P. compressa. These corals were maintained for 9 days under heat treatment or ambient conditions, with samples collected at the beginning and end of the experimental period from both conditions (see Methods). Finally, the third experiment was done in June 2021 on wild colonies of M. capitata, P. acuta, and P. compressa (Fig. 1A) from Kāneʻohe Bay, Oʻahu, Hawaiʻi (Fig. 1B) and analyzed to assess the extent of natural variability in the test strip results. It should be noted that the test strip results from each experiment are independent (i.e., analyzed separately to measure relative differences within the test population) due to differences in tissue freezing time, storage length, and handling, that affect the results.Figure 1Analysis of Hawaiian corals. (A) The three targeted species in Kāneʻohe Bay. Images created by D. Bhattacharya. (B) Sites of wild coral collection (marked by the yellow circles) in Kāneʻohe Bay in June 2021. Nubbins from two colonies (n = 3) of M. capitata were collected from six reefs: Reef 8, 9, Reef 13, Reef 27, Reef 41, and Reef 43 and from three sites near the Hawaiʻi Institute of Marine Biology (HIMB) on Coconut Island (Moku o Loʻe). This image was adapted from the Pacific Islands ocean observing system (https://www.pacioos.hawaii.edu/projects/coral/)33.Full size imagePortable device for test strip analysisTo facilitate field testing and allow accurate measurement of test strip results in the R, G, and B channels, a combination of 3D printing technology and computer vision was applied to the problem. A custom opaque phone holder was produced using 3D printing (Ultimaker Inc. [https://ultimaker.com]) to control light source quality and quantity during the procedure and for calculating test strip scores (Fig. 2A). This portable phone holder was designed to allow field personnel to conduct experiments and collect results in a convenient and user-friendly manner. Inside the black opaque phone holder, two LED diodes served as a controllable and consistent light source to eliminate noise introduced by any change in light levels (between or during) the test strip reactions, allowing for more accurate and consistent data collection19. Frames of the videos (taken using the smart phone positioned in the holder) were captured for different reagents according to the recommended reading time given by the user manual. For test strip RGB value measurement (Fig. 2B), an automated machine learning method, shown in Fig. 2C, was used as a replacement for traditional, manual methods. Our computer vision workflow, named TestStripDX, uses YOLOv418, which is one of the most mature, accurate, and popular computer vision models available, to isolate a target feature in an image and provide the metadata associated with the object (in this case, each test along the strip), such as position within the coordinates of a rectangular box. Sample pictures of strips were manually annotated and used for training the TestStripDX pipeline (Supplementary Data 1): i.e., these images were used to train a custom detector to identify the different reagent tests along a strip and thereafter, to provide the relevant RGB values for each test (Fig. 2C). Figure 2D shows the comparison between data collected using the manual (ImageJ) versus automated (TestStripDX) methods for a partial list of samples and reagent strips (data shown in Table S1). The very close association of the regression line (R2 = 0.996) with the artificial diagonal line validates the utility and accuracy of the computer vision method.Figure 2Portable instrument for test strip analysis. (A) The 3D-printed phone holder used to analyze the test strips. (B) Example test strip. (C) Flow chart for the training and application of the machine learning method (TestStripDX) used to analyze test strips. (D) Comparison of test strip colorimetric measurement values produced by TestStripDX and ImageJ, for a set of representative images (Table S1). The enforced diagonal matches very closely to the regression line.Full size imageCoral color score measurementsAs described above, color scores provide a proxy for coral health, and are generated by measuring the RGB channels of the bleached (unhealthy) and brown (healthy: i.e., pigmentation provided by the algal symbionts) areas of coral nubbins. For the CoralDX workflow, we trained a custom detector that can recognize coral nubbins, as well as red, green, and blue colored blocks (standards) which are used to normalize the R, G, and B channels in each image (see Supplementary Data 2). The images taken in this case were from a lab environment (Fig. 3A), but this approach could potentially be used in any location provided a background with a uniform color is used in the image. To achieve this goal, a background panel containing the red, blue, and green color blocks would be placed behind (smaller) coral nubbins to make the measurements. It is clear that for large colonies, this approach may prove challenging to apply but we expect that additional testing and modification to the method will allow us to design a better-suited tool for field use. To achieve our goal of automatically obtaining colorimetric measurements of coral nubbins, we needed to devise a method that provides YOLOv4 with areas that are not limited by the standard rectangular boxes used as input for this model. To accommodate irregular nubbin shapes, an additional step (training and testing steps are shown in Fig. 3B) was added to the automated method (Fig. 2C). Our approach uses computer vision-based edge detection to eliminate most of the background surrounding the edges of the coral nubbins, allowing for accurate quantification of the RGB values of the targeted piece of coral. After edge detection, we obtain a picture with a black background highlighting the coral shape as a “mask” (Fig. 3C; 2nd image from left). We then place the mask onto the original coral image, measuring non-zero R, G, and B values (Fig. 3C; 3rd and 4th images from left). This method is superior to selecting coral areas for manual analysis using tools such as the handsfree selection function in ImageJ, which are generally hard to manipulate, have low fault tolerance, and are time-consuming. The RGB values extracted from the coral nubbins and color blocks were investigated using principal component analysis (PCA) to generate Euclidean distances (color scores) among coral nubbins according to treatment, time point, and colony11,16. In Fig. 3D, the correlation between bleaching scores generated using CoralDX and ImageJ for M. capitata and P. acuta are presented for a representative set of nubbins cultured under ambient or thermal stress conditions in the 2019 bleaching experiment (Tables S2, S3). As is apparent, the regression lines show strong correlations (R2 = 0.968 for M. capitata and R2 = 0.991 for P. acuta), supporting the utility of the automated method. The data in both cases are very closely associated with the artificial diagonal line shown in the images, substantiating the strong positive correlation between the scores from the two approaches. The small differences in score values between the methods is primarily explained by difficulties in cropping the edges of coral nubbins that are heavily bleached; the contrast between the coral nubbin and the white background is lessened, creating discrepancies between the nubbin edges identified by CoralDX and ImageJ. This issue can be addressed by testing different color backgrounds to find the optimal set-up. Another potential contributing factor is the differences in the area selected for the color blocks between manual (ImageJ) versus automated methods, although this might only have a minor effect on the results. The CoralDX workflow is not computationally intensive and easily portable to other platforms, such as personal computers and smart phones, allowing for easy deployment in the field.Figure 3Analysis of coral color scores. (A) Example image of M. capitata coral nubbins used for automated color score analysis, with target areas marked with the black (coral nubbins) and yellow (color standards) boxes. The nubbin on the far right is used to demonstrate the masking procedure. (B) Flow chart for the training and application of the machine learning method (CoralDX) used to analyze coral color scores. (C) Example of image processing for one M. capitata coral nubbin showing (from left to right) the original, masked, segmented, and masked RGB nubbin, with the final (right) image used for color score measurements. (D) Comparison of color score values produced by CoralDX and ImageJ, for a set of representative coral nubbins from the 2019 bleaching experiment (Tables S2, S3). The enforced diagonal matches very closely to the regression line in both analyses.Full size imageTest strip resultsArmed with these new tools, we generated the test strip data for different coral samples, where larger “relative enzymatic activity” (REA, measured using the RGB values) values indicate increased reactivity (i.e., increased levels of products targeted by the test). The Accutest URS-14 100 strips test for ketones (using the sodium nitroprusside reaction) measures acetoacetate and assumes the presence of β-hydroxybutyrate and acetone. The former (β-hydroxybutyrate) acts as a signal to regulate metabolism and maintain energy homeostasis during nutrient deprivation. In this process, β-hydroxybutyrate is converted to acetoacetate. Ketone bodies are transported into tissues and converted into acetyl-CoA by thiolases, which then enters the TCA cycle and is oxidized in the mitochondria for energy. Bleaching in corals which are incapable of obtaining adequate energy stores through heterotrophy results in diminished growth rates, degraded reproductive capacity, amplified susceptibility to disease, and elevated mortality rates for the entire colony20. Although ketosis has not been explored in cnidarians, transcriptomic data generated from the M. capitata samples measured in this study12 demonstrate expression of the KEGG pathway for degradation of ketone bodies (Fig. S1). The combination of time-point and treatment (field, T1-AT, T1-HT, T3-HT, T5-HT) was the most significant factor impacting the ketone REA scores (p-value = 0.010) (see details of the PERMANOVA analysis in the “Methods” section). Given this framework, we find that the M. capitata samples remain steady throughout the bleaching period, except for a decrease in enzymatic activity at T3-HT (Fig. 4A). This result is supported by the transcriptomic data, which shows this pathway to be uniformly expressed at all timepoints, except T3-HT, when acetyl-CoA C-acetyltransferase is up-regulated in comparison to T1-HT (fold change [FC] = 1.52) and down-regulated at T5-HT (FC = − 1.84) (Table S4). A possible explanation for this result is that at T3-HT, when the first evidence of bleaching was present, the photosynthetic rate of the symbiotic algae was elevated due to the thermal stress, resulting in greater energy production and a decrease in the abundance of ketone bodies within the coral, but without sufficient stress to cause significant bleaching. During this time, acetyl-CoA C-acetyltransferase enzymatic activity favored the production of acetoacetyl-CoA. However, as dysbiosis continued and the corals no longer had access to the photosynthetic products provided by the symbionts, they produced more acetoacetate and ketosis was detected by the test strips. This response occurs despite the observation that M. capitata can persist heterotrophically and meet much of its energy needs in the absence of algal symbionts21. Interestingly, the three M. capitata field samples show similar amounts of ketone bodies to lab stressed corals, suggesting the presence of stressors in the natural environment.Figure 4Test strip results from the 2019 bleaching experiment. (A) Ketone test strip results for M. capitata, showing genotype-specific (see legend) differences in response. (B) Leukocytes test strip results for M. capitata, showing genotype-specific (see legend) differences in response. These are standard box plots, with the boxes representing the first (Q1) to third (Q3) quartiles. The lines in the boxes are the median (Q2) values and lines (“whiskers”) extending beyond the boxes are the minimum and maximum values, excluding outliers. (C) PCA of the ketone and leukocytes test strip data for the field, ambient, and T1-HT, T3-HT, and T5-HT timepoints.Full size imageThe leukocyte test measures the activity of leukocyte esterase (presence of white blood cells) and other signs of infection in human subjects. This test putatively assesses the coral innate immunity response, which includes the same phases in response to infection and loss of tissue integrity as other invertebrates: recognition, signaling, and effector response22. Corals contain multiple types of immune cells, such as amoebocytes and fibroblasts23. Amoebocytes are amoeboid cells residing in the mesoglea that remove necrotic tissue, encapsulate foreign particles, and generally display phagocytic activity to aid in organism defense against pathogens, which is the cnidarian principal mechanism of immunity24. Amoebocytes can be melanin-containing, agranular, or granular based cells, depending on the signaling pathway22. M. capitata shows an overall increase in enzymatic activity at T1-HT, signaling a heightened immune response (Fig. 4B). This relatively higher level of enzymatic activity decreases at T3 and T5. The field samples show levels that are comparable to the T3-HT and T5-HT thermal stress acclimated colonies and the T1 ambient (Amb) colonies. Again, a combination of time-point and treatment was the most significant factor impacting the leukocyte REA scores (p-value = 0.001). Additional research needs to be done to understand the cause of this response. Nonetheless, the results are consistent with the widely accepted hypothesis that M. capitata adapts well to bleaching conditions12. PCA of the M. capitata test strip results shows a clear separation of the T1-HT coral data from the T3-HT and T5-HT values, with the ambient and field samples intermixed among the latter two timepoints (Fig. 4C). This result again highlights the initial robust response to stress by M. capitata followed by acclimation to the heat treatment that is reflected in the field samples.A noteworthy aspect of the test strip results is the divergence in response to thermal stress among different colonies. This result has also been found for coral metabolomic data11. As described above, each holobiont integrates a complex set of biotic interactions between the host animal and microbiome, explaining the high variation in ketone and leukocytes test results, often between replicate nubbins from one colony and more frequently, between different colonies (Fig. 4A). Existing data using omics methods demonstrate that the stress response of the coral holobiont varies from colony to colony15. The metabolome is controlled by the coral animal genotype, microbial consortium, and environmental conditions, among other factors, and can fluctuate greatly based on individual metabolite turnover rates and the timing of sampling25,26. Therefore, accounting for natural variation in the stress response phenotype and its importance for effective testing methods is a crucial aspect of our work. Our results demonstrate that broad population level sampling (dozens to 100 s of colonies/genotypes) is likely needed to account for the inherent genetic and metabolic variation present in wild coral populations.We also did a more limited analysis of coral stress responses using the ketone and leukocyte test strips with three species (M. capitata, P. acuta, P. compressa) in a 2021 bleaching experiment in which we sampled multiple coral genotypes at time 0 and after 9 days of heat treatment (3ºC increase from ambient; see Methods). These results are based on analysis of 8–9 different coral genotypes (summarized in Fig. 5). Interpreted in the same way as described above, we see that there is substantial genotype-based variation in the stress response. Nonetheless, consistent with the 2019 data, M. capitata shows evidence of a thermal stress response in the ketone and leukocytes tests (Fig. 5A, B). P. acuta shows a more limited response, whereas P. compressa appears to have fully acclimated to the stress regime with lowered reactivity at the end of the experiment. Species identity was the most significant factor impacting the leukocyte and ketone REA scores (p-value = 0.002 and 0.033, respectively), but a combination of species and treatment (AT vs. HT) was found to be significant for leukocytes (p-value = 0.026). PCA of the M. capitata test results shows separation between the ambient samples and some of the high temperature treated samples along PC1; both the ketone and leukocytes tests contribute to the spread of samples along PC1. This reinforces the conclusion that there is evidence for a thermal stress response in the ketone and leukocytes tests of M. capitata however, the significant genotypic variability (particularly apparent when compared with Fig. 4) obscures the differences between the conditions for some of the samples. PCA of the P. acuta and P. compressa test results mirrors our conclusion that the ketone and leukocytes tests in these species show a limited response to stress, with no separation between the samples from the different conditions.Figure 5Test strip results from the 2021 bleaching experiment for three Hawaiian coral species. (A) Ketone test strip results for M. capitata, P. acuta, and P. compressa showing genotype-specific variation in response under the ambient (Amb) and high temperature (HT) treatments. (B) Leukocytes test strip results for M. capitata, P. acuta, and P. compressa showing genotype-specific variation in response under the ambient and high temperature (HT) treatments. These are standard box plots, with the boxes representing the first (Q1) to third (Q3) quartiles. The lines in the boxes are the median (Q2) values and lines (“whiskers”) extending beyond the boxes are the minimum and maximum values, excluding outliers. (C) PCA of the ketone and leukocytes test strip data for the Amb and HT treatments for the three Hawaiian coral species.Full size imageAnalysis of wild populationsTo assess natural variation, we collected apparently healthy M. capitata, P. acuta, and P. compressa nubbins (3 replicate nubbins per colony) from six reefs in Kāneʻohe Bay and from three sites near the Hawaiʻi Institute of Marine Biology on Coconut Island (Moku o Loʻe) (see Fig. 1) and analyzed these tissue extracts using the ketone and leukocytes tests. This analysis shows wide variation in the results with some interesting exceptions. Species identity was the most significant factor impacting the leukocytes and ketone REA scores (p-value = 0.002 and 0.001, respectively), but for leukocytes, colony identity, regardless of species, was also found to be a significant factor (p-value = 0.040). The M. capitata ketone test results are consistent among different reefs and within the same colony with the exception of some colonies (e.g., Colony 8 from Reef 9 and Colony 16 from Reef 43) that show wide intra-colony variation (Fig. 6A). Most of the ketone data for M. capitata fall between REAs of 10–20. In contrast, P. acuta shows more variation in the wild populations for the ketone test, suggesting that many of these coral colonies live under stressful conditions in the field (Fig. 6B). A similar situation to M. capitata, in terms of REAs, holds for the ketone test of P. compressa colonies that show more limited variation (Fig. 6C). The leukocytes test shows high variation for M. capitata (Fig. 6D), P. acuta (Fig. 6E), and P. compressa (Fig. 6F) colonies. These results again point out the complex nature of genome-environment interaction with respect to metabolic syndromes, both at the colony level and among different regions (replicates) of the same colony. For example, P. compressa Colonies 13–15 from Reef 9 show little to no intra-colony variation for the ketone test, yet another colony from this reef (Colony 16) shows high variation among replicates (Fig. 6C). In contrast, P. compressa Colonies 13–15 are far more variable when using the leukocytes test (Fig. 6F). On the basis of the more predictable lab-based results reported above, we interpret these “noisy” field data as evidence of the immense variation in the stress phenome of wild coral populations. Overall, the field results indicate that a starting set of test strip values, followed by repeated field sampling over time of wild colonies may be needed for accurate stress diagnosis, rather than the one-time measurement approach used here. Clearly, more work is needed with wild colony analysis, particularly under varying degrees of thermal stress and apparent bleaching, to fully realize the potential of the technique we present.Figure 6Test strip results from analysis of the 2021 collection of three Hawaiian coral species from the wild. (A and D) Ketone and leukocytes test strip results, respectively, for M. capitata showing intraindividual variation and among the different reefs (see legend) that were sampled (Fig. 1B). (B and E) Ketone and leukocytes test strip results, respectively, for P. acuta showing intraindividual variation and among the different reefs that were sampled. (C and F) Ketone and leukocytes test strip results, respectively, for P. compressa showing intraindividual variation and among the different reefs (see legend) that were sampled. These are standard box plots, with the boxes representing the first (Q1) to third (Q3) quartiles. The lines in the boxes are the median (Q2) values and lines (“whiskers”) extending beyond the boxes are the minimum and maximum values, excluding outliers.Full size imageAnalysis of transcriptomic dataTo identify pathways that may support the M. capitata leukocytes test strip results, which showed the most response in terms of change in REA at T1-HT (Fig. 4B), we analyzed existing transcriptomic (RNA-seq) data derived from the same coral nubbins. The RNA-seq and metabolomic data from these samples have been previously analyzed11,12. Here we searched for co-expression modules that contain genes that are up-regulated at the start of the thermal stress regime (T1-HT) when compared to the ambient treatment. It is at this timepoint that we find a strong cross-reaction with the leukocytes test, followed by loss of cross-reactivity at T3-HT and T5-HT (back down to T1-Amb levels), putatively indicating acclimation (Fig. 4B). As described above, the wound healing response in corals is complex and the (Urinalysis) leukocytes results need to be interpreted as a syndrome involving multiple pathways of stress and immune response. With these considerations in mind, we identified a module (Module 2; see Williams et al.12) of up-regulated genes that contains several markers associated with the coral stress response (Fig. 7). These include a tumor necrosis factor-activated receptor (TNFR)-Cys domain-containing protein (fold-change [FC] = 1.12) that is a well-known mediator of apoptosis and cell death that is functionally conserved in corals. Some members of the TNF family are associated with bleaching27. The most highly up-regulated gene in this module is C-type lysozyme 2 (FC = 2.42) that provides an anti-microbial function (e.g., digestion of peptidoglycan), and is likely expressed as a result of stress-induced dysbiosis in M. capitata. Other markers of stress that are up-regulated in Module 2 include E3 ubiquitin-protein ligase (FC = 1.03) involved in protein degradation, two protein disulfide-isomerase (FCs = 1.32, 1.09) involved in cellular defense against protein misfolding via chaperone activity28, and a metalloproteinase inhibitor 3 (FC = 1.42) which likely functions as a physiological anti-inflammatory molecule29. These data, although not directly substantiating the leukocytes test strip results, provide evidence that the wound healing and immune response pathways were up-regulated in Module 2 at T1-HT (albeit weakly, due to the stress resilience of M. capitata) as suggested by Fig. 4B.Figure 7Gene co-expression analysis. Module 2 representing significantly up-regulated genes in M. capitata from T1-HT in the 2019 differential gene expression analysis12. This module is enriched in genes involved in the wound healing and immune response. The legend for level of up-regulation is shown. Putative gene annotations are also shown.Full size image More

Erwin, D. H. & Tweedt, S. Ecological drivers of the Ediacaran-Cambrian diversification of Metazoa. Evol. Ecol. 26, 417–433 (2012).Article 

Google Scholar 
Laflamme, M., Darroch, S. A., Tweedt, S. M., Peterson, K. J. & Erwin, D. H. The end of the Ediacara biota: Extinction, biotic replacement, or Cheshire Cat?. Gondwana Res. 23, 558–573 (2013).ADS 
Article 

Google Scholar 
Mángano, M. G. & Buatois, L. A. Decoupling of body-plan diversification and ecological structuring during the Ediacaran-Cambrian transition: Evolutionary and geobiological feedbacks. Proc. R. Soc. B-Biol. Sci. 281, 20140038 (2014).Article 

Google Scholar 
Mángano, M. G. & Buatois, L. A. The Cambrian revolutions: Trace-fossil record, timing, links and geobiological impact. Earth-Sci. Rev. 173, 96–108 (2017).ADS 
Article 

Google Scholar 
Mángano, M. G. & Buatois, L. A. The rise and early evolution of animals: Where do we stand from a trace-fossil perspective?. Interface Focus 10, 20190103 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
Darroch, S. A. et al. Biotic replacement and mass extinction of the Ediacara biota. Proc. R. Soc. B 282, 20151003 (2015).PubMed 
PubMed Central 
Article 

Google Scholar 
Darroch, S. A., Smith, E. F., Laflamme, M. & Erwin, D. H. Ediacaran extinction and Cambrian explosion. Trends Ecol. Evol. 33, 653–663 (2018).PubMed 
Article 

Google Scholar 
Schiffbauer, J. D. et al. The latest Ediacaran Wormworld fauna: Setting the ecological stage for the Cambrian explosion. GSA Today 26, 4–11 (2016).Article 

Google Scholar 
Zamora, S., Deline, B., Javier Álvaro, J. & Rahman, I. A. The Cambrian Substrate Revolution and the early evolution of attachment in suspension-feeding echinoderms. Earth-Sci. Rev. 171, 478–491 (2017).ADS 
Article 

Google Scholar 
Hantsoo, K. G., Kaufman, A. J., Cui, H., Plummer, R. E. & Narbonne, G. M. Effects of bioturbation on carbon and sulfur cycling across the Ediacaran-Cambrian transition at the GSSP in Newfoundland, Canada. Can. J. Earth Sci. 55, 1240–1252 (2018).ADS 
CAS 
Article 

Google Scholar 
Boyle, R. A., Dahl, T. W., Bjerrum, C. J. & Canfield, D. E. Bioturbation and directionality in Earth’s carbon isotope record across the Neoproterozoic-Cambrian transition. Geobiology 16, 252–278 (2018).CAS 
PubMed 
Article 

Google Scholar 
van de Velde, S., Mills, B. J., Meysman, F. J., Lenton, T. M. & Poulton, S. W. Early Palaeozoic ocean anoxia and global warming driven by the evolution of shallow burrowing. Nat. Commun. 9, 2554 (2018).ADS 
PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Gougeon, R. C., Mangano, M. G., Buatois, L. A., Narbonne, G. M. & Laing, B. A. Early Cambrian origin of the shelf sediment mixed layer. Nat. Commun. 9, 1909 (2018).ADS 
PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Seilacher, A. & Pflüger, F. In Biostabilization of Sediments (eds. Krumbein W. E., Peterson D. M., & Stal L. J.) 97–105 (Bibliotheks und Informationsystem der Carl von Ossietzky Universität, 1994).Seilacher, A. Biomat-related lifestyles in the Precambrian. Palaios 14, 86–93 (1999).ADS 
Article 

Google Scholar 
Buatois, L. A. et al. Quantifying ecospace utilization and ecosystem engineering during the early Phanerozoic: The role of bioturbation and bioerosion. Sci. Adv. 6, eabb0618 (2020).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Hagadorn, J. W. & Bottjer, D. J. Restriction of a late Neoproterozoic biotope; suspect-microbial structures and trace fossils at the Vendian-Cambrian transition. Palaios 14, 73–85 (1999).ADS 
Article 

Google Scholar 
Buatois, L. A. & Mangano, M. G. Early colonization of the deep sea: Ichnologic evidence of deep-marine benthic ecology from the Early Cambrian of northwest Argentina. Palaios 18, 572–581 (2003).ADS 
Article 

Google Scholar 
Clausen, S., Álvaro, J. J. & Zamora, S. Replacement of benthic communities in two Neoproterozoic-Cambrian subtropical-to-temperate rift basins, High Atlas and Anti-Atlas, Morocco. J. Afr. Earth Sci. 98, 72–93 (2014).ADS 
Article 

Google Scholar 
Mángano, M. G. & Buatois, L. A. In The Trace-Fossil Record of Major Evolutionary Events: Volume 1: Precambrian and Paleozoic (eds. Mángano M. G. & Buatois L. A.) 73–126 (Springer Netherlands, 2016).Bayet-Goll, A., Buatois, L. A., Mangano, M. G. & Daraei, M. The interplay of environmental constraints and bioturbation on matground development along the marine depositional profile during the Ordovician Radiation. Geobiology 20, 33–270 (2022).Article 
CAS 

Google Scholar 
Bayet-Goll, A., Daraei, M., Geyer, G., Bahrami, N. & Bagheri, F. Environmental constraints on the distribution of matground and mixground ecosystems across the Cambrian Series 2–Miaolingian boundary interval in Iran: A case study for the central sector of northern Gondwana. J. Afr. Earth Sci. 176, 104120 (2021).Article 

Google Scholar 
Minter, N. J. et al. In The Trace-Fossil Record of Major Evolutionary Events: Volume 1: Precambrian and Paleozoic (eds. Mángano M. G. & Buatois L. A.) 157–204 (Springer Netherlands, 2016).Minter, N. J. et al. Early bursts of diversification defined the faunal colonization of land. Nat. Ecol. Evol. 1, 0175 (2017).Article 

Google Scholar 
Nemec, W. & Steel, R. J. In Fan Deltas: sedimentology and tectonic settings (eds. Nemec W. & Steel R. J.) 3–13 (Blackie and Son, 1988).Postma, G. An analysis of the variation in delta architecture. Terra Nova 2, 124–130 (1990).ADS 
Article 

Google Scholar 
Prior, D. B. & Bornhold, B. D. Submarine sedimentation on a developing Holocene fan delta. Sedimentology 36, 1053–1076 (1989).ADS 
Article 

Google Scholar 
Piper, D. J. W., Kontopoulos, N., Anagnostou, C., Chronis, G. & Panagos, A. G. Modern fan deltas in the western Gulf of Corinth, Greece. Geo-Mar. Lett. 10, 5–12 (1990).ADS 
Article 

Google Scholar 
Rasmussen, H. Nearshore and alluvial facies in the Sant Llorenç del Munt depositional system: Recognition and development. Sediment. Geol. 138, 71–98 (2000).ADS 
Article 

Google Scholar 
Steel, R., Rasmussen, H., Eide, S., Neuman, B. & Siggerud, E. Anatomy of high-sediment supply, transgressive tracts in the Vilomara composite sequence, Sant LlorencË del Munt, Ebro Basin, NE Spain. Sediment. Geol. 138, 125–142 (2000).ADS 
Article 

Google Scholar 
Zavala, C. et al. Deltas: A new classification expanding Bates’s concepts. J. Palaeogeogr. 10, 23 (2021).ADS 
Article 

Google Scholar 
Ekdale, A. A. & Lewis, D. W. Trace fossils and paleoenvironmental control of ichnofacies in a late Quaternary gravel and loess fan delta complex, New Zealand. Palaeogeogr. Palaeoclimatol. Palaeoecol. 81, 253–279 (1991).Article 

Google Scholar 
Buatois, L. A. & Mángano, M. G. Ichnology: Organism-substrate Interactions in Space and Time (Cambridge University Press, 2011).Book 

Google Scholar 
Hovikoski, J., Uchman, A., Alsen, P. & Ineson, J. Ichnological and sedimentological characteristics of submarine fan-delta deposits in a half-graben, Lower Cretaceous Palnatokes Bjerg Formation, NE Greenland. Ichnos 26, 28–57 (2019).Article 

Google Scholar 
Sendra, J., Reolid, M. & Reolid, J. Palaeoenvironmental interpretation of the Pliocene fan-delta system of the Vera Basin (SE Spain): Fossil assemblages, ichnology and taphonomy. Palaeoworld 29, 769–788 (2020).Article 

Google Scholar 
Kreis, L. K. et al. Lower Paleozoic map series: Saskatchewan. Miscellaneous Report 2004–8 (CD-ROM) (Saskatchewan Industry and Resources, 2004).
Google Scholar 
Marsh, A. & Love, M. In Saskatchewan Ministry of the Economy, Saskatchewan Geological Survey, Saskatchewan Geological Survey Vol. Open File 2014-1, set of 156 maps (2014).Sawatzky, H. B., Agarwal, R. G. & Wilson, W. Helium prospects of southwest Saskatchewan. 26 (Saskatchewan Department of Mineral Resources, 1960).Fyson, W. K. Deadwood and Winnipeg stratigraphy in southwestern Saskatchewan. Report 64, 37 (Saskatchewan Department of Mineral Resources, 1961).Kent, D. M. Paleotectonic controls on sedimentation in northern Williston Basin area, Saskatchewan. AAPG Bull. 67, 1345–1345 (1983).
Google Scholar 
Kent, D. M. & Haidl, F. M. The distribution of Ashern and Winnipegosis strata (Middle Devonian) on the Swift Current Platform, southern Saskatchewan. Summary of Investigations, Miscellaneous Report 93-4, 201–206 (Saskatchewan Geological Survey, Saskatchewan Energy and Mines, 1993).MacEachern, J. A., Zaitlin, B. A. & Pemberton, S. G. A sharp-based sandstone of the Viking Formation, Joffre Field, Alberta, Canada; criteria for recognition of transgressively incised shoreface complexes. J. Sediment. Res. 69, 876–892 (1999).ADS 
Article 

Google Scholar 
Pemberton, S. G., Frey, R. W. & Bromley, R. G. The ichnotaxonomy of Conostichus and other plug-shaped ichnofossils. Can. J. Earth Sci. 25, 866–892 (1988).ADS 
Article 

Google Scholar 
Hall, J. & Whitfield, R. P. Paleontology. In US Geol. Expl. 40th Par. Rept. 4, 197–302 (1877).Walcott, C. D. Cambrian and Lower Ozarkian trilobites. Smithson. Misc. Coll. 75, 53–60 (1924).
Google Scholar 
Meek, F. B. Preliminary paleontological report, consisting of lists and descriptions of fossils, with remarks on the ages of the rocks in which they were found. In U. S. Geol. Surv. Terr. 6th Ann. Rept., 429–518 (1873).Walcott, C. D. Cambrian geology and paleontology of Cambrian trilobites. Smithson. Misc. Coll. 64, 157–258 (1916).
Google Scholar 
Harding, S. C. & Ekdale, A. A. Trace fossils and glauconitic pellets provide insight into Cambrian siliciclastic marine environments. Palaios 33, 256–265 (2018).ADS 
Article 

Google Scholar 
Shillito, A. P. & Davies, N. S. The Tumblagooda Sandstone revisited: Exceptionally abundant trace fossils and geological outcrop provide a window onto Palaeozoic littoral habitats before invertebrate terrestrialization. Geol. Mag. 157, 1939–1970 (2020).ADS 
Article 

Google Scholar 
Shillito, A. P. & Davies, N. S. Archetypally Siluro-Devonian ichnofauna in the Cowie Formation, Scotland: Implications for the myriapod fossil record and Highland Boundary Fault Movement. Proc. Geol. Assoc. 128, 815–828 (2017).Article 

Google Scholar 
Buatois, L. A. et al. The invasion of the land in deep time: integrating Paleozoic records of paleobiology, ichnology, sedimentology, and geomorphology. Integr. Comp. Biol. 0, 1–35. https://doi.org/10.1093/icb/icac059 (2022).Davies, N. S. & Gibling, M. R. Paleozoic vegetation and the Siluro-Devonian rise of fluvial lateral accretion sets. Geology 38, 51–54 (2010).ADS 
Article 

Google Scholar 
McMahon, W. J. & Davies, N. S. Evolution of alluvial mudrock forced by early land plants. Science 359, 1022–1024 (2018).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Fedo, C. M. & Cooper, J. D. Braided fluvial to marine transition; the basal Lower Cambrian Wood Canyon Formation, southern Marble Mountains, Mojave Desert, California. J. Sediment. Res. 60, 220–234 (1990).
Google Scholar 
Eyre, B. Early Cambrian alluvial fan-deltas in the Georgina Basin, Australia. Aust. J. Earth Sci. 41, 27–36 (1994).ADS 
Article 

Google Scholar 
MacNaughton, R. B., Dalrymple, R. & Narbonne, G. M. Early Cambrian braid-delta deposits, MacKenzie Mountains, north-western Canada. Sedimentology 44, 587–609 (1997).ADS 
Article 

Google Scholar 
Muhlbauer, J. G. & Fedo, C. M. Architecture of a river-dominated, wave- and tide-influenced, pre-vegetation braid delta: Cambrian middle member of the Wood Canyon Formation, southern Marble Mountains, California, U.S.A. J. Sediment. Res. 90, 1011–1036 (2020).ADS 
Article 

Google Scholar 
McMahon, W. J., Davies, N. S. & Went, D. J. Negligible microbial matground influence on pre-vegetation river functioning: Evidence from the Ediacaran-Lower Cambrian Series Rouge, France. Precambrian Res. 292, 13–34 (2017).ADS 
CAS 
Article 

Google Scholar 
Mikuláš, R. Trace fossils from the Paseky Shale (Early Cambrian, Czech Republic). J. Czech Geol. Soc. 40, 37–44 (1995).
Google Scholar 
MacNaughton, R. B. & Narbonne, G. M. Evolution and ecology of Neoproterozoic-Lower Cambrian trace fossils, NW Canada. Palaios 14, 97–115 (1999).ADS 
Article 

Google Scholar 
Buatois, L. A. et al. Colonization of brackish-water systems through time: Evidence from the trace-fossil record. Palaios 20, 321–347 (2005).ADS 
Article 

Google Scholar 
Hofmann, R., Mángano, M. G., Elicki, O. & Shinaq, R. Paleoecologic and biostratigraphic significance of trace fossils from shallow- to marginal-marine environments from the Middle Cambrian (Stage 5) of Jordan. J. Paleontol. 86, 931–955 (2012).Article 

Google Scholar 
Mángano, M. G., Buatois, L. A., Hofmann, R., Elicki, O. & Shinaq, R. Exploring the aftermath of the Cambrian explosion: The evolutionary significance of marginal- to shallow-marine ichnofaunas of Jordan. Palaeogeogr. Palaeoclimatol. Palaeoecol. 374, 1–15 (2013).Article 

Google Scholar 
Mángano, M. G. et al. Were all trilobites fully marine? Trilobite expansion into brackish water during the early Palaeozoic. Proc. R. Soc. Lond. B Biol. Sci. 288, 20202263 (2021).
Google Scholar 
Siggerud, E. I. H. & Steel, R. J. Architecture and trace-fossil characteristics of a 10000–20000 Year, fluvial-to-marine sequence, SE Ebro Basin, Spain. J. Sediment. Res. 69, 365–383 (1999).ADS 
Article 

Google Scholar 
Lockley, M. G., Rindsberg, A. K. & Zeiler, R. M. The paleoenvironmental significance of the nearshore Curvolithus ichnofacies. Palaios 2, 255–262 (1987).ADS 
Article 

Google Scholar 
Folk, R. L. Petrology of Sedimentary Rocks (Hemphill Publishing Company, 1980).
Google Scholar 
Wentworth, C. K. A scale grade and class terms for clastic sediments. J. Geol. 30, 377–392 (1922).ADS 
Article 

Google Scholar 
Pettijohn, F. J., Potter, P. E. & Siever, R. Sand and Sandstone 2nd edn. (Springer, 1987).Book 

Google Scholar 
Ingram, R. L. Terminology for the thickness of stratification and parting units in sedimentary rocks. Geol. Soc. Am. Bull. 65, 937–938 (1954).ADS 
Article 

Google Scholar 
Reineck, H.-E. Sedimentgefüge im Bereich der südliche Nordsee. Abh. Senckb. Naturforsch. Ges. 505, 1–138 (1963).
Google Scholar  More

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.This is a summary of: Needham, D. M. et al. The microbiome of a bacterivorous marine choanoflagellate contains a resource-demanding obligate bacterial associate. Nat. Microbiol. https://doi.org/10.1038/s41564-022-01174-0 (2022). More

Ballantyne, G., Baldock, K. C. R., Rendell, L. & Willmer, P. G. Pollinator importance networks illustrate the crucial value of bees in a highly speciose plant community. Sci. Rep. 7, 8389 (2017).ADS 
PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Smith, M. R., Singh, G. M., Mozaffarian, D. & Myers, S. S. Effect of decreases of animal pollinators on human nutrition and global health: A modeling analysis. Lancet 386, 1964–1972 (2015).PubMed 
Article 

Google Scholar 
Jung, C. & Cho, S. Relationship between honey bee population and honey production in Korea: A historical trend analysis. J. Apic. 30(1), 7–12 (2015).
Google Scholar 
Abrol, D. P. Asiatic Honey Bee Apis cerana: Biodiversity Conservation and Agricultural Production (Springer, 2013).Book 

Google Scholar 
Chandel, Y. S., Kumar, A. & Srivastva, S. Comparative performance of Apis mellifera L. vis a vis Apis cerana Fab. on toria (Brassica campestris var Toria) in mid-hill zone of Himachal Pradesh, India. Indian J. Agric. Res. 34, 264–267 (2000).
Google Scholar 
Feng, M., Ramadan, H., Han, B., Yu, F. & Li, J. Hemolymph proteome changes during worker brood development match the biological divergences between western honey bees (Apis mellifera) and eastern honey bees (Apis cerana). BMC Genomics 15, 563–576 (2014).PubMed 
PubMed Central 
Article 

Google Scholar 
Peng, Y. S., Fang, Y., Xu, S. & Ge, L. The resistance mechanism of the Asian honey bee, Apis cerana Fabr., to an ectoparasitic mite, Varroa jacobsoni Oudemans. J. Invertebr. Pathol. 49, 54–60 (1987).Article 

Google Scholar 
McClenaghan, et al. Behavioral responses of honey bees, Apis cerana and Apis mellifera, to Vespa mandarinia marking and alarm pheromones. J. Apic. Res. 58(1), 141–148 (2018).Article 

Google Scholar 
Lin, Z. et al. Go east for better honey bee health: Apis cerana is faster at hygienic behavior than A. mellifera. PLoS ONE 11(9), e0162647 (2016).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Verma, L. R. & Dutta, P. C. Foraging behaviour of Apis cerana indicia and Apis mellifera in pollinating apple flowers. J. Apic. Res. 25, 197–201 (1986).Article 

Google Scholar 
Wang, Z. & Tan, K. Comparative analysis of olfactory learning of Apis cerana and Apis mellifera. Apidologie 45(1), 45–52 (2014).ADS 
Article 

Google Scholar 
Beekman, M. & Ratnieks, F. L. W. Long-range foraging by the honey-bee Apis mellifera L.. Funct. Ecol. 14, 490–496 (2000).Article 

Google Scholar 
Dyer, F. C. & Seeley, T. D. Dance dialects and foraging range in three Asian honey bee species. Behav. Ecol. Sociobiol. 28, 227–233 (1991).Article 

Google Scholar 
Koetz, A. H. Ecology, behaviour and control of Apis cerana with a focus on relevance to the Australian incursion. Insects 4(4), 558–592 (2013).PubMed 
PubMed Central 
Article 

Google Scholar 
Joshi, N. C. & Joshi, P. C. Foraging behavior of Apis spp. On apple flowers in a subtropical environment. N. Y. Sci. J. 3(3), 71–76 (2010).
Google Scholar 
Atwal, A. S. & Sharma, O. P. The dominance of Apis mellifera over Apis indica. Am. Bee J. 111, 343 (1971).
Google Scholar 
Kato, M., Shibata, A., Yasui, T. & Nagamasu, H. Impact of introduced honey bees, Apis mellifera, upon native bee communities in the Bonin (Ogasawara) Islands. Res. Popul. Ecol. 41, 217–228 (1999).Article 

Google Scholar 
Thorp, D. W., Wenner, A. M. & Barthell, J. F. Pollen and nectar resource overlap among bees on Santa Cruz Island. MBC Appl. Environ. Sci. 2020, 261–267 (2000).
Google Scholar 
Yang, G. Harm of introducing the western honey bee Apis mellifera L. to the Chinese honey bee Apis cerana F. and its ecological impact. Acta Entomol. Sin. 48, 401–406 (2005) ((in Chinese)).
Google Scholar 
Dubois, T., Pasquaretta, C., Barron, A. B., Gautrais, J. & Lihoreau, M. A model of resource partitioning between foraging bees based on learning. PLoS Comput. Biol. 17(7), e1009260 (2021).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Louveaux, J., Maurizio, A. & Vorwohl, G. Methods of melissopalynology. Bee World 59, 139–153 (1978).Article 

Google Scholar 
Hawkins, J., de Vere, N., Griffith, A. & Ford, C. R. Using DNA metabarcoding to Identify the floral composition of honey: A new tool for investigating honey bee foraging preferences. PLoS ONE 10(8), e0134735 (2015).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Smart, M. D. et al. A Comparison of honey bee-collected pollen from working agricultural lands using light microscopy and its metabarcoding. Environ. Entomol. 46(1), 38–49 (2016).
Google Scholar 
Keller, A. et al. Evaluating multiplexed next-generation sequencing as a method in palynology for mixed pollen samples. Plant Biol. 17, 558–566 (2015).CAS 
PubMed 
Article 

Google Scholar 
Richardson, R. T. et al. Rank-based characterization of pollen assemblages collected by honey bees using a multi-locus metabarcoding approach. Appl. Plant Sci. 3, 1500043 (2015).Article 

Google Scholar 
Sickel, W. et al. Increased efficiency in identifying mixed pollen samples by meta-barcoding with a dual-indexing approach. BMC Ecol. 15, 20 (2015).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Kress, W. J. & Erickson, D. L. DNA barcodes: Genes, genomics, and bioinformatics. PNAS 105, 2761–2762 (2008).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Hollingsworth, P. M., Graham, S. W. & Little, D. P. Choosing and using a plant DNA barcode. PLoS ONE 6, e19254 (2011).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Chen, S. et al. Validation of the ITS2 region as a novel DNA barcode for identifying medicinal plant species. PLoS ONE 5, e8613 (2010).ADS 
PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
CBOL Plant Working Group. A DNA barcode for land plants. PNAS 106, 12794–12797 (2009).ADS 
PubMed Central 
Article 

Google Scholar 
Pornon, A. et al. Using metabarcoding to reveal and quantify plant-pollinator interactions. Sci. Rep. 6, 27282 (2016).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Bell, K. L. et al. Quantitative and qualitative assessment of pollen DNA metabarcoding using constructed species mixtures. Mol. Ecol. 28(2), 431–455 (2018).PubMed 
Article 
CAS 

Google Scholar 
Baksay, S. et al. Experimental quantificarion of pollen with DNA metabarcoding using ITS1 and trnL. Sci. Rep. 10, 4202 (2020).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Ghosh, S. & Jung, C. Nutritional value of bee-collected pollens of hardy kiwi, Actinidia arguta (Actinidiaceae) and oak, Quercus sp. (Fagaceae). J. Asia Pac. Entomol. 20(1), 245–251 (2017).Article 

Google Scholar 
Brunet, J., Thairu, M. W., Henss, J. M., Link, R. I. & Kluevert, J. A. The effects of flower, floral display, and reward sizes on bumblebee foraging behaviour when pollen is the reward and plants are dichogamous. Int. J. Plant Sci. 176(9), 811–819 (2015).Article 

Google Scholar 
de Vere, N. et al. Using DNA metabarcoding to investigate honey bee foraging reveals limited flower use despite high floral availability. Sci. Rep. 7(1), 42838 (2017).ADS 
PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Barth, O. M. In O polen no mel brasileiro. Rio de Janeiro, Luxor 151 (1989).Lo, T., Gloag, R. S., Anderson, D. L. & Oldroyd, B. P. A molecular phylogeny of the genus Apis suggests that the giant honey bee of southern India, A. indica Fabricius, are valid species. Syst. Entomol. 35, 226–223 (2010).Article 

Google Scholar 
Pirk, C. W. W., Sole, C. L. & Crewe, R. M. Pheromones. In Honey Bees of Asia (eds Hepburn, H. R. & Radloff, S. E.) 207–214 (Springer, Berlin, 2011).Chapter 

Google Scholar 
Theisen-Jones, H. & Bienefeld, K. The Asian honey bee (Apis cerana) is significantly in decline. Bee World 93, 90–97 (2016).Article 

Google Scholar 
Sakagami, S. F. Some interspecific relations between Japanese and European honey bees. J. Anim. Ecol. 28, 51–68 (1959).Article 

Google Scholar 
Thomson, D. Competitive interactions between the invasive European honey bee and native bumble bees. Ecology 85, 458–470 (2004).Article 

Google Scholar 
Brosi, B. J. & Briggs, H. M. Single pollinator species losses reduce floral fidelity and plant reproductive function. PNAS 110(32), 13044–13048 (2013).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Iwasaki, J. M. et al. Floral usage partitioning and competition between social (Apis mellifera, Bombus terrestris) and solitary bees in New Zealand: Niche partitioning via floral preferences?. Austral Ecol. 43(8), 937–948 (2018).Article 

Google Scholar 
Rodrigues, C. S., Ferasso, D. C., Mossi, A. J. & Coelho, G. C. Pollen resources partitioning of stingless bees (Hymenoptera: Apidae) from the southern Atlantic forest Acta Scientiarum. Biol. Sci. 42, e48714 (2020).
Google Scholar 
Lucas, A. et al. Floral resource partitioning by individuals within generalised hoverfly pollination networks revealed by DNA metabarcoding. Sci. Rep. 8(1), 5133 (2018).ADS 
PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Lucek, K. et al. Metabarcoding of honey to assess differences in plant-pollinator interactions between urban and non-urban sites. Apidologie 50, 317–329 (2019).Article 

Google Scholar 
Tatsuno, M. & Osawa, N. Flower visitation patterns of the coexisting honey bees Apis cerana japonica and Apis mellifera (Hymentoptera: Apidae). Entomol. Sci. https://doi.org/10.1111/ens.12206 (2016).Article 

Google Scholar 
Kuang, B. Y. & Kuang, H. O. Biology of the Honey bee (Yunnan Science and Technology Press, 2002) ([In Chines]).
Google Scholar 
Ghorab, A. et al. Sensorial, melissopalynological and physico-chemical characteristics of honey from Babors Kabylia’s region (Algeria). Foods 10, 225 (2021).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Sasaki, M. In Bee’s eye view of flowering plants: Nectar- and pollen-source plants and related honey bee products. Kaiyusha, Tokyo, Japan (2010).Simpson, M. G. Diversity and classification of flowering plant: Eudicots. Plant Syst. 2010, 275–448 (2010).Article 

Google Scholar 
Wilms, W. & Wiechers, B. Floral resource partitioning between native Melipona bees and the introduced Africanized honey bee in the Brazilian Atlantic rain forest. Apidologie 28, 339–355 (1997).Article 

Google Scholar 
Klein, S. et al. Honey bees increase their foraging performance and frequency of pollen trips through experience. Sci. Rep. 9, 6778 (2019).ADS 
PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Bovo, S. et al. Shotgun metagenomics of honey DNA: Evaluation of a methodological approach to describe a multikingdom honey bee derived environmental DNA signature. PLoS ONE 13(10), e0205575 (2018).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Prosser, S. W. J. & Hebert, P. D. N. Rapid identification of the botanical and entomological sources of honey using DNA metabarcoding. Food Chem. 214, 183–191 (2017).CAS 
PubMed 
Article 

Google Scholar 
Dhaliwai, H. S. & Sharma, P. L. Foraging range of the Indian honey bee. J. Apic. Res. 13, 137–141 (1974).Article 

Google Scholar 
ESRI. ArcGis Pro (Version 10.6); ESRI Inc.: Redlands, CA, USA (2020).Palmieri, L., Bozza, E. & Giongo, L. Soft fruit traceability in food matrices using real-time PCR. Nutrients 1, 316–328 (2009).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Kress, W. J. et al. Plant DNA barcodes and a community phylogeny of a tropical forest dynamics plot in Panama. PNAS 106, 18621–18626 (2009).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Kress, J. W. & Erickson, L. D. A two-locus global DNA barcode for land plants: The coding rbcL gene complements the non-coding trnH-psbA spacer region. PLoS ONE 2(6), 1–10 (2007).
Google Scholar 
Bolger, A. M., Lohse, M. & Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 30(15), 2114–2120 (2014).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Rognes, T., Flouri, T., Nichols, B., Quince, C. & Mahé, F. VSEARCH: A versatile open source tool for metagenomics. PeerJ 4, e2584 (2016).PubMed 
PubMed Central 
Article 

Google Scholar 
Joshi, N. A. & Fass, J. N. Sickle: A sliding-window, adaptive, quality-based trimming tool for FastQ files (Version 1.33). https://github.com/najoshi/sickle (2011).Boyer, F. et al. OBITools: A Unix-inspired software package for DNA metabarcoding. Mol. Ecol. Res. 16, 176–182 (2016).CAS 
Article 

Google Scholar 
Altschul, S. F. et al. Gapped BLAST and PSI-BLAST: A new generation of protein database search programs. Nucleic Acids Res. 25, 3389–3402 (1997).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
R Core Team. R: A language and environment for statistical computing. https://www.gbif.org/tool/81287/r-a-language-andenvironment-for-statistical-computing (2022).Bosch, J., Retana, J. & Cerdá, X. Flowering phenology, floral traits and pollinator composition in a herbaceous Mediterranean plant community. Oecologia 109, 583–591 (1997).ADS 
PubMed 
Article 

Google Scholar 
Endress, P. K. Flower structure and trends of evolution in eudicots and their major subclades. Ann. Missouri Bot. Gard. 97(4), 541–583 (2010).Article 

Google Scholar 
Gómez, J. M., Torices, R., Lorite, J., Klingenberg, C. P. & Perfectti, F. The role of pollinators in the evolution of corolla shape variation, disparity and integration in a highly diversified plant family with a conserved floral bauplan. Ann. Bot. 117, 899–904 (2016).Article 

Google Scholar 
Watts, S., Dormann, C. F., González, M. M. & Ollerton, J. The influence of floral traits on specialization and modularity of plant-pollinator networks in a biodiversity hotspot in the Peruvian Andes. Ann. Bot. 118, 415–429 (2016).PubMed 
PubMed Central 
Article 

Google Scholar 
Lee, C. B. Coloured Flora of Korea. Hyangmunsa, Seoul, Republic of Korea (2003). More

DataWe used temperature data, rainfall data, and data on the incidence of malaria collected from 1901 to 2015 for 43 African countries to construct networks to determine the relationships between transmission of malaria and climate change elements, especially temperature and rainfall. Data resolution is given by the latitude and longitude of the capital city for every country in Africa. Temperature and rainfall data are provided in terms of monthly averages in the country wise. The nodes in the network represent the country, and the edges in the network represent the relationship between countries. We collected malaria data from Harvard Dataverse35 and the world malaria report from the WHO31. Data for temperature and rainfall were obtained from the Climate Change Knowledge Portal of the World Bank Group36.Network generation and analysisThe networks were constructed by using the threshold method where the network depends on the mean, standard deviation, and the real number ((n)) used to control the features of the network. Therefore, data for temperature, rainfall, and the incidence of malaria were divided into six groups mostly comprising ranges of 20 years (1900–1920, 1921–1940, 1941–1960, 1961–1980, 1981–2000) as well as the period from 2001 to 2015. The missing data in Malaria incidence data are filled by the average amount of malaria incidence collected per year.In Table S1, a malaria report from the World Health Organization shows that the rate of death is directly proportional to the incidence of malaria35. The death toll in Africa from malaria is about 98% of world deaths from malaria. Such deaths in African regions decrease thanks to efforts the WHO, governments, and the private sector have been conducting to prevent them. Weather and climate are among the factors that drive increases in malaria infections in different areas.We consider networks based on the threshold method (see the “Methods and Materials” section below). First, we fill the missing malaria incidence data, and we calculate normalized Pearson correlation coefficients of three-time series between two countries. Then, we obtain a correlation matrix for the countries. We estimate the average value of the correlation coefficients from the time intervals 1901–1920, 1921–1940, 1941–1960, 1961–1980, 1981–2000, and 2001–2015 for three time series: temperature, rainfall, and incidence of malaria. We summarize the averages and standard deviations of the correlation coefficients, as shown in Table S2. The mean values from the correlation in temperature are high, compared to those for rainfall and the incidence of malaria. The standard deviations in temperature and rainfall are large, but the standard deviation for the incidence of malaria is small.We chose an ad hoc threshold value of the correlation coefficients to generate sparse networks. The characteristic values for (n) of the threshold are given in Table S3. We consider three types of thresholds in order to observe changes in the networks according to the threshold.Let us define the normalized variance of each time series. We considered time series (T_{i} left( t right)), (M_{i} left( t right)), and (R_{i} left( t right)) in country (i) for temperature, the incidence of malaria, and rainfall, respectively. We define normalized variance as$$r_{ij} = frac{{x_{i} left( t right)x_{j} left( t right) – x_{i} left( t right)x_{j} left( t right)}}{{sigma_{i} sigma_{j} }}$$
(1)
where (x_{i} left( t right)) = (T_{i} left( t right)), (M_{i} left( t right)), (R_{i} left( t right)). We obtained a Pearson correlation matrix for each time series as follows:$$R_{S} = left[ {begin{array}{*{20}c} {r_{11} } & cdots & {r_{1N} } \ vdots & {r_{ij} } & vdots \ {r_{N1} } & cdots & {r_{NN} } \ end{array} } right]$$
(2)
where (S = T, M, R).We calculated the average value, (overline{r }), and the standard deviation, (sigma), for the correlation coefficients of the matrix. We applied the threshold method to generate a sparse network from the correlation matrix. Two countries are connected in the correlation network if and only if the value of the correlation coefficient is greater than, or equal to, the threshold value:$$r_{{ij}} = left{ {begin{array}{*{20}c} 1 & {{text{if}};r_{{ij}} ge bar{r}{text{ + n}}sigma } \ 0 & {{text{otherwise}}} \ end{array} } right.$$
(3)
where (r_{ij}) is the correlation coefficient between two countries, and (n) is an element of real numbers ((n in {mathbb{R}})). The value of (n) determines whether the network is sparsely or densely connected.We use Python programming language, packages, numpy for mathematical functions and random number generator, pandas for data analysis and manipulations, networkx for creation, manipulation, and studying the structure of the complex network, matplotlib for visualization and plotting graph and base map for map projection and visualization of geographic information. More

MacFarlane, G. R., Koller, C. E. & Blomberg, S. P. Accumulation and partitioning of heavy metals in mangroves: A synthesis of field-based studies. Chemosphere 69, 1454–1464 (2007).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Krauss, K. W. & Osland, M. J. Tropical cyclones and the organization of mangrove forests: A review. Ann. Bot. 125, 213–234 (2020).PubMed 

Google Scholar 
Kirk, G. J. D. & Krinzucker, H. J. The potential for nitrification and nitrate uptake in the rhizosphere of wetland plants: A modeling study. Ann. Bot. 96, 639–646 (2005).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Wu, C., Huang, L., Xue, S. G. & Pan, W. S. Oxic and anoxic conditions affect arsenic (As) accumulation and arsenite transporter expression in rice. Chemosphere 168, 969–975 (2017).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Tripathi, R. D. et al. Roles for root iron plaque in sequestration and uptake of heavy metals and metalloids in aquatic and wetland plants. Metallomics 6, 1789–1800 (2014).CAS 
PubMed 
Article 

Google Scholar 
Xiao, W. et al. Continuous flooding stimulates root iron plaque formation and reduces chromium accumulation in rice (Oryza sativa L.). Sci. Total Environ. 788, 147786 (2021).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Lee, C. H., Hsieh, Y. C., Lin, T. H. & Lee, D. Y. Iron plaque formation and its effect on arsenic uptake by different genotypes of paddy rice. Plant Soil 363, 231–241 (2013).CAS 
Article 

Google Scholar 
Dai, M. Y. et al. Phosphorus effects on radial oxygen loss, root porosity and iron plaque in two mangrove seedlings under cadmium stress. Mar. Pollut. Bull. 119, 262–269 (2017).CAS 
PubMed 
Article 

Google Scholar 
Liu, C. Y., Chen, C. L., Gong, X. F., Zhou, W. B. & Yang, J. Y. Progress in research of iron plaque on root surface of wetland plants. Acta Ecol. Sin. 34, 2470–2480 (2014).CAS 

Google Scholar 
Li, J., Liu, J. C., Yan, C. L., Du, D. L. & Li, H. L. The alleviation effect of iron on cadmium phytotoxicity on mangrove A. marina. Alleviation effect of rion on cadmium phytotoxicity in mangrove Avicennia marina (Forsk.) Vierh. Chemosphere 226, 413–420 (2019).ADS 
CAS 
Article 

Google Scholar 
Zhang, J. Y. et al. Effects of nano-Fe3O4-modified biochar on iron plaque formation and Cd accumulation in rice (Oryza sativa L.). Environ. Pollut. 260, 113970 (2020).CAS 
PubMed 
Article 

Google Scholar 
Farhat, Y. A., Kim, S. H., Seyfferth, A. L., Zhang, L. & Neumann, R. B. Altered arsenic availability, uptake, and allocation in rice under elevated temperature. Sci. Total Environ. 763, 143049 (2021).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Wu, X. Y. et al. A review of toxicity and mechanisms of individual and mixtures of heavy metals in the environment. Environ. Sci. Pollut. R. 23, 8244–8259 (2016).CAS 
Article 

Google Scholar 
Abubakar, U. S., Zulkifli, S. Z. & Ismail, A. Heavy metals bioavailability and pollution indices evaluation in the mangrove surface sediment of Sungai Puloh Malaysia. Environ. Earth. Sci. 77, 225 (2018).CAS 
Article 

Google Scholar 
Kulkarni, R., Deobagkar, D. & Zinjarde, S. Metals in mangrove ecosystems and associated biota: A global perspective. Ecotox. Environ. Safe. 153, 215–228 (2018).CAS 
Article 

Google Scholar 
Shi, C., Ding, H., Zan, Q. J. & Li, R. L. Spatial variation and ecological risk assessment of heavy metals in mangrove sediments across China. Mar. Pollut. Bull. 143, 115–124 (2019).CAS 
PubMed 
Article 

Google Scholar 
Cheng, H. et al. Mixture of Pb, Zn and Cu on root permeability and radial oxygen loss in the mangrove Bruguiera gymnorrhiza. Ecotoxicology 29, 691–697 (2020).CAS 
PubMed 
Article 

Google Scholar 
Cheng, S. S. et al. Temporal variations in physiological responses of kandelia obovata seedlings exposed to multiple heavy metals. Mar. Pollut. Bull. 124, 1089–1095 (2017).CAS 
PubMed 
Article 

Google Scholar 
Shen, X. X. et al. Does combined heavy metal stress enhance iron plaque formation and heavy metal bioaccumulation in Kandelia obovata?. Environ. Exp. Bot. 186, 104463 (2021).CAS 
Article 

Google Scholar 
Shen, X. X. et al. Interactive effects of single, binary and trinary trace metals (lead, zinc and copper) on the physiological responses of Kandelia obovata seedlings. Environ. Geochem. Hlth. 41, 135–148 (2019).CAS 
Article 

Google Scholar 
Youssef, T. & Saenger, P. Anatomical adaptive strategies to flooding and rhizosphere oxidation in mangrove seedlings. Aust. J. Bot. 44, 297–313 (1996).Article 

Google Scholar 
Cheng, H., Wang, Y. S., Fei, J., Jiang, Z. Y. & Ye, Z. H. Differences in root aeration, iron plaque formation and waterlogging tolerance in six mangroves along a continues tidal gradient. Ecotoxicology 24, 1659–1667 (2015).CAS 
PubMed 
Article 

Google Scholar 
Takahashi, H., Yamauchi, T., Colmer, T. D. & Nakazono, M. Aerenchyma formation in plants. In Low-Oxygen Stress in Plants (eds van Dongen, J. T. & Licausi, F.) 247–265 (Springer, 2014).Chapter 

Google Scholar 
Yamauchi, T., Colmer, T. D., Pedersen, O. & Nakazono, M. Regulation of root traits for internal aeration and tolerance to soil waterlogging-flooding stress. Plant Physiol. 176, 1118–1130 (2018).CAS 
PubMed 
Article 

Google Scholar 
Cheng, H. et al. The role of radial oxygen loss and root anatomy on zinc uptake and tolerance in mangrove seedlings. Environ. Pollut. 158, 1189–1196 (2010).CAS 
PubMed 
Article 

Google Scholar 
Liu, Y. et al. Mixed heavy metals tolerance and radial oxygen loss in mangrove seedlings. Mar. Pollut. Bull. 58, 1843–1849 (2009).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Wu, C., Li, H., Ye, Z., Wu, F. & Wong, M. H. Effects of As levels on radial oxygen loss and as speciation in rice. Environ. Sci. Pollut R. 20, 8334–8341 (2013).CAS 
Article 

Google Scholar 
Mendellshn, A., Kleiss, B. A. & Wakeley, J. S. Factors controlling the formation of oxidized root channels—a review. Wetlands 15, 37–46 (1995).Article 

Google Scholar 
Moller, C. L. & Sand-Jesen, K. Iron plaques improve the oxygen supply to root meristems of the freshwater plant Lobelia dortmanna. New Phytol. 179, 848–856 (2008).PubMed 
Article 
CAS 

Google Scholar 
Yang, J. X., Liu, Y. & Ye, Z. H. Root-induced changes of pH, eh, Fe (II) and fractions of Pb and Zn in rhizosphere soils of four wetland plants with different radial oxygen losses. Pedosphere 22, 518–527 (2012).CAS 
Article 

Google Scholar 
Hu, M., Li, F., Liu, C. & Wu, W. The diversity and abundance of as (III) oxidizers on root iron plaque is critical for arsenic bioavailability to rice. Sci. Rep. 5, 13611 (2015).ADS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Huang, G. X., Ding, C. F., Li, Y. S., Zhang, T. L. & Wang, X. X. Selenium enhances iron plaque formation by elevating the radial oxygen loss of roots to reduce cadmium accumulation in rice (Oryza sativa L.). J. Hazard. Mater. 398, 122860 (2020).CAS 
PubMed 
Article 

Google Scholar 
Thakur, S. et al. Plant-driven removal of heavy metals from soil: Uptake, translocation, tolerance mechanism, challenges, and future perspectives. Environ. Monit. Assess. 188, 206–212 (2016).PubMed 
Article 

Google Scholar 
Huang, H., Zhu, Y., Chen, Z., Yin, X. & Sun, G. Arsenic mobilization and speciation during iron plaque decomposition in a paddy soil. J. Soil. Sediment. 12, 402–410 (2012).CAS 
Article 

Google Scholar 
Zhong, S. Q. Effect of iron plaque on root growth and activity of two wetland plants. J. Hydroecol. 36, 74–79 (2015).
Google Scholar 
Khan, N. et al. Root iron plaque on wetland plants as a dynamic pool of nutrients and contaminants. Adv. Agron. 138, 1–96 (2016).Article 

Google Scholar 
Ma, H. H. et al. Formation of iron plaque on roots of Iris pseudacorus and its consequence for cadmium immobilization is impacted by zinc concentration. Ecotox. Environ. Safe. 193, 110306 (2020).CAS 
Article 

Google Scholar 
Martinez, S., Sáenz, M. E., Alberdi, J. L. & Di Marzio, W. D. Comparative ecotoxicity of single and binary mixtures exposures of cadmium and zinc on growth and biomarkers of Lemna gibba. Ecotoxicology 29, 571–583 (2020).CAS 
PubMed 
Article 

Google Scholar 
Huang, Y. Z., Hu, Y. & Liu, Y. X. Heavy metal accumulation in iron plaque and growth of rice plants upon exposure to single and combined contamination by copper, cadmium and lead. Acta Ecol. Sin. 29, 320–326 (2009).Article 

Google Scholar 
Deraison, H., Badenhausser, I., Börger, L. & Gross, N. Herbivore effect traits and their impact on plant community biomass: An experimental test using grasshoppers. Funct. Ecol. 29, 650–661 (2015).Article 

Google Scholar 
Wang, F., Wang, X. & Song, N. Polythylene microplastics increase cadmium uptake in lettuce (Lactuca sativa L.) by altering the soil microenvironment. Sci. Total Environ. 784, 147133 (2021).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Yu, H. et al. Microplastic residues in wetland ecosystems: Do they truly threaten the plant-microbe-soil system?. Environ. Int. 156, 106708 (2021).CAS 
PubMed 
Article 

Google Scholar 
He, B., Li, R. L., Chai, M. W. & Qiu, G. Y. Threat of heavy metal contamination in eight mangrove plants from the Futian mangrove forest, China. Environ. Geochem. Hlth. 36, 467–476 (2014).CAS 
Article 

Google Scholar 
Du, J. N., Yan, C. L. & Li, Z. D. Formation of iron plaque on mangrove Kandalar. Obovata (S.L.) Root surfaces and its role in cadmium uptake and translocation. Mar. Pollut. Bull. 74, 105–109 (2013).CAS 
PubMed 
Article 

Google Scholar 
Hao, Z. B., Cang, J. & Xu, Z. Plant Physiology Experiment (Harbin Institute of Technology Press, 2004).
Google Scholar 
Kludze, H. K., Delaune, R. D. & Patrick, W. H. A colorimetric method for assaying dissolved oxygen loss from container-grown rice roots. Agron. J. 86, 483–487 (1994).CAS 
Article 

Google Scholar 
Kludze, H. K., Delaune, R. D. & Patrick, W. H. Aerenchyma formation and methane and oxygen-exchange in rice. Soil Sci. Soc. Am. J. 57, 386–391 (1993).ADS 
CAS 
Article 

Google Scholar 
Mei, X. Q., Yang, Y., Tam, N. F. Y., Wang, Y. W. & Li, L. Roles of root porosity, radial oxygen loss, Fe plaque formation on nutrient removal and tolerance of wetland plants to domestic wastewater. Water Res. 50, 147–159 (2014).CAS 
PubMed 
Article 

Google Scholar 
Taylor, G. J. & Crowder, A. Use of the DCB technique for extraction of hydrous iron oxides from roots of wetland plants. Am. J. Bot. 70, 1254–1257 (1983).CAS 
Article 

Google Scholar 
USEPA (United States Environmental Protection Agency). Method 3052: microwave assisted acid digestion of siliceous and organically based matrices SW-846. DC: Washington (1996). More

Bromley, R. G. & Gale, A. S. The lithostratigraphy of the English Chalk Rock. Cretac. Res. 3, 273–306 (1982).Article 

Google Scholar 
Scholle, P. A., Arthur, M. A. & Ekdale, A. A. Pelagic environment. In Carbonate Depositional Environments (eds Scholle, P. A. et al.) 619–691 (Am. Ass. Petrol. Geol. Mem. 33, 1983).Chapter 

Google Scholar 
Gealy, E. L., Winterer, E. L. & Moberly, R. Methods, conventions, and general observations. Initial Rep. Deep Sea Drill. Proj. 7, 9–26 (1971).
Google Scholar 
Kroenke, L. W. et al. Ocean Drilling Program. Proc. ODP, Init. Repts. 130, College Station (1991).Dunham, R. L. Classification of carbonate rocks according to depositional texture. Mem. Am. Assoc. Petrol. Geol. 1, 108–121 (1962).
Google Scholar 
Quine, M. & Bosence, D. Stratal geometries, facies and sea-floor erosion in Upper Cretaceous chalk, Normandy, France. Sedimentology 38, 1113–1152 (1991).ADS 
Article 

Google Scholar 
Røgen, B., Gommesen, L. & Fabricius, I. L. Grain size distributions of Chalk from Image analysis of electron micrographs. Comput. Geosci. 27, 1071–1080 (2001).ADS 
Article 

Google Scholar 
Saïag, J. et al. Classifying chalk microtextures: Sedimentary versus diagenetic origin (Cenomanian–Santonian, Paris Basin, France). Sedimentology 66, 2976–3007 (2019).Article 
CAS 

Google Scholar 
Scholle, P. A. Chalk diagenesis and its relation to petroleum exploration: Oil from chalks, a modern miracle?. Bull. Am. Assoc. Petrol. Geol. 61, 982–1009 (1977).CAS 

Google Scholar 
Tagliavento, M., John, C. M., Anderskouv, K. & Stemmerik, L. Towards a new understanding of the genesis of chalk: Diagenetic origin of micarbs confirmed by clumped isotope analysis. Sedimentology 68, 513–530 (2021).CAS 
Article 

Google Scholar 
Bramlette, M. N. Significance of coccolithophorids in calcium-carbonate deposition. Bull. Geol. Soc. Am. 69, 121–126 (1958).Article 

Google Scholar 
Hattin, D. E. & Darko, D. A. Technique for determining coccolith abundance in shaly chalk of Greenhorn Limestone (Upper Cretaceous) of Kansas. Kansas Geol. Surv. Bull. 202, 1–11 (1971).
Google Scholar 
Houghton, S. D. Calcareous nannofossils. In Calcareous algae and Stromatolites (ed. Riding, R.) 217–266 (Springer, 1991).Chapter 

Google Scholar 
Bown, P. R., Lees, J. A. & Young, J. R. Calcareous nannoplankton evolution and diversity through time. In Coccolithophores—From Molecular Processes to Global Impact (eds Thierstein, H. R. & Young, J. R.) 481–508 (Springer, 2004).
Google Scholar 
Roth, P. H. Mesozoic paleoceanography of the North Atlantic and Tethys Oceans. In North Atlantic Paleoceanography (eds Summerhayes, C. P. & Shackleton, N. J.) 299–320 (Geological Society Special Publications, 1986).
Google Scholar 
Baumann, K.-H., Andruleit, H., Böckel, B., Geisen, M. & Kinkel, H. The significance of extant coccolithophores as indicators of ocean water masses, surface water temperature, and paleoproductivity: A review. Paläontol. Z. 79, 93–112 (2005).Article 

Google Scholar 
Miller, K. G. et al. The phanerozoic record of global sea-level change. Science 310, 1293–1298 (2005).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Ando, A., Huber, B. T., MacLeod, K. G. & Watkins, D. K. Early Cenomanian “hot greenhouse” revealed by oxygen isotope record of exceptionally well-preserved foraminifera from Tanzania. Paleoceanography 30, 1556–1572 (2015).ADS 
Article 

Google Scholar 
Ekdale, A. A. & Bromley, R. G. Comparative ichnology of shelf-sea and deep-sea chalk. J. Paleontol. 58, 322–332 (1984).
Google Scholar 
Savrda, C. E. Chalk and related deep-marine carbonates. In Trace Fossils as Indicators of Sedimentary Environments (eds Knaust, D. & Bromley, R. G.) 777–806 (Elsevier, 2012).Chapter 

Google Scholar 
Savrda, C. E., Foster, C. & Fluegeman, R. A unique Lower Paleocene shelf-sea chalk in the eastern U.S. Gulf coastal plain (Clayton Formation, western Alabama): Implications for depositional environment, sea-level dynamics and paleogeography. Palaeogeogr. Palaeoclimatol. Palaeoecol. 538, 109439 (2020).Article 

Google Scholar 
Erba, E., Watkins, D. & Mutterlose, J. Campanian dwarf calcareous nannofossils from Wodejebato Guyot. In Proc. Ocean Drill. Program Sci. Results (eds Haggerty, J. A. et al.) 141–155 (Ocean Drilling Program, 1995).
Google Scholar 
Hancock, J. M. The petrology of chalk. Proc. Geol. Assoc. 86, 499–535 (1975).Article 

Google Scholar 
Stanley, S. M. & Hardie, L. A. Secular oscillations in the carbonate mineralogy of reef-building and sediment-producing organisms driven by tectonically forced shifts in seawater chemistry. Palaeogeogr. Palaeoclimatol. Palaeoecol. 144, 3–19 (1998).Article 

Google Scholar 
Stanley, S. M., Ries, J. B. & Hardie, L. A. Seawater chemistry, coccolithophore population growth, and the origin of Cretaceous chalk. Geology 33, 593–596 (2005).ADS 
CAS 
Article 

Google Scholar 
Pemberton, S. G. et al. Ichnology and sedimentology of shallow to marginal marine systems: Ben Nevis and Avalon Reservoirs, Jeanne d’Arc Basin. Geol. Assoc. Can. Short Course Notes 15, 1–343 (2001).
Google Scholar 
Buatois, L. A. & Mángano, M. G. Ichnology: Organism–Substrate Interactions in Space and Time (Cambridge Press University, 2011).Book 

Google Scholar 
Frey, R. W. & Bromley, R. G. Ichnology of American chalks: The Selma Group (Upper Cretaceous), western Alabama. Can. J. Earth Sci. 22, 801–828 (1985).ADS 
Article 

Google Scholar 
Savrda, C. E. & Bottjer, D. Trace-fossil model for reconstructing oxygenation histories of ancient marine bottom waters: Application to Upper Cretaceous Niobrara Formation, Colorado. Palaeogeogr. Palaeoclimatol. Palaeoecol. 74, 49–74 (1989).Article 

Google Scholar 
Kennedy, W. J. Trace fossils in carbonate rocks. In The Study of Trace Fossils (ed. Frey, R. W.) 377–398 (Springer, 1975).Chapter 

Google Scholar 
Loucks, R. G., Gates, B. G. & Zahm, C. K. Depositional systems, lithofacies, nanopore to micropore matrix network, and reservoir quality of the Upper Cretaceous (Cenomanian) Buda Limestone in Dimmit County, southwestern Texas. Gulf Coast Assoc. Geol. Soc. 8, 281–300 (2019).
Google Scholar 
Valencia, F. L. et al. Depositional environments and controls on the stratigraphic architecture of the Cenomanian Buda Limestone in west Texas, U.S.A. Mar. Petrol. Geol. 133, 105275 (2021).Article 

Google Scholar 
Valencia, F. L., Laya, J. C., Buatois, L. A., Mángano, M. G. & Valencia, G. L. Sedimentology and stratigraphy of the Cenomanian Buda Limestone in central Texas, U.S.A.: Implications on regional and global depositional controls. Cretac. Res. 137, 105231 (2022).Article 

Google Scholar 
Martin, K. G. Stratigraphy of the Buda Limestone, south-central Texas. In Comanchean (Lower Cretaceous) Stratigraphy and Paleontology of Texas (ed. Hendricks, L.) 287–299 (Permian Basin Section SEPM 67 (8), 1967).
Google Scholar 
Mallon, A. J. & Swarbrick, R. E. Diagenetic characteristics of low permeability, non-reservoir chalks from the Central North Sea. Mar. Petrol. Geol. 25, 1097–1108 (2008).CAS 
Article 

Google Scholar 
Brasher, J. E. & Vagle, K. R. Influence of lithofacies and diagenesis on Norwegian North Sea chalk reservoirs. Am. Assoc. Petrol. Geol. Bull. 80, 746–769 (1996).CAS 

Google Scholar 
Hentz, T. F. & Ruppel, S. C. Regional stratigraphic and rock characteristics of eagle ford shale in its play area: Maverick Basin to East Texas Basin. Am. Ass. Petrol. Geol. Search and Discovery 10325 (2011).Robinson, W. C. Petrography and depositional environments of the Buda Limestone, northern Coahuila, Mexico. MS Thesis. The University of Texas, 156 (1982).Reaser, D. F. & Robinson, W. C. Cretaceous Buda Limestone in west Texas and northern Mexico. In Cretaceous Stratigraphy and Paleoecology, Texas and Mexico (ed. Scott, R. W.) 337–373 (Perkins Memorial volume, GCSSEPM Foundation, Special Publications in Geology 1, 2003).
Google Scholar 
Young, K. P. Cretaceous paleogeography: Implications of endemic ammonite faunas. Geol. Circ. (University of Texas at Austin, Bureau of Economic Geology) 72, 1–13 (1972).
Google Scholar 
Buatois, L. A. & Mángano, M. G. Ichnodiversity and ichnodisparity: Significance and caveats. Lethaia 46, 281–292 (2013).Article 

Google Scholar 
Buatois, L. A., Wisshak, M., Wilson, M. A. & Mángano, M. G. Categories of architectural designs in trace fossils: A measure of ichnodisparity. Earth Sci. Rev. 164, 102–181 (2017).ADS 
CAS 
Article 

Google Scholar 
Swinbanks, D. D. & Luternauer, J. L. Burrow distribution of thalassinidean shrimp on a Fraser Delta tidal flat, British Columbia. J. Paleontol. 61, 315–333 (1987).Article 

Google Scholar 
Carmona, N. B., Buatois, L. A. & Mángano, M. G. The trace fossil record of burrowing decapod crustaceans: Evaluating evolutionary radiations and behavioural convergence. In Trace Fossils in Evolutionary Palaeoecology (eds Webby, B. D. et al.) 141–153 (Wiley, 2004).
Google Scholar 
Baucon, A. et al. Ethology of the trace fossil Chondrites: Form, function and environment. Earth Sci. Rev. 202, 102989 (2020).CAS 
Article 

Google Scholar 
Pemberton, S. G. & Frey, R. W. Trace fossil nomenclature and the Planolites–Palaeophycus dilemma. J. Paleontol. 56, 843–881 (1982).
Google Scholar 
Rodríguez-Tovar, F. J. & Pérez-Valera, F. Trace fossil Rhizocorallium from the Middle Triassic of the Betic Cordillera, Southern Spain: Characterization and environmental implications. Palaios 23, 78–86 (2008).ADS 
Article 

Google Scholar 
Bown, T. M. & Kraus, M. J. Ichnofossils of the alluvial Willwood Formation (lower Eocene), Bighorn Basin, northwest Wyoming, USA. Palaeogeogr. Palaeoclimatol. Palaeoecol 43, 95–128 (1983).Article 

Google Scholar 
Uchman, A. Taxonomy and palaeoecology of flysch trace fossils: The Marnoso-arenacea Formation and associated facies (Miocene, Northern Apennines, Italy). Beringeria 15, 3–115 (1995).
Google Scholar 
Demírcan, H. & Uchman, A. The miniature trace fossil Bichordites kuzunensis isp. Nov., from early Oligocene prodelta sediments of the Mezardere Formation, Gökçeada Island, NW Turkey. Acta Geol. Pol. 62, 205–215 (2012).
Google Scholar 
Plaziat, J.-C. & Mahmoudi, M. Trace fossils attributed to burrowing echinoids: A revision including new ichnogenus and ichnospecies. Geobios 21, 209–233 (1988).Article 

Google Scholar 
Chamberlain, C. K. Morphology and ethology of trace fossils from the Ouachita Mountains, southeast Oklahoma. J. Paleontol. 45, 212–246 (1971).
Google Scholar 
Farrow, G. E. Bathymetric zonation of Jurassic trace fossils from the coast of Yorkshire, England. Palaeogeogr. Palaeoclimatol. Palaeoecol. 2, 103–151 (1966).Article 

Google Scholar 
Mángano, M. G., Buatois, L. A., West, R. R. & Maples, C. G. Contrasting behavioral and feeding strategies recorded by tidal-flat bivalve trace fossils from the upper carboniferous of eastern Kansas. Palaios 13, 335–351 (1998).ADS 
Article 

Google Scholar 
Pemberton, S. G., Frey, R. W. & Bromley, R. G. The ichnotaxonomy of Conostichus and other plug-shaped ichnofossils. Can. J. Earth Sci. 25, 866–892 (1988).ADS 
Article 

Google Scholar 
Nara, M. Rosselia socialis: A dwelling structure of a probable terebellid polychaete. Lethaia 28, 171–178 (1995).Article 

Google Scholar 
Wilson, M. A., Curran, H. A. & White, B. Paleontological evidence of a brief global sea-level event during the last interglacial. Lethaia 31, 241–250 (1998).Article 

Google Scholar 
Santos, A., Mayoral, E., Marques da Silva, C., Cachão, M. & Kullberg, J. C. Trypanites ichnofacies: Palaeoenvironmental and tectonic implications. A case study from the Miocene disconformity at Foz da Fonte (Lower Tagus Basin, Portugal). Palaeogeogr. Palaeoclimatol. Palaeoecol. 292, 35–43 (2010).Article 

Google Scholar 
Wilson, J. L. Carbonate Facies in Geological History (Springer, 1975).Book 

Google Scholar 
Tucker, M. E. & Wright, V. P. Carbonate Sedimentology (Blackwell Science, 1990).Book 

Google Scholar 
MacEachern, J. A. & Gingras, M. K. Recognition of brackish-water trace fossil assemblages in the Cretaceous western interior seaway of Alberta. In Sediment-Organism Interactions: A Multifaceted Ichnology (eds Bromley, R. G. et al.) 149–194 (Society for Sedimentary Geology Special Publication, 2007).
Google Scholar 
MacEachern, J. A., Zaitlin, B. A. & Pemberton, S. G. High-resolution sequence stratigraphy of early transgressive deposits, Viking Formation, Joffre Field, Alberta, Canada. Bull. Am. Assoc. Petrol. Geol. 82, 729–756 (1998).
Google Scholar 
Buatois, L. A., Netto, R. G. & Mángano, M. G. Ichnology of Permian marginal-marine to shallow-marine coal-bearing successions: Rio Bonito and Palermo formations, Parana Basin, Brazil. In Applied Ichnology (eds MacEachern, J. A. et al.) 167–177 (Society for Sedimentary Geology Short Course Notes, 2007).
Google Scholar 
Buatois, L. A. et al. Colonization of brackish-water systems through time: Evidence from the trace-fossil record. Palaios 20, 321–347 (2005).ADS 
Article 

Google Scholar 
Pemberton, S. G. & Wightman, D. M. Ichnological characteristics of brackish water deposits. In Applications of Ichnology to Petroleum Exploration: A Core Work-shop (ed. Pemberton, S. G.) 141–167 (Society of Economic Paleontologists and Mineralogists Core Workshop, 1992).Chapter 

Google Scholar 
Anderson, B. G. & Droser, M. L. Ichnofabrics and geometric configurations of Ophiomorpha within a sequence stratigraphic framework: An example from the Upper Cretaceous US western interior. Sedimentology 45, 379–396 (1998).ADS 
Article 

Google Scholar 
Buatois, L. A., Mángano, M. G. & Pattison, S. A. J. Ichnology of prodeltaic hyperpycnite–turbidite channel complexes and lobes from the Upper Cretaceous Prairie Canyon Member of the Mancos Shale, Book Cliffs, Utah, USA. Sedimentology 66, 1825–1860 (2019).Article 

Google Scholar 
Bhattacharya, J. P. & MacEachern, J. A. Hyperpycnal rivers and prodeltaic shelves in the Cretaceous seaway of North America. J. Sediment. Res. 79, 184–209 (2009).ADS 
Article 

Google Scholar 
Savrda, C. E. Ichnosedimentologic evidence for a noncatastrophic origin of Cretaceous-Tertiary boundary sand in Alabama. Geology 21, 1075–1078 (1993).ADS 
Article 

Google Scholar 
Schlager, W. Accommodation and supply-a dual control on stratigraphic sequences. Sediment. Geol. 86, 111–136 (1993).ADS 
Article 

Google Scholar 
Strasser, A. & Samankassou, E. Carbonate sedimentation rates today and in the past: Holocene of Florida Bay, Bahamas, and Bermuda vs. Upper Jurassic and Lower Cretaceous of the Jura Mountains (Switzerland and France). Geol. Croat. 56, 1–18 (2003).Article 

Google Scholar 
Moyano-Paz, D., Richiano, S., Varela, A. N., Gómez-Dacal, A. R. & Poire, D. G. Ichnological signatures from wave- and fluvial-dominated deltas: The La Anita Fromation, Upper Cretaceous, Austral-Magallanes Basin, Patagonia. Mar. Pet. Geol. 114, 104168 (2020).CAS 
Article 

Google Scholar 
De Gibert, J. M. & Ekdale, A. A. Trace fossil assemblages reflecting stressed environments in the Middle Jurassic Carmel Seaway of Central Utah. J. Paleontol. 73, 711–720 (1999).Article 

Google Scholar 
Gingras, M. K., MacEachern, J. A. & Dashtgard, S. E. Process ichnology and the elucidation of physico-chemical stress. Sediment. Geol. 237, 115–134 (2011).ADS 
CAS 
Article 

Google Scholar 
Smith, C. R., Levin, L. A., Hoover, D. J., McMurty, G. & Gage, J. D. Variations in bioturbation across the oxygen minimum zone in the northwest Arabian Sea. Deep-Sea Res. II 47, 227–257 (2000).ADS 
CAS 
Article 

Google Scholar 
Wignall, P. B., Newton, R. & Brookfield, M. E. Pyrite framboid evidence for oxygen-poor deposition during the Permian-Triassic crisis in Kashmir. Palaeogeogr. Palaeoclimatol. Palaeoecol. 216, 183–188 (2005).Article 

Google Scholar 
Kennedy, W. J. Burrows and surface traces from the Lower Chalk of southern England. Bull. Br. Mus. Nat. Hist. Geol. 15, 127–167 (1967).
Google Scholar 
Kennedy, W. J. & Garrison, R. E. Morphology and genesis of nodular chalks and hardgrounds in the Upper Cretaceous of southern England. Sedimentology 22, 311–386 (1975).ADS 
CAS 
Article 

Google Scholar 
Bromley, R. G. Some observations on burrows of thalassinidean Crustacea in chalk hardgrounds. Geol. Soc. Lond. Q. J. 123, 157–182 (1967).Article 

Google Scholar 
Bromley, R. G. Trace fossils at omission surfaces. In The Study of Trace Fossils (ed. Frey, R. W.) 399–428 (Springer, 1975).Chapter 

Google Scholar 
Hart, M. B., Harries, P. J. & Cárdenas, A. L. The Cretaceous/Paleogene boundary events in the Gulf Coast: Comparisons between Alabama and Texas. Gulf Coast Assoc. Geol. Trans. 63, 235–255 (2013).
Google Scholar 
Al Balushi, S. A. K. & Macquaker, J. H. S. Sedimentological evidence for bottom-water oxygenation during deposition of the Natih-B Member intrashelf-basinal sediments: Upper Cretaceous carbonate source rock, Natih Formation, North Sultanate of Oman. GeoArabia 16, 47–84 (2011).Article 

Google Scholar 
Lasseur, E. et al. A relative water-depth model for the Normandy Chalk (Cenomanian–Middle Coniacian, Paris Basin, France) based on facies patterns of metre-scale cycles. Sediment. Geol. 213, 1–26 (2009).ADS 
Article 

Google Scholar 
Dawson, W. C. & Reaser, D. F. Rhizocorallium in the upper Austin Chalk, Ellis County, Texas. Texas J. of Sci. 23, 207–214 (1980).
Google Scholar 
Dawson, W. C. & Reaser, D. F. Ichnology and paleoenvironments of the middle and upper Austin Chalk (Upper Cretaceous), northeastern Texas. Trans. Am. Assoc. Pet. Geol. Southwest Sec. 1985, 47–67 (1985).
Google Scholar 
Dawson, W. C. & Reaser, D. F. Trace fossils and paleoenvironments of lower and middle Austin Chalk (Upper Cretaceous), north-central Texas. Trans. Gulf Coast Assoc. Geol. Soc. 40, 161–173 (1990).
Google Scholar 
Dawson, W. C. & Reaser, D. F. Ichnology and Paleosubstrates of Austin Chalk (Cretaceous) Outcrops: Southern Dallas and Ellis Counties, Texas. Am. Assoc. Pet. Geol. Search Discovery Article #91004 (1991).Fürsich, F. T., Kennedy, W. J. & Palmer, T. J. Trace fossils at a regional discontinuity surface: The Austin/Taylor (Upper Cretaceous) contact in central Texas. J. Paleontol. 55, 537–551 (1981).
Google Scholar 
Morgan, R. F. A new ichnospecies of Gyrolithes from the Austin Chalk, Upper Cretaceous, Texas, USA. Ichnos 26, 1–7 (2018).Article 

Google Scholar 
Cooper, J. R., Godet, A. & Pope, M. C. Tectonic and eustatic impact on depositional features in the upper Cretaceous Austin Chalk Group of south-central Texas, USA. Sediment. Geol. 401, 105632 (2020).Article 

Google Scholar 
Loucks, R. G. et al. Geologic characterization of the type cored section for the Upper Cretaceous Austin Chalk Group in southern Texas: A combination fractured and unconventional reservoir. Am. Assoc. Pet. Geol. Bull. 104, 2209–2245 (2020).
Google Scholar 
Loucks, R. G., Reed, R. M., Ko, L. T., Zahm, C. K. & Larson, T. E. Micropetrographic characterization of a siliciclastic-rich chalk; Upper Cretaceous Austin Chalk Group along the onshore northern Gulf of Mexico, USA. Sediment. Geol. 412, 105821 (2021).CAS 
Article 

Google Scholar 
Bottjer, D. J. Paleoecology, Ichnology, and Depositional Environments of Upper Cretaceous Chalks (Annona Formation; chalk Member of Saratoga Formation), Southwestern Arkansas. PhD Dissertation, Indiana University, 424 (1978).Bottjer, D. J. Ichnology and depositional environments of Upper Cretaceous chalks, southwestern Arkansas (Annona Formation; chalk member, Saratoga Formation). Am. Assoc. Pet. Geol. Bull. 63, 422 (1979).
Google Scholar 
Bottjer, D. J. Trace fossils and paleoenvironments of two Arkansas Upper Cretaceous discontinuity surfaces. J. Paleontol. 59, 282–298 (1985).
Google Scholar 
Bottjer, D. J. Campanian-Maastrichtian chalks of southwestern Arkansas: Petrology, paleoenvironments and comparison with other North American and European chalks. Cretac. Res. 7, 161–196 (1986).Article 

Google Scholar 
Bayet-Goll, A., Neto de Carvalho, C., Monaco, P. & Sharafi, M. Sequence stratigraphic and sedimentologic significance of biogenic structures from chalky limestones of the Turonian-Campanian Abderaz Formation, Kopet-Dagh, Iran. In Cretaceous Period: Biotic Diversity and Biogeography (eds Khosla, A. & Lucas, S. G.) 19–43 (New Mex. Mus. Nat. His. Sci. Bull. 71, 2016).
Google Scholar 
Locklair, R. E. & Savrda, C. E. Ichnology of rhythmically bedded Demopolis Chalk (Upper Cretaceous, Alabama): Implications for paleoenvironment, depositional cycle origins, and tracemaker behavior. Palaios 13, 423–438 (1998).ADS 
Article 

Google Scholar 
Locklair, R. E. & Savrda, C. E. Ichnofossil tiering analysis of a rhythmically bedded chalk-marl sequence in the Upper Cretaceous of Alabama. Lethaia 31, 311–322 (1998).Article 

Google Scholar 
Kennedy, W. J. Trace fossils in the chalk environment. In Trace Fossils (eds Crimes, T. P. & Harper, J. C.) 263–282 (Geological Journal Special Issue 3, 1970).
Google Scholar 
Mortimore, R. N. & Pomerol, B. Stratigraphy and eustatic implications of trace fossil events in the Upper Cretaceous Chalk of northern Europe. Palaios 6, 216–231 (1991).ADS 
Article 

Google Scholar 
Foster, C. B. III. Geology of the Moscow Landing Section, Tombigbee River, Western Alabama, with Focus on Ichnologic Aspects of the Lower Paleocene Clayton Formation. M.Sc. Dissertation, Auburn University, 88 (2019).Gabdullin, R. R. Rhythmicity of the Upper Cretaceous Deposits of the East European Craton, Northwestern Caucasus and Southwestern Crimea: Structure, Classification, Formation Models (Mosk. Gos. Univ., 2002).
Google Scholar 
Baraboshkin, E. Y. & Zibrov, I. A. Characteristics of the Middle Cenomanian Rhythmic Sequence from Mount Selbukhra in Southwest Crimea. Moscow Univ. Geol. Bull. 67, 176–184 (2012).Article 

Google Scholar 
Blinkenberg, K. H., Lauridsen, B. W., Knaust, D. & Stemmerik, L. New ichnofabrics of the Cenomanian-Danian Chalk Group. J. Sediment. Res. 90, 701–712 (2020).ADS 
Article 

Google Scholar 
Ekdale, A. A. & Bromley, R. G. Trace fossils and ichnofabric in the Kjolby Gaard Marl, uppermost Cretaceous, Denmark. Bull. Geol. Soc. Denmark 31, 107–119 (1983).Article 

Google Scholar 
Ekdale, A. A. & Bromley, R. G. Cretaceous chalk ichnofacies in northern Europe. Geobios 8, 201–204 (1984).Article 

Google Scholar 
Ekdale, A. A. & Bromley, R. G. Analysis of composite ichnofabrics; An example in Uppermost Cretaceous chalk of Denmark. Palaios 6, 232–249 (1991).ADS 
Article 

Google Scholar 
Surlyk, F. et al. The cyclic Rørdal Member—A new lithostratigraphic unit of chronostratigraphic and palaeoclimatic importance in the upper Maastrichtian of Denmark. Bull. Geol. Soc. Denmark 58, 89–98 (2010).Article 

Google Scholar 
Lauridsen, B. W., Surlyk, F. & Bromley, R. G. Trace fossils of a cyclic chalk marl succession; the upper Maastrichtian Rørdal Member, Denamrk. Cretac. Res. 32, 194–211 (2011).Article 

Google Scholar 
Frey, R. W. Trace fossils of Fort Hays Limestone Member of Niobrara Chalk (Upper Cretaceous), west-central Kansas. Univ. Kansas Paleontol. Contrib. 53, 52 (1970).
Google Scholar 
Hattin, D. E. Stratigraphy and depositional environment of Smoky Hill Chalk Member, Niobrara Chalk (Upper Cretaceous) of the type area western Kansas. Kansas Geol. Surv. Bull. 225, 1–108 (1982).
Google Scholar 
Savrda, C. E. Ichnocoenoses in the Niobrara Formation: Implications for benthic oxygenation histories. In Stratigraphy and Paleoenvironments of the Cretaceous Western Interior Seaway, USA (eds Dean, W. E. & Arthur, M. A.) 137–151 (SEPM Society for Sedimentary Geology 6, 1998).Chapter 

Google Scholar 
Hattin, D. E. Widespread, synchronously deposited, burrow-mottled limestone beds in Greenhorn Limestone (Upper Cretaceous) of Kansas and southeastern Colorado. Am. Assoc. Pet. Geol. Bull. 55, 412–431 (1971).
Google Scholar 
Hattin, D. E. Stratigraphy and depositional environment of Greenhorn Limestone (Upper Cretaceous) of Kansas. Kansas Geol. Surv. Bull. 209, 128 (1975).
Google Scholar 
Savrda, C. E. Ichnology of the Bridge Creek Limestone: Evidence for temporal and spatial variations in paleo-oxygenation in the Western Interior Seaway. In Stratigraphy and Paleoenvironments of the Cretaceous Western Interior Seaway, USA (eds Dean, W. E. & Arthur, M. A.) 127–136 (SEPM Society for Sedimentary Geology 6, 1998).Chapter 

Google Scholar 
Rasmussen, S. L. & Surlyk, F. Facies and ichnology of an Upper Cretaceous chalk contourite drift complex, eastern Denmark, and the validity of contourite facies models. J. Geol. Soc. Lond. 169, 435–447 (2012).Article 

Google Scholar 
Surlyk, F. et al. Upper Campanian-Maastrichtian holostratigraphy of the eastern Danish Basin. Cretac. Res. 46, 232–256 (2013).Article 

Google Scholar 
Boussaha, M., Thibault, N., Anderskouv, K., Moreau, J. & Stemmerik, L. Controls on upper Campanian-Maastrichtian chalk deposition in the eastern Danish Basin. Sedimentology 64, 1998–2030 (2017).Article 

Google Scholar 
Reolid, J. & Betzler, C. The ichnology of carbonate drifts. Sedimentology 66, 1427–1448 (2019).Article 

Google Scholar 
Nygaard, E. Bathichnus and Its Significance in the Trace Fossil Association of Upper Cretaceous Chalk, Mors, Denmark 107–113 (Danm. Geol. Unders. Årbog, 1983).
Google Scholar 
Scholle, P. A., Albrechtsen, T. & Tirsgaard, H. Formation and diagenesis of bedding cycles in uppermost Cretaceous chalks of the Dan Field, Danish North Sea. Sedimentology 45, 223–243 (1998).ADS 
CAS 
Article 

Google Scholar 
Damholt, T. & Surlyk, F. Laminated–bioturbated cycles in Maastrichtian chalk of the North Sea: Oxygenation fluctuations within the Milankovitch frequency band. Sedimentology 51, 1323–1342 (2004).ADS 
Article 

Google Scholar 
Anderskouv, K. & Surlyk, F. Upper Cretaceous chalk facies and depositional history recorded in the Mona-1 core, Mona Ridge, Danish North Sea. Geol. Surv. Denmark Greenland Bull. 25, 1–60 (2011).Article 

Google Scholar 
Maliva, R. G. & Dickson, J. A. D. Microfacies and diagenetic controls of porosity in Cretaceous/Tertiary chalks, Eldfisk Field, Norwegian North Sea. Am. Assoc. Pet. Geol. Bull. 76, 1825–1838 (1992).
Google Scholar 
Knaust, D., Dorador, J. & Rodríguez-Tovar, F. J. Burrowed matrix powering dual porosity systems—A case study from the Maastrichtian chalk of the Gullfaks Field Norwegian North Sea. Mar. Petrol. Geol. 113, 104158 (2020).Article 

Google Scholar 
Phillips, C. & McIlroy, D. Ichnofabrics and biologically mediated changes in clay mineral assemblages from a deep-water, fine-grained, calcareous sedimentary succession: An example from the Upper Cretaceous Wyandot Formation, offshore Nova Scotia. Bull. Can. Petrol. Geol. 58, 203–218 (2010).Article 

Google Scholar 
Rodríguez-Tovar, F. J. & Hernández-Molina, F. J. Ichnological analysis of contourites: Past, present and future. Earth-Sci. Rev. 182, 28–41 (2018).ADS 
Article 

Google Scholar 
Miguez-Salas, O. & Rodríguez-Tovar, F. J. Ichnofacies distribution in the Eocene-Early Miocene Petra Tou Romiou outcrop, Cyprus: Sea level dynamics and palaeoenvironmental implications in a contourite environment. Int. J. Earth Sci. 108, 2531–2544 (2019).CAS 
Article 

Google Scholar 
Nelson, C. S. Bioturbation in middle bathyal, Cenozoic nannofossil oozes and chalks, southwest Pacific. In Initial Reports of the Deep Sea Drilling Project 90 (eds Kennett, J. P., von der Borch, C. C. et al.) 1189–1200 (Washington U.S. Government Printing Office, 1986).
Google Scholar 
Fütterer, D. K. Bioturbation and trace fossils in deep sea sediments of the Walvis Ridge, southeastern Atlantic, Leg 74. In Initial Reports of the Deep Sea Drilling Project 74 (eds Moore, T. C., Rabinowitz, P. D. et al.) 543–555 (Government Printing Office, 1984).
Google Scholar 
Wetzel, A. Ichnofabrics in Eocene to Maestrichtian sediments from Deep Sea Drilling Project Site 605, off the New Jersey coast. In Initial Reports of the Deep Sea Drilling Project 93 (eds. Hinte, J. E., Wise Jr., S. W. et al.) 825–835 (1987).Droser, M. L. & Bottjer, D. J. Trace fossils and ichnofabrics in Leg 119 cores. In Proceedings of the Ocean Drilling Program, Scientific Results 119 (eds. Barron, J., Larsen, B. et al.) 635–641 (1991).Desai, B. G. Ichnofabric analysis of bathyal chalks: The Miocene Inglis Formation of the Andaman and Nicobar Islands, India. J. Palaeogeogr. 10, 1–15 (2021).Article 

Google Scholar 
Warme, J. E., Kennedy, W. J. & Scheidermann, N. Biogenic sedimentary structures (trace fossils) in Leg 15 cores. In Initial Reports of the Deep Sea Drilling Project 15 (eds. Edgar, N. T., Saunders, J. B. et al.) 813–831 (1973).Maurrasse, F. Sedimentary structures of Caribbean Leg 15 sediments. In Initial Reports of the Deep-Sea Drilling Project 15 (eds. Edgar, T. et al.) (1974).Erba, E. & Premoli-Silva, I. Orbitally driven cycles in trace-fossil distribution from the Piobbico core (late Albian, central Italy). In Orbital Forcing and Cyclic Sequences, IAS Spec. Publ. 19 (eds De Boer, P. L. & Smith, D. G.) 211–225 (Blackwell Scientific, 1994).
Google Scholar 
Chamberlain, C. K. Trace fossils in DSDP cores of the Pacific. J. Paleontol. 49, 1074–1096 (1975).
Google Scholar 
Ekdale, A. A. Trace fossils in Deep Sea Drilling Project Leg 58 cores. In Initial Reports of the Deep Sea Drilling Project 58 (eds. de Vries Klein, G., Kobyashi, K. et al.) 601–605 (1980).Ekdale, A. A. Geologic history of the abyssal benthos: Evidence from trace fossils in Deep-Sea Drilling Project cores. PhD Dissertation, Rice University, 154 (1974).Ekdale, A. A. Abyssal trace fossils in worldwide Deep Sea Drilling Project cores. In Trace Fossils 2 (eds. Crimes, T. P. & Harper, J. C.) 163–182 (Geol. J., Spec. Iss. 9, 1977).Ekdale, A. A. & Berger, W. H. Deep-sea ichnofacies: Modern organism traces on and in pelagic carbonates of the western equatorial Pacific. Palaeogeogr. Palaeoclimatol. Palaeoecol. 23, 263–278 (1978).Article 

Google Scholar 
Ekdale, A. A., Muller, L. N. & Novak, M. T. Quantitative ichnology of modern pelagic deposits in the abyssal Atlantic. Palaeogeogr. Palaeoclimatol. Palaeoecol. 45, 189–223 (1984).CAS 
Article 

Google Scholar 
Savrda, C. E. Limited ichnologic fidelity and temporal resolution in pelagic sediments: Paleoenvironmental and paleoecologic implications. Palaios 29, 210–217 (2014).ADS 
Article 

Google Scholar 
Bromley, R. G. & Ekdale, A. A. Composite ichnofabrics and tiering of burrows. Geol. Mag. 123, 59–65 (1986).ADS 
Article 

Google Scholar 
Griffin, J. N. et al. Spatial heterogeneity increases the importance of species richness for an ecosystem process. Oikos 118, 1335–1342 (2009).Article 

Google Scholar 
Valentine, J. W. Overview of marine biodiversity. In Marine Macroecology (eds Witman, J. D. & Roy, K.) 3–28 (University of Chicago Press, 2009).Chapter 

Google Scholar 
Schlacher, T. A. et al. Soft-sediment benthic community structure in a coral reef lagoon—The prominence of spatial heterogeneity and “spot endemism”. Mar. Ecol. Prog. Ser. 174, 159–174 (1998).ADS 
Article 

Google Scholar 
Hummel, H. et al. Geographic patterns of biodiversity in European coastal marine benthos. J. Mar. Biol. Assoc. U.K. 97, 507–523 (2017).Article 

Google Scholar 
Harborne, A. R., Mumby, P. J., Żychaluk, K., Hedley, J. D. & Blackwell, P. G. Modeling the beta diversity of coral reefs. Ecology 87, 2871–2881 (2006).PubMed 
Article 

Google Scholar 
Christia, C., Giordani, G. & Papastergiadou, E. Environmental variability and macrophyte assemblages in coastal lagoon types of Western Greece (Mediterranean Sea). Water 10, 151 (2018).Article 
CAS 

Google Scholar 
Dorador, J., Rodríguez-Tovar, F. J., IODP Expedition 339 Scientists. Digital image treatment applied to ichnological analysis of marine core sediments. Facies 60, 39–44 (2014).Article 

Google Scholar 
Dorador, J. & Rodríguez-Tovar, F. J. High-resolution image treatment in ichnological core analysis: Initial steps, advances and prospects. Earth-Sci. Rev. 177, 226–237 (2018).ADS 
Article 

Google Scholar 
Taylor, A. M. & Goldring, R. Description and analysis of bioturbation and ichnofabric. J. Geol. Soc. 150, 141–148 (1993).ADS 
Article 

Google Scholar 
Cao, Y. M., Curran, A. H. & Glumac, B. Testing the use of photoshop and imageJ for evaluating ichnofabrics. 2015 GSA Annual Meeting in Baltimore, Maryland, USA, Paper No. 128-11 (The Geol. Soc. of Am., 2015). More