Effect of ETV5 knockout on gene expression profile in mouse muscle
-
摘要:目的
与野生型小鼠相比较,采用CRISPR/Cas9技术制备的ETV5基因纯合敲除小鼠在表现出内源性精原干细胞消融的同时伴随着体型和体质的明显弱势,本研究的目的是探究ETV5基因敲除对小鼠肌肉表达谱的影响。
方法采集6周龄的3只野生型公鼠和3只ETV5纯合敲除公鼠的肌肉组织并抽提RNA,进行转录组测序,并对测序结果进行生物信息学分析,对2组小鼠肌肉样品的差异表达基因进行聚类分析、GO和KEGG富集分析。
结果2组小鼠的肌肉组织共筛选出了574个差异表达基因,其中,上调基因292个,下调基因282个。发现了多个基因影响ETV5敲除小鼠的生长发育,其中,包括影响小鼠肌肉发育的基因Amd1和影响脂肪累积的基因Chrna2。GO和KEGG分析富集到的通路大多与脂肪代谢和生长发育相关。
结论本研究结果为ETV5敲除小鼠生长发育缓慢和延迟的分子机制提供了解释,为研究ETV5在生理和病理上的相关功能提供了参考。
Abstract:ObjectiveCompared with the wild type mice, the homozygous ETV5 knockout mice prepared by CRISPR/Cas9 technology showed significant weakness in body size and body weight along with endogenous spermatogonial stem cell ablation. The purpose of this study was to explore the effect of ETV5 knockout on the muscle expression profile of mice.
MethodThe muscle tissues of three wild-type male mice and three ETV5 homozygous knockout male mice aged six weeks were collected and total RNA was extracted for transcriptome sequencing. We analyzed the sequencing results using bioinformatic method. Cluster analysis, GO and KEGG enrichment analyses were performed on the differentially expressed genes in the muscle samples of two groups of mice.
ResultA total of 574 differentially expressed genes were screened out from the muscle tissues of two groups, including 292 up-regulated genes and 282 down-regulated genes. Several genes were found to affect the growth and development of ETV5 knockout mice. These genes included Amd1 affecting muscle development of mice, and Chrna2 affecting fat accumulation. Most of the pathways enriched by GO and KEGG analyses were related to fat metabolism and growth and development.
ConclusionThese results provide an explanation for the molecular mechanism of abnormal development of ETV5 knockout mice, and provide references for further study of in vivo function of ETV5 gene.
-
Keywords:
- ETV5 /
- Transcriptome sequencing /
- Mouse /
- Muscular development /
- Gene knockout /
- Gene expression profile
-
机械式排种器是播种机上使用最多的一种排种器,护种机构是排种器上的重要部件之一,其作业效果直接影响播种质量[1-3]。机械式排种器工作流程一般包括充种、清种、护种和投种4个过程[4-5]。护种机构安装于清种机构后,保护种子停留在型孔中并从清种区到达投种区完成投种[6-7]。护种机构在机械式排种器中被广泛应用,一般有固定板式、固定带式和随动柔性带式[8-13]等。罗锡文等[8,14]研究了一种可拆式弹性随动护种装置,试验结果表明,与固定护种板相比,显著降低了水稻种子的伤种率,种子破损率可降至0.2%以下。王冲等[9,15]研究了一种同步柔性皮带,并分析了同步理论和皮带对水稻种子最大支撑力,导出了同步的临界条件和影响同步效果的因素,与固定板式护种器相比,水稻种子破碎率可降至0.5%。廖庆喜等[16-17]研究的槽孔轮式排种器采用刷种轮与柔性带组合,降低了油菜籽粒的破碎率。张顺等[18-19]设计了一种具有气流清种与护种结构的气力滚筒式水稻精量排种器,提高了播种性能。李兆东等[20-21]设计了一种气流清种与气压护种结合的油菜精量气压式集排器,降低了传统护种方式下油菜种子的破损率。朱德泉等[22]设计了带有护种滑片的型孔轮式水稻精量排种器,播种性能有待进一步提高。
组合型孔排种器的随动柔性带式护种机构的护种带将种子保护在型孔中,护种带和排种轮之间的静摩擦力带动护种带同步转动,到达投种区时释放型孔中的种子,完成投种。柔性护种带与排种轮紧密贴合,防止进入型孔的种子在护种区由于重力和离心力的作用落出型孔而影响播种质量。本文针对排种器实际应用过程中容易出现护种带跑偏、打滑、磨损和断裂等现象,对护种机构进行了结构优化,并对护种带的参数进行了优选。柔性护种机构的结构和护种带的参数是保证护种带与排种轮之间同步的关键,也是保证护种机构可靠性和使用寿命的关键。本文在可拆式弹性随动护种装置的基础上,通过试验分析护种工作过程中引起成穴性下降、伤种率提高、护种带磨损和老化的原因,进一步优化护种机构的结构和参数,旨在提高型孔式水稻精量排种器的播种质量和可靠性,为机械式精密排种器的设计提供参考。
1. 工作原理与结构设计
图1为华南农业大学南方农业机械与装备关键技术教育部重点实验室设计的组合型孔排种器。其随动柔性护种机构的护种带将种子保护在型孔中,护种带和排种轮之间的静摩擦力带动护种带同步转动,到达投种区时释放型孔中的种子,完成投种。柔性护种带与排种轮紧密贴合,防止进入型孔的种子在护种区由于重力和离心力的作用落出型孔而影响播种质量。
图 1 排种器工作原理图1:限种机构;2:第1充种室;3:第2充种室;4:种子;5:排种轮;6:排种管;7:护种带架子;8:轴套;9:护种带;10:清种机构Figure 1. Schematic diagram of seed metering device1: Flow adjusting device; 2: The first filling chamber; 3: The second filling chamber; 4: Seed; 5: Hole-type roller; 6: Seed metering tube; 7: Protecting device; 8: Axle sleeve; 9: Guard belt; 10: Cleaning device1.1 护种机构结构设计
护种机构由护种带架子、轴套、轴和护种带组成(图2)。轴套和轴安装在护种带架子上,呈三角形分布。为了防止护种带在工作时跑偏,同时增加护种带的同步性能,在轴套上设计了3道防滑槽,护种带上设计了3道相应的防滑凸台;安装时,将护种带上的凸台与轴套上的防滑槽配合。
假设护种带和排种轮的接触面无相对运动(静摩擦力),同时护种带产生的压力对种子没有损伤,那么在护种区内伤种率为0,且护种带的磨损也最小。然而,在组合型孔排种器的实际应用过程中,护种带会出现不同步甚至停转的情况。进一步观察分析,护种带停止转动会降低护种带的寿命、影响伤种率和穴径。因此,应对护种机构的同步原理进行分析。
1.2 护种机构同步原理分析
排种器启动前期,排种轮在排种轴的带动下,线速度从0达到设定速度(v),同时护种带在与排种轮之间的摩擦力(f)的驱动下完成与排种轮的线速度同步,如图3所示。
图 3 护种机构同步力学分析M1~M3:对应转轴的摩擦力矩,N·m;F:排种轮作用在护种带上压力的合力,N;f1 :排种轮与护种带之间的最大静摩擦力,N;R:排种轮半径,m;r:转轴半径,m;ω:排种轮角速度,rad/sFigure 3. Mechanics analysis of synchronization for seed protecting deviceM1−M3: Friction torque of corresponding shaft, N·m; F: Resultant force of pressure exerted by the seed metering wheel on the seed guard belt, N; f1: Maximum static friction between seed metering wheel and guard belt, N; R: Radius of seed metering wheel, m; r: Radius of shaft, m; ω: Angular velocity of seed metering wheel, rad/s当护种带处于将要打滑的临界状态时,排种轮与护种带之间的最大静摩擦力(f1)作用在护种带上的力矩应等于3根转轴产生的力矩之和(M),如式(1)所示:
$${f_1}R = M = {M_1} + {M_2} + {M_3},$$ (1) $${f_1} = {F}{\mu _1},$$ (2) 式中, R为排种轮的半径,m;M为3根转轴的总摩擦力矩,N·m;M1、M2和M3分别为3根转轴的摩擦力矩,N·m;F为排种轮作用在护种带上压力的合力,N;µ1为排种轮与护种带的摩擦系数。
根据文献[9]对护种带的受力分析,护种带与排种轮同步瞬间,克服3 根转轴所产生的力矩相等,M1的大小与轴套和转轴的特性有关,即和两者之间的摩擦系数有关。
$${M_1} = {f_2}r = {F_2}{\mu _2}r,$$ (3) $${M_1} = {M_2} = {M_3},$$ (4) 式中,f2为轴套对转轴的摩擦力,N;F2为轴套对转轴压力,N;r为转轴半径,m;µ2为转轴与轴套之间的摩擦系数。
由式(1)~(4)可知,当护种带处于将要打滑的临界状态时,如式(5)所示:
$${F}{\mu _1}R = 3{F_2}{\mu _2}r,$$ (5) 当
${F}{\mu _1}R - 3{F_2}{\mu _2}r \geqslant 0$ 时,护种带能与排种轮实现同步。由于F和F2均与护种带形变产生的力有关,因此,增大µ1或减小µ2将直接影响护种带的同步性能。原轴套采用尼龙加工而成,如图4a所示,成本低且具有自润滑性,与不锈钢轴之间的摩擦系数约为0.4[23]。但是,一般材料在温度逐渐升高后,会加快磨损[24]。尼龙套在磨损之后的磨屑逐渐形成磨粒,在温度逐渐升高后磨屑或磨粒会粘在接触表面[25],使得轴套与转轴之间摩擦系数µ2增大,导致护种带打滑。
为此,对轴套进行了改进设计,采用同轴套和滚针轴承的结构,如图4b所示。滚针轴承与不锈钢轴之间的摩擦系数约为0.1[23]。滚针轴承更耐高温更耐磨。因此,采用铜轴套和滚针轴承结构可显著减小µ2。
2. 试验材料与方法
2.1 材料
试验在华南农业大学工程学院南方农业机械与装备关键技术教育部重点实验室进行。试验选用生产中常用的粳稻品种 ‘秀水134’种子和籼稻品种 ‘培杂泰丰’种子。
试验因素为护种带轴套结构(A)和不同硬度的护种带(B)。试验因素与水平如表1所示。
表 1 试验因素和水平Table 1. Factors and levels in experiments水平
Level轴套结构
Bush structure (A)护种带硬度/HA
Hardness of guard belt (B)1 尼龙轴套 40 2 滚针轴承&铜套结构 45 3 50 4 55 5 60 2.2 方法
2.2.1 护种带打滑率试验
以排种器工作0(刚开始工作)、25、50、75和100 h为测量时间点,用高速摄像机(型号:Photron FASTCAM Viewer)测定护种带的打滑率。在护种带和排种轴的边缘分别做好记号,高速摄像机设置像素为1 024×1 024,帧数为200帧,并确保护种带在高速摄像的取像范围内;排种器的工作转速固定在60 r/min,对应的水稻精量穴直播机大田作业参数:机具前进速度为1.2 m/s,穴距为15 cm(图5)。
通过高速摄像回放分析确定时间,t1和t2分别为记录护种带和排种轴转过3圈后所用的时间;利用SolidWorks分析三维模型,分别计算排种轮周长(l1)和护种带安装后的长度(l2)。打滑率(φ)的计算公式如下:
$$\varphi = \frac{{{l_1} {t_2} - {l_2} {t_1}}}{{{l_1} {t_2}}} \times 100{\text{%}}{\text{。}} $$ (6) 2.2.2 护种带磨损率试验
目前国际上一般采用三维形貌仪对橡胶磨损进行定量分析[26-31]。试验在华南理工大学测试中心进行,采用多功能三维形貌仪(Rtec-instruments UP)分别对使用前和工作100 h后的护种带表面三维形貌(磨损体积)测量。三维形貌仪的工作原理是利用三维共聚焦和白光干涉成像,通过原子力显微镜进行观察成像。将护种带粘贴至方形底座上,并将测试表面压平整。将三维形貌仪获取的图像导入Gwyddion软件中进行图像分析,计算磨损体积。
2.2.3 护种机构对播种性能影响试验
考虑直播机种植一季水稻的工作时长约为200 h,因此以直播机工作中期(即排种器工作100 h)作为参考测量时间点,不同处理的护种机构在多功能播种试验台上进行台架试验。采用人工清选的方式挑选出试验种子中的杂质和已破损种子,试验转速选择60 r/min,投种高度为20 cm,种床带速为1.5 m/s,型孔选择大型孔和小型孔。
伤种率试验参考GB/T6973—2005[32],记录传输带上排种器稳定工作时排出300穴种子中的每穴粒数和破碎种子的粒数(若有半粒种子情况,2个半粒算1粒),重复3次。
2个全新的组合型孔排种器分别安装全新护种机构1[轴套结构为尼龙轴套(A1)]和护种机构2[轴套结构为滚针轴承&铜套结构(A2)],并分别运转10 min,在排种口用容器分别接种。测量每次总播种量,并挑选出每组所排出的破碎种子,记录其数量(若有半粒种子情况,2个半粒算1粒),每组试验重复3次。2个护种机构工作100 h后,再重复上述试验过程。伤种率(K)计算公式如下:
$$K = \frac{{n m}}{{1\;000M}} \times 100{\text{%}} ,$$ (7) 式中:n为破损种子的数量,粒;m为湿种子的千粒质量,g;其中,‘秀水134’湿种子千粒质量为34.2 g,‘培杂泰丰’的湿种子千粒质量为33.4 g;M为试验10 min排种器的总排种量,g。
每一穴种子之间的最大距离作为穴径,用游标卡尺每组连续测量50成穴,重复3次。
3. 结果与分析
3.1 护种带打滑率
由表2的结果可知,试验刚开始时,试验组A2B1的打滑率为0.10%,为该时刻所有试验组中最优,且与其他组合有显著性差异(P<0.05);25和50 h后,试验组A2B3为同一时刻所有组别中最优,且与其他组别有显著性差异;75和100 h后,试验组A2B2、A2B3、A2B4和A2B5的打滑率相对较低,为同一时刻所有组别中最优,且A2B2、A2B3和A2B4之间没有显著性差异,与其他试验组之间均有显著性差异。双因素方差分析结果可知:试验因素 A和B对打滑率影响都极显著(P<0.01),且两者之间存在显著交互作用。A1与A2 之间差异极显著,且A2明显优于A1;B1与B2、B3、B4、B5之间差异极显著,而B2~B5之间差异不显著。
表 2 不同轴套结构不同工作时间对护种带打滑率的影响1)Table 2. Effects of different structure of axle sleeves on slip rate of guard belt at different working time% 试验组
Test group0 h 25 h 50 h 75 h 100 h A1B1 32.40±0.10f 21.28±0.15g 79.82±0.20h 68.81±0.2g 75.30±0.22e A1B2 2.30±0.01e 4.80±0.09e 5.41±0.10f 6.60±0.13e 6.21±0.14c A1B3 1.01±0.05d 3.40±0.11d 4.32±0.09e 5.52±0.14d 6.42±0.15c A1B4 0.97±0.02c 2.82±0.07c 4.27±0.01e 5.12±0.11c 6.35±0.12c A1B5 8.61±0.02c 9.80±0.10f 8.92±0.11g 10.21±0.18f 11.45±0.18d A2B1 0.10±0.01a 0.50±0.04b 1.20±0.03d 1.82±0.08b 1.81±0.08b A2B2 0.27±0.02b 0.45±0.03b 0.53±0.04b 0.62±0.05a 0.62±0.05a A2B3 0.16±0.01b 0.26±0.02a 0.25±0.02a 0.60±0.03a 0.52±0.03a A2B4 0.19±0.08b 0.57±0.03b 0.79±0.06c 0.71±0.02a 0.63±0.023a A2B5 0.25±0.01b 0.46±0.04b 0.43±0.04b 0.64±0.04a 0.64±0.035a 1) A1:尼龙轴套结构;A2:滚针轴承&铜套结构;B1~B5代表护种带硬度分别为40、45、50、55和60 HA。同列数据后不同小写字母表示差异显著(P<0.05,Duncan’s 法)
1) A1: Nylon bushing structure; A2: Needle bearing & copper bushing structure; B 1, B2, B3, B4 and B5 represented the hardness of 40, 45, 50, 55 and 60 HA, respectively. Different lowercase letters of data in the same column indicated significant differences (P<0.05, Duncan’s test)综合分析,试验组A2B3(轴套结构为滚针轴承&铜套、护种带硬度为50 HA)的护种带打滑率最小,效果最佳,且这种优势不随时间的延长而改变。
3.2 表面磨损试验结果
由表3可知,使用前所有护种带的表面形貌没有明显的差异,说明使用前所有护种带的表面较平整,且粗糙度的一致性较好。工作100 h后,比较同一硬度的护种带,采用A2轴套结构的磨损体积小于A1,且A1与A2差异显著。与轴套结构A1相比,A2与轴的摩擦系数较大,阻力也较大,而橡胶磨耗性能(磨耗量)随着温度和负荷的增大而增大[26-27],所以采用A1轴套的护种带磨损体积明显大于A2。
表 3 不同工作时间护种带磨损试验的磨损体积1)Table 3. The wear volume of guard belt in wear test at different working time×10−3 mm3 t/h A1B1 A1B2 A1B3 A1B4 A1B5 A2B1 A2B2 A2B3 A2B4 A2B5 0 2.5±0.2a 2.5±0.2a 2.6±0.3a 2.6±0.3a 2.4±0.3a 2.5±0.3a 2.5±0.2a 2.5±0.3a 2.4±0.2a 2.5±0.2a 100 205.6±14.6e 195.6±15.7e 170.1±9.7d 136.4±12.9c 98.7±8.9b 102.4±5.7b 98.7±4.7b 72.6±6.8a 70.5±5.4a 69.4±5.1a 1) A1:尼龙轴套结构;A2:滚针轴承&铜套结构;B1~B5代表护种带硬度分别为40、45、50、55和60 HA。同行数据的不同小写字母表示差异显著(P<0.05,Duncan’s 法)
1) A1: Nylon bushing structure; A2: Needle bearing & copper bushing structure; B 1, B2, B3, B4 and B5 represented the hardness of 40, 45, 50, 55 and 60 HA, respectively. Different lowercase letters of data in the same line indicated significant differences (P<0.05, Duncan’s test)比较同一轴套结构的护种带,磨损体积均随着的护种带硬度的增大而减小。单因素方差分析可知,试验组A2B3、A2B4、A2B5对护种带磨损体积影响最小,且差异不显著。
3.3 护种机构对播种性能的影响
图6a和6b的单因素方差分析结果表明,试验组A2B2、A2B3对水稻种子伤种率影响最小。由双因素方差分析可知,A和B对伤种率的影响均达极显著水平(P<0.01),且A与B之间存在显著交互作用;A1与A2之间差异极显著,且A2的伤种率明显小于A1,说明A2(滚针轴承&铜套结构)明显优于A1(尼龙轴套结构);试验组B2与B3之间无显著性差异,与其他试验组之间有显著性差异。
图 6 不同孔型护种机构对不同水稻品种伤种率和穴径的影响A1:尼龙轴套结构;A2:滚针轴承&铜套结构;B1~B5代表护种带硬度分别为40、45、50、55和60 HA。各图中,相同护种机构柱子上不同小写字母表示差异显著(P<0.05,Duncan’s 法)Figure 6. Effect of different protecting structure on hill diameter and seed damage rate of different rice varietyA1:Nylon bushing structure;A2:Needle bearing & copper bushing structure;B1, B2, B3, B4 and B5 represented the hardness of 40, 45, 50, 55 and 60 HA, respectively. Different lowercase letters on bars of the same seed protecting device indicated significant differences (P<0.05, Duncan’s test)试验组A2B2、A2B3、A2B4对穴径影响最小(图6c和6d)。由双因素方差分析可知,A对穴径影响显著,B对穴径无显著影响;A1、A2之间差异极显著,且A2明显优于A1。
综上所述,护种带的同步性对护种带的磨损和播种质量有显著影响。护种带的同步性越好(打滑率越小),护种带的磨损越小,伤种率越小,成穴性越好。
4. 结论
相同水稻品种播种时,小型孔的伤种率均高于大型孔的伤种率;在相同试验因素和水平下,‘培杂泰丰’(长粒型籼稻)的伤种率均高于‘秀水134’(圆粒形粳稻)的伤种率。原因可能有以下几个方面:小型孔总播种量小于大型孔,因此,小型孔的伤种率会高于大型孔;小型孔的容积相对大型孔较小,种子充入小型孔中时,种子(尤其是长粒型籼稻‘培杂泰丰’)可能会有部分露在型孔外面,由于护种带与排种轮之间的压力和打滑现象,种子会与排种轮、护种带发生相对移动,造成种子破损;种子在排种轮和护种带之间发生相对移动,会造成与型孔中的种子产生落种时间差,从而影响成穴性。
本研究根据护种带的同步原理,优化了轴承结构和护种带参数,提高了排种器的播种性能,延长了护种带使用寿命。在试验组内,轴套结构为滚针轴承&铜套、护种带硬度为50 HA(A2B3)的护种带磨损体积为72.6×10−3 mm−3,比优化前(A1B1)减少了约1/3,大、小型孔的伤种率分别为0.04%和0.15%,比优化前分别降低了约4/5和2/3;A2B3成穴性最好,播种性能最优。
-
图 9 qPCR检验差异表达基因在ETV5基因敲除和野生型小鼠肌肉组织的相对表达量
柱子上方的“*”“**”“***”分别表示2组差异达到0.05、0.01和0.001的显著水平(t检验)
Figure 9. Detecting the relative expression of DEGs in muscle tissue ofETV5 knockout and wild type mice by qPCR
“*”“**”“***”on bars indicate significant difference between two groups at the levels of 0.05, 0.01 and 0.001, respectively (t test)
表 1 部分差异表达基因
Table 1 Partial DEGs
上调基因
Up-regulated
geneFC1) 下调基因
Down-regulated
geneFC1) Chrna9 7.7151 Acaa2 −0.9966 Cdsn 6.9906 Ccdc85a −1.0183 Adrb3 6.8571 Homer2 −1.0468 Retn 4.8669 Crispld2 −1.0508 Slc7a10 4.2691 Rcan2 −1.0931 Fasn 3.9777 Rhobtb1 −1.1039 Acaa1b 3.8316 Khdrbs3 −1.1346 Cidec 3.5500 Spon1 −1.1809 Katnal2 3.0976 Nudt18 −1.2179 Elovl6 3.0859 Rps6ka1 −1.2612 Spon2 3.0029 Gdap1 −1.2711 Adig 2.8624 Nceh1 −1.3692 Rbp4 2.8034 Map2k6 −1.4616 Eda2r 2.5549 Grb14 −1.5417 Agpat2 2.4480 Dbp −1.7583 Smoc1 2.4085 Plin5 −1.7990 Igf2 2.2913 Hpn −1.9161 Ankrd1 2.1497 Smox −2.0291 Meg3 2.0777 Fabp3 −2.0493 B3galt2 2.0766 Mylk4 −2.0812 Slc1a3 2.0576 Mchr1 −2.1755 Perp 2.0052 Fabp3-ps1 −2.1760 Peg3 1.9896 Cacna2d4 −2.2257 Igfn1 1.9846 Amd2 −2.6730 Dapk1 1.8872 Bdh1 −2.9090 Pde3b 1.7084 Wisp2 −2.9357 Tbc1d31 1.6355 Amd-ps4 −2.9575 Fst 1.3919 Amd1 −3.0079 Thbs1 1.3821 Gapdh −3.7085 Il20rb 1.3807 Chrna2 −4.8778 Cd28 1.3191 Fam19a4 −5.1516 Apoe 0.9978 Etv5 −6.8862 1) FC表示log2 (差异倍数)
1) FC indicates log2 (Fold change) -
[1] LATCHMAN D S. Transcription-factor mutations and disease[J]. The New England Journal of Medicine, 1996, 334(1): 28-33. doi: 10.1056/NEJM199601043340108
[2] SHARROCKS A D. The ETS-domain transcription factor family[J]. Nature Reviews Molecular Cell Biology, 2001, 2(11): 827-837. doi: 10.1038/35099076
[3] OH S, SHIN S, JANKNECHT R. ETV1, 4 and 5: An oncogenic subfamily of ETS transcription factors[J]. Biochimica et Biophysica Acta, 2012, 1826(1): 1-12.
[4] O'BRYAN M K, GREALY A, STAHL P J, et al. Genetic variants in the ETV5 gene in fertile and infertile men with nonobstructive azoospermia associated with Sertoli cell-only syndrome[J]. Fertility and Sterility, 2012, 98(4): 827-835. doi: 10.1016/j.fertnstert.2012.06.013
[5] EO J, HAN K, MURPHY M K, et al. Etv5, an ETS transcription factor, is expressed in granulosa and cumulus cells and serves as a transcriptional regulator of the cyclooxygenase-2[J]. The Journal of Endocrinology, 2008, 198(2): 281-290. doi: 10.1677/JOE-08-0142
[6] BRINSTER R L, AVARBOCK M R. Germline transmission of donor haplotype following spermatogonial transplantation[J]. Proceedings of the National Academy of Sciences of the United States of America, 1994, 91(24): 11303-11307. doi: 10.1073/pnas.91.24.11303
[7] CHEN C, OUYANG W, GRIGURA V, et al. ERM is required for transcriptional control of the spermatogonial stem cell niche[J]. Nature, 2005, 436(7053): 1030-1034. doi: 10.1038/nature03894
[8] EO J, SHIN H, KWON S, et al. Complex ovarian defects lead to infertility in Etv5-/- female mice[J]. Molecular Human Reproduction, 2011, 17(9): 568-576. doi: 10.1093/molehr/gar021
[9] SCHLESSER H N, SIMON L, HOFMANN M C, et al. Effects of ETV5 (ets variant gene 5) on testis and body growth, time course of spermatogonial stem cell loss, and fertility in mice[J]. Biology of Reproduction, 2008, 78(3): 483-489. doi: 10.1095/biolreprod.107.062935
[10] JAMSAI D, CLARK B J, SMITH S J, et al. A missense mutation in the transcription factor ETV5 leads to sterility, increased embryonic and perinatal death, postnatal growth restriction, renal asymmetry and polydactyly in the mouse[J]. PLoS One, 2013, 8(10): e77311. doi: 10.1371/journal.pone.0077311
[11] ZHANG X, ZHAO X, LI G, et al. Establishment of Etv5 gene knockout mice as a recipient model for spermatogonial stem cell transplantation[J]. Biology Open, 2021, 10(1). doi: 10.1242/bio.056804.
[12] LOVE M I, HUBER W, ANDERS S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2[J]. Genome Biology, 2014, 15(12): 550. doi: 10.1186/s13059-014-0550-8
[13] KANEHISA M, GOTO S. KEGG: Kyoto encyclopedia of genes and genomes[J]. Nucleic Acids Research, 2000, 28(1): 27-30. doi: 10.1093/nar/28.1.27
[14] ARAKI S, SAIGA H, MAKABE K W, et al. Expression of AMD1, a gene for a MyoD1-related factor in the ascidian Halocynthia roretzi[J]. Roux’s Archives of Developmental Biology, 1994, 203(6): 320-327. doi: 10.1007/BF00457803
[15] JUN H, MA Y, CHEN Y, et al. Adrenergic-independent signaling via CHRNA2 regulates beige fat activation[J]. Developmental cell, 2020, 54(1): 106-116. doi: 10.1016/j.devcel.2020.05.017
[16] LUO J, XU Q, JIANG B, et al. Selectivity, cell permeability and oral availability studies of novel bromophenol derivative HPN as protein tyrosine phosphatase 1B inhibitor[J]. British Journal of Pharmacology, 2018, 175(1): 140-153. doi: 10.1111/bph.14080
[17] GRUNBERG J R, HOFFMANN J M, HEDJAZIFAR S, et al. Overexpressing the novel autocrine/endocrine adipokine WISP2 induces hyperplasia of the heart, white and brown adipose tissues and prevents insulin resistance[J]. Scientific Reports, 2017, 7: 43515. doi: 10.1038/srep43515
[18] CEREIJO R, GAVALDA-NAVARRO A, CAIRO M, et al. CXCL14, a brown adipokine that mediates brown-fat-to-macrophage communication in thermogenic adaptation[J]. Cell Metabolism, 2018, 28(5): 750-763. doi: 10.1016/j.cmet.2018.07.015
[19] SAKAKIBARA I, YANAGIHARA Y, HIMORI K, et al. Myofiber androgen receptor increases muscle strength mediated by a skeletal muscle splicing variant of Mylk4[J]. iScience, 2021, 24(4): 102303. doi: 10.1016/j.isci.2021.102303
[20] MASON R R, MOKHTAR R, MATZARIS M, et al. PLIN5 deletion remodels intracellular lipid composition and causes insulin resistance in muscle[J]. Molecular Metabolism, 2014, 3(6): 652-663. doi: 10.1016/j.molmet.2014.06.002
[21] BJURSELL M, GERDIN A K, PLOJ K, et al. Melanin-concentrating hormone receptor 1 deficiency increases insulin sensitivity in obese leptin-deficient mice without affecting body weight[J]. Diabetes, 2006, 55(3): 725-733. doi: 10.2337/diabetes.55.03.06.db05-1302
[22] WANG Y, DING Y, LI J. CRISPR-Cas9-mediated gene editing in mouse spermatogonial stem cells[J]. Methods in Molecular Biology, 2017, 1622: 293-305.
[23] GRABACKA M, PIERZCHALSKA M, DEAN M, et al. Regulation of ketone body metabolism and the role of PPARα[J]. International Journal of Molecular Sciences, 2016, 17(12): 2093.
[24] MAGADUM A, ENGEL F B. PPARβ/δ: Linking metabolism to regeneration[J]. International Journal of Molecular Sciences, 2018, 19(7): 2013.
[25] WANG S, DOUGHERTY E J, DANNER R L. PPARγ signaling and emerging opportunities for improved therapeutics[J]. Pharmacological Research, 2016, 111: 76-85. doi: 10.1016/j.phrs.2016.02.028