Genome-Wide Identification and Abiotic-Stress-Responsive Expression of CKX Gene Family in Liriodendron chinense

Liriodendron chinense is a tree species of the Magnoliaceae family, an ancient relict plant mainly used for landscaping and timber production due to its excellent material properties and ornamental value. The cytokinin oxidase/dehydrogenase (CKX) enzyme regulates cytokinin levels and plays an important role in plant growth, development, and resistance. However, too-high or too-low temperatures or soil drought can limit the growth of L. chinense, representing a key issue for research. Here, we identified the CKX gene family in the L. chinense genome and examined its transcriptional responses to cold, drought, and heat stresses. A total of five LcCKX genes, distributed on four chromosomes and divided into three phylogenetic groups, were identified across the whole L. chinense genome. Further analysis showed that multiple hormone- and stress-responsive cis-acting elements are located in the promoter regions of LcCKXs, indicating a potential role of these LcCKXs in plant growth, development, and response to environmental stresses. Based on existing transcriptome data, LcCKXs, especially LcCKX5, were found to transcriptionally respond to cold, heat, and drought stresses. Furthermore, quantitative reverse-transcription PCR (qRT-PCR) showed that LcCKX5 responds to drought stress in an ABA-dependent manner in stems and leaves and in an ABA-independent manner in roots. These results lay a foundation for functional research on LcCKX genes in the resistance breeding of the rare and endangered tree species of L. chinense.


Introduction
Cytokinins (CKs) are a group of plant hormones that play an important role in all aspects of plant growth and development, including apical dominance, stem or root branching, leaf spreading, lateral shoot growth, photosynthesis, seed germination, flower transition, and leaf senescence [1,2]. Cytokinins play opposite roles in shoot apical meristems (SAMs) and root apical meristems (RAMs). Plants with lowered cytokinin signal output or cytokinin content form larger RAMs and exhibit more rapidly growing roots [3,4]. Plant CKs promote developmental and physiological processes, drought tolerance, root architecture, and ultimately, crop productivity. CKs are degraded by cytokinin oxidases/dehydrogenases (CKXs), and CKs are precisely maintained in dynamic CK balance in nutritional tissues with the developmental regulation of the CKX-led irreversible degradation of catabolic CKs [5]. Studies have shown that the enhanced expression of CKX in roots to enhance cytokinin degradation leads to the formation of larger root systems in Arabidopsis, barley, oilseed rape, and rice [6].
The CKX gene family has been studied in a variety of plants, such as Arabidopsis, rice, wheat, tobacco, cotton, maize, soybean, chickpea, and others. For example, 7 CKXs (AtCKX1-AtCKX7) have been identified in Arabidopsis; a total of 13 CKXs

Phylogenetic Analysis of LcCKX Proteins
To further explore the phylogenetic relationships and evolutionary history of the LcCKX gene family, we constructed a phylogenetic tree containing CKX proteins from the basal angiosperm Amborella trichopoda, the magnoliophyte L. chinense, the monocotyledons Oryza sativa and Zea mays, and the dicotyledon Arabidopsis thaliana. A total of 44 CKXs were used to construct a phylogenetic evolutionary tree using the maximum-likelihood (ML) method, leading to three main phylogenetic groups (I, II, and III), two of which could be further divided into two subgroups (IIa and IIb, and IIIa and IIIb) ( Figure 1). Among them, group I comprises only three AmCKXs, while the rest of the groups contain CKXs from all five examined species. Specifically, LcCKX1 and -2 were grouped in IIIb, and LcCKX3, -4, and -5 were grouped in II ( Figure 1).

Analysis of Structure and Conserved Motifs of LcCKX Genes
Based on the analysis of the gene structure (Figure 2), we found that LcCKX1, LcCKX2, and LcCKX5 have four exons, while LcCKX3 and LcCKX4 have six and five exons, respectively. The identification of conserved motifs showed that all five LcCKXs protein sequences contain motifs 1, 3, 4, 6, 8, and 9 in the same order, hinting at the conserved function of these LcCKX proteins. Meanwhile, motifs 2 and 7 were only absent in LcCKX2 and -7, respectively, while motif 5 was present in LcCKX3-5, implying the functional differentiation of these LcCKXs.

Chromosome Localization and Tertiary Structure of LcCKX Genes
The chromosomal localization analysis showed that members of the LcCKX gene family are unevenly distributed on 13 L. chinense chromosomes (Figure 2a). Four CKX genes are distributed on three chromosomes, i.e., Chr2, -6, and -7, and the remaining one is located on a scaffold. The tertiary structure of a protein is the three-dimensional conformation of a protein molecule in its natural folded state. The tertiary structure is formed by further coiling and folding on the basis of the secondary structure, and the characterization of the tertiary structure of proteins is of great importance for the study  Figure 3b) show that LcCKX1 and LcCKX2 have a heart-shaped structure and a smaller mass. The LcCKX3, LcCKX4, and LcCKX5 models are more similar in structure and have a butterfly shape, and the protein structure morphologies show a high degree of similarity. These proteins all have one or several grooved structures. The grooved structures may provide the conditions for them to bind substrates to exert enzymatic activity, and the center of enzymatic activity may be located in these structures.

Analysis of Structure and Conserved Motifs of LcCKX Genes
Based on the analysis of the gene structure ( Figure 2), we found that LcCKX1 LcCKX2, and LcCKX5 have four exons, while LcCKX3 and LcCKX4 have six and five ex ons, respectively. The identification of conserved motifs showed that all five LcCKXs pro tein sequences contain motifs 1,3, 4, 6, 8, and 9 in the same order, hinting at the conserved function of these LcCKX proteins. Meanwhile, motifs 2 and 7 were only absent in LcCKX2

Chromosome Localization and Tertiary Structure of LcCKX Genes
The chromosomal localization analysis showed that members of the LcCKX gene family are unevenly distributed on 13 L. chinense chromosomes (Figure 2a). Four CKX genes are distributed on three chromosomes, i.e., Chr2, -6, and -7, and the remaining one is located on a scaffold. The tertiary structure of a protein is the three-dimensional conformation of a protein molecule in its natural folded state. The tertiary structure is formed by further coiling and folding on the basis of the secondary structure, and the characterization of the tertiary structure of proteins is of great importance for the study of the functional properties of proteins. The predictions of the 3D structures of LcCKXs (Figure 3b) show that LcCKX1 and LcCKX2 have a heart-shaped structure and a smaller mass. The LcCKX3, LcCKX4, and LcCKX5 models are more similar in structure and have a butterfly shape, and the protein structure morphologies show a high degree of similarity. These proteins all have one or several grooved structures. The grooved structures may provide the conditions for them to bind substrates to exert enzymatic activity, and the center of

Prediction of Cis-Acting Elements in LcCKX Promoter Regions
The cis-acting elements in the gene promoter region can be bound by specific transcription factors, thus regulating the expression of downstream genes. The prediction of cis-acting elements suggested that there are three main classes in the LcCKX family (Table 2), i.e., phytohormone signaling, environmental stress, and MYB binding sites. All LcCKX promoters contain abscisic acid (ABA) response elements, indicating that the LcCKX gene family might be sensitive to ABA. However, part of the LcCKX gene family also potentially responds to other phytohormones, such as ZT, MeJA, GA, and SA, indicating that these LcCKXs might be extensively involved in plant growth and development. Furthermore, most LcCKXs contain cis-elements that are related to cold and/or drought stress, implying a potential function of these LcCKXs in abiotic stress response.

Prediction of Cis-Acting Elements in LcCKX Promoter Regions
The cis-acting elements in the gene promoter region can be bound by specific t scription factors, thus regulating the expression of downstream genes. The predictio cis-acting elements suggested that there are three main classes in the LcCKX family (T 2), i.e., phytohormone signaling, environmental stress, and MYB binding sites. All Lc promoters contain abscisic acid (ABA) response elements, indicating that the LcCKX g family might be sensitive to ABA. However, part of the LcCKX gene family also po tially responds to other phytohormones, such as ZT, MeJA, GA, and SA, indicating these LcCKXs might be extensively involved in plant growth and development. Furt more, most LcCKXs contain cis-elements that are related to cold and/or drought st implying a potential function of these LcCKXs in abiotic stress response.

Gene Expression Pattern Analysis of LcCKX Family
The transcriptome sequence data of L. chinense were downloaded from the NCBI

Gene Expression Pattern Analysis of LcCKX Family
The transcriptome sequence data of L. chinense were downloaded from the NCBI SRA database, and the expression level of the genes was obtained using salmonid analysis. We determined the expression patterns of the LcCKX gene family in leaves under drought, cold (4 • C), and heat (40 • C) stresses ( Figure 4a). LcCKX5 showed a strong transcriptional response under all three stresses. Specifically, the expression of LcCKX5 was up-regulated and peaked after 12 h and 3 d in response to heat and drought stresses, respectively. In contrast, the expression level of LcCKX5 decreased under cold stress. In comparison, LcCKX1, -3, and -4 responded to these abiotic stresses, but less strongly, while LcCKX2 seemed not to be expressed at all. The above results indicate that LcCKX genes differentially responded to different stresses and that LcCKX5 might be a valuable target for further research on the resistance breeding of L. chinense. response under all three stresses. Specifically, the expression of LcCKX5 was up-regu and peaked after 12 h and 3 d in response to heat and drought stresses, respective contrast, the expression level of LcCKX5 decreased under cold stress. In compa LcCKX1, -3, and -4 responded to these abiotic stresses, but less strongly, while Lc seemed not to be expressed at all. The above results indicate that LcCKX genes diff tially responded to different stresses and that LcCKX5 might be a valuable target fo ther research on the resistance breeding of L. chinense. To elucidate the expression pattern of LcCKXs in the growth and developmen chinense, we constructed gene expression profiles for different stages of the somati bryogenesis of L. chinense (Figure 4b). LcCKX2, LcCKX3, and LcCKX5 were found involved in the process from embryonic callus to regeneration plantlet; LcCKX3 sh a transcriptional association during somatic embryo development, while LcCKX5 m be involved in somatic embryo maturation.
Meanwhile, when comparing the expression of the LcCKX gene family in diff organs of L. chinense (Figure 4c), the analysis showed that LcCKX1 was highly expr in different organs of the petals, especially in the shoot apex, while LcCKX2 was expressed in bracts and stamens than in other parts. In contrast, LcCKX3 and LcCKX low or no expression in different tissues. The above results suggest that LcCKX gen extensively involved in plant growth and development, as well as responses to a stresses, in L. chinense.

Expression Analysis of LcCKX5 under Abiotic Stress
To further investigate the expression pattern of LcCKX5 under drought stress PEG8000 was used to treat L. chinensis seedlings to simulate drought conditions, an To elucidate the expression pattern of LcCKXs in the growth and development of L. chinense, we constructed gene expression profiles for different stages of the somatic embryogenesis of L. chinense (Figure 4b). LcCKX2, LcCKX3, and LcCKX5 were found to be involved in the process from embryonic callus to regeneration plantlet; LcCKX3 showed a transcriptional association during somatic embryo development, while LcCKX5 might be involved in somatic embryo maturation.
Meanwhile, when comparing the expression of the LcCKX gene family in different organs of L. chinense (Figure 4c), the analysis showed that LcCKX1 was highly expressed in different organs of the petals, especially in the shoot apex, while LcCKX2 was more expressed in bracts and stamens than in other parts. In contrast, LcCKX3 and LcCKX4 had low or no expression in different tissues. The above results suggest that LcCKX genes are extensively involved in plant growth and development, as well as responses to abiotic stresses, in L. chinense.

Expression Analysis of LcCKX5 under Abiotic Stress
To further investigate the expression pattern of LcCKX5 under drought stress, 20% PEG8000 was used to treat L. chinensis seedlings to simulate drought conditions, and the transcriptional response of LcCKX5 was quantified using qRT-PCR. The internal reference genes and related primers are shown in the Table 3. The results show that the expression change of LcCKX5 differed between roots and shoots in response to drought stress. Specifically, under drought stress, the expression of LcCKX5 decreased and then increased over time in roots, while there were no significant changes in pairs of stems, and there was a gradual decrease in leaves (  Table 3. qRT-PCR primers used to quantify LcCKX5 gene expression.

Gene Name qRT-PCR Primer
LcCKX5 transcriptional response of LcCKX5 was quantified using qRT-PCR. The internal reference genes and related primers are shown in the Table 3. The results show that the expression change of LcCKX5 differed between roots and shoots in response to drought stress. Specifically, under drought stress, the expression of LcCKX5 decreased and then increased over time in roots, while there were no significant changes in pairs of stems, and there was a gradual decrease in leaves (  Table 3. qRT-PCR primers used to quantify LcCKX5 gene expression.

Discussion
Cytokinin dehydrogenase regulates the content of cytokinin and plays an important role in balancing the synthesis and degradation of cytokinin in plants [19,20]. CKX plays a very important role in balancing the synthesis and degradation of cytokinins in plants.
With the continuous development of plant genomics and reverse genetics, we can provide powerful tools to study CKX genes [21]. The decipherment of the genome of L. chinensis, a representative species of the Liriodendron genus in the Magnoliophyta family, provides new insights into the phylogenetic position of magnoliophytes in angiosperms [22]. CKX gene families have been identified or deduced in Arabidopsis, rice, wheat (Triticum aestivum L.), maize (Zea mays L.), cotton (Gossypium hirsutum L.), and alfalfa (Medicago sativa L.) [23][24][25]. The CKX gene family was identified or deduced in different plants, such as Arabidopsis, maize, cotton, and alfalfa. There are 7 AtCKX genes in Arabidopsis, 13 ZmCKXs in maize, 11 OsCKXs in rice, 8 VvCKXs in grape (Vitis vinifera L.), and 5 MnCKXs in mulberry (Morus notabilis) [7,[26][27][28][29]. In this study, five LcCKX genes were identified in L. chinensis. The number of genes in other species was small, which showed that CKX is a gene family with a small number of members. The five LcCKX genes identified fall into three subgroups, and the genes in each subgroup are highly conserved in terms of protein physicochemical properties, conserved motifs, and gene structure.
The identification of members of a species-specific gene family using the conserved structural domains of a particular gene family is an accurate and efficient method. Analysis of the conserved structural domains showed that most of them are relatively conserved in the family genes, and only some of them show some specificity; for example, motif 5 is only present in group II (LcCKX3~5) and motif 10 is only present in group IIa (LcCKX3~4). This specificity may be related to the functions assumed by different groups of CKX and to the differences in the higher structure of the proteins. Analysis of the 3D structure of the proteins showed that the butterfly-shaped protein structure corresponds to three members of group II, while the protein with a heart-shaped tertiary structure corresponds to two members of group I.
Cis-regulatory elements are specific motifs located in the promoter region of genes that act as binding sites for genes and play an important role in stress response by regulating the transcription of downstream genes [30]. In addition, some phytohormones, such as abscisic acid, salicylic acid, jasmonic acid, and ethylene, are also involved in regulating the adaptive response of plants to abiotic stresses [31]. According to the results of the predicted cis-acting elements, the LcCKX gene mainly responds to plant hormones, light, and abiotic stress. In other plant species, such as in maize, these cis-acting elements are also present in large numbers in the promoter regions of the CKX gene [32]. These results further suggest that the LcCKX gene may be important and involved in the plant response to abiotic stress. As the promoter regions of the five LcCKXs contain different types, numbers, and distributions of cis-acting elements, but also share the same characteristics, it is assumed that the different LcCKXs have different expression patterns.
The cytokinin group of plant hormones is involved in regulating several aspects of plant growth and development [33], many of which have direct effects on crop improvement, such as the regulation of the root crown structure [34], the regulation of inflorescence meristem tissue activity and seed yield [35,36], the regulation of leaf senescence and photosynthesis, and response to biotic and abiotic stresses [27,37,38]. According to the RNA-Seq data, under abiotic stress conditions such as drought, low-temperature treatments, and high-temperature treatments, the expression of some LcCKX genes can be changed, thus playing a defensive and protective role in plant growth and development processes, which may also be related to the resistance of goosefoot, but the exact mechanism of this effect needs to be studied in the future.
In other species, CKXs have diverse expression patterns, even when they are tandemly repeated. In maize, for example, ZmCKX1, ZmCKX6, and ZmCKX10 are expressed in all tissues, while ZmCKX7 and ZmCKX8 are only expressed in male ears and are clearly tissue-specific. ZmCKXs differ in number and expression in different tissues under the same stress conditions [32]. In Plasmodiophora brassicae, CKXs show different expression patterns in different tissues, indicating different functions [39]. Under normal conditions, three of these genes were highly expressed in all tissues studied, which is characteristic of the expression pattern of housekeeping genes [39].
Among tree species, CKXs have different expression patterns in different tissues. In studies related to CKX in Jatropha curcas, it was found that JcCKX1 is mainly expressed in flower buds, roots, and female flowers, while JcCKX2 shows very strong expression in female flowers and seeds. JcCKX3 is highly expressed in male flowers, and JcCKX4 shows high expression levels in mature leaves and female flowers, and extremely high expression in seeds. JcCKX5 is mainly expressed in stems, tender leaves, and fruits, while JcCKX6 is expressed in all tissues, and JcCKX7 is mainly expressed in roots [40]. In Malus domestica, the expression of CKX is significantly higher in leaves than in other tissues; MdCKX1, MdCKX7, MdCKX9, and MdCKX11/12 have high levels of expression in roots. In addition, compared with other tissues, all MdCKX levels are lower in the stem and axillary buds, indicating that CK accumulation is naturally limited to the stem and axillary buds. At the same time, Ming Tan et al. treated the experimental materials with 6-BA, and throughout the sampling time, MdCKX1, MdCKX2, MdCKX5, and MdCKX10 were significantly up-regulated, especially MdCKX10, whose transcription level increased by more than 100 times in 48 h. The MdCKX2, MdCKX4, MdCKX6, MdCKX7, and MdCKX8 transcripts did not respond to 6-BA treatment, and their expression levels were similar to those in untreated control buds after 24 and 48 h [41]. Combined with the relevant results of this study, it can be found that the expression pattern of CKX in different tissues of trees is similar and that there are also different expression locations. In general, the expression amount of CKX in roots is relative, while under non-biotic stress, the response of CKX in roots is faster and the response in leaves and stems is slower.
In summary, this research focused on LcCKX5 RNA-Seq data, namely the expression pattern of LcCKXs with respect to different tissues, time, and stress conditions, as well as functional predictions. The results could help to screen potential resistance genes and provide very detailed and reliable information for subsequent studies. 2 software was used to analyze and identify the protein sequences that contain both conserved structural domains. The candidate member sequences were compared using NCBI (https://www.ncbi.nlm.nih.gov/, accessed on 5 July 2022) to determine their CDS and position information, and the CDD database (https://www.ncbi.nlm.nih.gov/cdd/, accessed on 5 July 2022) was used to manually analyze and confirm that the candidate genes contained 2 conserved structural domains at the same time; finally, the members of the CKX gene family were obtained.

Protein Physicochemical Properties and Tertiary Structure Prediction of CKX Genes in L. chinense
According to the identified CKX gene ID, the chromosome location information was obtained from the GFF3 file. The physicochemical properties of the proteins, including relative molecular mass, theoretical isoelectric point, instability index, and hydrophilicity, were analyzed using the web-based online tool ExPASy (https://web.expasy.org/, accessed on 5 July 2022). The position information of the LcCKX genes on the chromosome was picked up from the L. chinense annotations using TBtools. The predictive analysis of the 3D structures of different members of the CKX gene family was performed according to the online website SWISS-MODEL (https://swissmodel.expasy.org/, accessed on 5 July 2022).

Phylogenetic Analysis of CKX Genes in L. chinense
MEGA v10.1.8 (Temple, Philadelphia, PA, USA) was used to examine CXKs from L. chinense, rice, and Arabidopsis to determine their phylogenetic relationship. We used MUSCLE implemented in MEGA to align the amino acid sequences and the maximumlikelihood estimation algorithm to create phylogenetic trees with a bootstrap value of 1000. DANMAN v9.0 (Lynnon Corporation, San Ramon, CA, USA) was used for the multi-fragment alignment of amino acid sequences.

Gene Structure and Conserved Motif Analysis of CKX Genes in L. chinense
The structure of the gene was produced based on the length of the CKX gene and the location information of the CDS. The conserved patterns of the aa sequences of all members of LcCKX were analyzed using MEME (http://memesuite.org/tools/meme/, accessed on 5 July 2022), and the phylogenetic tree, gene structure, and conserved patterns of the LcCKX gene were merged using TBtools.

Cis-Acting Elements of CKX Genes in L. chinensis
The genomic DNA sequences 2000 bp upstream of TAG in the sequences of different members of the LcCKX family were selected using whole-genome data of L. chinense, predicted and analyzed in cis using Plant Care (http://bioinformatics.psb.ugent.be/webtools/ plantcare/html/, accessed on 5 July 2022), and mapped and displayed using Origin-Pro9.06Bit software to predict gene function.

RNA-seq Analysis of LcCKX Gene Expression Levels in Different Organs and under Multiple Stresses
The drought stress transcriptome data of L. hybrid were annotated with accession number PRJNA679101 and can be downloaded from NCBI (https://www.ncbi.nlm.nih. gov/bioproject/PRJNA679101/; accessed on 5 September 2022). The data on somatic embryogenesis are unpublished. Heatmaps visualizing expressions were realized with TBtools software (v1.09). The log2 (TPM + 1) value was used for standardization and hierarchical cluster analysis.
Transcript data of different organs and hybrid data of L. chinense under high-temperature stress and drought stress were downloaded from NCBI. The transcript data on low temperature, L. chinense petal development, and hybrid L. chinense somatic embryogenesis are undisclosed data. The expression levels of related genes are listed in Table S2. All mRNA abundance values were measured in transcripts per million (TPM) based on the L. chinense genomic database.

Plant Materials and Abiotic Stress Treatment
Seedlings produced during the somatic embryogenesis of L. chinense were used as the starting material for this study. Prior to stress treatment, somatic embryogenic seedlings were removed from the culture vessels and domesticated in a greenhouse for 2 weeks (conditions: 22 • C, prolonged light for 16 h and darkness for 8 h, and relative humidity of 75%). Three biological replications were performed using 15% PEG8000 to simulate a natural drought environment. Roots, stems, and leaves were harvested after 6 h and 24 h of drought treatment, respectively. All experimental tissue samples were immediately frozen in liquid nitrogen and then stored at −80 • C.

RNA Extraction and Quantitative Real-Time PCR Analysis
The first-strand gene was synthesized from 1.0 mg of RNA using an Evo M-MLV RT kit with gDNA Clean (Changsha Precision Biotechnology (Hunan) Co.). An Equalbit 1× dsDNA HS Assay Kit (EQ121-01; Vazyme, Nanjing, China) was used to complete the quantification of all reverse cDNA. Polymerase chain reaction amplification was performed in 20 µL using a SYBR ® Green PreMix Pro Taq HS qPCR Kit (Precision Biotechnology (Hunan) Co., Ltd., Changsha, China) using a Roche LightCycler ® 480 Real-Time Polymerase Chain Reaction System. Three replicates were performed for each selected gene. The expression pattern of CKX-related genes under low-temperature stress was studied using qRT-PCR, and Primer 5.0 software was used to design primers in the non-conserved structural domain region of the genes. qPCR SYBR Green Master Mix (Vazyme) was used for real-time quantitative PCR, and GAPDH and 18s were used as the internal control genes. The real-time PCR cycling parameters were 95 • C for 30 s, followed by 45 cycles at 95 • C for 5 s and 60 • C for 30 s, with melting curve analysis. All reactions were performed in triplicate to ensure the repeatability of the results. Gene expression levels were calculated using 2 −∆∆Ct [42].

Conclusions
In this study, comprehensive analysis of the CKX gene family in the L. chinense genome was conducted, and five LcCKX genes were identified. Subsequently, gene structure analysis, phylogeny, chromosomal localization, gene duplication, and genome-wide identification and analysis of CKX family genes in L. chinense were carried out using bioinformatics and qRT-PCR. The differential expression of LcCKX5 genes in different tissues of L. chinense and the different expression trends under drought stress and ABA treatment may indicate that they play an important role in drought resistance and tissue development. The present study provides comprehensive information on the CKX genes in L. chinense, which could help to determine the functions of CKX genes.
Author Contributions: X.S. was the designer and performer of this study, completed data analysis, and wrote the first draft of the paper; L.Z. and W.W. assisted in some of the data analysis; L.X., Y.Y., Y.L., J.Z. and J.S. instructs the experimental operation; J.C. and Z.H. directed the experimental design, data analysis, paper writing, and full-text revision. All authors have read and agreed to the published version of the manuscript. Data Availability Statement: Transcriptome data of somatic embryogenesis and tissues have not yet been published. The drought stress transcriptome data of L. hybrid were annotated with accession number PRJNA679101 and can be downloaded from NCBI (https://www.ncbi.nlm.nih.gov/ bioproject/PRJNA679101/; accessed on 5 September 2022). The complete genome, transcript/protein sequences, and genome feature file of Lchi were downloaded from https://www.ncbi.nlm.nih.gov/ assembly/GCA_003013855.2 (accessed on 5 September 2022).