Donghua Li1
Xinlei Wang1
Yawei Fu1Chenxi Zhang1Yanfang Cao1Jie Wang1Yanhua Zhang1Yuanfang Li1Yi Chen1
Zhuanjian Li1
Wenting Li1
Ruirui Jiang1,2
Guirong Sun1,2Yadong Tian1,2
Guoxi Li1,2*
Xiangtao Kang1,2*- 1College of Animal Science and Veterinary Medicine, Henan Agricultural University, Zhengzhou, China
- 2Henan Innovative Engineering Research Center of Poultry Germplasm Resource, Henan Agricultural University, Zhengzhou, China
The Xichuan black-bone chicken, which is a rare local chicken species in China, is an important genetic resource of black-bone chickens. Tyrosine can affect melanin production, but the molecular mechanism underlying tyrosine-induced melanin deposition in Xichuan black-bone chickens is poorly understood. Here, the blackness degree and melanin content of the breast muscle of Xichuan black-bone chickens fed a basic diet with five levels of added tyrosine (i.e., 0.2, 0.4, 0.6, 0.8, and 1.0%; these groups were denoted test groups I-V, respectively) were assessed, and the results showed that 0.8% tyrosine was the optimal level of added tyrosine. Moreover, the effects of tyrosine supplementation on the proliferation and tyrosinase content of melanocytes in Xichuan black-bone chickens were evaluated. The results revealed a dose-dependent relationship between tyrosine supplementation and melanocyte proliferation. In addition, 417 differentially expressed genes (DEGs), including 160 upregulated genes and 257 downregulated genes, were identified in a comparative analysis of the transcriptome profiles constructed using the pooled total RNA from breast muscle tissues of the control group and test group IV, respectively (fold change ≥2.0, P < 0.05). These DEGs were mainly involved in melanogenesis, the calcium signaling pathway, the Wnt signaling pathway, the mTOR signaling pathway, and vascular smooth muscle contraction. The pathway analysis of the DEGs identified some key genes associated with pigmentation, such as DCT and EDNRB2. In summary, the melanin content of breast muscle could be markedly enhanced by adding an appropriate amount of tyrosine to the diet of Xichuan black-bone chickens, and the EDNRB2-mediated molecular regulatory network could play a key role in the biological process of tyrosine-induced melanin deposition. These results have deepened the understanding of the molecular regulatory mechanism of melanin deposition in black-bone chickens and provide a basis for the regulation of nutrition and genetic breeding associated with melanin deposition in Xichuan black-bone chickens.
Introduction
The black-bone chicken is an important resource in poultry production, and the fact that its body contains melanin significantly differentiates it from other chicken breeds (Tian et al., 2007; Tu et al., 2009b; Yu et al., 2018a). Many studies have shown that black-bone chickens have medical benefits, such as antioxidant activity (Xin et al., 2009) and the abilities to delay aging (Xu et al., 1999), treat anemia (Xie et al., 2009) and cure female menstrual abnormalities (Tu et al., 2009a). Due to their unique edible and medicinal value, black-bone chickens are highly preferable by consumers, and the exploitation and genetic breeding of black-bone chicken resources have become increasingly important. It is known that the melanin content in the meat, skin and bone of black-bone chickens is closely related to its product value and medicinal effects (Yu et al., 2018b). Therefore, the improvement of melanin deposition by nutritional regulation and genetic breeding is of great significance for the exploitation and utilization of black-bone chicken resources; however, the relevant research remains scarce.
Melanin is a type of high-protein molecule synthesized from melanin bodies in melanocytes (Diffey et al., 1995; Mackintosh, 2001), and studies have shown that nutritional factors can affect the deposition of melanin in the body (Ralli and Graef, 1943; Li et al., 2011). The ingestion of nutrients as a signal induces transcription, RNA processing and stability, protein synthesis and modification, which are processes that affect DNA replication, regulate gene expression, and maintain cell proliferation, differentiation and adaptation (Walker and Blackburn, 2004; Ruemmele and Garnier-Lengliné, 2012). Therefore, the addition of exogenous nutrition can effectively regulate the deposition of melanin in animals. In black-bone chickens, melanin is produced by tyrosinase through a series of biochemical reactions (Sato et al., 2007). Tyrosine, a type of aromatic amino acid, is a precursor of melanin synthesis (Scheiner et al., 2002), and black-bone chickens might have a precise synergistic regulatory mechanism between tyrosine metabolism and melanin synthesis that regulates the synthesis and deposition of melanin. The elucidation of the above mechanism is of great significance for improving the meat quality of black-bone chickens through the regulation of nutrition and for identifying the molecular markers related to melanin deposition in the breeding of black-bone chickens. However, the molecular mechanism underlying tyrosine-induced melanin deposition in black-bone chickens remains poorly understood.
The Xichuan black-bone chicken, which is a rare local chicken species, is mainly distributed in Xichuan County, Henan Province, China. This variety has five black parts (beak, skin, bones, legs, and meat) and is an important genetic resource of black-bone chickens. Because its meat is tender and delicious and has high nutritional and medicinal value, this variety has become an important object of production and development in recent years. Increased production performance and the production of high-quality poultry meat with a high content of melanin through nutritional regulation or genetic improvement are currently the main goal of the resource exploitation of Xichuan black-bone chickens. The breast muscle is an important part of skeletal muscle and the main source of poultry meat. Therefore, in this study, the breast muscle of Xichuan black-bone chicken was used as the research material to evaluate the effects of different concentrations of tyrosine on melanin deposition and to analyze the transcriptome profiles of breast muscle under tyrosine supplementation. We also verified the effect of tyrosine on the proliferation and tyrosinase content of melanocytes and determined the key genes and molecular regulatory network associated with tyrosine-induced melanin deposition. The main objectives of this study were to optimize the optimal supplementary amount of tyrosine in feed for the promotion of melanin deposition in the breast muscle of Xichuan black-bone chickens and to clarify the potential molecular mechanism underlying tyrosine-induced melanin deposition in black-bone chickens. The results obtained in this study will provide a better understanding of the molecular regulatory mechanisms of melanin deposition in black-bone chickens and contribute to nutrient regulation and genetic breeding related to melanin deposition in Xichuan black-bone chickens.
Materials and Methods
Ethics Statement
All sample collections and treatments were conducted strictly in accordance with the protocol approved by the Institutional Animal Care and Use Committee (IACUC) of Henan Agricultural University, China (11-0099).
Experimental Animals and Design
A total of 480 one-day-old healthy and weight-matched Xichuan black-bone chickens were randomly divided into six groups (the control group and test groups I-V), with five biological replicates of each group and 16 chickens per biological replicate. The control group was fed a basal diet (Supplementary Table S1). For test groups I-V, the basal diet was supplemented with 0.2, 0.4, 0.6, 0.8, and 1.0% tyrosine, respectively. Tyrosine (99%) was obtained from Swire Coca-Cola Drinks Co., Ltd., (Zhengzhou, Henan Province, China), and the trial period lasted 12 weeks. The chickens in all the treatment groups were fed and managed under the same conditions and were allowed free access to feed and water during the experimental period. The rooms were regularly cleaned and disinfected, and the chickens were routinely immunized.
Sample Collection
At the end of the feeding trial, all the experimental chickens were fasted for 12 h with continued access to water. Two Xichuan black-bone chickens from each replicate (a total of 60 chickens) were randomly selected. The live weight and dressing percentage of the 60 chickens were measured. These experimental chickens were euthanized by intravenous injection of KCl (1–2 mg/kg) under deep anesthesia, and 1–2 g of breast meat was frozen immediately in liquid nitrogen and stored at -80°C until RNA-seq analysis. The remaining breast meat was stripped from both sides for analysis of the melanin content.
Evaluation of Melanin Deposition in Breast Muscle
First, the breast muscle color of each experimental chicken was determined with a portable NR10QC colorimeter (3nh, China). The color parameter was measured based on the lightness index (L∗). All color values were obtained from three areas of each breast muscle. A higher reading indicated a lower degree of blackness, and vice versa. Second, the melanin content in the breast muscle of each experimental chicken was analyzed. Briefly, 50-g breast muscle samples from each chicken were weighed, and after the fascia, membrane and fat were removed, the samples were heated in an oven at 65°C for 24 h. The crushed samples were hydrolyzed with papain under neutral conditions for 3 h during heating at 55°C in a water bath, and after centrifugation, the precipitates were separated, soaked in 200 mL of 6 mol/L HCl and heated in an electric furnace for 30 min. The residue was collected, washed, wrapped with filter paper and placed in a Soxhlet extraction tube. Subsequently, the residue was degreased with ether in a water bath at 42°C, and the filter paper was washed repeatedly with distilled water and dried at 80°C in an oven. The melanin was removed and weighed using a microbalance. The melanin content (as the percentage of fresh tissue weight) was then calculated as follows: melanin weight (g)/sample weight (g) × 100%.
Exposure of Chicken Melanocytes to Tyrosine and Evaluation of the Effects of Tyrosine Exposure
To evaluate the effect of tyrosine on melanin deposition at the cellular level, we exposed melanocytes from Xichuan black-bone chickens to different concentrations of tyrosine. The initial density of the cells was 1 × 104 cells/mL. The control group was cultured with normal melanocyte medium, and the experimental groups were cultured with medium containing different concentrations of tyrosine (10-9, 10-8, 10-7, and 10-6). Fourth-generation melanocytes at the logarithmic phase of growth were digested with 0.25% trypsin (Gibco Company, United States) and seeded into cell culture plates for follow-up testing. On the one hand, the effect of tyrosine on the proliferation of melanocytes was evaluated by the Cell Counting Kit-8 (CCK-8) method on 96-well cell culture plates. Ten replicate wells were used for each group. Ten microliters of CCK-8 solution was added to each well, and the plates were incubated for 2 h at 37°C in an atmosphere with 5% CO2. The absorbance of each well at 450 nm was measured using a high-throughput multifunctional microplate test system. The proliferation rate was calculated as follows: [(A450 (with tyrosine) – A450 (medium only)]/[A450 (without tyrosine) – A450 (medium only)] × 100%. On the other hand, the effect of tyrosine on the tyrosinase content in melanocytes was detected using a chicken tyrosinase ELISA kit on six-well plates. Each group was tested in six replicate wells in triplicate. The plates were incubated for 72 h in a humidified 5% CO2 incubator at 37°C, and the intracellular components were then extracted using an ultrasonic cell disrupter. The content of tyrosinase in the cells was detected in strict accordance with the instructions provided with the tyrosinase ELISA kit. Moreover, the effects of tyrosine on endothelin receptor B subtype 2 (EDNRB2) gene expression in the above-treated melanocytes were also determined by qRT-PCR analysis.
Transcriptome Sequencing and Annotation of Breast Muscle
Based on the results of the melanin determination, the highest amount of melanin deposition in breast muscle tissue was detected in test group IV (0.8% dietary tyrosine level). Therefore, breast muscle samples from the control group and test group IV were selected for transcriptome sequencing. The total RNA from the breast muscle samples was extracted using an RNA extraction kit (Takara, Dalian, China) following the appropriate procedure and were processed using DNase I (RNase Free, QIAGEN). To characterize a general overview of genes expressed in the breast muscle tissues of the highest amount of melanin deposition, the pooled total RNA from three individuals of test group IV and the control group was used for cDNA library construction, respectively. The libraries were generated using the NEBNext® UltraTM Directional RNA Library Prep Kit for Illumina® (NEB, Ipswich, MA, United States). The processes necessary for library construction, such as mRNA isolation, fragmentation, first-strand cDNA synthesis, second-strand cDNA synthesis, terminal repair, 3’ end A-tailing, ligation, and enrichment, were completed according to the manufacturers’ instructions. Sequencing was performed according to the corresponding requirements with a paired-ending cDNA sequencing program. The sequencing process was controlled using the Illumina HiSeq 2000 platform data collection software program that was used for the real-time data analysis. The raw reads were processed to remove adaptor sequences, low-quality reads and reads containing poly-N sequences with Seqtk1. Using the spliced mapping algorithm in TopHat (version: 2.0.9) (Trapnell et al., 2012), the clean reads were mapped to the Ensembl Galgal5 reference genome. The number of fragments per gene was determined using StringTie (version: 1.3.0), and the values were normalized by the trimmed mean of M-values (TMM) method (Robinson and Oshlack, 2010). The fragments per kilobase of transcript per million mapped reads (FPKM) value of each gene was calculated using a Perl script (Poelchau et al., 2014).
Differential Expression and Functional Enrichment Analysis
EdgeR software was used to identify the differentially expressed genes (DEGs) between the control group and test group IV (Haydock et al., 2015). The threshold P-value was determined by controlling the false discovery rate (FDR), and the corrected p-value was the q-value (Kim and van de Wiel, 2008). The fold changes in the gene expression level were calculated based on the FPKM values. The screening conditions for DEGs were a q-value ≤0.05 and a fold change ≥2. The functional classifications of the DEGs were analyzed by Gene Ontology (GO) functional enrichment analysis (Young et al., 2010) and KOBAS software (Mao et al., 2005). The GO terms and KEGG pathways with P-values <0.05 were considered significantly enriched. The genes in the KEGG pathways related to pigmentation were submitted to STRINGv10.0 for protein-protein interaction (PPI) network analysis (Franceschini et al., 2013), and the PPI network was visualized using Cytoscape (version 3.6.0) (Shannon et al., 2003).
Quantitative Real-Time PCR (qRT-PCR) Analysis
To verify the accuracy of the RNA-seq results, seven genes were selected for qRT-PCR validation. The gene-specific primers were designed using the NCBI Primer BLAST online program (Supplementary Table S2) and synthesized by Nanjing Kingsray Biotechnology Co., Ltd. The qRT-PCR analysis was performed using a PrimeScriptTM RT Reagent Kit and SYBR® Premix Ex Taq II (Tli RNase H Plus; Takara, Dalian, China). The PCRs were conducted with an initial denaturation step at 95°C for 5 min followed by 34 cycles of 95°C for 15 s (denaturation), 60°C for 45 s (annealing), and 72°C for 40 s (extension) and a final extension at 72°C for 10 min. The reactions for each qRT-PCR were performed in triplicate. The relative expression levels of the genes were calculated using the 2-ΔΔCt method, and the levels were normalized to those of two housekeeping genes, β-actin and GAPDH (Vandesompele et al., 2002).
Statistical Analysis
All the data were statistically analyzed with SPSS 20.0 software to determine the significance of the differences (P < 0.05) between the test groups and the control group. The results are expressed as the mean ± standard deviation (SD).
Results
Concentration-Dependent Effect of Dietary Tyrosine on the Breast Muscle Melanin Content
The results indicate that dietary supplementation with different doses of tyrosine significantly improved the breast muscle blackness (Figure 1A) and increased the melanin content in the breast muscle of Xichuan black-bone chickens (Figure 1B). The lowest melanin content in breast muscle was detected in the control group, and a slightly increased content was found in test group I, but no significant differences were found between test group I and the control group (P > 0.05). However, the melanin contents of the test groups II-V were significantly higher than those of the control group and test group I (P < 0.05), and the maximal value was obtained with 0.8% tyrosine supplementation (Figure 1B). Although the melanin content was also elevated by 1.0% tyrosine supplementation, the content obtained with this level of tyrosine supplementation was lower than that obtained with 0.8% tyrosine supplementation. Thus, the 0.8% dose of dietary tyrosine was identified as the maximal effective dose for promoting the melanin content in breast muscle. In addition, a correlation analysis showed a significant correlation between the breast muscle darkness and the melanin content (P < 0.05) (Figure 1C).
Figure 1. Degree of blackness and melanin content of breast muscle tissue. (A) Breast muscle blackness. (B) Melanin content of breast muscle tissue. (C) Analysis of the correlation between the melanin content and the degree of blackness of breast muscle tissue. In x-axis, the C represents the control group, and the TI-TV represents test groups I-V, respectively. Data are expressed as mean ± SD (n = 60). The small letters a, b, etc., on the top of each bar represents statistical significance. Data with different small letters on each bar are statistically significant (P < 0.05), while data with the same letters are not statistically significant (P > 0.05).
Effect of Tyrosine Exposure on the Proliferation and Tyrosinase Content of Melanocytes
We also evaluated the effect of tyrosine in chicken melanocytes. Melanocytes were cultured with different concentrations of tyrosine for 72 h, and the melanocyte morphology was then observed under an inverted microscope. The test groups treated with different concentrations of tyrosine showed significantly higher numbers of cells and evident cytoplasm enlargement compared with the control group. In particular, the most obvious cytoplasmic change was detected in the group treated with 10-9 mol/L tyrosine (Figure 2). The analysis of the melanocyte proliferation rates also showed that the group treated with 10-9 mol/L tyrosine exhibited the maximal proliferation effect (P < 0.05) (Figure 3A). Moreover, the intracellular tyrosinase content was also measured using a microplate reader (Figure 3B), and the results showed that the tyrosinase content in the group treated with 10-9 mol/L tyrosine was significantly higher than that in the control group and that the maximal proliferation effect was obtained in the group treated with 10-6 mol/L tyrosine (P < 0.05). These results indicated that supplementation with an appropriate concentration of tyrosine can significantly increase the proliferation and tyrosinase content of melanocytes of Xichuan black-bone chickens.
Figure 2. The morphology of melanocytes treated with different concentrations of tyrosine for 72 h were observed by inverted phase contrast microscope. The A indicates the control melanocytes. The B, C, and D indicate the melanocytes cultured with 10-9, 10-8, and 10-6 mol/L tyrosine, respectively. The number of cells and cytoplasm in the experimental group increased significantly, especially in the 10-9mol/L tyrosine group.
Figure 3. Effects of treatment with different concentrations of tyrosine for 72 h on melanocytes. (A) Effects of different concentrations of tyrosine on melanocyte proliferation. (B) Effects of different concentrations of tyrosine on the tyrosinase content of melanocytes. Data are expressed as mean ± SD (n = 10). The small letters a, b, etc., on the top of each bar represents statistical significance. Data with different small letters on each bar are statistically significant (P < 0.05), while data with the same letters are not statistically significant (P > 0.05).
Characteristics of the Transcriptome Profiles of Breast Muscle With the Highest Content of Melanin
The above-described studies demonstrated that tyrosine can promote melanin deposition in the breast muscle of Xichuan black-bone chickens at the cellular level and in vivo. To further reveal the molecular regulatory mechanism underlying this biological process, we used the pooled total RNA from breast muscle tissues of the control group and test group IV (0.8% dietary tyrosine level and the highest melanin content) and performed transcriptome sequencing using the Illumina HiSeq 2000 platform. A total of 47,038,307 (95.34%) and 51,161,382 (95.46%) clean reads were obtained for the control group and test group IV, respectively (Table 1). After assembly, 15,074 (94.4%) reads were obtained for the control group, 15,173 (95.02%) reads were found for test group IV, and 14,277 reads were commonly expressed in both groups. These sequencing data have been submitted to the Genome Expression Omnibus (Accession Number GSE128028) of the National Center for Biotechnology Information (NCBI). The mapping rates for the reads that were uniquely mapped to the Galgal5 assembly of the chicken genome were 84.01 and 85.01% for the control group and test group IV, respectively. The comparison of the results aligned to their respective chromosomes showed no significant differences between the samples (Figure 4A). Most of these clean reads mapped to genes and coding regions (Figure 4B), consistent with the results of previous studies (Li H. et al., 2015; Zhang et al., 2017).
Table 1. Statistics for the cDNA library sequences.
Figure 4. Distribution of the identified mRNAs on each chromosome. (A) The outer ring represents the chicken genome and is labeled with the chromosome numbers and positions. The red circle shows the distribution of the mRNAs identified in the control group, and the orange circle shows the distribution of the mRNAs identified in test group IV. (B) Mapping region distribution.
The correlation analysis showed that the gene expression data were highly correlated between the two groups, with a correlation coefficient (R2) of 0.974 (Figure 5A), which indicated the rationality of the sample selection in this study. Among the 15,968 genes identified in the breast muscle of Xichuan black-bone chickens, 14,277 genes were coexpressed in the two groups, 798 genes were specifically expressed in the control group, and 896 genes were specifically expressed in test group IV (Figure 5B). Compared with the control group, 417 DEGs with a fold-change ≥2 and q-value ≤0.05 (Figure 5D) were detected in test group IV (Figure 5C and Supplementary Table S3), and these DEGs included 160 upregulated genes and 257 downregulated genes. Moreover, three downregulated (MITF, TYR, and EDNRB2) and four upregulated (ACACB, ELOVL6, ABRA, and GPX2) genes were selected for qRT-PCR verification. The log2 (fold change) values of the seven genes obtained from the qRT-PCR analysis were consistent with the RNA-seq results (Figure 6A). An analysis of the correlation between the RNA-seq and qRT-PCR results yielded a correlation coefficient of 0.95, which indicated the reliability of the sequencing results (Figure 6B).
Figure 5. RNA-seq transcriptome profiles distinguishing the two groups. (A) Correlations of the gene expression levels among the samples. (B) Number of genes identified in breast muscle tissue from chickens in the different groups. (C) Volcano plot of differentially expressed mRNAs between the control group and test group IV. (D) Total number of differentially expressed mRNAs between the control group and test group IV. Red represents gene upregulation, and blue represents gene downregulation.
Figure 6. Validation of the RNA-seq data by qRT-PCR. (A) Comparison of the fold change between RNA-seq and qRT-PCR. (B) Pearson correlation scattered plots based on the results of RNA-seq and qRT-PCR.
A functional enrichment analysis of these 417 DEGs was then performed to identify their physiological functions. The GO enrichment analysis revealed 24 annotations for biological process terms, 16 annotations for cellular component terms and 10 annotations for molecular function terms (Figure 7 and Supplementary Table S4). In addition, these DEGs were enriched in 91 KEGG pathways (Figure 8 and Supplementary Table S5), and these pathways can be divided into six types, namely, cellular processes, environmental information processing, genetic information processing, human diseases, metabolism, and organismal systems.
Figure 7. Gene Ontology analysis of differentially expressed mRNAs. (A) Comparison of GO annotations of functional genes. The horizontal axis shows the secondary nodes of three GO categories. The vertical axis displays the number of annotated differentially expressed genes. (B) Scatter plot of the top 30 enriched GO terms identified from the differentially expressed mRNAs in breast tissue from the control group and test group IV.
Figure 8. KEGG analysis of differentially expressed mRNAs. (A) Comparison of KEGG annotations of functional genes. (B) Scatter plot of the top 30 enriched pathways identified from the differentially expressed mRNAs in breast tissue from the control group and test group IV.
Key Genes in Tyrosine-Induced Melanin Deposition
Based on the results of the KEGG pathway enrichment analysis of the DEGs, some melanin synthesis-related pathways, including melanogenesis, the calcium signaling pathway, the Wnt signaling pathway, the mTOR signaling pathway, vascular smooth muscle contraction and adrenergic signaling in cardiomyocytes, were screened (Figure 9). Some key genes related to melanin synthesis were identified from the above pathways, and these were closely associated with pigmentation (Table 2). Among these genes, DCT and EDNRB2 were identified as genes that participate in many pathways, such as melanogenesis and the calcium signaling pathway. Thus, we focused on the effect of tyrosine on EDNRB2 gene expression in chicken melanocytes. Under the experimental treatment conditions with different concentrations of tyrosine, the expression of the EDNRB2 gene was significantly upregulated in a dose-dependent manner in melanocytes of Xichuan black-bone chickens (Figure 10). In particular, the highest EDNRB2 mRNA expression level was obtained in the group treated with 10-9 mol/L tyrosine, and this expression level was significantly different compared with the control group (P < 0.05). It is known that the EDNRB2 gene activates the downstream signaling molecule PKC, upregulates the downstream gene DCT and thereby affects melanin deposition.
Figure 9. Differentially expressed mRNAs involved in the melanin production pathway. The red nodes represent upregulated genes, and the closer to red indicate that the degree of the up-regulation is higher. The blue nodes represent downregulated genes, and the closer to blue indicate that the degree of the down-regulation is higher.
Table 2. Differentially expressed genes related to pigmentation.
Figure 10. Expression level of the EDNRB2 gene (determined by qRT-PCR analysis) in melanocytes of Xichuan black-bone chickens treated with different concentrations of tyrosine. Data are expressed as mean ± SD (n = 6). The small letters on the top of each bar represents statistical significance. Data with different small letters on each bar are statistically significant (P < 0.05), while data with the same letters are not statistically significant (P > 0.05).
Discussion
Animal phenotypes are influenced by both genetics and the environment. The most obvious difference between black-bone chickens and ordinary chickens is that the skin, muscle, periosteum and other tissues of black-bone chickens contain high levels of melanin (Hirano, 1990; Tian et al., 2007). Therefore, we first investigated the effects of six levels of tyrosine supplementation on melanin deposition in the breast muscles of Xichuan black-bone chickens. The results showed that the effects of dietary tyrosine on the breast muscle melanin content are concentration-dependent. In particular, the basal diet with 0.8% tyrosine significantly improved the blackness degree of breast muscle tissue and yielded the highest melanin content in breast muscle. Previous studies have shown that diets with low levels of tyrosine turn black cat hair a reddish-brown color (Yu et al., 2001), whereas dietary supplementation with sufficient tyrosine can prevent this change (Morris et al., 2002). Our study showed that appropriate tyrosine levels can promote melanin synthesis, which is consistent with the results of previous studies (Tsatmali et al., 2002).
Melanin is a type of high-protein molecule synthesized from melanin bodies in melanocytes (Diffey et al., 1995; Mackintosh, 2001). Black-bone chicken melanin is produced by TYR through a series of biochemical reactions (Sato et al., 2007). Thus, we also evaluated the effects of tyrosine exposure on the proliferation and tyrosinase content of melanocyte. The results indicated that treatment with different concentrations of tyrosine could induce the proliferation of chicken melanocytes cultured in vitro. The proliferation of chicken melanocytes increased with increases in the tyrosine concentrations within a certain concentration range, but tyrosine concentrations that exceed a certain range inhibited melanocyte proliferation to a certain degree. We also found that tyrosine dose-dependently affects the tyrosinase content of chicken melanocytes. Specifically, the tyrosinase content of melanocytes was significantly upregulated in the tyrosine-treated groups compared with the control group, and higher levels of tyrosine supplementation increased the intracellular tyrosinase content. Analogously, previous studies have also demonstrated that increases in the tyrosine concentrations can promote melanin synthesis in hamster melanoma cells, and a certain dose-dependent relationship has been found between tyrosine and tyrosinase (Słominski et al., 1988). These studies have indicated that the meat quality, particularly the blackness degree and melanin content, could be improved by the addition of an appropriate amount of tyrosine in Xichuan black-bone chicken farming.
Moreover, we highlighted the potential molecular mechanism underlying tyrosine-induced melanin deposition in the breast muscle of Xichuan black-bone chickens. Some key genes involved in tyrosine-induced melanin deposition, such as premelanosome (PMEL), melanophilin (MLPH), EDNRB2, DCT, and other pigment-related genes, were identified based on the transcriptome profiles constructed using the pooled total RNA from the breast muscle tissue of the control group and experimental group treated with 0.8% tyrosine supplementation. Therefore, PMEL, which is also known as gp100, ME20, PMEL17, and silver, is a melanoma-specific glycoprotein that plays an important role in the development of melanin bodies by forming a proteolytic fibrous matrix for melanin deposition (Leonhardt et al., 2013). MLPH, a small GTP-binding protein (Rab27a) and myosin-Va (Myosin-Va, Myo5a) form ternary complexes that play a key role in the transport of melanosomes (Hume et al., 2007). These complexes aggregate at the dendritic terminals of melanocytes and transfer melanosomes from melanocytes to adjacent keratinocytes. Mutations in the MLPH gene might affect the transport of melanosomes. A large number of studies have shown that variations in the MLPH gene affect animal hair color and skin color development (Hume et al., 2006; Bed’hom et al., 2012; Lehner et al., 2013; Bauer et al., 2018). Therefore, these genes might play a key role in tyrosine-induced melanin deposition in the breast muscle of Xichuan black-bone chickens.
Among the identified key genes, EDNRB2 is a paralog of the endothelin receptor B (EDNRB) gene and encodes a G protein-coupled receptor (GPCR) with seven transmembrane domains (Lecoin et al., 1998; Pla et al., 2005; Harris et al., 2008; Braasch and Schartl, 2014). Abnormalities in the EDNRB2 gene cause defects in pigment biosynthesis, which results in the production of pink skin and white fur (Kelsh et al., 2009; Fontanesi et al., 2014; Li L. et al., 2015). EDNRB2 is also involved in melanocyte differentiation and migration in vitro and in ovo (Pla et al., 2005). To date, EDNRB2 has been found to be associated with phenotypes in quail (Miwa et al., 2007), chickens (Kinoshita et al., 2014), frog (Kawasaki-Nishihara et al., 2011) and domestic duck (Li L. et al., 2015). In this study, EDNRB2 was also found to be upregulated in chicken melanocytes at the transcriptional level. Although the expression levels of EDNRB2 did not continue to increase with increases in the concentrations of tyrosine, its expression remained significantly higher than that in the group not treated with tyrosine. These studies indicate that EDNRB2 plays an important role in tyrosine-induced melanin deposition.
Based on the above results and those of a previous study, a network for the EDNRB2-mediated regulation of chicken melanin synthesis was constructed (Figure 11), and this network plays an important role in the regulation of melanocyte development (Hou et al., 2004; Saldana-Caboverde and Kos, 2010; Kawasaki-Nishihara et al., 2011). In this proposed network, ET3, an isomer of the endothelin receptor, has similar affinities for EDNRB and EDNRB2 and binds to these receptors to activate the downstream signal transduction molecule PKC, and the activation of PKC in turn triggers a physiological response in the cell (Pla et al., 2005; Kinoshita et al., 2014). This response then activates the expression of many specific genes, such as dopachrome tautomerase (DCT) (Guibert et al., 2004; Saternus et al., 2015; Salinas-Santander et al., 2018), tyrosinase (TYR) (Guibert et al., 2004; Xu et al., 2013), and tyrosinase-related protein 1 (TYRP1) (Alonso et al., 2008; Saternus et al., 2015), in melanocytes. As a dopamine isomerase downstream of tyrosine, DCT can catalyze the transformation of dopamine into 5,6-dihydroxyindole carboxylic acid, which can accelerate the formation of melanin and in turn affect melanin deposition in black-bone chicken breast muscle tissue.
Figure 11. Network of EDNRB2-mediated regulation of chicken melanin synthesis. In melanocytes, the ET3 gene binds to EDNRB2 gene to activates the downstream signal transduction molecule PKC, and the activation of PKC activates the expression of many specific genes such as TYR, TYRP1, and DCT. The DCT gene can catalyze the transformation of dopamine into 5,6-dihydroxyindole carboxylic acid to accelerate the formation of melanin. ET3, endothelin 3; EDNRB, endothelin receptor B; EDNRB2, endothelin receptor B subtype 2; PKC, protein kinase C; MITF, melanogenesis associated transcription factor; TYR, tyrosinase; TYRP1, tyrosinase-related protein 1; DCT, dopachrome tautomerase.
Conclusion
In summary, this study provides the first demonstration that tyrosine can promote melanin deposition in the breast muscle of Xichuan black-bone chickens at the cellular level and in vivo. For Xichuan black-bone chickens, the optimal tyrosine dietary supplementation level was found to be 0.8%. We also analyzed the characteristics of the transcriptome profiles of breast muscle tissue from the control group and test group IV (0.8% dietary tyrosine level and the highest melanin content). The EDNRB2-mediated regulatory network might play a key role in tyrosine-induced melanin deposition in the breast muscle of Xichuan black-bone chickens. These results provide insights into the molecular regulatory mechanisms for melanin deposition in black-bone chickens and are important for the nutrition regulation and the genetic improvement of the meat quality of Xichuan black-bone chickens.
Ethics Statement
The whole sample collection and treatment were conducted in strict accordance with the protocol approved by the Institutional Animal Care and Use Committee (IACUC) of Henan Agricultural University, China (11-0099).
Author Contributions
XK and GL conceived and designed the experiments. YC, JW, and YC contributed to animal husbandry. YZ, YL, and YF performed the experiments. YT and XW analyzed the data. CZ, RJ, and ZL contributed reagents and materials. GS and WL contributed analysis tools. DL wrote the manuscript. All authors reviewed and approved the final manuscript.
Funding
This work was supported by the Earmarked Fund for the National Natural Science Foundation of China (31572356), Modern Agro-Industry Technology Research System (CARS-40-K04), Agricultural Research Outstanding Talents and Innovation Team of Ministry of Agriculture, Program for Innovation Research Team of Ministry of Education (IRT1236), NSFC-Henan joint grant (U1804107), and Young Backbone Teacher Training Program of Henan Higher Educational Institution (2017GGJS036).
Conflict of Interest Statement
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.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2019.00457/full#supplementary-material
Footnote
References
Alonso, S., Izagirre, N., Smith-Zubiaga, I., Gardeazabal, J., Díaz-Ramón, J. L., Díaz-Pérez, J. L., et al. (2008). Complex signatures of selection for the melanogenic loci TYR, TYRP1 and DCT in humans. BMC Evol. Biol. 8:74. doi: 10.1186/1471-2148-8-74
Bauer, A., Kehl, A., Jagannathan, V., and Leeb, T. (2018). A novel MLPH variant in dogs with coat colour dilution. Anim. Genet. 49, 94–97. doi: 10.1111/age.12632
Bed’hom, B., Vaez, M., Coville, J. L., Gourichon, D., Chastel, O., Follett, S., et al. (2012). The lavender plumage colour in Japanese quail is associated with a complex mutation in the region of MLPH that is related to differences in growth, feed consumption and body temperature. BMC Genomics 13:442. doi: 10.1186/1471-2164-13-442
Braasch, I., and Schartl, M. (2014). Evolution of endothelin receptors in vertebrates. Gen. Comp. Endocrinol. 209, 21–34. doi: 10.1016/j.ygcen.2014.06.028
Diffey, B. L., Healy, E., Thody, A. J., and Rees, J. L. (1995). Melanin, melanocytes, and melanoma. Lancet 346:1713. doi: 10.1016/s0140-6736(95)92882-0
Fontanesi, L., Vargiolu, M., Scotti, E., Latorre, R., Faussone Pellegrini, M. S., Mazzoni, M., et al. (2014). The KIT gene is associated with the English spotting coat color locus and congenital megacolon in checkered giant rabbits (Oryctolagus cuniculus). PLoS One 9:e93750. doi: 10.1371/journal.pone.0093750
Franceschini, A., Szklarczyk, D., Frankild, S., Kuhn, M., Simonovic, M., Roth, A., et al. (2013). STRING v9.1: protein-protein interaction networks, with increased coverage and integration. Nucleic Acids Res. 41, D808–D815. doi: 10.1093/nar/gks1094
Guibert, S., Girardot, M., Leveziel, H., Julien, R., and Oulmouden, A. (2004). Pheomelanin coat colour dilution in french cattle breeds is not correlated with the TYR, TYRP1 and DCT transcription levels. Pigment Cell Res. 17, 337–345. doi: 10.1111/j.1600-0749.2004.00152.x
Harris, M. L., Hall, R., and Erickson, C. A. (2008). Directing pathfinding along the dorsolateral path - the role of EDNRB2 and EphB2 in overcoming inhibition. Development 135, 4113–4122. doi: 10.1242/dev.023119
Haydock, A., Terrao, M., Sekar, A., Ramasamy, G., Baugh, L., and Myler, P. J. (2015). RNA-Seq approaches for determining mRNA abundance in Leishmania. Methods Mol. Biol. 1201, 207–219. doi: 10.1007/978-1-4939-1438-8_12
Hirano, S. (1990). Observations on pigment granules in the bones of silky fowls. Arch. Histol. Cytol. 53, 89–93. doi: 10.1679/aohc.53.89
Hou, L., Pavan, W. J., Shin, M. K., and Arnheiter, H. (2004). Cell-autonomous and cell non-autonomous signaling through endothelin receptor B during melanocyte development. Development 131, 3239–3247. doi: 10.1242/dev.01193
Hume, A. N., Tarafder, A. K., Ramalho, J. S., Sviderskaya, E. V., and Seabra, M. C. (2006). A coiled-coil domain of melanophilin is essential for Myosin Va recruitment and melanosome transport in melanocytes. Mol. Biol. Cell 17, 4720–4735. doi: 10.1091/mbc.e06-05-0457
Hume, A. N., Ushakov, D. S., Tarafder, A. K., Ferenczi, M. A., and Seabra, M. C. (2007). Rab27a and MyoVa are the primary Mlph interactors regulating melanosome transport in melanocytes. J. Cell Sci. 120, 3111–3122. doi: 10.1242/jcs.010207
Kawasaki-Nishihara, A., Nishihara, D., Nakamura, H., and Yamamoto, H. (2011). ET3/Ednrb2 signaling is critically involved in regulating melanophore migration in Xenopus. Dev. Dyn. 240, 1454–1466. doi: 10.1002/dvdy.22649
Kelsh, R. N., Harris, M. L., Colanesi, S., and Erickson, C. A. (2009). Stripes and belly-spots – a review of pigment cell morphogenesis in vertebrates. Semin. Cell Dev. Biol. 20, 90–104. doi: 10.1016/j.semcdb.2008.10.001
Kim, K. I., and van de Wiel, M. A. (2008). Effects of dependence in high-dimensional multiple testing problems. BMC Bioinformatics 9:114. doi: 10.1186/1471-2105-9-114
Kinoshita, K., Akiyama, T., Mizutani, M., Shinomiya, A., Ishikawa, A., Younis, H. H., et al. (2014). Endothelin receptor B2 (EDNRB2) is responsible for the tyrosinase-independent recessive white (mo w) and mottled (mo) plumage phenotypes in the chicken. PLoS One 9:e86361. doi: 10.1371/journal.pone.0086361
Lecoin, L., Sakurai, T., Ngo, M. T., Abe, Y., Yanagisawa, M., and Le Douarin, N. M. (1998). Cloning and characterization of a novel endothelin receptor subtype in the avian class. Proc. Natl. Acad. Sci. U.S.A. 95, 3024–3029. doi: 10.1073/pnas.95.6.3024
Lehner, S., Gähle, M., Dierks, C., Stelter, R., Gerber, J., Brehm, R., et al. (2013). Two-exon skipping within MLPH Is associated with coat color dilution in rabbits. PLoS One 8:e84525. doi: 10.1371/journal.pone.0084525
Leonhardt, R. M., Vigneron, N., Jia, S. H., Graham, M., and Cresswell, P. (2013). Critical residues in the PMEL/Pmel17 N-terminus direct the hierarchical assembly of melanosomal fibrils. Mol. Biol. Cell 24, 964–981. doi: 10.1091/mbc.E12-10-0742
Li, G. H., Xu, H. Y., Xu, L. J., Qu, M. R., You, J. M., Yi, Z. H., et al. (2011). Effects of dietary selenium supplementation on growth performance and melanin content in tissues of taihe silky fowls. Sci. Agric. Sin. 44, 2777–2786. doi: 10.3864/j.issn.0578-1752.2011.13.016
Li, H., Wang, T., Xu, C., Wang, D., Ren, J., Li, Y., et al. (2015). Transcriptome profile of liver at different physiological stages reveals potential mode for lipid metabolism in laying hens. BMC Genomics 16:763. doi: 10.1186/s12864-015-1943-0
Li, L., Li, D., Liu, L., Li, S., Feng, Y., Peng, X., et al. (2015). Endothelin receptor B2 (EDNRB2) gene is associated with spot plumage pattern in domestic ducks (Anas platyrhynchos). PLoS One 10:e0125883. doi: 10.1371/journal.pone.0125883
Mackintosh, J. A. (2001). The antimicrobial properties of melanocytes, melanosomes and melanin and the evolution of black skin. J. Theor. Biol. 211, 101–113. doi: 10.1006/jtbi.2001.2331
Mao, X., Cai, T., Olyarchuk, J. G., and Wei, L. (2005). Automated genome annotation and pathway identification using the KEGG Orthology (KO) as a controlled vocabulary. Bioinformatics 21, 3787–3793. doi: 10.1093/bioinformatics/bti430
Miwa, M., Inoue-Murayama, M., Aoki, H., Kunisada, T., Hiragaki, T., Mizutani, M., et al. (2007). Endothelin receptor B2 (EDNRB2) is associated with the panda plumage colour mutation in Japanese quail. Anim. Genet. 38, 103–108. doi: 10.1111/j.1365-2052.2007.01568.x
Morris, J. G., Yu, S., and Rogers, Q. R. (2002). Red hair in black cats is reversed by addition of tyrosine to the diet. J. Nutr. 132, 1646S–1648S. doi: 10.1093/jn/132.6.1646S
Pla, P., Alberti, C., Solov’eva, O., Pasdar, M., Kunisada, T., and Larue, L. (2005). Ednrb2 orients cell migration towards the dorsolateral neural crest pathway and promotes melanocyte differentiation. Pigment Cell Res. 18, 181–187. doi: 10.1111/j.1600-0749.2005.00230.x
Poelchau, M. F., Huang, X., Goff, A., Reynolds, J., and Armbruster, P. (2014). An experimental and bioinformatics protocol for RNA-Seq analyses of photoperiodic diapause in the Asian tiger mosquito, Aedes albopictus. J. Vis. Exp. 93:e51961. doi: 10.3791/51961
Ralli, E. P., and Graef, I. (1943). Stimulating effect of adrenalectomy on hair growth and melanin deposition in black rats fed diets adequate and deficient in the filtrate factors of vitamin B. With a note on effect of diet on post-operative survival. Endocrinology 32, 1–12. doi: 10.1210/endo-32-1-1
Robinson, M. D., and Oshlack, A. (2010). A scaling normalization method for differential expression analysis of RNA-seq data. Genome Biol. 11:R25. doi: 10.1186/gb-2010-11-3-r25
Ruemmele, F. M., and Garnier-Lengliné, H. (2012). Why are genetics important for nutrition? Lessons from epigenetic research. Ann. Nutr. Metab. 60, 38–43. doi: 10.1159/000337363
Saldana-Caboverde, A., and Kos, L. (2010). Roles of endothelin signaling in melanocyte development and melanoma. Pigment Cell Melanoma Res. 23, 160–170. doi: 10.1111/j.1755-148X.2010.00678.x
Salinas-Santander, M., Trevino, V., De la Rosa-Moreno, E., Verduzco-Garza, B., Sánchez-Domínguez, C. N., Cantú-Salinas, C., et al. (2018). CAPN3, DCT, MLANA and TYRP1 are overexpressed in skin of vitiligo vulgaris Mexican patients. Exp. Ther. Med. 15, 2804–2811. doi: 10.3892/etm.2018.5764
Saternus, R., Pilz, S., Gräber, S., Kleber, M., März, W., Vogt, T., et al. (2015). A closer look at evolution: variants (SNPs) of genes involved in skin pigmentation, including EXOC2, TYR, TYRP1, and DCT, are associated with 25(OH)D serum concentration. Endocrinology 156, 39–47. doi: 10.1210/en.2014-1238
Sato, S., Otake, T., Suzuki, C., Saburi, J., and Kobayashi, E. (2007). Mapping of the recessive white locus and analysis of the tyrosinase gene in chickens. Poult. Sci. 86, 2126–2133. doi: 10.1093/ps/86.10.2126
Scheiner, S., Kar, T., and Pattanayak, J. (2002). Comparison of various types of hydrogen bonds involving aromatic amino acids. J. Am. Chem. Soc. 124, 13257–13264. doi: 10.1021/ja027200q
Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., et al. (2003). Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 13, 2498–2504. doi: 10.1101/gr.1239303
Słominski, A., Moellmann, G., Kuklinska, E., Bomirski, A., and Pawelek, J. (1988). Positive regulation of melanin pigmentation by two key substrates of the melanogenic pathway, L-tyrosine and L-dopa. J. Cell Sci. 89, 287–296.
Tian, Y., Xie, M., Wang, W., Wu, H., Fu, Z., Lin, L., et al. (2007). Determination of carnosine in black-bone silky fowl (Gallus gallus domesticus Brisson) and common chicken by HPLC. Eur. Food Res. Technol. 226, 311–314. doi: 10.1007/s00217-006-0528-1
Trapnell, C., Roberts, A., Goff, L., Pertea, G., Kim, D., Kelley, D. R., et al. (2012). Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat. Protoc. 7, 562–578. doi: 10.1038/nprot.2012.016
Tsatmali, M., Ancans, J., and Thody, A. J. (2002). Melanocyte function and its control by melanocortin peptides. J. Histochem. Cytochem. 50, 125–133. doi: 10.1177/002215540205000201
Tu, Y. G., Sun, Y. Z., Tian, Y. G., Xie, M. Y., and Chen, J. (2009a). Physicochemical characterisation and antioxidant activity of melanin from the muscles of Taihe Black-bone silky fowl (Gallus gallus domesticus Brisson). Food Chem. 114, 1345–1350. doi: 10.1016/j.foodchem.2008.11.015
Tu, Y. G., Xie, M. Y., Sun, Y. Z., and Tian, Y. G. (2009b). Structural characterization of melanin from black-bone silky fowl (Gallus gallus domesticus Brisson). Pigment Cell Melanoma Res. 22, 134–136. doi: 10.1111/j.1755-148X.2008.00529.x
Vandesompele, J., De Preter, K., Pattyn, F., Poppe, B., Van Roy, N., De Paepe, A., et al. (2002). Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multipleinternal control genes. Genome Biol. 3:research0034.1.
Walker, W. A., and Blackburn, G. (2004). Symposium introduction: nutrition and gene regulation. J. Nutr. 134, 2434S–2436S. doi: 10.1093/jn/134.9.2434s
Xie, M. Y., Tian, Y. G., and Yong-Gang, T. U. (2009). Review of the bioactive components of black-bone silky fowl and their functions. Modern Food Sci. Technol. 25, 461–465.
Xin, L. I., Tian, Y. G., Xie, M. Y., Sun, Y. Z., Liu, J. H., Tu, Y. G., et al. (2009). Extraction, isolation and antioxidant activity of melanin from black-bone silky fowl meat. Food Sci. 30, 117–120.
Xu, X. L., Zhang, S., and Chen, B. X. (1999). The effect of melanin from white silky fowl on antiageing in Drosophila melanogaster. J. Nanjing Agric. Univers. 22, 105–108.
Xu, Y., Zhang, X. H., and Pang, Y. Z. (2013). Association of tyrosinase (TYR) and tyrosinase-related protein 1 (TYRP1) with melanic plumage color in Korean quails (Coturnix coturnix). Asian-Australas. J. Anim. Sci. 26, 1518–1522. doi: 10.5713/ajas.2013.13162
Young, M. D., Wakefield, M. J., Smyth, G. K., and Oshlack, A. (2010). Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 11:R14. doi: 10.1186/gb-2010-11-2-r14
Yu, S., Rogers, Q. R., and Morris, J. G. (2001). Effect of low levels of dietary tyrosine on the hair colour of cats. J. Small Anim. Pract. 42, 176–180. doi: 10.1111/j.1748-5827.2001.tb01798.x
Yu, S., Wang, G., Liao, J., and Tang, M. (2018a). Transcriptome profile analysis identifies candidate genes for the melanin pigmentation of breast muscle in Muchuan black-boned chicken. Poult. Sci. 97, 3446–3455. doi: 10.3382/ps/pey238
Yu, S., Wang, G., Liao, J., Tang, M., and Sun, W. (2018b). Transcriptome profile analysis of mechanisms of black and white plumage determination in black-bone chicken. Cell Physiol. Biochem. 46, 2373–2384. doi: 10.1159/000489644
Keywords: breast muscle, chicken, melanin, melanocytes, tyrosine, transcriptome
Citation: Li D, Wang X, Fu Y, Zhang C, Cao Y, Wang J, Zhang Y, Li Y, Chen Y, Li Z, Li W, Jiang R, Sun G, Tian Y, Li G and Kang X (2019) Transcriptome Analysis of the Breast Muscle of Xichuan Black-Bone Chickens Under Tyrosine Supplementation Revealed the Mechanism of Tyrosine-Induced Melanin Deposition. Front. Genet. 10:457. doi: 10.3389/fgene.2019.00457
Received: 08 February 2019; Accepted: 30 April 2019;
Published: 15 May 2019.
Edited by:
Enrique Medina-Acosta, Universidade Estadual do Norte Fluminense Darcy Ribeiro, BrazilReviewed by:
Jong-Eun Park, National Institute of Animal Science, South KoreaKatia Cappelli, University of Perugia, Italy
Copyright © 2019 Li, Wang, Fu, Zhang, Cao, Wang, Zhang, Li, Chen, Li, Li, Jiang, Sun, Tian, Li and Kang. 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: Guoxi Li, bGlndW94aTA5MTRAMTI2LmNvbQ== Xiangtao Kang, eHRrYW5nMjAwMUAyNjMubmV0