Tao Feng
Hongxiang Ding
Jing Wang
Wei Xu
Yan Liu
Ákos Kenéz- 1Institute of Animal Husbandry and Veterinary Medicine (IAHVM), Beijing Academy of Agriculture and Forestry Sciences (BAAFS), Beijing, China
- 2Joint Laboratory of Animal Science Between IAHVM of BAAFS and Division of Agricultural Science and Natural Resource of Oklahoma State University, Beijing, China
- 3College of Animal Science and Technology, Henan University of Science and Technology, Luoyang, China
- 4College of Animal Science and Technology, Hebei North University, Zhangjiakou, China
- 5Department of Infectious Diseases and Public Health, City University of Hong Kong, Hong Kong, China
Introduction
Sheep industry is a major branch of animal husbandry throughout the north and central parts of China mainly for mutton and wool production (1, 2). Body weight gain during the fattening period is an important determinant for carcass weight. Sheep average daily gain (ADG) refers to the average weight gain of sheep during a certain period and is an important economic growth trait that improves production efficiency and economic benefits. Specifically, ADG has been reported to have a positive correlation coefficient (r = 0.53) with final body weight in growing lambs (3). The ADG of lamb or sheep is affected by genetic basis, nutrition level, growth stage, and management system (4–6). In sheep feedlots, lambs with the same genetic basis, same age, same management system, and very similar initial body weight often develop a large standard deviation in ADG, accounting for unwanted differences in final body weights (7, 8). To our knowledge, few studies reported the underlying metabolic mechanism of such an inter-individual difference.
Selecting those lambs expected to have high ADG, at the earliest as possible in their longevity, or at least before the fattening period, could obviously increase the profitability of sheep feedlots (9, 10). Recently, several serum components (hormones, metabolites, hematological, and biochemical parameters) have been identified as biomarkers to evaluate the residual feed intake of sheep (11, 12), an indicator of feed conversion efficiency. Metabolomics can qualitatively and quantitatively analyze hundreds of metabolites in diverse samples, which can be extensively used to study physiological and pathophysiological process such as starving and intrauterine growth restriction in sheep (13, 14). Blood is considered as an ideal sample in sheep metabolomics research and is potentially used to reflect the metabolic status on a whole body level (15). As we have known, alterations in the blood metabolome profiles of sheep during the fattening period, particularly in fattening sheep with the same genetic basis, same age, same management system, but high or low ADG, are still unclear. Therefore, the aims of the present study were (i) to reveal the metabolic characteristics of lambs with high or low ADG under the same management system and (ii) to investigate the potential metabolic pathways related to the growth performance of sheep.
Methods
Experiments were performed at the Experimental Station of Beijing Academy of Agriculture and Forestry Sciences in the Yangyuan county, Zhangjiakou city, Heibei province, northeast of China. A total of 200 crossbreed male lambs (Ovis aries) of Dorper rams and Mongolia ewes after weaning (45 days of age) were housed in eight sheltered outdoor paddocks and fed total mixed ration (TMR). Clean water and mineralized salt licks were available ad libitum. From 75 days of age, 50 lambs with similar body weight were selected and reared in individual pens indoors (0.7 × 1.0 m) until 120 days of age. Briefly, the lambs were acclimatized, lasting 15 days before formal assessment. At 90 days, 40 lambs with similar body weight were selected to do ADG research lasting 30 days. Lamb TMR was compounded based on the recommendations of sheep feeding standard in China (NY/T816-2004) and contained digestible energy of 11.83 MJ·kg−1, metabolic energy of 9.73 MJ·kg−1, 14.61% crude protein, 0.39% calcium, and 0.25% phosphorus. Body weight of lambs was accurately measured in the morning before feeding, and at 75, 90, and 120 days of age, using calibrated electronic scales. ADG was calculated based on body weight. Differences in ADG between the high ADG and low ADG lamb group were analyzed using a t-test. A P < 0.05 indicated statistical significance.
At 120 days of age, after weighting, blood samples were drawn from the jugular vein of the top seven lambs with the greatest ADG and the bottom seven lambs with the lowest ADG using needles and vacutainers covered with anti-coagulant (BD Vacutainer, USA) for a minimum of 6 ml. The blood was placed at room temperature for 4 h and then centrifuged at 2,000 g for 30 min at 4°C. Serum separation was carefully proceeded. The serum was aliquoted and rapidly frozen by dry ice. Frozen serum samples were stored at −80°C until metabolomics analyses.
Using 100 μl of serum, metabolites were extracted using methanol. Extracts were sonicated, and after centrifugation, the supernatant was gently added to sample vials for LC-MS/MS analysis. A pooled quality control sample (QC) was performed for system conditioning and quality control. Chromatographic separation of the metabolites was operated on a Thermo UHPLC system equipped with an ACQUITY UPLC HSS T3 (100 mm × 2.1 mm i.d., 1.8 μm; Waters, Milford, USA).
Following LC-MS/MS analyses, the raw data were inputted into the Progenesis QI 2.3 (Non-linear Dynamics, Waters, USA) for peak picking and alignment. Mass spectra of these metabolic characteristics were discerned through the accurate mass, MS/MS fragments spectra, and isotope ratio difference, by scanning in public available biochemical databases such as the Human Metabolome Database (HMDB) (http://www.hmdb.ca/) and the Metlin database (https://metlin.scripps.edu/). A multivariate statistical analysis was conducted using “ropls” (Version 1.6.2, http://bioconductor.org/packages/release/bioc/html/ropls.html) R package from Bioconductor on Majorbio Cloud Platform (https://cloud.majorbio.com). Principal component analysis (PCA) was applied to check outliers and present trends. Partial least squares discriminate analysis (PLS-DA) was used to identify the general metabolic changes in serum of sheep between high and low ADG. Variable importance in the projection (VIP) was computed by an orthogonal partial least squares discriminate analysis (OPLS-DA) model. Differential metabolites among ESI groups were summarized and annotated into their biochemical pathways through metabolic enrichment and pathway analysis based on database matching (KEGG, http://www.genome.jp/kegg/). Furthermore, Volcano plot was used to compare the size of the fold change to statistical significance.
Results
At the beginning of the lamb fattening trial, when lambs were 90 days of age, the mean (SD in parentheses) body weight of seven original lambs corresponding to the high ADG was 26.7 (0.9) kg, while the mean body weight of seven original lambs corresponding to the low ADG was 26.8 (0.8) kg. At the end of the lamb trial, the body weight of the seven lambs with the highest weight gain was 35.67 (0.8) kg and the body weight of the seven lambs with the lower weight gain was 32.56 (1.0) kg. The ADG of high weight gain lambs was 298.1 (15.9) g·day−1, which differed (P < 0.01) from the ADG of low weight gain lambs 191.9 (23.6) g·day−1, while the average ADG of all lambs tested was 239.5 (34.8) g·day−1 (n = 40).
Variation of ADG depends on sheep breeds and ages. Previous studies showed that crossbreed of specialized mutton breeds and local sheep breeds had greater ADG than local sheep. For local sheep breeds under a barn feeding fattening system, 6-month Altay and Hu lambs presented ADG from 100 to 200 g·day−1 (16), while ADG of 3-month Ningxia Tan sheep lambs was between 90 and 130 g·day−1 (17) and ADG of Small Tail Han sheep lambs was between 140 and 180 g·day−1 (18). Furthermore, the ADG of crossbreed lambs of Dorper and Small Tailed Han sheep was between 265 and 322 g·day−1 in pens (19). In our study, crossbreed lambs of Dorper rams and Mongolia ewes at 4-month age exhibited excellent growth performance with an average ADG of 240 g·day−1. Crossbreed lambs are recommended to produce lamb meat in the north part of China.
In the LC-MS spectra of lamb serum with high or low ADG, 10,231 metabolites were initially found. After quality control and discernment, 462 compounds were reliably detected. The PCA score plot presented that the first and second principal components (PCs) clarified 22.3 and 12.6% of the variation, respectively (Figure 1A). As expected, the separated plot representing high and low ADG can be observed in the PCA plot. Next, PLS-DA was executed to exhibit the variations between the high and low ADG lambs. As shown in Figure 1B, the PLS-DA analysis demonstrated that the serum metabolites of the low ADG lambs distinctly differed from those of the high ADG lambs. Correspondingly, the values of R2Y and Q2 were 0.998 and 0.796, respectively (Figure 1), indicating good interpretability and predictability by this PLS-DA model. A value of Q2 = 1 indicates a perfect discrimination of metabolites profiles between groups.
Figure 1. Serum metabolite profile of lambs with different average daily gain (ADG). (A) PCA score plot for serum metabolites of high ADG lambs (H) and low ADG lambs (L). (B) PLS-DA score plot for serum metabolites of high ADG lambs (H) and low ADG lambs (L), and model overview showing high R2Y and Q2. (C) Enhanced Volcano plots of OPLS-DA showing fold changes (log2FC) and the negative logarithm (base 10) of the P-values of 57 differential serum metabolites in high ADG lambs compared with low ADG lambs. (D) Expression profile, VIP score, and P-value of the top 30 differential serum metabolites in high ADG lambs compared with low ADG lambs.
In comparison, a total 57 differential serum metabolites were found according to the Volcano plot (P < 0.05, VIP > 1.0, and log2FC > 1 or < 1.0), of which 35 metabolites showed up-regulation and 22 showed down-regulation (Figure 1C and Supplementary Table 1). A total of 50 differential serum metabolites were annotated in six superclasses according to the HMDB database, of which 36 belonged to lipids and lipid-like molecules, 5 belonged to organic nitrogen compounds, 4 belonged to organic acids and derivatives, 2 belonged to benzenoids, 2 belonged to organic oxygen compounds, and 1 belonged to organoheterocyclic compounds. Expression profile and VIP of the top 30 metabolites based on the OPLS-DA model are shown in Figure 1D. Regarding the KEGG pathway, the four pathways including at least two differential serum metabolites annotated were metabolic pathways [L-Histidinol, Myristic acid, D-Sedoheptulose 7-phosphate, L-Arginine, and PC(14:1(9Z)/20:2(11Z,14Z))], biosynthesis of amino acids (L-Histidinol, D-Sedoheptulose 7-phosphate, and L-Arginine), glycerophospholipid metabolism [PC(14:1(9Z)/20:2(11Z,14Z)), LysoPC(20:5(5Z,8Z,11Z,14Z,17Z)), and LysoPC(18:3(6Z,9Z,12Z))], and histidine metabolism (L-Histidinol and 3-Methyl-L-histidine), respectively. Based on our results, the serum metabolome prolife of lambs was affected by high or low ADG.
As shown in Figure 1D and Supplementary Table 1, most of the lipids and lipid-like molecules were accumulated in serum samples of high ADG lambs compared with low ADG ones, including LysoPC(20:5 (5Z,8Z,11Z,14Z,17Z)) and PC(14:1(9Z)/20:2 (11Z,14Z)), which are metabolites involved in glycerophospholipid metabolism. In detail, LysoPC(20:5(5Z,8Z,11Z,14Z,17Z)) and PC(14:1(9Z)/20:2(11Z,14Z)) had a high concentration in serum of high ADG lambs [log2(FC) = 1.04 and 1.02]. Lysophospholipids (LPL) mainly include lysophosphatidylcholine (LPC), lysophosphatidic acid (LPA), lysophosphatidylethanolamine (LPE), and lysophosphatidylinositol (LPI), which are derivatives of phospholipid with absence of a fatty acid chain by hydrolysis (20). LPL could be a potent feed additive to improve production and feed efficiency according to studies in non-ruminant animals (21, 22), as well as in ruminants (23, 24). Recently, research showed that LPL supplementation could increase ADG in lambs, potentially through altering feed digestion (25), and intermediating bacterial phospholipid turnover as one of the cellular growth factor or potent lipid mediator in bacteria (26). Moreover, LPC could alter enterocyte monolayer permeability via protein kinase C (27).
Amino acids in tissue and in serum seemed to change under various physiological status, such as starvation, fasting, grazing, and stress in sheep (14) and dairy cows (28). During starvation in sheep, circulating amino acids had a general trend to increase mainly because muscle proteins were mobilized to improve gluconeogenesis in the livers by enhancing amino acids supply (14). Three metabolites (L-Arginine, L-Histidinol, and D-Sedoheptulose 7-phosphate) were identified to increase in circulation of low ADG lambs (Supplementary Table 1). Arginine is a conditionally essential amino acid in livestock and has a potential role on regulating energy partitioning between fat and lean deposition (29). Intriguingly, in lower body weight suckling lambs affected by intrauterine growth restriction, arginine supplementation increased ADG and decreased feed conversion rate (30). In the growing period in Dorper and Damara sheep, seasonal weight loss resulting from dietary restriction resulted in an increased arginine level in liver, but a decreased level in Australia Merino, while histidine level increased in all three sheep studied during dietary restriction (14). Plasma arginine was reported to increase in feed-restricted dairy cows, too (28). Collectively, these results showed that circulating amino acid concentrations changed to satisfy requirements of growth needs and normal metabolism in sheep.
Histidine can be a substrate for gluconeogenesis and protein synthesis; however, it can also affect the active site of enzymes (31). L-Histidinol and 3-Methyl-L-histidine (annotated in histidine metabolism pathway) were down-regulated in lambs with high ADG (Figure 1D and Supplementary Table 1), indicating that histidine metabolism had a trend of less activity, potentially to support faster growth and development in lambs in an intensive fattening system. Moreover, histidine metabolism seems to be up-regulated under nutritional restriction. It was reported that plasma L-histidine level decreased in barn confinement sheep compared with free grazing sheep due to high body weight gain (32).
Collectively, the purpose of this report was to reveal serum metabolome profiles of fattening lambs in a barn feeding fattening system, with particular attention to unique metabolites in lambs with high or low ADG. Our findings showed that differential metabolites affected by ADG belonged to lipids and lipid-like molecules, organic nitrogen compounds, organic acids and derivatives, benzenoids, organic oxygen compounds, and organoheterocyclic compounds. The identified metabolites have an effect on regulating metabolic pathways, biosynthesis of amino acids, glycerophospholipid metabolism, and histidine metabolism. These results indicate that selected serum metabolites could have potential application to estimate sheep with different ADG. Further larger-size studies with more various cohorts of sheep are desired to validate our finding.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.
Ethics Statement
The animal study was reviewed and approved by Beijing Academy of Agriculture and Forestry Sciences.
Author Contributions
TF and YL conceived the study. JW and TF obtained funding. HD performed animal trials and data collection. WX and TF performed data interpretation. TF and HD wrote the manuscript. WX and ÁK performed manuscript revision. All authors read and approved the final manuscript content.
Funding
This work was supported by the Science Foundation of Institute of Animal Husbandry and Veterinary Medicine, Beijing Academy of Agriculture and Forestry Sciences (XMSSYJJ202101), the Achievement Transformation Project of Beijing Academy of Agriculture and Forestry Sciences (2018017 and 2020607), and the Key Research and Development Plan Program of Hebei Province (20326629D).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We thank all staff of the Experimental Station of Beijing Academy of Agriculture and Forestry Sciences (located in the Yangyuan county, Zhangjiakou city, Heibei province) for their assistance during the sample collection.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2021.662536/full#supplementary-material
References
1. Li Q, Lu Z, Jin M, Fei X, Quan K, Liu Y, et al. verification and analysis of sheep tail type-associated gene polymorphisms. Animals. (2020) 10:89. doi: 10.3390/ani10010089
2. Quan K, Li J, Han H, Wei H, Zhao J, Si HA, et al. Review of Huang-Huai sheep, a new multiparous mutton sheep breed first identified in China. Trop Anim Health Prod. (2020) 53:35. doi: 10.1007/s11250-020-02453-w
3. Zhang X, Wang W, Mo F, La Y, Li C, Li F. Association of residual feed intake with growth and slaughtering performance, blood metabolism, and body composition in growing lambs. Sci Rep. (2017) 7:12681. doi: 10.1038/s41598-017-13042-7
4. Arthy V, Venkataramanan R, Sivaselvam SN, Sreekumar C, Balasubramanyam D. Genetic evaluation of growth in farmers' flocks of Madras Red sheep under long-term selection in a group breeding scheme. Trop Anim Health Prod. (2018) 50:1463–71. doi: 10.1007/s11250-018-1581-z
5. Santos A, Giraldez FJ, Mateo J, Frutos J, Andrés S. Programming Merino lambs by early feed restriction reduces growth rates and increases fat accretion during the fattening period with no effect on meat quality traits. Meat Sci. (2018) 135:20–6. doi: 10.1016/j.meatsci.2017.08.007
6. Atti N, Mahouachi M. The effects of diet, slaughter weight and docking on growth, carcass composition and meat quality of fat-tailed Barbarine lambs. A review. Trop Anim Health Prod. (2011) 43:1371–8. doi: 10.1007/s11250-011-9865-6
7. Beltrão ES, de Azevedo Silva AM, Filho JMP, de Moura JFP, de Oliveira JPF, Oliveira RL, et al. Effect of different blend levels of spineless cactus and Mombasa hay as roughage on intake, digestibility, ingestive behavior, and performance of lambs. Trop Anim Health Prod. (2021) 53:140. doi: 10.1007/s11250-021-02585-7
8. Chen X, Mi H, Cui K, Zhou R, Tian S, Zhang L. Effects of diets containing finger millet straw and corn straw on growth performance, plasma metabolites, immune capacity, and carcass traits in fattening lambs. Animals. (2020) 10:1285. doi: 10.3390/ani10081285
9. Muir SK, Linden N, Knight M, Behrendt R, Kearney G. Sheep residual feed intake and feeding behaviour: are ‘nibblers' or ‘binge eaters' more efficient. Anim Prod Sci. (2018) 58:1459–64. doi: 10.1071/AN17770
10. Jackson T, Heard J, Malcolm B. System changes to a lamb farm in south-west victoria: some pre-experimental modelling. AFBM J. (2014) 11:1–18. doi: 10.22004/ag.econ.198725
11. Goldansaz SA, Markus S, Berjanskii M, Rout M, Guo AC, Wang Z, et al. Candidate serum metabolite biomarkers of residual feed intake and carcass merit in sheep. J Anim Sci. (2020) 98:skaa298. doi: 10.1093/jas/skaa298
12. Rincon-Delgado RM, Gutierrez H, Perez-Vszq ED, Muro-Reyes A, Diaz-Garci LH, Banuelos VR, et al. Relationship of residual feed intake on specific hematological and biochemical parameters in Rambouillet sheep. Agri J. (2011) 6:87–91. doi: 10.3923/aj.2011.87.91
13. Chang EI, Wesolowski SR, Gilje EA, Baker PR, Reisz JA, D'Alessandro A, et al. Skeletal muscle amino acid uptake is lower and alanine production is greater in late gestation intrauterine growth-restricted fetal sheep hindlimb. Am J Physiol Regul Integr Comp Physiol. (2019) 317:R615–29. doi: 10.1152/ajpregu.00115.2019
14. Ribeiro DM, Madeira MS, Kilminster T, Scanlon T, Oldham C, Greeff J, et al. Amino acid profiles of muscle and liver tissues of Australian Merino, Damara and Dorper lambs under restricted feeding. J Anim Physiol Anim Nutr. (2019) 103:1295–302. doi: 10.1111/jpn.13148
15. Psychogios N, Hau DD, Peng J, Guo AC, Mandal R, Bouatra S, et al. The human serum metabolome. PLoS ONE. (2011) 6:e16957. doi: 10.1371/journal.pone.0016957
16. Zhou J, Ji K, Liu H, Zhang Y, Degen AA, Jiao D, et al. Effect of air temperature on growth performance, apparent digestibilities, rumen fermentation and serum metabolites in Altay and Hu lambs. J Anim Physiol Anim Nutr. (2020) 104:1023–33. doi: 10.1111/jpn.13318
17. Jiang B, Wang T, Zhou Y, Li F. Effects of enzyme + bacteria treatment on growth performance, rumen bacterial diversity, KEGG pathways, and the CAZy spectrum of Tan sheep. Bioengineered. (2020) 11:1221–32. doi: 10.1080/21655979.2020.1837459
18. Du H, Erdene K, Chen S, Qi S, Bao Z, Zhao Y, et al. Correlation of the rumen fluid microbiome and the average daily gain with a dietary supplementation of Allium mongolicum Regel extracts in sheep. J Anim Sci. (2019) 97:2831–43. doi: 10.1093/jas/skz139
19. Jia P, Cui K, Ma T, Wan F, Wang W, Yang D, et al. Influence of dietary supplementation with Bacillus licheniformis and Saccharomyces cerevisiae as alternatives to monensin on growth performance, antioxidant, immunity, ruminal fermentation and microbial diversity of fattening lambs. Sci Rep. (2018) 8:16712. doi: 10.1038/s41598-018-35081-4
20. Mnasri T, Hérault J, Gauvry L, Loiseau C, Poisson L, Ergan F, et al. Lipase-catalyzed production of lysophospholipids. OCL. (2017) 24:D405. doi: 10.1051/ocl/2017011
21. Polycarpo GV, Burbarelli MF, CarÃo AC, Merseguel CE, Dadalt JC, Maganha SR, et al. Effects of lipid sources, lysophospholipids and organic acids in maize-based broiler diets on nutrient balance, liver concentration of fat-soluble vitamins, jejunal microbiota and performance. Br Poult Sci. (2016) 57:788–98. doi: 10.1080/00071668.2016.1219019
22. Wang QQ, Long SF, Hu JX, Li M, Pan L, Piao XS. Effects of dietary lysophospholipid complex supplementation on lactation performance, and nutrient digestibility in lactating sows. Anim Feed Sci Technol. (2019) 251:56–63. doi: 10.1016/j.anifeedsci.2018.12.009
23. Lee C, Morris DL, Copelin JE, Hettick JM, Kwon IH. Effects of lysophospholipids on short-term production, nitrogen utilization, and rumen fermentation and bacterial population in lactating dairy cows. J Dairy Sci. (2019) 102:3110–20. doi: 10.3168/jds.2018-15777
24. He Y, Zhong R, Cheng L, You P, Li Y, Sun X. Effects of the supplementation of lysophospholipids through pelleted total mixed rations on blood biochemical parameters and milk production and composition of mid-lactation dairy cows. Animals. (2020) 10:E215. doi: 10.3390/ani10020215
25. Huo Q, Li B, Cheng L, Wu T, You P, Shen S, et al. Dietary supplementation of lysophospholipids affects feed digestion in lambs. Animals. (2019) 9:E805. doi: 10.3390/ani9100805
26. Zheng L, Lin Y, Lu S, Zhang J, Bogdanov M. Biogenesis, transport and remodeling of lysophospholipids in Gram-negative bacteria. Biochim Biophys Acta Mol Cell Biol Lipids. (2017) 1862:1404–13. doi: 10.1016/j.bbalip.2016.11.015
27. Sawai T, Lampman R, Hua Y, Segura B, Drongowski RA, Coran AG, et al. Lysophosphatidylcholine alters enterocyte monolayer permeability via a protein kinase C/Ca2+ mechanism. Pediatr Surg Int. (2002) 187:591–4. doi: 10.1007/s00383-002-0860-x
28. Laeger T, Görs S, Metges CC, Kuhla B. Effect of feed restriction on metabolites in cerebrospinal fluid and plasma of dairy cows. J Dairy Sci. (2012) 95:1198–208. doi: 10.3168/jds.2011-4506
29. Wu G, Bazer FW, Dai Z, Li D, Wang J, Wu Z. Amino acid nutrition in animals: protein synthesis and beyond. Annu Rev Anim Biosci. (2014) 2:387–417. doi: 10.1146/annurev-animal-022513-114113
30. Zhang H, Peng A, Guo S, Wang M, Loor JJ, Wang H. Dietary N-carbamylglutamate and l-arginine supplementation improves intestinal energy status in intrauterine-growth-retarded suckling lambs. Food Funct. (2019) 10:1903–14. doi: 10.1039/c8fo01618f
31. Brosnan ME, Brosnan JT. Histidine metabolism and function. J Nutr. (2020) 150(Suppl. 1):2570–5S. doi: 10.1093/jn/nxaa079
Keywords: metabolomics, sheep, average daily gain, metabolic pathway, serum
Citation: Feng T, Ding H, Wang J, Xu W, Liu Y and Kenéz Á (2021) Metabolite Profile of Sheep Serum With High or Low Average Daily Gain. Front. Vet. Sci. 8:662536. doi: 10.3389/fvets.2021.662536
Received: 01 February 2021; Accepted: 02 March 2021;
Published: 05 May 2021.
Edited by:
Haidong Yao, Karolinska Institutet (KI), SwedenReviewed by:
Shengguo Zhao, Chinese Academy of Agricultural Sciences, ChinaFugui Fang, Anhui Agricultural University, China
Copyright © 2021 Feng, Ding, Wang, Xu, Liu and Kenéz. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yan Liu, liuyanxms@163.com; Tao Feng, fengtao_gs@163.com