|Home | About | Journals | Submit | Contact Us | Français|
Phenol oxidases (POs) catalyze the oxidation of dopa and dopamine to melanin, which is crucial for cuticle formation and innate immune maintenance in insects. Although, Laccase 2, a member of the PO family, has been reported to be a requirement for melanin-mediated cuticle tanning in the development stages of some insects, whether it participates in cuticle construction and other physiological processes during the metamorphosis of mosquito pupae is unclear.
The association between the phenotype and the expression profile of Anopheles sinensis Laccase 2 (AsLac2) was assessed from pupation to adult eclosion. Individuals showing an expression deficiency of AsLac2 that was produced by RNAi and their phenotypic defects and physiological characterizations were compared in detail with the controls.
During the dominant expression period, knockdown of AsLac2 in pupae caused the cuticle to be unpigmented, and produced thin and very soft cuticles, which further impeded the eclosion rate of adults as well as their fitness. Moreover, melanization immune responses in the pupae were sharply decreased, leading to poor resistance to microorganism infection. Both the high conservation among Laccase 2 homologs and a very similar genomic synteny of the neighborhood in Anopheles genus implies a conservative function in the pupal stage.
To our knowledge, this is the first study to report the serious phenotypic defects in mosquito pupae caused by the dysfunction of Laccase 2. Our findings strongly suggest that Laccase 2 is crucial for Anopheles cuticle construction and melanization immune responses to pathogen infections during pupal metamorphosis. This irreplaceability provides valuable information on the application of Lacccase 2 and/or other key genes in the melanin metabolism pathway for developing mosquito control strategies.
The online version of this article (doi:10.1186/s13071-017-2118-4) contains supplementary material, which is available to authorized users.
Melanin is not only the substrate used for cuticle tanning in insects, but it is also involved in innate immune responses against exogenous pathogen infections through melanotic encapsulation [1–8]. Thus, melanin metabolism determines whether these two essential physiological processes can be normally activated and affect the typical development of insects. For the Mosquitoes, the adults pose a threat to human health because they transmit malignant diseases by biting humans [9–11]. During the pre-developmental stage of the adult, the development status of the mosquito pupae directly affects the eclosion rate, and further affects the population growth rate of adults.
During pupal development, the cuticle gradually darkens and sclerotizes, which provides enough support for pupae to break out of the puparium, protects them from mechanical injuries, and facilitates the emergence of adults; in addition, some intermediates of melanin metabolism are known to be involved in encapsulation, which is helpful for the survival of the pupae in natural ecosystems [7, 11–16]. Therefore, melanin metabolism during the pupal stage is crucial for mosquito development, survival, and reproduction.
The metabolism of melanin with the hydroxylation of tyrosine to dihydroxy phenylalanine (dopa) by the rate-limiting enzyme tyrosine hydroxylase (TH), followed by the decarboxylation of dopa to dopamine by the dopa decarboxylase (DDC) [5, 17–19]. These melanin precursors must be further oxidized by phenol oxidases (POs) to quinones and quinone methides [19–25]. Subsequently, the quinones and quinone methides conjugate with cuticular proteins to construct and tan the cuticle, and are also involved in melanotic encapsulation [3, 5, 14, 15, 17, 18, 22–25]. Laccase, a member of the phenol oxidase family, conservatively contains three cupredoxin-like domains and four copper ions that reside in a T1 copper site and a T2/T3 tricopper center [5, 26]. Laccases oxidize a broad range of substrates, including polyphenol, methoxy-substituted phenol, aminophenol, and phenylenediamine [5, 27–29]. There are two major types of Laccase genes, Laccase 1 and Laccase 2, that have been identified in many insects [30–44]. Laccase 1 has been reported to likely be involved in cuticular sclerotization [33, 35, 37], while Laccase 2 is mainly expressed in the cuticle or egg shell and is directly involved in melanin-mediated cuticle tanning [38–44]. In the red flour beetle, stinkbug, and honey bee, the pupal cuticle or newly-molted adult becomes white and more flexible due to the depletion of Laccase 2 [38, 40, 42]. Additionally, the dysfunction of Laccase 2 in mosquito eggs results in pale and fragile eggshells, which further causes the eventual collapse of the eggs . This evidence suggests that Laccase 2 is an important regulator in melanin synthesis and deposition. Its dysfunction results in melanin precursors that fail to be oxidized, and impairs the pigmentation and sclerotization of the cuticle or exoskeleton. Although the expression of Laccase 2 is enriched in the epidermis, we speculated that it is also expressed in immune tissues, such as hemolymph or fat body, and may be involved in melanin synthesis, which participates in melanotic encapsulation immune responses. However, knowledge of Laccase 2 functions in cuticle tanning and pathogen resistance during mosquito pupal development is still limited.
Knockdown of the Laccase 2 gene (AsLac2) in the pupal stage was performed by RNAi to better understand its function in the pupal development of Anopheles sinensis. Our results revealed that the pupae with an unpigmented cuticle developed abnormally, had a severe mortality rate, and were susceptible to exogenous microorganism infection. These findings support that AsLac2 is not only required for cuticle tanning, but also for melanotic immune responses of the malaria vector mosquitoes, Anopheles sinensis. Moreover, the high genomic synteny of the neighborhood of Laccase 2 in species of the genus Anopheles may be a strong indicator of the gene functional conservation in the specific developmental stages. The present study deepens our understanding of the function of the AsLac2 gene in mosquitoes and provides a new reference for mosquito control.
The Anopheles sinensis LS-WX strain was reared at 27 °C with 80% humidity under a 12 h/12 h (light/dark) photoperiod. The larvae in the different developmental stages were fed fry food in clean water, and the adults were provided 10% glucose solution.
BLAST analysis was performed to search homologous Laccase from the Anopheles sinensis genome and transcriptome databases using three insect Laccase 2 proteins (Protein ID: AGAP006176-PB, NP_001034487, and BAG70891) as queries. Fragments from the transcriptome data  were assembled using the SeqMan program (https://www.dnastar.com/). The signal peptide was predicted by SignalP 4.1 (http://www.cbs.dtu.dk/services/SignalP) and conserved domains were analyzed by SMART (http://smart.embl-heidelberg.de/). Total RNA was extracted from the pupae (after pupation 32 h) using the TRizol reagent (Invitrogen, Shanghai, China) according to the manufacturer’s instructions. Total RNA (1 μg) was reverse-transcribed with random primer using the First-Strand cDNA Synthesis Kit (Takara, Dalian, China). Five pairs of primers were designed based on the putative AsLac2 sequence to obtain the full open reading frame (ORF). The sequence of 3′UTR was obtained by the rapid amplification of cDNA ends (RACE) technique using the GeneRacer kit (Invitrogen, Shanghai, China). PCR products were isolated and subcloned into PMD-19 vector (Takara, Dalian, China) for sequencing. The primers are listed in Additional file 1: Table S1.
Homologous Laccase were searched in the NCBI database (http://www.ncbi.nlm.nih.gov/) and Vectorbase (http://www.vectorbase.org/) using the BLASTP program. Amino acid sequences of the divergence part in the C-terminus were aligned with MUSCLE (http://www.ebi.ac.uk/Tools/msa/muscle/). The best-fit evolutionary model (WAG+G) and the genetic distance were estimated by MEGA 5.0 (http://www.megasoftware.net/) . Maximum likelihood phylogenetic analysis was conducted by MEGA 5.0 and bootstrap values were obtained based on 1,000 bootstrap replications. Additionally, the An. gambiae Laccase 2 gene and its adjacent genes were used as the templates to search for Laccase 2 homologs in other insect genomes, and their genome locations and distributions were compared in detail. We selected 7.5 kb upstream DNA sequences of the Laccase 2 gene to predict cis-acting regulatory elements would respond to hormonal signals during pupal development using JASPAR (http://jaspar.genereg.net/) with 90% confidence settings.
Three individuals at each sampling point were collected at different development stages (pupae were collected at 0, 8, 16, 24, 32 and 40 h after pupation and the adults were collected at 0, 3, 6, 9, 12, 24, 48 and 72 h after the eclosion) for phenotype observations (Olympus, Shanghai, China) and gene expression analysis of AsLac2. The cuticle, fat body, and hemolymph of the pupae were dissected at 32 h after pupation for tissue expression pattern analysis. Gene expression was determined by qRT-PCR as previously described . Each sample was used to perform three biological repeats. The Ribosomal protein L49 (RPL49) gene was used as the internal control (see primers listed in Additional file 1: Table S1).
Double-stranded RNA (dsRNA) of AsLac2 (dsLac2 and dsLac2-2, designed to avoid the off-target effect) were synthesized with a T7 RiboMAX™ Express RNAi System (Promega, Chongqing, China) according to the manufacturer’s instructions. For microinjection, the dsRNA was dissolved in RNase-free water and the concentration was estimated. Based on AsLac2 temporal expression patterns, each of the pupae was injected with 800 ng of dsRNA (volume: 180–200 nl) into the thorax through the back of the dorsal plate within 2 h after the pupation according to the temporal expression pattern of AsLac2. A red fluorescent protein gene (dsRed) was used as the control. The tanning degree of the pupae was checked at 38 h after dsLac2 or dsRed injection, while the checking time for adults was 3 h after the eclosion. Three pupae or adults in the dsLac2- and dsRed-injected groups at each sampling point were also collected for qRT-PCR analyses . Each sample was used to perform three biological repeats. The Ribosomal protein L49 (RPL49) gene was used as the internal control. The primers used for dsRNA synthesis are listed in Additional file 1: Table S1.
The dorsal plates of the adults from dsLac2- and dsRed-injected groups were dissected under the microscope (Olympus, Japan). All of the tissue treatments and slice preparations were carried out as previously described . The thickness of each cuticle layer was measured with Image J software (https://imagej.net/Welcome) as previously described . Four equal spaced measurements for each dorsal plate section, and the average value was calculated as the mean thickness in one biological repeat. Four biological replicates were performed in the dsLac2- and the dsRed-injected groups, respectively.
The pathogens, Serratia marcescens (Sm) and Bacillus bombyseptie (Bb), were incubated to the logarithmic phase (OD600=0.6–0.8) in LB medium at 37 °C. Individuals in the dsLac2- (failed to pigment) and the dsRed-injected groups were sampled according to the melaninization degree, and then injected with 0.12 μl of bacteria solution at 26 h after pupation, respectively (the colour of the pupa cuticle in the dsRed group started to darken, while the cuticle of Lac2-silenced pupae failed to tan). The pupal survival rate was recorded every 2 h after being infected. Six pupae at the same developmental stage of the same size were homogenized at 4 °C with 400 μl of PBS (pH=7.0) and centrifuged (500× g for 5 min at 4 °C). Total protein concentrations (8 mg/ml) were determined using the Bradford method (BBI, Hong Kong, China). Melanization reactions were incubated at 30 °C for 3 h, followed by adding 1 mM of phenylthiourea (PTU) to terminate the reactions and A490 nm values were measured to estimate the amount of melanin. Each sample was used to perform three biological repeats.
The Anopheles sinensis Laccase 2 (AsLac2) is 4,391 bp (without poly A) with a 2,265 bp ORF encoding 754 amino acids, and a 2,126 bp 3′UTR (GenBank Number: KY132102); a 26-residue signal peptide was predicted, suggesting it is a secreted protein. Conserved domain prediction revealed that AsLac2 contains three domains: Cu-oxidase3, Cu-oxidase, and Cu-oxidase 2 (Fig. 1a).
Phylogenetic analysis showed that representative insect Laccases were grouped into two distinct clades (Fig. 1b), Laccase 1 and Laccase 2. The newly cloned Laccase gene was identified as a Laccase 2 homologue. Three Cu-oxidase domains for typical Laccase 2 are highly conserved in different insects and exhibit 81.6 to 100% identities (Additional file 2: Table S2). In particularly, within the genus Anopheles, the identities range from 92.3% to 100% for three Cu-oxidase domains (Additional file 2: Table S2). These features suggest that their functions are likely conserved in different insect species. Moreover, chromosomal or scaffold synteny analyses revealed the neighborhood genes around Laccase 2 are more conservative across the genus Anopheles than in other insect species (Fig. 1c), which may imply the conservation of this genome fragment during the evolution and divergence of mosquitoes. Further analyses found that upstream regulatory elements of Lac2, like ‘BR-C’, ‘E74’, ‘FTZ’, and ‘ECR::USP’, occurs in all analyzed Lac2 (Additional file 3: Table S3), suggesting that Lac2 is regulated by ecdysone during development [47–49].
Both male and female Anopheles sinensis reared in our laboratory required about 40–42 h to complete the eclosion. According to our observations of pupae phenotype in different developmental stages, the cuticle tanning pattern is the same in the male and female pupae . The cuticle colour of the female pupae is a little transparent with a slight yellowish tint within 24 h after pupation (Fig. 2a) and subsequently, strong melanization and hardening were observed until the late pupal stage, especially at 32 h after pupation (Fig. 2a). Therefore, we concluded that the expression levels of key genes involved in melanin synthesis and cuticle tanning should be significantly increased with the tanning degrees. The expression of AsLac2 was almost undetectable at 0 h after pupation, while its expression started to be upregulated at 8 h, was remarkably increased at 16 h, and was maintained to 40 h after pupation (Fig. 2b). This expression pattern positively correlated well with the cuticle melanization and sclerotization process, but negatively correlated with the derived hemolymph ecdysone titer during the mosquito pupal stage  (Fig. 2b). Interestingly, a retarded expression pattern of AsLac2 was observed when compared with the upstream, rate-limited gene AsTH , indicating that the catalytic function of AsLac2 occurred after the accumulation of catecholamine. Furthermore, we found that the expression intensity of AsLac2 in the pupae (from the early-middle to the last stage) was significantly stronger than that in the adults (Fig. 2a), implying its importance in normal pupal development (Fig. 2). Therefore, we concluded that Aslac2 is closely linked to pupae cuticle tanning.
Normally, the pupal cuticle darkens and gradually hardens within 24 h of development, and becomes much darker at 38 h after pupation. In this study, no pigmentation blocking was observed in the dsRed-injected individuals as they exhibited a normal development and tanning process (Fig. 3a). In contrast, dsLac2-injected pupae exhibited very little or even no melanism at 38 h after pupation, and were accompanied by very soft cuticles (Fig. 3a). At this time point, the expression level of AsLac2 was much lower than that in the dsRed-injected group (Fig. 3a). To avoid the off-target effect, we designed another pair of primers (dsLas2-2) without overlapping with dsLac2. There were 80 and 66% of individuals whose cuticle was not tanned and flexible that were obtained in the dsLac2 and dsLac2-2 groups, respectively (Fig. 3a, Additional file 4: Table S4). In the dsRed-injected group, the cephalothorax was almost completely tanned at 3 h after adult emergence (Fig. 3b); however, the dsLac2-injected pupae were extremely struggling to emerge and the adults displayed a very light body colour and soft cephalothorax (Fig. 3b). Moreover, we found that the thickness of the dorsal plate section of the cephalothorax in the dsLac2-injected pupae (very fragile and easy to break) (Fig. 4b) was only about 54% (Fig. 4c) of that in the control group (Fig. 4a). The expression level of dsLac2 was also significantly lower than in the control group (Fig. 3b). Our present data indicates that AsLac2 dysfunction can lead to cuticle defects in the melanizaiton and sclerotization processes, suggesting that it is crucial for cuticle tanning (Fig. 3, Fig. 4a, ,bb and andc).c). In the control group, the pupae started to emerge at 37 h and completed the pupation at 41 h with the eclosion rate reaching 79% (Fig. 4d). However, only a 50% eclosion rate was obtained from 45 h to 49 h after pupation in the dsLac2-injected group (Fig. 4d). Obviously, there was an 8-h eclosion retardation in the dsLac2-injected groups (Fig. 4d). Some of the remaining pupae barely emerged, were physically weak and had wing vibration defects, which distinctly impaired the adults’ ability to survive. The obvious time lag, low emergence rate, and physical weakness definitely affected the mosquito population.
In the dsLac2-injected group, pupa death started at 2 h after Bb challenge, then the mortality increased with the time going from 2 to 10 h and all individuals died at 10 h after the injection (Fig. 5a). In contract, large-scale mortality in the dsRed group occurred at 20–24 h after the bacteria challenge, which was significantly delayed compared to those in the dsLac2-injected group (P<0.01) (Fig. 5a). Although, the mortality rate in the dsRed- and dsLac2-injected groups caused by Sm was not significantly more severe that those infected by Bb, the survival and life span of the dsLac2-injected pupae were still obviously decreased by Sm (P<0.01) (Fig. 5a). The in vivo assay results revealed that the melanization degree of dsLac2-injected pupal humor was significantly lower than those in the controls after Bb or Sm infection, suggesting that the melanization immune response was severely attenuated in the dsLac2-injected pupae (Fig. 5b). In conclusion, it is worth noting that dsLac2-injected pupae died rapidly and severely after Bb or Sm infections (Fig. 5). These results verified that the AsLac2 gene is involved in Anopheles sinensis’ innate immunity. AsLac2 dysfunction greatly suppressed the melanization immune responses to bacterial infections, which further caused quick and severe death.
Laccase 2 is mainly synthesized in epidermis cells and secreted to the procuticle [5, 34, 38, 40–44]. Subsequently, it catalyzes the oxidation of catecholamine to quinones and/or quinone methides, which are involved in cuticle pigmentation, as well as stabilizes and reinforces the cuticle structure via catalyzing the cross-linking among various cuticular components, such as proteins-proteins, proteins-melanin, or proteins-other quinones and/or quinone methides [3–5, 44, 50]. This implies that its normal function controls the supply of important materials for cuticle construction and maintenance. Therefore, the expression pattern of Laccase 2 in a specific developmental stage should reflect the characteristics and changes of the cuticle. In the present study, the expression of AsLac2 was significantly upregulated in the middle and the late pupal stage (Fig. 2). The good agreement of this expression pattern and rapid cuticle tanning strongly supports this view. Serious defects in the pupal cuticle caused by AsLac2 silencing in the high expression period undoubtedly imply that it plays a decisive role in the maintenance of cuticle integrity during mosquito pupal metamorphosis.
Moreover, in most insects, Laccase 2 is irreplaceable and dominantly expressed in pupae, while other melanin metabolism genes, such as TH or Yellow [38, 40, 51–53], are not always dominantly expressed in the pupal stage. Therefore, we thought that the pupae could protect itself from external adverse factors mainly through passive defenses during pupal and/or eclosion periods. The high expression of Laccase 2 in this stage can promote quinone and/or quinonemethide accumulation to enhance, harden, and steady the pupal cuticle (Figs. 2, ,33 and and4).4). It is also beneficial for adult cuticle formation to adapt to the environment. On the other hand, ecdysone titer peaks during the early and the middle pupa developmental stage further trigger apolysis and promote the formation of the new pupal cuticle [47–49]. The upregulation of AsLac2 in this period verifies the general rule that the formation of the new cuticle and the expression patterns of Laccase 2 are modulated by the ecdysone titer and their functional conservations in the insect evolutionary processes.
Two alternative isoforms of Laccase 2, Lac2A and Lac2B, were identified in some insects [34, 38]. They share the same N-terminal structure and three Cu-oxidase domains, but the C-terminals are quite divergent. Previous studies showed that Lac2A was definitely up-regulated in the cuticle, but not Lac2B [34, 38]. Although the expression pattern and abundance were different in vivo, no significant differences in substrate preferences were observed between the recombinant Lac2A and Lac2B . In this study, two alternative splicing forms of Laccase 2 were also predicted in the Anopheles sinensis genome (Additional file 5: Figure S1; Additional file 6: Dataset 1). The phylogenetic analysis revealed that insect LAC2s were grouped into two distinct clads (Additional file 5: Figure S1). Lac2As are highly conserved with high bootstrap supports, but Lac2Bs are quite divergent (Additional file 5: Figure S1). This result implies that Lac2B may be more diverse in functions. However, no transcripts of AsLac2B were detected in our analysis; this may be because their expression is undetectable or they are dominantly expressed in other development stages or tissues. Although the dsRNA fragment was designed according to the common region of Lac2A and predicted Lac2B, no expression of Lac2B was detected in the high expression period and the tissues selected for Lac2A silencing. Moreover, the defective phenotype caused by RNAi was due to the dramatic downregulation of Lac2A expression, which is in good agreement with previous reports [38, 40–42, 44]. Therefore, we concluded that the defective phenotype was caused by Lac2A, but not Lac2B. Further studies will be conducted to investigate the expression pattern of Lac2B in other development stages (such as the adult stages) and more fine tissues, which can further help to understand its functions.
Melanin generated from the precursors of dopa and dopamine by the activated phenol oxidases (POs) is important for immune responses [5, 7–9, 12, 14–16, 19, 23, 24, 55, 56]. Genes (TH, DDC etc.) required for melanin synthesis are distributed in immune tissues including fat body and hemolymph [12, 57, 58]. Thus, immune tissues can utilize these precursors in melanization immune reactions catalyzed by phenol oxidases, such as tyrosinase or prophenoloxidases. These reactions are crucial for pathogen resistance. Laccase 2, a member of the phenol oxidase family, can catalyze the precursors to form melanin. If this gene can be expressed in some other immune tissues besides the cuticle of invertebrates [31, 34], theoretically, it can catalyze the terminal oxidation reaction for melanin synthesis. Our results revealed that AsLac2 was expressed in the pupal immune tissues (Additional file 7: Figure S2); this is similar to Laccase 2 expression being slightly increased in Anopheles gambiae after bacterial infection . Combined with our results showing that the melanization in the AsLac2-silenced pupae was significantly impaired and the resistance to bacteria challenges was dramatically decreased (Fig. 5), we notice the close relationship between the expression of AsLac2 and melanization immune responses. We speculate that AsLac2 may exert prophenoloxidase activity and involve melanotic immune responses. The dysfunction of AsLac2 can result in insufficient melanin synthesis and further impair immune responses to bacterial pathogens. Additionally, it is unclear whether immune responses are impaired by the feedback inhibition to catecholamine synthase caused by the temporarily accumulation of catecholamine, which is due to the silencing of AsLac2. This feedback inhibition may reduce the efficiency of melanin metabolism, in which the formation of the melanin precursors is further affected. To test this hypothesis, in vitro assays will be conducted to test the activity of total phenoloxidases in melanization and to determine the catecholamine content in dsLac2-injected pupae in subsequent research.
To our knowledge, the present results are the first to illustrate that the suppression of Laccase 2 results in serious adverse effects, which are almost lethal for the wild mosquitoes during the pupal development of Anopheles sinensis. As a rate-limiting enzyme in the final step of melanin metabolism, Laccase 2 has a similar expression pattern as the initial rate-limiting enzyme gene TH in Anopheles sinensis . Silencing these two genes produced similar physiological defective phenotypes, further suggesting that Laccase 2 and TH are crucial for normal pupal development. Our studies further corroborated that melanin metabolism is indispensable for the normal development of mosquito pupae. Therefore, the key genes and/or their regulatory elements in the melanin metabolism pathway may be valuable targets for mosquito prevention and control.
Primers used in this study. (XLSX 11 kb)
Amino acid sequence identity of Cu-oxidase domains of LAC2 orthologs. (PDF 109 kb)
Prediction of cis-acting transcriptional regulatory elements upstream. 7.5 kb of Laccase 2. (XLSX 9 kb)
Statistical analysis of pupal cuticle tanning degree in the RNAi experiment at 38 h after pupation. (DOC 28 kb)
Alternative splicing and amino acid sequence analyses of Laccase 2 in representative insect species. a Predicted alternative splicing forms of AsLac2. Red and green boxes represent the special exon of Laccase 2A and Laccase 2B, respectively. b Maximum Likelihood phylogenetic tree of two Laccase 2 forms in different insect species. c Genetic distance estimations among LAC2As (red) and LAC2Bs (yellow). (PDF 97 kb)
Dataset 1. The identified Laccase 2A and the predicted Laccase 2B isoform in Anopheles sinensis. (DOCX 15 kb)
Gene expression patterns of AsLac2 in pupal cuticle, fat body, and hemolymph. AsRPL49 was used as the internal control. (PDF 72 kb)
We would like to thank Dr. Gao Xiong (State Key Laboratory of Silkworm Genome Biology, Southwest University, Chongqing, China) for providing guidance in the frozen cuticle section.
This work was supported by National Natural Science Foundation of China (No. 31302038, 31672363, and 31372265), Science and Technology Research Project of Chongqing Municipal Education Commission (No. KJ1500327), Par-Eu Scholars Program (20136666), Coordinated Research Project of the International Atomic Energy Agency (18268), and National Key Program of Science and Technology Foundation Work of China (2015FY210300).
The data supporting the conclusions of this article are included within the article and the additional files. The sequence was submitted to the GenBank database under accession number KY132102.
LQ and BC conceived and designed the experiments; BC supplied the mosquito; MHD reared the mosquito; LQ, MHD, and FLS performed the experiments; LQ, MHD, ZWY, YJH, and ZTY analyzed the data; LQ and YJH wrote and revised the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Ming-Hui Du, Email: moc.kooltuo@hmdsumiced.
Zheng-Wen Yan, Email: moc.621@newgnehznay.
You-Jin Hao, Email: moc.liamtoh@nijuoyoah.
Zhen-Tian Yan, Email: moc.361@813tzy.
Feng-Ling Si, Email: moc.liamg@7121gnilgnefis.
Bin Chen, Email: email@example.com.
Liang Qiao, Email: moc.361@uwsgnailoaiq.