Protoplast-based engineering has become an important tool for basic plant molecular biology research and developing genome-edited crops. Uncaria rhynchophylla is a traditional Chinese medicinal plant with a variety of pharmaceutically important indole alkaloids. In this study, an optimized protocol for U. rhynchophylla protoplast isolation, purification, and transient gene expression was developed. The best protoplast separation protocol was found to be 0.8 M D-mannitol, 1.25% Cellulase R-10, and 0.6% Macerozyme R-10 enzymolysis for 5 h at 26 °C in the dark with constant oscillation at 40 rpm/min. The protoplast yield was as high as 1.5 × 10⁷ protoplasts/g fresh weight, and the survival rate of protoplasts was greater than 90%. Furthermore, polyethylene glycol (PEG)-mediated transient transformation of U. rhynchophylla protoplasts was investigated by optimizing different crucial factors affecting transfection efficiency, including plasmid DNA amount, PEG concentration, and transfection duration. The U. rhynchophylla protoplast transfection rate was highest (71%) when protoplasts were transfected overnight at 24 °C with the 40 µg of plasmid DNA for 40 min in a solution containing 40% PEG. This highly efficient protoplast-based transient expression system was used for subcellular localization of transcription factor UrWRKY37. Finally, a dual-luciferase assay was used to detect a transcription factor promoter interaction by co-expressing UrWRKY37 with a UrTDC-promoter reporter plasmid. Taken together, our optimized protocols provide a foundation for future molecular studies of gene function and expression in U. rhynchophylla.

E. rutaecarpa var. officinalis is a traditional Chinese medicinal plant known for its therapeutic effects, which encompass the promotion of digestion, the dispelling of cold, the alleviation of pain, and the exhibition of anti-inflammatory and antibacterial properties. The principal active component of this plant, limonin, is a potent triterpene compound with notable pharmacological activities. Despite its significance, the complete biosynthesis pathway of limonin in E. rutaecarpa var. officinalis remains incompletely understood, and the underlying molecular mechanisms remain unexplored. The main purpose of this study was to screen the reference genes suitable for expression analysis in E. rutaecarpa var. officinalis, calculate the expression patterns of the genes in the limonin biosynthesis pathway, and identify the relevant enzyme genes related to limonin biosynthesis. The reference genes play a pivotal role in establishing reliable reference standards for normalizing the gene expression data, thereby ensuring precision and credibility in the biological research outcomes. In order to identify the optimal reference genes and gene expression patterns across the diverse tissues (e.g., roots, stems, leaves, and flower buds) and developmental stages (i.e., 17 July, 24 August, 1 September, and 24 October) of E. rutaecarpa var. officinalis, LC-MS was used to analyze the limonin contents in distinct tissue samples and developmental stages, and qRT-PCR technology was employed to investigate the expression patterns of the ten reference genes and eighteen genes involved in limonin biosynthesis. Utilizing a comprehensive analysis that integrated three software tools (GeNorm ver. 3.5, NormFinder ver. 0.953 and BestKeeper ver. 1.0) and Delta Ct method alongside the RefFinder website, the best reference genes were selected. Through the research, we determined that Act1 and UBQ served as the preferred reference genes for normalizing gene expression during various fruit developmental stages, while Act1 and His3 were optimal for different tissues. Using Act1 and UBQ as the reference genes, and based on the different fruit developmental stages, qRT-PCR analysis was performed on the pathway genes selected from the “full-length transcriptome + expression profile + metabolome” data in the limonin biosynthesis pathway of E. rutaecarpa var. officinalis. The findings indicated that there were consistent expression patterns of HMGCR, SQE, and CYP450 with fluctuations in the limonin contents, suggesting their potential involvement in the limonin biosynthesis of E. rutaecarpa var. officinalis. This study lays the foundation for further research on the metabolic pathway of limonin in E. rutaecarpa var. officinalis and provides reliable reference genes for other researchers to use for conducting expression analyses.

Loranthus (Taxillus chinensis) is a facultative, hemiparasite and stem parasitic plant that attacks other plants for living. Transcriptome sequencing and bioinformatics analysis were applied in this study to identify the gene expression profiles of fresh seeds (CK), baby (FB), and adult haustoria tissues (FD). We assembled 160,571 loranthus genes, of which 64,926, 35,417, and 47,249 were aligned to NR, GO, and KEGG pathway databases, respectively. We identified 14,295, 15,921, and 16,402 genes in CK, FB, and FD, respectively. We next identified 5,480 differentially expressed genes (DEGs) in the process, of which 258, 174, 81, and 94 were encoding ribosomal proteins (RP), transcription factors (TF), ubiquitin, and disease resistance proteins, respectively. Some DEGs were identified to be upregulated along with the haustoria development (e.g., 68 RP and 26 ubiquitin genes). Notably, 36 RP DEGs peak at FB; 10 ER, 5 WRKY, 6 bHLH, and 4 MYB TF genes upregulated only in FD. Further, we identified 4 out of 32 microRNA genes dysregulated in the loranthus haustoria development. This is the first haustoria transcriptome of loranthus, and our findings will improve our understanding of the molecular mechanism of haustoria.

The hemiparasitic Taxillus chinensis (DC.) Danser is a root-parasitizing medicinal plant with photosynthetic ability, which is lost in other parasitic plants. However, the cultivation and medical application of the species are limited by the recalcitrant seeds of the species, and even though the molecular mechanisms underlying this recalcitrance have been investigated using transcriptomic and proteomic methods, genome resources for T. chinensis have yet to be reported. Accordingly, the aim of the present study was to use nanopore, short-read, and high-throughput chromosome conformation capture sequencing to construct a chromosome-level assembly of the T. chinensis genome. The final genome assembly was 521.90 Mb in length, and 496.43 Mb (95.12%) could be grouped into nine chromosomes with contig and scaffold N50 values of 3.80 and 56.90 Mb, respectively. In addition, a total of 33,894 protein-coding genes were predicted, and gene family clustering identified 11 photosystem-related gene families, thereby indicating photosynthetic ability, which is a characteristic of hemiparasitic plants. This chromosome-level genome assembly of T. chinensis provides a valuable genomic resource for elucidating the genetic basis underlying the recalcitrant characteristics of T. chinensis seeds and the evolution of photosynthesis loss in parasitic plants.

Uncaria rhynchophylla (Miq.) Miq. ex Havil, a traditional medicinal herb, is enriched with several pharmacologically active terpenoid indole alkaloids (TIAs). At present, no method has been reported that can comprehensively select and evaluate the appropriate reference genes for gene expression analysis, especially the transcription factors and key enzyme genes involved in the biosynthesis pathway of TIAs in U. rhynchophylla. Reverse transcription quantitative PCR (RT-qPCR) is currently the most common method for detecting gene expression levels due to its high sensitivity, specificity, reproducibility, and ease of use. However, this methodology is dependent on selecting an optimal reference gene to accurately normalize the RT-qPCR results. Ten candidate reference genes, which are homologues of genes used in other plant species and are common reference genes, were used to evaluate the expression stability under three stress-related experimental treatments (methyl jasmonate, ethylene, and low temperature) using multiple stability analysis methodologies. The results showed that, among the candidate reference genes, S-adenosylmethionine decarboxylase (SAM) exhibited a higher expression stability under the experimental conditions tested. Using SAM as a reference gene, the expression profiles of 14 genes for key TIA enzymes and a WRKY1 transcription factor were examined under three experimental stress treatments that affect the accumulation of TIAs in U. rhynchophylla. The expression pattern of WRKY1 was similar to that of tryptophan decarboxylase (TDC) under ETH treatment. This research is the first to report the stability of reference genes in U. rhynchophylla and provides an important foundation for future gene expression analyses in U. rhynchophylla. The RT-qPCR results indicate that the expression of WRKY1 is similar to that of TDC under ETH treatment. It may coordinate the expression of TDC, providing a possible method to enhance alkaloid production in the future through synthetic biology.

Taxillus chinensis (DC.) Danser, a parasitic plant that belongs to the Loranthaceae family, has a long history of being used in the Chinese medicine. We observed that the loranthus seeds were sensitive to temperature and could lose viability below 0°C quickly. Thus, we performed small RNA sequencing to study the microRNA (miRNA) regulation in the loranthus seeds under cold stress. In total, we identified 600 miRNAs, for the first time, in the loranthus seeds under cold stress. Then, we detected 224, 229, and 223 miRNAs (TPM > 1) in A0 (control), A1 (cold treatment for 12 h at 0°C), and A2 (cold treatment for 36 h at 0°C), respectively. We next identified 103 differentially expressed miRNAs (DEmiRs) in the loranthus seeds in response to cold. Notably, miR408 was upregulated during the cold treatment, which can regulate genes encoding phytocyanin family proteins and phytophenol oxidases. Some DEmiRs were specific to A1 and may function in early response to cold, such as gma‐miR393b‐3p, miR946, ath‐miR779.2‐3p, miR398, and miR9662. It is interesting that ICE3, IAA13, and multiple transcription factors (e.g., WRKY and CRF4 and TCP4) regulated by the DEmiRs have been reported to respond cold in other plants. We further identified 4, 3, and 4 DEmiRs involved in the pathways “responding to cold,” “responding to abiotic stimulus,” and “seed development/germination,” respectively. qRT‐PCR was used to confirm the expression changes of DEmiRs and their targets in the loranthus seeds during the cold treatment. This is the first time to study cold‐responsive miRNAs in loranthus, and our findings provide a valuable resource for future studies.