- Open Access
Mapping the potential distribution of high artemisinin-yielding Artemisia annua L. (Qinghao) in China with a geographic information system
Chinese Medicine volume 5, Article number: 18 (2010)
Artemisia annua L. is an important source for artemisinin, a potent drug for treating malaria. This study aims to map and predict the potential geographic distribution of A. annua L. in China.
The Geographic Information System for traditional Chinese medicine (TCM-GIS) was developed and used to map the potential geographic distribution of A. annua L.
Climatic, edaphic and topographic characteristics of A. annua L. microhabitats in Youyang County were mapped to find distribution patterns. The maps identified that certain habitats in the Chongqing region and some potential regions, especially in Guizhou Province, possess similarity indices of ≥98%. In particular, high quality microhabitats A. annua L. were found in the Wuling mountains region.
The present study demonstrates a GIS approach to predict potential habitats for A. annua L. TCM-GIS is a powerful tool for assessing bioclimatic suitability for medicinal plants.
Artemisia annua L. (Qinghao, Annual Wormwood) is a strongly fragrant, annual herbaceous plant used in Chinese medicine . A. annua L. is the only natural botanical source for artemisinin (Qinghaosu) [2, 3] and a potential source for essential oils for the perfume industry . A. annua L. is now cultivated in China, Vietnam, India, Romania, Kenya and Tanzania . Artemisinin, an endoperoxide sesquiterpene lactone in the aerial parts of A. annua L., is more efficacious, faster and less toxic than chloroquine in treating malaria. In addition, artemisinin is a potent anti-cancer agent, a possible antibacterial agent as well as a natural pesticide [6, 7]. Chemical and biological synthesis of artemisinin is still under development due to poor yields [8–11]. Therefore, wild or cultivated A. annua L. is a major source for artemisinin [2, 3, 12].
The artemisinin content is highly dependent on plant ecotypes, ecological interactions, seasonal and geographic variations [13–18]. In fact, artemisinin is absent in some A. annua L. Artemisinin was first isolated in China and some Chinese germplasm has relatively higher artemisinin levels than those of Europe, North America, East Africa and Australia [2, 13, 16, 17, 19, 20]. In Youyang County, Chongqing, China, the hometown of A. annua L., the plants have high (0.9%) levels of artemisinin. In 2006 the county became a national protected geographic area recognized by the General Administration of Quality Supervision, Inspection and Quarantine of China . As the demand for artemisinin remains high around the world, finding suitable geographic regions for A. annua L. is a critical research area for the World Health Organization .
The geographic information system (GIS) technology manages geographic information with applications for various fields such as natural resources, transportation planning, environmental studies and vegetation distribution studies [23–26]. Recently updated, the geographic information system for traditional Chinese medicine (TCM-GIS) captures, stores, analyzes and displays geographically referenced information to analyze genetic, ecological and geographic patterns of the spatial distribution of a target species. Using the TCM-GIS, our previous studies analyzed the potential habitats and distributions of Chinese medicinal plants such as Glycyrrhiza uralensis Fisch., Panax quinquefolium and Panax ginseng[27–29]. The present study aims to characterize the eco-environmental conditions in the A. annua L. production areas in Youyang County and predict the potential distributions of A. annua L. with a high artemisinin-yielding potentials.
The spatial distribution of A. annua L. was based on the following four sources: (1) the flora of China ; (2) scientific literature concerning the geographic distribution of A. annua L. in China ; (3) the Chinese Virtual Herbarium (CVH) , (4) germplasm accessions from the Sharing Information System for Chinese Medicinal Plant Germplasm Resources ; (5) field data of wild A. annua L. and interviews in Youyang County in 2008. Due to the excellent quality of A. annua L. from the habitats in Youyang County [31, 34–36], a total of 180 accessions of A. annua L. germplasm were collected and used in the present study.
The potential distribution mapping program TCM-GIS and geo-referenced datasets were used to develop eco-adaptation models. The TCM-GIS package included three databases, namely (1) a basic geographic information database including digital line graphics and a digital elevation model (scale: 1:1,000,000), (2) a soil database (scale: 1:4,000,000), (3) and a climate database (mean values between 1971 and 2000). All three databases were used for spatial analysis and model calibration.
Raster and vector are two main data models in the TCM-GIS. Raster layers (1 × 1 km2 resolution) were used for the eco-environmental analysis and cluster analysis. Vector layers were used to derive and identify the spatial extent and location of suitable habitats through overlay analysis. Moreover, global positioning system data on the locations of the 180 accessions were obtained for villages such as Banqiao, Zhongduo, Mawang and Nanmu and used in the TCM-GIS analysis (Figure 1).
In the present study, 14 eco-environmental variables were chosen for the predication of spatial distribution in Youyang County. These variables, namely (1) average temperature in January (ATJA), (2) average temperature in February (ATF), (3) average temperature in March (ATM), (4) average temperature in April (ATAP), (5) average temperature in May (ATMA), (6) average temperature in June (ATJ), (7) average temperature in July (ATJU), (8) average temperature in August (ATA), (9) average annual temperature (AAT), (10) annual sunshine time (AST), (11) total annual precipitation (TAP), (12) relative humidity (RH), (13) altitude (AL), (14) and soil properties (SP), were classified into three categories: topography, climate and edaphology (Table 1).
An optimal range was established by identifying minima and maxima for eco-environmental variables (e.g. elevation and temperature) at sample collection sites. The A. annua L. macro-habitats were characterized by examining the mean, minimal and maximal values, standard deviation (SD), standard error (SE), and coefficient of variation (CV) of these variables (Table 2). Prior to distance analysis, we normalized the raster grid data representing each variable. We derived the mean absolute deviation using the following equation:
where xkf was the measured values of the variable f and mf is the mean for the variable f. For the determination of similarity between grid data and eco-factor ranges, the statistical distance was calculated with the Minkowski distance equation :
which is a generalization of the Euclidean distance and Manhattan distance; in general the shorter the distance, the greater the similarity. The comprehensive similarity index (SI) of each factor layer was calculated with an overlay analysis with various weighting values. Finally, maps with two ranks of predictive distributions were generated, followed by a grid-based spatial cluster analysis, vector-based overlaying, intersection analysis and an area calculation (Figures 2, 3, 4, Table 3).
The most favorable region for A. annua L. growth is one that has an SI range of 99%-100%, while the second-most favorable region is one that has an SI range of 98%-99%.
Results and Discussion
The climatic, edaphic and topographic characteristics of known A. annua L. habitats are listed in Table 2. While low CV values for RH (CV: 0.33), TAP (1.28), AST (3.33), ATJU (4.60), AAT (4.69), ATA (6.23), ATJ (6.77) and ATMA (6.81) suggested that these could be the major limiting factors affecting the distribution of high quality A. annua L., high CV values for AL (29.79), ATJA (21.46) and ATF (21.43) suggested otherwise. According to the CV values, weighting value for each parameter was divided into levels I (0.15), II (0.08), III (0.06) and IV (0.03) and weighting values should add up to one. In addition, datasets of eco-factors from known habitats in Youyang County were as follows: ATJA = 1.2-5.6°C, ATF = 2.0-6.0°C, ATM = 4.0-10.0°C, ATAP = 10.0-16.0°C, ATMA = 14.0-20.0°C, ATJ = 18.0-24.0°C, ATJU = 21.6-27.3°C, ATA = 20.0-26.0°C, AAT = 15.9-21.0°C, AST = 1048-1200 h, TAP = 1169-1267 mm, RH = 79.2-80.6%, AL = 498-1010 mm. Soil types were mainly yellow soil, yellow sandy soil, limestone soil, paddy soil and brown soil with pH value at 6-7 and organic matter content ≥1.3%. Thus, we assumed that these conditions were optimal for the growth of high artemisinin-yielding A. annua L.
A. annua L. is a short-day plant. Non-juvenile plants are very responsive to short photoperiodic stimuli and flower about two weeks after induction. They require about 1000 hours of sunlight per year. Our results suggest that annual sunlight time is a critical factor for the growth of A. annua L., which is consistent with previous studies [5, 38]. Previous findings that A. annua L. requires a strict watering regime during the preliminary growth stages [5, 39] are also consistent with our results.
Figures 2 and 3 are the maps derived from the TCM-GIS analyses. The predicted areas were primarily located in the Wuling Mountain region in central China, covering Guizhou, Chongqing, Hunan, Hubei and Sichuan (25°14'-31°38' N to 104°31'-111°51'E). The predicted habitat density was high in northeastern Guizhou, southeastern Chongqing, northwestern Hunan, southwestern Hubei and parts of southern Sichuan.
The total favorable regions (SI 98%-99%) made up 1.60% of China's total land area covering 162 counties and cities (a total of 60,292 km2), among which Guizhou took the lead with 31,150 km2 including 68 counties and cities. The most favorable region for A. annua L. (SI 99%-100%) was in the 58 counties and cities in Guizhou Province with a predicted area of 54,350 km2. The second largest predicted area (14,330 km2) was in the 12 counties and cities in Chongqing, followed by Hunan, Hubei and Sichuan (Figure 4). The counties and cities with significant areas of potential habitat are listed in Table 3. The data indicated that Youyang County contained the largest favorable area with more than 4000 km2. Unexpectedly, the total predicted areas in Wuchuan and Zunyi Counties in Guizhou exceeded 2000 km2.
One of the world's largest artemisinin manufacturers and its affiliates operate A. annua L. farms in the Chongqing Wulingshan Mountain Range [40, 22]. Apart from this, Guizhou may be another important region for A. annua L. cultivation, particularly in the northeastern part of the province. Our model predicted that 13% of this area is potential A. annua L. habitats [41, 42]. Our model did not predict Guangxi Province, known for its habitats of A. annua L. of relatively low quality, as a region for A. annua L. cultivation possibly due to the subtropical climate, low altitude and red soil in Guangxi which are very different from those in other A. annua L. regions in China .
Interviews with the locals suggest that the Guizhou region and Youyang County have comparative advantages for A. annua L. growth with a high-yield variety and minimal pests. Furthermore, the northeastern Guizhou is home to wild populations of A. annua L. which may be an alternative source for artemisinin.
Using the TCM-GIS, we aimed to determine the optimal ecological factors from known habitats and the results showed that RH, TAP, AST, STJU, AAT and SP were important limiting factors. We also aimed to map the distribution of potential regions for the development of A. annua L. in China based on selected climatic, soil and topographical values. Using bioclimatic similarity theory and the TCM-GIS, we predicted the potential growing areas at the county level, particularly in northeastern Guizhou Province. The TCM-GIS is adequate for predicting and identifying potential areas for A. annua L. cultivation.
Using a higher resolution raster and vector spatial databases, we improved the resolution of species distribution considerably on the national surveys conducted in the 1960s, 1970s and 1980s. While most of the survey data were based largely on personal experiences and rough estimates, the model used in the present study is relatively objective.
The present study demonstrates a GIS approach to predict the potential habitats for A. annua L. TCM-GIS is a powerful tool for assessing bioclimatic suitability for medicinal plants.
traditional Chinese medicine geographic information system
geographic information system
coefficient of variation
average temperature in January
average temperature in February
average temperature in March
average temperature in April
average temperature in May
average temperature in June
average temperature in July
average temperature in August
average annual temperature
annual sunshine time
total annual precipitation
Chinese Virtual Herbarium.
Wang CW: The forests of China: with a survey of grassland and desert vegetation. Maria Moors Cabot Foundation Publication. 1961, Cambridge Massachusetts: Harvard University, 5: 155-164.
Jain DC, Mathur AK, Gupta MM, Singh AK, Verma RK, Gupta AP, Kumar S: Isolation of high artemisinin-yielding clones of Artemisia annua. Phytochemistry. 1996, 5: 993-1001. 10.1016/S0031-9422(96)00369-X.
Klayman DL: Qinghaosu (artemisinin): an antimalarial drug from China. Science. 1985, 228: 1049-1055. 10.1126/science.3887571.
Simon JE, Charles E, Cebert L, Grant J, Janick J, Whipkey A: Artemisia annua L: a promising aromatic and medicinal. Advances in New Crops: Proceeding of the First National Symposium New Crops: Research, Development, Economics. 1990, Portland: Timber Press, Incorporated, 522-526.
World Health Organization: WHO monograph on good agricultural and collection practices (GACP) for Artemisia annua L. Geneva. 2006
Efferth T, Dunstan H, Sauerbrey A, Miyachi H, Chitambar CR: The anti-malarial artesunate is also active against cancer. Int J Oncol. 2001, 18: 767-773.
Kumar S, Khanuja SPS, Kumar TRS, Jain DC, Srivastava S, Bhattacharya AK, Saikia D, Shasany AK, Darokar MP, Sharma RP: Method for the use of alpha arteether as an anti-bacterial and anti-fungal agent. US Patent 6127405. 2000
Avery MA, Chong WKM, Jennings-White C: The total synthesis of (+)-Artemisinin and (+)-9-Desmethylartemisinin. Tetrahedron Lett. 1987, 28: 4629-4632. 10.1016/S0040-4039(00)96582-1.
Ravindranathan T, Kumar MA, Menon RB, Hiremath SV: Stereoselective synthesis of artemisinin. Tetrahedron Lett. 1990, 31: 755-758. 10.1016/S0040-4039(00)94621-5.
Van Geldre E, Vergauwe A, Eeckhout Van den E: State of art of the production of the antimalarial compound artemisinin in plants. Plant Mol Biol. 1997, 33: 199-209. 10.1023/A:1005716600612.
Chen DH, Li GF, Ye HC: Expression of a chimeric farnesyl diphosphate synthase gene in Artemisia annua L. transgenic plants via Agrobacterium tumefaciens-mediated transformation. Plant Sci. 2000, 155: 179-185. 10.1016/S0168-9452(00)00217-X.
Kumar S, Gupta SK, Singh P, Bajpai P, Gupta MM, Singh D, Gupta AK, Ram G, Shasany AK, Sharma S: High yields of artemisinin by multi-harvest of Artemisia annua crops. Ind Crops Prod. 2004, 1: 77-90. 10.1016/j.indcrop.2003.07.003.
Delabays N, Benakis A, Collet G: Selection and breeding for high artemisinin (Qinghaosu) yielding strains of Artemisia annua. Acta Hort (ISHS). 1993, 330: 203-208.
Liersch R, Soicke H, Stehr C, Tullner HU: Formation of artemisinin in Artemisia annua during one vegetation period. Planta Med. 1986, 52: 387-390. 10.1055/s-2007-969193.
Singh A, Vishwakarma RA, Husain A: Evaluation of Artemisia annua strains for higher artemisinin production. Planta Med. 1988, 64: 475-476. 10.1055/s-2006-962515.
Charles DJ, Simon JE, Wood KV, Heinstein P: Germplasm variation in artemisinin content of Artemisia annua using an alternative method of artemisinin analysis from crude plant extracts. J Nat Prod. 1990, 53: 157-159. 10.1021/np50067a021.
Woerdenbag HJ, Pras N, Chan NG, Bang BT, Bos R, van Uden W, Van YP, van Boi NV, Batterman S, Lugt CB: Artemisinin related sesquiterpenes and essential oil in Artemisia annua during a vegetation period in Vietnam. Planta Med. 1994, 60: 272-275. 10.1055/s-2006-959474.
Gupta SK, Singh P, Bajpai P, Ram G, Singh D, Gupta MM, Jain DC, Khanuja SPS, Kumar S: Morphogenetic variation for artemisinin and volatile oil in Artemisia annua. Ind Crops Prod. 2002, 16: 217-224. 10.1016/S0926-6690(02)00049-3.
Trigg EI: Qinghaosu (artemisinin) as an antimalarial drug. Econ Med Plant Res. 1989, 3: 19-55.
Klayman DL, Lin AJ, Acton N, Scovill JP, Hoch JM, Milhous WK, Theoharides AD, Dobek AS: Isolation of artemisinin (qinghaosu) from Artemisia annua growing in the United States. J Nat Prod. 1984, 47: 715-717. 10.1021/np50034a027.
Genearal Administration of Quality Supervision, Inspection and Quarantine of the People's Republic of China. [http://kjs.aqsiq.gov.cn/xxgkml/ywxx/dlbzcpbh/dlbzcpzybzzcdj/bsjg/200710/t20071012_44016.htm]
World Health Organization Global Malaria Programme: Proceedings of the Meeting on the production of artemisinin and artemisinin-based combination therapies: 6-7 June 2005; Arusha. 2006, Global Malaria Programme
Goodchild MF: Geographic information systems and science: today and tomorrow. Procedia: Earth and Planetary Science. 2009, 1: 1037-1043. 10.1016/j.proeps.2009.09.160.
Segura S, Coppens d'Eeckenbrugge G, Lòpez L, Grum M, Guarino L: Mapping the potential distribution of five species of Passiflora in Andean countries. Genet Resour Crop Evol. 2003, 12: 555-566. 10.1023/A:1024424012134.
Chen SL, Suo FM, Han JP, Xie CX, Yao H, Li XW, Li Y, Wei JH: Analysis on ecological suitability and regionalization of traditional Chinese medicinal materials. Chin Tradit Herb Drugs. 2007, 38: 481-487.
Smith-Ramirez C, Diaz I, Pliscoff P, Valdovinos C, Mendez MA, Larrain J, Samaniego H: Distribution patterns of flora and fauna in southern Chilean coastal rain forests: Integrating natural history and GIS. Biodivers and Conserv. 2007, 16: 2627-2648. 10.1007/s10531-006-9073-2.
Wang JY, Zhao RH, Sun CZ, Chen SL Liu SQ, Wei JH, Xiao XX, Zhou YQ: Suitability evaluation of Glycyrrhiza Uralensis Fisch.'s distributive area based on TCMGIS. Mod Chin Med. 2006, 8: 4-8.
Chen SL, Zhou YQ, Xie CX, Zhao RH, Sun CZ, Wei JH, Liu SQ, Gao WW: Suitability evaluation of Panax quinquefolium producing area based on TCMGIS. Zhongguo Zhong Yao Za Zhi. 2008, 33: 741-745.
Wang Y, Wei JH, Chen SL, Sun CZ, Zhao RH, Liu SQ, Xiao XX, Wang JY, Zhou YQ: Analysis of Panax Ginseng' producing area based on TCMGIS. Asia-Pacific Trad Med. 2006, 6: 73-78.
Academiae Sinicae Edita: Flora Reipublicae Popularis Sinicae. 1989, Beijing: Science Press, 76 (2): 62-65.
Zhong GY, Zhou HR, Ling Y, Hu M, Zhao PP: Investigation on ecological environment and quantitative analysis of Artemisinin of Sweet Wormwood (Artemisia annua). Chin Med Mat. 1998, 29: 264-267.
The Chinese Virtual Herbarium (CVH). [http://www.cvh.org.cn]
Germplasm Accessions From the Sharing Information System for Chinese Medicinal Plant Germplasm Resources. [http://www.tcm-resources.com/IS/querymain.asp]
Yang SP, Yang X, Huang JG, Ding DR: Advances in researches on Artemisinin production. J Trop Subtrop Bot. 2004, 12: 189-194.
Xu XY, Zheng YM, Fu SQ, Han YM, Yang YH: Studies on HPLC Fingerprint of A rtemisia annua L. from Chongqing. Lishizhen Med Mater Med Res. 2009, 20: 1188-1189.
Xia ZZ: Studies on the methods for assay of artemisinin in Artemisia annual L. and its semi-synthetic derivatives. Master thesis. 2006, Chongqing Medical University, School of Pharmaceutical Sinence
He YQ: Multivariate Statistical Analysis. 2004, Beijing: China Renmin University Press
Purdue University Annual Wormwood (Artemisia annua L.). [http://www.hort.purdue.edu/newcrop/cropfactsheets/artemisia.pdf]
Wei X, Li F, Xu CQ, Fu XH, Huang ZF: Study on biological characteristics of Artemisia annua L. Guihaia. 1997, 17: 166-168.
Holley Pharmaceuticals Company. [http://holleypharma.com/page/meuf/Holley_Pharmaceuticals.html]
Zhou AA, Zhen WX, Ge FH: Determination of artemisinin in Artemisia annuna by HPLC-ELSD. Zhong Yao Cai. 2006, 29: 242-245.
Wang YL, Zhang L, Zhang HM: Determination of artemisinin in Artemisia annua from different habitats in Changde region Hunan province. Pharmaceut care Res. 2007, 5: 381-382.
The authors would like to thank the National Natural Science Foundation of China for its support through project No. 3050081.
The authors declare that they have no competing interests.
LFH, SLC and CXX designed the study and performed the analyses. BZD helped with data analysis. All authors wrote the manuscript. All authors read and approved the final version of the manuscript.