Fossil evidence reveals how plants responded to cooling during the Cretaceous-Paleogene transition
BMC Plant Biology volume 19, Article number: 402 (2019)
Around the Cretaceous-Paleogene (K-Pg) boundary, an obvious global cooling occurred, which resulted in dramatic changes in terrestrial ecosystems and the evolutionary trends of numerous organisms. However, how plant lineages responded to the cooling has remained unknown until now. Between ca. 70–60 Ma Mesocyparis McIver & Basinger (Cupressaceae), an extinct conifer genus, was distributed from eastern Asia to western North America and provides an excellent opportunity to solve this riddle.
Here we report a new species, Mesocyparis sinica from the early Paleocene of Jiayin, Heilongjiang, northeastern China. By integrating lines of evidence from phylogeny and comparative morphology of Mesocyparis, we found that during ca.70–60 Ma, the size of seed cone of Mesocyparis more than doubled, probably driven by the cooling during the K-Pg transition, which might be an effective adaptation for seed dispersal by animals. More importantly, we discovered that the northern limit of this genus, as well as those of two other arboreal taxa Metasequoia Miki ex Hu et Cheng (gymnosperm) and Nordenskioldia Heer (angiosperm), migrated ca.4–5° southward in paleolatitude during this time interval.
Our results suggest that the cooling during the K-Pg transition may have been responsible for the increase in size of the seed cone of Mesocyparis and have driven the migration of plants southwards.
The cooling from the Maastrichtian to the Paleocene that occurred during the Cretaceous-Paleogene (K-Pg) transition (ca. 70–60 Ma), has been recorded by land plant leaf data  and marine oxygen isotope records . These dramatic environmental and climate changes [3,4,5] altered the composition of terrestrial ecosystems and the evolutionary trends of numerous organisms. Examples include the extinction of non-avian dinosaurs and the rise of mammals, the decline of gymnosperms, and the upsurge of angiosperms [6,7,8]. Some zoological cases indicated how animal taxa changed their morphological features during the K-Pg transition [8,9,10]. For instance, a previous study found that the maximum body size and dental complexity disparity of multituberculate mammals across the K-Pg boundary was associated with dietary expansion . On the other hand, botanical evidence revealed regional vegetation change during this time interval [11, 12] and provided estimates of paleotemperature in both hemispheres [2, 12,13,14]. Leaf physiognomic analysis showed that the 13 °C mean annual temperature isotherm (mesothermal-microthemal boundary) of high latitudinal western North America moved southwards from ca. 65–75 °N in the Maastrichtian to ca. 55 °N in the Danian [15, 16]. Besides, a study focused on the range changes of Metasequoia around the world since the Late Cretaceous has shown the knowledge of the distributional changes of this genus across the K-Pg transition . Nevertheless, it has remained relatively unknown how specific plant taxa responded to the cooling during the K-Pg transition.
The extinct genus Mesocyparis McIver & Basinger (Cupressaceae) was distributed in eastern Asia and western North America from the Late Cretaceous to the Paleocene, ca. 70–60 Ma [18,19,20,21,22]. Here, we report new fossil material of this genus from the upper part of the Wuyun Formation, Jiayin County, Heilongjiang Province. The new Mesocyparis fossils from Northeast China, together with the other four species of this genus in the Northern Hemisphere, provide us with an opportunity to evaluate adaptions within the genus.
Genus: Mesocyparis McIver et Basinger.
Species: Mesocyparis sinica Y. M. Cui, W. Wang, D. K. Ferguson, J. Yang et Y. F. Wang, sp. nov.
Diagnosis: Seed cones ovate to orbicular, borne in opposite pairs (Fig. 1a and Fig.2a). Cone scales 4 in number, woody, obovate to ovate, decussate, approximately equal in size (Fig. 1b, c and Fig. 2b). Umbo near apex of scale, leaflike, erect, with apex tapered (Fig. 1b, c: arrow 2, 2a and 2b). Branches pinnate and opposite, forming flat frond-like sprays (Fig. 1d). Leaves on both vegetative and fertile shoots decussate, scaly and dimorphic with facial and lateral leaves (Fig. 1a, d, Fig. 2a and c).
Holotype: compression fertile shoot with two mature seed cones and two immature seed cones, Palaeobotanical Museum of China, Institute of Botany, Chinese Academy of Sciences (IBCAS), specimen no. WY0801 (Fig. 1)a, b and c.
Etymology: The specific name, sinica, refers to the fossil discovery in China.
Type locality: Wuyun coalmine of Jiayin County, Heilongjiang Province, NE China (49° 14′ N, 129° 28′ E; Fig. 3).
Horizon and age: Wuyun Formation, Danian, Paleocene.
Description: Seed cones borne in opposite pairs, ovate, 4.4–5.2 mm long, 4.0–5.1 mm wide (Fig. 1a). Bracts 4 in number, woody, decussate, obovate to ovate, 3.7–4.3 mm long, 2.8–3.1 mm wide, and sizes of upper pair and basal pair nearly identical (Fig. 1b and c). Umbo leaflike, erect, at top of scale, with apex acuminate (Fig. 1b and c: arrow 2). Pair of immature seed cones arising from axis of fertile branch measure, 1.4–1.5 mm long, 0.9–1.3 mm wide (Fig. 1a: arrow 1). Immature cone scales peltate, 0.7–1.5 mm long, 0.8–0.9 mm wide, apex abrupt to blunt.
Fertile shoots flattened, leaves decussate and dimorphic, facial and lateral leaves morphologically different (Fig. 1a). Facial leaves triangular to rhombic, 0.9–2.7 mm long, 0.7–1.2 mm wide, thickened along midrib. Tips of leaves are acute to blunt, rising above basal part of subsequent leaves. Lateral leaves folded, falcate to triangular, 1.2–2.2 mm long, 0.4–0.8 mm wide, apex acute, rising above basal part of subsequent lateral leaves (Fig. 1a).
Branching pinnate, flattened, branchlets plagiotropic, borne in opposite pairs, forming flat frond–like sprays, about 3.2 cm long. Basal pair of facial leaves arising from the axils of branch lateral leaves. Branchlets on same side parallel, 1.9–2.8 mm apart (Fig. 1d). Scale-like leaves dimorphic, facial leaves adpressed, lateral leaves folded along midrib. There are some morphological differences between leaves of branches and branchlets (Fig. 1e and f). Facial leaves of branchlets ovate to rhombic, 1.1–2.1 mm long, 0.8–1.6 mm wide, adpressed, thickened along midrib, apex acuminate to acute, rising above basal part of subsequent leaves. Lateral leaves of branchlets folded, falcate, 0.9–2.4 mm long, 0.4–0.8 mm wide, apex acute; most of them overlap basal part of subsequent lateral leaves (Fig. 1e). Facial leaves of branches 1.8–3.8 mm long, 0.4–0.6 mm wide, apex acuminate to acute, thickened slightly, rising slightly above basal part of subsequent leaves. Lateral leaves of branches folded, falcate, 1.7–2.7 mm long, 0.4–0.5 mm wide, with apex acute, and basal part covered by earlier facial leaves (Fig. 1f).
Comparison between M. sinica and other species of Mesocyparis
The newly-found fossils from Jiayin and other Mesocyparis records (Figs. 1 and 2) display significant differences (Fig. 2 and Table 1): The seed cones in the Jiayin fossil were erect, with an acuminate apex, while those of M. borealis from Saskatchewan, Canada were reflexed, with an acute apex . Seed cones from Jiayin were borne in opposite pairs, with cone scales ovate to obovate, and umbo near the apex of the cone scale, while the leaves of fertile shoots were scaly and dimorphic. In contrast, the seed cones of M. umbonata from Alberta, Canada were decussate, with peltate cone scales bearing a central umbo, and leaves of fertile shoots monotypic . The cone scales of the Jiayin specimens are nearly equal in size, while the upper pair of cone scales of M. beringiana from Koryak, Russia was larger than the basal pair, while the apex of the umbo was acute. At the same time, the seed cones and foliage of M. beringiana were smaller than those of the other three species and the Jiayin fossils [19, 20, 22]. Most cone scales of M. rosanovii from the Amur Region, Russia were arranged in decussate pairs, only a few being borne in opposite pairs. Moreover, the cone scales were peltate, with an acute umbo apex, and the leaves of the fertile shoots monotypic.
Here we provide a key to the fossil species of Mesocyparis based on cone characters as below:
Phylogenetic position of Mesocyparis
Here we treat the Wuyun fossils as a new species based on a detailed morphological comparison, which shows two or more significant morphological differences between the Wuyun fossils and other Mesocyparis species (Fig. 2, Table 1 and the Key in “Comparison between M. sinica and other species of Mesocyparis” section). The full fossil description and a morphological comparison with the other four species in Mesocyparis are provided in supplementary text of the electronic supplementary material.
Our family-wide phylogenetic analysis indicates that Mesocyparis is a member of the Cupressoideae and clusters with Hesperocyparis, Callitropsis–Xanthocyparis, and Juniperus–Cupressus (BS = 51%) (Fig. 4). Mesocyparis shares whorled cone scales equal (Additional file 1: Char. 30:1 and Table S1), and ovule shape obpyriform (Additional file 1: Char. 34:2 and Table S1), with these five genera. The phylogenetic tree of Mesocyparis is shown in Fig. 5, Additional file 2: Figure. S1 and Additional file 3: Table S2. In Mesocyparis, the two western North American species, M. borealis and M. umbonata, formed a clade (WNA clade, BS = 70%), and the three Asian species, M. beringiana, M. rosanovii and M. sinica, formed another (EA clade, BS = 65%). Mesocyparis sinica is sister to M. rosanovii (BS = 62%).
Distributional changes during the K-Pg transition
Paleomaps  indicate that the landmasses of Eastern Asia and western North America were connected from the Late Cretaceous to the Paleocene. However, North America in the Late Cretaceous was divided by the Western Interior Seaway, while in the Paleocene western and eastern North America became contiguous due to the retreat of this seaway. Fossil records show Mesocyparis ranged from 59°N–76°N paleolatitude (Horseshoe Canyon, Alberta, Canada  and Kyoyak, Russia [18,19,20] during the Late Cretaceous (Fig. 6a) while in the Paleocene it ranged from 52 to 72°N paleolatitude (Amur, Russia  and Heilongjiang, China (this study) in Asia and Wyoming & Alaska, USA and Alberta, Canada  of western North America) (Fig. 6b). This Mesocyparis distribution (Fig. 6), combined with phylogenetic analysis (Fig. 5), reveals that the genus probably first occurred at high latitudes in the Northern Hemisphere and then divided into 2 clades —— the East Asian (EA) clade and the Western North American (WNA) clade.
Morphological changes during the K-Pg transition
When plotting the scored 22 morphological characters against the species-level tree of Mesocyparis, we found some significant morphological changes in the genus from Maastrichtian to Paleocene: 1) seed cones became larger (average cone size increased about 130%, detail estimates in Methods), 2) the umbo shape on the scale changed from erect (Fig. 2h and k) to reflexed (Fig. 2e), and 3) the umbo position on the cone scales moved from the middle (Fig. 2h and k) to the apex (Fig. 2e).
In general, the increasing size of the seed cone in Mesocyparis may be associated with the climate cooling from Maastrichtian to Paleocene across the K-Pg boundary. A negative correlation between the cone size and mean annual temperature (MAT) is apparent in modern Cupressus, one of the nearest living relatives of extinct Mesocyparis (Fig. 7 and Additional file 4: Table S3). In addition, it has been observed that fruits and seeds of angiosperms in North America and Europe increased during the K-Pg transition .
Specific morphological features can be linked to a specific biological function, and may allow us to speculate about the reason behind Mesocyparis seed cones more than doubling in size across the K-Pg transition. The increasing size of seed cones might augment their visibility and attractiveness for herbivorous and/or omnivorous animals, thereby enhancing seed dispersal, which could expand the plant distribution and increase their ability to survive temporal and spatial environmental change. Today some bird species (such as azure-winged magpie and light-vented bulbul) feed on intact seed cones and we envisage similar feeding strategies existed in early birds during the K-Pg transition. The increasing size of the seed cone might be an adaptive mechanism for seed dispersal (i.e., conspicuously larger seed cones were more attractive for those cone-eating/removal animals). An animal’s preference for larger seed cones might select and even amplify the increasing size of seed cones in Mesocyparis. Evidence has shown that the rapid radiation and diversification of birds [25, 26] and mammals  that, occurred during this period produced numerous potential vectors for seed dispersal of both gymnosperms and angiosperms.
Cooling during the K-Pg transition drove Mesocyparis southwards
In general, the geographical distribution of plants at high latitudes is strongly constrained by the terrestrial thermal regime, particularly the mean coldest monthly temperature (MCMT) and/or the mean annual temperature (MAT) [28, 29]. Thus, the northern limit of a specific terrestrial plant taxon can reflect temperature changes with relatively high sensitivity .
The cooling during the K-Pg transition (ca. 70–60 Ma, global average temperature fell by ca. 4 °C)  could have driven the northern limit of Mesocyparis south by ca. 4° N from the Maastrichtian to the Paleocene (Fig. 8 and Additional file 5: Table S4), which means that the northern limit of survival for the plant retreated southwards by ca. 444 km. Similar responses are found in Metasequoia (Cupressaceae) and Nordenskioldia (Trochodendraceae). Metasequoia’s northern limit moved south by ca. 4 ° (Fig. 8 and Additional file 6: Table S5), while Nordenskioldia’s moved almost 5 ° (Fig. 8 and Additional file 7: Table S6). These findings indicate that the southward migration of the northern limit of Mesocyparis as driven by the cooling during the K-Pg transition is by no means unique.
We report a new fossil species, Mesocyparis sinica from the early Paleocene of Jiayin, Heilongjiang, northeastern China. By integrating lines of evidence from phylogeny and the comparative morphology of Mesocyparis, we found that during ca. 70–60 Ma, the size of the seed cone of Mesocyparis increased by ca. 130%, probably driven by the cooling during the K-Pg transition, which might be an effective adaptation for seed dispersal by animals. Meanwhile, our analyses indicate that the northern limits of this genus, as well as those of two other arboreal taxa Metasequoia (gymnosperm) and Nordenskioldia (angiosperm), have migrated ca. 4–5° southwards in paleolatitude during this interval, in association with the cooling. This is the first report to link plant evolutionary and ecological responses to the climate cooling during the K-Pg transition. This study not only enriches our understanding of how specific plant taxa responded to the climate changes during the K-Pg transition, but also provides a reversed case study for the scientific community to speculate about plant responses and distributional changes under global warming today.
The specimens of Mesocyparis sinica described here were collected from the upper part of the Wuyun Formation, at Wuyun coalmine (49° 14′ N, 129° 28′ E; Fig. 3) in Jiayin County, Heilongjiang Province, NE China. The enclosing sediments of the Jiayin fossil are grey sandy mudstone. Although some authors [32, 33] suggested it was of Maastrichtian (Late Cretaceous) age, the Wuyun Formation is now considered to belong to the Danian stage of the Early Paleocene based on palynotology [34,35,36], vegetation composition [37,38,39,40] and stratigraphic correlation [41,42,43]. In addition, Single zircon LA-ICP-MS dating shows the age of the lower part of the Wuyun Formation to be 64.1 ± 0.7 Ma [44, 45].
Materials in this study
The fossil specimens of Mesocyparis sinica examined in this article are deposited in the National Museum of Plant History of China, Institute of Botany, Chinese Academy of Sciences: WY0801 (Figs. 1a, b and c), WY0802 (Figs. 1d, e and f), WY0803. These specimens were exposed and observed under a stereomicroscope, which identified by Dr. Yi-Ming Cui and Prof. Yu-Fei Wang. The characters were described following the terminology of . Attempts were made to extract structure of seed cones under micro-CT and cuticles of leaf, but these efforts were unsuccessful.
Morphological and distributional information for the other four species of extinct Mesocyparis were collected from literature: M. borealis from , M. umbonata from , M. beringiana from [19, 20, 22], and M. rosanovii from .
The morphological and distributional information of extant Cupressus were measured and collected from specimens in the Herbarium of the Institute of Botany, Chinese Academy of Sciences, Beijing (PE) (Additional file 1: Table S1).
We first used the “backbone constraint” tree approach  to determine the phylogenetic position of Mesocyparis in the Cupressaceae s.l. The constrained tree is modified from the two trees of Mao et al.  and Yang et al. , in which significant conflicting nodes (BS ≥ 85%) between the two trees were collapsed. The morphological matrix was modified from Farjon  and included 53 characters (Additional file 1: Table S1). Following the family-wide analysis, we selected Juniperus monticola, Xanthocyparis vietnamensis, and Hesperocyparis arizonica as outgroups and built the species-level relationships of Mesocyparis. All five species of Mesocyparis were included. The 22 morphological characters were coded based on the literature [18,19,20,21,22] and our observations (Additional file 3: Table S2). Phylogenetic analysis for each matrix was conducted using maximum parsimony method in PAUP* version 4.0b10 . Heuristic searches were performed with 1000 random sequence addition replicates, tree-bisection-reconnection branch swapping, MulTrees in effect, and steepest descent off. Bootstrapping was carried out with 1000 replicates, using a heuristic search strategy.
Estimates of average cone size changes during the K-Pg transition
Here we use a volume formula of an ellipsoid to estimate the average cone size of each species in Mesocyparis, for the shapes of the seed cones in this genus are ovoid to ellipsoidal. Therefore, we set the seed cone length (scl) as the polar-radius of the ellipsoid and seed cone width (scw) as the two equatorial-radii. The estimated volume (V) of Mesocyparis is expressed by the following ellipsoid volume equation:
where scl = median value of seed cone length, scw = median value of seed cone width. The scls of five species in Mesocyparis are 3.9 mm (M. umbonata), 3 mm (M. beringiana), 4 mm (M. borealis), 4.25 mm (M. rosanovii) and 4.8 mm (M. sinica), while their scws are 2.95 mm (M. umbonata), 3 mm (M. beringiana), 4 mm (M. borealis), 3.5 mm (M. rosanovii) and 4.55 mm (M. sinica). Based on the volume equation, the estimated volumes are 5.657 mm3 (M. umbonata), 4.5 mm3 (M. beringiana), 10.667 mm3 (M. borealis), 8.677 mm3 (M. rosanovii) and 16.562 mm3 (M. sinica). Thus, we can estimate that the average volumes of Mesocyparis in the Maastrichtian and the Paleocene were 5.079 mm3 and 11.969 mm3 respectively, so the increase in volume of Mesocyparis from the Maastrichtian to the Paleocene was about 136%.
Correlational analysis of seed cone sizes of Cupressus and temperature
To explore the influence of climate change on the morphology of seed cones in Mesocyparis, we calculated the potential relationship between the seed cone sizes of extant Cupressus and temperatures in their habitats. We examined 45 samples of five Cupressus species (C. chengiana S. Y. Hu, C. funebris Endl., C. duclouxiana Hickel, C. torulosa D. Don and C. gigantea Cheng et L. K. Fu) naturally occurring in China from PE Herbarium, which have female cones and entire information on locality and elevation. The diameters of the female cones were measured to represent their seed cone sizes. Mean annual temperature (MAT) used the climate data set of New et al.  and then corrected for the altitude of the plant locality using a temperature lapse rate of 5 °C /km (Additional file 4: Table S3).
The coordinates of paleolatitude and paleolongitude were compiled from previous reports of Mesocyparis [18,19,20,21,22], Metasequoia (Cupressaceae, gymnosperm)  and Nordenskioldia (Trochodendraceae, angiosperm)  in the literature, and converted from today’s latitude and longitude to those inferred for the Cretaceous and Paleocene using “Point Tracker for windows” software, to plot the sites on paleogeographic maps representing two intervals of geologic time. The base maps with paleocoastlines were modified from paleogeographic maps of the late Cretaceous (ca. 94 Ma) and the Paleocene (ca. 60 Ma) from Scotese .
Availability of data and materials
All data generated or analysed during this study are included in this published article and its additional files.
Institute of Botany, Chinese Academy of Sciences
Cretaceous – Paleogene
Megaannus = Million years ago
Mean annual temperature
Mean coldest monthly temperature
Herbarium of the Institute of Botany, Chinese Academy of Sciences, Beijing
- s.l. :
- scl :
Seed cone length
- scw :
Seed cone width
- V :
Western North American
Wilf P, Johnson KR, Huber BT. Correlated terrestrial and marine evidence for global climate changes before mass extinction at the cretaceous-Paleogene boundary. Proc Natl Acad Sci U S A. 2003;100:599–604.
Li L, Keller G. Variability in late cretaceous climate and deep waters: evidence from stable isotopes. Mar Geol. 1999;161:171–90.
Schulte P, Alegret L, Arenillas I, Arz JA, Barton PJ, Bown PR, et al. The Chicxulub asteroid impact and mass extinction at the cretaceous-Paleogene boundary. Science. 2010;327:1214–8.
Renne PR, Sprain CJ, Richards MA, Self S, Vanderkluysen L, Pande K. State shift in Deccan volcanism at the cretaceous-Paleogene boundary, possibly induced by impact. Science. 2015;350:76–8.
Chiarenza AA, Mannion PD, Lunt DJ, Farnsworth A, Jones LA, Kelland SJ, et al. Ecological niche modelling does not support climatically-driven dinosaur diversity decline before the cretaceous/Paleogene mass extinction. Nat Commun. 2019;10:1091.
Chen M, Strömberg CAE, Wilson GP. Assembly of modern mammal community structure driven by late cretaceous dental evolution, rise of flowering plants, and dinosaur demise. Proc Natl Acad Sci. 2019;116:9931–40.
Magallón S, Hilu KW, Quandt D. Land plant evolutionary timeline: gene effects are secondary to fossil constraints in relaxed clock estimation of age and substitution rates. Am J Bot. 2013;100:556–73.
Wilson GP, Evans AR, Corfe IJ, Smits PD, Fortelius M, Jernvall J. Adaptive radiation of multituberculate mammals before the extinction of dinosaurs. Nature. 2012;483:457–60.
Sibert EC, Norris RD. New age of fishes initiated by the cretaceous-Paleogene mass extinction. Proceedings of the National Academy of Sciences, USA. 2015;112:8537–42.
Grossnickle DM, Newham E. Therian mammals experience an ecomorphological radiation during the late cretaceous and selective extinction at the K-Pg boundary. Proc R Soc B Biol Sci. 2016;283:20160256.
Tschudy RH, Tschudy BD. Extinction and survival of plant life following the cretaceous/tertiary boundary event, Western interior, North America. Geology. 1986;14:667–70.
Spicer RA, Collinson ME. Plants and floral change at the cretaceous-Paleogene boundary: three decades on. Special Paper of the Geological Society of America. 2014;505:117–32.
Wolfe JA, Upchurch GR. Leaf assemblages across the cretaceous-tertiary boundary in the Raton Basin, New Mexico and Colorado. Proceedings of the National Academy of Sciences, USA. 1987;84:5096–100.
Bowman VC, Francis JE, Askin RA, Riding JB, Swindles GT. Latest cretaceous–earliest Paleogene vegetation and climate change at the high southern latitudes: palynological evidence from Seymour Island, Antarctic peninsula. Palaeogeogr Palaeoclimatol Palaeoecol. 2014;408:26–47.
Wolfe JA, Upchurch GR. North American nonmarine climates and vegetation during the late cretaceous. Palaeogeogr Palaeoclimatol Palaeoecol. 1987;61:33–77.
Davies-Vollum KS. Early Palaeocene palaeoclimatic inferences from fossil floras of the western interior, USA. Palaeogeogr Palaeoclimatol Palaeoecol. 1997;136:145–64.
Liu YJ, Arens NC, Li CS. Range change in Metasequoia: relationship to palaeoclimate. Bot J Linn Soc. 2007;154:115–27.
McIver EE, Basinger JF. Mesocyparis borealis gen. Et sp. nov.: fossil Cupressaceae from the early tertiary of Saskatchewan, Canada. Can J Bot. 1987;65:2338–51.
Golovneva LB. New genus Microconium (Cupressaceae) from the late cretaceous deposits of the north-east of the USSR. Botanicheskii Zhurnal. 1988;73:1179–83 (in Russian).
Moiseeva MG. The Maastrichtian flora of the Amaam lagoon area (northeastern Russia). Stratigr Geol Correl. 2012;20:579–679.
Kodrul TM, Tekleva MV, Krassilov VA. A new conifer species, Mesocyparis rosanovii sp. nov. (Cupressaceae, Coniferales), and Transberingian floristic connections. Paleontol J. 2006;40:328–38.
McIver EE, Aulenback KR. Morphology and relationships of Mesocyparis umbonata sp.nov.: fossil Cupressaceae from the late cretaceous of Alberta, Canada. Can J Bot. 1994;72:273–95.
Scotese CR. Atlas of earth history. Arlington: PALEOMAP Project; 2001. ISBN:0–9700020–0–9
Eriksson O, Friis EM, Löfgren P. Seed size, fruit size, and dispersal systems in angiosperms from the early cretaceous to the late tertiary. Am Nat. 2000;156:47–58.
O'Leary MA, Bloch JI, Flynn JJ, Gaudin TJ, Giallombardo A, Giannini NP, et al. The placental mammal ancestor and the post-K-Pg radiation of placentals. Science. 2013;341:662–7.
Jarvis ED, Mirarab S, Aberer AJ, Li B, Houde P, Li C, et al. Whole-genome analyses resolve early branches in the tree of life of modern birds. Science. 2014;346:1320–31.
Ksepka DT, Stidham TA, Williamson TE. Early Paleocene landbird supports rapid phylogenetic and morphological diversification of crown birds after the K-Pg mass extinction. Proc Natl Acad Sci U S A. 2017;114:8047–52.
Woodward FI. Climate and Plant Distribution. Cambridge: Cambridge University Press; 1987.
Criddle RS, Hopkin MS, Mcarthur ED, Hansen LD. Plant distribution and the temperature coefficient of metabolism. Plant Cell Environ. 1994;17:233–43.
Beerling DJ. The impact of temperature on the northern distribution limits of the introduced species Fallopia japonica and Impatiens glandulifera in north-West Europe. J Biogeogr. 1993;20:45–53.
Kergoat GJ, Bouchard P, Clamens AL, Abbate JL, Jourdan H, Jabbour-Zahab R, et al. Cretaceous environmental changes led to high extinction rates in a hyperdiverse beetle family. BMC Evol Biol. 2014;14:220.
Hsu J. Late cretaceous and Cenozoic vegetation in China, emphasizing their connections with North America. Ann Mo Bot Gard. 1983;70:490–508.
Tao JR, Xiong XZ. The latest cretaceous flora of Heilongjiang Province and the floristic relationship between East Asia and North America. Acta Phytotaxonomica Sinica. 1986;24:1–15.
Markevich VS, Bugdaeva EV, Ashraf AR, Sun G. Boundary of cretaceous and Paleogene continental deposits in Zeya-Bureya Basin, Amur (Heilongjiang) river region. Global Geology. 2011;14:144–59.
Markevich V, Bugdaeva E, Nichols D, Sun G. Paleogene coal-forming plants of the Zeya-Bureya Basin (Amur River region). Paleontol J. 2010;44:1321–31.
Hao H, Ferguson DK, Feng GP, Ablaev A, Wang YF, Li CS. Early Paleocene vegetation and climate in Jiayin. NE China Climatic Change. 2010;99:547–66.
Sun G, Dong ZM, Akhmetiev M, Markevich V, Ashraf AR, Godefroit P, et al. Late Cretaceou- Paleocene biota and the K-Pg boundary, with discussion on extinction of dinosaurs. Shanghai: Shanghai Science and Technology Press; 2014. ISBN:978–7–5428–6111–5.
Sun G, Akhmetiev M, Markevich V, Ashraf AR, Bugdaeva E, Godefroit P, et al. Late cretaceous biota and the cretaceous-Paleocene (K-Pg) boundary in Jiayin of Heilongjiang, China. Global Geology. 2011;14:115–43.
Sun G, Quan C, Sun CL, Sun YW, Luo KL, Lu JS. Some new knowledge on subdivision and age of Wuyun Formation in Jiayin of Heilongjiang, China. J Jilin Univ (Earth Sci Ed). 2005;35:137–42.
Quan C, Sun G, Zhou ZY. A new tertiary Ginkgo (Ginkgoaceae) from the Wuyuan formation of Jiayin, Heilongjiang, northeastern China and its paleoenvironmental implications. Am J Bot. 2010;97:446–57.
Wang YF, Li CS, Li ZY, Fu DZ. Wuyunanthus gen. Nova flower of Celastraceae from the Palaeocene of north-East China. Bot J Linn Soc. 2001;136:323–7.
Feng GP, Ablaev AG, Wang YF, Li CS. Paleocene Wuyun Flora in Northeast China: Ulmus furcinervis of Ulmaceae. Acta Phytotaxonomica Sinica. 2003;45:146–51.
Wang Q, Ablaev AG, Wang YF, Li CS. Paleocene Wuyun flora in Northeast China: Woodwardia bureiensis, Dryopteris sp. and Osmunda sachalinensis. Acta Phytotaxonomica Sinica. 2006;44:712–20.
Shigeyuki S, Sun G, Knittel U, Tohru D, Hideki I. Radiometric zircon ages of a tuff sample from the Baishantou member of Wuyun formation, Jiayin: a contribution to the search for the K-T boundary in Heilongjian River area. China Acta Geologica Sinica. 2011;85:1351–8.
Knittel U, Suzuki S, Akhmetiev MA, Sun G, Walia M, Hung CH. 66±1 ma single zircon U-Pb date confirms the location of the non-marine K-Pg boundary in the Amur/Heilongjiang River area, (Russia, China). Neues Jahrb Geol Palaontol Abh. 2013;270:1–11.
Farjon A. A monograph of Cupressaceae and Sciadopitys. London: RBG Kew; 2000.
Endress PK, Doyle JA. Reconstructing the ancestral angiosperm flower and its initial specializations. Am J Bot. 2009;96:22–66.
Mao K, Milne RI, Zhang L, Peng Y, Liu J, Thomas P, et al. Distribution of living Cupressaceae reflects the breakup of Pangea. Proc Natl Acad Sci U S A. 2012;109:7793–8.
Yang ZY, Ran JH, Wang XQ. Three genome-based phylogeny of Cupressaceae s.l. : further evidence for the evolution of gymnosperms and southern hemisphere biogeography. Mol Phylogenet Evol. 2012;64:452–70.
Swofford DL, PAUP*. Phylogenetic analysis using parsimony (*and other methods), version 4.0b10. Sunderland, MA: Sinauer Associates; 2003.
New M, Lister D, Hulme M, Makin I. A high-resolution data set of surface climate over global land areas. Clim Res. 2002;21:1–25.
Wang YH, Ferguson DK, Feng GP, Wang YF, Zhilin SG, Li CS, Svetlana PT, Yang J, Ablaev AG. The phytogeography of the extinct angiosperm Nordenskioeldia (Trochodendraceae) and its response to climate changes. Palaeogeogr Palaeoclimatol Palaeoecol. 2009;280:183–92.
Wang M, Lloyd GT. Rates of morphological evolution are heterogeneous in early cretaceous birds. Proc R Soc B Biol Sci. 2016;283:20160214.
We thank Prof. Robert A. Spicer for correcting the manuscript, Prof. Dianne Edwards for reading the previous versions of the manuscript and Ying-Bao Sun for help in drawing the seed cones of Mesocyparis.
This research was supported by the Strategic Priority Research Program of Chinese Academy of Sciences (No. XDB26000000, XDB310301) for collecting materials, the National Natural Science Foundation of China (No. 31570223, 41790454, 31770231, and 31470315) for doing experiments and analyzing data, the Basic Project of the State Key Laboratory of Palaeobiology and Stratigraphy (20192101) and the Chinese Academy of Sciences President’s International Fellowship Initiative (No. 2018VBA0016) for writing and revising the manuscript.
Ethics approval and consent to participate
The fossil materials of Mesocyparis sinica used in this study were collected by Prof. Yu-Fei Wang in Jiayin County, Heilongjiang, China. According to both national and local guidelines, there’re no specific permissions required for collecting and using the fossils in the study.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Morphological characters matrix including both extant and extinct species of the Cupressaceae s.l. (DOCX 34 kb)
Fig. S1. The strict consensus tree of Mesocyparis inferred from the 22 morphological characters. (TIF 337 kb)
Table S2. Morphological characters matrix of Mesocyparis and its outgroups. (DOCX 17 kb)
Table S3. The seed cones sizes of Cupressus, their geographic information and mean annual temperature (MAT) of their habitats. (DOCX 19 kb)
Table S4. Fossil localities of Mesocyparis and the estimates of their paleo-latitudes and palaeo-longitudes. (DOCX 15 kb)
Table S5. Fossil localities of Metasequoia and the estimates of their paleo-latitudes and palaeo-longitudes. (DOCX 20 kb)
Table S6. Fossil localities of Nordenskioldia and the estimates of their paleo-latitudes and paleo-longitudes. (DOCX 16 kb)
About this article
Cite this article
Cui, YM., Wang, W., Ferguson, D.K. et al. Fossil evidence reveals how plants responded to cooling during the Cretaceous-Paleogene transition. BMC Plant Biol 19, 402 (2019). https://doi.org/10.1186/s12870-019-1980-y