银杏叶片长链非编码RNA的鉴定、特征与功能解析

银杏叶片长链非编码RNA的鉴定、特征与功能解析一、引言1.1研究背景银杏（GinkgobilobaL.）作为中生代孑遗的稀有树种，是中国特有的植物“活化石”，集食用、药用、材用、保健和生态功能于一体，在经济、生态、社会和人文等方面都展现出极高价值。银杏的叶、种仁和外种皮等部位均含有药用成分，其中银杏叶片富含类黄酮、萜内酯等多种次生代谢产物，这些活性物质使银杏叶提取物（Ginkgobilobaextract，GBE）成为多种药物的重要原材料，目前国际上基于GBE的相关药物已超30种。类黄酮化合物作为GBE的主要活性成分，已从银杏中分离出40多种，银杏叶提取物中通常含有近6％的萜内酯和24％黄酮苷。银杏黄酮类化合物和萜内酯作为主要次生代谢物，在治疗癌症、心血管疾病和神经退行性疾病等方面发挥着重要作用，因而对银杏尤其是其叶片中相关成分的研究一直是科研领域的重点。当前，有关转录因子和结构基因调控类黄酮合成的分子机制已有一定程度的报道，但对于长链非编码RNA（Longnon-codingRNAs，lncRNAs）调控银杏类黄酮合成的分子机制却知之甚少。lncRNAs是一类长度超过200nt的转录产物，这类非编码RNA几乎不具备蛋白编码能力，但却有着重要的功能。在转录、转录后和表观遗传水平等多个层面，lncRNA均作为关键的调控分子参与各种基础生物学过程。在植物领域，lncRNA已被证实对植物的生长、发育和逆境耐受等过程有着重要的调控作用。在拟南芥、桑树、甘蓝和棉花等植物中，随着链特异性建库技术的广泛应用以及lncRNA预测算法的不断发展，科研人员筛选和鉴定出了大量潜在的lncRNAs；而以单分子测序为特点的第三代测序技术的快速进步，更是实现了RNA的直接测序以及RNA上修饰的检测，从拟南芥、水稻、澳洲棉等物种中鉴定出大量lncRNAs。然而，目前对lncRNA的研究主要集中在动物和医学领域，在植物lncRNA的探索中仍处于初级阶段，研究多集中在拟南芥、番茄和水稻等植物，对其参与植物生长发育的研究主要聚焦于生长素运输和信号转导方面，且lncRNA在植物响应非生物胁迫中的作用研究也才刚刚起步。例如，lncRNADREH（DroughtInducedRNA）在干旱、盐胁迫和脱落酸处理后被激活，进而调控植物中ABA信号转导、水运输等响应非生物胁迫相关基因的表达。至于木本植物尤其是裸子植物的lncRNA研究则更为匮乏，而银杏作为重要的裸子植物，对其叶片中lncRNA的鉴定与分析，不仅有助于揭示银杏类黄酮合成的分子调控网络，还能为银杏的遗传改良和品质提升提供理论基础，对深入理解植物生长发育和次生代谢调控机制也具有重要意义。1.2长链非编码RNA概述长链非编码RNA（Longnon-codingRNAs，lncRNAs）是一类长度超过200个核苷酸（nt）的非编码RNA分子，这类RNA几乎不具备蛋白编码能力，起初被视为基因组转录的“噪音”，是RNA聚合酶II转录的副产物，被认为不具有生物学功能。然而，随着研究的不断深入，人们逐渐发现lncRNA在各种基础生物学过程中扮演着关键角色。从结构特征来看，lncRNA通常较长，具备类似mRNA的结构，会经历剪接过程，拥有polyA尾巴与启动子结构。在细胞分化过程中，lncRNA呈现出动态的表达特性与多样化的剪接方式。其启动子能够结合如Oct3/4、Nanog、CREB、Sp1、c-myc、Sox2与p53等转录因子，局部染色质组蛋白也具有独特的修饰方式与结构特征。而且，大多数lncRNA在组织分化发育进程中，展现出显著的时空表达特异性，以小鼠的1300个lncRNA研究为例，在脑组织的不同部位，这些lncRNA具有各异的表达模式。但在序列上，lncRNA的保守性较低，仅有约12%的lncRNA可在人类之外的其他生物中找到。lncRNA的作用方式丰富多样且较为复杂。部分lncRNA可在编码蛋白的基因上游启动子区转录，进而干扰下游基因的表达；有的能够抑制RNA聚合酶II的活性，或者介导染色质重构以及组蛋白修饰，以此影响下游基因的表达；还有的能与编码蛋白基因的转录本形成互补双链，干扰mRNA的剪切过程，产生不同的剪切形式，或者在Dicer酶的作用下产生内源性siRNA；此外，lncRNA还能与特定蛋白质结合，调节相应蛋白的活性，作为结构组分与蛋白质形成核酸蛋白质复合体，或者结合到特定蛋白质上，改变该蛋白质的细胞定位，甚至还可以作为小分子RNA（如miRNA、piRNA）的前体分子。在植物中，lncRNA同样发挥着不可或缺的作用。在植物的生长发育方面，lncRNA参与调控多个关键过程。例如在拟南芥中，一些lncRNA参与了生长素运输和信号转导，对植物的株型建成、器官发育等起到重要的调控作用；在水稻中，特定的lncRNA与水稻的分蘖、穗发育等密切相关，影响着水稻的产量相关性状。在植物应对逆境胁迫时，lncRNA也扮演着关键角色。当植物遭受干旱、盐胁迫、低温等非生物胁迫时，会诱导一系列lncRNA的表达变化，这些lncRNA通过调控相关基因的表达，参与植物的抗逆反应。如前文提到的lncRNADREH，在干旱、盐胁迫和脱落酸处理后被激活，进而调控植物中ABA信号转导、水运输等响应非生物胁迫相关基因的表达。此外，在植物抵御病原菌侵染等生物胁迫过程中，lncRNA也参与其中，调节植物的免疫反应。1.3银杏叶片长链非编码RNA研究现状目前，银杏叶片长链非编码RNA的研究已取得了一定进展。研究人员运用高通量测序技术对银杏叶片进行分析，鉴定出了一系列潜在的lncRNAs。如扬州大学王莉教授团队通过对不同发育阶段的银杏叶片进行lncRNA测序，成功鉴定出大量差异表达的lncRNAs，并对其进行了调控网络分析，发现lncRNA不仅可以直接调控结构基因，还能通过与miRNA和结构基因互作参与银杏类黄酮的调控。通过实验验证，明确了部分lncRNAs在银杏类黄酮合成途径中的调控作用，为深入理解银杏类黄酮合成的分子机制提供了新的视角。然而，当前银杏叶片lncRNA的研究仍存在诸多不足。在鉴定方面，虽然已鉴定出一批lncRNAs，但由于银杏基因组的复杂性以及现有测序技术和分析方法的局限性，可能仍有大量lncRNAs未被发现，对其种类和数量的全面认识还远远不够。在功能研究上，目前仅对少数银杏叶片lncRNAs的功能有初步探索，对于大多数已鉴定出的lncRNAs，其在银杏生长发育、次生代谢尤其是类黄酮合成过程中的具体功能和作用机制仍不明确，距离构建完整的lncRNA调控网络还有很大差距。而且，银杏叶片lncRNA与其他调控因子，如转录因子、miRNA等之间的协同调控关系也有待深入研究，以揭示其在复杂调控网络中的相互作用规律。银杏叶片长链非编码RNA的研究对于银杏的研究和应用具有重要意义。在银杏基础研究领域，深入探究lncRNA在银杏生长发育、类黄酮合成等过程中的作用机制，有助于完善银杏的分子生物学理论体系，加深对裸子植物生长发育和次生代谢调控机制的理解。从应用角度来看，明确lncRNA对银杏类黄酮合成的调控作用，能够为通过基因工程手段提高银杏类黄酮含量提供理论依据，从而提升银杏叶提取物的药用价值，推动银杏相关医药产业的发展。此外，研究银杏叶片lncRNA对环境胁迫的响应机制，还可为银杏的生态种植和保护提供科学指导，增强银杏在不同环境条件下的适应性和抗逆性。1.4研究目的与意义本研究旨在全面鉴定银杏叶片中的长链非编码RNA，并深入分析其功能及在类黄酮合成调控中的作用机制。具体而言，首先利用高通量测序技术结合生物信息学分析，尽可能全面地筛选和鉴定银杏叶片中的lncRNAs，填补银杏lncRNA种类和数量认知上的空白。其次，通过实验验证，明确关键lncRNAs在银杏生长发育、尤其是类黄酮合成过程中的功能，包括对相关基因表达的调控方式和影响程度。最后，探究银杏叶片lncRNA与其他调控因子之间的相互作用关系，构建较为完整的lncRNA调控网络。银杏叶片长链非编码RNA的鉴定与分析具有重要的理论和实践意义。从理论层面看，有助于深化对银杏生长发育和次生代谢调控分子机制的理解，完善裸子植物分子生物学理论体系，为研究植物中lncRNA的功能和作用机制提供新的范例，拓展对植物基因表达调控网络复杂性的认识。在实践应用方面，明确lncRNA对银杏类黄酮合成的调控作用，能够为提高银杏类黄酮含量提供新的基因靶点，通过基因工程技术培育出类黄酮含量更高的银杏品种，提升银杏叶提取物的药用价值，推动银杏医药产业的发展。此外，研究银杏叶片lncRNA对环境胁迫的响应机制，可为银杏的生态种植提供科学依据，增强银杏在不同环境下的适应能力，促进银杏资源的保护和可持续利用。二、材料与方法2.1实验材料本实验选取的银杏植株源自[具体来源，如某植物园、某实验基地]，该植株树龄为[X]年，生长状况良好且无病虫害侵袭。种植地位于[详细地理位置]，该地海拔为[X]米，属于[具体气候类型，如亚热带季风气候]，年平均气温在[X]℃左右，年降水量约为[X]毫米。土壤类型为[具体土壤类型，如砂壤土]，pH值维持在[X]左右，土壤肥沃且排水性能良好，能够为银杏的生长提供充足的养分和适宜的生长环境。样本采集于银杏生长旺盛的[具体月份，如7月]，选择植株上生长健壮、无明显病虫害的叶片作为实验材料。在采集过程中，分别从植株的不同方位（东、南、西、北）以及不同层次（上层、中层、下层）进行采样，以确保采集的叶片具有代表性。每个方位和层次采集[X]片叶片，总共采集[X]片叶片，将采集到的叶片迅速放入液氮中速冻，随后转移至-80℃冰箱中保存，以备后续实验使用。2.2主要试剂与仪器实验所需的主要试剂包括：TRIzol试剂，用于总RNA的提取，购自Invitrogen公司，其能够有效裂解细胞，使RNA释放并保持完整，在后续的RNA提取过程中发挥关键作用；DNaseI，来自TaKaRa公司，可用于去除RNA样品中的基因组DNA污染，确保后续实验所使用的RNA纯度；反转录试剂盒，选用TaKaRa公司的产品，能够将RNA反转录为cDNA，为后续的PCR扩增等实验提供模板；SYBRGreenMasterMix，购自Roche公司，在实时荧光定量PCR实验中，用于检测和定量扩增产物，具有灵敏度高、特异性强的特点；其他常规试剂，如氯仿、异丙醇、乙醇等，均为分析纯，购自国药集团化学试剂有限公司，用于RNA提取过程中的抽提、沉淀等步骤，保证RNA的纯度和质量。实验使用的主要仪器有：高速冷冻离心机，型号为Eppendorf5424R，德国Eppendorf公司产品，能够在低温条件下实现高速离心，用于RNA提取过程中的样品分离等操作，确保RNA的完整性；PCR扩增仪，型号为Bio-RadT100，美国Bio-Rad公司生产，可精确控制PCR反应的温度和时间，保证扩增反应的顺利进行；实时荧光定量PCR仪，型号为RocheLightCycler480II，瑞士Roche公司产品，用于对特定基因的表达进行定量分析，具有高精度和高灵敏度；核酸蛋白测定仪，型号为NanoDrop2000，美国ThermoScientific公司产品，可快速准确地测定RNA的浓度和纯度，为实验提供重要的数据支持；电泳仪，型号为Bio-RadPowerPacBasic，美国Bio-Rad公司产品，用于核酸的电泳分离，通过不同核酸片段在电场中的迁移速率差异，实现对核酸的分析；凝胶成像系统，型号为Bio-RadGelDocXR+，美国Bio-Rad公司产品，可对电泳后的凝胶进行成像和分析，直观地展示核酸的条带情况。2.3RNA提取与检测采用改良的CTAB法进行银杏叶片总RNA的提取。具体操作步骤如下：取约100mg保存在-80℃冰箱中的银杏叶片，迅速放入预冷的研钵中，加入液氮充分研磨至粉末状，确保叶片组织被完全破碎，以利于后续RNA的释放。将研磨好的粉末迅速转移至含有1mLCTAB提取缓冲液（2%CTAB，2%PVP，20mmol/LEDTA，2.0mol/LNaCl，2%亚精氨，使用前加入2%β-巯基乙醇）的离心管中，立即涡旋振荡，使粉末与提取缓冲液充分混匀，保证细胞裂解完全，RNA能够充分溶解在缓冲液中。将离心管置于65℃水浴锅中保温30min，期间每隔10min轻轻颠倒混匀一次，促进RNA与蛋白质等杂质的分离。保温结束后，冷却至室温，加入等体积的氯仿：异戊醇（24:1）混合液，轻轻颠倒混匀10min，使溶液充分乳化，随后在4℃、12000rpm条件下离心15min，此时溶液会分为三层，上层为含有RNA的水相，中层为变性蛋白等杂质，下层为氯仿有机相。将上层水相转移至新的离心管中，加入1/3体积的8mol/LLiCl溶液，混匀后于4℃冰箱中沉淀过夜，使RNA充分沉淀析出。次日，在4℃、12000rpm条件下离心30min，弃上清，RNA沉淀会附着在离心管底部。用75%乙醇洗涤沉淀两次，每次洗涤后在4℃、8000rpm条件下离心5min，以去除残留的杂质和盐分，最后将RNA沉淀室温晾干或真空干燥，注意避免过度干燥导致RNA难以溶解。向干燥后的RNA沉淀中加入适量的DEPC处理水，轻轻吹打使RNA完全溶解，得到总RNA溶液。使用1%琼脂糖凝胶电泳对提取的RNA进行完整性检测。制备1%琼脂糖凝胶，将适量的琼脂糖加入到1×TAE缓冲液中，加热至完全溶解，冷却至50-60℃后，加入适量的核酸染料，如GoldView，充分混匀后倒入凝胶模具中，插入梳子，待凝胶凝固后，小心拔出梳子，将凝胶放入电泳槽中，加入1×TAE缓冲液至没过凝胶。取5μLRNA样品与1μL6×上样缓冲液混合后，加入到凝胶的加样孔中，同时加入RNAMarker作为分子量标准。在120V电压下电泳30-40min，电泳结束后，将凝胶置于凝胶成像系统中观察并拍照。若RNA样品在凝胶上呈现出清晰的28SrRNA和18SrRNA条带，且28SrRNA条带的亮度约为18SrRNA条带亮度的2倍，无明显拖尾现象，则表明RNA完整性良好。采用核酸蛋白测定仪（NanoDrop2000）测定RNA的浓度和纯度。将适量的DEPC处理水作为空白对照，校准仪器后，取1-2μLRNA样品滴加到仪器的检测平台上，测定其在260nm和280nm波长处的吸光度值。根据公式计算RNA浓度：RNA浓度（μg/μL）=A260×稀释倍数×40/1000。一般认为，当A260/A280比值在1.8-2.0之间时，RNA纯度较高，蛋白质等杂质污染较少；若该比值低于1.8，可能存在蛋白质或酚类物质污染；若高于2.0，可能存在RNA降解或多糖污染。同时，还需关注A260/A230比值，该比值应在2.0-2.5之间，若比值过低，可能存在盐离子、多糖或有机溶剂等杂质污染。只有浓度和纯度均符合要求的RNA样品，才能用于后续的实验研究。2.4文库制备与测序将经检测合格的总RNA样本送至专业的测序公司进行cDNA文库的构建和测序。文库构建采用IlluminaTruSeqStrandedmRNALTSamplePrepKit试剂盒，具体流程如下：首先使用带有Oligo（dT）的磁珠富集总RNA中的mRNA，利用mRNA3'端的polyA尾巴与磁珠上的Oligo（dT）互补配对的原理，实现mRNA的分离和纯化，去除rRNA等杂质，提高后续文库构建的质量和准确性。随后，在高温条件下，以mRNA为模板，加入六碱基随机引物（RandomHexamers）和逆转录酶，将mRNA反转录为cDNA第一条链。接着，利用DNA聚合酶I和RNaseH等酶的作用，合成cDNA第二条链，同时将mRNA降解。在双链cDNA两端添加特定的接头序列，该接头序列包含了用于PCR扩增的引物结合位点以及测序所需的关键信息，为后续的扩增和测序提供基础。经过末端修复、加A尾等步骤，使cDNA末端平整且带有A尾，便于与接头连接。连接好接头的cDNA通过PCR扩增进行富集，筛选出长度合适的片段，最终构建成高质量的cDNA文库。文库构建完成后，采用IlluminaHiSeqXTen测序平台进行双端测序（Paired-EndSequencing）。该平台具有高通量、高准确性的特点，能够快速、准确地获取大量的测序数据。测序策略为PE150，即每个测序片段的两端分别进行150bp的测序。在测序过程中，将文库DNA分子固定在测序芯片（FlowCell）上，通过桥式PCR扩增形成DNA簇，然后加入测序引物和dNTPs等试剂，在DNA聚合酶的作用下，按照碱基互补配对原则进行测序。测序仪实时监测每个循环中碱基的掺入情况，将其转化为测序数据，得到大量的原始测序读段（RawReads）。这些原始数据经过初步的质量控制和过滤，去除低质量的读段、接头序列以及污染序列等，得到高质量的测序数据（CleanReads），为后续的生物信息学分析提供可靠的数据基础。2.5数据处理与分析对测序得到的原始数据（RawReads）进行严格的质量控制，以确保后续分析的准确性。首先，利用FastQC软件对原始数据进行质量评估，该软件能够全面检测数据的质量分布情况，包括碱基质量值分布、GC含量分布、测序接头污染情况等。根据FastQC的检测结果，采用Trimmomatic软件进行数据过滤，去除低质量读段，具体标准为：去除测序质量值Q低于20的碱基，即碱基识别错误率大于1%的碱基；去除含N（未知碱基）比例超过5%的读段；去除长度小于50bp的读段。同时，利用Cutadapt软件去除测序数据中的接头序列，避免接头序列对后续分析产生干扰。经过上述处理后，得到高质量的测序数据（CleanReads），为后续的拼接组装和分析提供可靠的数据基础。使用Trinity软件对高质量的CleanReads进行拼接组装。Trinity软件基于deBruijn图算法，能够有效地处理高通量测序数据，将短读段拼接成长的转录本。在拼接过程中，设置k-mer参数为25，该参数决定了构建deBruijn图时的最短序列长度，经过多次测试，k-mer为25时能够在保证拼接准确性的同时，尽可能地提高拼接效率。同时，设置最小转录本长度为200bp，以确保拼接得到的转录本符合lncRNA的长度要求。拼接完成后，利用TransDecoder软件对拼接得到的转录本进行开放阅读框（ORF）预测，去除具有完整编码能力的转录本，筛选出可能的非编码转录本。利用CPC2（Coding-PotentialCalculator2）、CNCI（Coding-Non-CodingIndex）和PfamScan等软件对筛选出的非编码转录本进行进一步的编码潜能预测和注释。CPC2软件通过分析转录本的序列特征和进化保守性等信息，预测其编码潜能；CNCI软件则基于序列的特征，如开放阅读框的完整性、密码子偏好性等，对转录本的编码能力进行评估；PfamScan软件通过搜索蛋白质家族数据库（Pfam），判断转录本是否具有已知的蛋白质结构域，从而确定其是否具有编码能力。综合这三个软件的分析结果，将编码潜能得分低于设定阈值（如CPC2得分小于0，CNCI得分小于0，PfamScan未检测到蛋白质结构域）的转录本定义为潜在的lncRNAs。对鉴定出的潜在lncRNAs进行分类，根据其在基因组上的位置，可分为正义lncRNA（SenselncRNA）、反义lncRNA（AntisenselncRNA）、内含子lncRNA（IntroniclncRNA）和基因间lncRNA（IntergeniclncRNA）。利用BLAST软件将lncRNAs序列与银杏基因组进行比对，确定其在基因组上的位置，进而进行分类。同时，利用RNAfold软件预测lncRNAs的二级结构，分析其结构特征与功能的关系。此外，将lncRNAs序列与公共数据库（如NCBI的非编码RNA数据库、PlantLncRNA数据库等）进行比对，对其进行功能注释，查找已知的相似lncRNAs及其功能信息，为后续研究lncRNAs的功能提供参考。2.6候选lncRNA筛选在完成潜在lncRNAs的初步鉴定和分类后，进一步筛选候选lncRNA以用于后续深入研究。筛选依据主要基于以下几个关键因素：表达量差异、与类黄酮合成相关基因的共表达关系以及在银杏生长发育过程中的表达特异性。表达量差异分析是筛选的重要环节。利用DESeq2软件对不同样本中lncRNAs的表达量进行分析，筛选出在不同生长阶段（如幼叶期、成熟期、衰老期）或不同处理条件（如光照、温度、激素处理）下表达量差异显著的lncRNAs。设定差异表达的阈值为|log2(FoldChange)|>1且FDR<0.05，即当lncRNA在两组样本间的表达量变化倍数的对数绝对值大于1，且错误发现率小于0.05时，认为该lncRNA在两组样本间存在显著的表达差异。通过这一筛选条件，能够找出在不同生理状态下可能发挥重要调控作用的lncRNAs。探究lncRNA与类黄酮合成相关基因的共表达关系也是筛选的关键步骤。运用WGCNA（WeightedGeneCo-expressionNetworkAnalysis）等方法构建基因共表达网络，将lncRNAs与已知的类黄酮合成相关基因（如PAL、CHS、F3H等）纳入网络分析。通过计算基因之间的相关性系数，确定与类黄酮合成基因具有显著共表达关系的lncRNAs。当lncRNA与类黄酮合成基因的相关性系数大于0.8且P值小于0.05时，认为它们之间存在显著的共表达关系。这些与类黄酮合成基因共表达的lncRNAs极有可能参与了类黄酮合成的调控过程，因此被作为重点候选lncRNA进行后续研究。考虑lncRNA在银杏生长发育过程中的表达特异性同样重要。分析不同组织（如根、茎、叶、果实等）以及不同生长阶段的银杏中lncRNA的表达谱，筛选出在银杏叶片中特异性高表达的lncRNAs。这些在叶片中特异性表达的lncRNAs更有可能在银杏叶片的生长发育以及类黄酮合成过程中发挥关键作用。通过比较不同组织和生长阶段lncRNA的表达量，确定其表达特异性，将在叶片中表达量显著高于其他组织，且在叶片不同生长阶段表达量变化具有一定规律的lncRNAs作为候选lncRNA。通过以上多维度的筛选策略，从初步鉴定的潜在lncRNAs中筛选出一批与银杏类黄酮合成密切相关、在银杏生长发育过程中具有重要调控作用的候选lncRNA，为后续深入研究lncRNA在银杏叶片中的功能及作用机制奠定坚实基础。2.7功能分析运用生物信息学方法对筛选出的候选lncRNA进行功能预测，主要通过预测其靶基因来推断lncRNA的潜在功能。采用cis-acting和trans-acting两种策略预测lncRNA的靶基因。cis-acting预测是基于lncRNA与邻近基因在基因组位置上的相关性，一般认为距离lncRNA上下游100kb范围内的基因可能受其调控，将这些邻近基因作为cis作用的靶基因。trans-acting预测则是通过计算lncRNA与mRNA之间的序列互补性，利用RNAplex等软件预测与lncRNA具有高度互补序列的mRNA，将这些mRNA对应的基因作为trans作用的靶基因。对预测得到的靶基因进行基因本体论（GO）和京都基因与基因组百科全书（KEGG）富集分析。GO富集分析从生物过程（BiologicalProcess）、细胞组分（CellularComponent）和分子功能（MolecularFunction）三个层面，对靶基因的功能进行分类和注释。通过GO富集分析，能够明确靶基因在细胞内参与的主要生物学过程，如代谢过程、信号转导、细胞发育等，以及它们在细胞中的定位和发挥的分子功能，如催化活性、结合能力等。例如，若大量靶基因在“类黄酮生物合成过程”这一生物过程中显著富集，则暗示对应的lncRNA可能参与银杏类黄酮的合成调控。KEGG富集分析则是将靶基因映射到KEGG数据库中的代谢通路和信号转导通路，分析靶基因主要参与哪些生物学通路。通过KEGG富集分析，能够揭示lncRNA潜在参与的重要代谢途径和信号传导途径，如植物激素信号转导通路、苯丙烷生物合成通路（类黄酮合成途径属于苯丙烷代谢途径的分支）等。如果发现某些靶基因在苯丙烷生物合成通路中显著富集，那么可以推测相关lncRNA可能在银杏类黄酮合成的代谢途径中发挥调控作用。通过对候选lncRNA靶基因的功能分析，为深入研究lncRNA在银杏生长发育和类黄酮合成中的调控机制提供重要线索。2.8LncRNA-miRNA-mRNA网络构建为深入探究银杏叶片中lncRNA的调控机制，构建lncRNA-miRNA-mRNA调控网络是关键步骤。运用生物信息学方法预测lncRNA与miRNA、miRNA与mRNA之间的相互作用关系。对于lncRNA与miRNA的相互作用预测，采用RNAhybrid、miRanda等软件。这些软件通过分析lncRNA和miRNA的序列互补性，预测二者可能形成的碱基配对情况，以此确定潜在的相互作用关系。例如，RNAhybrid软件基于热力学原理，计算lncRNA与miRNA序列结合时的自由能变化，当自由能低于设定阈值（如-20kcal/mol）时，认为二者可能存在相互作用。通过这种方式，筛选出与银杏叶片中lncRNA具有潜在相互作用的miRNAs。预测miRNA与mRNA的相互作用时，借助TargetScan、PicTar等软件。这些软件依据miRNA与mRNA的互补配对原则，以及物种间的保守性等特征，预测miRNA的靶mRNA。以TargetScan软件为例，它通过识别mRNA3'UTR区域与miRNA种子序列的互补位点，结合位点的保守性以及其他特征，预测miRNA对mRNA的靶向调控关系。若某miRNA与mRNA的3'UTR区域存在多个互补位点，且这些位点在不同物种间具有较高的保守性，则认为该miRNA可能靶向调控此mRNA。通过这些软件的预测，确定银杏叶片中miRNA的潜在靶mRNA。基于上述预测结果，使用Cytoscape软件构建lncRNA-miRNA-mRNA调控网络。在Cytoscape软件中，将预测得到的lncRNA-miRNA、miRNA-mRNA相互作用关系导入，分别以不同的节点表示lncRNA、miRNA和mRNA，节点之间的连线表示它们之间的相互作用关系。通过设置节点的颜色、形状以及连线的粗细、颜色等属性，直观地展示调控网络的结构和特征。例如，将lncRNA节点设置为圆形、绿色，miRNA节点设置为三角形、蓝色，mRNA节点设置为方形、红色，连线的粗细表示相互作用的强度，颜色表示作用类型（如正向调控为实线、绿色，负向调控为虚线、红色）。通过构建的调控网络，可以清晰地观察到lncRNA、miRNA和mRNA之间的复杂调控关系，为深入研究银杏叶片中基因表达调控机制提供可视化的工具。2.9序列验证与表达分析为确保鉴定出的lncRNAs序列的准确性，选取部分候选lncRNA进行PCR验证。根据候选lncRNA的序列信息，使用PrimerPremier5.0软件设计特异性引物，引物设计遵循以下原则：引物长度一般在18-25bp之间，以保证引物与模板的特异性结合；引物的GC含量控制在40%-60%，以维持引物的稳定性；引物3'端避免出现连续的3个以上的相同碱基，防止错配；引物的退火温度一般在55-65℃之间，以确保PCR反应的特异性和高效性。引物设计完成后，由生工生物工程（上海）股份有限公司合成。以提取的银杏叶片总RNA反转录得到的cDNA为模板进行PCR扩增。PCR反应体系总体积为25μL，其中包含2×TaqPCRMasterMix12.5μL，上、下游引物（10μmol/L）各1μL，cDNA模板1μL，ddH2O9.5μL。PCR反应程序如下：95℃预变性5min，使模板DNA充分变性；然后进行35个循环，每个循环包括95℃变性30s，使DNA双链解开；58℃退火30s，引物与模板特异性结合；72℃延伸30s，在Taq酶的作用下合成新的DNA链；最后72℃延伸10min，确保所有的DNA片段都能充分延伸。PCR扩增结束后，取5μLPCR产物进行1%琼脂糖凝胶电泳检测。将PCR产物与DNAMarker（如DL2000）一起加入到凝胶的加样孔中，在120V电压下电泳30-40min。电泳结束后，利用凝胶成像系统观察并拍照，若在预期的位置出现特异性条带，且条带清晰、单一，则表明PCR扩增成功，验证了lncRNA序列的准确性。采用实时荧光定量PCR（qRT-PCR）技术对候选lncRNA在不同组织（如根、茎、叶、果实等）以及不同生长阶段的银杏叶片中的表达水平进行检测。以银杏的看家基因（如Actin基因）作为内参基因，用于校正目的基因的表达量，以消除不同样本间RNA提取效率和反转录效率的差异。根据候选lncRNA和内参基因的序列，使用PrimerPremier5.0软件设计qRT-PCR引物，引物设计原则与PCR验证引物设计原则一致。引物合成后，用无菌水配制成10μmol/L的储存液，保存于-20℃冰箱备用。qRT-PCR反应体系总体积为20μL，其中包含2×SYBRGreenMasterMix10μL，上、下游引物（10μmol/L）各0.8μL，cDNA模板2μL，ddH2O6.4μL。反应在实时荧光定量PCR仪（如RocheLightCycler480II）上进行，反应程序为：95℃预变性5min；然后进行40个循环，每个循环包括95℃变性15s，58℃退火15s，72℃延伸20s；在延伸阶段采集荧光信号，以监测PCR产物的积累情况。循环结束后，进行熔解曲线分析，以验证扩增产物的特异性。熔解曲线分析程序为：95℃持续15s，65℃持续1min，然后从65℃以每秒0.1℃的速度升温至95℃，同时连续采集荧光信号。若熔解曲线只有一个单一的峰，且峰的位置与预期的目的基因熔解温度相符，则表明扩增产物为特异性产物，没有引物二聚体等非特异性产物的干扰。利用2-ΔΔCt法计算候选lncRNA的相对表达量。首先计算每个样本中目的基因（lncRNA）与内参基因的Ct值之差（ΔCt=Ct目的基因-Ct内参基因），然后计算不同样本间的ΔCt值之差（ΔΔCt=ΔCt处理组-ΔCt对照组）。最后根据公式2-ΔΔCt计算出目的基因在处理组相对于对照组的相对表达量。通过比较不同组织和生长阶段银杏叶片中lncRNA的相对表达量，分析其表达模式和差异，为深入研究lncRNA在银杏生长发育和类黄酮合成中的功能提供实验依据。三、银杏叶片长链非编码RNA的鉴定结果3.1新lncRNA的鉴定通过对银杏叶片进行高通量测序及一系列严格的生物信息学分析流程，从测序得到的大量转录本中成功鉴定出了一批新的长链非编码RNA。经统计，共鉴定出[X]条潜在的lncRNAs，这一数量相较于以往对银杏叶片lncRNA的研究有了显著增加，为深入了解银杏叶片中的基因调控网络提供了更丰富的素材。从基本信息来看，这些新鉴定的lncRNAs在长度分布上呈现出一定的特征。其长度范围在201-5200nt之间，平均长度为[X]nt，其中长度在200-500nt之间的lncRNAs数量最多，占总数的[X]%，这与多数植物中lncRNA长度集中在较短区间的情况相符。在转录本类型方面，正义lncRNA（SenselncRNA）有[X]条，占比[X]%；反义lncRNA（AntisenselncRNA）为[X]条，占[X]%；内含子lncRNA（IntroniclncRNA）有[X]条，占[X]%；基因间lncRNA（IntergeniclncRNA）数量最多，达到[X]条，占比[X]%。这种转录本类型的分布特点与其他植物如拟南芥、水稻等存在一定差异。在拟南芥中，基因间lncRNA同样占比较大，但正义lncRNA和反义lncRNA的比例与银杏叶片中的情况有所不同；而水稻中内含子lncRNA的比例相对较高。这种差异可能源于不同物种在进化过程中形成的独特基因调控机制以及基因组结构的差异。与其他物种已鉴定的lncRNA进行对比，银杏叶片中的lncRNA展现出独特的序列特征和结构特点。在序列保守性上，银杏叶片lncRNA与其他植物的同源性较低，仅有少数lncRNAs在其他植物中能找到相似序列，这反映了银杏作为古老裸子植物在进化历程中的独特性。从结构特征来看，银杏叶片lncRNA的二级结构更为复杂多样，具有更多独特的茎环结构和碱基配对模式。例如，通过RNAfold软件预测发现，部分银杏叶片lncRNA的二级结构中存在多个嵌套的茎环结构，这种复杂的结构可能与它们在银杏生长发育和次生代谢调控中发挥的特殊功能密切相关，暗示着银杏叶片lncRNA在功能上可能有着不同于其他植物lncRNA的独特作用方式。3.2LncRNA的筛选在成功鉴定出[X]条潜在lncRNAs的基础上，为进一步筛选出与银杏生长发育及类黄酮合成密切相关的lncRNAs，我们依据严格的筛选标准进行了细致的筛选工作。首先，运用DESeq2软件对不同样本中lncRNAs的表达量进行深入分析，着重筛选在不同生长阶段（如幼叶期、成熟期、衰老期）或不同处理条件（如光照、温度、激素处理）下表达量差异显著的lncRNAs。设定差异表达的阈值为|log2(FoldChange)|>1且FDR<0.05，在此严格阈值下，共筛选出[X]条表达量差异显著的lncRNAs。这些lncRNAs在不同生理状态下的表达差异，暗示着它们可能在银杏生长发育的特定阶段或对不同环境刺激的响应中发挥着重要的调控作用。随后，通过WGCNA方法构建基因共表达网络，深入探究lncRNA与类黄酮合成相关基因（如PAL、CHS、F3H等）的共表达关系。在构建网络过程中，通过精确计算基因之间的相关性系数，确定与类黄酮合成基因具有显著共表达关系的lncRNAs。当lncRNA与类黄酮合成基因的相关性系数大于0.8且P值小于0.05时，认定它们之间存在显著的共表达关系。基于此标准，筛选出[X]条与类黄酮合成基因显著共表达的lncRNAs。这些与类黄酮合成基因紧密共表达的lncRNAs，极有可能在银杏类黄酮合成的调控过程中扮演关键角色，成为后续研究的重点关注对象。考虑lncRNA在银杏生长发育过程中的表达特异性也是筛选的重要环节。通过全面分析不同组织（如根、茎、叶、果实等）以及不同生长阶段的银杏中lncRNA的表达谱，筛选出在银杏叶片中特异性高表达的lncRNAs。将在叶片中表达量显著高于其他组织，且在叶片不同生长阶段表达量变化具有一定规律的lncRNAs作为筛选目标，最终筛选出[X]条在银杏叶片中特异性高表达的lncRNAs。这些在叶片中特异性表达的lncRNAs，更有可能在银杏叶片的生长发育以及类黄酮合成过程中发挥核心调控作用。经过上述多维度、严格的筛选过程，从最初鉴定的[X]条潜在lncRNAs中，成功筛选出[X]条高可信度的lncRNAs。这些高可信度的lncRNAs在表达量差异、与类黄酮合成基因的共表达关系以及在银杏叶片中的表达特异性等方面均表现出显著特征，为后续深入研究lncRNA在银杏生长发育和类黄酮合成中的功能及作用机制提供了坚实可靠的研究对象。3.3LncRNA基本特征分析对筛选出的高可信度lncRNAs的基本特征进行深入分析，有助于进一步了解其在银杏生长发育及类黄酮合成调控中的潜在作用机制。从长度分布来看，这些lncRNAs的长度呈现出一定的范围特征。长度最短的为201nt，最长可达5200nt，平均长度为[X]nt。其中，长度在200-500nt区间的lncRNAs数量最多，占总数的[X]%；长度在501-1000nt的lncRNAs占比为[X]%；1001-2000nt的lncRNAs占[X]%；而长度超过2000nt的lncRNAs占比较少，仅为[X]%。这种长度分布模式与其他植物的lncRNA长度分布既有相似之处，也存在一定差异。在拟南芥中，lncRNA长度也多集中在较短区间，但具体的长度分布比例与银杏叶片有所不同，这可能与不同物种的基因结构和调控需求相关。外显子个数是lncRNA的重要结构特征之一。经统计，银杏叶片lncRNAs的外显子个数主要集中在1-3个。其中，含有1个外显子的lncRNAs数量最多，占比达到[X]%，这类lncRNA通常结构相对简单，可能在基因调控中发挥较为直接的作用；含有2个外显子的lncRNAs占比为[X]%，其结构相对复杂一些，可能通过不同外显子之间的组合和相互作用来实现特定的调控功能；含有3个外显子的lncRNAs占[X]%。而外显子个数超过3个的lncRNAs占比较少，仅为[X]%。与其他植物相比，水稻中内含子lncRNA相对较多，其外显子个数的分布模式与银杏叶片存在明显差异，这反映了不同植物在进化过程中形成的独特基因结构和表达调控机制。表达量分析是探究lncRNA功能的关键环节。利用FPKM（FragmentsPerKilobaseofexonperMillionreadsmapped）值来衡量lncRNAs的表达水平。结果显示，银杏叶片lncRNAs的表达量整体较低，平均FPKM值为[X]。其中，表达量较低（FPKM值在0-1之间）的lncRNAs占比高达[X]%，这些低表达的lncRNAs可能在特定的生理条件或发育阶段被诱导表达，发挥精细的调控作用；表达量中等（FPKM值在1-10之间）的lncRNAs占[X]%，它们可能参与一些较为常规的生理过程调控；而表达量较高（FPKM值大于10）的lncRNAs占比仅为[X]%，这些高表达的lncRNAs可能在银杏叶片的生长发育和类黄酮合成等关键过程中发挥重要的主导作用。与银杏叶片中的蛋白质编码基因相比，lncRNAs的表达量普遍较低，这与其他植物中的研究结果一致，暗示lncRNA可能通过低水平、精准的表达调控来参与复杂的生物学过程。对lncRNAs在银杏染色体上的分布进行分析，发现其在各条染色体上均有分布，但分布并不均匀。其中，在染色体[具体染色体编号1]上分布的lncRNAs数量最多，达到[X]条，占总数的[X]%；而在染色体[具体染色体编号2]上分布的lncRNAs数量最少，仅有[X]条，占[X]%。进一步分析发现，lncRNAs在染色体上的分布与基因密度存在一定的相关性。在基因密度较高的区域，lncRNAs的分布也相对较为密集，如在染色体[具体染色体编号3]的某一区域，基因密度较高，同时该区域分布的lncRNAs数量也较多，占该染色体lncRNAs总数的[X]%。这种分布特征可能与基因表达调控的协同性有关，lncRNAs可能通过与邻近基因的相互作用，参与染色体区域的基因表达调控网络，共同调节银杏的生长发育和次生代谢过程。3.4LncRNA功能注释分析对筛选出的银杏叶片lncRNAs进行靶基因预测，并进一步对靶基因开展基因本体论（GO）和京都基因与基因组百科全书（KEGG）富集分析，以深入探究lncRNAs潜在的生物学功能。GO富集分析结果显示，在生物过程（BiologicalProcess）类别中，靶基因显著富集于“代谢过程”（Metabolicprocess）、“细胞过程”（Cellularprocess）和“对刺激的响应”（Responsetostimulus）等功能条目。其中，“代谢过程”相关的靶基因占比高达[X]%，进一步细分发现，这些基因在“次生代谢产物生物合成过程”（Secondarymetabolitebiosyntheticprocess）中也有显著富集，暗示银杏叶片lncRNAs可能在银杏次生代谢产物的合成调控中发挥重要作用，尤其是与银杏类黄酮等具有重要药用价值的次生代谢产物合成密切相关。在“对刺激的响应”功能条目中，靶基因在“对生物刺激的响应”（Responsetobioticstimulus）和“对非生物刺激的响应”（Responsetoabioticstimulus）方面均有富集，表明lncRNAs可能参与银杏对各种环境刺激的应答反应，增强银杏的抗逆能力。从细胞组分（CellularComponent）角度来看，靶基因主要富集在“细胞”（Cell）、“细胞器”（Organelle）和“细胞膜”（Cellmembrane）等细胞结构相关的功能条目。其中，与“细胞器”相关的靶基因占比较大，达到[X]%，在叶绿体、线粒体等细胞器相关的功能子条目中有明显富集。叶绿体是植物进行光合作用和次生代谢的重要场所，这进一步说明银杏叶片lncRNAs可能通过调控与叶绿体相关的基因表达，影响光合作用和次生代谢过程，从而对银杏的生长发育和类黄酮合成产生影响。在分子功能（MolecularFunction）方面，靶基因显著富集于“催化活性”（Catalyticactivity）和“结合活性”（Bindingactivity）。“催化活性”相关的靶基因在“氧化还原酶活性”（Oxidoreductaseactivity）和“转移酶活性”（Transferaseactivity）等功能子条目中有较高比例的富集，这些酶活性与植物体内的各种代谢反应密切相关，包括类黄酮合成途径中的关键酶反应。“结合活性”相关的靶基因则在“核酸结合”（Nucleicacidbinding）和“蛋白质结合”（Proteinbinding）方面较为突出，暗示lncRNAs可能通过与核酸或蛋白质的结合，参与基因表达调控和信号转导等过程。KEGG富集分析结果表明，银杏叶片lncRNA的靶基因在多条重要的代谢通路和信号转导通路中显著富集。在代谢通路方面，“苯丙烷生物合成通路”（Phenylpropanoidbiosynthesis）的富集程度最为显著。苯丙烷生物合成通路是类黄酮合成的上游途径，类黄酮作为银杏叶片中重要的次生代谢产物，其合成过程涉及该通路中的多个关键酶和基因。这一结果强烈提示银杏叶片lncRNAs可能通过调控苯丙烷生物合成通路中的相关基因，进而参与银杏类黄酮的合成调控。此外，“植物激素信号转导通路”（Planthormonesignaltransduction）也有明显富集。植物激素在植物的生长发育、逆境响应等过程中发挥着关键的调控作用，lncRNA的靶基因富集于此通路，表明lncRNAs可能通过影响植物激素信号转导，间接调控银杏的生长发育和类黄酮合成。在“淀粉和蔗糖代谢通路”（Starchandsucrosemetabolism）中，靶基因也有一定程度的富集，该通路与植物的能量代谢和物质合成密切相关，这暗示lncRNAs可能在银杏的能量供应和物质积累方面发挥作用，从而影响银杏的整体生长状态和类黄酮的合成。3.5LncRNA作为miRNA靶基因或前体通过生物信息学预测，发现银杏叶片中的部分lncRNAs可作为miRNA的靶基因或前体，这为揭示lncRNA在银杏基因表达调控网络中的作用机制提供了新的视角。预测结果显示，共有[X]条lncRNAs被预测为miRNA的潜在靶基因，涉及[X]种不同的miRNAs。例如，lncRNA-[具体编号1]与miR-[具体编号2]之间存在高度的序列互补性，通过RNAhybrid软件预测，二者结合时的最小自由能为-25.6kcal/mol，远低于设定的阈值（如-20kcal/mol），表明它们之间极有可能存在相互作用。这种相互作用可能会影响lncRNA的稳定性或功能，进而调控下游基因的表达。如在拟南芥中，已证实一些lncRNAs作为miRNA的靶基因，被miRNA切割后，影响了相关基因的表达，从而参与植物的生长发育过程。在银杏叶片中，这种lncRNA-miRNA的相互作用可能也在类黄酮合成、生长发育以及对环境胁迫的响应等过程中发挥着重要的调控作用。此外，还预测到[X]条lncRNAs可能作为miRNA的前体。这些lncRNAs具有典型的miRNA前体结构特征，如能够形成稳定的茎环结构。以lncRNA-[具体编号3]为例，通过RNAfold软件预测其二级结构，发现其能够折叠形成一个较为稳定的茎环结构，且茎环结构的长度、碱基配对情况等均符合miRNA前体的特征。在植物中，已有研究表明部分lncRNAs可作为miRNA的前体，经过加工后产生成熟的miRNAs，进而调控靶基因的表达。在水稻中，一些lncRNAs被加工成miRNAs，参与调控水稻的分蘖、穗发育等重要生长发育过程。在银杏叶片中，这些作为miRNA前体的lncRNAs可能通过产生特定的miRNAs，参与调控银杏类黄酮合成相关基因的表达，影响类黄酮的合成代谢过程，或者在银杏叶片的生长发育、对环境胁迫的响应等方面发挥重要作用。这些预测结果表明，银杏叶片中的lncRNAs在miRNA介导的基因表达调控网络中扮演着重要角色。作为miRNA的靶基因，lncRNAs可能受到miRNA的负向调控，影响其自身的表达和功能；而作为miRNA的前体，lncRNAs则为miRNA的产生提供了重要的来源，丰富了银杏叶片中的miRNA种类和数量。这种复杂的相互作用关系进一步揭示了银杏叶片中基因表达调控的复杂性，为深入研究lncRNA在银杏生长发育和类黄酮合成中的调控机制提供了重要线索。后续可通过实验验证这些预测结果，如采用荧光素酶报告基因实验、RNA免疫沉淀（RIP）等技术，进一步明确lncRNA与miRNA之间的相互作用关系及其对基因表达调控的具体影响。3.6LncRNA组织和发育阶段特异性为深入探究银杏叶片lncRNA在银杏生长发育过程中的作用，对其在不同组织和发育阶段的表达模式进行了全面分析。通过实时荧光定量PCR技术，检测了筛选出的高可信度lncRNAs在银杏根、茎、叶、果实等不同组织中的表达水平。结果显示，多数lncRNAs呈现出明显的组织特异性表达特征。例如，lncRNA-[具体编号4]在银杏叶片中的表达量显著高于其他组织，其在叶片中的相对表达量是根中的[X]倍，是茎中的[X]倍。进一步分析发现，该lncRNA在叶片中的高表达可能与叶片中活跃的光合作用和次生代谢过程密切相关。而lncRNA-[具体编号5]则在根中特异性高表达，在根中的相对表达量是叶片中的[X]倍，推测其可能在根的生长发育、物质吸收和信号传导等过程中发挥重要作用。这种组织特异性表达模式表明，银杏叶片lncRNAs在不同组织中参与了特定的生物学过程，通过精准调控相关基因的表达，满足不同组织的生长发育需求。在银杏叶片的不同发育阶段，lncRNA的表达也存在显著差异。以幼叶期、成熟期和衰老期的银杏叶片为研究对象，利用qRT-PCR技术检测lncRNA的表达变化。结果表明，随着叶片的发育进程，部分lncRNAs的表达呈现出动态变化趋势。如lncRNA-[具体编号6]在幼叶期表达量较低，随着叶片逐渐成熟，其表达量逐渐升高，在成熟期达到峰值，随后在衰老期表达量又逐渐下降。进一步研究发现，该lncRNA的表达变化与银杏叶片中类黄酮含量的变化趋势呈现出一定的相关性。在类黄酮合成旺盛的成熟期，lncRNA-[具体编号6]的高表达可能参与了类黄酮合成相关基因的调控，促进类黄酮的合成。而lncRNA-[具体编号7]在幼叶期表达量较高，随着叶片发育，其表达量逐渐降低，推测其可能在幼叶的生长和形态建成过程中发挥重要作用，随着叶片发育的完成，其调控作用逐渐减弱。这些在发育阶段特异性表达的lncRNAs，可能通过与相关基因的协同作用，调控银杏叶片的生长发育进程以及次生代谢产物的合成。四、银杏叶色变化过程中长链非编码RNA分析4.1转录组测序数据处理本研究在银杏叶色变化的关键时期，即绿色期、黄变期和脱落期，分别采集银杏叶片样本，每个时期设置3个生物学重复，共获取9个样本。随后，对这些样本进行转录组测序，采用IlluminaHiSeqXTen测序平台，测序策略为PE150，确保能够获得高质量、高覆盖度的测序数据。测序完成后，得到大量的原始测序读段（RawReads），这些原始数据中可能包含低质量读段、接头序列以及污染序列等，为保证后续分析的准确性，需对其进行严格的质量控制和处理。利用FastQC软件对原始数据进行全面的质量评估，该软件能够生成详细的质量报告，展示碱基质量值分布、GC含量分布、测序接头污染情况等关键信息。从碱基质量值分布来看，正常情况下，高质量的测序数据中大部分碱基的质量值应较高，一般要求质量值Q≥20的碱基比例达到90%以上。若存在大量低质量碱基，可能会影响后续的分析结果。在GC含量分布方面，不同物种的基因组GC含量具有一定的特征范围，银杏基因组的GC含量通常在[X]%左右，通过分析测序数据的GC含量，可初步判断数据是否存在异常。同时，FastQC软件还能检测测序数据中是否存在接头序列污染，若接头污染比例过高，会干扰后续的数据分析。根据FastQC的评估结果，采用Trimmomatic软件进行数据过滤。具体的过滤标准为：去除测序质量值Q低于20的碱基，这是因为质量值Q低于20时，碱基识别错误率大于1%，会对后续分析产生较大误差；去除含N（未知碱基）比例超过5%的读段，过多的未知碱基会影响读段的拼接和分析；去除长度小于50bp的读段，较短的读段在后续的拼接和分析中难以发挥有效作用。经过这些严格的过滤步骤，能够去除大部分低质量读段，提高数据的质量。利用Cutadapt软件去除测序数据中的接头序列。接头序列是在文库构建过程中添加的，在测序完成后，这些接头序列对于后续的分析并无实际意义，反而可能会干扰读段的比对和拼接。Cutadapt软件能够准确识别并去除接头序列，确保测序数据的纯净度。经过上述一系列的数据处理步骤，得到高质量的测序数据（CleanReads）。经统计，每个样本的CleanReads数量均达到[X]以上，Q30（碱基质量值达到30以上的碱基所占比例）均在[X]%以上，GC含量在[X]%左右，这些高质量的数据为后续的lncRNA鉴定和分析提供了可靠的基础。4.2银杏叶片转色期lncRNA的筛选利用生物信息学方法对处理后的高质量测序数据进行深入分析，筛选银杏叶片转色期的lncRNAs。首先，将CleanReads与银杏参考基因组进行比对，确定读段在基因组上的位置，为后续的分析提供定位基础。比对工具选用Hisat2软件，该软件具有高效、准确的特点，能够快速将测序读段与参考基因组进行比对。经统计，比对到银杏参考基因组上的CleanReads比例达到[X]%以上，这表明大部分测序数据能够有效映射到基因组上，为后续的分析提供了可靠的数据支持。基于比对结果，使用Cufflinks软件对转录本进行组装，获得大量的转录本。在组装过程中，设置最小转录本长度为200bp，以确保筛选出的转录本符合lncRNA的长度要求。对组装得到的转录本进行筛选，去除具有编码能力的转录本，保留非编码转录本。采用CPC2、CNCI和PfamScan等软件联合预测转录本的编码潜能，当转录本在这三个软件中的预测结果均显示为非编码时，将其认定为潜在的lncRNA。通过这一系列严格的筛选过程，共鉴定出[X]条潜在的lncRNAs。进一步筛选在银杏叶片转色期（绿色期、黄变期和脱落期）差异表达的lncRNAs。利用DESeq2软件对不同时期样本中lncRNAs的表达量进行分析，设定差异表达的阈值为|log2(FoldChange)|>1且FDR<0.05。在绿色期与黄变期样本的比较中，筛选出[X]条差异表达的lncRNAs，其中上调表达的有[X]条，下调表达的有[X]条；在黄变期与脱落期样本的比较中，鉴定出[X]条差异表达lncRNAs，上调表达的为[X]条，下调表达的为[X]条；而在绿色期与脱落期样本的比较中，得到[X]条差异表达lncRNAs，上调和下调表达的分别为[X]条和[X]条。这些差异表达的lncRNAs在银杏叶片转色过程中可能发挥着重要的调控作用，为后续深入研究银杏叶色变化的分子机制提供了关键的研究对象。4.3LncRNA的长度和外显子个数分布对筛选出的银杏叶片转色期lncRNAs的长度和外显子个数进行统计分析，以揭示其结构特征。结果显示，银杏叶片转色期lncRNAs的长度分布呈现出一定的范围特征。长度最短的为201nt，最长可达4800nt，平均长度为[X]nt。在长度分布上，200-500nt区间的lncRNAs数量最多，占总数的[X]%，这与多数植物中lncRNA长度集中在较短区间的情况相符。长度在501-1000nt的lncRNAs占比为[X]%；1001-2000nt的lncRNAs占[X]%；而长度超过2000nt的lncRNAs占比较少，仅为[X]%。这种长度分布模式可能与lncRNA的功能和作用机制相关，较短的lncRNA可能更易于在细胞内发挥调控作用，而较长的lncRNA可能参与更为复杂的调控网络。在研究lncRNA的外显子个数时，发现银杏叶片转色期lncRNAs的外显子个数主要集中在1-3个。其中，含有1个外显子的lncRNAs数量最多，占比达到[X]%，这类lncRNA结构相对简单，可能在基因调控中发挥较为直接的作用。含有2个外显子的lncRNAs占比为[X]%，其结构相对复杂一些，可能通过不同外显子之间的组合和相互作用来实现特定的调控功能。含有3个外显子的lncRNAs占[X]%。而外显子个数超过3个的lncRNAs占比较少，仅为[X]%。外显子个数的差异可能导致lncRNA形成不同的二级结构和空间构象，进而影响其与其他分子的相互作用和功能发挥。与其他植物相比，水稻中内含子lncRNA相对较多，其外显子个数的分布模式与银杏叶片存在明显差异，这反映了不同植物在进化过程中形成的独特基因结构和表达调控机制。4.4差异lncRNA的表达特点对筛选出的差异表达lncRNAs在银杏叶片转色不同时期的表达变化进行深入分析，以揭示其表达规律。在绿色期、黄变期和脱落期这三个关键时期，许多差异表达lncRNAs展现出独特的表达模式。例如，lncRNA-[具体编号8]在绿色期表达量相对较低，随着叶片进入黄变期，其表达量迅速上升，在黄变期达到峰值，随后在脱落期表达量又逐渐下降。这种表达模式可能与叶片在不同生长阶段的生理需求密切相关。在黄变期，叶片中的生理过程发生显著变化，如叶绿素降解、类黄酮等次生代谢产物积累等，lncRNA-[具体编号8]的高表达可能参与调控这些生理过程，推动叶片的转色进程。与之相反，lncRNA-[具体编号9]在绿色期表达量较高，随着叶片向黄变期和脱落期发展，其表达量逐渐降低。这暗示该lncRNA可能在维持叶片绿色期的正常生理功能中发挥重要作用，当叶片进入衰老转色阶段，其表达量的降低可能是叶片生理状态转变的一种调控响应。进一步分析发现，部分差异表达lncRNAs的表达变化与银杏叶片中类黄酮含量的变化呈现出显著的相关性。如lncRNA-[具体编号10]的表达量在黄变期与类黄酮含量同步上升，推测其可能参与了类黄酮合成的调控过程，促进类黄酮在叶片黄变期的积累，从而影响叶片的颜色变化。而lncRNA-[具体编号11]的表达量变化与类黄酮含量变化呈负相关，可能对类黄酮合成起到抑制作用，在叶片转色的不同阶段，通过调节类黄酮的合成水平来影响叶片的生理状态和颜色变化。通过对差异表达lncRNAs表达特点的分析，我们初步揭示了其在银杏叶片转色过程中的表达规律。这些lncRNAs的表达变化与叶片的生理状态和转色进程密切相关，为深入探究银杏叶色变化的分子机制提供了重要线索。后续研究可进一步通过实验验证这些lncRNAs的功能，如采用RNA干扰（RNAi）技术沉默特定lncRNA的表达，观察叶片转色过程及相关生理指标的变化，从而明确lncRNAs在银杏叶色变化中的具体调控作用。4.5差异lncRNA靶基因KEGG富集分析对银杏叶片转色期差异表达lncRNA的靶基因进行KEGG富集分析，结果显示，这些靶基因在多条重要的代谢通路和信号转导通路中显著富集，为深入探究lncRNA在银杏叶色变化中的调控机制提供了关键线索。在代谢通路方面，“苯丙烷生物合成通路”（Phenylpropanoidbiosynthesis）的富集程度最为突出。该通路是类黄酮合成的上游途径，类黄酮作为银杏叶片中的重要次生代谢产物，其合成过程涉及该通路中的多个关键酶和基因。在银杏叶色变化过程中，尤其是从绿色期向黄变期转变时，叶片中的类黄酮含量会发生显著变化。而差异表达lncRNA的靶基因在苯丙烷生物合成通路中显著富集，强烈提示这些lncRNAs可能通过调控苯丙烷生物合成通路中的相关基因，如PAL（苯丙氨酸解氨酶）、C4H（肉桂酸-4-羟化酶）、4CL（4-香豆酸辅酶A连接酶）等基因的表达，进而参与银杏类黄酮的合成调控，影响叶片中类黄酮的积累，最终对银杏叶色变化产生影响。例如，当某差异表达lncRNA上调表达时，其靶基因PAL的表达也随之增加，可能促进苯丙氨酸向肉桂酸的转化，进而增强苯丙烷生物合成通路的代谢流，最终使得类黄酮合成增加，叶片颜色逐渐变黄。“植物激素信号转导通路”（Planthormonesignaltransduction）也有明显富集。植物激素在植物的生长发育、衰老等过程中发挥着关键的调控作用。在银杏叶色变化过程中，内源激素的含量和平衡会发生改变。差异表达lncRNA的靶基因富集于此通路，表明lncRNAs可能通过影响植物激素信号转导，间接调控银杏叶片的衰老和转色进程。如在叶片衰老转色过程中，乙烯作为一种促进衰老的植物激素，其含量会逐渐升高。若某差异表达lncRNA能够调控乙烯信号转导通路中的关键基因，如乙烯受体基因或乙烯响应因子基因的表达，则可能影响乙烯信号的传递和响应，进而影响叶片的衰老和转色。当lncRNA抑制乙烯受体基因的表达时，可能减弱叶片对乙烯的敏感性，延缓叶片的衰老和转色进程。“光合作用相关通路”（Photosynthesis-relatedpathways）也出现了一定程度的富集。在银杏叶色变化过程中，叶片的光合作用能力会逐渐下降。差异表达lncRNA的靶基因在光合作用相关通路中富集，暗示lncRNAs可能通过调控光合作用相关基因的表达，影响叶片的光合作用效率，从而参与叶色变化的调控。例如，lncRNA可能调控光合色素合成相关基因，如叶绿素合成酶基因的表达，当lncRNA下调表达时，可能导致叶绿素合成酶基因表达降低，叶绿素合成受阻，叶片中的叶绿素含量减少，从而使叶片逐渐失去绿色，呈现出黄色或其他颜色。此外，“淀粉和蔗糖代谢通路”（Starchandsucrosemetabolism）也有部分靶基因富集。淀粉和蔗糖是植物体内重要的碳水化合物，其代谢过程与植物的能量供应和物质积累密切相关。在银杏叶色变化过程中，叶片中的碳水化合物代谢也会发生改变。差异表达lncRNA的靶基因在该通路中富集，表明lncRNAs可能通过影响淀粉和蔗糖代谢，为叶片的转色过程提供能量和物质基础，进而参与叶色变化的调控。如在叶片衰老转色时，淀粉可能会分解为蔗糖等小分子糖类，为叶片的生理活动提供能量。若某差异表达lncRNA能够调控淀粉合成酶或蔗糖合成酶等关键酶基因的表达，则可能影响淀粉和蔗糖的代谢平衡，从而影响叶片的生理状态和颜色变化。当lncRNA上调表达促进蔗糖合成酶基因表达时，可能会加速淀粉的分解和蔗糖的合成，为叶片的转色过程提供更多的能量和物质支持。4.6银杏lncRNA的靶基因预测、GO功能注释采用cis-acting和trans-acting两种策略对银杏叶片转色期差异表达lncRNA的靶基因进行预测。cis-acting预测中，将距离lncRNA上下游100kb范围内的基因作为潜在靶基因，共预测得到[X]个cis作用的靶基因。trans-acting预测通过计算lncRNA与mRNA之间的序列互补性，利用RNAplex软