MicroRNA (miRNA)是新近发现并引起广泛关注的一类非编码小RNA分子,它在生物体内发挥重要的基因表达调控作用,是一套全新的表达调控机制。它广泛存在于动物、植物、细菌、病毒等多种生物物种中,在最新的miRBase Virsion16.0中共收录miRNA 15172,成熟1miRNA序列超过17000条,涉及142个物种,其中大多数miRNA在生物进化过程中高度保守。目前发现,1niRNA是生物体内最丰富最重要的一类基因调控因子,在内源基因的调控以及基因组对抗外源核酸侵入中起着重要的作用。在真核生物中,miRNA基因约占整个基因组基因总数的2%,其可能调控的mRNA占全部mRNA的比例高达30%。研究发现miRNAs在细胞凋亡、癌症发生、胚胎发育、干细胞分化、炎症、生长发育、病毒感染等过程中发挥重要作用。研究表明miRNA通过与靶基因mRNA特异性的碱基配对,引起靶mRNA的降解或者抑制其翻译,从而对基因表达进行转录后调控。每个miRNA可以有多个靶基因,而几个miRNAs也可以调节同一个基因。这种复杂的调节网络既可以通过一个miRNA来调控多个基因的表达,也可以通过几个1miRNAs的组合来精细调控某个基因的表达。在线虫、果蝇、小鼠和人等物种已经发现的miRNAs功能研究中,多数具有和其他参与调控基因表达的分子一样的特征,即在不同组织、不同发育阶段中miRNA的水平有显著差异,这种miRNA表达模式具有的位相性和时序性,提示miRNA有可能参与基因表达的调控。正是由于miRNA具有的基因表达调控功能和特点,使其对基因功能研究、生物进化探索和人类疾病防治具有重大意义,因而成为继siRNA (small interfering RNA)之后的新的研究热点。在蚊虫的miRNA研究领域,2005年,Wang等首次运用生物信息学的技术方法,鉴定了38个冈比亚按蚊的miRNAs。2007年,Winter等首次实验验证冈比亚按蚊miRNAs,比较了冈比亚按蚊不同时期的miRNA表达谱,并通过对敲除Dicer1和Agol mRNA(这两种酶与miRNA的生成有关)前后冈比亚按蚊感染疟原虫的生长发育情况表明miRNA作为一类新型的调节因子可能参与按蚊对疟原虫的防御反应。Mead等和Li等人通过实验分别鉴定了斯氏按蚊埃及伊蚊miRNAs,并揭示miRNA可能在生长发育及传播疾病过程中起到重要调节作用。2010年Scalsky等分析致卷库蚊成蚊和白纹伊蚊C7/10细胞感染西尼罗病毒前后miRNA表达差异,提示miRNA可能与病毒感染有关。白纹伊蚊作为登革病毒传播重要的传播媒介之一,主要分布在我国南部和东南亚地区。近年来登革热病例正在迅猛增加,其流行已威胁到全球二分之一人口的健康安全。在我国,登革热的流行以广东省最为严峻。自1978年广东佛山首次暴发登革热以来,广东省已发生过12次大小不等的登革热流行,近年来更常有暴发流行。登革热在广东的流行,不仅给广大人民群众的生命健康带来了严重危害,而且还经常转化为突发的公共卫生事件,给社会带来不安定,并严重影响经济的发展。由于目前还没有有效的疫苗可用于登革热的预防,预防登革热只能从控制传染源和传播媒介方面着手,而控制传染源难度较大,故此,而今预防登革热的唯一对策就是控制其传播媒介----埃及伊蚊和白纹伊蚊,但目前这方面的基础研究,如蚊虫对登革病毒的易感性和传播以及现代生物技术在蚊媒及其传染病防治中的应用,还很欠缺、亟待加强。因此,弄清楚白纹伊蚊对登革病毒的易感性基础对于制定有针对性的防制措施非常重要。白纹伊蚊miRNA是否存在?其序列构成如何?有何表达特征?在登革病毒的感染过程中起什么作用?能否从中找到阻断登革热传播的靶标呢?在本研究开始之前,有关白纹伊蚊miRNAs研究尚未见报道。本研究利用新一代Solexa高通量测序技术,通过对白纹伊蚊miRNAs的组学特征以及miRNAs在白纹伊蚊不同发育阶段的表达特征进行深入分析。Solexa测序是Illumina公司推出的新一代高通量、大规模测序技术,具有高准确性、高通量、高灵敏度、和低运行成本等突出优势,目前已广泛应用于基因组测序、转录组测序以及小RNA测序领域。其中在miRNA研究领域,通过小分子RNA大规模测序分析,可以从中获得物种全基因组水平的miRNA图谱,实现包括新miRNA分子的挖掘,其作用靶基因的预测和鉴定、样品间差异表达分析、miRNAs聚类、进化分析以及表达谱分析等科学应用。因为这一测序技术不需要依赖于基因组序列,因此可以对目前尚无基因组序列的物种进行测序分析,拓宽了miRNA的研究领域,既极大的丰富了miRNA的种类,又为miRNA的结构特点以及功能的研究提供了便利。在本课题开展过程中,人们通过Solexa高通量测序技术对鸡胚、家蚕、蜜蜂、蝗虫、库蚊和人等多个物种的小RNA进行了深度测序,极大的丰富了miRNA的种类,促进了miRNA的深入研究。本课题利用新一代高通量测序和生物信息学的技术方法,首次对白纹伊蚊miRNA的组学进行系统性分析,并将其组学特征与果蝇、按蚊的miRNA进行比较分析；应用Solexa和Northern blot分析miRNA在白纹伊蚊不同生长发育阶段的表达特征,为进一步探索miRNA在蚊虫生长发育过程中的调控机制提供依据；应用Solexa测序技术与生物信息学方法分析经登革病毒感染前后的miRNA表达差异及其靶基因,为登革病毒调控机制的研究提供新的思路和方法。研究目的：1.通过新一代高通量测序技术Solexa测序和生物信息学的技术方法,对白纹伊蚊miRNAs的组学进行系统性分析,丰富蚊虫miRNA并揭示新的蚊虫特异miRNA。2.应用Solexa测序和Northern blot对miRNA在白纹伊蚊不同生长发育阶段的表达特征及种系进化分析,为进一步探索miRNA在蚊虫生长发育过程中的调控机制提供依据。3.应用Solexa测序技术与生物信息学的方法分析经登革病毒感染前后的miRNA表达差异及预测其靶基因,初探调控登革病毒靶标,为研究登革病毒感染调控机制奠定基础研究方法：1高通量测序及其信息学分析白纹伊蚊miRNA方法的建立：1.1样品收集：分别收集12-24h虫卵、早期幼虫、晚期幼虫、虫蛹、羽化3-5d雄蚊、雌蚊等6个不同阶段样本。1.2 Total RNA提取和鉴定：用Trizol提取各个不同阶段样本的Total RNA,核酸测量仪检测浓度及纯度,15%尿素变性聚丙烯酰胺凝胶电泳鉴定所提Total RNA的质量,将6个样品total RNA等量混匀后用Agilent生物分析仪分析RNA完整性。1.3传统方法建库及验证：采用TA克隆和传统Sanger测序对所提Total RNA鉴定是否可以用于Solexa测序。1.4 Solexa高通量测序：首先将18-30nt范围的small RNA从Total RNA中分离出来,两端分别加上特定接头后体外反转录成cDNA,经RT-PCR扩增后用Solexa测序仪对DNA片段进行单向末端直接测序。1.5原始序列的处理与统计：通过去接头,去低质量,去污染、统计序列长度分布等过程,获得去杂质的“干净”序列(Clean序列)。1.6小RNA序列的比对分析：将上述clean序列通过与NCBI Genbank (http:www.ncbi.nlm. nih.gov/)中rRNA、tRNA、snRNA、snoRNA、repeat、exon、intron、miRNA等库进行比对和分类注释。1.7已知1miRNA鉴定和新miRNA预测：用SOAP软件选择没有比对任何注释的小RNA与埃及伊蚊基因组序列比对,Ensembl法(http://www.ensembl.org)获得了与埃及伊蚊吻合序列及其100nt侧翼基因序列,用MIREAP软件鉴定和预测miRNA,最后用RNAfold(http://ma.tbi.univie.ac.at/cgi-bin/ RNAfold.cgi)分析其最低能量构象。1.8用Northern blot方法验证Solexa测序技术分析鉴定miRNAs2白纹伊蚊IniRNAs组学的鉴定：2.1用上述建立的Solexa高通量测序和生物信息学分析虫卵、幼虫、蛹、雌雄成虫及吸血后的6个发育阶段miRNAs。2.2将鉴定miRNAs与miRBase中库蚊、按蚊、果蝇、意蜂、家蚕、原鸡、斑马鱼、小鼠和人等物种miRNA进行种系进化分析,并将获得的miRNAs与三个蚊种基因组序列进行比对。2.3 miRNA基因簇分析：以Altuvia等的miRNA成簇分析方法为基础加以改善进行miRNA成簇分析。以1000nt作为相邻miRNA间距离的最大阈值来寻找染色体上同一方向可成簇miRNA。3白纹伊蚊miRNAs表达谱分析及Northern blot验证：3.1 miRNA差异表达分析：根据Solexa测序结果,统计白纹伊蚊miRNA在各个阶段表达量,并归一化。然后以未吸血雌蚊为参照,制1og2-rato、Scatter plot图,比较2个样品中miRNA表达差异。3.2 miRNA表达谱分析：将6个发育阶段各个miRNA求平均值,用MeV软件聚类分析各个发育阶段miRNA表达。3.3 Northtern blot验证对不同发育阶段样本及吸血前后白纹伊蚊miRNA的表达谱分析。4登革病毒感染后白纹伊蚊miRNA的Solexa分析：4.1 Solexa测序方法对白纹伊蚊感染登革病毒前后表达差异分析,寻找特异性表达miRNA。4.2 RNAhybrid软件预测针对登革病毒基因组的miRNA。研究结果：1通过对白纹伊蚊5个发育阶段混合样品的Solexa测序以及生物信息学分析,共测得序列12663936条。按照rRNAetc->repeat -> exon -> intron>known miRNA-的优先级顺序将小RNA遍历,分别是5094481、986、51020、307750、2368768。进一步鉴定79种已知miRNAs,6种新的miRNAs。2 Northern blot验证miR-184在各个期均有表达,而miR-2941仅在虫卵中表达,证明Solexa高通量测序方法的可靠性。3用上述建立的Solexa高通量测序方法对白纹伊蚊虫卵、幼虫、蛹、雌雄成虫及吸血后的雌蚊6个发育阶段miRNAs组学分析,共获得已知91种miRNAs,37个新miRNAs。4将获得91种已知miRNAs进行种系进化分析,发现白纹伊蚊与昆虫类miRNA在进化上比脊椎动物更为保守,其中与库蚊miRNA同源性最高。5 37个新miRNA是蚊虫特有的,进一步与冈比亚蚊、致卷库蚊基因组序列比对,其中miR-B-m0097, miR-F-m0156和miR-M-m0019在按蚊基因组序列中发现,在库蚊中除了上述miRNA外,miR-2941、miR-2943、miR-M-m0036等10种miRNA在库蚊发现,其余26个均为伊蚊特异miRNA.6 128个miRNA中19个pre-miRNAs是高度相似或相同的,如miR-137-1和miR-137-2、miR-263a和miR-263b、miR-2941-1和miR-2941-2等。7共发现10个白纹伊蚊基因簇,其中8个基因簇各个簇的miRNA距离小于lkb。基因簇aal-miR-306b、aal-miR-79、aal-miR-9b共同位于染色体CH477970.1,已经明确编码基因功能即丝氨酸苏氨酸激酶的内含子上miRNA簇。8 128miRNA在白纹伊蚊各个生长发育阶段聚类分析显示未成虫和成虫2个大类,其中未成虫中包括虫卵、幼虫和蛹,其中蛹和幼虫阶段的miRNA表达特征上更为接近。从纵向miRNA在各个阶段表达分析,共有虫卵特异性、幼虫特异性、蛹特异性、雄虫特异性、雌虫特异性及吸血后雌虫+幼虫特异性6个簇。其中虫卵和雄虫显示明显期特异性表达。31miRNAs在虫卵中呈高表达,如miR-315、miR-2941、miR-2943、miR-2944、miR-2945、miR-E-m0165等。54miRNAs在雄虫中呈高表达,其中某些miRNA在雌虫中也呈高表达如let-7、miR-210等。9 Northern blot对6种保守性1niRNAs与6种蚊虫特异性miRNAs进行验证。所有的白纹伊蚊miRNAs至少在一个发育阶段出现20nt左右的条带,并且除了个别miRNAs表达略有不同外,其与miRNAs表达谱与Solexa测序结果高度相似。10 Solexa高通量分析吸血前后的miRNAs表达差异分析,大部分miRNAs在吸血前后表达量保持一致,miR-263、miR-283、miR-308、miR、miR-92、miR-998、miR-989、miR-2941、miR-2944、miR-184吸血后表达量增加,miR-B-m003、miR-7、miR-993、miR-957在吸血后表达量降低。进一步northern blot分析4种miRNAs (miR-184、miR-998、miR-989、miR-2941)均在吸血后表达水平增加,与Solexa结果高度吻合。11经Solexa测序,感染和对照样本中均可检测到91种保守1niRNAs和37种新miRNAs。白纹伊蚊经登革病毒感染后miRNA表达谱差异显著,其中38个白纹伊蚊miRNAs经登革病毒感染后的表达下调,21个miRNAs经登革病毒感染后表达上调。12用RNAhybrid软件和SMVLight软件筛选登革病毒靶基因,结合表达谱差异,共获得23种登革病毒感染后差异表达的miRNAs,其靶基因分别针对登革病毒基因组非结构蛋白(1、2、4、5)和包膜蛋白基因等。研究结论：1首次用Solexa高通量测序技术并结合生物信息学方法对白纹伊蚊虫卵、幼虫、蛹、雌雄成虫及吸血后的雌蚊6个生长发育阶段miRNAs组学进行系统性分析,丰富蚊虫miRNA,并揭示新的蚊虫特异的miRNAs。2白纹伊蚊系统性表达谱分析表明白纹伊蚊miRNAs在不同发育阶段表达具有期特异性。二者结果高度相似性,进一步验证了Solexa高通量测序技术可作为miRNA表达谱分析的可靠性。3雌蚊吸血前后的miRNAs表达差异表明miRNA与雌蚊的生长发育和吸血有关,对其功能及作用机制的研究将有助于我们对蚊虫繁殖、发育及疾病传播的进一步了解。4经登革病毒感染后白纹伊蚊miRNA表达谱差异显著,提示miRNA可能再在抗病毒、抗感染中起着重要的调控作用。为研究登革病毒感染调控机制奠定基础。