法医学杂志文章

法医学杂志文章

分子层面对生物的研究，在个体水平上主要是看单个基因的变化以及全转录本的变化（RNA-seq）；在对个体的研究的基础上，开始了群体水平的研究。如果说常规的遗传学主要的研究对象是个体或者个体家系的话，那么群体遗传学则是主要研究由不同个体组成的群体的遗传规律。 在测序技术大力发展之前，对群体主要是依靠表型进行研究，如加拉巴哥群岛的13中鸟雀有着不同的喙，达尔文认为这是自然选择造成的后果 。达尔文的进化论对应的观点可以简单概括为“物竞天择，适者生存”，这也是最为大众所接受的一种进化学说。直到1968年，日本遗传学家提出了中性进化理论[2]，也叫中性演化理论。中性理论的提出很大程度上是基于分子生物化学的发展。可以这样理解中性理论：一群人抽奖，在没有内幕的情况下，每个人抽到一等奖的概率是相等的，这个可能性和参与抽奖的人的身高、年龄、爱好等因素都没有关系。中性理论常作为群体遗传研究中的假设理论（CK）来计算其他各种统计指标。 群体遗传学，研究的单位是群体，比如粳稻、籼稻、野生稻，就能够构成不同的群体；我们国内的各省份的水稻也可以作为一个个群体。 群体遗传学大概可以分为群体内的研究和群体间的研究。比如研究云南元阳的水稻的遗传多样性；如果研究是的云南元阳的水稻和东北的水稻，那就可以算成是群体间的研究。群体间和群体内的研究是相互的。 测序价格的急剧下降[3]使得大规模的群体测序得以实现。

常见的变异类型有SNP、IdDel、SV、CNV等。重测序中最关注的是SNP，其次是InDel。其他的几种结构变异的研究不是太多。

有参考基因组的物种的全基因组测序叫做重测序，没有参考基因组的物种的全基因组测序则需要从头组装。随着测序价格的降低，越来越多物种的参考基因组都已经测序组装完成。 plant genomes [4]网站实时显示全基因组测序已经完成的植物，其中2012年以后爆发式增长。在群体遗传学研究中更多的是有参考基因组的物种，尤其是模式物种，植物中常见的是拟南芥、水稻和玉米。

主要的分析流程见下图。现在的测序公司基本上都会帮客户完成整个的分析流程，因为主要耗费的资源是计算资源。我认为在整个分析的流程中最重要的是Linux目录的构建，混乱的目录会导致后续的分析频频出问题，重测序分析会生成很多的中间文件，良好的目录管理会使得项目分析流程井然有序。 该部分涉及到的软件的安装和基础的Linux基础知识就不详细说明了。

正选择似乎可以更好地用自然选择来解释。就是一个基因or位点能够使个体有着更强的生存力或者是育性，这样就会使得这个个体的后代更多，如此一来，这个基因or位点在群体中就越来越多。

正选择能够使有利的突变基因or位点在群体中得到传播，但是与此同时却降低了群体的多态性水平。也就是说原先该位点周围的核苷酸组成是多样性的，在经过正选择之后，这个位点周围核苷酸的多样性就渐渐的趋于同质化了。这就好比一块田，里面本来有水稻和稗草及其他杂草，由于稗草的适应性增强，稗草在逐渐增多，水稻慢慢变少，最后甚至是只剩下了稗草。 我们将这种选择之后多态性降低的情况叫做选择扫荡（Selective Sweep)。检测选择扫荡的软件有SweeD[7]。选择扫荡有可能是人工选择的结果，如2014年 Nature Genetics关于非洲栽培稻的文章就使用了SweeD来检测非洲栽培稻基因组上受人工选择的区域[8]。

负选择和正选择刚好是相反的。简单理解成群体中的某个个体出现了一个致命的突变，从而自己或者是后代从群体中被淘汰。这也导致群体中该位点的多态性的降低。就好比我有10株水稻，其中一株在成长过程中突然不见了，那么对我的这个小的水稻群体来说，这个消失的水稻的独有的位点在群体中就不见了，整体的多态性就降低了。

平衡选择指多个等位基因在一个群体的基因库中以高于遗传漂变预期的频率被保留，如杂合子优势。

平衡选择检测的算法有BetaScan2[10]，这是个Python脚本，输入文件只需要过滤好的SNP数据即可。

计算公式为： 其中 是有效群体大小， 是每个位点的突变速率。 但是群体大小往往是无法精确知道的，需要对其进行估计。

分离位点数 是 的估计值，表示相关基因在多序列比对中表现出多态性的位置。计算公式为： 其中 为分离位点数量，比如SNP数量。 为个体数量的倒数和：

指的是核苷酸多样性，值越大说明核苷酸多样性越高。通常用于衡量群体内的核苷酸多样性，也可以用来推演进化关系[11]。计算公式为： 可以理解成现在群体内两两求 ，再计算群体的均值。计算的软件最常见的是 vcftools ，也有对应的R包 PopGenome 。通常是选定有一定的基因组区域，设定好窗口大小，然后滑动窗口进行计算。 3KRGP文章就计算了水稻不同亚群间4号染色体部分区域上的 值[12]，能够看出控制水稻籽粒落粒性的基因 Sh4 位置多态性在所有的亚群中都降低了。说明这个基因在所有的亚群中都是受到选择的，这可能是人工选择的结果。

Tajima's D是日本学者Tajima Fumio 1989年提出的一种统计检验方法，用于检验DNA序列在演化过程中是否遵循中性演化模型[14]。计算公式为： D值大小有如下三种生物学意义：

叫固定分化指数，用于估计亚群间平均多态性大小与整个种群平均多态性大小的差异，反映的是群体结构的变化。其简单估计的计算公式为： 的取值范围是[0,1]。当 时，表明亚群间有着明显的种群分化。 在中性进化条件下， 的大小主要取决于遗传漂变和迁移等因素的影响。假设种群中的某个等位基因因为对特定的生境的适应度较高而经历适应性选择，那该基因的频率在种群中会升高，种群的分化水平增大，使得种群有着较高的 值。 值可以和GWAS的结果一起进行分析， 超过一定阈值的区域往往和GWAS筛选到的位点是一致的，如2018年棉花重测序的文章[15]：

ROD可以基于野生群体和驯化群体间核苷酸多态性参数 的差异识别选择型号，也可以测量驯化群体和野生型群体相比损失的多态性。计算公式为： 和 一样，ROD也可以和GWAS结合起来：

群体结构分析可以简单理解成采样测序的这些个体可以分成几个小组，以及给每个个体之间的远近关系是怎么样的。群体结构分析三剑客， 分别是 进化树 、 PCA 和 群体结构图 。

进化树就是将个体按照远近关系分别连接起来的图。

常用的绘图软件是 Phylip 和 Snpphylo 。进化树修饰的软件有 MEGA ， ggtree 等，推荐网页版工具 iTOL ，无比强大。 外群定根法：当群体的个体的差异很小时，可以引入其他物种作为根。如在对三叶草建树时可以引入水稻的序列作为根进行建树。

PCA是很常见的降维方法，如微生物研究中常用来检验样品分群情况。PCA计算的软件很多，plink可以直接用vcf文件计算PCA，R语言也可以进行PCA计算。

PCA图在群体重测序中有如下几种作用：

进化树和PCA能够看出来群体是不是分层的，但是无法知道群体分成几个群合适，也无法看出群体间的基因交流，更无法看出个体的混血程度。这时候就需要群体分层图了。

可以将进化树和群体分层图结合进行展示，如下图：

先了解下概念，此处借鉴基迪奥生物网站的解释[22]。 要理解 LD 衰减图，我们就必须先理解连锁不平衡（Linkage disequilibrium，LD）的概念。连锁不平衡是由两个名词构成，连锁 + 不平衡。前者，很容易让我们产生概念混淆；后者，让这个概念变得愈加晦涩。因此从一个类似的概念入手，大家可能更容易理解 LD 的概念，那就是基因的共表达。 基因的共表达，通常指的是两个基因的表达量呈现相关性。比较常见的例子就是：转录组因子和靶基因间的关系。因为转录因子对它的靶基因有正调控作用，所以转录因子的表达量提高会导致靶基因的表达量也上调，两者往往存在正相关关系。这个正相关关系，可以使用相关系数 来度量，这个数值在 - 1~1 之间。总而言之，相关性可以理解为两个元素共同变化，步调一致。 类似的，连锁不平衡（LD）就是度量两个分子标记的基因型变化是否步调一致，存在相关性的指标。如果两个 SNP 标记位置相邻，那么在群体中也会呈现基因型步调一致的情况。比如有两个基因座，分别对应 A/a 和 B/b 两种等位基因。如果两个基因座是相关的，我们将会看到某些基因型往往共同遗传，即某些单倍型的频率会高于期望值。 参照王荣焕等[23]的方法进行LD参数计算：

随着标记间的距离增加，平均的LD程度将降低，呈现出衰减状态，这种情况叫LD衰减。LD衰减分析的作用：

GWAS(genome-wide association study)，全基因组关联分析，常用在医学和农学领域。简单理解成将SNP等遗传标记和表型数据进行关联分析，检测和表型相关的位点，然后再倒回去找到对应的基因，研究其对表型的影响。这些被研究的表型在医学上常常是疾病的表型；在农学上常常是受关注的农艺性状，比如水稻的株高、产量、穗粒数等。GWAS思想首次提出是在心肌梗塞的治疗上[24]，首次应用是在2005年的文章上[25]。

目前使用最广泛的模型是混合线性模型[26]：

所有的参数软件（如Emmax）会自动完成计算。

GWAS结果文件通常只有两个图，一个是曼哈顿图，另外一个是Q-Q图。一般是先看Q-Q图，如果Q-Q正常，曼哈顿图的结果才有意义。

MSMC（multiple sequentially Markovian coalescent）[27]，底层算法很复杂，类似于PSMC。MSMC的主要功能是推断有效群体大小和群体分离历史。

这样看起来更直观：

LAMP(Local Ancestry in Admixed Populations，混杂群体的局部族源推断)，用于推断采用聚类的方法假设同时检测的位点间不存在重组情况，对每组相邻的 SNP 进行检测分析[28]，在运算速度和推断准确度上都有了质的飞跃。

用于推断群体分离和混合[29]。图是这样的：

测序方案关系到后续的分析，不同的样本量对应不同的测序方法和分析方法。

[1]. 自然选择(维基百科) [2]. Kimura, Motoo. "Evolutionary rate at the molecular level." Nature . 217.5129 (1968): 624-626 . [3]. 测序价格变化趋势 [4]. plant genomes [5]. DePristo, Mark A., et al. "A framework for variation discovery and genotyping using next-generation DNA sequencing data." Nature Genetics . 43.5 (2011): 491. [6]. Biswas, Shameek, and Joshua M. Akey. "Genomic insights into positive selection." ** TRENDS in Genetics . 22.8 (2006): 437-446. [7]. Pavlidis, Pavlos, et al. "Sweed: likelihood-based detection of selective sweeps in thousands of genomes." Molecular biology and evolution 30.9 (2013): 2224-2234. [8]. Wang, Muhua, et al. "The genome sequence of African rice (Oryza glaberrima) and evidence for independent domestication." Nature Genetics 46.9 (2014): 982. [9]. Bamshad, Michael, and Stephen P. Wooding. "Signatures of natural selection in the human genome." Nature Reviews Genetics 4.2 (2003): 99. [10]. Siewert, Katherine M., and Benjamin F. Voight. "BetaScan2: Standardized statistics to detect balancing selection utilizing substitution data." BioRxiv (2018): 497255. [11]. Yu, N.; Jensen-Seaman MI; Chemnick L; Ryder O; Li WH (March 2004). Genetics . 166 (3): 1375–83. [12]. Wang, Wensheng, et al. "Genomic variation in 3,010 diverse accessions of Asian cultivated rice." Nature 557.7703 (2018): 43. [13]. Li, C., Zhou, A. & Sang, T. Rice domestication by reducing shattering. Science 311, 1936–1939 (2006). [14]. Tajima, Fumio. "Statistical method for testing the neutral mutation hypothesis by DNA polymorphism." Genetics 123.3 (1989): 585-595. [15]. Du, Xiongming, et al. "Resequencing of 243 diploid cotton accessions based on an updated A genome identifies the genetic basis of key agronomic traits." Nature Genetics 50.6 (2018): 796. [16]. Lu, Kun, et al. "Whole-genome resequencing reveals Brassica napus origin and genetic loci involved in its improvement." Nature communications . 10.1 (2019): 1154. [17]. Zhou, Z., Jiang, Y., Wang, Z. et al. Resequencing 302 wild and cultivated accessions identifies genes related to domestication and improvement in soybean. Nat Biotechnol 33, 408–414 (2015). [18]. Liang, Z., Duan, S., Sheng, J. et al. Whole-genome resequencing of 472 Vitis accessions for grapevine diversity and demographic history analyses. Nat Commun 10, 1190 (2019). [19]. Alexander, D.H., Lange, K. Enhancements to the ADMIXTURE algorithm for individual ancestry estimation. BMC Bioinformatics 12, 246 (2011). [20]. Francis, Roy M. "pophelper: an R package and web app to analyse and visualize population structure." Molecular ecology resources 17.1 (2017): 27-32. [21]. . [22]. . [23]. WANG Rong-Huan, WANG Tian-Yu, LI Yu. Linkage disequilibrium in plant genomes[J]. HEREDITAS , 2007, 29(11): 1317-1323. [24]. Ozaki, K., Ohnishi, Y., Iida, A. et al. Functional SNPs in the lymphotoxin-α gene that are associated with susceptibility to myocardial infarction. Nat Genet 32, 650–654 (2002). [25]. Klein, Robert J., et al. "Complement factor H polymorphism in age-related macular degeneration." Science 308.5720 (2005): 385-389. [26]. Yu, Jianming, et al. "A unified mixed-model method for association mapping that accounts for multiple levels of relatedness." Nature genetics 38.2 (2006): 203. [27]. Schiffels, Stephan, and Richard Durbin. "Inferring human population size and separation history from multiple genome sequences." Nature genetics 46.8 (2014): 919. [28]. Sankararaman, Sriram, et al. "Estimating local ancestry in admixed populations." The American Journal of Human Genetics 82.2 (2008): 290-303. [29]. Pickrell, Joseph K., and Jonathan K. Pritchard. "Inference of population splits and mixtures from genome-wide allele frequency data." PLoS genetics 8.11 (2012): e1002967. [30]. Chen, Jia-hui, et al. "Genome-wide analysis of Cushion willow provides insights into alpine plant divergence in a biodiversity hotspot." Nature communications 10.1 (2019): 1-12. [31]. 孙宽，侯一平。法医族源推断的分子生物学进展 [J]. 法医学杂志 ，2018,34 (03):286-293. [32].

中国农业大学的于静娟有哪些论文？

[1]智联腾,赵倩,敖光明,于静娟. 天然彩色棉纤维特异表达启动子LTP3的克隆及其在烟草中的表达特异性[J]. 热带生物学报,2011,(2).
[2]裴黎,牛慧媛,涂政,彭建雄,马原,从斌,毛泽善,于静娟. 改良高盐低pH方法提取陈旧大麻DNA初探[J]. 中国法医学杂志,2009,(3).
[3]李博,于静娟,赵倩,朱登云,敖光明. 基因pf40在玉米中的遗传转化[J]. 中国农业科学,2009,(9).
[4]PROMMEE Wittaya,王梅珍,朱登云,赵倩,敖光明,于静娟. 水稻悬浮细胞再生及根癌农杆菌介导转化的影响因素(英文)[J]. 中国生物化学与分子生物学报,2009,(11).
[5]梁玉玲,于静娟. 新型白化型除草剂靶标酶对羟苯基丙酮酸双加氧酶及其耐性转基因植物研究进展[J]. 中国生物工程杂志,2009,(12).
[6]夏玉凤,赵倩,于静娟,敖光明. 融合基因gfp-lacZ在克隆载体pGEM3Z-f(+)中高效表达[J]. 生物技术,2006,(4).
[7]王小娟,王海彬,于静娟,杨俊兴,黄志伟. 利用哺乳动物细胞双杂交体系研究仙灵骨葆对雌激素受体的作用[J]. 中国新药杂志,2006,(23).
[8]杨龙峰,于静娟,陈明勇,敖光明. 传染性法氏囊病病毒VP2基因在转基因烟草中的初步表达[J]. 动物医学进展,2007,(2).
[9]杨翅春,于静娟,赵倩,朱登云,敖光明. 玉米花青素调节基因Lc对转基因烟草和矮牵牛花色的影响[J]. 农业生物技术学报,2007,(1).
[10]刘颖慧,于静娟,敖光明,赵倩. 根癌农杆菌介导转化谷子的影响因素(英文)[J]. 中国生物化学与分子生物学报,2007,(7).
[11]李邱华,洪波,仝征,杨英杰,马超,于静娟,高俊平. 新铁炮百合遗传转化体系的建立及Zm401基因的导入[J]. 农业生物技术学报,2008,(1).
[12]黄大军,于静娟,文建成,顾晓明,房毅,李伟华. 滇Ⅰ型不育系繁殖父本花粉密度与母本异交结实的关系[J]. 云南农业大学学报,2008,(3).
[13]裴黎,牛慧媛,毛泽善,于静娟,凃政,彭建雄,张颖,张贵芹,郭佳. 云南大麻DNA的提取及检测初步研究[J]. 刑事技术,2008,(5).
[14]夏玉凤,赵倩,于静娟,敖光明. PF40的亚细胞定位研究[J]. 生物化学与生物物理进展,2005,(11).
[15]冯晓燕,于静娟,赵倩,敖光明. 谷子脂转移蛋白cDNA的克隆及特性研究[J]. 农业生物技术学报,2003,(1).
[16]张文河,赵倩,于静娟,朱登云,敖光明. 转兔防御素基因(NP-1)玉米植株的获得及其抗病性分析[J]. 农业生物技术学报,2003,(4).
[17]赵倩,于静娟,朱登云,敖光明. 番茄花粉特异表达的高赖氨酸蛋白基因(tsb)的克隆与特性分析[J]. 中国生物化学与分子生物学报,2004,(2).
[18]赵倩,赖凡,于静娟,朱登云,敖光明. 谷子类胡萝卜素生物合成途径中番茄红素β-环化酶基因的克隆[J]. 中国农业大学学报,2004,(1).
[19]薛静,于静娟,赵倩,朱登云,敖光明. 谷子f103基因的克隆及其特性分析[J]. 中国农业大学学报,2004,(1).
[20]王为民,赵倩,于静娟,朱登云. 水稻营养品质的改良[J]. 中国生物工程杂志,2004,(5).
[21]王为民,赵倩,朱登云,于静娟,敖光明. 高赖氨酸蛋白质基因和bar基因导入水稻及转基因植株的检测[J]. 农业生物技术学报,2004,(4).
[22]赵倩,于静娟,朱登云,敖光明. 改变天冬氨酸激酶和二氢吡啶二羧酸合成酶调控活性提高玉米中游离赖氨酸含量[J]. 农业生物技术学报,2004,(3).
[23]张革平,于静娟,朱登云,敖光明,赵倩. 外源高赖氨酸蛋白质基因对转基因玉米种子中蛋白质组分的影响[J]. 中国农业科学,2004,(12).
[24]郎志宏,于静娟,朱登云,赵倩,敖光明. 高赖氨酸蛋白基因SBgLR的克隆及其对提高玉米种子中蛋白质和赖氨酸含量的作用[J]. 农业生物技术学报,2004,(5).
[25]薛静,于静娟,赵倩,朱登云,敖光明. 谷子种子特异表达f128基因的克隆及其特性[J]. 农业生物技术学报,2004,(5).
[26]刘颖慧,于静娟,赵倩,朱登云,敖光明. 根癌农杆菌介导谷子的遗传转化[J]. 农业生物技术学报,2005,(1).
[27]王为民,赵倩,于静娟,朱登云,敖光明. 水稻转高赖氨酸蛋白质基因(sb401)植株的获得及种子中蛋白质和氨基酸的含量分析[J]. 作物学报,2005,(5).
[28]于静娟,国凤利,赵德刚,傅永福,韩玉珍,敖光明,孟繁静. 矮牵牛花同源异型基因fbp2的克隆及其对烟草花形态的影响[J]. 植物学报,1999,(1).
[29]张秀君,刘俊起,赵倩,于静娟,敖光明. 用基因枪将高赖氨酸基因导入玉米及转基因植株的检测[J]. 农业生物技术学报,1999,(4).
[30]傅永福,赵德刚,韩玉珍,于静娟,孟繁静. 花同源异型基因FBP2调节叶片中过氧化物酶的表达[J]. 植物生理学报,2000,(4).
[31]张七仙,刘俊起,李成霞,张秀君,赵倩,于静娟,敖光明. 一种简捷的同时提取植物总DNA和总RNA的方法[J]. 中国农业科学,2001,(4).
[32]敖光明,于静娟,赵倩. 天然彩色棉的研究进展及其有待解决的问题[J]. 中国农业科技导报,2001,(1).
[33]孙学辉,敖光明,于静娟,赵倩. 高赖氨酸基因导入玉米自交系的研究[J]. 农业生物技术学报,2001,(2).
[34]李成霞,刘俊起,于静娟,赵倩,敖光明. 玉米花粉特异性基因ZM401 cDNA片段的克隆及其表达研究[J]. 农业生物技术学报,2001,(4).
[35]于静娟,徐南方,孟繁静. 大麻雄蕊特异蛋白质的初步研究[J]. 中国农业大学学报,1997,(3).
[36]于静娟,敖光明. 水稻10kD富硫醇溶蛋白基因在马铃薯中的表达[J]. Acta Botanica Sinica,1997,(4).
[37]于静娟,韩玉珍,赵德刚,傅永福,孟繁静. 玉米赤霉烯酮与大麻的性别表达[J]. 中国农业大学学报,1998,(5).
[38]李秀菊,于静娟,姚槐应,孟繁静. 玉米赤霉烯酮浸种对玉米幼苗抗冷性的影响[J]. 北京农业大学学报,1995,(3).
[39]于静娟,敖光明. 水稻10kd富硫醇溶蛋白基因的分离及序列分析[J]. 农业生物技术学报,1993,(1).
[40]王珲,于静娟,孟繁静. 冬小麦不同发育时期内源玉米赤霉烯酮含量变化对抽穗的影响(简报)[J]. 北京农业大学学报,1992,(1).
[41]郎志宏,于静娟,朱登云,赵倩,敖光明. 高赖氨酸蛋白基因SBgLR的克隆及其对提高玉米种子中蛋白质和赖氨酸的作用[A]. 林栖凤.科学出版社[C].: 科学出版社,2004:.
[42]于静娟,梁翰文,朱登云,薛静,赵倩,敖光明. 一种种子特异性高效启动子及其应用[Z]. CN101063139: ,2007.
[43]敖光明,于静娟,赵倩,郎志宏. 高赖氨酸蛋白基因以及提高禾本科作物种子中赖氨酸和蛋白质含量的方法[Z]. CN1317570: ,2001.
[44]敖光明,赵倩,于静娟,朱登云,彭鹏,张秀君,孙学辉. 改良禾本科作物品质的方法[Z]. CN1317572: ,2001.
[45]敖光明,于静娟,赵倩,戴晓燕,,李成霞. 与花粉发育相关的非编码RNA的cDNA以及制备雄性不育植物的方法[Z]. CN1523108: ,2004.
[46]敖光明,赵倩,于静娟,冯晓燕,,刘颖慧. 植物顶端发育相关的cDNA以及减弱植物顶端优势增加分枝数目的新方法[Z]. CN1523109: ,2004.
[47]苏震,徐文英,于静娟,赵琳娜,刘凤霞,周少霞. 一种植物耐低温蛋白及其编码基因与应用[Z]. CN101289502: ,2008.
[48]苏震,徐文英,于静娟,赵琳娜,刘凤霞,周少霞. 植物耐低温蛋白及其编码基因与应用[Z]. CN101289503: ,2008.
[49]于静娟,敖光明,赵倩,刘俊起,李成霞,王冬雪. 一种花药绒毡层和花粉特异性高效启动子及其应用[Z]. CN101532015: ,2009.
[50]于静娟,赵琳娜,赵倩,敖光明. 利用植物铝诱导表达基因增加植物抗铝性的方法[Z]. CN101532029: ,2009.

高鹏的发表的文章：

(*为通讯作者)1. Gao, P., Chak-Lui Wong, C., Kwok-Kwan Tung, E., Man-Fong Lee, J., Wong, C-M., Oi-Lin Ng, I*. Deregulation of microRNA expression occurs early and accumulates in early stages of HBV-associated multistep hepatocarcinogenesis, Journal of Hepatology 2011;54(6):1177-84.（SCI收录，IF=9.35）2. Gao, P*，Wei JM, Li PY, Zhang CJ, Jian WC, Zhang YH, Xing AY, Zhou GY*.Screening of deoxyribozyme with high reversal efficiency against multidrug resistance in breast carcinoma cells. J Cell Mol Med. 2011;15(10):2130-8.（SCI收录，IF=5.22）3. Gao P*, Zhou GY, Zhang QH, Su ZX, Zhang TG, Xiang L, Wang Y, Zhang SL, Mu K. Lymphangiogenesis in gastric carcinoma correlates with prognosis. Journal of pathology, 2009;218(2):192-200. （SCI收录，IF=7.42）4. Gao P*, Yang X, Xue YW, Zhang XF, Wang Y, Liu WJ, Wu XJ. Promoter methylation of glutathione S-transferase pi1 and multidrug resistance gene 1 in bronchioloalveolar carcinoma and its correlation with DNA methyltransferase 1 expression. Cancer, 2009 115(14):3222-32. （SCI收录，IF=5.43）5. Gao P*, Zhou GY, Yin G, Liu Y, Liu ZY, Zhang J, Hao CY. Lymphatic Vessel density as a prognostic indicator for patients with stage I cervical carcinoma. Human Pathology, 2006;37(6)：719-725. （SCI收录，IF=3.03）6. Gao P*, Zhou GY, Guo LL, Zhang QH, Zhen JH, Fang AJ, Lin XY. Reversal of drug resistance in breast carcinoma cells by anti-mdr1 ribozyme regulated by the tumor-specific MUC-1 promoter. Cancer Letters. 2007;256(1):81-9. （SCI收录，IF=4.86）7. Gao P*, Zhou GY*, Zhang QH, Xiang L, Zhang SL, Li C, Sun YL. Clinicopathological significance of peritumoral lymphatic vessel density in gastric carcinoma. Cancer Letters. 2008 18;263(2):223-30. （SCI收录，IF=4.86）8. Gao P*, Zhou GY*, Lei DP, Zhang XF, Li L, Xu JW, Lin XY. Selection of antisense oligonucleotides for reversal of multidrug resistance in breast carcinoma cells. Cytotherapy. 2007;9(8) 795-801. （SCI收录，IF=3.47）9. Gao P, Zhou GY*, Zhang QH, Li H, Mu K, Yuan YP, Zhang J, Wang BH. Reversal MDR in breast carcinoma cells by ribozyme designed according the secondary structure of mdr1 mRNA. Chinese Journal of Physiology, 2006;49(2) :96-103. （SCI收录，IF=0.68）10. Guo L L, Gao P* (通讯作者), Wu YG, Jian WC, Hao CY, Li H, Lin XY. Alteration of Cyclin D1 in Chinese Patients with Breast Carcinoma and its Correlation with Ki-67, pRb, and p53. Archives of Medical Research, 2007;38(8):846-52. （SCI收录，IF=1.98）11. Gao P*, Zhou GY, Liu Y, Li JS, Zhen JH, Yuan YP. Alteration of cyclin D1 in gastric carcinoma and its clinicopathologic significance. World Journal of Gastroenterology,2004，10（20）2936-2939.（ IF=2.02）12. Gao P*, Yang X, Zhao H. Reply to Promoter methylation of glutathione S-transferase pi1 and multidrug resistance gene 1 in bronchioloalveolar carcinoma and its correlation with DNA methyltransferase 1 expression. Cancer. 2010 Apr 1;116(7):1839.（author reply, SCI收录，IF=5.43）13. Gao P*, Ai-Yan Xing, Geng-Yin Zhou, Ting-Guo Zhang, Chao Gao. The molecular mechanism of microRNA-145 to suppress invasion-metastasis cascade in gastric cancer. Oncogene.2011. Accepted. (SCI收录，IF=7.46）14. Ai-Yan Xing, Xiao-Fang Zhang, Xiu-Quan He, Bin Wang , Wen-Jun Liu, Chao Gao, Duan-Bo Shi, Gao P* (通讯作者). miR-145 is underexpressed in manifold human cancers and might serve as a tumor suppressor. Ann Surg Oncol. ed. (SCI收录，IF=4.19）15. 高鹏， 周庚寅，张庆慧等。乳腺癌组织中多形上皮粘蛋白1表达与肿瘤侵袭力的研究。中华医学杂志，2005；85（6）：381-384。16. 高鹏，周庚寅等。转染抗mdr-1核酶基因逆转肿瘤细胞耐药性的研究，中华病理学杂志,2004，33(3):251-254.17. 高鹏，周庚寅等。负向调节多药耐药基因1转录逆转肿瘤细胞多药耐药。中华病理学杂志，2003；32（6）：563-566。18. 高鹏，周庚寅等.硫化型反义寡核苷酸逆转乳腺癌多药耐药的研究。中华实验外科杂志，2004,21(2),79-80.19. 高鹏，周庚寅等。乳腺癌C-erbB2改变与内分泌治疗敏感性及预后的关系。中国普通外科杂志，2003，10：735-738。20. Gao P, Zhou Geng Yin, Wu Ya Guang, et al. Expression of cyclin D1 protein in breast tumors and its relationship with genetic alteration. US-chinese journal of lymphology and oncology,2003，2：7-10.21. 高鹏，周庚寅等.细胞表面糖蛋白表达异常与涎腺多形性腺瘤复发的关系.华西口腔医学杂志，2005，23（2）：164-166。22. 高鹏，周庚寅等。乳腺癌细胞周期素D1基因改变和表转录异常对其蛋白表达的影响。中国现代普通外科进展，2004，7：31-34。23. 高鹏，王美清等.卵巢上皮性肿瘤多种癌基因改变及蛋白表达的研究。现代妇产科进展，2001，10：346-348。24. 高鹏，王美清等.乳腺良恶性增生中Cyclin D1表达的研究。山东医科大学学报，2001，39：423-425.25. 高鹏，周庚寅. 肿瘤多药耐药及其基因治疗。中国现代普通外科进展，2005，8（2）：72-75。26. 高鹏，周庚寅.免疫组织化学的原理及应用. 山东卫生，2003，27（3）：52-53. 杨希,高鹏*,薛玉文,王妍,刘文君,张晓芳.肺腺癌中DNA 甲基转移酶1 的表达.山东大学学报(医学版),2009:47(1)30-33.27. Sun Y, Zhang T, Gao P, Meng B, Gao Y, Wang X, Zhang J, Wang H, Wu X, Zheng W, Zhou G. Targeting glucosylceramide synthase downregulates expression of the multidrug resistance gene MDR1 and sensitizes breast carcinoma cells to anticancer drugs. Breast Cancer Res Treat. 2010;121(3):591-9.（SCI收录）28. Yin G, Liu YQ, Gao P, Wang XH. Male urethritis glandularis: case report. Chin Med J (Engl). 2007;120(16):1460-1. （SCI收录）29. Zhang X, Li J, Qiu Z, Gao P, Wu X, Zhou G.. Co-suppression of MDR1 (multidrug resistence 1) and GCS (glucosylceramide synthase) restores sensitivity to mutidrug resistance breast cancer cells by RNA interference (RNAi). Cancer Biol Ther. 2009;8(12):1117-21. （SCI收录）30. Sun YL, Zhou GY, Li KN, Gao P, Zhang QH, Zhen JH, Bai YH, Zhang XF. Suppression of glucosylceramide synthase by RNA interference reverses multidrug resistance in human breast cancer cells. Neoplasma. 2006;53(1):1-8. （SCI收录）31. Zhang Y, Wang H, Wei L, Li G, Yu J, Gao Y, Gao P,et al. Transcriptional modulation of BCRP gene to reverse multidrug resistance by toremifene in breast adenocarcinoma cells. Breast Cancer Res Treat. 2010;123(3):679-89. （SCI收录）32. Zhang X, Wu X, Li J, Sun Y, Gao P, Zhang C, Zhang H, Zhou G. MDR1 (multidrug resistence 1) can regulate GCS (glucosylceramide synthase) in breast cancer cells. J Surg Oncol. 2011 Oct;104(5):466-71.（SCI收录）33. Mu K, Li L, Yang Q, Zhang T, Gao P, Meng B, Liu Z, Wang Y, Zhou G. Detection of CHK1 and CCND1 gene copy number changes in breast cancer with dual-colour fluorescence in-situ hybridization. Histopathology. 2011;58(4):601-7. （SCI收录）34. ABCG2 is associated with HER-2 Expression, lymph node metastasis and clinical stage in breast invasive ductal carcinoma. Xiang L, Su P, Xia S, Liu Z, Wang Y, Gao P, Zhou G. Diagn Pathol. 2011 Sep 27;6:90. （SCI收录）35. 刘媛，高鹏等.胎盘细胞凋亡与胎儿生长受限的研究.中华妇产科杂志，2002，12：721-722.36. 周庚寅，高鹏等.疑似星形细胞瘤的脱髓鞘假瘤2例报道及文献复习. 临床与实验病理学杂志.2004.1237. 孙妍琳，高鹏等.肿瘤细胞多药耐药与神经酰胺的糖基化.中国现代普通外科进展.2004，7（3）：139-141.38. 第4作者.PTEN过表达增加人乳腺癌MCF_7细胞对阿霉素的药敏性。中华普通外科杂志，2005，20（4）：240-242.39. 第4作者.葡萄籽多酚逆转人乳腺癌多药耐药性及其机制的研究. 中华普通外科杂志，2004，19（8）：488-490.40. 第4作者.RNA干扰沉默缺氧诱导因子-1α逆转乳腺癌的耐药性. 中华病理学杂志，2006，35（6）：357-360 。41. 第4作者. 吸烟对呼吸道上皮损伤机制的研究.中国公共卫生，2005.21（8）：897-899。42. 第4作者. HIF-1α P-gp在腋淋巴结阴性乳腺癌组织中的表达及意义.中国肿瘤临床.2007,34（2）10-13.43. 第4作者. 葡萄籽多酚对多药耐药的逆转作用及其机制. 山东大学学报（医学版）.2004，42（4）387-389.44. 第4作者.大鼠急性心肌梗死心肌ER表达的变化规律.中国法医学杂志.2005,20(2):84-85.45. 第4作者. 体外瞬时转染野生型PTEN过表达对乳腺癌MCF_7细胞的影响. 山东大学学报（医学版）.2005,43(3):203-207.46. 第4作者. 原发性肝癌中HIF-1α、P-gp的表达及相关性的研究山东大学学报（医学版）.2007,45(3):246-249.47. 第4作者. 脑皮质星形胶质细胞的纯化培养. 山东大学学报（医学版）.2006,44(8):857-859.48. 第5作者. 葡萄籽多酚对逆转人乳腺癌细胞MCF-7/ZAIZ在裸鼠的多药耐药逆转作用. 中华外科杂志，2004，42（13）795-798.49. 第5作者.葡萄糖神经酰胺合成酶在人乳腺癌细胞的表达及其与多药耐药的关系。中华病理学杂志，2005，34（2）：109-110 。50. 第5作者. 腋淋巴结阴性乳腺癌中基因表达PTEN和微血管密度的关系及意义.中国肿瘤临床，2005，32（5）：248-251。51. 第5作者.人乳腺癌细胞中葡萄糖神经酰胺合成酶基因的表达和意义. 中国现代普通外科进展，2005，8（3）:141-143.52. 第6作者.GCS特异性小干扰RNA表达载体的构建及其逆转乳腺癌细胞耐药的研究。中华医学杂志，2005，85（8）：518-521 。53. 第5作者. 小分子干扰RNA联合反义脱氧寡核苷酸靶向逆转人乳腺癌细胞MCF-7/ADR多药耐药的研究.中华乳腺病杂志(电子版),2008，2（4）：436-442 。