收藏切换
Genetic parameters and response to selection for body weight in turbot (Scophthalmus maximus, Linnaeus)
收藏切换
PDF
Ding LYU1, 2, , Weiji WANG2, , Sheng LUAN2, Yulong HU2, Jiantao GUAN2, Zhixiang LI2, Huanhuan WU2, Jie KONG1, 2, *, Shoutang LIU3
Acta Oceanologica Sinica | 2018, 37(6) : 47 - 51
Less
收藏切换
Acta Oceanologica Sinica | 2018, 37(6): 47-51
Marine Biology
Genetic parameters and response to selection for body weight in turbot (Scophthalmus maximus, Linnaeus)
Full
Ding LYU1, 2, , Weiji WANG2, , Sheng LUAN2, Yulong HU2, Jiantao GUAN2, Zhixiang LI2, Huanhuan WU2, Jie KONG1, 2, *, Shoutang LIU3
Affiliations
  • 1 Wuxi Fisheries College, Nanjing Agricultural University, Wuxi 214081, China
  • 2 Key Laboratory for Sustainable Utilization of Marine Fisheries Resources of Ministry of Agriculture, Yellow Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences, Qingdao 266071, China
  • 3 Yellow Sea Fisheries Limited Company, Haiyang 265122, China
Published: 2018-06-25 doi: 10.1007/s13131-018-1150-3
Outline
收藏切换

Genetic parameters and response to selection were estimated for harvest body weight in turbot. The data consisted of 10 952 individuals of 508 full-sib families from three generations (G0, G1, and G2). The heritability estimates for G0, G1, and G2 were 0.11±0.08, 0.18±0.09, and 0.17±0.07, respectively. Over three generations, the heritability estimate was 0.19±0.04. Maternal and common environmental effects were 0.10±0.04, 0.14±0.04, and 0.13±0.03 within each generation and 0.12±0.01 across generations. The selection differential in growth was 18.24 g in G0 and 21.19 g in G1 corresponding to an average of 19.72 g per generation. The genetic gains were also calculated, they were 22.06 g in G1 and 11.93 g in G2, corresponding to 6.36% and 3.52% body weight. The total genetic gain after two generations was 10.10% body weight, which indicated that the selective breeding program for the body weight trait in turbot was successful.

turbot  /  genetic parameter  /  selective breeding  /  body weight  /  Scophthalmus maximus
Ding LYU, Weiji WANG, Sheng LUAN, Yulong HU, Jiantao GUAN, Zhixiang LI, Huanhuan WU, Jie KONG, Shoutang LIU. Genetic parameters and response to selection for body weight in turbot (Scophthalmus maximus, Linnaeus)[J]. Acta Oceanologica Sinica, 2018 , 37 (6) : 47 -51 . DOI: 10.1007/s13131-018-1150-3
Turbot (Scophthalmus maximus, Linnaeus, 1758) is a marine fish distributed on the Atlantic coasts of Europe, including the Baltic Sea, Black Sea, and Mediterranean Sea (Blanquer et al., 1992; Lei and Liu, 1995). With the advantages of fast growth, high cold resistance, rich taste, and good nutrition, turbot is the most widely cultured commercial flatfish globally. Since its introduction into China from Britain as juvenile by the Yellow Sea Fisheries Research Institute in 1992, turbot aquaculture has developed into an important industry with an annual production of 56 300 t in 2013 (National Technology Research and Development Center of Flatfish Culture Industry, 2014). To promote the development of the turbot aquaculture industry and increase turbot production, there is a major effort in genetic improvement programs aimed at body weight traits to breed fast-growing varieties, because the commercial value of a turbot depends primarily on its body weight.
Selective breeding is a basic approach for genetic improvement that can offer the opportunity of continuous genetic gain; any gain in a core breeding group can be multiplied and expressed in millions of progeny (Lind et al., 2012). Selective breeding has been conducted for many aquatic species to improve the body weight trait (O’Flynn et al., 1999; Bentsen et al., 2012; Luan et al., 2012; Sui et al., 2015). In China, turbot breeding began more than a decade ago and two cross varieties, Danfa Turbot (Registration No. GS-02-001-2010) and Duobao No.1 (Registration No. GS-02-001-2014), have been ratified by the National Certification Committee for Aquatic Varieties of China. However, there have been no reports, to date, on genetic gains for any traits in either turbot or the newly selected turbot varieties.
A selective breeding program based on best linear unbiased prediction (BLUP) to improve harvest body weight for cultured turbot carried out by the Yellow Sea Fisheries Research Institute was initiated in 2006. The selective breeding program for turbot has now been underway for over ten years and phenotypic data for harvest body weight has been obtained for two selective generations (G1 and G2). Here, the heritability for harvest body weight was estimated based on within and across generation datasets using the restricted maximum likelihood (REML) method; the response to selection after two generation selection was calculated in the form of selection differential and genetic gains. The results provide a reference point with directive significance for further turbot selection breeding work.
The breeding program was performed at the Genetic Breeding Center of Seawater Flatfish (36°40′24″N, 121°09′00″E), Yellow Sea Fisheries Research Institute, Haiyang City, Shandong Province, China. The progenies were cultured in two farms: one in the Genetic Breeding Center of Seawater Flatfish and the other one (36°54′14.10″N, 121°48′37″E) in Rushan City, Shandong Province, China.
The base turbot generation (G0) was bred by introduced populations from five countries (Denmark, France, Spain, Britain, and Chile), all of which were introduced as about 5 cm long juvenile in 2003 and 2004. Healthy individuals that fed well and exhibited normal morphology were chosen as the brood stock. The brood stock were reared in 5 m×5 m×0.6 m (L×W×H) ponds at a density of 2–3 kg/m2. Three months before breeding, illumination and water temperature in the culture ponds were controlled to stimulate gonad development. Illumination intensity on the water was 200–600 lx and the time increased gradually from 8 h to 10 h per day. Water temperature increased from 8°C to 14°C. The pond water was exchanged by continuous running water.
The base generation (G0) was established in 2006 and 2007 using hybridization between populations and artificial insemination. Lively individuals, which had no trauma and body weight of more than 1 kg, could be selected as parent fish. Full- and half-sib families were produced via an unbalanced nested mating design by randomly mating selected males and females from different populations. The numbers of turbot from different populations used to produce the base generation (G0) are show in Table 1, and the numbers of family in each hybridized combination were shown in Table 2.
Fertilized eggs were incubated in mini flow water at 13.5–15°C until hatching (~3–5 d). All of the families were then reared separately in fiberglass-reinforced plastic tanks (0.5 m3) at a density of thirty thousand larvae (~50 mL fertilized eggs). The temperature of the larvae culture was maintained at 16–20°C. Salinity was 25–30 and dissolved oxygen was maintained at 4.9–8.0 mg/L. The rearing procedure was kept as consistent as possible for the different families. At the 35th and 70th days after hatching, 1 000 larvae and 400 juvenile fish per family were selected randomly and reared in new fiberglass-reinforced plastic tanks (0.5 m3). At the 3rd months after hatching, 80 individuals with a greater body weight than the mean from each family were selected randomly, weighed, and marked with a VIE tag. The fish were then randomly distributed into either one or several 5 m×5 m×0.6 m (L×W× H) test ponds in two farms. The VIE tag was updated once every three months. Approximately 15 months post-hatching, all of the surviving individuals were landed, toweled dry and measured harvest body weight using electronic scale in turn. The progenies were genetically evaluated for the harvest body weight trait.
The subsequent generations (G1 and G2) were established by family and within-family selection strategy. From G0 to G2, all the full-sib families were ranked according their average family breeding values and about 40 families with high average family breeding values were selected as the parents of the next generation. Selected brood individuals were injected with passive integrated transponder (PIT) markers to distinguish them. Full- and half-sib families in G1 and G2 were produced via an unbalanced nested mating design by mating selected males and females from different families. Families with high average family breeding values were mated first and the candidate parents with higher breeding values were mated first. Brood individuals was mated randomly and inbreeding coefficient was maintained at <1%. Rearing and tagging of families in G1 and G2 were the same as G0.
Because we used PIT markers and VIE tags, a complete pedigree containing all phenotypically measured individuals was available and used in this study. The mixed model should include all effects that potentially played a role on the harvest body weight. The fixed effects considered in this study were year and farm. Each effect contains different factors that impact growth. The year effect contained annual breeding environment factors (e.g., water temperature and quality) and breeding operation. The farm effect contained raising environment, management level, stocking density, etc. A factor at different levels may have different effects on another factor, so their interaction was also included. The animal model was used to estimate the variance component using REML method with ASREML software (Gilmour et al., 2009) as follows:
${y_{ijkm}} = \mu + {a_i} + {F_j} + {Y_k} + {F_j}\times {Y_k} + {A_i}\left( {{F_j}\times {Y_k}} \right) + {d_m} + {e_{ijkm}},$
where yijkm is the phenotypic observation for harvest body weight, μ is the overall mean, ai is the random additive genetic effect of animal i, Fj is the fixed effect of farm j (two levels), Yk is the fixed effect of year (seven levels), Fj×Yk is the interaction between farm and year, Ai (Fj×Yk) is a linear covariate nested within the interaction between year and farm, dm is the random maternal and common environmental effects on full-sib m, and eijkm is the random error term.
The distribution of the random effects a, d, and e were assumed to be normal, with means of zero. The variance-covariance matrix was as follows:
$\text{V}\left[ {\begin{array}{*{20}{c}}a\\d\\e\end{array}} \right] = \left[ {\begin{array}{*{20}{c}}{\text{A}\sigma _a^2} & 0 & 0\\0 & {{\text{I}_d}\sigma _d^2} & 0\\0 & 0 & {{\text{I}_e}\sigma _e^2}\end{array}} \right],$
where ${\sigma _a^2}$, ${\sigma _d^2}$, and ${\sigma _e^2}$ are the variance of the random effects a, d, and e, respectively. A is the n×n matrix of additive genetic relationship, and Id and Ie are n×1 identity matrices. n is the number of individuals.
Total phenotypic variance (${\sigma _p^2}$) was calculated as the sum of the additive genetic variance (${\sigma _a^2}$), the maternal and common environmental effects variance (${\sigma _d^2}$), and the error variance (${\sigma _e^2}$). The heritability (h2) was computed as $\sigma _a^2/\sigma _p^2$, and maternal and common environment effects (c2) as $\sigma _d^2/\sigma _p^2$. All of the calculations were carried out in ASREML software (Gilmour et al., 2009).
The z-score was used to test whether the heritability significantly varied from zero:
$z = \frac{{{h^2}}}{{se}},$
where h2 is the estimates of heritability, and se is the standard errors of heritability.
The selection differential for growth was estimated for each generation by comparing the mean estimated breeding value of the brooders that actually produced progeny and the mean of all of the fish before the brooders were selected. Estimated breeding values were calculated based on phenotype information available at the time of selection. The genetic gain in each generation was also calculated by comparing the mean estimated breeding value of the current and previous generation. The cumulative genetic gain after two selections was expressed as a percentage based on the following formula:
${P_c} = \mathop \prod \limits_i^n \left( {1 + {P_i}} \right) - 1,$
where Pc is the cumulative genetic gain (%), Pi is the genetic gain (%) for the ith generation, and i is the number of the selection generation (i=1, 2).
The main statistical results of data set are shown in Table 3 and Table 4. The least square means of harvest body weight in Rushan farm were larger than Haiyang in 2009 and 2011 and progenies in the other years were only cultured in one farm. In total, 508 full-sib turbot families were constructed from 2006 to 2014. There were no and very few mature female fish in 2008 and 2012, respectively, therefore, no families were produced in either year. Because the selected parent individuals develop at different rates, each generation contained families in either two or three years. The number of full-sib families for each year ranged from 54 to 97. Because some parent individuals were mated with more than one partner, there were also some maternal half-sib families and paternal half-sib families in each year. The data set consisted of 10 952 individuals in total. The average number of individuals in each full-sib family ranged from 12.63 to 31.14 in different years.
The variance components, heritability, and maternal and common environmental effects are listed in Table 5. Four heritability (h2) estimators for the harvest body weight were obtained. Heritability within G0 (0.11±0.08) was lower than that within G1, G2, and across generations. The other estimates were very similar, 0.18±0.09 within G1, 0.17±0.07 within G2, and 0.19±0.04 across generations. The heritability estimate for body weight in G0 was not significantly different from zero and in G1 and G2 were significantly different from zero. The heritability estimate across generations was very significantly different from zero. Maternal and common environmental effects (c2) within and across generations closely corresponded to heritability. The heritability and maternal and common environmental effects standard errors across generations were smaller than those within generations.
The average estimated breeding value of all individuals and brooders in each generation are listed in Table 6. The growth selection differential was 18.24 g in G0 and 21.19 g in G1 corresponding to an average of 19.72 g per generation. The genetic gains were also obtained, they were 22.06 g in G1 and 11.93 g in G2, corresponding to 6.36% and 3.52% body weight. The cumulative genetic gain after two generations was 10.10% body weight.
In this study, sex was not included in the mixed effect model analysis, because it is difficult to determine turbot gender morphologically at 15 months old and, moreover, dissecting to observe the gonad inevitably results in the death of the individual. Furthermore, Lyu et al. (2016) reported that variance analysis revealed no significant fixed effects of sex on harvest body weight. Wang et al. (2014) also obtained similar results, although female turbot had a growth advantage compared with males the difference was not significant in the early stage.
Heritability estimates for body weight within and across generations were all low to medium in 15 month-old turbot, based on the following categorization: low (0.05–0.15), medium (0.20–0.40), high (0.45–0.60), and very high (>0.65) (Cardellino and Rovira, 1987). There is a limited number of publications regarding the genetic parameters for body weight traits in turbot. Most previous estimates for body weight of adult turbot were medium to high (Gjerde et al., 1997; Ma et al., 2008; Lyu et al., 2016). Although genetic parameter estimates are only applicable to specific populations and environments, different populations, environmental conditions, and statistical models all result in different estimators (Gall et al., 1993); the heritability estimated in this study was lower than previous reports, which suggests that our heritability values may be underestimates. The most likely reason is that the sum of fixed environment effects was underestimated and phenotypic variance was overestimated because some possible subsistent fixed effect of pond was ignored. Furthermore, selecting individuals 3 months post-hatching may have resulted in a slight bias in the heritability estimate, because it might increase resemblance in family members. The heritability values for harvest body weight remained stable across generations, which suggested the potential for continued genetic improvement through selective turbot breeding in China.
A set of relatively large maternal and common environmental effects was estimated in this study. They were 0.10±0.04, 0.14± 0.04, and 0.13±0.03 in G0, G1, and G2, respectively, and 0.12±0.01 across generations, which were similar to the heritability estimate values. Several studies have shown that maternal effects are high during early growth but have no obvious effect on adult fish (Liu et al., 2011; Furutsuka-Uozumi and Tabata, 1999; Shimada et al., 2007); therefore, these relatively large maternal and common environmental effects were likely the result of rearing full-sib families separately before they had reached a suitable size for tagging. Although the number of fish in each fiberglass-reinforced plastic tank was homogenized twice, at 35 and 70 days after hatching, the large variation in death rate in different families was associated with a large variation in breeding density, which resulted in a pronounced environmental effect. These results suggested better synchronized breeding and standardized family production is required.
Notably, the standard error of heritability across generations were much smaller than those within generations and the heritability estimate for body weight in G0 was not significantly different from zero. These indicate that a larger dataset (larger numbers of individuals, families, as well as generations) is necessary to accurately estimate genetic parameters.
In this study, we found a substantial selection response for growth. The selection differential calculated was 18.24 and 21.19 in G0 and G1, respectively. As a consequence of selection differential, the genetic gains were 22.06 and 11.93, corresponding to 6.36% and 3.52% body weight in subsequent generations. The small difference between selection differential and genetic gain may have been because the parents selected produced different amounts of progeny. Moreover random genetic drift was also a possible reason (Falconer and Mackay, 1996). The least squares mean of harvest body weight of G1 was 347.07 g, much larger than that of G0 (299.25 g). The above results indicate that the selection of G0 was successful. However, the actual phenotypic changes in G2 seemed incompatible with average breeding value changes. The most likely explanation was that the environment in G1 was far better than in G2. The changes in the phenotypic value might have been the result of a combined effect of selection differential and environmental changes. Therefore, the actual phenotypic change should only correspond to the expected response to selection when the environmental effects on the parental and progeny generations are identical (Gall et al., 1993).
After two generation selection, both average estimated breeding and phenotypic values were improved, which indicated that the BLUP selective breeding program for the body weight trait in turbot was successful. Nonetheless, the average genetic gains in this study were relatively small (<10% in one generation), compared with earlier reports on other aquatic species (Charo-Karisa et al., 2006; Liu et al., 2014; Luan et al., 2012; Maluwa and Gjerde, 2007; Rezk et al., 2009; Thodesen et al., 2012). Gjedrem (2000) also reported that genetic gain in growth per generation for aquatic animal species ranged from 10% to 20%. The main reason for low genetic gain in this study is that the proportion of turbot with well-developed gonads was low, which resulted in a low selection intensity.
In this study, the number of full-sib families for each year had a large variation, ranging from 54 to 97, and also there were significant difference among least square means of harvest body weight by each year. These results suggested aquaculture environment and management situations were far from the best condition, especially in gonad development process. Number of cultivated families was restricted by a low proportion of turbot with well-developed gonads (especially in female). A sufficient number of full-sib families form the basis of increasing selection intensity and avoiding inbreeding, which results in phenotypic depression (Thodesen et al., 2005).
  • The Taishan Scholar Program for Seed Industry under contract No. ZR2014CQ001; the Accurate Identification and Selection Breeding Creative Utilization of Turbot Germplasm Resources under contract No. 2016LZGC031-2.
Bentsen H B, Gjerde B, Nguyen N H, et al. 2012. Genetic improvement of farmed tilapias: genetic parameters for body weight at harvest in Nile tilapia (Oreochromis niloticus) during five generations of testing in multiple environments. Aquaculture, 338–341: 56–65
Blanquer A, Alayse J P, Berrada-Rkhami O, et al. 1992. Allozyme variation in turbot (Psetta maxima) and brill (Scophthalmus rhombus) (Osteichthyes, Pleuronectoformes, Scophthalmidae) throughout their range in europe. Journal of Fish Biology, 41(5): 725–736
Cardellino R, Rovira J. 1987. Mejoramiento Genético Animal. Buenos Aires, Argentina: Hemisferio Sur
Charo-Karisa H, Komen H, Rezk M A, et al. 2006. Heritability estimates and response to selection for growth of nile tilapia (Oreochromis niloticus) in low-input earthen ponds. Aquaculture, 261(2): 479–486
Falconer D S, Mackay T F C. 1996. Introduction to Quantitative Genetics. 4th ed. Longman Group Limited, Harlow, Essex, UK: Pearson
Furutsuka-Uozumi K, Tabata K. 1999. Estimation of the heritability of quantitative traits in fry and juvenile hirame, Paralichthys olivaceus. Aquaculture Science, 47(1): 49–54
Gall G A E, Bakar Y B, Famula T. 1993. Estimating genetic change from selection. Aquaculture, 111(1–4): 75–88
Gilmour A R, Gogel B J, Cullis B R, et al. 2009. ASReml User Guide Release 3.0. Hemel Hempstead, UK: VSN International Ltd
Gjerde B, Roer J E, Lein I, et al. 1997. Heritability for body weight in farmed turbot. Aquaculture International, 5(2): 175–178
Gjedrem T. 2000. Genetic improvement of cold-water fish species. Aquaculture Research, 31(1): 25–33
Lei Qilin, Liu Xinfu. 1995. An primary study on culture of turbot, Scophthalmus maeoticus L. Modern Fisheries Information (in Chinese), 10(11): 1–3
Lind C E, Ponzoni R W, Nguyen N H, et al. 2012. Selective breeding in fish and conservation of genetic resources for aquaculture. Reproduction in Domestic Animals, 47(S4): 255–263
Liu Baosuo, Zhang Tianshi, Kong Jie, et al. 2011. Estimation of genetic parameters for growth and upper thermal tolerance traits in turbot Scophthalmus maximus. Journal of Fisheries of China, 35(11): 1601–1606
Liu Jianyong, Lai Zhifu, Fu Xueli, et al. 2014. Genetic parameters and selection responses for growth and survival of the small abalone Haliotis diversicolor after four generations of successive selection. Aquaculture, 436: 58–64
Luan Sheng, Yang Guoliang, Wang Junyi, et al. 2012. Genetic parameters and response to selection for harvest body weight of the giant freshwater prawn Macrobrachium rosenbergii. Aquaculture, 362–363: 88–96
Lyu Ding, Wang Weiji, Luan Sheng, et al. 2016. Estimating genetic parameters for growth traits with molecular relatedness in turbot (Scophthalmus maximus, Linnaeus). Aquaculture, 468: 149–155
Ma Aijun, Wang Xinan, Yang Zhi, et al. 2008. The growth traits and their heritability of young turbot (Scophthalmus maximus L.). Oceanologia et Limnologia Sinica (in Chinese), 39(5): 499–504
Maluwa A O, Gjerde B. 2007. Response to selection for harvest body weight of Oreochromis shiranus. Aquaculture, 273(1): 33–41
National Technology Research and Development Center of Flatfish Culture Industry. 2014. Annual Report 2013 of National Technology System for Flatfish Culture Industry (in Chinese). Qingdao: China Ocean University Press
O’Flynn F M, Bailey J K, Friars G W. 1999. Responses to two generations of index selection in Atlantic salmon (Salmo salar). Aquaculture, 173(1–4): 143–147
Rezk M A, Ponzoni R W, Khaw H L, et al. 2009. Selective breeding for increased body weight in a synthetic breed of Egyptian Nile tilapia, Oreochromis niloticus: Response to selection and genetic parameters. Aquaculture, 293(3–4): 187–194
Shimada Y, Shikano T, Murakami N, et al. 2007. Maternal and genetic effects on individual variation during early development in Japanese flounder Paralichthys olivaceus. Fisheries Science, 73(2): 244–249
Sui Juan, Luan Sheng, Luo Kun, et al. 2015. Genetic parameters and response to selection of harvest body weight of the Chinese shrimp Fenneropenaeus chinensis after five generations of multi-trait selection. Aquaculture, 542: 134–141
Thodesen J, Hu Honglang, Kong Jie. 2005. Inbreeding and its impact on aquaculture. Journal of Fisheries of China (in Chinese), 29(6): 849–856
Thodesen J, Rye M, Wang Yuxiang, et al. 2012. Genetic improvement of tilapias in China: Genetic parameters and selection responses in growth of Nile tilapia (Oreochromis niloticus) after six generations of multi-trait selection for growth and fillet yield. Aquaculture, 322–323: 51–64
Wang Xinan, Ma Aijun, Huang Zhihui, et al. 2014. Developmental differences between female and male groups in turbot (Scophthalmus maximus) breeding families. Journal of Fisheries of China (in Chinese), 38(4): 465–470
Year 2018 volume 37 Issue 6
PDF
35
19
Cite this Article
BibTeX
Article Info
doi: 10.1007/s13131-018-1150-3
  • Receive Date:2017-05-17
  • Online Date:2026-04-14
  • Published:2018-06-25
Article Data
Affiliations
History
  • Received:2017-05-17
  • Accepted:2017-09-05
Funding
The Taishan Scholar Program for Seed Industry under contract No. ZR2014CQ001; the Accurate Identification and Selection Breeding Creative Utilization of Turbot Germplasm Resources under contract No. 2016LZGC031-2.
Affiliations
    1 Wuxi Fisheries College, Nanjing Agricultural University, Wuxi 214081, China
    2 Key Laboratory for Sustainable Utilization of Marine Fisheries Resources of Ministry of Agriculture, Yellow Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences, Qingdao 266071, China
    3 Yellow Sea Fisheries Limited Company, Haiyang 265122, China

Corresponding:

References
Share
https://castjournals.cast.org.cn/joweb/aos/EN/10.1007/s13131-018-1150-3
Share to
QR

Scan QR to access full text

Cite this article
BibTeX
Citations
表12种不同金属材料的力学参数

Family
属数
Number of
genus
种数
Number of
species
占总种数比例
Percentage of
total species (%)

Genus
种数
Number of
species
占总种数比例
Percentage of total
species (%)
鹅膏菌科Amanitaceae 2 11 5.26 鹅膏菌属 Amanita 10 4.78
小菇科 Mycenaceae 2 12 5.74 丝盖伞属 Inocybe 5 2.39
多孔菌科 Polyporaceae 8 14 6.70 蜡蘑属 Laccaria 5 2.39
红菇科 Russulaceae 3 23 11.00 小皮伞属 Marasmius 6 2.87
小菇属 Mycena 11 5.26
光柄菇属 Pluteus 5 2.39
红菇属 Russula 17 8.13
栓菌属 Trametes 5 2.39
关闭全屏
  • BibTeX
  • EndNote
  • RefWorks
  • TxT