Ecology

1.Fernandez, R. J. Do humans create deserts?. Trends Ecol. Evol. 17, 6–7 (2002).
Google Scholar 
2.Geist, H. J. & Lambin, E. F. Dynamic causal patterns of desertification. Bioscience 54(9), 817–829 (2004).
Google Scholar 
3.Imeson, A. Desertification, Land Degradation and Sustainability (Routledge, 2012).
Google Scholar 
4.Romm, J. Desertification: The next dust bowl. Nature 478, 450–451 (2011).ADS 
CAS 
PubMed 

Google Scholar 
5.Portnov, B. A. & Safriel, U. N. Combating desertification in the Negev: Dryland agriculture vs. dryland urbanization. J. Arid Environ. 56, 659–680 (2004).ADS 

Google Scholar 
6.Salvati, L., Bajocco, S., Ceccarelli, T., Zitti, M. & Perini, L. Towards a process-based evaluation of land vulnerability to soil degradation in Italy. Ecol. Ind. 11(5), 1216–1227 (2011).CAS 

Google Scholar 
7.Santini, M., Caccamo, G., Laurenti, A., Noce, S. & Valentini, R. A multi-model GIS framework for desertification risk assessment. Appl. Geogr. 30(3), 394–415 (2010).
Google Scholar 
8.Bajocco, S., Salvati, L. & Ricotta, C. Land degradation vs. Fire: A spiral process?. Prog. Phys. Geogr. 35(1), 3–18 (2011).
Google Scholar 
9.Incerti, G., Feoli, E., Salvati, L., Brunetti, A. & Giovacchini, A. Analysis of bioclimatic time series and their neural network-based classification to characterise drought risk patterns in South Italy. Int. J. Biometeorol. 51(4), 253–263 (2007).ADS 
CAS 
PubMed 

Google Scholar 
10.Salvati, L., Perini, L., Sabbi, A. & Bajocco, S. Climate Aridity and land use changes: A regional-scale analysis. Geogr. Res. 50(2), 193–203 (2012).
Google Scholar 
11.Coluzzi, R. et al. Investigating climate variability and long-term vegetation activity across heterogeneous Basilicata agroecosystems. Geomat. Nat. Haz. Risk 10(1), 168–180 (2019).
Google Scholar 
12.Imbrenda, V. et al. Analysis of landscape evolution in a vulnerable coastal area under natural and human pressure. Geomat. Nat. Haz. Risk 9(1), 1249–1279 (2018).
Google Scholar 
13.Imbrenda, V. et al. Land degradation and metropolitan expansion in a peri-urban environment. Geomat. Nat. Haz. Risk 12(1), 1797–1818 (2021).
Google Scholar 
14.Kairis, O., Karavitis, C., Kounalaki, A., Salvati, L. & Kosmas, C. The effect of land management practices on soil erosion and land desertification in an olive grove. Soil Use Manag. 29(4), 597–606 (2013).
Google Scholar 
15.Kairis, O., Karavitis, C., Salvati, L., Kounalaki, A. & Kosmas, K. Exploring the impact of overgrazing on soil erosion and land degradation in a dry Mediterranean agro-forest landscape (Crete, Greece). Arid Land Res. Manag. 29(3), 360–374 (2015).
Google Scholar 
16.Karamesouti, M. et al. Land-use and land degradation processes affecting soil resources: Evidence from a traditional Mediterranean cropland (Greece). CATENA 132, 45–55 (2015).
Google Scholar 
17.Kosmas, C. et al. Land degradation and long-term changes in agro-pastoral systems: An empirical analysis of ecological resilience in Asteroussia-Crete (Greece). CATENA 147, 196–204 (2016).
Google Scholar 
18.Jongman, R. H. G. Homogenisation and fragmentation of the European landscape: Ecological consequences and solutions. Landsc. Urban Plan. 58(2), 211–221 (2002).
Google Scholar 
19.Lavado Contador, J. F., Schnabel, S., Gomez Gutierrez, A. & Pulido Fernandez, M. Mapping sensitivity to land degradation in Extremadura, SW Spain. Land Degrad. Dev. 20(2), 129–144 (2009).
Google Scholar 
20.Otto, R., Krusi, B. O. & Kienast, F. Degradation of an arid coastal landscape in relation to land use changes in southern Tenerife (Canary Islands). J. Arid Environ. 70, 527–539 (2007).ADS 

Google Scholar 
21.Braje, T. J., Leppard, T. P., Fitzpatrick, S. M. & Erlandson, J. M. Archaeology, historical ecology and anthropogenic island ecosystems. Environ. Conserv. 44(3), 286–297 (2017).
Google Scholar 
22.Rick, T., Ontiveros, M. Á. C., Jerardino, A., Mariotti, A., Méndez, C. & Williams, A. N. Human-environmental interactions in Mediterranean climate regions from the Pleistocene to the Anthropocene. Anthropocene, 100253 (2020).23.Bajocco, S., Ceccarelli, T., Smiraglia, D., Salvati, L. & Ricotta, C. Modeling the ecological niche of long-term land use changes: The role of biophysical factors. Ecol. Ind. 60, 231–236 (2016).
Google Scholar 
24.Antrop, M. Landscape change and the urbanization process in Europe. Landsc. Urban Plan. 67(1), 9–26 (2004).
Google Scholar 
25.Bakra, N., Weindorf, D. C., Bahnassy, M. H. & El-Badawi, M. M. Multi-temporal assessment of land sensitivity to desertification in a fragile agro-ecosystem: Environmental indicators. Ecol. Ind. 15(1), 271–280 (2012).
Google Scholar 
26.Pacheco, F. A. L., Fernandes, L. F. S., Junior, R. F. V., Valera, C. A. & Pissarra, T. C. T. Land degradation: Multiple environmental consequences and routes to neutrality. Curr. Opin. Environ. Sci. Health 5, 79–86 (2018).
Google Scholar 
27.Gomes, E. et al. Agricultural land fragmentation analysis in a peri-urban context: From the past into the future. Ecol. Ind. 97, 380–388 (2019).
Google Scholar 
28.Gulcin, D. & Yilmaz, K. T. The assessment of landscape fragmentation in an agricultural environment: Degradation or contribution to ecosystem services?. Fresenius Environ. Bull. 25(12), 7941–7950 (2017).
Google Scholar 
29.Pili, S., Grigoriadis, E., Carlucci, M., Clemente, M. & Salvati, L. Towards sustainable growth? A multi-criteria assessment of (changing) urban forms. Ecol. Ind. 76, 71–80 (2017).
Google Scholar 
30.Haddad, N. M., Brudvig, L. A. & Clobert, J. Habitat fragmentation and its lasting impact on Earth’s ecosystems. Sci. Adv. 1(2), 1–9 (2015).
Google Scholar 
31.Kouba, Y., Gartzia, M., El Aich, A. & Alados, C. L. Deserts do not advance, they are created: Land degradation and desertification in semiarid environments in the Middle Atlas, Morocco. J. Arid Environ. 158, 1–8 (2018).ADS 

Google Scholar 
32.Nagendra, H., Munroe, D. K. & Southworth, J. From pattern to process: Landscape fragmentation and the analysis of land use/land cover change. Agric. Ecosyst. Environ. 101(2), 111–115 (2004).
Google Scholar 
33.Lin, Y., Han, G., Zhao, M. & Chang, S. X. Spatial vegetation patterns as early signs of desertification: A case study of a desert steppe in Inner Mongolia, China. Landsc. Ecol. 25(10), 1519–1527 (2010).
Google Scholar 
34.Kéfi, S. et al. Spatial vegetation patterns and imminent desertification in Mediterranean arid ecosystems. Nature 449(7159), 213–217 (2007).ADS 
PubMed 

Google Scholar 
35.Girvetz, E. H., Thorne, J. H., Berry, A. M. & Jaeger, J. A. Integration of landscape fragmentation analysis into regional planning: A statewide multi-scale case study from California, USA. Landsc. Urban Plan. 86(3), 205–218 (2008).
Google Scholar 
36.Hargis, C. D., Bissonette, J. A. & David, J. L. The behavior of landscape metrics commonly used in the study of habitat fragmentation. Landsc. Ecol. 13(3), 167–186 (1998).
Google Scholar 
37.Llausàs, A. & Nogué, J. Indicators of landscape fragmentation: The case for combining ecological indices and the perceptive approach. Ecol. Ind. 15(1), 85–91 (2012).
Google Scholar 
38.Salvati, L. & Zitti, M. The environmental “risky” region: Identifying land degradation processes through integration of socio-economic and ecological indicators in a multivariate regionalization model. Environ. Manage. 44(5), 888 (2009).ADS 
PubMed 

Google Scholar 
39.Ferrara, A. et al. Updating the MEDALUS-ESA framework for worldwide land degradation and desertification assessment. Land Degrad. Dev. 31(12), 1593–1607 (2020).
Google Scholar 
40.Delfanti, L. et al. Solar plants, environmental degradation and local socioeconomic contexts: A case study in a Mediterranean country. Environ. Impact Assess. Rev. 61, 88–93 (2016).
Google Scholar 
41.Cowie, A. L. et al. Land in balance: The scientific conceptual framework for Land Degradation Neutrality. Environ. Sci. Policy 79, 25–35 (2018).
Google Scholar 
42.Lanfredi, M. et al. A geostatistics-assisted approach to the deterministic approximation of climate data. Environ. Model. Softw. 66, 69–77 (2015).
Google Scholar 
43.Coluzzi, R. et al. Density matters? Settlement expansion and land degradation in Peri-urban and rural districts of Italy. Environ. Impact Assess. Rev. 92, 106703 (2022).
Google Scholar 
44.Xie, H., Zhang, Y., Wu, Z. & Lv, T. A bibliometric analysis on land degradation: Current status, development, and future directions. Land 9(1), 28 (2020).
Google Scholar 
45.Ferrara, A., Salvati, L., Sateriano, A. & Nolè, A. Performance evaluation and costs assessment of a key indicator system to monitor desertification vulnerability. Ecol. Ind. 23, 123–129 (2012).
Google Scholar 
46.Salvati, L. From simplicity to complexity: The changing geography of land vulnerability to degradation in Italy. Geogr. Res. 51(3), 318–328 (2013).
Google Scholar 
47.Salvati, L. et al. Italy’s renewable water resources as estimated on the basis of the monthly water balance. Irrig. Drain. J. Int. Commiss. Irrig. Drain. 57(5), 507–515 (2008).
Google Scholar 
48.Salvati, L. et al. Assessing the effectiveness of sustainable land management policies for combating desertification: A data mining approach. J. Environ. Manage. 183, 754–762 (2016).CAS 
PubMed 

Google Scholar 
49.Recanatesi, F. et al. A fifty-year sustainability assessment of Italian agro-forest districts. Sustainability 8(1), 32 (2016).
Google Scholar 
50.Bajocco, S., De Angelis, A. & Salvati, L. A satellite-based green index as a proxy for vegetation cover quality in a Mediterranean region. Ecol. Ind. 23, 578–587 (2012).
Google Scholar 
51.Smiraglia, D. et al. The latent relationship between soil vulnerability to degradation and land fragmentation: A statistical analysis of landscape metrics in Italy, 1960–2010. Environ. Manage. 64(2), 154–165 (2019).ADS 
PubMed 

Google Scholar 
52.Smiraglia, D., Ceccarelli, T., Bajocco, S., Salvati, L. & Perini, L. Linking trajectories of land change, land degradation processes and ecosystem services. Environ. Res. 147, 590–600 (2016).CAS 
PubMed 

Google Scholar 
53.Zambon, I., Benedetti, A., Ferrara, C. & Salvati, L. Soil matters? A multivariate analysis of socioeconomic constraints to urban expansion in Mediterranean Europe. Ecol. Econ. 146, 173–183 (2018).
Google Scholar 
54.Zambon, I. et al. Land quality, sustainable development and environmental degradation in agricultural districts: A computational approach based on entropy indexes. Environ. Impact Assess. Rev. 64, 37–46 (2017).
Google Scholar 
55.Basso, B. et al. Evaluating responses to land degradation mitigation measures in Southern Italy. Int. J. Environ. Res. 6(2), 367–380 (2012).
Google Scholar 
56.Salvati, L. & Zitti, M. Land degradation in the Mediterranean Basin: Linking bio-physical and economic factors into an ecological perspective. Biota 6, 67–77 (2005).
Google Scholar 
57.Qi, Y. et al. Temporal-spatial variability of desertification in an agro-pastoral transitional zone of northern Shaanxi Province, China. CATENA 88(1), 37–45 (2012).
Google Scholar 
58.Sklenicka, P. Classification of farmland ownership fragmentation as a cause of land degradation: A review on typology, consequences, and remedies. Land Use Policy 57, 694–701 (2016).
Google Scholar 
59.Vos, W. & Meekes, H. Trends in European cultural landscape development: Perspectives for a sustainable future. Landsc. Urban Plan. 46(1), 3–14 (1999).
Google Scholar 
60.Mao, D. et al. Land degradation and restoration in the arid and semiarid zones of China: Quantified evidence and implications from satellites. Land Degrad. Dev. 29(11), 3841–3851 (2018).
Google Scholar 
61.Symeonakis, E., Calvo-Cases, A. & Arnau-Rosalen, E. Land use change and land degradation in southeastern Mediterranean Spain. Environ. Manage. 40(1), 80–94 (2007).ADS 
PubMed 

Google Scholar 
62.Ibanez, J., Martinez Valderrama, J. & Puigdefabregas, J. Assessing desertification risk using system stability condition analysis. Ecol. Model. 213, 180–190 (2008).
Google Scholar 
63.Hill, J., Stellmes, M., Udelhoven, T., Röder, A. & Sommer, S. Mediterranean desertification and land degradation: Mapping related land use change syndromes based on satellite observations. Global Planet. Change 64(3), 146–157 (2008).ADS 

Google Scholar 
64.Sommer, S. et al. Application of indicator systems for monitoring and assessment of desertification from national to global scales. Land Degrad. Dev. 22(2), 184–197 (2011).
Google Scholar 
65.Vogt, J. V. et al. Monitoring and Assessment of Land Degradation and Desertification: Towards new conceptual and integrated approaches. Land Degrad. Dev. 22(2), 150–165 (2011).
Google Scholar 
66.Scarascia, M. V., Battista, F. D. & Salvati, L. Water resources in Italy: Availability and agricultural uses. Irrig. Drain. J. Int. Commiss. Irrig. Drain. 55(2), 115–127 (2006).
Google Scholar 
67.Wang, H., Yuan, H., Xu, X. & Liu, S. Landscape structure of desertification grassland in source region of Yellow River. J. Appl. Ecol. 17(9), 1665–1670 (2006).
Google Scholar 
68.Wang, J. et al. Spatio-temporal pattern of land degradation from 1990 to 2015 in Mongolia. Environ. Dev. 34, 100497 (2020).
Google Scholar 
69.Briassoulis, H. Governing desertification in Mediterranean Europe: The challenge of environmental policy integration in multi-level governance contexts. Land Degrad. Dev. 22(3), 313–325 (2011).
Google Scholar 
70.Juntti, M. & Wilson, G. A. Conceptualising desertification in Southern Europe: Stakeholder interpretations and multiple policy agendas. Eur. Environ. 15, 228–249 (2005).
Google Scholar 
71.Gisladottir, G. & Stocking, M. Land degradation control and its global environmental benefits. Land Degrad. Dev. 16, 99–112 (2005).
Google Scholar  More

NATURE PODCAST
14 January 2022

Recreating the lost sounds of spring

How citizen science is helping us hear lost soundscapes.

Geoff Marsh

Geoff Marsh

View author publications

You can also search for this author in PubMed
 Google Scholar

Twitter

Facebook

Email

Subscribe
Subscribe

iTunes
Google Podcast
acast
RSS

The researcher resurrecting our declining soundscapes.

Your browser does not support the audio element.

Download MP3

As our environments change, so too do the sounds they make — and this change in soundscape can effect us in a whole host of ways, from our wellbeing to the way we think about conservation. In this Podcast Extra we hear from one researcher, Simon Butler, who is combining citizen science data with technology to recreate soundscapes lost to the past. Butler hopes to better understand how soundscapes change in response to changes in the environment, and use this to look forward to the soundscapes of the future.Nature Communications: Bird population declines and species turnover are changing the acoustic properties of spring soundscapesNever miss an episode: Subscribe to the Nature Podcast on Apple Podcasts, Google Podcasts, Spotify or your favourite podcast app. Head here for the Nature Podcast RSS feed

doi: https://doi.org/10.1038/d41586-022-00023-8

Related Articles

Using sound to explore events of the Universe

Pioneers of sound

Healing sound: the use of ultrasound in drug delivery and other therapeutic applications

Subjects

Ecology

Conservation biology

Latest on:

Ecology

Wind power versus wildlife: root mitigation in evidence
Correspondence 11 JAN 22

Two million species catalogued by 500 experts
Correspondence 11 JAN 22

EU Nature Restoration Law needs ambitious and binding targets
Correspondence 11 JAN 22

Jobs

Postdoctoral Research Assistant

Cancer Research UK Beatson Institute
glasgow, United Kingdom

Postdoctoral Training Fellow in Data Science (ID:53)

Institute of Cancer Research (ICR)
London, United Kingdom

Post Doctoral Research Fellow in Bioinformatics (ID:54)

Institute of Cancer Research (ICR)
London, United Kingdom

Postdoctoral Training Fellow – Systems and Precision Medicine team (ID:87)

Institute of Cancer Research (ICR)
London, United Kingdom More

Comparison of spectral water index methodsFigure 3 shows the results of different spectral water index methods. Through overlay analysis with the original image and detailed visual analysis, it was found that AWEIsh had the best extraction performance and could accurately identify the boundary of the water body. NDWI, MNDWI, and EWI had different degrees of leakage extraction; NDWI and EWI had an evident leakage extraction in the northwest corner of the image, and the water leakage extraction of MNDWI was mainly concentrated in the middle of the image. There was a lot of water body misidentified in WI, especially in the southeast corner of the image.Figure 3Results of water extraction from different spectral water index methods. (a) The original image. The threshold values and extracted water bodies from (b) NDWI, (c) MNDWI, (d) EWI, and (f) AWEIsh methods determined by the OTSU algorithm. (e) The extraction result of WI.Full size imageBased on the visual interpretation of the water boundary, 100 test samples were selected and the confusion matrix32 was calculated to obtain the extraction accuracy of the water body from three aspects: commission error, omission error, and overall accuracy (Table 1). As seen from the table, the overall accuracy of AWEIsh was the highest, attaining a value of 99.14%, with extremely low commission and omission errors. WI had the lowest overall accuracy and the highest commission error, and could not distinguish water bodies and low reflectivity features effectively. The overall accuracies of NDWI, MNDWI, and EWI were similar. Comparing the results of the visual interpretation and quantitative analysis, the rapid extraction model of surface water based on the Google Earth Engine utilizing AWEIsh index was used for assessing the spatiotemporal changes of water bodies.Table 1 Accuracy comparison of different spectral water index methods.Full size tableTime series analysis of typical reservoir areaTo understand the inter-annual variation trend and intra-annual variation of the reservoir area in the dry zone of Sri Lanka, time series analysis was conducted with the Maduru Oya Reservoir as the case study area. The Maduru Oya Reservoir is the second largest reservoir in Sri Lanka, located in the east-central region, which is the main water source for irrigation and drinking, and has a high incidence of chronic kidney disease of unknown aetiology (CKDu). Figure 4 shows the inter-annual and intra-annual variations of Maduru Oya Reservoir area.Figure 4Observed area change in the Maduru Oya Reservoir. (a) Inter-annual variation of the Maduru Oya Reservoir area; (b) Intra-annual variation of the Maduru Oya Reservoir area in 2017.Full size imageFigure 4 shows that the inter-annual fluctuation of Maduru Oya Reservoir area is slight, while the intra-annual fluctuation is significant. From 1988 to 2018, the reservoir area showed an overall increasing trend with slight float; the smallest area was recorded in 1992 (27.43 km2) and the largest area in 2013 (42.97 km2) (Fig. 4a). The rainy season in the dry zone of Sri Lanka occurs from October to February, and the dry season occurs from March to September. In 2017, the maximum area of the Maduru Oya Reservoir was noted in February, and the minimum area was noted in September. The area in February was 2.24 times bigger than that of September, with a difference of 31.58 km2. The maximum area of reservoirs or lakes generally occurs at the end of the wet season (February), and the minimum area occurs at the end of the dry season (September)2, which is consistent with the occurrence of maximum and minimum area in the Maduru Oya Reservoir in 2017(Fig. 4b). The area of the reservoir increased significantly in May during the dry season. According to meteorological data33, there were persistent strong winds and torrential rains in Sri Lanka in May 2017, resulting in an abnormal increase in the reservoir area. Generally, the period in which the area increased was from October to February (rainy season), while March to September (dry season) was the period in which the area decreased regardless of the influence of abnormal weather factors. The intra-annual fluctuation of the reservoir was severe, and there was a risk of drought and flooding at the same time. This observation implied that the seasonal regulation of water resources must be focussed in the future.Analysis of spatiotemporal change of inland lakes and reservoirsTo systematically analyze the spatiotemporal variation characteristics of inland water in Sri Lanka in recent years, and considering the cloud cover of Landsat-5/8 images, 1995, 2005 and 2015 were selected as the study year with an interval of 10 years. The distribution information of surface water in three stages was obtained by running the rapid extraction model of surface water in the Google Earth Engine. According to statistics, the surface water areas of Sri Lanka in 1995, 2005, and 2015 were 1654.18 km2, 1964.86 km2, and 2136.81 km2, respectively. In the past 20 years, the water area of Sri Lanka has increased significantly. To further analyse the spatiotemporal changes of inland lakes and reservoirs, a 5-m buffer data of rivers in 2015 were produced in ArcGIS10.3 software; further, the area corresponding to the river channels were removed from the three images and only the lagoon areas were preserved. Lagoons are ubiquitous in the coastal areas of Sri Lanka, with flood discharge, aquaculture, coastal protection, and other functions34. The results consisting of the extracted lakes, reservoirs, and lagoons are shown in Fig. 5.Figure 5Water extraction results for Sri Lanka in 1995, 2005, and 2015. The administrative boundary data of Sri Lanka comes from the Humanitarian Data Exchange (HDX) open platform (https://data.humdata.org). The maps were generated by geospatial analysis of ArcGIS software (version ArcGIS 10.3; http://www.esri.com/software/arcgis/arcgis-for-desktop).Full size imageThe overall water area of lakes and reservoirs in Sri Lanka showed an increasing trend from 1995 to 2015, and the lagoon area increased over these 20 years (Fig. 5). Because the lagoon does not belong to inland freshwater sensu stricto, the corresponding statistical analysis was not included in the following step. According to statistics, the total area covered of lakes and reservoirs in Sri Lanka were 1020.41 km2, 1270.53 km2, and 1417.68 km2 in 1995, 2005, and 2015 respectively. In the past 20 years, the area of lakes and reservoirs in Sri Lanka has increased by a considerable margin, attaining a value of 397.27 km2. To further analyse the spatiotemporal variation of inland lakes and reservoirs, they were divided into four grades according to their area: I ( More

Field studiesPhotosynthetic carbon fixation was investigated at eight different stations in the Pearl River estuary of the South China Sea (Fig. 1a and Supplementary Table 1), where the phytoplankton assemblages were dominated by diatoms45 during the time of our investigation (June 2015). Samples were collected from 10 to 20 m depths and transferred immediately into 50 mL quartz tubes and sealed to prevent gas exchange. The samples were inoculated with 100 μL of 5 μCi (0.185 MBq) NaH14CO3 solution for 2.15 h. All the incubations were carried out under incident solar radiation, attenuated with neutral density filters to simulate light intensities at the sampling depths, and the temperature was controlled with flow-through surface seawater.After incubation, the cells were filtered onto glass-fiber filters (25 mm, Whatman GF/F, USA) and stored at −20 ° C until measurement, during which the filters were exposed to HCl fumes overnight and dried (20 °C, 6 h) to remove unincorporated NaH14CO3 as CO2. The incorporated radioactivity was measured by liquid scintillation counting (LS 6500, Beckman Coulter, USA), and photosynthetic carbon fixation rates were estimated as previously reported46. Since the measurements were carried out under varying and low light levels similar to in situ levels at depths of 10 and 20 m, we normalized the photosynthetic rates to light intensity (μmol C (μg Chl a)−1 h−1 (μmol photons m−2 s−1)−1) to obtain the light use efficiency of photosynthesis (PLUE). This was done to allow for a meaningful comparison among different stations according to the linear relationship of photosynthetic carbon fixation under low solar irradiance levels46, which lies within the range of sunlight levels used in the present fieldwork ( More

1.Bird, C. S. et al. A global perspective on the trophic geography of sharks. Nat. Ecol. Evol. 2(2), 299–305. https://doi.org/10.1038/s41559-017-0432-z (2018).Article 
PubMed 
PubMed Central 

Google Scholar 
2.Ferretti, F., Worm, B., Britten, G. L., Heithaus, M. R. & Lotze, H. K. Patterns and ecosystem consequences of shark declines in the ocean. Ecol. Lett. 13(8), 1055–1071. https://doi.org/10.1111/j.1461-0248.2010.01489.x (2010).Article 
PubMed 
PubMed Central 

Google Scholar 
3.Hammerschlag, N., Schmitz, O. J., Flecker, A. S., Lafferty, K. D., Sih, A., Atwood, T. B., Gallagher, A. J., Irschick, D. J., Skubel, R., & Cooke, S. J. Ecosystem function and services of aquatic predators in the anthropocene. In Trends in Ecology and Evolution Vol. 34, Issue 4, 369–383. (Elsevier Ltd, 2019). https://doi.org/10.1016/j.tree.2019.01.0054.Heithaus, M. R., Frid, A., Wirsing, A. J. & Worm, B. Predicting ecological consequences of marine top predator declines. Trends Ecol. Evol. 23(4), 202–210. https://doi.org/10.1016/j.tree.2008.01.003 (2008).Article 
PubMed 
PubMed Central 

Google Scholar 
5.Williams, J. J., Papastamatiou, Y. P., Caselle, J. E., Bradley, D. & Jacoby, D. M. P. Mobile marine predators: An understudied source of nutrients to coral reefs in an unfished atoll. Proc. R. Soc. B Biol. Sci. 285(1875), 20172456. https://doi.org/10.1098/rspb.2017.2456 (2018).Article 

Google Scholar 
6.Dulvy, N. K., Simpfendorfer, C. A., Davidson, L. N. K., Fordham, S. V., Bräutigam, A., Sant, G., & Welch, D. J. Challenges and priorities in shark and ray conservation. In Current Biology, Vol. 27, Issue 11, R565–R572. (Cell Press, 2017). https://doi.org/10.1016/j.cub.2017.04.038.7.MacNeil, M. A. et al. Global status and conservation potential of reef sharks. Nature 583(7818), 801–806. https://doi.org/10.1038/s41586-020-2519-y (2020).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
8.MacKeracher, T., Diedrich, A. & Simpfendorfer, C. A. Sharks, rays and marine protected areas: A critical evaluation of current perspectives. Fish Fish. 20(2), 255–267. https://doi.org/10.1111/faf.12337 (2019).Article 

Google Scholar 
9.Albano, P. S. et al. Successful parks for sharks: No-take marine reserve provides conservation benefits to endemic and threatened sharks off South Africa. Biol. Conserv. 261, 109302 (2021).Article 

Google Scholar 
10.Bond, M. E. et al. Reef sharks exhibit site-fidelity and higher relative abundance in marine reserves on the Mesoamerican Barrier reef. PLOS ONE 7(3), e32983. https://doi.org/10.1371/journal.pone.0032983 (2012).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
11.Ruppert, J. L. W. et al. Human activities as a driver of spatial variation in the trophic structure of fish communities on Pacific coral reefs. Glob. Change Biol. 24(1), e67–e79. https://doi.org/10.1111/gcb.13882 (2018).Article 

Google Scholar 
12.Valdivia, A., Cox, C. E. & Bruno, J. F. Predatory fish depletion and recovery potential on Caribbean reefs. Sci. Adv. 3, e1601303 (2017).ADS 
Article 

Google Scholar 
13.Dwyer, R. G. et al. Individual and population benefits of marine reserves for reef sharks. Curr. Biol. 30(3), 480–489. https://doi.org/10.1016/j.cub.2019.12.005 (2020).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
14.R Core Team. R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria. https://www.R-project.org/. (2021).15.Wickham, H, ggplot2: Elegant Graphics for Data Analysis. Springer, New York. ISBN 978-3-319-24277-4 (2016).16.Kahle, D. & Wickham, H. ggmap: spatial visualization with ggplot2. R J. 5(1), 144–161 (2013).Article 

Google Scholar 
17.Desbiens, A. A. et al. Revisiting the paradigm of shark-driven trophic cascades in coral reef ecosystems. Ecology 102(4), e03303. https://doi.org/10.1002/ecy.3303 (2021).Article 
PubMed 
PubMed Central 

Google Scholar 
18.Morrissey, J. E. & Gruber, S. H. Habitat selection by juvenile lemon sharks, Negaprion brevirostris. Environ. Biol. Fishes 38, 311–319 (1993).Article 

Google Scholar 
19.Clementi, G. et al. Anthropogenic pressures on reef-associated sharks in jurisdictions with and without directed shark fishing. Mar. Ecol. Prog. Ser. 661, 175–186. https://doi.org/10.3354/meps13607 (2021).ADS 
Article 

Google Scholar 
20.Juhel, J. B. et al. Isolation and no-entry marine reserves mitigate anthropogenic impacts on grey reef shark behavior. Sci. Rep. 9(1), 1–11. https://doi.org/10.1038/s41598-018-37145-x (2019).CAS 
Article 

Google Scholar 
21.Goetze, J. S. et al. Fish wariness is a more sensitive indicator to changes in fishing pressure than abundance, length or biomass. Ecol. Appl. 27, 1178–1189 (2017).Article 

Google Scholar 
22.Mitchell, J. D. et al. Quantifying shark depredation in a recreational fishery in the Ningaloo Marine Park and Exmouth Gulf, Western Australia. Mar. Ecol. Prog. Ser. 587, 141–157. https://doi.org/10.3354/meps12412 (2018).ADS 
Article 

Google Scholar 
23.Mitchell, J. D. et al. A novel experimental approach to investigate the potential for behavioural change in sharks in the context of depredation. J. Exp. Mar. Biol. Ecol. 530–531, 151440. https://doi.org/10.1016/j.jembe.2020.151440 (2020).Article 

Google Scholar 
24.Speed, C. W., Cappo, M. & Meekan, M. G. Evidence for rapid recovery of shark populations within a coral reef marine protected area. Biol. Cons. 220, 308–319. https://doi.org/10.1016/j.biocon.2018.01.010 (2018).Article 

Google Scholar 
25.Bond, M. E., Albanese, J. V., Heithaus, E. A. B. M. R. & Cerrato, R. D. G. R. Top predators induce habitat shifts in prey within marine protected areas. Oecologia 190(2), 375–385. https://doi.org/10.1007/s00442-019-04421-0 (2019).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
26.Lester, E. K. et al. Relative influence of predators, competitors and seascape heterogeneity on behaviour and abundance of coral reef mesopredators. Oikos 130, 2239–2249. https://doi.org/10.1111/oik.08463 (2021).Article 

Google Scholar 
27.Phenix, L. et al. Evaluating the effects of large marine predators on mobile prey behavior across subtropical reef systems. Ecol. Evol. 9, 13740–13751 (2019).Article 

Google Scholar 
28.Shea, B. D. et al. Effects of exposure to large sharks on the abundance and behavior of mobile prey fishes along a temperate coastal gradient. PLOS ONE 15(3), e0230308. https://doi.org/10.1371/journal.pone.0230308 (2020).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
29.Sherman, C. S., Heupel, M. R., Moore, S. K., Chin, A. & Simpfendorfer, C. A. When sharks are away, rays will play: Effects of top predator removal in coral reef ecosystems. Mar. Ecol. Prog. Ser. 641, 145–157. https://doi.org/10.3354/meps13307 (2020).ADS 
Article 

Google Scholar 
30.Ryan, K. L., Hall, N. G., Lai, E. K., Smallwood, C. B., Tate, A., Taylor, S. M., & Wise, B. S. Statewide Survey of Boat-Based Recreational Fishing in Western Australia 2017/18, 8. Fisheries Research Report No. 297 (2019).31.Cresswell, A. K. et al. Disentangling the response of fishes to recreational fishing over 30 years within a fringing coral reef reserve network. Biol. Cons. 237, 514–524. https://doi.org/10.1016/j.biocon.2019.06.023 (2019).Article 

Google Scholar 
32.Strydom, S. et al. Too hot to handle: Unprecedented seagrass death driven by marine heatwave in a World Heritage Area. Glob. Change Biol. 26(6), 3525–3538. https://doi.org/10.1111/gcb.15065 (2020).ADS 
Article 

Google Scholar 
33.Goetze, J. S., & Fullwood, L. A. F. Fiji’s largest marine reserve benefits reef sharks. In Coral Reefs Vol. 32, Issue 1, 121–125. (Springer, 2013). https://doi.org/10.1007/s00338-012-0970-4.34.Juhel, J. B. et al. Reef accessibility impairs the protection of sharks. J. Appl. Ecol. 55(2), 673–683. https://doi.org/10.1111/1365-2664.13007 (2018).Article 

Google Scholar 
35.Birt, M. J. et al. Isolated reefs support stable fish communities with high abundances of regionally fished species. Ecol. Evol. 11(9), 4701–4718. https://doi.org/10.1002/ece3.7370 (2021).Article 
PubMed 
PubMed Central 

Google Scholar 
36.Fitzpatrick, R. et al. A comparison of the seasonal movements of tiger sharks and green turtles provides insight into their predator-prey relationship. PLOS ONE 7(12), e51927. https://doi.org/10.1371/journal.pone.0051927 (2012).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
37.Mourier, J. et al. Extreme inverted trophic pyramid of reef sharks supported by spawning groupers. Curr. Biol. 26(15), 2011–2016. https://doi.org/10.1016/j.cub.2016.05.058 (2016).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
38.Braccini, M., Molony, B. & Blay, N. Patterns in abundance and size of sharks in northwestern Australia: Cause for optimism. ICES J. Mar. Sci. 77(1), 72–82. https://doi.org/10.1093/icesjms/fsz187 (2020).Article 

Google Scholar 
39.Holmes, T., Rule, M., Bancroft, K., Shedrawi, G., Murray, K., Wilson, S., & Kendrick, A. Ecological Monitoring in the Ningaloo Marine Reserves 2017 (2017).40.Martín, G., Espinoza, M., Heupel, M. & Simpfendorfer, C. A. Estimating marine protected area network benefits for reef sharks. J. Appl. Ecol. 57(10), 1969–1980. https://doi.org/10.1111/1365-2664.13706 (2020).Article 

Google Scholar 
41.Ferreira, L. C. et al. Crossing latitudes-long-distance tracking of an apex predator. PLOS ONE 10(2), e0116916. https://doi.org/10.1371/journal.pone.0116916 (2015).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
42.Priede, I. G., Bagley, P. M., Smith, A., Creasey, S. & Merrett, N. R. Scavenging deep demersal fishes of the Porcupine Seabight, north-east Atlantic: Observations by baited camera, trap and trawl. J. Mar. Biol. Assoc. 74(3), 481–498. https://doi.org/10.1017/S0025315400047615 (1994).Article 

Google Scholar 
43.Stobart, B. et al. Performance of baited underwater video: Does it underestimate abundance at high population densities?. PLOS ONE 10(5), e0127559. https://doi.org/10.1371/journal.pone.0127559 (2015).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
44.Papastamatiou, Y. P., Lowe, C. G., Caselle, J. E. & Friedlander, A. M. Scale-dependent effects of habitat on movements and path structure of reef sharks at a predator-dominated atoll. Ecology 90(4), 996–1008 (2009).Article 

Google Scholar 
45.Rizzari, J. R., Frisch, A. J. & Magnenat, K. A. Diversity, abundance, and distribution of reef sharks on outer-shelf reefs of the Great Barrier Reef Australia. Mar. Biol. 161(12), 2847–2855. https://doi.org/10.1007/s00227-014-2550-3 (2014).Article 

Google Scholar 
46.Speed, C., Field, I., Meekan, M. & Bradshaw, C. Complexities of coastal shark movements and their implications for management. Mar. Ecol. Prog. Ser. 408, 275–293. https://doi.org/10.3354/meps08581 (2010).ADS 
Article 

Google Scholar 
47.Espinoza, M., Cappo, M., Heupel, M. R., Tobin, A. J. & Simpfendorfer, C. A. Quantifying shark distribution patterns and species-habitat associations: Implications of marine park zoning. PLOS ONE 9(9), e106885. https://doi.org/10.1371/journal.pone.0106885 (2014).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
48.Mourier, J., Planes, S. & Buray, N. Trophic interactions at the top of the coral reef food chain. Coral Reefs 32(1), 285. https://doi.org/10.1007/s00338-012-0976-y (2013).ADS 
Article 

Google Scholar 
49.Raoult, V., Broadhurst, M. K., Peddemors, V. M., Williamson, J. E. & Gaston, T. F. Resource use of great hammerhead sharks (Sphyrna mokarran) off eastern Australia. J. Fish Biol. 95(6), 1430–1440. https://doi.org/10.1111/jfb.14160 (2019).Article 
PubMed 
PubMed Central 

Google Scholar 
50.Andrzejaczek, S. et al. Biologging tags reveal links between fine-scale horizontal and vertical movement behaviors in tiger sharks (Galeocerdo cuvier). Front. Mar. Sci. 6(May), 1–13. https://doi.org/10.3389/fmars.2019.00229 (2019).ADS 
Article 

Google Scholar 
51.Andrzejaczek, S. et al. Depth-dependent dive kinematics suggest cost-efficient foraging strategies by tiger sharks. R. Soc. Open Sci. 7(8), 200789. https://doi.org/10.1098/rsos.200789 (2020).ADS 
Article 
PubMed 
PubMed Central 

Google Scholar 
52.Brooks, E. J., Sloman, K. A., Sims, D. W. & Danylchuk, A. J. Validating the use of baited remote underwater video surveys for assessing the diversity, distribution and abundance of sharks in the Bahamas. Endang. Species Res. 13(3), 231–243. https://doi.org/10.3354/esr00331 (2011).Article 

Google Scholar 
53.Santana-Garcon, J. et al. Calibration of pelagic stereo-BRUVs and scientific longline surveys for sampling sharks. Methods Ecol. Evol. 5(8), 824–833. https://doi.org/10.1111/2041-210X.12216 (2014).Article 

Google Scholar 
54.Barnett, A., Abrantes, K. G., Seymour, J. & Fitzpatrick, R. Residency and spatial use by reef sharks of an isolated seamount and its implications for conservation. PLOS ONE 7(5), e36574. https://doi.org/10.1371/journal.pone.0036574 (2012).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
55.Papastamatiou, Y. P. et al. Activity seascapes highlight central place foraging strategies in marine predators that never stop swimming. Mov. Ecol. 6(1), 1–15. https://doi.org/10.1186/s40462-018-0127-3 (2018).Article 

Google Scholar 
56.Vianna, G. M. S., Meekan, M. G., Meeuwig, J. J. & Speed, C. W. Environmental influences on patterns of vertical movement and site fidelity of grey reef sharks (Carcharhinus amblyrhynchos) at aggregation sites. PLOS ONE 8(4), e60331. https://doi.org/10.1371/journal.pone.0060331 (2013).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
57.Lear, K. O., Whitney, N. M., Morris, J. J. & Gleiss, A. C. Temporal niche partitioning as a novel mechanism promoting co-existence of sympatric predators in marine systems. Proc. R. Soc. B: Biol. Sci. 288(1954), 20210816. https://doi.org/10.1098/rspb.2021.0816 (2021).Article 

Google Scholar 
58.Queiroz, N. et al. Global spatial risk assessment of sharks under the footprint of fisheries. Nature 572(7770), 461–466. https://doi.org/10.1038/s41586-019-1444-4 (2019).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
59.Langlois, T. et al. A field and video annotation guide for baited remote underwater stereo-video surveys of demersal fish assemblages. Methods Ecol. Evol. 11(11), 1401–1409. https://doi.org/10.1111/2041-210X.13470 (2020).Article 

Google Scholar 
60.Lin, X. & Zhang, D. Inference in generalized additive mixed models by using smoothing splines. J. R. Stat. Soc. 61(2), 381–400 (1999).MathSciNet 
Article 

Google Scholar 
61.Fisher, R., Wilson, S. K., Sin, T. M., Lee, A. C. & Langlois, T. J. A simple function for full-subsets multiple regression in ecology with R. Ecol. Evol. 8(12), 6104–6113. https://doi.org/10.1002/ece3.4134 (2018).Article 
PubMed 
PubMed Central 

Google Scholar 
62.Mullahy, J. Specification and testing of some modified count data models. J. Econom. 33, 341–365 (1986).MathSciNet 
Article 

Google Scholar 
63.Tweedie, M. An index which distinguishes between some important exponential families. In Statistics: Applications and New Directions: Proceedings of the Indian Statistical Institute Golden Jubelee International Conference Vol. 604 (1984).64.Wood, S. N. Generalized Additive Models: An Introduction with R 2nd edn. (Chapman and Hall/CRC, 2017).Book 

Google Scholar 
65.Burnham, K. P. & Anderson, D. R. Multimodel inference: Understanding AIC and BIC in model selection. Sociol. Methods Res. 33(2), 261–304. https://doi.org/10.1177/0049124104268644 (2004).MathSciNet 
Article 

Google Scholar 
66.Burnham, K. P. & Anderson, D. R. Model Selection and Multimodel Inference; A Practical Information-Theoretic Approach 2nd edn. (Springer, 2002).MATH 

Google Scholar 
67.Ward‐Paige, C. A., Keith, D. M., Worm, B. & Lotze, H. K. Recovery potential and conservation options for elasmobranchs. J. Fish Biol. 80(5), 1844–1869 (2012).68.Graham, F et al. Use of marine protected areas and exclusive economic zones in the subtropical western North Atlantic Ocean by large highly mobile sharks. Divers. Distrib. 22(5), 534–546 (2016).69.Morgan, A., Calich, H., Sulikowski, J. & Hammerschlag, N. Evaluating spatial management options for tiger shark (Galeocerdo cuvier) conservation in US Atlantic Waters. ICES J. Mar. Sci. 77(7–8), 3095–3109 (2020).70.Harvey, E. S. & Shortis, M. R. A system for stereo-video measurement of sub-tidal organisms. Mar. Technol. Soc. J. 29(4), 10–22 (1995).71.R Core Team (2020). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. URL https://www.R-project.org/72.McLean, D. L. et al. Distribution, abundance, diversity and habitat associations of fishes across a bioregion experiencing rapid coastal development. Estuar. Coast Shelf S. 178, 36–47 (2016).73.Althaus, F.et al. A standardised vocabulary for identifying benthic biota and substrata from underwater imagery: the CATAMI classification scheme. PloS one 10(10), e0141039 (2015).74.Wilson, S. K., Graham, N. A. J. & Polunin, N. V. C. Appraisal of visual assessments of habitat complexity and benthic composition on coral reefs. Mar. Biol. 151(3), 1069–1076 (2007).75.Roff, G. et al. The ecological role of sharks on coral reefs. Trends Ecol. Evol. 31(5), 395–407 (2016). More

Environmental isotopic baseline of the Yingpan cemeteryMultiple body tissues of Yingpan Man, including his serial dentine, bone collagen, hair keratin and muscle, were isotopically analyzed. The δ13C measurements show a mix of C3 and C4 dietary inputs with values that range from − 17.3‰ to − 13.7‰, while the δ15N results are remarkably elevated and range from 14.9‰ to 19.9‰. According to the isotopic literature, plant δ15N values above 20‰ are relatively rare in terrestrial environments22,26,27,28. However, the δ15N values of the plants, animals and Yingpan Man are remarkably elevated (up to 23.3‰) compared to isotopic results from other parts of China and Europe (Fig. 3)12. This suggests that the entire ecosystem of this portion of Xinjiang is 15N-enriched. In addition, seeds of Nitraria pamirica having elevated δ15N results (up to 27.9‰) were found in the western region of the Pamir Plateau near the Afghanistan border23. Grains of wheat and barley from Shichengzi (ca. 202 BC–220 AD; ~ 400 km from Yingpan) yielded δ15N values that range from 14.6‰ to 19.8‰ (mean ± SD = 17.2 ± 1.5‰; n = 10) and 13.4‰ to 19.8‰ (mean ± SD = 17.3 ± 1.9‰; n = 12), respectively29. This evidence suggests that the entire Tarim Basin and greater Xinjiang has some of the most elevated terrestrial δ15N values in Eurasia. This is also supported by past isotopic work in China, which demonstrated a correlation between human δ15N values and annual mean precipitation, with individuals from Xinjiang having the highest δ15N values of all regions studied12. The cause of this 15N-enrichment must be at least partially environmental based on past isotopic studies27,30. The extreme aridity of the Yingpan cemetery site, which is located in the Taklamakan Desert, the driest region of China and characterized by little rainfall of 0 to 100 mm/yr and high evapotranspiration rates  > 2500 mm/year (Fig. 1)31, results in a large plant 15N-enrichment by intensive evaporation of 15N-depleted ammonia (NH4+) from the soil. These elevated δ15N results are then translated up the food chain to the domestic animals and humans. This is compelling evidence that Yingpan man was born and raised in a 15N-enriched environment that was extremely arid, and that he consumed wheat/barley, broomcorn millet and grape that were grown locally.Breastfeeding, weaning and childhood dietary patterns of Yingpan ManIsotopic analysis of dentine serial sections in human teeth permit an investigation of individual dietary patterns over the time period when the teeth were developing24. For Yingpan Man, his first molar (M1), represents the period of his life from birth to approximately 10 years old32. The first five 1 mm dentine sections from the M1 crown, corresponding to the first ~ 2.2 years of his life, show a steady decrease in δ15N of 2.5‰ (Fig. 4). This is evidence that Yingpan Man was fully weaned off breastmilk at or before the age of ~ 2.2 years old (Fig. 4)33,34, and is similar to findings of the Zhou Dynasty (1122–771 BC) sites of Boyangcheng35, Xiyasi and Changxinyuan36.Archaeological evidence for infant feeding practices is rarely preserved in Xinjiang. However, two extraordinary feeding vessels (one made of goat breast skin, the other made from an ox horn) were previously discovered with the burial of a 10-month old infant mummy in the southern Tarim Basin4,7. This infant dates to ~ 1000 BC and was found at the Zhagunluke site which is ~ 460 km from the Yingpan cemetery (Fig. 1). In addition, historical documents also provide supporting evidence that children were fed with animal milk in ancient Xinjiang. For example, the Kharosthī scripts (written in a version of the north Indian Prakrit language) are a collection of contracts, letters and other documents (e.g. wood tablets, leather, silk, paper, etc.) that detail life, events, trading, taxes and agricultural practices during the third to fifth century AD of the Shanshan Kingdom (~ 250 km from Yingpan, see Fig. 6)37. It is recorded in these documents that a “milk fee” of cattle or camel was paid to the adoptee when children were adopted in ancient Shanshan, and this exchange was not only protected by the law but also needed to be ensured by witnesses37,38. Thus, there may have been a common tradition of feeding infants with milk from domestic animals (e.g. goat, cattle, camel, etc.) during the weaning process in ancient Xinjiang.Figure 6Yingpan Man’s grave goods mirroring the cultural contact across the Eurasian continent (~ 300 AD). (a) Decorative patterns of vase and flowers on the sideboards of Yingpan Man’s coffin. (b) Decorative patterns of nude putties, goats, bulls and pomegranate trees on Yingpan Man’s woolen caftan. (c) Yingpan Man’s white hemp mask with golden diadem (possibly related to the death masks from Tashtyk culture). (d) Woolen trousers with decorative patterns of double quatrefoil florals surrounded by lozenges made up of solid circles and flowers. (e) Miniature funerary summer coat. (f) Miniature funerary winter robe. (g) Embroidered armband. (h) Yellow painted wood coffin with decorative patterns. (i) A broken brocade that was decorated with Chinese characters of “Shou” and “You”. (j) Silk fish. (k) Crowing cockerel pillow. (l) Perfume sachet. (m) A tufted carpet with decorative patterns of a male lion. (Map was generated using GMT 5.2.1. The original pictures of artifacts were previously published14 and provided by Wenying Li. Pictures were modified using Adobe Photoshop CC 2015 V.1.2. The final layout was created in Adobe Illustrator CC 2019 V.23.1.1).Full size imageIn contrast, the δ13C results of the first four dentine sections show little variability, which indicates that Yingpan Man consumed other foods/liquids from an early age and was not exclusively breastfed for a significant amount of time after birth (Fig. 4; Fig. S4)34. Between 4 and 9 mm, there is an increase in the δ13C values of 1.7‰, which corresponds to a late weaning and childhood diet from approximately 1.8 to 4.2 years old, where more millet was consumed (Fig. 4; Fig. S4). This millet possibly took the form of a gruel as there is little evidence that Yingpan Man consumed increased amounts of animal protein during this period. Similar isotopic patterns in dentine representing individuals that consumed increasing amounts of millet during weaning and early childhood were found at the Late Neolithic (4500 BP) Gaoshan site in Sichuan Province39. This suggests that millet may have had a long history of being used as a weaning or childhood food in China. From 9 to 14 mm, the dentine serial sections show a decline in δ13C values of 1.4‰ which indicates a period of increasing C3 foods in the diet from approximately 4.2 to 6.6 years old (Fig. 4; Fig. S4). Then from 14 to 21 mm the δ13C values increase again by 1.3‰, evidence of a change back to millet consumption from approximately 6.6 to 10 years old. Thus, over the first 10 years of life, Yingpan Man had frequent dietary shifts between C3-based and C4-based foods (Fig. 4; Fig. S4).Seasonal diet of Yingpan ManIn agreement with his dentine serial sections, the δ13C results of Yingpan Man’s hair display fluctuations between C3-based (e.g. wheat, barley) and C4-based (millet) foods (Fig. 5a). These hair δ13C results appear to follow a periodic trend, and suggest that his diet changed during different seasons of the year. In contrast, apart from the last 6–7 months of life, Yingpan Man’s hair δ15N values show little variation (Fig. 5b). However, a strong correlation was found between the variability of both the hair δ13C and δ15N values (Sperman’s r = − 0.534; sig = 0; n = 46; Fig. S5a), but not between the δ13C and δ34S values (Sperman’s r = 0.100; sig = 0.539; n = 40; Fig. 5; Fig. S5b). The lack of significant changes in the δ15N values is evidence that Yingpan Man’s protein consumption was relatively constant for approximately the last 4 years before death. This constant protein consumption, but variable intake of C3 and C4 plants, could suggest that in addition to directly consuming different amounts of wheat/barley and millet, Yingpan Man may have actually consumed domestic animals (goats, sheep, cattle) that were foddered on these crops at different times of the year. This possibility could explain the periodic variability in the δ13C values as well as the lack of variability in the δ15N values. In addition, the δ34S hair results also show little variation in the serial sections (~ 1.5‰). Yingpan Man’s δ34S results have a terrestrial range between 10.4‰ and 11.9‰40,41,42 and are similar to past δ34S values from the Proto-Shang site of Nancheng in Hebei Province43. While few δ34S studies have been reported for China for comparison, this lack of isotopic variability in sulfur suggests that Yingpan Man was likely not a Silk Road traveler, but stayed close to the local area during the last years of his life42. Future research involving strontium analysis on Yingpan Man’s hair serial sections and teeth will hopefully support or refute these findings44.Regional evidence in support of seasonal dietary changes comes from apatite δ13C and δ18O results of serial sections of tooth enamel of domesticates from the pastoral sites of Dali, Begash and Tasbas in Kazakhstan (~ 750 km from Yingpan)45. Specifically, data from early and middle Bronze Age sheep, goats and cattle displayed periodic patterns in their δ13Capa and δ18Oapa results that reflect the consumption of more millet during winter and more C3 plants during summer months, based on environmental inputs of body water45. A combination of radiocarbon dates and the application of the dietary mixing model (MixSIAR) identified early sheep and goats from this region (e.g. Dali, 2705–2545 cal BC) to be winter foddered with up to 44–50% of millet intake45. During later periods, this reliance on millet fodder during the winter months increased with some goats from Begash having 50–60% millet in their diets. Further, the relative contribution of millet to the diet of sheep and goats from the early phase (2345–2080 cal BC) of Begash reached up to 68–74% of the whole diet during winter months, especially from November to December45. This key study indicates that there is a long history and precedent for seasonal C3 and C4 feeding of domestic animals in Central Asia, and supports the isotopic findings of a seasonal diet in the hair of Yingpan Man.In addition, isotopic analysis of sequential hair samples from mummies recovered from the Oglakhty cemetery in the Minusinsk Basin of southern Siberia, Russia (~ 900 km north of Yingpan) also show seasonal dietary variation with millet and fish consumed during the late summer and autumn46. These mummies of the Tashtyk culture date to the same period as Yingpan Man (third to fourth centuries AD), and interestingly were also buried with white painted gypsum funerary masks that are similar to the one that Yingpan Man was wearing46,47. This unique burial tradition could suggest links such as trade or an association between Yingpan Man and the Tashtyk culture and additional research is necessary to explore this possibility in more detail (Fig. 6).Historical documents such as the Kharosthī scripts provide additional valuable information about seasonal diets and the foods consumed by the inhabitants of ancient Xinjiang37. These texts describe how the people of the Tarim Basin cultivated mainly wheat, millet and barley as their main cereal crops and that grapes were carefully managed for the production of wine37. Autumn (around the 10th month of the year), was mentioned as the time for harvesting crops and trading crops as well as wine and animals, and paying debts37. Thus, autumn would have been the time of the year with the most abundant amount of food resources, especially C3 foods like fruit, vegetables and wine. This is still true today, as the harvesting of agricultural products in the Tarim Basin mainly occurs in the month of October48. In contrast, during the winter months food resources would have been scarce with only non-perishable crops like millets and wheat available to guarantee the food supply “in the harsh winters of Inner Asia”45. Thus, more C4 foods were likely consumed during the winter and spring months while more C3 foods were consumed during the summer and autumn months.If this information is applied to Yingpan Man’s hair δ13C results, it would suggest that the δ13C values decreased during the summer and autumn months (June to October) and likely reach a nadir during the middle of autumn or October. Therefore, Yingpan Man’s hair sections which are 8 to 10 cm, 18 to 25 cm and 36 to 43 cm from his scalp reflect the period of middle autumn (between August to October) while the hair sections which are 4 to 6 cm, 13 to 15 cm and 29 to 33 cm from his scalp represent late winter (December to February). This would suggest that Yingpan Man died ~ 4 months after the last drop in his δ13C hair values or in spring, possibly March or April (Fig. 5). In addition, the clothes in which Yingpan Man was buried, as well as his associated wardrobe, also provide information about the timing of his death14,15. The miniature robe placed on his chest was designed for winter, as it was long, double layered and the interior was lined with sheep’s wool (Fig. 2f). Whereas the miniature coat placed near his wrist was designed for summer as it was shorter and made only of a single layer of silk (Fig. 2p). Yingpan Man’s caftan, in which he was buried, was double layered with the outer layer made of wool and the inner lining made of silk, and this was likely designed either for the spring or autumn (Fig. 2k)14. Thus, Yingpan Man’s burial clothes combined with the isotopic and historical evidence indicate that Yingpan Man died sometime during the spring months14,15.To better illustrate Yingpan Man’s seasonal dietary variation, two hair sections respectively representing the highest (Sample A: − 14.0 ± 0.1‰, 15.2 ± 0.1‰) and lowest δ13C values (Sample B: − 17.3 ± 0.2‰, 15.8 ± 0.1‰) were analyzed by a Bayesian mixing model with the application of FRUITS (Food Reconstruction Using Isotopic Transferred signals)49,50. As displayed in Fig. S6, isotopic data for millet, wheat/barley, grape and sheep/goat from the Yingpan cemetery were incorporated into the mixing model as feasible dietary sources for Yingpan Man in both scenarios (Sample A and Sample B). The isotopic fractionation between human hair and diet is corrected with an offset of 4.0 ± 0.5‰ for δ13C values51,52,53 and 4.5 ± 0.5‰ for δ15N values51,54. The relative contribution of different macronutrients is defined in the mixing model according to published records with 74 ± 4% of the carbon originating from protein and 26 ± 4% originating from carbohydrates and lipids and all of the nitrogen originating from protein55,56. The high hair δ13C value of − 14.0‰ (Sample A) represents a heavy reliance on millet consumption (39–67%; median = 53%) and a low amount of dietary input from wheat/barley (0–51%; median = 12%), grape (1–46%; median = 25%) and sheep/goat (0–30%; median = 4%), which likely reflects Yingpan Man’s diet during the winter and spring months (Fig. S6 and Table S13). In contrast, the low hair δ13C value of − 17.3‰ (Sample B) likely represents a summer/autumn diet with a decline in C4 foods, like millet (8–44%; median = 28%), as well as increased reliance on C3 foods like wheat/barley (1–71%; median = 21%) and sheep/goat (0–81%; median = 9%) (Fig. S6 and Table S13). However, the importance of sheep/goat and wheat/barley is likely under-estimated here given their small sample size and highly elevated δ15N values. Nonetheless, the FRUITS mixing model indicates that the varying consumption of plant foods, especially millet, is clearly responsible for the δ13C shifts of Yingpan Man’s diet.The last months of Yingpan ManThe hair δ13C results provide an estimate for the time of year Yingpan Man likely died. In addition, the hair δ15N and δ34S values can provide evidence for how Yingpan Man may have died. The last ~ 6 cm (closest to the scalp) of Yingpan Man’s hair show a general rise in δ15N by ~ 1‰ (Fig. 5b). This pattern is uncharacteristic compared to the other hair δ15N results that display little variation. This unique 15N-enrichment could represent a period of catabolic wasting due to the recycling of tissue proteins as a result of a prolonged illness25,57,58. Additional support for some sort of disease or period of illness comes from the fact that there is little change in the δ13C values but a slight decrease in the δ34S values during the last ~ 6 cm of Yingpan Man’s hair. Tissue catabolism is known to cause an increase in δ15N but little change in δ13C values in human hair25. Further, δ34S results are known to decrease in the red blood cells and serum (by ~ 1.5‰) of patients suffering from liver cancer59. As the last ~ 6 cm of hair displayed a slight decrease in δ34S by ~ 1‰, with the largest decline during the last month of life, this evidence in conjunction with the δ13C and δ15N values suggests that Yingpan Man did not die suddenly but likely suffered some type of debilitating sickness over the last months of his life before he succumbed. However, as most wasting diseases or illnesses leave no traces on human skeletons, it is difficult to define the specific disease that caused Yingpan Man’s death, and a detailed paleo-pathological study of Yingpan Man is needed in the future.The grave goods buried with Yingpan Man also suggest he may have suffered a compromised health status before death. In particular, a piece of tattered yellow brocade decorated with brown and blue images of vines, animals, birds, flowers, as well as the Chinese characters of “Shou” and “You” was found placed at a prominent position on the right side of Yingpan Man’s head (Fig. 7)14. According to Chinese historical literature sources, “You” means “blessing” or “blessed”. While, “Shou” means “long live” or “healthy” (Fig. 7)60. Similar grave goods with Chinese characters e.g. “Yan Nian Yi Shou Da Yi Zi Sun” (meaning “live longer and benefit the descendants”), “Chang Le Ming Guang Cheng Fu Shou You” (meaning “always be happy, be bright, be lucky and be blessed”), “Yong Chang” (meaning “always be prosperous”), were also frequently unearthed from contemporary and later cemeteries of the Tarim Basin, especially the nearby sites of the Loulan culture (third century BC to 448 AD, renamed as Shanshan in 77 BC, ~ 250 km from Yingpan), as grave goods that carry good wishes for the dead15. However, Yingpan Man’s brocade shows significant traces of wear or “rubbing” (Fig. 7). This is interesting as this brocade was not complete or new but well-worn to the point of being tattered and frayed, yet it was still placed at a very important position in Yingpan Man burial—just beside his head14. This suggests that this brocade may have been an important “lucky” health charm for Yingpan Man that was frequently used either by himself or by those who cared for him before burial. Thus, both isotopic and archaeological evidence suggest Yingpan Man suffered some type of illness during the last ~ 6 months of his life, likely in winter, and that he succumbed to this illness in the following spring.Figure 7Photo and sketch showing details of the embroidered brocade that was found with Yingpan Man burial. (The original picture and sketch were previously published14 and provided by Wenying Li. Pictures were modified using Adobe Photoshop CC 2015 V.1.2. The final layout was created in Adobe Illustrator CC 2019 V.23.1.1).Full size imageWho was Yingpan Man?The social identity and status of Yingpan Man is enigmatic given the various cultural components of his grave goods and since the collection of his funerary objects are unique compared to all of the other burials found in the same cemetery14. This has created much controversy and debate about the social identity of Yingpan Man7,14,15,16. Thus, “Who was Yingpan Man?” and “Why was he buried here, in a normal unmarked grave with such lavish and exotic grave goods?” is an active topic of debate.The physical anthropology of Yingpan Man was investigated but not formally published. According to which, the metric and nonmetric index of traits of his skeletal remains suggest a mix of both European and Mongolian features (Dong Wei, personal communication). However, facial reconstruction conducted by Dong Wei (personal communication) indicates that Yingpan Man’s facial structure is more characteristic of the features from Western Eurasia. This is consistent with the image of the face painted on his death mask and the fact that he wears a golden diadem across his forehead, which is more traditionally associated with Greece4,7. However, Yingpan Man’s white hemp mask is similar in style to the white painted gypsum masks of the Tashtyk culture from the Minusinsk Basin of Russia46. In addition, Yingpan Man and the mummies of the Oglakhty cemetery of Tashtyk date to nearly the same period (third to fourth centuries AD), and the polychrome silk cloth from the Tarim Basin has also been discovered in Oglakhty46. This evidence, while circumstantial, could suggest some form of association or that links with the Tashtyk may have taken place during the lifetime of Yingpan Man, possibly through trade or familial relationships (Fig. 6).Other components of Yingpan Man’s burial provide important clues about his social identity and cultural affiliations in life61. The styles and types of grave goods of Yingpan Man display an unique mix of both Eastern and Western cultures and traditions that were likely common to inhabitants of Silk Road trading towns in Xinjiang during the third to fourth centuries AD (Fig. 6)14,15,16. However, there appears to be a strong eastern component in some of his funerary arrangements. For example, the burial practices associated with Yingpan Man: covering his face, filling his nose, burying his body fully clothed, covering it with a silken burial shroud, as well as using miniature funeral objects as grave goods are in accordance with the suggested burial rites of “Yan” (meaning “covering”), “Zhen” (meaning “filling”), “She Min Mu” (meaning “covering the eyes”), “Qin” (meaning “quilt” or “burial shroud”) mentioned in the Confucian literature of Yili (meaning “Rites”, formed during the Zhou Dynasties (1046 BC to 256 BC))15,62. In addition, the styles and designs of some of his grave goods are indicative of Chinese spiritual beliefs. For instance, the diamond and circle-shaped designs of decorative patterns on the cover and sideboards of Yingpan Man’s coffin are argued by some scholars to be the traditional “Lianbi” pattern (meaning “linked jades”) which symbolizes the jade burial suits (“Jinlü Yuyi”) that were popular in Han Dynasty burials of high status nobles in central and southern China (e.g. Nanyuewang tomb, Mancheng tomb)15. In particular, as jades are believed to be the ideal material for embalming in ancient China63, the “Lianbi” pattern on Yingpan Man’s coffin carries the symbolic meaning of preserving his body forever so that his soul and spirt will reach heaven (Fig. 6; Fig. S1)15.The brocade “health charm” found to the right side of Yingpan Man also shows a clear affiliation with China, as it was decorated with the Chinese characters of “Shou” and “You”. This is important as some historical documents suggest that Kharosthī was the common language in the Tarim Basin at this time37,38, but Yingpan Man was likely more accustomed to Chinese spiritual beliefs according to this brocade (Figs. 6, 7). Moreover, the crowing cockerel pillow is also interesting as it was recorded in the Han Dynasty literature of Lunheng (meaning “On Balance”, compiled by Chong Wang during the Eastern Han Dynasty in 88 AD), that the dead “turns into ghosts” and “hurt the living ones”64, and that a crowing cockerel was believed to be able to expel evil spirts and ghosts65. Thus, a crowing cockerel pillow was frequently used in funerary practices in China since the Han Dynasty66, and it is still a common grave good today in some areas of modern China (e.g. Shandong and Guizhou Provinces)65.The eight decorative pearls attached to the crowing cockerel pillow are also important status markers regarding the identity of Yingpan Man (Fig. 8). In Xinjiang, pearls would have been long distance imported products from the coastal regions which are  > 3000 km from the Yingpan cemetery. Ancient China, Egypt, Persia, Greece and India were known to have produced and prized pearls67. However, according to Chinese historic literature sources: Shangshu (meaning “the Book of Documents”; written by pre-Qin philosophers during the Zhou Dynasty ~ 1000 BC)68, Hanshu (meaning “the Book of Han”; written by Gu Ban during the Eastern Han Dynasty in 105 AD)69 and Hou Hanshu (meaning “the Book of Later Han”; written by Ye Fan during the Southern Dynasties from 432 to 445 AD)70, pearls were known as royal tributes in China since the Xia Dynasty (~ 2000 BC), and were mainly produced in the coastal cities of southern China, such as Panyu (modern Guangzhou City in Guangdong Province), Hepu (modern Hainan Province), Zhuya (modern Guangxi and Guangdong Province), etc.67 These coastal regions of China were connected to the Tarim Basin area via the Silk Road trading routes and were possibly the source for the pearls found on Yingpan Man’s pillow (Fig. 6). In particular, a silk pillow fully covered with unpolished natural pearls (placed beneath the tomb owner’s head, weighing ~ 470 g) and a lacquer box of pearls (weighing 4,117 g) were unearthed from the tomb of the King of the Nanyue Kingdom (a vassal state of Han) in modern Guangzhou71. This demonstrates a similar preference for pearls as decorations on pillows or as grave goods for high status individuals.Figure 8Crowing cockerel pillow. (a–d) Photo and sketch of the crowing cockerel pillow from Yingpan Man burial. (e–h) Decorative patterns of the four mythical beasts on Han Dynasty eaves tiles. (The original pictures and sketches were previously published14 and provided by Wenying Li. Pictures were modified using Adobe Photoshop CC 2015 V.1.2. The final layout was created in Adobe Illustrator CC 2019 V.23.1.1).Full size imageMoreover, the decorative patterns found on Yingpan Man’s crowing cockerel pillow are remarkable. These include the images of: a monkey-shaped face, a griffin-shaped beast and a net-shaped pattern as well as the images of the mythical beasts of the “blue dragon”, “red sparrow” and “white tiger” (Fig. 8d)17,47,61,72. Together with another mythical beast known as the “black turtle-snake”, these four mythical beasts were believed to be the guardians of the four cardinal directions and also represents four different colors and elements according to traditional Chinese cosmology (Fig. 8e–h). The dragon guarding the East and representing the color blue and the element of wood, the sparrow guarding the South and representing the color red as well as the element of fire, the tiger guarding the West and representing the color white and the element of metal, and the turtle-snake guarding the North and representing the color black and the element of water62,73. Thus, the images of these four mythical beasts were frequently integrated into the designs of cities, buildings, cemeteries and objects of ancient China for their symbolic function of protection63. In this context, the net-shaped pattern on Yingpan Man’s crowing cockerel pillow is likely a symbol of the four cardinal directions, while the monkey shaped face possibly symbolizes a human, that is to be protected. Notably, here, in the design of Yingpan Man’s crowing cockerel pillow, the image of the northern guardian of the turtle-snake is replaced by the Greek mythical beast known as Griffin (Fig. 8d)16,47,74. This change in imagery and symbolism could suggest that contact between the East and the West was via a northern route and that this was associated with some type of Greek influence or the Griffin. According to historical literature sources, a large number of Greeks migrated into Central Asia with Alexander the Great’s eastward expedition3,5,6. Hellenistic kingdoms such as Bactria (a.k.a. Daxia in the Chinese literature, located in the area of Pamir Plateau in modern Afghanistan, ~ 1900 km from Yingpan) subsequently became centers of Greek influence in this area and exported elements of Hellenistic culture to surrounding kingdoms that were located along the main routes of the Silk Road, such as Dayuan (located in the Fergana Valley, ~ 1300 km from Yingpan)2,6,75. Some ancient cities in the Tarim Basin area, such as Yingpan, Maideke and Yuansha, are also argued to have been influenced by Hellenistic elements as their city walls were circular in shape and this was clearly a unique architecture compared to the square-shaped traditional Chinese city walls63,76. It is highly probably that these Hellenistic elements were transported to the ancient city of Yingpan by the Silk Road trading routes which followed a north–south direction in this region, and therefore would support that Greek elements would be associated with the northern direction in this region of Xinjiang. In conclusion, Yingpan Man’s crowing cockerel pillow displays an unprecedented combination of Chinese spiritual beliefs incorporated with western motifs, and symbolizes the intertwining of Eastern and Western cultures in this part of Central Asia (Figs. 6, 8).Even more significant uses of western components in Yingpan Man’s burial is also visualized in other grave goods, e.g. his: carpet, coffin, caftan. The decorative pattern of the male lion on Yingpan Man’s tufted woolen carpet is clearly an imported element as the lion was originally from Africa, Southern Europe, West Asia and India. Lions were not introduced into China until 87 AD, being sent to Emperor Zhang of the Eastern Han Dynasty as a gift from Pacorus II, the King of the Anxi Empire (a.k.a. Parthian Empire (247 BC to 224 AD, replaced by the Persian in 226 AD, located in modern Iran))70. In particular, this exquisite lion-decorated carpet suggests a very high social status for Yingpan Man as lions were deified in ancient China since the Han Dynasty and were used as symbols of power, authority and royalty as it was in other areas of the Eurasia (Fig. 2a)15,61. Moreover, it is also argued by some scholars that the motif of the lion on this carpet is from the Buddhist art of India, as lions are frequently mentioned in Buddhist stories and a Buddhist temple was also found at the site of Yingpan (Fig. 6)17,77.In addition, though the main decorative patterns on Yingpan Man’s coffin are the “Lianbi” pattern, images of vines, leaves, pomegranate flowers and vases were also depicted inside of the diamond patterns (Fig. S1)15. Among which, the image of pomegranate flowers is clearly an imported element from the West as pomegranates were originally domesticated in the middle East ~ 5000 years ago and were not introduced into China until the Han to Jin Dynasties, or with the flow of goods along the Silk Road78. In particular, the pomegranate was used as a symbol of health, fertility and rebirth as mentioned in many ancient cultures, especially in Greek and Turkish myths78. Thus, this is particularly interesting as the image of pomegranate appears not only on Yingpan Man’s coffin as flowers on the headboard and footboard (Fig. 6; Fig. S1), but also on his caftan as decorative patterns of fruited-trees (Figs. 6, 9). Specifically, the decorative patterns on Yingpan Man’s woolen caftan consists of six sets of nude puttis and animals with fruited pomegranate trees standing in between (Fig. 9)14. In particular, each of these six sets of images is composed of a symmetrical pair of confronting muscular puttis or animals (goats or bulls) that are either leaning away from or toward each other14. The nude puttis are holding either a spear, sword or a shield with capes swirling from their shoulders, while the animals of goats and bulls are in the pose of jumping and the bulls have laurel wreaths around their waists (Fig. 9)14. A similar design of stance and composition of figures was also discovered in a mosaic floor from the Villa of Good Fortune at Olynthos, Greece (paired female figures with weapons, fourth century BC)79 as well as another mosaic floor from Pella, the Macedonian capital (two nude youths with capes flying on their shoulder and weapons held in their hand, about to attack an animal in between of them, 325 to 300 BC)79. Thus, it is suggested that this caftan was the work of a weaver familiar with both Western and Eastern motifs as the character and poses of the nude puttis are clearly Western in style and appearance, the fruited pomegranate tree is believed to be a Persian motif, while the paired facing animals of goats and bulls are similar to the animal art of Central Asia8. However, this caftan was completed using the technique of double-weaving, the lining and belt were made of silk14, while the design of the side slits are indicative of a localized adaptation for horse-riding15. Moreover, analysis using high performance liquid chromatography (HPLC) on the red and yellow threads of this caftan suggest that they were dyed in a local workshop with indigenous materials of Rubia Tinctorum and Populus Pruinosa Schrenk, respectively80,81. In conclusion, the design and style of Yingpan Man’s caftan is unprecedented, and it is a masterpiece that combines both Greek/Roman, Persian, Central Asian, Chinese and local elements (Fig. 6).Figure 9Photo and sketch of the decorative patterns on Yingpan Man’s caftan. (The original pictures and sketches were previously published14 and provided by Wenying Li. Pictures were modified using Adobe Photoshop CC 2015 V.1.2. The final layout was created in Adobe Illustrator CC 2019 V.23.1.1).Full size imageThe opulence and fine quality of the objects buried with Yingpan Man indicate that he must have had a high social status before death14,47. Given the importance of the town of Yingpan as a trading center on the Silk Road, the excavators who discovered Yingpan Man suggested that he was a wealthy merchant from the West14. Others have proposed that Yingpan Man might have been a Sogdian merchant since the Sogdians (an Iranian-speaking people whose homeland lay near Samarkand in what is now Uzbekistan) were the richest traders along the route7. However, given the relatively young age of Yingpan Man before his death (~ 30 years old), it is unlikely that he amassed all his fortune and high social status only through trade, and possibly by inheritance or military feats. The government of the Jin Dynasty established the administrative organization of “Xiyu Zhangshi Fu” (meaning “Chief Governor of the Western Regions”) in ancient Xinjiang, and the capital city of “Xiyu Zhangshi Fu” was located nearby Lop Nur and is very close to the ancient city of Yingpan (~ 185 km away)15,16. In addition, comparison of Yingpan Man’s burial to other contemporary burials from Gansu also suggest that Yingpan Man was possibly a military official from the government of Central China16. More supporting evidence comes from the embroidered armband that was buried with Yingpan Man as colorful armbands were suggested to be used by soldiers for the protection from evil forces in ancient Xinjiang (Fig. 6)15. An additional explanation is that Yingpan Man was a noble or even a king of the nearby state named Shan (a.k.a. Moshan)76, and the ancient city of Yingpan was suggested to be the capital city of this state15. An alternative explanation is that Yingpan Man belonged to a local noble family who were displaced from Bactria to the southern Tarim Basin after civil strife in the Kushan Empire at the end of the second century AD, given the popularity of Kushan arts in this area during the Han to Jin Dynasties82,83. The isotopic evidence presented here, in particular the hair δ34S results, add additional information to Yingpan Man’s identity. The lack of δ34S isotopic variability (10.4‰ to 11.9‰) over the last ~ 3–4 years of life indicates that Yingpan Man was not a Silk Road traveler or merchant, at least during this period of his life. Thus, Yingpan Man appears to have been a local, possibly a governmental official or royal to this region of the Tarim Basin, perhaps from the nearby state of Shan. This might suggest why he was buried in the Yingpan cemetery as it was purported to be capital of this ancient state. More

1.Sicard, A., Michalakis, Y., Gutiérrez, S. & Blanc, S. The strange lifestyle of multipartite viruses. PLoS Pathog. 12, 1–19 (2016).
Google Scholar 
2.Lucía-sanz, A. & Manrubia, S. Multipartite viruses: Adaptive trick or evolutionary treat ?. NPJ Syst. Biol. Appl. 34, 1–11 (2017).
Google Scholar 
3.Rojas, M. R., Hagen, C., Lucas, W. J. & Gilbertson, R. L. Exploiting chinks in the plant’s armor: Evolution and emergence of geminiviruses. Annu. Rev. Phytopathol. 43, 361–394 (2005).CAS 
PubMed 

Google Scholar 
4.Zerbini, F. M. et al. ICTV virus taxonomy profile: Geminiviridae. J. Gen. Virol. 98, 131–133 (2017).CAS 
PubMed 
PubMed Central 

Google Scholar 
5.Varsani, A. et al. Establishment of three new genera in the family Geminiviridae: Becurtovirus, Eragrovirus and Turncurtovirus. Arch. Virol. 159, 1873–1882 (2014).CAS 
PubMed 

Google Scholar 
6.Zhou, X. Advances in understanding begomovirus satellites. Annu. Rev. Phytopathol. 51, 357–381 (2013).CAS 
PubMed 

Google Scholar 
7.Lozano, G. et al. Characterization of non-coding DNA satellites associated with sweepoviruses (Genus Begomovirus, Geminiviridae)—Definition of a distinct class of Begomovirus-associated satellites. Front. Microbiol. 7, 1–13 (2016).
Google Scholar 
8.Fondong, V. N. Geminivirus protein structure and function. Mol. Plant Pathol. 14, 635–649 (2013).CAS 
PubMed 
PubMed Central 

Google Scholar 
9.Shafiq, M., Asad, S., Zafar, Y., Briddon, R. W. & Mansoor, S. Pepper leaf curl Lahore virus requires the DNA B component of tomato leaf curl New Delhi virus to cause leaf curl symptoms. Virol. J. 7, 367 (2010).CAS 
PubMed 
PubMed Central 

Google Scholar 
10.Hanley-Bowdoin, L. et al. Geminiviruses : models for plant DNA replication, transcription and cell cycle regulation. Crit. Rev. Plant Sci. 18, 71–106 (1999).CAS 

Google Scholar 
11.De Bruyn, A. et al. East African cassava mosaic-like viruses from Africa to Indian ocean islands: molecular diversity, evolutionary history and geographical dissemination of a bipartite begomovirus. BMC Evol. Biol. 12, 1–18 (2012).
Google Scholar 
12.Navas-Castillo, J., Fiallo-Olivé, E. & Sanchez-Campos, S. Emerging virus diseases transmitted by whiteflies. Annu. Rev. Phytopathol. 49, 219–248 (2011).CAS 
PubMed 

Google Scholar 
13.Rey, M. E. C. et al. Diversity of dicotyledenous-infecting geminiviruses and their associated DNA molecules in southern Africa, including the South-west Indian ocean islands. Viruses 4, 1753–1791 (2012).CAS 
PubMed 
PubMed Central 

Google Scholar 
14.Brown, J. K. et al. Revision of Begomovirus taxonomy based on pairwise sequence comparisons. Arch. Virol. 160, 1593–1619 (2015).CAS 
PubMed 

Google Scholar 
15.Ouattara, A. et al. Diversity, distribution and prevalence of vegetable-infecting geminiviruses in Burkina Faso. Plant Pathol. 69, 379–392 (2019).
Google Scholar 
16.Tiendrébéogo, F. et al. Characterization of pepper yellow vein Mali virus in Capsicum sp. Burkina Faso. Plant Pathol. J. 7, 155–161 (2008).
Google Scholar 
17.Zhou, Y. C. et al. Evidence of local evolution of tomato-infecting begomovirus species in West Africa: Characterization of tomato leaf curl Mali virus and tomato yellow leaf crumple virus from Mali. Arch. Virol. 153, 693–706 (2008).CAS 
PubMed 

Google Scholar 
18.Hamilton, W. D. O., Bisaro, D. M., Coutts, R. H. A. & Buck, K. W. Demonstration of the bipartite nature of the genome of a single-stranded DNA plant virus by infection with the doned DNA component. Nucleic Acids Res. 11, 7387–7396 (1983).CAS 
PubMed 
PubMed Central 

Google Scholar 
19.Padidam, R., Beachy, R. N. & Fauquet, C. M. Tomato leaf curl geminivirus from India has a bipartite genome and coat protein is not essential for infectivity. J. Gen. Virol. 76, 25–35 (1995).CAS 
PubMed 

Google Scholar 
20.Rochester, D. E., DePaulo, J. J., Fauquet, C. M. & Beachy, R. N. Complete nucleotide sequence of the geminivirus tomato yellow leaf curl virus Thailand isolate. J. Gen. Virol. 75, 477–485 (1994).CAS 
PubMed 

Google Scholar 
21.Chakraborty, S., Pandey, P. K., Banerjee, M. K., Kalloo, G. & Fauquet, C. M. Tomato leaf curl Gujarat virus a new begomovirus species causing a severe leaf curl disease of tomato in Varanasi India. Phytopathology 93, 1485–1495 (2003).CAS 
PubMed 

Google Scholar 
22.Sattar, M. N. et al. First identification of begomoviruses infecting tomato with leaf curl disease in Burkina Faso. Plant Dis. 99, 732–732 (2015).
Google Scholar 
23.Ouattara, A. et al. Tomato leaf curl Burkina Faso virus: a novel tomato-infecting monopartite begomovirus from Burkina Faso. Arch. Virol. 162, 1427–1429 (2017).CAS 
PubMed 

Google Scholar 
24.Tiendrébéogo, F. et al. Molecular and biological characterization of pepper yellow vein Mali virus (PepYVMV) isolates associated with pepper yellow vein disease in Burkina Faso. Arch. Virol. 156, 483–487 (2011).PubMed 

Google Scholar 
25.Chen, L.-F. et al. A severe symptom phenotype in tomato in Mali is caused by a reassortant between a novel recombinant begomovirus (Tomato yellow leaf curl Mali virus ) and a betasatellite. Mol. Plant Pathol. 10, 415–430 (2009).CAS 
PubMed 
PubMed Central 

Google Scholar 
26.Rojas, M. R. et al. Functional analysis of proteins involved in movement of the monopartite begomovirus, tomato yellow leaf curl virus. Virology 291, 110–125 (2001).CAS 
PubMed 

Google Scholar 
27.Ranjan, P., Kumar, R. V. & Chakraborty, S. Differential pathogenicity among tomato leaf curl Gujarat virus isolates from India. Virus Genes 47, 524–531 (2013).CAS 
PubMed 

Google Scholar 
28.Jyothsna, P. et al. Infection of tomato leaf curl New Delhi virus (ToLCNDV), a bipartite begomovirus with betasatellites, results in enhanced level of helper virus components and antagonistic interaction between DNA B and betasatellites. Appl. Microbiol. Biotechnol. 97, 5457–5471 (2013).CAS 
PubMed 

Google Scholar 
29.Duan, Y. P., Powell, C. A., Purcifull, D. E., Broglio, P. & Hiebert, E. Phenotypic variation in transgenic tobacco expressing mutated geminivirus movement/pathogenicity (BC1) proteins. Mol. Plant- Microbe Interact. 10, 1065–1074 (1997).CAS 
PubMed 

Google Scholar 
30.Hussain, M., Mansoor, S., Iram, S., Fatima, A. N. & Zafar, Y. The nuclear shuttle protein of tomato leaf curl New Delhi virus is a pathogenicity determinant. J. Virol. 79, 4434–4439 (2005).CAS 
PubMed 
PubMed Central 

Google Scholar 
31.Geoghegan, J. L. & Holmes, E. C. The phylogenomics of evolving virus virulence. Nat. Rev. Genet. 19, 756–769 (2018).CAS 
PubMed 
PubMed Central 

Google Scholar 
32.Péréfarres, F. et al. Frequency-dependent assistance as a way out of competitive exclusion between two strains of an emerging virus. Proc. R. Soc. B Biol. Sci. 281, 1–9 (2014).
Google Scholar 
33.Wang, H. L., Gilbertson, R. L. & Lucas, W. J. Spatial and temporal distribution of bean dwarf mosaic geminivirus in Phaseolus vulgaris and Nicotiana benthamiana. Phytopathology 86, 1204–1214 (1996).
Google Scholar 
34.Hanley-Bowdoin, L., Bejarano, E. R., Robertson, D. & Mansoor, S. Geminiviruses: Masters at redirecting and reprogramming plant processes. Nat. Rev. Microbiol. 11, 777–788 (2013).CAS 
PubMed 

Google Scholar 
35.Londono, A., Riego-Ruiz, L. & Arguello-Astorga, G. R. DNA-binding specificity determinants of replication proteins encoded by eukaryotic ssDNA viruses are adjacent to widely separated RCR conserved motifs. Arch. Virol. 155, 1033–1046 (2010).CAS 
PubMed 

Google Scholar 
36.Gutiérrez, S., Michalakis, Y., Van Munster, M. & Blanc, S. Plant feeding by insect vectors can affect life cycle, population genetics and evolution of plant viruses. Funct. Ecol. 27, 610–622 (2013).
Google Scholar 
37.Lee, C. H. et al. A single amino acid substitution in the movement protein enables the mechanical transmission of a geminivirus. Mol. Plant Pathol. 00, 1–18 (2020).
Google Scholar 
38.Froissart, R., Doumayrou, J., Vuillaume, F., Alizon, S. & Michalakis, Y. The virulence-transmission trade-off in vector-borne plant viruses: A review of (non-)existing studies. Philos. Trans. R. Soc. B 365, 1907–1918 (2010).CAS 

Google Scholar 
39.Ditta, G., Stanfield, S. & Corbin, D. Broad host range DNA cloning system for Gram-negative bacteria: Construction of a gene bank of Rhizobium meliloti. Proc. Natl. Acad. Sci. 77, 7347–7351 (1980).CAS 
PubMed 
PubMed Central 

Google Scholar 
40.Lapidot, M., Cohen, L., Machbash, Z. & Levy, D. Development of a scale for evaluation of tomato yellow leaf curl virus resistance level in tomato plants. Phytopathology 96, 1404–1408 (2006).CAS 
PubMed 

Google Scholar 
41.Vernerey, M. S., Pirolles, E., Blanc, S. & Sicard, A. Localizing genome segments and protein products of a multipartite virus in host plant cells. Bio-Protoc. 9, 1–14 (2019).
Google Scholar 
42.Lefeuvre, P., Hoareau, M., Delatte, H., Reynaud, B. & Lett, J. M. A multiplex PCR method discriminating between the TYLCV and TYLCV-Mld clades of tomato yellow leaf curl virus. J. Virol. Methods 144, 165–168 (2007).CAS 
PubMed 

Google Scholar 
43.Urbino, C. et al. Within-host dynamics of the emergence of tomato yellow leaf curl virus recombinants. PLoS ONE 8, 1–14 (2013).
Google Scholar 
44.Conflon, D. et al. Accumulation and transmission of alphasatellite, betasatellite and tomato yellow leaf curl virus in susceptible and Ty-1 resistant tomato plants. Virus Res. 253, 124–134 (2018).CAS 
PubMed 

Google Scholar 
45.R Development Core Team. R: A language and environment for statistical computing. (2017).46.Pinheiro, J. C., Bates, D. M., DebRoy, S., Sarkar, D. & Team, R. C. nlme: Linear and nonlinear mixed effects models. (2016).47.Hothorn, T., Bretz, F. & Westfall, P. Simultaneous inference in general parametric models. Biometrical J. 50, 346–363 (2008).MathSciNet 
MATH 

Google Scholar  More

Natural diets of tropical birds vary within speciesWe collected 62 regurgitated samples (using the tartar emetic method 22) from multiple tropical bird species representing four bird orders (Columbiformes–Pigeons, Coraciiformes–Kingfishers, Psittaciformes–Parrots, and Passeriformes–Passerines). First, we characterized diet components visually and then through metabarcoding of 52 of these samples using universal primers targeting invertebrates (Cytochrome c oxidase subunit I: COI gene) and plants (Internal transcribed spacer 2: ITS2 gene) (Table S1 and Fig. 2). Through visual identification, we identified plant material in 26 samples. The most common visually identified invertebrate orders were Araneae (spiders—27 samples), and Coleoptera (beetles—27 samples) (Table S2). Metabarcoding sequences were analysed using the OBITools software25. Overall, we found 47 plant operational taxonomic units (OTUs—97% sequence similarity threshold) and 180 invertebrate OTUs (Table S3). Plant items were dominated by the orders Rosales (27.7% OTUs), Fabales (8.5% OTUs), and Sapindales (8.5% OTUs). Except for four OTUs, all plants were identified to the genus level. Of the invertebrate OTUs, 54 belonged to feather mites (known feather symbionts), endoparasites, and rotifers (likely due to accidental consumption along with drinking water), and these OTUs were removed from further analyses, leaving 126 potential dietary invertebrate OTUs. Invertebrate samples were dominated by the classes Insecta (67.5% OTUs) and Arachnida (28.6% OTUs). At the order-level, dietary items were mainly represented by Araneae (spiders—28.6% OTUs), Hemiptera (true bugs—15.9% OTUs), Diptera (flies—14.3% OTUs), and Lepidoptera (moths and butterflies—10.3% OTUs). However, 77% of the invertebrate OTUs could not be identified to genus level, highlighting the limited research on genotyping invertebrate communities in Papua New Guinea.Figure 2Natural diets of wild birds vary between individuals of the same species and the results of the two identification methods of dietary components (visual identification and metabarcoding). Relative abundances based on the presence/absence of data of different dietary components are indicated in colours. Only invertebrates are separated into taxonomic orders as visual identification is unable to identify plant orders. Individuals depicted with asterisks had both crop microbiome and diet samples (dataset 1), while black font represents individuals with both cloacal microbiomes and diet samples (dataset 2). Individuals are clustered according to the species (each species is given a six-letter code name) and their literature-based dietary guilds. The order of the species is indicated with illustrations (Columbiformes–Pigeons, Coraciiformes–Kingfishers, Passeriformes–Passerines and Psittaciformes–Parrots), while ‡ represents diet samples with a complete consensus between the two identification methods.Full size imageDiet item identification differed markedly between visual and metabarcoding methods (Fig. 2, Tables S2 and S3). The diet components of individuals also varied notably within species (Figs. 2 and S1). Only diets of 12 out of 52 individuals were fully congruent between the two methods (Fig. 2). Of these 12 samples, eight had only plant material. Identification of invertebrate orders also differed between the two methods (Fig. 2, Table 1). Both methods identified the arthropod orders Hemiptera, Diptera, Orthoptera (crickets and locusts), and Araneae in the same samples (Fig. 2 and Table 1), while metabarcoding detected lower proportions of Coleoptera than the visual identification (Table 1).Table 1 Comparison between diet items identified in the regurgitated samples from the two approaches (visual identification and metabarcoding).Full size tableComparison of microbiomes and consumed diet itemsFor subsequent comparisons of diets and microbiomes, we utilised individual datasets from both visual identification (diet components identified at the order level) and metabarcoding (both OTU and order level), and a combination (order level) of both approaches (for details see “Methods” section on identifying prey items). Due to differences between the diet identification methods, a combination of the results was used to circumscribe the full diversity of consumed diets and to account for inherent biases associated with the two methods (i.e., the inability to identify plant material and smaller body parts of invertebrates visually, and extraction and sequencing biases associated with metabarcoding). We separated the microbiome dataset into three datasets due to sequencing limitations: dataset 1 included 12 birds with successfully sequenced crop microbiomes and diets identified using both methods, dataset 2 included 27 birds with successfully sequenced cloacal microbiomes and diets, and dataset 3 included 17 birds for which we obtained successfully sequenced crop and cloacal microbiomes (Table S1). Prior to subsequent analyses, each microbiome dataset was rarefied to even sequencing depths using the sample with the lowest number of sequences26 (Fig. S2).Crop microbiome similarity did not align with the consumed diet similarity (dataset 1)Out of the collected crop samples (N = 62), samples from only 19 individuals were successfully sequenced for their microbiomes. Of these individuals, we acquired diet samples for 12 individuals. Bacterial 16S rRNA MiSeq sequences were analysed using the DADA2 pipeline27 within QIIME228. There were 351,867 bacterial sequences (mean ± SD: 29,322 ± 33,009) in the crop microbiomes prior to rarefaction (Table S4). After rarefaction, bacterial sequences were identified to 615 amplicon sequence variants (ASVs—100% sequence similarity). Crop microbiomes were dominated by Proteobacteria (53.6%), Actinobacteria (18.9%), and Firmicutes (17.9%). Alpha diversities of individual microbiomes were calculated using the diversity function in the microbiome package29 and they did not differ significantly between host orders [Chao1 richness: Kruskal Wallis (KW) χ2 = 4.559, df = 3, p = 0.2271; Shannon’s diversity index: χ2 = 2.853, df = 3, p = 0.4149], or literature-based dietary guilds (Chao1 richness: KW χ2 = 4.317, df = 2, p = 0.1155; Shannon’s diversity index: KW χ2 = 2.852, df = 2, p = 0.2403) (Fig. S3).The compositional differences of crop microbiomes were investigated with the adonis2 function in the vegan package30 using permutational multivariate analyses of variance tests (PERMANOVA). These analyses revealed that the bird host order did not influence the crop microbiome composition (PERMANOVA10,000 permutations: Bray–Curtis: F = 1.251, R2 = 0.0993, p = 0.1911; Jaccard: F = 1.154, R2 = 0.0962, p = 0.2191) (Fig. S1). The effect of feeding guild was masked by host order as they are strongly correlated in this dataset. Furthermore, the lack of an effect of host taxa on crop microbiomes may be a result of the small sample sizes.We further investigated whether alpha diversity of the crop microbiomes was influenced by the diet item diversity of individuals. The Chao1 richness estimates of the microbiomes and the richness of the consumed diet items (number of different diet items based on the combined results) of individuals were not significantly correlated (Table S5), suggesting that the diet richness does not impact crop microbiome richness. However, Shannon’s diversity index of crop microbiomes and diet diversity were marginally significantly negatively associated (Table S5). This suggests that despite the lack of an association between diet and microbiome richness, crop microbiome evenness could be influenced by diet diversity.We then explored the association between the crop microbiome composition and the consumed diets, investigating correlations between Bray–Curtis and Jaccard dissimilarities of microbiomes, and Jaccard dissimilarity of diets using Mantel tests in the vegan package30. The compositional similarity of the diets based on any of the methods (visual, metabarcoding—both OTU and order-level separately, and combined) did not correlate significantly with crop microbiome compositions (Table 2 and Fig. S4). We observed similar non-significant associations between diets and microbiomes when investigating host orders separately (Table S6). This suggests that overall crop microbiomes of individuals are not completely modelled by the composition of the consumed diets.Table 2 Results of Mantel tests between the crop (dataset 1) and the cloacal (dataset 2) microbiome similarities (measured with both Bray–Curtis and Jaccard distances) and the consumed diet similarities (measured with Jaccard distances).Full size tableHost-taxon specific cloacal microbes are associated with different diet items (dataset 2)We obtained 27 individuals from 15 bird species with successfully sequenced cloacal microbiomes and diet samples (based on both metabarcoding and visual identification). Prior to rarefying, we acquired 818,272 bacterial sequences from the cloacal swab samples (mean ± SD: 30,306 ± 20,903) (Table S7). After rarefaction, bacterial sequences were assigned to 1,324 ASVs that belonged to Actinobacteria (35.9%), Proteobacteria (32.6%), Firmicutes (21.2%) and Tenericutes (5.0%). Cloacal microbiome alpha diversity did not differ significantly between different bird orders (Chao1 richness: KW χ2 = 2.624, df = 3, p = 0.4532; Shannon’s diversity: χ2 = 6.595, df = 3, p = 0.0861) or literature-based dietary guilds (Chao1 richness: KW χ2 = 1.128, df = 3, p = 0.7703; Shannon’s diversity: KW χ2 = 1.673, df = 3, p = 0.6429) (Fig. S5).However, cloacal microbiome beta diversity was significantly influenced by host bird order (PERMANOVA10,000 permutations: Bray–Curtis: F = 2.159, R2 = 0.2055, p  More