Skip to main content

ORIGINAL RESEARCH article

Front. Vet. Sci., 10 July 2023
Sec. Livestock Genomics

Melatonin regulates the periodic growth of secondary hair follicles through the nuclear receptor RORα

Zeyu Lu&#x;Zeyu Lu1Jing Wu&#x;Jing Wu2Jing WuJing Wu3Tiejia ZhangTiejia Zhang4Junyang LiuJunyang Liu1Qing MuQing Mu1TerigeleTerigele1Zixian WuZixian Wu1Yanjun ZhangYanjun Zhang1Rui SuRui Su1Zhihong LiuZhihong Liu1Zhiying WangZhiying Wang1Ruijun WangRuijun Wang1Lv QiLv Qi1Yanhong Zhao
Yanhong Zhao1*
  • 1State Key Laboratory of Animal Genetics and Breeding and Reproduction, College of Animal Science, Inner Mongolia Agricultural University, Hohhot, Inner Mongolia, China
  • 2College of Food Science and Engineering, Inner Mongolia Agricultural University, Hohhot, Inner Mongolia, China
  • 3Shangdu County Vocational and Technical School, Ulanqab, Inner Mongolia, China
  • 4Zhangbei Liang Mianjing People's Government, Zhangjiakou, Hebei, China

Cashmere is the fine bottom hair produced by the secondary hair follicles of the skin. This hair is economically important. Previous studies by our research group have shown that exogenous melatonin (MT) can regulate the periodic growth of secondary hair follicles, induce the secondary development of villi, and alter the expression of some genes related to hair follicle development. Few studies on the regulation of villus growth by MT binding receptors have been published. In this study, MT was implanted subcutaneously behind the ear of Inner Mongolia cashmere goats. RT-qPCR, in situ hybridization, Western blot analysis, immunofluorescence and RNAi techniques were used to investigate the receptors and functions of MT in regulating the development of secondary hair follicles in Inner Mongolia cashmere goats. The results showed that MT binds to the nuclear receptor RORα on dermal papilla stimulates hair follicle development and promotes villus growth. The RORα mRNA expression in the skin of Inner Mongolia cashmere goats was periodic and showed a trend of first increasing and then decreasing. The expression began to increase in February, peaked in April, and reached the lowest level in May. RORα significantly affected the mRNA expression of β-catenin gene, a key gene in hair follicle development, in the presence of MT. It will lay a solid molecular foundation for further research on the regulation mechanism between MT receptor and villus growth and development and to achieve artificial regulation of villus growth time and yield to improve the effect of villus production.

1. Introduction

As an important textile material, cashmere is known as “fiber gem” and “soft gold” because of its light, soft and warm fiber and has economic value in animal husbandry. The skin of the Inner Mongolia cashmere goat is composed of cashmere and wool. The primary hair follicles develop into wool, and the secondary hair follicles develop into cashmere. The primary hair follicles generally appear at 55 days of gestational age and develop completely at approximately 135 days. Secondary hair follicles begin at 75 days of gestational age and are fully developed at 6 months after birth (1, 2). Hair follicle development is divided into three periods: the growth phase, catagen and telogen. The growth phase starts in April every year, and growth accelerates from August to September; during catagen, which starts in October, the hair stem stops growing, and epithelial cells begin to undergo apoptosis; during telogen, which occurs from December to March of the following year, hair follicle growth enters a relative resting phase (2).

The growth and development and periodic changes of hair follicles are complex processes regulated by many factors, including endocrine hormones and active molecules such as growth factors and cytokines, including melatonin (MT). MT is an indole hormone secreted mainly by the pineal gland and by some plants. Many studies have shown that MT has a wide range of functions, including immune regulation (3), treatment of animal diseases (4), antioxidant effects (5, 6), treatment of psychiatric diseases (7), use as a sleep aid (8, 9), resistance to radiation damage (10), antitumor effects (11), animal reproduction (1215), induction of hair follicle development prior to the growth phase and other functions (1619). In addition to directly acting on its receptors to alter circadian activity, MT can bind with highly chimeric sites in the pituitary to affect the secretion of nervous system hormones and then affect circadian activity. Ross et al. (20) found that the PER1 mRNA level in the pituitary of hamsters increased with the decrease in MT content in the blood.

Light intensity and sunshine duration had a major influence on MT secretion in animals, and MT content in the blood of animals changed periodically with the season. When the light duration was short, the content of MT in the blood increased and otherwise decreased. Studies have shown that MT can promote the villus growth of cashmere goats in the non-cashmere period. Since the nights are long and the days are short in winter, MT can promote the villus growth of cashmere goats and stimulate hair follicles to increase cashmere production (21, 22). MT was shown to have a regulatory effect on the periodic growth of villi of seasonal molting animals (23). Duan Chunhui (24) implanted MT in seasonal long-haired sheep in April and June and the results showed that MT could improve the wool yield and length of half-sibling sheep and non-half-sibling sheep and had a more obvious effect on promoting wool growth of half-sibling sheep. Wang Yan (25) found through qPCR that MT could activate Wnt signaling in catagen and telogen and improve the expression of Wnt10b. Through immunohistochemistry and HE staining, Mu Qing (26) found that after MT implantation, the number of secondary hair follicles in May and June was significantly higher than that in the control group, and the distance to the epidermis was closer than that in the control group, indicating that MT may activate related genes to restart and promote hair follicle development in advance and play an important role in the transition stage from telogen to the growth phase. Chang Zili (27) also showed that MT can induce the early occurrence of hair follicles, resulting in the secondary development of hair follicles, and MT downregulated the expression of Delta1, Notch1 and hairless in the Notch pathway during the growth phase, promoting villus growth.

Hormones function by binding to the corresponding receptor. MT is a highly lipophilic substance that easily crosses the cell membrane and nuclear membrane and can also be dissolved in aqueous solution; thus, it has various roles in tissue cells. Two classes of MT receptors have been identified thus far: membrane receptors (MT1, MT2, and possibly MT3) and nuclear receptors (RZR/ROR, a member of the orphan receptor subfamily). Mammalian MT receptors consist of MT1 (Mel1a), MT2 (Mel1b) and their related receptor GPR50 (28). Mel1c is the MT receptor of nonmammalian animals. Mel1a is mostly distributed in the hypothalamus, pituitary and suprachiasmatic nuclei (SCN), while Mel1b is mostly distributed in the retina (29, 30). MT secretion is increased at night, and MT binds to the G-protein-coupled receptor and activates MTNR1a, which enables the cells in the animal to detect radiant light energy and coordinate the circadian activities of cells within 24 h through the light reception of the retina, consistent with the natural photocycle (29). Mel1b has an antagonistic effect on MT and can resist the inhibitory effect due to the increase in MT concentration, which is opposite to the effect of Mel1a (3036).

The MT nuclear receptor RORα, also known as the MT ligand, shares a set of DNA domains and ligand domains with RORβ and RORγ. A previous study showed that RORα1 could be activated by MT under specific conditions, and MT was identified as the ligand of RORα. Later, Wang et al. (37) found that 7A-hydroxyl and 24S-hydroxyl cholesterol could inhibit the transcriptional activity of RORα and RORγ in liver cells. LBD binding of RORγt small molecule inhibitors was shown to inhibit the function of Th17, which further confirmed that the natural ligands of RORα and RORγt are cholesterol and its derivatives (38). Studies have shown that CS mRNA binding to SR1001 in the liver can inhibit the circadian activity of RORα and CS mRNA (39) expression, suggesting that RORα plays a key role in circadian rhythm.

Breast (40), melanoma (41), hepatocellular carcinoma (42), and colon cancers (43), etc. have been reported to be closely related to RORα. Whether RORα is involved in MT regulation of villus growth has not yet been reported. In this study, real-time fluorescence quantitative PCR, in situ hybridization, RNAi and other methods were used to preliminarily identified RORα as the receptor of MT and the role of RORα in MT regulation of villus cycle growth.

2. Materials and methods

2.1. Animals

One-year-old Inner Mongolia Albasi cashmere goats were obtained from Jinlai Animal Husbandry Technology Co., Ltd., Inner Mongolia Autonomous Region and were fed a mixture of coarse feed and granular concentrate according to the sheep farm feeding standard (400 g/d, coarse feed 50%, granular concentrate 50%). The animal management and operation procedures adopted in this study were in accordance with the guidelines of the Animal Feeding and Use Committee of Inner Mongolia Agricultural University and were carried out in strict accordance with animal welfare and ethical guidelines. The animals are shown in Figure 1.

FIGURE 1
www.frontiersin.org

Figure 1. Images of 1-year-old Inner Mongolian cashmere goats. Photo taken of Inner Mongolia cashmere goats from Inner Mongolia Jinlai Animal Husbandry Technology Co., Ltd. (A) Inner Mongolia cashmere goats, (B) Inner Mongolia cashmere goat.

2.2. Implantation of MT

The growth cycle of cashmere was 12 months. MT capsules were implanted in the subcutaneous tissue behind the ears of the cashmere goats every 2 months (provided by the Animal Science Laboratory of Northeast Forestry University), and the amount implanted was 2 mg/kg for each sheep.

2.3. Sample collection

All methods were implemented in accordance with the relevant guidelines and regulations of the European Commission (1997). The method was also approved by the Ethics Committee of Inner Mongolia Agricultural University. The animals used in this test were from Inner Mongolia Jinlai Animal Husbandry Technology Co., Ltd. Disinfect and sterilize the sampling environment and sampling tools, and the sampling personnel wear experimental clothing, masks, gloves, etc. to prevent sample contamination. One-year-old female Inner Mongolia cashmere goats with similar body weights were randomly divided into two groups: the experimental group (MT) and the control group (no preparation). Local anesthesia with procaine was used to minimize pain. Skin tissue was collected from the experimental group and the control group from the posterior edge of the shoulder every month for a year. The goats with skin injury were observed for 1 week to confirm their recovery. Each skin sample was approximately 1 cm2, part of which was transported back to the laboratory in liquid nitrogen to extract total RNA (to synthesize cDNA) and total protein from the sample. The other part was placed in 4% paraformaldehyde and transported to the laboratory for preparation of paraffin sections. A cashmere goat was randomly selected and slaughtered 4 h after sunset on the day of skin tissue collection. The head was collected, the midline skin was cut lengthwise along the nasal bone and peeled off laterally until the skull was completely exposed, and a cross cut was made in the middle of the nasal bone until the brain tissue inside the head was observed. To prevent damage to the cerebral cortex, we cut the bone on either side through a cross incision, flipped the skull cap 180° to expose brain tissue, and moved the cerebellum aside to expose the junction between brain stem and brain, allowing visualization of the hypothalamus of the cashmere goat. The hypothalamus was clipped along the root, and the hypothalamus tissue was rapidly transported back to the laboratory in liquid nitrogen to extract total RNA (for synthesis of cDNA) from the sample. Mouse skin tissue cDNA was provided by the State Key Laboratory of Animal Genetics, Breeding and Reproduction.

2.4. Primer design

Based on the MTNR1a, MTNR1b, RORα, SFRP1, TCHHL1, and β-catenin DNA sequences of other species and goat sequences in GenBank (Nut Cloud link1), the conserved cDNA sequences of other species were obtained by DNAStar homology comparison. Primers were designed according to the transitron principle, and β-actin was used to correct the expression level (Table 1).

TABLE 1
www.frontiersin.org

Table 1. Primer sequences.

2.5. In situ hybridization

In situ hybridization was performed on skin tissue sections to detect the position of complementary sequences in eukaryotic cells using RNA probes of known sequences. Based on the 550 bp MTNR1a cDNA, a 364 bp MTNR1a probe PCR primer was designed. The RORα probe template was obtained by PCR amplification and purification using the RORα plasmid as a template. The localization of MTNR1a and RORα in Inner Mongolia cashmere goat skin was studied by specifically combining a DIG-labeled nucleic acid probe with cashmere goat skin tissue according to the principle of base complementation.

2.6. Immunofluorescence staining detection of RORα protein expression in hair follicles

The skin tissue was cut into 0.15 cm3 pieces and placed in PBS containing 30% sucrose for dehydration. After the slicer was precooled, the tissue pieces were cut into 8 μm thick frozen sections and sealed with PBS buffer of 10% BSA (containing 0.4% Triton to increase permeability) at room temperature and dark conditions for 1 h. The slices were incubated with 0.02 mol/L primary antibody (Anti-ROR alpha, ab60134, Abcam, America) containing 0.4% Triton prepared with PBS at 4°C overnight and then incubated with 0.01 mol/L secondary antibody [Goat Anti-Rabbit IgG H&L (Alexa Fluor®) preabsorbed, ab150081, Abcam, America] containing 0.4% Triton prepared with PBS at room temperature for 1 h. Rinse with 0.01 mol/L PBS after each step. The sections were dried naturally and sealed for fluorescence microscopy to detect the location of RORα protein in skin hair follicles.

2.7. Western blot

Total protein was extracted and isolated from the skin tissues of goats with and without MT implantation in February, September and December and electrically imprinted onto a PVDF (polyvinylidene fluoride) membrane. After being sealed at room temperature for 1 h, the membrane was incubated with primary antibody overnight. After incubation, the membrane was washed with 15 ml of TBS/T. After incubation with secondary antibody for 1.5 h, the membrane was rinsed again, and the protein expression of RORα in skin hair follicles was detected.

2.8. Cell culture

Fibroblasts were cultured by the tissue block method. The sample was washed 3–5 times with PBS containing 3% double antibody. Tissue blocks were cut into 2 mm2, placed in a 25 mm2 culture flask, flipped until tissue block faced up, and placed in an incubator at 37°C and 5% CO2 for 4 h. After removal, the samples were turned over, complete medium was added to soak the tissue blocks, and the samples were returned to the incubator and continued to grow. The medium was changed every 3 days according to the growth state of the cells. Fibroblasts were cultured as donor cells, and the cells were inoculated in 24-well plates at a density of 5 × 103 cells/ml. The cells were counted with a blood cell counter for 8 consecutive days. The cell number in 3 wells at each time point was calculated, and the average value was used to generate the cell growth curve.

2.9. The effect of RNAi on the expression of RORα and other genes in fibroblasts

A RORα-specific shRNA vector was constructed, and the RORα target sequence was synthesized according to RORα sequence and interference fragment design principles (Table 2). An shRNA template was used for the annealing reaction, restriction enzyme digestion of the vector and target fragment was performed, the reaction product and linearized plasmid were connected and transformed into receptive cells, and positive clones were identified and resequenced. GFP-labeled RORα shRNA cells were transfected, and RORα mRNA expression was detected by RT-qPCR.

TABLE 2
www.frontiersin.org

Table 2. RORα target sequence.

2.10. Statistical analysis

CT values of the number of cycles in each group were determined by comparative CT value method in relative quantification method. mRNA expression levels were analyzed and compared by 2-(ΔΔCT) quantitative analysis method, and expressed in the form of mean ± standard error. SAS software was used for one-way ANOVA and significance analysis, P < 0.05 was considered statistically significant.

3. Results

3.1. Expression of MTNR1a, MTNR1b, and RORα in different tissues of cashmere goats

According to the above primer design principle, the specific MTNR1a cDNA band from the cashmere goat hypothalamus was amplified. No target fragment of MTNR1a cDNA was amplified from cashmere goat skin (Figure 2A). The MTNR1b fragment could not be amplified from either goat skin or hypothalamus cDNA, but specific bands were amplified from mouse skin cDNA (Figure 2B). The RORα fragment was amplified from the cDNA of both tissues (Figure 2C). From the above results, we preliminarily concluded that MTNR1a is expressed in the hypothalamus of cashmere goats but not the skin. MTNR1b was not expressed in the skin and hypothalamus of cashmere goats. RORα was expressed in the skin and hypothalamus of cashmere goats in each month.

FIGURE 2
www.frontiersin.org

Figure 2. RT-PCR electrophoretic map. (A) MTNR1a M: DL2000 DNA marker. 1: Skin amplification results. 2, 3: Results of hypothalamus amplification. 4: Negative control. (B) MTNR1b M: DL2000 DNA marker. 1: Negative control. 2: Results of mouse skin amplification. 3: Results of hypothalamus amplification in cashmere goats. 4: Results of skin amplification in cashmere goats. (C) RORα M: DL2000 DNA marker. 1–4: Results of skin amplification in cashmere goats. 5: Negative control. 6: Results of amplification of the hypothalamus from cashmere goats.

3.2. In situ hybridization was used to detect the expression of MTNR1a in different tissues of cashmere goats

The amplified sequences of the hypothalamus and skin were transcribed into MTNR1a cDNA probes, and FISH detection was performed. The results showed a strong positive signal in the hypothalamus, but no positive signal was detected in the skin (Figure 3), indicating that MTNR1a was widely expressed in the hypothalamus tissue but not in the skin tissue.

FIGURE 3
www.frontiersin.org

Figure 3. Results of in situ hybridization of the MTNR1a cDNA probe in different tissues. (A) An antisense strand probe was hybridized in situ in hypothalamus tissue. (B) The sense strand probe was hybridized in situ in hypothalamus tissue. (D) The antisense strand chain probe was hybridized in situ in skin tissue. (E) The sense strand probe was hybridized in situ in skin tissue. (C,F) Blank control. Magnification is 100×.

3.3. RORα mRNA expression changed dynamically during the villus development cycle

In skin samples from Inner Mongolia cashmere goats, the target fragments of RORα and β-actin were 417 and 208 bp, respectively, and no bands were found in the negative control. The changes in RORα mRNA expression in different months showed that RORα mRNA expression changed dynamically during the villus development cycle (Figure 4), with the highest expression in April, followed by February, and the lowest expression in May. In general, the expression level was low in the growing stage and increased in catagen and telogen, with a significant difference.

FIGURE 4
www.frontiersin.org

Figure 4. The relative expression of RORα mRNA in the skin of cashmere goats in different months. The expression of RORα mRNA changed dynamically throughout the villus development cycle. qPCR results are the control (n = 3). a, b, and c represent significant differences (p < 0.05).

3.4. In situ hybridization was used to detect RORα expression in different parts of cashmere goat skin

In situ hybridization of hair follicles from cashmere goat skin in February, April, August, October, and December and hair follicles from cashmere goat skin implanted with MT was carried out by using a RORα cDNA probe. The distribution of RORα in skin hair follicles in February during the regression period, August during the growth period and the implantation group was randomly selected (Figure 5). Strong positive expression signals were found in the two different months and the implanted group. In February, strong positive signals were detected in the hair stem. In August, strong expression was observed in the hair stem cortex and outer root sheath. In the MT group, the positive signals of hair follicles were mainly found in the outer root sheath of the hair shaft and hair shaft.

FIGURE 5
www.frontiersin.org

Figure 5. In situ hybridization of the RORα mRNA probe in skin tissue. (A) Results of in situ hybridization of the antisense strand probe in skin tissue of the MT-implanted group in February. (B) Results of in situ hybridization of the sense strand probe in skin tissue of the MT-implanted group in February. (D) Results of in situ hybridization of the antisense strand probe in skin tissue of the MT-implanted group in August. (E) Results of in situ hybridization of the sense strand probe in skin tissue of the MT-implanted group in August. (C,F) Blank control. Magnification is 40×.

3.5. Expression distribution of RORα in skin hair follicles

Immunofluorescence results showed that RORα was expressed in February, September and December of the hair follicle development cycle, and in February, September, and December RORα was expressed in telogen, the growth stage and catagen, respectively. In the MT-implanted group, RORα protein of primary hair follicles and secondary hair follicles was expressed in the external root sheath and dermal papilla of the hair follicle at three stages, while RORα protein of primary and secondary hair follicles in the control group was only expressed in the external root sheath of the hair follicle at three stages (Figures 6, 7). We concluded that the binding of MT and RORα on dermal papilla stimulates the development of hair follicles and promotes villus growth.

FIGURE 6
www.frontiersin.org

Figure 6. The expression and distribution of RORα protein in primary hair follicles in February, September and December were detected by immunofluorescence. (A–C) The expression of RORα protein in primary hair follicles of the experimental group, control group and negative control group in February. (D–F) The expression of RORα protein in primary hair follicles of the experimental group, control group and negative control group in September. (G–I) The expression of RORα protein in primary hair follicles of the experimental group, control group and negative control group in December. The magnification of (A–E,G,H) is 400×, (F) is 200×, (I) is 100×. The white arrow points to a fluorescent signal.

FIGURE 7
www.frontiersin.org

Figure 7. The expression and distribution of RORα protein in secondary hair follicles in February, September and December were detected by immunofluorescence. (A–C) The expression of RORα protein in secondary hair follicles of the experimental group, control group and negative control group in February. (D–F) The expression of RORα protein in secondary hair follicles of the experimental group, control group and negative control group in September. (G–I) The expression of RORα protein in secondary hair follicles of the experimental group, control group and negative control group in December. The magnification of (A,D,G) is 400×, (E,H) is 200×, (B,C,F,I) is 100×. The white arrow points to a fluorescent signal.

3.6. MT upregulated RORα protein expression in skin hair follicles

To further verify the effect of MT on RORα protein, we used Western blotting to detect the changes in RORα protein in February, September, and December (Figures 8A,B). The expression of RORα in the implanted group was significantly higher than that in the control group in February and showed a downward trend from September to December, with no significant difference compared with the control group (Figure 8C). This result indicated that MT can upregulate RORα expression in February.

FIGURE 8
www.frontiersin.org

Figure 8. The protein expression of RORα and β-actin in the hair follicle development period of Inner Mongolia cashmere goats was detected by Western blotting. (A) 1, 2, and 3 are RORα proteins in the implanted MT groups of February, September, and December, respectively; 4, 5, and 6 are RORα proteins of the control group in February, September, and December, respectively. (B) 1, 2, and 3 are β-actin proteins in the implanted MT groups of February, September, and December, respectively; 4, 5, and 6 are β-actin proteins of the control group in February, September, and December, respectively. (C) Relative protein expression level of RORα in hair follicles in February, September, and December. The data are presented as the means ± SEMs of triplicates from three independent experiments, using fold change compared to the control group (n = 3). The measure of β-actin was used as Intrinsic reference protein for Western blot. a, b represent significant differences (p < 0.05).

3.7. RNAi inhibited RORα expression in fibroblasts

shRNA targeting the RORα gene was successfully constructed by a chemical synthesis method. After transfection with liposomes into fibroblasts, the transfection efficiency of the RORα shRNA was detected by RT-qPCR. The relative RORα mRNA expression was determined by the 2-(ΔΔCT) method. However, the expression level of RORα in the transfected control group and the embedded fibroblasts was much lower than before transfection, and there was no significant difference between the negative control group and the implanted group (Figure 9). The results showed that MT significantly inhibited RORα.

FIGURE 9
www.frontiersin.org

Figure 9. The effect of RNA interference on RORα gene expression in Inner Mongolia cashmere goat skin fibroblasts was detected by qPCR. shRNA was loaded into Inner Mongolia cashmere goat skin fibroblasts by RNA interference technology, and the expression of RORα was detected by qPCR before and after transfection in different groups. qPCR results are presented as means ± SEMs of triplicates from three independent experiments, using fold change compared to the control (n = 3). a, b, c represent significant differences (p < 0.05).

3.8. The effect of RORα on SFRP1, TCHHL1, and β-catenin expression

The β-catenin mRNA level was not significantly different before and after transfection in the control group but was significantly different before and after transfection in the MT group. Before transfection, the β-catenin mRNA level in the MT group was significantly higher than that in the control group (Figure 10A). There was a significant difference in TCHHL1 mRNA expression between the control group and the other groups, but there was no significant difference in TCHHL1 mRNA expression between the MT group and the MT-transfected group (Figure 10B). The SFRP1 mRNA expression level did not change before and after MT treatment but decreased before and after transfection, but the difference was not significant (Figure 10C).

FIGURE 10
www.frontiersin.org

Figure 10. Effect of downregulation RORα expression on its related genes. (A) RNA interference inhibited RORα expression, and qPCR showed that MT could increase the expression of β-catenin. (B) RNA interference inhibited RORα expression, and qPCR showed that MT could lower the expression of TCHHL1. (C) RNA interference inhibited RORα expression, and qPCR showed that MT had no significant effect on SFRP1. qPCR results are presented as means ± SEMs of triplicates from three independent experiments, using fold change compared to the control (n = 3). a, b, c represent significant differences (p < 0.05).

4. Discussion

MTNR1a is expressed rhythmically in the hypothalamus, and MTNR1a mRNA expression is higher in the evening than at dawn (44). MTNR1b mRNA is expressed in different tissues and organs of different animals. Migaud et al. (45) detected the expression of MTNR1a in the hypocolliculus papilla precursor and pituitary nodules of sheep, but MTNR1b mRNA was not detected, indicating that MTNR1b does not exist in sheep. At present, most studies on MTNR1b are related to polycystic ovary syndrome (46) and gestational diabetes mellitus (47), and few studies on mammals have been published. The results of Expression of MTNR1a was detected in the hypothalamus and MTNR1b was detected in mouse skin in this study were consistent with the results of Chen et al. (48). In this study, the MTNR1a, MTNR1b and RORα mRNA expression in Inner Mongolia cashmere goat skin tissue in February, April, August, October and December was detected by RT-PCR and in situ hybridization. The positive control was the hypothalamus of Inner Mongolia cashmere goat, and MTNR1a mRNA expression was only detected in the hypothalamus tissue. In addition, mouse skin tissue was used as the positive control, and MTNR1b mRNA expression was detected only in mouse skin tissue. MTNR1b mRNA was not expressed in the skin and hypothalamus of Inner Mongolia cashmere goats. However, RORα mRNA expression was detected in the skin tissue at all 5 months.

MT plays a role in immune system regulation, lymphocyte regulation and antitumor activity by binding with RORα, and the expression of RORα mRNA was detected only in the skin tissue of cashmere goats, indicating that skin tissue is the target organ of RORα and RORα is an MT nuclear receptor. However, the expression of RORα in hair follicles is not clear. Therefore, we verified the localization of RORα in skin tissue by in situ hybridization. In the control group, RORα had a strong signal in different parts of the hair shaft, hair shaft and hair follicle of cashmere goat skin in February, April, August, October, and December. RORα expression in the implanted group was mainly found in the outer root sheath of the hair shaft and the hair shaft, indicating that RORα may be involved in periodic hair follicle growth. However, whether the change in RORα has a significant effect on the genes related to hair follicle development needs further verification. Therefore, real-time quantitative PCR was used to assess this issue.

August and October are in the long day period, and MT secretion gradually decreases; December, February and April are in the short day period, and MT secretion gradually increases. RORα mRNA expression changes dynamically with periodic changes in MT secretion during the hair follicle development cycle (49). The MT concentration has an antagonistic effect on RORα mRNA expression; that is, when exogenous MT concentration increases, RORα mRNA expression decreases (50), thus regulating the periodic growth of hair follicles. This conclusion is consistent with the results of this study. However, the RORα mRNA expression in February was not consistent with this pattern, which may be caused by RNA degradation. MT implantation promoted RORα expression at the protein level more directly than at the mRNA level. Therefore, further verification of the function of RORα protein in cashmere goat skin hair follicles is needed.

Western blotting and immunofluorescence were used to further verify the effect of MT on RORα expression. Univariate analysis of the experimental results showed that the RORα protein level in the experimental group was higher than that in the control group in February. Because the hair follicles are dormant in February, high levels of RORα protein could be stored for hair follicle redevelopment, indicating that MT implantation may promote RORα protein expression. There was no significant difference in September and December. Immunofluorescence analysis showed that the RORα protein was expressed in February, September and December of the hair follicle development cycle. The RORα protein was expressed in the outer root sheath and dermal papilla of the primary hair follicles and secondary hair follicles in the experimental group but was only expressed in the outer root sheath of the primary and secondary hair follicles in the control group. Therefore, MT could promote villus growth. RORα receptors may bind to the dermal papilla to promote villus growth and development. Hair stems were observed under a fluorescence microscope before and after immunofluorescence, and whether the RORα protein is expressed in hair stems needs further verification.

Zhao Xiaohong (51) and Kang Jian (52) downregulated the expression of RORα in MGC803 cells and fibroblasts by the RNAi technique and found that RORα expression affected the related proteins β-catenin, MMP-9, and TIMP-3 and the related gene FGF5. Xi Wenhui (53) transfected MTNR1a into fibroblasts to explore the expression of genes related to hair follicle development and found that the MT membrane receptor was not expressed in skin hair follicles, and only the nuclear receptor RORα was significantly expressed. Therefore, shRNA was linked to an interference vector in this study, and liposome transient transfection technology was used to transfect fibroblasts. The real-time fluorescence quantitative PCR results showed that the expression level of RORα increased with the decrease in MT secretion with exposure to long sunlight, whereas the expression level of RORα decreased with the increase in MT secretion, indicating that MT has a negative regulatory effect on its nuclear receptor RORα, which further verified the research of Chang Zili (27). MT regulates the periodic development of hair follicles by binding its own nuclear receptor RORα. Does RORα expression affect villus development? We screened the β-catenin, SFRP1 and TCHHL1 genes, which are related to hair follicle development, for verification.

The β-catenin is a key component of the Wnt/β-catenin signaling pathway, SFRP1 mainly acts upstream of Wnt signaling and inhibits β-catenin phosphorylation (5456). Liu Wei (57) showed that the expression of β-catenin was proportional to the hair follicle growth cycle in a study of exogenous MT implantation in cashmere goats. Loss of β-catenin expression resulted in a block of hair lamella formation and abnormal development of hair follicles, while high expression of β-catenin promoted the differentiation and development of hair follicle stem cells (58, 59). TCHHL1 is expressed in the root sheath of hair follicles, which determines the formation of hair follicle morphology and hair shape (60). In this experiment, the β-catenin mRNA level in the experimental group was significantly higher than that in the control group, and the difference was significant before and after transfection, indicating that RORα expression has an effect on β-catenin and that MT positively regulates β-catenin expression through RORα. MT did not affect the expression of SFRP1, but the SFRP1 mRNA level was not significantly decreased after transfection, indicating that the downregulation of RORα expression had almost no effect on the expression of SFRP1. The TCHHL1 mRNA level in the experimental group was significantly lower than that in the control group, and the difference was significant before and after transfection, indicating that downregulation of RORα expression has an impact on TCHHL1 and that MT regulates TCHHL1 expression through RORα. RORα regulates the coexpression of SFRP1, β-catenin and TCHHL1, We hypothesized that there may be some regulatory relationship between the three genes, which needs to be investigated.

5. Conclusion

In this study, the expression patterns and distribution of MT receptors MTNR1a, MTNR1b, and RORα in hair follicles of Inner Mongolia cashmere goats in different periods and different tissues were analyzed. The results showed that MTNRla and MTNRIb were absent in the skin tissues. In the control group, the expression signals of RORα were mainly evident in the hair shaft and outer root sheath of hair shaft, while in the MT group, the expression signals were concentrated in the dermal papilla and outer root sheath of hair shaft. It was confirmed that MT regulates the periodic growth of hair follicles by promoting the expression of nuclear receptor RORα in the dermal papilla, and also found that MT can up-regulates β-catenin and down-regulates TCHHL1 related to hair follicle development, and the expression of both genes is down-regulated after inhibiting the expression of RORα. These results provide theoretical basis and data basis for further elucidate the mechanism of MT regulation of secondary hair follicle development.

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.

Ethics statement

The animal study was reviewed and approved by the Scientific Research and Academic Ethics Committee of Inner Mongolia Agricultural University and the Biomedical Research Ethics of Inner Mongolia Agricultural University [Approval No. (2020) 056].

Author contributions

JW (third author) and TZ conceived and designed the study, and obtained, analyzed, and interpreted the data. ZeL drafted the manuscript and revised it. ZeL, JW (second author), QM, T, JL, and ZXW analyzed and interpreted the data. YHZ conceived and designed the study, received funding and oversaw the study. All authors were involved in the design and implementation of the experiment and in the writing of the manuscript.

Funding

This research was supported by grants from the National Key Research and Development Program of China (2022YED1300200); National Natural Science Foundation of China (32160772 and 31860628); Special Project for Cultivating High Level Achievements in the College of Animal Science (BZX202204); and Innovative Research Team in Universities of Inner Mongolia Autonomous Region (NMGIRT2322).

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.

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.

Supplementary material

The Supplementary material for this article can be found online at:https://www.frontiersin.org/articles/10.3389/fvets.2023.1203302/full#supplementary-material

Footnotes

References

1. Luo, JY. Transcriptional regulation mechanism of hair follicle development in yaks with seasonal changes. Master Thesis: Lanzhou University (2020).

Google Scholar

2. Ma, LN. Effects of TNFRSF1A on hair follicle development cycle of Inner Mongolia cashmere goats. Master Thesis: Inner Mongolia Agricultural University (2019).

Google Scholar

3. Xiao, ZP, Shi, Q, and Fan, ZY. Physiological function of melatonin and its effect on body rhythm (2016) 8:9–11.

Google Scholar

4. Wei, L, Wang, X, Li, QF, Liu, CG, and Bai, XJ. Analysis of the correlation between the breed, sex, age, exogenous melatonin and the incidence of trichomophagism in minks. J Northeast Agric Univ. (2011) 42:49–53.

Google Scholar

5. Duan, JF, and Cai, W. Research progress of melatonin and schizophrenia. Int J Psychiatry. (2013) 40:151–5.

Google Scholar

6. Pieri, C, Marra, M, Moroni, F, Recchioni, R, and Marcheselli, F. Melatonin: A peroxyl radical scavenger more effective than vitamin E. Life Sci. (1994) 55:PL271–6. doi: 10.1016/0024-3205(94)00666-0

CrossRef Full Text | Google Scholar

7. Huang, W, and Jiang, RH. Research Progress on the relationship between melatonin and mental illness. Journal of Practical Cardio-cerebral pulmonary vascular diseases. (2017) 25:8–11.

Google Scholar

8. Skene, DJ, Vivien-Roels, B, Sparks, DL, Hunsaker, JC, Pévet, P, Ravid, D, et al. Daily variation in the concentration of melatonin and 5-methoxytryptophol in the human pineal gland: effect of age and Alzheimer's disease. Brain Res. (1990) 528:170–4. doi: 10.1016/0006-8993(90)90214-V

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Zhong, WJ, Zheng, ZQ, Dong, MZ, Shi, SL, Liu, YL, and Zhang, RQ. Animal sleep test of melatonin. Chinese J Food Hygiene. (1998) 4

Google Scholar

10. Canyilmaz, E, Uslu, GH, Bahat, Z, Kandaz, M, Mungan, S, Haciislamoglu, E, et al. Comparison of the effects of melatonin and genistein on radiation-induced nephrotoxicity: results of an experimental study. Biomed Rep. (2016) 4:45–50. doi: 10.3892/br.2015.547

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Bai, YH, Wang, ZQ, Li, JJ, Liu, HL, and Xu, Z. Study on the antitumor activity of melatonin northwest. Pharm J. (2003) 1:13–4.

Google Scholar

12. Xiao, F. Effects of exogenous melatonin on seasonal reproduction and winter hair growth of foxes and semen quality of imported Finnish foxes. Master Thesis. China Agricultural University; (2004).

Google Scholar

13. Shi, JM, Li, HL, Shi, WQ, Zhang, HH, and Liu, GS. Effects of exogenous melatonin on growth and reproductive performance of Kunming mice. Chin J Anim Husbandry. (2008) 16:4098–20. doi: 10.11569/wcjd.v16.i36.4098

CrossRef Full Text | Google Scholar

14. Xiong, LR, Huo, LJ, Zhang, P, and Xiong, JJ. Effects of melatonin on antler production and reproduction of Cervidae. The first China deer industry development conference and the establishment conference of the deer industry branch of the China animal husbandry association. Baotou, Inner Mongolia, China (2010).

Google Scholar

15. Feng, Q, Jing, LZ, Sui, CS, Xu, W, Wang, G, and Wang, LH. Research progress on the application of melatonin in fur animal breeding. J Feed Industry. (2008) 7:52–4.

Google Scholar

16. Liu, JC, Yin, XZ, and Fang, TQ. Effects of melatonin on growth promotion and cashmere yield of Chinese cashmere goats during non-cashmere period. J Zool. (1998) 3:9–13.

Google Scholar

17. Zheng, JT, Ren, KL, Cao, L, Li, YP, Niu, XY, and Feng, GL. Effects of melatonin on growth performance and fur quality of Rex rabbits. J. Anim. Nutrition. (2017) 29:443–8.

Google Scholar

18. Niu, XY, Zheng, JT, Cao, L, Li, YP, Tang, YP, Huang, SF, et al. Effects of melatonin on growth, fur quality and slaughter performance of Rex rabbits. Agric Sci. (2015) 43:1679–83.

Google Scholar

19. Liu, NH, and Xie, QE. Effects of melatonin implantation on early maturation of mink. J Anim Husbandry Vet Commun. (2009):33.

Google Scholar

20. Ross, AW, and Morgan, PJ. The pars tuberalis as a target of the central clock. Cell Tissue Res. (2002) 309:163–71. doi: 10.1007/s00441-002-0582-3

CrossRef Full Text | Google Scholar

21. Du, H Y. Effects of exogenous melatonin on the performance of spring wool production of alpaca. Master thesis. Shanxi Agricultural University (2005)

Google Scholar

22. Cong, Y, Deng, H, Feng, Y, Chen, Q, and Sun, Y. Melatonin implantation from winter solstice could extend the cashmere growth phase effectively. Small Rumin Res. (2011) 99:48–53. doi: 10.1016/j.smallrumres.2011.03.055

CrossRef Full Text | Google Scholar

23. Wang, JM, Fu, SY, Li, C, Bai, X, Liu, B, Ji, XY, et al. The effect of melatonin on the let-7 family and its target genes in the skin of Inner Mongolia cashmere goats. J Acta Vet Zootech Sinica. (2019) 50:1367–79.

Google Scholar

24. Duan, CH. Study on the mode and mechanism of melatonin regulating cashmere growth. Doctor thesis. China Agricultural University (2016).

Google Scholar

25. Wang, Y. Effects of melatonin on the expression of genes related to hair follicle development in Inner Mongolia cashmere goats based on RNA-seq. Master Thesis, Inner Mongolia Agricultural University (2017).

Google Scholar

26. Mu, Q. Effects of melatonin on hair follicle growth cycle and Wnt10B and β-catenin gene expression in Inner Mongolia cashmere goats. Master thesis. Inner Mongolia Agricultural University (2018).

Google Scholar

27. Chang, ZL. Effects of continuous implantation of melatonin on villus growth characteristics and related genes of cashmere goats. Doctor thesis. Inner Mongolia Agricultural University (2010).

Google Scholar

28. Liu, C, Yu, XF, Teng, B, Zhao, PY, and Ren, PG. Research progress of the role of GPR50 in obesity and related diseases. J. Life Sci. (2021) 33:61–7.

Google Scholar

29. Zhao, YH. Expression of melatonin receptor gene in skin hair follicle of Inner Mongolia cashmere goat. Doctor thesis. Inner Mongolia Agricultural University; (2009).

Google Scholar

30. Jockers, R, Delagrange, P, Dubocovich, ML, Markus, RP, Renault, N, Tosini, G, et al. Update on melatonin receptors: IUPHAR review 20: melatonin receptors. Br J Pharmacol. (2016) 173:2702–25. doi: 10.1111/bph.13536

CrossRef Full Text | Google Scholar

31. Frungieri, MB, Mayerhofer, A, Zitta, K, Pignataro, OP, Calandra, RS, and Gonzalez-Calvar, SI. Direct effect of melatonin on Syrian Hamster testes: melatonin subtype 1a receptors, inhibition of androgen production, and interaction with the local corticotropin-releasing hormone system. Endocrinology. (2005) 146:1541–52. doi: 10.1210/en.2004-0990

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Soares, JM, Masana, MI, Erşahin, Ç, and Dubocovich, ML. Functional melatonin receptors in rat ovaries at various stages of the estrous cycle. J Pharmacol Exp Ther. (2003) 306:694–702. doi: 10.1124/jpet.103.049916

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Roy, D, and Belsham, DD. Melatonin receptor activation regulates GnRH gene expression and secretion in GT1–7 GnRH neurons. J Biol Chem. (2002) 277:251–8. doi: 10.1074/jbc.M108890200

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Masana, MI, Doolen, S, Ersahin, C, Al-Ghoul, WM, Duckles, SP, Dubocovich, ML, et al. MT 2 melatonin receptors are present and functional in rat caudal artery. J Pharmacol Exp Ther. (2002) 302:1295–302. doi: 10.1124/jpet.302.3.1295

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Liu, C, Weaver, DR, Jin, X, Shearman, LP, Pieschl, RL, Gribkoff, VK, et al. Molecular dissection of two distinct actions of melatonin on the suprachiasmatic circadian clock. Neuron. (1997) 19:91–102. doi: 10.1016/S0896-6273(00)80350-5

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Roberts, JEF, Wiechmann, A, and Hu, D-N. Melatonin receptors in human uveal melanocytes and melanoma cells: melatonin receptors and melanoma cells. J Pineal Res. (2002) 28:165–71. doi: 10.1034/j.1600-079X.2001.280306.x

CrossRef Full Text | Google Scholar

37. Wang, Y, Kumar, N, Solt, LA, Richardson, TI, Helvering, LM, Crumbley, C, et al. Modulation of retinoic acid receptor-related orphan receptor α and γ activity by 7-oxygenated sterol ligands. J Biol Chem. (2010) 285:5013–25. doi: 10.1074/jbc.M109.080614

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Cheng, H, Li, XF, Zhao, WP, and Li, XF. A new study on the differentiation of Th17 cells mediated by RORα and RORγ. Chin. J. Clin. Physic. (2015) 9:1413–7.

Google Scholar

39. Solt, LA, Kumar, N, Nuhant, P, Wang, Y, Lauer, JL, Liu, J, et al. Suppression of TH17 differentiation and autoimmunity by a synthetic ROR ligand. Nature. (2011) 472:491–4. doi: 10.1038/nature10075

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Xiong, G, Wang, C, Evers, BM, Zhou, BP, and Xu, R. RORα suppresses breast tumor invasion by inducing SEMA3F expression. Cancer Res. (2012) 72:1728–39. doi: 10.1158/0008-5472.CAN-11-2762

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Slominski, AT, Kim, TK, Takeda, Y, Janjetovic, Z, Brożyna, AA, Skobowiat, C, et al. RORα and ROR γ are expressed in human skin and serve as receptors for endogenously produced noncalcemic 20-hydroxy- and 20,23-dihydroxyvitamin D. FASEB J. (2014) 28:2775–89. doi: 10.1096/fj.13-242040

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Ou, Z, Shi, X, Gilroy, RK, Kirisci, L, Romkes, M, Lynch, C, et al. Regulation of the human Hydroxysteroid Sulfotransferase (SULT2A1) by RORα and RORγ and its potential relevance to human liver diseases. Mol Endocrinol. (2013) 27:106–15. doi: 10.1210/me.2012-1145

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Chen, Q, Chen, L, Zhao, R, Yang, X-D, Imran, K, and Xing, C-G. Microarray analyses reveal liver metastasis-related genes in metastatic colorectal cancer cell model. J Cancer Res Clin Oncol. (2013) 139:1169–78. doi: 10.1007/s00432-013-1424-2

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Shao, SC, Yan, BB, Zhang, TW, Liang, P, Zhou, SF, Li, XH, et al. Association analysis of MTNR1A gene polymorphism and birth number in five sheep populations. J Inner Mongolia Agric Univ. (2021) 42:1–7.

Google Scholar

45. Migaud, M, Daveau, A, and Malpaux, B. MTNR1A melatonin receptors in the ovine premammillary hypothalamus: day-night variation in the expression of the Transcripts1. Biol Reprod. (2005) 72:393–8. doi: 10.1095/biolreprod.104.030064

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Yang, JF, Yang, HM, and Cheng, ZJ. The relationship between single nucleotide polymorphism of melatonin receptor 1B gene and polycystic ovary syndrome. J Appl Gynecol Endocrinol. (2016) 3:145–8.

Google Scholar

47. Guo, QL, Wang, R, Li, WL, and Suo, DM. The relationship between mTNR1B rs10830962 and CDKAL1 rs7754840 gene polymorphism and the incidence of gestational diabetes mellitus. J Appl Med. (2019) 35:2848–52.

Google Scholar

48. Chen, JN, Zhang, QW, Wang, Q, Wang, YY, Ma, YJ, Zhang, Y, et al. Expression of MTNR1a and MTNR1b in different tissues of yak and its localization in testis. J Vet Sci. (2020) 40:64–70.

Google Scholar

49. Wang, LF, Lu, DX, Sun, HZ, Zhao, XY, and Shan, D. Effects of light and melatonin on nitrogen allocation and cashmere production performance of Inner Mongolia cashmere goats. J Scientia Agric Sinica. (2006) 5:1004–10.

Google Scholar

50. Kang, J, and Feng, Y. Effect of maternal melatonin on expression of RORα gene in visceral tissue of offspring Rex rabbits %J. Chin J Anim Husbandry. (2020) 56:86–94.

Google Scholar

51. Zhao, XH. DADS up-regulation of RORA inhiBITS proliferation, migration and invasion of human gastric cancer MGC803 cells. Master thesis. University of South China; (2012)

Google Scholar

52. Kang, J. Study on specific expression and RNA interference of FGF5 gene in cashmere goats. Master thesis. Northwest A&F University; (2014)

Google Scholar

53. Xi, WH. Study on transfection of MT1 gene into immortalized skin fibroblasts of Inner Mongolia cashmere goat. Master thesis. Inner Mongolia Agricultural University; (2011).

Google Scholar

54. Zhu, KR, Sun, QF, and Zhang, YQ. Long non-coding RNA LINP1 induces tumorigenesis of Wilms’ tumor by affecting Wnt/β-catenin signaling pathway. Eur Rev Med Pharmacol Sci. (2019) 23:5691–8. doi: 10.26355/eurrev_201907_18306

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Lin, L, Liu, GL, Yang, LN, and Gao, L. Crocin improves spatial memory in Alzheimer's disease rats based on Wnt/β-catenin signaling pathway. Chin J Gerontol. (2021) 41:153–7.

Google Scholar

56. Ma, JW, Wang, J, and Song, W. Role and mechanism of SFRP1 and β-catenin in esophageal cancer %. Clinical laboratory medicine. (2022) 19:1326–9.

Google Scholar

57. Liu, W. Study on expression of key genes of Wnt pathway in one cycle of cashmere goat villus growth. Master thesis. Inner Mongolia Agricultural University (2010).

Google Scholar

58. Eastman, Q, and Grosschedl, R. Regulation of LEF-1/TCF transcription factors by Wnt and other signals. Curr Opin Cell Biol. (1999) 11:233–40. doi: 10.1016/S0955-0674(99)80031-3

PubMed Abstract | CrossRef Full Text | Google Scholar

59. Liu, GY, Bai, LY, Li, FC, Sun, HT, Yang, LP, Jiang, WX, et al. Research progress in the regulatory signaling pathways of hair follicle development and periodic growth. J Anim Husb Vet Med. (2021) 53:125–9.

Google Scholar

60. Yamakoshi, T, Makino, T, Ur Rehman, M, Yoshihisa, Y, Sugimori, M, and Shimizu, T. Trichohyalin-like 1 protein, a member of fused S100 proteins, is expressed in normal and pathologic human skin. Biochem Biophys Res Commun. (2013) 432:66–72. doi: 10.1016/j.bbrc.2013.01.084

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: melatonin, RORα, secondary hair follicle, hair follicles grow periodically, Inner Mongolia cashmere goat

Citation: Lu Z, Wu J, Wu J, Zhang T, Liu J, Mu Q,, Wu Z, Zhang Y, Su R, Liu Z, Wang Z, Wang R, Qi L and Zhao Y (2023) Melatonin regulates the periodic growth of secondary hair follicles through the nuclear receptor RORα. Front. Vet. Sci. 10:1203302. doi: 10.3389/fvets.2023.1203302

Received: 10 April 2023; Accepted: 22 June 2023;
Published: 10 July 2023.

Edited by:

Ran Di, Chinese Academy of Agricultural Sciences, China

Reviewed by:

Wenlin Bai, Shenyang Agricultural University, China
Cecilia Dall'Aglio, University of Perugia, Italy

Copyright © 2023 Lu, Wu, Wu, Zhang, Liu, Mu, Terigele, Wu, Zhang, Su, Liu, Wang, Wang, Qi and Zhao. 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: Yanhong Zhao, MTM5NDcxOTY0MzJAMTYzLmNvbQ==

These authors have contributed equally to this work and share first authorship

Disclaimer: 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.