Sika deer are known to prefer oak leaves, which are rich in tannins and toxic to most mammals; however, the genetic mechanisms underlying their unique ability to adapt to living in the jungle are still unclear. In identifying the mechanism responsible for the tolerance of a highly toxic diet, we have made a major advancement by explaining the genome of sika deer. We generated the first high-quality, chromosome-level genome assembly of sika deer and measured the correlation between tannin intake and RNA expression in 15 tissues through 180 experiments. Comparative genome analyses showed that the UGT and CYP gene families are functionally involved in the adaptation of sika deer to high-tannin food, especially the expansion of the UGT family 2 subfamily B of UGT genes. The first chromosome-level assembly and genetic characterization of the tolerance to a highly toxic diet suggest that the sika deer genome may serve as an essential resource for understanding evolutionary events and tannin adaptation. Our study provides a paradigm of comparative expressive genomics that can be applied to the study of unique biological features in non-model animals.

The conservation of species diversity and improvement of forest structure are essential roles of the Natural Reserve Policy and the Natural Forest Protection Program (NFPP) in China. However, the long-term effects of NFPP are still not well-defined, and a natural reserve (Liangshui) and surrounding region were surveyed as a proxy of NFPP for approaching the protection effects. Our results showed that long-term conservation significantly altered the dominant species in the herb layer (80% of species), followed by shrub (58%) and tree layers (50%); there was a 1.6-8.0-fold increase in abundance in Corylus shrubs, Acer trees and Carex grass, but a 1.3–10.0-fold abundance decrease in larch trees, Athyrium herbs and Lonicera shrubs. In contrast, tree species diversity and distribution evenness increased by 31% and 23.4% in the reserve, respectively. Forest protection in the reserve also led to the forest structural alteration with the observation of larger-sized trees and shorter herbs, but relatively sparse forests (smaller tree density). Structural equation modeling manifested that the reserve directly altered forest structure, at a coefficient of 0.854, nearly two-fold higher than its impact on diversity (0.459) and dominant species (−0.445). The most affected parameters were plant size (trees and herbs) and tree density related to forest structure, tree diversity, herb richness and evenness for diversity traits, and Oxalidaceae and Rosaceae for dominant species. This study provides basic data that can be used to evaluate the impact of the nature reserve in NE China, and these findings can be used to guide the implementation of NFPP in the long-term in the future.

To explore how climatic factors influence tree growth within the context of global climate changes, we used a dendroclimatological analysis to understand the response of Larix gmelinii to climatic variations along a broad latitudinal gradient from 47.27° to 52.66° N in the Greater Khingan Mountains of Northeastern China. The growth-limiting climate factors and a detailed association between radial growth and climate were identified along the gradient using redundancy analysis (RDA) and standard correlation function analysis over the period 1960–2013. The results showed that temperatures during current June to July represented the most important factor affecting tree radial growth in the study area. Across all studied latitudes, Larix gmelinii growth might be decreasing in radial growth due to higher monthly maximum temperature (Tmax) and monthly mean temperatures (Tm) in the current June, especially for the stands at low and middle latitudes. With continued warming, Larix gmelinii radial growth at high latitudes (e.g., Mangui (MG) and Mohe (MH)) might be reduced by warmer temperatures in July. In addition, Larix gmelinii might be decreasing in radial growth from decreasing precipitation. Our results show that there is a decreasing trend in Larix gmelinii radial growth under the observed general increase of temperatures in the Greater Khingan Mountains in recent years.

Association is the basic unit of plant community classification. Exploring the distribution of plant associations can help improve our understanding of biodiversity conservation. Different associations depend on different habitats and studying the association level is important for ecological restoration, regional ecological protection, regulating the ecological balance, and maintaining biodiversity. However, previous studies have only focused on suitable distribution areas for species and not on the distribution of plant associations. Larix gmelinii is a sensitive and abundant species that occurs along the southern margin of the Eurasian boreal forests, and its distribution is closely related to permafrost. In this study, 420 original plots of L. gmelinii forests were investigated. We used a Maxent model and the ArcGIS software to project the potential geographical distribution of L. gmelinii associations in the future (by 2050 and 2070) according to the climate scenarios RCP 2.6, RCP 4.5, and RCP 8.5. We used the multi‐classification logistic regression analysis method to obtain the response of the suitable area change for the L. gmelinii alliance and associations to climate change under different climate scenarios. Results revealed that temperature is the most crucial factor affecting the distribution of L. gmelinii forests and most of its associations under different climate scenarios. Suitable areas for each association type are shrinking by varying degrees, especially due to habitat loss at high altitudes in special terrains. Different L. gmelinii associations should have different management measures based on the site conditions, composition structure, growth, development, and renewal succession trends. Subsequent research should consider data on biological factors to obtain more accurate prediction results.