A large-scale sampling was undertaken during a research cruise across the South China Sea in August 2016, covering an area of about 100,000 km2 to investigate the molecular diversity and distributions of micro-eukaryotic protists, with a focus on the potentially harmful microalgal (HAB) species along the east coast of Peninsular Malaysia. Environmental DNAs from 30 stations were extracted and DNA metabarcoding targeting the V4 and V9 markers in the 18S rDNA was performed. Many protistan molecular units, including previously unreported HAB taxa, were discovered for the first time in the water. Our findings also revealed interesting spatial distribution patterns, with a marked signal of compositional turnover between latitudinal regimes of water masses, where dinophytes and diatom compositions were among the most strongly enhanced at the fronts, leading to distinct niches. Our results further confirmed the widespread distribution of HAB species, such as the toxigenic Alexandrium tamiyavaichii and Pseudo-nitzschia species, and the fish-killing Margalefidinium polykrikoides and Karlodinium veneficum. The molecular information obtained from this study provides an updated HAB species inventory and a toolset that could facilitate existing HAB monitoring schemes in the region to better inform management decisions.
This study describes two novel species of marine dinophytes in the genus Alexandrium. Morphological characteristics and phylogenetic analyses support the placement of the new taxa, herein designated as Alexandrium limii sp. nov. and A. ogatae sp. nov. Alexandrium limii, a species closely related to A. taylorii, is distinguished by having a shorter 2'/4' suture length, narrower plates 1' and 6'', with larger length: width ratios, and by the position of the ventral pore (Vp). Alexandrium ogatae is distinguishable with its metasert plate 1' having almost parallel lateral margins, and by lacking a Vp. Production of paralytic shellfish toxins (PSTs), cycloimines, and goniodomins (GDs) in clonal cultures of A. ogatae, A. limii, and A. taylorii were examined analytically and the results showed that all strains contained GDs, with GDA as major variants (6-14 pg cell-1) for all strains except the Japanese strain of A. limii, which exclusively had a desmethyl variant of GDA (1.4-7.3 pg cell-1). None of the strains contained detectable levels of PSTs and cycloimines.
Coastal eutrophication is one of the pivotal factors driving occurrence of harmful algal blooms (HABs), whose underlying mechanism remained unclear. To better understand the nutrient regime triggering HABs and their formation process, the phytoplankton composition and its response to varying nitrogen (N) and phosphorus (P), physio-chemical parameters in water and sediment in Johor Strait in March 2019 were analyzed. Surface and sub-surface HABs were observed with the main causative species of Skeletonema, Chaetoceros and Karlodinium. The ecophysiological responses of Skeletonema to the low ambient N/P ratio such as secreting alkaline phosphatase, regulating cell morphology (volume; surface area/volume ratio) might play an important role in dominating the community. Anaerobic sediment iron-bound P release and simultaneous N removal by denitrification and anammox, shaped the stoichiometry of N and P in water column. The decrease of N/P ratio might shift the phytoplankton community into the dominance of HABs causative diatoms and dinoflagellates.
Microplastic pollution is a global issue that has a detrimental impact on food safety. In marine environments, microplastics are a threat to marine organisms, as they are often the same size range as prey and are mistaken as food. Consumption of microplastics has led to the damage of digestive organs and a reduction in growth and reproductive output. In this study, microplastic pollution was assessed across three commercially available shrimp species that were obtained from the supermarkets of Singapore. A total of 93 individuals were studied from the Pacific white leg shrimp, Litopenaeus vannamei, the Argentine red shrimp Pleoticus muelleri and the Indian white shrimp Fenneropenaeus indicus. Microplastic fibers, fragments, film and spheres were identified from the digestive tract of these organisms. Microplastic abundance ranged from 13.4 to 7050 items. F. indicus exhibited the highest number of microplastics. Microplastic film was the most abundant in L. vannamei individuals (93-97%) and spheres were the most abundant in P. muelleri (70%) and F. indicus (61%) individuals. This study demonstrates that microplastic contamination is evident in commonly consumed shrimp and highlights the role of shrimp in the trophic transfer and accumulation of microplastics in seafood. The consumption of microplastic-containing seafood is a route of exposure to humans and has implications on human health and food security. Capsule: Microplastics were examined in three shrimp species from the supermarkets of Singapore. Microplastics ranged from 13.4 to 7050 items of shrimp.
Recent molecular phylogenetic studies of Gambierdiscus species flagged several new species and genotypes, thus leading to revitalizing its systematics. The inter-relationships of clades revealed by the primary sequence information of nuclear ribosomal genes (rDNA), however, can sometimes be equivocal, and therefore, in this study, the taxonomic status of a ribotype, Gambierdiscus sp. type 6, was evaluated using specimens collected from the original locality, Marakei Island, Republic of Kiribati; and specimens found in Rawa Island, Peninsular Malaysia, were further used for comparison. Morphologically, the ribotype cells resembled G. scabrosus, G. belizeanus, G. balechii, G. cheloniae and G. lapillus in thecal ornamentation, where the thecal surfaces are reticulate-foveated, but differed from G. scabrosus by its hatchet-shaped Plate 2', and G. belizeanus by the asymmetrical Plate 3'. To identify the phylogenetic relationship of this ribotype, a large dataset of the large subunit (LSU) and small subunit (SSU) rDNAs were compiled, and performed comprehensive analyses, using Bayesian-inference, maximum-parsimony, and maximum-likelihood, for the latter two incorporating the sequence-structure information of the SSU rDNA. Both the LSU and SSU rDNA phylogenetic trees displayed an identical topology and supported the hypothesis that the relationship between Gambierdiscus sp. type 6 and G. balechii was monophyletic. As a result, the taxonomic status of Gambierdiscus sp. type 6 was revised, and assigned as Gambierdiscus balechii. Toxicity analysis using neuroblastoma N2A assay confirmed that the Central Pacific strains were toxic, ranging from 1.1 to 19.9 fg P-CTX-1 eq cell-1, but no toxicity was detected in a Western Pacific strain. This suggested that the species might be one of the species contributing to the high incidence rate of ciguatera fish poisoning in Marakei Island.
In the coastal countries of Southeast Asia, fish is a staple diet and certain fish species are food delicacies to local populations or commercially important to individual communities. Although there have been several suspected cases of ciguatera fish poisoning (CFP) in Southeast Asian countries, few have been confirmed by ciguatoxins identification, resulting in limited information for the correct diagnosis of this food-borne disease. In the present study, ciguatoxin-1B (CTX-1B) in red snapper (Lutjanus bohar) implicated in a CFP case in Sabah, Malaysia, in December 2017 was determined by single-quadrupole selected ion monitoring (SIM) liquid chromatography/mass spectrometry (LC/MS). Continuous consumption of the toxic fish likely resulted in CFP, even when the toxin concentration in the fish consumed was low. The identification of the fish species was performed using the molecular characterization of the mitochondrial cytochrome c oxidase subunit I gene marker, with a phylogenetic analysis of the genus Lutjanus. This is the first report identifying the causative toxin in fish-implicated CFP in Malaysia.
Effects of aquaculture activities on the environmental parameters and phytoplankton community structure were investigated in a semi-enclosed lagoon located at Semerak River, Malaysia. Elevated concentrations of phosphate and ammonia were observed at the aquaculture area and the inner lagoon. Relatively low dissolved oxygen, high total chlorophyll a, and high phytoplankton abundances but low species richness were recorded. Chaetoceros, Pseudo-nitzschia brasiliana, Blixaea quinquecornis, and Skeletonema blooms were observed, and some were associated with anoxia condition. Eutrophication level assessed by UNTRIX suggests that the water quality in the lagoon is deteriorating. Dissolved inorganic phosphorus and nitrogen at the impacted area were 15 and 12 times higher than the reference sites, respectively. Such trophic status indices could provide a useful guideline for optimal aquaculture management plan to reduce the environmental impact caused by aquaculture.
Information on the diversity and distribution of harmful microalgae in the Gulf of Thailand is very limited and mainly based on microscopic observations. Here, we collected 44 water samples from the Gulf of Thailand and its adjacent water (Perhentian Island, Malaysia) for comparison in 2018. DNA metabarcoding was performed targeting the partial large subunit ribosomal RNA gene (LSU rDNA D1-D3) and the internal transcribed spacers (ITS1 and ITS2). A total of 50 dinoflagellate genera (made up of 72 species) were identified based on the LSU rDNA dataset, while the results of ITS1 and ITS2 datasets revealed 33 and 32 dinoflagellate genera comprising 69 and 64 species, respectively. Five potentially toxic Pseudo-nitzschia (Bacillariophyceae) species were detected, with four as newly recorded species in the water (Pseudo-nitzschia americana/brasilliana, Pseudo-nitzschia simulans/delicatissima, P. galaxiae and P. multistriata). The highest relative abundances of P. galaxiae and P. multistriata were found in Trat Bay and Chumphon (accounting for 0.20% and 0.06% of total ASVs abundance, respectively). Three paralytic shellfish toxin producing dinoflagellate species were detected: Alexandrium tamiyavanichii, Alexandrium fragae, and Gymnodinium catenatum. The highest abundance of A. tamiyavanichii was found in the surface sample of Chumphon (CHO7 station), accounting for 1.95% of total ASVs abundance. Two azaspiracid producing dinoflagellate species, Azadinium poporum ribotype B, Azadinium spinosum ribotype A, and a pinnatoxin producing dinoflagellate species Vulcanodinium rugosum, with two ribotypes B and C, were revealed from the datasets although with very low abundances. Six fish killing dinoflagellate species, including Margalefidinium polykrikoides group IV, Margalefidinium fulvescens, Karenia mikimotoi, Karenia selliformis ribotype B, Karlodinium australe, and Karlodinium digitatum were detected and all representing new records in this area. The findings of numerous harmful microalgal species in the Gulf of Thailand highlight the potential risk of human intoxication and fish killing events.
Thirty-four strains of Heterocapsa were established from Malaysian waters and their morphologies were examined by light, scanning, and transmission electron microscopy. Three species, H. bohaiensis, H. huensis, and H. rotundata, and three new species, H. borneoensis sp. nov., H. limii sp. nov., and H. iwatakii sp. nov. were described in this study. The three species were differentiated morphologically by unique characteristics of cell size, shape, displacement of the cingulum, shape and position of nucleus, the number and position of pyrenoids, and body scale ultrastructure. The species delimitations were robustly supported by the molecular data. A light-microscopy-based key to species of Heterocapsa is established, with two major groups, i.e., species with a single pyrenoid, and species with multiple pyrenoids. Bioassays were conducted by exposing Artemia nauplii to Heterocapsa densities of 1-5 × 105 cells mL-1, and treatments exposed to H. borneoensis showed naupliar mortality, while no naupliar death was observed in the treatments exposed to cells of H. bohaiensis, H. huensis, H. limii, and H. iwatakii. Naupliar death was observed during the initial 24 h for both tested H. borneoensis strains, and mortality rates increased up to 50% after 72-h exposure. This study documented for the first time the diversity and cytotoxic potency of Heterocapsa species from Malaysian waters.
Coastal ecosystems are often subjected to anthropogenic disturbances that lead to water quality deterioration and an increase in harmful algal bloom (HAB) events. Using the next-generation molecular tool of 18S rDNA metabarcoding, we examined the community assemblages of HAB species in the Johor Strait, Malaysia between May 2018 and September 2019, covering 19 stations across the strait. The molecular operational taxonomic units (OTUs) of HAB taxa retrieved from the dataset (n = 194) revealed a much higher number of HAB taxa (26 OTUs) than before, with 12 taxa belong to new records in the strait. As revealed in the findings of this study, the diversity and community structure of HAB taxa varied significantly over time and space. The most common and abundant HAB taxa in the strait (frequency of occurrence >70%) comprised Heterosigma akashiwo, Fibrocapsa japonica, Pseudo-nitzschia pungens, Dinophysis spp., Gymnodinium catenatum, Alexandrium leei, and A. tamiyavanichii. Also, our results demonstrated that the HAB community assemblages in the strait were dependent on the interplay of environmental variables that influence by the monsoonal effects. Different HAB taxa, constitute various functional types, occupied and prevailed in different environmental niches across space and time, leading to diverse community assemblages and population density. This study adds to the current understandings of HAB dynamics and provides a robust overview of temporal-spatial changes in HAB community assemblages along the environmental gradients in a tropical eutrophic coastal ecosystem.
Saxitoxins (STXs) constitute a family of potent sodium channel blocking toxins, causative agents of paralytic shellfish poisoning (PSP), and are produced by several species of marine dinoflagellates and cyanobacteria. Two STX-core genes, sxtA and sxtG, have been well elucidated in Alexandrium but the expression of these genes under various nutritional modes in tropical species remains unclear. This study investigates the physiological responses of a tropical Pacific strain of Alexandrium minutum growing with nitrate or ammonium, and with various nitrogen to phosphorus (N:P) supply ratios. The transcriptional responses of the sxt genes were observed. Likewise, a putative sxtI encoding O-carbamoyltransferase (herein designated as AmsxtI) was recovered from the transcriptomic data, and its expression was investigated. The results revealed that the cellular toxin quota (Qt) was higher in P-depleted, nitrate-grown cultures. With cultures at similar N:P (<16), cells grown with excess ammonium showed a higher Qt than those grown with nitrate. sxtA1 was not expressed under any culture conditions, suggesting that this gene might not be involved in STX biosynthesis by this strain. Conversely, sxtA4 and sxtG showed positive correlations with Qt, and were up-regulated in P-depleted, nitrate-grown cultures and with excess ambient ammonium. On the other hand, AmsxtI was expressed only when induced by P-depletion, suggesting that this gene may play an important role in P-recycling metabolism, while simultaneously enhancing toxin production.
The marine dinoflagellate Alexandrium minutum is known to produce saxitoxins that cause paralytic shellfish poisoning in human worldwide through consumption of the contaminated shellfish mollusks. Despite numerous studies on the growth physiology and saxitoxin production of this species, the knowledge on the molecular basis of nutrient uptakes in relation to toxin production in this species is limited. In this study, relative expressions of the high-affinity transporter genes of nitrate, ammonium, and phosphate (AmNrt2, AmAmt1 and AmPiPT1) and the assimilation genes, nitrate reductase (AmNas), glutamine synthase (AmGSIII) and carbamoyl phosphate synthase (AmCPSII) from A. minutum were studied in batch clonal culture condition with two nitrogen sources (nitrate: NO3- or ammonium: NH4+) under different N:P ratios (high-P: N:P of 14 and 16, and low-P: N:P of 155). The expression of AmAmt1 was suppressed in excess NH4+-grown condition but was not observed in AmNrt2 and AmNas. Expressions of AmAmt1, AmNrt2, AmNas, AmGSIII, AmCPSII, and AmPiPT1 were high in P-deficient condition, showing that A. minutum is likely to take up nutrients for growth under P-stress condition. Conversely, relative expression of AmCPSII was incongruent with cell growth, but was well correlated with toxin quota, suggesting that the gene might involve in arginine metabolism and related toxin production pathway. The expression of AmGSIII is found coincided with higher toxin production and is believed to involve in mechanism to detoxify the cells from excess ammonium stress. The gene regulation observed in this study has provided better insights into the ecophysiology of A. minutum in relation to its adaptive strategies in unfavorable environments.
To assess the effects of environmental changes on phytoplankton community structure in a mangrove ecosystem, phytoplankton distribution in Matang mangrove, Malaysia was examined. Phytoplankton and water samples, and in situ environmental parameters from three estuaries with differing levels of disturbance were examined monthly for one year. Two species, Cyclotella choctawhatcheeana and Skeletonema costatum, were dominant in the least disturbed and moderately disturbed areas, respectively. Skeletonema costatum was also the most dominant in the most disturbed area. Significant differences in phytoplankton density and biodiversity between the least and most disturbed areas were also observed. Principle component 1 (salinity, conductivity, total solids/water transparency and nitrogenous compounds) and PC2 (dissolved oxygen, pH and temperature) explained 60.4% of the total variance. This study illustrated that changes in phytoplankton community structure in Matang mangrove estuaries were significantly correlated with environmental parameters which were in turn influenced by ecosystem disturbance levels as well as seasonal changes.
In a field survey in the Taiwan Strait during April 2016, the species composition and the domoic acid production of the diatom genus Pseudo-nitzschia were investigated. A total of 80 strains of Pseudo-nitzschia were established, and species identification was determined based on a combination of morphological and molecular data. Fourteen taxa were recognized, i.e., P. americana, P. brasiliana, P. calliantha, P. cuspidata, P. galaxiae, P. lundholmiae, P. multiseries, P. multistriata, P. pseudodelicatissima, P. pungens var. aveirensis, P. pungenus var. pungens and P. sabit, as well as two novel species P. chiniana C.X. Huang & Yang Li and P. qiana C.X. Huang & Yang Li. Morphologically, P. chiniana is characterized by striae comprising one or two rows of poroids, and valve ends that are normally dominated by two rows of poroids within each stria. Whereas P. qiana is unique by having a narrow valve width (1.3-1.5 μm) and sharply pointed valve ends. Both taxa constitute their own monophyletic lineage in the phylogenetic analyses inferred from LSU and ITS2 rDNA, and are well differentiated from other Pseudo-nitzschia species. Pseudo-nitzschia chiniana forms a group with P. abrensis and P. batesiana in LSU and ITS trees, whereas P. qiana is sister to P. lineola. When comparing ITS2 secondary structure, five CBCs and seven HCBCs are recognized between P. chiniana and P. abrensis, and four CBCs and ten HCBCs between P. chiniana and P. batesiana. Two CBCs and eight HCBCs are found between P. qiana with P. lineola. The ability of the strains to produce domoic acid was assessed, including a potential toxin induction by the presence of brine shrimps. Results revealed production of domoic acid in six strains belonging to three species. Without presence of brine shrimps, cellular DA (pDA) was detected in four P. multiseries strains (1.6 ± 0.3, 26.6 ± 2.7, 68.3 ± 4.2 and 56.9 ± 4.7 fg cell-1, separately), one strain of P. pseudodelicatissima (0.8 ± 0.2 fg cell-1) and one strain of P. lundholmiae (2.5 ± 0.4 fg cell-1). In the presence of brine shrimps, pDA contents increased significantly (p
Different clades belonging to the cosmopolitan marine diatom Pseudo-nitzschia pungens appear to be present in different oceanic environments, however, a 'hybrid zone', where populations of different clades interbreed, has also been reported. Many studies have investigated the sexual reproduction of P. pungens, focused on morphology and life cycle, rather than the role of sexual reproduction in mixing the genomes of their parents. We carried out crossing experiments to determine the sexual compatibility/incompatibility between different clades of P. pungens, and examined the genetic polymorphism in the ITS2 region. Sexual reproduction did not occur only between clades II and III under any of experimental temperature conditions. Four offspring strains were established between clade I and III successfully. Strains established from offspring were found interbreed with other offspring strains as well as viable with their parental strains. We confirmed the hybrid sequence patterns between clades I and III and found novel sequence types including polymorphic single nucleotide polymorphisms (SNPs) in the offspring strains. Our results implicate that gene exchange and mixing between different clades are still possible, and that sexual reproduction is a significant ecological strategy to maintain the genetic diversity within this diatom species.
Pseudo-nitzschia pungens is a widely distributed marine pennate diatom. Hybrid zones, regions in which two different genotypes may interbreed, are important areas for speciation and ecology, and have been reported across the globe for this species. However, sexual reproduction between differing clades in the natural environment is yet to be observed and is difficult to predict. Here we carried out experiments using two mono-clonal cultures of P. pungens from different genotypes to measure the frequency and timing of sexual reproduction across varying biotic (growth phases and cell activity potential) and abiotic conditions (nutrients, light, turbulence). We found the mating rates and number of zygotes gradually decreased from exponential to late stationary growth phases. The maximum zygote abundance observed was 1,390 cells mL-1 and the maximum mating rate was 7.1%, both which occurred during the exponential growth phase. Conversely, only 9 cells mL-1 and a maximum mating rate of 0.1% was observed during the late stationary phase. We also found the higher the relative potential cell activity (rPCA) in parent cells, as determined by the concentration of chlorophyll a per cell and the ratio of colony formation during parent cultivations, revealed higher mating rates. Furthermore, sexual events were reduced under nutrient enrichment conditions, and mating pairs and zygotes were not formed under aphotic (dark) or shaking culture conditions (150 rpm). In order to understand the sexual reproduction of Pseudo-nitzschia in the natural environment, our results highlight that it is most likely the combination of both biotic (growth phase, Chl. a content) and abiotic factors (nutrients, light, turbulence) that will determine the successful union of intraspecific populations of P. pungens in any given region.
In 2015, a remarkably high density bloom of Alexandrium minutum occurred in Sungai Geting, a semi-enclosed lagoon situated in the northeast of Peninsular Malaysia, causing severe discoloration and contaminated the benthic clams (Polymesoda). Plankton and water samples were collected to investigate the mechanisms of bloom development of this toxic species. Analysis of bloom samples using flow cytometry indicated that the bloom was initiated by the process of active excystment, as planomycetes (>4C cells) were observed in the early stage of the bloom. Increase in planozygotes (2C cells) was evident during the middle stage of the bloom, coinciding with an abrupt decrease in salinity and increase of temperature. The bloom was sustained through the combination of binary division of vegetative cells, division of planozygotes, and cyst germination through continuous excystment. Nutrient depletion followed by precipitation subsequently caused the bloom to terminate. This study provides the first continuous record of in situ life-cycle stages of a natural bloom population of A. minutum through a complete bloom cycle. The event has provided a fundamental understanding of the pelagic life-cycle stages of this tropical dinoflagellate, and demonstrated a unique bloom development characteristic shared among toxic Alexandrium species in coastal embayments.
Species of the genus Chattonella (Raphidophyceae) are a group of marine protists that are commonly found in coastal waters. Some are known as harmful microalgae that form noxious blooms and cause massive fish mortality in finfish aquaculture. In Malaysia, blooms of Chattonella have been recorded since the 1980s in the Johor Strait. In this study, two strains of Chattonella were established from the strait, and morphological examination revealed characteristics resembling Chattonella subsalsa. The molecular characterization further confirmed the species' identity as C. subsalsa. To precisely detect the cells of C. subsalsa in the environment, a whole-cell fluorescence in-situ hybridisation (FISH) assay was developed. The species-specific oligonucleotide probes were designed in silico based on the nucleotide sequences of the large subunit (LSU) and internal transcribed spacer 2 (ITS2) of the ribosomal DNA (rDNA). The best candidate signature regions in the LSU-rRNA and ITS2-rDNA were selected based on hybridisation efficiency and probe parameters. The probes were synthesised as biotinylated probes and tested by tyramide signal amplification with FISH (FISH-TSA). The results showed the specificity of the probes toward the target cells. FISH-TSA has been proven to be a potential tool in the detection of harmful algae in the environment and could be applied to the harmful algal monitoring program.
Oral Squamous Cell Carcinoma (OSCC) remains a cancer with poor prognosis and high recurrence rate. Even with multimodal treatment options available for OSCC, tumor drug resistance is still a persistent problem, leading to increased tumor invasiveness among OSCC patients. An emerging trend of thought proposes that extracellular vesicles (EVs) play a role in facilitating tumor progression and chemoresistance via signaling between tumor cells. In particular, exosomes and microvesicles are heavily implicated in this process by various studies. Where primary studies into a particular EV-mediated chemoresistance mechanism in OSCC are limited, similar studies on other cancer cell types will be used in the discussion below to provide ideas for a new line of investigation into OSCC chemoresistance. By understanding how EVs are or may be involved in OSCC chemoresistance, novel targeted therapies such as EV inhibition may be an effective alternative to current treatment options in the near future. In this review, the current understandings on OSCC drug mechanisms under the novel context of exosomes and microvesicles were reviewed, including shuttling of miRNA content, drug efflux, alteration of vesicular pH, anti-apoptotic signaling, modulation of DNA damage repair, immunomodulation, epithelial-to-mesenchymal transition and maintenance of tumor by cancer stem cells.
In this study, inter- and intraspecific genetic diversity within the marine harmful dinoflagellate genus Coolia Meunier was evaluated using isolates obtained from the tropics to subtropics in both Pacific and Atlantic Ocean basins. The aim was to assess the phylogeographic history of the genus and to clarify the validity of established species including Coolia malayensis. Phylogenetic analysis of the D1-D2 LSU rDNA sequences identified six major lineages (L1-L6) corresponding to the morphospecies Coolia malayensis (L1), C. monotis (L2), C. santacroce (L3), C. palmyrensis (L4), C. tropicalis (L5), and C. canariensis (L6). A median joining network (MJN) of C. malayensis ITS2 rDNA sequences revealed a total of 16 haplotypes; however, no spatial genetic differentiation among populations was observed. These MJN results in conjunction with CBC analysis, rDNA phylogenies and geographical distribution analyses confirm C. malayensis as a distinct species which is globally distributed in the tropical to warm-temperate regions. A molecular clock analysis using ITS2 rDNA revealed the evolutionary history of Coolia dated back to the Mesozoic, and supports the hypothesis that historical vicariant events in the early Cenozoic drove the allopatric differentiation of C. malayensis and C. monotis.