基于细胞色素b基因的中国岩羊不同地理种群遗传差异分析

为了揭示中国岩羊不同地理种群的遗传差异,探讨岩羊亚种分化的分子机制,采用中国岩羊不同地理种群的细胞色素b(Cyt b)基因的全序列,分析了碱基变异情况、遗传距离以及核苷酸序列差异。用最大似然法和贝叶斯法构建分子系统树并对获得的拓扑结构进行分析。结果发现,西藏亚种与四川亚种Cyt b基因平均序列差异为4.2%(±0.007),处于偶蹄目亚种的序列差异范围内,支持了目前对岩羊西藏亚种的分类地位。四川亚种内部各地理种群之间的遗传距离(0.033±0.0 111)与它们分别到西藏种群的遗传距离(0.042±0.007)差异不显著(t=1.824,P=0.084),说明四川亚种内部各地理种群间已经发生较显著的遗传分化。其中,四川、甘肃和青海种群亲缘关系较近,并与四川亚种内部的其它种群已产生了显著的遗传分化。因此认为四川亚种内部各地理种群的种下分化需要深入研究。     

四纹豆象不同地理种群的遗传分化

【目的】通过对四纹豆象Callosobruchus maculatus Fabricius不同地理种群mtDNA-Cytb和COⅠ基因部分序列进行比较, 分析其不同地理种群间的遗传分化情况, 为揭示其与生物入侵的关系及入侵过程中种群系统发育地理格局与演变机制提供依据。【方法】用PCR产物直接测序法对分别来自中国海南、喀麦隆、韩国和泰国的四纹豆象4个地理种群的mtDNA-Cytb和COⅠ序列进行测序, 运用软件MEGA3.1对四纹豆象不同地理种群mtDNA-Cytb和COⅠ序列进行序列分析, 以绿豆象C. chinensis为外群构建了不同单倍型的分子系统树。【结果】34条420 bp Cytb序列中共检测到14个多态位点和5种单倍型, 33条822 bp COⅠ序列中检测到28个多态位点和9种单倍型, 其中4种单倍型为独享单倍型, 其余为全部或部分种群的共享单倍型。AMOVA分析结果显示, 四纹豆象4个地理种群间的遗传结构差异并不明显, 遗传差异主要发生在地理种群内。对4个地理种群进行了Fst值和基因流动统计, 结果表明4个地理种群间既存在着一定数量的基因交流, 也存在一定程度的遗传分化。【结论】根据单倍型分布格局初步推测, 中国不可能是四纹豆象的原产地, 而喀麦隆有可能是原产地之一, 并且喀麦隆种群与泰国种群之间的基因交流比较充分, 而中国种群与其他种群之间的遗传分化相对较大。     

福建东张水库陆封型香鱼的遗传变异

测定了福建东张水库18尾陆封型香鱼的ND4-tRNA^Ser基因533 bp序列,发现2个可变位点和2个简约信息位点,共检测到4种单倍型。基于Kimura 2-parameter模型的遗传距离在0.001 9～0.003 9之间,在NJ系统树上无明显谱系结构;单倍型多样性(Hd)和核苷酸多样性(Pi)分别为0.627和0.001 37,平均核苷酸差异数为0.732,远低于两侧洄游型日本香鱼指名亚种(Hd为0.793,Pi为0.002 94),与被列为IUCN濒危物种名录的日本琉球香鱼亚种(Hd=0.519,Pi=0.001 11)相当。东张水库香鱼遗传多样性如此贫乏到底是历史原因还是近期陆封的结果,仍需进一步研究。目前亟待采取有效措施保护东张水库陆封香鱼。     

中国舞毒蛾六个地理种群的RAPD分析及SCAR标记构建

舞毒蛾Lymantria dispar L.是世界性农林害虫, 包含不同的亚种, 其中亚洲舞毒蛾的雌蛾具有较强的飞行能力, 已成为国际性的重要检疫性有害生物。然而, 不同舞毒蛾亚种及种群间形态难辨, 因此采用传统的手段鉴别舞毒蛾亚种种群是很困难的。本研究首先采用RAPD标记分析了中国舞毒蛾6个地理种群的遗传多态性。结果表明, 所检测的舞毒蛾种群的遗传分化系数G_st为0.7571, 由此推算出的平均有效迁移数（基因流参数）N_em为0.1604, 说明不同舞毒蛾种群间的遗传分化程度较高, 缺乏广泛的基因流动。本研究在RAPD遗传分析基础之上, 筛选出了4个舞毒蛾种群的特异性遗传位点, 然后对这些特异性位点进行了克隆测序、 序列分析和位点特异性引物设计。结果表明, 其中2个舞毒蛾种群的位点特异性引物可产生序列特征性扩增区域（SCAR）标记。经验证, 这些标记可被用来鉴别特定的舞毒蛾地理种群, 因此有助于对这些舞毒蛾地理种群的分布与扩散进行监测。     

大石鸡兰州亚种的mtDNA种群遗传结构、单倍型分布和遗传多样性

大石鸡(Alectorismagna)是中国西北部的特有种。我们测定了大石鸡兰州亚种(A.m.lanzhouensis)8个地理种群106个样本的线粒体DNA控制区5′端458bp序列,研究其种群遗传结构和遗传多样性。27个变异位点共确定25种单倍型,其中单倍型M2广泛分布,而许多单倍型为一些地方种群特有。单倍型分布沿着南北方向变化,存在明显的地理结构。8个种群中核苷酸多样性最高的是定西种群,0.0069,最低的是海原种群,0.0028;单倍型多样性最高的是武山种群,0.86,最低的是北道种群,0.52。北方种群比南方种群具有更高的遗传多样性。系统发生树和单倍型分布表明,大石鸡兰州亚种存在两个明显的分支。溯祖理论、更新世冰期和花粉支持兰州亚种起源于兰州盆地,这个盆地是其遗传多样性的中心。     

中国迁粉蝶属(鳞翅目,粉蝶科)基于NDⅠ和COⅠ基因部分序列的系统发育关系

测定了中国迁粉蝶属3种4亚种以及迁粉蝶指名亚种5个生态型的线粒体COⅠ基因和NDⅠ基因部分序列,以灵奇尖粉蝶为外群,分别采用最大简约法、最大似然法和贝叶斯法构建分子系统树,探讨了中国迁粉蝶属3种4亚种间的系统发育关系,并首次明确了迁粉蝶指名亚种5个生态型的亲缘关系。结果显示在测得的COⅠ基因661bp序列和NDⅠ基因472bp序列中,有变异位点227个,简约位点119个,A+T含量(75.5%)明显偏高。迁粉蝶属4亚种分子系统关系为:(镉黄迁粉蝶指名亚种+((梨花迁粉蝶指名亚种+梨花迁粉蝶海南亚种)+迁粉蝶指名亚种)),迁粉蝶指名亚种5个生态型分子系统关系为:(无纹型+((有纹型+红角型)+(银斑型+血斑型)))。     

线粒体DNA序列变异显示鸭绿江茴鱼为黑龙江茴鱼同物异名

经测序获得北极茴鱼鸭绿江亚种Thymallus arcticus yaluensis和黑龙江茴鱼T.grubii共26尾的线粒体DNA控制区基因全序列,以下游黑龙江茴鱼T.tugarinae、北极茴鱼指名亚种T.a.arcticus和欧洲茴鱼T.thymallus作外群,构建了茴鱼属部分种分子系统发育树,以探讨北极茴鱼鸭绿江亚种分类地位的有效性。结果显示,鸭绿江茴鱼与黑龙江茴鱼的4个地理种群(呼玛河、甘河、牡丹江、松花江)间的Kimura 2-parameter遗传距离仅为0·007~0·015,处于黑龙江茴鱼不同地理种群间的变异范围之内(0·004~0·015),明显小于属内种间的遗传距离(0·046~0·062)。以邻接法(NJ)和最大简约法(MP)构建的茴鱼属部分种系统发育树表明,鸭绿江茴鱼未构成独立的单系群,也未与北极茴鱼T.a.arcticus形成亚种关系,而是聚到黑龙江茴鱼种群内部,共同构为一个单系群。研究表明,鸭绿江茴鱼T.a.yaluensis并非有效种,应为黑龙江茴鱼T.grubii的同物异名,可视为黑龙江茴鱼的一个地理生态种群。     

不同地理种群美洲斑潜蝇及近缘种的rDNA-ITS1序列分析和比较

【目的】通过对美洲斑潜蝇Liriomyza sativae Blanchard 不同地理种群及近缘种间的核糖体DNA第一内转录间隔区（rDNA-ITS1）进行比较,分析美洲斑潜蝇不同地理种群间的遗传分化情况,并为美洲斑潜蝇与近缘种间提供分子鉴别标记。【方法】用PCR产物直接测序法及克隆测序法对我国美洲斑潜蝇8个地理种群的rDNA-ITS1序列进行测序,并调用GenBank中3个近缘种的rDNA-ITS序列,运用软件MEGA3.1对美洲斑潜蝇不同地理种群及近缘种间的rDNA-ITS1序列进行分析。【结果】美洲斑潜蝇8个地理种群间的分化程度较低,只有8个变异位点,遗传距离都在0.02以下,但4个近缘种间的碱基差异显著,遗传距离为0.149～0.390,有126个变异位点,12个美洲斑潜蝇特异性识别位点。【结论】虽然基于rDNA-ITS1序列所显示的美洲斑潜蝇各地理种群之间的遗传分化很小,但是其分化趋势与地理分布基本相吻合;得到的12个特异性识别位点不仅可以作为美洲斑潜蝇与其近缘种间鉴别的分子标记,而且可为今后设计鉴别性PCR引物提供重要的参考依据。     

暗褐蝈螽不同地理种群间的遗传分化

利用线粒体COI基因片段研究了我国吉林、辽宁、内蒙古、河北和四川的暗褐蝈螽(Gampsocleis sedakovii)12个地理种群间的遗传分化。结果表明: 36条626 bp的mtDNA-COI基因片段中共检测到单倍型29种, 多态位点71个, 其中简约信号位点37个, 单一多态位点34个。分子变异等级分析(AMOVA)的计算结果显示, 群体内变异仅占总变异的37.23%, 明显小于群体间变异, Fst值为0.62770, 各群体间呈现明显的遗传分化。最大简约法(MP)系统发育分析结果显示, 暗褐蝈螽的12个地理种群以极高的自举支持度(100%)聚类, 形成两个独立的分枝。由于两个分枝的聚类结果与基于形态学特征划分的两个亚种(Gampsocleis sedakovii sedakovii和Gampsocleis sedakovii obscura)并不对应。综合采集地的生境, 初步推测暗褐蝈螽两个亚种间的形态学差异可能是由于它们所处生境不同所引起。在12个地理种群中, 只有内蒙古通辽(NTL)的个体在两个分枝中均有出现。单倍型分布格局研究发现, 单倍型H10是内蒙古通辽(NTL)、鄂温克(NEWK)和吉林吉林(JJL)3个地理种群的共享单倍型, 说明它们来自于共同的祖先。研究结果支持我国的东北地区是暗褐蝈螽遗传分化中心的观点, 但不支持基于形态学特征划分的两个亚种。     

高原裸裂尻鱼头部轮廓形状及其变异的几何形态测量分析

为了明确描述高原裸裂尻鱼的表型多样性,运用现代几何形态测量法,基于高原裸裂尻鱼头部轮廓的形态坐标数据,对其5个地理种群的样本做了定量分析。结果主要如下:①主成分分析（PCA）和典型变量分析（CVA）显示,5个地理种群的样本的变异主要位于口型特征、下吻部形状、下颌骨形状、下颌锐利角质程度和头部形状等性状。②基于主成分分析、典型变量分析和聚类分析,高原裸裂尻鱼三个亚种（班公湖亚种、玛法木湖亚种和指名亚种）相互分开,在PC1轴上从负值到正值依次排列,从而确定了高原裸裂尻鱼亚种划分的有效性。③通过形态差异的显著性检验,基于马氏距离和普氏距离分析了高原裸裂尻鱼各地理种群之间的形态变异距离,根据差异显著性将该5个地理种群的高原裸裂尻鱼分为湖泊型亚种和河流型亚种,随后的聚类分析同样支持该结果。总之,通过几何形态测量法,准确描述了高原裸裂尻鱼头部轮廓的形态多样性,按照外形的差异,将5个地理种群分为两种生态型,并推测该差异主要是由于所生活的水体环境和食物差异所造成。对高原裸裂尻鱼形态差异的识别和研究,能为其生态学研究和物种多样性保护提供相应指导,丰富微进化理论。     

Mitochondrial phylogeography of the Dalmatian wall lizard, Podarcis melisellensis (Lacertidae)

A 903 bp section of the mitochondrial cytochrome b gene was sequenced from 73 specimens of Podarcis melisellensis collected at 52 localities distributed over the major part of the species’ range. In addition, parts of the 12S (about 470 bp) and 16S rRNA (about 500 bp) genes were analysed for 11 representative samples leading to a congruent phylogeny. Our study includes representatives of all 20 subspecies recognized today. The phylogenetic analysis of the sequence data revealed three main clades: mainland with nearby islands, Vis archipelago, and Lastovo archipelago. The degree of mitochondrial DNA divergence among these clades suggests a separation of the respective population groups during the earliest Pleistocene. The phylogenetic pattern observed within the species is in sharp contrast to the actual taxonomic division into subspecies. A correlation between genetic diversity of P. melisellensis populations and paleogeography of the regions they inhabit is discussed.     

Molecular evidence for the subspecific differentiation of blue sheep (Pseudois nayaur) and polyphyletic origin of dwarf blue sheep (Pseudois schaeferi)

Blue sheep (Pseudois nayaur), a Central Asian ungulate with restricted geographic distribution, exhibits unclear variation in morphology and phylogeographic structure. The composition of species and subspecies in the genus Pseudois is controversial, particularly with respect to the taxonomic designation of geographically restricted populations. Here, 26 specimens including 5 dwarf blue sheep (Pseudois schaeferi), which were collected from a broad geographic region in China, were analyzed for 2 mitochondrial DNA fragments (cytochrome b and control region sequences). In a pattern consistent with geographically defined subspecies, we found three deeply divergent mitochondrial lineages restricted to different geographic regions. The currently designated two subspecies of blue sheep, Pseudois nayaur nayaur and Pseudois nayaur szechuanensis, were recognized in the phylogenetic trees. In addition, the Helan Mountain population showed distinct genetic characteristics from other geographic populations, and thus should be classified as a new subspecies. In contrast, dwarf blue sheep clustered closely with some blue sheep from Sichuan Province in the phylogenetic trees. Therefore, dwarf blue sheep appear to be a subset of Pseudois nayaur szechuanensis. After considering both population genetic information and molecular clock analysis, we obtained some relevant molecular phylogeographic information concerning the historical biogeography of blue sheep. These results also indicate that western Sichuan was a potential refugium for blue sheep during the Quaternary period.     

韩国与中国东北地区松鼠线粒体DNA控制区的序列变异

为了研究松鼠东北亚种(Sciurus vulgaris manchuricusThomas)不同种群的序列变异水平并进一步确定分类地位,我们分析了韩国5个地点和中国东北2个地点的松鼠标本的线粒体DNA控制区的全序列(1 058 bp)。39个韩国松鼠标本显示出21种单倍型,这些单倍型间的平均Tamura-Nei距离为1·0%; 24个中国松鼠标本显示21种单倍型,单倍型间的平均Tamura-Nei距离为1·4% (1 058 bp的全序列中发生变异的位点有119个,占11·2%)。韩国松鼠和中国松鼠间的平均距离为1·3%。并且韩国和中国松鼠的所有42个单倍型形成了一个单系分支,Fst值为0·04,表明在两个国家的松鼠间没有发生遗传分化。因此,序列分析的分子生物学的结果支持现行的分类,即来自韩国的朝鲜亚种(S·v·coreae)是中国北部地区松鼠东北亚种(S·v·manchuricus)的同物异名。这还需要进一步对北朝鲜和中国东北其它地区更多标本的分子和形态学分析来验证这一结论。     

Population genetic structure of the melon fly, <Emphasis Type="Italic">Bactrocera cucurbitae</Emphasis> (Diptera: Tephritidae), from China and Southeast Asia

The melon fly, Bactrocera cucurbitae Coquillett, is a species of fruit flies of significant agricultural interest. Of supposed Indian origin, the melon fly is now widely distributed throughout South East Asia up to China, while it has been recently eradicated from Japan. The population structure of seven geographic populations from coastal China, as well as samples from other regions of South East Asia and Japan, including lab colonies, have been studied using a 782 bp fragment of mitochondrial cytochrome oxidase I (COI) gene sequence. The observed genetic diversity was exceedingly low, considering the geographic scale of the sampling, and one single haplotype was found to be predominant from Sri Lanka to China. We confirm that Bactrocera cucurbitae exists in South East Asia as a single phyletic lineage, that Chinese populations are genetically uniform, and that no apparent genetic differentiation exists between these and three available Japanese melon fly sequences.     

基于线粒体基因组控制区序列的中国冰清绢蝶遗传分化和谱系生物地理研究

【目的】分析中国冰清绢蝶Parnassius glacialis遗传多样性和遗传分化情况及其系统发生关系,推测其起源及分歧时间,并探讨其历史生物地理分布格局的成因。【方法】以2008-2016年采集于中国13个地区的冰清绢蝶的325头及绢蝶属其他11个种的11头个体样品,对其线粒体基因组的AT富集区(控制区)序列进行PCR扩增和序列测定,采用MEGA 6.0, Dna SP 5.1和Arlequin 3.5等软件分析其遗传多样性和遗传分化情况;以其他近缘绢蝶种类作为外类群,采用PhyML3.0, MrBayes 3.2和BEAST V1.8.3等软件重建冰清绢蝶的系统发生树,并利用宽松分子钟和校准方法推测中国冰清绢蝶的起源和分歧时间;同时结合现今冰清绢蝶的地理分布特点和第四纪以来的地球环境事件,分析其扩散路径,探讨其历史生物地理分布格局及成因。【结果】冰清绢蝶AT富集区序列长度介于487～495 bp之间,平均为491 bp,其长度的差异主要体现在Poly(T)或Poly(A)的数目差异,A+T平均含量为95.76%。供试冰清绢蝶13个地理居群325头个体中,共检测出基于线粒体基因组AT富集区序列的239个单倍型,单倍型多样性Hd为0.9971;核苷酸多态性指数π值为0.02948, Theta (per site)Eta值为0.06594。系统发生分析和分子定年显示,冰清绢蝶和白绢蝶P.stubbendorfii为近缘姊妹种,二者在距今7.52百万年(Ma)的中新世晚期开始分化。冰清绢蝶的祖先可能起源于我国西南地区(现今横断山脉-喜马拉雅附近的青藏高原东北缘)距今1.53 Ma的更新世时期,先由起源地先期扩散到达小陇山、秦岭一带;后来,由于第四纪的冰期-间冰期轮回事件和扩散路径的不同而分化出2大支系,但两个支系之间存在少量的单倍型相互混杂现象;随后,2大支系继续向低海拔、低纬度地区扩散,即通过伏牛山脉向东北方向扩散到嵩山、泰山、昆嵛山一带,经大别山系向东南方向扩散到琅琊山、紫金山和天目山一带。【结论】中国13个居群的冰清绢蝶没有形成明显的地理分化格局,其遗传变异主要来自居群内部,居群间的遗传分化现象尚不明显;冰清绢蝶起源于更新世时期的较高海拔地区,由于第四纪的冰期-间冰期轮回事件以及扩散路径的不同而形成2大支系并向中国东部和南部低海拔和低纬度地区扩散。     

新疆草兔的种群遗传结构和亚种分化

新疆草兔 (Lepus capensis) 的群体遗传结构至今无系统的研究报道,亚种水平的分类也长期存在争议.该文测定了形态分类上的新疆草兔3个亚种共87个个体的线粒体DNA (mtDNA)控制区(control region,D-Loop)592 bp的序列,经分析发现148个多态性位点,共定义了44个单倍型.新疆草兔的单倍型多样度(h,0.977 ± 0.005)和核苷酸多样度(π,0.064 ± 0.031)都较高,显示了较高的遗传多样性.分子变异分析(AMOVA)结果显示,4个地理群体间的显著分化可能是由地理隔离造成的.群体遗传结构分析显示,新疆草兔包含4个进化枝,并且每个进化枝都对应特定的分布区域,显示了明显的系统地理结构.该研究的结果支持形态分类上草兔西域亚种(L.c.lehmanni)的分类地位; 但中亚亚种(L.c.centrasiaticus)被分为两个独立的进化枝,提示可能存在两个亚种; 帕米尔亚种(L.c.pamirensis)与其他亚种间的遗传距离在13%以上,提示其可能已达到种的分化水平.     

Large sequence divergence of mitochondrial DNA genotypes of the control region within populations of the african antelope, kob (Kobus kob)

The genetic differentiation among kob populations (Kobus kob) representing two recognized subspecies was examined using mitochondrial control region sequences. Two distinct lineages (estimated sequence divergence of 9.8%) exhibited different geographical distributions and do not coincide with previously recognized ranges of subspecies. The presence of the two lineages was further supported with sequences of mitochondrial cytochrome b gene. One lineage was predominant in the west and southern ranges of the populations studied and the other was commonly found in a more northern distribution (Murchison populations) in Uganda. Murchison and the geographically intermediate Toro populations (Uganda) represented the area of overlap. The existence of the two lineages in the area of overlap is hypothesized to have resulted from a range expansion and secondary contact of the two lineages of kob that evolved in allopatry. The existence of the kob during the Pleistocene offers a plausible explanation for the observed biogeographic pattern. Our mitochondrial data reveal two examples of discordance between a gene tree and presumed species tree as: (i) the two lineages co-occur in the kob subspecies, Kobus kob thomasi (Uganda kob); and (ii) the puku, which was included in the analysis because of its controversial taxonomic status (currently recognized as a distinct species from the kob), is paraphyletic with respect to the kob. Significant degrees of heterogeneity were detected between populations. Relatively high genetic variation was observed in the populations, however, the inclusion of distinct lineages influences the population structure and nucleotide diversity of the kob populations.     

Using dense locality sampling resolves the subtle genetic population structure of the dispersive fish species Plecoglossus altivelis

In dispersive species with continuous distributions, genetic differentiation between local populations is often absent or subtle and thus difficult to detect. To incorporate such subtle differentiation into management plans, it may be essential to analyse many samples from many localities using adequate numbers of high‐resolution genetic markers. Here, we evaluated the usefulness of dense locality sampling in resolving genetic population structure in the ayu (Plecoglossus altivelis), a dispersive fish important in Japanese inland fisheries. Genetic variability in, and differentiation between, ayu populations around the Japan–Ryukyu Archipelago were investigated in 4746 individuals collected from 120 localities by genotyping 12 microsatellite markers. These individuals represented the two subspecies of ayu, namely the Ryukyuan subspecies (Plecoglossus altivelis ryukyuensis) and both amphidromous and landlocked forms of the nominotypical subspecies (P. a. altivelis) along the archipelago. We successfully detected an absence of genetic differentiation within the landlocked form and subtle but significant differentiation and clear geographic patterns of genetic variation among populations of the amphidromous form, which had been considered genetically homogeneous. This suggests that dense locality sampling effectively resolves subtle differences in genetic population structure, reducing stochastic deviation in the detection of genetic differentiation and geographic patterns in local populations of this dispersive species. Resampling analyses based on empirical data sets clearly demonstrate the effectiveness of increasing the number of locality samples for stable and reliable estimations of genetic fixation indices. The genetic population structure observed within the amphidromous form provides useful information for identifying management or conservation units in ayu.     

Investigating Concordance among Genetic Data, Subspecies Circumscriptions and Hostplant Use in the Nymphalid Butterfly Polygonia faunus

Subspecies are commonly used taxonomic units to formally describe intraspecific geographic variation in morphological traits. However, the concept of subspecies is not clearly defined, and there is little agreement about what they represent in terms of evolutionary units, and whether they can be used as reliably useful units in conservation, evolutionary theory and taxonomy. We here investigate whether the morphologically well-characterized subspecies in the North American butterfly Polygonia faunus are supported by genetic data from mitochondrial sequences and eight microsatellite loci. We also investigate the phylogeographic structure of P. faunus and test whether similarities in host-plant use among populations are related to genetic similarity. Neither the nuclear nor the mitochondrial data corroborated subspecies groupings. We found three well defined genetic clusters corresponding to California, Arizona and (New Mexico+Colorado). There was little structuring among the remaining populations, probably due to gene flow across populations. We found no support for the hypothesis that similarities in host use are related to genetic proximity. The results indicate that the species underwent a recent rapid expansion, probably from two glacial refugia in western North America. The mitochondrial haplotype network indicates at least two independent expansion phases into eastern North America. Our results clearly demonstrate that subspecies in P. faunus do not conform to the structuring of genetic variation. More studies on insects and other invertebrates are needed to better understand the scope of this phenomenon. The results of this study will be crucial in designing further experiments to understand the evolution of hostplant utilization in this species.