介绍

Ducrosia是开花的属药材属于家庭吗伞形科d . anethifolia是三个Ducrosia物种广泛分布在南部埃及、伊朗、伊拉克、阿富汗和巴基斯坦1。这种草药在伊朗通常称为Moshkbu和Moshgak2。这芳香的天线部分药用植物应用镇痛缓解背痛,头痛,感冒,绞痛3。此外,Moshgak添加到许多食品调味料在一些地区的伊朗2

芳香精油(EOs)是微妙的,和挥发性液体从花等植物材料中提取根,树皮、种子、水果皮,和木头4。大约90%的全球光电生产消耗的风味和芳香产业,大部分是在化妆品、香水、饮料和食物。此外,超过25%的全球药品处方来源于植物来源。因此,plant-originated产品将会越来越需要研究者和生产者继续寻找新的植物来源5。单萜是最常见的精油在医学草药的代谢物6。然而,Sefidkon Javidtash7正癸醇癸醛,α-thujene,香紫苏醇,柠檬烯为主要成分的干Moshgak天线的部分。

今天,许多生物活性成分的使用,如酚醛树脂和黄酮类化合物对人类健康的伟大价值6。最近,许多研究已经报道精油的组成、酚类物质和类黄酮抗氧化活动在不同医学草药。化合物的数量和他们的生物活性的能力可能会受到遗传和环境因素的影响8,9。干旱的一个主要环境因素影响代谢产物的浓度10。干旱或渗透压力导致过度生成活性氧(ROS),可以在高浓度有害细胞的组件11。植物展示一个广泛的抗压力或压力避免机制生存压力条件10。植物在干旱胁迫条件下主要产生较高含量的次生代谢物与有条件12。有很多研究,突出了水分胁迫对伞形科植物的次生代谢物的影响,包括Cuminum植物1314。总酚类(TPC)和总类黄酮(交通)和抗氧化活性也可以提高在不同的药用植物水分胁迫的反应1,8,14,15。同时,Mottaghipisheh et al。16显示的抗氧化活性d . anethifoliaethanolic和乙酸乙酯提取物小于丁羟甲苯(二叔丁基对甲酚)作为合成抗氧化剂,是作为积极的控制。评估各年人口的大量的繁殖计划。此外,种植一年可以高度影响药用植物的次生代谢产物等9,17,18lavandulifolia士达古19

有有限的研究d . anethifolia人群。Sabbaghi和Rahimmalek20.评估伊朗Moshgak种群的形态和分子变异。Arbabi et al。2还评价了形态和油组件的叶子在这个物种在一个地理区域在伊朗东南部地区。然而,没有研究进行水分胁迫的影响在精油的成分,酚醛树脂,黄酮类化合物,不同人群的抗氧化能力d . anethifolia在连续两年气相色谱分析结果提供有价值的数据有许多信息是很难的过程基于单变量方法,多元分析提供新的工具,用于探索和色谱数据的解释21,22,23

本研究的目的是:(1)评估的不同人群D.anethifolia实验场条件下选择人口高植物化学的内容;(2)评估化学成分和水分胁迫条件和多酚含量(3)使用多元分析选择和引入精英人群进行进一步培养和处理。

结果

水压力和人口对精油的影响的内容

之间的方差分析结果显示显著差异Moshgak人口精油的研究内容。此外,今年的影响、灌溉、交互和灌溉××的人口人口EOs(表的相互作用是重要的1)。提取的精油含量的颜色是淡黄色的。然而,精油在连续两年的趋势(表2)表明,更高的EOs发生在第二年。EOs内容也增加了在干旱胁迫下(表1)。

表1方差分析(意思是正方形)的植物化学成分特征24 Moshgak人群(P)三个层次的评价灌溉(I)和2年(Y)的2019年和2020年。
表2比较的灌溉制度,年,人口(交互)8 Moshgak人口的植物化学的特征。

在第一个农艺年(2019),最高和最低观察精油含量的增加的数量Abadeh(83%)和Ferydunshahr(52%),分别为(表2)。

在农艺第二年(2020年),水分亏缺造成显著的精油在所有Moshgak人群。Qazvin市最高的EOs增强观察和Abarkuh2人群(83%),而Abarkuh1最低(56%)(表2)。

干旱胁迫和人口对精油成分的影响

精油成分表中列出3。的气相色谱峰对应d . anethifolia和主要化合物的质谱见图S1。结果表明化学多态性Moshgak人群高,和26个EOs组件通过原产地灌溉制度和人口的影响。EOs的主要成分检测独联体-chrysanthenyl醋酸和α-pinene(表3)。结果还显示,独联体-chrysanthenyl醋酸是最丰富的组件。在控制条件下,这种化合物的最高和最低数量在Semirom(66.4%)和Abadeh(39.1%),分别为。此外,精油成分的浓度在所有Moshgak种群水分胁迫下显著升高(表3)。最高的和最低的增加独联体-chrysanthenyl醋酸在Abarkuh1(21.32%)和Hormozgan(6.42%),分别为(表3)。水分亏缺下,α-pinene科曼地毯数量最高的增长(6.1%),最低的在Qazvin(0.61%)(表3)。的最大增加观察trans-verbenol Abadeh人口。桧烯没有出现在Abadeh, Ferydunshahr1, Abarkuh1人群。然而,在水果从干旱强调植物桧烯注意到(表的出现3)。

表3的精油成分的八个Moshgak种群在三个级别的灌溉系统。

干旱胁迫和人口对酚类含量的影响

方差分析的结果显示显著差异Moshgak人口研究的总酚含量(TPC)和类黄酮含量(交通)。在目前的研究中,TPC和交通的影响,灌溉(表1)。之间的交互作用是显著的人口今年灌溉×××灌溉和人口特征。

在实验的第一个赛季,在控制条件下,“Ferydunshahr”表现出最高的TPC(95.87毫克的TAE DW / g),而“科曼地毯”显示最高的交通(4.13毫克的QE / g (DW)。总酚和类黄酮含量升高(表模拟干旱胁迫下所有人群2)。

第二年,在控制条件下,最高的观察TPC Ferydunshahr(99.34毫克的TAE DW / g),而最大的交通是记录在科曼地毯(4.68毫克的QE / g (DW)。干旱强调所有样本人口总酚和总黄酮含量有了显著提高。最高的TPC增加和交通的数量Abarkuh(52.87%)和科曼地毯(1.37%),分别为(表2)。

干旱胁迫和人口对DPPH自由基清除的影响分析

方差分析的结果表明Moshgak人口显著差异对DPPH自由基清除试验(集成电路研究50)。此外,今年的影响,灌溉,人口年灌溉××人口人口互动是重要的集成电路50(表1)。第一年,在控制条件下,抗氧化活性最高(最低的IC50值)获得了在Hormozgan(166.17µg /毫升),而Abarkuh1拥有最低的抗氧化活性(IC最低50值)(726.81µg /毫升)。在干旱胁迫下抗氧化活性增加(减少集成电路50在所有Moshgak数量(表)2)。

第二年,在控制条件下,最低和最高的集成电路50观察在Hormozgan(216.63µg /毫升)和Semirom(730.13µg /毫升),分别(表吗3)。此外,活动增加以应对干旱胁迫(表2)。

精油的主要组件之间的相关性和TPC和交通

精油的主要成分之间的相关系数和总酚类物质和类黄酮是在表的内容6。在控制条件下,最高coefficeintsα-pinene至β-pinene (r = 90 * *),以及myrcane和柠檬烯(r = 84 * *)。独联体-chrysanthenyl醋酸与β-pinene弱正相关(r = 0.66 *),但与柠檬烯略负相关(r =−0.46 *)(表6)。此外,p-cymene之间的关系和总酚醛树脂(r = 66 * *)和集成电路50(r = 0.48 *)是积极的,但与精油负面内容的相关性(r =−0.54 *)。然而,没有明显相关性独联体-chrysanthenyl醋酸和类黄酮或酚醛树脂(表的内容4)。

表4之间的皮尔逊相关系数主要精油成分的水果精油、总酚、总flavnoid八Moshgak人口控制和DPPH(较低的三角形),中度水分胁迫(上三角形)条件。

在中度水分胁迫的条件下,皮尔森相关矩阵表示和α-pinene之间显著正相关反式马鞭草烯醇(r = 0.95 * *), myrcane (r = 0.90 * *),β-pinene (r = 0.81 * *)和柠檬烯(r = 0.79 * *)(表5)。此外,之间的关系独联体-chrysanthenyl醋酸和α-pinene (r =−0.55 *),柠檬烯(r =−0.59 *),γ-terpinene (r =−66 * *),和trans-verbenol (r =−0.53 *) -(表4)。

院长的精油呈负相关(r =−0.74 * *)、交通(r =−0.47 *),但它呈正相关,TPC (r = 0.69 * *)。此外,之间不存在显著相关独联体-chrysanthenyl醋酸和植物化学的特征(表4)。

在严重水分胁迫条件下,α-pinene与trans-verbenol显著正相关(r = 0.97 * *),月桂烯(r = 0.87 * *)和桧烯(r = 0.83 * *),但相关负面了独联体-chrysanthenyl醋酸(r =−0.58 *)(表6)。同时,独联体-chrysanthenyl醋酸有负相关反式马鞭草烯醇(r =−0.57 *)和柠檬烯(r =−0.54 *),但它与IC略呈正相关50(r = 0.61 * *)(表5)。

表5之间的皮尔逊相关系数主要精油成分的水果精油、总酚、总类黄酮和DPPH八Moshgak人口严重的水分胁迫的条件下。

层次聚类分析和主成分分析

探讨可能的相似之处和人口之间的关系,进行分层聚类分析(HCA)基于2年的平均值(无花果。1)。HCA的结果可以应用于选择优越的父母基于合适的遗传距离和独联体-chrysanthenyl醋酸为目标穿越未来Moshgak育种程序。根据HCA, Moshgak人口研究被分为三组。集群由12正常和适度的水分胁迫条件的人群。第二组包括Abarkuh数量由中度和重度水分胁迫的条件下最高独联体p-cymene -chrysanthenyl乙酸,乙酸香叶酯、总酚和集成电路50。八的数量由中度和重度水分胁迫的条件下被安置在第三组。他们被高α-pinene杰出、桧烯β-pinene, myrcane,柠檬烯,γ-terpinene, trans-verbenol,癸醛、总类黄酮(无花果。1)。

图1
图1

8 Moshgak种群进行聚类分析的系统树图生成三个灌溉制度基于植物化学的字符使用基于平方计算欧几里得不同病房。

主成分分析(PCA)是基于主要的植物化学成分特征平均2年了。主成分分析证实了HCA(图的结果。2)。主成分分析结果表明,前两个组件解释总变异的39%和19%,分别(无花果。2)。第一个PC (PC1)与α-pinene正相关,桧烯,β-pinene, myrcane,柠檬烯,γ-terpinene,反式马鞭草烯醇、癸醛和类黄酮。第二个电脑(PC2)有关联的精油含量高,独联体p-cymene -chrysanthenyl乙酸,乙酸香叶酯和总酚和集成电路50。Abarkuh1和Abarkuh2种群在中度和重度水分胁迫条件下在同一组分类。Abadeh和科曼地毯人口在中度和重度水分胁迫的条件下互相分类附近。此外,Ferydunshahr Semirom Hormozgan,和Qazvin市人口严重水分胁迫条件下在同一组分类。其他人分为另一组。biplot分析的基础上,根据高PC1选择和PC2值会导致高的人群独联体-chrysanthenyl乙酸、乙酸香叶酯和酚醛树脂(无花果。2)。

图2
图2

PCA分析基于植物化学的字符在八Moshgak人口三个灌溉制度。

讨论

草药已经使用了几千年,许多世界各地的人们。由于过度捕杀的自然栖息地,保护程序的扩张,草药的评估和应用对遗传资源的保护和利用至关重要18。狭窄的遗传多样性信息和干旱胁迫的影响在Moshgak可用。在这项研究中,我们分析各种Moshgak人口来自不同地理区域的伊朗植物化学的特点和在干旱胁迫条件下植物精油成分。这项研究的结果提出了种群之间的显著差异和灌溉制度对所有特征研究,代表人口中多样性种质对干旱胁迫的响应。以前,许多重要的形态特征的变化和石油当地Moshgak人口是报道的内容2,7。这项研究的结果表明,人口显著不同的EOs和响应干旱胁迫的反应不同。EOs的增加产量以应对水赤字可能是由于碳重吸收,作为植物生长减缓,与压力相关的代谢产物的生物合成24。以前的研究报告不同的精油含量范围的叶子Moshgak人群16,25,而在目前的研究中,在果实含油量范围从1.31到3.68%。EOs内容可以受到环境因素的影响4,19。Alinian和Razmjoo26观察到干旱胁迫会导致增加孜然EO (Cuminum植物l .)。此外,Bettaieb et al。13报道说,水果孜然精油减少为了应对严重的干旱胁迫,虽然在温和的反应条件下的水赤字增加。这些发现表明,EOs取决于基因型,和干旱的水平。gc - ms分析结果表明化学Moshgak种群之间的多态性。在目前的研究中,这种植物的石油成分的水果是第一次决定。Arbabi et al。2报道称,独联体n-decanal -chrysanthyl乙酸,正癸醇,α-pinene叶子的主要化合物的数量d . anethifolia。此外,Karami Bahloli27dodecanal n-decanal报道,独联体-chrysanthenyl醋酸,α-pinene精油的主要成分d . anethifolia叶子。在目前的研究中,独联体-chrysanthyl醋酸和α-pinene是主要的化合物。看来的生产独联体醋酸-chrysanthenyl之后才开始花的外观。的生产独联体-chrysanthenyl醋酸在树叶很低,这种化合物的生物合成进行了主要在成年的花序28。抑制了前列腺素合成酶chrysanthenyl醋酸和可以镇痛属性29日。此外,醋酸chrysanthenyl可能导致antimigraine活动,由于其抑制前列腺素30.独联体-chrysanthyl醋酸是一种单萜烃2。先前的研究显示,单萜可以高度受环境压力的影响31日。环境压力条件下植物产生高浓度的单萜,因为低碳增长的分配,暗示国防和经济增长之间的竞争32

在目前的研究中,一些精油化合物的增加以应对干旱胁迫。符合我们的发现,Bettaieb et al。13报道称,灌溉制度的应用可能会导致重大的变化孜然EOs质量。例如,γ-terpinene显著增加的百分比来响应干旱胁迫。各级EO化合物的变化的水赤字是由于生物合成有关的酶的活性变化以应对水分亏缺13。单萜的形成是由萜烯合成酶催化的,由发展和激活stress-linked因素33

在这项研究中,总酚、总类黄酮含量增加水资源短缺,尽管他们追求不同的趋势,表明Moshgak干旱压力的生化反应可能是人口依赖。同样,孜然在水分胁迫下的酚醛内容大大增加13。相比之下,Neffati et al。34报道称,胡荽的多酚含量降低盐度胁迫下(c一L)。不同的趋势酚醛分数可以估算模拟干旱胁迫下的合成各种组件TPC和薄膜电路各种人群努力适应不同土壤水分水平。许多研究人员报道,TPC积累的组件和薄膜电路的模拟干旱胁迫下可以增加提高作物质量在药用植物7。哈里森和被35报道,改善TPC和薄膜电路响应干旱胁迫可以降解的结果更大的酚醛成更小的组件。此外,酚类物质有很强的抗氧化活性,因此参与宣传和活性氧(ROS)时不可避免地产生光合作用代谢是被环境限制13。本研究的结果表明,集成电路的数量是明显不同的50和对干旱胁迫的反应不同的反应。此外,抗氧化活性的变化观察根据他们的TPC drought-stressed植物。在确凿的这一发现,Bettaieb et al。36表明,孜然的抗氧化活性的变化是反映TPC变异。其他研究人员报道,孜然抗氧化活动表现出与中度和重度水分亏缺处理种子孜然的最大抗氧化活性13。暴露的植物对不良环境条件如盐度、极端干旱、或温度应力可以增强活性氧的生产13。低分子量的抗氧化剂能有效清除有害自由基11

在目前的研究中,积极的集成电路之间的关系也被观察到50在干旱条件下,TPC。相比,我们发现,先前的研究显示抗氧化活性和TPC之间显著相关c .植物13c一34

本研究的结果导致了人口高独联体-chrysanthenyl醋酸、精油从种子和总酚含量使用PCA(高PC1和PC2)。Abarkuh人口压力下被分为另一组由于p-cymene含量高。Hosseini et al。18报道称,选择基于biplot导致高的基因型鉴定茴香精油产量和主要组件。

各种关系中得到特征。正相关的化合物显示选择根据化合物的可以改进的另一个原因18。本研究的结果显示α-pinene有负相关独联体-chrysanthenyl醋酸。此外,观察之间的负相关独联体-chrysanthenyl醋酸与反式马鞭草烯醇,柠檬烯。负一些异构体化合物之间的联系也在先前的研究报告4,6

水赤字引发部分气孔关闭和大量减少了有限公司2涌入到叶子,这直接导致一个被动增加的所有进程消耗NADPH+H+降低高度,包括生物合成次生植物产品。另外,酶负责天然产物的生物合成可能积极向上调节。相应的增加NADPH+H+消费明显导致与压力相关的能量的耗散,从而防止氧自由基的生成。这些衔接手段揭示压力机制,导致增加自然产品,如类异戊二烯生物合成,生物碱,或酚醛树脂37

结论

Moshgak种质用于当前的研究表明各种植物化学的变化特征和精油的化合物可以用于进一步育种程序。总酚、总类黄酮含量,抗氧化活性和EOs增加以应对水资源短缺。同时,精油成分增加应对干旱胁迫。此外,独联体-chrysanthenyl醋酸是最丰富的精油成分这种植物的果实。最高的增加独联体-chrysanthenyl醋酸(21.32%),以应对水分胁迫是Abarkuh1观察。显著差异对干旱胁迫的反应显示,每个人可以利用某些目的。

优势种群中引入当前研究合适的候选人的进一步发展新的药用品种和制药项目。

材料和方法

道德的声明。本研究进行了按照伊斯法罕科技大学的伦理委员会,伊朗。写,知情同意是获得所有的志愿者。所有实验协议在伊斯法罕大学伦理委员会批准的技术,伊朗。

植物材料和实验条件

8个种群的种子通过药用基因银行的核心研究机构(CRF),伊斯法罕医学大学,伊斯法罕,伊朗。收集样本,允许从研究所的森林和牧场和符合当地和国家指导方针和立法。植物的决心是由Mehdi Rahimmalek博士和样本存放在伊斯法罕科技大学的植物标本,伊斯法罕,伊朗。地理位置、基因银行数量和种群的特征如表所示6和无花果。3。IUT的实验进行了研究农场,伊斯法罕(北纬32°38ʹ,东经51 39°ʹ东)。所有人群的种子种植在随机完全区组设计(RCBD三个复制和三个灌溉制度(3治疗)。每个情节包括两行2米长70厘米。土壤质地粘壤土,pH值7.5和有机质含量的1%。化肥使用100公斤N /公顷和播种前P 100公斤/公顷。艾伦灌溉制度是根据指定的程序等。37。进行了实验研究和现场研究基于先前的文献描述的标准方法18,32,36

表6 8 Moshgak种群的地理位置及其特征。
图3
图3

网站收集Moshgak人口来自不同地理区域的伊朗,在这项研究中使用。人口基因银行数量如表所示1

水分胁迫处理

在灌溉、土壤样本获得0-15土壤深度,15 - 30、30 - 45厘米三种环境来估计土壤含水量和计算灌溉用水内容基于45厘米深度加油。灌溉深度估计根据以下公式:

$ ${\文本我}{}\,= \ \离开[{\离开({{\文本{FC}} - \ uptheta} \右)/ 100}\右]{文本\ D{}} \, \ \,{\文本{B}});$ $

在那里,我是灌溉(cm)的深度,FC (−0.03 MPa)土壤水分重量在田间持水量(22%)、θ(−1.5 MPa)在灌溉时土壤重量含水率(10%),D是神经根带深度(60厘米),B是土壤容重的土壤根区(1.4 g厘米吗−3)38。种子被播种在三个实验条件,包括控制条件(C),中等(MS),严重水分胁迫(SS)。的植物灌溉播种时间,之后每周进行灌溉,直到开花初期阶段(20%开花)。此外,灌溉直到40、60和80%,分别(FC)的能力。获取应用程序每个水的体积,Parshall水槽公式(Id =我×p;p是我的分数可能会耗尽从根区)应用。最后,搞笑= (Id / Ea)×100,与Ea代表灌溉效率(%),应该是平均65%。

精油的提取

成熟的水果从d . anethifolia植物生长在每一农艺年度收集和干在树荫下4天的25°C。十个样品从每个人口的水果被混合,然后应用于测量精油。石油是通过使用Clevenger-type hydro-distillation装置4 h。为此,每个样本应用的50克。蒸馏精油是采购使用乙醚作为收集溶剂(v / v),收集玻璃容器,并将其保持在4°C到分析。最后,精油内容提出了基于风干的重量百分比种子标本。

气相色谱分析-质谱法(gc - ms)分析

精油成分进行了分析使用一个安捷伦7890气相色谱仪(美国安捷伦科技,帕洛阿尔托,CA)和一个女士HP-5毛细管柱(30 m×0.25毫米,0.25μm膜厚度)。烤箱温度设定在70°C 2分钟紧随其后的是一个增加到200°C(持有= 1分钟)。注入器和检测器温度设定在210和220°C,分别。分析是通过氦载气(进气压力17.7 psi;分裂1:10)。气相色谱仪是耦合的安捷伦5975 c质谱仪作为选择性检测器的电离电压70电动汽车;离子源、质量四极温度230,和150°C, 35 - 500 m / z的质量范围。GC -化学站操作软件应用程序。保留指数是通过使用获得的保留时间的正烷烃(C8-C21)注射前精油分析在相同的条件下,使用Kovats指数公式39。精油成分是比较确定的保留指数(RIs)与文献中报道39以及比较他们与那些记录在NIST质谱08年(国家标准与技术研究院),Wiley275。L (ChemStation数据系统)。成分的百分比获得总油成分的峰面积。

Methanolic提取和总酚含量

总酚含量与Folin-Ciocalteu试剂得到解释Hodaei et al。8有一些变化。短暂,提供样品提取、25毫升的80%甲醇添加和动摇25°C 24 h。提取过滤,和2.5毫升的Folin-Ciocalteu试剂(1:10稀释蒸馏水)和2毫升7.5%的碳酸钠添加到0.5毫升的methanolic提取和动摇。混合物在45°C加热15分钟。最后,提取液的吸光度是通过分光光度计在765海里。水和试剂应用作为一个空白样品。丹宁酸应用参考标准和总酚含量被报道为毫克的鞣酸等价物(TAE)在干燥的基础上每克样品。

总类黄酮提取的内容

氯化铝比色过程被一用于总类黄酮的测定等。40小的改变。methanolic提取(125µL)是与75年混合µL 5%的NaNO2解决方案,在黑暗中柯坪6分钟后,150µL AlCl3(10%)添加和孵化5分钟。然后,750μL氢氧化钠(1米)和2500μL蒸馏水补充道。15分钟的孵化后,提取了粉色和吸光度是读510海里。最后,交通在毫克的槲皮素等价物(QE)每克的水果干重。

DPPH自由基清除实验

的自由基清除活性d . anethifolia确定使用的程序由一et al。40小的改变。短暂,5毫升的0.1毫米与0.1毫升的甲醇溶液DPPH提取混合在选定的浓度(100,和300 ppm)。混合物被动摇,保持30分钟25°C。样品的吸光度是阅读使用紫外可见分光光度计在517 nm甲醇(80%)作为空白和叔丁基羟基hdroxytoluene(二叔丁基对甲酚)应用标准抗氧化剂在1 - 100μg /毫升比较。吸光度测量在517 nm和空白。抑制的百分比的样本估计的基础上,下面的公式:

$ $ \ mathrm抑制{\ %}= \压裂{{(}_ {B} -{一}_{一})}{{一}_ {B}} \乘以100 $ $

在那里,\ ({}_ {}\)\ ({}_ {B} \)的吸光度值百分比的DPPH自由基水果提取样本和控制,分别。随后,抑制比例对提取浓度绘制和50%的抑制浓度(IC50计算(图)的DPPH值。S1)。

统计分析

数值计算数据的正常使用proc正常测试的SAS软件9.2版本(SAS研究所卡里、数控、美国)41。方差分析(方差分析)是使用一个随机区组设计执行2年灌溉和三个层次,用SAS软件9.2版本的GLM过程(SAS研究所卡里、数控、美国)37。意思是比较用最少的显著差异(LSD)测试使用SAS软件发布9.2 (p < 0.05) (SAS研究所卡里、数控、美国)37

聚类分析是根据病房的方法使用的平方欧氏通过统计图形软件版本不同。17.2 (http://www.statgraphics.com)。主成分分析(PCA)进行识别之间的关联d . anethifolia人口和所有测量特征使用统计图形软件版本。17.2 (http://www.statgraphics.com)。