Next Article in Journal
Ecological and Landscape Factors Affecting the Spread of European Mistletoe (Viscum album L.) in Urban Areas (A Case Study of the Kaliningrad City, Russia)
Previous Article in Journal
Mutations of OsPLDa1 Increase Lysophospholipid Content and Enhance Cooking and Eating Quality in Rice
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Taxonomic Notes on the ‘Mahat’ (Artocarpus lacucha and A. thailandicus, Moraceae) Species Complex in Thailand

by
Chaiwat Aneklaphakij
1,
Somnuk Bunsupa
1,
Yotsawate Sirichamorn
2,
Bhanubong Bongcheewin
3 and
Veena Satitpatipan
1,*
1
Department of Pharmacognosy, Faculty of Pharmacy, Mahidol University, Bangkok 10400, Thailand
2
Department of Biology, Faculty of Science, Silpakorn University, Sanam Chandra Palace Campus, Nakhon Pathom 73000, Thailand
3
Department of Pharmaceutical Botany, Faculty of Pharmacy, Mahidol University, Bangkok 10400, Thailand
*
Author to whom correspondence should be addressed.
Plants 2020, 9(3), 391; https://doi.org/10.3390/plants9030391
Submission received: 15 February 2020 / Revised: 13 March 2020 / Accepted: 20 March 2020 / Published: 22 March 2020
(This article belongs to the Section Plant Systematics, Taxonomy, Nomenclature and Classification)

Abstract

:
‘Mahat’ is a well-known medicinal plant utilized in Thailand. The Thai name ‘Mahat’ has been used in many scientific articles for years. However, it is, unpredictably, a homonym of two scientific names in Flora of Thailand, i.e., A. lacucha and A. thailandicus. Additionally, both species are complex due to their high morphological variation. This causes difficulties in species identification especially when this Thai name is referred to as the scientific name for research publication, quality control of pharmaceutical raw materials, and registration of pharmaceutical products. In this study, we scrutinized the taxonomy of ‘Mahat’ by detailed examination of its morphology and distribution, including molecular and qualitative phytochemical studies. Leaf surfaces were inspected using scanning electron microscopy. The phylogeny of both species was studied using DNA sequences of nuclear and plastid regions. Chromatographic fingerprints, focusing on the major active compound oxyresveratrol, were identified using high-performance liquid chromatography. According to our current study, phylogenetic evidence showed that some samples of both species were clustered together in the same clade and phytochemical fingerprints were almost identical. These results are valuable data for taxonomic revision in the near future and reveal the possible utilization of A. thailandicus as a new material source of oxyresveratrol in the pharmaceutical industry.

1. Introduction

Moraceae is one of the most well-known plant families of the world, comprising about forty genera and over a thousand species [1]. The family is mainly distributed in the tropics [2]. Molecular and morphological evidences [3,4,5] indicate that the family is monophyletic; however, the classification at the tribal or infra-tribal level is still obscure due to the extraordinary diversity of complex inflorescence structures, pollination syndromes, and breeding systems [2]. The genus Artocarpus J.R. Forst. & G. Forst. is the third-largest genus of the family and the largest genus in the tribe Artocarpeae [2]. The center of distribution of the genus is in Borneo [6]. In 1847, Trécul divided the genus into two subgenera, Jaca Trécul and Pseudojaca Trécul [7]. The former subgenus has spirally arranged leaves, amplexicaul stipules, and annulate stipule scars, whereas the latter has alternate-distichous leaves, nonamplexicaul stipules, and lateral stipule scars. In 1959, Jarrett renamed the subgenus Jaca to Artocarpus [8,9]. Thereafter, Zerega and colleagues reduced the genus Prainea King to Artocarpus subgen. Prainea (King) Zerega, and the series Cauliflori Jarrett raised to Artocarpus subgen. Cauliflori (Jarrett) Zerega [2]. Thus, Artocarpus includes four subgenera, Artocarpus, Pseudojaca, Prainea, and Cauliflori [2].
In Thailand, the genus Artocarpus contains fourteen species with one endemic species, A. thailandicus C.C. Berg [10,11]. Several species are used for medicinal purposes, including A. heterophyllus Lam., A. lacucha Roxb. ex Buch.-Ham., A. thailandicus, and A. altilis (Parkinson) Fosberg. In Thai traditional medicine, the plant with the vernacular name ‘Mahat’ (common name: monkey jack) is widely used as an anthelmintic drug [12]. Furthermore, ‘Mahat’ is a major source of oxyresveratrol (trans-2,3′,4,5′-tetrahydroxystilbene), which is an active ingredient in skin whitening products, with antityrosinase, antioxidant, and anti-inflammatory activity [13,14,15,16,17]. Due to the morphological similarities, A. lacucha and A. thailandicus share the vernacular name ‘Mahat’. Although the Moraceae has already been revised in Flora of Thailand [11] with a great deal of taxonomic information, the key to species of the genus Artocarpus is based mainly on leaf and inflorescence characteristics which remain troublesome for the identification of ‘Mahat’. For example, presence of indumentum in the areoles are not constant characteristics in all examined specimens, or even in all leaves of the same specimen. This variation was also mentioned in Flora of Thailand [11]. Furthermore, the use of inflorescence characters can be problematic in non-flowering collections. This is affecting quality control of the raw materials in the pharmaceutical industry and the registration of pharmaceutical products with the Thai Food and Drug Administration (Thai FDA).
In this study, we attempted to verify the taxonomy of A. lacucha and A. thailandicus, which are used as the crude drug named ‘Mahat’ in Thai traditional medicine and herbal pharmaceutical products, using macroscopic and microscopic morphology, molecular phylogeny, and qualitative phytochemistry. We believed that our results provide valuable data for the taxonomic revision of ‘Mahat’ in the future as well as being helpful in the pharmaceutical product industry, as A. thailandicus is possibly a new source of plant raw material for oxyresveratrol consumption.

2. Results

2.1. Macroscopic Morphology and Distribution

The macroscopic morphology of collected samples and herbarium specimens were examined. The morphological characters of both species conformed to the plant description in Flora of Thailand [11]. The main characters differentiating the species are summarized in Figure 1 and Table 1. The distribution of both species in Thailand is presented in Figure 2. Furthermore, the typification details of A. lacucha were not mentioned in Flora of Thailand. However, during our herbarium observations we found the type specimen as shown below.
Type:—BANGLADESH. Rangpur: Chilmari, 16 December 1808, Collector Unknown 4655.B [holotype K (http://specimens.kew.org/herbarium/K001039600, male, and http://specimens.kew.org/herbarium/K001039601, female, last acc. 1 September 2018)].
In addition, there was the confusion of the nomenclature, i.e., A. lacucha and A. lakoocha Roxb., and its typification. Mabberley confirmed that A. lacucha was the correct scientific name and the type specimen was collected in K-W with Hamilton’s label [18]. The species name A. lakoocha, however, was mentioned in the samples’ DNA sequences obtained from GenBank database and it still was used in phylogenetic analyses in order to show that they are GenBank sequences. Moreover, in the plant description of A. lacucha, A. dadah Miq. was presented as one of the synonyms. Williams et al., in their phylogenetic and biogeographic investigation of Artocarpus showed that A. dadah showed monophyly with A. lacucha [6]. Hence, A. dadah should be discarded from the description in Flora of Thailand [11].

2.2. Leaf Surface Morphology

Fourteen samples were selected, as representatives of the species and covering all floristic regions in Thailand, for scanning electron microscopy (SEM) of leaf surface characteristics. The upper and lower leaf surfaces, including the indumentum of individual species, are shown in Figure 3 and Figure 4, respectively. The upper leaf surfaces of both species were similar. They could be glabrous (Figure 3A,C) or have some straight and uncinate smooth hairs on the veins and lamina (Figure 3B,D). The lower leaf surfaces of specimens identified as A. lacucha had rather varied indumentum, from gnarled papillae to dense long straight-gnarled hairs on the veinlets, with smooth straight and/or uncinate hairs located sparsely on the lateral veins, but with almost glabrous areoles (Figure 4A–D). Similar to A. lacucha, the lower leaf surface of specimens identified as A. thailandicus had long straight-gnarled hairs on the veinlets. However, the areoles were densely covered with tomentose hairs and dense, smooth, straight, or uncinate hairs on the lateral veins (Figure 4E,F). In addition, lower surface of some A. thailandicus specimens (PBM05203, PBM05204, and PBM05205) showed gnarled papillae on the veinlets and lack of the indumentum in the areoles (Figure 4G,H).

2.3. Phylogenetics

Fifty-seven taxa were included in the phylogenetic analyses. All taxa provided ITS and trnH-psbA data, and thirty-five and fifty-six taxa provided ETS and trnL-F sequences, respectively. Maximum parsimony (MP) and maximum likelihood (ML) analyses of the separate markers, nuclear (ITS and ETS) and chloroplast (trnL-F and trnH-psbA) loci (Figures S1–S24) did not show any major conflicts in topology and thus all were combined (Figures S25–S28). The combined four aligned regions were 1726–1829 bp long. MP analysis of the four combined regions showed 741 potentially informative characters and resulted in three most parsimonious trees (MPTs) with a length = 282, the consistency index = 0.707, the retention index = 0.903 and the composite index = 0.779 (0.639) for all sites and parsimony-informative sites (in parentheses). For ML analysis, the discrete gamma distribution (5 categories (+G, parameter)) was 0.454. MP analysis conformed to ML analysis. The MP and ML consensus tree (Figure 5) revealed a major clade containing specimens of A. lacucha, A. thailandicus and some other Artocarpus species included in the analyses (sequences obtained from GenBank) with several minor clades. Most of them have low supports and some parts of the cladogram had an unresolved topology. Two minor clades, indicated in the blue and red square dashed line, showed that A. lacucha and A. thailandicus clustered together, with low bootstrap support. Both species were likely to be polyphyletic as they presented poor resolution of relationships. Only some samples of A. lacucha formed only one well-supported clade with 99% bootstrap support (indicated by green square dashed line). Likewise, A. thailandicus seems to be presented in more than two clades, all with low support values. Additionally, the leaf indumentum from SEM results were mapped and showed some correlation with cladogram, as shown in Figure 5.

2.4. Qualitative Phytochemistry

In this study, we focused mainly on the qualitative analysis of the major active compound, oxyresveratrol, in the extract of both species by high-performance liquid chromatography (HPLC). Based on the HPLC chromatograms produced by both species were almost identical (Figure 6). Oxyresveratrol was detected at the retention time of 17.3 min in the extract of both species.

3. Discussion

There have been many publications reporting the use of ‘Mahat’ from Thailand as a material source in various fields such as biogeography, phylogeny, phytochemistry, and biological activity testing. Nevertheless, A. lacucha [6,12,15,16,19,20,21,22,23] has been cited more frequently than A. thailandicus [6,24]. This is likely because the name A. lacucha is better-known and has been used for longer time (since 1814) [25] than A. thailandicus (since 2005) [10]. The confusion over the naming of these plants has been exacerbated by the high variation in their vegetative characters. Furthermore, the heartwood, which has no distinguishing characteristics, has also been used for research and pharmaceutical production. It is thus more difficult to distinguish the species origin of ‘Mahat’ raw materials, and consequently their adulteration. Therefore, further taxonomic investigation of these two species in Thailand is needed. If there are other evidences that reject the monophyly of A. thailandicus, its taxonomic status at species level should also be discarded. According to our studies, however, samples of A. lacucha were collected only in Thailand, while this species is distributed from South, East and Southeast Asia to Oceania [11]. More samples of A. lacucha from the areas outside Thailand should be included into the analyses to clarify the complexity of both species.
The species circumscriptions of Artocarpus subgen. Pseudojaca is one of the most difficult [6]. The alteration of taxonomic rank at species level could occur due to variations in leaf and inflorescence morphology [6]. Although the morphology of A. lacucha and A. thailandicus are different and appear sufficient to be used for efficient identification, there are some limitations to those characters. According to the taxonomic key of the genus Artocarpus in Flora of Thailand [11] and the results from our study, the presence of hairs in the areoles may be used for the identification of A. thailandicus. This character was dominant, particularly when examined under a stereo microscope, but the variation in hair density was observed in some specimens. Nevertheless, in some examined specimens, this character was absent (Pooma et al. 5749, Aneklaphakij 10, Aneklaphakij 11, Aneklaphakij 13). The specimen “Wongprasert s.n.” showed a presence of hairs only in some leaves. Additionally, this variation was also mentioned in Flora of Thailand [11]. Moreover, the inflorescence is another significant part used to differentiate each species and can be problematic in non-flowering specimens. Furthermore, we hypothesized that Berg et al. [11] examined mostly the ancient herbaria specimens using only stereo microscope with less magnification than SEM. Thus, some characteristics could not be observed; for example, the gnarled papillae on the veinlets could be seen only under SEM. Therefore, our leaf morphology result was probably different from the previous work [11] and could imply that the species identification using leaf indumentum might not be suitable criterion, although, the characteristic is not affected by environmental factors (the first author’s personal observation). The second character was indumentum color, which is white and brown, the interpretation of which may be subjective. The third character was the length of the peduncle of staminate and pistillate inflorescences, with A. lacucha having longer peduncles than A. thailandicus. However, this character could vary during different stages of maturity. Moreover, it is suggested that hybridization between these two species is likely to occur [11]. Evidences of interspecific hybridization are still unclear and hybrid determination is required. However, this biological phenomenon can be common in the genus Artocarpus, as found among species of commercial breadfruits, A. altilis (Park.) Fosberg, A. mariannensis Trécul and A. camansi Blanco, whose hybrids shared morphological and cytological characters [26]. Cross-pollination by wind with the small, powdery pollen grains spread were also reported for some species of breadfruit trees [27,28]. Therefore, it can be hypothesized that interspecific hybridization between A. lacucha and A. thailandicus can also be possible. The MP & ML strict consensus trees showed that some of A. lacucha and A. thailandicus occurred in the same clade, with low bootstrap support, and tree topology was unresolved. This phenomenon could occur when studying the population scale as found in in the genus Secale L. (Poaceae) by Maraci and colleagues [29]. It may be caused by the high similarity in the sequences of both species and the other Artocarpus spp. included in the study or the speciation of them might be based on other evidences rather than molecular data such as geographical origin [29]. Even though both species evolved from the same common ancestor, as in the Williams and colleagues’ study [6], the incongruence between morphological characters and molecular evidence, which is a common phenomenon in Moraceae such as Ficus krishnae C. DC. [30], could result from asymmetric evolution rates in DNA sequences and phenotypic morphogenesis. In addition, from this phenomenon, the connected morphological variation and molecular data probably could not be observed. Nevertheless, there was a well-defined clade (indicated by green square dashed line, Figure 5) of A. lacucha. We discovered that the SEM results of samples PBM05199, PBM05206, and PBM05214 shared the synapomorphic characters which the upper leaf surface was glabrous and lower leaf surface was gnarled papillae on the veinlets, sparsely smooth straight or uncinate hairs on lateral veins. Therefore, we attempted to correlate the micromorphology of leaf surface with molecular data as shown in Figure 5. Only a specimen in the early diverging clade (encircled by red square dashed line) was SEM analyzed and it shows straight and uncinate smooth hair on upper surface and gnarled long straight hairs mostly on veinlets of leaf lower surface. Those characteristics of leaf indumentum were hypothesized to be also present among other samples in this clade and also conformed to the description in Flora of Thailand [11], even though Berg et al. did not use SEM for leaf morphology examination. However, those features have evolved more than once in phylogeny of ‘Mahat’, as shown in Figure 5. For the later diverging clades (encircled by green and blue square dashed line), they could be noticed that later evolved three traits: 1. glabrous upper surface (shown by light blue square), 2. gnarled papillae on the veinlets with sparsely smooth straight or uncinate hairs on lateral veins (shown by yellow star), and 3. gnarled long straight hairs mostly on the veinlets with tomentose hair in the areoles and densely smooth straight or uncinate hairs on the lateral veins on lower surface (shown by grey star). Leaf indumentum is beneficial for plant protection against herbivory, pathogens or strong sunlight and also relates to drought tolerance by facilitating condensation of air moisture onto the plant surface [31,32,33]. The upper surface of A. lacucha, however, seemingly become more glabrous and the reasons of this evolutionary trend are still unclear. In order to compensate the absence of leaf indumentum on the upper surface, A. lacucha might develop a thicker cuticle. More detailed studies of leaf anatomy, especially the cuticle thickness, are required to understand trend of evolution in A. lacucha’s leaf indumentum. Moreover, leaf anatomical characteristics apart from indumentum are potentially useful for taxonomic purposes in Moraceae. For example, in Ficus subsection Urostigma, the epidermis and lithocysts characteristics could be used for identification of Ficus subsection Conosycea [34]. Since the relationship between these two species are still not well-supported to ensure species circumscription, this result needs to be confirmed. Sampling from a wider range of distribution, including areas outside Thailand and the use of additional DNA markers, might help to resolve the taxonomy of ‘Mahat’. Most of the morphological differences observed are relevant to the expression of important regulatory genes [35,36]; however, the studied DNA regions might be housekeeping regions that are not involve in the expression of the relevant characters or not involved in their regulation [30]. Thus, genes related to the expression of morphological characters that differ between these two species should be investigated in future.
Chromatography is one tool for classifying and identifying plants, especially cases where plants produce unique chromatogram patterns for specific chemical markers [37,38]. We investigated the use of high-performance liquid chromatography to differentiate the phytochemical fingerprint of the major active compound, oxyresveratrol, in the extract of both species. Oxyresveratrol is commonly collected from the heartwood; however, we tested twigs in this study due to the limitations of sample collection. Oxyresveratrol was detected in both species’ extracts; however, there is likely variation caused by differences in season, age and geography of the samples. The aligned chromatograms showed high similarity among species. Thus, phytochemical analysis conformed to the molecular studies, that the species could not be differentiated. Moreover, these species can be used interchangeably in traditional medicine and pharmaceutical products as the active compound is identical. However, future studies should investigate metabolites using high-performance liquid chromatography coupled with high-resolution mass spectrometry.
According to our study, both molecular technique and phytochemical study focusing on major active compound could be useful tools for discrimination of plant species, especially in the quality control of raw materials for pharmaceutical product production. Generally, the raw material is in a powder form which is very difficult for plant species identification. Therefore, the molecular technique and qualitative phytochemical analysis can be used for raw material quality examination prior using for production. These two techniques are very suitable because a small amount of sample is used, the process consumes a moderate amount of time and could reduce production costs in cases where undesired raw material is found before mass production.

4. Materials and Methods

4.1. Taxonomic Study

The plant samples were collected from the fieldwork between October 2015 and November 2017. For representativeness, plants from all Thai floristic regions were sampled [11]. All the samples were first identified using the identification key presented in Flora of Thailand. Each sample was examined under the stereo and light microscopes supplemented with an ocular micrometer. All voucher specimens were deposited at the Herbarium of Department of Pharmaceutical Botany (PBM), Mahidol University, Thailand. There were 35 collections in this study. In total, 18 samples were identified to A. lacucha, collected from Northern, North-eastern, Central and Peninsular regions of Thailand, whereas 17 samples identified as A. thailandicus were collected from all floristic regions except Peninsular region. The list of studied samples is shown in Table 2. Moreover, 72 and 23 specimens of A. lacucha and A. thailandicus, respectively, from BK, BKF, E, K, L, and PBM (Thiers 2019 [continuously updated]) were also studied. All specimens were seen by the first author. The GPS coordinates of the collections were retrieved from the herbarium labels, the fuzzy gazetteer (http://isodp.hof-university.de/fuzzyg/query) and Google Maps. The map of distribution in Thailand based on the localities of all collections was constructed using SimpleMappr, http://www.simplemappr.net/ [39].

4.2. Scanning Electron Microscopy (SEM)

Leaves were collected at the third to the fourth node from the branch apex, with at least three leaves for each accession. Samples were cleaned and dried in a hot air oven at 60 °C for a day, then preserved in silica gel for the subsequent scanning electron microscope (SEM) study. The upper and lower leaf surfaces of seven samples of A. lacucha (voucher PBM05199, PBM05206, PBM05214, PBM05218, PBM05221, PBM05225, PBM05231) and A. thailandicus (voucher PBM05202, PBM05203, PBM05204, PBM05205, PBM05210, PBM05219, PBM05229) were selected as representatives of the species and all Thai floristic regions. Dried leaves were cut to the size of 5 × 5 mm from the median area between the margin and midrib, and mounted directly on stubs using double-sided carbon tape, then coated with gold by a Polaron Range SC7620 Sputter Coater & CA7625 Carbon Accessary and platinum palladium by Hitachi E102 sputter coater, Japan. The samples were photographed by FEI Quanta 450, Hillsboro, OR, USA (accelerating voltage at 10.0 kV) at the Scientific Equipment Center, Faculty of Science, Kasetsart University and by Hitachi SU8010, Japan (accelerating voltage at 10.0 kV) at the Center of Nanoimaging, Faculty of Science, Mahidol University, Bangkok, Thailand.

4.3. DNA Extraction, Amplification and Sequencing

Young leaves of collected samples were dried and stored in silica gel for DNA extraction. They were milled into the fine powder using a mortar and pestle. The DNA extraction process followed the protocol of the plant genomic DNA extraction kit (RBC Bioscience, Taiwan). Purified DNA was kept at −20 °C before being used as a template in a polymerase chain reaction (PCR). In this study, specific primer sequences in Table 3 (synthesized by Macrogen Inc., Seoul, South Korea) were employed for amplifying the nuclear genome (ITS and ETS region) and chloroplast genome (trnL-F and trnH-psbA intergenic spacer). PCR buffer, dNTPs and Taq polymerase reagents were purchased from RBC Bioscience, Taiwan. DNase/RNase free distilled water was bought from Invitrogen, USA. The PCR condition of each gene was based on previous studies (See in the reference of each primer sequences). After the PCR reaction, the product was checked using gel electrophoresis. The PCR products were purified using a HiYield™ Gel/PCR DNA Fragments Extraction Kit, then the concentration was determined using a SmartSpecTM Plus Spectrophotometer (Bio-rad, Philadelphia, PA, USA). The samples containing a concentration of more than 50 ng/µL were delivered for sequencing. For the sequencing procedure, the capillary electrophoresis sequencing method (standard sequencing) was performed by Macrogen sequencing services (Macrogen Inc., Seoul, South Korea). Apart from our samples, we retrieved four samples of A. lacucha and three samples of A. thailandicus from the NCBI GenBank database to be included in the phylogenetic analysis. However, these samples contain only ITS region, trnL-F and trnH-psbA intergenic spacer.

4.4. Sequence Alignment and Phylogenetic Analysis

Subsequent to the sequencing process being accomplished, each sequence was aligned by ClustalW [44] using the biological sequence alignment editor (BioEdit). The positions containing more than 75% missing or undetermined characters were eliminated. The analyses of each separate locus, nuclear and chloroplast loci, as well as four combined molecular DNA markers were analyzed using the maximum parsimony (MP) and maximum likelihood (ML) method for phylogenetic tree reconstruction by molecular evolutionary genetics analysis version 7.0 (MEGA7) [45]. The MP tree was obtained using the Tree-Bisection-Regrafting (TBR) algorithm [46] (p. 126) with search level 1 in which the initial trees were obtained by the random addition of sequences (10 replicates). Model selection for ML analysis was conducted using the model selection tool as implemented in MEGA7. From the model test result, the ML tree was conducted based on the Tamura 3-parameter model [47] for four combined regions; chloroplast loci; trnL-F and trnH-psbA intergenic spacer, while Kimura 2-parameter model [48] was used for the nuclear loci, ITS, and ETS regions. The initial tree for the heuristic search was obtained automatically by applying Neighbor-Join and BioNJ algorithms to a matrix of pairwise distances estimated using the Maximum Composite Likelihood (MCL) approach and then selecting the topology with a superior log-likelihood value. For four combined regions and nuclear loci, a discrete gamma distribution was used to model evolutionary rate differences among the sites (5 categories (+G, parameter)). In both methods, the bootstrap strict consensus tree inferred from 1000 replicates was taken to represent the evolutionary history of the taxa analyzed and the percentage of replicate trees in which the associated taxa clustered together in the bootstrap test (1000 replicates) were shown next to the branches [49]. Branches corresponding to partitions reproduced in less than 50% bootstrap replicates were collapsed. The sequences of each region from each sample were submitted into the NCBI GenBank database. All of the accession numbers were shown in Table 4. The sequences of A. lacucha, A. thailandicus, A. dadah, A. heterophyllus Lam., A. nitidus subsp. griffithii (King) F.M. Jarrett, A. nitidus subsp. lingnanensis (Merr.) F.M. Jarrett, A. fulvicortex F.M. Jarrett, A. integer (Thunb.) Merr., and Ficus racemosa L. in the NCBI GenBank database were retrieved and also included in the phylogenetic analysis. In addition, A. integer and A. heterophyllus were selected as outgroup. Bootstrap percentages were described by Chantarasuwan and colleagues [40] as strong (85–100%), moderate (75–84%), low (50–74%), and no support (≤50%).

4.5. Qualitative Phytochemical Analysis

The twigs of all collected samples were chopped, dried in a hot air oven at 60 °C for 1–2 days and ground using a hammer miller into the moderated powder (through sieve No.14), then kept in a dry place. The high-performance liquid chromatographic (HPLC) method was used for qualitative phytochemical analysis. The twig powder of each sample was macerated with methanol AR grade (Fisher Scientific, Pittsburgh, PA, USA) in the ratio of 1:30 at room temperature for 72 h on the shaker (100 rpm). The extract was filtered by Whatman No.1 filter paper, followed by evaporating to dryness under reduced pressure at 40 °C using a rotary evaporator (Buchi, Flawil, Switzerland). All of the extracts were diluted with methanol HPLC grade (Honeywell Burdick & Jackson™, Muskegon, MI, USA) to the concentration of 1 mg/mL before injection. The chromatographic analysis was performed on a Shimadzu, LC-20AD liquid chromatography system containing a quaternary solvent delivery system, autosampler (SIL-20AHT), degasser (DGU-20A5R), column oven (CTO-20A) and photodiode array detector (PAD, SPD-M20A). Each chromatographic data acquisition was performed using Lab Solutions software. The completely validated chromatographic condition was adapted from the previous study [50]. Chromatographic separation was achieved on a C18 reverse phase column (Phenomenex, Torrance, CA, USA, Gemini-NX C18, 250 mm × 4.6 mm i.d., 5 µm) using a 1 mL/min mobile phase gradient system. The mobile phase comprises acetonitrile (Sigma, St. Louis, MO, USA) and 0.5% aqueous acetic acid (v/v). The gradient elution profile was 18%–25% acetonitrile at 0–10 min, and 25%–40% acetonitrile at 10–25 min, which was then returned to equilibration at the initial conditions for 10 minutes before the injection of the next sample. The column temperature was maintained at 30 °C and the injection volume was 20 µL. The detection was examined at 320 nm. The oxyresveratrol standard was prepared at a concentration of 6.25 µg/mL and injected with the same amount and under the same conditions as the sample. Carbamazepine (Himedia, Mumbai, India) in the concentration of 100 µg/mL was augmented as the internal standard.

5. Conclusions

‘Mahat’ is a well-known medicinal plant in Thailand. The ambiguity of the taxonomic classification of the source of ‘Mahat’ is pivotal for quality control of its pharmaceutical product yields. We revised the plant morphology and distribution and investigated the molecular evidence and phytochemical fingerprints. Morphological characters were highly variable especially the indumentum in the areoles, and phylogenetic evidence demonstrated that A. lacucha is likely to be polyphyletic, with only one well-supported clade. In addition, A. thailandicus seems to form more than one clade. Moreover, qualitative phytochemical analysis by HPLC showed that the phytochemical fingerprints of the major active compound, oxyresveratrol, were almost identical. Hence, A. thailandicus might be utilized in pharmaceutical industry as another new plant material source of oxyresveratrol in the future. According to our present results, we still keep both A. lacucha and A. thailandicus as separate species even if the molecular and phytochemical evidences show these two species could be the same. In order to reduce A. thailandicus as a synonym of A. lacucha, further deeply detailed taxonomic researches of both species should be re-investigated.

Supplementary Materials

The following are available online at https://www.mdpi.com/2223-7747/9/3/391/s1, Figure S1: ITS region maximum parsimony consensus tree of A. lacucha and A. thailandicus with other Artocarpus species. A. integer and A. heterophyllus were chosen as an outgroup, Figure S2: ETS region maximum parsimony consensus tree of A. lacucha and A. thailandicus, Figure S3: trnL-F intergenic spacer maximum parsimony consensus tree of A. lacucha and A. thailandicus with other Artocarpus species. A. integer and A. heterophyllus were chosen as an outgroup, Figure S4: trnH-psbA intergenic spacer maximum parsimony consensus tree of A. lacucha and A. thailandicus with other Artocarpus species. F. racemosa was chosen as an outgroup, Figure S5: Nuclear loci maximum parsimony consensus tree of A. lacucha and A. thailandicus with other Artocarpus species. A. integer and A. heterophyllus were chosen as an outgroup, Figure S6: Chloroplast loci maximum parsimony consensus tree of A. lacucha and A. thailandicus with other Artocarpus species. F. racemosa was chosen as an outgroup, Figure S7: ITS region maximum likelihood consensus tree of A. lacucha and A. thailandicus with other Artocarpus species. A. integer and A. heterophyllus were chosen as an outgroup, Figure S8: ETS region maximum likelihood consensus tree of A. lacucha and A. thailandicus, Figure S9: trnL-F intergenic spacer maximum likelihood consensus tree of A. lacucha and A. thailandicus with other Artocarpus species. A. integer and A. heterophyllus were chosen as an outgroup, Figure S10: trnH-psbA intergenic spacer maximum likelihood consensus tree of A. lacucha and A. thailandicus with other Artocarpus species. F. racemosa was chosen as an outgroup, Figure S11: Nuclear loci maximum likelihood consensus tree of A. lacucha and A. thailandicus with other Artocarpus species. A. integer and A. heterophyllus were chosen as an outgroup, Figure S12: Chloroplast loci maximum likelihood consensus tree of A. lacucha and A. thailandicus with other Artocarpus species. F. racemosa was chosen as an outgroup, Figure S13: ITS region maximum parsimony (without branches collapsed) tree of A. lacucha and A. thailandicus with other Artocarpus species. A. integer and A. heterophyllus were chosen as an outgroup, Figure S14: ETS region maximum parsimony (without branches collapsed) tree of A. lacucha and A. thailandicus, Figure S15: trnL-F intergenic spacer maximum parsimony (without branches collapsed) tree of A. lacucha and A. thailandicus with other Artocarpus species. A. integer and A. heterophyllus were chosen as an outgroup, Figure S16: trnH-psbA intergenic spacer maximum parsimony (without branches collapsed) tree of A. lacucha and A. thailandicus with other Artocarpus species. F. racemosa was chosen as an outgroup, Figure S17: Nuclear loci maximum parsimony (without branches collapsed) tree of A. lacucha and A. thailandicus with other Artocarpus species. A. integer and A. heterophyllus were chosen as an outgroup, Figure S18: Chloroplast loci maximum parsimony (without branches collapsed) tree of A. lacucha and A. thailandicus with other Artocarpus species. F. racemosa was chosen as an outgroup, Figure S19: ITS region maximum likelihood (without branches collapsed) tree of A. lacucha and A. thailandicus with other Artocarpus species. A. integer and A. heterophyllus were chosen as an outgroup, Figure S20: ETS region maximum likelihood (without branches collapsed) tree of A. lacucha and A. thailandicus, Figure S21: trnL-F intergenic spacer maximum likelihood (without branches collapsed) tree of A. lacucha and A. thailandicus with other Artocarpus species. A. integer and A. heterophyllus were chosen as an outgroup, Figure S22: trnH-psbA intergenic spacer maximum likelihood (without branches collapsed) tree of A. lacucha and A. thailandicus with other Artocarpus species. F. racemosa was chosen as an outgroup, Figure S23: Nuclear loci maximum likelihood (without branches collapsed) tree of A. lacucha and A. thailandicus with other Artocarpus species. A. integer and A. heterophyllus were chosen as an outgroup, Figure S24: Chloroplast loci maximum likelihood (without branches collapsed) tree of A. lacucha and A. thailandicus with other Artocarpus species. F. racemosa was chosen as an outgroup, Figure S25: Four combined regions maximum parsimony consensus tree of A. lacucha and A. thailandicus with other Artocarpus species. A. integer and A. heterophyllus were chosen as an outgroup. Two minor clades, the blue and red square dashed line, showed that both species were clustered together in the same clade with low bootstrap support. Some parts of the tree showed unresolved topology, Figure S26: Four combined regions maximum likelihood consensus tree of A. lacucha and A. thailandicus with other Artocarpus species. A. integer and A. heterophyllus were chosen as an outgroup. Two minor clades, the blue and red square dashed line, showed that both species were clustered together in the same clade with low bootstrap support. Some parts of the tree showed unresolved topology, Figure S27: Four combined regions maximum parsimony (without branches collapsed) tree of A. lacucha and A. thailandicus with other Artocarpus species. A. integer and A. heterophyllus were chosen as an outgroup, Figure S28: Four combined regions maximum likelihood (without branches collapsed) tree of A. lacucha and A. thailandicus with other Artocarpus species. A. integer and A. heterophyllus were chosen as an outgroup, Supplementary file S1: Four loci multiple alignment.

Author Contributions

Conceptualization, C.A., S.B., Y.S., B.B. and V.S.; Data curation, C.A.; Formal analysis, C.A.; Investigation, C.A., S.B., B.B. and V.S.; Methodology, C.A., S.B., B.B. and V.S.; Project administration, V.S.; Resources, V.S.; Supervision, S.B., Y.S., B.B. and V.S.; Validation, S.B., Y.S., B.B. and V.S.; Visualization, C.A.; Writing—original draft, C.A.; Writing—review & editing, Y.S., B.B. and V.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Mahidol Medical Scholars Program (MSP), Mahidol University, Thailand and Young Researcher Development Program 2019 from National Research Council of Thailand.

Acknowledgments

The authors are very grateful to thank Assoc. Prof. Suppachai Tiyaworanant from Faculty of Pharmaceutical Sciences, Khon Kaen University, Thailand and Asst. Prof. Charan Leeratiwong from Faculty of Science, Prince of Songkla University, Thailand for their help in plant sample collection. The authors wish to also thank Department of Pharmacognosy and Department of Pharmaceutical Botany, Faculty of Pharmacy, Mahidol University for the laboratory devices, facilities as well as the help of the scientist.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Zhang, Q.; Onstein, R.E.; Little, S.A.; Sauquet, H. Estimating divergence times and ancestral breeding systems in Ficus and Moraceae. Ann. Bot. 2018, 123, 191–204. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Zerega, N.J.C.; Nur Supardi, M.N.; Motley, T.J. Phylogeny and Recircumscription of Artocarpeae (Moraceae) with a Focus on Artocarpus. Syst. Bot. 2010, 35, 766–782. [Google Scholar] [CrossRef] [Green Version]
  3. Datwyler, S.L.; Weiblen, G.D. On the origin of the fig: Phylogenetic relationships of Moraceae from ndhF sequences. Am. J. Bot. 2004, 91, 767–777. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Zerega, N.J.C.; Clement, W.L.; Datwyler, S.L.; Weiblen, G.D. Biogeography and divergence times in the mulberry family (Moraceae). Mol. Phylogenet. Evol. 2005, 37, 402–416. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Clement, W.L.; Weiblen, G.D. Morphological evolution in the mulberry family (Moraceae). Syst. Bot. 2009, 34, 530–552. [Google Scholar] [CrossRef]
  6. Williams, E.W.; Gardner, E.M.; Harris, R., III; Chaveerach, A.; Pereira, J.T.; Zerega, N.J.C. Out of Borneo: Biogeography, phylogeny and divergence date estimates of Artocarpus (Moraceae). Ann. Bot. 2017, 119, 611–627. [Google Scholar] [CrossRef] [Green Version]
  7. Trécul, A. Mémoire sur la famille des Artocarpées. Ann. Sci. Nat. 1847, 3, 6–383. [Google Scholar]
  8. Jarrett, F.M. Studies in Artocarpus and allied genera, I. General considerations. J. Arnold Arbor. 1959, 40, 1–29. [Google Scholar]
  9. Jarrett, F.M. Studies in Artocarpus and allied genera, II. A revision of Prainea. J. Arnold Arbor. 1959, 40, 30–37. [Google Scholar]
  10. Berg, C.C. A New Species of Artocarpus (Moraceae) from Thailand. Blumea 2005, 50, 531–533. [Google Scholar] [CrossRef] [Green Version]
  11. Berg, C.C.; Pattharahirantricin, N.; Chantarasuwan, B. Moraceae. In Flora of Thailand; Santisuk, T., Larsen, K., Eds.; The Forest Herbarium, Department of National Parks, Wildlife and Plant Conservation: Bangkok, Thailand, 2011; Volume 10, pp. 475–675. [Google Scholar]
  12. Maneechai, S.; Likhitwitayawuid, K.; Sritularak, B.; Palanuvej, C.; Ruangrungsi, N.; Sirisa-ard, P. Quantitative Analysis of Oxyresveratrol Content in Artocarpus lakoocha and ‘Puag-Haad’. Med. Princ. Pract. 2009, 18, 223–227. [Google Scholar] [CrossRef] [PubMed]
  13. Lorenz, P.; Roychowdhury, S.; Engelmann, M.; Wolf, G.; Horn, T.F.W. Oxyresveratrol and resveratrol are potent antioxidants and free radical scavengers: Effect on nitrosative and oxidative stress derived from microglial cells. Nitric Oxide 2003, 9, 64–76. [Google Scholar] [CrossRef] [PubMed]
  14. Chung, K.O.; Kim, B.Y.; Lee, M.H.; Kim, Y.R.; Chung, H.Y.; Park, J.H.; Moon, J.O. In-vitro and in-vivo anti-inflammatory effect of oxyresveratrol from Morus alba L. J. Pharm. Pharmacol. 2003, 55, 1695–1700. [Google Scholar] [CrossRef] [PubMed]
  15. Likhitwitayawuid, K.; Sornsute, A.; Sritularak, B.; Ploypradith, P. Chemical transformations of oxyresveratrol (trans-2,4,3′,5′-tetrahydroxystilbene) into a potent tyrosinase inhibitor and a strong cytotoxic agent. Bioorg. Med. Chem. Lett. 2006, 16, 5650–5653. [Google Scholar] [CrossRef]
  16. Tengamnuay, P.; Pengrungruangwong, K.; Pheansri, I.; Likhitwitayawuid, K. Artocarpus lakoocha heartwood extract as a novel cosmetic ingredient: Evaluation of the in vitro anti-tyrosinase and in vivo skin whitening activities. Int. J. Cosmet. Sci. 2006, 28, 269–276. [Google Scholar] [CrossRef]
  17. Likhitwitayawuid, K. Stilbenes with tyrosinase inhibitory activity. Curr. Sci. 2008, 94, 44–52. [Google Scholar]
  18. Mabberley, D.J. Francis Hamilton’s Commentaries with Particular Reference to Meliaceae. Taxon 1977, 26, 529. [Google Scholar] [CrossRef]
  19. Charoenlarp, P.; Radomyos, P.; Harinasuta, T. Treatment of taeniasis with Puag-Haad: A crude extract of Artocarpus lakoocha wood. Southeast Asian J. Trop. Med. Public Health 1981, 12, 568–570. [Google Scholar]
  20. Puntumchai, A.; Kittakoop, P.; Rajviroongit, S.; Vimuttipong, S.; Likhitwitayawuid, K.; Thebtaranonth, Y. Lakoochins A and B, New Antimycobacterial Stilbene Derivatives from Artocarpus lakoocha. J. Nat. Prod. 2004, 67, 485–486. [Google Scholar] [CrossRef]
  21. Likhitwitayawuid, K.; Sritulalak, B.; Benchanak, K.; Lipipun, V.; Mathew, J.; Schinazi, R.F. Phenolics with antiviral activity from Millettia erythrocalyx and Artocarpus lakoocha. Nat. Prod. Res. 2005, 19, 177–182. [Google Scholar] [CrossRef]
  22. Chuanasa, T.; Phromjai, J.; Lipipun, V.; Likhitwitayawuid, K.; Suzuki, M.; Pramyothin, P.; Hattori, M.; Shiraki, K. Anti-herpes simplex virus (HSV-1) activity of oxyresveratrol derived from Thai medicinal plant: Mechanism of action and therapeutic efficacy on cutaneous HSV-1 infection in mice. Antivir. Res. 2008, 80, 62–70. [Google Scholar] [CrossRef] [PubMed]
  23. Sasivimolphan, P.; Lipipun, V.; Likhitwitayawuid, K.; Takemoto, M.; Pramyothin, P.; Hattori, M.; Shiraki, K. Inhibitory activity of oxyresveratrol on wild-type and drug-resistant varicella-zoster virus replication in vitro. Antivir. Res. 2009, 84, 95–97. [Google Scholar] [CrossRef] [PubMed]
  24. Boonkla, W.; Pitchuanchom, S.; Meepowpan, P.; Thaisuchat, H.; Nuntasaen, N.; Punyanitya, S.; Pompimon, W. Aromatic compound glucopyranoside from new species Artocarpus thailandicus. Am. J. Appl. Sci. 2011, 8, 1093–1097. [Google Scholar] [CrossRef]
  25. Royal Botanic Garden, C.; Carey, W.; Roxburgh, W. Classis XXI. MONOECIA MONANDRIA. Hort. Bengal. 1814, 1, 65–70. Available online: https://www.biodiversitylibrary.org/bibliography/95337#/summary (accessed on 22 March 2020).
  26. Ragone, D. Chromosome numbers and pollen stainability of three species of Pacific Island breadfruit (Artocarpus, Moraceae). Am. J. Bot. 2001, 88, 693–696. [Google Scholar] [CrossRef]
  27. Jarrett, F.M. Studies in Artocarpus and allied genera. III. A revision of Artocarpus subgenus Artocarpus. J. Arnold Arbor. 1959, 40, 113–155. [Google Scholar]
  28. Brantjes, N.B.M. Nectar and pollination of breadfruit, Artocarpus altilis (Moraceae). Acta Bot. Neerl. 1981, 30, 345–352. [Google Scholar] [CrossRef]
  29. Maraci, Ö.; Özkan, H.; Bilgin, R. Phylogeny and genetic structure in the genus Secale. PLoS ONE 2018, 13, e0200825. [Google Scholar] [CrossRef]
  30. Anand, K.K.; Jena, S.N.; Chaudhary, L.B.; Singh, M. Conflict between morphological and molecular data: A case study of Ficus krishnae (Moraceae). Phytotaxa 2016, 247, 143–147. [Google Scholar] [CrossRef] [Green Version]
  31. Dalin, P.; Ågren, J.; Björkman, C.; Huttunen, P.; Kärkkäinen, K. Leaf trichome formation and plant resistance to herbivory. In Induced Plant Resistance to Herbivory; Schaller, A., Ed.; Springer: Dordrecht, The Netherlands, 2008; pp. 89–105. [Google Scholar]
  32. War, A.R.; Paulraj, M.G.; Sharma, H.C. Mechanisms of plant defense against insect herbivores. Plant Signal Behav. 2012, 7, 1306–1320. [Google Scholar] [CrossRef] [Green Version]
  33. Jhonson, H.B. Plant pubescence: An ecological perspective. Bot. Rev. 1975, 41, 233–258. [Google Scholar] [CrossRef]
  34. Chantarasuwan, B.; Baas, P.; van Heuven, B.J.; Baider, C.; van Welzen, P.C. Leaf anatomy of Ficus subsection Urostigma (Moraceae). Bot. J. 2014, 175, 259–281. [Google Scholar] [CrossRef] [Green Version]
  35. Rebeiz, M.; Tsiantis, M. Enhancer evolution and the origins of morphological novelty. Curr. Opin. Genet. Dev. 2017, 45, 115–123. [Google Scholar] [CrossRef] [PubMed]
  36. Das Gupta, M.; Tsiantis, M. Gene networks and the evolution of plant morphology. Curr. Opin. Plant Biol. 2018, 45, 82–87. [Google Scholar] [CrossRef] [PubMed]
  37. Lazarowych, N.J.; Pekos, P. Use of fingerprinting and marker compounds for identification and standardization of botanical drugs: Strategies for Applying Pharmaceutical HPLC Analysis to Herbal Products. Ther. Innov. Regul. Sci. 1998, 32, 497–512. [Google Scholar] [CrossRef]
  38. Giri, L.; Andola, H.; Purohit, V.K.; Rawat, M.S.M.; Rawal, R.; Bhatt, I.D. Chromatographic and spectral fingerprinting standardization of traditional medicines: An overview as modern tools. Res. J. Phytochem. 2010, 4, 234–241. [Google Scholar] [CrossRef] [Green Version]
  39. SimpleMappr, an Online Tool to Produce Publication-Quality Point Maps. Available online: http://www.simplemappr.net (accessed on 2 September 2019).
  40. Chantarasuwan, B.; Berg, C.C.; Kjellberg, F.; Rønsted, N.; Garcia, M.; Baider, C.; van Welzen, P.C. A New Classification of Ficus Subsection Urostigma (Moraceae) Based on Four Nuclear DNA Markers (ITS, ETS, G3pdh, and ncpGS), Morphology and Leaf Anatomy. PLoS ONE 2015, 10, e0128289. [Google Scholar] [CrossRef]
  41. Baldwin, B.G.; Markos, S. Phylogenetic utility of the external transcribed spacer (ETS) of 18S-26S rDNA: Congruence of ETS and ITS trees of Calycadenia (Compositae). Mol. Phylogenet Evol. 1998, 10, 449–463. [Google Scholar] [CrossRef]
  42. Taberlet, P.; Gielly, L.; Pautou, G.; Bouvet, J. Universal primers for amplification of three non-coding regions of chloroplast DNA. Plant Mol. Biol. 1991, 17, 1105–1109. [Google Scholar] [CrossRef]
  43. Sang, T.; Crawford, D.J.; Stuessy, T.F. Chloroplast DNA phylogeny, reticulate evolution and biogeography of Paeonia (Paeoniaceae). Am. J. Bot. 1997, 84, 1120–1136. [Google Scholar] [CrossRef] [Green Version]
  44. Chenna, R.; Sugawara, H.; Tadashi, K.; Lopez, R.; Gibson, T.J.; Higgins, D.G.; Thompson, J.D. Multiple sequence alignment with the Clustal series of programs. Nucleic Acids Res. 2003, 31, 3497–3500. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Nei, M.; Kumar, S. Molecular Evolution and Phylogenetics; Oxford University Press: New York, NY, USA, 2000; pp. 1–333. [Google Scholar]
  47. Tamura, K. Estimation of the number of nucleotide substitutions when there are strong transition-transversion and G + C-content biases. Mol. Biol. Evol. 1992, 9, 678–687. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Kimura, M. A simple method for estimating evolutionary rate of base substitutions through comparative studies of nucleotide sequences. J. Mol. Evol. 1980, 16, 111–120. [Google Scholar] [CrossRef] [PubMed]
  49. Felsenstein, J. Confidence limits on phylogenies: An approach using the bootstrap. Evolution 1985, 39, 783–791. [Google Scholar] [CrossRef] [PubMed]
  50. Huang, H.; Zhang, J.; Chen, G.; Lu, Z.; Wang, X.; Sha, N.; Shao, B.; Li, P.; Guo, D.A. High performance liquid chromatographic method for the determination and pharmacokinetic studies of oxyresveratrol and resveratrol in rat plasma after oral administration of Smilax china extract. Biomed. Chromatogr. 2008, 22, 421–427. [Google Scholar] [CrossRef]
Figure 1. Inflorescences of A. lacucha and A. thailandicus. (A) staminate inflorescence of A. lacucha, subglobose (B) staminate inflorescence of A. thailandicus, obovoid (C) pistillate inflorescence of A. lacucha, obovoid (D) pistillate inflorescence of A. thailandicus, obovoid.
Figure 1. Inflorescences of A. lacucha and A. thailandicus. (A) staminate inflorescence of A. lacucha, subglobose (B) staminate inflorescence of A. thailandicus, obovoid (C) pistillate inflorescence of A. lacucha, obovoid (D) pistillate inflorescence of A. thailandicus, obovoid.
Plants 09 00391 g001
Figure 2. Distribution map of A. lacucha (yellow-circle) and A. thailandicus (pink-star) in Thailand.
Figure 2. Distribution map of A. lacucha (yellow-circle) and A. thailandicus (pink-star) in Thailand.
Plants 09 00391 g002
Figure 3. Upper surface of leaf micromorphology study by SEM of A. lacucha and A. thailandicus. (A) Upper surface of A. lacucha (voucher PBM05214), glabrous (B) Upper surface of A. lacucha (voucher PBM05218), with straight and uncinate smooth hair (C) Upper surface of A. thailandicus, glabrous (voucher PBM05203) (D) Upper surface of A. thailandicus (voucher PBM05229), with straight and uncinate smooth hair.
Figure 3. Upper surface of leaf micromorphology study by SEM of A. lacucha and A. thailandicus. (A) Upper surface of A. lacucha (voucher PBM05214), glabrous (B) Upper surface of A. lacucha (voucher PBM05218), with straight and uncinate smooth hair (C) Upper surface of A. thailandicus, glabrous (voucher PBM05203) (D) Upper surface of A. thailandicus (voucher PBM05229), with straight and uncinate smooth hair.
Plants 09 00391 g003
Figure 4. Lower surface of leaf micromorphology study by SEM of A. lacucha and A. thailandicus. (A,B) Lower surface of A. lacucha (PBM05214), gnarled papillae on the veinlets, sparsely smooth straight or uncinate hairs on lateral veins. (C,D) Lower surface of A. lacucha (PBM05225), gnarled long straight hairs mostly on veinlets. (E,F) Lower surface of A. thailandicus (PBM05229), gnarled long straight hairs mostly on the veinlets with tomentose hair in the areoles and densely smooth straight or uncinate hairs on the lateral veins. (G,H) Lower surface of A. thailandicus (PBM05205), gnarled papillae on the veinlets, lack of tomentose hair in the areoles.
Figure 4. Lower surface of leaf micromorphology study by SEM of A. lacucha and A. thailandicus. (A,B) Lower surface of A. lacucha (PBM05214), gnarled papillae on the veinlets, sparsely smooth straight or uncinate hairs on lateral veins. (C,D) Lower surface of A. lacucha (PBM05225), gnarled long straight hairs mostly on veinlets. (E,F) Lower surface of A. thailandicus (PBM05229), gnarled long straight hairs mostly on the veinlets with tomentose hair in the areoles and densely smooth straight or uncinate hairs on the lateral veins. (G,H) Lower surface of A. thailandicus (PBM05205), gnarled papillae on the veinlets, lack of tomentose hair in the areoles.
Plants 09 00391 g004
Figure 5. Maximum parsimony (MP) and maximum likelihood (ML) analysis of A. lacucha and A. thailandicus with other Artocarpus species. A. integer and A. heterophyllus were chosen as outgroup. The bootstrap percentage of MP and ML were shown as above and below the branches, respectively. Two minor clades, the blue and red square dashed line, indicated that both species were clustered together in the same clade with low bootstrap support. Some parts of the trees showed unresolved topology. Some samples of A. lacucha formed only one well-supported clade with 99% bootstrap support (indicated by green square dashed line). The SEM results were correlated. Light-blue square: upper surface glabrous; purple square: upper surface with straight and uncinate smooth hair; yellow star: lower surface with gnarled papillae on the veinlets, sparsely smooth straight or uncinate hairs on lateral veins; pink star: lower surface with gnarled long straight hairs mostly on veinlets; grey star: lower surface with gnarled long straight hairs mostly on the veinlets with tomentose hair in the areoles and densely smooth straight or uncinate hairs on the lateral veins.
Figure 5. Maximum parsimony (MP) and maximum likelihood (ML) analysis of A. lacucha and A. thailandicus with other Artocarpus species. A. integer and A. heterophyllus were chosen as outgroup. The bootstrap percentage of MP and ML were shown as above and below the branches, respectively. Two minor clades, the blue and red square dashed line, indicated that both species were clustered together in the same clade with low bootstrap support. Some parts of the trees showed unresolved topology. Some samples of A. lacucha formed only one well-supported clade with 99% bootstrap support (indicated by green square dashed line). The SEM results were correlated. Light-blue square: upper surface glabrous; purple square: upper surface with straight and uncinate smooth hair; yellow star: lower surface with gnarled papillae on the veinlets, sparsely smooth straight or uncinate hairs on lateral veins; pink star: lower surface with gnarled long straight hairs mostly on veinlets; grey star: lower surface with gnarled long straight hairs mostly on the veinlets with tomentose hair in the areoles and densely smooth straight or uncinate hairs on the lateral veins.
Plants 09 00391 g005
Figure 6. HPLC chromatograms of A. lacucha and A. thailandicus twig extract alignment. Oxyresveratrol standard and carbamazepine (internal standard) were detected at the retention time of 17.3 (Figure 6A) and 28.3 min, respectively. The HPLC chromatograms of oxyresveratrol (black line), A. lacucha voucher PBM05224 (blue line) and A. thailandicus voucher PBM05227 (red line) were selected as representative of each species and aligned in order to comparison.
Figure 6. HPLC chromatograms of A. lacucha and A. thailandicus twig extract alignment. Oxyresveratrol standard and carbamazepine (internal standard) were detected at the retention time of 17.3 (Figure 6A) and 28.3 min, respectively. The HPLC chromatograms of oxyresveratrol (black line), A. lacucha voucher PBM05224 (blue line) and A. thailandicus voucher PBM05227 (red line) were selected as representative of each species and aligned in order to comparison.
Plants 09 00391 g006
Table 1. The main differentiating characters between A. lacucha and A. thailandicus.
Table 1. The main differentiating characters between A. lacucha and A. thailandicus.
CharactersA. lacuchaA. thailandicus
Indumentum of the areoles on the lower surface of laminaAbsentPresent
Shape and size of staminate inflorescenceMostly subglobose
0.6–2 cm diam.
Mostly obovoid, cylindrical, or clavate
0.2–0.4 cm diam.
Size of pistillate inflorescence1.5 cm diam.0.7–1 cm diam.
The length of peduncle of staminate inflorescence0.2–2 cmSessile or up to 0.1 cm
The length of peduncle of pistillate inflorescence0.2–4 cmUp to 0.7 cm
Table 2. List of the studied materials.
Table 2. List of the studied materials.
TaxonFloristic RegionLocalityPBM NumberCollector Number
A. lacuchaCentralSaraburi, Chaloem Phra Kiat, lat 14.67633° N lon 100.8868° E05197Aneklaphakij 2
CentralSaraburi, Chaloem Phra Kiat, lat 14.67217° N lon 100.887° E05198Aneklaphakij 4
CentralSaraburi, Chaloem Phra Kiat, lat 14.66953° N lon 100.8898° E05199Aneklaphakij 5
CentralNakhon Pathom, Phutthamonthon, lat 13.79072° N lon 100.317° E05200Aneklaphakij 6
CentralNakhon Pathom, Phutthamonthon, lat 13.79072° N lon 100.317° E05201Aneklaphakij 7
North-easternKhon Kaen, Wiang Kao, lat 16.79369° N lon 102.2577° E05206Aneklaphakij 14
North-easternKhon Kaen, Wiang Kao, lat 16.79369° N lon 102.2577° E05207Aneklaphakij 15
North-easternKhon Kaen, Phon, lat 15.88319° N lon 102.5332° E05214Aneklaphakij 23
North-easternKhon Kaen, Phon, lat 15.88319° N lon 102.5332° E05215Aneklaphakij 24
North-easternKhon Kaen, Phon, lat 15.88319° N lon 102.5332° E05216Aneklaphakij 25
North-easternKhon Kaen, Phon, lat 15.88319° N lon 102.5332° E05217Aneklaphakij 26
North-easternKhon Kaen, Phon, lat 15.88319° N lon 102.5332° E05218Aneklaphakij 27
CentralBangkok, Prawet, lat 13.68917° N lon 100.6596° E05221Aneklaphakij 31
NorthernChiang Mai, Mae Rim, lat 18.89464° N lon 98.85864° E05223Aneklaphakij 36
NorthernChiang Mai, Mae Rim, lat 18.89456° N lon 98.85869° E05224Aneklaphakij 37
NorthernChiang Mai, Mae Rim, lat 18.89133° N lon 98.86214° E05225Aneklaphakij 38
PeninsularSongkhla, Hat Yai, lat 7.008417° N lon 100.5059° E05230Aneklaphakij 44
PeninsularSongkhla, Hat Yai, lat 7.008528° N lon 100.5062° E05231Aneklaphakij 47
A. thailandicusSouth-westernPhetchaburi, Kaeng Krachan, lat 12.9025° N lon 99.64325° E05202Aneklaphakij 9
EasternUbon Ratchathani, Samrong, lat 15.01044° N lon 104.7846° E05203Aneklaphakij 10
EasternUbon Ratchathani, Samrong, lat 15.01044° N lon 104.7846° E05204Aneklaphakij 11
South-easternPrachin Buri, Mueang Prachin Buri, lat 14.10967° N lon 101.4367° E05205Aneklaphakij 13
North-easternKhon Kaen, Mueang Khon Kaen, lat 16.4695° N lon 102.8271° E05208Aneklaphakij 17
North-easternKhon Kaen, Mueang Khon Kaen, lat 16.47194° N lon 102.8218° E05209Aneklaphakij 18
North-easternKhon Kaen, Mueang Khon Kaen, lat 16.47706° N lon 102.8258° E05210Aneklaphakij 19
North-easternKhon Kaen, Mueang Khon Kaen, lat 16.47706° N lon 102.8258° E05211Aneklaphakij 20
North-easternKhon Kaen, Mueang Khon Kaen, lat 16.47975° N lon 102.8276° E05212Aneklaphakij 21
North-easternKhon Kaen, Mueang Khon Kaen, lat 16.47447° N lon 102.8254° E05213Aneklaphakij 22
CentralBangkok, Prawet, lat 13.68853° N lon 100.6606° E05219Aneklaphakij 28
CentralBangkok, Prawet, lat 13.68853° N lon 100.6606° E05220Aneklaphakij 30
NorthernChiang Mai, Mae Rim, lat 18.89742° N lon 98.85928° E05222Aneklaphakij 33
NorthernChiang Mai, Mae Rim, lat 18.89136° N lon 98.8625° E05226Aneklaphakij 39
NorthernChiang Mai, Mae Rim, lat 18.88869° N lon 98.86169° E05227Aneklaphakij 40
NorthernChiang Mai, Mae Rim, lat 18.88869° N lon 98.86169° E05228Aneklaphakij 41
NorthernChiang Mai, Mae Rim, lat 18.89° N lon 98.85842° E05229Aneklaphakij 42
Table 3. Primer sequences used in the study.
Table 3. Primer sequences used in the study.
RegionPrimer SequenceReferences
ITSForward: 5’–AACAAGGTTTCCGTAGGTGA–3’
Reverse: 5’–TATGCTTAAAYTCAGCGGGT–3’
[2]
ETSForward: 5’–GACCCTTGGTTCCTGTGTTGC–3’
Reverse: 5’–ACTTACACATGCATGGCTTAATCT–3’
[40,41]
trnL-Fprimer c: 5’–CGAAATCGGTAGACGCTACG–3’
primer d: 5’–GGGGATAGAGGGACTTGAAC–3’
primer e: 5’–GGTTCAAGTCCCTCTATCCC–3’
primer f: 5’–ATTTGAACTGGTGACACGAG–3’
[42]
trnH-psbAForward: 5’–CGCGCATGGTGGATTCACAAATC–3’
Reverse: 5’–GTTATGCATGAACGTAATGCTC–3’
[43]
Table 4. Accession number of the sequences deposited in NCBI GenBank database.
Table 4. Accession number of the sequences deposited in NCBI GenBank database.
PBM NumberITSETStrnL-FtrnH-psbA
05197MK850217MK873914MK873937MK873960
05198MK850218MK873915MK873938MK873961
05199MH801139MH807231MH817052MH817064
05200MK850219MK873916MK873939MK873962
05201MK850220MK873917MK873940MK873963
05202MH801140MH807232MH817053MH817065
05203MH801141MH807233MH817054MH817066
05204MK850221MK873918MK873941MK873964
05205MH801142MH807234MH817055MH817067
05206MH801143MH807235MH817056MH817068
05207MK850222MK873919MK873942MK873965
05208MK850223MK873920MK873943MK873966
05209MK850224MK873921MK873944MK873967
05210MH801144MH807236MH817057MH817069
05211MK850225MK873922MK873945MK873968
05212MK850226MK873923MK873946MK873969
05213MK850227MK873924MK873947MK873970
05214MK850228MK873925MK873948MK873971
05215MK850229MK873926MK873949MK873972
05216MK850230MK873927MK873950MK873973
05217MK850231MK873928MK873951MK873974
05218MH801145MH807237MH817058MH817070
05219MH801146MH807238MH817059MH817071
05220MK850232MK873929MK873952MK873975
05221MH801147MH807239MH817060MH817072
05222MK850233MK873930MK873953MK873976
05223MK850234MK873931MK873954MK873977
05224MK850235MK873932MK873955MK873978
05225MH801148MH807240MH817061MH817073
05226MK850236MK873933MK873956MK873979
05227MK850237MK873934MK873957MK873980
05228MK850238MK873935MK873958MK873981
05229MH801149MH807241MH817062MH817074
05230MK850239MK873936MK873959MK873982
05231MH801150MH807242MH817063MH817075

Share and Cite

MDPI and ACS Style

Aneklaphakij, C.; Bunsupa, S.; Sirichamorn, Y.; Bongcheewin, B.; Satitpatipan, V. Taxonomic Notes on the ‘Mahat’ (Artocarpus lacucha and A. thailandicus, Moraceae) Species Complex in Thailand. Plants 2020, 9, 391. https://doi.org/10.3390/plants9030391

AMA Style

Aneklaphakij C, Bunsupa S, Sirichamorn Y, Bongcheewin B, Satitpatipan V. Taxonomic Notes on the ‘Mahat’ (Artocarpus lacucha and A. thailandicus, Moraceae) Species Complex in Thailand. Plants. 2020; 9(3):391. https://doi.org/10.3390/plants9030391

Chicago/Turabian Style

Aneklaphakij, Chaiwat, Somnuk Bunsupa, Yotsawate Sirichamorn, Bhanubong Bongcheewin, and Veena Satitpatipan. 2020. "Taxonomic Notes on the ‘Mahat’ (Artocarpus lacucha and A. thailandicus, Moraceae) Species Complex in Thailand" Plants 9, no. 3: 391. https://doi.org/10.3390/plants9030391

APA Style

Aneklaphakij, C., Bunsupa, S., Sirichamorn, Y., Bongcheewin, B., & Satitpatipan, V. (2020). Taxonomic Notes on the ‘Mahat’ (Artocarpus lacucha and A. thailandicus, Moraceae) Species Complex in Thailand. Plants, 9(3), 391. https://doi.org/10.3390/plants9030391

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop