Benefits and Limitations of DNA Barcoding and Metabarcoding in Herbal Product Authentication

Introduction

Herbal medicines play an important role globally in the health care sector and in industrialised countries they are often considered as an alternative to mono‐substance medicines. Current quality and authentication assessment methods rely mainly on morphology and analytical phytochemistry‐based methods detailed in pharmacopoeias. Herbal products however are often highly processed with numerous ingredients, and even if these analytical methods are accurate for quality control of specific lead or marker compounds, they are of limited suitability for the authentication of biological ingredients.

Objective

To review the benefits and limitations of DNA barcoding and metabarcoding in complementing current herbal product authentication.

Method

Recent literature relating to DNA based authentication of medicinal plants, herbal medicines and products are summarised to provide a basic understanding of how DNA barcoding and metabarcoding can be applied to this field.

Results

Different methods of quality control and authentication have varying resolution and usefulness along the value chain of these products. DNA barcoding can be used for authenticating products based on single herbal ingredients and DNA metabarcoding for assessment of species diversity in processed products, and both methods should be used in combination with appropriate hyphenated chemical methods for quality control.

Conclusions

DNA barcoding and metabarcoding have potential in the context of quality control of both well and poorly regulated supply systems. Standardisation of protocols for DNA barcoding and DNA sequence‐based identification are necessary before DNA‐based biological methods can be implemented as routine analytical approaches and approved by the competent authorities for use in regulated procedures. © 2017 The Authors. Phytochemical Analysis Published by John Wiley & Sons Ltd.     

Identification of herbal medicinal materials using DNA barcodes

Herbal medicinal materials have been used worldwide for centuries to maintain health and to treat disease. However, adulteration of herbal medicines remains a major concern of users and industry for reasons of safety and efficacy. Identification of herbal medicinal materials by DNA technology has been widely applied,started from the mid-1990s. In recent years, DNA barcoding of global plant species using four standard barcodes (rbcL, matK, trnH-psbA and ITS) has been a major focus in the fields of biodiversity and conservation. These DNA barcodes can also be used as reliable tools to facilitate the identification of herbal medicinal materials for the safe use of herbs, quality control, and forensic investigation. Many studies have applied these DNA barcodes for the identification of herbal medicinal species and their adulterants. The present article reviews efforts in the identification of herbal medicinal materials using the standard DNA barcodes and other DNA sequence-based markers.     

Development of an activation tagging system for the basidiomycetous medicinal fungus Antrodia cinnamomea

This study describes the development of an efficient and reliable activation tagging system for the medicinal fungus Antrodia cinnamomea. For successful Agrobacterium tumefaciens-mediated transformation, different parameters were considered. The Agrobacterium concentration of 5 × 10⁸ cfu ml⁻¹, 1 mm acetosyringone, 25-d-old mycelia at 0.2 g ml⁻¹, and co-culture period of 6 d were found to be the most optimal conditions for enhancing the transformation efficiency. The mitotic stability of transferred DNA (T-DNA) was demonstrated by growing eight randomly selected putative transformants in malt extract agar medium for five subcultures. Insertion of T-DNA into the genome of transformants was confirmed by PCR and Southern hybridization. Results showed that 88 % of the mutants contained a single T-DNA insertion. Two of the mutants were observed with different triterpenoid profiles compared with the untransformed cultures. Our results suggest a new functional genomics approach to tag the triterpenoid biosynthesis genes in A. cinnamomea.     

Botryorhodines A-D, antifungal and cytotoxic depsidones from Botryosphaeria rhodina, an endophyte of the medicinal plant Bidens pilosa

An endophytic fungus (Botryosphaeria rhodina) was isolated from the stems of the medicinal plant Bidens pilosa (Asteraceae) that is known for its anti-inflammatory, antiseptic and antifungal effects. The ethyl acetate extract of the fungal isolate exhibits significant antifungal activity as well as potent cytotoxic and antiproliferative effects against several cancer cell lines. Activity-guided fractionation resulted in the isolation of a complex of four depsidones, botryorhodines A-D and the auxin indole carboxylic acid. Botryorhodine A and B show moderate to weak cytotoxic activities against HeLa cell lines with a CC₅₀ of 96.97 μM and 36.41 μM, respectively. In addition, they also show antifungal activity against a range of pathogenic fungi such as Aspergillus terreus (MIC 26.03 μM for botryorhodine A and 49.70 μM for B) and the plant pathogen Fusarium oxysporum (MIC 191.60 μM for botryorhodine A and 238.80 μM for B). A potential role of the endophyte in modulating fungal populations living within or attacking the host plant is discussed.     

A simple bacterial transformation method using magnesium- and calcium-aminoclays

An efficient and user-friendly bacterial transformation method by simple spreading cells with aminoclays was demonstrated. Compared to the reported transformation approaches using DNA adsorption or wrapping onto (in)organic fibers, the spontaneously generated clay-coated DNA suprastructures by mixing DNA with aminoclay resulted in transformants in both Gram-negative (Escherichia coli) and Gram-positive cells (Streptococcus mutans). Notably, the wild type S. mutans showed comparable transformation efficiency to that of the E. coli host for recombinant DNA cloning. This is a potentially promising result because other trials such as heat-shock, electroporation, and treatment with sepiolite for introducing DNA into the wild type S. mutans failed. Under defined conditions, the transformation efficiency of E. coli XL1-Blue and S. mutans exhibited ~ 2 × 10⁵ and ~ 6 × 10³ CFU/μg of plasmid DNA using magnesium-aminoclay. In contrast, transformation efficiency was higher in S. mutans than that in E. coli XL1-Blue for calcium-aminoclay. It was also confirmed that each plasmid transformed into E. coli and S. mutans was stably maintained and that they expressed the inserted gene encoding the green fluorescent protein during prolonged growth of up to 80 generations.     

DNA barcoding to assess species adulteration in raw drug trade of “Bala” (genus: Sida L.) herbal products in South India

The genus Sida L. (family: Malvaceae) is widely used in India and many other countries including China, South East Asia, Africa and South America for treating various neurological disorders and for improving general health and vigour. However, as with many other herbal medicines, it is believed that the Sida products sold in the market may be adulterated with other related or unrelated plant species. In this study, we investigate species adulteration in the raw herbal trade of Sida natural health products (NHPs) in southern India. DNA barcoding was used as a tool to identify the ingredients in the NHPs. A biological reference material (BRM) library for Sida and closely related species was developed using taxonomically authenticated species. DNA barcodes for the species were developed using one nuclear (ITS) and two chloroplast regions (matK and psbA-trnH). The psbA-trnH and ITS region were found to effectively discriminate all species with an interspecific distance of 0.133 and 0.149 and intra-species distance of 0.007 and 0.015 respectively. These DNA barcodes were used to identify the ingredients in raw Sida herbal products obtained from 10 markets in Southern India. Our study indicated that species adulteration in the market samples is rampant especially in case of Sida cordifolia, where all the market samples analyzed were Sida acuta. We discuss the results and the need for a robust herbal drug authentication system to regulate the quality in raw herbal trade market.     

Modulation of ethanol toxicity by Asian ginseng (Panax ginseng) in Japanese ricefish (Oryzias latipes) embryogenesis

Alcohol consumption by women during pregnancy often induces fetal alcohol spectrum disorder (FASD) in children who have serious central nervous system (CNS), cardiovascular, and craniofacial defects. Prevention of FASD, other than women abstaining from alcohol drinking during pregnancy, is not known. A limitation of the use of synthetic anti-alcoholic drugs during pregnancy led us to investigate herbal products. In particular, many plants including Asian ginseng (Panax ginseng) have therapeutic potential for the treatment of alcoholism. We used Japanese ricefish (medaka) (Oryzias latipes), an animal model of FASD, for identifying herbal medicines that can attenuate ethanol toxicity. Fertilized eggs in standard laboratory conditions were exposed to ginseng (PG) root extract (0–2 mg/mL) either 0–2 (group A) or 1–3 (group B) day post fertilization (dpf) followed by maintenance in a clean hatching solution. The calculated IC₅₀ as determined 10 dpf in A and B groups were 355.3 ± 1.12 and 679.7 ± 1.6 μg/mL, respectively. Simultaneous exposure of embryos in sub-lethal concentrations of PG (50–200 μg/mL) and ethanol (300 mM) for 48 h disrupted vessel circulation and enhanced mortality. However, PG (100 μg/mL) may partially protect trabecular cartilage (TC) deformities in the neurocranium in B group embryos induced by ethanol (300 mM). To understand the mechanism, embryonic ethanol concentration was measured at 2 dpf and adh5, adh8, aldh2, aldh9a, catalase, GST, and GR mRNAs were analyzed at 6 dpf. It was observed that although ethanol is able to reduce adh8 and GST mRNA contents, the simultaneous addition of PG was unable to alter ethanol level as well as mRNA contents in these embryos. Therefore, antagonistic effects of PG on ethanol toxicity are mediated by a mechanism which is different from those regulating ethanol metabolism and oxidative stress.     

Non-target impacts of soybean rust fungicides on the fungal entomopathogens of soybean aphid

Soybean aphid, Aphis glycines, has caused serious economic damage to soybean across the North Central US since its introduction to North America in 2000. The management of another invasive soybean pest, Asian soybean rust, Phakopsora pachyrhizi, using foliar fungicide applications has the potential to impact soybean aphid populations by suppressing beneficial fungal entomopathogens. In 2005 and 2006, we applied recommended soybean rust fungicide treatments, consisting of strobilurin and triazole fungicides, to small soybean plots in two locations to assess if such applications might suppress aphid fungal epizootics. In Lamberton, MN, in 2005, during the epizootic, fungicide-treated plots averaged 2.0 ± 0.7% (mean ± SE) disease prevalence while untreated plots averaged 14.2 ± 5.6%. In 2007, we applied strobilurin and strobilurin-triazole mix fungicides to single-plant microplots either before or after release of Pandora neoaphidis, the most commonly observed aphid pathogen in 2005 and 2006. Treatments that contained a mixture of two active ingredients significantly lowered peak and cumulative aphid disease prevalence in both early and late reproductive stage soybeans indicating that fungicide mixtures used to manage soybean rust can negatively impact an aphid-specific fungal pathogen. However, no consistent soybean aphid population response was observed in these studies of low levels of aphid fungal infection.     

The presence of exotic Hypnea flexicaulis (Rhodophyta) in the Mediterranean Sea as indicated by morphology, rbcL and cox1 analyses

Here we report the first finding of Hypnea flexicaulis Yamagishi and Masuda in the Mediterranean Sea (Lagoon of Venice, Italy), identified through molecular analyses using the plastid ribulose-1,5-bisphosphate carboxylase/oxygenase (rbcL) and the mitochondrial protein-coding cytochrome c oxidase subunit I (cox1) genes. The phylogenetic reconstruction, based on rbcL + cox1 multiple alignment, showed that all specimens of H. flexicaulis from Venice, Korea, Philippines and Taiwan were included in a monophyletic group supported by a bootstrap value of 100%.It is highly probable that H. flexicaulis has been introduced from Indo-Pacific populations, in particular the Korean one, probably via ship traffic or shellfish transfers.The use of DNA barcoding combined with morphological observations was, in this case, a rapid way to identify this allochthonous species.     

A duplex-PCR method for species- and pathovar-level identification and detection of the quarantine plant pathogen Xanthomonas arboricola pv. pruni

A PCR-based method was developed for the stone fruit quarantine pathogen Xanthomonas arboricola pv. pruni (Xap), which provides rapid, sensitive and specific in planta detection and isolate identification. Primers specific for Xap were identified using random amplified polymorphic DNA (RAPD). Simplex PCR with these primers had a limit of detection per PCR reaction of approximately 10 CFU for isolate cultures and 50 CFU for plant material when used on tenfold dilutions of isolate culture or genomic DNA extracted from spiked samples, respectively. The primers were adapted as a high-throughput single-step screening based on a digoxigenin-labeled DNA probe assay with a detection limit of 4 × 10² CFU from isolate cultures. A duplex-PCR method was designed that includes the pathovar-level with species-level primers based on species-specific regions of the quinate metabolic gene qumA, increasing diagnostic confidence and offering the first molecular test for all X. arboricola pathovars.     

Efficient shoots regeneration and genetic transformation of Bacopa monniera

Bacopa monniera is an important source of metabolites with pharmaceutical value. It has been regarded as a valuable medicinal plant and its entire commercial requirement is met from wild natural population. Recently, metabolic engineering has emerged as an important solution for sustained supply of assured and quality raw material for the production of active principles. Present report describes efficient in vitro multiplication and transformation method for genetic manipulation of this species. MS medium supplemented with 2 mgl⁻¹ BA and 0.2 mgl⁻¹ IAA was found optimum for maximum shoot regeneration (98.33 %) from in vitro leaves with 2–3 longitudinal cuts. Agrobacterium tumefaciens-mediated transformation method was used for generating transgenic B. monniera plants. Putative transformants were confirmed by GUS assay and PCR based confirmation of hptII gene. DNA blot analysis showed single copy insertion of transgene cassette. An average of 87.5 % of the regenerated shoots were found PCR positive for hptII gene and GUS activity was detected in leaves of transgenic shoots at a frequency of 82.5 % The efficient multiple shoots regeneration system described herein may help in mass production of B. monniera plant. Also, the high frequency transformation protocol described here can be used for genetic engineering of B. monniera for enhancement of its pharmaceutically important metabolites.     

Nutritional value and intake of aquatic ferns (Azolla filiculoides Lam. and Salvinia molesta Mitchell.) in sows

Aquatic ferns (AFs) such as Azolla filiculoides and Salvinia molesta are grown on swine lagoons in the tropics and used in diets for pigs. The present work is aimed at evaluating their potential as feed ingredients for sows. When presented with ad libitum AFs, gilts weighing 110 ± 14 kg (mean ± SD), were able to ingest 9.1–9.7 kg fresh AF per day (from 597 to 630 g dry matter (DM) per day) and from 1240 to 1428 g DM per day when presented in a dry, ground form. A digestibility study was conducted, using sows weighing 213 ± 9 kg (mean ± SD), which were fed diets containing maize, soybean meal and 0, 150 or 300 g AF kg⁻¹ diet. The presence of AFs had a negative impact on the faecal digestibility of the crude protein, NDF and energy content of the whole diet (P<0.001) and on the ileal protein digestibility, especially with 300 g AFs kg⁻¹ diet. The level of AFs in the diet had no effect on stomach weight (P>0.05) but increased the weight of the rest of the gastrointestinal tract (P<0.001). The rate of AF fibre fermentation in the pig large intestine was measured using an in vitro gas test. The rates were much lower than tropical tree foliage, which can also be used in pig diets in the tropics. This could partly explain the low apparent digestibility of AFs in pigs. In conclusion, the inclusion level of AFs in rations for sows should be limited to 150 g AFs kg⁻¹ diet due to the low digestibility and energy density, as well as the negative impact on the digestibility of the whole diet.     

Using a starch-rich mutant of Arabidopsis thaliana as feedstock for fermentative hydrogen production

A mutant plant (Arabidopsis thaliana), sex1-1 (starch excess 1-1), accumulating high starch content in leaves was created to serve as better biomass feedstock for a H₂-producing strain Clostridium butyricum CGS2, which efficiently utilizes starch for H₂ production but cannot assimilate cellulosic materials. The starch content of the mutant plant increased to 10.67 mg/fresh weight, which is four times higher than that of wild type plant. Using sex1-1 mutant plant as feedstock, C. butyricum CGS2 could produce 490.4 ml/l of H₂ with a H₂ production rate of 32.9 ml/h/l. The H₂ production performance appeared to increase with the increase in the concentration of mutant plant from 2.5 to 10 g/l. The highest H₂ to plant biomass yield was nearly 49 ml/g for the mutant plant. This study successfully demonstrated the feasibility of using a starch-rich mutant plant for more effective bioH₂ production with C. butyricum CGS2.     

Amounts of nuclear DNA in marine halophytes

The amount of nuclear DNA, expressed as the C-value, was estimated for 13 marine halophytic plant species from six families. Plant material was collected in the nature reserve of the Strunjan saltpan in the Northern Adriatic and comprised all halophytic species inside the investigated area. Reproductive region of the shoot or root tips of halophytes were dissected, nuclei were Feulgen stained and 2C-values were measured by DNA image cytometry as follows: Crithmum maritimum (4.38 pg DNA), Artemisia caerulescens (6.43 pg), Aster tripolium (21.43 pg), Inula crithmoides (3.63 pg), Atriplex portulacoides (1.83 pg), A. prostrata (1.51 pg), Salicornia europaea (2.75 pg), Salsola soda (2.62 pg), Sarcocornia fruticosa (5.91 pg), Suaeda maritima (2.11 pg), Limonium angustifolium (5.06 pg), Puccinellia palustris (8.15 pg) and Ruppia cirrhosa (4.65 pg). With the exception of the C-value estimate for A. caerulescens, which has been listed in the Plant DNA C-values Database, the C-values represent the first estimates for all the examined species. In addition, the C-value for R. cirrhosa is also the first report for the family Ruppiaceae. The investigated halophytes had a smaller genome size compared to other known C-values for species within a particular family and also when compared to the mean values of dicots and monocots. The study also showed that halophylic annuals have a smaller genome size (2.49 pg) than perennial ones (7.45 pg DNA).     

Sorgoleone,a sorghum root exudate: Algicidal activity and acute toxicity to the ricefish Oryzias latipes

The discovery of natural and natural-based compounds has resulted in its application as an alternative to synthetic algicides to control harmful algae in aquatic systems. Of the many natural-product-based algicides, sorgoleone, a natural plant product from Sorghum bicolor root exudates has been investigated for its controlling effect on different algal species and its acute fish toxicity. Growth of the blue green algal species Microcystis aeruginosa Kützing was completely inhibited by the crude methanol extract of sorghum root at 20 μg mL⁻¹. The most noticeable inhibition was observed in the bioassay of n-hexane soluble extract, where 98% growth inhibition occurred in M. aeruginosa at the concentration of 1.25 μg mL⁻¹. Sorgoleone very effectively controlled blue green algae inhibiting 97% of M. aeruginosa at 0.5 μg mL⁻¹ and 99% of Anabaena affinis Lemmermann at 4 μg mL⁻¹. In contrast, inhibition of the green algae species Chlorella vulgaris Beijerinck and Scenedensmus spp. at 16 μg mL⁻¹ sorgoleone was 87 and 68%, respectively. There were no mortalities or adverse effects observed in any of the fish exposed to water control, solvent control, and a nominal concentration of 1 μg mL⁻¹ during the test period. The no observed effect concentration (NOEC) value was 1.5 μg mL⁻¹ for the tested fish (O. latipes). Sorgoleone can be considered as an effective and an ecologically and environmentally sustainable approach to controlling harmful algae.     

Solitary bees – Potential vectors for gene flow from cultivated to wild strawberries

The genus Fragaria (Rosaceae) contains 24 plant species, including hybrid species such as the octoploid garden strawberry (F. × ananassa). Natural hybridization between Fragaria species has repeatedly been reported, and the potential future cultivation of genetically modified strawberries has made the study of hybridization potential between F. × ananassa and its wild relatives increasingly important. In Europe, F. × ananassa is the only octoploid species present, and the most likely candidate for hybridization is the common diploid woodland strawberry (F. vesca). To date, it is unknown whether pollinator spectra of the two Fragaria species overlap and thus might promote interspecific gene flow. We carried out a survey of flower visitors in northwestern Switzerland to identify major flower visitors of F. vesca and F. × ananassa. This survey indicated that wild bees are the most important shared flower visitors of F. × ananassa and F. vesca. Therefore, we studied flower choice behavior of the common wild bee Osmia bicornis in a greenhouse experiment. Osmia bicornis did not discriminate between F. × ananassa and F. vesca flowers. We conclude that wild bees are important shared flower visitors of both F. × ananassa and F. vesca and are potential vectors for gene flow between cultivated and wild strawberries.     

Relaxant effects of Schisandra chinensis and its major lignans on agonists-induced contraction in guinea pig ileum

In this study, the herbal extracts of Schisandra chinensis were demonstrated to inhibit the contractions induced by acetylcholine (ACh) and serotonin (5-HT) in guinea pig ileum, and the 95% ethanol extract was more effective than the aqueous extract. Analysis with High Performance Liquid Chromatography (HPLC) indicated that schisandrin, schisandrol B, schisandrin A and schisandrin B were the major lignans of Schisandra chinensis, and the ethanol extract contained higher amount of these lignans than the aqueous extract. All four lignans inhibited the contractile responses to ACh, with EC₂₀ values ranging from 2.2 ± 0.4 μM (schisandrin A) to 13.2 ± 4.7 μM (schisandrin). The effectiveness of these compounds in relaxing the 5-HT-induced contraction was observed with a similar magnitude. Receptor binding assay indicated that Schisandra lignans did not show significant antagonistic effect on muscarinic M3 receptor. In Ca²⁺-free preparations primed with ACh or KCl, schisandrin A (50 μM) attenuated the contractile responses to cumulative addition of CaCl₂ by 37%. In addition, schisandrin A also concentration-dependently inhibited ACh-induced contractions in Ca²⁺-free buffer. This study demonstrates that Schisandra chinensis exhibited relaxant effects on agonist-induced contraction in guinea pig ileum, with schisandrin, schisandrol B, schisandrin A and schisandrin B being the major active ingredients. The antispasmodic action of schisandrin A involved inhibitions on both Ca²⁺ influx through L-type Ca²⁺ channels and intracellular Ca²⁺ mobilization, rather than specific antagonism of cholinergic muscarinic receptors.     

Nuclear and chloroplast DNA sequences data support the origin of Potamogeton intortusifolius J.B. He in China as a hybrid between P. perfoliatus Linn. and P. wrightii Morong

Internal transcribed spacers (ITS) of nuclear ribosomal DNA and chloroplast rbcL gene sequence data were used to test the hypothesis that natural populations of Potamogeton intortusifolius J.B. He in China originated from hybridization between P. perfoliatus Linn. and P. wrightii Morong. Based on ITS sequences data, P. intortusifolius possessed heterozygous rDNA genotypes which confirmed the hybrid origin of P. intortusifolius. Chloroplast rbcL gene sequences of P. intortusifolius from Yichang population revealed the same chloroplast haplotype as P. perfoliatus and the samples of P. intortusifolius from Weinan population had the same chloroplast haplotype as P. wrightii, which indicated that both putative parental species had been the maternal parent and that the two populations of P. intortusifolius had independent origins. This study confirms P. intortusifolius as a reciprocal hybrid. Because P. × intortusifolius in China has the same hybrid origin as P. × anguillanus Koidz. in Japan, it is suggested that P. × intortusifolius should be a synonym of P. × anguillanus.     

Association between osteoporosis and polymorphisms of the bone Gla protein, estrogen receptor 1, collagen 1-A1 and calcitonin receptor genes in Turkish postmenopausal women

In this study, we have investigated the association between osteoporosis and osteocalcin (BGLAP) − 298 C>T, estrogen receptor 1 (ER1) 397 T>C, collagen type1 alpha 1 (Col1A1) 2046 G>T and calcitonin receptor (CALCR) 1340 T>C polymorphisms. Genomic DNA was obtained from 266 persons (158 osteoporotic and 108 healthy controls). Genomic DNA was extracted from EDTA-preserved peripheral venous blood of patients and controls by a salting-out method and analyzed by PCR-RFLP. As a result, there was no statistically significant difference in the genotype and allele frequencies of patients and controls for BGLAP − 298 C>T, Col1A1 2046 G>T, ER1 397 T>C and CALCR 1340 T>C polymorphisms. However, ER1 CC genotype compared with TT + TC genotypes was found to increase the two fold the risk of osteoporosis [p = 0.039, OR = 2.156, 95% CI (1.083–4.293)] and CALCR CC genotype compared with TT + TC genotypes was found to have protective effect against osteoporosis [p = 0.045, OR = 0.471, 95% CI (0.237–0.9372)]. In the combined genotype analysis, ER1/CALCR TCCC combined genotype was estimated to have protective effect against osteoporosis [p = 0.0125, OR = 0.323, 95% CI (0.1383–0.755)] whereas BGLAP/Col1A1 CCTT and ER1/CALCR CCTT combined genotypes were estimated as risk factors for osteoporosis in Turkish population (p = 0.027, p = 0.009 respectively).     

Comparative growth kinetics and virulence of four different isolates of entomopathogenic fungi in the house fly (Muscadomestica L.)

Virulence (speed of kill) of a fungal entomopathogen against a particular host insect depends on biological properties of the specific isolate-host combination, together with factors such as fungal dose. How these intrinsic and extrinsic factors affect the actual pattern and extent of fungal growth invivo is poorly understood. In this study we exposed adult house flies (Muscadomestica L.) to surfaces treated with high and low doses of Beauveriabassiana (isolates BbGHA and Bb5344), Metarhiziumanisopliae (strain MaF52) and M.anisopliae var. acridum (isolate Ma189) and used quantitative real-time PCR with species-specific primers to examine the relationship between fungal growth kinetics and virulence. At the highest dose, all fungal isolates killed flies significantly faster than controls, with BbGHA, Bb5344 and MaF52 roughly equivalent in virulence (median survival time (±SE) = 5.0 ± 0.10, 5.0 ± 0.08 and 5.0 ± 0.12 days, respectively) and Ma189 killing more slowly (MST = 8.0 ± 0.20 days). At the lower dose, effective virulence was reduced and only flies exposed to isolates BbGHA and Bb5344 died significantly faster than controls (MST = 12 ± 1.36, 15 ± 0.64, 18 ± 0.86 and 21.0 ± 0.0 days for BbGHA, Bb5344, MaF52 and Ma189, respectively). Real-time PCR assays revealed that flies exposed to surfaces treated with the high dose of spores had greater spore pickup than flies exposed to the low dose for each isolate. After pickup, a general pattern emerged for all isolates in which there was a significant reduction of recovered fungal DNA 48 h after exposure followed by a brief recovery phase, a stable period of little net change in fungal sequence counts, and then a dramatic increase in sequence counts of up to three orders of magnitude around the time of host death. However, while the patterns of growth were similar, there were quantitative differences such that higher final sequence counts were recovered in insects infected with the most lethal isolates and with the higher dose. These results suggest that variation in virulence between isolates, species and doses is determined more by quantitative rather than qualitative differences in fungal growth kinetics.