基于COⅠ基因的我国悬铃木方翅网蝽(半翅目: 网蝽科)种群遗传多样性和遗传结构分析

doi:10.11707/j.1001-7488.20210810

摘要/Abstract

摘要：

目的: 研究悬铃木方翅网蝽在我国不同地理种群间的遗传多样性和遗传结构，探讨其种群在我国的入侵趋势。方法: 共测定12个悬铃木方翅网蝽地理种群中240个个体COⅠ基因的部分核苷酸序列，利用DnaSP V6.12.03等软件对方翅网蝽不同地理种群间的遗传分化程度、基因流水平、分子变异情况、种群历史动态等进行分析。结果: 悬铃木方翅网蝽12个地理种群的COⅠ核苷酸序列含有13个变异位点和9个单倍型，其中有5个共享单倍型，Hap1为祖先单倍型。方翅网蝽总种群单倍型多样性指数H_d为0.584，地理种群间H_d值在0~0.758之间。种群总固定系数F_ST值为0.48，各地理种群间F_ST值在-0.05~0.89之间，种群间基因分化程度差异较大；种群总基因流N_m值为0.55，且各地理种群间N_m值在-85.25~47.58之间，种群间基因交流程度不同。悬铃木方翅网蝽种群内个体间变异占总变异的52.19%，同时华中组群经历了方翅网蝽种群扩张事件。结论: 12个悬铃木方翅网蝽地理种群在不同组群间和地理种群间具有显著分化的遗传结构。悬铃木方翅网蝽的遗传分化主要来自种群内个体间变异，各地理种群间遗传距离与地理距离无明显的相关性，地理距离未明显阻碍地理种群间的基因交流。同时，推测在中国主要分布区内悬铃木方翅网蝽以上海为中心向全国进行扩散。

关键词: 悬铃木方翅网蝽, COⅠ基因, 固定系数, 遗传分化, 基因流

Abstract:

Objective: Corythucha ciliata(Hemiptera: Tingidae), an important invasive pest in China, mainly damages on Platanus plants. It has broken out in multiple regions of China. This research tries to explore the genetic differentiation among the geographic population of C.ciliata, and probe its spread path. Method: The sectional nucleotide sequences of COⅠ gene of 240 individuals from 12 geographic populations were determined. Based on their nucleotide sequences, the genetic differentiation, gene flow level, molecular variance, and population demographic history were analyzed by using DnaSP V6.12.03 and other analysis software. Result: There were 13 mutation sites and 9 haplotypes in the COⅠ gene among 12 geographic populations. Among them, 5 haplotypes were shared by all the geographic populations, and the haplotype H1 was an ancestor haplotype. Haplotype diversity index (H_d) of the total population was 0.584, and the H_d value among geographic populations ranged from 0 to 0.758. The fixed coefficient(F_ST) of the total population was 0.48, and the value between geographic populations ranged from-0.05 to 0.89, suggesting that the degree of differentiation between different populations of C. ciliata was obviously different. Gene flow value(N_m) was 0.55 within the total population and ranged from-85.25 to 47.58 between geographic populations, indicating that the degree of gene exchange was distinct difference among populations. The individual variation within the populations accounted for 52.19% of the total variation, while populations of C. ciliata from the Central China experienced a population expansion event. Conclusion: There are significant differentiation of genetic structures among different groups and different geographical populations of C.ciliata. The genetic differentiation of geographic populations primarily derives from variation of intraspecies. The genetic distance of the geographical populations is no significant association with geographical distance, and the geographical distance has no obvious effect on the gene flow of populations. Furthermore, we speculate that the insect spread from Shanghai as the center to other regions in China.

Key words: Corythucha ciliata, mt DNA COⅠ gene, fixed coefficient, genetic differentiation, gene flows

中图分类号:

S763.35

张元臣,卢绍辉,龚东风,马兴莉. 基于COⅠ基因的我国悬铃木方翅网蝽(半翅目: 网蝽科)种群遗传多样性和遗传结构分析[J]. 林业科学, 2021, 57(8): 102-111.

Yuanchen Zhang,Shaohui Lu,Dongfeng Gong,Xingli Ma. Analysis of Genetic Diversity and Genetic Structure in Geographic Populations of Corythucha ciliata(Hemiptera: Tingidae) from China Based on Mitochondrial DNA COⅠ Gene Sequences[J]. Scientia Silvae Sinicae, 2021, 57(8): 102-111.

图/表 11

表1

图1

表2

图2

表3

表4

表5

图3

表6

图4

表7

参考文献 0

	宾淑英, 吴仲真, 张鹤, 等. 分子标记技术在昆虫入侵研究中的应用. 昆虫学报, 2014, 57 (9): 1094- 1104.
	Bin S Y , Wu Z Z , Zhang H , et al. Application of molecular markers in studies of insect invasion. Acta Entomologica Sinica, 2014, 57 (9): 1094- 1104.
	陈圣宾, 欧阳志云, 徐卫华, 等. Beta多样性研究进展. 生物多样性, 2010, 18 (4): 323- 335.
	Chen S B , Ouyang Z Y , Xu W H , et al. A review of beta diversity studies. Biodiversity Science, 2010, 18 (4): 323- 335.
	付宁宁, 刘佳, 渠成, 等. 悬铃木方翅网蝽化学感受蛋白ccilcsp1的结构及其结合寄主挥发物的预测分析. 林业科学, 2017, 53 (10): 109- 117. doi: 10.11707/j.1001-7488.20171012
	Fu N N , Liu J , Qu C , et al. Analysis of Corythucha ciliata CcilCSP1 structure and prediction of its binding to host-plant volatiles. Scientia Silvae Sinicae, 2017, 53 (10): 109- 117. doi: 10.11707/j.1001-7488.20171012
	郝德君, 陆佳伟, 殷丽. 3种生物农药对悬铃木方翅网蝽的室内毒力及解毒酶活性的影响. 林业科学, 2012, 48 (7): 86- 91.
	Hao D J , Lu J W , Yin L . Toxicity of three biopesticides and their effect on detoxification enzyme activity of Corythucha ciliata. Scientia Silvae Sinicae, 2012, 48 (7): 86- 91.
	呼晓庆, 杨兆福. 基于线粒体COI, Cytb和COⅡ基因的中国草地螟不同地理种群遗传分化分析. 昆虫学报, 2019, 62 (6): 720- 733.
	Hu X Q , Yang Z F . Analysis of the genetic differentiation among different geographic populations of Loxostege sticticalis(Lepidoptera: Crambidae) in China based on mitochondrial COI, Cytb and COII genes. Acta Entomologica Sinica, 2019, 62 (6): 720- 733.
	鞠瑞亭, 李博. 悬铃木方翅网蝽: 一种正在迅速扩张的城市外来入侵害虫. 生物多样性, 2011, 18 (6): 638- 646.
	Ju R T , Li B . Sycamore lace bug, Corythucha ciliata, an invasive alien pest rapidly spreading in urban China. Biodiversity Science, 2011, 18 (6): 638- 646.
	李传仁, 夏文胜, 王福莲. 悬铃木方翅网蝽在中国的首次发现. 动物分类学, 2007, 32 (4): 944- 946. doi: 10.3969/j.issn.1000-0739.2007.04.038
	Li C R , Xia W S , Wang F L . First records of Corythucha ciliata(Say) in China(Hemiptera, Tingidae). Acta Zootaxonomica Sinica, 2007, 32 (4): 944- 946. doi: 10.3969/j.issn.1000-0739.2007.04.038
	李慧, 周顺武, 王亚非. 西太平洋副热带高压异常与中国长江中下游夏季降水关系研究综述. 气象与环境学报, 2013, 29 (1): 93- 102. doi: 10.3969/j.issn.1673-503X.2013.01.016
	Li H , Zhou S W , Wang Y F . A review on relationship between subtropical high anomaly over West Pacific and summer precipitation in the middle-lower reaches of the Yangtze river. Journal of Meteorology and Environment, 2013, 29 (1): 93- 102. doi: 10.3969/j.issn.1673-503X.2013.01.016
	李峰奇, 付宁宁, 张连忠, 等. 悬铃木方翅网蝽生物学、化学生态学及防治研究进展. 昆虫学报, 2018, 61 (9): 1076- 1086.
	Li F Q , Fu N N , Zhang L Z , et al. Advances in biology, chemical ecology and control of the sycamore lace bug, Corythucha ciliata(Hemiptera: Tingidae). Acta Entomologica Sinica, 2018, 61 (9): 1076- 1086.
	李菁, 张颖, 王振营, 等. 基于线粒体DNA COⅡ基因的亚洲玉米螟中国不同地理种群遗传分化及基因流研究. 昆虫学报, 2010, 53 (10): 1135- 1143.
	Li J , Zhang Y , Wang Z Y , et al. Genetic differentiation and gene flow among different geographical populations of the Asian corn borer, Ostrinia furnacalis(Guenée)(Lepidoptera: Crambidae) in China estimated by mitochondrial COⅡ gene sequences. Acta Entomologica Sinica, 2010, 53 (10): 1135- 1143.
	刘蕾, 周晶, 刘俊杰, 等. 长江中下游地区持续性暴雨年代际变化特征及环流形势. 气象与环境学报, 2018, 34 (3): 30- 38.
	Liu L , Zhou J , Liu J J , et al. Decadal characteristics of persistent heavy precipitation in the middle and lower reaches of the Yangtze River and its circulation pattern. Journal of Meteorology and Environment, 2018, 34 (3): 30- 38.
	卢绍辉. 2020. 悬铃木方翅网蝽种群遗传结构及扩散传播机制研究. 长沙: 中南林业科技大学博士学位论文.
	Lu S H. 2020. Study on the population genetic structure and diffusion mechanism of Sycamore lace bug, Corythucha ciliata(Say). Changsha: PhD thesis of Central South University of Forestry and Technology. [in Chinese]
	卢绍辉, 魏美才, 袁国军, 等. 悬铃木方翅网蝽飞行能力测试及远距离扩散轨迹分析. 昆虫学报, 2018, 61 (9): 97- 105.
	Lu S H , Wei M C , Yuan G J , et al. Flight capacity test and trajectory analysis of long-distance dispersal of Corythucha ciliata (Hemiptera: Tingidae). Acta Entomologica Sinica, 2018, 69 (9): 97- 105.
	王福莲, 李传仁, 刘万学, 等. 新入侵物种悬铃木方翅网蝽的生物学特性与防治技术研究进展. 林业科学, 2008, 44 (6): 137- 142. doi: 10.3321/j.issn:1001-7488.2008.06.023
	Wang F L , Li C R , Liu W X , et al. Advance in biological characteristics and control techniques of the new invasive sycamore lace bug(Corythucha ciliata). Scientia Silvae Sinicae, 2008, 44 (6): 137- 142. doi: 10.3321/j.issn:1001-7488.2008.06.023
	王克勤, 刘兴龙, 邵天玉. 中国东北地区亚洲玉米螟遗传多样性及寄主专化性分析. 昆虫学报, 2018, 61 (9): 1054- 1066.
	Wang K Q , Liu X L , Shao T Y . Analysis of the genetic diversity and host specialization of Ostrinia furnacalis(Lepidoptera: Crambidae) in Northeast China. Acta Entomologica Sinica, 2018, 61 (9): 1054- 1066.
	王蒙, 徐浪, 张润杰, 等. 基于线粒体COI基因的桔小实蝇种群遗传分化研究. 昆虫学报, 2014, 57 (12): 1424- 1438.
	Wang M , Xu L , Zhang R J , et al. Population genetic differentiation analysis of the oriental fruit fly, Bactrocera dorsalis(Diptera: Tephritidae) based on mtDNA COI gene. Acta Entomologica Sinica, 2014, 57 (12): 1424- 1438.
	王戎疆, 万宏, 龙玉, 等. 利用线粒体COII基因序列对中国尾蛱蝶属系统分化的研究(鳞翅目: 蛱蝶科). 昆虫学报, 2004, 47 (2): 243- 247. doi: 10.3321/j.issn:0454-6296.2004.02.018
	Wang R J , Wan H , Long Y , et al. Phylogenetic analysis of Polyura in China inferred from mitochondrial COⅡ sequences(Lepidoptera: Nymphalidae). Acta Entomologica Sinica, 2004, 47 (2): 243- 247. doi: 10.3321/j.issn:0454-6296.2004.02.018
	肖娱玉, 王凤, 鞠瑞亭, 等. 上海地区悬铃木方翅网蝽的生活史及发生情况. 昆虫知识, 2010, 47 (2): 404- 408.
	Xiao W Y , Wang F , Ju R T . Life history and occurrence of Corythucha ciliata in Shanghai. Chinese Bulletin of Entomology, 2010, 47 (2): 404- 408.
	许有鹏, 于瑞宏, 马宗伟. 长江中下游洪水灾害成因及洪水特征模拟分析. 长江流域资源与环境, 2005, 14 (5): 638- 643. doi: 10.3969/j.issn.1004-8227.2005.05.020
	Xu Y P , Yu R H , Ma Z W . Analysis on the cause of formation of flood disaster and flood characteristics in the middle and lower reaches of the Yangtze River. Resources and Environment in the Yangtze Basin, 2005, 14 (5): 638- 643. doi: 10.3969/j.issn.1004-8227.2005.05.020
	杨顺义, 周兴隆, 陈露露, 等. 甘肃省二斑叶螨地理种群的遗传分析. 植物保护学报, 2018, 45 (6): 1328- 1334.
	Yang S Y , Zhou X L , Chen L L , et al. The genetic analysis of two-spotted spider mite Tetranychus urticae geographical populations in Gansu Province. Journal of Plant Protection, 2018, 45 (6): 1328- 1334.
	姚瑞华, 赵越, 王东, 等. 长江中下游流域水环境现状及污染防治对策. 人民长江, 2014, 45 (S1): 45- 47.
	Yao R H , Zhao Y , Wang D , et al. Status of water environment in middle and lower reaches of the Yangtze River and pollution prevention measures. Yangze River, 2014, 45 (S1): 45- 47.
	虞国跃, 王合, 朱晓清, 等. 北京发现悬铃木方翅网蝽为害. 植物保护, 2014, 40 (5): 200- 202. doi: 10.3969/j.issn.0529-1542.2014.05.039
	Yu G Y , Wang H , Zhu X Q , et al. The first discovery of the sycamore lace bug Corythucha ciliata in Beijing. China Plant Protection, 2014, 40 (5): 200- 202. doi: 10.3969/j.issn.0529-1542.2014.05.039
	Amouroux P , Normand F , Nibouche S , et al. Invasive mango blossom gall midge, Procontarinia mangiferae (Felt)(Diptera: Cecidomyiidae) in Reunion Island: ecological plasticity, permanent and structured populations. Biological Invasions, 2013, 15 (8): 1677- 1693. doi: 10.1007/s10530-012-0400-0
	Arias M B , Elfekih S , Vogler A P . Population genetics and migration pathways of the Mediterranean fruit fly Ceratitis capitata inferred with coalescent methods. PeerJ, 2018, 6, e5340. doi: 10.7717/peerj.5340
	Aslam A F M , Sultana S , Rain F F , et al. Molecular characterization and identification of three stored grain pests based on mitochondrial cytochrome C oxidase subunit I(COI) gene sequences. Bangladesh Journal of Zoology, 2019, 47 (1): 1- 11. doi: 10.3329/bjz.v47i1.42016
	Barthlott W , Hostert A , Kier G , et al. Patterns of vascular plant diversity at continental to global scales. Erdkunde, 2007, 61 (4): 305- 315. doi: 10.3112/erdkunde.2007.04.01
	Beerli P . Comparison of Bayesian and maximum-likelihood inference of population genetic parameters. Bioinformatics, 2006, 22 (3): 341- 345. doi: 10.1093/bioinformatics/bti803
	Bieńkowski A O , Orlova-Bienkowskaja M J . Invasive agricultural pest Drosophila suzukii(Diptera, Drosophilidae) appeared in the Russian Caucasus. Insects, 2020, 11 (11): 826. doi: 10.3390/insects11110826
	Bures G , Czencz K , Geza F , et al. Results of studies on leaf pests of Sycamore maple(Platanus sp.) trees. Novenyvedelem, 1997, 33 (1): 23- 27.
	Cameron S L . Insect mitochondrial genomics: implications for evolution and phylogeny. Annual Review of Entomology, 2014, 59 (1): 95- 117. doi: 10.1146/annurev-ento-011613-162007
	Chin N L , Kek S P . Chemical and genetic markers for identification of honey origin from its bee speciation. International Journal of Food Engineering, 2018, 4 (4): 304- 307.
	Cravedi P . Organisms harmful to plants in the city. Informatore Fitopato-Logico, 2000, 50 (12): 7- 8.
	Halbert S E , Meeker J R . The sycamore lace bug, Corythucha ciliata (Say)(Hemiptera: Tingidae). Annals of the Entomological Society of America, 1983, 76 (2): 262- 265. doi: 10.1093/aesa/76.2.262
	Jermiin L S , Crozier R H . The cytochrome b region in the mitochondrial DNA of the ant Tetraponera rufoniger: Sequence divergence in hymenoptera may be associated with nucleotide content. Journal of Molecular Evolution, 1994, 38 (3): 282- 294.
	Ju R T , Li Y Z , Wang F , et al. Spread of an exotic species, Corythucha ciliata(Say, 1832), in China(Hemiptera: Tingidae). Entomological News, 2009, 120 (4): 409- 414. doi: 10.3157/021.120.0410
	Ju R T , Wang F , Li B . Effects of temperature on the development and population growth of the sycamore lace bug, Corythucha ciliata. Journal of Insect Science, 2011, 11 (16): 16.
	Kimura M , Weiss G . The stepping stone model of population structure and the decease of genetic correlations with distance. Genetics, 1964, 49, 561- 576. doi: 10.1093/genetics/49.4.561
	Lanzavecchia S B , Cladera J L , Faccio P , et al. Origin and distribution of Ceratitis capitata mitochondrial DNA haplotypes in Argentina. Annals of the Entomological Society of America, 2008, 101 (3): 627- 638. doi: 10.1603/0013-8746(2008)101[627:OADOCC]2.0.CO;2
	Lau J A . Genetic variation in invasive species response to direct and indirect species interactions. Biological Invasions, 2015, 17 (2): 651- 659. doi: 10.1007/s10530-014-0756-4
	Maceljski M , Balarin I . Studies of Corythuca ciliata, a new pest in Europe. Anzeiger fur Schadlingskunde Pflanzen-und Umweltschutz, 1974, 47, 165- 170. doi: 10.1007/BF01868932
	Men Q , Xue G , Mu D , et al. Mitochondrial DNA markers reveal high genetic diversity and strong genetic differentiation in populations of Dendrolimus kikuchii Matsumura(Lepidoptera: Lasiocampidae). PLoS ONE, 2017, 12 (6): e0179706. doi: 10.1371/journal.pone.0179706
	Palumbo J C , Perring T M , Millar J G , et al. Biology, ecology, and management of an invasive stink bug, Bagrada hilaris, in North America. Annual Review of Entomology, 2016, 61 (1): 323- 333.
	Petit R J , Aguinagalde I , de Beaulieu J L , et al. Glacial refugia: hotspots but not melting pots of genetic diversity. Science, 2003, 300 (5625): 1563- 1565. doi: 10.1126/science.1083264
	Qin Y J , Buahom N , Krosch M N , et al. Genetic diversity and population structure in Bactrocera correcta(Diptera: Tephritidae) inferred from mtDNA cox1 and microsatellite markers. Scientific Reports, 2016, 6, 38476. doi: 10.1038/srep38476
	Seraphim N , Barreto M A , Almeida G S S , et al. Genetic diversity of Parides ascanius(Lepidoptera: Papilionidae: Troidini): implications for the conservation of Brazil's most iconic endangered invertebrate species. Conservation Genetics, 2016, 17 (3): 533- 546. doi: 10.1007/s10592-015-0802-5
	Siegert N W , McCullough D G , Liebhold A M , et al. Dendrochronological reconstruction of the epicentre and early spread of emerald ash borer in North America. Divers Distrib, 2014, 20, 847- 858. doi: 10.1111/ddi.12212
	Slatkin M . Isolation by distance in equilibrium and non-equilibrium populations. Evolution, 1993, 47 (1): 264- 279. doi: 10.1111/j.1558-5646.1993.tb01215.x
	Streito J C . Note sur quelques especes envahissantes de Tingidae: (Corythucha ciliata Say, 1932), Stephanitis pyrioides(Scott, 1874) et Stephanitis takeeyai Drake & Maa, 1995(Hemiptera: Tingidae). L'Entomologiste, 2006, 62 (1/2): 31- 36.
	Tajima F . DNA polymorphism in a subdivided population: the expected number of segregating sites in the two-subpopulation model. Genetics, 1989, 123 (1): 229- 240. doi: 10.1093/genetics/123.1.229
	Wang X , Jiang L , Qian E , et al. Genetic polymorphisms and haplotypic structure analysis of the Guizhou Gelao ethnic group based on 35 Y-STR loci. Legal Medicine, 2020, 43, 101666. doi: 10.1016/j.legalmed.2019.101666
	Xie D F , Li M J , Tan J B , et al. Phylogeography and genetic effects of habitat fragmentation on Endemic urophysa(ranunculaceae) in Yungui plateau and adjacent regions. PLoS ONE, 2017, 12 (10): e0186378. doi: 10.1371/journal.pone.0186378
	Yang W Y , Tang X T , Cai L , et al. Isolation and characterization of nine microsatelite loci from the sycamore lace bug Corythucha ciliata(Hemiptera: Tingidae). Florida Entomologist, 2014, 97 (3): 1070- 1074. doi: 10.1653/024.097.0310
	Yang W Y , Tang X T , Ju R T , et al. The population genetic structure of Corythucha ciliata(Say)(Hemiptera: Tingidae) provides insights into its distribution and invasiveness. Scientific Reports, 2017, 7 (1): 635. doi: 10.1038/s41598-017-00279-5
	Ye D , Hu Y , Song M , et al. Clonality-climate relationships along latitudinal gradient across China: adaptation of clonality to environments. PLoS ONE, 2014, 9 (4): e94009. doi: 10.1371/journal.pone.0094009

单倍型 Haplotype	共享种群数量 Number of shared populations	共享个体数量 Number of shared individuals	种群中出现的频率 Occurrence frequency in population
Hap1	9	123	CS(20) SH(20) WH(19) BJ(12) CD(4) CQ(11) HZ(15) GY(6) ZZ(17)
Hap2	1	1	WH(1)
Hap3	9	94	BJ(8) SJZ(15) CD(16) CQ(9) HZ(4) JN(14) TY(17) GY(8) ZZ(3)
Hap4	3	5	SJZ(1) JN(1) GY(3)
Hap5	1	8	SJZ(5) JN(2) GY(1)
Hap6	1	1	HZ(1)
Hap7	2	6	JN(3) TY(3)
Hap8	1	1	GY(1)
Hap9	1	1	GY(1)

代码Code	Hap	N	H_d	π	k
CS	Hap1(20)	1	0.000	0.000	0.000
SH	Hap1(20)	1	0.000	0.000	0.000
WH	Hap1(19) Hap2(1)	2	0.100	0.001	0.100
BJ	Hap1(11) Hap3(9)	2	0.521	0.003	2.084
SJZ	Hap3(14) Hap4(1) Hap5(5)	3	0.468	0.002	1.195
CD	Hap1(4) Hap3(16)	2	0.337	0.002	1.347
CQ	Hap1(11) Hap3(9)	2	0.521	0.003	2.084
HZ	Hap1(15) Hap3(4)Hap6(1)	3	0.416	0.002	1.447
JN	Hap3(14) Hap4(1) Hap5(2) Hap7(3)	4	0.500	0.002	1.700
TY	Hap3(17) Hap7(3)	2	0.268	0.001	0.805
GY	Hap1(6) Hap3(8)Hap4(3) Hap5(1) Hap8(1) Hap9(1)	6	0.758	0.005	3.853
ZZ	Hap1(17) Hap7(3)	2	0.268	0.001	1.074
Z		9	0.584	0.003	2.276

变异 Variation	自由度 Degree of freedom	离差平方和 Sum of squares deviations	变异率 Variation rate	固定系数 Fixed parameter
组群之间 Among groups	3	74.34	26.03	F_CT=0.260^△
组群内种群之间 Among populations within groups	11	49.25	21.78	F_SC=0.294^△
种群内个体之间 Among individuals within populations	228	150.20	52.19	F_ST = 0.478^△
总变异Total variation	239	273.79

所属区域 Region	华中地区 Central China	华东地区 East China	西南地区 Southwest China	华北地区 North China
华中地区Central China		10.12	0.49	0.32
华东地区East China	0.05		0.92	0.59
西南地区Southwest China	0.51	0.35		361.82
华北地区North China	0.61	0.46	0.01

	CS	SH	WH	BJ	SJZ	CD	CQ	HZ	JN	TY	GY	ZZ
CS		0.00	0.00	0.69	0.08	0.13	0.69	2.87	0.15	0.06	0.91	4.25
SH	0.00		0.00	0.69	0.08	0.13	0.69	2.87	0.15	0.06	0.91	4.25
WH	0.00	0.00		0.72	0.09	0.14	0.72	3.12	0.16	0.07	0.93	4.64
BJ	0.42	0.42	0.49		0.83	2.12	-10.01	5.38	1.75	1.02	47.58	2.81
SJZ	0.86	0.86	0.85	0.38		6.48	0.83	0.29	42.24	9.23	1.54	0.22
CD	0.79	0.79	0.78	0.19	0.07		2.12	0.51	-22.43	-85.25	3.23	0.37
CQ	0.42	0.42	0.41	-0.05	0.38	0.19		5.38	1.75	1.02	47.58	2.81
HZ	0.150	0.15	0.14	0.09	0.63	0.50	0.09		0.50	0.30	3.35	-12.46
JN	0.77	0.77	0.76	0.22	0.01	-0.02	0.22	0.50		-19.16	3.19	0.37
TY	0.89	0.89	0.88	0.33	0.05	-0.01	0.33	0.62	-0.03		1.73	0.22
GY	0.36	0.36	0.35	0.01	0.25	0.13	0.01	0.13	0.14	0.22		2.32
ZZ	0.11	0.11	0.10	0.15	0.69	0.57	0.15	-0.042	0.57	0.70	0.18