Ecology

Study area and period
Our study area is the tropics between 15°N–35°S6. We divided the study area into three continents and studied them separately: South America, Africa, and Australasia. Australasia includes Australia and southeast Asia, but excludes southern India. Our results are generated on 0.25° spatial resolution. We classify a cell as forest if it contained at least 50% tree cover (‘forest cover’ in this manuscript) in 1999 according to the dataset from ref. 41. The moisture recycling simulations were carried out for 2003–2014 (‘recent climate’), for which a consistent set of input data was available (see also ref. 11). ‘Late 21st century’ refers to 2071–2100.
Local-scale forest hysteresis
Previous research has shown that tropical forests may have local-scale tipping points at certain mean annual rainfall levels, but are also affected by the seasonality of that rainfall5,6,8. Local-scale tipping points for forest were determined using tree cover data following a method from ref. 6. Using potential analysis42, an empirical stability landscape (as in Fig. 1a) is constructed based on spatial distributions of tree cover against environmental variables such as mean annual rainfall for each continent separately. For each value of the environmental variable, the probability density of tree cover was determined using the MATLAB function ksdensity with a bandwidth of 5%. We applied Gaussian weights to the environmental variable with a standard deviation of 0.05 times the length of the axis of the environmental variable. Local maxima of the resulting probability density function are interpreted as stable states, where we ignored local maxima below a threshold value of 0.004. We used Landsat tree cover data for 2000 on 30 m resolution downloaded for every 0.01°43. We masked out human-used areas, water bodies, and bare ground using the ESA GlobCover land cover dataset for 2009 on 300 m resolution (values 11–30 and ≥190). From the resulting dataset we randomly sampled one million locations for each continent and used them to construct the stability landscapes6 against mean annual rainfall and average MCWD. MCWD is the cumulative difference between evapotranspiration and rainfall using monthly averages of those fluxes calculated for each calendar year44. It is set to zero when monthly rainfall exceeds monthly evapotranspiration and becomes more negative with an increasing water deficit. Following ref. 11, for both mean annual rainfall and MCWD, we took monthly data from the GLDAS 2.0 dataset45 for 1970–1999 so the 30-year period leading up to the land-cover sample (for the year 2000) was used.
Forest evapotranspiration
To estimate the fraction of evapotranspiration attributable to forest cover we used the large-scale hydrological model PCR-GLOBWB, run at 0.5° resolution46. Per grid cell, the model simulates evapotranspiration for a range of land-cover types. Here, we are specifically interested in the evapotranspiration of forests, or ‘tall natural vegetation’ in PCR-GLOBWB. Note that we here account for both forest transpiration and canopy interception evaporation instead of, as in ref. 11, only transpiration.
PCR-GLOBWB computes the water balance in two soil layers and a groundwater layer. Soil type, fractional area of saturated soil, and the spatiotemporal distribution of groundwater depth are accounted for (see refs. 46,47). It includes six land-cover types, with spatially varying parameters46: tall and short natural vegetation, pasture, rainfed crops, and paddy and non-paddy irrigated crops. The model was forced with WATCH Forcing Data ERA-Interim precipitation, temperature, and reference potential evapotranspiration for 1979–201448. We used monthly evapotranspiration output of PCR-GLOBWB, implying that we assume that forest component of evapotranspiration remains equal within each month. For detailed model descriptions and validation, we refer to earlier studies11,46,49,50.
Atmospheric moisture tracking
As an essential step in estimating the forest–rainfall feedback, we determined where the moisture from enhanced evapotranspiration precipitates again by using atmospheric moisture tracking. The method for atmospheric moisture tracking resembles that in ref. 11. Apart from the expansion of the study area to the entire tropics, a notable difference is that we here used ERA5 reanalysis data rather than ERA-Interim, meaning that the simulations were based on finer resolution input data (i.e. on 0.25° instead of 0.75° for spatial resolution, and 1 h instead of 3 h for temporal resolution). ERA5 has better performance than ERA-Interim regarding wind fields and rainfall, especially in the tropics51,52,53. Below we summarize the method (see also ref. 11).
We used a Lagrangian method of moisture tracking that is based on previous studies11,54,55,56 that track parcels of evaporated moisture forward through the atmosphere to their subsequent precipitation location. Moisture particles that enter the atmosphere are assigned a random location within the 0.25° grid cell and random starting height in the atmosphere scaled with the humidity profile, and their trajectories are then tracked through the atmosphere. The trajectories are forced with the three-dimensional ERA5 reanalysis estimates of wind speed and direction, which were linearly interpolated at every time step of 0.25 h. Water particles in the atmosphere have an equal probability of raining out, regardless of vertical position. Rainfall A (mm per time step) at location x,y and time t that has evaporated from any location of release in any cell is ref. 56

$$A_{x,y,t} = P_{x,y,t}frac{{W_{{mathrm{parcel}},t}E_{{mathrm{source}},t}}}{{{mathrm{TPW}}_{x,y,t}}},$$
(1)

where P is rainfall in mm per time step, Wparcel is the water in the tracked parcel in mm, Esource is its fraction of water that evapotranspired from the source, and TPW is the precipitable water in the atmospheric water column in mm. Every time step, the amount of water in the parcel is updated based on evapotranspiration ET into the parcel and rainfall P from it:

$$W_{{mathrm{parcel}},t} = W_{{mathrm{parcel}},t – 1} + ({mathrm{ET}}_{x,y,t} – P_{x,y,t})frac{{W_{{mathrm{parcel}},t – 1}}}{{{mathrm{TPW}}_{x,y,t}}}.$$
(2)

The fraction of water in the parcel that has evapotranspired from the source then becomes

$$E_{{mathrm{source}},t} = frac{{E_{{mathrm{source}},t – 1}W_{{mathrm{parcel}},t – 1} – A_{x,y,t}}}{{W_{{mathrm{parcel}},t}}}.$$
(3)

Thus, the amount of water that was tracked from the source location decreases with precipitation along its trajectory. Parcels were followed until either less than 5% of its original amount was left in the atmosphere, or the tracking time was 30 days. Any moisture remaining in the parcel when the trajectories end is assumed to rain out over non-land areas, thus not contributing to our analysis. We analysed each continent separately for reasons of computability. However, small moisture flows between forests in different continents can be expected, as has been simulated for flows from Africa to the Amazon57. Over all land points, ERA5 hourly evapotranspiration is linearly interpolated to every 0.25-h time step. The moisture flow mij in mm per month linking evapotranspiration in cell i to rainfall in cell j where (left[ {x,y} right] , {it{epsilon }} , j) over the course of a given month becomes

$$m_{ij} = mathop {sum}limits_{t = 0}^{{mathrm{month}}} {A_{j,t}} cdot frac{{{mathrm{ET}}_{i,t}}}{{W_{i,t}}},$$
(4)

where ETi,t is the evapotranspiration in mm per time step, and Wi,t is the tracked amount of water from source cell i at time step t.
At continental scales, evaporated moisture can rain down and re-evaporate multiple times. This also means that forest evapotranspiration can enhance rainfall multiple times. We accounted for this ‘cascading moisture recycling’ following refs. 11,14, in which the rainfall attributed to an upwind forest source is subsequently tracked upon re-evaporation. After six re-evapotranspiration cycles, cascading moisture recycling has decreased to practically zero11. Therefore, following ref. 11, seven iterations of cascading moisture recycling were performed.
Hysteresis experiments
We determined the hysteresis of tropical forests through a series of iterative runs; each one started either from a fully forested continent or from a fully deforested continent. We simulated the hypothetical mean annual rainfall levels across the (tropical part of the) continent given this initial condition, that is, rainfall without any forest evapotranspiration or rainfall including evapotranspiration from an entirely forested continent. Next, based on the empirical bifurcation diagrams (i.e. nonforest, bistable forests, and stable forests in each continent are determined based on the bistability ranges shown in Supplementary Figs. 1−3), we determined either the minimal distribution of tropical forest (in case of a no-forest initial condition, based on the higher end of the bistability range) or the maximal distribution (in case of a fully forested initial condition, based on the lower end of the bistability range) at these rainfall levels. Thus, in the simulations with an empty initial condition, only stable forests (green in Fig. 1) would regrow; in those with a full initial condition, both stable and bistable forests (green and yellow in Fig. 1) would disappear. Because the resulting new distribution of forest would generate different levels of rainfall, the procedure was repeated with the respective forest distribution as initial condition. This occurred until rainfall levels had (practically) stabilized between iterative runs (up to three iterations).
We assumed that no other vegetation type replaces the forest in order to show the theoretically possible distributions of tropical forests. This may lead to an overestimation of the actual effects of forest on rainfall, especially if forests would be replaced by highly transpiring crops58. Furthermore, land-cover changes will alter wind patterns and therefore the expected coupling between forests through evapotranspiration and rainfall59. Fossil fuel emissions not only alter the climate, but the emitted CO2 also fertilizes plants. This increases trees’ water-use efficiency, reducing their water demand, but it also increases biomass production60. The effects of this CO2 fertilization on the water cycle might be small61, but its net effects on tropical forest hysteresis remains uncertain.
For display of Fig. 2, the resolution of rainfall values was increased by a factor of 2, to 0.125° and smoothed using a two-dimensional Gaussian kernel with a standard deviation of 0.5°.
Climate-change scenario
As the estimate of late 21st-century rainfall conditions, we used the rainfall output from the SSP5-8.5 scenario simulations by seven available CMIP6 models62. These models are BCC-CSM2-MR, CanESM5, CNRM-CM6-1, CNRM-ESM2, IPSL-CM6A-LR, MRI-ESM2.0, and UKESM1.0-LL. We took the mean across the model outputs for the annual rainfall values for 2071–2100. The scenario is considered a severe climate-change scenario. Because the moisture tracking model is forced with atmospheric reanalysis data, we assumed that (forest-induced) moisture flows in the scenario are the same as in the period of our atmospheric simulations (2003–2014).
We assumed that a tipping point from a forested to a nonforested state occurs when mean annual rainfall in a forested cell (forest cover ≥ 50%) is currently (2003–2014) above the lower tipping point (Supplementary Figs. 1–3), but is reduced to below the lower tipping point in the climate-change scenario. Similarly, a tipping point from a nonforested to a forested state occurs when mean annual rainfall in a nonforested cell (forest cover  More

Scopes and databases
The CEAP dataset comprises all the thermal power plants operating in China, totalling 2,714 plants (or 6,267 units), from 2014 to 2017 in 26 provinces and 4 municipalities (except Hong Kong, Macao, Taiwan and Tibet; Table 1). The thermal power plants produce electricity by combusting a variety of fossil energies, which fall into 4 categories: coal, gas plus oil, biomass and others (detailed in Table 2).
Table 1 China’s thermal power plants in CEAP.
Full size table

Table 2 Fuel type descriptions.
Full size table

The CEAP dataset integrates two databases, i.e., the CEMS data and unit-specific information. The CEMS data—the direct, real-time measurements of stack gas concentrations of PM, SO2 and NOX from China’s power plant stacks—are monitored by China’s CEMS network and reported to the China Ministry of Ecology and Environment (MEE; http://www.envsc.cn/). The CEMS data are recorded on a source and hourly basis. In total, the CEMS dataset covers 4,622 emission sources (i.e., power plant stacks) associated with 5,606 units (accounting for 98% of China’s thermal power capacity), 35,064 hours from 2014 to 2017, and 3 air pollutants (i.e., PM, SO2 and NOX) for each source-hour sample (Table 3). The MEE has also provided stack-specific information (regarding latitude and longitude, heights, temperature, diameter, etc.; http://permit.mee.gov.cn/).
Table 3 CEMS coverage of China’s thermal power plant units or stacks in CEAP.
Full size table

Unit-specific information is also derived from the MEE, involving activity levels (energy consumption and power generation), operating capacities, geographic allocations and pollution control equipment (particularly the types and removal efficiencies) at a yearly frequency. Due to data availability, the unit information is available only until 2016, and the activity levels for 2017 are projected following the overall trends in provincial thermal power generation between 2016 and 2017 (which are available in the China Energy Statistical Yearbooks26), under the assumption that new units constructed in 2017 have the same structures of installed capacities, energy uses and regions as those of the existing units in 2016.
With a combination of the two datasets, the CEAP dataset provides nationwide, plant-level, dynamic PM, SO2 and NOX emissions from China’s thermal power plants from 2014 to 2017. Relative to existing inventories, the CEAP dataset is innovative in that it incorporates comprehensive real CEMS-measured emission data, avoiding the use of average emission factors and the associated operational assumptions and uncertain parameters.
Pre-processing of CEMS data
We have been exclusively granted access to the data from China’s CEMS network. Generally, the CEMS consists of a sampling system (for filtering and sampling flue gas), an online analytical component (for monitoring flue gas parameters, particularly emission concentrations) and a data processing system (for collecting, processing and reporting monitoring data)27,28. According to the GB13223-2003 regulation29, the CEMS network should cover all power plant furnaces that burn coal (except stoker and spreader stoker) and oil and generate >65 tons of steam each hour, as well as those that burn pulverized coal and gas. Thus, some power plants have not yet been incorporated into the CEMS network (accounting for 3–4% of the total thermal power capacity from 2014 to 2017) because their furnaces did not meet the requirements necessary to install a CEMS. For the power plants outside the CEMS network, we assume their stack concentrations are similar to the averages of the units with similar fuel types and similar regions within the CEMS network.
To guarantee the reliability of CEMS data, China’s government has made great efforts in developing specific regulations and technical guidelines for power plants and local entities to follow and supervise, respectively24,28,30,31,32. These official documents elaborate on all the processes required to regulate the CEMS network, including not only CEMS installation, operation, inspection, maintenance and repair but also CEMS data collection, processing, reporting, analysis and storage28,32,33. Since 2014, all state-monitored companies have been mandated to report their CEMS data to the local governments through a series of online platforms for different provinces (listed in Supplementary Table 1). Local entities have random onsite inspections to check the truthfulness of the reported results on at least a quarterly basis23,24,28,32,34; this system enables a comparison of CEMS data across different firms to explore potential outliers and abnormalities and prevent data manipulation28,35. Then, the governments release the inspection results to the public through the same online platforms (listed in Supplementary Table 1)24,36,37. Severe financial penalties and criminal punishments can be imposed on firms that adopt data manipulation (in terms of deleting, distorting and forging CEMS data, for example)38,39.
The malfunction of CEMS monitors may also introduce large uncertainty to CEMS data during the processes of operation (indication errors, span drift, zero drift, etc.), maintenance (particularly the failure to perform calibration and reference tests) and data reporting (invalid data communication, data missing, etc.)24,28. Accordingly, each power plant is required to make at least one A-, B- and C-grade overhaul for 32–80, 14–50 and 9–30 days per 4–6, 2–3 and 1 year(s), respectively, as well as one D-grade overhaul (if needed) for 5–15 days per year, to check, maintain and upgrade its technologies, thereby reducing measurement uncertainty40. During these overhauls, CEMS operators conduct CEMS calibration (i.e., zero and span calibration), maintenance procedures (e.g., examining and cleaning major CEMS components and replacing or upgrading parts, if necessary, such as optical lens, filter and sampling meter) and a reference test (i.e., relative accuracy test audit). Furthermore, third-party operators examine CEMS operation and maintenance routines, to guarantee standardized CEMS operation and facilitate improvement in CEMS data accuracy27,28,31. All the related activities should be documented according to standardized requirement contents27,28. Even with the aforementioned efforts, there is still a small proportion of nulls and outliers in the CEMS database, which represent 1% and 0.1% of the total operating hours, respectively, from 2014 to 2017. We treat these samples seriously by following the relevant official documents, which have been released by China’s government. Table 4 provides the treatment methods for nulls or zeros, which can be divided into 3 types based on duration. On the one hand, we consider nulls and/or zeros that span at least 5 successive days as a downtime or overhaul and omit them in the estimation, according to the regulation27. On the other hand, missing data lasting  24 hours are assumed around the valid values near the time and set to the monthly averages27:

$${hat{C}}_{s,i,y,m,h}={bar{C}}_{s,i,y,m,bullet }$$
(1)

where ({C}_{s,i,y,m,h}) denotes the stack gas concentrations of pollutant s emitted by unit i for year y, month m and hour h (i.e., the actual measurements monitored by the CEMS network), defined as the amount of pollutants per unit of emitted stack gas (g m−3)41,42; ({widehat{C}}_{s,i,y,m,h}) is the imputation for the missing data ({C}_{s,i,y,m,h}); ({bar{C}}_{s,i,y,m,.}) is the mean of the hourly valid values for the same pollutant, unit, year and month as ({C}_{s,i,y,m,h}). In contrast, the missing data for 1–24 hour(s) are interpolated with the arithmetic averages of the two nearest valid points before and after them27,43:

$${widehat{C}}_{s,i,y,m,h}=frac{{C}_{s,i,y,m,h-l}+{C}_{s,i,y,m,h+q}}{2}$$
(2)

where ({C}_{s,i,y,m,h-l}) and ({C}_{s,i,y,m,h{rm{+}}q}) represent the nearest last known measurements (l hour(s) before) and next known measurements (q hour(s) after), respectively, for the missing data ({C}_{s,i,y,m,h}), namely, the series data ({C}_{s,i,y,m,h-l+1}),…, ({C}_{s,i,y,m,h}),…, ({C}_{s,i,y,m,h+q-1}) are all missing values. Furthermore, we treat the measurements that are out of the measurement ranges of CEMS instruments (outside of which the data are unreliable30,44; detailed in Supplementary Table 2) as abnormal data and process them in a similar way to nulls according to the official regulation27.
Table 4 Treatment methods for nulls and the relevant official documents.
Full size table

CEMS-based estimation of emission factors and absolute emissions
The introduction of real CEMS-monitored measurements provides a direct estimation for emission factors on a source and hourly basis, avoiding the use of average emission factors with many assumptions and uncertain parameters17,42,44.

$$E{F}_{s,i,y,m,h}={C}_{s,i,y,m,h}{V}_{i,y}$$
(3)

In Eq. (3), (E{F}_{s,i,y,m,h}) indicates the emission factor, defined as the amount of emissions per unit of fuel use (in g kg−1 for solid or liquid fuel and in g m−3 for gas fuel), and ({V}_{i,y}) is the theoretical flue gas rate, defined as the expected volume of flue gas per unit of fuel use under standard production conditions (m3 kg−1 for solid or liquid fuel and m3 m−3 for gas fuel)42, which was estimated by the China Pollution Source Census (2011)45 based on sufficient field measurements (detailed in Table 5). Based on Eq. (3), abated emission factors can be directly obtained even without the use of removal efficiencies and the relevant parameters, because CEMS monitors the gas concentrations at stacks after the effect of control equipment (if any).
Table 5 Theoretical flue gas rate.
Full size table

Notably, recent clean air policies (particularly different emissions standards) target stack concentrations, such that a large proportion of missing data exist regarding other measurements (particularly flue gas rates, with missing data accounting for 34.62%, 31.91%, 29.97% and 42.96% of the total samples in 2014, 2015, 2016 and 2017, respectively). Accordingly, we introduce theoretical flue gas rates into the estimation to avoid significant underestimation of the actual volume when there are too many missing data values46. In addition, the adoption of theoretical flue gas rates can address flue gas leakage, a common problem in power plants that greatly distorts the real flue gas volume46. The theoretical flue gas rates are derived from the China Pollution Source Census, with values varying across operating capacities, fuel types and boiler types42,45. Thus, the actual volume of flue gas is computed in terms of the theoretical flue gas rate times actual fuel consumption.
The absolute emissions of PM, SO2 and NOX from individual power plants can be estimated in terms of the emission factors times the activity levels21:

$${E}_{s,i,y,m}=E{F}_{s,i,y,m}{A}_{i,y,m}$$
(4)

where ({E}_{s,i,y,m}) represents the air pollution emissions (g); and ({A}_{i,y,m}) is the activity data, i.e., the amount of fuel use (kg for solid or liquid fuel and m3 for gas fuel). In the CEAP dataset, power plant emissions are estimated on a monthly basis (the smallest scale for activity data), in which the yearly unit-level activity data are allocated at the monthly scale using the monthly province-level thermal power generation as weights16:

$${A}_{i,y,m}=frac{{F}_{{p}_{i},y,m}}{{sum }_{m=1}^{12}{F}_{{p}_{i},y,m}}{A}_{i,y}$$
(5)

where ({F}_{{p}_{i},y,m}) denotes the thermal power generation by province Pi, which is obtained from the Chinese Energy Statistics Yearbook26, and ({p}_{i}) indicates the province where unit i is located. More

1.
Daws, M. I., Mullins, C. E., Burslem, D. F. R. P., Paton, S. R. & Dalling, J. W. Topographic position affects the water regime in a semideciduous tropical forest in Panamá. Plant Soil 238, 79–89. https://doi.org/10.1023/A:1014289930621 (2002).
CAS  Article  Google Scholar 
2.
Sundqvist, M. K., Sanders, N. J. & Wardle, D. A. Community and ecosystem responses to elevational gradients: processes, mechanisms, and insights for global change. Annu. Rev. Ecol. Evol. Syst. 44, 261–280. https://doi.org/10.1146/annurev-ecolsys-110512-135750 (2013).
Article  Google Scholar 

3.
Jucker, T. et al. Topography shapes the structure, composition and function of tropical forest landscapes. Ecol. Lett. 21, 989–1000. https://doi.org/10.1111/ele.12964 (2018).
Article  PubMed  PubMed Central  Google Scholar 

4.
Moeslund, J. E., Arge, L., Bøcher, P. K., Dalgaard, T. & Svenning, J.-C. Topography as a driver of local terrestrial vascular plant diversity patterns. Nord. J. Bot. 31, 129–144. https://doi.org/10.1111/j.1756-1051.2013.00082.x (2013).
Article  Google Scholar 

5.
Holland, P. G. & Steyn, D. G. Vegetational responses to latitudinal variations in slope angle and aspect. J. Biogeogr. 2, 179–183. https://doi.org/10.2307/3037989 (1975).
Article  Google Scholar 

6.
Yetemen, O., Istanbulluoglu, E. & Duvall, A. R. Solar radiation as a global driver of hillslope asymmetry: Insights from an ecogeomorphic landscape evolution model. Water Resour. Res. 51, 9843–9861. https://doi.org/10.1002/2015wr017103 (2015).
ADS  Article  Google Scholar 

7.
Bennie, J., Hill, M. O., Baxter, R. & Huntley, B. Influence of slope and aspect on long-term vegetation change in British chalk grasslands. J. Ecol. 94, 355–368. https://doi.org/10.1111/j.1365-2745.2006.01104.x (2006).
Article  Google Scholar 

8.
Cantlon, J. E. Vegetation and microclimates on north and south slopes of Cushetunk Mountain, New Jersey. Ecol. Monogr. 23, 241–270. https://doi.org/10.2307/1943593 (1953).
Article  Google Scholar 

9.
Warren, R. J. Mechanisms driving understory evergreen herb distributions across slope aspects: as derived from landscape position. Plant Ecol. 198, 297–308. https://doi.org/10.1007/s11258-008-9406-1 (2008).
Article  Google Scholar 

10.
Bennie, J., Huntley, B., Wiltshire, A., Hill, M. O. & Baxter, R. Slope, aspect and climate: Spatially explicit and implicit models of topographic microclimate in chalk grassland. Ecol. Model. 216, 47–59. https://doi.org/10.1016/j.ecolmodel.2008.04.010 (2008).
Article  Google Scholar 

11.
Burnett, B. N., Meyer, G. A. & McFadden, L. D. Aspect-related microclimatic influences on slope forms and processes, northeastern Arizona. J. Geophys. Res. Earth Surf. https://doi.org/10.1029/2007jf000789 (2008).
Article  Google Scholar 

12.
Geroy, I. J. et al. Aspect influences on soil water retention and storage. Hydrol. Process. 25, 3836–3842. https://doi.org/10.1002/hyp.8281 (2011).
ADS  Article  Google Scholar 

13.
Huang, Y.-M., Liu, D. & An, S.-S. Effects of slope aspect on soil nitrogen and microbial properties in the Chinese Loess region. CATENA 125, 135–145. https://doi.org/10.1016/j.catena.2014.09.010 (2015).
CAS  Article  Google Scholar 

14.
Lozano-García, B., Parras-Alcántara, L. & Brevik, E. C. Impact of topographic aspect and vegetation (native and reforested areas) on soil organic carbon and nitrogen budgets in Mediterranean natural areas. Sci. Total Environ. 544, 963–970. https://doi.org/10.1016/j.scitotenv.2015.12.022 (2016).
ADS  CAS  Article  PubMed  Google Scholar 

15.
Rasmussen, C. & Tabor, N. J. Applying a quantitative pedogenic energy model across a range of environmental gradients. Soil Sci. Soc. Am. J. 71, 1719–1729. https://doi.org/10.2136/sssaj2007.0051 (2007).
ADS  CAS  Article  Google Scholar 

16.
Broxton, P. D., Troch, P. A. & Lyon, S. W. On the role of aspect to quantify water transit times in small mountainous catchments. Water Resour. Res. https://doi.org/10.1029/2008wr007438 (2009).
Article  Google Scholar 

17.
Casanova, M., Messing, I. & Joel, A. Influence of aspect and slope gradient on hydraulic conductivity measured by tension infiltrometer. Hydrol. Process. 14, 155–164. https://doi.org/10.1002/(sici)1099-1085(200001)14:1%3c155::aid-hyp917%3e3.0.co;2-j (2000).
ADS  Article  Google Scholar 

18.
Gutiérrez-Jurado, H. A., Vivoni, E. R., Istanbulluoglu, E. & Bras, R. L. Ecohydrological response to a geomorphically significant flood event in a semiarid catchment with contrasting ecosystems. Geophys. Res. Lett. https://doi.org/10.1029/2007gl030994 (2007).
Article  Google Scholar 

19.
Wang, L., Wei, S., Horton, R. & Shao, M. A. Effects of vegetation and slope aspect on water budget in the hill and gully region of the Loess Plateau of China. CATENA 87, 90–100. https://doi.org/10.1016/j.catena.2011.05.010 (2011).
Article  Google Scholar 

20.
Pook, E. & Moore, C. The influence of aspect on the composition and structure of dry sclerophyll forest on Black Mountain, Canberra. ACT. Aust. J. Bot. 14, 223–242. https://doi.org/10.1071/BT9660223 (1966).
Article  Google Scholar 

21.
Armesto, J. J. & Martίnez, J. A. Relations between vegetation structure and slope aspect in the Mediterranean region of Chile. J. Ecol. 66, 881–889. https://doi.org/10.2307/2259301 (1978).
Article  Google Scholar 

22.
Badano, E. I., Cavieres, L. A., Molina-Montenegro, M. A. & Quiroz, C. L. Slope aspect influences plant association patterns in the Mediterranean matorral of central Chile. J. Arid Environ. 62, 93–108. https://doi.org/10.1016/j.jaridenv.2004.10.012 (2005).
ADS  Article  Google Scholar 

23.
Zapata-Rios, X., Brooks, P. D., Troch, P. A., McIntosh, J. & Guo, Q. Influence of terrain aspect on water partitioning, vegetation structure and vegetation greening in high-elevation catchments in northern New Mexico. Ecohydrology 9, 782–795. https://doi.org/10.1002/eco.1674 (2016).
Article  Google Scholar 

24.
Poulos, H. M. & Camp, A. E. Topographic influences on vegetation mosaics and tree diversity in the Chihuahuan Desert Borderlands. Ecology 91, 1140–1151. https://doi.org/10.1890/08-1808.1 (2010).
Article  PubMed  Google Scholar 

25.
Kutiel, P. & Lavee, H. Efffect of slope aspect on soil and vegetation properties along an aridity transect. Isr. J. Plant Sci. 47, 169. https://doi.org/10.1080/07929978.1999.10676770 (1999).
Article  Google Scholar 

26.
Sternberg, M. & Shoshany, M. Influence of slope aspect on Mediterranean woody formations: comparison of a semiarid and an arid site in Israel. Ecol. Res. 16, 335–345. https://doi.org/10.1046/j.1440-1703.2001.00393.x (2001).
Article  Google Scholar 

27.
Méndez-Toribio, M., Meave, J. A., Zermeño-Hernández, I. & Ibarra-Manríquez, G. Effects of slope aspect and topographic position on environmental variables, disturbance regime and tree community attributes in a seasonal tropical dry forest. J. Veg. Sci. 27, 1094–1103. https://doi.org/10.1111/jvs.12455 (2016).
Article  Google Scholar 

28.
Gallardo-Cruz, J. A., Pérez-García, E. A. & Meave, J. A. β-Diversity and vegetation structure as influenced by slope aspect and altitude in a seasonally dry tropical landscape. Landsc. Ecol. 24, 473–482. https://doi.org/10.1007/s10980-009-9332-1 (2009).
Article  Google Scholar 

29.
Paudel, S. & Vetaas, O. R. Effects of topography and land use on woody plant species composition and beta diversity in an arid Trans-Himalayan landscape, Nepal. J. Mt. Sci. 11, 1112–1122. https://doi.org/10.1007/s11629-013-2858-3 (2014).
Article  Google Scholar 

30.
Zhang, R. The Dry Valley of the Hengduan Mountains Regions (Science Press, Beijing, 1992).
Google Scholar 

31.
Vegetation, E. B. O. S. Sichuan Vegetation (People’s Publishing House of Sichuan, Beijing, 1980).
Google Scholar 

32.
Guan, W. et al. Vegetation classification and the main types of vegetation of the dry valley of Minjiang River. J. Mt. Res. 22, 679–686 (2004).
Google Scholar 

33.
Ma, K.-M. et al. Multiple-scale soil moisture distribution and its implications for ecosystem restoration in an arid river valley, China. Land Degrad. Dev. 15, 75–85. https://doi.org/10.1002/ldr.584 (2004).
CAS  Article  Google Scholar 

34.
Lu, T., Ma, K. M., Zhang, W. H. & Fu, B. J. Differential responses of shrubs and herbs present at the Upper Minjiang River basin (Tibetan Plateau) to several soil variables. J. Arid Environ. 67, 373–390. https://doi.org/10.1016/j.jaridenv.2006.03.011 (2006).
ADS  Article  Google Scholar 

35.
Xu, X.-L., Ma, K.-M., Fu, B.-J., Song, C.-J. & Liu, W. Relationships between vegetation and soil and topography in a dry warm river valley, SW China. CATENA 75, 138–145. https://doi.org/10.1016/j.catena.2008.04.016 (2008).
Article  Google Scholar 

36.
Solon, J., Degórski, M. & Roo-Zielińska, E. Vegetation response to a topographical-soil gradient. CATENA 71, 309–320. https://doi.org/10.1016/j.catena.2007.01.006 (2007).
Article  Google Scholar 

37.
Birkeland, P. W. Soils and Geomorphology (Oxford University Press, Oxford, 1984).
Google Scholar 

38.
Loik, M. E., Breshears, D. D., Lauenroth, W. K. & Belnap, J. A multi-scale perspective of water pulses in dryland ecosystems: climatology and ecohydrology of the western USA. Oecologia 141, 269–281. https://doi.org/10.1007/s00442-004-1570-y (2004).
ADS  Article  PubMed  Google Scholar 

39.
Schwinning, S. & Sala, O. E. Hierarchy of responses to resource pulses in arid and semi-arid ecosystems. Oecologia 141, 211–220. https://doi.org/10.1007/s00442-004-1520-8 (2004).
ADS  Article  PubMed  Google Scholar 

40.
Fernandez-Going, B. M., Harrison, S., Anacker, B. & Safford, H. Climate interacts with soil to produce beta diversity in Californian plant communities. Ecology 94, 2007–2018. https://doi.org/10.1890/12-2011.1 (2013).
CAS  Article  PubMed  Google Scholar 

41.
Araya, Y. N., Gowing, D. J. & Dise, N. Does soil nitrogen availability mediate the response of grassland composition to water regime?. J. Veg. Sci. 24, 506–517. https://doi.org/10.1111/j.1654-1103.2012.01481.x (2013).
Article  Google Scholar 

42.
Araya, Y. N. et al. A fundamental, eco-hydrological basis for niche segregation in plant communities. New Phytol. 189, 253–258. https://doi.org/10.1111/j.1469-8137.2010.03475.x (2011).
Article  PubMed  Google Scholar 

43.
Li, Y. J., Bao, W. K. & Wu, N. Spatial patterns of the soil seed bank and extant vegetation across the dry Minjiang River valley in southwest China. J. Arid Environ. 75, 1083–1089. https://doi.org/10.1016/j.jaridenv.2011.05.012 (2011).
ADS  Article  Google Scholar 

44.
Li, F., Bao, W., Liu, J. & Wu, N. Eco-anatomical characteristics of Sophora davidii leaves along an elevation gradient in upper Minjiang River dry valley. Chin. J. Appl. Ecol. 17, 5–10 (2006).
Google Scholar 

45.
Li, F., Bao, W. & Zhu, L. Species diversity and spatial distribution of legumes in the dry valley of Minjiang River, SW China. J. Mt. Sci. 1, 76–84 (2010).
Google Scholar 

46.
Song, C. J. et al. Distribution patterns of shrubby N-fixers and non-N fixers in an arid valley in Southwest China: implications for ecological restoration. Ecol. Res. 25, 553–564. https://doi.org/10.1007/s11284-009-0685-3 (2010).
Article  Google Scholar 

47.
Liu, G. et al. Aboveground biomass of main shrubs in dry valley of Minjiang River. Acta Ecol. Sin. 23, 1757–1764. https://doi.org/10.1023/A:1022289509702 (2003).
Article  Google Scholar 

48.
Ellenberg, H. et al. Zeigerwerte von pflanzen in Mitteleuropa (1992).

49.
Shreve, F. Soil temperature as influenced by altitude and slope exposure. Ecology 5, 128–136. https://doi.org/10.2307/1929010 (1924).
Article  Google Scholar 

50.
Gutiérrez-Jurado, H. A. et al. On the observed ecohydrologic dynamics of a semiarid basin with aspect-delimited ecosystems. Water Resour. Res. 49, 8263–8284. https://doi.org/10.1002/2013wr014364 (2013).
Article  Google Scholar 

51.
Liu, G. et al. Distribution regulation of aboveground biomass of three main shrub types in the dry valley of Minjiang River. J. Mt. Sci. 21, 24–32 (2003).
Google Scholar 

52.
Pugnaire, F. I. & Luque, M. T. Changes in plant interactions along a gradient of environmental stress. Oikos 93, 42–49. https://doi.org/10.1034/j.1600-0706.2001.930104.x (2001).
Article  Google Scholar 

53.
Chase, J. M. & Leibold, M. A. Spatial scale dictates the productivity–biodiversity relationship. Nature 416, 427–430. https://doi.org/10.1038/416427a (2002).
ADS  CAS  Article  Google Scholar 

54.
Sarr, D. A., Hibbs, D. E. & Huston, M. A. A hierarchical perspective of plant diversity. Q. Rev. Biol. 80, 187–212. https://doi.org/10.1086/433058 (2005).
Article  PubMed  Google Scholar 

55.
Myers, N., Mittermeier, R. A., Mittermeier, C. G., Da Fonseca, G. A. & Kent, J. Biodiversity hotspots for conservation priorities. Nature 403, 853. https://doi.org/10.1038/35002501 (2000).
ADS  CAS  Article  PubMed  Google Scholar 

56.
Ricklefs, R. E. Community diversity: relative roles of local and regional processes. Science 235, 167–171. https://doi.org/10.1126/science.235.4785.167 (1987).
ADS  CAS  Article  PubMed  Google Scholar 

57.
Pang, X. Y., Bao, W. K. & Ning, W. U. Reasons of dry valley climate characteristic and its formation reason in upstream of Minjiang River. Resour. Environ. Yangtze Basin 17, 46–53 (2008).
Google Scholar 

58.
Minchin, P. R. An evaluation of the relative robustness of techniques for ecological ordination. Vegetatio 69, 89–107. https://doi.org/10.1007/BF00038690 (1987).
Article  Google Scholar 

59.
Clarke, K. R. Non-parametric multivariate analyses of changes in community structure. Aust. J. Ecol. 18, 117–143. https://doi.org/10.1111/j.1442-9993.1993.tb00438.x (1993).
Article  Google Scholar 

60.
Oksanen, J. et al.vegan: Community Ecology Package. R package version 2.5-2. 2018 (2018).

61.
Wang, Y. J., Huang, C. D., Zhang, J., Yang, W. Q. & Wang, X. S. Species Diversity, biomass and their relationship of shrubberies in an arid valley of the Minjiang River. Arid Zone Res. 27, 567–572. https://doi.org/10.3724/SP.J.1077.2010.01263 (2010).
CAS  Article  Google Scholar 

62.
Gotelli, N. J. & Colwell, R. K. Estimating species richness. In: Biological Diversity: Frontiers in Measurement and Assessment, Vol. 12 (eds Magurran, A. & McGill, B.) 39–54 (Oxford University Press, Oxford, 2011).
Google Scholar 

63.
Colwell, R. K. & Coddington, J. A. Estimating terrestrial biodiversity through extrapolation. Philos. Trans. R. Soc. Lond. B Biol. Sci. 345, 101–118 (1994).
ADS  CAS  Article  Google Scholar 

64.
Smith, E. P. & van Belle, G. Nonparametric estimation of species richness. Biometrics 40, 119–129 (1984).
Article  Google Scholar 

65.
Maurer, B. A. & McGill, B. J. In Biological Diversity: Frontiers in Measurement and Assessment (eds Magurran, A. E. & McGill, B. J.) 55–65 (Oxford University Press, Oxford, 2011).

66.
Fisher, R. A., Corbet, A. S. & Williams, C. B. The relation between the number of species and the number of individuals in a random sample of an animal population. J. Anim. Ecol. 12, 42–58 (1943).
Article  Google Scholar 

67.
Sørensen, T. A. A method of establishing groups of equal amplitude in plant sociology based on similarity of species content and its application to analyses of the vegetation on Danish commons. K. Dan. Vidensk. Selsk. Biol. Skr. 5, 1–34 (1948).
Google Scholar 

68.
Koleff, P., Gaston, K. J. & Lennon, J. J. Measuring beta diversity for presence-absence data. J. Anim. Ecol. 72, 367–382. https://doi.org/10.1046/j.1365-2656.2003.00710.x (2003).
Article  Google Scholar 

69.
Anderson, M. J. Distance-based tests for homogeneity of multivariate dispersions. Biometrics 62, 245–253. https://doi.org/10.1111/j.1541-0420.2005.00440.x (2006).
MathSciNet  Article  PubMed  MATH  Google Scholar 

70.
Anderson, M. J., Ellingsen, K. E. & McArdle, B. H. Multivariate dispersion as a measure of beta diversity. Ecol. Lett. 9, 683–693. https://doi.org/10.1111/j.1461-0248.2006.00926.x (2006).
Article  PubMed  Google Scholar 

71.
R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, Vienna, 2018). More

1.
Reynolds, M. P. (eds) Climate Change and Crop Production (CABI Publishing, 2010).
2.
Morales-Castilla, I. et al. Diversity buffers winegrowing regions from climate change losses. Proc. Natl Acad. Sci. USA 117, 2864–2869 (2020).
CAS  Article  Google Scholar 

3.
Nicotra, A. B. et al. Plant phenotypic plasticity in a changing climate. Trends Plant Sci. 15, 684–692 (2010).
CAS  Article  Google Scholar 

4.
Bellon, M. R., Hodson, D. & Hellin, J. Assessing the vulnerability of traditional maize seed systems in Mexico to climate change. Proc. Natl Acad. Sci. USA 108, 13432–13437 (2011).
CAS  Article  Google Scholar 

5.
Driedonks, N., Rieu, I. & Vriezen, W. H. Breeding for plant heat tolerance at vegetative and reproductive stages. Plant Reprod. 29, 67–79 (2016).
CAS  Article  Google Scholar 

6.
Mercer, K. L. & Perales, H. Evolutionary response of landraces to climate change in centers of crop diversity. Evol. Appl. 3, 480–493 (2010).
Article  Google Scholar 

7.
FAOSTAT (FAO, 2020); http://www.fao.org/faostat/en/#home

8.
Gibson, R., Mwanga, R. O. M., Namanda, S., Jeremiah, S. C. & Barker, I. Review of Sweetpotato Seed Systems in East and Southern Africa Working Paper 2009-1 (International Potato Center, 2009).

9.
Laurie, R. N., Laurie, S. M., Du Plooy, C. P., Finnie, J. F. & Van Staden, J. Yield of drought-stressed sweet potato in relation to canopy cover, stem length and stomatal conductance. J. Agric. Sci. 7, 201–214 (2015).
Google Scholar 

10.
Yang et al. High-throughput deep sequencing reveals the important role that microRNAs play in the salt response in sweet potato (Ipomoea batatas L.). BMC Genomics 21, 164 (2020).
CAS  Article  Google Scholar 

11.
Warren, J. F. Typhoons and droughts: food shortages and famine in the Philippines since the seventeenth century. Int. Rev. Environ. Hist. 4, 27–44 (2018).
Article  Google Scholar 

12.
Challinor, A. J. et al. A meta-analysis of crop yield under climate change and adaptation. Nat. Clim. Change 4, 287–291 (2014).
Article  Google Scholar 

13.
Roullier, C. et al. Disentangling the origins of cultivated sweet potato (Ipomoea batatas (L.) Lam.). PLoS ONE 8, e62707 (2013).
CAS  Article  Google Scholar 

14.
Kassali, R. Economics of sweet potato production. Int. J. Veg. Sci. 17, 313–321 (2011).
Article  Google Scholar 

15.
Omotobora, B. O., Adebola, P. O., Modise, D. M., Laurie, S. M. & Gerrano, A. S. Greenhouse and field evaluation of selected sweetpotato (Ipomoea batatas (L.) Lam.) accessions for drought tolerance in South Africa. Am. J. Plant Sci. 5, 3328–3339 (2014).
Article  Google Scholar 

16.
Woolfe, J. A. Sweetpotato: An Untapped Food Resource (Cambridge Univ. Press, 1992).

17.
Jayne, T. S., Villareal, M., Pingali, P. & Hemrich, G. Interactions between the Agricultural Sector and the HIV/AIDS Pandemic: Implications for Agricultural Policy ESA Working Paper No. 04-46 (FAO, 2004).

18.
Lebot, V. in Root and Tuber Crops (ed. Bradshaw, J. E.) 97–125 (Springer, 2010).

19.
Zhao et al. Temperature increase reduces global yields of major crops in four independent estimates. Proc. Natl Acad. Sci. USA 114, 9326–9331 (2017).
CAS  Article  Google Scholar 

20.
Lobell, D. B., Schlenker, W. & Costa-Roberts, J. Climate trends and global crop production since 1980. Science 333, 616–620 (2011).
CAS  Article  Google Scholar 

21.
Nangombe et al. Record-breaking climate extremes in Africa under stabilized 1.5 °C and 2 °C global warming scenarios. Nat. Clim. Change 8, 375–380 (2018).
Article  Google Scholar 

22.
Seneviratne, S. I. et al. in Special Report on Managing the Risks of Extreme Events and Disasters to Advance Climate Change Adaptation (eds Field, C. B. et al.) 109–230 (Cambridge Univ. Press, 2012).

23.
O’sullivan et al. Thermal limits of leaf metabolism across biomes. Glob. Change Biol. 23, 209–223 (2017).
Article  Google Scholar 

24.
Tack, J., Lingenfelser, J. & Jagadish, S. K. Disaggregating sorghum yield reductions under warming scenarios exposes narrow genetic diversity in US breeding programs. Proc. Natl Acad. Sci. USA 114, 9296–9301 (2017).
CAS  Article  Google Scholar 

25.
Araújo et al. Heat freezes niche evolution. Ecol. Lett. 16, 1206–1219 (2013).
Article  Google Scholar 

26.
Singh, K. D. & Mandal, R. C. Performance of coleus and sweet potato in relation to seasonal variations, time of planting. J. Root Crops 2, 17–22 (1976).
Google Scholar 

27.
Gajanayake, B. R., Reddy, K., Shankle, M. W., Arancibia, R. A. & Villordon, A. Quantifying storage root initiation, growth, and developmental responses of sweetpotato to early season temperature. Agr. J. 106, 1795–1804 (2014).
Article  Google Scholar 

28.
Boeck, H. J. D., Velde, H. V. D., Groote, T. D. & Nijs, I. Ideas and perspectives: heat stress: more than hot air. Biogeosciences 13, 5821–5825 (2016).
Article  Google Scholar 

29.
Wahid, A., Gelani, S., Ashraf, M. & Foolad, M. Heat tolerance in plants: an overview. Environ. Exp. Bot. 61, 199–223 (2007).
Article  Google Scholar 

30.
Bita, C. & Gerats, T. Plant tolerance to high temperature in a changing environment: scientific fundamentals and production of heat stress-tolerant crops. Front. Plant Sci. 4, 273 (2013).
Article  Google Scholar 

31.
Jones, H. G. Plants and Microclimate: A Quantitative Approach to Environmental Plant Physiology (Cambridge Univ. Press, 2013).

32.
Dong, N., Prentice, I. C., Harrison, S. P., Song, Q. H. & Zhang, Y. P. Biophysical homoeostasis of leaf temperature: a neglected process for vegetation and land-surface modelling. Glob. Ecol. Biogeogr. 26, 998–1007 (2017).
Article  Google Scholar 

33.
Reynolds, M. P. et al. Evaluating physiological traits to complement empirical selection for wheat in warm environments. Euphytica 100, 84–95 (1998).
Article  Google Scholar 

34.
Park, S. et al. Orange protein has a role in phytoene synthase stabilization in sweetpotato. Sci. Rep. 6, 33563 (2016).
CAS  Article  Google Scholar 

35.
Pironon, S. et al. Potential adaptive strategies for 29 sub-Saharan crops under future climate change. Nat. Clim. Change 9, 758–763 (2019).
Article  Google Scholar 

36.
Herrera, J. M. et al. Lessons from 20 years of studies of wheat genotypes in multiple environments and under contrasting production systems. Front. Plant Sci. 10, 1745 (2020).
Article  Google Scholar 

37.
Khoury, C. K. et al. Distributions, ex situ conservation priorities, and genetic resource potential of crop wild relatives of sweetpotato [Ipomoea batatas (L.) Lam., I. series Batatas]. Front. Plant Sci. 6, 251 (2015).
Article  Google Scholar 

38.
Alwang, J. et al. Pathways from research on improved staple crop germplasm to poverty reduction for smallholder farmers. Agric. Sys. 172, 16–27 (2019).
Article  Google Scholar 

39.
Pilling, D., Bélanger, J. & Hoffmann, I. Declining biodiversity for food and agriculture needs urgent global action. Nat. Food 1, 144–147 (2020).
Article  Google Scholar 

40.
Renard, D. & Tilman, D. National food production stabilized by crop diversity. Nature 571, 257–260 (2019).
CAS  Article  Google Scholar 

41.
Beck, H. et al. Present and future Köppen-Geiger climate classification maps at 1-km resolution. Sci. Data 5, 180214 (2018).
Article  Google Scholar 

42.
Patterson, H. D. & Williams, E. R. A new class of resolvable incomplete block designs. Biometrika 63, 83–92 (1976).
Article  Google Scholar 

43.
Kumar, A. et al. Improving the efficiency of wheat breeding experiments using alpha lattice design over randomized complete block design. Cereal Res. Commun. 48, 95–101 (2020).
Article  Google Scholar 

44.
Khan, M. et al. Comparative efficiency of alpha lattice design and complete randomized block design in wheat, maize and potato field trials. J. Res. Dev. Manag. 11, 115–118 (2015).
Google Scholar 

45.
Faye, E., Rebaudo, F., Yánez-Cajo, D., Cauvy-Fraunié, S. & Dangles, O. A toolbox for studying thermal heterogeneity across spatial scales: from unmanned aerial vehicle imagery to landscape metrics. Methods Ecol. Evol. 7, 437–446 (2016).
Article  Google Scholar 

46.
Liaw, A. & Wiener, M. Classification and regression by randomForest. R News 2, 18–22 (2002).
Google Scholar 

47.
Fick, S. E. & Hijmans, R. J. WorldClim 2: new 1-km spatial resolution climate surfaces for global land areas. Int. J. Clim. 37, 4302–4315 (2017).
Article  Google Scholar 

48.
Knox, J., Hess, T., Daccache, A. & Wheeler, T. Climate change impacts on crop productivity in Africa and South Asia. Environ. Res. Lett. 7, 034032 (2012).
Article  Google Scholar 

49.
Adoption of the Paris Agreement FCC/CP/2015/L.9/Rev.1 (UNFCCC, 2015). More

1.
Klicka, J. & Zink, R. M. Pleistocene effects on North American songbird evolution. Proc. R. Soc. B Biol. Sci. 266, 695–700 (1999).
Article  Google Scholar 
2.
Johnson, N. K. & Cicero, C. New mitochondrial DNA data affirm the importance of pleistocene speciation in North American birds. Evolution 58, 1122–1130 (2004).
PubMed  Article  Google Scholar 

3.
Avise, J. C. & Walker, D. Pleistocene phylogeographic effects on avian populations and the speciation process. Proc. R. Soc. B Biol. Sci. 265, 457–463 (1998).
CAS  Article  Google Scholar 

4.
April, J., Hanner, R. H., Dion-Côté, A.-M. & Bernatchez, L. Glacial cycles as an allopatric speciation pump in north-eastern American freshwater fishes. Mol. Ecol. 22, 409–422 (2013).
CAS  PubMed  Article  Google Scholar 

5.
Ribera, I. & Vogler, A. P. Speciation of Iberian diving beetles in Pleistocene refugia (Coleoptera, Dytiscidae). Mol. Ecol. 13, 179–193 (2004).
CAS  PubMed  Article  Google Scholar 

6.
Hewitt, G. Some genetic consequences of ice ages, and their role in divergence and speciation. Biol. J. Linn. Soc. 58, 247–276 (1996).
Article  Google Scholar 

7.
Hewitt, G. M. Genetic consequences of climatic oscillations in the Quaternary. . Philos. Trans. R. Soc. Lond., B Biol. Sci. 359, 183–195 (2004).
CAS  Article  Google Scholar 

8.
Hewitt, G. Post-glacial re-colonization of European biota. Biol. J. Linn. Soc. 68, 87–112 (1999).
Article  Google Scholar 

9.
Skrzecz, I. & Bulka, M. Insect assemblages in Norway spruce (Picea abies (L.) Karst.) stumps in the Eastern Sudetes. Folia For. Pol. Ser. A. 52, 98–107 (2010).
Google Scholar 

10.
Röder, J. et al. Arthropod species richness in the Norway Spruce (Picea abies (L.) Karst.) canopy along an elevation gradient. For. Ecol. Manage. 259, 1513–1521 (2010).
Article  Google Scholar 

11.
Tollefsrud, M. M. et al. Genetic consequences of glacial survival and postglacial colonization in Norway spruce: combined analysis of mitochondrial DNA and fossil pollen. Mol. Ecol. 17, 4134–4150 (2008).
CAS  PubMed  Article  Google Scholar 

12.
Latałowa, M. & van der Knaap, W. O. Late quaternary expansion of Norway spruce Picea abies (L.) Karst. in Europe according to pollen data. Quat. Sci. Rev. 25, 2780–2805 (2006).
ADS  Article  Google Scholar 

13.
Gałka, M. & Tobolski, K. Macrofossil evidence of early Holocene presence of Picea abies (Norway spruce) in NE Poland. Ann. Bot. Fenn. 3847, 129–141 (2013).
Article  Google Scholar 

14.
Dering, M. & Lewandowski, A. Finding the meeting zone: Where have the northern and southern ranges of Norway spruce overlapped?. For. Ecol. Manage. 259, 229–235 (2009).
Article  Google Scholar 

15.
Giesecke, T. & Bennett, K. D. The Holocene spread of Picea abies (L.) Karst. in Fennoscandia and adjacent areas. J. Biogeogr. 31, 1523–1548 (2004).
Article  Google Scholar 

16.
Taberlet, P., Fumagalli, L., Wust-Saucy, A. G. & Cosson, J. F. Comparative phylogeography and postglacial colonization routes in Europe. Mol. Ecol. 7, 453–464 (1998).
CAS  PubMed  Article  Google Scholar 

17.
Zając, A. & Zając, M. Atlas rozmieszczenia roślin naczyniowych w Polsce. (ed. Zając, A. & Zając, M) (Pracownia Chorologii Komputerowej Instytutu Botaniki Uniwersytetu Jagiellońskiego, 2001).

18.
Sallé, A., Arthofer, W., Lieutier, F., Stauffer, C. & Kerdelhué, C. Phylogeography of a host-specific insect: Genetic structure of Ips typographus in Europe does not reflect past fragmentation of its host. Biol. J. Linn. Soc. 90, 239–246 (2007).
Article  Google Scholar 

19.
Stauffer, C., Lakatos, F. & Hewitt, G. M. Phylogeography and postglacial colonization routes of Ips typographus L. (Coleoptera, Scolytidae). Mol. Ecol. 8, 763–773 (1999).
Article  Google Scholar 

20.
Bertheau, C. et al. Divergent evolutionary histories of two sympatric spruce bark beetle species. Mol. Ecol. 22, 3318–3332 (2013).
CAS  PubMed  Article  Google Scholar 

21.
Mayer, F. et al. Comparative multilocus phylogeography of two Palaearctic spruce bark beetles: influence of contrasting ecological strategies on genetic variation. Mol. Ecol. 24, 1292–1310 (2015).
PubMed  Article  Google Scholar 

22.
Schebeck, M. et al. Pleistocene climate cycling and host plant association shaped the demographic history of the bark beetle Pityogenes chalcographus. Sci. Rep. 8, 14207 (2018).
ADS  PubMed  PubMed Central  Article  CAS  Google Scholar 

23.
Forsse, E. & Solbreck, C. Migration in the bark beetle Ips typographus L: duration, timing and height of flight. Zeitschrift für Angewandte Entomologie 100, 47–57 (1985).
Article  Google Scholar 

24.
Putz, J., Vorwagner, E. M. & Hoch, G. Flight performance of Monochamus sartor and Monochamus sutor, potential vectors of the pine wood nematode. For. J. 62, 195–201 (2016).
Google Scholar 

25.
Kawai, M. et al. Genetic Structure of Pine Sawyer Monochamus alternatus (Coleoptera: Cerambycidae) Populations in Northeast Asia: Consequences of the Spread of Pine Wilt Disease. Environ. Entomol. 35, 569–579 (2006).
Article  Google Scholar 

26.
Shoda-Kagaya, E. Genetic differentiation of the pine wilt disease vector Monochamus alternatus (Coleoptera: Cerambycidae) over a mountain range – revealed from microsatellite DNA markers. Bull. Entomol. Res. 97, 167–174 (2007).
CAS  PubMed  Article  Google Scholar 

27.
Haran, J., Rousselet, J., Tellez, D., Roques, A. & Roux, G. Phylogeography of Monochamus galloprovincialis, the European vector of the pinewood nematode. J. Pest Sci. 2004(91), 247–257 (2018).
Article  Google Scholar 

28.
Koutroumpa, F. A., Rougon, D., Bertheau, C., Lieutier, F. & Roux-Morabito, G. Evolutionary relationships within European Monochamus (Coleoptera: Cerambycidae) highlight the role of altitude in species delineation. Biol. J. Linn. Soc. 109, 354–376 (2013).
Article  Google Scholar 

29.
Plewa, R. et al. Morphology, genetics and Wolbachia endosymbionts support distinctiveness of Monochamus sartor sartor and M. s. urussovii (Coleoptera: Cerambycidae). Arthropod Syst. Phylogeny 76, 123–135 (2018).
Google Scholar 

30.
Ballard, J. W. O. & Whitlock, M. C. The incomplete natural history of mitochondria. Mol. Ecol. 13, 729–744 (2004).
PubMed  Article  Google Scholar 

31.
Toews, D. P. L. & Brelsford, A. The biogeography of mitochondrial and nuclear discordance in animals. Mol. Ecol. 21, 3907–3930 (2012).
CAS  PubMed  Article  Google Scholar 

32.
Bazin, E. Population size does not influence mitochondrial genetic diversity in animals. Science 312, 570–572 (2006).
ADS  CAS  PubMed  Article  Google Scholar 

33.
Plavil’shhikov, N. Zhuki-drovoseki, chast’3. Podsemejstvo Lamiinae, ch. 1. (ed. Plavil’shhikov, N.) 510–514 (Fauna SSSR. Zhestkokrylye, 1958).

34.
Sama, G. Atlas of the Cerambycidae of Europe and the Mediterranean Area (ed. Sama, G.) 98–99 (Vit Kabourek, 2002).

35.
Sama, G. & Löbl, L. Cerambycidae, Western Palaearctic taxa.In:Catalogue of Palaearctic Coleoptera. Volume 6. Chrysomeloidea (ed. Löbl, I. & Smetana, A.) 281–283 (Apollo Books, 2010).

36.
Wallin, H., Schroeder, M. & Kvamme, T. A review of the European species of Monochamus Dejean, 1821 (Coleoptera, Cerambycidae) -with a description of the genitalia characters. Nor. J. Entomol. 60, 11–38 (2013).
Google Scholar 

37.
Rossa, R., Goczał, J. & Tofilski, A. Within and between-species variation of wing venation in genus Monochamus (Coleoptera: Cerambycidae). J. Insect Sci. 16, 5 (2016).
PubMed  PubMed Central  Article  Google Scholar 

38.
Ravazzi, C. Late quaternary history of spruce in southern Europe. Rev. Palaeobot. Palynol. 120, 131–177 (2002).
Article  Google Scholar 

39.
Terhürne-Berson, R. Changing distribution patterns of selected conifers in the Quaternary of Europe caused by climatic variations.https://d-nb.info/975944576/34 (2005).

40.
Ralska-Jasiewiczowa, M Latałowa, M. et al. Late Glacial and Holocene history of vegetation in Poland based on isopollen maps(ed. Ralska-Jasiewiczowa, M Latałowa, M. et al.) 147–159 (W. Szafer Institute of Botany, Polish Academy of Sciences, 2004).

41.
Adams, D. C., Rohlf, F. J. & Slice, D. E. Geometric morphometrics: ten years of progress following the ‘revolution’. Ital. J. Zool. 71, 5–16 (2004).
Article  Google Scholar 

42.
Zelditch, M. L., Swiderski, D. L., Sheets, H. D. & Fink, W. L. Geometric Morphometrics for Biologist. A Primer(ed. Zelditch, M. L., Swiderski, D. L., Sheets, H. D. & Fink, W. L.). 263–490 (2004).

43.
Oleksa, A. & Tofilski, A. Wing geometric morphometrics and microsatellite analysis provide similar discrimination of honey bee subspecies. Apidologie 46, 49–60 (2015).
Article  Google Scholar 

44.
Hurtado-Burillo, M. et al. A geometric morphometric and microsatellite analyses of Scaptotrigona mexicana and S. pectoralis (Apidae: Meliponini) sheds light on the biodiversity of Mesoamerican stingless bees. J. Insect Conserv. 5, 753–763 (2016).
Article  Google Scholar 

45.
Haack, R. A. Exotic bark- and wood-boring Coleoptera in the United States: recent establishments and interceptions. Can. J. For. Res. 36, 269–288 (2006).
Article  Google Scholar 

46.
Shoda, E., Kubota, K. & Makihara, H. Geographical structuring of mitochondrial DNA in Semanotus japonicus (Coleoptera: Cerambycidae). Appl. Entomol. Zool. 38, 339–345 (2003).
CAS  Article  Google Scholar 

47.
Shoda, E., Kubota, K. & Makihara, H. Geographical structure of morphological characters in Semanotus japonicus (Coleoptera: Cerambycidae) in Japan. Appl. Entomol. Zool. 38, 369–377 (2003).
Article  Google Scholar 

48.
Rossa, R., Goczał, J., Pawliczek, B. & Ohbayashi, N. Hind wing variation in Leptura annularis complex among European and Asiatic populations (Coleoptera, Cerambycidae). Zookeys 724, 31–42 (2017).
Article  Google Scholar 

49.
Kolk, A. & Starzyk, J. Atlas szkodliwych owadów lesnych (The atlas of harmful forest insects). 268–488 (Multico,1996).

50.
Cherepanov, A. Usachi Severnoj Azii (Lamiinae) (ed. Cherepanov, A.) 97–102 (Nauka Publishers, 1983).

51.
Danilevsky, M. L. A check list of the longicorn beetles (Cerambycoidea) of Russia.https://www.zin.ru/Animalia/Coleoptera/doc/List_of_Russia.doc (2015).

52.
Haran, J. & Roux-Morabito, G. Development of 12 microsatellites loci for the longhorn beetle Monochamus galloprovincialis (Coleoptera Cerambycidae), vector of the Pine Wood Nematode in Europe. Conserv. Genet. Res. 6, 975–977 (2014).
Article  Google Scholar 

53.
Chybicki, I. J. & Burczyk, J. Simultaneous estimation of null alleles and inbreeding coefficients. J. Hered. 100, 106–113 (2009).
CAS  PubMed  Article  Google Scholar 

54.
Guo, S. W. & Thompson, E. A. Performing the exact test of Hardy-Weinberg proportion for multiple alleles. Biometrics 1, 361–372 (1992).
MATH  Article  Google Scholar 

55.
Paradis, E. Pegas: An R package for population genetics with an integrated-modular approach. Bioinformatics 26, 419–420 (2010).
CAS  PubMed  Article  Google Scholar 

56.
Agapow, P. M. & Burt, A. Indices of multilocus linkage disequilibrium. Mol. Ecol. Notes 1, 101–102 (2001).
CAS  Article  Google Scholar 

57.
Kamvar, Z. N., Tabima, J. F. & Gr̈unwald, N. J. Poppr, ,. An R package for genetic analysis of populations with clonal, partially clonal, and/or sexual reproduction. PeerJ e281, 1–14 (2014).
Google Scholar 

58.
Storey, J. D., Bass, A. J., Dabney, A. & Robinson, D. qvalue: Q-value estimation for false discovery rate control. R package version 2.14.1. https://qvalue.princeton.edu (2019).

59.
Goudet, J. & Jombart, T. hierfstat: estimation and tests of hierarchical F-statistics. https://github.com/jgx65/hierfstat (2015).

60.
Hsieh, T. C., Ma, K. H. & Chao, A. iNEXT: iNterpolation and EXTrapolation for species diversity. R package version 2.0.19. https://chao.stat.nthu.edu.tw/blog/software-download/ (2019).

61.
Chao, A. et al. Rarefaction and extrapolation with Hill numbers: a framework for sampling and estimation in species diversity studies. Ecol. Monogr. 84, 45–67 (2014).
Article  Google Scholar 

62.
Evanno, G., Regnaut, S. & Goudet, J. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol. Ecol. 14, 2611–2620 (2005).
CAS  PubMed  PubMed Central  Article  Google Scholar 

63.
Jombart, T. et al. Discriminant analysis of principal components: a new method for the analysis of genetically structured populations. BMC Genet. 11, 94 (2010).
PubMed  PubMed Central  Article  Google Scholar 

64.
Jombart, T. Adegenet: a R package for the multivariate analysis of genetic markers. Bioinformatics 24, 1403–1405 (2008).
CAS  PubMed  Article  Google Scholar 

65.
Pastore, M. Overlapping: a R package for estimating overlapping in empirical distributions. J. Open Source Softw. 3, 1023 (2018).
ADS  Article  Google Scholar 

66.
Cornuet, J.-M. et al. DIYABC v20: a software to make approximate Bayesian computation inferences about population history using single nucleotide polymorphism, DNA sequence and microsatellite data. Bioinformatics 30, 1187–1189 (2014).
CAS  PubMed  Article  Google Scholar 

67.
Goczał, J., Rossa, R. & Tofilski, A. Intersexual and intrasexual patterns of horn size and shape variation in the European rhinoceros beetle: quantifying the shape of weapons. Biol. J. Linn. Soc. 127, 34–43 (2019).
Article  Google Scholar 

68.
Rohlf, F. & Slice, D. Extensions of the procrustes method for the optimal superimposition of landmarks. Syst. Zool. 39, 40–59 (1990).
Article  Google Scholar 

69.
Cattell, R. B. The scree test for the number of factors. Multivariate Behav. Res. 1, 245–276 (1966).
CAS  PubMed  Article  Google Scholar 

70.
TIBCO Software Inc. Statistica (data analysis software system) version 13. https://www.statsoft.pl/statistica-i-tibco-software/ (2017).

71.
Schliep, K. P. phangorn: phylogenetic analysis in R. Bioinformatics 27, 592–593 (2011).
CAS  PubMed  PubMed Central  Article  Google Scholar 

72.
R Core Team. R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing. https://www.r-project.org/ (2015). More

1.
Mante, V., Sussillo, D., Shenoy, K. V. & Newsome, W. T. Context-dependent computation by recurrent dynamics in prefrontal cortex. Nature 503, 78–84 (2013).
CAS  PubMed  PubMed Central  Article  Google Scholar 
2.
Raposo, D., Kaufman, M. T. & Churchland, A. K. A category-free neural population supports evolving demands during decision-making. Nat. Neurosci. 17, 1784–1792 (2014).
CAS  PubMed  PubMed Central  Article  Google Scholar 

3.
Roy, J. E., Buschman, T. J. & Miller, E. K. PFC neurons reflect categorical decisions about ambiguous stimuli. J. Cogn. Neurosci. 26, 1283–1291 (2014).
PubMed  PubMed Central  Article  Google Scholar 

4.
Siegel, M., Buschman, T. J. & Miller, E. K. Cortical information flow during flexible sensorimotor decisions. Science 348, 1352–1355 (2015).
CAS  PubMed  PubMed Central  Article  Google Scholar 

5.
Katz, L. N., Yates, J. L., Pillow, J. W. & Huk, A. C. Dissociated functional significance of decision-related activity in the primate dorsal stream. Nature 535, 285–288 (2016).
CAS  PubMed  PubMed Central  Article  Google Scholar 

6.
Hanks, T. D. et al. Distinct relationships of parietal and prefrontal cortices to evidence accumulation. Nature 520, 220–223 (2015).
CAS  PubMed  PubMed Central  Article  Google Scholar 

7.
Erlich, J. C., Brunton, B. W., Duan, C. A., Hanks, T. D. & Brody, C. D. Distinct effects of prefrontal and parietal cortex inactivations on an accumulation of evidence task in the rat. eLife 4, e05457 (2015).
PubMed Central  Article  PubMed  Google Scholar 

8.
Goard, M. J., Pho, G. N., Woodson, J. & Sur, M. Distinct roles of visual, parietal, and frontal motor cortices in memory-guided sensorimotor decisions. eLife 5, e13764 (2016).
PubMed  PubMed Central  Article  Google Scholar 

9.
Bruce, C. J. & Goldberg, M. E. Primate frontal eye fields. I. Single neurons discharging before saccades. J. Neurophysiol. 53, 603–635 (1985).
CAS  PubMed  Article  Google Scholar 

10.
Kim, J. N. & Shadlen, M. N. Neural correlates of a decision in the dorsolateral prefrontal cortex of the macaque. Nat. Neurosci. 2, 176–185 (1999).
PubMed  Article  Google Scholar 

11.
Ding, L. & Gold, J. I. Neural correlates of perceptual decision making before, during, and after decision commitment in monkey frontal eye field. Cereb. Cortex 22, 1052–1067 (2011).
PubMed  PubMed Central  Article  Google Scholar 

12.
Stokes, M. G. et al. Dynamic coding for cognitive control in prefrontal cortex. Neuron 78, 364–375 (2013).
CAS  PubMed  PubMed Central  Article  Google Scholar 

13.
Rigotti, M. et al. The importance of mixed selectivity in complex cognitive tasks. Nature 497, 585–590 (2013).
CAS  PubMed  PubMed Central  Article  Google Scholar 

14.
Purcell, B. A. et al. Neurally constrained modeling of perceptual decision making. Psychol. Rev. 117, 1113–1143 (2010).
PubMed Central  Article  PubMed  Google Scholar 

15.
Heitz, R. P. & Schall, J. D. Neural mechanisms of speed-accuracy tradeoff. Neuron 76, 616–628 (2012).
CAS  PubMed  PubMed Central  Article  Google Scholar 

16.
Parthasarathy, A. et al. Mixed selectivity morphs population codes in prefrontal cortex. Nat. Neurosci. 20, 1770–1779 (2017).
CAS  PubMed  Article  Google Scholar 

17.
Beck, J. M. et al. Probabilistic population codes for Bayesian decision making. Neuron 60, 1142–1152 (2008).
CAS  PubMed  PubMed Central  Article  Google Scholar 

18.
Churchland, M. M. et al. Neural population dynamics during reaching. Nature 487, 51–56 (2012).
CAS  PubMed  PubMed Central  Article  Google Scholar 

19.
Kaufman, M. T., Churchland, M. M., Ryu, S. I. & Shenoy, K. V. Cortical activity in the null space: permitting preparation without movement. Nat. Neurosci. 17, 440–448 (2014).
CAS  PubMed  PubMed Central  Article  Google Scholar 

20.
Lara, A. H., Cunningham, J. P. & Churchland, M. M. Different population dynamics in the supplementary motor area and motor cortex during reaching. Nat. Commun. 9, 2754 (2018).
CAS  PubMed  PubMed Central  Article  Google Scholar 

21.
Akaike, H. A new look at the statistical model identification. IEEE Transactions on Automatic Control 19, 716–723 (1974).
Article  Google Scholar 

22.
Aoi, M. & Pillow, J. W. Model-based targeted dimensionality reduction for neuronal population data. Adv. Neural Inform. Process. Syst. 31, 6690–6699 (2018).
Google Scholar 

23.
Elsayed, G. F. & Cunningham, J. P. Structure in neural population recordings: an expected byproduct of simpler phenomena? Nat. Neurosci. 20, 1310–1318 (2017).
PubMed  PubMed Central  Article  CAS  Google Scholar 

24.
Rossi-Pool, R. et al. Decoding a decision process in the neuronal population of dorsal premotor cortex. Neuron 96, 1432–1446 (2017).
CAS  PubMed  Article  Google Scholar 

25.
Williamson, R. C. et al. Scaling properties of dimensionality reduction for neural populations and network models. PLoS Comput. Biol. 12, e1005141 (2016).
PubMed  PubMed Central  Article  CAS  Google Scholar 

26.
Cunningham, J. P. & Byron, M. Y. Dimensionality reduction for large-scale neural recordings. Nature Neurosci. 17, 1500–1509 (2014).
CAS  PubMed  Article  Google Scholar 

27.
Kobak, D. et al. Demixed principal component analysis of neural population data. eLife 5, e10989 (2016).
PubMed  PubMed Central  Article  CAS  Google Scholar 

28.
Machens, C. K. Demixing population activity in higher cortical areas. Front. Comput. Neurosci. 4, 126 (2010).
PubMed  PubMed Central  Article  Google Scholar 

29.
Machens, C. K., Romo, R. & Brody, C. D. Functional, but not anatomical, separation of ‘what’ and ‘when’ in prefrontal cortex. J. Neurosci. 30, 350–360 (2010).
CAS  PubMed  PubMed Central  Article  Google Scholar 

30.
Romo, R., Brody, C. D., Hernández, A. & Lemus, L. Neuronal correlates of parametric working memory in the prefrontal cortex. Nature 399, 470–473 (1999).
CAS  PubMed  Article  Google Scholar 

31.
Hernández, A., Zainos, A. & Romo, R. Temporal evolution of a decision-making process in medial premotor cortex. Neuron 33, 959–972 (2002).
PubMed  Article  Google Scholar 

32.
Romo, R., Hernández, A., Zainos, A., Lemus, L. & Brody, C. D. Neuronal correlates of decision-making in secondary somatosensory cortex. Nat. Neurosci. 5, 1217–1225 (2002).
CAS  PubMed  Article  Google Scholar 

33.
Romo, R., Hernández, A. & Zainos, A. Neuronal correlates of a perceptual decision in ventral premotor cortex. Neuron 41, 165–173 (2004).
CAS  PubMed  Article  Google Scholar 

34.
Elsayed, G. F., Lara, A. H., Kaufman, M. T., Churchland, M. M. & Cunningham, J. P. Reorganization between preparatory and movement population responses in motor cortex. Nat. Commun. 7, 13239 (2016).
CAS  PubMed  PubMed Central  Article  Google Scholar 

35.
Gold, J. I. & Shadlen, M. N. The influence of behavioral context on the representation of a perceptual decision in developing oculomotor commands. J. Neurosci. 23, 632–651 (2003).
CAS  PubMed  PubMed Central  Article  Google Scholar 

36.
Kiani, R., Hanks, T. D. & Shadlen, M. N. Bounded integration in parietal cortex underlies decisions even when viewing duration is dictated by the environment. J. Neurosci. 28, 3017–3029 (2008).
CAS  PubMed  PubMed Central  Article  Google Scholar 

37.
Meister, M. L. R., Hennig, J. A. & Huk, A. C. Signal multiplexing and single-neuron computations in lateral intraparietal area during decision-making. J. Neurosci. 33, 2254–2267 (2013).
CAS  PubMed  PubMed Central  Article  Google Scholar 

38.
Latimer, K. W., Yates, J. L., Meister, M. L. R., Huk, A. C. & Pillow, J. W. Single-trial spike trains in parietal cortex reveal discrete steps during decision-making. Science 349, 184–187 (2015).
CAS  PubMed  PubMed Central  Article  Google Scholar 

39.
Yates, J. L., Park, I. M., Katz, L. N., Pillow, J. W. & Huk, A. C. Functional dissection of signal and noise in mt and lip during decision-making. Nat. Neurosci. 20, 1285–1292 (2017).
PubMed  PubMed Central  Article  CAS  Google Scholar 

40.
Vernet, M., Quentin, R., Chanes, L., Mitsumasu, A. & Valero-Cabré, A. Frontal eye field, where art thou? anatomy, function, and non-invasive manipulation of frontal regions involved in eye movements and associated cognitive operations. Front. Integr. Neurosci. 8, 1–24 (2014).
Google Scholar 

41.
Bruce, C. J., Goldberg, M. E., Bushnell, M. C. & Stanton, G. B. Primate frontal eye fields. II. Physiological and anatomical correlates of electrically evoked eye movements. J. Neurophysiol. 54, 714–734 (1985).
CAS  PubMed  Article  Google Scholar 

42.
Rizzolatti, G., Riggio, L., Dascola, I. & Umiltá, C. Reorienting attention across the horizontal and vertical meridians: evidence in favor of a premotor theory of attention. Neuropsychologia 25, 31–40 (1987).
CAS  PubMed  Article  Google Scholar 

43.
Thompson, K. G., Biscoe, K. L. & Sato, T. R. Neuronal basis of covert spatial attention in the frontal eye field. J. Neurosci. 25, 9479–9487 (2005).
CAS  PubMed  PubMed Central  Article  Google Scholar 

44.
Schall, J. D. On the role of frontal eye field in guiding attention and saccades. Vision Res. 44, 1453–1467 (2004).
PubMed  Article  Google Scholar 

45.
Storey, J. D. A direct approach to false discovery rates. J. R. Stat. Soc. Series B Stat. Methodol. 64, 479–498 (2002).
Article  Google Scholar 

46.
Brody, C. D., Hernández, A., Zainos, A. & Romo, R. Timing and neural encoding of somatosensory parametric working memory in macaque prefrontal cortex. Cereb. Cortex 13, 1196–1207 (2003).
PubMed  Article  Google Scholar 

47.
Roitman, J. D. & Shadlen, M. N. Response of neurons in the lateral intraparietal area during a combined visual discrimination reaction time task. J. Neurosci. 22, 9475–9489 (2002).
PubMed Central  Article  PubMed  Google Scholar 

48.
Lawrence, N. Probabilistic non-linear principal component analysis with Gaussian process latent variable models. J. Mach. Learn. Res. 6, 1783–1816 (2005).
Google Scholar 

49.
Cunningham, J. P. & Ghahramani, Z. Linear dimensionality reduction: survey, insights, and generalizations. J. Mach. Learn. Res. 16, 2859–2900 (2015).

50.
Mackevicius, E. L. et al. Unsupervised discovery of temporal sequences in high-dimensional datasets, with applications to neuroscience. eLife 8, e38471 (2019).
PubMed  PubMed Central  Article  Google Scholar  More

1.
van Belkum A, Soriaga LB, LaFave MC, Akella S, Veyrieras JB, Barbu EM, et al. Phylogenetic distribution of CRISPR-Cas systems in antibiotic-resistant Pseudomonas aeruginosa. MBio. 2015;6:e01796–15.
2.
Makarova KS, Wolf YI, Iranzo J, Shmakov SA, Alkhnbashi OS, Brouns SJ, et al. Evolutionary classification of CRISPR-Cas systems: a burst of class 2 and derived variants. Nat Rev Microbiol. 2019;18:1–17.

3.
Portillo MC, Gonzalez JM. CRISPR elements in the Thermococcales: evidence for associated horizontal gene transfer in Pyrococcus furiosus. J Appl Genet. 2009;50:421–30.
CAS  PubMed  Article  Google Scholar 

4.
Watson BN, Staals RH, Fineran, PC. CRISPR-Cas-mediated phage resistance enhances horizontal gene transfer by transduction. MBio. 2018;9:e02406–17.

5.
Varble A, Meaden S, Barrangou R, Westra ER, Marraffini LA. Recombination between phages and CRISPR-cas loci facilitates horizontal gene transfer in staphylococci. Nat Microbiol. 2019;4:956–63.
CAS  PubMed  PubMed Central  Article  Google Scholar 

6.
Levin BR. Nasty viruses, costly plasmids, population dynamics, and the conditions for establishing and maintaining CRISPR-mediated adaptive immunity in bacteria. PLoS Genet. 2010;6:e1001171.
PubMed  PubMed Central  Article  CAS  Google Scholar 

7.
van Houte S, Buckling A, Westra ER. Evolutionary ecology of prokaryotic immune mechanisms. Microbiol Mol Biol Rev. 2016;80:745–63.
PubMed  PubMed Central  Article  Google Scholar 

8.
Rollie C, Chevallereau A, Watson BN, Chyou TY, Fradet O, McLeod I, et al. Targeting of temperate phages drives loss of type I CRISPR-Cas systems. Nature. 2020;578:149–53.
CAS  PubMed  PubMed Central  Article  Google Scholar 

9.
Westra ER, van Houte S, Oyesiku-Blakemore S, Makin B, Broniewski JM, Best A, et al. Parasite exposure drives selective evolution of constitutive versus inducible defense. Curr Biol. 2015;25:1043–9.
CAS  PubMed  Article  Google Scholar 

10.
Vale PF, Lafforgue G, Gatchitch F, Gardan R, Moineau S, Gandon S. Costs of CRISPR-Cas-mediated resistance in Streptococcus thermophilus. Proc R Soc B: Biol Sci. 2015;282:20151270.
Article  CAS  Google Scholar 

11.
van Houte S, Ekroth AK, Broniewski JM, Chabas H, Ashby B, Bondy-Denomy J, et al. The diversity-generating benefits of a prokaryotic adaptive immune system. Nature. 2016;532:385–8.
PubMed  PubMed Central  Article  CAS  Google Scholar 

12.
Stern A, Keren L, Wurtzel O, Amitai G, Sorek R. Self-targeting by CRISPR: gene regulation or autoimmunity? Trends Genet. 2010;26:335–40.
CAS  PubMed  PubMed Central  Article  Google Scholar 

13.
Levy A, Goren MG, Yosef I, Auster O, Manor M, Amitai G, et al. CRISPR adaptation biases explain preference for acquisition of foreign DNA. Nature. 2015;520:505–10.
CAS  PubMed  PubMed Central  Article  Google Scholar 

14.
Vercoe RB, Chang JT, Dy RL, Taylor C, Gristwood T, Clulow JS, et al. Cytotoxic chromosomal targeting by CRISPR/Cas systems can reshape bacterial genomes and expel or remodel pathogenicity islands. PLoS Genet. 2013;9:e1003454.
CAS  PubMed  PubMed Central  Article  Google Scholar 

15.
Staals RH, Jackson SA, Biswas A, Brouns SJ, Brown CM, Fineran PC. Interference-driven spacer acquisition is dominant over naive and primed adaptation in a native CRISPR–Cas system. Nat Commun. 2016;7:1–13.
Article  CAS  Google Scholar 

16.
Weissman JL, Stoltzfus A, Westra ER, Johnson PL. Avoidance of Self during CRISPR immunization. Trends Microbiol. 2020;28:543–53.

17.
Koonin EV. Open questions: CRISPR biology. BMC Biol. 2018;16:1–3.
Article  CAS  Google Scholar 

18.
Patterson AG, Yevstigneyeva MS, Fineran PC. Regulation of CRISPR–Cas adaptive immune systems. Curr Opin Microbiol. 2017;37:1–7.
CAS  PubMed  Article  Google Scholar 

19.
Quax TE, Voet M, Sismeiro O, Dillies MA, Jagla B, Coppee JY, et al. Massive activation of archaeal defense genes during viral infection. J Virol. 2013;87:8419–28.
CAS  PubMed  PubMed Central  Article  Google Scholar 

20.
Young JC, Dill BD, Pan C, Hettich RL, Banfield JF, Shah M, et al. Phage-induced expression of CRISPR-associated proteins is revealed by shotgun proteomics in Streptococcus thermophilus. PloS ONE. 2012;7:e38077.
CAS  PubMed  PubMed Central  Article  Google Scholar 

21.
Agari Y, Sakamoto K, Tamakoshi M, Oshima T, Kuramitsu S, Shinkai A. Transcription profile of Thermus thermophilus CRISPR systems after phage infection. J Mol Biol. 2010;395:270–81.
CAS  PubMed  Article  Google Scholar 

22.
Chabas H, van Houte S, Høyland-Kroghsbo NM, Buckling A, Westra ER. Immigration of susceptible hosts triggers the evolution of alternative parasite defence strategies. Proc R Soc B: Biol Sci. 2016;283:20160721.
Article  Google Scholar 

23.
Pawluk A, Staals RH, Taylor C, Watson BN, Saha S, Fineran PC, et al. Inactivation of CRISPR-Cas systems by anti-CRISPR proteins in diverse bacterial species. Nat Microbiol. 2016;1:1–6.
Article  CAS  Google Scholar 

24.
Bondy-Denomy J, Garcia B, Strum S, Du M, Rollins MF, Hidalgo-Reyes Y, et al. Multiple mechanisms for CRISPR-Cas inhibition by anti-CRISPR proteins. Nature. 2015;526:136–9.
CAS  PubMed  PubMed Central  Article  Google Scholar 

25.
Stanley SY, Borges AL, Chen KH, Swaney DL, Krogan NJ, Bondy-Denomy J, et al. Anti-CRISPR-associated proteins are crucial repressors of anti-CRISPR transcription. Cell. 2019;178:1452–64.
CAS  PubMed  PubMed Central  Article  Google Scholar 

26.
Cady KC, Bondy-Denomy J, Heussler GE, Davidson AR, O’Toole GA. The CRISPR/Cas adaptive immune system of Pseudomonas aeruginosa mediates resistance to naturally occurring and engineered phages. J Bacteriol. 2012;194:5728–38.
CAS  PubMed  PubMed Central  Article  Google Scholar 

27.
Chevallereau A, Meaden S, Fradet O, Landsberger M, Maestri A, Biswas A, et al. Exploitation of the cooperative behaviors of anti-CRISPR phages. Cell Host Microbe. 2020;27:189–98.
CAS  PubMed  PubMed Central  Article  Google Scholar 

28.
Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 2011;17:10–12.
Article  Google Scholar 

29.
Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol. 2012;19:455–77.
CAS  PubMed  PubMed Central  Article  Google Scholar 

30.
Zhang J, Kobert K, Flouri T, Stamatakis A. PEAR: a fast and accurate Illumina Paired-End reAd mergeR. Bioinformatics. 2014;30:614–20.
CAS  PubMed  Article  Google Scholar 

31.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10.
CAS  PubMed  Article  Google Scholar 

32.
Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7:335–6.
CAS  PubMed  PubMed Central  Article  Google Scholar 

33.
Mahé F, Rognes T, Quince C, de Vargas C, Dunthorn M. Swarm v2: highly-scalable and high-resolution amplicon clustering. PeerJ. 2015;3:e1420.
PubMed  PubMed Central  Article  Google Scholar 

34.
Chao A. Nonparametric estimation of the number of classes in a population. Scand J Stat. 1984;11:265–70.

35.
Biswas A, Staals RH, Morales SE, Fineran PC, Brown CM. CRISPRDetect: a flexible algorithm to define CRISPR arrays. BMC Genom. 2016;17:1–14.
Article  CAS  Google Scholar 

36.
Biswas A, Gagnon JN, Brouns SJ, Fineran PC, Brown CM. CRISPRTarget: bioinformatic prediction and analysis of crRNA targets. RNA Biol. 2013;10:817–27.
CAS  PubMed  PubMed Central  Article  Google Scholar 

37.
Joshi N. & Fass J. Sickle: A sliding-window, adaptive, quality-based trimming tool for fastQ files (version 1.33) [software]. 2011. https://github.com/najoshi/sickle.

38.
Li H, Durbin R. Fast and accurate short read alignment with Burrows–Wheeler transform. Bioinformatics. 2009;25:1754–60.
CAS  PubMed  PubMed Central  Article  Google Scholar 

39.
Anders S, Pyl PT, Huber W. HTSeq—a Python framework to work with high-throughput sequencing data. Bioinformatics. 2015;31:166–9.
CAS  PubMed  PubMed Central  Article  Google Scholar 

40.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.
PubMed  PubMed Central  Article  CAS  Google Scholar 

41.
Galili T. dendextend: an R package for visualizing, adjusting and comparing trees of hierarchical clustering. Bioinformatics. 2015;31:3718–20.
CAS  PubMed  PubMed Central  Article  Google Scholar 

42.
Richter C, Dy RL, McKenzie RE, Watson BN, Taylor C, Chang JT, et al. Priming in the Type IF CRISPR-Cas system triggers strand-independent spacer acquisition, bi-directionally from the primed protospacer. Nucleic Acids Res. 2014;42:8516–26.
CAS  PubMed  PubMed Central  Article  Google Scholar 

43.
Fineran PC, Gerritzen MJ, Suárez-Diez M, Künne T, Boekhorst J, van Hijum SA, et al. Degenerate target sites mediate rapid primed CRISPR adaptation. Proc Natl Acad Sci USA. 2014;111:E1629–38.
CAS  PubMed  Article  Google Scholar 

44.
Xue C, Seetharam AS, Musharova O, Severinov K, Brouns J, Severin SJ, et al. CRISPR interference and priming varies with individual spacer sequences. Nucleic Acids Res. 2015;43:10831–47.
CAS  PubMed  PubMed Central  Article  Google Scholar 

45.
Wortham BW, Oliveira MA, Patel CN. Polyamines in bacteria: pleiotropic effects yet specific mechanisms. In: (eds Perry RD, Fetherston JD). The genus Yersinia. New York, NY: Springer; 2007. p. 106–15.

46.
Bru JL, Rawson B, Trinh C, Whiteson K, Høyland-Kroghsbo NM, Siryaporn A. PQS produced by the Pseudomonas aeruginosa stress response repels swarms away from bacteriophage and antibiotics. J Bacteriol. 2019;201:e00383–19.

47.
Doron S, Fedida A, Hernández-Prieto MA, Sabehi G, Karunker I, Stazic D, et al. Transcriptome dynamics of a broad host-range cyanophage and its hosts. ISME J. 2016;10:1437–55.
CAS  PubMed  Article  Google Scholar 

48.
Landsberger M, Gandon S, Meaden S, Rollie C, Chevallereau A, Chabas H, et al. Anti-CRISPR phages cooperate to overcome CRISPR-Cas immunity. Cell. 2018;174:908–16.
CAS  PubMed  PubMed Central  Article  Google Scholar 

49.
Wei Y, Terns RM, Terns MP. Cas9 function and host genome sampling in Type II-A CRISPR-Cas adaptation. Genes Dev. 2015;29:356–61.
CAS  PubMed  PubMed Central  Article  Google Scholar 

50.
Fokine A, Rossmann MG. Common evolutionary origin of procapsid proteases, phage tail tubes, and tubes of bacterial type VI secretion systems. Structure. 2016;24:1928–35.
CAS  PubMed  PubMed Central  Article  Google Scholar 

51.
Modell JW, Jiang W, Marraffini LA. CRISPR-Cas systems exploit viral DNA injection to establish and maintain adaptive immunity. Nature. 2017;544:101–4.
CAS  PubMed  PubMed Central  Article  Google Scholar 

52.
van Sluijs L, van Houte S, van Der Oost J, Brouns SJ, Buckling A, Westra ER. Addiction systems antagonize bacterial adaptive immunity. FEMS Microbiol Lett. 2019;366:fnz047.
CAS  PubMed  PubMed Central  Google Scholar 

53.
Wagemans J, Blasdel BG, Van den Bossche A, Uytterhoeven B, De Smet J, Paeshuyse J, et al. Functional elucidation of antibacterial phage ORFans targeting Pseudomonas aeruginosa. Cell Microbiol. 2014;16:1822–35.
CAS  PubMed  Article  Google Scholar 

54.
van Houte S, Ros VI, van Oers MM. Walking with insects: molecular mechanisms behind parasitic manipulation of host behaviour. Mol Ecol. 2013;22:3458–75.
PubMed  Article  Google Scholar 

55.
Taylor TB, Buckling A. Bacterial motility confers fitness advantage in the presence of phages. J Evolut Biol. 2013;26:2154–60.
CAS  Article  Google Scholar 

56.
Heussler GE, Cady KC, Koeppen K, Bhuju S, Stanton BA, O’Toole GA. Clustered regularly interspaced short palindromic repeat-dependent, biofilm-specific death of Pseudomonas aeruginosa mediated by increased expression of phage-related genes. MBio. 2015;6:e00129–15.

57.
Hernandez CA, Koskella B. Phage resistance evolution in vitro is not reflective of in vivo outcome in a plant‐bacteria‐phage system. Evolution. 2019;73:2461–75.
CAS  PubMed  Article  Google Scholar 

58.
Alseth EO, Pursey E, Luján AM, McLeod I, Rollie C, Westra ER. Bacterial biodiversity drives the evolution of CRISPR-based phage resistance. Nature. 2019;574:549–52.
CAS  PubMed  PubMed Central  Article  Google Scholar 

59.
Westra E, Levin BR. How important is CRISPR-Cas for protecting natural populations of bacteria against infections with badass DNAs? BioRxiv. 2020. https://www.biorxiv.org/content/10.1101/2020.02.05.935965v2. More

1.
Bergquist PR. Sponges. London, United Kingdom: Hutchinson and Co. Ltd; 1978.
2.
Brain CKB, Prave AR, Hoffmann KH, Fallick AE, Botha A, Herd DA, et al. The first animals: Ca. 760-million-year-old sponge-like fossils from Namibia. S Afr J Sci. 2012;108:1–8.
Article  Google Scholar 

3.
Hooper JNA, van Soest RWM. Systema Porifera: a guide to the classification of sponges. New York: Kluwer Academic/Plenum Publishers; 2002.
Google Scholar 

4.
Webster NS. Sponge disease: a global threat? Environ Microbiol. 2007;9:1363–75.
CAS  PubMed  Article  PubMed Central  Google Scholar 

5.
Bell JJ. The functional roles of marine sponges. Estuar Coast Shelf Sci. 2008;79:341–53.
Article  Google Scholar 

6.
de Goeij JM, van Oevelen D, Vermeij MJA, Middelburg JJ, Osinga R, de Goeij AFPM, et al. Surviving in a Marine Desert: the sponge loop retains resources within coral reefs. Science. 2013;342:108–10.
PubMed  Article  CAS  PubMed Central  Google Scholar 

7.
Taylor MW, Radax R, Steger D, Wagner M. Sponge-associated microorganisms: evolution, ecology, and biotechnological potential. Microbiol Mol Biol Rev. 2007;71:295–347.
CAS  PubMed  PubMed Central  Article  Google Scholar 

8.
Wehrl M, Steinert M, Hentschel U. Bacterial uptake by the marine sponge Aplysina aerophoba. Microb Ecol. 2007;53:355–65.
PubMed  Article  PubMed Central  Google Scholar 

9.
Thomas T, Rusch D, DeMaere MZ, Yung PY, Lewis M, Halpern A, et al. Functional genomic signatures of sponge bacteria reveal unique and shared features of symbiosis. ISME J. 2010;4:1557–67.
CAS  PubMed  Article  PubMed Central  Google Scholar 

10.
Nguyen MTHD, Liu M, Thomas T. Ankyrin-repeat proteins from sponge symbionts modulate amoebal phagocytosis. Mol Ecol. 2013;23:1635–45.
PubMed  Article  CAS  Google Scholar 

11.
Reynolds D, Thomas T. Evolution and function of eukaryotic-like proteins from sponge symbionts. Mol Ecol. 2016;25:5242–53.
CAS  PubMed  Article  PubMed Central  Google Scholar 

12.
Jahn MT, Arkhipova K, Markert SM, Stigloher C, Lachnit T, Pita L, et al. A phage protein aids bacterial symbionts in eukaryote immune evasion. Cell Host Microbe. 2019;26:542–50.e5.
CAS  PubMed  PubMed Central  Article  Google Scholar 

13.
Degnan SM. The surprisingly complex immune gene repertoire of a simple sponge, exemplified by the NLR genes: a capacity for specificity? Dev Comp Immunol. 2015;48:269–74.
CAS  PubMed  Article  PubMed Central  Google Scholar 

14.
Hentschel U, Usher KM, Taylor MW. Marine sponges as microbial fermenters. FEMS Microbiol Ecol. 2006;55:167–77.
CAS  PubMed  Article  PubMed Central  Google Scholar 

15.
Moitinho-Silva L, Steinert G, Nielsen S, Hardoim CCP, Wu YC, McCormack GP, et al. Predicting the HMA-LMA status in marine sponges by machine learning. Front Microbiol. 2017;8:1–14.
Article  Google Scholar 

16.
Moitinho-Silva L, Nielsen S, Amir A, Gonzalez A, Ackermann GL, Cerrano C, et al. The sponge microbiome project. Gigascience. 2017;6:1–7.
CAS  PubMed  Article  PubMed Central  Google Scholar 

17.
Hentschel U, Hopke J, Horn M, Friedrich AB, Wagner M, Hacker J, et al. Molecular evidence for a uniform microbial community in sponges from different oceans. Appl Environ Microbiol. 2002;68:4431–40.
CAS  PubMed  PubMed Central  Article  Google Scholar 

18.
Simister RL, Deines P, Botté ES, Webster NS, Taylor MW. Sponge-specific clusters revisited: a comprehensive phylogeny of sponge-associated microorganisms. Environ Microbiol. 2012;14:517–24.
CAS  PubMed  Article  PubMed Central  Google Scholar 

19.
Taylor MW, Tsai P, Simister RL, Deines P, Botte E, Ericson G, et al. “Sponge-specific” bacteria are widespread (but rare) in diverse marine environments. ISME J. 2013;7:438–43.
CAS  PubMed  Article  PubMed Central  Google Scholar 

20.
Croué J, West NJ, Escande M-L, Intertaglia L, Lebaron P, Suzuki MT. A single betaproteobacterium dominates the microbial community of the crambescidine-containing sponge Crambe crambe. Sci Rep. 2013;3:2583.
PubMed  PubMed Central  Article  Google Scholar 

21.
Fan L, Reynolds D, Liu M, Stark M, Kjelleberg S, Webster NS, et al. Functional equivalence and evolutionary convergence in complex communities of microbial sponge symbionts. PNAS. 2012;109:E1878–87.
CAS  PubMed  Article  Google Scholar 

22.
Simister RL, Taylor MW, Rogers KM, Schupp PJ, Deines P. Temporal molecular and isotopic analysis of active bacterial communities in two New Zealand sponges. FEMS Microbiol Ecol. 2013;85:195–205.
PubMed  Article  Google Scholar 

23.
Gantt SE, López-Legentil S, Erwin PM. Stable microbial communities in the sponge Crambe crambe from inside and outside a polluted Mediterranean harbor. FEMS Microbiol Lett. 2017;364:1–7.
Article  CAS  Google Scholar 

24.
Bin JeongJ, Kim KH, Park JS. Sponge-specific unknown bacterial groups detected in marine sponges collected from Korea through barcoded pyrosequencing. J Microbiol Biotechnol. 2015;25:1–10.
Article  Google Scholar 

25.
Thiel V, Neulinger SC, Staufenberger T, Schmaljohann R, Imhoff JF. Spatial distribution of sponge-associated bacteria in the Mediterranean sponge Tethya aurantium. FEMS Microbiol Ecol. 2006;59:47–63.
PubMed  Article  CAS  Google Scholar 

26.
Fieth RA, Gauthier M-EA, Bayes J, Green KM, Degnan SM. Ontogenetic changes in the bacterial symbiont community of the tropical demosponge Amphimedon queenslandica: metamorphosis is a new beginning. Front Mar Sci. 2016;3:1–20.
Article  Google Scholar 

27.
Batani G. Fluorescence in situ hybridisation for the localisation and culturing of marine bacteria: co-localisation of symbionts in sponges (Unpublished PhD thesis chapter). UNSW Sydney. 2018.

28.
Waterworth SC, Jiwaji M, Kalinski JCJ, Parker-Nance S, Dorrington RA. A place to call home: an analysis of the bacterial communities in two Tethya rubra Samaai and Gibbons 2005 populations in algoa bay, South Africa. Mar Drugs. 2017;15:95.
PubMed Central  Article  CAS  Google Scholar 

29.
Wu S, Ou H, Liu T, Wang D, Zhao J. Structure and dynamics of microbiomes associated with the marine sponge Tedania sp. during its life cycle. FEMS Microbiol Ecol. 2018;94:1–9.
Article  CAS  Google Scholar 

30.
Matcher GF, Waterworth SC, Walmsley TA, Matsatsa T, Parker-Nance S, Davies-Coleman MT, et al. Keeping it in the family: coevolution of latrunculid sponges and their dominant bacterial symbionts. Microbiologyopen. 2017;6:1–13.
Article  CAS  Google Scholar 

31.
Jackson SA, Flemer B, McCann A, Kennedy J, Morrissey JP, O’Gara F, et al. Archaea appear to dominate the microbiome of Inflatella pellicula deep sea sponges. PLoS ONE. 2013;8:1–8.
Google Scholar 

32.
Steinert G, Taylor MW, Deines P, Simister RL, de Voogd NJ, Hoggard M, et al. In four shallow and mesophotic tropical reef sponges from Guam the microbial community largely depends on host identity. PeerJ. 2016;4:e1936.
PubMed  PubMed Central  Article  CAS  Google Scholar 

33.
Gauthier M-EA, Watson JR, Degnan SM. Draft genomes shed light on the dual bacterial symbiosis that dominates the microbiome of the coral reef sponge Amphimedon queenslandica. Front Mar Sci. 2016;3:1–18.
Article  Google Scholar 

34.
Parks DH, Imelfort M, Skennerton CT, Hugenholtz P, Tyson GW. CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Cold Spring Harb Lab Press. 2014;25:1043–55.
Google Scholar 

35.
Mukherjee S, Stamatis D, Bertsch J, Ovchinnikova G, Verezemska O, Isbandi M, et al. Genomes OnLine database (GOLD) v.6: data updates and feature enhancements. Nucleic Acids Res. 2017;45:D446–56.
CAS  PubMed  Article  Google Scholar 

36.
Huntemann M, Ivanova NN, Mavromatis K, James Tripp H, Paez-Espino D, Palaniappan K, et al. The standard operating procedure of the DOE-JGI microbial genome annotation pipeline (MGAP v.4). Stand Genom Sci. 2015;10:4–9.
Article  CAS  Google Scholar 

37.
Eddy S. Profile hidden Markov models. Bioinformatics. 1998;14:755–63.
CAS  PubMed  Article  Google Scholar 

38.
Pruesse E, Peplies J, Glöckner FO. SINA: accurate high-throughput multiple sequence alignment of ribosomal RNA genes. Bioinformatics. 2012;28:1823–9.
CAS  PubMed  PubMed Central  Article  Google Scholar 

39.
Ludwig W, Strunk O, Westram R, Richter L, Meier H, Yadhukumar, et al. ARB: a software environment for sequence data. Nucleic Acids Res. 2004;32:1363–71.
CAS  PubMed  PubMed Central  Article  Google Scholar 

40.
Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, et al. The SILVA ribosomal RNA gene database project: Improved data processing and web-based tools. Nucleic Acids Res. 2013;41:590–6.
Article  CAS  Google Scholar 

41.
Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30:1312–3.
CAS  PubMed  PubMed Central  Article  Google Scholar 

42.
Letunic I, Bork P. Interactive tree of life (iTOL) v3: an online tool for the display and annotation of phylogenetic and other trees. Nucleic Acids Res. 2016;44:W242–5.
CAS  PubMed  PubMed Central  Article  Google Scholar 

43.
Chaumeil P-A, Mussig AJ, Hugenholtz P, Parks DH. GTDB-Tk: a toolkit to classify genomes with the Genome Taxonomy Database. Bioinformatics. 2020;36:1925–7.
Google Scholar 

44.
Parks DH, Chuvochina M, Waite DW, Rinke C, Skarshewski A, Chaumeil PA, et al. A standardized bacterial taxonomy based on genome phylogeny substantially revises the tree of life. Nat Biotechnol. 2018;36:996–1004.
CAS  PubMed  Article  PubMed Central  Google Scholar 

45.
Konstantinidis KT, Tiedje JM. Towards a genome-based taxonomy for prokaryotes. J Bacteriol. 2005;187:6258–64.
CAS  PubMed  PubMed Central  Article  Google Scholar 

46.
Rodriguez-R LM, Konstantinidis KT. The enveomics collection: a toolbox for specialized analyses of microbial genomes and metagenomes. Peer J Prepr. 2016;4:e1900v1.
Google Scholar 

47.
Thomas T, Moitinho-Silva L, Lurgi M, Björk JR, Easson C, Astudillo-García C, et al. Diversity, structure and convergent evolution of the global sponge microbiome. Nat Commun. 2016;7:11870.
CAS  PubMed  PubMed Central  Article  Google Scholar 

48.
Díez-Vives C, Moitinho-Silva L, Nielsen S, Reynolds D, Thomas T. Expression of eukaryotic-like protein in the microbiome of sponges. Mol Ecol. 2017;26:1432–51.
Article  CAS  Google Scholar 

49.
Öztürk B, De Jaeger L, Smidt H, Sipkema D. Culture-dependent and independent approaches for identifying novel halogenases encoded by Crambe crambe (marine sponge) microbiota. Sci Rep. 2013;3:1–9.
Article  Google Scholar 

50.
Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics. 2011;12:1–16.
Article  Google Scholar 

51.
Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10:R25.
PubMed  PubMed Central  Article  CAS  Google Scholar 

52.
Li B, Ruotti V, Stewart RM, Thomson JA, Dewey CN. RNA-Seq gene expression estimation with read mapping uncertainty. Bioinformatics. 2010;26:493–500.
PubMed  Article  CAS  PubMed Central  Google Scholar 

53.
Wagner GP, Kin K, Lynch VJ. Measurement of mRNA abundance using RNA-seq data: RPKM measure is inconsistent among samples. Theory Biosci. 2012;131:281–5.
CAS  PubMed  Article  PubMed Central  Google Scholar 

54.
Webster NS, Hill RT. The culturable microbial community of the great barrier reef sponge Rhopaloeides odorabile is dominated by an α-Proteobacterium. Mar Biol. 2001;138:843–51.
CAS  Article  Google Scholar 

55.
Engelberts JP, Robbins SJ, de Goeij JM, Aranda M, Bell SC, Webster NS. Characterization of a sponge microbiome using an integrative genome-centric approach. ISME J. 2020;14:1100–10.
CAS  PubMed  Article  PubMed Central  Google Scholar 

56.
Moitinho-Silva L, Díez-Vives C, Batani G, Esteves AIS, Jahn MT, Thomas T. Integrated metabolism in sponge-microbe symbiosis revealed by genome-centered metatranscriptomics. ISME J. 2017;11:1651–66.
CAS  PubMed  PubMed Central  Article  Google Scholar 

57.
Lobb B, Tremblay BJM, Moreno-Hagelsieb G, Doxey AC. An assessment of genome annotation coverage across the bacterial tree of life. Microb Genom. 2020;6:e000341.
PubMed Central  Google Scholar 

58.
Lopez JV, Ranzer LK, Ledger A, Schoch B, Duckworth A, Mccarthy PJ, et al. Comparison of bacterial diversity within the coral reef sponge, axinella corrugata, the encrusting coral erythropodium caribaeorum. Proc. 11 Int. Coral Reef Symp. 2008;2:1362–6.
Google Scholar 

59.
Yarza P, Yilmaz P, Pruesse E, Glöckner FO, Ludwig W, Schleifer KH, et al. Uniting the classification of cultured and uncultured bacteria and archaea using 16S rRNA gene sequences. Nat Rev Microbiol. 2014;12:635–45.
CAS  PubMed  Article  PubMed Central  Google Scholar 

60.
Konstantinidis KT, Rosselló-Móra R, Amann R. Uncultivated microbes in need of their own taxonomy. ISME J. 2017;11:2399–406.
PubMed  PubMed Central  Article  Google Scholar 

61.
Bowers RM, Kyrpides NC, Stepanauskas R, Harmon-smith M, Doud D, Jarett J, et al. Minimum information about a single amplified genome (MISAG) and a metagenome-assembled genome (MIMAG) of bacteria and archaea. Nat Biotechnol. 2017;35:725–31.
CAS  PubMed  PubMed Central  Article  Google Scholar 

62.
Chuvochina M, Rinke C, Parks DH, Rappé MS, Tyson GW, Yilmaz P, et al. The importance of designating type material for uncultured taxa. Syst Appl Microbiol. 2019;42:15–21.
PubMed  Article  PubMed Central  Google Scholar 

63.
Tripp E. Crowell’s handbook of classical mythology. New York: Thomas Y. Crowell Company; 1970.

64.
Athanassakis AN. Hesiod: Theogony, Works and Days, Shield. 2nd ed. Baltimore and London: The John Hopkins University Press; 2005.

65.
Webster NS, Negri AP, Munro MMHG, Battershill CN. Diverse microbial communities inhabit Antarctic sponges. Environ Microbiol. 2004;6:288–300.
PubMed  Article  PubMed Central  Google Scholar 

66.
Coelho FJRC, Cleary DFR, Gomes NCM, Pólonia ARM, Huang YM, Liu LL, et al. Sponge prokaryote communities in Taiwanese coral reef and shallow hydrothermal vent ecosystems. Microb Ecol. 2018;75:239–54.
CAS  PubMed  Article  PubMed Central  Google Scholar 

67.
Sipkema D, de Caralt S, Morillo JA, Al-Soud WAB, Sørensen SJ, Smidt H, et al. Similar sponge-associated bacteria can be acquired via both vertical and horizontal transmission. Environ Microbiol. 2015;17:3807–21.
CAS  PubMed  Article  PubMed Central  Google Scholar 

68.
Bjork JR, Diez-Vives C, Astudillo-Garcia C, Archie E, Montoya JM. Vertical transmission of sponge microbiota is inconsistent and unfaithful. Nat Ecol Evol. 2019;3:1172–83.
PubMed  PubMed Central  Article  Google Scholar 

69.
Gonzalez-Zapata FL, Bongaerts P, Ramírez-Portilla C, Adu-Oppong B, Walljasper G, Reyes A, et al. Holobiont diversity in a reef-building coral over its entire depth range in the mesophotic zone. Front Mar Sci. 2018;5:1–13.
Article  Google Scholar 

70.
Yang S, Sun W, Zhang F, Li Z. Phylogenetically diverse denitrifying and ammonia-oxidizing bacteria in corals Alcyonium gracillimum and Tubastraea coccinea. Mar Biotechnol. 2013;15:540–51.
CAS  PubMed  Article  Google Scholar 

71.
Høj L, Levy N, Baillie BK, Clode PL, Strohmaier RC, Siboni N, et al. Crown-of-thorns sea star Acanthaster cf. solaris has tissue—characteristic microbiomes with potential roles in health and reproduction. Appl Environ Microbiol. 2018;84:1–18.
Article  Google Scholar 

72.
Schmitt S, Deines P, Behnam F, Wagner M, Taylor MW. Chloroflexi bacteria are more diverse, abundant, and similar in high than in low microbial abundance sponges. FEMS Microbiol Ecol. 2011;78:497–510.
CAS  PubMed  Article  PubMed Central  Google Scholar 

73.
Bergquist PR, Kelly-Borges M. An evaluation of the genus Tethya (Porifera: Demospongiae: Hadromerjda) with descriptions of new species from the southwest Pacific. Beagle Rec Mus Art Galleries North Territ. 1991;8:37–72.
Google Scholar 

74.
Hoshino T, Yilmaz LS, Noguera DR, Daims H, Wagner M. Quantification of target molecules needed to detect microorganisms by fluorescence in situ hybridization (FISH) and catalyzed reporter deposition-FISH. Appl Environ Microbiol. 2008;74:5068–77.
CAS  PubMed  PubMed Central  Article  Google Scholar 

75.
Moran NA, Wernegreen JJ. Lifestyle evolution in symbiotic bacteria: insights from genomics. Trends Ecol Evol. 2000;15:321–6.
CAS  PubMed  Article  PubMed Central  Google Scholar 

76.
Vacelet J, Donadey C. Electron microscope study of the association between some sponges and bacteria. J Exp Mar Bio Ecol. 1977;30:301–14.
Article  Google Scholar 

77.
Maldonado M. Intergenerational transmission of symbiotic bacteria in oviparous and viviparous demosponges, with emphasis on intracytoplasmically-compartmented bacterial types. J Mar Biol Assoc UK. 2007;87:1701–13.
Article  Google Scholar 

78.
Ilan M, Abelson A. The life of a sponge in a sandy lagoon. Biol Bull. 1995;189:363–9.
CAS  PubMed  Article  PubMed Central  Google Scholar 

79.
Borisov VB, Gennis RB, Hemp J, Verkhovsky MI. The cytochrome bd respiratory oxygen reductases. Biochim Biophys Acta Bioenerg. 2011;1807:1398–413.
CAS  Article  Google Scholar 

80.
Lavy A, Keren R, Yahel G, Ilan M. Intermittent hypoxia and prolonged suboxia measured in situ in a marine sponge. Front Mar Sci. 2016;3:1–11.
Article  Google Scholar 

81.
Hoffmann F, Larsen O, Thiel V, Rapp HT, Pape T, Michaelis W, et al. An anaerobic world in sponges. Geomicrobiol J. 2005;22:1–10.
Article  Google Scholar 

82.
Moreno-Vivian C, Cabello C, Martinez-Luque M, Blasco R, Castillo F. Prokaryotic nitrate reduction: molecular properties and functional distinction among bacterial nitrate reductases. J Bacteriol. 1999;181:6573–84.
CAS  PubMed  PubMed Central  Article  Google Scholar 

83.
Toffanin A, Wu Q, Maskus M, Casella S, Abruña HD, Shapleigh JP. Characterization of the gene encoding nitrite reductase and the physiological consequences of its expression in the nondenitrifying Rhizobium “hedysari” strain HCNT1. Appl Environ Microbiol. 1996;62:4019–25.
CAS  PubMed  PubMed Central  Article  Google Scholar 

84.
Jones PM, George AM. The ABC transporter structure and mechanism: perspectives on recent research. Cell Mol Life Sci. 2004;61:682–99.
CAS  PubMed  Article  PubMed Central  Google Scholar 

85.
Mulligan C, Fischer M, Thomas GH. Tripartite ATP-independent periplasmic (TRAP) transporters in bacteria and archaea. FEMS Microbiol Rev. 2011;35:68–86.
CAS  PubMed  Article  PubMed Central  Google Scholar 

86.
Vastermark A, Wollwage S, Houle ME, Rio R, Saier MH. Expansion of the APC superfamily of secondary carriers. Proteins Struct Funct Bioinform. 2014;82:2797–811.
CAS  Article  Google Scholar 

87.
Pao SS, Paulsen IT, Saier MH. Major facilitator superfamily. Microbiol Mol Biol Rev. 1998;62:1–34.
CAS  PubMed  PubMed Central  Article  Google Scholar 

88.
Ziegler C, Bremer E, Krämer R. The BCCT family of carriers: from physiology to crystal structure. Mol Microbiol. 2010;78:13–34.
CAS  PubMed  PubMed Central  Google Scholar 

89.
Andrade SLA, Einsle O. The Amt/Mep/Rh family of ammonium transport proteins (Review). Mol Membr Biol. 2007;24:357–65.
CAS  PubMed  Article  PubMed Central  Google Scholar 

90.
Kamke J, Sczyrba A, Ivanova N, Schwientek P, Rinke C, Mavromatis K, et al. Single-cell genomics reveals complex carbohydrate degradation patterns in poribacterial symbionts of marine sponges. ISME J. 2013;7:2287–300.
CAS  PubMed  PubMed Central  Article  Google Scholar 

91.
Clifford EL, Hansell DA, Varela MM, Nieto-Cid M, Herndl GJ, Sintes E. Crustacean zooplankton release copious amounts of dissolved organic matter as taurine in the ocean. Limnol Oceanogr. 2017;62:2745–58.
CAS  PubMed  PubMed Central  Article  Google Scholar 

92.
Karimi E, Keller-Costa T, Slaby BM, Cox CJ, da Rocha UN, Hentschel U, et al. Genomic blueprints of sponge-prokaryote symbiosis are shared by low abundant and cultivatable Alphaproteobacteria. Sci Rep. 2019;9:1–15.
Article  CAS  Google Scholar 

93.
Dahl C, Franz B, Hensen D, Kesselheim A, Zigann R. Sulfite oxidation in the purple sulfur bacterium Allochromatium vinosum: Identification of SoeABC as a major player and relevance of SoxYZ in the process. Microbiology. 2013;159:2626–38.
CAS  PubMed  Article  PubMed Central  Google Scholar 

94.
Bardischewsky F, Quentmeier A, Friedrich CG. The flavoprotein SoxF functions in chemotrophic thiosulfate oxidation of Paracoccus pantotrophus in vivo and in vitro. FEMS Microbiol Lett. 2006;258:121–6.
CAS  PubMed  Article  PubMed Central  Google Scholar 

95.
Gregersen LH, Bryant DA, Frigaard NU. Mechanisms and evolution of oxidative sulfur metabolism in green sulfur bacteria. Front Microbiol. 2011;2:1–14.
Article  CAS  Google Scholar 

96.
Lavy A, Keren R, Yu K, Thomas BC, Alvarez-Cohen L, Banfield JF, et al. A novel Chromatiales bacterium is a potential sulfide oxidizer in multiple orders of marine sponges. Environ Microbiol. 2018;20:800–14.
CAS  PubMed  Article  PubMed Central  Google Scholar 

97.
Tian R-M, Wang Y, Bougouffa S, Gao ZM, Cai L, Bajic V, et al. Genomic analysis reveals versatile heterotrophic capacity of a potentially symbiotic sulfur-oxidizing bacterium in sponge. Environ Microbiol. 2014;16:3548–61.
CAS  PubMed  Article  PubMed Central  Google Scholar 

98.
Keller MD, Kiene RP, Matrai PA, Bellows WK. Production of glycine betaine and dimethylsulfoniopropionate in marine phytoplankton. II. N-limited chemostat cultures. Mar Biol. 1999;135:249–57.
CAS  Article  Google Scholar 

99.
Diaz MR, Visscher PT, Taylor BF. Metabolism of dimethylsulfoniopropionate and glycine betaine by a marine bacterium. FEMS Microbiol Lett. 1992;96:61–5.
CAS  Article  Google Scholar 

100.
Sun J, Steindler L, Thrash JC, Halsey KH, Smith DP, Carter AE, et al. One carbon metabolism in SAR11 pelagic marine bacteria. PLoS ONE. 2011;6:e23973.
CAS  PubMed  PubMed Central  Article  Google Scholar 

101.
Equar MY, Tani Y, Mihara H. Purification and properties of glycine oxidase from pseudomonas putida KT2440. J Nutr Sci Vitaminol. 2015;61:506–10.
CAS  PubMed  Article  Google Scholar 

102.
Karimi E, Slaby BM, Soares AR, Blom J, Hentschel U, Costa R. Metagenomic binning reveals versatile nutrient cycling and distinct adaptive features in alphaproteobacterial symbionts of marine sponges. FEMS Microbiol Ecol. 2018;94:1–18.
Article  CAS  Google Scholar 

103.
Moreno-Pino M, Cristi A, Gillooly JF, Trefault N. Characterizing the microbiomes of Antarctic sponges: a functional metagenomic approach. Sci Rep. 2020;10:1–12.
Article  CAS  Google Scholar 

104.
Mayer J, Huhn T, Habeck M, Denger K, Hollemeyer K, Cook AM. 2,3-Dihydroxypropane-1-sulfonate degraded by Cupriavidus pinatubonensis JMP134: purification of dihydroxypropanesulfonate 3-dehydrogenase. Microbiology. 2010;156:1556–64.
CAS  PubMed  Article  Google Scholar 

105.
Michael AJ. Polyamines in eukaryotes, bacteria, and archaea. J Biol Chem. 2016;291:14896–903.
CAS  PubMed  PubMed Central  Article  Google Scholar 

106.
Tsukamoto S, Kato H, Hirota H, Fusetani N. Pseudoceratidine: a new antifouling spermidine derivative from the marine sponge Pseudoceratina purpurea. Tetrahedron Lett. 1996;37:1439–40.
CAS  Article  Google Scholar 

107.
Igarashi K, Kashiwagi K. Polyamine transport in bacteria and yeast. Biochem J. 1999;344:633–42.
CAS  PubMed  PubMed Central  Article  Google Scholar 

108.
Shah P, Swiatlo E. A multifaceted role for polyamines in bacterial pathogens. Mol Microbiol. 2008;68:4–16.
CAS  PubMed  Article  Google Scholar 

109.
Schneider BL, Reitzer L. Pathway and enzyme redundancy in putrescine catabolism in Escherichia coli. J Bacteriol. 2012;194:4080–8.
CAS  PubMed  PubMed Central  Article  Google Scholar 

110.
Mou X, Sun S, Rayapati P, Moran MA. Genes for transport and metabolism of spermidine in Ruegeria pomeroyi DSS-3 and other marine bacteria. Aquat Microb Ecol. 2010;58:311–21.
Article  Google Scholar 

111.
Dasu VV, Nakada Y, Ohnishi-Kameyama M, Kimura K, Itoh Y. Characterization and a role of Pseudomonas aeruginosa spermidine dehydrogenase in polyamine catabolism. Microbiology. 2006;152:2265–72.
CAS  PubMed  Article  PubMed Central  Google Scholar 

112.
Tofalo R, Cocchi S, Suzzi G. Polyamines and gut microbiota. Front Nutr. 2019;6:1–5.
Article  CAS  Google Scholar 

113.
Sharfstein ST, Keasling JD. Polyphosphate metabolism in Escherichia coli. Ann N Y Acad Sci. 1994;745:77–91.
CAS  PubMed  Article  PubMed Central  Google Scholar 

114.
Zhang F, Blasiak LC, Karolin JO, Powell RJ, Geddes CD, Hill RT. Phosphorus sequestration in the form of polyphosphate by microbial symbionts in marine sponges. Proc Natl Acad Sci USA. 2015;112:4381–6.
CAS  PubMed  Article  PubMed Central  Google Scholar 

115.
Wang L, Yan J, Wise MJ, Liu Q, Asenso J, Huang Y, et al. Distribution patterns of polyphosphate metabolism pathway and its relationships with bacterial durability and virulence. Front Microbiol. 2018;9:1–10.
Article  Google Scholar 

116.
Romano S, Schulz-Vogt HN, González JM, Bondarev V. Phosphate limitation induces drastic physiological changes, virulence-related gene expression, and secondary metabolite production in Pseudovibrio sp. strain FO-BEG1. Appl Environ Microbiol. 2015;81:3518–28.
CAS  PubMed  PubMed Central  Article  Google Scholar 

117.
Bergquist PR. In: eLS (ed) Porifera (Sponges). John Wiley & Sons, Ltd. 2001.

118.
Kovacs-Simon A, Titball RW, Michell SL. Lipoproteins of bacterial pathogens. Infect Immun. 2011;79:548–61.
CAS  PubMed  Article  PubMed Central  Google Scholar 

119.
Aoki S, Ito M, Iwasaki W. From β- To α-proteobacteria: The origin and evolution of rhizobial nodulation genes nodij. Mol Biol Evol. 2013;30:2494–508.
CAS  PubMed  Article  PubMed Central  Google Scholar 

120.
Darby AC, Cho NH, Fuxelius HH, Westberg J, Andersson SGE. Intracellular pathogens go extreme: genome evolution in the Rickettsiales. Trends Genet. 2007;23:511–20.
CAS  PubMed  Article  PubMed Central  Google Scholar 

121.
Kuo CH, Ochman H. Inferring clocks when lacking rocks: the variable rates of molecular evolution in bacteria. Biol Direct. 2009;4:35.
PubMed  PubMed Central  Article  CAS  Google Scholar  More