Impacts of different ecological types of earthworm on aggregate distribution and stability in typical latosolic red and red soils
-
摘要:目的
综合分析华南地区不同生态类型蚯蚓对土壤团聚体分布、团聚体水稳定性以及土壤团聚体破坏率的影响，以期为蚯蚓改善土壤结构提供科学依据和理论支持。
方法采用室内盆栽培养试验，在华南地区具有代表性的赤红壤和红壤中分别接种生态类型不同的4种蚯蚓：赤子爱胜蚓Eisenia fetida、南美岸蚓Pontoscolex corethrurus、壮伟远盲蚓Amynthas robustus和参状远盲蚓A. aspergillum，研究不同生态类型蚯蚓对土壤结构的影响。
结果与空白对照相比，南美岸蚓、壮伟远盲蚓和参状远盲蚓的添加均显著增加了2种土壤中的大团聚体(d>2 000 μm)比例(P<0.05)，其中，赤红壤中其分别提升了35.20%、44.81%和37.88%，红壤中其分别提升了14.92%、25.31%和20.18%；与空白对照相比，壮伟远盲蚓的添加均显著降低了赤红壤水稳性大团聚体和水稳性小团聚体(250 μm≤d≤2 000 μm)的比例，却显著提升了赤红壤水稳性微团聚体(d<250 μm)的比例(37.84%，P<0.05)；4种蚯蚓的添加均显著提升了2种土壤的团聚体破坏率(P<0.05)，其中添加参状远盲蚓的2种土壤的团聚体破坏率均最低。主成分分析结果显示：在2种不同的土壤中，不同蚯蚓作用后土壤的团聚体分布和稳定性之间存在显著差异(P<0.05)。
结论华南不同生态类型的蚯蚓对赤红壤和红壤的团聚体分布和稳定性的影响不同；内栖型蚯蚓(南美岸蚓和壮伟远盲蚓)对土壤团聚体的形成和土壤结构的改善效果最好；深栖型蚯蚓(参状远盲蚓)对土壤水稳性团聚体的破坏作用最小。
Abstract:ObjectiveTo comprehensively analyze the effects of different ecological types of earthworm on soil aggregate distribution, water stability of aggregate and percentage of aggregate disruption (PAD) in South China, and provide scientific basis and theoretical support for utilizing earthworms to improve soil structure.
MethodThe laboratory culture experiment was carried out in this study. Eisenia fetida, Pontoscolex corethrurus, Amynthas robustus and A. aspergillum were added in typical latosolic red soil and red soil in South China to study the impacts of different ecological species of earthworm on soil structure.
ResultCompared with the blank control, the addition of P. corethrurus, A. robustus and A. aspergillum significantly promoted the formation of large-aggregates (d > 2 000 μm) in two soil types ( P < 0.05), The large-aggregate proportion in latosolic red soil increased respectively by 35.20%, 44.81% and 37.88%, and increased respectively by 14.92%, 25.31% and 20.18% in red soil. A. robustus significantly reduced the contents of water stable large-aggregates (d > 2 000 μm) and small-aggregates (250 μm≤ d≤2 000 μm) in latosolic red soil, but significantly increased the proportion of water stable micro-aggregates (d < 250 μm, P < 0.05), which was 37.84% higher than that of the blank control. The addition of four types of earthworms significantly increased PAD values in both two soil types ( P < 0.05), and the PAD value in adding A. aspergillum treatment was the lowest. The results of principal component analysis showed that there were significant differences in the distributions and stabilities of soil aggregates under the action of different earthworms in two different soils (P < 0.05).
ConclusionThe impacts of different earthworm types in South China on latosolic red soil and red soil aggregate distribution and stability are different. The endogenic earthworm (P. corethrurus and A. robustus) has the best impact on the improvement of soil structure. The endo-anecic (A. aspergillum) had the least destructive impact on water stable soil aggregate.
-
Keywords:
- earthworm /
- soil structure /
- soil aggregate /
- distribution /
- stability /
- percentage of aggregate disruption
-
染色体作为生物遗传物质的载体，对当前生命遗传研究开展尤为重要。染色体消除技术(Chromosome elimination technology)近年来随着基因编辑技术的发展不断进步，在动物模型构建、家畜育种改良中的性别控制、特定染色体的功能研究以及人类非整倍体疾病治疗等领域具有重要的应用价值。目前通过该技术已成功在人、小鼠、斑马鱼等物种上实现了特定染色体的消除[1-3]，然而打靶效率低、非靶点的额外失活及操作过程的繁琐性大大制约了该技术的发展[3-5]。但随着新一代CRISPR/Cas9技术的出现，灵活的靶点选择、较高的编辑效率以及对于靶点多重切割的可能性使该技术的有效运用成为可能。
目前，研究者们已通过锌指核酸酶(Zinc finger nuclease，ZFN)技术、Cre/loxP系统、TKNEO基因敲入等多种方法实现了目标染色体的敲除[5-7]。ZFN技术通过人工改造的核酸内切酶(锌指蛋白单元)与特定的DNA序列结合赋予其精确的打靶功能[8]，但其脱靶效应、毒性及锌指结构的低适用性造成其基因操作效率低下[9]；Cre/loxP系统主要由Cre重组酶与2个loxP位点组成，由于loxP的位点不同会产生DNA链的插入、删除、易位等不同重组方式[10]，但该系统的使用需要引入特定表达元件且本身会对细胞产生毒性[4, 11]；TKNEO是胸苷激酶(Thymidine kinase，TK)和新霉素抗性(Neomycin resistance，NEO)的融合编码基因，利用传统基因打靶技术敲入目标染色体进行药筛可产生自发性染色体缺失[3]，但其操作的繁琐性和不确定性使其并不具有广泛的适用性。
CRISPR/Cas9技术近年来已经成为了基因编辑的主要手段，其由大肠埃希菌Ⅱ型CRISPR/Cas系统改造而来，由于其引导RNA(Small guide RNA，sgRNA)序列的灵活性而能对多个位点进行灵活编辑，从而实现基因的定点敲入、敲除或其他多种修饰[12-16]。该系统主要由2个部分组成，Cas9蛋白和sgRNA，其中后者负责识别目标靶序列，前者依靠其自身含有的2个切割域完成对目标序列的切割，之后依靠细胞本身的2种修复机制(同源重组修复和非同源末端连接)实现目标序列的精确编辑[17-21]。最新研究表明，CRISPR/Cas9系统可通过产生多重DNA双链断裂(Double-strand breakage，DSB)引起Y染色体降解[22-23]，对这种方法，虽然部分细胞会因自身修复机制保持基因组的稳定性，但是仍有相当的几率可产生特定染色体消除的细胞，通过筛选出此种染色体消除的细胞系作为供体应用于动物克隆，可制备出特定染色体消除的动物，因此达到控制动物本身及其后代的性别进行偏移的目的。
本研究拟通过CRISPR/Cas9技术，在猪Y染色体内寻找能够被特异sgRNA识别的重复序列，进行多位点切割以实现目标染色体的降解。本研究旨在为特定染色体及其携带基因功能的研究提供一种新的工具，同时也为实现大型家畜性别控制提供一种新的技术方案。
1. 材料与方法
1.1 材料
本试验所用的CRISPR/Cas9载体HP180_px330_GFP及猪胎儿成纤维细胞(Porcine embryonic fibroblast，PEF)由国家生猪种业工程技术研究中心保存。设计的sgRNA及相关引物由华大基因合成。质粒与基因组DNA抽提试剂盒、限制性内切酶BpiI、gRNA体外转录试剂盒GeneArtTM Precision gRNA Synthesis Kit、荧光定量试剂盒PowerUpTMSYBRTMGreen Master Mix购自赛默飞公司；T4 DNA连接酶购自TaKaRa公司；Cas9内切酶购自NEB公司；感受态细胞Trans5α、6×DNA Loading buffer、RNA loading buffer购自全式金生物。高糖DMEM 培养液、胰酶、胎牛血清均购自Gibco公司。
1.2 方法
1.2.1 寻找Y染色体特异靶向重复序列及sgRNA设计
在Ensembl获取整个Y染色体注释文件，在Linux下用awk命令提取注释文件中的所有基因，并在Ensembl上将每个基因与全基因组序列进行逐一比对，最终筛选出符合条件的多拷贝基因。
利用在线网站CCTop - CRISPR/Cas9 target online predictor( https://cctop.cos.uni-heidelberg.de/)对所选择靶基因进行sgRNA设计，确定最终序列与位点，并送至华大基因合成。
1.2.2 sgRNA体外切割验证
对设计好的sgRNA体外转录进行切割，初步筛选出效率较高的进行载体构建。首先进行靶片段的PCR扩增。根据本试验靶基因的选择，使用Primer Premier 6进行相应靶点引物设计，合成目标靶片段。然后根据赛默飞公司的GeneArtTM Precision gRNA Synthesis Kit说明书设计sgRNA体外转录模板引物，并最终进行Gene_ID为ENSSSCG00000043058、ENSSSCG00000046447的转录模板合成(以下简称为43058、46447)。其合成原理(图1)及相关引物(表1、表2)如下所示。
![]()
图 1 sgRNA合成原理图体外合成的gRNA由2条34~38 bp寡核苷酸与tracrRNA片段退火形成，其中Target-F由T7启动子和靶序列的16~20 bp组成，Target-R由tracrRNA片段的部分序列和与目标序列互补的19~20 bp组成Figure 1. Schematic diagram of sgRNA synthesisTwo 34−38 bp oligonucleotides are needed to assemble the synthetic gRNA with tracrRNA fragment, among them the Target-F is composed of the T7 promoter and 16−20 bp sequence of the target sequence, and the Target-R is composed of partial sequence of tracrRNA fragment and 19−20 bp sequence complementary to the target sequence表 1 靶片段扩增引物Table 1. Primers of target fragment amplification引物名称 Primer name 引物序列(5′→3′) Primer sequence 产物长度/bp Product length 43058target_F GTCCCATTTTCCAGGCCTTA 616 43058target_R GTTCAGAGGCACACAATGTG 46447target_F AGTTGAAGGCCACTTGGTCA 657 46447target_R ATGACAATTTGTCCTGGTGGCC 表 2 sgRNA体外转录模板合成引物Table 2. Primers for synthesis of sgRNA in vitro transcription template引物名称 Primer name 引物序列(5'→3') Primer sequence T43058-F1 TAATACGACTCACTATAGTCACTGGGGACTTGGACTT T43058-R1 TTCTAGCTCTAAAACCAAGTCCAAGTCCCCAGTG T43058-F2 TAATACGACTCACTATAGGTGCATCTTTAGGAGACAC T43058-R2 TTCTAGCTCTAAAACGTGTCTCCTAAAGATGCACC T46447-F1 TAATACGACTCACTATAGATGTGGCTGCAGAAACTCG T46447-R1 TTCTAGCTCTAAAACCGAGTTTCTGCAGCCACAT T46447-F2 TAATACGACTCACTATAGAAGAAGCGCTCTAGAACAGG T46447-R2 TTCTAGCTCTAAAACCCTGTTCTAGAGCGCTTCT 体外Cas9切割验证主要包含3个过程：Cas9-sgRNA的自发组装、Cas9-sgRNA复合体与靶片段的特异结合、Cas9-sgRNA对靶点进行特异切割。采用20 μL反应体系：体外转录sgRNA模板100 ng，靶片段250 ng，Cas9蛋白1 μL，10×Cas9 Reaction buffer 2 μL，最后加Nuclease-free water补足至20 μL。
在进行Cas9体外切割时，反应体系首先加入除靶片段外其他成分进行Cas9与sgRNA组装(组装时间30 min)，再加入相应靶片段孵育1 h，切割完毕后在反应体系中加入2 μL的蛋白酶K再次孵育20 min。最终切割产物加入6×DNA Loading Buffer进行20 g/L的凝胶电泳，分析电泳条带灰度计算切割效率。
1.2.3 CRISPR/Cas9载体构建
确定有效的sgRNA后进行真核细胞表达载体构建。首先使用限制性内切酶BpiI对载体HP180_px330_GFP(图2)双酶切使其线性化，合成的sgRNA经退火形成双链由T4连接酶在16 ℃条件下连接1 h，然后热应激转入Trans5α感受态细胞，于220 r/min、37 ℃培养1 h后涂布含氨苄青霉素的LB平板，隔夜培养后挑取单克隆，进行菌液PCR并送测序，测序正确菌液进行质粒抽提，用于PEF电转。
1.2.4 PEF培养与转染
试验所用细胞经鉴定为公猪细胞(使用引物F：TCATAGCTCAAACGATGGACGTG，R：ACACAATGAAAGCGTTCATGGGTC，扩增片段大小为92 bp)。PEF复苏后加入含体积分数为12%胎牛血清的高糖DMEM培养液，置于38 ℃、CO2体积分数为5%的培养箱中培养。当复苏后的细胞融合度达到80%~90%时，用质量分数为0.05%的胰酶消化后进行细胞计数。取1×106 mL−1细胞悬液加入1.5 mL离心管，加入100 μL电转液及目的质粒以Lonza核转仪进行电转，电转程序A003。转染24 h后，荧光显微镜下拍照并观察转染细胞荧光信号强弱，判断转染效果。
1.2.5 流式分选及单细胞克隆团鉴定
细胞转染2~3 d后进行流式分选(Fluorescence activating cell sorter，FACS)。细胞分选前经质量分数为0.05%的胰酶消化，加入含体积分数为2%胎牛血清的 DPBS重悬，并在上机前经40 µm细胞筛过滤防止细胞团块的干扰。根据GFP标记激发波长488 nm分选出阳性细胞(5×104~10×104个)并进行极限稀释(梯度稀释降低细胞密度)培养，细胞培养形成单个细胞团后取部分细胞抽提DNA进行Y染色体性别决定区(Sex determination region of Y chromosome, SRY)基因鉴定。引物信息如表3所示。
表 3 靶基因及SRY基因PCR引物Table 3. PCR primers of target gene and SRY gene引物名称 Primer name 引物序列(5′→3′) Primer sequence 产物长度/bp Product length 43058_F CATTGTGATGCATGGGCGATTT 238 43058_R GAGACCTCAAGTCCAAGTCCC 46447_F ATGTGGCTGCAGAAACTCGT 169 46447_R ATTTGTCCTGGTGGCCTGAC SRY_F GTGGCTGGGATGCAAGTGGAAA 228 SRY_R GCCTTGGCGACTGTGTATGTGAA 1.2.6 实时定量荧光PCR检测Y染色体敲除效率
分选后待细胞融合度达到80%~90%提取基因组DNA，并设计相应靶基因及SRY定量引物(表3)观察其DNA含量是否显著降低。以基因组为模版检测目的基因在DNA水平的含量是否变化，并以基因组水平的β-actin为内参。反应体系与反应产物均按照相应试剂盒(PrimeScriptTM RT reagent Kit with gDNA Eraser与PowerUpTM SYBRTM Green Master Mix)说明进行。
1.2.7 核型鉴定
对SRY鉴定为阴性细胞继续扩繁培养，并传代至T25培养瓶，待细胞生长至融合度为70%~80%时进行核型分析，判断染色体核型是否缺失。核型分析由杭州卡优太谱生物科技公司进行。
2. 结果与分析
2.1 多拷贝基因查询及sgRNA设计结果
2.1.1 多拷贝基因查询结果
通过Ensembl的blast筛选出了35个位于Y染色体上的多拷贝基因，统计结果如表4所示。
表 4 猪Y染色体多拷贝基因信息1)Table 4. Muti-copy genes in Y chromosome of pig基因编号 Gene_ID 起始位置 Begin position 结束位置 End position 长度/bp Length 拷贝数 Copy number ENSSSCG00000041799 20570323 20570724 402 90 ENSSSCG00000045041 23748599 23749006 408 42 ENSSSCG00000043614 22488485 22488925 441 39 ENSSSCG00000043058 21558744 21559151 408 39 ENSSSCG00000050250 20115553 20117373 1821 32 ENSSSCG00000049276 20156946 20158766 1821 32 ENSSSCG00000041211 20259563 20261383 1821 32 ENSSSCG00000041323 20851299 20853119 1821 32 ENSSSCG00000051247 23675355 23677155 1821 32 ENSSSCG00000044850 22938437 22940148 1721 32 ENSSSCG00000035359 19637403 19639087 1685 32 ENSSSCG00000038218 20710734 20712418 1685 32 ENSSSCG00000031994 20902449 20904133 1685 32 ENSSSCG00000032210 21694455 21696139 1685 32 ENSSSCG00000033039 19929688 19931370 1683 32 ENSSSCG00000044818 21140962 21142262 1301 32 ENSSSCG00000039954 22859918 22861602 1685 31 ENSSSCG00000049743 23367269 23368953 1685 31 ENSSSCG00000042334 19885774 19887455 1682 31 ENSSSCG00000042971 21982938 21984619 1682 31 ENSSSCG00000045651 23897691 23899366 1676 31 ENSSSCG00000051215 20442628 20444448 1821 30 ENSSSCG00000041507 22640161 22641981 1821 30 ENSSSCG00000045761 24007552 24009372 1821 29 ENSSSCG00000041750 23944189 23946009 1821 29 ENSSSCG00000050201 21301627 21303492 1821 29 ENSSSCG00000048908 19541113 19542933 1821 27 ENSSSCG00000045241 21443949 21445655 1707 26 ENSSSCG00000046033 23455676 23457401 1726 24 ENSSSCG00000036112 22537239 22541180 3942 3 ENSSSCG00000043960 23320436 23322365 1930 2 ENSSSCG00000031258 21822265 21828538 6274 * ENSSSCG00000051413 21205077 21216797 11721 * ENSSSCG00000046447 22138546 22177691 39146 * ENSSSCG00000042484 23810416 23850459 40244 * 1) “*”代表基因序列太长，不存在完全拷贝片段，拷贝数统计不完整，但存在部分序列多拷贝，因此也可作为多拷贝基因的一种；加粗Gene_ID为靶基因 1) “*” represents that the gene sequence is too long, there is no complete copy fragment, and the copy number statistics are incomplete, but there are multiple copies of some sequences, so it can also be used as a kind of multi copy gene；Bold gene ID represents the target gene 通过再次与NCBI数据库进行比对，以及基于靶点基因组测序的真实结果与sgRNA设计的实际情况，我们最终选择了网站预测评分较高、符合测序结果且在Y染色体上重复拷贝数靠前(多位点DSB)的Gene_ID为ENSSSCG00000043058、ENSSSCG00000046447的基因进行了sgRNA设计。这2个基因在全基因组的比对结果如图3所示。
![]()
图 3 ENSSSCG00000043058和ENSSSCG00000046447在猪染色体上的位置示意图所选择的靶基因仅存在于Y染色体上并有多个拷贝，有利于实现Y的多位点切割Figure 3. Locations of ENSSSCG00000043058 and ENSSSCG00000046447 on pig chromosomeThe selected target genes only exist on the Y chromosome and have multiple copies, which is conducive to the realization of Y multipoint cutting2.1.2 sgRNA设计结果
对ENSSSCG00000043058和ENSSSCG00000046447位点进行sgRNA设计(表5)。在sgRNA及其互补序列末端分别加上与BpiI酶切HP180_px330_GFP质粒互补的黏性末端，由华大基因合成相应的Oligo。
表 5 sgRNA靶点位置和序列Table 5. Location and sequence of sgRNA targetsgRNA Y染色体位置 Position on Y chromosome 靶点序列 Target sequence 43058sg1 21558841~21558860 TCACTGGGGACTTGGACTTG 43058sg2 21558795~21558814 ATGGAGGATTCTGAACCAGT 46447sg1 22254580~22254599 GATGTGGCTGCAGAAACTCG 46447sg2 22254625~22254644 AAGAAGCGCTCTAGAACAGG 2.2 sgRNA体外切割结果
以野生型PEF DNA为模板进行PCR扩增，产物进行凝胶电泳分析后切胶回收得到sgRNA体外转录模板。根据GeneArtTM Precision gRNA Synthesis Kit试剂盒的体系与方法进行sgRNA的转录模板合成(图4A)与体外sgRNA转录(图4B)，转录产物经纯化后电泳验证条带完整性(长度100 bp)，电泳结果与预期一致。
![]()
图 4 靶片段模板DNA扩增与sgRNA体外转录结果M：DL2000 DNA marker；1和2分别代表扩增靶片段43058和46447，长度分别为616和657 bp；3、4、5和6代表转录产物43058sg1、43058sg2、46447sg1和46447sg2，长度均为100 bpFigure 4. Amplification of target fragment template DNA and in vitro transcription of sgRNAM:DL2000 DNA marker; 1 and 2 represent the amplified target fragment 43058 and 46447 respectively, and the lengths are 616 and 657 bp respectively; 3, 4, 5 and 6 represent transcripts 43058sg1, 43058sg2, 46447sg1 and 46447sg2 respectively, and the lengths all are 100 bp将扩增好的靶片段、sgRNA按照一定比例与Cas9混合，按照方法中的步骤进行反应得到切割产物，43058sg1的理论切割产物长度为525和91 bp，43058sg2的理论切割产物长度为444和172 bp，46447sg1的理论切割产物长度为498和159 bp，46447sg2的理论切割产物长度为453和204 bp。凝胶电泳分析表明43058sg2与46447sg2均达到了预期的切割效果，并且43058sg2实现了接近完全的靶基因切割(图5)。灰度分析表明，46447sg2的体外切割效率也达到了58%(图6)，最终我们选择 43058sg2和46447sg2进行载体构建。
![]()
图 5 设计的sgRNA体外切割电泳图M：DL7000 DNA marker；1、4：未转染细胞对照；2、3、5、6为设计的43058sg1、43058sg2、46447sg1以及46447sg2的实际切割产物Figure 5. Electrophoregram of in vitro cleavage of designed sgRNAM: DL7000 DNA marker; 1, 4: Control of untransfected cells; 2, 3, 5, 6 are the actual cutting products of designed 43058sg1, 43058sg2, 46447sg1 and 46447sg2 respectively2.3 HP180_px330_GFP载体构建结果
HP180_px330_GFP骨架上有预留的BpiI双酶切位点，其线性化后与退火形成具有黏性末端的双链sgRNA，在连接酶作用下形成闭合环状结构。其测序结果如图7所示，结果与设计的sgRNA结果一致，表明载体构建成功。
![]()
图 7 HP180_px330_GFP质粒测序结果绿色下划线部分为BpiI双酶切位点，红色箭头标注部分为针对各靶位点构建的识别序列Figure 7. HP180_ px330_ GFP plasmid sequencing resultsThe underlined part in green is BpiI double digestion site, and the marked part in red arrow is the recognition sequence constructed for each target site2.4 细胞性别鉴定及细胞转染检测
2.4.1 性别鉴定
对试验所用PEF进行SRY鉴定，电泳结果显示目的细胞为雄性(图8)。
![]()
图 8 SRY基因检测结果M：DL700 DNA marker；1：空白水对照；2：阴性对照母猪DNA；3：试验所用公猪细胞Figure 8. SRY gene assay resultM：DL700 DNA marker;1: Blank water control; 2: DNA of negative control sow; 3: Boar cells used in the experiment2.4.2 细胞转染
将构建的CRISPR/Cas9载体通过细胞核转染导入雄性PEF中。由于载体带有增强绿色荧光蛋白(Enhanced green fluorescent protein，EGFP)表达单元，可通过观察EGFP的表达评估载体的转染效率。对已电转的细胞培养24 h后观察荧光信号，确定已构建好的载体成功电转进入细胞。结果(图9)表明质粒已成功表达。
![]()
图 9 PEF电转后荧光表达效果图A、B分别为2次重复试验的转染效果示意图；A1、B1为细胞明场图片，A2、B2分别为A1、B1同视野的EGFP表达图Figure 9. EGFP expression after PEF electroporationA and B are cell transfection results of two repeated experiments；A1 and B1 are transfected cells imaged in bright field, and A2 and B2 are EGFP expression images of the same vision as A1 and B1, respectively2.5 FACS分选阳性细胞及单克隆鉴定结果
通过FACS分选EGFP阳性细胞，我们发现对照组与试验组GFP信号存在显著差异，其中对照组未检测到GFP信号，而阳性对照的细胞转染效率达到了71%，这与荧光观察结果类似。通过培养EGFP阳性细胞，获得单细胞克隆团。接着我们对单克隆细胞团提取微量DNA，进行相应的SRY检测，发现1/3的细胞克隆未检测到SRY信号(图10)，显示这些细胞的Y染色体存在缺失的可能。
![]()
图 10 单克隆细胞团SRY基因鉴定示意图M：DL700 DNA marker；1~24代表单克隆细胞团编号Figure 10. SRY gene identification of monoclonal cell clustersM: DL700 DNA marker; 1−24 represent the number of monoclonal cell clusters2.6 实时定量荧光PCR检测目的基因及SRY基因
由于选用的目的基因在PEF上并不表达mRNA，因此采用qPCR检测目的基因(内参基因选择β-actin)和SRY基因的DNA含量。我们抽提EGFP阳性细胞的基因组DNA，进行qPCR分析目的基因(43058和46477)和SRY在DNA水平的含量。qPCR结果显示在进行靶基因敲除后，靶基因含量均显著降低。其中靶基因43058和46447分别下降了约20%和30%(P<0.05)，对应的SRY也下降了7%和11%。46447敲除效果相对于43058的基因敲除效果更为显著，不仅靶基因敲除效率较高，SRY基因DNA水平也下降较多(图11)。
![]()
图 11 EGFP阳性细胞中靶基因和SRY基因相较于正常雄性细胞的DNA水平检测“*”“**”“***”分别表示差异达到0.05、0.01和0.001的显著水平(t检验)Figure 11. Detection of DNA levels of target genes and SRY gene in EGFP-positive cells compared with normal male cells“*” “**” “***” indicate significant differences at 0.05, 0.01 and 0.001 levels respectively (t test)2.7 核型分析结果
核型鉴定结果显示上述流式分选得到的EGFP阳性细胞存在染色体缺失现象，表明设计的sgRNA在细胞水平上能够实现整个Y染色体的删除。在鉴定的12个单细胞核型中，有1个完全缺失Y染色体(图12)。
![]()
图 12 EGFP阳性细胞核型分析结果示意图Figure 12. Diagram of karyotype assay result of EGFP-positive cells3. 讨论与结论
CRISPR/Cas9近年来由于其基因编辑的灵活性和高效率而被广泛应用。其通过DSB导致的非保守非同源末端连接修复途径造成目的基因的敲除效果，但是该种途径常常造成异常的插入缺失，甚至破坏染色体和基因组的完整性[24-26]。而以天然的(或定向的)DNA模板的同源重组修复(Homologous recombination repair，HDR)途径则可以对目的基因位点进行精确修饰而产生意向的目的基因型[14, 27]，本研究利用CRISPR系统的高效切割能力对Y染色体上多个重复基因位点进行切割，通过产生多重DSB损伤基因组和DNA修复的稳定性，引起Y染色体降解。实现整条染色体删除是建立在CRISPR/Cas9系统能对多个位点进行同时编辑的基础上的，虽然早期的研究已经证明能够产生大片段的缺失，但这对于整条染色体的消融还远远不够。2015年，Yang等[28]首先通过靶向内源性逆转率病毒(PERVs)的62个重复位点在猪上实现了PERVs的全基因组降解；随后，研究人员进一步利用CRISPR/Cas9对靶点切割的多重性，成功在ES细胞、体内细胞和受精卵中删除了Y染色体，并对其他常染色体(7号、14号)及X染色体进行了特异性消除[22-23]。实现该技术的关键是在Y染色体上寻找能够被sgRNA特异识别的靶向序列，这些序列在其他常染色体或X染色体上不能被识别，通过对其进行多重切割从而实现整条染色体的消融。因此在选择靶点时，单一数据库的比对结果不能很好地反映个体染色体组的真实情况，最终选择靶基因是在Ensembl比对的基础上重新与NCBI、Sanger测序结果进行比对，选择多个数据库共有序列并与个体基因组序列结合，增加靶点的可靠性。
不同于一般的sgRNA效率检测采用T7E1酶切和测序验证，由于本试验靶点为Y上的多拷贝重复序列，这些片段在Y上广泛分布且保守性不强，可能会存在一些SNP或者突变，导致扩增未编辑前的靶基因产物时出现套峰现象，和编辑后的测序结果无法区分。因此在这里只能通过体外切割进行sgRNA的效率检测。最终以ID为ENSSSCG00000043058、ENSSSCG00000046447的基因为靶点设计sgRNA，通过对Y上不同靶点，即散在重复和簇状聚集的靶基因的打靶，探索了不同打靶位点对编辑效率的影响，结果表明簇状聚集的打靶效果要明显优于散在重复序列，这与Zuo等[22]和Adikusuma等[23]的研究结果一致，其原因可能是长臂偶然断裂造成的染色体截断或易位导致的。另外由于我们选择的基因为非蛋白编码基因，所以在定量分析时最终选择以DNA为模板。在检测靶点DNA含量变化的同时我们也对SRY含量进行了定量，结果均表明其表达量随靶基因降低而降低，虽然这种变化并不明显，但从侧面说明了整条Y染色体降解的可能性。核型分析的结果也表明部分单细胞克隆在染色体水平确实存在着Y染色体的丢失，表明通过多位点断裂造成整条染色体的缺失是可行的。
本研究通过寻找能够特异性识别Y染色体的多拷贝重复序列设计sgRNA，并通过体外切割、qPCR、细胞转染、核型分析等手段证明了该系统在体外体内都有较高的切割效率，可成功实现整条染色体的消除，为未来染色体功能及动物性别控制研究提供了新的方向。
-
![]()
图 1 不同生态类型蚯蚓处理的赤红壤各粒级(d)团聚体的分布情况
LCK：赤红壤不接种蚯蚓；LEf：赤红壤+赤子爱胜蚓；LPc：赤红壤+南美岸蚓；LAr：赤红壤+壮伟远盲蚓；LAa：赤红壤+参状远盲蚓；图中相同粒级柱子上方的不同小写字母表示处理间差异显著(P<0.05，Duncan’ s法)
Figure 1. Distributions of latosolic red soil aggregates with different size (d) processed by different ecological types of earthworm
LCK: Latosolic red soil without earthworm; LEf: Latosolic red soil+E. fetida; LPc: Latosolic red soil+P. corethrurus; LAr: Latosolic red soil+A. robustus; LAa: Latosolic red soil+A. aspergillum. Different lowercase letters on the bars of the same graded column indicated significant differences among different treatments(P<0.05, Duncan’ s test)
![]()
图 2 不同生态类型蚯蚓处理的红壤各粒级(d)团聚体的分布情况
RCK：红壤不接种蚯蚓；REf：红壤+赤子爱胜蚓；RPc：红壤+南美岸蚓；RAr：红壤+壮伟远盲蚓；RAa：红壤+参状远盲蚓. 图中相同粒级柱子上方的不同小写字母表示处理间差异显著(P<0.05，Duncan’ s法)
Figure 2. Distributions of red soil aggregates with different size (d) processed by different ecological types of earthworm
RCK: Red soil without earthworm; REf: Red soil+E. fetida; RPc: Red soil+P. corethrurus; RAr: Red soil+A. robustus; RAa: Red soil+A. aspergillum. Different lowercase letters on the bars of the same graded column indicated significant differences among different treatments(P<0.05, Duncan’ s test)
![]()
图 3 不同生态类型蚯蚓处理的赤红壤各粒级(d)水稳性团聚体分布情况
LCK：赤红壤不接种蚯蚓；LEf：赤红壤+赤子爱胜蚓；LPc：赤红壤+南美岸蚓；LAr：赤红壤+壮伟远盲蚓；LAa：赤红壤+参状远盲蚓；图中相同粒级柱子上方的不同小写字母表示处理间差异显著(P<0.05，Duncan’ s法)
Figure 3. Distributions of water stable aggregates with different size (d) in latosolic red soil processed by different ecological types of earthworm
LCK: Latosolic red soil without earthworm; LEf: Latosolic red soil+E. fetida; LPc: Latosolic red soil+P. corethrurus; LAr: Latosolic red soil+A. robustus; LAa: Latosolic red soil+A. aspergillum. Different lowercase letters on the bars of the same graded column indicated significant differences among different treatments (P<0.05, Duncan’ s test)
![]()
图 4 不同生态类型蚯蚓处理的红壤各粒级(d)水稳性团聚体分布情况
RCK：红壤不接种蚯蚓；REf：红壤+赤子爱胜蚓；RPc：红壤+南美岸蚓；RAr：红壤+壮伟远盲蚓；RAa：红壤+参状远盲蚓. 图中相同粒级柱子上方的不同小写字母表示处理间差异显著(P<0.05，Duncan’ s法)
Figure 4. Distributions of water stable aggregates with different size (d) in red soil processed by different ecological types of earthworm
RCK: Red soil without earthworm; REf: Red soil+E. fetida; RPc: Red soil+P. corethrurus; RAr: Red soil+A. robustus; RAa: Red soil+A. aspergillum. Different lowercase letters on the bars of the same graded column indicated significant differences among different treatments(P<0.05, Duncan’ s test)
![]()
图 5 不同生态类型蚯蚓处理的土壤团聚体破坏率(PAD)
CK：不接种蚯蚓；Ef：接种赤子爱胜蚓；Pc：接种南美岸蚓；Ar：接种壮伟远盲蚓；Aa：接种参状远盲蚓. 图中相同土壤类型柱子上方的不同小写字母表示不同生态类型蚯蚓处理间差异显著(P<0.05，Duncan’ s法)
Figure 5. Percentage of aggregate disruption(PAD) of soil aggregates processed by different ecological types of earthworm
CK: Without earthworm; Ef: Adding E. fetida; Pc: Adding P. corethrurus; Ar: Adding A. robustus; Aa: Adding A. aspergillum. Different lowercase letters on the bars of the same soil type indicated significant differences among different treatments by different ecological types of earthworm(P<0.05, Duncan’ s test)
![]()
图 6 不同生态类型蚯蚓对土壤团聚体分布及稳定性影响的主成分分析
DL：大团聚体含量；DS：小团聚体含量；DM：微团聚体含量；WL：水稳性大团聚体含量；WS：水稳性小团聚体含量；WM：水稳性微团聚体含量；PAD：土壤团聚体破坏率；LCK：赤红壤不接种蚯蚓；LEF：赤红壤+赤子爱胜蚓；LPc：赤红壤+南美岸蚓；LAr：赤红壤+壮伟远盲蚓；LAa：赤红壤+参状远盲蚓；RCK：红壤不接种蚯蚓；REf：红壤+赤子爱胜蚓；RPc：红壤+南美岸蚓；RAr：红壤+壮伟远盲蚓；RAa：红壤+参状远盲蚓
Figure 6. Principal component analyses of the effects of different ecological types of earthworm on the distributions and stabilities of soil aggregates
DL: Large-aggregate content; DS: Small-aggregate content; DM: Micro-aggregate content; WL: Water stable large-aggregate content; WS: Water stable small-aggregate content; WM: Water stable micro-aggregate content; PAD: Percentage of aggregate disruption of soil aggregate; LCK: Latosolic red soil without earthworm; LEf: Latosolic red soil+E. fetida; LPc: Latosolic red soil+P. corethrurus; LAr: Latosolic red soil+A. robustus; LAa: Latosolic red soil+A. aspergillum; RCK: Red soil without earthworms; REf: Red soil+E. fetida; RPc: Red soil+P. corethrurus; RAr: Red soil+A. robustus; RAa: Red soil+A. aspergillum
表 1 供试土壤的基本理化性质
Table 1 Physichemical properties of the soil sample
土壤类型
Soil typepH(H2O) w/(g·kg−1) 碳/氮质量比
Carbon/nitrogen mass ratiow(黏粒1))/%
Clay content阳离子交换量/(cmol·kg−1)
Cation exchange capacity有机碳
Organic carbon全氮
Total nitrogen赤红壤 Latosolic red soil 4.25 27.00 1.85 14.60 22.00 12.3 红壤 Red soil 6.77 21.30 2.92 7.30 32.10 13.3 1) d<2 000 μm 
下载: 导出CSV
-
[1] 赵其国. 中国东部红壤地区土壤退化的时空变化、机理及调控[M]. 北京: 科学出版社, 2002. [2] 张凤荣. 土壤地理学[M]. 北京: 中国农业出版社, 2001: 126-130. [3] 冯瑞章, 周万海, 龙瑞军, 等. 江河源区不同退化程度高寒草地土壤物理、化学及生物学特征研究[J]. 土壤通报, 2010, 41(2): 263-269. [4] 江仁涛, 李富程, 沈凇涛. 川西北高寒草地退化对土壤团聚体组成及稳定性的影响[J]. 水土保持研究, 2018, 25(4): 36-42. [5] 黄福珍. 论蚯蚓对土壤结构形成及性态的影响[J]. 土壤学报, 1979, 16(3): 211-217. [6] 张卫信, 陈迪马, 赵灿灿. 蚯蚓在生态系统中的作用[J]. 生物多样性, 2007, 15(2): 142-153. doi: 10.3321/j.issn:1005-0094.2007.02.005 [7] 王笑, 王帅, 滕明姣, 等. 两种代表性蚯蚓对设施菜地土壤微生物群落结构及理化性质的影响[J]. 生态学报, 2017, 37(15): 5146-5156. [8] KARACA A. Biology of earthworms[M]. London: Chapman & Hall, 2011: 39-50.
[9] 郝桂玉, 黄民生, 徐亚同. 蚯蚓及其在生态环境保护中的应用[J]. 环境科学研究, 2004, 17(3): 75-77. doi: 10.3321/j.issn:1001-6929.2004.03.020 [10] 袁新田, 焦加国, 朱玲, 等. 不同秸秆施用方式下接种蚯蚓对土壤团聚体及其中碳分布的影响[J]. 土壤, 2011, 43(6): 968-974. [11] 于建光, 胡锋, 李辉信, 等. 接种蚯蚓对土壤团聚体分布、稳定性及有机碳赋存的影响[J]. 水土保持学报, 2010, 24(3): 175-179. [12] BOSSUYT H, SIX J, HENDRIX P F. Protection of soil carbon by micro-aggregates within earthworm casts[J]. Soil Biol Biochem, 2005, 37(2): 251-258. doi: 10.1016/j.soilbio.2004.07.035
[13] SCHRADER S, ZHANG H. Earthworm casting: Stabilization or destabilization of soil structure?[J]. Soil Biol Biochem, 1997, 29: 469-475. doi: 10.1016/S0038-0717(96)00103-4
[14] 赵长巍. 接种蚯蚓对黑土水稳性团聚体含量及有机碳蓄积的影响[D]. 长春: 吉林农业大学, 2007. [15] BOSSUYT H, SIX J, HENDRIX P F. Rapid incorporation of carbon from fresh residues into newly formed stable micro-aggregates within earthworm casts[J]. Eur J Soil Sci, 2004, 55(2): 393-399. doi: 10.1111/j.1351-0754.2004.00603.x
[16] JONGMANS A G, PULLEMAN M M, MARINISSEN J C Y. Soil structure and earthworm activity in a marine silt loam under pasture versus arable land[J]. Biol Fertil Soil, 2001, 33(4): 279-285. doi: 10.1007/s003740000318
[17] WINSOME T, MCCOLL J G. Changes in chemistry and aggregation of a California forest soil worked by the earthworm Argilophilus papillifer eisen (megascolecidae)[J]. Soil Biol Biochem, 1998, 30(13): 1677-1687. doi: 10.1016/S0038-0717(97)00262-9
[18] SHIPITALO M J, PROTZ R. Chemistry and micro-morphology of aggregation in earthworm casts[J]. Geoderma, 1989, 45(3): 357-374.
[19] 徐芹, 肖能文. 中国陆栖蚯蚓[M]. 北京: 中国农业出版社, 2011: 281-293. [20] 黄锦华. 蚯蚓对南亚热带森林根系来源碳转化及土壤生态过程的影响[D]. 广州: 中国科学院华南植物园, 2015. [21] ZHANG B G, LI G T, SHEN T S, et al. Changes in microbial biomass C, N, and P and enzyme activities in soil incubated with the earthworms Metaphire guillelmi or Eisenia fetida[J]. Soil Biol Biochem, 2000, 32(14): 2055-2062. doi: 10.1016/S0038-0717(00)00111-5
[22] ZHANG C, MORA P, DAI J, et al. Earthworm and organic amendment effects on microbial activities and metal availability in a contaminated soil from China[J]. Appl Soil Ecol, 2016, 104: 54-66. doi: 10.1016/j.apsoil.2016.03.006
[23] 吕贻忠, 李保国. 土壤学实验[M]. 北京: 中国农业出版社, 2010: 88-92. [24] SIX J, ELLIOTT E T, PAUSTIAN K. Soil macro-aggregate turnover and micro-aggregate formation: a mechanism for C sequestration under no-tillage agriculture[J]. Soil Biol Biochem, 2000, 32(14): 2099-2103. doi: 10.1016/S0038-0717(00)00179-6
[25] 史奕, 陈欣, 沈善敏. 土壤团聚体的稳定机制及人类活动的影响[J]. 应用生态学报, 2002, 13(11): 1491-1494. doi: 10.3321/j.issn:1001-9332.2002.11.032 [26] THIOULOUSE J, CHESSEL D, DOLEDEC S, et al. ADE-4: A multivariate analysis and graphical display software[J]. Stat Comput, 1997(7): 75-83.
[27] SCULLION J, MALIK A. Earthworm activity affecting organic matter, aggregation and microbial activity in soils restored after opencast mining for coal[J]. Soil Biol Biochem, 2000, 32(1): 119-126. doi: 10.1016/S0038-0717(99)00142-X
[28] SPURGEON D J, KEITH A M, Schmidt O, et al. Land-use and land-management change: Relationships with earthworm and fungi communities and soil structural properties[J]. BMC Ecol, 2013, 13(1): 1-13. doi: 10.1186/1472-6785-13-1
[29] BLANCHART E, LAVELLE P, BRAUDEAU E, et al. Regulation of soil structure by geophagous earthworm activities in humid savannas of cote d'ivoire[J]. Soil Biol Biochem, 1997, 29(3/4): 431-439. doi: 10.1016/S0038-0717(96)00042-9
[30] 王清奎, 汪思龙. 土壤团聚体形成与稳定机制及影响因素[J]. 土壤通报, 2005, 36(3): 415-421. doi: 10.3321/j.issn:0564-3945.2005.03.031 [31] PULLEMAN M M, SIX J, BREEMEN N V, et al. Soil organic matter distribution and micro-aggregate characteristics as affected by agricultural management and earthworm activity[J]. Eur J Soil Sci, 2005, 56(4): 453-467. doi: 10.1111/j.1365-2389.2004.00696.x
[32] MAGID J, KJAERGAARD C, GORISSEN A, et al. Drying and re-wetting of a loamy sand soil did not increase the turnover of native organic matter, but retarded the decomposition of added 14C-labelled plant material[J]. Soil Biol Biochem, 1999, 31(4): 595-602. doi: 10.1016/S0038-0717(98)00164-3
[33] 章明奎, 何振立, 陈国潮, 等. 利用方式对红壤水稳定性团聚体形成的影响[J]. 土壤学报, 1997(4): 359-366. doi: 10.3321/j.issn:0564-3929.1997.04.002 [34] 刘艳, 马茂华, 吴胜军, 等. 干湿交替下土壤团聚体稳定性研究进展与展望[J]. 土壤, 2018, 50(5): 4-16. [35] 张池, 周波, 吴家龙, 等. 蚯蚓在我国南方土壤修复中的应用[J]. 生物多样性, 2018, 26(10): 1091-1102. doi: 10.17520/biods.2018151 [36] 张池, 陈旭飞, 周波, 等. 华南地区壮伟环毛蚓(Amynthas robustus)和皮质远盲蚓(Amynthas corticis)对土壤酶活性和微生物学特征的影响[J]. 中国农业科学, 2012, 45(13): 2658-2667. [37] LAVELLE P, SPAIN A V. Soil ecology[M]. Berlin: Springer Netherlands, 2001: 44-57.
-
期刊类型引用(1)
1. 解雅茹，金昊延，孔辰，蔡蓓，张令锴. CRISPR/Cas9系统在家畜生殖细胞中的研究进展. 畜牧兽医学报. 2025(02): 479-491 . 
百度学术
其他类型引用(0)
计量
- 文章访问数: 1083
- HTML全文浏览量: 8
- PDF下载量: 1515
- 被引次数: 1
