Macaques in China: Evolutionary dispersion and subsequent development

Abstract Depicting a taxonomic group's evolutionary trajectory as a function of changes in the geographical landscape and its historical distribution is critical for constructing informed conservation strategies. Based on fossil sites from the Pliocene to the Holocene, and historical records since 1175 AD, we established macaques’ dispersal pathways into and through China. These routes include internal pathways starting from the southeast corner of the Qinghai‐Tibet Plateau and Mts. Hengduan in western China, and the routes through the estuaries of the three major rivers (Yangtze, Yellow, and Pearl). Our results indicate that macaques used the three rivers and avoided the higher elevation of the plateaus to promote their radiation. They occupied the whole mainland and islands in the Pleistocene and experienced shrunken distribution in the Holocene due to climate changes and human‐induced activities. A prominent China‐wide reduction occurred between 1817 and 1917; and a remarkable retraction from central China happened between 1918 and 2018 following further eco‐social development and human expansion in central China, particularly since the second half of the last century. Starting in 1175 there was a restriction of range to higher altitudes, so that macaques have contracted their range to the west, and the Qinghai‐Tibet Plateau and Mts. Hengduan have become an important sanctuary. We predict that if the current climate and human‐induced changes are not reversed by decisive conservation actions, macaques in east and central China will likely be extinct in the near future.

pathways into and through China. These routes include internal pathways starting from the southeast corner of the Qinghai-Tibet Plateau and Mts. Hengduan in western China, and the routes through the estuaries of the three major rivers (Yangtze, Yellow, and Pearl). Our results indicate that macaques used the three rivers and avoided the higher elevation of the plateaus to promote their radiation.
They occupied the whole mainland and islands in the Pleistocene and experienced shrunken distribution in the Holocene due to climate changes and human-induced  (Li, Pan, & Oxnard, 2002;Lotze & Worm, 2009;Turvey, Crees, & Di Fonzo, 2015;Zhao et al., 2018). Historical records of the distribution of sufficient duration are only available for a few civilizations. One of these is that of China with its written history dating back to 1,250 years (Li & Hong, 2007).
There are currently 28 species of nonhuman primates in China.
Such hypothesized pathways are, unfortunately, not schematically backgrounded with geographic landscapes, such as plateaus, and major waterways in China, so that they cannot tell us specifically how a group selected dispersal ways according to the structure of the landscapes.
Climate changes after the Pleistocene and extensive human activities in the Holocene significantly impacted macaque distributions in China Pan et al., 2016). Currently, M. mulatta and M.  (Fan et al., 2017;Zhao et al, 2015). Macaque populations in China are generally more negatively impacted than those outsides because of China's remarkable development that has been dependent principally on the extraction of the natural resources, including extensive deforestation, particularly since the 1950s Pan et al., 2016).
Here we examine fossil and extant macaque taxa distribution in China to explore the dispersal patterns of Macaca during the Pliocene to Holocene, and to predict future distributions following their current changing trajectories. We discuss possible external dispersal routes used by the sinica and fascicularis groups, and M. leonina from the silenus groups, into China; propose their within China dispersal routes relying on mainland fossil locations; and display the changes of spatiotemporal distribution trajectories of the macaques in recent history.

| MATERIALS AND METHODS
Our project is based on fossil records and a broad literature review from academic journals, governmental annals, magazines, and books, many written in Chinese and not easily accessible to the international community. We particularly relied on the two by books written in Chinese by Rongsheng Wen (Wen, 2009(Wen, , 2013Wen, 2018) that review archives and annals at the county level since 1175 AD, and indicate where large mammals occur. Although different names were given to the same animal groups at different places and historical periods, in most cases this does not influence primate identification (Li et al., 2002;Wen, 2009Wen, , 2013. There were four primate types in China: (a) Gibbons called "Yuan" in Chinese, often described by their long upper limbs, slim body structure, and unique calls (Zhang, 2015).  (Liu & Li, 2013), which was translated into English in 1943 (Arthur, 1943). In this novel "Sun Wukong," the Monkey King, is an incarnation of the macaques and is also mentioned in Confucianism, Buddhism, and Taoism-"Three-Religion-in-One Ideology" (Liu & Li, 2013;Lou, 2016). We (Li and Pan) validated historical name variations with the methods used in Li et al. (2002).
Although government archives and annals started in China more than 2,000 years ago in the Qin Dynasty (Dai, Liu & Xiao, 2007), the detailed documentation of animal distribution started about 400 years ago (Li et al., 2002). Thus, we used the fossil records, and historical records since 1716 to generate distribution profiles for each century (see Supporting Information Database). Evaluation of the fossil records of Chinese macaques was based mainly on book chapters (Jablonski, 2008), book publications (Lucas, 2001;Wu & Poirier, 1995), and journal papers (e.g., Wu, 2004). Longitudes, latitudes, and altitudes of the locations were determined from Google Earth. Commonly, the same records for a given location in each century were documented in different years and by different publications or archives, only one record was used for mapping. Spatiotemporal distribution maps are generated with ESRI ArcGIS Desktop 10.6, and SPSS Version 20 was used to illustrate altitudinal changes from 1175 AD to the present.  Delson (1980) referring to morphology and geographic distribution: Macaca sylvanus found only in North Africa was separated from the sylvanus-silenus group, and M. arctoides was combined with the sinica group, which is supported by the evidence from Y Chromosome loci (Tosi et al., 2000), Alu elements for genetic markers (Li et al., 2009), and homologous regions of DNA (Deinard & Smith, 2001); and (c) Groves (2000) separated M. leonina from M. nemestrina based on morphology, which is supported by molecular study (Roos, Thanh, Walter, & Nadler, 2007). Macaca pagensis also was separated from M. nemestrina (Brandon-Jones et al., 2004;Tosi et al., 2000). Macaca siberu, exclusively on the island of Siberut, the west coast of Sumatra, was separated from M. pagensis by Kitchener and Groves (2002) referring to its isolated distribution and morphology. M. leucogenys was reported in southeastern Tibet (Li, Zhao, & Fan, 2015), and M. munzala was recorded in North Eastern India and Tibet (Sinha et al., 2005). The combination of M. leucogenys, and M. munzala with the sinca group is based on morphology (Li et al., 2015;Sinha et al., 2005), and genetic and molecular evidence (Chakraborty et al., 2007;Fan et al., 2017). Red-colored species are the taxa found in China, and bolded species are the taxa newly categorized 3 | RESULTS  (Abegg & Thierry, 2002;Delson, 1980;Fa, 1989)

| External dispersion
The ancestor of Asian macaques separated from their African counterparts in the Miocene between 5.3 and 5.9 Ma (Alba et al., 2014;Delson, 1975;Fa, 1989;Tosi, Morales, & Melnick, 2003). This has been validated based on the evidence from chromosomes and molecular markers (Tosi et al., 2003), gross morphology (Fooden, 1976), fossils (Delson, 1980;Takai et al., 2014), geographic distribution pattern (Abegg & Thierry, 2002), and genetics (Li & Zhang, 2005 F I G U R E 2 Fossil-bearing locations of the macaques in China discovered in the Pliocene, Pleistocene, and Holocene (see Supporting Information Database for more details). The dispersal patterns of the sinica and fascicularis groups, and M. leonina are illustrated referring to the hypotheses by Delson (1980), Fa (1989), and Tosi et al. (2003), except for entrances of the three major rivers. The squared area is the southeast corner of the Qinghai-Tibet Plateau and Mts. Hengduan in Yunnan (SQPMH) F I G U R E 3 Age differentiation of the fossils among three major Chinese geographic regions. Western China (Chongqing, Gansu, Guizhou, Ningxia, Qinghai, Shaanxi, Sichuan, Xinjiang, Yunnan, and Xizang), Central China (Anhui, Heilongjiang, Henan, Hubei, Hunan, Inner Mongolia, Jiangxi, Jilin, and Shanxi), and Eastern China (Beijing, Fujian, Guangdong, Guangxi, Hainan, Hebei, Jiangsu, Liaoning, Shandong, Shanghai, Tianjin and Zhejiang (Dong & Liang, 2014;Huang, Ma, & Huang, 2017;Wei, 1998). SD: standard deviation Macaques entered China using two routes; one for the sinica group, and M. leonina from the silenus group through the SQPMH, and the other for the fascicularis group along the coastlines of South and East China Seas (Figure 2).
After having reached south India along the coast and spread to Southeast Asia, where the silenus group was originated. The sinica group, however, dispersed to northern India and approached the Himalayas and Qinghai-Tibet Plateau (Sakai et al., 2006), which formed a barrier to its further northern dispersion, so that this ancestral group moved eastward into the SQPMH, and subsequently dispersed eastward China (orange arrows in Figure 2), while some moved southward diverging into M. arctoides in Cambodia; Laos, Malaysia, Myanmar, Thailand, and Viet Nam (Delson, 1980;Fa, 1989;Fooden, 1982). M. leonina from the silenus group in the Southeast Asia moved northward (black arrows in Figure 2), and entered into the SQPMH. Interestingly, different from other taxa, this species still remains there, mostly in Western Yunnan Province .
The dispersal procedure of the fascicularis group was closely related to climate change; the formation of the glaciation led to the lower sea levels during the Pliocene and Pleistocene, and land bridges formed in East Asia at different times (Voris, 2000). At that time Hainan, Taiwan, and Japan were frequently linked to the mainland until approximately 18,000 years ago (Adams & Faure, 1997). This allowed the diversification and led to M. mulatta in Hainan, M. cyclopis in Taiwan, and M. fuscata in Japan (Fa, 1989). A puzzling question is how M. mulatta got into the mainland. That more fossil localities are near the estuaries of the Yangtze, Yellow, and Pearl Rivers in Figure 2 suggests they used the lower parts of the rivers for mainland radiation. The dispersal procedure for the fascicularis group is thought to have occurred between 2 and 0.7 Ma (Abegg & Thierry, 2002;Delson, 1980) from Sundaland. On the other hand, the fact that M. mulatta is widely distributed in East and South Asia suggests that their ancestors also dispersed northward from the Sundaland (red arrows in Figure 2).
That fossils to the east of China are younger than those in the west suggests that the external dispersal for the sinica group may start earlier than the fascicularis group (Figure 3). This suggestion is supported by nucleotide research indicating that Taiwan and Japan were colonized around only 0.38-0.44 Ma ago (Chu, Lin, & Wu, 1999). in the Late Miocene-Early Pliocene due to the orogenic events (Ding & Liu, 1989;Lin et al., 2001). Macaques likely used what is now the Yellow River drainage to cross the river and reached the most northeastern corner of China (Liaoning and Jilin). (b) They also likely used the area between the Yangtze and Pearl Rivers that originate in the SOPMH.
During the glaciation in the Late Pleistocene, the water level of the Yangtze River dropped by more than 20 m in some sections (Yang, 1986), thus macaques were likely able to cross the Yangtze. (c) South of the Pearl River contains many macaque fossils suggesting this was also another dispersal route. After having gotten into the mainland through estuaries of the three rivers, populations of the M. mulatta started their radiation in China.

| Distribution alteration in recent history
Historical records of the last 300 years reveal many changes in macaque distribution in China since 1716 (Figure 4). This was mainly caused by human-induced activities Pan et al., 2016) (Figure 4b,c). In addition to climate change (Jablonski & Whitfort, 1999;Zhang et al., 2007), distribution shifts were also associated with economic and social development, particularly its central China (Zhu et al., 2018), such as the cultivated land area increased most quickly in the early Qing Dynasty (Ge et al., 2004). This resulted from technology development and population growth (Lu, 2015), from 388,150,057 in 1799 (Wang, 1997) to about 600 million in 1950 (http://www.china-profile.com/data/fig_WPP2008_TFR_1. htm). The most remarkable change as shown in Figure 5 was, however, the increased altitudinal scale, resulting in the SQPMH becoming an important sanctuary for them.
China experienced dramatic social and environmental change that is well documented starting from the first dynasty (Shang) more than 3,000 years ago (Wakeman, 1977). Since 1780, China suffered 25 major domestic and international conflicts (Elleman, 2001), including the two Opium Wars in 1834-1842 and 1856-1860 (Feige & Miron, 2008;Hanes & Sanello, 2002). Archives and historical records were likely destroyed or lost during the wars, which may explain some locations not being represented in 1817-1917, but in the previous century (1716-1816). West China, particularly SQPMH, was a sanctuary for plants and animals since the Pleistocene (Jablonski, 1993) and is now a biodiversity spot (Tang et al., 2006), and contrast east and central China, it is a relatively lower level of human occupancy, so that some macaque taxa moved back to the Qinghai-Tibet Plateaus and the SOPMH.

| CONCLUSION
Based on the records of fossil-bearing sites from the Pliocene to Holocene, and historical distribution records since 1175 AD, we evaluated macaque dispersal pathways into and throughout China.