Application of DNA Barcoding to Identify Medicinal Plants

Summary

The present protocol describes the use of DNA barcoding technology to authenticate plant-based medicinal materials, and the medicinal plant Angelica sinensis (Oliv.) Diels was identified as an example.

Abstract

Medicinal plants are valuable resources globally and are used worldwide to maintain health and treat disease; however, the presence of adulteration obstructs their development. DNA barcoding, a technique for species identification by standard DNA regions, facilitates prompt and accurate identification of traditional medicinal plants. The process of DNA barcoding entails six basic steps: 1) processing the medicinal plants, 2) extracting high-quality total DNA from the medicinal plants using centrifugal column method, 3) amplifying target DNA region internal transcribed spacer 2 (ITS2) with universal primers of plants and performing Sanger sequencing, 4) splicing and aligning sequence to obtain the target sequence, 5) matching the barcode sequence against the barcode library for identification, 6) aligning sequence, comparing intraspecific and interspecific variation, constructing phylogenetic neighbor-joining tree. As shown in the results, the universal primer can amplify the target region. Basic Local Alignment Search Tool (BLAST) demonstrates the percentage identified was 100%, and the neighbor-joining tree demonstrates that the splicing sequences were clustered with the A. sinensis OR879715.1 clade, and the clade support value is 100. This protocol provides a reference for applying DNA barcoding technology as an effective method to identify medicinal plants and adulterants.

Introduction

Medicinal plants have a wide range of pharmacological effects and are important substances for the treatment and prevention of diseases. The market demand for medicinal plants in herbal medicines and pharmaceutical products is huge and still growing. With the increasing market for medicinal plants, the problem of adulteration has hindered the development of medicinal plants. Currently, medicinal plant adulteration suffers from these reasons: 1) the similar morphology of medicinal plants makes it difficult to identify and use them correctly^1,2,3,4,5, 2) the increasing demand for medicinal plants has led to insufficient supply in the market^6,7, 3) medicinal plants are expensive and have fluctuating prices, and the use of cheap herbs instead of economically valuable materials has led to adulteration and market profiteering^8,9. To solve the problems of proper identification and use of medicinal plants, there is a need for a technology that allows non-specialists to identify the source¹⁰.

There are limitations to properly identifying and using medicinal plants by appearance and odor alone¹¹. For more accurate identification and quality control, physicochemical methods are employed¹². Thin-layer chromatography (TLC), for instance, is a rapid method for identifying medicinal herbs and is included in the Chinese Pharmacopoeia. Nevertheless, it requires reference standards to identify the plants¹³. High-performance liquid chromatography (HPLC) can perform both qualitative and quantitative tests, making it suitable for medicine quality testing. However, this instrument is expensive and requires specialized knowledge to operate^14,15. DNA barcoding technology is a rapid and accurate species identification technology that identifies species by using stable DNA sequences to differentiate between plants with similar morphology. Hebert et al. proposed DNA barcoding technology to identify species, which uses a recognized and relatively short DNA sequence within the genome¹⁶, Chen et al. proposed ITS2 as a universal sequence for the molecular identification of medicinal plants and identified more than 6600 plants and found a high identification ability¹⁷. The study of medicinal plants based on ITS2 and chloroplast gene fragments has demonstrated the ability to distinguish species at the species level, with applications in the fields of resource protection, market regulation, and international trade^17,18,19. The application of DNA barcoding can overcome the limitations of traditional identification methods, which is helpful in ensuring the quality and safety of medicinal plants, preventing resource abuse and confusion²⁰. It has proven beneficial for studying medicinal plants, supporting the growth of the herbal industry in a sustainable manner, and providing a guarantee that the consumer uses the correct medication^21,22,23,24.

The paper outlines protocols for how to apply DNA barcoding to identify medicinal plants using A. sinensis as an example. DNA barcoding utilizes one or more standardized short genetic markers in an organism's DNA to recognize it as belonging to a particular species. Current methods for identifying medicinal plants have certain limitations. Traditional morphological identification relies heavily on experts' extensive experience, which can be influenced by human factors, while chemical identification is prone to issues like adulteration of key compounds. In contrast, DNA barcoding provides a more accurate means of species identification, offering advantages such as speed, high reproducibility, and stability. Furthermore, this technology enables centralized management and sharing of existing species sequence data through the internet and information platforms, significantly improving the efficiency and reliability of species identification^11,12. Due to their similar morphology and medicinal parts, A. sinensis is often confused with other species and is frequently misused or intentionally substituted on the market²⁵. Yang et al. identified A. anomala using DNA barcoding to distinguish it from false drugs²⁶. Yuan et al. collected 23 species of Angelica and used DNA barcoding to differentiate between the various species²⁷. By following the procedures described in this protocol, users will be able to authenticate plant-based medicinal materials from pre-processing to ultimately matching the barcoding sequence with the barcode database (Figure 1).

Protocol

1. Sample preparation

1. The present study utilized two-year-old A. sinensis with low temperature and long sunshine at an altitude of 2800 m from Min County, Gansu Province, as experimental materials.Three plants were selected from which more than three roots were selected for experimentation. Rinse the roots first with sterilized water and then clean them with an alcohol sprayer or alcohol swabs to prevent contamination that could affect the accuracy of the experimental results. Cut the roots into pieces smaller than 5 mm using a sterilized scalpel and mix for further analysis. (Figure 2).
   NOTE: Alternatively, place the roots in liquid nitrogen to crisp up, making them easier to break. Depending on the shape, hardness, and other characteristics of the medicinal plants, use instruments to deal with them.
2. Prepare three tubes of experimental samples for parallel experiments. Label each sample as AS-1, AS-2, and AS-3. Take about 100 mg fresh or 20 mg dry-weight tissue (air-dried) and place it in a 2 mL centrifuge tube.
3. Place two 5 mm stainless steel beads into the 2 mL centrifuge tubes. Transfer the tubes to a high-throughput tissue grinder operating at 60 Hz for 60 s.

2 DNA extraction from sample

NOTE: This study uses a Plant Genomic DNA Extraction Kit, which is based on the CTAB method. DNA extraction is performed following the instruction manual with some modifications, including the addition of a unique Nuclear Isolation Buffer (NIB). The protocol improves DNA purity by first separating the contaminants from the tissue and then adding an initial nuclei isolation step^28,29,30.

1. Add 800 µL of pre-cooled NIB and place the sample on a vortex mixer, vibrating for 5 min.
2. Place in a centrifuge and process at 13,780 x g for 10 min. Remove the supernatant.
3. Repeat steps 2.1-2.2 for an additional 3 rounds.
   NOTE: NIB extracts higher-quality DNA, and its preparation contains some hazardous chemicals that must be handled in a fume hood, making it more difficult for novices to handle. However, this step is not mandatory. Users may opt to follow this recommendation or skip these steps. Its detailed formulation is in Table 1.
4. Add 600 µL of Buffer P1 and 5 µL of RNase A (10 mg/mL). Place the sample on a vortex mixer and vortex to mix. Incubate the samples in a water bath at 65 °C for 1.5 h, inverting the tubes 2x-3x during incubation.
5. Add 200 µL of Buffer P2, mix thoroughly, place on ice for 15 min, and then centrifuge at 13,780 x g for 10 min.
6. Carefully aspirate the supernatant onto a filter column AF, centrifuge at 13,780 x g for 2 min, and transfer the lower filtrate collected in the collection tube to a new 1.5 mL centrifuge tube.
7. Calculate the volume of filtrate, add 1.5x the volume of P3, and shake immediately to get a mixed solution.
8. Add 650 µL of the prepared mixture to an adsorbent column adsorption column (AC) and incubate for 5 min, then centrifuge at 13,780 x g for 30-60 s.
9. Add the filtrate obtained from step 2.8 to the AC again, centrifuge, and discard the waste liquid. Repeat the procedure with the remaining mixed solution from step 2.7.
10. Add 600 µL of rinse solution WB, centrifuge at 13,780 x g for 30 s, discard the waste solution, and repeat the procedure 1x.
11. Place the column AC back into the empty collection tube, centrifuge at 13,780 x g for 2 min, and leave at room temperature to remove residual ethanol.
12. Take the adsorbent column AC and place it in a clean 1.5 mL centrifuge tube. Add 45 µL of sterilized deionized water (ddH₂O) to the middle of the adsorbent membrane, which has been heated to 65 °C. Leave at room temperature for 5 min. Centrifuge at 13,780 x g for 2 min.
13. Add the resulting solution back to the column and centrifuge at 13,780 x g for 2 min at room temperature.
14. Collect the filtered solution. Use a spectrophotometer to determine DNA concentration; OD260 represents the absorbance of nucleic acids. Use the following criteria for DNA quality determination standard: OD260/OD280=1.80-2.00, the normal range; OD260/OD280>2.00; there is a large amount of RNA contamination; OD260/OD280<1.80, there are proteins, phenols, and other contaminants.

3 Target fragment amplification and detection

NOTE: Select appropriate amplification primers based on the plant's proposed DNA barcoding regions; reaction conditions are shown in Table 2.

1. Set up the reaction conditions corresponding to primer amplification. Use 0.2 mL centrifuge tubes filled with 9.5 µL of ddH₂O, 12.5 µL of 2x PCR Master Mix, and 1 µL each of forward primer, reverse primer, and template DNA for a total volume of 25 µL. Prepare a negative control, replacing the template with water.
2. Add 20 mL of 1x Tris Acetate-EDTA (TEA) buffer and 0.2 g of agarose to heat and prepare a 1% agarose gel, then add 2 µL of nucleic acid gel stain (10,000x) and mix thoroughly. Pour the mixture into a gel casting tray and allow it to solidify to obtain the agarose gel. Add 3 µL of the PCR product and 2 µL of DNA marker 5000 bp to the gel.
3. Set the voltage at 120 V, the current at 400 mA, and electrophoresis for 40 min. Use a versatile gel image analysis system to view gels and check the results.

4. Data collection and analysis

NOTE: There is a wide range of sequence assembly and analysis software; we use Codon code Aligner V11.0.1 and MEGA11 as examples for illustration.

1. Use the aligner to open the sequencing file to observe the peak map and splice the bidirectional sequencing peak map files.
   1. Select Create a New Project and click OK to go to the software homepage. Select Add Samples under the File menu, select the File to be processed as sequence trace format, and import it directly into the software as the imported file at Unassembled Samples.
   2. Choose Contig, then, in the section, choose Advanced Assembly, and then choose Assemble in Groups. First-time users, after selecting Assemble in Groups, will see this dialog: Choose Define names. Modify the Define sample name parts according to the file name and click on Assemble. The aligned sequence will be displayed.
   3. Retrieve the complete ITS2 sequences according to hidden Markov model (HMM)-based model analysis. Use the 5.8s (CGAGGGCACGCCTGCCTGGGCGTCA) and 28s (CGACCCCAGGTCAGGCGGGGCCACC) sequences³¹ to search the aligned sequence from the beginning to remove the 5.8s and 28s regions.
   4. Select Go, then, in the section, select Search Sequence and click OK to get the matching sequence. Choose Contig, choose Delete in the section to remove the 5.8s and 28s regions and the low-quality sequences, and then the software will output the splicing sequences for matching.
2. Use the splicing sequences obtained from the three parallel experimental steps for identification, select BLAST search using the nucleotide database at The National Center for Biotechnology Information (NCBI), and select the Species Identification Tool and the Specific Primer ITS2 corresponding in the Global Pharmacopoeia Genome Database (GPGD) database.
   NOTE: NCBI³² and GPGD³³ include the nucleotide sequences of medicinal plants and animals and they are recognized as authoritative databases. Their major methods for identification are the methods based on sequence comparison³⁴ (BLAST). The percentage identity from the identity results is used to determine whether the species with the greatest similarity is the target genus.
3. Use MEGA11 for sequence alignment, analyze and compare intraspecific and interspecific sequence variation, and construct a phylogenetic neighbor-joining tree.
   1. Download the sequences of three species of Angelica and one species of Peucedanum praeruptorum as an out-group from NCBI. Analyze these sequences together with the result sequence obtained.
   2. Click File, choose Open a File/Session to import the file, and a dialog will appear. Choose Align > DNA, and select Align by ClustalW for sequence comparison. Choose Phylogenetic Analysis in Data, then return to the home page. Click Construct/Test Neighbor-Joining Tree in PHYLOGENY. The software generates a phylogenetic tree.

Results

Sample DNA quality

The range of OD260/OD280 absorbance ratios was 1.80-1.84. The DNA quantity measured by spectrophotometer for each sample was greater than 100 ng/µL (Table 3), indicating a good quality of sample DNA extraction. This indicates that the samples were not contaminated by protein, RNA, or reagents during the extraction process and that they were of good quality and suitable for downstream PCR amplification.

Discussion

Molecular identification technology is easier to learn and master than traditional identification methods. It overcomes the limitations of traditional medicinal herb identification, as it is not affected by the plant's growth stage and does not rely on subjective judgment or the accumulation of specialized expertise³⁵. Other species are confused with A. sinensis due to its similar morphology and medicinal parts, but with the use of DNA barcoding, it can be definitively identified...

Materials

Name	Company	Catalog Number	Comments
CTAB Rapid Plant Genomic DNA Extraction Kit	Shanghai Huiling Biotechnology Co.	NG411M	Suitable for rapid extraction of high quality genomic DNA from different tissues of a wide range of plants
DL5000 DNA Marker	Nanjing vazyme Bio-technology Co.	MD102-02	Ready-to-use product, take an appropriate amount of this product directly for electrophoresis when running the gel.
Electrophoresis	Beijing Liuyi Biotechnology Co.	DYY-6C	Adopt touch screen design, can display set voltage, set current at the same time, parallel output
Ethylenediaminetetraacetic acid	BeijingpsaitongBiotechnologyCo.,Ltd	E70015-100G	Nuclear Isolation Buffer formulation reagents.
Goldview Nucleic Acid Gel Stain(10,000×)	Yisheng Biotechnology (Shanghai) Co., Ltd	10201ES03	When using agarose gel electrophoresis to detect DNA, it binds to DNA and produces a strong fluorescent signal.
High-Speed Tabletop Centrifuge	Changsha High-tech Industrial Development Zone Xiangyi Centrifuge Instrument Co.	H1650	For fast and efficient separation of samples, this compact and lightweight centrifuge offers reliable safety
High-Throughput Tissue Grinder	Shanghai Jingxin Industrial Development Co.	Tiss-48	A high-frequency vibration instrument for grinding samples
http://www.gpgenome.com/	Institute of Chinese Materia Medica, China Academy of Chinese Medical Sciences	-	HerbGenomics database includes genome sequences, gene sets, organelle genomes, low coverage genome data and DNA barcode sequences.
https://www.ncbi.nlm.nih.gov/	U.S. National Library of Medicine	-	The National Center for Biotechnology Information advances science and health by providing access to biomedical and genomic information.
Kylin-Bell	Haimen Qilinbel Instrument Manufacturing Co.	VORTEX-5	Rapid mixing in the form of a high-speed vortex, mixing speed, uniformity, thoroughness
LifeTouch	Hangzhou BORI Technology Co.	TC-96/G/H(b)B	Adoption of advanced thermoelectric refrigeration technology and newly created TAS technology to enhance its overall performance
Multi-functional Gel Image Analysis System	Southwest Operation Center of Shanghai Tianneng Life Science Co.	Tanon-Mini Space 2000	Performs rapid gene amplification experiments with a gradient function for mapping amplification conditions and a gradient temperature range of up to 30°C
NaCl	Beijing Solarbio Science&Technology Co.,Ltd.	S8210-100	Nuclear Isolation Buffer formulation reagents.
Nanodrop One	Genes Ltd.	ND ONE	Quantify DNA, RNA and protein samples in seconds with just 1-2 µL of sample
Polyvinyl pyrrolidone	Shanghai yuanye Bio-Technology Co., Ltd	S30268-500g	Nuclear Isolation Buffer formulation reagents.
Snowflake Ice Maker	Shanghai Zhixin Experimental Instrument Technology Co.	ZX-60X	Adopting rotary extrusion ice making method, fast ice making speed and high efficiency of ice production
Stainless Steel Beads for Tissue Homogenizer	Beyotime Biotechnology.	F6623	Equipment for grinding and mixing of tissue and other samples by vibration
Tris Acetate-EDTA buffer	Beyotime Biotech Inc	ST716	TAE is a commonly used buffer for DNA electrophoresis, frequently employed in agarose gel electrophoresis.
Tris-HCl	Beijing Solarbio Science&Technology Co.,Ltd.	T8230	Nuclear Isolation Buffer formulation reagents.
Water bath Kettle	Shanghai Senxin Experimental Instrument Co.	DK-8D	Precise thermostat and temperature regulation, accurate and reliable temperature control
β-mercaptoethanol	Shanghai Eon Chemical Technology Co.	R054186-100ml	Nuclear Isolation Buffer formulation reagents.