Change and correlation analysis of pigment contents and color value during growth of Chlorella vulgaris
-
摘要:
该文基于CIE-L*a*b*颜色分析系统,采用全自动分光测色仪对热带普通小球藻 (Chlorella vulgaris) 生长过程中的颜色值变化进行监测,同时测定其细胞密度及色素的质量浓度变化并进行相关分析。结果显示,热带普通小球藻生长过程中细胞密度分为接种期、快速生长期、稳定期和衰亡期4个阶段,绿色素和黄色素质量浓度和颜色参数值黄度值(b*)、色差值(△E)和饱和度(C)均呈增加趋势,而颜色参数值明度值(L*)、红度值(a*)、相角(H°)呈降低趋势;颜色参数值与细胞密度、色素质量浓度具有较强的相关性,线性回归模型拟合度较高,在快速生长期尤为显著,其中b*和△E可作为较佳的指示指标。通过量化热带普通小球藻生长过程中的颜色参数变化,建立非接触性数量及质量动态变化预测模型,可快速准确预测,为提高养殖过程中精细操作及效率提供理论基础。
Abstract:Automatic color spectrometer was used based on CIE-L*a*b* color system to monitor the color value change in the growth process of Chlorella vulgaris. The cell density and pigment contents, as well as their correlational relationship were analyzed. The results show that C. vulgaris growth was divided into four phases (seeding, rapidly growing, stable and declining). The general trends of green and yellow pigments were increasing, as well as yellow value (b*), color difference (ΔE) and chroma (C), but those of lightness value (L*), red value (a*) and hue angle (H°) were decreasing. Strong correlation had been observed between color values and cell density, pigment content. The linear regression models had a high fit degree especially in rapid growing phase. The b* and ΔE values could be used as good indicators. The establishment of non-contact quantity and quality dynamic prediction model can make the prediction quick and accurate, and provide a theoretical basis for improving the precision operation and efficiency in the process of aquaculture.
-
Keywords:
- Chlorella vulgaris /
- Growth process /
- Color parameters /
- Cell density /
- Pigment content /
- Correlation analysis
-
三角帆蚌 (Hyriopsis cumingii) 所产珍珠占中国珍珠产量的80%以上[1],且品质良好、经济价值较高。但三角帆蚌培育珍珠的大小、光泽度等均与性别有关,在产珠性能方面雄性优于雌性[2],而人工雌雄分养情况下能促进雌性三角帆蚌的个体发育及产珠性能[3]。三角帆蚌属于雌雄异体,不存在异型染色体,目前性别决定无法从染色体层面进行研究,所以性别分化与决定相关基因的研究成为性别调控的热点,对揭示复杂的性别决定与分化机制研究及创造经济效益具有重大意义[4]。
丝裂原活化蛋白激酶 (MAPK) 也称为细胞外信号调节激酶,受各种细胞外生长因子及受体的相互作用而被激活,会影响许多组织特异性生物活性,如细胞增殖、存活和分化[5-6]。然而,MAPK在雌性生殖细胞中不遵循这种传统的模式[7]。在卵母细胞中,MAPK独立于生长因子和酪氨酸激酶受体而被激活,独立于转录调控发挥作用[8],受特殊上游调节剂激酶 MOS的控制,在卵母细胞减数分裂中起关键作用[9-10]。Mos/MAPK通路在脊椎动物中启动卵母细胞的成熟发育[10],C-MOS基因在Mos/MAPK通路中位于MAP2K1基因的上游。MAP2K1基因作为MAPK家族的一分子,参与细胞增殖、分化、基因转录调控和发育过程[11-13]。在脊椎动物中,MAP2K1基因参与小鼠 (Mus musculus) 的卵巢发育[14],促进水牛 (Bubalus) 卵母细胞的表达与发育[15],对鹅 (Anser cygnoides orientalis) 的卵泡发育起到支持与促进作用[5]。在无脊椎动物中,MAP2K1基因对斑节对虾 (Penaeus monodon) 卵母细胞的生长发育起促进作用[14];能促进果蝇 (Drosophilid) 的卵母细胞成熟[16]。
在贝类中,紫贻贝 (Mytilus galloprovincialis)、虾夷扇贝 (Mizuhopecten yessoensis) 的MAP2K1基因均有报道[17-18]。本文对三角帆蚌MAP2K1基因进行了克隆,分析了该基因在性腺不同阶段和雌雄性腺组织中的表达模式,探讨了该基因在三角帆蚌性腺发育中的作用。
1. 材料与方法
1.1 实验材料
从浙江省金华市武义实验基地选取健康的1—8月龄和1、2、3龄三角帆蚌,暂养1周后对1—8月龄个体的性腺组织、1龄和2龄个体各组织 (性腺、闭壳肌、肝胰腺、鳃、外套膜、斧足)、3龄个体的性腺组织进行取样,样本立刻放入液氮,速冻后用镊子放入样本盒,−80 ℃保存。
1.2 总RNA提取与cDNA合成
按TRIzol法提取总RNA,利用NanoDrop 2000分度计进行核酸检测,通过1%琼脂糖凝胶电泳[19] (电压180 V,电流200 mA) 检查RNA完整性情况 (3条条带5.8S、18S、28S)。检查合格后将RNA保存于−80 ℃。将RNA进行反转录 (Prime ScriptTM RT reagent Kit with gDNA Eraser试剂盒,TaKaRa, 日本) 得到cDNA模板,cDNA保存于−40 ℃。
1.3 基因序列的获取
通过转录组库[20]中得到的序列设计3'RACE的inner和outer引物 (表1),并参照SMART 3'RACE试剂盒 (Clontech, America) 说明书进行MAP2K1基因的3'端PCR扩增。电泳 (电压180 V,电流200 mA) 检测产物,割胶回收纯化。用PMD19-T (TaKaRa,大连) 载体连接纯化产物,16 ℃下连接4 h。准备冰盒,在冰上将连接产物转化入大肠杆菌DH5α (TaKaRa,大连) 中 30 min。准备含氨苄的固体培养基 (4 ℃冰箱保存),将所得菌液均匀地涂在上面,37 ℃培养10 h。挑选白色菌株送至生工 (上海) 测序。
表 1 本实验所需引物Table 1. Primers for this study引物名
Primer's name引物序列 (5'—3')
Primer sequence (5'−3')用途
PurposeK1-O TATTCCTACAGCAGACTCGTGATTG 3'RACE outer K1-I GACCTGGATTTGATAATAGATGTGGG 3'RACE inner K1-YG-F AACAAGCTGAATCTGACGCTG 荧光定量正向引物 K1-YG-R GCTCCTTCAGTTTCTTGGTCAGT 荧光定量反向引物 EFl-αF GGAACTTCCCAGGCAGACTGTGC 内参正向引物 EFl-αR TCAAAACGGGCCGCAGAGAAT 内参反向引物 K1-YW-F AACAAGCTGAATCTGACGCTG 原位杂交正向引物 K1-YW-R TAATACGACTCACTATAGGGGCTCCTTCAGTTTCTTGGTCAGT 原位杂交反向引物 M-GR-1F TTGGGTATGGTTCGCAGTT 干扰链1正向引物 M-GR-1R GCCCTTGTAAGCTCATGTTT 干扰链1正向引物 T7+M-GR-1F TAATACGACTCACTATAGGGTTGGGTATGGTTCGCAGTT T7+干扰链1正向引物 T7+M-GR-1R TAATACGACTCACTATAGGGGCCCTTGTAAGCTCATGTTT T7+干扰链1反向引物 M-GR-2F TATGAACGAACATCGGAGGC 干扰链2正向引物 M-GR-2R CCCATACGCAACAACCC 干扰链2反向引物 T7+M-GR-2F TAATACGACTCACTATAGGGTATGAACGAACATCGGAGGC T7+干扰链2正向引物 T7+M-GR-2R TAATACGACTCACTATAGGGCCCATACGCAACAACCC T7+干扰链2反向引物 M-GR-3F TAACGCCTCACGACGAC 干扰链3正向引物 M-GR-3R CCGGCTCTTGCATTCC 干扰链3反向引物 T7+M-GR-3F TAATACGACTCACTATAGGGTAACGCCTCACGACGAC T7+干扰链3正向引物 T7+M-GR-3R TAATACGACTCACTATAGGGCCGGCTCTTGCATTCC T7+干扰链3反向引物 1.4 基因生物信息学分析
使用NCBI (https://www.ncbi.nlm.nih.gov/) 中的ORF finder预测MAP2K1基因的开放阅读框 (ORF)、Blast得到不同物种MAP2K1基因序列、Premier设计引物;蛋白基本理化性质利用ProtParam tool (https://web.expasy.org/aprotparam/) 分析;二级结构利用SOPMA (https://npsa-prabi.ibcp.fr/NPSA/npsa_sopma.html) 预测;三级结构利用I-TASSER (https://zhanglab.ccmb.med.umich.edu/I-TASSER/) 预测;跨膜结构域利用TMHMM Serverv 2.0程序 (http://www.cbs.dtu.dk/services/TMHMM/) 预测并用 (http://www.genedenovo.com/news/582.html) 解读结果;信号肽利用SignalP-5.0 (http://www.cbs.dtu.dk/services/SignalP/) 预测并用 (http://www.cbs.dtu.dk/services/SignalP-3.0/output.php) 解读结果,基因表达量差别利用SPSS 18.0软件分析;表达量差别条形图利用SigmaPlot 12.5绘制;NJ (邻接法) 系统进化树利用MEGA 7.0构建,Bootstrap分析其可靠性,重复1 000次。
1.5 荧光定量分析
选取EF-1α为内参基因,使用CFX96 TouchTM Real-Time PCR Detection System (Bio-Rad) 仪器进行定量分析。定量体系 (20 μL) 为cDNA (使用前混合均匀) 1.6 μL,2×TB Green Premix ExTaq 10 μL,F (上游引物) 0.8 μL,R (下游引物) 0.8 μL,RNase Free Water 6.8 μL。定量程序为95 ℃ 3 min;95 ℃ 5 s,60 ℃ 30 s,循环40次。
1.6 原位杂交
原位杂交的引物具有特异性,上游引物保持不变,将T7启动子序列 (TAATACGACTCACTATAGGG) 加入到下游引物5'端前。三角帆蚌性腺cDNA作为模板进行扩增,得到目的条带割胶回收,纯化后作为模板,使用T7 (T7 High Efficiency Transcription) 试剂盒和DIG RNA Labeling Mix得到标记探针,cDNA 1 µL,4 μL 5×T7 Transcription Reaction、8 μL 10 mmol·L−1 NTP Mix,2 μL T7 Transcription Enzyme Mix,RNase-free water 5 μL。反应程序为37 ℃,2 h。将所得探针纯化后保存于−80 ℃。
取2龄健康雌、雄三角帆蚌的性腺组织 (组织上连带少许表皮),将其在4%多聚甲醛与DEPC中固定2 h,转移至70%乙醇中,在乙醇中可保存约10 d。进而进行组织梯度脱水、组织透明、浸蜡、包埋、切片 (厚度控制在6~8 μm)、粘片。利用DIG nucleic acid detection kit进行原位杂交,并进行封片处理,在显微镜下观察杂交信号并拍照。
1.7 RNA干扰
1.7.1 dsRNA的合成
使用Primer Premier 5.0软件在C-MOS基因的ORF区设计3对引物,在上下游引物前加上T7序列,将每组中正常引物与加有T7序列的引物交叉重组,共形成6对引物,分别进行PCR扩增,扩增后回收纯化cDNA序列 (质量浓度≥167 ng·μL−1)。利用T7 High Efficiency Transcription 试剂盒 (全式金,北京) 将cDNA逆转录成单链RNA,混匀,37 ℃ PCR孵育2 h。将同一片段的两管单链RNA在冰上等量混匀,将混匀物放置在70 ℃ PCR仪中反应10 min,反应结束立刻取出后室温静置20 min。向每20 μL的反应液加入1 μL RNase A (用无酶水稀释200倍) 和1 μL DNase I,放置在37 ℃ PCR仪中反应30 min。这一步骤是为了去除DNA,使之退火合成dsRNA。将所得溶液中加入1/10体积的3 mol·L−1醋酸铵和等量异丙醇,冰中静置5 min;放置离心机中12 000 r·min−1、4 ℃、10 min。弃上清,70%乙醇洗涤沉淀;12 000 r·min−1、4 ℃离心10 min。超净工作台内风干15 min;加入100 μL RNase-free water溶解。检测dsRNA的质量和浓度,−80 ℃保存。
1.7.2 dsRNA干扰链注射
取1龄三角帆蚌按照每组15只分为4组 (干扰链1、干扰链2、干扰链3和对照组),实验组每只注射50 μL (200 ng·μL−1) 干扰链,阴性对照组每只注射50 μL (200 ng·μL−1) 生理盐水 (Normal saline, NS)。在第12、第24和第48小时提取性腺组织并立刻放置于液氮中。利用TRIzol法提取总RNA,试剂盒将其反转为cDNA,将反转后的cDNA作为模板进行荧光定量实验,检测直接干扰C-MOS基因后在性腺中的表达量情况及下游基因MAP2K1的表达量情况。
2. 结果
2.1 MAP2K1基因cDNA的克隆及序列特征
本实验克隆得到三角帆蚌MAP2K1基因cDNA的序列,其中5' 非编码区 (UTR) 为145 bp,3'UTR为2 070 bp,开放阅读框 (ORF) 1 194 bp,编码397个氨基酸 (图1)。预测相对分子质量44.24 kD,等电点为6.81,平均亲水系数为−0.379,推测为亲水蛋白。根据MAP2K1蛋白既不具备跨膜结构、也不具备信号肽的特点,推测其不属于膜蛋白,而属于胞内蛋白。通过SOPMA分析得到MAP2K1基因具有S-TKC结构域 (72—372 aa,图2-a),通过I-TASSER预测MAP2K1基因的蛋白质三级结构 (图2-b)。
2.2 MAP2K1基因系统进化
利用GeneDoc进行不同物种之间的序列比对,结果显示来自多种物种的MAP2K1基因结构域重合度很高 (图3),说明在多种物种之间,MAP2K1基因进化和结构具有保守性,推断其具有功能相似性。
根据同源性比对结果分析,三角帆蚌MAP2K1基因与紫贻贝、虾夷扇贝、泥蚶 (Tegillarca granosa) 的同源性较高,同源性分别为85.52%、82.65%、80.59%。与智人 (Homo sapiens)、北美属兔 (Ochotona princeps) 等同源性较低,同源性分别为72.27%、72.53%。系统进化树显示三角帆蚌与紫贻贝、虾夷扇贝、泥蚶聚为一支,亲缘关系较近 (图4)。
2.3 MAP2K1基因在各组织和各时期的表达量差异
在雌雄三角帆蚌的不同组织中,MAP2K1基因均有一定的表达量,但是MAP2K1基因在雌性性腺中的相对表达量高于另5个组织,雌雄差异也比较显著 (P<0.05,图5)。
在三角帆蚌幼龄 (1—8月) 性腺发育阶段过程中,MAP2K1基因均有一定程度的表达 (图6)。说明MAP2K1基因在三角帆蚌幼龄性腺发育阶段起到一定作用。
三角帆蚌1龄性腺中,雌性相对表达量显著高于雄性 (P<0.01),但在雄性中也有一定表达 (图7);等到2龄左右,性腺相对表达量雌雄差异极显著 (P<0.01),雄性中几乎不表达;在3龄性腺中,雌雄表达量均处于非常低的水平。结果表明,MAP2K1基因在三角帆蚌2龄卵巢发育过程中起到非常重要的作用。
2.4 性腺组织原位杂交
结果显示,实验组中三角帆蚌雌性性腺的卵母细胞、卵子上,蓝紫色杂交信号非常显著。实验组中三角帆蚌雄性性腺的精原细胞、精母细胞上,蓝紫色杂交信号存在极少,雌性及雄性阴性对照组中均未出现杂交信号 (图8)。根据结果推测,MAP2K1基因可能参与三角帆蚌的卵巢发育过程。
2.5 RNA干扰结果
三条干扰链的干扰结果为,干扰链1在雌性和雄性中的干扰率分别为88.43%和77.90%;干扰链2在雌性和雄性中的干扰率分别为93.10%和78.10%;干扰链3在雌性和雄性中的干扰率分别为96.39%和86.26% (图9)。干扰率结果显示干扰链3对C-MOS基因的干扰效果最佳,为此选择干扰链3进行实验,检测C-MOS基因在12、24和48 h干扰后在性腺中的表达情况,结果显示干扰链3在第24小时的干扰效率最佳 (图10)。
干扰C-MOS基因后,荧光定量检测其下游基因MAP2K1的表达量情况 (图11)。结果显示,MAP2K1基因在雌性中的表达量下降了82.31%,在雄性中的表达量下降了73.60%。说明在三角帆蚌中,C-MOS基因对MAP2K1基因有显著的调控作用。
3. 讨论
本研究克隆了三角帆蚌MAP2K1基因的cDNA序列,编码一个长为1 194 bp的ORF。蛋白质结构预测分析显示,MAP2K1基因中含有S-TKC催化结构域。蛋白苏氨酸/丝氨酸激酶是十分常见的一类蛋白激酶[21],主要作用是蛋白磷酸化、去磷酸化,有利于细胞内活性[22]。S-TKC是这类激酶中一个位于核心部位重要的催化结构域[23]。在S-TKC中有一个极其保守的区域,在催化结构域它已被证明是涉及ATP的结合,即赖氨酸残基附近N-末端的一个富含甘氨酸的残基。在S-TKC的中间区域,存在一个天冬氨酸残基,在催化结构域它已被证明对酶的催化活性起显著作用[24-26]。拥有S-TKC催化结构域的基因通常参与许多生理生化过程[27-30],若其构造折叠或变化,会对丝氨酸/苏氨酸蛋白激酶产生影响,从而对个体的生殖发育过程造成极大损害[31],这说明MAP2K1基因可能参与了三角帆蚌个体的生殖发育过程。
不同物种之间的序列比对结果显示,MAP2K1基因在三角帆蚌中的氨基酸序列与其他物种结构域重合度较高,说明与其他物种相似度较高,尤其与虾夷扇贝等软体动物相似度更高,由此推断在进化过程中,MAP2K1基因呈现高度保守状态。各组织、早期发育阶段、不同发育时期的定量结果说明,MAP2K1基因在卵巢中特异性高表达,在1—8月龄时期均有不同程度的表达,在1、2、3龄卵巢中的表达量均高于睾丸,对小鼠[14]、水牛[15]、鹅[5]、斑节对虾[14]及果蝇[16]MAP2K1基因的研究中,均发现雌性表达量高于雄性,与本实验结果相近,由此推测MAP2K1基因是三角帆蚌中的偏雌性基因,既参与了三角帆蚌的早期发育阶段,也参与了卵巢发育过程。根据原位杂交的探针定位显示,MAP2K1基因在卵母细胞及卵子上存在明显的信号,推测MAP2K1基因可能参与了三角帆蚌的卵巢发育过程。
RNA干扰结果显示,3条干扰链中干扰链2和干扰链3对C-MOS基因起到抑制作用,但是干扰链3的干扰效果最佳,在雌性和雄性中的干扰率分别为84.21%和96.87%。在选取干扰链3观察干扰C-MOS基因后MAP2K1基因的表达情况中,MAP2K1基因在雌性中的表达量下降了82.31%,在雄性中的表达量下降了73.60%。结果显示在三角帆蚌中,C-MOS基因对MAP2K1基因有显著的调控作用。在脊椎动物卵母细胞MAPK信号通路中,MAPK独立于生长因子和酪氨酸激酶受体而被激活,受特殊上游调节剂激酶MOS的控制,在卵母细胞减数分裂中起关键作用[9]。综上,推测在三角帆蚌中,MAP2K1基因可能参与了三角帆蚌的卵巢发育过程,是偏雌性基因。
-
表 1 热带普通小球藻生长过程中颜色参数与细胞密度、色素参数质量浓度的简单相关分析
Table 1 Simple correlation analysis among color parameters and cell density, pigment contents of C. vulgaris in process of growth
参数
Parameter细胞密度
Cell density叶绿素a
Chl-a叶绿素b
Chl-b总叶绿素
Chlorophyll类胡萝卜素
Carotemoid catabolite类叶比值
Ratio of carotenoid and
chlorophyll明度值L* Lightness value −0.734* −0.912** −0.919** −0.918** −0.981** −0.913** 红度值a* Red value −0.918** −0.940** −0.957** −0.949** −0.939** −0.948** 黄度值b* Yellow value 0.857** 0.943** 0.954** 0.950** 0.970** 0.948** 色差值△E Color difference 0.845** 0.941** 0.952** 0.948** 0.974** 0.946** 相角H° Hue angle −0.766** −0.689* −0.703* −0.696* −0.639* −0.774** 饱和度C Chroma 0.859** 0.943** 0.955** 0.950** 0.970** 0.948** 注:*. 在0.05水平 (双侧) 上显著相关;**. 在0.01水平 (双侧) 上显著相关 Note: *. Significant correlation (Bilateral) at 0.05 level; **. Very significant difference (Bilateral) at 0.01 level. 表 2 热带普通小球藻生长过程中细胞密度、色素参数质量浓度与颜色参数的回归分析
Table 2 Analysis of regression among cell density, pigment contents and color parameters of C. vulgaris in process of growth
参数
Parameter回归方程
Regression equation调整R2
Adjusted R2F P 细胞密度 Cell density y细胞密度=−282.126 3+2.170 2xL−198.395 5xa−26.132 0xb+2.023 4x△E 0.971 6 43.093 2 0.000 2 叶绿素a Chl-a y叶绿素a=11.568 4−0.113 5xL−0.263 4 xa+0.0940 1xb−0.159 5x△E 0.914 0 13.693 9 0.003 6 叶绿素b Chl-b y叶绿素b=5.740 5−0.056 4xL−0.174 5xa−0.004 4xb−0.040 8x△E 0.953 1 25.768 9 0.000 6 总叶绿素 Chlorophyll y总叶绿素=17.297 0−0.169 7xL−0.437 4xa+0.089 7xb−0.200 3x△E 0.929 9 16.990 1 0.002 0 类胡萝卜素 Carotemoid catabolite y类胡萝卜素=9.873 3−0.096 7xL−0.183 9xa+0.019 94xb−0.069 78x△E 0.973 1 45.680 3 0.000 1 总色素 Plamochromic pigment y总色素=27.170 4−0.266 4xL−0.621 4xa+0.109 7xb−0.270 1x△E 0.945 0 25.760 0 0.000 6 -
[1] 鲍亦璐. 微藻培养过程的营养优化与控制研究[D]. 广州: 华南理工大学, 2012: 1-2. [2] MOUHIM R F, CORNET J F, FONTANE T, et al. Production, isolation and preliminary characterization of the exopolysaccharide of the cyanobacterium Spirulina platensis[J]. Biotechnol Lett, 1993, 15(6): 567-572. doi: 10.1007/BF00138541
[3] OTERO A, DIGNA G, MORALES E D, et al. Manipulation of the biochemical composition of the eicosapentaenoic acid-rich microalga Isochrysis galbana in semicontinuous cultures[J]. Biotechnol Appl Bioc, 2010, 26(3): 171-177.
[4] PULZ O, GROSS W. Valuable products from biotechnology of microalgae[J]. Appl Microbiol Biot, 2004, 65(6): 635-648. doi: 10.1007/s00253-004-1647-x
[5] JIANG L, LUO S, FAN X, et al. Biomass and lipid production of marine microalgae using municipal wastewater and high concentration of CO2[J]. Appl Energ, 2011, 88(10): 3336-3341. doi: 10.1016/j.apenergy.2011.03.043
[6] TAKAGI M, KARSENO Y T. Effect of salt concentration on intracellular accumulation of lipids and triacylglyceride in marine microalgae Dunaliella cells[J]. J Biosci Bioeng, 2006, 101(3): 223-226. doi: 10.1263/jbb.101.223
[7] 潘冰峰, 李祖义. 利用海洋微藻生产富含DHA的单细胞油脂[J]. 生物工程进展, 2000(6): 43-45. doi: 10.3969/j.issn.1671-8135.2000.06.010 [8] 丁丽欢, 王珺, 严虹羽, 等. 热带普通小球藻培养模式的筛选及其培养基的优化[J]. 微生物学通报, 2018, 45(5): 981-989. [9] SAFI C, ZEBIB B, MERAH O, et al. Morphology, composition, production, processing and applications of Chlorella vulgaris: a review[J]. Renew Sust Energ Rev, 2014, 35: 265-278. doi: 10.1016/j.rser.2014.04.007
[10] 常杰. 小球藻在动物营养中的应用研究进展[J]. 饲料研究, 2016(2): 7-10. [11] 过伟民, 魏春阳, 张艳玲, 等. 烤烟表面颜色的量化及其与胡萝卜素类物质的关系[J]. 烟草科技, 2012(1): 62-68. doi: 10.3969/j.issn.1002-0861.2012.01.013 [12] 丁根胜, 张庆明, 巴金莎, 等. 烟叶颜色色度学指标与烤烟品质的关系分析[J]. 中国烟草科学, 2011, 32(4): 14-18. doi: 10.3969/j.issn.1007-5119.2011.04.004 [13] 谢鹏飞, 邓小华, 周清明, 等. 密集烘烤过程中烟叶颜色、形态和水分变化及相互关系[J]. 作物研究, 2012, 26(5): 486-490. doi: 10.3969/j.issn.1001-5280.2012.05.19 [14] 曹云涛, 孔维宝, 葸玉琴, 等. 不同营养方式对普通小球藻生长特性和细胞组成的影响[J]. 食品与发酵工业, 2011, 37(10): 45-51. [15] 葸玉琴, 崇梅, 朱巧巧, 等. 不同氮源对混养小球藻生长和部分生化组成的影响[J]. 西北师范大学学报(自然科学版), 2015, 51(1): 82-86. [16] 李小嫄, 王洪伟, 童华荣. 食品感官评价技术在茶叶品质评价中的应用研究进展[J]. 食品安全质量检测学报, 2015, 6(5): 1542-1547. [17] 孟凡娟, 许向阳, 李景富. 番茄果实色素含量和表面颜色相关性研究[J]. 东北农业大学学报, 2006(4): 459-462. doi: 10.3969/j.issn.1005-9369.2006.04.007 [18] TIMOTHY J, GRAEME E, HOBSON A. General method for the instrumental assessment of the color tomato fruit during ripening[J]. J Sci Food Agric, 1984, 35(12): 1277-1281. doi: 10.1002/jsfa.2740351202
[19] 蔡鸿昌, 崔海信, 宋卫堂, 等. 黄瓜初花期叶片光合色素含量与颜色特征的初步研究[J]. 农业工程学报, 2006(9): 34-38. doi: 10.3321/j.issn:1002-6819.2006.09.007 [20] 于国锋, 冯媛, 郑宏伟, 等. CIE-L*a*b*颜色模型在烟草中的应用研究进展[J]. 河南农业科学, 2019, 48(2): 1-7. [21] TAIZ L, ZEIGER E, MOLLER I, et al. Plant physiology and development[M]. 6th ed. Massachusetts: Sinauer Associates Inc, 2016: 10-12.
[22] 田喜, 何绍兰, 吕强, 等. 高光谱图像信息的柑橘叶片光合色素含量分析技术研究[J]. 光谱学与光谱分析, 2014, 34(9): 2506-2512. doi: 10.3964/j.issn.1000-0593(2014)09-2506-07 -
期刊类型引用(1)
1. 路俊怡,姜立妍,龙凯,王韬,吴正理,李艳红. HcTLR1通过MyD88-NF-κB信号通路参与三角帆蚌抗菌免疫应答. 南方水产科学. 2024(06): 19-30 . 本站查看
其他类型引用(0)