The Fritillaria species ranked as a well-known traditional medicine in China and has become rare due to excessive harvesting. To find reasonable strategy for conservation and cultivation, identification of new ecological distribution of Fritillaria species together with prediction of those responses to climate change are necessary. In terms of current occurrence records and bioclimatic variables, the suitable habitats for Fritillaria delavayi, Fritillaria taipaiensis, and Fritillaria wabuensis were predicted. In comparison with Maxent and GARP, Biomod2 obtained the best AUC, KAPPA and TSS values of larger than 0.926 and was chosen to construct model. Temperature seasonality was indicated to put the greatest influence on Fritillaria taipaiensis and Fritillaria wabuensis, while isothermality was of most importance for Fritillaria delavayi. The current suitable areas for three Fritillaria species were distributed in south-west China, accounting for approximately 17.72%, 23.06% and 20.60% of China's total area, respectively. During 2021–2100 period, the suitable habitats of F. delavayi and F. wabuensis reached the maximum under SSP585 scenario, while that of F. taipaiensis reached the maximum under SSP126 scenario. The high niche overlap among three Fritillaria species showed correlation with the chemical composition (P ≤ 0.05), while no correlation was observed between niche overlap and DNA barcodes, indicating that spatial distribution had a major influence on chemical composition in the Fritillaria species. Finally, the acquisition of species-specific habitats would contribute to decrease in habitat competition, and future conservation and cultivation of Fritillaria species.

Fritillaria species, a well-known Chinese traditional medicine for more than 2,000 years, have become rare resources due to excessive harvesting. In order to balance the economical requirement and ecological protection of Fritillaria species, it is necessary to determine (1) the important environmental variables that were responsible for the spatial distribution, (2) distribution change in response to climate change in the future, (3) ecological niche overlap between various Fritillaria species, and (4) the correlation between spatial distribution and phylogenies as well. In this study, the areas with potential ecological suitability for Fritillaria cirrhosa, Fritillaria unibracteata, and Fritillaria przewalskii were predicted using MaxEnt based on the current occurrence records and bioclimatic variables. The result indicated that precipitation and elevation were the most important environmental variables for the three species. Moreover, the current suitable habitats of F. cirrhosa, F. unibracteata, and F. przewalskii encompassed 681,951, 481,607, and 349,199 km², respectively. Under the scenario of the highest concentration of greenhouse gas emission (SSP585), the whole suitable habitats of F. cirrhosa and F. przewalskii reach the maximum from 2021 to 2100, while those of F. unibracteata reach the maximum from 2021 to 2100 under the scenario of moderate emission (SSP370) from 2021 to 2100. The MaxEnt data were also used to predict the ecological niche overlap, and thus high overlap occurring among three Fritillaria species was observed. The niche overlap of three Fritillaria species was related to the phylogenetic analysis despite the non-significance (P > 0.05), indicating that spatial distribution was one of the factors that contributed to the speciation diversification. Additionally, we predicted species-specific habitats to decrease habitat competition. Overall, the information obtained in this study provided new insight into the potential distribution and ecological niche of three species for the conservation and management in the future.

Drought stress is a major constraint on plant growth and agricultural productivity. Caffeic acid O-methyltransferase (COMT), an enzyme involved in the methylation of various substrates, plays a pivotal role in plant responses to abiotic stress. The involvement of COMTs in drought response, particularly through the enhancement of lignin and melatonin biosynthesis, remains poorly understood. In this study, LcCOMT was firstly proposed to be associated with the biosynthesis of both lignin and melatonin, as demonstrated through sequence comparison, phylogenetic analysis, and conserved motif identification. In vitro enzymatic assays revealed that LcCOMT effectively methylates N-acetylserotonin to melatonin, albeit with a higher Km value compared to caffeic acid. Site-directed mutagenesis of residues Phe171 and Asp269 resulted in a significant reduction in catalytic activity for caffeic acid, with minimal impact on N-acetylserotonin, underscoring the specificity of these residues in substrate binding and catalysis. Under drought conditions, LcCOMT expression was significantly upregulated. Overexpression of LcCOMT gene in Arabidopsis plants conferred enhanced drought tolerance, characterized by elevated lignin and melatonin levels, increased chlorophyll and carotenoid content, heightened activities of antioxidant enzymes peroxidase (POD), catalase (CAT), and superoxide dismutase (SOD), and reduced malondialdehyde (MDA) and hydrogen peroxide (H₂O₂) accumulation. This study is among the few to demonstrate that COMT-mediated drought tolerance is achieved through the simultaneous promotion of lignin and melatonin biosynthesis. LcCOMT represents the first functionally characterized COMT in Apiaceae family, and it holds potential as a target for genetic enhancement of drought tolerance in future crop improvement strategies.

The Plant-specific AP2/ERF gene family encodes proteins involved in various biological and physiological processes. Although the genome of Coptis chinensis Franch, a plant producing benzylisoquinoline alkaloids (BIAs), has been sequenced at the chromosome level, studies on the AP2/ERF gene family in C. chinensis are lacking. Thus, a genome-wide identification of AP2/ERF gene family in C. chinensis was conducted to explore its role in BIAs biosynthesis. A total of 96 CcAP2/ERF genes were identified and categorized into five subfamilies, including 43 ERFs, 32 DREBs, 17 AP2s, 3 RAVs, and 1 Soloist, based on their structural domains. These CcAP2/ERF genes were unevenly distributed across nine chromosomes. Analysis of gene duplication events identified 17 CcAP2/ERF gene pairs in the genome, with 7 involved in tandem duplication events and 10 involved in segmental duplicate events, indicating that both types of duplications contributed to the expansion of the AP2/ERF gene family. The Ka/Ks ratio analysis suggested that the CcAP2/ERF gene family underwent strong purifying selection. Two phytohormones, methyl jasmonate and abscisic acid, were identified as potential key inducers of BIAs biosynthesis due to the cis-acting element prediction. Analysis of the spatial transcriptomic data revealed that 28 differentially expressed AP2/ERF genes had the highest or relatively higher expression levels in the rhizome, 17 of which positively correlated with the tissue-specific accumulation of BIAs. Further real-time PCR verification and protein-protein interaction analysis indicated that DREB1B might be one of the central regulators in a highly complex BIAs biosynthesis network. These findings provide significant insight into the function of AP2/ERF genes in C. chinensis, particularly in the regulatory network of BIAs biosynthesis in C. chinensis. This study also identifies candidate genes for metabolic engineering to increase BIAs content in C. chinensis.