Estimation of genetic parameters for ammonia nitrogen and freshwater tolerance traits in Penaeus monodon
-
摘要:
该研究以高氨氮和淡水胁迫致死时间为衡量指标,估计了斑节对虾(Penaeus monodon)耐氨氮和淡水应激性状的遗传参数。以斑节对虾“南海1号”和非洲品系为亲本,建立了27个全同胞家系,含5个半同胞家系。利用单性状动物模型和ASReml软件估计斑节对虾幼体阶段的耐氨氮和淡水应激性状的方差组分和遗传参数。通过最佳线性无偏预测法估计所有个体和家系耐氨氮和淡水应激性状的育种值。斑节对虾耐氨氮和耐低盐的遗传力分别为0.11±0.04和0.29±0.08,且统计检验达到显著水平。斑节对虾耐氨氮和淡水应激性状的家系表型值相关系数为0.15,表现为低度线性正相关。斑节对虾耐氨氮和淡水应激性状的家系育种值相关系数为0.57,表现为中度线性正相关。因此在进行斑节对虾耐氨氮性状选育时,耐淡水应激性状也能得到一定的改良。
Abstract:Taking lethal time of high ammonia nitrogen and freshwater tolerance as indices, we estimated the heritability of ammonia nitrogen tolerance and freshwater tolerance of Penaeus monodon. We established 27 full-sib families including five half-sib families of P. monodon based on two breeding populations of P. monodon from Nanhai No.1 and African strain. Besides, we estimated the variance components and genetic parameters for ammonia nitrogen tolerance and freshwater tolerance by using one trait animal model and ASReml software. Then we calculated the estimated breeding values of all animals and families with best linear unbiased prediction. The heritability for ammonia nitrogen tolerance and freshwater tolerance were 0.11±0.04 and 0.29±0.08, respectively (P<0.05). The family phenotypic value correlation coefficients between ammonia nitrogen tolerance and freshwater tolerance was 0.15, showing low linear positive correlation and statistically insignificant difference. The family breeding value correlation coefficients between ammonia nitrogen tolerance and freshwater tolerance was 0.57, showing moderate linear positive correlation and statistically significant difference (P<0.01). Thus, the resistance to freshwater stress can also be improved to some extent with selection of characteristics of ammonia tolerance of P. monodon.
-
微卫星(microsatellite)DNA是由1~6个核苷酸为单位多次串联重复的简单序列,又称简短串联重复序列(short tandem repeat,STR)、简单序列重复(simple sequence repeat,SSR) 或简单序列长度多态性(simple sequence length polymorphism,SSLP),是20世纪80年代末发展起来的一种分子标记。微卫星具有高度的多态性[1]和共显性,在真核生物基因组中随机分布,每代变异超过2%[2]。其分析可以分辨一个碱基对的差异[3],且等位基因的条带易于识别和解释,比其它分子标记带来更多的信息量[4]。因此在种群遗传结构分析、种群遗传多样性检测、遗传图谱的构建及生产性状位点的连锁分析与QTL定位分析中得到了广泛的应用[5-7]。
由于不同种类之间的微卫星引物序列通用性较差。因此必须首先从实验生物基因组中获得微卫星DNA序列,设计引物,筛选多态性微卫星标记。然而筛选微卫星座位的工作繁重,其使用受到一定的限制。随着分子生物学技术的发展,相继产生了一些新的微卫星DNA分离方法。本文对几种微卫星位点分离技术进行介绍并对其进行分析比较,为选择适合的方法提供参考。
1. 微卫星位点的分离方法
1.1 生物信息学方法(Data mining)与近缘种交叉扩增(cross amplification)
获得微卫星最简捷的途径就是通过互连网从公共数据库(如EMBL、GenBank、EST数据库等)查找微卫星DNA,如河豚(Fugu rubripes)[8]、中国明对虾(Fenneropenaeus orientalis)[9]和栉孔扇贝(Chlamys farreri)[10]等。徐鹏等[9]利用生物信息学方法从10 446个中国明对虾ESTs序列中筛选微卫星DNA,共发现微卫星序列229个,占整个ESTs数据的2.19%,其中含双碱基重复序列146个和3碱基重复序列58个,分别占在ESTs微卫星序列总数的63.76%和25.33%,大部分为完美型(perfect)的重复序列。根据筛选的微卫星序列设计19对引物并进行多态性检测,在有扩增产物的16对引物中,首次筛选得到8个中国明对虾微卫星标记。其次,利用近缘种的已知微卫星引物进行交叉扩增(cross amplification)筛选也可获得一定的有效微卫星引物。鲁双庆等[11]研究了鲤(Cyprinus carpio)微卫星引物对远缘种黄鳝(Monopteras albus)的适用性。结果显示,31对鲤微卫星引物中有11对引物能对黄鳝DNA模板扩增出特异性带谱。每对引物扩增的等位基因数3~13个,平均每个位点5.6个,显示了较高的多态性。林凯东等[12]首次将鲤的微卫星引物用于草鱼(Ctenopharyngodon idellus)基因组分析。
1.2 小片段DNA克隆法(small DNA fragments cloning)
此方法的基本原理是构建目标生物基因组小片段DNA文库,通过杂交筛选出含有微卫星序列的阳性克隆。首先用限制性内切酶充分消化DNA,琼脂糖凝胶电泳,选取400~900 bp片段克隆,构建基因组文库。然后将重组克隆转移到杂交膜上,采用32P标记的重复序列探针如(AT)n进行杂交,放射自显影,获得阳性克隆、测序,设计PCR引物并进行扩增筛选(图 1)。
孙效文等[13]利用上述方法从鲤鱼文库中筛选2 000个菌落,获得阳性克隆45个,有22个含有微卫星,其中完美型的占63.6%,非完美型的占22.7%,混合型的占13.7%。陈微等[14]用此方法从牙鲆(Paralichthys olivaceus)文库中筛选阳性克隆45个,共得到20个微卫星序列,其中完美型13个,非完美型5个,混合型2个。PANIEGO等[15]也用此方法分离向日葵(Helianthus annnus)微卫星标记。测序503个微卫星克隆,设计了271对微卫星引物。
利用微卫星重复序列作探针筛选基因组文库是一种耗时费力的方法[16]。获得阳性克隆的比例通常较少(0.04%~12%之间)。对于那些基因组中微卫星DNA含量不是很丰富的种类(比如鸟类和植物),或需要大量的微卫星标记用于研究种群遗传结构[17-18]或构建遗传图谱[19]时,用此方法分离微卫星难以满足需要。
1.3 以RAPD为基础的PCR分离微卫星法(RAPD-based PIMA approach)
为了避免基因组文库的构建和筛选的繁琐性,一些国外学者先后报道用随机扩增多态DNA(RAPD)方法去扩增未知的微卫星序列[20-21]。其中用PCR分离微卫星(PCR isolation of microsatellite approach,PIMA)[22]就是用RAPD引物从目标生物基因组中获得随机扩增片段,这些扩增产物克隆到T-载体上,然后用含有重复序列引物和载体引物筛选阳性克隆、测序(图 1,PIMA支路)。RAPD片段比随机基因组克隆含有较多的微卫星重复序列[23-25]。与传统的方法相比,PIMA省略了DNA酶切、片段大小选择以及接头连接这些步骤,虽然实验操作简单但研究报道较少。
1.4 引物伸长法(Primer extension approach)
基因组DNA酶切后,选一定大小的片段插入到质粒(pBluescript)或噬菌体(M13mp18)载体上,然后转化、洗脱获得含有单链环状DNA(ssDNA)并构建文库。以ssDNA为模板,与含有重复序列的寡聚核苷酸或带有生物素标记的重复序列进行伸长反应,获得富含微卫星位点双链DNA的次级文库。最后Southern杂交或PCR筛选阳性克隆并测序(图 2)。
OSTRANDER等[26]和PAETKAU[27]运用引物伸长法分别对犬齿类和鲍鱼基因组文库进行实验,并分别获得40%~50%和最高达100%阳性克隆。OSTRANDER等构建了含60 000克隆的单链DNA文库,对于基因组频率低于1%的特定重复序列有600座位(含所需的重复序列)在富集文库中表现出来。并且该方法对分离富含2个寡核苷酸重复序列如(AC)n比较有效,对分离3个和4个寡核苷酸重复序列是否有效还不肯定。此方法实验步骤繁琐使其适用性受到一定限制。
1.5 选择性杂交法(selective hybridization method)
选择性杂交法第一步与传统方法一样,也是基因组DNA经过酶切后与接头或载体连接,用连有接头的DNA片段与标记的重复序列探针选择性杂交,捕获含有微卫星序列的DNA片段。最后PCR扩增、克隆。对重组克隆可直接测序、PCR筛选或Southern blot筛选(图 3)。
此方法实验步骤简单,报道较多[28-32]。其中酶切片段与接头或载体连接是非常重要的一步,因为连接插入片段的量不足和多联体的形成都会限制下一步的实验。酶切片段大小的选择可以在酶切后[30]或在连接后[31]进行,在连接后选择片段大小有利于去除非特异性连接。重复序列探针与DNA的杂交既可在尼龙膜上[28, 32]也可在包被有链霉亲和素的磁珠上[30-31]进行。LI等[33-34]采用选择杂交法已成功分离出8个皱纹盘鲍(Haliotis discus hannai)微卫星标记。
1.6 磁珠富集法(magnetic beads enrichment method)
磁珠富集法分离微卫星分子标记是一种简单高效的方法[31, 35], 已经应用于一些植物和动物微卫星分子标记的分离[26, 36-37]。实验生物基因组DNA经酶切,连上接头,PCR扩增,与生物素标记的重复序列探针杂交,磁珠富集,构建PCR富集文库(图 4)。然后再进行筛选即可得到微卫星分子标记。可用不同的酶和接头进行酶切、连接。如果按AFLP的方法进行酶切、连接、扩增,再富集,即为FIASCO法(fast isolation by AFLP of sequences containing repeats)[38]。包括单酶切扩增产物富集[38]和双酶切预扩增产物富集[36]。
孙效文等[13]用此方法筛选鲤鱼微卫星标记。获得微卫星314个,完美型占79.0 %,非完美型占14.3 %,混合型占6.7 %,重复次数超过10的有293个,占93.3 %。孙效文等[38]用磁珠富集法分离草鱼微卫星DNA标记。筛选获得阳性克隆132个,86.36 %含有微卫星序列。高国庆等[36]用FIASCO法从AFLP片段中分离花生(Arachis hypogaeal)微卫星DNA标记,结果回收纯化14个片段,测序后发现都含有简单重复序列。从预扩增的AFLP片段中富集SSR可获得较多的含简单重复序列的片段,而选择性扩增的AFLP片段不经富集直接测序的方法效率较低[39]。此方法在其他种类(比如一些鸟类、鱼类、甲壳类和红珊瑚)的富集率在50%~90%[38]。这些研究表明,用生物素-磁珠富集法克隆微卫星效率高,成本低,所获微卫星质量高。
2. 结语
直接克隆法耗时费力且筛选的效率也比较低。引物伸长法虽有所报道,其实验操作繁琐使其适用性受到限制。磁珠富集法是一种高效而简单快速的分离方法,已经应用于一些植物和动物微卫星分子标记的分离。整个过程可在一星期完成。如果利用篮白斑筛选克隆,从理论上讲,每一个白色菌斑都应当含有微卫星序列。但是操作过程中一些因素会影响筛选微卫星效率。最主要的影响因素是磁珠的平衡及洗液和洗涤温度的严格控制。从国内发表的文献上看,目前从基因组分离微卫星的主要方法是小片段克隆法和磁珠富集法以及稍作改进的一些方法。还有一些从基因组中克隆微卫星的方法,但最有效的还是基于PCR扩增的磁珠富集法。
从整体上看,微卫星分子标记在陆生动植物的分离和应用比较早,水产动物起步较晚,大多还处于分离筛选的初级阶段。在国内,只有鲤[41],对虾[9, 42],栉孔扇贝[10]等少数种类已分离出微卫星分子标记。网上基因库可利用的资源也还非常有限,远不能满足研究和应用的需要。而我国是一个渔业大国,水产动物种类多,遗传差异大,在遗传多样性分析、连锁图谱构建、数量性状基因(quantitative trait loci, QTL)定位分析、分子标记辅助育种等方面需要大量的微卫星标记。因此,选择有效的分离方法可以加速水产动物尤其主要养殖种类的微卫星分子标记的筛选及其应用进程。随着水产养殖生物的基因组测序和大规模cDNA测序,生物信息学方法分离微卫星DNA及其在近缘种的扩增筛选将越来越受到重视。
-
图 1 斑节对虾家系耐氨氮 (a) 和淡水应激 (b) 存活时间的箱线图
盒子包含了中位数、第一四分位数和第三四分位数;×、−、△和○依次表示平均值、最小值、最大值和异常值
Figure 1. Box plots of ammonia nitrogen tolerance (a) and freshwater tolerance (b) survival time for P. monodon families
The boxes contain median, 25th and 75th percentiles; ×, −, △ and ○ represent mean, minimum, maximum and outliers, respectively.
表 1 斑节对虾27个家系耐氨氮耐和淡水应激存活时间的表型参数
Table 1 Phenotypic parameters of ammonia nitrogen tolerance and freshwater tolerance for 27 families of P. monodon
性状
trait平均值
mean标准差
standard deviation变异系数/%
coefficient of variation最大值
maximum最小值
minimum耐氨氮存活/h
survival for ammonia nitrogen tolerance35.66 24.34 68.26 98 2 耐淡水应激存活/min
survival for freshwater tolerance122.56 53 43.63 330 15 表 2 斑节对虾耐氨氮和淡水应激性状的方差组分和遗传力
Table 2 Variance components, heritabilities of ammonia nitrogen tolerance and freshwater tolerance for P. monodon
性状
trait方差 variance 表型方差 (σp 2)
phenotypic variance加性方差 (σα 2)
additive genetic variance残差方差 (σe 2)
residual variance遗传力 (h2)
heritability遗传力显著性
significantly different heritability耐氨氮/h
ammonia nitrogen tolerance588.27 66.30 521.97 0.11±0.04 * (P<0.05) 耐淡水应激/min
freshwater tolerance2 663.50 783.48 1 880.00 0.29±0.08 * (P<0.05) -
[1] 江世贵, 杨丛海, 周发林, 等. 斑节对虾种虾繁育技术[M]. 北京: 海洋出版社, 2013: 1-5. [2] 周凯敏, 江世贵, 黄建华, 等. 斑节对虾Chitinase-2基因的克隆及其在蜕皮和幼体发育过程中的表达分析[J]. 南方水产科学, 2017, 13(4): 59-68. doi: 10.3969/j.issn.2095-0780.2017.04.008 [3] 孙苗苗, 黄建华, 杨其彬, 等. 13个斑节对虾家系的生长及抗氨氮特性比较[J]. 上海海洋大学学报, 2011, 20(4): 510-516. [4] 李永, 杨其彬, 苏天凤, 等. 氨氮对斑节对虾的毒性及免疫指标的影响[J]. 上海海洋大学学报, 2012, 21(3): 358-362. [5] 黄建华, 李永, 杨其彬, 等. 斑节对虾家系氨氮耐受性的比较[J]. 南方水产科学, 2012, 8(6): 37-43. doi: 10.3969/j.issn.2095-0780.2012.06.006 [6] 盛志廉, 陈瑶生. 数量遗传学[M]. 北京: 科学出版社, 2001: 42-74. [7] 姜鹏, 韩林强, 白俊杰, 等. 草鱼生长性状的遗传参数和育种值估计[J]. 中国水产科学, 2018, 25(1): 18-25. [8] 卢薛, 孙际佳, 王海芳, 等. 鳜鱼生长性状遗传参数的估计[J]. 中国水产科学, 2016, 23(6): 1268-1278. [9] 郑静静, 刘建勇, 蒋湘, 等. 日本囊对虾早期生长性状遗传参数估计[J]. 中国水产科学, 2017, 24(4): 710-717. [10] 栾生, 边文冀, 邓伟, 等. 斑点叉尾鮰基础群体生长和存活性状遗传参数估计[J]. 水产学报, 2012, 36(9): 1313-1321. [11] 梁冰冰, 王鸿霞. 文蛤早期发育阶段生长性状遗传参数的估计[J]. 海洋科学, 2016, 40(3): 47-52. doi: 10.11759/hykx20150215001 [12] 郑静静, 刘建勇, 刘加慧, 等. 凡纳滨对虾(Litopenaeus vannamei)不同生长阶段体重的遗传参数和育种值估计[J]. 海洋与湖沼, 2016, 47(5): 1005-1012. doi: 10.11693/hyhz20160400075 [13] 李吉涛, 李健, 刘萍, 等. 脊尾白虾(Exopalaemon carinicauda)体长和体重遗传力的估计[J]. 海洋与湖沼, 2013, 44(4): 968-972. [14] 李浩, 薛宝宝, 李炼星, 等. 缢蛏快长选育系早期生长性状的遗传参数估计[J]. 海洋渔业, 2018, 40(3): 342-349. doi: 10.3969/j.issn.1004-2490.2018.03.010 [15] 韩奋杰, 张伟杰, 秦宇博, 等. 中间球海胆幼体及稚海胆生长性状的遗传参数估计[J]. 大连海洋大学学报, 2017, 32(2): 145-149. [16] 王雪磊, 李琪, 于红, 等. 长牡蛎壳金性状遗传参数评估及与生长性状的关联分析[J]. 水产学报, 2016, 40(12): 1889-1896. [17] HYUMA K, NOBUYUKI I, AKIHIRO K. Estimation of heritability of tolerance to low-salinity by factorial mating system in rainbow trout (Oncorhynchus mykiss)[J]. Suisan Zoshoku, 2001b, 49(2): 253-260.
[18] HYUMA K, NOBUYUKI I, AKIHIRO K. Estimation of heritability of tolerance to low-oxygen water in rainbow trout (Oncorhynchus mykiss)[J]. Suisan Zoshoku, 2002, 50(3): 369-374.
[19] 张嘉晨, 曹伏君, 刘建勇, 等. 凡纳滨对虾(Litopenaeus vannmei)生长和耐低溶氧性状的遗传参数估计和遗传获得评估[J]. 海洋与湖沼, 2016, 47(4): 869-875. doi: 10.11693/hyhz20160400084 [20] 郝登春, 栾生, 曹宝祥, 等. 凡纳滨对虾家系幼虾淡化和养殖阶段存活性状遗传参数估计[J]. 渔业科学进展, 2018, 39(1): 90-97. [21] 王晓清, 王志勇, 何湘蓉. 大黄鱼(Larimichthys crocea)耐环境因子试验及其遗传力的估计[J]. 海洋与湖沼, 2009, 40(6): 781-785. doi: 10.3321/j.issn:0029-814X.2009.06.017 [22] 王正, 高保全, 刘萍, 等. 三疣梭子蟹(Portunus trituberculatus)耐低盐的遗传力估计[J]. 渔业科学进展, 2015, 36(3): 74-78. doi: 10.11758/yykxjz.20150312 [23] 黄付友. “黄海1号”中国对虾生长性状和对高pH、高氨氮抗性遗传力的估计[D]. 青岛: 中国海洋大学, 2007: 33-41. [24] SUN M M, HUANG J H, JIANG S G, et al. Estimates of heritability and genetic correlations for growth-related traits in the tiger prawn Penaeus monodon[J]. Aquacult Res, 2015, 46(6): 1363-1368. doi: 10.1111/are.2015.46.issue-6
[25] MACBETH M, KENWAY M, SALMON M, et al. Heritability of reproductive traits and genetic correlations with growth in the black tiger drawn Penaeus monodon reared in tanks[J]. Aquaculture, 2007, 270(1/2/3/4): 51-56.
[26] KENWAY M, MACBETH M, SALMON M, et al. Heritability and genetic correlations of growth and survival in black tiger prawn Penaeus monodon reared in tanks[J]. Aquaculture, 2006, 259(1/2/3/4): 138-145.
[27] BENZIE J H, KENWAY M, TROTT L. Estimates for the heritability of size in juvenile Penaeus monodon prawns from half-sib matings[J]. Aquaculture, 1997, 152(1/2/3/4): 49-53.
[28] KRISHNA G, GOPIKRISHNA G, GOPAL C A, et al. Genetic parameters for growth and survival in Penaeus monodon cultured in India[J]. Aquaculture, 2011, 318(1/2): 74-78.
[29] 杨其彬, 温为庚, 黄建华, 等. 斑节对虾4个不同群体建立家系的生长及成活[J]. 南方水产, 2010, 6(3): 36-40. doi: 10.3969/j.issn.1673-2227.2010.03.007 [30] MARTINEZ V, NEIRA R, GALL G E. Estimation of genetic parameters from pedigreed populations: lessons from analysis of alevin weight in coho salmon Oncorhynchus kisutch[J]. Aquaculture, 1999, 180(3/4): 223-236.
[31] CASTILLO-JUAREZ H, CASARES J C, VILLELA C C, et al. Heritability for body weight at harvest size in the Pacific white shrimp, Penaeus (Litopenaeus) vannamei, from a multi-environment experiment using univariate and multivariate animal models[J]. Aquaculture, 2007, 273(1): 42-49. doi: 10.1016/j.aquaculture.2007.09.023
[32] 刘宝锁, 张天时, 孔杰, 等. 大菱鲆生长和耐高温性状的遗传参数估计[J]. 水产学报, 2011, 35(11): 1601-1606. [33] ARGUE B J, ARCE S M, LOTZ J M, et al. Selective breeding of Pacific white shrimp (Litopenaeus vannamei) for growth and resistance to Taura Syndrome Virus[J]. Aquaculture, 2002, 204(3/4): 447-460.
[34] GITTERLE T, SALTE R, GJERDE B, et al. Genetic (co)variation in resistance to White Spot Syndrome Virus (WSSV) and harvest weight in Penaeus Litopenaeus vannamei[J]. Aquaculture, 2005, 246(1/2/3/4): 139-149.
[35] 栾生, 孔杰, 张天时, 等. 基于表型值和育种值的中国对虾生长、抗逆性状相关分析[J]. 海洋水产研究, 2008, 29(3): 14-20. [36] CHARO K H, KOMEN H, REZK M A, et al. Heritability estimates and response to selection for growth of Nile tilapia (Oreochromis niloticus) in lowinput earthen ponds[J]. Aquaculture, 2006, 261(2): 479-486. doi: 10.1016/j.aquaculture.2006.07.007
[37] GITTERLE T, RYE M, SALTE R, et al. Genetic (co)variation in harvest body weight and survival in Penaeus (Litopenaeus) vannamei under standard commercial conditions[J]. Aquaculture, 2005, 243(1/2/3/4): 83-92.
[38] 蒋湘, 刘建勇, 赖志服. 九孔鲍 (Haliotis diversicolor supertexta)耐低盐与生长性状的遗传参数评估[J]. 海洋与湖沼, 2014, 45(3): 542-547. doi: 10.11693/hyhz20130400034 [39] 殷宗俊, 张勤. 利用GLMM方法估计家畜阈性状的遗传力[J]. 中国农业大学学报, 2005, 10(6): 33-36. doi: 10.3321/j.issn:1007-4333.2005.06.008 [40] ØDEGÅRDA J, OLESEN I, GJERDE B, et al. Evaluation of statistical models for genetic analysis of challenge test data on ISA resistance in Atlantic salmon (Salmo salar): prediction of progeny survival[J]. Aquaculture, 2007, 266(1/2/3/4): 70-76.
[41] 赵存发, 高佃平, 李金泉, 等. 内蒙古白绒山羊体重性状遗传力的估计[J]. 畜牧与饲料科学, 2010, 20(6): 73-74. doi: 10.3969/j.issn.1672-5190.2010.06.029 [42] 何玉英, 王清印, 谭乐义, 等. 中国对虾生长性状的遗传力和遗传相关估计[J]. 安徽农业科学, 2011, 39(17): 10499-10502. doi: 10.3969/j.issn.1009-4229-B.2011.04.037 -
期刊类型引用(4)
1. 毛嘉烜,汪振华,章守宇,林军,邹沁东,陈怡卉. 大陈岛海域潮下带岩礁生境鱼类群聚特征. 上海海洋大学学报. 2025(03): 604-617 . 百度学术
2. 廖东荣,黎傲雪,陈宁,吕少梁,王学锋. 湛江港皮氏叫姑鱼种群生物学特征和资源开发状态. 南方水产科学. 2024(03): 27-35 . 本站查看
3. Wen Yang,Wenjia Hu,Bin Chen,Hongjian Tan,Shangke Su,Like Ding,Peng Dong,Weiwei Yu,Jianguo Du. Impact of climate change on potential habitat distribution of Sciaenidae in the coastal waters of China. Acta Oceanologica Sinica. 2023(04): 59-71 . 必应学术
4. 杨雯,胡文佳,陈彬,谭红建,董鹏,俞炜炜,杜建国. 基于MaxEnt模型的中国近海主要石首鱼科鱼类潜在适生区. 生态学杂志. 2022(09): 1825-1834 . 百度学术
其他类型引用(6)