- 1National Engineering Laboratory for Animal Breeding, Key Laboratory of Animal Genetics and Breeding of the Ministry of Agricultural, Beijing Key Laboratory for Animal Genetic Improvement, College of Animal Science and Technology, China Agricultural University, Beijing, China
- 2Inner Mongolia Golden Grassland Ecological Technology Group Co., LTD., Inner Mongolia, China
- 3College of Animal Science, Xinjiang Agricultural University, Urumqi, China
Melatonin is an indole-like neuroendocrine hormone. A large number of studies have shown that melatonin can improve production performance of ewes, but it is not clear in lambs. In this study, the growth and development of the 2-month-old lambs implanted with melatonin were monitored for 60 days. The results showed that the growth rate of body weight and body skew length of lambs with melatonin treatment were significantly improved compared to the controls. The similar results were also observed in red blood cell count, hematocrit, red blood cell volume distribution width, the levels of growth hormone, testosterone, immunoglobulin A, immunoglobulin M and albumin. In addition, the cross sectional area of muscle fibers and adipose cells of lambs with melatonin implantation were also significantly increased compared to the controls (P<0.05). To further explore the potential mechanisms, the muscle and adipose tissue were selected for transcriptome sequencing. KEGG enrichment results showed that melatonin regulated the expression of genes related to apoptotic signaling pathway in muscle and adipocytes. Since the intestinal microbiota are involved in the nutritional balance and animal growth, the 16SrRNA sequencing related to the intestinal microbiota was also performed. The data indicated that the structural differences of fecal microflora mainly occur in the pathways of Cardiovascular disease, Excretory system and Signaling molecules and interaction. In brief, melatonin promotes the growth and development of lambs. The potential mechanisms may be that melatonin increased the growth hormone and testosterone mediated apoptosis signaling pathway and regulated intestinal microbial flora. Our results provide valuable information for melatonin to improve the production of sheep husbandry in the future.
Introduction
Melatonin is a tryptophan derivative, mainly secreted by the pineal gland of vertebrates at night and it is also widely distributed in bacteria and plants (1). Melatonin has a variety of functions in mammals, such as regulating biological rhythm, participating in immune regulation, inhibiting tumorigenesis, reducing oxidative stress and delaying aging (2).. Many physiological functions of melatonin are mediated primarily through the activation of receptors MT1 and MT2 (3). As to its functions in reproductive system, melatonin regulates seasonal reproductive activities of sheep such as modification of estrus cycle through the hypothalamic-pituitary-gonad axis (4). For example, in vitro, melatonin supplementation improves the maturation rate and cleavage rate of sheep oocytes by exerting its antioxidant function (5). In vivo melatonin subcutaneously injection during estrus increased the level of estradiol in the blood of ewe, and promoted embryo development as well as the pregnancy rate in the awes with natural mating or embryo transferring (6). Melatonin also improved the reproductive capacity of rams by modulation of testicular function (7). Evidence has been accumulated to show melatonin’s great application value in the sheep husbandry (8). The mechanisms of melatonin on mammalian physiology are multiples. These include its antioxidant and anti-inflammatory activities, energy metabolism, biorhythm regulation and mitochondrial biogenesis (9). Some of these activities are mediated by melatonin receptors and others are receptor-independent (10). In addition, some actions vary based on the dosages and delivery routes.
There are several melatonin delivery routes to improve the growth and reproductive performance of sheep, but a predominant route is subcutaneous implantation. By use of this method during cashmere production cycle, melatonin improved regeneration rate of skin secondary hair follicles and cashmere production (11). In the non-breeding season, melatonin regulated the estrus time of ewes and advances their estrus (12). Combined the controlled internal drug release devices and equine chorionic gonadotrophin (CIDR-ECG) with melatonin implantation in Anglo-Nubian goats significantly increased their luteal number, progesterone concentration, pregnancy rate and fertility (13). In the non-breeding season (at the beginning of May) melatonin implantation significantly improved the endocrine and exocrine function of testicles and scrotal circumference as well as quality parameters of the ejaculate in the Hungarian native breed Black Racka rams (14). Lamb weaning stress has always been an important issue which directly limits sheep production. As we know that melatonin implantation in pregnant ewes promoted fetal growth and reduced mortality of newborn lambs (15). However, the effects of melatonin on growth and development of lambs after weaning has not been investigated.
The growth and development of lambs involve in muscle and fat growth, immune activities and intestinal microflora regulation. Melatonin can regulate the release of inflammatory factors by interfering T cell differentiation and humoral and cellular immune activities (16). Melatonin can inhibit TLR4-mediated IL-6 expression in ewes choroid plexus (17). Melatonin is mainly synthesized in mitochondria (18). Mitochondria is the main site of energy metabolism (19). Melatonin maintains mitochondrial homeostasis and the ATP synthesis (20). Therefore, melatonin affects the metabolic activities of the body through its function on mitochondria. Evolutionally, mitochondria are the decedents of the bacteria of purple nonsulfur bacteria (R. rubrum) which have the capacity to synthesize melatonin (21). This principal is also applied to the intestinal microflora. Intestinal microflora not only affects nutrient absorption, immune regulation, but also energy metabolism in mammals (22, 23). Melatonin can modify the distribution and structures of intestinal microflora in suckling piglets (24). We hypothesize that the growth and development of lambs will also be influenced by their intestinal microflora which can be regulated with melatonin supplementation.
In the study, we will investigate whether melatonin administration will impact the growth and development of lambs. If it does, what are the potential mechanisms including the signal transduction pathways? Whether the crosstalk between melatonin and intestinal microbiota is involved in this process? To clarify these fundamental issues will provide valuable information not only for basic research but also for increasing the production in sheep husbandry.
Materials and methods
Experimental animals
The lambs of the Hu Sheep (Ovis aries) were from Golden Grassland Ecological Technology Co., LTD., Autonomous Region of Inner Mongolia. A total of 120 healthy 2-month-old male lambs with the body weight of about 19kg were selected. All the animal experimental design was in accordance with the regulations of animal Welfare Committee of China Agricultural University (No : AW01502202-1-2).
Chemicals and regents
Follicle-Stimulating Hormone, FSH (B03PZB), Growth Hormone, GH (B12PZB), Luteinizing Hormone, LH(B04PZB), Androgens,T (B10PZB) and cortisol, COR (COR-D10PZB) were purchased from Beijing North Biotechnology Research Institute Co., LTD.; Estradiol, E2(B05PZB) was purchased from Tianjin Union Medical Science and Technology Co., LTD.; Superoxide dismutase, SOD ((A001-1-2), Malonaldehyde, MDA(A003-1-2), Glutathione peroxidase, GSH-PX (A005-1-2), Catalase, CAT (A007-1-1), Immunoglobulin m, IGM (Jh-00015), Immunoglobulin g, IGG (Jh-00013) and Immunoglobulin a, IGA(Jh-00014) were purchased from Nanjing Jiancheng Biological Engineering Research Institute Co., LTD.; Glutamic-pyruvic transaminase, ALT (Alanine substrate method, B2008),Glutamic oxalacetic transaminase, AST (aspartic acid substrate method, B2009), Total protein, TP (biuret method, B2010), Albumin, ALB (bromocresol green method, B2011), alkaline phosphatase, ALP (NPP substrate-AMP buffer method, B2014), Urea nitrogen, UREA(Urease-glutamate dehydrogenase method, B2017), Serum creatinine, CRE (Creatine oxidase method, B2018), Uric acid, UA (uricase method, B2019), Glucose, GLU (Glucose oxidase method,. B2007) were purchased from Beijing Beijian•xinchuangyuan Biotechnology Co., LTD.
Experimental design
Lambs were randomly divided into 3 groups which was with 40 lambs for each group. Melatonin was implanted subcutaneously in the neck (Institute of Special Products, Chinese Academy of Agricultural Sciences, No: 2019/4/18) with 0mg, 3mg/kg and 4.5mg/kg, respectively for each lamb in each group. The lambs were fed with the total mixed diet at 9:00 and 17:00, but accessed to the water freely. The ratio of concentrate to coarse feed was 1:3. The basic components of concentrate feed include corn,soybean meal,corn dry distiller’s grains,calcium carbonate,vitamins and trace elements. Nutrients include crude protein ≥18%, crude fiber ≥15%, ash ≤10%, calcium 0.5%-20%, total phosphorus 0.3%, sodium chloride 0.3%-20%, moisture ≤13%, lysine ≥0.5%. The body height, oblique length, chest circumference, chest width, chest depth, anterior duct and posterior duct were measured morning on day 0d, 30d and 60d after an overnight fasting.
Routine blood test
5mL blood was collected from the jugular vein of the lambs, on 60d of the study. The samples were stored in EDTA tubes and transported at 4°C. Blood routine tests were performed with blood cell analyzer (Shenzhen Mindray Bio-medical Electronics Co., LTD., BC2800Vet, Shenzhen, China).
Detection of hormones and physiological and biochemical parameters
5mL blood was collected from the jugular vein of the lambs on day 0d, 30d and 60d. The samples were kept at room temperature for 2h, centrifuged at 3000r/min for 8min, and the serum was stored at -20°C. All samples were tested according to the instructions of the kit as above 1.1. Chemicals and regents.
(1) Reproductive hormones [FSH, GH, LH, T, E2] and Antioxidant factors [SOD, MDA, GSH-PX, CAT, COR] were detected by radioimmunoassay (Xi’an Nuclear Instrument Factory, XH6080; Xian, China) in Beijing North Biotechnology Research Institute Co., LTD.
(2) Physiological and biochemical parameters: ALT, TP, ALP, UREA, CRE, UA, GLU were measured following the manufacturer’s instructions.
(3) Immune factors: IGM, IGG, IGA were detected by 721 spectrophotometer and colorimetry.
(4) Melatonin: The serum samples were mixed with methanol in a ratio of 1:4 and then oscillated in a vortex. After centrifugation (12000 r/min for 10 min), the supernatant was collected and filtered into the membrane for use. Melatonin detection was performed using high performance liquid mass spectrometry (AGilent1290-G6470, Santa Clara, CA, USA) at the Central Laboratory of The Institute of Animal Science, Chinese Academy of Agricultural Sciences (Beijing, China).
Tissue staining
The sheep were stabilized, the surface of the longissimus muscle of the back was sheared, disinfected with iodophor, and then deiodized with 75% alcohol. Using a scalpel, slice the sterilized skin and fat in turn, using forceps and direct scissors to obtain the longest muscle and fat in the back. The longissimus muscle and back adipose tissue were fixed with 4% paraformaldehyde (Beijing Regen Biotechnology Co., LTD., No. DF0135) for 24h, then the tissue and label were placed in the dehydration box. The dehydration box is put into the dehydrator (Wuhan Junjie Electronics Co., LTD., JJ-12J, Wuhan, China) and then the gradient alcohol (Sinopagic Chemical Reagent Co., LTD., 100092683) is used for dehydration, and then the tissue is dipped in wax. The wax-soaked tissue is embedded in the embedding machine (Wuhan Junjie Electronics Co., LTD., JB-P5, Wuhan, China). The trimmed wax blocks were cooled at -20°C (Wuhan Junjie Electronics Co., LTD., JB-P5, Wuhan, China) and sliced in a paraffin slicer (Shanghai Leica Instrument Co., LTD., RM2016, Shanghai, China) with a thickness of 4μm. The slices were floated in the spreading machine (Jinhua Kedi Instrument Equipment Co., LTD., KD-P, Zhejiang, China) and the tissues were flatten in water at 40°C. The slides (Servicebio, G6004) were picked up and the slices were baked in the oven at 60°C (Tianjin Laiborey Instrument Equipment Co., LTD., GFL-230, Tianjin, China). Water drying wax baking out of room temperature storage for later. The slices were then placed in xylene I (Sinopagic Chemical Reagent Co., LTD., 10023418) for 20min, xylene II for 20min, anhydrous ethanol I for 5min, anhydrous ethanol II for 5min, 75% alcohol for 5min, and washed with water. The slices were stained with hematoxane dye (HE dye kit purchased from Servicebio, G1003, Wuhan, China) and sealed with neutral gum (Sinopharming Chemical Reagents Co., LTD., 10004160). Microscopic examination and image analysis were performed using Nikon Eclipse E100 (Tokyo, Japan).
RNA-seq and qPCR
Total RNA was extracted from muscle and adipose tissue by TRIzol (Ambion,15596026), and genomic DNA was removed by DNase I (TaKara). The quality of RNA samples was detected by 2100 Bioanalyser (Agilent) and ND-2000(NanoDrop Technologies). After the samples were qualified (OD260/280 = 1.8-2.2, OD260/230≥2.0, RIN≥6.5), 28S:18S≥1.0, >2μg) were sequenced.
TruSeqTM RNA Sample Preparation Kit (Illumina, San Diego, CA) was used for the construction of RNA library. SuperScript double-stranded cDNA synthesis Kit (Invitrogen, CA) was used to inversely synthesize cDNA and form a stable double-stranded structure. After cDNA enrichment by PCR (sample Preparation Kit (Illumina, San Diego, CA) Kit), DNA clean Beads (DNA Clean Beads) are screened for 200-300 bp bands. After quantified by TBS380 (Picogreen), Illumina HiSeq XTEN/NovaSeq 6000 sequencing platform was used for high-throughput sequencing with a read length of PE150.
QPCR procedure: Incubation at 95°C for 3 min (activation of FastStart DNA polymerase), 95°C for 10 seconds, 60°C for 30 seconds. 39 cycles were obtained with appropriate extension time for a single fluorescence acquisition. FAS, P53, GADD45B, PARP1, CytC, BAx, Bcl-2, APAF1, FAS-L, TNF-α, IL1B, Bcl-XL, ATM, IL3, CASP3, CASP8 were selected for verification, and GAPDH was used as control. All gene validation was repeated for 3 times, and gene expression was evaluated by 2-△△CT method. The primers were verified as shown in Table 1. Additional Requirements
16S intestinal microbiome sequencing
Rectal feces were collected in a sterile environment, frozen with liquid nitrogen and stored at -80°C. Stool microbiome DNA extraction method refer to PowerSoil DNA Isolation Kit (MoBio Laboratories, Carlsbad, CA) [Omega DNA Kit]. The samples are stored at -20°C after passing quality inspection. Primers F (5’-ACTCCTACGGGAGgCAGCAGcag-3’) and R (5’-GGACTACNNGGG TATctaat-3’) were used to amplify V3-V4 region of bacterial 16SrRNA gene. PCR products were used to construct microbial diversity sequencing library. Illumina Miseq PE300 high-throughput sequencing platform was used for Paired end sequencing in Beijing Alvesen Gene Technology Co., LTD. QIIME1 (V1.8.0) software was used to split the samples according to Barcode sequence, and Pear (V0.9.6) software was used to filter and splice the data. The uPARSE algorithm of Vsearch (V2.7.1) software was used to perform OTU clustering (Operational Taxonomic Units) for high-quality sequences, and the similarity threshold was 97%. The RDP Classifier algorithm was used for comparison with Silva128 database, and the 70% confidence threshold was set to obtain the species classification information corresponding to each OTU. Venn plots, rank-abundance curves and species composition histogram were analyzed using R (V3.6.0) software based on species annotation and relative Abundance results. R (V3.6.0) software was used for PLS-DA analysis and partial least square discriminant analysis. QIIME1 (V1.8.0) was used to calculate the beta diversity distance matrix, R (V3.6.0) software was used for clustering heat map analysis based on Weighted Unifrace distance, and Phython (V2.7) software was used for LEfSe analysis.
Statistical analysis
The data were sorted by Excel. The results were expressed as mean ± SEM. One-way ANOVA was used to analyze the normality of the samples and followed by Duncan test to determine the statistically significant difference between the groups with the help of SPSS 26.0(IBM SPSS Statistics, Armonk, NY, USA). P<0.05 was taken as the criterion of significance of difference.
Results
Effects of melatonin on growth performance of lambs
The body weight (Figure 1A), body oblique length (Figure 1C), chest circumference (Figure 1D), chest width (Figure 1E) and chest depth (Figure 1F) in melatonin implanted group were significantly greater than those in the control group (P<0.05) with the better results in 4.5 than 3 mg/kg melatonin treated group. There were no significant differences in body height (Figure 1B), anterior canal circumference (Figure 1G) and posterior canal circumference (Figure 1H) among all groups (P>0.05).
Figure 1 Effects of melatonin on growth performance of lambs. (A) body weight; (B) body height; (C) body oblique length; (D) chest circumference; (E) chest width; (F) chest width; (G) anterior tube circumference; (H) posterior tube circumference Different letters indicate P<0.05 vs control.
Effects of melatonin on the parameters of physiology and biochemistry in lambs
After 30 days of melatonin implantation, the albumin content of lambs increased gradually and was significantly higher than that of control group (P<0.05). No significant differences on globulin and white bulb ratio were detected among the groups (P>0.05) (Table 2). Interestingly, the RBC count, hematocrit and RBC volume distribution width were significantly higher in lambs with 4.5mg/kg melatonin implantation than those in the 3mg/kg group (P<0.05) while these differences were not observed between the lambs with melatonin implantation and controls (P>0.05) (Table 2). WBC count, mean RBC volume, mean RBC hemoglobin content, mean PLATELET volume, platelet distribution width and platelet hematocrit had no significant differences among all groups (P>0.05) (Table 3) and these were the case on levels of glucose, urea and uric acid among groups (P>0.05) (Table 4).
Effects of melatonin on reproductive hormones in lambs
Melatonin implantation significantly increased the blood melatonin levels of the lambs compared to the controls (P<0.05) (Figure 2A). The growth hormone (Figure 2B), estradiol (Figure 2C) and androgen (Figure 2D) in lambs with 4.5mg/kg of melatonin implantation were significantly increased compared to the controls (P<0.05). After 1 month melatonin implantation, androgen level in 4.5mg/kg group was significantly higher than that of 3mg/kg group (P<0.05) and at 2 months, the growth hormone level had the similar change as the androgen (P<0.05). The levels of luteinizing hormone releasing hormone (Figure 2E) and follicle-stimulating hormone releasing hormone (Figure 2F) had no significant difference among groups (P>0.05).
Figure 2 Effects of melatonin on blood melatonin levels and reproductive hormones in lambs. (A) Melatonin levels; (B) Growth hormone; (C) Estradiol; (D) Androgens; (E) Follicle-stimulating hormone; (F) Luteinizing hormone-releasing hormone. Note: Different lowercase letters indicate P<0.05 vs control.
Effects of melatonin on antioxidant and immune factors in lambs
Cortisol (Figure 3A) and malondialdehyde (Figure 3C) were significantly lower and superoxide dismutase (Figure 3B) was significantly higher in lambs with 4.5mg/kg melatonin implantation than those in other groups (P<0.05). Glutathione peroxidase (Figure 3D) and catalase (Figure 3E) were not significantly differences among groups (P>0.05). The levels of immunoglobulin A (Figure 4A) and immunoglobulin M (Figure 4C) in lambs with 4.5mg/kg melatonin implantation were significantly higher than those in other groups (P<0.05). But, there was no significant difference in the immunogiobulin G (Figure 4B) among the groups.
Figure 3 Effects of melatonin on antioxidant factors of lambs. (A) cortisol; (B) superoxide dismutase; (C) malondialdehyde; (D) glutathione peroxidase; (E) catalase. Different lowercase letters indicate P<0.05, and the same lowercase letters indicate P>0.05 vs contrpl.
Figure 4 Effects of melatonin on immunoglobulins of lambs. (A) IgA; (B) IgG; (C) IgM Different lowercase letters indicate significant differences (P<0.05), and the same lowercase letters indicate insignificant differences (P>0.05).
Effects of melatonin on longissimus muscle and back adipose tissue in lambs
The cross-sectional area of longismus muscle tissue was significantly larger (Figure 5C) in lambs with 90 mg melatonin implantation than that in the other groups (Figures 5A, B) (P<0.05). The cross-sectional area of adipose tissue implanted with 3mg/kg (Figures 6B) and 4.5mg/kg (Figures 6C) melatonin was significantly higher than that of the control group (Figures 6A). The cross-sectional areas of adipose tissue were significantly larger in lamb with 60 and 90 mg melatonin implantations than that of the control group (Figures 5D, 6D).
Figure 5 Cross-sectional muscle fiber diagram of the back muscles. (A) The cross section of the longissimus muscle in control group, 200×, 50×, and 20×, respectively; (B) The transverse section of the longissimus muscle in 3mg/kg group, 200×, 50×, and 20×, respectively; (C) The cross section of the longissimus muscle in 4.5mg/kg group, 200×, 50×, and 20×, respectively; (D) Statistical analysis of the cross-sectional area of the longissimus dorsi. Different lowercase letters indicate P<0.05 vs control.
Figure 6 Cross-sectional cytogram of dorsal adipose tissue. (A) Transverse sections of the back adipose tissue in the control group, 200×, 50×, and 20×, respectively; (B) Transverse sections of the back adipose tissue in the 3mg/kg group, 200×, 50×, and 20×, respectively; (C) The cross section of the back adipose tissue in the 4.5mg/kg group, 200×, 50×, and 20×, respectively; (D) Statistical analysis of the cross-sectional area of dorsal adipose tissue. Different lowercase letters indicate P<0.05 vs control.
Transcriptome analysis of back longissimus muscle and adipose tissue of lambs
The transcriptome analyses were performed from the back longissimus muscles and adipose tissues of lambs with 4.5mg/kg melatonin implantation and controls. The significant differences in gene expressions in back longissimus muscle and adipose tissue have been observed between the melatonin and the control groups (Figures 7A, 8A). A total of 346 differentially expressed genes were detected in muscle tissue, of which 237 were up and 109 were down-regulated (Figure 7B). A total of 492 differentially expressed genes were detected in adipose tissue, of which 284 were up and 208 were down-regulated (Figure 8B). GO enrichment analysis showed that the differentially expressed genes in muscle tissue were mainly enriched in extracellular matrix and positive regulation of metabolic process pathway (Figure 7C). KEGG enrichment analysis showed that these differentially expressed genes were mainly enriched in Apoptosis and Base excision repair pathway (Figure 7D). GO enrichment analysis showed that the differentially expressed genes in adipose tissue were mainly enriched in positive regulation of cell death and response to lipopolysaccharide pathway (Figure 8C). KEGG enrichment analysis showed that these genes were mainly enriched in Apoptosis and Measles pathways (Figure 8D). Based on the KEGG enrichment analysis, differentially expressed genes in muscle and adipose tissues both were convergently enriched in the Apoptosis pathway. Therefore, sixteen genes in Apoptosis pathway were selected for validation. The expression levels of CytC, GADD45B, IL-3, PARP1 and TNF-ɑ in muscles were significantly higher and IL1-β were significantly lower in lambs with melatonin implantation than those in the controls (Figure 9). The expression levels of CytC, FAS, P53 and PARP1 in adipose tissue were significantly higher and IL1-β was significantly lower in lambs with melatonin implantation than those in controls (Figure 10). The expression of CytC, PARP1 and IL1-β genes in both muscle tissue and adipose tissue showed the similar changes after melatonin intervention in the Apoptosis pathway.,
Figure 7 Transcriptome sequencing of back muscle tissue. (A) Cluster analysis; (B) Expression difference analysis; (C) GO enrichment analysis; (D) KEGG enrichment analysis.
Figure 8 Transcriptome sequencing of back muscle tissue. (A) Cluster analysis; (B) Expression difference analysis; (C) GO enrichment analysis; (D) KEGG enrichment analysis.
Figure 9 Quantitative analysis of apoptosis and collateral genes in back muscle tissue (A) APAF1; (B) BAX; (C) BCL-XL; (D) CASP8; (E) ATM; (F) BCL-2; (G) CASP3; (H) CytC; (I) FAS; (J) FAS-L; (K) GADD45B; (L) IL-1β; (M) IL-3; (N) P53; (O) PARP1; (P) TNF-ɑ; Different lowercase letters indicate P<0.05 vs control.
Figure 10 Quantitative analysis of apoptosis dredging genes in dorsal adipose tissue. (A) APAF1; (B) BAX; (C) BCL-XL; (D) CASP8; (E) ATM; (F) BCL-2; (G) CASP3; (H) CytC;(I) FAS; (J) FAS-L; (K) GADD45B; (L) IL-1β; (M) IL-3; (N) P53; (O) PARP1; (P) TNF-ɑ; Different lowercase letters indicate P<0.05 vs control.
Effects of melatonin on intestinal microflora of lambs
In order to explore the effect of melatonin on the intestinal microflora structure in the lambs, 6 lambs with 4.5mg/kg melatonin implantation and 6 lambs as untreated controls were selected for 16SrDNA sequencing, respectively. A total of 1902 Operational Taxonomic Units (OTUs) were identified by sequencing, among which 369 differences OTUs (Figure 11A) were identified. The values of Rank curve decreased steadily, indicating the good species uniformity and richness (Figure 11B). The analysis of species composition showed that there was no significant difference in the structure of microflora at the genus level among samples. But the proportion of different microflora among individuals was different (Figure 11C). Partial Least Squares Discrimination Analysis (PLS-DA) showed the significant differences of microbial composition between the melatonin treated lambs and the controls, the PC1 and PC2 were 13.54% and 11.87%, respectively (Figure 11D). LDA Effect Size analysis was used to evaluate the structure of microbial flora. When LDA score was set up as 3, melatonin treated and control groups were divided into 5 and 6 different microflora (Figure 11E) respectively, among them o_Clostridia_vadiaBB60_group had the most obvious influence. In addition, PICRUSt was used to analyze microbial community function. Based on the data of KEGG analysis, the differences in structural and functional pathways of fecal microflora are mainly concentrated in the pathways of Cardiovascular disease, Excretory system and Signaling molecules and interaction (Figure 11F). Melatonin treatment may accelerate the rate of Cardiovascular and cerebrovascular metabolic cycle. Therefore, the differences in Cardiovascular disease pathway are the most obvious.
Figure 11 Sequencing analysis of gut microbiota. (A) Venn diagram; (B) Rank abundance diagram; (C) Histogram of genus-level species composition analysis; (D) OTU-based partial least squares discriminant analysis; (E) LDA distribution of LEfSe analysis based on taxonomic information Histogram; Species with differences in LDA score greater than 3, species with significant differences in abundance in different groups, the length of the histogram represents the impact size of significantly different species; (F) PICRUSt predicted differential pathway heat map A.
Discussion
Sheep is one of the most important species in animal husbandry, and mutton is one of the main meats in human diet. With the rapid development of social economy, sheep farms and associated industry are upsized to extremely large-scale. Traditional feeding methods no longer meet the demand of consumer for high-quality mutton. To produce high-quality mutton at low cost and meet consumer’s increased demand are huge challenge for scientists. Melatonin is a naturally occurring molecule present in almost all organisms with a variety of biological functions (25). A large number of studies have shown that melatonin can improve the efficiency of sheep reproductive performance (26, 27). The growth retardation of lambs after weaning is a frequently faced problem in large-scale breeding farms. To explore this problem, 120 post-weaning lambs were selected in this study to evaluate whether melatonin implantation will impact the growth retardation of lambs related to the weaning. Based on the pilot study, to achieve the obvious effects on the growth and development of lambs required around 60 days after melatonin implantation. Thus the 60 days or experimental period was selected.
The results showed that the body weight, body length, chest circumference, chest width and chest depth of lambs with melatonin implantation were significantly increased compared to the controls, and the effect was better in the 4.5mg/kg than that in 3mg/kg group. Many previous studies have reported the methods to improve growth and development of lambs. However, majority of these methods focused on the feed additives and antibiotics (28, 29). For example, dietary supplementation of oregano and cobalt mixture significantly improved growth performance, feed conversion, carcass weight, slaughter rate and immunity of lambs (30). However, these xenobiotics may have some unfavorable metabolic consequences for human body (31). In current study, melatonin was selected since melatonin is a naturally occurring molecule found in human body and with minimal side effects (32). Compared with other synthetic growth promoting products, melatonin has huge safety margin to animals and humans. In addition, previous studies have reported the regulatory effects of melatonin on sheep reproductive activities (33). Here, we reported that melatonin implantation at a certain dosage significantly improved the growth performance of lambs.
Firstly, we evaluated the safety parameters of the different dosages of melatonin on lambs. The results showed that alanine aminotransferase and albumin in 4.5mg/kg group were significantly higher than those in control group (P<0.05), but there were no significant differences in glucose, urea and uric acid level compared to the control (P>0.05). When we look at physiological markers in the blood, the results are variable in different melatonin dose groups. The RBC count, hematocrit and RBC distribution width in 4.5mg/kg group were significantly lower than those of the control group. However, no significant differences were detected between 3mg/kg group and control. Alanine aminotransferase (34), albumin (35) and RBC number (36) are all important indicators of the safety of a treatment., Melatonin treatment did not significantly altered pathophysiological states of lambs, As a result, melatonin treatment will not affect the normal growth of lambs.
The growth hormone, estradiol and testosterone in lambs with melatonin implantation were significantly increased compared to the controls. These hormones are often associated with muscle growth and body structural formations (37, 38). The effect of melatonin on the increased plasma growth hormone (GH) concentration was reported in Japanese quails (39). Melatonin is reported to reduce androgen level in male Siberian hamsters (40). These results suggest that melatonin can affect hormonal levels in animals. But, the observation is different with our findings and this may relate to the species specific. Melatonin is a potent antioxidant (41). As we observed in the study that it reduced the level of malondialdehyde and increased the activity of SOD. These actions of melatonin may prevent the muscle degradation under oxidative stress of lambs after weaning (42). Most importantly, it was observed that melatonin treatment increased the cross-sectional area of muscle fiber and adipose tissue cells. This is the direct evidence to show how melatonin improved the growth performance of lambs. Genetic information showed that the differentially expressed genes in muscle and adipose tissues are enriched in many functional pathways after melatonin treatment. KEGG analysis found that the most significant differentially expressed genes in muscle and adipose tissues impacted by melatonin was Apoptosis pathway. Apoptosis, which removes aged and damaged cells during development and protects cells from external interference and maintains homeostasis (43). Several studies have reported the effects of melatonin on apoptosis (44, 45). It is possible that the effects of melatonin on muscle growth and development might involve in gene expression of apoptotic signaling pathway.
The muscle and fat transcriptome transcription analysis further demonstrated that melatonin regulated the Apoptosis pathway function. Sixteen genes in Apoptosis pathway were selected for expression verification. Three genes showed the same expression pattern in muscle tissue and adipose tissue with melatonin treatment. These included CytC, PARP1 and IL1-β. CytC is an indispensable element of the electron transport chain in mitochondrial for energy metabolism (46). PARP1 gene is the hub of intracellular inflammatory response, oxidative stress and metabolic response (47). IL1-β is one of the important factors in the inflammatory regulatory network, initiating and controlling the inflammatory response (48). CytC, PARP1 and IL1-β are key genes to maintain metabolism and immunity. These results provide further evidence that melatonin may affect the growth and development of lambs by regulating the expression of these genes in the Apoptosis pathway. In addition to these internal pathways, some external pathways may also participate in the growth performance of lambs. This includes the intestinal microbiota of the lambs. It is well known that these microbiota affect the nutritional and physiological processes of hosts (49). The gastrointestinal tract of mammals inhabits trillions of microorganisms, and different types of intestinal microorganisms maintain a relatively stable number to maintain nutrient absorption (50). It has been reported that melatonin regulates the distribution and metabolism of intestinal microorganisms (51, 52). The crosstalk between the microbiota and the tissues of the host will influence the growth and development status of the host (53). It has been showed that melatonin treatment increased the relative abundance of families Prevotellaceae, Muribaculaceae while it reduced the relative abundance of Succinivibrionaceae, Veillonellaceae in ruminal of lactating cows (54). In this study, it was found that melatonin also affected rumen microflora of lambs. PICRUSt analysis showed that melatonin mainly affected the Cardiovascular disease, Excretory system and Signaling molecules and interaction pathways related to the microbiota. These pathways significantly affect the body health and nutrient absorption of lambs. This suggests that melatonin may also enhance nutrient absorption and promote muscle and fat development in lambs by partially influencing gut microbes.
In conclusion, in this study, we first reported that melatonin implantation im-proved the growth rate of lambs by increasing the sizes of their muscle and fat cells. The potential mechanisms are multiples. These include melatonin increases the growth hormone and androgen secretion, reduces the oxidative stress, enriches the expressions of genes in Apoptosis pathway, such as CytC, PARP1 and IL1-β genes which are the key genes to maintain metabolism and immunity. Finally, melatonin regulates the intestinal microbiota to target their Cardiovascular disease, Excretory system and Signaling molecules and interaction pathways.
In conclusion, melatonin implantation improved the growth rate of lambs by increasing the sizes of their muscle and fat cells. The underlying mechanisms are illustrated in Figure 12. Melatonin increases the growth hormone and testosterone secretion, reduces the oxi-dative stress in lambs. Melatonin enriches the expressions of genes in Apoptosis pathway, especially in CytC, PARP1 and IL1-β genes. Melatonin regulates the intestinal microbiota to target their Cardiovascular disease, Excretory system and Signaling molecules and interaction pathways. These findings provide novel information for the application of melatonin to improve the production of sheep husbandry in the future.
Data availability statement
The 16S rRNA microbial sequencing data and Transcriptome sequencing analysis data have been successfully submitted to the National Center for Biotechnology Information. SRA data: PRJNA865106 and PRJNA866689.
Ethics statement
The animal study was reviewed and approved by Animal Welfare and Experimental Animal Ethics Committee of China Agricultural University (No:AW01502202-1-2). Written informed consent was obtained from the owners for the participation of their animals in this study.
Author contributions
WM and HW designed the experiments and prepared the original manuscript. GDL and LY performed muscle and fat samples were collected and transcriptome sequencing analysis. LW, MZ, and SG collected and analyzed blood samples. SX and XG are responsible for feeding and management of lamps. FL participates in experimental design and provides technical services. PJ participated in the data analysis of physiological and biochemical indicators. AW participates in experimental design and editing the manuscript. GSL participates in experimental design, supervised the project and revised the article. All authors contributed to the article and approved the submitted version.
Funding
National Key Research and development Program (2021YFD1300903), the Science and Technology Major Project of Inner Mongolia (2021ZD0023-1); Key Research and Development Project of Hainan Province (ZDYF2021XDNY174); National Transgenic Key Project of the Ministry of Agriculture of China(2018ZX0800801B).
Acknowledgments
This research thank for National Key Research and development Program (2021YFD1300903), the Science and Technology Major Project of Inner Mongolia (2021ZD0023-1); Key Research and Development Project of Hainan Province (ZDYF2021XDNY174); National Transgenic Key Project of the Ministry of Agriculture of China(2018ZX0800801B).
Conflict of interest
Authors SX, XG, and FL were employed by Inner Mongolia Golden Grassland Ecological Technology Group Co., LTD., China.
The remaining 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.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1. Calvo JR, Gonzalez-Yanes C, Maldonado MD. The role of melatonin in the cells of the innate immunity: a review. J Pineal Res (2013) 55:103–20. doi: 10.1111/jpi.12075
2. Tordjman S, Chokron S, Delorme R, Charrier A, Bellissant E, Jaafari N, et al. Melatonin: Pharmacology, functions and therapeutic benefits. Curr Neuropharmacol (2017) 15(3):434–43. doi: 10.2174/1570159X14666161228122115
3. Zlotos DP, Jockers R, Cecon E, Rivara S, Witt-Enderby PA. MT1 and MT2 melatonin receptors: ligands, models, oligomers, and therapeutic potential. J Medicinal Chem (2014) 57:3161–85. doi: 10.1021/jm401343c
4. Pool KR, Rickard JP, de Graaf SP. Overcoming neuroendocrine and metabolic barriers to puberty: the role of melatonin in advancing puberty in ewe lambs. Domest Anim Endocrinol (2020) 72:106457. doi: 10.1016/j.domaniend.2020.106457
5. Casao A, Abecia J, Pérez JC, Blanco TM, Vázquez M, Forcada F. The effects of melatonin on in vitro oocyte competence and embryo development in sheep. Spanish J Agric Res (2010) 8:35–41. doi: 10.5424/sjar/2010081-1141
6. Song Y, Wu H, Wang X, Haire A, Zhang X, Zhang J, et al. Melatonin improves the efficiency of super-ovulation and timed artificial insemination in sheep. PeerJ (2019) 7:e6750. doi: 10.7717/peerj.6750
7. Pool KR, Rickard JP, Pini T, de Graaf SP. Exogenous melatonin advances the ram breeding season and increases testicular function. Sci Rep (2020) 10:9711. doi: 10.1038/s41598-020-66594-6
8. Palacin I, Forcada F, Abecia J. Meta-analysis of the efficacy of melatonin implants for improving reproductive performance in sheep. Spanish J Agric Res (2011) 9:730–43. doi: 10.5424/sjar/20110903-348-10
9. Cipolla-Neto J, Amaral FG, Afeche SC, Tan DX, Reiter RJ. Melatonin, energy metabolism, and obesity: a review. J Pineal Res (2014) 56:371–81. doi: 10.1111/jpi.12137
10. Emet M, Ozcan H, Ozel L, Yayla M, Halici Z. Hacimuftuoglu a. a review of melatonin, its receptors and drugs. Eurasian J Med (2016) 48:135–41. doi: 10.5152/eurasianjmed.2015.0267
11. Yang CH, Duan CH, Wu ZY, Li Y, Luan YY, Fu XJ, et al. Effects of melatonin administration to cashmere goats on cashmere production and hair follicle characteristics in two consecutive cashmere growth cycles. Domest Anim Endocrinol (2021) 74:106534. doi: 10.1016/j.domaniend
12. Cevik M, Yilmazer C. Kocyigit. a effects of melatonin implantation on the fertility potentials of kivircik and charollais ewes and rams during the non-breeding season. Pol J Vet Sci (2017) 20:501–6. doi: 10.1515/pjvs-2017-0060
13. El-Mokadem MY, El-Din A, Ramadan TA, Rashad AMA, Taha TA, Samak MA. Manipulation of reproductive seasonality using melatonin implantation in Anglo-Nubian does treated with controlled internal drug release and equine chorionic gonadotropin during the nonbreeding season. J Dairy Sci (2017) 100:5028–39. doi: 10.3168/jds.2016-12240
14. Egerszegi I, Sarlós P, Rátky J, Solti L, Faigl V, Kulcsár M, et al. Effect of melatonin treatment on semen parameters and endocrine function in black racka rams out of the breeding season. Small Ruminant Res (2014) 116:192–8. doi: 10.1016/j.smallrumres.2013.11.001
15. Flinn T, Gunn JR, Kind KL, Swinbourne AM, Weaver AC, Kelly JM, et al. Maternal melatonin implants improve twin merino lamb survival. J Anim Sci (2020) 98:skaa344. doi: 10.1093/jas/skaa344
16. Zhao CN, Wang P, Mao YM, Dan YL, Wu Q, Li XM, et al. Potential role of melatonin in autoimmune diseases. Cytokine Growth Factor Rev (2019) 48:1–10. doi: 10.1016/j.cytogfr.2019.07.002
17. Kowalewska M, Herman AP, Szczepkowska A, Skipor J. The effect of melatonin from slow-release implants on basic and TLR-4-mediated gene expression of inflammatory cytokines and their receptors in the choroid plexus in ewes. Res Vet Sci (2017) 113:50–5. doi: 10.1016/j.rvsc.2017.09.003
18. Tan D-X, Manchester LC, Qin L, Reiter RJ. Melatonin: a mitochondrial targeting molecule involving mitochondrial protection and dynamics. Int J Mol Sci (2016) 17:2124. doi: 10.3390/ijms17122124
19. Chandel NS. Mitochondria as signaling organelles. BMC Biol (2014) 12:1–7. doi: 10.1186/1741-7007-12-34
20. de Almeida Chuffa LG, Seiva FRF, Cucielo MS, Silveira HS, Reiter RJ, Lupi LA. Mitochondrial functions and melatonin: a tour of the reproductive cancers. Cell Mol Life Sci (2019) 76:837–63. doi: 10.1007/s00018-018-2963-0
21. Tan DX, Manchester LC, Liu X, Rosales-Corral SA, Acuna-Castroviejo D, Reiter RJ. Mitochondria and chloroplasts as the original sites of melatonin synthesis: a hypothesis related to melatonin’s primary function and evolution in eukaryotes. J Pineal Res (2013) 54:127–38. doi: 10.1111/jpi.12026
22. Krajmalnik-Brown R, Ilhan ZE, Kang DW, DiBaise JK. Effects of gut microbes on nutrient absorption and energy regulation. Nutr Clin Pract (2012) 27:201–14. doi: 10.1177/0884533611436116
23. Leser TD, Mølbak L. Better living through microbial action: the benefits of the mammalian gastrointestinal microbiota on the host. Environ Microbiol (2009) 11:2194–206. doi: 10.1111/j.1462-2920.2009.01941.x
24. Xia S, Gao W, Li Y, Ma J, Gong S, Gao Z, et al. Effects of melatonin on intestinal function and bacterial compositions in sucking piglets. J Anim Physiol Anim Nutr (Berl) (2022) 00:1–10. doi: 10.1111/jpn.13675
25. Szczepanik M. Melatonin and its influence on immune system. J Physiol Pharmacol (2007) 58:115–24. doi: 10.2170/physiolsci.SC008107
26. Zarazaga LA, Gatica MC, Celi I, Guzman JL. Reproductive performance is improved during seasonal anoestrus when female and male murciano-granadina goats receive melatonin implants and in payoya goats when females are thus treated. Reprod Domest Anim (2012) 47:436–42. doi: 10.1111/j.1439-0531.2011.01899.x
27. Palacín I, Abecia J-A, Forcada F, Casao A, Cebrián J-Á, Muiño T, et al. Effects of exogenous melatonin treatment on out-of-season ram fertility. Ital J Anim Sci (2016) 7:199–206. doi: 10.4081/ijas.2008.199
28. Direkvandi E, Mohammadabadi T, Salem AZ. Effect of microbial feed additives on growth performance, microbial protein synthesis, and rumen microbial population in growing lambs. Transl Anim Sci (2020) 4:txaa203. doi: 10.1093/tas/txaa203
29. Garcia-Galicia IA, Arras-Acosta JA, Huerta-Jimenez M, Rentería-Monterrubio AL, Loya-Olguin JL, Carrillo-Lopez LM, et al. Natural oregano essential oil may replace antibiotics in lamb diets: Effects on meat quality. Antibiotics (2020) 9:248. doi: 10.3390/antibiotics9050248
30. Wu JP, Zhou R, Liu LS, Casper DP, Lang X, Wang CL, et al. Growth performance, nutrient digestibility, blood parameters, and carcass characteristics by lambs fed an oregano and cobalt blend. Animal (2021) 15:100365. doi: 10.1016/j.animal.2021.100365
31. Dutta T, Yadav S, Chatterjee A. Antibiotics as feed additives for livestock: human health concerns. Indian J Anim. Heal (2019) 58:121. doi: 10.36062/ijah.58.2SPL.2019.121-136
32. Acuña-Castroviejo D, Escames G, Venegas C, Díaz-Casado ME, Lima-Cabello E, López LC, et al. Extrapineal melatonin: sources, regulation, and potential functions. Cell Mol Life Sci (2014) 71:2997–3025. doi: 10.1007/s00018-014-1579-2
33. Talpur H, Chandio I, Brohi R, Worku T, Rehman Z, Bhattarai D, et al. Research progress on the role of melatonin and its receptors in animal reproduction: A comprehensive review. Reprod Domest Anim (2018) 53:831–49. doi: 10.1111/rda.13188
34. Valenti L, Pelusi S, Bianco C, Ceriotti F, Berzuini A, Iogna Prat L, et al. Definition of healthy ranges for alanine aminotransferase levels: a 2021 update. Hepatol Commun (2021) 5:1824–32. doi: 10.1002/hep4.1794
35. Morsch CM, Gonçalves LF and Barros E. Health-related quality of life among haemodialysis patients–relationship with clinical indicators, morbidity and mortality. J Clin Nurs (2006) 15:498–504. doi: 10.1111/j.1365-2702.2006.01349.x
36. Carson JL, Grossman BJ, Kleinman S, Tinmouth AT, Marques MB, Fung MK, et al. Red blood cell transfusion: a clinical practice guideline from the AABB. Ann Internal Med (2012) 157:49–58. doi: 10.7326/0003-4819-157-1-201206190-00429
37. Mavalli MD, DiGirolamo DJ, Fan Y, Riddle RC, Campbell KS, Van Groen T, et al. Distinct growth hormone receptor signaling modes regulate skeletal muscle development and insulin sensitivity in mice. J Clin Invest (2010) 120:4007–20. doi: 10.1172/JCI42447
38. Singh D, Sanyal S, Chattopadhyay N. The role of estrogen in bone growth and formation: changes at puberty. Cell Health Cytoskeleton (2011) 3:1. doi: 10.2147/CHC.S8916
39. Zeman M, Buyse J, Lamosova D, Herichova I, Decuypere E. Role of melatonin in the control of growth and growth hormone secretion in poultry. Domest Anim Endocrinol (1999) 17:199–207. doi: 10.1016/S0739-7240(99)00037-5
40. Munley KM, Deyoe JE, Ren CC, Demas GE. Melatonin mediates seasonal transitions in aggressive behavior and circulating androgen profiles in male Siberian hamsters. Horm Behav (2020) 117:104608. doi: 10.1016/j.yhbeh.2019.104608
41. Tan D-X, Manchester LC, Esteban-Zubero E, Zhou Z, Reiter RJ. Melatonin as a potent and inducible endogenous antioxidant: synthesis and metabolism. Molecules (2015) 20:18886–906. doi: 10.3390/molecules201018886
42. Bantounou M, Plascevic J, Galley HF. Melatonin and related compounds: Antioxidant and anti-inflammatory actions. Antioxidants (2022) 11:532. doi: 10.3390/antiox11030532
43. Vicencio JM, Galluzzi L, Tajeddine N, Ortiz C, Criollo A, Tasdemir E, et al. Senescence, apoptosis or autophagy? Gerontology (2008) 54:92–9. doi: 10.1159/000129697
44. Fernández A, Ordóñez R, Reiter RJ, González-Gallego J, Mauriz JL. Melatonin and endoplasmic reticulum stress: relation to autophagy and apoptosis. J Pineal Res (2015) 59:292–307. doi: 10.1111/jpi.12264
45. Hosseinzadeh A, Kamrava SK, Joghataei MT, Darabi R, Shakeri-Zadeh A, Shahriari M, et al. Apoptosis signaling pathways in osteoarthritis and possible protective role of melatonin. J Pineal Res (2016) 61:411–25. doi: 10.1111/jpi.12362
46. Zhao RZ, Jiang S, Zhang L, Yu ZB. Mitochondrial electron transport chain, ROS generation and uncoupling. Int J Mol Med (2019) 44:3–15. doi: 10.3892/ijmm.2019.4188
47. Pazzaglia S, Pioli C. Multifaceted role of PARP-1 in DNA repair and inflammation: pathological and therapeutic implications in cancer and non-cancer diseases. Cells (2019) 9:41. doi: 10.3390/cells9010041
48. Shamash S, Reichert F, Rotshenker S. The cytokine network of wallerian degeneration: tumor necrosis factor-α, interleukin-1α, and interleukin-1β. J Neurosci (2002) 22:3052–60. doi: 10.1523/JNEUROSCI.22-08-03052.2002
49. Zoetendal EG, Collier CT, Koike S, Mackie RI, Gaskins HR. Molecular ecological analysis of the gastrointestinal microbiota: a review. J Nutr (2004) 134:465–72. doi: 10.1093/jn/134.2.465
50. Ma N, Zhang J, Reiter RJ, Ma X. Melatonin mediates mucosal immune cells, microbial metabolism, and rhythm crosstalk: A therapeutic target to reduce intestinal inflammation. Med Res Rev (2020) 40:606–32. doi: 10.1002/med.21628
51. Yin J, Li Y, Han H, Chen S, Gao J, Liu G, et al. Melatonin reprogramming of gut microbiota improves lipid dysmetabolism in high-fat diet-fed mice. J Pineal Res (2018) 65:e12524. doi: 10.1111/jpi.12524
52. Burcelin R. Regulation of metabolism: a cross talk between gut microbiota and its human host. Physiology (2012) 27:300–7. doi: 10.1152/physiol.00023.2012
53. Pandeya DR, D’Souza R, Rahman MM, Akhter S, Kim H-J, Hong S-T. Host-microbial interaction in the mammalian intestine and their metabolic role inside. BioMed Res-India (2012) 23:9–21. doi: 10.1080/14763141.2012.671354
Keywords: melatonin, lamb, muscle, fat, intestinal microbe, apoptosis
Citation: Ma W, Wu H, Li G, Yan L, Wang L, Zhao M, Guan S, Xu S, Guo X, Liu F, Ji P, Wusiman A and Liu G (2022) Melatonin promotes the growth and development of lambs by increasing growth hormone and testosterone, targeting on apoptosis signaling pathway and intestinal microflora. Front. Endocrinol. 13:966120. doi: 10.3389/fendo.2022.966120
Received: 10 June 2022; Accepted: 25 July 2022;
Published: 19 August 2022.
Edited by:
Yanli Zhang, Nanjing Agricultural University, ChinaReviewed by:
Roland Eghoghosoa Akhigbe, Ladoke Akintola University of Technology, NigeriaJun Wang, Jilin Agriculture University, China
Copyright © 2022 Ma, Wu, Li, Yan, Wang, Zhao, Guan, Xu, Guo, Liu, Ji, Wusiman and Liu. 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: Guoshi Liu, Z3NobGl1QGNhdS5lZHUuY24=
†These authors share first authorship