Can Wang
Houyuan Lu
Jianping Zhang
Limi Mao
Yong Ge- 1Department of Archaeology, School of History and Culture, Shandong University, Jinan, China
- 2Key Laboratory of Cenozoic Geology and Environment, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing, China
- 3Center for Excellence in Tibetan Plateau Earth Science, Chinese Academy of Sciences, Beijing, China
- 4College of Earth and Planetary Sciences, University of Chinese Academy of Sciences, Beijing, China
- 5Nanjing Institute of Geology and Palaeontology, Chinese Academy of Sciences, Nanjing, China
- 6Key Laboratory of Vertebrate Evolution and Human Origins, Institute of Vertebrate Paleontology and Paleoanthropology, Chinese Academy of Sciences, Beijing, China
- 7Center for Excellence in Life and Paleoenvironment, Chinese Academy of Sciences, Beijing, China
In the last decade, our understanding of rice domestication has improved by new archaeological findings using advanced analytical techniques such as morphological and morphometric analyses on rice grains, spikelet bases and phytoliths, and ancient DNA analysis on rice remains. Previous studies have considered the size of rice bulliform phytoliths as a proxy for tracking the domestication process. These phytoliths are often abundant and well preserved in sediments, and their shape is under the control of numerous genes, which may shift toward larger sizes by genetic mutation in domestication. Therefore, it has been assumed that the bulliforms of domesticated rice are usually larger than those of wild ones; however, morphometric data supporting this assumption are lacking in the literature, thereby requiring additional evidence to test its veracity. In this study, the vertical and horizonal lengths of bulliform phytoliths were measured in four rice species (domesticated Oryza sativa and wild Oryza rufipogon, Oryza officinalis, and Oryza meyeriana) from different regions of southern China. We found that the bulliform morphometric data of wild and domesticated rice overlapped and that there was no statistically significant difference between them. Therefore, bulliform size could not be used as a diagnostic indicator to distinguish domesticated rice from wild species and is a supporting rather than conclusive proxy for determining the domesticated status of rice in archaeological research. We further found that larger rice bulliform sizes likely occurred at the locations with higher temperature, precipitation, and water levels, indicating hydrothermal environment is an alternative factor influencing the size of rice bulliform phytoliths. For further archaeological use of an increasing size trend of bulliform phytoliths to reveal the process of rice domestication, we present some suggestions for controlling the influence of hydrothermal factors. Even so, the combination of bulliform phytolith size with other established criteria is strongly suggested to provide precise identification of wild and domesticated rice in future research.
Introduction
Asian rice (Oryza sativa L.) is one of the most important crops and forms a staple food for more than half of the global population (Nayar, 2014). Understanding its origins and domestication from wild rice (Oryza rufipogon Griff.) is thus an important aim for researchers. There are two major subspecies of domesticated O. sativa, Oryza sativa indica, which is thought to have originated in the Himalayan region, and Oryza sativa japonica, which is thought to have originated in China (Londo et al., 2006). Oryza rufipogon is generally recognized as the ancestor of Oryza sativa (Wei et al., 2012), with two distinct domestication events leading to the two subspecies of domesticated rice (Londo et al., 2006). However, this hypothesis is still being vigorously debated as other evidence supports a single origin of Asian rice (Huang et al., 2012). There are more than 20 wild rice species recognized in the genus Oryza (Nayar, 2014; Stein et al., 2018), which belongs to the family: Poaceae and tribe: Oryzeae. Three of these species are found in China, Oryza rufipogon, Oryza officinalis, and Oryza meyeriana (Fan et al., 2000).
In the last decade, scientific understanding of rice domestication has greatly improved by new archaeological and genetic evidence using advanced analytical techniques such as flotation combined with morphometric analysis on rice grains and spikelet bases (Liu et al., 2007; Fuller et al., 2009; Gross and Zhao, 2014; Zheng et al., 2016), phytolith analysis (Wu et al., 2014; Huan et al., 2015; Zuo et al., 2017), pan-genome analysis (Wang et al., 2018a), and genome-wide association studies (Huang et al., 2012; Civáň and Brown, 2017; Choi and Purugganan, 2018). According to recent findings, a general consensus among scholars has been reached that rice was first domesticated in the middle and lower Yangtze River regions of southern China (Molina et al., 2011; Gross and Zhao, 2014; Silva et al., 2015; Choi et al., 2017; Zuo et al., 2017), although some consider the Pearl River region to also be a part of the original area where rice domestication occurred (Huang et al., 2012; Wei et al., 2012).
The estimated time of origin of rice domestication is before 13,000 BP, based on molecular clock analysis (Molina et al., 2011; Choi et al., 2017; Choi and Purugganan, 2018), which is much older than the earliest archaeological date of domestication (<10,000 BP) (Fuller et al., 2014; Larson et al., 2014). This disagreement may result from genetic studies that mainly focused on identifying the origins of alleles associated with domestication (e.g., sh4, rc, laba1, prog1), which likely emerged in wild rice prior to domestication (Choi et al., 2017; Civáň and Brown, 2017). On the other hand, archaeological studies attempted to detect the first appearance of morphological traits associated with domestication in archaeobotanical remains (Fuller et al., 2009; Jones and Liu, 2009). The molecular and archaeological chronologies may date two main phases in the macroevolutionary process: the emergence of a trait and the success of that trait (the trait becomes quantitatively significant within a population) (Katz, 2019). However, the most notable chronological dispute over the rice domestication is between two archaeological opinions: one suggests that the process of rice domestication may have begun around 10,000–9,000 BP (Liu et al., 2007; Zheng et al., 2007; Wu et al., 2014; Zheng et al., 2016; Zuo et al., 2017), while the other suggests that domestication of rice did not occur until around 8,000–6,000 BP (Fuller et al., 2007; Fuller et al., 2008; Fuller et al., 2009; Fuller et al., 2010; Larson et al., 2014). This debate is largely attributable to the differences in the methods employed, and the criteria used by various authors to identify domestication in rice remain from early sites in the Yangtze River region such as Shangshan, Kuahuqiao, and Xiaohuangshan. Establishing accurate and feasible criteria for distinguishing between domesticated and wild rice is thus of prime importance.
Three important lines of archaeological evidence have frequently been used in China, including grain size and morphological characteristics, spikelet bases, and phytoliths (Liu and Chen, 2012; Fuller and Castillo, 2014; Fuller, 2018). Grain morphometrics are considered to be semi-domestication traits and may not be diagnostic indicators of early domestication (Fuller, 2007; Fuller et al., 2007), mostly due to the considerable variation and overlap in length between domesticated and wild populations which leads to some proportion of false assignments in ancient rice grains (Fuller, 2007). Many studies have focused on the form of spikelet base, which is thought to be the most diagnostic trait in rice remains in terms of identifying domestication status (Fuller, 2018). However, these are insufficient, because the distinctive characteristics of immaturity, shattering, and nonshattering states and/or wild, japonica, and indica rice based on spikelet bases are divergent among the criteria provided by different researchers (Zheng et al., 2007; Pan, 2008; Fuller et al., 2009; Fuller et al., 2010; Pan, 2011; Gross and Zhao, 2014; Zheng et al., 2016), and their diagnostic power for domestication has the potential to be more reliable. More importantly, these macrobotanical remains do not preserve well in early sediments with acid soil; therefore, very few have been recovered from sites dated earlier than 9,000 BP (Zhao, 2011; Qin, 2012; Zhao and Jiang, 2016).
Phytoliths have played an important role in the identification of rice remains recovered from early archaeological sites, due to their high resistance to decomposition (Piperno, 2006; Ball et al., 2016). Double-peaked phytoliths from husks and bulliform phytoliths from rice leaves are both certainly diagnostic indicators of Oryza and show variation within and between species. A number of identification criteria based on these phytoliths have been suggested and widely used in the last 20 years (Zhao et al., 1998; Lu et al., 2002; Gu et al., 2013; Wu et al., 2014; Huan et al., 2015; Hilbert et al., 2017; Zuo et al., 2017; Wang et al., 2018b). Although the utility of these methods for distinguishing domesticated from wild rice is still under debate (Fuller et al., 2010; Fuller, 2018), they are recognized as key alternative methodologies besides using morphological domestication data of rice macroremains. Three-dimensional measurements and discriminant function analysis of double-peaked phytoliths are useful in determining the wild/domesticated nature of rice remains; however, double-peaked phytoliths usually present their side and top view under the microscope, which does not meet the requirements of morphometric analysis (present in front view) (Zhao et al., 1998; Gu, 2009), making the work arduous in most cases.
Rice bulliform phytoliths are abundant in rice leaves and are often well preserved and represented in archaeological sediments (Fujiwara, 1993; Wang and Lu, 1993). Generally, bulliform phytoliths in Oryza have a distinctive fan shape with numerous scale-like decorations on the half round side (lateral side) (Lu et al., 2002; Wang and Lu, 2012; Gu et al., 2013). Morphological measurements and number of scale-like decorations along the scalloped edge have been employed to distinguish wild Oryza species from domesticated ones. Studies on modern rice plants and paddy surface soils have suggested that bulliform phytoliths with ≥9 scale-like decorations were likely domesticated, while those with <9 were generally wild (Lu et al., 2002; Huan et al., 2015). Whether this feature is a useful domestication indicator remains inconclusive and requires further validation; in addition, genetic explanatory mechanisms of bulliform scale-like decoration variation between species remain unclear and deserve further study.
Bulliform shape of rice appears to be under the control of 16 genes (QTLs) (Zheng et al., 2003a) and phytoliths may shift toward larger sizes as a result of genetic mutation during the domestication process (Zheng et al., 2003b; Piperno, 2006; Luo et al., 2016). Some researchers have, therefore, assumed that the bulliforms of domesticated rice are usually larger than wild ones and that the trend of increasing size in rice bulliform phytoliths could reflect domestication of rice (Zheng et al., 2003b, Zheng et al., 2004; Fuller et al., 2007). In recent years, vertical and horizontal lengths (i.e., sizes) of rice bulliforms have increasingly been used as a proxy for tracking the domestication process or determining the degree of domestication at different sites, such as the Tanghu (Zhang et al., 2012), Zhuzhai (Wang et al., 2018b), Shunshanji (Luo et al., 2016), Shangshan, Hehuanshan, and Huxi sites (Zuo et al., 2017; Qiu et al., 2019). However, to date, morphometric data from modern rice plants supporting this method and their assumptions are missing (Pearsall et al., 1995; Zhang and Wang, 1998; Ma and Fang, 2007; Gu et al., 2013), and thus additional evidence is required to test its veracity.
Moreover, some studies have argued that changes in phytolith size may not only be triggered by domestication but also influenced by environmental factors, such as CO2 concentrations (Ge et al., 2010), evapotranspiration rates (Issaharou-Matchi et al., 2016), and water levels in the growing habitat (Fuller, 2018). Fuller (2018) indicated that the 16 genes suggested by Zheng et al. (2003a) only explained between 37 and 54% of bulliform variation, suggesting that the environment or growing conditions also play an essential role. These authors further speculated that if shifting bulliform morphology was merely a phenotypic response to environmental conditions, it would be a less useful indicator of domestication. Thus, without data supporting the exclusion of environmental factors, bulliform phytolith measurements alone may not be an accurate identification tool for distinguishing between domesticated and wild Oryza species.
In the present study, we tested whether the size of bulliform phytoliths is an effective statistical indicator for distinguishing between wild rice and domesticated rice, based on the comparative analysis of morphometric data from different rice species. In addition, we attempted to examine how growing conditions, especially climate and water levels, influence bulliform size. Finally, we discussed how rice bulliform phytolith morphometry can be used in archaeological research.
Materials and Methods
In the present study, a total of 24 specimens of Oryza were sampled. The samples consisted of six specimens of the domestic O. sativa, and for the wild species, 16 specimens of O. rufipogon, 1 specimen of O. officinalis, and 1 specimen of O. meyeriana all collected in southern China (Table 1; Figure 1). All 6 specimens of O. sativa and 8 of the 16 specimens of O. rufipogon were sampled from the test paddy field belonging to Wuhan Botanical Garden, Chinese Academy of Sciences (CAS), at Huazhong Agricultural University, Wuhan, Hubei Province. Another eight specimens of O. rufipogon were sampled from Hainan, Yunnan, Hunan, and Jiangxi Provinces. The specimens of O. officinalis and O. meyeriana were sampled from Hainan Province.
Table 1 Information on the rice plants studied and measured data of vertical length and horizontal length of bulliform phytoliths from the studied samples.
Figure 1 Locations of sample collection sites. Test paddy field in Huazhong Agricultural University, Wuhan, Hubei Province (1); Shuitaoshuxia, Dongxiang, Jiangxi Province (2); Anjiashan, Dongxiang, Jiangxi Province (3); Huli Marsh, Chaling, Hunan Province (4); Yuanjiang protection area of wild rice, Yunnan Province (5); Xishuangbanna Tropical Botanical Garden, CAS, Yunnan Province (6); Hulu village, Wenchang, Hainan Province (7); Tanshen village, Wenchang, Hainan Province (8); Mingxing village, Wanning, Hainan Province (9); Zhangxian village, Lingshui, Hainan Province (10); Luhuitou Park, Sanya, Hainan Province (11).
Field collection of rice plants was assisted by a team of investigators from Wuhan Botanical Garden and Nanjing Institute of Geology and Palaeontology, CAS, with the permission of the owner or regulatory body for each location. The criteria for categorizing and naming the species collected were accepted from the classification scheme of the genus Oryza in the Flora of China (Liu and Phillips, 2006) (available at http://flora.huh.harvard.edu/china/PDF/PDF22/Oryza.pdf). The specimens of O. sativa and O. rufipogon in Wuhan were both mature when we collected them between September 21 and 23, 2011, but the ripening rate of O. rufipogon was very low due to low temperatures. In Hainan Province, the specimens of O. rufipogon were in anthesis and immature while the specimens of O. officinalis and O. meyeriana were mature and had begun shattering when we collected them between December 1 and 7, 2012. The specimens of O. rufipogon in Chaling, Hunan Province, were immature when we collected them between September 18 and 28, 2010. The specimens of O. rufipogon in Jiangxi and Yunnan Province were already mature when we collected them during October and November 2010, respectively. All plant samples were preserved at the Institute of Geology and Geophysics, CAS, Beijing.
There are different hydrological environments among our sampling sites. In Wuhan, specimens of O. sativa and O. rufipogon were simultaneously cultivated at paddy fields with shallow water (Figures 2a, b), which sometimes needs draining to achieve moderate draught during the pustulation and fruiting stage; O. rufipogon from Xishuangbanna Tropical Botanical Garden, CAS, Jinghong County, had a similar habitat. In Chaling, O. rufipogon grows in the Huli Marsh where there is perennial stagnant water (Figure 2c); O. rufipogon from Yuanjiang has a similar habitat. In Anjiashan and Shuitaoshuxia, Dongxiang County, O. rufipogon grows in a seasonal wetland with 0–100-cm water depth. In Hainan Province, specimens of O. rufipogon from Wenchang and Wanning Cities grow in permanent wetlands where there are ponds filled with deep water with 30–150-cm depth (Figures 2d–g, j); O. officinalis in Lingshui County grows in a ravine stream in the valley and prefers shady and wet habitats (Figures 2h, k); O. meyeriana in Sanya City grows in an understory bush on a hill and prefers shady and dry habitats (Figures 2i, l).
Figure 2 Photos of parts of sampling locations and rice plants. Domesticated rice paddy in Huazhong Agricultural University, Wuhan, Hubei Province (site 1 in Figure 1) (a); wild rice paddy in Huazhong Agricultural University, Wuhan, Hubei Province (site 1 in Figure 1) (b); site of Oryza rufipogon population in Huli Marsh, Chaling, Hunan Province (site 4 in Figure 1) (c); site of O. rufipogon population in Hulu village, Wenchang, Hainan Province (site 7 in Figure 1) (d); site and plants of O. rufipogon population in Mingxing village, Wanning, Hainan Province (site 9 in Figure 1) (e, f); site of O. rufipogon population in Tanshen village, Wenchang, Hainan Province (site 8 in Figure 1) (g); site of O. officinalis population in Zhangxian village, Lingshui, Hainan Province (site 10 in Figure 1) (h); site of O. meyeriana population in Luhuitou Park, Sanya, Hainan Province (site 11 in Figure 1) (i); plant of O. rufipogon in Hainan Province (j); plant of O. officinalis in Hainan Province (k); plant of O. meyeriana in Hainan Province (l); the pictures of rice plants were taken by Dr. Limi Mao. The individual in Figure 2c was Dr. Jianping Zhang who had approved the publication of this image.
For each rice specimen, we selected all leaf blades from the bottom to the top of a single plant, making sure to keep the leaf blade intact for phytolith extraction. This is because there is significant difference in bulliform phytolith size among different leaf blades of the same plant and different parts of the same leaf blade, with the smaller bulliform phytoliths from the lower leaves (Fuller and Qin, 2009). Bulliform size tends to decrease from leaf base to leaf apex of the same leaf blade (Wang et al., 1997). Therefore, variation in bulliform phytoliths from a few randomly selected rice leaves does not reflect the overall data, and only the selection of all intact leaves can guarantee the representativeness and reliability of the data.
Leaf blades were cleaned with distilled water in an ultrasonic cleaner, and then prepared for wet oxidation: 1) all samples were cut into 1–3-cm pieces and placed in 20 ml of 65% saturated nitric acid for over 12 h, then heated in a water bath for 20 min to oxidize organic materials completely. 2) The solutions were centrifuged at 3,000 rpm for 6 min, decanted and rinsed three times with distilled water, and then rinsed with 95% ethanol until the supernatants were clear. 3) The extracted phytoliths were mounted onto microscopic slides in neutral resin, and the residual samples were transferred to storage vials.
A Leica DM750 light microscope at 600× magnification was used for photomicrography and phytolith counting. One hundred or more bulliform phytoliths, including asymmetric types, were counted in each sample, except for the sample of O. meyeriana which produced only a few bulliform phytoliths. Two morphometric parameters, vertical and horizonal lengths (VL and HL), were measured to describe the size of bulliform phytoliths. The measurements were taken from images using the ImageJ software (version 1.48r.). Descriptive statistics of morphometric data were performed using Excel software, and the mean ± SD of each sample was plotted using the Grapher software to perform a comparative analysis. Discriminant function analysis in SPSS 24.0 software was then used to statistically determine the differences in bulliform sizes in different species.
In order to understand whether bulliform phytolith size correlated with environmental conditions and to what degree, a Pearson correlation analysis was performed. Eight environmental variables were chosen: altitude, mean annual precipitation (MAP), July precipitation (MP7), January precipitation (MP1), mean annual temperature (MAT), July temperature (MT7), January temperature (MT1), and relative humidity (HHH). Modern climatic data for the 11 sampling sites (Figure 1) were obtained from the nearest meteorological station to each site, since the spatial variation of climatic parameters exhibits a clear gradient across these locations. These data can be collected from the databases (1981–2010) of the National Meteorological Information Center, China (http://data.cma.cn/). Origin 8.5.1 software was used to conduct the correlation analysis of bulliform morphometrics and environmental variables, which were plotted into scatter plots. A linear regression was inserted into these scatter plots, and then the Pearson correlation coefficients (r) and significances (P) were taken to statistically evaluate the correlation.
Results
Morphological Contrast of Bulliform Phytoliths in the Four Oryza Species
Overall, the bulliform phytoliths in O. sativa and O. rufipogon have similar shapes with an intact circular part of the fan, round arc of the scalloped edge, and ridge-like tubercle on the lateral side (Figure 3). Significant intraspecific morphological variation, however, was also found in both O. sativa (Figures 3a–h) and O. rufipogon (Figures 3i–p). Moreover, we noted that the bulliform phytoliths from the Xishuangbanna O. rufipogon had a very specific shape with a small circular part of the fan, angular scalloped edge, large deep decorations, and without the ridge-like tubercle on the lateral side (Figures 3q–t). This shape was not only different from that of O. sativa but also distinct among other O. rufipogon specimens.
Figure 3 Bulliform phytoliths in some studied Oryza sativa and O. rufipogon species. Bulliform phytoliths from O. sativa (AA31) (a–d); bulliform phytoliths from O. sativa (AA48) (e–h); bulliform phytoliths from O. rufipogon (HL-3) (i–l); bulliform phytoliths from O. rufipogon (Z13) (m–p); bulliform phytoliths from O. rufipogon (XSBN) (q–t); scale bar = 30 μm. The red arrows point to the ridge-like tubercle.
Bulliform phytoliths of O. officinalis are mostly large and full in shape with a rounded arc of the scalloped edge, irregular large deep decorations, longer handles, and a shorter intact circular part of the fan but without the ridge-like tubercle on the lateral side (Figure S1). From an overall perspective, this shape is similar to that of O. sativa and O. rufipogon.
Bulliform phytoliths of O. meyeriana are generally long and thin and not full in shape like the other varieties. These phytoliths are very small and appear similar to a teardrop or nail, with an angular scalloped edge, small irregular, but deep decorations, a longer handle, and shorter circular part of the fan not intact, without the ridge-like tubercle on the lateral side (Figure S2). This shape is different from that of the other three Oryza species.
Morphometric Analysis of Bulliform Phytoliths in the Four Oryza Species
Overall, our morphometric data demonstrated that the bulliforms of domesticated rice were not always larger than wild ones, and there was no significant difference in size. Table 1 shows the mean values of vertical length (VL) and horizontal length (HL) of bulliform phytoliths from the studied samples. For original measured data, see the Datasheet S1.
In the six O. sativa specimens, the maximum mean of VL of bulliform phytoliths was 44.52 ± 6.22 μm (AA37), while the minimum was 37.77 ± 4.74 μm (AA38); the maximum mean of HL of bulliform phytoliths was 38.31 ± 6.26 μm (AA31), while the minimum was 32.78 ± 4.47 μm (AA38). In all 609 bulliform phytoliths from O. sativa, the maximum VL was 60.27 μm occurring in sample AA36, while the minimum was 25.25 μm occurring in sample AA38; the maximum HL was 61.96 μm occurring in sample AA31, while the minimum was 20.94 μm occurring in sample AA48.
In the 16 O. rufipogon specimens, the maximum mean of VL of bulliform phytoliths was 52.23 ± 6.84 μm (WN-7), while the minimum was 29.89 ± 4.01 μm (Z13); the maximum mean of HL of bulliform phytoliths was 43.05 ± 6.74 μm (WN-7), while the minimum was 24.19 ± 3.99 μm (Z13). In all 1,649 bulliform phytoliths, the maximum VL was 70.11 μm occurring in sample WN-7, while the minimum was 21.30 μm occurring in sample Z13; the maximum HL was 59.81 μm occurring in sample WN-7, while the minimum was 15.54 μm occurring in sample Z13.
In the only specimen of O. officinalis, the mean values of VL and HL of bulliform phytoliths were 51.13 ± 8.77 and 40.43 ± 7.34 μm, respectively. In all 100 bulliform phytoliths, the maximum VL was 80.00 μm, while the minimum was 32.94 μm; the maximum HL was 62.76 μm, while the minimum was 25.17 μm.
In the only specimen of O. meyeriana, the mean values of VL and HL of bulliform phytoliths were 38.02 ± 5.35 and 29.37 ± 6.29 μm, respectively. In all 34 bulliform phytoliths, the maximum VL was 50.59 μm, while the minimum was 30.49 μm; the maximum HL was 44.26 μm, while the minimum was 18.18 μm.
Figure 4 presents a comparison of the sizes of bulliform phytoliths from the studied rice species. The results revealed that the values of VL and HL from O. rufipogon were scattered and widely overlapped with the other three Oryza species. Although the bulliform sizes of O. sativa were larger compared with O. meyeriana, they partly overlapped with O. rufipogon (mean VL: 37–45 μm; HL: 32–40 μm) and were significantly smaller than those of O. officinalis. The bulliform size of O. meyeriana was smaller compared with other Oryza species and only larger than a few specimens of O. rufipogon. The bulliform size of the specimen of O. officinalis was larger, exceeding most of the studied samples.
Figure 4 Oryza bulliform phytolith measurements from the studied species. VL, vertical length; HL, horizontal length of rice bulliform; gray error bar represents ± 1 SD.
Parameters VL and HL were important and used in the discriminant function analysis (Table S1). The following two canonical discriminate functions were used in the analysis: function 1 explained 89.5% of the variance, and function 2 explained 10.5% of the variance. Parameter VL had the largest absolute correlation with function 1, indicating that it contributed most to function 1; parameters HL had the largest absolute correlation with function 2, indicating that it contributed more to function 2 (Table S1). These two functions were used to plot the data (Figure 5). Four groups without distinct centroids were obtained; O. sativa had a clear intersection with O. rufipogon, whereas there was a slight distinction for O. officinalis and O. meyeriana with regard to the other species. The accuracy of the classification was ascertained by cross validating the results (Table 2). Only 36.8% of the original data and 36.7% of the cross-validated data were correctly classified, suggesting that the discriminant functions obtained using parameters VL and HL could not be successfully used to discriminate O. sativa, O. rufipogon, O. meyeriana, and O. officinalis. Thus, it supported the conclusion that there were no significant differences in bulliform sizes between wild and domesticated rice.
Figure 5 Discriminant function analyses of Oryza sativa, O. rufipogon, O. meyeriana, and O. officinalis.
Table 2 Classification results of the discriminant function analysis.
The 24 specimens of rice were then divided into two groups: mature and immature (Figure 6). As can be seen, the sizes of bulliform phytoliths from mature rice were scattered and partly overlapped with the immature rice. In contrast, bulliform phytoliths of immature rice were slightly larger than most mature species. There was no significant difference in bulliform size between mature and immature rice species.
Figure 6 Contrast of Oryza bulliform phytolith sizes from the studied specimens depending on whether they were mature or immature. The sample codes in the figure refer to the field numbers: specimens WN-7, TS-1, and HL-3 were O. rufipogon from Wanning and Wenchang, Hainan Province, and specimen P4 was O. rufipogon from Chaling, Hunan Province.
We further compared the bulliform size of rice species in terms of their growing region (Figure 7). It was found that the bulliforms of O. rufipogon and O. officinalis growing in the tropical Hainan were the largest with mean vertical and horizontal lengths greater than 45 and 40 μm, respectively. The bulliforms of O. rufipogon growing in Chaling, Dongxiang, and Yuanjiang and O. sativa growing in Wuhan were similar in size, with mean VL and HL ranges of 40–45 and 33–40 μm, respectively. O. rufipogon growing in Wuhan and Xishuangbanna and O. meyeriana growing in Hainan, have the smallest bulliform phytoliths with mean VL and HL ranges of 33–39 and 29–34 μm, respectively. It is noted that O. rufipogon and O. officinalis in the Hainan population have the most favorable habitat in terms of water availability, with permanent deep water or perennial ravine streams; O. rufipogon in Chaling, Dongxiang, and Yuanjiang grow in marshes and seasonal wetlands where stagnant water converges at the roots of rice; O. rufipogon in Wuhan and Xishuangbanna grow in paddy fields with shallow water, which is occasionally drained to maintain relatively dry habitats; O. meyeriana in Hainan prefers a drier environment.
Figure 7 Contrast of Oryza bulliform phytolith sizes from the studied specimens in terms of their growing regions. WH, Wuhan; HN, Hainan; YN, Yunnan; CL, Chaling; DX, Dongxiang; XSBN, O. rufipogon from Xishuangbanna; YNY-1, O. rufipogon from Yuanjiang.
In addition, the effect of habitat wetness on bulliform phytolith size was investigated. We found that the bulliform phytolith size of specimens of O. rufipogon native to the warmer and wetter sites of Wenchang, Hainan (field no. Z23), and Dongxiang, Jiangxi (field no. Z1), and cultivated in the paddy fields in Wuhan, was smaller than that of native species (Figure 8). Similarly, in the Wuhan paddy field, the species of O. rufipogon native to the warmer and wetter sites of Guangdong (field nos. Z13, Z15), Guangxi (field nos. Z12, Z14, Z27), and India (field no. Z41) also have smaller bulliforms (Figure 8). In the tropical Hainan, under parallel climate conditions, O. officinalis in the aquatic environment has a significantly larger bulliform size than those of O. meyeriana in a dry habitat (Figure 8). Therefore, the sizes of bulliform phytoliths from wild rice with preferable water habitats were mostly larger than those of wild rice under relatively dry conditions.
Figure 8 Contrast of bulliform phytolith sizes from some studied wild rice in terms of water levels in their habitats. Red represents wet habitats; blue represents relatively dry habitats. The sample codes in the figure refer to the field numbers. WH, Wuhan; HN, Hainan; DX, Dongxiang; GX, Guangxi; GD, Guangdong.
Finally, Figure 9 shows the comparison of bulliform phytolith sizes between O. sativa and O. rufipogon growing in the adjacent test paddy field in Wuhan. For O. sativa, the ranges of mean VL and HL of bulliforms fell into 37–45 and 32–48 μm, respectively. For O. rufipogon, the mean VL and HL ranges were 29–41 and 24–34 μm, respectively. Thus, our data also indicated that the bulliform phytoliths of O. rufipogon may be generally smaller than those of O. sativa if they were artificially grown in the same environment.
Figure 9 Contrast of bulliform phytolith sizes from domesticated and wild rice growing at the adjacent test paddy field in Wuhan. Demonstrating that the bulliform phytoliths of these domesticated rice specimens were generally larger compared with wild specimens (Oryza rufipogon) in the same environment. WH, Wuhan.
Correlation Analysis of Bulliform Morphometrics and Environmental Variables
Summary statistics for the eight environmental variables of different sites are given in Table 3. The VL and HL values of bulliform phytoliths plotted against the different environmental variables and the results are shown in Figure 10 and Figure 11. MAP, MAT, MT1, and HHH had a positive correlation (r = 0.438–0.610) with VL and HL, and the linear regression analysis revealed that these correlations were significant (p < 0.01 or 0.05) (Figure 10), demonstrating that the bulliform phytolith sizes are affected by changes in these climatic parameters. In contrast, other environmental characteristics, including MP7, MT7, MP1, and altitude, were not significantly correlated with changes in bulliform sizes (Figure 11).
Table 3 Summary of the environmental variables in the 11 sampling sites used for correlation analysis.
Figure 10 Scatter plots of bulliform phytolith sizes as defined by VL and HL for studied species versus different climatic variables. Mean VL values vs. observed MAP values (A); mean HL values vs. observed MAP values (B); mean VL values vs. observed MAT values (C); mean HL values vs. observed MAT values (D); mean VL values vs. observed MT1 values (E); mean HL values vs. observed MT1 values (F); mean VL values vs. observed HHH values (G); mean HL values vs. observed HHH values (H). Linear regression analysis of these data indicates significant correlation.
Figure 11 Scatter plots of bulliform phytolith sizes as defined by VL and HL for studied species versus different environmental variables. Mean VL values vs. observed MP7 values (A); mean HL values vs. observed MP7 values (B); mean VL values vs. observed MT7 values (C); mean HL values vs. observed MT7 values (D); mean VL values vs. observed MP1 values (E); mean HL values vs. observed MP1 values (F); mean VL values vs. altitude values (G); mean HL values vs. altitude values (H). Linear regression analysis of these data indicates no significant correlation.
Discussion and Conclusions
Causes of Variations in Rice Bulliform Phytolith Morphometry
The results of the present study indicate that morphometric measurements of bulliform phytoliths from wild and domesticated rice widely overlap (Figure 4), exhibiting little diagnostic potential for taxonomic identification at the species level. These results thus support the conclusion that the morphometry of bulliform phytoliths is not as informative for distinguishing between domesticated rice and wild rice, as suggested by previous studies (e.g., Pearsall et al., 1995; Wang and Lu, 2012; Gu et al., 2013). However, the above morphometric data have often been overlooked, and an increasing number of studies have recently used bulliform phytolith size as a proxy to track the rice domestication process (Zhang et al., 2012; Luo et al., 2016; Zuo et al., 2017; Qiu et al., 2019). The bulliform phytoliths from domesticated rice really were larger than those from wild ones in the same test paddy field in Wuhan city (Figure 9), possibly indicating that domestication may result in an increase in bulliform phytolith size. The genetic and phylogenetic signal for this bulliform size variation has not been well revealed to date. Furthermore, this increase did not necessarily result from domestication and may be caused by other factors.
Previous studies have suggested that two factors, plant maturity and environmental conditions, may affect rice bulliform phytolith size. It is suggested that the bulliform phytoliths from mature rice leaves are usually larger than those from immature leaves (Zheng et al., 2003a; Qin et al., 2006; Fuller et al., 2007). This hypothesis may be real for the plants growing in the same location, but is not supported when comparing bulliform size of rice species from different sites (Figure 6). Because environmental conditions are seemingly much more important, and the effect of degree of maturity could be ignored when multiple environmental factors are considered. According to the results of the present study (Figures 10 and 11), from a statistical point of view, we can conclude that the larger rice bulliform phytolith sizes, as defined by the higher VL and HL values, likely occurred at the locations with higher temperature and precipitation. Therefore, the increasing trend in rice bulliform phytolith size in some archaeological records (Zheng et al., 2003b; Zheng et al., 2004; Luo et al., 2016; Zuo et al., 2017; Qiu et al., 2019) may also be caused by climatic changes during the early and middle Holocene when temperature and precipitation were gradually rising.
The present study also revealed that the growing microenvironment, such as water environment, can also influence the size of rice bulliform phytoliths. Rice growing under wetter conditions usually produced larger bulliform phytoliths than those growing under drier conditions (Figures 7 and 8). Therefore, except climate regimes, changes in wet/dry habitat for rice should be considered for the use of bulliform size to track the process of rice domestication.
It should be pointed out that the present study just revealed the hydrothermal condition as one of the environmental factors influencing rice bulliform size. Some factors such as plant genotypes, soil fertility, light period length on photosynthesis, and other abiotic factors may also cause these variations, which were not controlled or excluded for this study. Conditional plantation experiment under controlling environments and genotypes in test paddy field is needed to further test if and how water levels and temperature can affect the bulliform phytolith size.
Implications for Archaeology of Rice Domestication
Changes in bulliform phytolith size of rice are regulated not only by domestication, which possibly represents genetic changes, but also by environmental factors. Given that environmental factors influence bulliform phytolith size of rice and that the role of genetic background has not yet been firmly established similar to established domestication traits such as non-shattering and increased seed size, bulliform measurement was considered as a semi-domestication trait (Fuller and Qin, 2009). Therefore, the use of rice bulliform phytolith size as an index for determining domesticated plants from their wild ancestors should be conditional. In other words, if the increasing size trend of bulliform phytoliths is used to reveal the process of rice domestication, the influence of hydrothermal conditions should be excluded first.
For further archaeological use of this index, we suggest that: 1) if the time series of rice bulliform phytolith size from a region is long, then the climatic changes (fluctuations in temperature and precipitation) through time should be considered, and the results from quantitative reconstructions of paleoclimate could be used as an independent variable to explain bulliform size variation; 2) the spatial scale of studied regions should be small and without a clear climatic gradient; 3) parallel comparison of rice domestication processes in different regions using bulliform size should consider climatic differences between the regions; and 4) the changes in rice arable systems (wet/dry growing conditions) in any studied archaeological sites should be first revealed using the promising sensitive/fixed phytolith morphotype model defined by Weisskopf et al. (2015).
Even though the influence of environmental factors has been controlled or excluded, and rice bulliform phytoliths shifting toward larger sizes are interpreted as reflecting the domestication process, it is still not possible to provide a determinate range of bulliform size for identifying domesticated rice, due to the wide overlap observed in the bulliform morphometric data between modern wild and domesticated rice (Figure 4). Thus, rice bulliform phytolith size is a supporting rather than conclusive proxy for determining the domesticated status of rice in archaeological research. Combination of bulliform phytolith size with other established criteria can provide precise identification of wild and domesticated rice.
Finally, notably, frequent gene exchange occurs between domesticated and wild rice, and there is a co-evolutionary relationship between them (Song et al., 2006; Zhao et al., 2010; Ge and Sang, 2011; Choi et al., 2017). Recent large-scale genomic analysis showed that O. rufipogon populations are widely affected by the gene flow of domesticated rice, and that the existing O. rufipogon are actually a hybrid swarm (Wang et al., 2017). This indicates that it is difficult to rule out the interference of the domesticated rice gene flow when using the existing O. rufipogon species as a reference for phytolith morphological analysis. Wild rice populations in the region with higher rice farming intensity are more affected by the introgression of the domesticated rice gene, and the genetic relationship with domesticated rice is closer (Song et al., 2003; Song et al., 2006), leading to the possibility of bias in the morphometric measurements of their bulliform phytoliths. Further research utilizing archaeological rice remains combined with ancient DNA analysis (e.g., Tanaka et al., 2010; Castillo et al., 2016), may reduce the interference of domesticated rice gene flow, and thus generate credible results to establish suitable criteria for distinguishing between domesticated rice and wild rice.
Data Availability
The datasets generated during the current study are available from the corresponding author on reasonable request.
Author Contributions
HL and CW designed research. JZ, CW, and LM collected samples. CW and JZ performed experiments. CW and YG analyzed data. CW and HL wrote the paper. All authors read and approved the final manuscript.
Funding
This work is supported by the National Natural Science Foundation of China (Grant Nos. 41701233, 41830322 and 41430103), the Strategic Priority Research Program of Chinese Academy of Sciences (Grant No. XDB26000000), and the Youth Top Talent Program of the Education Department of Hebei Province (Grant No. BJ2018118).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We are grateful to Dr. Wen Zhou and Yumei Li for their assistance in the field work on rice sampling; Dr. Yajie Dong and Deke Xu for kind help with data processing. We appreciate the constructive comments from the editor and three reviewers. We also would like to thank Editage (www.editage.cn) for English language editing.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2019.01037/full#supplementary-material
Datasheet S1 | Morphometric data of bulliform phytoliths from the 24 studied samples.
References
Ball, T., Chandler-Ezell, K., Dickau, R., Duncan, N., Hart, T. C., Iriarte, J., et al. (2016). Phytoliths as a tool for investigations of agricultural origins and dispersals around the world. J. Archaeol. Sci. 68, 32–45. doi: 10.1016/j.jas.2015.08.010
Castillo, C. C., Tanaka, K., Sato, Y.-I., Ishikawa, R., Bellina, B., Higham, C., et al. (2016). Archaeogenetic study of prehistoric rice remains from Thailand and India: evidence of early japonica in South and Southeast Asia. Archaeol. Anthropol. Sci. 8 (3), 523–543. doi: 10.1007/s12520-015-0236-5
Choi, J. Y., Platts, A. E., Fuller, D. Q., Hsing, Y. I., Wing, R. A., Purugganan, M. D. (2017). The rice paradox: multiple origins but single domestication in Asian rice. Mol. Biol. Evol. 34 (4), 969–979. doi: 10.1093/molbev/msx049
Choi, J. Y., Purugganan, M. D. (2018). Multiple origin but single domestication led to Oryza sativa. G3 (Bethesda). 8 (3), 797–803. doi: 10.1534/g3.117.300334
Civáň, P., Brown, T. A. (2017). Origin of rice (Oryza sativa L.) domestication genes. Genet. Resour. Crop. Evol. 64 (6), 1125–1132. doi: 10.1007/s10722-017-0518-0
Fan, S. G., Zhang, Z. J., Liu, L., Liu, H. X., Liang, C. Y. (2000). The species, geographical distribution of wild rice and their characteristics in China. J. Wuhan Bot. Res. 18 (5), 417–425 (in Chinese).
Fujiwara, H. (1993). “Research into history of rice cultivation using plant opal analysis,” in Current Research in Phytolith Analysis: Application in Achaeology and Paleoecology. Eds. Pearsall, D. M., Piperno, D. R. (Philadelphia: University of Pennsylvania), 147–158.
Fuller, D. Q. (2007). Contrasting patterns in crop domestication and domestication rates: recent archaeobotanical insights from the Old World. Ann. Bot. 100 (5), 903–924. doi: 10.1093/aob/mcm048
Fuller, D. Q. (2018). Rice: a user guide for archaeologists (version 1.0). Available at: https://www.ucl.ac.uk/rice/resources.
Fuller, D. Q., Castillo, C. (2014). “Rice: origins and development,” in Encyclopedia of Global Archaeology, vol. 6339–6343. Ed. Smith, C. (Berlin: Springer). doi: 10.1007/978-1-4419-0465-2_2185
Fuller, D. Q., Denham, T., Arroyo-Kalin, M., Lucas, L., Stevens, C. J., Qin, L., et al. (2014). Convergent evolution and parallelism in plant domestication revealed by an expanding archaeological record. Proc. Natl. Acad. Sci. U.S.A. 111 (17), 6147–6152. doi: 10.1073/pnas.1308937110
Fuller, D. Q., Harvey, E., Qin, L. (2007). Presumed domestication? Evidence for wild rice cultivation and domestication in the fifth millennium BC of the Lower Yangtze region. Antiquity 81 (312), 316–331. doi: 10.1017/S0003598X0009520X
Fuller, D. Q., Qin, L. (2009). Archaeobotany in the origin of the rice agriculture. Relics From South 3, 38–45 (in Chinese). doi: 10.3969/j.issn.1004-6275.2009.03.006
Fuller, D. Q., Qin, L., Harvey, E. (2008). Rice archaeobotany revisited: comments on Liu et al. (2007). Antiquity 82 (315). Project Gallery, online (www.antiquity.ac.uk/ProjGall/fuller1/).
Fuller, D. Q., Qin, L., Zheng, Y., Zhao, Z., Chen, X., Hosoya, L. A., et al. (2009). The domestication process and domestication rate in rice: spikelet bases from the Lower Yangtze. Science 323 (5921), 1607–1610. doi: 10.1126/science.1166605
Fuller, D. Q., Sato, Y. I., Castillo, C., Qin, L., Weisskopf, A. R., Kingwell-Banham, E. J., et al. (2010). Consilience of genetics and archaeobotany in the entangled history of rice. Archaeol. Anthropol. Sci. 2 (2), 115–131. doi: 10.1007/s12520-010-0035-y
Ge, S., Sang, T. (2011). Inappropriate model rejects independent domestications of indica and japonica rice. Proc. Natl. Acad. Sci. U.S.A. 108 (39), E755–E755. doi: 10.1073/pnas.1111601108
Ge, Y., Jie, D. M., Guo, J. X., Liu, H. M., Shi, L. X. (2010). Response of phytoliths in Leymus chinensis to the simulation of elevated global CO2 concentrations in Songnen Grassland, China. Chin. Sci. Bull. 55 (32), 3703–3708. doi: 10.1007/s11434-010-4123-2
Gross, B. L., Zhao, Z. J. (2014). Archaeological and genetic insights into the origins of domesticated rice. Proc. Natl. Acad. Sci. U.S.A. 111 (17), 6190–6197. doi: 10.1073/pnas.1308942110
Gu, H. B. (2009). “An overview of the methods distinguishing the rice phytolith between Oryza sativa subsp. Hsien and Oryza sativa subsp. Keng. Journal of Human Archaeology, The Institute of Archaeology of Hunan Province (Changsha: Yuelu Publishing House), 268–276 (in Chinese).
Gu, Y. S., Zhao, Z. J., Pearsall, D. M. (2013). Phytolith morphology research on wild and domesticated rice species in East Asia. Quat. Int. 287, 141–148. doi: 10.1016/j.quaint.2012.02.013
Hilbert, L., Neves, E. G., Pugliese, F., Whitney, B. S., Shock, M., Veasey, E., et al. (2017). Evidence for mid-Holocene rice domestication in the Americas. Nat. Ecol. Evol. 1 (11), 1693–1698. doi: 10.1038/s41559-017-0322-4
Huan, X., Lu, H., Wang, C., Tang, X., Zuo, X., Ge, Y., et al. (2015). Bulliform phytolith research in wild and domesticated rice paddy soil in south China. PLoS One 10 (10), e0141255. doi: 10.1371/journal.pone.0141255
Huang, X., Kurata, N., Wei, X., Wang, Z. X., Wang, A., Zhao, Q., et al. (2012). A map of rice genome variation reveals the origin of cultivated rice. Nature 490 (7421), 497–501. doi: 10.1038/nature11532
Issaharou-Matchi, I., Barboni, D., Meunier, J. D., Saadou, M., Dussouillez, P., Contoux, C., et al. (2016). Intraspecific biogenic silica variations in the grass species Pennisetum pedicellatum along an evapotranspiration gradient in South Niger. Flora 220, 84–93. doi: 10.1016/j.flora.2016.02.008
Jones, M. K., Liu, X. (2009). Origins of agriculture in East Asia. Science 324 (5928), 730–731. doi: 10.1126/science.1172082
Katz, O. (2019). Conflict and complementarity of paleontological and molecular chronologies? Paleobiology 45 (1), 7–20. doi: 10.1017/pab.2018.44
Larson, G., Piperno, D. R., Allaby, R. G., Purugganan, M. D., Andersson, L., Arroyo-Kalin, M., et al. (2014). Current perspectives and the future of domestication studies. Proc. Natl. Acad. Sci. U.S.A. 111 (17), 6139–6146. doi: 10.1073/pnas.1323964111
Liu, L., Chen, X. (2012). The Archaeology of China: From the Late Paleolithic to the Early Bronze Age. Cambridge: Cambridge University Press. doi: 10.1017/CBO9781139015301
Liu, L., Lee, G. A., Jiang, L. P., Zhang, J. Z. (2007). Evidence for the early beginning (c. 9000 cal. BP) of rice domestication in China: a response. Holocene 17 (8), 1059–1068. doi: 10.1177/0959683607085121
Liu, L., Phillips, S. M. (2006). ORYZA. @ in Flora of China Vol. 22. Wu, ZY, Raven, PH, Hong, DY, editors. Beijing: Science Press, 182–184.
Londo, J. P., Chiang, Y. C., Hung, K. H., Chiang, T. Y., Schaal, B. A. (2006). Phylogeography of Asian wild rice, Oryza rufipogon, reveals multiple independent domestications of cultivated rice, Oryza sativa. Proc. Natl. Acad. Sci. U.S.A. 103 (25), 9578–9583. doi: 10.1073/pnas.0603152103
Lu, H. Y., Liu, Z. X., Wu, N. Q., Berne, S., Saito, Y., Liu, B. Z., et al. (2002). Rice domestication and climatic change: phytolith evidence from East China. Boreas 31 (4), 378–385. doi: 10.1111/j.1502-3885.2002.tb01081.x
Luo, W. H., Yang, Y. Z., Yao, L., Chen, Z. J., Li, J., Yin, C. L., et al. (2016). Phytolith records of rice agriculture during the Middle Neolithic in the middle reaches of Huai River region, China. Quat. Int. 426, 133–140. doi: 10.1016/j.quaint.2016.03.010
Ma, X. L., Fang, J. Y. (2007). Silicas in leaves of eight wild rice species. Acta. Bot. Boreal.-Occident. Sin. 27 (8), 1531–1536 (in Chinese).
Molina, J., Sikora, M., Garud, N., Flowers, J. M., Rubinstein, S., Reynolds, A., et al. (2011). Molecular evidence for a single evolutionary origin of domesticated rice. Proc. Natl. Acad. Sci. U.S.A. 108 (20), 8351–8356. doi: 10.1073/pnas.1104686108
Pan, Y. (2008). Immature wild rice harvesting at Kuahuqiao, China? Antiquity 82 (316). Project Gallery, online (www.antiquity.ac.uk/ProjGall/pan/).
Pan, Y. (2011). The process of early agricultural emergence in the Yangzi Delta and Qiantang Drainage during 10000-6000 BP in terms of human ecology and archaeobotany. [dissertation/doctor’s thesis]. [Shanghai]: Fudan University (in Chinese).
Pearsall, D., Piperno, D., Dinan, E., Umlauf, M., Zhao, Z., Benfer, R. (1995). Distinguishing rice (Oryza sativa Poaceae) from wild Oryza species through phytolith analysis: results of preliminary research. Econ. Bot. 49 (2), 183–196. doi: 10.1007/BF02862923
Piperno, D. R. (2006). Phytoliths: a Comprehensive Guide for Archaeologists and Paleoecologists. New York: AltaMira Press.
Qin, L. (2012). “Archaeobotanical research and expectations on Chinese argiculture origin,” in Archaeological Research (9), eds. School of Archeaology and Museology of Peking University, Center of the Study of Chinese Archaeology, Peking University (Beijing: Cultural Relics Press), 260–315 (in Chinese).
Qin, L., Fuller, D. Q., Harvey, E. (2006). Subsistence of Hemudu site, and reconsideration of issues in the study of early rice from lower Yangtze. Dongfang Kaogu 3, 307–350 (in Chinese).
Qiu, Z. W., Jiang, L. P., Wang, C. S., David, V., Wu, Y. (2019). New evidence for rice cultivation from the Early Neolithic Hehuashan site. Archaeol. Anthropol. Sci. 11 (4), 1259–1272. doi: 10.1007/s12520-018-0602-1
Stein, J. C., Yu, Y., Copetti, D., Zwickl, D. J., Zhang, L., Zhang, C., et al. (2018). Genomes of 13 domesticated and wild rice relatives highlight genetic conservation, turnover and innovation across the genus Oryza. Nat. Genet. 50 (2), 285–296. doi: 10.1038/s41588-018-0040-0
Silva, F., Stevens, C. J., Weisskopf, A., Castillo, C., Qin, L., Bevan, A., et al. (2015). Modelling the geographical origin of rice cultivation in Asia using the rice archaeological database. PLoS One 10 (9), e0137024. doi: 10.1371/journal.pone.0137024
Song, Z., Lu, B., Zhu, Y., Chen, J. (2003). Gene flow from cultivated rice to the wild species Oryza rufipogon under experimental field conditions. New Phytol. 157 (3), 657–665. doi: 10.1046/j.1469-8137.2003.00699.x
Song, Z., Zhu, W., Rong, J., Xu, X., Chen, J., Lu, B.-R. (2006). Evidences of introgression from cultivated rice to Oryza rufipogon (Poaceae) populations based on SSR fingerprinting: implications for wild rice differentiation and conservation. Evol. Ecol. 20 (6), 501–522. doi: 10.1007/s10682-006-9113-0
Tanaka, K., Honda, T., Ishikawa, R. (2010). Rice archaeological remains and the possibility of DNA archaeology: examples from Yayoi and Heian periods of Northern Japan. Archaeol. Anthropol. Sci. 2 (2), 69–78. doi: 10.1007/s12520-010-0036-x
Wang, C., Lu, H. Y. (2012). Research progress of fan-shaped phytolith of rice and relative issues. Quat. Sci. 32 (2), 269–281 (in Chinese). doi: 10.3969/j.issn.1001-7410.2012.02.13
Wang, C., Lu, H. Y., Gu, W. F., Zuo, X. X., Zhang, J. P., Liu, Y. F., et al. (2018b). Temporal changes of mixed millet and rice agriculture in Neolithic-Bronze Age Central Plain, China: Archaeobotanical evidence from the Zhuzhai site. Holocene 28 (5), 738–754. doi: 10.1177/0959683617744269
Wang, C. L., Udatsu, T., Fujiwara, H., Tang, L. H., Zou, J. S. (1997). Morphological characteristics of silica bodies from motor cells in rice (Oryza sativa L.) and their differences between indica and japonica. Jiangsu J. Agr. Sci. 13 (3), 129–138 (in Chinese).
Wang, H., Vieira, F. G., Crawford, J. E., Chu, C., Nielsen, R. (2017). Asian wild rice is a hybrid swarm with extensive gene flow and feralization from domesticated rice. Genome Res. 27, 1029–1038. doi: 10.1101/gr.204800.116
Wang, W., Mauleon, R., Hu, Z., Chebotarov, D., Tai, S., Wu, Z., et al. (2018a). Genomic variation in 3,010 diverse accessions of Asian cultivated rice. Nature 557 (7703), 43–49. doi: 10.1038/s41586-018-0063-9
Wang, Y. J., Lu, H. Y. (1993). The Study of Phytolith and Its Application. Beijing: China Ocean Press (in Chinese).
Wei, X., Qiao, W. H., Chen, Y. T., Wang, R. S., Cao, L. R., Zhang, W. X., et al. (2012). Domestication and geographic origin of Oryza sativa in China: insights from multilocus analysis of nucleotide variation of O. sativa and O. rufipogon. Mol. Ecol. 21 (20), 5073–5087. doi: 10.1111/j.1365-294X.2012.05748.x
Weisskopf, A., Qin, L., Ding, J., Sun, G., Fuller, D. (2015). Phytoliths and rice: from wet to dry and back again in the Neolithic Lower Yangtze. Antiquity 89 (347), 1051–1063. doi: 10.15184/aqy.2015.94
Wu, Y., Jiang, L. P., Zheng, Y. F., Wang, C. S., Zhao, Z. J. (2014). Morphological trend analysis of rice phytolith during the early Neolithic in the Lower Yangtze. J. Archaeol. Sci. 49, 326–331. doi: 10.1016/j.jas.2014.06.001
Zhang, J., Lu, H., Gu, W., Wu, N., Zhou, K., Hu, Y., et al. (2012). Early mixed farming of millet and rice 7800 years ago in the Middle Yellow River region, China. PLoS One 7 (12), e52146. doi: 10.1371/journal.pone.0052146
Zhang, W. X., Wang, L. L. (1998). Phytoliths in leaves of 7 Oryza species. J. China Agr. Univ. 3 (3), 21–25 (in Chinese).
Zhao, K., Wright, M., Kimball, J., Eizenga, G., McClung, A., Kovach, M., et al. (2010). Genomic diversity and introgression in O. sativa reveal the impact of domestication and breeding on the rice genome. PLoS One 5 (5), e10780. doi: 10.1371/journal.pone.0010780
Zhao, Z. J., Pearsall, D., Benfer, R., Piperno, D. (1998). Distinguishing rice (Oryza sativa poaceae) from wild Oryza species through phytolith analysis, II finalized method. Econ. Bot. 52 (2), 134–145. doi: 10.1007/BF02861201
Zhao, Z. J. (2011). New archaeobotanic data for the study of the origins of agriculture in China. Curr. Anthropol. 52, S295–S306. doi: 10.1086/659308
Zhao, Z. J., Jiang, L. P. (2016). Analyses of plant remains from floatation work carried out at the Shangshan site, Pujiang County, Zhejiang Province. Relics From South 3, 109–116 (in Chinese). doi: 10.3969/j.issn.1004-6275.2016.03.013
Zheng, Y. F., Crawford, G. W., Jiang, L., Chen, X. (2016). Rice domestication revealed by reduced shattering of archaeological rice from the Lower Yangtze valley. Sci. Rep. 6, 28136. doi: 10.1038/srep28136
Zheng, Y. F., Dong, Y. J., Matsui, A., Udatsu, T., Fujiwara, H. (2003a). Molecular genetic basis of determining subspecies of ancient rice using the shape of phytoliths. J. Archaeol. Sci. 30 (10), 1215–1221. doi: 10.1016/S0305-4403(02)00248-0
Zheng, Y. F., Jiang, L. P., Zheng, J. M. (2004). Study of the remains of ancient rice from the Kuahuqiao site in Zhejiang Province. Chin. J. Rice. Sci. 18 (2), 119–124 (in Chinese). doi: 10.3321/j.issn.1001-7216.2004.02.007
Zheng, Y. F., Matsui, A., Fujiwar, H. (2003b). Phytoliths of rice detected in the Neolithic sites in the valley of the Taihu Lake in China. Environ. Archaeol. 8, 177–183. doi: 10.1179/env.2003.8.2.177
Zheng, Y. F., Sun, G. P., Chen, X. G. (2007). Characteristics of the short rachillae of rice from archaeological sites dating to 7000 years ago. Chin. Sci. Bull. 52 (12), 1654–1660. doi: 10.1007/s11434-007-0258-1
Keywords: rice, bulliform phytolith, Oryza sativa, Oryza rufipogon, domestication, morphometric analysis
Citation: Wang C, Lu H, Zhang J, Mao L and Ge Y (2019) Bulliform Phytolith Size of Rice and Its Correlation With Hydrothermal Environment: A Preliminary Morphological Study on Species in Southern China. Front. Plant Sci. 10:1037. doi: 10.3389/fpls.2019.01037
Received: 02 May 2019; Accepted: 24 July 2019;
Published: 22 August 2019.
Edited by:
Rivka Elbaum, The Hebrew University of Jerusalem, IsraelReviewed by:
Wricha Tyagi, Central Agricultural University, IndiaDongying Gao, University of Georgia, United States
Ofir Katz, Dead Sea and Arava Science Center, Israel
Copyright © 2019 Wang, Lu, Zhang, Mao and Ge. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Can Wang, c2hhbmRhd2FuZ2NhbkAxNjMuY29t; Houyuan Lu, aG91eXVhbmx1QG1haWwuaWdnY2FzLmFjLmNu