ADMA代谢酶基因多态性与中国汉族人群动脉硬化性心脑血管疾病易感性及机制研究

【摘要】第一章ADMA代谢酶基因多态性与中国汉族人群冠心病易感性的关联研究研究背景冠心病(Coronaryheartdisease,CHID)是最主要的心血管疾病,也是造成死亡的主要疾病之一。冠心病的风险因素包括遗传因素和环境因素两个方面。动脉粥样硬化(Atherosclerosis,AS)是冠心病发生的最根本病理改变。非对称二甲基精氨酸(Asymmetricdimethylarginine,ADMA)一种内源性一氧化氮合酶(Nitricoxidesynthase,NOS)抑制剂,通过抑制NO的合成,导致内皮功能紊乱。同时,ADMA通过抑制NO的生物合成和促进炎症过程,在动脉硬化的发生发展中起重要作用。人体内大部分ADMA经ADMA代谢酶代谢灭活,ADMA代谢酶主要包括二甲基精氨酸二甲胺水解酶1(Dimethylargininedimethylaminohydrolase1,DDAH1)和丙氨酸乙醛酸转氨酶2(Alanine-glyoxylateaminotransferase2,AGXT2)代谢灭活。研究表明,ADMA灭活酶基因多态性与心血管疾病(Cardiovasculardisease,CVD)易感性和血浆中ADMA水平均相关,本研究的目的旨在阐明DDAH1和AGXT2基因多态性是否与中国汉族人冠心病易感性相关,同时探究这两个基因多态对血浆中ADMA的影响,为CHD寻找新的易感位点和新的防治策略。方法1.病例对照研究DDAH1和AGXT2基因多态性与CHD易感性的关系收集1103例正常对照和942例CHD患者血液标本,提取基因组DNA。利用聚合酶链式反应-限制性片段长度多态性(PCR-RFLP)的分析方法,对正常对照组和CHD样本进行AGXT2rs37369和1DDAH1rs233113位点基因分型。应用卡方检验和Logistic回归分析基因型与CHD发病风险的关系。2.DDAH1和AGXT2基因多态性对健康志愿者血浆中ADMA水平影响利用酶联免疫吸附测定法(Enzyme-linkedimmunosorbentassay,ELISA)方法测定311例健康志愿者血浆中ADMA浓度。单因素方差分析(One-wayANOVA)和独立样本t检验分析DDAH1和AGXT2基因型与血浆中ADMA水平间关系。结果1.Logistic回归分析对CHD潜在混杂因素进行校正后发现,DDAH1rs233113AT基因型和T等位(AT+TT)可以显著降低总人群中CHD的发病风险(OR=0.79,95%CI:0.65-0.95,P=0.013和OR=0.82,95%CI:0.69-0.98,P=0.028);在非吸烟人群中,rs233113AT基因型(OR=0.73,95%CI-0.58-0.92,P=0.009)和T等位(OR=0.78,95%CI:0.63-0.98,P=0.030)可以显著降低CHD的发病风险；对糖尿病进行分层分析发现,AT基因型(OR=0.77,95%CI:0.63-0.94,P=0.011)和T等位(OR=0.79,95%CI:0.65-0.96,P=0.015)可降低无糖尿病个体CHD的发病风险；同时AT基因型(OR=0.73,95%CI:0.56-0.93,P=0.013)和T等位(OR=0.76,95%CI：0.60-0.96,P=0.023)可以显著降低不吸烟人群中无糖尿病个体CHD发病风险。2.CHD患者中AGXT2rs37369GG基因型的频率显著高于对照组(18.5%vs14.8%,P=0.025);与rs37369A(AA+AG)基因型相比,rs37369GG基因型可显著升高吸烟人群中CHD的发病风险(OR=2.21,95%CI:1.24-3.92,P=0.007),以及显著升高吸烟且同时患有糖尿病个体的CHD发病风险(OR=3.32,95%CI：1.14-9.67,P=0.028)。3.在吸烟个体中,DDAH1rs233113AT基因型个体血浆中ADMA水平显著高于其它两种基因型(P=0.038和P=0.036)；在不吸烟人群中,AGXT2rs37369GG基因型个体血浆中ADMA水平显著低于A等位(AA+AG)携带者(0.45±0.05μM,n=43vs0.65±0.04μM,n=184,P=0.003),但是在吸烟人群中rs37369GG基因型个体血浆中ADMA水平显著高于rs37369A等位基因携带者(0.92±0.16μM,n=14vs0.51±0.05μM,n=70,P=0.004),吸烟人群中AGXT2rs37369GG基因型个体血浆ADMA水平也是显著高于不吸烟人群rs37369GG基因型个体(P=-0.012)。4.携带DDAHlrs233113和AGXT2rs37369两个SNP位点任一风险基因型均可显著升高冠心病的发病风险,且健康对照吸烟个体中同时携带两个SNP位点风险基因型的血浆中ADMA水平显著升高(P=0.003)。结论1.DDAHlrs233113多态可降低中国汉族人群冠心病的发病风险,尤其可以降低不吸烟和不伴糖尿病的个体冠心病的发病风险。2.AGXT2rs37369多态可显著升高中国汉族吸烟人群冠心病的发病风险,尤其可升高吸烟且同时伴有糖尿病的个体冠心病的发病风险,其机制可能与吸烟状态下血浆中ADMA水平升高有关。3.DDAH1rs233113和AGXT2rs37369多态性之间存在联合作用,共同影响中国汉族人群冠心病的发病风险。第二章ADMA代谢酶基因多态性与中国汉族人群缺血性脑卒中易感性和预后的关联研究研究背景脑卒中(Stroke)是全世界范围内最主要的致残原因,在所有死因中位居第二,缺血性脑卒中(Ischemicstroke,IS)是最常见的脑卒中类型,占全部脑卒中的60%-80%。我国人群脑卒中发病率和患病率总体呈上升趋势,并且其致残率和致死率较高,是严重危害我国老年人健康与生命的主要疾病。IS病因也包括环境因素和遗传因素。颈动脉粥样硬化是IS最重要的病因和危险因素,它与梗塞的发生、发展、复发以及梗死的部位密切相关。非对称二甲基精氨酸(Asymmetricdimethylarginine,ADMA)通过抑制一氧化氮合酶(Nitricoxidesynthase,NOS)活性抑制NO的合成,从而导致内皮功能紊乱。ADMA造成的内皮功能紊乱是动脉硬化发生的起始步骤。人体内绝大部分ADMA由二甲基精氨酸二甲胺水解酶1(Dimethylargininedimethylaminohydrolase1,DDAH1)和丙氨酸乙醛酸转氨酶2(Alanine-glyoxylateaminotransferase2,AGXT2)代谢灭活,研究表明ADMA代谢酶基因多态性与循环中ADMA水平相关。本研目的旨在探讨AGXT2和DDAH1基因的遗传多态性是否与IS易感性和预后相关,为IS寻找新的易感位点和新的防治策略。方法1.病例对照研究AGXT2和DDAHl基因多态性与IS易感性的关系收集1103例正常对照和1012例IS患者血液标本,提取基因组DNA。利用聚合酶链式反应-限制性片段长度多态性(PCR-RFLP)方法,对正常对照组和IS样本进行AGXT2rs37369、DDAH1rs233113和rs233112位点的基因分型。应用卡方检验和Logistic回归分析基因型与IS发病风险的关系。2.AGXT2和DDAH1基因多态性与IS预后关系对部分IS患者进行电话随访,详细记录患者是否有事件发生(包括是否死亡、是否复发脑梗塞和心梗梗塞),事件发生时间等信息,应用卡方检验和Kaplan-Meier分析AGXT2和DDAH1基因多态性对IS患者预后的影响。结果1.AGXT2rs37369GG基因型和A等位(AA+AG)分布在病例组和对种组分布无显著差异(P=0.163),同时与IS易感性也不相关,。2.DDAH1rs233113与总人群、吸烟人群和不吸烟人群中IS易感性不相关。以高血压和糖尿病分层分析发现,与rs233113AA基因型相比,AT基因型可以显著降低糖尿病患者IS发病风险(OR=0.73,95%CI=0.53-0.99,P=0.046)。3.对照组中DDAH1rs233112G等位基因(AG+GG)频率显著高于病例组(71.6%vs66.7%,P=0.020)。Logistic回归分析发现,与rs233112AA基因型个体相比,总人群中AG基因型和GG基因型个体以及G等位(AG+GG)携带者IS的发病风险均显著降低(OR=0.78,95%CI=0.63-0.96,P=0.019;OR=0.76,95%CI=0.58-0.99,P=0.045和OR＝0.77,95%CI=0.63-0.94,P=0.011);rs233112AG基因型、GG基因型和G等位携带(AG+GG)还可显著降低高血压个体IS的发病风险(OR=0.73,95%CI=0.57-0.92,P=0.007；OR=0.71,95％CI=O.52-0.96,P=0.026和OR=0.72,95%CI=0.58-0.90,P=0.003)；同时rs233112AG基因型和G等位(AG+GG)可以显著降低糖尿病患者IS的发病风险(OR=0.59,95%CI=0.42-0.84,P=0.003和OR=0.65,95%CI=0.47-0.90.P=0.010)。4.与同时携带DDAH1rs233113和rs233112保护基因型个体比,同时携带两个位点风险基因型(rs233113AA,rs233112AA)个体IS风险显著升高(OR=1.30,95%CI=1.05-1.62,P=0.017),且IS患者风险单倍型(rs233113A/rs233112A)的频率显著高于rs233113A/rs233112G和rs23113T/rs233112A单倍型(P=0.017和P=0.026)。5.AGXT2rs37369、DDAH1rs233113和rs233112与IS预后不相关,但是在发现人群和合并人群中有事件发生的IS患者中rs233113T等位(AT+TT)频率有升高趋势(P=0.087和P=0.097)。结论1.DDAH1rs233113和rs233112可以降低IS的发病风险,尤其可以降低伴有糖尿和高血压的个体IS的发病风险。2.DDAH1rs233113和rs233112可以通过联合作用升高IS的发病风险。3.AGXT2rs37369多态与IS易感性不相关。4.AGXT2rs37369、DDAH1rs233113和rs233112与IS的预后无显著相关。第三章miR-455-5p对DDAHl表达调节作用及DDAHl遗传多态性的影响研究背景通过本论文第一、二章的研究,我们发现DDAH1rs233113多态与冠心病(Coronaryheartdisease,CHD)的易感性相关,同时DDAH1rs233113和rs233112多态还与缺血性脑卒中(Ischemicstroke,IS)易感性相关。本课题组的前期研究发现,rs233113和rs233112多态可显著升高原代培养的脐静脉内皮细胞(HUVEC)内DDAH1的表达,升高细胞裂解液对ADMA的代谢活性,降低细胞内ADMA的浓度,但这两个SNP是否为功能性SNP,以及两个SNP发挥作用的机制仍不清楚。微小RNA(MicroRNA,miRNA)一类内源性的长约22个碱基的非编码RNA,通过与靶基因3’非编码区Jntranslatedregions,UTR)以碱基互补配对方式结合调控靶基因的表达。DDAH1的表达受微小RNA(MicroRNA,miRNA)调控。生物信息学分析发现,miR-455-5p在人DDAH13’UTR区有两个的miRNA识别元件(MicroRNArecognitionelements,MRE),其中MRE1位于rs233113上游100bp左右处,MRE2位于rs233112处。通过查阅文献,我们发现miR-455-5p是一种与炎症反应密切相关的miRNA,炎症状态下其表达显著升高。因此,基于本学位论文前两章的研究基础和以上背景,我们提出DDAH1rs233112和rs233113多态可能通过影响miR-455-5p与DDAHl的结合,从而影响miR-455-5p对DDAH1的表达调控,最终影响动脉硬化性心脑血管疾病的发生发展。本研究旨在查明miR-455-5p是否对DDAH1产生直接抑制作用,查明rs233112和rs233113多态在miR-455-5p介导DDAH1表达调节中的作用,以明确rs233112和rs233113多态性影响心脑血管疾病发病风险和预后的分子机制,为动脉硬化性心脑血管疾病寻找新的遗传标记物和防治策略。方法1.miR-455-5p模拟物(mimic)转染HUVEC-C观察对miR-455-5p表达、DDAH1mRNA和蛋白表达的影响使用不同浓度miR-455-5pmimic(20.40.80nM)转染HUVEC-C24h,real-timePCR检测niR-455-5p和DDAH1mRNA表达,western-blot检测DDAH1蛋白表达。2.健康志愿者外周血单个核细胞(PBMC)中观察miR-455-5p和DDAH1mRNA表达水平间相关性及rs233113和rs233112多态性的影响招募健康志愿者并抽取外周血,分离PBMC,提取PBMC中的DNA和RNA,利用聚合酶链式反应-限制性片段长度多态性(PCR-RFLP)方法对每个个体进行DDAH1rs233113和rs233112基因分型,应用real-timePCR检测miR-455-5p和DDAH1mRNA表达,Pearson相关性分析miR-455-5p与DDAH1mRNA表达水平间的相关性。3.原代HUVEC中观察miR-455-5p与DDAH1mRNA表达、ADMA代谢活性和细胞内外ADMA水平间相关性及rs233113和rs233112多态性的影响分离和培养原代HUVEC,PCR-RFLP对脐静脉组织进行基因分型,第4代内皮细胞于6孔板中培养24小时后,real-timePCR检测DDAH1mRNA表达,western-blot检测DDAH1蛋白表达,ELISA检测细胞内外ADMA浓度以及ADMA代谢活性,Pearson相关性分析进行结果分析。4.肿瘤坏死因子α(TNF-α)处理HUVEC-C观察miR-455-5P、DDAH1mRNA和蛋白表达变化不同浓度TNF-α(25、50、100ng/ml)处理HUVEC-C24h后,real-timePCR检测miR-455-5p和DDAH1mRNA表达,western-blot检测DDAH1蛋白表达。5.HUVEC-C和HEK293细胞中miR-455-5p与不同报告基因质粒共转染对报告基因活性的影响PCR扩增DDAH1基因3’-UTR区含miR-455-5p两个MRE的序列,连接到双报告基因质粒pmirGLO上,定点突变两个位点分别包含rs233113和rs233112多态位点不同单倍型的质粒(rs233113A-rs233112A、rs233113A-rs233112G、rs233113T-rs233112A和rs233113T-rs233112G),以及两个MRE分别为野生型和突变体的质粒[MRE1-A-A(野生型质粒)、MRE1-A-MRE2mut、MRE1-T-MRE2mut、MRE1mut-A-A、MRE1mut-A-G、MRE1mut-A-MRE2mut],将报告基因质粒与不同浓度miR-455-5pmimic(10、20、40nM)或阴性对照(20nM)共同转染HUVEC-C和HEK293细胞24小时后,收集细胞并取裂解液,双荧光素酶报告实验测定荧光素酶活性,查明miR-455-5p模拟物对不同质粒报告基因活性的影响。结果1.HUVEC-C转染miR-455-5pmimic24h呈浓度依赖性地上调细胞内miR-455-5p表达(P<0.05),同时DDAH1三个转录本的mRNA表达及DDAH1蛋白表达显著下降。2.不同浓度TNF-α(25、50、100ng/ml)处理HUVEC24h均可显著上调HUVEC细胞内miR-455-5p表达(P<0.05),100ng/mlTNF-α可显著下调DDAH1三个转录本mRNA表达(P<0.05),50ng/ml和100ng/mlTNF-α处理可显著下调DDAH1蛋白表达,但25ng/mlTNF-α对DDAH1的蛋白表达无显著影响。3.在携带rs233113T等位基因的个体来源的PBMC中,miR-455-5p和DDAH1-V1mRNA表达呈显著负相关(R=-0.636,P=0.026,n=12)；在原代培养的HUVEC中,]miR-455-5p和DDAH1-V1mRNA表达也呈显著的负相关(R=-0.512,P=0.018,n=21),且这种显著的相关性主要出现在携带rs233113T等位和rs233112G等位的个体来源的HUVEC中(rs233113T等位基因携带者：R=-0.764,P=0.017,n=9;rs233112G等位基因携带者：R=-0.660,P=-0.010,n=14)。4.在原代培养的HUVEC中,miR-455-5p表达水平与细胞裂解液ADMA的代谢活性间呈显著的负相关(R=-0.608,P=0.036,n=12),与细胞内ADMA水平呈显著正相关(R=0.628,P=0.029,n=12)。5.在HUVEC-C和HEK293细胞中,10、20、40nM的miR-455-5pmimic均可显著下调rs233113A-rs233112A、rs233113A-rs233112G和rs233113T-rs233112A报告基因质粒活性；但对于rs233113T-rs233112G报告基因质粒,在HUVEC-C中只有20和40nMmiR-455-5pmimic处理可显著下调该质粒的报告基因活性,而在HEK293细胞中,只有40nMmiR-455-5pmimic才可显著下调该报告基因质粒的报告基因活性。6.HUVEC-C中,10、20、40nM的miR-455-5pmimic均可以显著下调MRE1-A-A(野生型质粒)和MRE1-A-MRE2mut(MRE2识别位点突变质粒)质粒报告基因活性,20和40nM的1niR-455-5pmimic可显著下调MRE1mut-A-A(MRE1识别位点突变质粒)质粒报告基因活性,但是各浓度miR-455-5pmimic均不影响MRE1-T-MRE2mut.MRE1mut-A-G和MRE1mut-A-MRE2mut质粒的报告基因活性。结论1.DDAHl是miR-455-5p的靶基因,同时通过DDAH13’-UTR区的两个MRE发挥抑制DDAH1表达的作用；2.TNF-α可通过上调miR-455-5p的表达,从而下调DDAH1的表达。3.DDAH1rs233113和rs233112多态性分别通过MRE1和MRE2影响miR-455-5p下调DDAH1表达的作用。4.DDAH1rs233113和rs233112构成的单倍型对DDAH1报告基因质粒活性的影响高于单个SNP位点的影响。
【作者】周继朋；
【导师】陈小平；
【作者基本信息】中南大学，药学，2014，博士
【关键词】冠心病；二甲基精氨酸二甲胺水解酶1；丙氨酸乙醛酸转氨酶2；非对称二甲基精氨酸；单核苷酸多态性；缺血性脑卒中；微小RNA；微小RNA识别元件；报告基因；3’非编码区；

【参考文献】
[1]王伟.江苏省人口分布合理性评价[D].南京财经大学,统计学,2013,硕士.
[2]周炳强.辽东湾地区吹填土电渗加固实验及模拟研究[D].吉林大学,地质工程,2014,硕士.
[3]吕乐华.高速列车摩擦接地装置材料的制备及其载流摩擦磨损性能的研究[D].河南科技大学,材料加工工程,2013,硕士.
[4]刘艳伟.E1A基因对人肺腺癌细胞增殖和转移的影响[D].山西医科大学,病理学,2002,硕士.
[5]徐国芹.基于AGLET移动代理技术的分布式入侵检测系统的研究[D].辽宁工程技术大学,计算机应用技术,2004,硕士.
[6]滕青.英语专业大学生二外法语学习动机调查研究[D].山东大学,英语语言文学,2012,硕士.
[7]刘秀红.基于模糊PID控制的电阻炉炉温系统的研究[D].燕山大学,测试计量技术及仪器,2004,硕士.
[8]陈芳.桂林民俗文化旅游开发研究[D].中南林学院,2003.
[9]周铁军.具三段常数不连续信号传递函数的二元离散神经网络模型的动力学研究[D].湖南大学,应用数学,2003,硕士.
[10]高兵.XML在关系数据库中存储和查询方法的研究[D].哈尔滨工程大学,2004.
[11]刘倩,王坤,吴向阳,关庆猛,李赫楠,吴君.CdS量子点/石墨烯复合物的制备及电致化学发光性能研究[A].中国化学会.中国化学会第27届学术年会第09分会场摘要集[C].中国化学会:,2010:1.
[12]张桂华.大学生公民意识教育研究[D].山东师范大学,思想政治教育,2013,硕士.
[13]陈英杰.会计信息质量与投资者利益保护问题研究[D].广西师范大学,企业管理,2013,硕士.
[14]李鸣.药品生产企业质量受权人制度构建研究[D].复旦大学,药剂学,2012,硕士.
[15]张晓东.太阳能电动车模型设计[D].沈阳理工大学,车辆工程,2012,硕士.
[16]邓浛笑.上市公司股份代持法律问题研究[D].华中师范大学,民商法学,2014,硕士.
[17]杜美.鱼精蛋白包覆PLGA纳米粒用于乙肝疫苗载体的初步研究[D].华中科技大学,生物化学与分子生物学,2013,硕士.
[18]王广生.内置式永磁同步电机设计及弱磁性能研究[D].湖南大学,电气工程,2012,硕士.
[19]余盛君.从绘画语言的找寻看创作过程[D].新疆师范大学,美术（专业学位）,2013,硕士.
[20]王曼君.质量体系文件的编制[J].建筑技术.1997(08)
[21]周意,李峰峰.基于CRM的银行业务应用模块分析[J].中国金融电脑,2004,02:30-33.
[22]孔娟.姜黄素对心肌梗死后大鼠心肌NF-κB、PPAR-γ、BCL-2表达的影响[D].新乡医学院,内科学,2012,硕士.
[23]马杰.城市轨道交通车站通行设施通过能力研究[D].北京交通大学,2010.
[24]窦婷婷.农村人力资本与农民收入增长的关系研究[D].河北经贸大学,行政管理,2012,硕士.
[25]廖鹏.磷酸酯酶SCP1调控c-Myc蛋白去磷酸化及相关生物学功能的研究和SCP1细胞膜定位的发现及其棕榈酰化修饰的研究[D].华东师范大学,生物医学,2012,博士.
[26]许再香.从《喧哗与骚动》看福克纳的女性观[D].安徽大学,英语语言文学,2014,硕士.
[27]郑文光.对在任者投机行为可能引致经济增长的机制的研究[D].东北财经大学,西方经济学,2012,硕士.
[28]李然.粉末压延技术制备高硅铁硅合金[D].武汉理工大学,材料加工工程,2004,硕士.
[29]杨新国.初中生成就目标、学习策略与学业成绩的研究[D].广西师范大学,课程与教学论,2003,硕士.
[30]王静.全球价值链中的中国汽车产业升级研究[D].吉林大学,世界经济,2013,硕士.
[31]王绪海.基于特定编码方式的RFID防碰撞算法研究[D].南京邮电大学,电子与通信工程（专业学位）,2013,硕士.
[32]王英杰.南瓜籽仁粉冲调饮品的开发及造粒工艺的研究[D].浙江工商大学,食品工程,2014,硕士.
[33]孙嘉轶.基于再制造的闭环供应链结构选择及协调模型研究[D].哈尔滨理工大学,管理科学与工程,2013,博士.
[34]陈伟,于泳,徐殿国,徐壮.基于自适应补偿的异步电机静止参数辨识方法[J].中国电机工程学报,2012,06:156-162.
[35]马丽卿.收入结构对我国上市商业银行绩效影响的实证分析[D].山西财经大学,金融学,2013,硕士.
[36]卢阳阳.生活污水同步硝化反硝化脱氮研究[D].北京交通大学,2014.
[37]周婷.基于改进Hopfield网络的注塑模具制造车间生产调度的研究[D].华南理工大学,机械工程,2012,硕士.
[38]肖留军.基于中医传承辅助系统的治疗泄泻（久泻）组方规律研究[D].广州中医药大学,中医内科学（专业学位）,2013,硕士.
[39]张兵.一种高效的序列模式增量挖掘算法（NPSP）[D].广西师范大学,计算机软件与理论,2003,硕士.
[40]郑涛.基于虚拟仪器的高速水洞数据采集系统开发[D].西北工业大学,机械电子工程,2004,硕士.
[41]郭娇龙.中国内地省域建筑业劳动生产率水平及影响因素研究[D].长安大学,技术经济及管理,2014,硕士.
[42]于化龙.西北地区FDI对国内投资的挤入挤出效应研究[D].西北师范大学,统计学,2013,硕士.
[43]汪琼.子空间分解类算法在参数估计中的应用[D].扬州大学,信号与信息处理,2011,硕士.
[44]王剑.CF公司库存管理优化研究[D].浙江理工大学,工商管理（专业学位）,2014,硕士.
[45]刘海.课堂模糊限制语的语用学研究[D].苏州大学,外国语言学及应用语言学,2013,硕士.
[46]王占中.基于CORBA的入侵检测系统研究[D].广西大学,控制理论与控制工程,2004,硕士.
[47]吕婕.大型球磨机筒体应力有限元分析[D].沈阳工业大学,机械电子工程,2012,硕士.
[48]刘健辰,章兢,张红强,何敏.时变时滞离散广义Markov跳变系统的鲁棒稳定性[J].控制与决策,2011,05:667-672+678.
[49]李丽娜.长江口滨岸潮滩大型底栖动物重金属的分布累积及其生态毒理效应[D].华东师范大学,自然地理学,2004,硕士.
[50]刘敏.基于收益共享的农超对接模式及利益分配问题研究[D].浙江工业大学,2013.