*4.1. Identification and Sequence Analysis of NAC Genes from L. indica*

As a reference sequence, the NAC protein sequence was acquired from the *Arabidopsis* TAIR database (https://www.arabidopsis.org/) (accessed on 17 March 2022) and members of the putative *LiNAC* genes were sought for using the *L. indica* transcriptome database (the raw sequence data were obtained from the website https://bigd.big.ac.cn/gsa) (accessed on 17 March 2022) and a BLASTP search. Then, we used the local HMMER3.0 software(Robert, D.F.; Ashburn, VA, USA) and the Hidden Markov (HMM) profile of the NAC protein (PF02365), which was downloaded from the Pfam database (http://pfam.xfam.org) (accessed on 17 March 2022). The potential gene members discovered using the two search techniques were pooled. The conserved NAM domain was discovered in all potential *LiNAC* genes using the Batch Web CD-Search Tool (https://www.ncbi.nlm.nih.gov/cdd) (accessed on 17 March 2022) and SMART (http://smart.embl.de/) (accessed on 17 March 2022), ensuring the correctness of the results.
