More stories

  • in

    A Tunguska sized airburst destroyed Tall el-Hammam a Middle Bronze Age city in the Jordan Valley near the Dead Sea

    Melted quartz grainsCrystalline quartz melts between 1670 °C (tridymite) and 1713 °C (cristobalite), and because quartz is pervasive and easily identified, melted grains serve as an important temperature indicator. At TeH, we observed that unmelted potsherds displayed no melted quartz grains, indicating exposure to low temperatures. On the other hand, most quartz grains on the surfaces of pottery, mudbricks, and roofing clay exhibited some degree of melting, and unmelted quartz grains were rare. Nearly all quartz grains found on broken, unmelted surfaces of potsherds were also unmelted. On melted pottery and mudbricks, melted quartz has an estimated density of 1 grain per 5 mm2.Melted quartz grains at TeH exhibit a wide range of morphologies. Some show evidence of partial melting that only melted grain edges and not the rest of the grain (Figs. 22, 23). Others displayed nearly complete melting with diffusion into the melted Ca–Al–Si matrix of pottery or mudbrick (Fig. 22). Melted quartz grains commonly exhibit vesiculation caused by outgassing (Figs. 22, 23), suggesting that those grains rose above quartz’s melting point of ~ 1713 °C.Figure 22SEM images of melted quartz grains on melted potsherd from the palace. (a) Highly melted quartz grain from the upper surface of melted pottery; shows flow lines of molten quartz in darker ‘neck’ at upper right; (b) manually constructed EDS-based phase map showing 100% quartz grain (green) embedded in Ca–Al–Si matrix of melted pottery (red); blue marks mixing zone between SiO2 and matrix at approximately  > 1713 °C, the melting point of quartz. Yellow arrow points to area depleted in oxygen, indicating high-temperature transformation to elemental Si mixed with melted SiO2. (c) Highly melted quartz grain; (d) manually constructed EDS-based phase map showing diffusion/mixing zone in blue with arrow pointing to bubble, indicating outgassing as grain reached temperatures above its melting point. (e) Quartz grain that has almost completely melted; (f) manually constructed EDS-based phase map showing the small remnant of a melted quartz grain (green) with a wide mixing zone (blue).Full size imageFigure 23SEM images of melted quartz grains on melted mudbrick from the palace. (a) Highly melted quartz grain; (b) manually constructed EDS-based phase map indicates center is pure SiO2 surrounded by melted mudbrick. Arrow points to vesicles indicating outgassing as grain temperature rose above ~ 1713 °C, the melting point of quartz. (c) The surface of a flattened quartz grain showing flow marks toward the upper right. High temperatures are required to lower the viscosity sufficiently for quartz to flow. (d) Manually constructed EDS-based phase map with an arrow pointing to vesicles indicating outgassing at high temperatures. (e) Close up of grain in panel ‘c’ showing flow marks (schlieren) at arrows. (f) Shattered, melted quartz splattered onto mudbrick meltglass; (g) manually constructed EDS-based phase map indicating that the blue area is SiO2; the yellow area is a shattered, thermally altered Fe-oxide grain.Full size imageAn SEM–EDS elemental map of one melted grain showed that the quartz had begun to dissociate into elemental Si (Fig. 22b). Another grain (Fig. 23c–e) displayed flow marks consistent with exposure to temperatures above 1713 °C where the viscosity of quartz falls low enough for it to flow easily. Another SEM–EDS analysis confirmed that one agglutinated mass of material is 100 wt.% SiO2 (Fig. 23f, g), suggesting that this polycrystalline quartz grain shattered, melted, and partially fused again.Discussion of melted quartzMoore et al.17 reported that during heating experiments, many quartz grains  50-µm-wide remained visually unaltered up to ~ 1700 °C. By 1850 °C, all quartz grains fully melted. These experiments establish a particle-size dependency and confirm confirmed the melting point for  > 50-µm-wide TeH quartz grains between ~ 1700–1850 °C. Melted  > 50-µm-wide quartz grains on the surfaces of melted pottery and mudbrick from the TeH destruction layer indicate exposure to these unusually high temperatures  > 1700 °C.Previously, Thy et al.70 proposed that glass at Abu Hureyra did not form during a cosmic impact, but rather, formed in biomass slag that resulted from thatched hut fires. However, Thy et al. did not determine whether or not high-temperature grains existed in the biomass slag. To test that claim, Moore et al.17 analyzed biomass slag from Africa and found only low-temperature melted grains with melting points of ~ 1200 °C, consistent with a temperature range for biomass slag of 1155–1290 °C, as reported by Thy et al.71. Upon testing the purported impact glass from Abu Hureyra, Moore et al.17 discovered high-temperature mineral grains that melt in the range of 1713° to  > 2000 °C, as are also found in TeH glass. These test results suggest that the melted glass from Abu Hureyra must have been exposed to higher temperatures than those associated with fires in thatched huts. Because of the presence of high-temperature minerals at TeH, we conclude that, as at Abu Hureyra, the meltglass could not have formed simply by burning thatched huts or wood-roofed, mudbrick buildings.Melted Fe- and Si-rich spherulesThe presence of melted spherulitic objects (“spherules”) has commonly been used to help identify and investigate high-temperature airburst/impact events in the sedimentary record. Although these objects are referred to here as “spherules,” they display a wide range of other impact-related morphologies that include rounded, sub-rounded, ovate, oblate, elongated, teardrop, dumbbell, and/or broken forms17,72,73,74,75,76,77,78,79,80,81,82. Optical microscopy and SEM–EDS are commonly used to identify and analyze spherules and the processes by which they are formed. Care is needed to conclusively distinguish high-temperature spherules produced by cosmic impacts from other superficially similar forms. Other such objects that frequently occur in sediments include anthropogenic spherules (typically from modern coal-fired power plants), authigenic framboids (Supporting Information, Fig. S7), rounded detrital magnetite, and volcanic spherules.Spherules in TeH sediment were investigated from stratigraphic sequences that include the MB II destruction layer at four locations: palace, temple, ring road, and wadi (Fig. 24). For the palace (Field UA, Square 7GG), the sequence spanned 28 cm with 5 contiguous samples of sediment ranging from 3-cm thick for the MB II destruction layer to 13-cm thick for some outlying samples. In the palace, 310 spherules/kg (Fig. 24d) were observed in the destruction layer with none found in samples above and below this layer. For the temple (Field LS, Square 42J), 5 continuous samples spanned 43 cm and ranged in thickness from 6 to 16 cm; the MB II layer contained ~ 2345 Fe- and Si-rich spherules/kg with 782/kg in the sample immediately below and none at other levels (Fig. 24c). Six contiguous samples from the ring road (Field LA, Square 28 M) spanned 30 cm with all 5 cm thick; the MB II destruction layer at this location contains 2150 spherules/kg with none detected in younger or older samples (Fig. 24b). Five discontinuous samples from the wadi spanned 170 cm, ranging from 10-cm thick for the destruction layer up to 20-cm thick for other samples; the MB II destruction layer at this location contained 2780 spherules/kg with none in samples from other levels (Fig. 24a, Supporting Information, Table S3). Notably, when melted mudbrick from the ring road was being mounted for SEM analysis, numerous loose spherules were observed within vesicles of the sample, confirming a close association between the spherules and meltglass. At all four locations, the peaks in high-temperature spherule abundances occur in the MB II destruction layer dating to ~ 1650 BCE.Figure 24Spherule abundances. (a)–(d). Number per kg for Fe- and S-rich spherules from 4 locations. Depths are in cm above or below the bottom of the destruction layer.Full size imageSEM images of spherules are shown in Figs. 25, 26, 27 and 28, and compositions are listed in Supporting Information, Table S4. The average spherule diameter was 40.5 µm with a range of 7 to 72 µm. The dominant minerals were Fe oxides averaging 40.2 wt.%, with a range of up to 84.1 wt.%; elemental Fe with a range of up to 80.3 wt.%; SiO2 averaging 20.9 wt.%, ranging from 1.0 to 45.2 wt.%; Al2O3 averaging 7.8 wt.% with a range of up to 15.6 wt.%; and TiO2 averaging 7.1 wt.% with a range of up to 53.1 wt.%. Fourteen spherules had compositions  > 48 wt.% of oxidized Fe, elemental Fe, and TiO2; five spherules contained  75 wt.% Fe with no Ti. Eight of 23 spherules analyzed contained detectable levels of Ti at up to 53.1 wt.%.Figure 25SEM images of mostly silica-rich spherules from TeH. (a)–(d) Representative spherules from the ring road on the lower tall. SEM images of iron-rich spherules. (e)–(f) Fe-rich spherules from the temple complex. (g) temple spherule containing ~ 3.7 wt.% Cr. (h) Broken, vesicular spherule from temple containing 1.4 wt.% Ni and 3.7 wt.% Cr. SEM images of titanium-rich spherules. Ti content of these ranges from 18.9 to 1.2 wt.%, averaging 10.7 wt.%. (i)–(k) Spherules from the ring road. (l) Spherule from the wadi site.Full size imageFigure 26SEM image of rare-earth (REE) spherule. (a) REE-rich 72-µm-wide spherule from the palace, dominantly composed of Fe, La, Ce, and O. (b) Close up of REE blebs found on the spherule. (c)–(f) SEM–EDS elemental maps showing composition. La = 15.6 wt.% and Ce = 21.0 wt.%. Ce is enriched over Fe and La in the middle part of the spherule, as seen in panels ‘d’ through ‘f’.Full size imageFigure 27SEM images of a spherule mainly composed of Fe and Si. (a) Fe–Ti-rich 54-µm-wide spherule from the palace. Spherule displays a protrusion to the left, suggesting aerodynamic shaping when molten, after which the tail detached. (b) A focused ion beam (FIB) was used to section the spherule, revealing inclusions of wassonite or titanium sulfide (TiS; yellow arrows) that are lighter-colored than the matrix. (c)–(f) Color-coded SEM–EDS elemental maps, showing the distribution of Ti, S, Si, and Fe and the location of the TiS grains. The spherule is dominantly composed of Fe and Si with minor amounts of Ti and S found in TiS inclusions.Full size imageFigure 28Fe-rich spherules embedded in meltglass. (a) Optical photomicrograph of a 167-µm-wide piece of meltglass with embedded Fe-rich spherules. (b) SEM image of same grain as in panel ‘a’. Melted quartz grain (Qtz) is embedded in Ca–Al–Si-rich matrix, which has the same composition as melted mudbrick. (c) SEM close-up image of the boxed area and panel ‘b’, showing splattered Fe-rich spherule.Full size imageTwo unusual spherules from the palace contain anomalously high percentages of rare-earth elements (REEs) at  > 37 wt.% of combined lanthanum (La), and cerium (Ce) (Fig. 26), as determined by preliminary measurements using SEM–EDS. Minor oxides account for the rest of the spherules’ bulk composition (Table S1).One 54-µm-wide sectioned spherule contains titanium sulfide (TiS) with a melting point of ~ 1780° C. TiS, known as wassonite, was first identified in meteorites (Fig. 27) and has been reported in impact-related material17,81,83. However, TiS sometimes occurs as an exsolution product forming fine networks in magnetite and ilmenite and can be of terrestrial origin.One unusual piece of 167-µm-wide Ca–Al–Si meltglass contains nearly two dozen iron oxide spherules on its surface (Fig. 28). The meltglass contains a completely melted quartz grain as part of the matrix (Fig. 28b). Most of the spherules appear to have been flattened or crushed by collision with the meltglass while they were still partially molten (Fig. 28c).Discussion of spherules and meltglassMelted materials from non-impact-related combustion have been reported in multiple studies. Consequently, we investigated whether Ca-, Fe-, and Si-rich spherules and meltglass (mudbrick, pottery, plaster, and roofing clay) may have formed normally, rather than from a cosmic impact event. For example, (i) glassy spherules and meltglass are known to form when carbon-rich biomass smolders below ground at ~ 1000° to 1300 °C, such as in midden mounds71. They also form in buried peat deposits84, underground coal seams85, burned haystacks86, and in large bonfires, such as at the Native American site at Cahokia, Illinois, in the USA87. (ii) Also, ancient fortifications (hillforts) in Scotland and Sweden, dating from ~ 1000 BCE to 1400 AD, have artificially vitrified walls that melted at temperatures of ~ 850° to 1000 °C88. (iii) Partially vitrified pottery and meltglass derived from the melting of wattle and daub (thatch and clay) with estimated temperatures of ~ 1000 °C have been reported in burned houses of the Trypillia culture in Ukraine89,90. (iv) Vitrified mudbricks and pottery that melted at 17 investigated biomass glass from midden mounds in Africa and found no high-temperature minerals. For this contribution, we used SEM–EDS to examine aluminosilicate meltglass from an underground peat fire in South Carolina, USA; meltglass in coal-fired fly ash from New Jersey, USA; and mining slag from a copper mine in Arizona, USA. All these meltglass examples display unmelted quartz and contain no other high-temperature melted grains, consistent with low-temperature melting at  97% wt.% FeO, as are found at TeH. Nor can these low temperatures produce meltglass and spherules embedded with melted zircon (melting point = 1687 °C), chromite (2190 °C), quartz (1713 °C), platinum (1768 °C), and iridium (2466 °C). Moore et al.17 confirmed that the melting of these high-temperature minerals requires minimum temperatures of ~ 1500° to 2500 °C.This evidence demonstrates that although the matrix of the spherules and meltglass at TeH likely experienced incipient melting at temperatures lower than ~ 1300 °C, this value represents only the minimum temperature of exposure, because the high-temperature minerals embedded in them do not melt at such low temperatures. Instead, the spherules and meltglass at TeH must have reached temperatures greater than ~ 1300 °C, most likely involving brief exposure to ambient temperatures of ~ 2500 °C, the melting point of iridium. These temperatures far exceed those characteristic of city fires and other types of biomass burning. In summary, all of this evidence is consistent with very high temperatures known during cosmic impacts but inconsistent with other known natural causes.Calcium carbonate spherules and plasterIn sediments of the destruction layer, we observed amber-to-off-white-colored spherules (Fig. 29) at high concentrations of ~ 240,000/kg in the palace, ~ 420/kg in the temple, ~ 60/kg on the ring road, and ~ 910/kg in the wadi (Supporting information, Table S2). In all four profiles, the spherules peak in the destruction layer with few to none above or below. Peak abundances of calcium carbonate spherules are closely associated with peak abundances of plaster fragments, which are the same color. By far the most spherules (~ 250× more) occurred in the destruction layer of the palace, where excavations showed that nearly every room and ceiling was surfaced with off-white lime-based plaster. Excavators uncovered high-quality lime plaster fragments still adhering to mudbricks inside the MB II palace complex, and in one palace room, we uncovered fragments of melted plaster (Fig. 29e). In contrast, lime plaster was very rarely used in buildings on the lower tall, including those near the temple.Figure 29Images of calcium carbonate spherules and melted plaster from TeH. (a) Photomicrographs of translucent, amber-colored CaCO3 spherules from the destruction layer in the palace. (b) SEM image of 83-µm carbonate spherule with impact or outgassing crater at arrow. (c) Photomicrograph of ~ 2-mm-wide piece of partially melted palace plaster from oxygen/propylene torch test, showing incipient melting at 1500 °C. Arrows point to hemispheric droplets emerging as spherules. (d) 142-µm cluster of 8 carbonate spherules with apparent impact or outgassing crater at arrow. (e) 64 × 30 mm piece of melted plaster that broke off the palace wall and became melted. It is composed only of calcium, carbon, and oxygen.Full size imageTo explore a potential connection between plaster and spherules, we performed SEM–EDS on samples of the palace plaster. Comparison of SEM–EDS analyses shows that the plaster composition has a  > 96% similarity to the spherule composition: CaCO3 = 71.4 wt.% in plaster versus 68.7 wt.% in the spherules; elemental C = 23.6 versus 26.3 wt.%; SiO2 = 2.4 versus 1.8 wt.%; MgO = 1.7 versus 2.0 wt.%; and SO3 = 0.94 versus 1.2 wt.%. The high carbon percentage and low sulfur content indicate that the plaster was made from calcium carbonate and not gypsum (CaSO4·2H2O). SEM imaging revealed that the plaster contains small plant parts, commonly used in plaster as a binder, and is likely the source of the high abundance of elemental C in the plaster. Inspection showed no evidence of microfossils, such as coccoliths, brachiopods, and foraminifera. The morphology of the spherules indicates that they are not authigenic or biological in origin.Discussion of carbonate plaster and spherulesOne of the earliest known uses of CaCO3-based plaster was in ~ 6750 BCE at Ayn Ghazal, ~ 35 km from TeH in modern-day Amman, Jordan97. At that site, multi-purpose lime plaster was used to make statues and figurines and to coat the interior walls of buildings. Because the production of lime-based plaster occurred at least 3000 years before TeH was destroyed, the inhabitants of TeH undoubtedly were familiar with the process. Typically, lime powder was produced in ancient times by stacking wood/combustibles interspersed with limestone rocks and then setting the stack on fire. Temperatures of ~ 800–1100 °C were required to transform the rocks into crumbly chalk, which was then mixed with water to make hydrated lime and plastered onto mudbrick walls97.At TeH, fragments of CaCO3-based plaster are intermixed in covarying abundances with CaCO3-based spherules with both compositions matching to within 96%. This similarity suggests that the carbonate spherules are derived from the plaster. We infer that the high-temperature blast wave from the impact event stripped some plaster from the interior walls of the palace and melted some into spherules. However, it is difficult to directly melt CaCO3, which gives off CO2 at high temperatures and decomposes into lime powder. We investigated this cycle in a heating experiment with an oxygen/propylene torch and found that we could decompose the plaster at ~ 1500 °C, the upper limit of the heating test, and begin incipient melting of the plaster. The heated plaster produced emergent droplets at that temperature but did not transform into free spherules (Supporting Information, Text S2).Similar spherules have been reported from Meteor Crater, where spherules up to ~ 200 μm in diameter are composed entirely of CaCO3 formed from a cosmic impact into limestone98,99. One of several possible hypotheses for TeH is that during the impact event, the limestone plaster converted to CaO with an equilibrium melting point of 2572 °C. However, it is highly likely that airborne contaminants, such as sodium and water vapor, reacted with the CaO and significantly lowered the melting point, allowing spherule formation at ≥ 1500 °C.The proposed chemical sequence of events of plaster formation and the later impact are as follows:

    1.

    Limestone was heated to ~ 800–1100 °C, decomposing to quicklime:

    $${text{CaCO}}_{{3}} to {text{ CaO }} + {text{ CO}}_{{2}}$$

    2.

    Quicklime was mixed with water to make a wet plaster:$${text{CaO }} + {text{ H}}_{{2}} {text{O }} to {text{ Ca}}left( {{text{OH}}} right)_{{2}}$$

    3.

    The plaster hardened and slowly absorbed CO2 to revert to CaCO3:$${text{Ca}}left( {{text{OH}}} right)_{{2}} + {text{ CO}}_{{2}} to {text{ H}}_{{2}} {text{O }} + {text{ CaCO}}_{{3}}$$

    4.

    The high-temperature impact event melted some plaster into spherules:$${text{CaCO}}_{{3}} to {text{ CaO }}left( {{text{spherules}}} right) , + {text{ CO}}_{{2}} left( { > {15}00^circ {text{C}}} right)$$

    5.

    CaO spherules slowly absorbed CO2 to revert to CaCO3:$${text{Ca }} + {text{ CO}}_{{2}} to {text{ CaCO}}_{{3}} left( {text{as spherules}} right)$$

    General discussion of all spherulesAccording to the previous investigations17,72,81,82, Fe-rich spherules such as those found at TeH typically melt at  > 1538 °C, the melting point of iron (Table 1). Because of the presence of magnetite (Fe3O4) in the REE spherule, its melting point is inferred to be  > 1590 °C (Table 1). The Si-rich spherules are similar in composition to TeH sediment and mudbrick, and thus, we propose that they were derived from the melting of these materials at  > 1250 °C. The carbonate-rich spherules likely formed at  > 1500 °C.Several studies describe a mechanism by which spherules could form during a low-altitude cosmic airburst100,101. When a bolide enters Earth’s atmosphere, it is subjected to immense aerodynamic drag and ablation, causing most of the object to fragment into a high-temperature fireball, after which its remaining mass is converted into a high-temperature vapor jet that continues at hypervelocity down to the Earth’s surface. Depending on the altitude of the bolide’s disruption, this jet is capable of excavating unconsolidated surficial sediments, melting them, and ejecting the molten material into the air as Si- and Fe-rich spherules and meltglass. This melted material typically contains a very low percentage (17.Melted zircons in pottery and mudbricks were observed (Fig. 30) at an estimated density of 1 grain per 20 mm2. On highly melted surfaces, nearly all zircons showed some degree of melting. In contrast, nearly all zircons found on broken interior surfaces were unmelted (Fig. 30d), except those within ~ 1 mm of melted surfaces. This implies that the temperature of the surrounding atmosphere was higher than the internal temperatures of the melting objects. Unmelted potsherds displayed only unmelted minerals.Figure 30SEM images of melted zircon grains. (a) Melted TeH zircon grain with bubbles at yellow arrow due to high-temperature dissociation and/or entrapped porosity. (b) Melted TeH zircon grain decorated with bubbles along the fracture line at upper arrow; arrows labeled “Bd” point to bright granular baddeleyite, ZrO2, formed during the high-temperature dissociation of zircon. (c) Almost fully melted TeH zircon grain mixing into the Ca–Al–Si matrix. (d) A typical unmelted zircon grain from TeH with straight, euhedral edges. Grain shows cracks on the top surface from possible thermal or mechanical damage. (e) For comparison, from cosmic airburst/impact at Dakhleh Oasis in Egypt: melted zircon decorated with lines of bubbles (arrow).Full size imageThe melted zircons in TeH materials exhibit a wide range of morphologies. Most showed evidence of sufficient melting to alter or destroy the original distinctive, euhedral shape of the grains. Also, the grains were often decorated with vesicles that were associated with fractures (Fig. 30a, c).Stoichiometric zircon contains 67.2 wt.% and 32.8 wt.% ZrO2 and SiO2 respectively, but in several TeH samples, we observed a reduction in the SiO2 concentration due to a loss of volatile SiO from the dissociation of SiO2. This alteration has been found to occur at 1676 °C, slightly below zircon’s melting point of 1687 °C103. This zircon dissociation leads to varying ZrO2:SiO2 ratios and to the formation of distinctive granular textures of pure ZrO2, also known as baddeleyite104 (Figs. 30, 31, 32). With increasing time at temperature, zircon will eventually convert partially or completely to ZrO2. Nearly all zircons observed on the surfaces of melted materials were either melted or showed some conversion to baddeleyite. We observed one zircon grain (Fig. 32d–e) displaying granular ZrO2 associated with three phases that span a wide range of SiO2 concentrations, likely formed at temperatures above 1687 °C. This extreme temperature and competing loss of SiO over an inferred duration of only several seconds led to complex microstructures, where grains melted, outgassed, and diffused into the surrounding matrix.Figure 31SEM images of other melted zircon grains in palace potsherd. (a) Two melted zircon grains adjacent to a previously discussed melted quartz grain; (b) close-up of same zircon grains; (c) manually constructed EDS-based phase map showing baddeleyite grains in green. The blue area represents melted zircon, while the red background represents the Ca–Al–Si matrix of the melted pottery. (d) Manually constructed EDS-based phase map of zircon grain showing small baddeleyite grains in green at the top.Full size imageFigure 32SEM images of melted zircon grains in mudbrick meltglass from the palace. (a) Thermally distorted zircon grain with a “hook” that resulted from the flow of molten material at  > 1687 °C; the darker area represents unrelated debris on top of zircon. (b) Manually constructed EDS-based phase map showing baddeleyite grains (Bd = ZrO2) in green, zircon in blue, and melted mudbrick in red. (c) Zircon grain showing limited thermal alteration, yet sufficient to cause dissociation into bright baddeleyite grains at ~ 1676 °C. (d) Zircon grain exhibiting three phases of thermal alteration, as shown in detail in (e), where a manually constructed EDS-based phase map demonstrates that high temperatures caused bubbling in the center band of zircon (purple = Hi) producing sub-micron-sized grains of baddeleyite (e.g., at arrow). Medium temperatures caused zircon to melt and flow (blue = Lo), and lower temperatures at the left end of grain produced thermal cracks (medium blue = Med). The green area marks the high-Si diffusion zone resulting from the dissociation of zircon. (f) Zircon grain from TeH has been fully converted to granular baddeleyite.Full size imageDiscussion of melted zirconZircon grains have a theoretical, equilibrium melting point of ~ 1687 °C. Under laboratory heating17, zircon grains showed no detectable alteration in shape at ~ 1300 °C but displayed incipient melting of grain edges and dissociation to baddeleyite beginning at ~ 1400 °C with increasing dissociation to 1500 °C17. Most zircon grains  120 µm were still recognizable but displayed considerable melting17. These experiments establish a lower melting range for TeH zircon grains of ~ 1400° to 1500 °C.Patterson105 showed that zircon dissociation becomes favorable above 1538 °C and particles between 1 and 100 µm in size melted and dissociated when passing through a plasma, forming spherules with various amounts of SiO2 glass containing ZrO2 crystallites ranging in size from 5 nm to 1 µm. The majority of zircon crystals were monoclinic, but tetragonal ZrO2 was observed for the smaller crystallite sizes. Residence times were in the order of 100 ms, and the specific ZrO2 to SiO2 ratio within each spherule depended on the particle’s time at temperature106.Bohor et al.104 presented images of impact-shocked zircons from the K-Pg impact event at 66 Ma that are morphologically indistinguishable from those at TeH. Decorated zircon grains are uncommon in nature but commonly associated with cosmic impact events, as evidenced by two partially melted zircons from the known airburst/impact at Dakhleh Oasis, Egypt (Fig. 30e). The presence of bubbles indicates that temperatures reached at least 1676 °C, where the zircon began to dissociate and outgas. Similar dissociated zircon grains also have been found in tektite glass and distal fallback ejecta (deposited from hot vapor clouds). Granular baddeleyite-zircon has been found in the ~ 150-km-wide K-Pg impact crater107 and the 28-km-wide Mistatin Lake crater in Canada107. The dissociation of zircon requires high temperatures of ~ 1676 °C104, implying that TeH was exposed to similar extreme conditions.Melted chromite grainsExamples of melted chromite, another mineral that melts at high temperatures, were also observed. Thermally-altered chromite grains were observed in melted pottery, melted mudbricks, and melted roofing clay from the palace. Their estimated density was 1 grain per 100 mm2, making them rarer than melted zircon grains. The morphologies of chromite grains range from thermally altered (Fig. 33a) to fully melted (Fig. 33b, d). One chromite grain from the palace displays unusual octahedral cleavage or shock-induced planar fractures (Fig. 33b). The typical chemical composition for chromite is 25.0 wt.% Fe, 28.6 wt.% O, and 46.5 wt.% Cr, although the Cr content can vary from low values to ~ 68 wt.%. SEM images reveal that, as chromite grains melted, some Cr-rich molten material migrated into and mixed with the host melt, causing an increase in Cr and Fe, and corresponding depletion of Si. The ratio of Cr to Fe in chromite affects its equilibrium melting point, which varies from ~ 1590 °C for a negligible amount of Cr up to ~ 2265 °C for ~ 46.5 wt.% Cr as in chromite or chromian magnetite ((Fe)Cr2O4), placing the melting point of TeH chromite at close to 2265 °C.Figure 33SEM images of melted chromite grains found on a melted potsherd from the palace. (a) Shattered, polycrystalline chromite grain that appears to have become agglutinated while molten. (b) Melted chromite grain, displaying cleavage (lamellae) suggestive of thermal and/or mechanical shock metamorphism at ~ 12 GPa; (c) close-up image showing angles between three sets of crystalline cleavage; (d) manually constructed EDS-based phase map showing chromite (purple) embedded in Ca–Al–Si matrix. The lines mark three sets of cleavage extending across the entire grain. A melt tail merging with the matrix is observed to trail off to the upper right of the grain at arrow.Full size imageDiscussion of melted chromiteChromite grains theoretically melt at ~ 2190 °C. Moore et al.17 reported the results of heating experiments in which chromite grains in bulk sediment showed almost no thermal alteration up to ~ 1500 °C (Supporting Information, Fig. S8). At temperatures of ~ 1600 °C and ~ 1700 °C, the shapes of chromite grains were intact but exhibited limited melting of grain edges. These results establish a range of ~ 1600° to 1700 °C for melting chromite grains.Because chromite typically does not exhibit cleavage, the grain exhibiting this feature is highly unusual. Its origin is unclear but there are several possibilities. The cleavage may have resulted from exsolution while cooling in the source magma. Alternately, the lamellae may have resulted from mechanical shock during a cosmic impact, under the same conditions that produced the shocked quartz, as reported by Chen et al.108 for meteorites shocked at pressures of ~ 12 GPa. Or they may have been formed by thermal shock, i.e., rapid thermal loading followed by rapid quenching. This latter suggestion is supported by the observation that the outside glass coating on the potsherd does not exhibit any quench crystals, implying that the cooling progressed very rapidly from liquid state to solid state (glass). This is rare in terrestrial events except for some varieties of obsidian, but common in melted material produced by atomic detonations (trinitite), lightning strikes (fulgurites), and cosmic airburst/impacts (meltglass)81. More investigations are needed to determine the origin of the potentially shocked chromite.Nuggets of Ir, Pt, Ru, Ni, Ag, Au, Cr, and Cu in meltglassUsing SEM–EDS, we investigated abundances and potential origins (terrestrial versus extraterrestrial) of platinum-group elements (PGEs) embedded in TeH meltglass, in addition to Ni, Au, and Ag. Samples studied include melted pottery (n = 3); melted mudbrick (n = 6); melted roofing clay (n = 1), and melted lime-based building plaster (n = 1). On the surfaces of all four types of meltglass, we observed melted metal-rich nuggets and irregularly shaped metallic splatter, some with high concentrations of PGEs (ruthenium (Ru), rhodium (Rh), palladium (Pd), osmium (Os), iridium (Ir), and platinum (Pt)) and some nuggets enriched in silver (Ag), gold (Au), chromium (Cr), copper (Cu), and nickel (Ni) with no PGEs (Figs. 34, 35). Importantly, these metal-rich nuggets were observed only on the top surfaces of meltglass and not inside vesicles or on broken interior surfaces.Figure 34SEM images of nuggets of melted metals in mudbrick meltglass from the palace. (a)–(c) Pt-dominant TeH nuggets enriched in ruthenium (Ru), rhodium (Rh), palladium (Pd), osmium (Os), iridium (Ir), and platinum (Pt). (d)–(f) Fe-dominant TeH splatter is also enriched in PGEs. (g)–(i) Nuggets enriched in varying percentages and combinations of nickel (Ni), chromium (Cr), copper (Cu), and silver (Ag).Full size imageFigure 35Average composition of selected metal-rich nuggets from the palace. (a-h) Silver (Ag), gold (Au), chromium (Cr), copper (Cu), iridium (Ir), nickel (Ni), platinum (Pt), and ruthenium (Ru), showing wt.% in selected nuggets from the destruction layer of the palace (7GG).Full size imageUsing SEM–EDS, we identified variable concentrations and assemblages of PGEs. The metallic particles appear to have melted at high temperatures based on the minimum melting points of the elements: iridium at 2466 °C; platinum = 1768 °C; and ruthenium = 2334 °C, indicating a temperature range of between approximately 1768° and 2466 °C. Our investigations also identified two PGE groups, one with nuggets in which Pt dominates Fe and the other with metallic splatter in which Fe dominates Pt.Pt-dominant nuggetsWe conducted 21 measurements on Pt-dominant TeH nuggets on meltglass (Fig. 34a–c). The nuggets average ~ 5 µm in length (range 1–12 µm) with an estimated concentration of 1 nugget per 10 mm2. For these nuggets, Fe concentrations average 1.0 wt.%, Ir = 6.0 wt.%, and Pt = 44.9 wt.% (Supporting Information, Tables S6, S7). The presence of PGEs was confirmed by two SEM–EDS instruments that verified the accurate identification of PGEs through analyses of several blanks that showed no PGE content. Some concentrations are low ( Pt or Pt  > Fe were found to be consistent between the two instruments.To determine the source of TeH nuggets and splatter, we constructed ternary diagrams. Terrestrial PGE nuggets are commonly found in ore bodies that when eroded, can become concentrated in riverine placer deposits, including those of the Jordan River floodplain. To compare Fe–Ir–Pt relationships among the TeH nuggets, we compiled data from nearby placer deposits in Greece109, Turkey110,111, and Iraq112, along with distant placers in Russia113,114,115, Canada116, and Alaska, USA117,118. The compilation of 109 Pt-dominant placer nuggets indicates that the average Fe concentration is 8.2 wt.%, Ir = 2.9 wt.%, and Pt = 80.3 wt.%. For the Ir-dominant placer nuggets (n = 104), Fe = 0.4 wt.%, Ir = 47.8 wt.%, and Pt = 5.3 wt.% (Supporting Information, Tables S6, S7). The ternary diagrams reveal that the values for Pt-dominant TeH nuggets overlap with Pt-dominant terrestrial placer nuggets but the Fe-dominant splatter is dissimilar (Fig. 36a).Figure 36Ternary diagrams for PGE-rich grains. Comparison of Fe–Ir–Pt ratios of PGE-rich nuggets fused into the surfaces of TeH meltglass. There are two populations of TeH nuggets (red diamonds): Pt-dominant at #1 (top) and Fe-dominant at #2 (bottom left). (a) TeH Pt-dominant nugget group #1 (red diamonds) overlaps Pt-dominant but not Ir-dominant nuggets (blue circles) from placers and ophiolite deposits in Greece, Turkey, Iraq, Russia, Canada, and the USA. The Fe-dominant TeH nugget group #2 is geochemically dissimilar to all known placer nuggets, suggesting that these nuggets are not placer-derived. (b) TeH nuggets (red diamonds) compared to nuggets in carbonaceous chondrites (light gray circles) and nuggets in cosmic spherules (dark gray circles). Pt-dominant TeH nuggets in group #1 are a poor match, but Fe-dominant TeH splatter is an excellent match with chondritic meteorites and cosmic spherules, suggesting that they may be extraterrestrial in origin and that the impactor may have been a chondrite. (c) TeH nuggets (red diamonds) are a poor match for most nuggets in iron meteorites (purple circles), but an excellent match for nuggets found in comets (green circles). These data suggest that Fe-dominant PGE nuggets at TeH may have originated from cometary material. (d) Semi-log comparison of PGEs ruthenium (Ru), rhodium (Rh), palladium (Pd), osmium (Os), iridium (Ir), and platinum (Pt), normalized to CI chondrites. TeH Fe-dominant splatter (red line) is an excellent match for PGE nuggets in carbonaceous chondrites (blue line), cosmic spherules (purple line), micrometeorites (dark blue line), and iron meteorites (gray line). In contrast, TeH PGE nuggets are a poor match for bulk material from CI-normalized CV-type chondrites (e.g., Allende; orange line) and CM-type chondrites (e.g., Murchison; brown line).Full size imageFe-dominant splatterWe made 8 measurements on TeH Fe-dominant PGE splatter (Fig. 34d–f). The metal-rich areas average ~ 318 µm in length (range 20–825 µm) with an estimated concentration of 1 PGE-rich bleb per mm2, 100× more common than the TeH nuggets. Average concentrations are Fe = 17.5 wt.%, Ir = 4.7 wt.%, and Pt = 1.5 wt.%.We explored a potential extraterrestrial origin by constructing ternary diagrams for comparison of TeH Fe-dominant splatter with known meteorites and comets (Fig. 36b, c). We compiled data for 164 nuggets extracted from carbonaceous chondritic meteorites (e.g., Allende, Murchison, Leoville, and Adelaide)119,120,121,122, seafloor cosmic spherules123,124, iron meteorites122,125, Comet Wild 2126, and cometary dust particles126. For average weight percentages, see Supporting Information, Tables S6, S7. The Fe-dominant TeH splatter (Fig. 36b) closely matches nuggets from carbonaceous chondrites and cosmic spherules but is a weak match for most iron meteorites (Fig. 36c). In addition, the TeH nuggets are similar to four cometary particles, two of which were collected during the Stardust flyby mission of Comet Wild 2 in 2004126. For average weight percentages, see Supporting Information, Tables S6, S7.To further explore an extraterrestrial connection for TeH Fe-dominant splatter, we compiled wt.% data for TeH PGEs (Rh, Ru, Pd, Os, Ir, and Pt) and normalized them to CI chondrites using values from Anders and Grevasse127. We compared those values to CI-normalized nuggets in carbonaceous chondrites, including CV-type chondrites (e.g., Allende) and CM types (e.g., Murchison)119,120,122,128,129,130,131, seafloor cosmic spherules124, micrometeorites123, and iron meteorites122,125. These results are shown in Fig. 36d.The TeH Fe-dominant splatter closely matches all types of extraterrestrial material with a similar pattern among all data sets: Pd has the lowest normalized values and Os and/or Ir have the highest, closely followed by Pt. The TeH splatter was also compared to the CI-normalized wt.% of bulk meteoritic material from CV- and CM-type chondrites (Fig. 36d). The composition of TeH splatter shows poor correlation with bulk chondritic materials, although the splatter is an excellent geochemical match with the PGE nuggets inside them. In summary, the CI normalization of PGEs suggests an extraterrestrial origin for the Fe-dominant TeH splatter, just as the ternary diagrams also suggest an extraterrestrial source. The correspondence of these two independent results suggests that the quantification of PGEs is sufficiently accurate in this study.Another unusually abundant element, Mo, is also associated with Fe-dominant splatter but not with Pt-dominant nuggets. Mo averages 0.3 wt.% with up to 1.1 wt.% detected in Fe-dominant splatter but with none detected in TeH Pt-dominant nuggets. Mo also is not reported in any terrestrial placer nuggets and occurs in low concentrations (less than ~ 0.02 wt.%) in iron meteorites. In contrast, Mo is reported at high concentrations in PGE nuggets from carbonaceous chondrites (~ 11.5 wt.%), cosmic spherules (0.6 wt.%), and cometary material (5.8 wt.%). Thus, the Mo content of TeH splatter appears dissimilar to terrestrial material but overlaps values of known cosmic material, suggesting an extraterrestrial origin.Based on the volume and weight of the meltglass, we estimate that the extraterrestrial-like metallic TeH Fe-dominant splatter represents  More

  • in

    Confirmation of Oryctes rhinoceros nudivirus infections in G-haplotype coconut rhinoceros beetles (Oryctes rhinoceros) from Palauan PCR-positive populations

    Insects and virusOryctes rhinoceros was collected from Amami, Kagoshima, Japan in 2017 and Ishigaki, Okinawa, Japan in 2018. The insects were brought back to the lab in Tokyo and maintained in a moisture mushroom mat substrate (Mushroom Mat, Tsukiyono Kinokoen, Japan) which was also served as food for larvae. The temperature was held at 25–30 °C with a 16-h light / 8-h dark photoperiod. To collect eggs, 2 or 3 female adults were put in a plastic case containing a moisture mushroom mat substrate with a male adult beetle. The insect jelly (Dorcus Jelly, Fujikon, Japan) was provided ad libitum as food for adults. After 2 weeks, we collected eggs, and about 10 eggs were placed in a plastic cup with a moisture mushroom mat substrate until hatched larvae developed to the second instar. This strain was used in all bioassays in this study. All Japanese O. rhinoceros were confirmed as CRB-G.The OrNV-X2B isolate used in this study was originally isolated from Philippine CRB and obtained from AgResearch in New Zealand.Cell culturesFRI-AnCu-35 (AnCu35) cells were obtained from Genebank of NARO (Tsukuba, Japan)27. This continuous cell line was developed from embryos of the cupreous chafer, Anomala cuprea (Coleptera: Scarabaeidae). The cells were maintained as adherent cultures in 25 cm2 tissue culture flasks (Falcon, Corning, USA) at 25 °C in 5 ml of 10% Fetal Bovine Serum (Gibco, Thermo Fisher Scientific, USA) supplemented Grace’s insect medium (Gibco). Cells were passaged in the above culture medium until the cell monolayer reached 70% confluence.DNA extraction and identification of haplotypes in Palauan populationCRB specimens were collected in Palau using pheromone traps containing ethyl 4-methyloctanoate (ChemTica Internacional, Costa Rica). Adults were dissected to collect midgut and gut tissues to avoid cross contamination between dissection of individuals, which were immediately soaked into 0.1 μg/ml gentamicin solution to prevent bacterial contamination during transportation at room temperature. Specimens were stored at − 30 °C after arrival to Tokyo. The tissues were homogenized in cell lysis solution (10 mM Tris–HCl, 100 mM EDTA, 1% SDS, pH 8.0) using pestles in 1.5 ml microcentrifuge tubes. Homogenates were centrifuged at 12,000× g for 5 min at 4 °C. Proteinase K (200 µg/ml final concentration) (Nippon Gene Co. Ltd., Japan) was added to the supernatant and incubated at 50 °C for 5 h. To remove contaminating RNA, RNase A solution (100 µg/ml final concentration) (Nippon Gene Co. Ltd.) was added. After a 30 min incubation at 37 °C, the mixture was placed on ice and supplemented with 200 μl of Protein Precipitation Solution (Qiagen, Germany), and then centrifuged at 17,000× g for 15 min at 4 °C. The supernatant was isopropanol-precipitated, pelleted by centrifugation, and washed with 70% ethanol. Finally, precipitated DNA was dissolved in distilled MilliQ water. The concentrations of each DNA solution were measured by using NanoVue Plus (GE Healthcare, Buckinghamshire, England, UK). The sample DNA was diluted to 10 ng/μl and used for PCR. The following primer pair was used to amplify a 523 bp fragment of the COI gene: C1-J-1718Oryctes (5′-GGAGGTTTCGGAAATTGACTTGTTCC-3′) and C1-N-2191Oryctes (5′-CCAGGTAGAATTAAAATRTATACCTC-3′)9. Each 10 μl PCR reaction contained: 5 μl Emerald Amp (Takara, Japan), 0.3 μl forward primer (10 μM), 0.3 μl reverse primer (10 μM), 3.4 μl Milli-Q water (Merck Millipore, USA), and 1 µl template DNA. PCR amplifications were performed in a Life ECO thermocycler (Bioer Technology, China) with a cycling profile of 35 cycles of 94 °C denaturation (30 s), 50 °C annealing (45 s), 72 °C extension (1 min) with an initial denaturation of 3 min at 94 °C and a final extension of 5 min at 72 °C. A 5 μl aliquot of each PCR amplicon was checked by agarose gel electrophoresis (1.5%, 1 × TBE), stained with Midori green (Nippon Genetics, Japan) and fluorescence visualized over UV light. Photographs were recorded using an E-BOX-VX2 /20 M (E & M, Japan).For direct sequencing, the PCR products were purified using a QIAquick PCR Purification Kit (Qiagen). The purified DNA was sequenced using BigDye Terminator Kit ver. 3.1 (Applied Biosystems, USA) and performed by the 3700 DNA analyzer (Applied Biosystems). The obtained sequences were analyzed using MEGA X software28 and the G haplotype was identified by the presence of the (A→G) point mutation in the COI region as previously described9.Virus detection in Palauan populationUsing the same samples as above, virus detection was carried out by PCR. The following primer pair was used to amplify a 944 bp fragment of the OrNV-gp054 gene (GrBNV-gp83-like protein): OrNV15a (5′-ATTACGTCGTAGAGGCAATC-3′) and OrNV15b (5′-ATGATCGATTCGTCTATGG-3′)29. PCR amplifications were performed as above.Transmission electron microscopy (TEM) was also used for detection of OrNV within a subset of PCR positive CRB tissue samples. After washing in phosphate-buffered saline (PBS), midgut and fat body samples of Palauan CRB adults from Melekeok and Aimeliik (respectively; two each), were subjected to following resin fixation as described previously30: tissues were fixed in 5% glutaraldehyde for 1 h, rinsed 4 times with Millonig’s phosphate buffer (0.18% NaH2PO4 × H2O, 2.33% Na2HPO4 × 7H2O, 0.5% NaCl, pH 7.4), post-fixed and stained in 1% OsO4 for 2 h and dehydrated in an ethanol series. Following the final dehydration step, the ethanol was replaced by QY-1 (Nisshin EM, Tokyo), and the tissues were embedded in epoxy resin comprising 47% TAAB EPON812, 19% DDSA, 32% MNA and 2% DMP30 (Nisshin EM, Tokyo). Then, they were cut into 70 nm thick sections with a diamond knife on an Ultracut N ultramicrotome (Leica, Vienna, Austria), attached to grids and observed using TEM (JEM-1400Plus, JEOL, Japan).Isolation of OrNV from Palauan samples and infectivity to Japanese CRB larvaeVirus isolation was carried out using a modification of a method previously described23. The frozen tissues of two virus positive CRB-G from Melekeok were washed with PBS twice, and after grounding with 1 ml PBS by pestles, centrifuged at 6,000 g × 5 min at 4 °C. The supernatant was filtered by 0.45 µm pore sized filter (Merck, USA) and transferred to a 1.5 ml ultracentrifuge tube in a clean bench. Virus was pelleted by centrifugation at 4 °C, 98,600 g for 30 min using a TLA55 rotor. After separation, the supernatant was discarded and the pellet was suspended in 500 μl of PBS and designated as “virus solution”. A portion of this solution (30 µl/larva) was intrahemocoelically injected into 82nd instar CRB to evaluate its infectivity. This experiment had no biological replicates due to the very small amount of inoculum available. Intrahemocoelically injected larvae were reared in the insect rearing mat at 25 °C for two weeks. Following death, larval cadavers were immediately dissected to collect midgut for following RNA extraction to detect expression of a viral gene, and electron microscopy observation. Total RNA was extracted from larval tissue samples using ISOGEN (Nippon Gene Co. Ltd., Tokyo, Japan), as described in the manufactural protocol. The total RNA samples were treated with RNAse-free recombinant DNAse I (TaKaRa, Japan) to remove the contaminating DNAs. The DNAse I treated total RNA samples (approximately 100 ng/µl) were used as templates for cDNA synthesis using a TaKaRa RNA PCR Kit (AMV) ver. 3.0 (TaKaRa, Japan). PCR reactions were conducted as above using OrNV15a and b primers (detects gene GrBNV-gp83-like gene). This experiment was conducted in triplicate.Inoculum preparation using FRI-AnCu-35 cellsOrNV isolates were propagated using the FRI-AnCu-35 (AnCu35) cell line for further analyses following methods previously described for the DSIR-Ha-1179 cell line system9,12. AnCu35 was a Coleopteran cell line readily available in Japan, and was inoculated with the Palau OrNV solution prepared above and the OrNV-X2B isolate which was provided by AgResearch, New Zealand. When the cell culture reached 25% confluency, a 100 µl aliquot of virus solution was inoculated and incubated at 25 °C. The virus-treated cells were observed by optical microscope.Quantification of viral copy number using qPCR was conducted as follows. To measure the amount of OrNV virus produced by the AnCu35 cell line, DNA was extracted as described above for tissue samples from 1.5 ml of the virus treated cell’s suspension at 10 dpi (3 suspensions per each virus isolate). The extracted DNA was subjected to quantitative PCR (qPCR) following previously described methods31. The primer pair for qPCR was designed from the genome sequence of the P74 homolog of OrNV, a viral structural protein that is conserved widely among nudiviruses, polydnaviruses and baculoviruses32, to amplify a region of 82 bp of OrNV-X2B-gp120 (OrNV-p74_f2026: 5′-ATCGCCGGTGTGTTTATGG-3′, OrNV-p74_r2107: 5′-AGAGGGCTAACGCTACGAC-3′). The qPCR reaction was performed by using Step One Plus Real-Time PCR System (Life Technologies, USA). The reaction mixture contained 10 ng of template DNA, 5 µl of FastStart Universal SYBR Green Master Mix (ROX) (Roche, Switzerland), 0.3 µl forward primer (10 µM), 0.3 µl reverse primer (10 µM), and 3.4 µl Milli-Q water. The qPCR cycle condition was as follows: 95 °C 10 min; 40 cycle of 95 °C 15 s, 60 °C 1 min. At the end of the cycles, a dissociation curve analysis of the amplified product was performed as follows: 95 °C 15 s, 60 °C 1 min, 95 °C 15 s. The Ct value of each sample DNA was measured twice using two wells as technical replicates. The quantity of the viral genome (ng) in each sample was calculated from a standard curve generated from 29.7 to 29.7 × 10–5 ng of purified PCR amplicon from the OrNV P74 gene. The viral copies in 1 ng of sample DNA was estimated from the molecular weight of qPCR target region (p74). The virus titer was determined from average copy numbers of three virus suspensions as follows. The p74 qPCR amplicon was 83 bp, and the molecular weight of the amplicon was calculated as the length of dsDNA (83 bp) × 330 daltons × 2 nt/bp = 54,780 daltons (g/mol). DNA weight of 1 copy of virus genome was calculated as 54,780 g/mol/Avogadro constant (6.023 × 1023 molecules/mol) = 9.095 × 10–20 g/ molecule. Amplicons of the above region was purified by QIA quick PCR purification kit (Qiagen) and 29.7 ng/ul of DNA was obtained for use as a quantification standard. This is equivalent to 3.266 × 1011 copies of p74 gene (because the amplicon is 9.095 × 10–20 g/copy). Based on qPCR using the serial dilutions (× 10 – 105) of the standard DNA prepared above, Ct values were examined by each concentration of viral DNA. Ct-value = − 3.3112x – 1.4219 (x: diluton factor of 10x). Accordingly, copy number of p74 = 3.266 × 1011+x. Viral copy number (copy number of p74 genes) was calculated from Ct-value from the above formula.Viral replication in CRB larvae by time course and killing speedField collected CRB-G larvae from Japan were inoculated with the OrNV-Palau1 and -X2B isolates to examine establishment of infection over time using qPCR. The inoculum was prepared from supernatant collected from OrNV infected AnCu35 cell cultures at 10 dpi, passed through a 0.45 µm filter, and preserved at 4 °C until use.Second instar CRB was inoculated intrahemocoelically with 30 μl of the virus solution prepared from cell-culture per larva using a microinjector (Kiya Kogyo Seisakusho, Japan) fitted with a micro-syringe (Ito Seisakusho, Japan). The virus doses of OrNV-Palau1 and -X2B strains used for inoculation were confirmed to be comparable by absolute quantification using the above qPCR method (Palau1: 3.1 × 105 copies/ng, X2B: 3.3 × 105copies/ng; the mean titer of 3 DNA templates, respectively). As a mock treatment, CRB was injected with 30 µl PBS. The inoculated larvae were kept individually in plastic containers with a rearing mat in a 25 °C incubator. The samples were collected at 3, 6, and 9 dpi (25–30 larvae per time point) into 15 ml tubes and stored at − 30 °C until the DNA was extracted as above. Total DNA was extracted from whole, individual larvae which were dissected to remove midgut contents to prevent interference to Taq polymerase, and subjected to qPCR as above. Changes in viral copy number within the same virus strain over time were analyzed by one-way, nonparametric Steel–Dwass tests using JMP@ 9.0.0 software (SAS Institute, Cary, NC). Differences in virus copy number between strains were analyzed in the same way, but to correct for errors in the test values due to multiple comparisons, Bonferroni’s correction was used to set the α-value for the test at 0.008333. Ten larvae were inoculated and examined per each treatment-time point with three replications.To estimate killing speed, CRB-G larvae from Japan were inoculated with the OrNV-Palau1 and -X2B isolates as described previously. Intrahemocoelically inoculated larvae were reared individually in plastic containers with a rearing mat in a 25 °C incubator. Mortality of inoculated larvae were observed every day. Forty larvae were examined in a replicate with three replications carried out for virus treatments (total 120 larvae). The mock PBS inoculation treatment was done only once (total 37 larvae).Genome sequencingGenome sequencing of the OrNV-Palau1 isolate and X2B isolate was conducted. For obtaining high quality DNA, virus particles were purified, from 3 mL of AnCu35 culture supernatant collected six days after inoculation with OrNV. Virus containing supernatant was transferred to Ultra-Clear polyallomer tubes (Beckman Coulter, USA) with a 20–50% (w/w) sucrose density gradient and subjected to ultracentrifugation at 72,100 g, 4 °C, for 1 h. After ultracentrifugation, the white virus band was collected in a 1.5 ml tube. The solution was then subjected to ultracentrifugation at 110,000 g, 4 °C for 1 h to precipitate the viral particles33. Then, DNA was extracted from purified OrNV virions as described above. For the sequencing analysis, DNA libraries were prepared using the Nextera XT DNA Library Prep Kit (Illumina, USA). Amplified libraries were sequenced on Illumina HiSeq 2500 instrument using paired-end 2 × 150 bp chemistry which was performed by Novogene (Beijing, China). Contigs of each strain from NGS reads were generated by assembly using Unicycler (version 0.4.8)34. The gaps between contigs were further closed with Sanger sequences obtained by PCR direct sequencing using appropriate specific primers, and the sequence was aligned by minimap2 (version 2.17)35. The assembly and sequences of contigs were also confirmed by mapping to the OrNV isolate Solomon Islands genome sequence (GenBank accession no. MN623374.1) with NGS reads and Sanger sequences using minimap2. The mapped reads (SAM files) were converted to BAM format using SAMtools (version 1.10)36. After the sorting and indexing of BAM files, the consensus sequences were generated using bcftools (version 1.10.2)37.ORFs of at least 50 codons in size that possessed significant amino acid sequence similarity with ORFs from OrNV-Ma07 were identified with Lasergene GeneQuest (DNAStar, v. 17) and BLASTp. ORFs with no significant matches to other sequences also were selected for annotation if (a) they did not overlap a larger ORF by  > 75 bp, and (b) they were predicted to be protein-encoding by both the fgenesV0 (http://www.softberry.com/berry.phtml?topic=index&group=programs&subgroup=gfindv) and Vgas38 programs.OrNV genome sequences were compared by pairwise alignment using the Martinez/Needleman-Wunsch method as implemented in Lasergene MegAlign 15. Pairwise sequence identities were determined from these alignments as previously described39. Differences in ORF content and distribution of selected OrNV genomic regions were visualized with Mauve version 2015022640.Phylogenetic inferenceTo infer the relationships among OrNV isolates on the basis of nucleotide sequence alignments, the DNA polymerase ORFs of completely sequenced isolates (Table 2), OrNV-PV50516, and a set of nine isolates from Indonesia17 were aligned by MUSCLE as implemented in Lasergene MegAlign Pro v. 17 (DNAStar). Phylogeny was inferred by maximum likelihood using MEGA X28 with the Tamura-Nei (TN93) model41, with ambiguous data eliminated prior to analysis. Tree reliability was evaluated by bootstrap with 500 replicates. More

  • in

    “Indirect development” increases reproductive plasticity and contributes to the success of scyphozoan jellyfish in the oceans

    1.Cartwright, P. et al. Exceptionally preserved jellyfishes from the middle Cambrian. PLoS One 2, e1121 (2007).ADS 
    Article 

    Google Scholar 
    2.Walcott, C. D. Cambrian Geology and Paleontology II: No. 3—Middle Cambrian Holothurians and Medusae Vol. 3 (Smithsonian Institution, 1911).
    Google Scholar 
    3.Willoughby, R. H. & Robison, R. A. Medusoids from the Middle Cambrian of Utah. J. Paleontol. 53, 494–500 (1979).
    Google Scholar 
    4.Rigby, S. & Milsom, C. V. Origins, evolution, and diversification of zooplankton. Annu. Rev. Ecol. Syst. 31, 293–313 (2000).Article 

    Google Scholar 
    5.Young, G. A. & Hagadorn, J. W. The fossil record of cnidarian medusae. Palaeoworld 19, 212–221 (2010).Article 

    Google Scholar 
    6.Technau, U. & Steele, R. E. Evolutionary crossroads in developmental biology: Cnidaria. Development 138, 1447 (2012).Article 

    Google Scholar 
    7.Hoegh-Guldberg, O., Poloczanska, E. S., Skirving, W. & Dove, S. Coral reef ecosystems under climate change and ocean acidification. Front. Mar. Sci. https://doi.org/10.3389/fmars.2017.00158 (2017).Article 

    Google Scholar 
    8.Hagadorn, J. W., Dott, R. H. & Damrow, D. Stranded on a Late Cambrian shoreline: Medusae from central Wisconsin. Geology 30, 147–150 (2002).ADS 
    Article 

    Google Scholar 
    9.Boero, F. Review of jellyfish blooms in the Mediterranean and Black Sea. Studies and Reviews. General Fisheries Commission for the Mediterranean, Vol. 92 (FAO, Rome, 2013).10.Brotz, L., Cheung, W., Kleisner, K., Pakhomov, E. & Pauly, D. Increasing jellyfish populations: Trends in large marine ecosystems. Hydrobiologia 690, 3–20 (2012).Article 

    Google Scholar 
    11.Condon, R. H. et al. Recurrent jellyfish blooms are a consequence of global oscillations. Proc. Natl. Acad. Sci. 110, 1000–1005. https://doi.org/10.1073/pnas.1210920110 (2013).ADS 
    Article 
    PubMed 

    Google Scholar 
    12.Arai, M. Pelagic coelenterates and eutrophication: A review. Hydrobiologia 451, 69–87. https://doi.org/10.1023/A:1011840123140 (2001).Article 

    Google Scholar 
    13.Purcell, J. E., Malej, A. & Benović, A. in Ecosystems at the Land-Sea Margin: Drainage Basin to Coastal Seas Vol. 55 Coastal and Estuarine Studies Ch. 8, 241–263 (American Geophysical Union, 1999).14.Lynam, C. P. et al. Have jellyfish in the Irish Sea benefited from climate change and overfishing?. Glob. Change Biol. 17, 767–782. https://doi.org/10.1111/j.1365-2486.2010.02352.x (2011).ADS 
    Article 

    Google Scholar 
    15.Richardson, A. J., Bakun, A., Hays, G. C. & Gibbons, M. J. The jellyfish joyride: Causes, consequences and management responses to a more gelatinous future. Trends Ecol. Evol. 24, 312–322 (2009).Article 

    Google Scholar 
    16.Lucas, C. H., Graham, W. M. & Widmer, C. Jellyfish life histories: Role of polyps in forming and maintaining scyphomedusa populations. Adv. Mar. Biol. 63, 133–196 (2012).Article 

    Google Scholar 
    17.Helm, R. R. Evolution and development of scyphozoan jellyfish. Biol. Rev. 93, 1228–1250 (2018).Article 

    Google Scholar 
    18.Jarms, G. & Morandini, A. C. World Atlas of Jellyfish (Dölling und Galitz Verlag, Germany, 2019).
    Google Scholar 
    19.Piraino, S., Boero, F., Aeschbach, B. & Schmid, V. Reversing the life cycle: medusae transforming into polyps and cell transdifferentiation in Turritopsis nutricula (Cnidaria, Hydrozoa). Biol. Bull. 180, 302–312 (1996).Article 

    Google Scholar 
    20.De Vito, D., Piraino, S., Schmich, J., Bouillon, J. & Boero, F. Evidence of reverse development in Leptomedusae (Cnidaria, Hydrozoa): the case of Laodicea undulata (Forbes and Goodsir 1851). Mar. Biol. 149, 339–346 (2006).Article 

    Google Scholar 
    21.He, J., Zheng, L., Zhang, W. & Lin, Y. Life cycle reversal in Aurelia sp.1 (Cnidaria, Scyphozoa). PLoS One 10, e0145314 (2015).Article 

    Google Scholar 
    22.Sandrini, L. R. & Avian, M. Biological cycle of Pelagia noctiluca: Morphological aspects of the development from planula to ephyra. Mar. Biol. 74, 169–174. https://doi.org/10.1007/BF00413920 (1983).Article 

    Google Scholar 
    23.Jarms, G., Båmstedt, U., Tiemann, H., Martinussen, M. B. & Fosså, J. H. The holopelagic life cycle of the deep-sea medusa Periphylla periphylla (Scyphozoa, Coronatae). Sarsia 84, 55–65 (1999).Article 

    Google Scholar 
    24.Dawson, M. N. & Hamner, W. M. A character-based analysis of the evolution of jellyfish blooms: Adaptation and exaptation. Hydrobiologia 616, 193–215. https://doi.org/10.1007/s10750-008-9591-x (2009).Article 

    Google Scholar 
    25.Ceh, J., Gonzalez, J., Pacheco, A. S. & Riascos, J. M. The elusive life cycle of scyphozoan jellyfish—Metagenesis revisited. Sci. Rep. 5, 12037. https://doi.org/10.1038/srep12037. http://www.nature.com/srep/2015/150708/srep12037/abs/srep12037.html#supplementary-information (2015).26.Campos, L., Gonzállez, K. & Ceh, J. First report of a precocious form of strobilation in a jellyfish, the South American Pacific sea nettle Chrysaora plocamia. Mar. Biodivers. 50, 85 (2020).Article 

    Google Scholar 
    27.Henroth, L. & Grondähl, F. On the biology of Aurelia aurita (L.) 1. Release and growth of Aurelia aurita (L.) ephyrae in the Gullmar Fjiord, western Sweden, 1982–83. Ophelia 22, 189–199 (1983).Article 

    Google Scholar 
    28.Hirai, E. On the developmental cycles of Aurelia aurita and Dactylometra pacifica. Bull. Mar. Biol. Stn Asamushi IX, 81 (1958).
    Google Scholar 
    29.Kakinuma, Y. An experimental study of the life cycle and organ differentiation of Aurelia aurita Lamarck. Bull. Mar. Biol. Stn. Asamushi XV, 101–113 (1975).
    Google Scholar 
    30.Yasuda, T. Ecological studies on the jelly-fish, Aurelia aurita, in Urazoko Bay, Fukui Prefecture-XI. An observation on ephyra formation. Publ. Seto Mar. Biol. Lab. XXII, 75–80 (1975).Article 

    Google Scholar 
    31.Suzuki, K. S. et al. Seasonal alternation of the ontogenetic development of the moon jellyfish Aurelia coerulea in Maizuru Bay, Japan. PLoS One 14, e0225513. https://doi.org/10.1371/journal.pone.0225513 (2019).CAS 
    Article 
    PubMed 
    PubMed Central 

    Google Scholar 
    32.Avian, M. In Workshop on Jellyfish in the Mediterranean Sea Vol. 2 (eds Rottini Sandrini, L. & Avian, M.) 47–59 (Nova Thalassia, 1986).
    Google Scholar 
    33.Costello, J. et al. Project Meduza in the context of its historical time. Ann. Ser. Hist. Nat. 19, 1–18 (2009).
    Google Scholar 
    34.Margiotta, F. et al. Do plankton reflect the environmental quality status? The case of a post-industrial Mediterranean Bay. Mar. Environ. Res. 160, 104980 (2020).CAS 
    Article 

    Google Scholar 
    35.Schiariti, A. et al. Asexual reproduction strategies and blooming potential in Scyphozoa. Mar. Ecol. Prog. Ser. 510, 241–253 (2014).ADS 
    Article 

    Google Scholar 
    36.Yasuda, T. Ecological studies on the jelly-fish, Aurelia aurita, in Urazoko Bay, Fukui Prefecture-IV. Monthly change in the bell-length composition and breeding season. Bull. Jpn. Soc. Sci. Fish. 37, 364–370 (1971).Article 

    Google Scholar 
    37.Suryan, R. M. et al. Environmental forcing on life history strategies: Evidence for multi-trophic level responses at ocean basin scales. Prog. Oceanogr. 81, 214–222 (2009).ADS 
    Article 

    Google Scholar 
    38.Dawson, M. N. Macro-morphological variation among cryptic species of the moon jellyfish, Aurelia (Cnidaria: Scyphozoa). Mar. Biol. 143, 369–379 (2003).Article 

    Google Scholar 
    39.Benović, A. et al. Ecological characteristics of the Mljet Island seawater lakes (South Adriatic Sea) with special reference to their resident population of medusae. Sci. Mar. 64, 197–206 (2000).Article 

    Google Scholar 
    40.Prieto, L., Astorga, D., Navarro, G. & Ruiz, J. Environmental control of phase transition and polyp survival of a massive-outbreaker jellyfish. PLoS One 5, e13793. https://doi.org/10.1371/journal.pone.0013793 (2010).ADS 
    CAS 
    Article 
    PubMed 
    PubMed Central 

    Google Scholar 
    41.Purcell, J. et al. Temperature effects on asexual reproduction rates of scyphozoan polyps from the NW Mediterranean Sea. Hydrobiologia 690, 169–180 (2012).CAS 
    Article 

    Google Scholar 
    42.Kikinger, R. Cotylorhiza tuberculata (Cnidaria: Scyphozoa)—Life history of a stationary population. PSZN Mar. Ecol. 13, 333–362 (1992).Article 

    Google Scholar 
    43.Djeghri, N., Pondaven, P., Stibor, H. & Dawson, M. N. Review of the diversity, traits, and ecology of zooxanthellate jellyfishes. Mar. Biol. 166, 147 (2019).Article 

    Google Scholar 
    44.Glynn, P. W. & Colgan, M. W. Sporadic disturbances in fluctuating coral reef environments: El Niño and coral reef development in the Eastern Pacific. Am. Zool. 32, 707–718. https://doi.org/10.1093/icb/32.6.707 (1999).Article 

    Google Scholar  More

  • in

    Study on hyperspectral estimation model of soil organic carbon content in the wheat field under different water treatments

    1.Aryal, D. R., De Jong, B. H., Ochoa-Gaona, S., Esparza-Olguin, L. & Mendoza-Vega, J. Carbon stocks and changes in tropical secondary forests of southern Mexico. Agr. Ecosyst. Environ. 195, 220–230 (2014).Article 

    Google Scholar 
    2.Aryal, D. R., De Jong, B. H., Ochoa-Gaona, S., Mendoza-Vega, J. & Esparza-Olguin, L. Successional and seasonal variation in litterfall and associated nutrient transfer in semi-evergreen tropical forests of SE Mexico. Nutr. Cycl. Agroecosys. 103(1), 45–60 (2015).CAS 
    Article 

    Google Scholar 
    3.Aryal, D. R. et al. Soil organic carbon depletion from forests to grasslands conversion in Mexico: A review. Trop. Agric. 8, 181 (2018).CAS 

    Google Scholar 
    4.Gao, W., Yang, J., Ren, S. R. & Liu, H. L. The trend of soil organic carbon, total nitrogen, and wheat and maize productivity under different long-term fertilizations in the upland fluvo-aquic soil of North China. Nutr. Cycl. Agroecosys. 103, 61–73 (2015).CAS 
    Article 

    Google Scholar 
    5.Qi, H., Paz-Kagan, T., Karnieli, A., Jin, X. & Li, S. Evaluating calibration methods for predicting soil available nutrients using hyperspectral VNIR data. Soil Till Res. 175, 267–275 (2018).Article 

    Google Scholar 
    6.Dong, X., Tian, J., Zhang, R., He, D. & Chen, Q. Study on the relationship between soil emissivity spectra and content of soil element. Spectrosc. Spect. Anal. 37(02), 557–564 (2017).CAS 

    Google Scholar 
    7.Kemper, T. & Sommer, S. Estimate of heavy metal contamination in soils after a mining accident using reflectance spectroscopy. Environ. Sci. Technol. 36(12), 2742–2747 (2002).ADS 
    CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    8.Panigrahi, N. & Das, B. S. Canopy spectral reflectance as a predictor of soil water potential in rice. Water Resour. Res. 54(4), 2544–2560 (2018).ADS 
    Article 

    Google Scholar 
    9.Peddle, D. R., White, H. P., Soffer, R. J., Miller, J. R. & LeDrew, E. F. Reflectance processing of remote sensing spectroradiometer data. Comput. Geoences. 27(2), 203–213 (2001).ADS 

    Google Scholar 
    10.Ben-Dor, E. et al. Using imaging spectroscopy to study soil properties. Remote Sens. Environ. 113, S38–S55 (2009).Article 

    Google Scholar 
    11.Rossel, R. A., Walvoort, D. J., Mcbratney, A. B., Janik, L. J. & Skjemstad, J. O. Visible, near infrared, mid infrared or combined diffuse reflectance spectroscopy for simultaneous assessment of various soil properties. Geoderma 131(1), 59–75 (2006).ADS 
    Article 
    CAS 

    Google Scholar 
    12.Cheng, H. et al. Estimating heavy metal concentrations in suburban soils with reflectance spectroscopy. Geoderma 336, 59–67 (2019).ADS 
    CAS 
    Article 

    Google Scholar 
    13.Ding, J., Yang, A., Wang, J., Sagan, V. & Yu, D. Machine-learning-based quantitative estimation of soil organic carbon content by VIS/NIR spectroscopy. PeerJ 6(3), e5714 (2018).PubMed 
    PubMed Central 
    Article 
    CAS 

    Google Scholar 
    14.Gobrecht, A., Bendoula, R., Roger, J.-M. & Bellon-Maurel, V. A new optical method coupling light polarization and vis–NIR spectroscopy to improve the measurement of soil carbon content. Soil Till Res. 155, 461–470 (2016).Article 

    Google Scholar 
    15.Gu, X., Wang, Y., Song, X. & Xu, X. The Inversion Model of Soil Organic Matter of Cultivated Land Based on Hyperspectral Technology. Remote Sensing for Agriculture, Ecosystems, and Hydrology XVII (International Society for Optics and Photonics, 2015).
    Google Scholar 
    16.Nawar, S., Buddenbaum, H., Hill, J., Kozak, J. & Mouazen, A. M. Estimating the soil clay content and organic matter by means of different calibration methods of vis-NIR diffuse reflectance spectroscopy. Soil Till Res. 155, 510–522 (2016).Article 

    Google Scholar 
    17.Yu, X., Liu, Q., Wang, Y., Liu, X. & Liu, X. Evaluation of MLSR and PLSR for estimating soil element contents using visible/near-infrared spectroscopy in apple orchards on the Jiaodong peninsula. CATENA 137, 340–349 (2016).CAS 
    Article 

    Google Scholar 
    18.Ji, W. J., Li, X., Li, C. X., Zhou, Y. & Shi, Z. Using different data mining algorithes to predict soil organic matter based on visible-near infrared spectroscopy. Spectrosc. Spect. Anal. 32(09), 2393–2397 (2012).CAS 

    Google Scholar 
    19.Douglas, R. K., Nawar, S., Alamar, M. C., Mouazen, A. M. & Coulon, F. Rapid prediction of total petroleum hydrocarbons concentration in contaminated soil using vis-NIR spectroscopy and regression techniques. Sci. Total Environ. 616, 147–155 (2018).ADS 
    PubMed 
    Article 
    CAS 
    PubMed Central 

    Google Scholar 
    20.Mouazen, A. M. & Al-Asadi, R. A. Influence of soil moisture content on assessment of bulk density with combined frequency domain reflectometry and visible and near infrared spectroscopy under semi field conditions. Soil Till Res. 176, 95–103 (2018).Article 

    Google Scholar 
    21.Rossel, R. A. & Behrens, T. Using data mining to model and interpret soil diffuse reflectance spectra. Geoderma 158(1), 46–54 (2010).ADS 
    CAS 
    Article 

    Google Scholar 
    22.Nawar, S. & Mouazen, A. M. Comparison between random forests, artificial neural networks and gradient boosted machines methods of on-line vis-NIR spectroscopy measurements of soil total nitrogen and total carbon. Sensors. 17(10), 2428 (2017).ADS 
    PubMed Central 
    Article 
    CAS 

    Google Scholar 
    23.Wang, J., Chen, Y., Chen, F., Shi, T. & Wu, G. Wavelet-based coupling of leaf and canopy reflectance spectra to improve the estimation accuracy of foliar nitrogen concentration. Agr. Forest Meteorol. 248, 306–315 (2018).ADS 
    Article 

    Google Scholar 
    24.Hong, Y. et al. Combining fractional order derivative and spectral variable selection for organic matter estimation of homogeneous soil samples by vis–NIR spectroscopy. Remote Sens. 10(3), 479 (2018).ADS 
    Article 

    Google Scholar 
    25.Sorenson, P. T. et al. Monitoring organic carbon, total nitrogen, and pH for reclaimed soils using field reflectance spectroscopy. Can J Soil Sci. 97(2), 241–248 (2017).CAS 
    Article 

    Google Scholar 
    26.Gomez, C., Rossel, R. A. V. & Mcbratney, A. B. Soil organic carbon prediction by hyperspectral remote sensing and field vis-NIR spectroscopy: An Australian case study. Geoderma 146(3–4), 403–411 (2008).ADS 
    CAS 
    Article 

    Google Scholar 
    27.Shi, T. Z. et al. Comparison of multivariate methods for estimating soil total nitrogen with visible/near-infrared spectroscopy. Plant Soil. 366(1–2), 363–375 (2013).CAS 
    Article 

    Google Scholar 
    28.Stenberg, B., Rossel, R. A. V., Mouazen, A. M. & Wetterlind, J. Chapter five-visible and near infrared spectroscopy in soil science. Adv. Agron. 107, 163–215 (2010).CAS 
    Article 

    Google Scholar 
    29.Uddin, M. P., Mamun, M. A. & Hossain, M. A. PCA-based feature reduction for hyperspectral remote sensing image classification. IETE Tech. Rev. 5, 1–21 (2020).
    Google Scholar 
    30.Cambule, A. H., Rossiter, D. G., Stoorvogel, J. J. & Smaling, E. M. A. Building a near infrared spectral library for soil organic carbon estimation in the Limpopo National Park Mozambique. Geoderma 183, 41–48 (2012).ADS 
    Article 
    CAS 

    Google Scholar 
    31.Kawamura, K. et al. Vis-NIR spectroscopy and PLS regression with waveband selection for estimating the total C and N of paddy soils in Madagascar. Remote Sens. 9(10), 1081 (2017).ADS 
    Article 

    Google Scholar 
    32.Leone, A. P., Viscarra-Rossel, R. A., Amenta, P. & Buondonno, A. Prediction of soil properties with PLSR and vis-NIR spectroscopy: Application to mediterranean soils from southern Italy. Curr. Anal. Chem. 8(2), 283–299 (2012).CAS 
    Article 

    Google Scholar 
    33.Wang, S., Chen, Y., Wang, M., Zhao, Y. & Li, J. SPA-based methods for the quantitative estimation of the soil salt content in saline-alkali land from field spectroscopy data: A case study from the Yellow River irrigation regions. Remote Sens. 11(8), 967 (2019).ADS 
    CAS 
    Article 

    Google Scholar 
    34.Barnes, E. M. et al. Remote- and ground-based sensor techniques to map soil properties. Photogramm. Eng Rem S. 69(6), 619–630 (2003).Article 

    Google Scholar 
    35.Priori, S. et al. Field-scale mapping of soil carbon stock with limited sampling by coupling gamma-ray and vis-NIR spectroscopy. Soil Sci Soc Am J. 80(4), 954–964 (2016).ADS 
    CAS 
    Article 

    Google Scholar 
    36.Amin, I., Fikrat, F., Mammadov, E. & Babayev, M. Soil organic carbon prediction by vis-NIR spectroscopy: Case study the Kur-Aras plain Azerbaijan. Commun. Soil Sci. Plan. 51(6), 726–734 (2020).CAS 
    Article 

    Google Scholar 
    37.Yu, L. et al. Hyperspectral estimation of soil organic matter content based on partial least squares regression. Trans. CSAE. 31(14), 103–109 (2015).
    Google Scholar 
    38.Liu, Y. F., Lu, Y. N., Guo, L., Xiao, F. T. & Chen, Y. Y. Construction of calibration set based on the land use types in visible and near-infrared (VIS-NIR)model for soil organic matter estimation. Acta Pedol. Sin. 53, 332–341 (2016).
    Google Scholar 
    39.Zhou, X. M. & Zhang, T. Analysis of the April 2019 atmospheric circulation and weather. Meteor. Mon. 45(7), 1028–1036 (2019).
    Google Scholar 
    40.Guan, L. & Zhang, T. Analysis of the May 2019 atmospheric circulation and weather. Meteor. Mon. 45(8), 1181–1188 (2019).
    Google Scholar 
    41.Li, X., He, Y. & Wu, C. Non-destructive discrimination of paddy seeds of different storage age based on Vis/NIR spectroscopy. J. Stored Prod. Res. 44(3), 264–268 (2008).Article 

    Google Scholar 
    42.Boško, M. & Bensa, A. Prediction of soil organic carbon using VIS-NIR spectroscopy: Application to Red Mediterranean soils from Croatia. Eurasian J. Soil Sci. 6(4), 365–373 (2017).
    Google Scholar 
    43.McCarty, G. W., Reeves, J. B. III., Reeves, V. B., Follett, R. F. & Kimble, J. M. Mid-infrared and near-infrared diffuse reflectance spectroscopy for soil carbon measurement. Soil Sci. Soc. Am. J. 66(2), 640–646 (2002).ADS 
    CAS 
    Article 

    Google Scholar 
    44.Gholizadeh, A. et al. Comparing different data preprocessing methods for monitoring soil heavy metals based on soil spectral features. Soil Water Res. 10(4), 218–227 (2015).CAS 
    Article 

    Google Scholar 
    45.Wang, X., Xue, L., He, X. W. & Liu, M. H. Vitamin C content estimation of chilies using Vis/NIR spectroscopy. Int. Conf. Electr. Inf. Control Eng. 2011, 1894–1897 (2011).
    Google Scholar 
    46.Lee, K. S. et al. Wavelength identification and diffuse reflectance estimation for surface and profile soil properties. Am. Soc. Agric. Biol. Eng. 52(3), 683–695 (2009).CAS 

    Google Scholar  More

  • in

    Genome-wide analysis reveals associations between climate and regional patterns of adaptive divergence and dispersal in American pikas

    Alexander DH, Novembre J, Lange K (2009) Fast model-based estimation of ancestry in unrelated individuals. Genome Res 19:1655–1664CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Alexander DH, Shringarpure SS, Novembre J, Lange K (2015) Admixture 1.3 software manual. UCLA Hum Genet Softw Distrib, Los Angeles
    Google Scholar 
    Angert AL, Bontrager MG, Ågren J (2020) What do we really know about adaptation at range edges? Annu Rev Ecol Evol Syst 51:341–361Article 

    Google Scholar 
    Araújo MB, Pearson RG, Thuiller W, Erhard M (2005) Validation of species–climate impact models under climate change. Glob Change Biol 11:1504–1513Article 

    Google Scholar 
    Astle W, Balding DJ (2009) Population structure and cryptic relatedness in genetic association studies. Stat Sci 24:451–471Article 

    Google Scholar 
    Attard CRM, Beheregaray LB, Möller LM (2018) Genotyping-by-sequencing for estimating relatedness in nonmodel organisms: Avoiding the trap of precise bias. Mol Ecol Resour 18:381–390CAS 
    PubMed 
    Article 

    Google Scholar 
    Baird NA, Etter PD, Atwood TS, Currey MC, Shiver AL, Lewis ZA et al. (2008) Rapid SNP discovery and genetic mapping using sequenced RAD markers. PLoS One 3:1–7Article 
    CAS 

    Google Scholar 
    Barbosa S, Mestre F, White TA, Paupério J, Alves PC, Searle JB (2018) Integrative approaches to guide conservation decisions: Using genomics to define conservation units and functional corridors. Mol Ecol 27:3452–3465PubMed 
    Article 

    Google Scholar 
    Beever EA, Brussard PF, Berger J (2003) Patterns of apparent extirpation among isolated populations of pikas (Ochotona princeps) in the Great Basin. J Mammal 84:37–54Article 

    Google Scholar 
    Beever EA, Ray C, Mote PW, Wilkening JL (2010) Testing alternative models of climate-mediated extirpations. Ecol Appl 20:164–178PubMed 
    Article 

    Google Scholar 
    Beever EA, Ray C, Wilkening JL, Brussard PF, Mote PW (2011) Contemporary climate change alters the pace and drivers of extinction. Glob Change Biol 17:2054–2070Article 

    Google Scholar 
    Beever EA, Perrine JD, Rickman T, Flores M, Clark JP, Waters C et al. (2016) Pika (Ochotona princeps) losses from two isolated regions reflect temperature and water balance, but reflect habitat area in a mainland region. J Mammal 97:1495–1511Article 

    Google Scholar 
    Bellard C, Bertelsmeier C, Leadley P, Thuiller W, Courchamp F (2012) Impacts of climate change on the future of biodiversity. Ecol Lett 15:365–377PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Blois JL, Williams JW, Fitzpatrick MC, Jackson ST, Ferrier S (2013) Space can substitute for time in predicting climate-change effects on biodiversity. Proc Natl Acad Sci USA 110:9374–9379CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Browning SR (2008) Missing data imputation and haplotype phase inference for genome-wide association studies. Hum Genet 124:439–450CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Calkins MT, Beever EA, Boykin KG, Frey JK, Andersen MC (2012) Not-so-splendid isolation: modeling climate-mediated range collapse of a montane mammal Ochotona princeps across numerous ecoregions. Ecography 35:780–791Article 

    Google Scholar 
    Carlson SM, Cunningham CJ, Westley PAH (2014) Evolutionary rescue in a changing world. Trends Ecol Evol 29:521–530PubMed 
    Article 

    Google Scholar 
    Castillo JA, Epps CW, Davis AR, Cushman SA (2014) Landscape effects on gene flow for a climate-sensitive montane species, the American pika. Mol Ecol 23:843–856PubMed 
    Article 

    Google Scholar 
    Castillo JA, Epps CW, Jeffress MR, Ray C, Rodhouse TJ, Schwalm D (2016) Replicated landscape genetic and network analyses reveal wide variation in functional connectivity for American pikas. Ecol Appl 26:1660–1676PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    Ceballos G, Ehrlich PR, Raven PH (2020) Vertebrates on the brink as indicators of biological annihilation and the sixth mass extinction. Proc Natl Acad Sci USA 117:13596–13602CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Ceballos G, Ehrlich PR, Barnosky AD, García A, Pringle RM, Palmer TM (2015) Accelerated modern human–induced species losses: entering the sixth mass extinction. Sci Adv 1:e1400253PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Chapman JA, Flux JE (1990) Rabbits, hares and pikas: status survey and conservation action plan. IUCN.Chypre M, Zaidi N, Smans K (2012) ATP-citrate lyase: a mini-review. Biochem Biophys Res Commun 422:1–4CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    Cornuet JM, Luikart G (1996) Description and power analysis of two tests for detecting recent population bottlenecks from allele frequency data. Genetics 144:2001–2014CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Danecek P, Auton A, Abecasis G, Albers CA, Banks E, DePristo MA et al. (2011) The variant call format and VCFtools. Bioinformatics 27:2156–2158CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Diaz HF, Grosjean M, Graumlich L (2003) Climate variability and change in high elevation regions: past, present and future. Clim Change 59:1–4Article 

    Google Scholar 
    Eckert CG, Samis EK, Lougheed SC (2008) Genetic variation across species’ geographical ranges: the central-marginal hypothesis and beyond. Mol Ecol 17:1170–1188CAS 
    PubMed 
    Article 

    Google Scholar 
    Erb LP, Ray C, Guralnick R (2011) On the generality of a climate-mediated shift in the distribution of the American pika (Ochotona princeps). Ecology 92:1730–1735PubMed 
    Article 

    Google Scholar 
    Excoffier L, Foll M, Petit RJ (2009) Genetic consequences of range expansions. Annu Rev Ecol Evol Syst 40:481–501Article 

    Google Scholar 
    Flanagan SP, Forester BR, Latch EK, Aitken SN, Hoban S (2018) Guidelines for planning genomic assessment and monitoring of locally adaptive variation to inform species conservation. Evol Appl 11:1035–1052PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    Foll M, Gaggiotti O (2008) A genome-scan method to identify selected loci appropriate for both dominant and codominant markers: a Bayesian perspective. Genetics 180:977–993PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Franks SJ, Hoffmann AA (2012) Genetics of climate change adaptation. Annu Rev Genet 46:185–208CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    Frichot E, François O (2015) LEA: an R package for landscape and ecological association studies. Methods Ecol Evol 6:925–929Article 

    Google Scholar 
    Frichot E, Schoville SD, Bouchard G, François O (2013) Testing for associations between loci and environmental gradients using latent factor mixed models. Mol Biol Evol 30:1687–1699CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Funk WC, McKay JK, Hohenlohe PA, Allendorf FW (2012) Harnessing genomics for delineating conservation units. Trends Ecol Evol 27:489–496PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Galbreath KE, Hafner DJ, Zamudio KR (2009) When cold is better: climate-driven elevation shifts yield complex patterns of diversification and demography in an Alpine specialist (American Pika, Ochotona Princeps). Evolution 63:2848–2863CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    Gautier M (2015) Genome-wide scan for adaptive divergence and association with population-specific covariates. Genetics 201:1555–1579CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Götz S, García-Gómez JM, Terol J, Williams TD, Nagaraj SH, Nueda MJ et al. (2008) High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res 36:3420–3435PubMed 
    PubMed Central 
    Article 
    CAS 

    Google Scholar 
    Gradogna A, Gavazzo P, Boccaccio A, Pusch M (2017) Subunit‐dependent oxidative stress sensitivity of LRRC8 volume‐regulated anion channels. J Physiol 595:6719–6733CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Hafner DJ, Sullivan RM (1995) Historical and ecological biogeography of Nearctic pikas (Lagomorpha: Ochotonidae). J Mammal 76:302–321Article 

    Google Scholar 
    Hafner DJ, Smith AT (2010) Revision of the subspecies of the American pika, Ochotona princeps (Lagomorpha: Ochotonidae). J Mammal 91:401–417Article 

    Google Scholar 
    Hanson JO, Marques A, Veríssimo A, Camacho-Sanchez M, Velo-Antón G, Martínez-Solano Í et al. (2020) Conservation planning for adaptive and neutral evolutionary processes. J Appl Ecol 57:2159–2169Article 

    Google Scholar 
    Harrisson KA, Pavlova A, Telonis-Scott M, Sunnucks P (2014) Using genomics to characterize evolutionary potential for conservation of wild populations. Evol Appl 7:1008–1025PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Heled J, Drummond AJ (2008) Bayesian inference of population size history from multiple loci. BMC Evol Biol 8:289PubMed 
    PubMed Central 
    Article 
    CAS 

    Google Scholar 
    Henry P, Russello MA (2013) Adaptive divergence along environmental gradients in a climate-change-sensitive mammal. Ecol Evol 3:3906–3917CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Henry P, Sim Z, Russello MA (2012) Genetic evidence for restricted dispersal along continuous altitudinal gradients in a climate change-sensitive mammal: the American pika. PLoS One 7:e39077CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Hewitt G (2000) The genetic legacy of the Quaternary ice ages. Nature 405:907–913CAS 
    PubMed 
    Article 
    PubMed Central 

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

    Google Scholar 
    Hinzpeter A, Lipecka J, Brouillard F, Baudoin-Legros M, Dadlez M, Edelman A et al. (2006) Association between Hsp90 and the ClC-2 chloride channel upregulates channel function. Am J Physiol Cell Physiol 290:C45–56CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    Hoban S (2018) Integrative conservation genetics: Prioritizing populations using climate predictions, adaptive potential and habitat connectivity. Mol Ecol Resour 18:14–17PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    Hoffmann AA, Sgrò CM (2011) Climate change and evolutionary adaptation. Nature 470:479–485CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    Holbrook JD, DeYoung RW, Janecka JE, Tewes ME, Honeycutt RL, Young JH (2012) Genetic diversity, population structure, and movements of mountain lions (Puma concolor) in Texas. J Mammal 93:989–1000Article 

    Google Scholar 
    IPCC (2014) AR5 Synthesis Report: Climate Change 2014. IPCC.Jentsch TJ, Lutter D, Planells-Cases R, Ullrich F, Voss FK (2016) VRAC: molecular identification as LRRC8 heteromers with differential functions. Pflüg Arch Eur J Physiol 468:385–393CAS 
    Article 

    Google Scholar 
    Johnston AN, Bruggeman JE, Beers AT, Beever EA, Christophersen RG, Ransom JI (2019) Ecological consequences of anomalies in atmospheric moisture and snowpack. Ecology 100:e02638PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    Johnston KM, Freund KA, Schmitz OJ (2012) Projected range shifting by montane mammals under climate change: implications for Cascadia’s National Parks. Ecosphere 3:art97Article 

    Google Scholar 
    Kilham L (1958) Territorial behavior in pikas. J Mammal 39:307–307Article 

    Google Scholar 
    Kimura M (1971) Theoretical foundation of population genetics at the molecular level. Theor Popul Biol 2:174–208CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    Klingler KB, Jahner JP, Parchman TL, Ray C, Peacock MM (2021) Genomic variation in the American pika: signatures of geographic isolation and implications for conservation. BMC Ecol Evol 21:2CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Lambers JHR (2015) Extinction risks from climate change. Science 348:501–502CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    Latch EK, Dharmarajan G, Glaubitz JC, Rhodes OE (2006) Relative performance of Bayesian clustering software for inferring population substructure and individual assignment at low levels of population differentiation. Conserv Genet 7:295–302Article 

    Google Scholar 
    Latch EK, Scognamillo DG, Fike JA, Chamberlain MJ, Rhodes Jr OE (2008) Deciphering ecological barriers to North American River Otter (Lontra canadensis) gene flow in the Louisiana landscape. J Hered 99:265–274CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    Lee KM, Coop G (2017) Distinguishing among modes of convergent adaptation using population genomic data. Genetics 207:1591–1619PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Lee KM, Coop G (2019) Population genomics perspectives on convergent adaptation. Philos Trans R Soc B Biol Sci 374:20180236Article 

    Google Scholar 
    Lemay MA, Russello MA (2015) Genetic evidence for ecological divergence in kokanee salmon. Mol Ecol 24:798–811CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    Li H (2013) Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. ArXiv13033997 Q-Bio.MacArthur RA, Wang LCH (1974) Behavioral thermoregulation in the pika Ochotona princeps: a field study using radiotelemetry. Can J Zool 52:353–358CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    McCain CM (2019) Assessing the risks to United States and Canadian mammals caused by climate change using a trait-mediated model. J Mammal 100:1808–1817
    Google Scholar 
    Meirmans PG, Van Tienderen PH (2004) GENOTYPE and GENODIVE: two programs for the analysis of genetic diversity of asexual organisms. Mol Ecol Notes 4:792–794Article 

    Google Scholar 
    Morin PA, Luikart G, Wayne RK, the SNP Workshop Group (2004) SNPs in ecology, evolution and conservation. Trends Ecol Evol 19:208–216Article 

    Google Scholar 
    Moritz C, Agudo R (2013) The future of species under climate change: resilience or decline? Science 341:504–508CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    Moritz C, Patton JL, Conroy CJ, Parra JL, White GC, Beissinger SR (2008) Impact of a century of climate change on small-mammal communities in Yosemite National Park, USA. Science 322:261–264CAS 
    PubMed 
    Article 

    Google Scholar 
    Morrison SF, Hik DS (2007) Demographic analysis of a declining pika Ochotona collaris population: linking survival to broad-scale climate patterns via spring snowmelt patterns. J Anim Ecol 76:899–907PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    Moyer‐Horner L, Mathewson PD, Jones GM, Kearney MR, Porter WP (2015) Modeling behavioral thermoregulation in a climate change sentinel. Ecol Evol 5:5810–5822PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Mussmann SM, Douglas MR, Chafin TK, Douglas ME (2019) BA3-SNPs: contemporary migration reconfigured in BayesAss for next-generation sequence data. Methods Ecol Evol 10:1808–1813Article 

    Google Scholar 
    Parmesan C (2006) Ecological and evolutionary responses to recent climate change. Annu Rev Ecol Evol Syst 37:637–669Article 

    Google Scholar 
    Paz-Vinas I, Loot G, Hermoso V, Veyssière C, Poulet N, Grenouillet G et al. (2018) Systematic conservation planning for intraspecific genetic diversity. Proc R Soc B Biol Sci 285:20172746Article 

    Google Scholar 
    Peacock MM (1997) Determining natal dispersal patterns in a population of North American pikas (Ochotona princeps) using direct mark-resight and indirect genetic methods. Behav Ecol 8:340–350Article 

    Google Scholar 
    Peacock MM, Smith AT (1997) Nonrandom mating in pikas Ochotona princeps: evidence for inbreeding between individuals of intermediate relatedness. Mol Ecol 6:801–811CAS 
    PubMed 
    Article 

    Google Scholar 
    Pew J, Muir PH, Wang J, Frasier TR (2015) related: an R package for analysing pairwise relatedness from codominant molecular markers. Mol Ecol Resour 15:557–561PubMed 
    Article 

    Google Scholar 
    Pickett STA (1989) Space-for-time substitution as an alternative to long-term studies. In: Likens GE (ed) Long-term studies in ecology. Springer New York, New York, NY, p 110–135Chapter 

    Google Scholar 
    Piry S, Alapetite A, Cornuet J-M, Paetkau D, Baudouin L, Estoup A (2004) GENECLASS2: a software for genetic assignment and first-generation migrant detection. J Hered 95:536–539CAS 
    PubMed 
    Article 

    Google Scholar 
    R Core Team (2019) R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria, http://www.R-project.org/Rankin AM, Galbreath KE, Teeter KC (2017) Signatures of adaptive molecular evolution in American pikas (Ochotona princeps). J Mammal 98:1156–1167Article 

    Google Scholar 
    Rannala B, Mountain JL (1997) Detecting immigration by using multilocus genotypes. Proc Natl Acad Sci USA 94:9197–9201CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Razgour O, Taggart JB, Manel S, Juste J, Ibáñez C, Rebelo H et al. (2018) An integrated framework to identify wildlife populations under threat from climate change. Mol Ecol Resour 18:18–31PubMed 
    Article 

    Google Scholar 
    Rellstab C, Gugerli F, Eckert AJ, Hancock AM, Holderegger R (2015) A practical guide to environmental association analysis in landscape genomics. Mol Ecol 24:4348–4370PubMed 
    Article 

    Google Scholar 
    Ritland K (1996) Estimators for pairwise relatedness and individual inbreeding coefficients. Genet Res 67:175–185Article 

    Google Scholar 
    Robson KM, Lamb CT, Russello MA (2016) Low genetic diversity, restricted dispersal, and elevation-specific patterns of population decline in American pikas in an atypical environment. J Mammal 97:464–472Article 

    Google Scholar 
    Rochette NC, Rivera‐Colón AG, Catchen JM (2019) Stacks 2: analytical methods for paired-end sequencing improve RADseq-based population genomics. Mol Ecol 28:4737–4754CAS 
    PubMed 
    Article 

    Google Scholar 
    Rousset F (2008) GENEPOP’007: a complete re-implementation of the GENEPOP software for Windows and Linux. Mol Ecol Resour 8:103–106PubMed 
    Article 

    Google Scholar 
    Rubidge EM, Patton JL, Lim M, Burton AC, Brashares JS, Moritz C (2012) Climate-induced range contraction drives genetic erosion in an alpine mammal. Nat Clim Change 2:285–288Article 

    Google Scholar 
    Russello MA, Waterhouse MD, Etter PD, Johnson EA (2015) From promise to practice: pairing non-invasive sampling with genomics in conservation. PeerJ 3:e1106PubMed 
    PubMed Central 
    Article 
    CAS 

    Google Scholar 
    Savolainen O, Lascoux M, Merilä J (2013) Ecological genomics of local adaptation. Nat Rev Genet 14:807–820CAS 
    PubMed 
    Article 

    Google Scholar 
    Sjodin BMF, Galbreath KE, Lanier HC, Russello MA (2021) Chromosome-level reference genome assembly for the American pika (Ochotona princeps). J Hered https://doi.org/10.1093/jhered/esab031Smith AT (1974b) The distribution and dispersal of pikas: influences of behavior and climate. Ecology 55:1368–1376Article 

    Google Scholar 
    Smith AT (1974a) The distribution and dispersal of pikas: consequences of insular population structure. Ecology 55:1112–1119Article 

    Google Scholar 
    Smith AT (2020) Conservation status of American pikas (Ochotona princeps). J Mammal 101:1466–1488Article 

    Google Scholar 
    Smith AT, Ivins BL (1983) Colonization in a pika population: dispersal vs philopatry. Behav Ecol Sociobiol 13:37–47Article 

    Google Scholar 
    Smith AT, Weston ML (1990) Ochotona princeps. Mamm Species 352:1–8.Article 

    Google Scholar 
    Smith AT, Millar CI (2018) American pika (Ochotona princeps) population survival in winters with low or no snowpack. West North Am Nat 78:126–132Article 

    Google Scholar 
    La Sorte FA, Jetz W (2010) Projected range contractions of montane biodiversity under global warming. Proc R Soc B Biol Sci 277:3401–3410Article 

    Google Scholar 
    Stern DL (2013) The genetic causes of convergent evolution. Nat Rev Genet 14:751–764CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    Stewart JAE, Wright DH, Heckman KA (2017) Apparent climate-mediated loss and fragmentation of core habitat of the American pika in the Northern Sierra Nevada, California, USA. PLoS One 12:e0181834PubMed 
    PubMed Central 
    Article 
    CAS 

    Google Scholar 
    Stewart JAE, Perrine JD, Nichols LB, Thorne JH, Millar CI, Goehring KE et al. (2015) Revisiting the past to foretell the future: summer temperature and habitat area predict pika extirpations in California. J Biogeogr 42:880–890Article 

    Google Scholar 
    Varner J, Dearing MD (2014) The importance of biologically relevant microclimates in habitat suitability assessments. PLoS One 9:e104648PubMed 
    PubMed Central 
    Article 
    CAS 

    Google Scholar 
    Waldvogel A-M, Feldmeyer B, Rolshausen G, Exposito‐Alonso M, Rellstab C, Kofler R et al. (2020) Evolutionary genomics can improve prediction of species’ responses to climate change. Evol Lett 4:4–18PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Wang T, Hamann A, Spittlehouse DL, Murdock TQ (2012) ClimateWNA—high-resolution spatial climate data for Western North America. J Appl Meteorol Climatol 51:16–29Article 

    Google Scholar 
    Waterhouse MD, Erb LP, Beever EA, Russello MA (2018) Adaptive population divergence and directional gene flow across steep elevational gradients in a climate-sensitive mammal. Mol Ecol 27:2512–2528PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    Weir BS, Cockerham CC (1984) Estimating F-statistics for the analysis of population structure. Evolution 38:1358–1370CAS 
    PubMed 

    Google Scholar 
    Wiens JJ (2016) Climate-related local extinctions are already widespread among plant and animal species. PLoS Biol 14:e2001104PubMed 
    PubMed Central 
    Article 
    CAS 

    Google Scholar 
    Wilkening JL, Ray C (2016) Characterizing predictors of survival in the American pika (Ochotona princeps). J Mammal 97:1366–1375Article 

    Google Scholar 
    Wilkening JL, Ray C, Beever EA, Brussard PF (2011) Modeling contemporary range retraction in Great Basin pikas (Ochotona princeps) using data on microclimate and microhabitat. Quat Int 235:77–88Article 

    Google Scholar 
    Wilson GA, Rannala B (2003) Bayesian inference of recent migration rates using multilocus genotypes. Genetics 163:1177–1191PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Wogan GOU, Wang IJ (2018) The value of space-for-time substitution for studying fine-scale microevolutionary processes. Ecography 41:1456–1468Article 

    Google Scholar 
    Yandow LH, Chalfoun AD, Doak DF (2015) Climate tolerances and habitat requirements jointly shape the elevational distribution of the American pika (Ochotona princeps), with implications for climate change effects. PLoS One 10:e0131082PubMed 
    PubMed Central 
    Article 
    CAS 

    Google Scholar 
    Zgurski JM, Hik DS (2012) Polygynandry and even-sexed dispersal in a population of collared pikas, Ochotona collaris. Anim Behav 83:1075–1082Article 

    Google Scholar 
    Zhang D, Pan J, Cao J, Cao Y, Zhou H (2020) Screening of drought-resistance related genes and analysis of promising regulatory pathway in camel renal medulla. Genomics 112:2633–2639CAS 
    PubMed 
    Article 

    Google Scholar 
    Zhang Z, Schwartz S, Wagner L, Miller W (2000) A greedy algorithm for aligning DNA sequences. J Comput Biol 7:203–214CAS 
    PubMed 
    Article 

    Google Scholar  More

  • in

    Lethal microbial blooms delayed freshwater ecosystem recovery following the end-Permian extinction

    1.Paerl, H. W. & Otten, T. G. Harmful cyanobacterial blooms: causes, consequences, and controls. Microb. Ecol. 65, 995–1010 (2013).CAS 
    PubMed 
    Article 

    Google Scholar 
    2.Visser, P. M. et al. How rising CO2 and global warming may stimulate harmful cyanobacterial blooms. Harmful Algae 54, 145–159 (2016).CAS 
    PubMed 
    Article 

    Google Scholar 
    3.Lürling, M., Mendes e Mello, M., van Oosterhout, F., de Senerpont Domis, L. & Marinho, M. M. Response of natural cyanobacteria and algae assemblages to a nutrient pulse and elevated temperature. Front. Microbiol. 9, 1851 (2018).PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    4.Low-Décarie, E., Fussmann, G. F. & Bell, G. Aquatic primary production in a high-CO2 world. Trends Ecol. Evol. 29, 223–232 (2014).PubMed 
    Article 

    Google Scholar 
    5.Stanley, S. M. Estimates of the magnitudes of major marine mass extinctions in earth history. PNAS 113, E6325–E6334 (2016).ADS 
    CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    6.Sun, Y. D. et al. Lethally hot temperatures during the Early Triassic Greenhouse. Science 338, 366–370 (2012).ADS 
    CAS 
    PubMed 
    Article 

    Google Scholar 
    7.Frank, T. D. et al. Pace, magnitude, and nature of terrestrial climate change through the end Permian extinction in southeastern Gondwana. Geology 49, https://doi.org/10.1130/G48795.1 (2021).8.Wu, Y. et al. Six-fold increase of atmospheric pCO2 during the Permian–Triassic mass extinction. Nat. Commun. 12, 2137 (2021).ADS 
    CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    9.Burgess, S. D., Muirhead, J. D. & Bowring, S. A. Initial pulse of Siberian Traps sills as the trigger of the end-Permian mass extinction. Nat. Commun. 8, 164 (2017).ADS 
    CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    10.Mays, C. et al. Refined Permian–Triassic floristic timeline reveals early collapse and delayed recovery of south polar terrestrial ecosystems. GSA Bull. 132, 1489–1513 (2020).CAS 
    Article 

    Google Scholar 
    11.Chu, D. et al. Ecological disturbance in tropical peatlands prior to marine Permian-Triassic mass extinction. Geology 48, 288–292 (2020).ADS 
    Article 

    Google Scholar 
    12.Retallack, G. J., Veevers, J. J. & Morante, R. Global coal gap between Permian–Triassic extinction and Middle Triassic recovery of peat-forming plants. GSA Bull. 108, 195–207 (1996).CAS 
    Article 

    Google Scholar 
    13.Fielding, C. R. et al. Age and pattern of the southern high-latitude continental end-Permian extinction constrained by multiproxy analysis. Nat. Commun. 10, 385 (2019).ADS 
    CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    14.Fielding, C. R. et al. Sedimentology of the continental end-Permian extinction event in the Sydney Basin, eastern Australia. Sedimentology 68, 30–62 (2021).CAS 
    Article 

    Google Scholar 
    15.Metcalfe, I., Crowley, J. L., Nicoll, R. S. & Schmitz, M. High-precision U-Pb CA-TIMS calibration of Middle Permian to Lower Triassic sequences, mass extinction and extreme climate-change in eastern Australian Gondwana. Gondwana Res. 28, 61–81 (2015).ADS 
    CAS 
    Article 

    Google Scholar 
    16.Vajda, V. et al. End-Permian (252 Mya) deforestation, wildfires and flooding—an ancient biotic crisis with lessons for the present. Earth Planet. Sci. Lett. 529, 115875 (2020).CAS 
    Article 

    Google Scholar 
    17.McLoughlin, S. et al. Dwelling in the dead zone—vertebrate burrows immediately succeeding the end-Permian extinction event in Australia. Palaios 35, 342–357 (2020).ADS 
    Article 

    Google Scholar 
    18.Lamb, A. L., Wilson, G. P. & Leng, M. J. A review of coastal palaeoclimate and relative sea-level reconstructions using δ13C and C/N ratios in organic material. Earth-Sci. Rev. 75, 29–57 (2006).ADS 
    CAS 
    Article 

    Google Scholar 
    19.Mays, C., Vajda, V. & McLoughlin, S. Permian–Triassic non-marine algae of Gondwana—distributions, natural affinities and ecological implications. Earth-Sci. Rev. 212, 103382 (2021).CAS 
    Article 

    Google Scholar 
    20.McLoughlin, S. et al. Age and paleoenvironmental significance of the Frazer Beach Member—a new lithostratigraphic unit overlying the end-Permian extinction horizon in the Sydney Basin, Australia. Front. Earth Sci. 8, 600976 (2021).Article 

    Google Scholar 
    21.Huber, J. K. A postglacial pollen and nonsiliceous algae record from Gegoka Lake, Lake County, Minnesota. J. Paleolimnol. 16, 23–35 (1996).ADS 
    Article 

    Google Scholar 
    22.Woodward, C. A. & Shulmeister, J. A Holocene record of human induced and natural environmental change from Lake Forsyth (Te Wairewa), New Zealand. J. Paleolimnol. 34, 481–501 (2005).ADS 
    Article 

    Google Scholar 
    23.Pacton, M., Gorin, G. & Fiet, N. Occurrence of photosynthetic microbial mats in a Lower Cretaceous black shale (central Italy): a shallow-water deposit. Facies 55, 401–419 (2009).Article 

    Google Scholar 
    24.Pacton, M., Gorin, G. E. & Vasconcelos, C. Amorphous organic matter—Experimental data on formation and the role of microbes. Rev. Palaeobot. Palynol. 166, 253–267 (2011).Article 

    Google Scholar 
    25.Tyson, R. V. Sedimentary Organic Matter: Organic Facies and Palynofacies (Chapman & Hall, 1995).26.Retallack, G. J. Earliest Triassic claystone breccias and soil-erosion crisis. J. Sediment. Res. 75, 679–695 (2005).ADS 
    Article 

    Google Scholar 
    27.Augland, L. E. et al. The main pulse of the Siberian Traps expanded in size and composition. Sci. Rep. 9, 18723 (2019).ADS 
    CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    28.Retallack, G. J. Post-apocalyptic greenhouse paleoclimate revealed by earliest Triassic paleosols in the Sydney Basin. Aust. GSA Bull. 111, 52–70 (1999).CAS 
    Article 

    Google Scholar 
    29.Woodward, C., Shulmeister, J., Larsen, J., Jacobsen, G. E. & Zawadzki, A. The hydrological legacy of deforestation on global wetlands. Science 346, 844–847 (2014).ADS 
    CAS 
    PubMed 
    Article 

    Google Scholar 
    30.Stevenson, R. J. & Smol, J. P. In Freshwater Algae of North America: Ecology and Classification (eds Wehr, J. D., Sheath, R. G. & Kociolek, P.) Ch. 21, 921–962 (Academic Press, 2015).31.Lindström, S., Bjerager, M., Alsen, P., Sanei, H. & Bojesen-Koefoed, J. The Smithian–Spathian boundary in North Greenland: implications for extreme global climate changes. Geol. Mag. 157, 1547–1567 (2020).ADS 
    Article 
    CAS 

    Google Scholar 
    32.de Leeuw, J. W., Versteegh, G. J. M. & van Bergen, P. F. in Plants and Climate Change, Plant Ecology (eds Rozema, J., Aerts, R. & Cornelissen, H.) Vol. 182, 209–233 (Springer, 2006).33.Baudelet, P.-H., Ricochon, G., Linder, M. & Muniglia, L. A new insight into cell walls of Chlorophyta. Algal Res 25, 333–371 (2017).Article 

    Google Scholar 
    34.Graham, L. E. & Gray, J. In Plants Invade the Land: Evolutionary and Environmental Perspectives (eds Gensel, P. G. & Edwards, D.) 140–158 (Columbia University Press, 2001).35.Demura, M., Ioki, M., Kawachi, M., Nakajima, N. & Watanabe, M. M. Desiccation tolerance of Botryococcus braunii (Trebouxiophyceae, Chlorophyta) and extreme temperature tolerance of dehydrated cells. J. Appl. Phycol. 26, 49–53 (2014).CAS 
    PubMed 
    Article 

    Google Scholar 
    36.Del Cortona, A. et al. Neoproterozoic origin and multiple transitions to macroscopic growth in green seaweeds. PNAS 117, 2551–2559 (2020).PubMed 
    PubMed Central 
    Article 
    CAS 

    Google Scholar 
    37.Wheeler, A., Van de Wetering, N., Esterle, J. S. & Götz, A. E. Palaeoenvironmental changes recorded in the palynology and palynofacies of a Late Permian Marker Mudstone (Galilee Basin, Australia). Palaeoworld 29, 439–452 (2020).Article 

    Google Scholar 
    38.Reynolds, C. S., Huszar, V., Kruk, C., Naselli-Flores, L. & Melo, S. Towards a functional classification of the freshwater phytoplankton. J. Plankton Res. 24, 417–428 (2002).Article 

    Google Scholar 
    39.Low-Décarie, E., Fussmann, G. F. & Bell, G. The effect of elevated CO2 on growth and competition in experimental phytoplankton communities. Glob. Change Biol. 17, 2525–2535 (2011).ADS 
    Article 

    Google Scholar 
    40.von Alvensleben, N., Magnusson, M. & Heimann, K. Salinity tolerance of four freshwater microalgal species and the effects of salinity and nutrient limitation on biochemical profiles. J. Appl. Phycol. 28, 861–876 (2016).Article 
    CAS 

    Google Scholar 
    41.Chu, D. et al. Microbial mats in the terrestrial Lower Triassic of North China and implications for the Permian–Triassic mass extinction. Palaeogeog. Palaeoclimatol. Palaeoecol. 474, 214–231 (2017).ADS 
    Article 

    Google Scholar 
    42.Guo, W. et al. Secular variations of ichnofossils from the terrestrial Late Permian–Middle Triassic succession at the Shichuanhe section in Shaanxi Province, North China. Glob. Planet. Change 181, 102978 (2019).Article 

    Google Scholar 
    43.Lee, J. Y. et al. Future global climate: Scenario-based projections and near-term information. In Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change (eds Masson-Delmotte, V., et al.) 195 pp. (Cambridge University Press, 2021).44.de Jersey, N. J. Palynology of the Permian-Triassic transition in the western Bowen Basin. Geol. Surv. Qld. Publ. 374, 1–39 (1979).
    Google Scholar 
    45.Lindström, S. & McLoughlin, S. Synchronous palynofloristic extinction and recovery after the end-Permian event in the Prince Charles Mountains, Antarctica: Implications for palynofloristic turnover across Gondwana. Rev. Palaeobot. Palynol. 145, 89–122 (2007).Article 

    Google Scholar 
    46.Grebe, H. Permian plant microfossils from the Newcastle Coal Measures/Narrabeen Group Boundary, Lake Munmorah, New South Wales. Rec. Geol. Surv. NSW 12, 125–136 (1970).
    Google Scholar 
    47.Mishra, S. et al. A new acritarch spike of Leiosphaeridia dessicata comb. nov. emend. from the Upper Permian and Lower Triassic sequence of India (Pranhita-Godavari Basin): its origin and palaeoecological significance. Palaeogeog. Palaeoclimatol. Palaeoecol. 567, 110274 (2021).ADS 
    Article 

    Google Scholar 
    48.Grice, K. et al. Photic zone euxinia during the Permian-Triassic superanoxic event. Science 307, 706–709 (2005).ADS 
    CAS 
    PubMed 
    Article 

    Google Scholar 
    49.Kershaw, S. et al. Microbialites and global environmental change across the Permian–Triassic boundary: a synthesis. Geobiology 10, 25–47 (2012).CAS 
    PubMed 
    Article 

    Google Scholar 
    50.Schneebeli-Hermann, E. et al. Palynofacies analysis of the Permian–Triassic transition in the Amb section (Salt Range, Pakistan): implications for the anoxia on the South Tethyan Margin. J. Asian Earth Sci. 60, 225–234 (2012).ADS 
    Article 

    Google Scholar 
    51.van Soelen, E. E. & Kürschner, W. M. Late Permian to Early Triassic changes in acritarch assemblages and morphology in the Boreal Arctic: new data from the Finnmark Platform. Palaeogeog. Palaeoclimatol. Palaeoecol. 505, 120–127 (2018).ADS 
    Article 

    Google Scholar 
    52.Spina, A., Cirilli, S., Utting, J. & Jansonius, J. Palynology of the Permian and Triassic of the Tesero and Bulla sections (Western Dolomites, Italy) and consideration about the enigmatic species Reduviasporonites chalastus. Rev. Palaeobot. Palynol. 218, 3–14 (2015).Article 

    Google Scholar 
    53.Thomas, B. M. et al. Unique marine Permian‐Triassic boundary section from Western Australia. Aust. J. Earth Sci. 51, 423–430 (2004).ADS 
    Article 

    Google Scholar 
    54.Schneebeli-Hermann, E. & Bucher, H. Palynostratigraphy at the Permian-Triassic boundary of the Amb section, Salt Range, Pakistan. Palynology 39, 1–18 (2015).Article 

    Google Scholar 
    55.Lei, Y. et al. Phytoplankton (acritarch) community changes during the Permian-Triassic transition in South China. Palaeogeog. Palaeoclimatol. Palaeoecol. 519, 84–94 (2019).ADS 
    Article 

    Google Scholar 
    56.Algeo, T. J. et al. Plankton and productivity during the Permian–Triassic boundary crisis: An analysis of organic carbon fluxes. Glob. Planet. Change 105, 52–67 (2013).ADS 
    Article 

    Google Scholar 
    57.van Soelen, E. E., Twitchett, R. J. & Kürschner, W. M. Salinity changes and anoxia resulting from enhanced run-off during the late Permian global warming and mass extinction event. Climate 14, 441–453 (2018).
    Google Scholar 
    58.Kaiho, K. et al. Effects of soil erosion and anoxic–euxinic ocean in the Permian–Triassic marine crisis. Heliyon 2, e00137 (2016).PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    59.Bond, D. P. G. & Grasby, S. E. On the causes of mass extinctions. Palaeogeog. Palaeoclimatol. Palaeoecol. 478, 3–29 (2017).ADS 
    Article 

    Google Scholar 
    60.Lindström, S., Erlström, M., Piasecki, S., Nielsen, L. H. & Mathiesen, A. Palynology and terrestrial ecosystem change of the Middle Triassic to lowermost Jurassic succession of the eastern Danish Basin. Rev. Palaeobot. Palynol. 244, 65–95 (2017).Article 

    Google Scholar 
    61.Garel, S. et al. Paleohydrological and paleoenvironmental changes recorded in terrestrial sediments of the Paleocene–Eocene boundary (Normandy, France). Palaeogeog. Palaeoclimatol. Palaeoecol. 376, 184–199 (2013).ADS 
    Article 

    Google Scholar 
    62.van de Schootbrugge, B. & Gollner, S. In Ecosystem Paleobiology and Geobiology, The Paleontological Society Papers (eds Bush, A. M., Pruss, S. B. & Payne, J. L.) 19, 87–114 (The Paleontological Society, 2013).63.Mata, S. A. & Bottjer, D. J. Microbes and mass extinctions: paleoenvironmental distribution of microbialites during times of biotic crisis. Geobiology 10, 3–24 (2012).CAS 
    PubMed 
    Article 

    Google Scholar 
    64.Peterffy, O., Calner, M. & Vajda, V. Early Jurassic microbial mats—a potential response to reduced biotic activity in the aftermath of the end-Triassic mass extinction event. Palaeogeog. Palaeoclimatol. Palaeoecol. 464, 76–85 (2016).ADS 
    Article 

    Google Scholar 
    65.Schoene, B. et al. U-Pb constraints on pulsed eruption of the Deccan Traps across the end-Cretaceous mass extinction. Science 3636, 862–866 (2019).ADS 
    Article 
    CAS 

    Google Scholar 
    66.Hull, P. M. et al. On impact and volcanism across the Cretaceous-Paleogene boundary. Science 367, 266–272 (2020).ADS 
    CAS 
    PubMed 
    Article 

    Google Scholar 
    67.Vajda, V., Ocampo, A., Ferrow, E. & Bender Koch, C. Nano particles as the primary cause for long-term sunlight suppression at high southern latitudes following the Chicxulub impact—evidence from ejecta deposits in Belize and Mexico. Gondwana Res. 27, 1079–1088 (2015).ADS 
    CAS 
    Article 

    Google Scholar 
    68.Sepúlveda, J., Wendler, J. E., Summons, R. E. & Hinrichs, K.-U. Rapid resurgence of marine productivity after the Cretaceous-Paleogene mass extinction. Science 326, 129–132 (2009).ADS 
    PubMed 
    Article 
    CAS 

    Google Scholar 
    69.Bralower, T. J. et al. Origin of a global carbonate layer deposited in the aftermath of the Cretaceous-Paleogene boundary impact. Earth Planet. Sci. Lett. 548, 116476 (2020).CAS 
    Article 

    Google Scholar 
    70.Schaefer, B. et al. Microbial life in the nascent Chicxulub crater. Geology 48, 328–332 (2020).ADS 
    CAS 
    Article 

    Google Scholar 
    71.Milligan, J. N., Royer, D. L., Franks, P. J., Upchurch, G. R. & McKee, M. L. No evidence for a large atmospheric CO2 spike across the Cretaceous‐Paleogene boundary. Geophys. Res. Lett. 46, 3462–3472 (2019).ADS 
    CAS 
    Article 

    Google Scholar 
    72.Strother, P. K. & Wellman, C. H. The Nonesuch Formation Lagerstätte: a rare window into freshwater life one billion years ago. J. Geol. Soc. 178, jgs2020–jgs2133 (2021).Article 

    Google Scholar 
    73.Sepkoski, J. J., Bambach, R. K. & Droser, M. L. In Cycles and Events in Stratigraphy (eds Einsele, G., Ricken, W. & Seilacher, A.) 298–312 (Springer-Verlag, 1991).74.Tyson, R. V. Calibration of hydrogen indices with microscopy: a review, reanalysis and new results using the fluorescence scale. Org. Geochem. 37, 45–63 (2006).CAS 
    Article 

    Google Scholar 
    75.Benninghoff, W. S. Calculation of pollen and spore density in sediments by addition of exotic pollen in known quantities. Pollen et. Spores 4, 332–333 (1962).
    Google Scholar 
    76.Maher, L. J. Statistics for microfossil concentration measurements employing samples spiked with marker grains. Rev. Palaeobot. Palynol. 32, 153–191 (1981).Article 

    Google Scholar 
    77.Simpson, M. G. Plant Systematics (Academic Press, 2019).78.Evitt, W. R. A discussion and proposals concerning fossil dinoflagellates, hystrichospheres, and acritarchs, II. PNAS 49, 298–302 (1963).ADS 
    CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    79.Rampino, M. R. & Eshet, Y. The fungal and acritarch events as time markers for the latest Permian mass extinction: an update. Geosci. Front. 9, 147–154 (2018).Article 

    Google Scholar 
    80.Combaz, A. Les palynofaciès. Rev. Micropaléontol. 7, 205–218 (1964).
    Google Scholar 
    81.Hammer, Ø., Harper, D. A. T. & Ryan, P. D. PAST: paleontological statistics software package for education and data analysis. Palaeontol. Electron. 4, 4 (2001).
    Google Scholar 
    82.Wei, W. & Algeo, T. J. Elemental proxies for paleosalinity analysis of ancient shales and mudrocks. Geochim. Cosmochim. Acta 287, 341–366 (2020).ADS 
    CAS 
    Article 

    Google Scholar 
    83.Rowe, H., Hughes, N. & Robinson, K. The quantification and application of handheld energy-dispersive x-ray fluorescence (ED-XRF) in mudrock chemostratigraphy and geochemistry. Chem. Geol. 324–325, 122–131 (2012).ADS 
    Article 
    CAS 

    Google Scholar 
    84.Blakey, R. C. Global paleogeography and tectonics in deep time. https://deeptimemaps.com/global-series-details/. Accessed 16 June 2020 (2016).85.Zhuravlev, A. Y. & Wood, R. A. Anoxia as the cause of the mid-Early Cambrian (Botomian) extinction event. Geology 24, 311–314 (1996).ADS 
    CAS 
    Article 

    Google Scholar 
    86.Zhang, W., Shi, X., Jiang, G., Tang, D. & Wang, X. Mass-occurrence of oncoids at the Cambrian Series 2–Series 3 transition: Implications for microbial resurgence following an Early Cambrian extinction. Gondwana Res. 28, 432–450 (2015).ADS 
    Article 
    CAS 

    Google Scholar 
    87.Vecoli, M. Fossil microphytoplankton dynamics across the Ordovician–Silurian boundary. Rev. Palaeobot. Palynol. 148, 91–107 (2008).Article 

    Google Scholar 
    88.Xie, S. et al. Contrasting microbial community changes during mass extinctions at the Middle/Late Permian and Permian/Triassic boundaries. Earth Planet. Sci. Lett. 460, 180–191 (2017).ADS 
    CAS 
    Article 

    Google Scholar 
    89.Eshet, Y., Rampino, M. R. & Visscher, H. Fungal event and palynological record of ecological crisis and recovery across the Permian-Triassic boundary. Geology 23, 967–970 (1995).ADS 
    Article 

    Google Scholar 
    90.Richoz, S. et al. Hydrogen sulphide poisoning of shallow seas following the end-Triassic extinction. Nat. Geosci. 5, 662–667 (2012).ADS 
    CAS 
    Article 

    Google Scholar 
    91.Lindström, S. et al. No causal link between terrestrial ecosystem change and methane release during the end-Triassic mass extinction. Geology 40, 531–534 (2012).ADS 
    Article 
    CAS 

    Google Scholar 
    92.van de Schootbrugge, B. et al. End-Triassic calcification crisis and blooms of organic-walled “disaster species”. Palaeogeog. Palaeoclimatol. Palaeoecol. 244, 126–141 (2007).ADS 
    Article 

    Google Scholar 
    93.Slater, S. M., Twitchett, R. J., Danise, S. & Vajda, V. Substantial vegetation response to Early Jurassic global warming with impacts on oceanic anoxia. Nat. Geosci. 12, 462–467 (2019).ADS 
    CAS 
    Article 

    Google Scholar 
    94.Polgári, M. et al. Mineral and chemostratigraphy of a Toarcian black shale hosting Mn-carbonate microbialites (Úrkút, Hungary). Palaeogeog. Palaeoclimatol. Palaeoecol. 459, 99–120 (2016).ADS 
    Article 

    Google Scholar 
    95.Xu, W. et al. Carbon sequestration in an expanded lake system during the Toarcian oceanic anoxic event. Nat. Geosci. 10, 129–134 (2017).ADS 
    CAS 
    Article 

    Google Scholar 
    96.Kashiyama, Y. et al. Diazotrophic cyanobacteria as the major photoautotrophs during mid-Cretaceous oceanic anoxic events: nitrogen and carbon isotopic evidence from sedimentary porphyrin. Org. Geochem. 39, 532–549 (2008).CAS 
    Article 

    Google Scholar 
    97.Jarvis, I. et al. Microfossil assemblages and the Cenomanian-Turonian (late Cretaceous) oceanic anoxic event. Cretac. Res 9, 3–103 (1988).Article 

    Google Scholar 
    98.Layeb, M., Ben Fadhel, M., Layeb-Tounsi, Y. & Ben Youssef, M. First microbialites associated to organic-rich facies of the Oceanic Anoxic Event 2 (Northern Tunisia, Cenomanian–Turonian transition). Arab. J. Geosci. 7, 3349–3363 (2014).CAS 
    Article 

    Google Scholar 
    99.Pearce, M. A., Jarvis, I. & Tocher, B. A. The Cenomanian–Turonian boundary event, OAE2 and palaeoenvironmental change in epicontinental seas: new insights from the dinocyst and geochemical records. Palaeogeog. Palaeoclimatol. Palaeoecol. 280, 207–234 (2009).ADS 
    Article 

    Google Scholar 
    100.Kuypers, M. M. M., Pancost, R. D., Nijenhuis, I. A. & Sinninghe Damsté, J. S. Enhanced productivity led to increased organic carbon burial in the euxinic North Atlantic basin during the late Cenomanian oceanic anoxic event. Paleoceanography 17, 1051 (2002).ADS 
    Article 

    Google Scholar 
    101.Dodsworth, P., Eldrett, J. S. & Hart, M. B. Cretaceous Oceanic Anoxic Event 2 in eastern England: further palynological and geochemical data from Melton Ross. figshare https://doi.org/10.6084/m9.figshare.c.4987205.v3 (2020).102.Schwab, K. W., Bayliss, G. S., Smith, M. A. & Yoder, N. B. Mushroom and broccoli-head shaped algal fragments from the Eagle Ford Shale of south Texas and Coahuila, Mexico. Search and Discovery 70134 (2013).103.Lyson, T. R. et al. Exceptional continental record of biotic recovery after the Cretaceous–Paleogene mass extinction. Science 366, 977–983 (2019).ADS 
    CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    104.Scasso, R. A. et al. A high-resolution record of environmental changes from a Cretaceous-Paleogene section of Seymour Island. Antarctica. Palaeogeog. Palaeoclimatol. Palaeoecol. 555, 109844 (2020).ADS 
    Article 

    Google Scholar 
    105.Sosa-Montes de Oca, C. et al. Minor changes in biomarker assemblages in the aftermath of the Cretaceous-Paleogene mass extinction event at the Agost distal section (Spain). Palaeogeog. Palaeoclimatol. Palaeoecol. 569, 110310 (2021).ADS 
    Article 

    Google Scholar 
    106.Sluijs, A. et al. Environmental precursors to rapid light carbon injection at the Palaeocene/Eocene boundary. Nature 450, 1218–1222 (2007).ADS 
    CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    107.Junium, C. K., Dickson, A. J. & Uveges, B. T. Perturbation to the nitrogen cycle during rapid Early Eocene global warming. Nat. Commun. 9, 3186 (2018).ADS 
    PubMed 
    PubMed Central 
    Article 
    CAS 

    Google Scholar 
    108.Pagani, M. et al. Arctic hydrology during global warming at the Palaeocene/Eocene thermal maximum. Nature 442, 671–675 (2006).ADS 
    CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    109.Kender, S. et al. Marine and terrestrial environmental changes in NW Europe preceding carbon release at the Paleocene–Eocene transition. Earth Planet. Sci. Lett. 353–354, 108–120 (2012).ADS 
    Article 
    CAS 

    Google Scholar 
    110.Huisman, J. et al. Cyanobacterial blooms. Nat. Rev. Microbiol. 16, 471–483 (2018).CAS 
    PubMed 
    Article 

    Google Scholar 
    111.Paerl, H. W., Hall, N. S. & Calandrino, E. S. Controlling harmful cyanobacterial blooms in a world experiencing anthropogenic and climatic-induced change. Sci. Tot. Environ. 409, 1739–1745 (2011).CAS 
    Article 

    Google Scholar  More

  • in

    A symbiotic aphid selfishly manipulates attending ants via dopamine in honeydew

    1.Darwin, C. On the origin of species. (D. Appleton and Co., 1871). https://doi.org/10.5962/bhl.title.28875.2.Thompson, J. N. Mutualistic webs of species. Science (80–) 312, 372–373 (2006).CAS 
    Article 

    Google Scholar 
    3.Bronstein, J. L. The exploitation of mutualisms. Ecol. Lett. 4, 277–287 (2001).Article 

    Google Scholar 
    4.Bshary, R. & Grutter, A. S. Experimental evidence that partner choice is a driving force in the payoff distribution among cooperators or mutualists: The cleaner fish case. Ecol. Lett. 5, 130–136 (2002).Article 

    Google Scholar 
    5.Kiers, E. T., Rousseau, R. A., West, S. A. & Denison, R. F. Host sanctions and the legume–rhizobium mutualism. Nature 425, 78–81 (2003).ADS 
    CAS 
    Article 

    Google Scholar 
    6.Heil, M., Barajas-Barron, A., Orona-Tamayo, D., Wielsch, N. & Svatos, A. Partner manipulation stabilises a horizontally transmitted mutualism. Ecol. Lett. 17, 185–192 (2014).Article 

    Google Scholar 
    7.Hindsbo, O. Effects of Polymorphus (Acanthocephala) on colour and behaviour of Gammarus lacustris. Nature 238, 333 (1972).ADS 
    Article 

    Google Scholar 
    8.Thomas, F., Renaud, F., de Meeus, T. & Poulin, R. Manipulation of host behaviour by parasites: Ecosystem engineering in the intertidal zone?. Proc. R. Soc. B Biol. Sci. 265, 1091–1096 (1998).Article 

    Google Scholar 
    9.Thomas, F. et al. Do hairworms (Nematomorpha) manipulate the water seeking behaviour of their terrestrial hosts?. J. Evol. Biol. 15, 356–361 (2002).Article 

    Google Scholar 
    10.Kadoya, E. Z., Ishii, H. S. & Williams, N. M. Host manipulation of bumble bee queens by Sphaerularia nematodes indirectly affects foraging of non-host workers. Ecology 96, 1361–1370 (2015).Article 

    Google Scholar 
    11.Hojo, M. K., Pierce, N. E. & Tsuji, K. Lycaenid caterpillar secretions manipulate attendant ant behavior. Curr. Biol. 25, 2260–2264 (2015).CAS 
    Article 

    Google Scholar 
    12.Poulin, R., Brodeur, J. & Moore, J. Parasite manipulation of host behaviour: Should hosts always lose?. Oikos 70, 479 (1994).Article 

    Google Scholar 
    13.Heil, M. et al. Divergent investment strategies of Acacia myrmecophytes and the coexistence of mutualists and exploiters. Proc. Natl. Acad. Sci. USA. 106, 18091–18096 (2009).ADS 
    CAS 
    Article 

    Google Scholar 
    14.Watanabe, S., Murakami, T., Yoshimura, J. & Hasegawa, E. Color polymorphism in an aphid is maintained by attending ants. Sci. Adv. 2, (2016).15.Watanabe, S., Yoshimura, J. & Hasegawa, E. Ants improve the reproduction of inferior morphs to maintain a polymorphism in symbiont aphids. Sci. Rep. 8, (2018).16.Sakata, H. Density-dependent predation of the ant Lasius niger (Hymenoptera: Formicidae) on two attended aphids Lachnus tropicalis and Myzocallis kuricola (Homoptera: Aphididae). Res. Popul. Ecol. (Kyoto) 37, 159–164 (1995).Article 

    Google Scholar 
    17.Evans, P. D. Biogenic Amines in the insect nervous system. Adv. In Insect Phys. 15, 317–473 (1980).CAS 
    Article 

    Google Scholar 
    18.Aonuma, H. & Watanabe, T. Octopaminergic system in the brain controls aggressive motivation in the ant Formica japonica. Acta Biol. Hung. 63, 63–68 (2012).Article 

    Google Scholar 
    19.Stevenson, P. A., Dyakonova, V., Rillich, J. & Schildberger, K. Octopamine and experience-dependent modulation of aggression in crickets. J. Neurosci. 25, 1431–1441 (2005).CAS 
    Article 

    Google Scholar 
    20.Kostowski, W. & Tarchalska, B. The effects of some drugs affecting brain 5-HT on the aggressive behaviour and spontaneous electrical activity of the central nervous system of the ant Formica rufa. Brain Res. 38, 143–149 (1972).CAS 
    Article 

    Google Scholar 
    21.Szczuka, A. et al. The effects of serotonin, dopamine, octopamine and tyramine on behavior of workers of the ant Formica polyctena during dyadic aggression tests. Acta Neurobiol. Exp. (Wars) 73, 495–520 (2013).
    Google Scholar 
    22.Way, M. J. Mutualism between ants and honeydew-producing homoptera. Annu. Rev. Entomol. 8, 307–344 (1963).Article 

    Google Scholar 
    23.Hafer-Hahmann, N. Behavior out of control: Experimental evolution of resistance to host manipulation. Ecol. Evol. 9, 7237–7245 (2019).Article 

    Google Scholar 
    24.Martinez, J., Fleury, F. & Varaldi, J. Heritable variation in an extended phenotype: The case of a parasitoid manipulated by a virus. J. Evol. Biol. 25, 54–65 (2012).Article 

    Google Scholar 
    25.Engelstädter, J. & Hurst, G. D. D. The ecology and evolution of microbes that manipulate host reproduction. Annu. Rev. Ecol. Evol. Syst. 40, 127–149 (2009).Article 

    Google Scholar 
    26.Rosenthal, G. G. & Servedio, M. R. Chase-away sexual selection: Resistance to ‘resistance’. Evolution (N.Y.) 53, 296 (1999).
    Google Scholar 
    27.Woodring, J., Wiedemann, R., Fischer, M. K., Hoffmann, K. H. & Völkl, W. Honeydew amino acids in relation to sugars and their role in the establishment of ant-attendance hierarchy in eight species of aphids feeding on tansy (Tanacetum vulgare). Physiol. Entomol. 29, 311–319 (2004).CAS 
    Article 

    Google Scholar 
    28.Stadler, B. & Dixon, A. F. G. Ecology and evolution of aphid-ant interactions. Annu. Rev. Ecol. Evol. Syst. 36, 345–372 (2005).Article 

    Google Scholar 
    29.Tsuji, K. & Dobata, S. Social cancer and the biology of the clonal ant Pristomyrmex punctatus (Hymenoptera: Formicidae). Myrmecological News 15, 91–99 (2011).
    Google Scholar 
    30.Vellend, M. Conceptual synthesis in community ecology. Q. Rev. Biol. 85, 183–206 (2010).Article 

    Google Scholar 
    31.Agawa, H. & Kawata, M. The effect of color polymorphism on mortality in the aphid Macrosiphoniella yomogicola. Ecol. Res. 10, 301–306 (1995).Article 

    Google Scholar 
    32.Watanabe, S., Murakami, Y. & Hasegawa, E. Effects of attending ant species on the fate of colonies of an aphid, Macrosiphoniella yomogicola (Matsumura) (Homoptera: Aphididae), in an ant-aphid symbiosis. Entomol. News 128, 325 (2019).Article 

    Google Scholar 
    33.Wada-Katsumata, A., Yamaoka, R. & Aonuma, H. Social interactions influence dopamine and octopamine homeostasis in the brain of the ant Formica japonica. J. Exp. Biol. 214, 1707–1713 (2011).CAS 
    Article 

    Google Scholar 
    34.Aonuma, H. & Watanabe, T. Changes in the content of brain biogenic amine associated with early colony establishment in the queen of the ant, formica japonica. PLoS One 7, (2012).35.Aonuma, H. Serotonergic control in initiating defensive responses to unexpected tactile stimuli in the trap-jaw ant Odontomachus kuroiwae. J. Exp. Biol. 223, jeb228874 (2020).Article 

    Google Scholar 
    36.R Core Team. R: A Language and Environment for Statistical Computing. (2020).37.Bates, D., Mächler, M., Bolker, B. & Walker, S. Fitting linear mixed-effects models using {lme4}. J. Stat. Softw. 67, 1–48 (2015).Article 

    Google Scholar 
    38.Venables, W. N. & Ripley, B. D. Modern Applied Statistics with S. (Springer, 2002).39.Hothorn, T. & Hornik, K. exactRankTests: Exact Distributions for Rank and Permutation Tests. (2019). More

  • in

    First tracks of newborn straight-tusked elephants (Palaeoloxodon antiquus)

    The MTS proboscidean tracks and trackmakersRounded-to-elliptical tracks, with an axial length range from 9.6 to 54.5 cm (pes), were found mostly isolated and as manus-pes couples, or associated forming at least eight short trackways (see Table 1). They reveal good preservation in one 6-footprint trackway (see below), two converging trackways and some couples, showing anteriorly directed, wide, short and blunt toe impressions (Figs. 2, 3 and 4). Toe impressions are not commonly visible in elephant footprints9,13, (but see27), which attests to cases of exceptional preservation in Matalascañas tracks. Preservation as true tracks is identified through expulsion marginal rims (e.g., Fig. 4a, g) and possible ejecta (Fig. 3b,e). Large and flat sole surfaces sometimes show evidence of pockmarks23 (Fig. 4f).Table 1 Measurements of Proboscipeda tracks, ordered from smallest to largest in length.Full size tableFigure 2Proboscidean tracks (Proboscipeda panfamilia) attributed in the MTS to straight-tusked elephants. (a–h) Morphological features of small-sized tracks produced by calves and juveniles. Examples of manus impressions in (a) PAT/MTS/011a, (b) PAT/MTS/016 and (f) PAT/MTS/015x, and for further interpretation of (a) see Fig. 3; the latter two with drag marks made during the foot-off event. (c) and (g) PAT/MTS/002a,b: Manus-pes couple found isolated showing heteropody and different number of toe impressions (interpretation as left-side tracks by peak pressure deformation in the left side of the track according to27); interpretation in (c). (d) PAT/MTS/014 and (e) PAT/MTS/007a: Calf-sized pes with three toe impressions. (h) PAT/MTS/011 h: Badly preserved manus of a calf. Scale bar = 5 cm.Full size imageFigure 3Photograph, outline, high-resolution 3D and false-coloured 3D images of the PAT/MTS/0011a track representing the best preserved manus of a juvenile-sized Proboscipeda track. (a) and (c) From the photograph and high-resolution images, five toe impressions in the anterior part of the rounded track are clear (especially toes I–IV). (b) and (f) The false coloured images in orthogonal (b) and oblique angle views (f) highlight the deepening of the track fore- and outwards, thus revealing a peak pressure pattern typical of left forefoot (toes III–IV), as well as a possible ejecta mound in front of the track. The poorly evident and narrow expulsion rim developed around the track is the result of the high cohesiveness and plasticity of the clayey fine-sand substrate. (d) Contour map supporting previous interpretation. (e) The cross-section of the track details the anterior migration of the foot pressure during its rotation, creating a peak pressure in the foot-off event that is represented in the deepest part of the track. Scale bars are 10 cm.Full size imageFigure 4Large-sized Proboscipeda tracks attributed to P. antiquus adults. (a) to (d) PAT/MTS/001: Right manus showing clearly 5 toe impressions and the frontal and lateral displacement rims (morphological interpretation based on the orthogonal (b) and oblique (d) depth and contour (c) maps). (e) and (f) PAT/MTS/010e: Deeper manus with pockmarks; toe pad impressions indicated (I–III). (g) PAT/MTS/004a,b: large manus-pes couple where the hind foot deformed the fore foot during overstepping, and revealing a typical elephantine gait; the toe impressions in both tracks indicate the direction of movement. Scale bar = 10 cm.Full size imageIrrespective of the track size, pes are elliptical to sub-rounded, with the length axis larger than the width and manus are circular or elliptical, with the width axis larger than the length (Figs. 2c and 4d, g for small and large size tracks, respectively). The safest way to differentiate between pes and manus is through the orientation of the track provided by the toe impressions, or by the orientation of the longer axis in trackways. When arranged in trackways, manus-pes couples show the typical elephantine gait, showing a short pace resulting from the fore- and hind feet on the same side swinging forward simultaneously below the body, as it is known from modern elephant gait28. In some cases, the partial impression of a pes overstepping the proximal part of a manus can be seen (Fig. 2c, g). Based on similar preservational style and opposing directions of movement without overlapping at the meeting point, a converging pair of trackways was apparently produced contemporaneously by an adult and a rather small juvenile. Sharp edges of the toe impressions indicate the presence of nails. These are found mostly in well preserved, smaller-sized tracks (Fig. 2a, d, e) because nails are commonly worn down in adult elephants and not always shown in their tracks13. These morphological features allow us to attribute the MTS trackways to the ichnospecies Proboscipeda panfamilia used previously for describing, among other tracksites, those tracks attributed confidently to the straight-tusked elephant Palaeoloxodon antiquus in the paleogeographical context of southern Europe11,14 (see supplementary Table S1).Manus-pes couples, when showing overstepping, were not considered in Table 1 (Fig. 2c, g). Overstepping depends on the speed of walking; at faster speeds the overstepping is only partial or there is no overstepping; elephants maintain the footfall pattern at all speeds, shifting toward a calculated 25% phase offset between limbs as they increase speed28 (Fig. 2g). The smallest tracks usually do not show overstepping possibly because of the greater activity, with longer pace and stride lengths, demonstrated by calves and juveniles when compared to adults. Manus or pes showing a large width-length ratio (below 0.80–0.96 sensu25) were not considered for the estimates since they represent slippage.Younger elephants have more pliable skin and musculature than adults. Also, the greater expansion and distribution of the weight in heavier adult animals is enough to reduce or negate toe impressions in some types of sediments, such as compacted substrates24,29. Interpreting the sedimentological data for the paleosol where MTS was developed15,17,30, suggests a drying clayey-sandy substrate14 that was still plastic enough to absorb the impact of the limbs during the locomotion of the elephants (presence of expulsion rims and absence of radial pressure cracks), and preserving, in many cases, the morphological details of the feet in good condition (Figs. 2a, 3, 4a; see Fig. 2h for a badly preserved example).Ichnological inference about the height, body mass and age of Palaeoloxodon antiquus in the MTSSeveral methods have been proposed for estimating the height at the shoulders for proboscideans, and the relationship between body mass and age with shoulder height 1,31,32. A linear relationship between foot length and shoulder height was confirmed by Lee and Moss33 from extant elephants and compared with fossil examples by Pasenko24. Pes length has been especially used in studies as an indicator of shoulder height21,34,35,36. Among Asian elephants, manus circumference has been shown to have a similar predictive relationship with shoulder height33. These parameters were determined for each isolated track (or representative track in a trackway), including manus and pes (Table 1), using equations previously proposed31,33 (see Methods). A similar approach has been applied to mammoth track studies in North America21,27, where modern ontogenetic and body-mass data has been used to provide age and size estimates from fossil tracks.From the skeletal record, sexual dimorphism of P. antiquus was observed to be more accentuated than in extant elephants, especially in terms of size differences1. During the first 10 years of life, both male and female African bush elephant foot lengths increase rapidly, with the fastest growth shown in the first two years for calves33,37. In P. antiquus, males would have continued to grow until their fifties according to bone data1, while females would have been much smaller as result of energy expenditure with reproduction, flattening the growth curve just after puberty. That is why the equations of Lee and Moss33 that discriminates the shoulder height from tracks for males and females have been applied. However, by comparison with the study of Marano and Palombo32 (based on the progress of eruption and degree of wear of teeth compared to extant elephants), and the body mass correlation of Larramendi et al.1 for calculating the age of P. antiquus, our MTS ages obtained from the application of the regression curve of Lee and Moss33 are underestimated and must be analysed as minimum age approximations for track lengths corresponding to adolescent and adult animals, especially for males. The obtained estimations from tracks are subject to a level of uncertainty related to biotic and abiotic factors that can distort the data (i.e., taphonomy) as it happens also with the calculations taken from skeletal proportions. Therefore, McNeil et al.21 even included data from frozen mammoth carcasses on the growth curve of Lee and Moss33 for correcting size discrepancies along ontogeny. For P. antiquus, our best data for comparison comes, however, from the flesh reconstructions1.Ontogenetic implicationsBased on the best fossil site found for this species in Europe, corresponding to 70 individual Palaeoloxodon antiquus specimens recovered in Geiseltal, Germany, Larramendi et al.1 developed the best reconstruction, so far, of the life appearance of this species and discussed size, body mass, ontogeny and sexual dimorphism. The Neumark-Nord bone site may be contemporary or slightly older than MTS, corresponding to late Middle Pleistocene-to-Eemian interglacial period1. The authors found that the body mass of P. antiquus males was up to three times more that of male Asian elephants and twice that of extant male African bush elephants. The large size determined for straight-tusked elephants (with an estimated  > 400 cm shoulder height in the flesh and body mass of 13 tonnes) and a later complete epiphyseal-diaphyseal fusion of limb bones (not yet totally fused at an estimated age of 47 years), in comparison with extant elephants, suggests that this species had a longer lifespan of 80 years or more1. Sexual dimorphism of P. antiquus was observed to be more accentuated than in extant elephants, with females generally not exceeding 300 cm at the shoulders with an estimated weight of not more than 5.5 tonnes, while males continued to grow until their fifties1. Males in extant elephant species grow more rapidly than females after puberty (i.e., around 7 years in age), which are affected by a trade-off between growth and reproduction. Under normal nutritional conditions, the growth rate is generally higher in males than females leading to a marked difference in size between sexes at already around 10 years in age33,37,38,39.The ontogenetic variation in growth projected for the MTS, when compared to what we known from extant proboscideans, is expressed in the track size distribution plot, with the definition of five age classes (Fig. 5; see also Table 1): calves under 2 years in age (when extant elephants experience fastest growth rates in both sexes), juveniles between 2 and 7 years in age (up to when elephant females reach their sexual maturity and therefore experience a strong reduction of growth rate in comparison to males), 7–15 years in age which include pre-puberty males and young female adults, over 15 years in age and  More