红鳍东方鲀结合CCL25b核心启动子区转录因子的鉴定

Identification of Transcription Factors Binding the Core Promoter Region of CCL25b in Takifugu rubripes

摘要: 基因转录调控一直是生物信息学研究的一个重要内容，转录因子与启动子区域上的特异性调控元件的结合在基因表达和调控中具有重要意义，是构建基因调控网络的一个核心问题。趋化因子配体25 (chemokine 25, CCL25)是一种属于CC趋化因子家族的小细胞因子，在免疫细胞调节和炎症过程如T细胞归巢和慢性组织炎症中均发挥重要作用。为了探究CCL25在表达过程中受到的转录调控机制，本实验我们初次采用聚合酶链式反应(PCR)扩增出红鳍东方鲀CCL25b的核心启动子区，并利用DNA pull down方法对结合在该序列的蛋白质进行富集，最后结合质谱鉴定数据及GO富集分析初步筛选出6个转录因子(C1QBP、PURA、ARHGAP35、NME2、RAP2C和DRG1)，这些初步结果将有助于推进对与CCL25表达调节机制有关的潜在生物学过程的理解。

Abstract: Gene transcriptional regulation has been an important element of bioinformatics research. The binding of transcription factors to specific regulatory elements on promoter regions is important in gene expression and regulation, and is a central issue in the construction of gene regulatory networks. Chemokine ligand 25 (CCL25) is a small cytokine belonging to the CC chemokine family that plays an important role in immune cell regulation and inflammatory processes such as T cell homing and chronic tissue inflammation. In order to investigate the transcriptional regulation mechanism of CCL25 in its expression, Polymerase chain reaction (PCR) was used to amplify the core promoter region of CCL25b in Takifugu rubripes, and then DNA pull down was used to enrich the proteins bound to this sequence. Finally, we combined the mass spectrometry identification data with GO enrichment analysis to initially identify six transcription factors (C1QBP, PURA, ARHGAP35, NME2, RAP2C and DRG1), and these results will help advance the understanding of the underlying biological processes associated with the regulatory mechanism of CCL25 expression.

1. 引言

趋化因子配体25 (chemokine 25, CCL25)是一种属于CC趋化因子家族的小细胞因子，也称为胸腺表达的趋化因子(TECK)。研究表明CCL25在T细胞的发育中起作用 [1]。在豚鼠骨关节炎模型中也被证实，关节内注射CCL25联合透明质酸可保护软骨免受蛋白多糖丢失和与膝关节软骨破坏相关的形态恶化 [2]。多项研究表明，趋化因子受体9 (chemokine receptor 9, CCR9)及其专用配体CCL25在多种恶性肿瘤中过表达，并与肿瘤的增殖、凋亡、侵袭、迁移和耐药密切相关 [3]。综上，CCL25在免疫细胞调节和炎症过程如T细胞归巢和慢性组织炎症中发挥重要作用。因此，全面了解CCL25表达及其活性的调控机制具有重要意义。

基因的转录是由被称为启动子的上游调控DNA的结构决定的。理解真核生物和原核生物的转录调控，需要对与感兴趣的启动子结合的调控蛋白进行鉴定。转录因子与启动子区域上的特异性调控元件的结合在基因表达和调控中至关重要 [4]。激活转录这一活动需要与DNA结合的反式激活子募集转录共激活因子或抑制子，以及众多与基础转录因子相互作用的相关蛋白，因而分析这些与DNA结合的大型蛋白质复合物是阐明调节基因表达机制的重要步骤 [5]。转录因子一般由两个功能结构域组成，一是DNA结合结构域，二是转录活化结构域。根据作用效果的差异，转录因子可分为转录抑制因子和转录激活因子两大类。当细胞中没有诱导物时，转录抑制因子可以与调控序列结合并使附近的启动子失去转录活性。一旦诱导物出现，抑制因子即从靶DNA区域脱离，同时启动子的转录活性恢复。而转录激活因子一般需要与效应配体相互作用才能实现其激活功能。根据启动子DNA结合位点的不同，抑制因子和激活因子之间是可以转化的 [6]。

目前，已有多种方法可用于DNA结合蛋白的检测和分离，比如：凝胶迁移实验(EMSA)、足迹实验(Foot-printing Assay)、干扰实验(Interference Assays)、染色质免疫共沉淀(ChIP)、DNA-蛋白质印迹杂交(Southwestern Blot)和DNA pull down等。而相较于以上所有的分离方法，DNA pull down的优点在于可以同时分析复合物中的大量蛋白质。Wu-Guo Deng等人将环氧合酶-2启动子作为探针，与核提取物孵育，最后通过免疫印迹评估复合物中的蛋白质，检测到一系列转录因子 [7]。为了探究蛋白混合物中的结合蛋白，质谱技术可以对蛋白质复合物进行更加深入地定量分析。Diem Hong Tran等人将DNA pull down与二维凝胶电泳和质谱分析相结合，鉴定了六种与人D-氨基酸氧化酶基因5’侧翼区域结合的蛋白质 [8]。该方法在确定反激活因子p300/CBP和相关蛋白与环氧化酶-2 (COX-2)启动子结合的调控方面同样被证明是有用的。

本实验中，我们使用PCR技术扩增获得了CCL25b的启动子序列，并成功制备了红鳍东方鲀脑和鳃组织总蛋白。通过DNA pull down技术获得了与CCL25b基因的启动子序列相互作用的蛋白，这些蛋白经聚丙烯酰胺凝胶电泳检测并进行质谱鉴定分析。我们的结果首次报道了CCL25b基因启动子序列相互作用蛋白，这可能成为进一步研究CCL25b基因表达调控因子的潜在靶点之一。

2. 材料与方法

2.1. 引物设计

本实验以红鳍东方鲀基因组DNA为模板，在NCBI数据库中检索CCL25b基因序列，确定出CCL25b的核心启动子区域，选取基因转录起始位点上游2000 bp及下游500 bp为靶序列进行巢式PCR反应以克隆启动子片段，利用Primer Primer5.0设计两对引物(表1)，其中CCL25bF1和CCL25bR1为CCL25b巢式PCR的外引物，CCL25bF2b和CCL25bR2U为CCL25b巢式PCR的内引物，UPL为具有生物素标记的特异性通用引物。CCL25bF1和CCL25bR1引物委托上海捷瑞生物工程有限公司合成并纯化，CCL25bF2b、CCL25bR2U和生物素标记的特异性通用引物委托华大基因合成并纯化。不同引物在基因CCL25b核心启动子区的定位如图1所示，其中TSS为转录起始位点。

Table 1. Primer sequences

表1. 引物序列

Figure 1. Cloning strategy of promoters

图1. 启动子的克隆策略

2.2. 基因组DNA的提取及启动子片段的标记

取适量红鳍东方鲀的肌肉组织，使用天根生化科技(北京)有限公司的DNA提取试剂盒(目录号：DP324)提取基因组DNA，提取的DNA置于−20℃冰箱保存备用。

以红鳍东方鲀基因组DNA为模板，使用外引物对进行普通PCR反应以克隆启动子片段，反应体系(20 μL)如表2。

Table 2. PCR reaction system

表2. PCR反应体系

反应条件：94℃预变性5 min；94℃变性30 s，59℃退火30 s，72℃延伸2 min 10 s，共进行34个循环；72℃延伸10 min，4℃保存。

用1%琼脂糖凝胶电泳检测PCR产物条带大小是否正确，PCR产物条带检测合格后，使用天根生化科技(北京)有限公司的DNA凝胶回收试剂盒(目录号：DP209)回收纯化正确条带，实验步骤参照试剂盒所提供方法。

将回收到的基因片段作为二次PCR的模板，进行降落PCR反应以克隆带有生物素标记的启动子片段，反应体系(20 μL)如表3。

Table 3. Touchdown PCR reaction system

表3. 降落PCR反应体系

反应条件：94℃预变性5 min；94℃变性30 s，72℃退火与延伸1 min 30 s，共进行5个循环；94℃变性30 s，68℃退火30 s，72℃延伸1 min 10 s，共进行5个循环；94℃变性30 s，60℃退火30 s，72℃延伸1 min 10 s，共进行25个循环；72℃延伸10 min，4℃保存。

用2%琼脂糖凝胶电泳检测PCR产物条带大小是否正确，在检测产物条带大小正常的情况下，收集带有生物素标记的启动子片段并用天根生化科技(北京)有限公司的普通DNA产物纯化试剂盒(目录号：DP204)进行纯化，实验步骤参照试剂盒所提供方法。将纯化后带有生物素标记的启动子片段送至北京六合华大基因科技有限公司进行测序，使用NanoDrop2000 (Thermo Scientific, USA)进行浓度测定并置于−20℃保存以备下游实验。

2.3. 活性组织总蛋白的制备

使用上海生工的《一步法动物组织活性蛋白提取试剂盒》(产品编号：C500006)提取红鳍东方鲀脑、鳃组织总蛋白。将提取得到的组织活性总蛋白用biosharp的BCA蛋白浓度测定试剂盒(产品编号：BL521A)及紫外分光光度计进行蛋白标准曲线的绘制并计算所提蛋白浓度，实验步骤参照试剂盒所提供方法。最后将所提活性蛋白分装并置于−80℃保存备用。

2.4. 启动子–磁珠的结合

按照Streptavidin Magnetic Beads (New England BioLabs, USA. 目录号：S1420S)说明书中的磁珠与核酸的结合容量，每1 mg磁珠可结合500 pmol以上的单链25 bp生物素标记的寡核苷酸，结合2.2中测得的带有生物素标记的启动子片段浓度，将磁珠与生物素标记的启动子片段共孵育。首先取出于4℃存放的磁珠并涡旋振荡使磁珠彻底悬浮，然后取30 μL磁珠置于1.5 mL无菌无酶的离心管中，向其中加入200 μL Wash/Binding Buffer，并涡旋以悬浮磁珠，置于磁力架上静置2 min，去上清，重复三次；接着向洗涤后的磁珠中加入520 μL生物素标记的CCL25b启动子片段(按照磁珠结合核苷酸标准的1.5倍原则加入)，涡旋以悬浮磁珠，4℃过夜，然后将磁珠与启动子片段的复合物置于磁力架上静置2 min，去上清；向所得沉淀中加入500 μL Wash/Binding Buffer，涡旋以悬浮磁珠，置于磁力架上静置2 min，去上清，重复三次；向上述所得沉淀中加入100 μL TE并用移液器轻轻吹吸混匀，得启动子–磁珠复合物。

2.5. DNA pull down

根据2.3中所测得的红鳍东方鲀鳃、脑组织总蛋白的浓度，按照每30 μL磁珠加入500 μg组织总蛋白的标准，制备启动子–磁珠–蛋白复合物。首先将活性组织总蛋白从−80℃取出并置于冰盒内缓慢复温，取110 μL脑组织蛋白样品、170 μL鳃组织蛋白样品分别向其中加入2 μL蛋白酶抑制剂、2 μL DTT，混匀后加入到过夜孵育后的启动子–磁珠复合物中并置于4℃孵育5 h。然后将孵育后的复合物置于磁力架上静置2 min，去除上清液，收集磁珠。用冰冷PBS (500 μL)洗涤上述收集到的磁珠，重复三次，尽可能去除上清，收集沉淀，得到启动子–磁珠–蛋白复合物。向上述所得复合物中加72 μL PBS悬浮，取出15 μL用于SDS-PAGE垂直电泳检测，剩余57 μL置于磁力架上加磁力2分钟，弃上清，沉淀置于−80℃保存用于质谱鉴定。

2.6. 银染

SDS-PAGE垂直电泳完毕后，我们使用上海生工的低背景银染试剂盒(产品编号：C510027)对凝胶进行染色，实验步骤参照试剂盒所提供方法，染色完毕后进行拍照。

2.7. 质谱检测

为了分析DNA pull down得到的CCL25b启动子结合蛋白的种类，我们将步骤2.5中得到的启动子–磁珠–蛋白复合物送到上海拜谱生物科技有限公司进行质谱分析(QE鉴定)。质谱测试数据用MaxQuant1.6.1.0软件检索Uniprot Protein Database数据库。基于强度定量法(Intensity based absolute quantification, iBAQ)对鉴定的蛋白质进行定量 [9]，其基本原理为每个蛋白质鉴定到的全部肽段曲线下的面积之和(AUCs)与此蛋白的理论酶切肽段数(N)之比为此蛋白的表达量(AUCs/N)。

2.8. GO富集分析

将质谱结果分析后，通过R语言包clusterProfiler v3.18.1 [10] 对鉴定得到的蛋白进行GO富集分析，以此确定质谱检测到的蛋白质行使的主要生物学功能。其中adjusted P-values < 0.05的GO术语即为质谱检测到的蛋白质显著富集的GO term。

2.9. 蛋白质相互作用分析

结合GO富集分析结果以及UniProt (https://www.uniprot.org/)和GeneCards (https://www.genecards.org/)网站对蛋白基本功能的注释，我们将富集到“核酸结合”和“核苷结合”术语上的关键蛋白导入到STRING v11.0数据库(https://string-db.org/)中，筛选可信度(combined score ≥ 0.15)较高的蛋白质互作网络(Protein-protein interaction network, PPI)，之后将数据导入到Cytoscape v3.7.1 (https://cytoscape.org/)进行可视化输出。

3. 结果与分析

3.1. 生物素标记启动子片段的制备

使用外引物对通过PCR扩增技术、琼脂糖凝胶电泳、切胶纯化并回收，获得一次PCR产物(图2(A))。为保证所回收启动子片段的纯度及特异性，再以一次PCR产物为模板，使用内引物对进行降落PCR扩增技术、琼脂糖凝胶电泳、普通DNA产物纯化并回收，即获得目标启动子片段(图2(B))，图2表明目标启动子片段纯度较高，无杂带污染。测序结果同样显示启动子片段能够很好的比对到红鳍东方鲀基因组上。

(A) M. DL2000 DNA Marker; M’. λ DNA-HindIII Digest DNA Marker; 1, 2. CCL25b启动子片段(2013 bp)；(B) M. DL2000 DNA Marker; 1, 4. 生物素标记的CCL25b启动子片段(1145 bp); 2, 5. 未被生物素标记的CCL25b启动子片段(1120 bp); 3, 6. 阴性对照(无条带)

Figure 2. Preparation of biotin-labeled promoter fragments

图2. 生物素标记的启动子片段的制备

3.2. DNA pull down

根据步骤2.5中的实验操作，将制备的启动子–磁珠复合物分别与脑、鳃组织总蛋白孵育，纯化后将特异结合的蛋白进行SDS-PAGE垂直电泳及银染，银染结果如图3所示。

M. 蛋白质分子量标准；1. CCL25b启动子在脑中钓取的蛋白(转录因子)；2. CCL25b启动子在鳃中钓取的蛋白(转录因子)；3. 生物素标记的CCL25b启动子片段；M’. DL2000 DNA Marker

Figure 3. Silver staining results of DNA pull down products

图3. DNA pull down产物SDS-PAGE后银染

3.3. 质谱结果分析与转录因子筛选

完成样品胰酶酶解等前处理、二维色谱上机分析以及单次数据检索分析，质谱结果显示在脑中共鉴定得到1705个蛋白质，在鳃中共鉴定得到2458个蛋白质。将这些鉴定得到的蛋白通过R语言包clusterProfiler v3.18.1进行GO富集分析，GO富集分析结果显示有56个编码生物学功能与核酸结合相关的差异性蛋白质的基因ID，有144个编码生物学功能与核苷结合相关的差异性蛋白质的基因ID (图4)。在UniProt和GeneCards网站依据蛋白的基本功能分别对上述200个差异性蛋白质的基因ID进行检索，结果显示富集到“DNA结合”这一术语上的蛋白包括RNA聚合酶、转录调控蛋白和其它结合蛋白，我们从中筛选出6个生物学功能与转录过程相关的蛋白质，其在脑和鳃中的iBAQ值如表4。

3.4. 蛋白质相互作用分析

经STRING数据库及Cytoscape软件构建的可视化PPI网络如图所示(图5)。整个网络包含两种类型的互作：DNA pull down实验验证和来自文献的互作关系。主网络中共包含22种蛋白质，经实验验证得到的潜在转录因子有6个，包括PURA、ARHGAP35、NME2、RAP2C、C1QBP和DRG1。其中ARHGAP35 (9)的度最高。与上述转录因子以及靶基因CCL25b存在直接或间接相互作用的蛋白有15个。其中CDH1 (13)、RASA1 (10)和RHOA (10)的度相对较高。

Figure 4. GO enrichment analysis

图4. GO富集分析

Table 4. Transcriptional regulators obtained by mass spectrometry identification (6)

表4. 质谱鉴定获得的转录调控因子(6个)

注：不同的颜色代表不同类型的互作关系。每个节点都是一个蛋白。节点的大小代表节点的度值。红色表示靶基因CCL25b，绿色表示经实验验证得到的潜在转录因子，蓝色表示与靶基因以及潜在转录因子相互作用的蛋白。

Figure 5. Protein interaction network diagram

图5. 蛋白质相互作用网络图

4. 讨论

遗传信息传递的过程较为复杂，主要涉及染色质激活、转录起始、转录后修饰、转录产物在细胞内的转运、翻译起始、翻译后修饰等诸多步骤，而针对基因表达的调控可发生于上述的任一环节。转录作为基因组表达的第一步，是真核细胞所有生命活动的基础，参与细胞的生长和分裂、基因组结构和修复，与炎症和癌症等疾病发生发展也有密切关系，因而发生在转录水平的调节，尤其是转录起始水平的调节，对基因表达有着至关重要的作用 [4]。真核生物转录起始需要RNA聚合酶对起始点上游DNA序列进行识别和结合，生成转录前起始复合体，而这种识别和结合依靠于转录因子。转录因子是一种可以调控基因表达功能的蛋白质分子，也称为反式作用因子。

为了研究红鳍东方鲀CCL25b基因在细胞中的表达调控机制，本研究中，我们首次通过DNA pull down技术对结合在红鳍东方鲀CCL25b启动子区域的转录因子进行了富集。将组织活性总蛋白、5’末端被生物素标记的双链DNA启动子片段和带有链霉亲和素标记的磁珠共同孵育后，由于生物素与链霉亲和素的极高亲和力，复合物可被下拉 [11]。随后我们将下拉产物进行银染，观察到有特异性条带之后将下拉产物送质谱检测。将质谱鉴定得到的蛋白进行GO富集分析，富集分析结果表明有56个编码蛋白的基因富集到“translation regulator activity, nucleic acid binding”这一术语上，有144个编码蛋白的基因富集到“nucleoside binding”这一术语上。此外，富集结果分析表明许多经质谱鉴定到的基因还富集到了与翻译、酶结合、核糖体形成等相关的GO术语上，我们猜测造成这种现象的原因可能有以下两点：其一可能是由于在整个pull down过程中特别是在提取活性组织总蛋白时没有进行相应的RNaseA的处理，导致大量的mRNA会与启动子片段下游(能转录产生mRNA的部分)产生杂交，继而将与mRNA相结合的核糖体钓取出来。其二可能是由于本次实验中我们采用的蛋白是组织总蛋白，胞质中的蛋白可能会对实验结果产生一定的影响。

本实验中我们基于iBAQ对鉴定的蛋白质进行定量。根据蛋白质iBAQ值的高低以及在UniProt和NCBI等数据库中进行的文献调研和生物信息学分析结果，本文选择了DNA结合蛋白C1QBP、PURA、ARHGAP35、NME2、RAP2C、DRG1进行后续的分析。

研究表明基因PURA在人类中编码Pur-alpha蛋白 [12]，Pur-alpha是高等生物中发现的第一个序列特异性单链DNA结合蛋白 [12]。它与单链和双链DNA都有结合，在其结合部位的富含嘌呤的链上与G残基进行接触。Pur-alpha属于由四个成员组成的Pur蛋白家族，其中还包括Pur-beta (GenBank AY039216.1; GI:14906267) [12] 和两种形式Pur-gamma (变体A，GenBank AF195513.2；变体B，GenBank AY077841) [13]。从细菌到人类的Pur蛋白质序列都包含一个高度保守的氨基酸片段 [12] [14] (参见NCBI smart00712)。人类的Pur-alpha包含这个Pur结构域的三个重复，细菌的Pur-alpha只包含一个 [12] [15]。这种进化上的保守性意味着这个结构域的特定序列对整个生物体光谱中的大多数物种的生存都很重要。由于Pur-alpha在低等生物和人类中的功能差别很大，因此Pur结构域的这种基本性质引起了人们的关注。研究表明Pur-alpha对哺乳动物的大脑和血细胞发育至关重要 [16]，在我们的研究中也发现转录因子PURA在红鳍东方鲀脑中的iBAQ值要远远高于鳃中的，这与之前的研究结果保持一致。在人类中，Pur-alpha的功能是激活细胞核中的转录，促进细胞质中的RNA运输，并在细胞周期中调节DNA复制 [14]。在某些功能中，Pur-alpha与家族成员Pur-beta相互作用 [17] [18]。一些细胞周期调控功能可能是通过Pur-alpha与Cyclin/Cdk蛋白激酶的结合来实现的，Cyclin/Cdk蛋白激酶对调节细胞周期转换点的蛋白质进行磷酸化 [19] [20]。所有生物体对Pur-alpha的要求都是通过Pur-alpha与核酸结合的能力以及与调节和运输蛋白的相互作用来实现的。

ARHGAP35，也称为GRLF1，其在人体中编码糖皮质激素受体的DNA结合因子1。研究表明人糖皮质激素受体DNA结合因子与糖皮质激素受体基因(hGR基因)的启动子区相关联，是糖皮质激素受体转录的阻遏物。此外，生物化学分析表明，ARHGAP35的相互作用具有序列特异性，ARHGAP35的转录功效是通过其与特定序列基序的相互作用来调节的，其表达水平受糖皮质激素调节。研究发现糖皮质激素是大脑正常发育所必需的，它启动终端成熟，重塑轴突和树突，并影响细胞的生存 [21]，也可能在海马发育中发挥作用。这与我们的研究结果保持一致，我们的质谱鉴定数据表明在红鳍东方鲀的脑中鉴定到了ARHGAP35，而在红鳍东方鲀的鳃中没有发现ARHGAP35。此外，糖皮质激素通过与糖皮质激素受体结合影响细胞，被激活的糖皮质激素受体–糖皮质激素复合物在细胞核中上调抗炎蛋白的表达(这一过程被称为转录激活)，并通过阻止其他转录因子从细胞核中转录到细胞内，抑制细胞内促炎蛋白的表达(转录抑制) [22]。糖皮质激素也被证明在T淋巴细胞的发育和稳态中起作用，这在T细胞系对糖皮质激素的敏感性增加或降低的转基因小鼠中得到了证明 [23]。有趣的是，目标基因CCL25也被认为在T 细胞的发育中发挥作用 [1]。鉴于功能上的相关性，我们推测ARHGAP35很有可能作为CCL25的转录因子而发挥作用。

研究表明NME2作为MYC基因的转录激活剂，能够非特异性结合DNA [24] [25]。与MYC基因启动子的核酸酶超敏元件(NHE) III (1)区域内的单链鸟嘌呤和富胞苷链结合而不与双链NHE III (1)结合 [25]。而Myc蛋白作为转录因子，又可以通过结合增强子盒序列(E-boxes)和招募组蛋白乙酰转移酶(HATs)来激活许多促增殖基因的表达。研究表明Myc的功能是通过招募转录延伸因子来调控活跃转录基因的转录延伸 [26]。它也可以作为转录抑制因子，通过结合Miz-1转录因子并置换p300共同激活因子，抑制Miz-1目标基因的表达。此外，myc在控制DNA复制方面有直接作用 [27]，这种活动可能有助于癌细胞的DNA扩增。

RAP2C基因编码的蛋白质是Ras GTPase超家族的Ras相关蛋白亚家族的成员。该家族的成员是小的GTP酶，它们充当分子开关来调节细胞增殖、分化和凋亡。据报道，该蛋白质可激活血清反应元件的体外转录活性。可能在SRE介导的基因转录中起作用 [28]。研究表明基因C1QBP可能参与调节FOXC1转录活性和NFY/CCAAT结合因子复合物介导的转录 [29] [30]。有趣的是，我们的质谱数据显示RAP2C蛋白和C1QBP蛋白只在红鳍东方鲀的的脑组织中鉴定到，这可能与转录因子的组织特异性相关。

Drg1是一种人类GTPase，属于Obg家族，在古生菌和真核生物中发现，但在真细菌中未发现 [31] [32]。DRG1通过水解GTP中的γ-磷酸键，催化GTP向GDP的转化，具有内在的GTP酶活性 [33] [34]。我们猜测这可能与转录因子的活性相关。

尽管转录领域的研究在过去的半个世纪里取得了重大的进展，但理解转录的机制和调控对生命科学来说仍然是一个巨大的挑战，具有核心的重要性。本实验的质谱鉴定获得了多个候选的DNA结合蛋白，其中不乏非特异结合蛋白。虽然这些蛋白在调控CCL25b基因表达中的实际功能尚不清楚，但我们为寻找CCL25b的DNA相互作用蛋白提供了新的思路。这些初步结果将有助于我们进一步了解与CCL25b表达调控机制相关的潜在生物学过程。

基金项目

国家重点研发计划资助(编号：2018YFD0900301-10)。

NOTES

^*通讯作者。