Open Access
How to translate text using browser tools
Allen L. Szalanski, James W. Austin, J. Mckern, Matthew T. Messenger
Author Affiliations +

Genetic evidence for a new subterranean termite species herein named Reticulitermes okanaganensis is provided based on DNA sequence analysis. Partial sequences of the mitochondrial DNA rRNA 16S gene were obtained from 27 samples of R. okanaganensis from British Columbia, Idaho, Oregon, Nevada and California. Five nucleotide sites were variable among the four observed haplotypes. One haplotype occurred only once, while the most common haplotype, O3, occurred in 37% of the samples. Molecular phylogenetic analysis of R. okanaganensis relative to five other North American Reticulitermes species has clarified its distinct position within the genus.

There is a general consensus that the genus Reticulitermes is in need of revision (Weesner 1970; Nutting 1990; Scheffrahn & Su 1994; Forschler & Jenkins 1999). This is an especially difficult because of the lack of discrete morphological characters which provide for accurate identification of specimens to the specifies level. For this reason, the application of non-morphological identification methods such as the evaluation of cuticular hydrocarbon analysis and mtDNA markers have been used.

Most recently, the application of the mitochondrial rRNA 16S gene has been applied to identify Reticulitermes populations from the south central United States (Austin et al. 2004a, 2004b, 2004c) for species across North America (Austin et al. 2005b), and for clarification of exotic introductions of nearctic Reticulitermes around the world (Austin et al. 2005a; Su et al. 2006). The application of this marker has tremendous potential for molecular diagnostics of Reticulitermes, with increased accuracy of positive species identifications (Szalanski et al. 2003), and clarifying the identities of exotic introductions around the world (Austin et al. 2005a) and from North America (Austin et al. 2005b). Recently, Copren et al. (2005) found evidence for as many as seven new species of Reticulitermes from the western United States based upon cuticular hydrocarbon phenotypes, but resolved to designate them as putative new species after attempting to corroborate their relationship with molecular phylogenetics and reproductive flight dates. These discrepancies likely are attributed to the environmental plasticity of cuticular hydrocarbons and stresses the need for fixed character states for species identification such as mtDNA sequences.

We provide the first genetic evidence of a new species, named herein as R. okanaganensis after its first collection site in the Okanagan region of British Columbia, Canada, and provide a phylogenetic analysis of Reticulitermes by applying the 16S mtDNA gene. In addition, the geographic distribution of Reticulitermes species and haplotypes of species throughout the region are discussed.

Materials and Methods

Termites were collected from various locations in British Columbia, Idaho, Oregon, Nevada, and California, both from our own collecting efforts and from the 2002 National Termite Survey (Table 1). Samples were preserved in 100% ethanol.

Alcohol-preserved specimens were allowed to dry on filter paper, and DNA was extracted according to Liu & Beckenbach (1992) on individual whole worker termites with the Puregene DNA isolation kit D-5000A (Gentra, Minneapolis, MN). Extracted DNA was resuspended in 50 μl of Tris:EDTA and stored at -20°C. Polymerase chain reaction was conducted with the primers LR-J-13007 (5’-TTACGCTGTTATCCCTAA-3’) (Kambhampati & Smith 1995) and LR-N-13398 (5’-CGCCTGTTTATCAAAAACAT-3’) (Simon et al. 1994). These PCR primers amplify an approximately 428-bp-region of the mtDNA 16S rRNA gene. The PCR reactions were conducted with 2 μl of the extracted DNA (Szalanski et al. 2000), with a profile consisting of 35 cycles of 94°C for 45 s, 46°C for 45 s, and 72°C for 60 s. Amplified DNA from individual termites was purified and concentrated with minicolumns (Wizard PCRpreps, Promega) according to the manufacturer’s instructions. Samples were sent to The University of Arkansas Medical Center DNA Sequencing Facility (Little Rock, AR) for direct sequencing in both directions. GenBank accession numbers were DQ389178-DQ389180 and DQ438936 for the Reticulitermes okanaganensis haplotypes found in this study. Consensus sequences for each sample were obtained with Bioedit 5.09 (Hall 1999). Mitochondrial DNA haplotypes were aligned with MacClade v4 (Sinauer Associates, Sunderland, MA).

The distance matrix option of PAUP* 4.0b10 (Swofford 2001) was used to calculate genetic distances according to the Kimura 2-parameter model of sequence evolution (Kimura 1980). Mitochondrial 16S sequences from R. flavipes, R. hesperus, R. tibialis, R. hageni, R. flavipes, and R. virginicus (Szalanski et al. 2003; Austin et al. 2004a,b,c) were added to the DNA sequence dataset for comparison. DNA sequences from the Formosan termite, Coptotermes formosanus Shiraki (GenBank AY558910), and Heterotermes aureus (Snyder) (GenBank AY280399), were added to act as outgroup taxa. DNA sequences were aligned by CLUSTAL W (Thompson et al. 1994). Maximum likelihood and unweighted parsimony analysis on the alignments were conducted by PAUP* 4.0b10 (Swofford 2001). Gaps were treated as missing characters for all analysis. The reliability of trees was tested with a bootstrap test (Felsenstein 1985). Parsimony bootstrap analysis included 1,000 resamplings with the Branch and Bound algorithm of PAUP*. For maximum likelihood analysis, the default likelihood parameters were used (HKY85 six-parameter model of nucleotide substitution, empirical base frequencies with the exception of the transition/transversion ratio, which was set to 2.596954:1). These parameters were used to carry out a heuristic search by PAUP* with a neighbor joining tree as the starting tree.


DNA sequencing of the 16S rRNA amplicon from R. okanaganensis revealed an average size of 428 bp. The average base frequencies were A = 0.41, C = 0.23, G = 0.13, and T = 0.23. Among the 27 R. okanaganensis mtDNA rRNA 16S DNA sequences, a total of 5 nucleotide sites were variable (Table 2). Four distinct haplotypes (lineages) were observed (Table 1), and genetic divergence among these haplotypes ranged from 0.2 to 0.5 percent. One haplotype, O4, occurred only once, while the most common haplotype, O3, occurred in 37% of the R. okanaganensis samples. While haplotype O4 was only found in Nevada, haplotype O1 was found over the largest geographical area (Table 1).

We conducted a phylogenetic analysis on R. okanaganensis relative to all described North American Reticulitermes species to clarify the phylogenetic relationships of R. okanaganensis within the genus. Parsimony analysis of the aligned Reticulitermes spp. and the outgroup taxa used 436 characters, of which 90 were variable (20%) and 59 (14%) were parsimony informative. This analysis had a single consensus tree with a length = 154, and a CI value of 0.695 and verified the distinct monophyly of ‘R. okanaganensis’ in the genus (Fig. 1). Maximum likelihood analysis (-ln L = 1363.74728) also supported a distinct clade for R. okanaganensis and was identical to the maximum parsimony analysis. Reticulitermes okanaganensis formed a sister group to R. flavipes, while R. hesperus grouped with R. tibialis, and R. hageni formed a distinct clade with R. virginicus.


Phylogenetic analysis of North American Reticulitermes revealed a distinct clade for R. okanaganensis. There is clear genetic evidence for R. okanaganensis being identified as a discrete species, possessing multiple haplotypes over a broad geographical range including California, Idaho, Oregon, Nevada, and British Columbia (Fig. 2). The genetic isolation observed demonstrates that this group represents a distinct taxonomic entity, albeit cryptic due to overlapping morphological characters with other congeners from the region that merits further investigation.

The identification of new species from northern California has been reported by Haverty & Nelson (1997), Haverty et al. (1999) applying cuticular hydrocarbons, and further investigated with ethological data (Getty et al. 2000a, 2000b). It has been proposed that extensive collecting from western states may produce as many as 6 new Reticulitermes species, and perhaps at least three Reticulitermes species in Mexico (Myles 2000). We have already identified other morphologically and genetically distinct nearctic Reticulitermes from this region and we are attempting to collect additional material in the near future so that subsequent investigations can be performed and the overall relationships of western nearctic Reticulitermes can be better understood.

Reticulitermes okanaganensis has become a problem in British Columbia, where there is an unprecedented number of attacks to structures, particularly around Kelowna, BC (Hugh Philip unpublished). Reports of difficulties controlling Reticulitermes infestations from urban structures in northern California (Kistner & Sbragia 2001) prompts us to consider whether these difficult control scenarios are due to the structures themselves, or the possibility of infestations from this newly identified species, that may have slightly different nutritional, ecological, and behavioral requirements.

The western subterranean termite, Reticulitermes hesperus, is the most destructive termite found in California (Lewis 2001), but it is highly probable that R. okanaganensis accounts for a significant portion of damage previously thought to be associated with R. hesperus. Snyder (1949) originally described R. hesperus based on material from Little Bear Lake, San Bernadino mountains, California and samples previously collected from Vancouver, British Columbia, and suggests the range to extend from British Columbia through Washington, Oregon, California, and Nevada, extending down to Baja, California (Spencer 1937). The northern limit of its range occurs in the upper Frazer valley from Lytton to Kamloops, British Columbia (Spencer 1937, 1945). Our genetic evaluation of populations from along the Pacific and throughout the western states suggest that R. hesperus does not generally occur away from the pacific coast region, as is described by Weesner (1965), and is isolated west of the Sierra Nevada/Cascades mountain ranges (unpublished). Reticulitermes tibialis overlaps entirely with the range of R. okanaganensis, and extends east as far as Indiana, but is not commonly found east of Texas and Oklahoma (Banks & Snyder 1920; Austin et al. 2004a). Light and Pickens (1934) indicate collections of R. tibialis from the northern tip of Idaho, while the U.S. Dept. of Agriculture (1959a, 1959b) have reports from Twin Falls and Lewiston to Lapwari, Idaho. Weesner (1970) concedes that collections from western California to as far east as Elko and Reno, Nevada appear to be R. tibialis, but are too difficult to assign to either R. hesperus or R. tibialis.

Reticulitermes okanaganensis haplotype O1, which has been reported in Osoyoos, BC, Canada, is also found in El Dorado county California, Washoe county Nevada, and Ada and Nez Perce counties in Idaho. Osoyoos, B.C. is located in the southern interior ecoprovince of Canada, a region characterized as a desert climatic zone. We are currently investigating the distribution of Reticulitermes species throughout the western United States. From this study and ongoing research, we expect that R. okanaganensis will be found in the drier regions north of the Cascades and East of the Sierra Nevada Ranges. While the possible importation of this species either to the United States or to Canada should be considered, in all probability we are dealing with a cryptic endemic species. Establishment of this species in Canada may be attributed to trade between the US and Canada, and may account for the unprecedented number of structures being attacked in British Columbia. It may be equally plausible that urban sprawl and encroachment into previously undeveloped areas for housing needs has contributed to the increase frequency of attacks to structures there.

An important criterion for determining the extent of genetic variation for a species lies in the ability to sample from populations evenly distributed within the species range (Mayr & Ashlock 1991). For this reason, future studies of this unknown species, including a proper species description, are required. This study represents an important first step towards this endeavor.


We thank M. Rust for critical review of this manuscript. Thanks are extended to R. Saran Univ. of CA Riverside, P. Pachamuthu, and S. Vega of Western Exterminating (Anaheim, CA), numerous pest management professionals, and especially to H. Philip of the Canadian Ministry of Agriculture, Food and Fisheries, Food Safety and Quality Branch Plant Health Unit, Kelowna, B. C. Canada, for providing samples. This research was supported in part by the University of Arkansas, Arkansas Agricultural Experiment Station, and the City of New Orleans Mosquito and Termite Control Board.

Literature Cited


J. W. Austin, A. L. Szalanski, and M. T. Messenger . 2004a. Genetic variation and distribution of the subterranean termite genus Reticulitermes (Isoptera: Rhinotermitidae) in Arkansas and Louisiana. Florida Entomol 87:473–480. Google Scholar


J. W. Austin, A. L. Szalanski, R. E. Gold, and B. T. Foster . 2004b. Genetic variation and geographical distribution of the subterranean termite genus Reticulitermes in Texas. Southwest. Entomol 29:1–11. Google Scholar


J. W. Austin, A. L. Szalanski, and B. M. Kard . 2004c. Distribution and genetic variation of Reticulitermes (Isoptera: Rhinotermitidae) in Oklahoma. Florida Entomol 87:145–151. Google Scholar


J. W. Austin, A. L. Szalanski, R. H. Scheffrahn, M. T. Messenger, S. Dronnet, and A-G. Bagnères . 2005a. Genetic evidence for the synonymy of two Reticulitermes Species: Reticulitermes flavipes (Kollar) and Reticulitermes santonensis (Feytaud). Ann. Entomol. Soc. Amer 98:395–401. Google Scholar


J. W. Austin, A. L. Szalanski, R. H. Scheffrahn, and M. T. Messenger . 2005b. Genetic variation of Reticulitermes flavipes (Isoptera: Rhinotermitidae) in North America using the mitochondrial rRNA 16S gene. Ann. Entomol. Soc. Amer. 98: (in press). Google Scholar


N. A. Banks and T. E. Snyder . 1920. A Revision of the Nearctic Termites. U.S. Nat. Mus., Bull. 108: 228 pp. Google Scholar


N. A. Banks 1946. Species Distinction in Reticulitermes (Isoptera: Rhinotermitidae). MS thesis. University of Chicago, Chicago, Illiniois. 29 pp. Google Scholar


K. A. Copren, L. J. Nelson, E. L. Vargo, and M. I. Haverty . 2005. Phylogenetic analyses of mtDNA sequences corroborate taxonomic designations based on cuticular hydrocarbons in subterranean termites. Mol. Phylogen. Evol 35:689–700. Google Scholar


S. E. Donovan, D. T. Jones, W. A. Sands, and P. Eggleton . 2000. Morphological phylogenetics of termites (Isoptera). Biol. J. Linnean Soc 70:467–513. Google Scholar


A. E. Emerson and E. M. Miller . 1943. A key to the termites of Florida. Entomol. News 54:184–187. Google Scholar


J. Felsenstein 1985. Confidence limits on phylogenies: An approach using the bootstrap. Evolution 39:783–791. Google Scholar


B. T. Forschler and T. M. Jenkins . 1999. Evaluation of subterranean termite biology using genetic, chemotaxonomic, and morphometric markers and ecological data: a testimonial for multi-disciplinary efforts. Trends Entomol 2:71–79. Google Scholar


G. M. Getty, M. I. Haverty, K. A. Copren, and V. R. Lewis . 2000a. Response of Reticulitermes spp. (Isoptera: Rhinotermitidae) in northern California to baiting with hexaflumuron with Sentricon Termite Colony Elimination System. J. Econ. Entomol 93:1498–1507. Google Scholar


G. M. Getty, M. I. Haverty, and V. R. Lewis . 2000b. Agonistic behavior between recently collected and laboratory cultured Reticulitermes spp. (Isoptera: Rhinotermitidae) from northern California. Pan-Pacific Entomol 76 243–250. Google Scholar


T. A. Hall 1999. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl. Acids. Symp. Ser 41:95–98. Google Scholar


M. I. Haverty and L. J. Nelson . 1997. Cuticular hydrocarbons of Reticulitermes (Isoptera: Rhinotermitidae) from northern California indicate undescribed species. Comp. Biochem. Physiol B118:868–880. Google Scholar


M. I. Haverty, L. J. Nelson, and B. T. Forschler . 1999. New cuticular hydrocarbon phenotypes of Reticulitermes (Isoptera: Rhinotermitidae) from the United States. Sociobiol 34:1–21. Google Scholar


N. C. Hostettler, D. W. Hall, and R. H. Scheffrahn . 1995. Intracolony morphometric variation and labral shape in Florida Reticulitermes (Isoptera: Rhinotermitidae) soldiers significance for identification. Florida Entomol 78:119–129. Google Scholar


S. Kambhampati and P. T. Smith . 1995. PCR primers for the amplification of four insect mitochondrial gene fragments. Insect Molec. Biol 4:233–236. Google Scholar


M. Kimura 1980. A simple method for estimating evolutionary rate of base substitutions through comparative study of nucleotide sequences. J. Molec. Evol 16:111–120. Google Scholar


D. H. Kistner and R. J. Sbragia . 2001. The use of the Sentricon Termite Colony Elimination System for controlling termites in difficult control sites in northern California. Sociobiol 37:265–280. Google Scholar


V. R. Lewis July 2001. Pest Notes: Termites. UC ANR Publication 7415. 6 pp. Also available online at  www.ipm.ucdavis.eduGoogle Scholar


H. Liu and A. T. Beckenbach . 1992. Evolution of the mitochondrial cytochrome oxidase II gene among 10 orders of insects. Mol. Phylogen. Evol 41:41–52. Google Scholar


S. F. Light 1934. The termite fauna of North America with special reference to the United States. pp. 127-135 In C. A. Kofoid, S. F. Light, A. C. Horner, M. Randall, W. B. Herms, and E. E. Bowe [eds.], Termites and Termite Control, 2nd ed., University of California Press, Berkeley, CA. Google Scholar


S. F. Light and A. L. Pickens . 1934. American subterranean termites and their distribution. pp. 150-156 In C. A. Kofoid [ed.], Termites and Termite Control, Univ. of California Press, Berkeley, California. Google Scholar


E. Mayr and P. D. Ashlock . 1991. Principles of Systematic Zoology, 2nd ed. McGraw-Hill Book Company, New York. Google Scholar


E. M. Miller 1949. A Handbook on Florida Termites. University of Miami Press, Coral Gables, Florida. 30 pp. Google Scholar


T. G. Myles and M. M. Hooten . 2000. Evaluation of the Occurrence of Termite Species and Potential Termite Burrowing Depths in the Area of the Nevada Test Site: Present Day and for the Next 10,000 Years. Neptune and Co., Inc. Document Number 05100-01. 17 pp. Google Scholar


W. L. Nutting 1990. Insecta, Isoptera. pp. 997-1032 In D. L. Dindal, D. L. [ed.], Soil Biology Guide. Wiley and Sons, New York. Google Scholar


A. L. Pickens 1934a. The biology and economic significance of the western subterranean termite, Reticulitermes hesperus. pp. 157-183 In C. A. Kofoid, S. F. Light, A. C. Horner, M. Randall, W. B. Herms, and E. E. Bowe [eds.], Termites and Termite Control, 2nd ed. University of California Press, Berkeley, CA. Google Scholar


A. L. Pickens 1934b. The barren-lands subterranean termite, Reticulitermes tibialis. pp. 184-186 In C. A. Kofoid, S. F. Light, A. C. Horner, M. Randall, W. B. Herms, and E. E. Bowe [eds.], Termites and Termite Control, 2nd ed. University of California Press, Berkeley, CA. Google Scholar


R. H. Scheffrahn and N-Y. Su . 1994. Keys to soldier and winged adult termites (Isoptera) of Florida. Florida Entomol 77:460–474. Google Scholar


C. Simon, F. Frati, A. Beckenbach, B. Crespi, H. Liu, and P. Flook . 1994. Evolution, weighting, and phylogenetic utility of mitochondrial gene sequences and a compilation of conserved polymerase chain reaction primers. Ann. Entomol. Soc. Amer 87:651–701. Google Scholar


T. E. Snyder 1949. Catalog of the Termites (Isoptera) of the World. Smithsonian Misc. Coll 112:1–490. Google Scholar


R. R. Sokal and T. J. Crovello . 1970. The biological species concept: a critical evaluation. Amer. Natur 104:127–153. Google Scholar


G. J. Spencer 1937. The termite situation in British Columbia in 1936. Proc. Entomol. Soc. Brit. Columbia 33:42–43. Google Scholar


SpencerG. J. 1945. On the incidence, density and decline of certain insects in British Columbia. Proc. Entomol. Soc. Brit. Columbia 42:21. Google Scholar


N-Y. Su, W. Ye, R. Ripa, R. H. Scheffrahn, and R. M. Giblin-Davis . 2006. Identification of Chilean Reticulitermes (Isoptera: Rhinotermitidae) Inferred from Three Mitochondrial Gene DNA Sequences and Soldier Morphology. Ann. Entomol. Soc. Am 99:352–363. Google Scholar


D. L. Swofford 2001. PAUP*: Phylogenetic analysis using parsimony (*and other methods), ver. 4.0b10. Sinauer, Sunderland, Massachusetts. Google Scholar


A. L. Szalanski, D. S. Sikes, R. Bischof, and M. Fritz . 2000. Population genetics and phylogenetics of the endangered American burying beetle, Nicrophorus americanus (Coleoptera: Silphidae). Ann. Entomol. Soc. Amer 93:589–594. Google Scholar


A. L. Szalanski, J. W. Austin, and C. B. Owens . 2003. Identification of Reticulitermes spp. (Isoptera: Rhinotermitidae) from the South Central United States by PCR-RFLP. J. Econ. Entomol 96:1514–1519. Google Scholar


J. D. Thompson, D. G. Higgins, and T. J. Gibson . 1994. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignments through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res 22:4673–4680. Google Scholar


U.S. Dept. Agr., Plant Pest Control Division 1959a. Coop. Econ. Insect Rept. 9, 362.  Google Scholar


U.S. Dept. Agr., Plant Pest Control Division 1959b. Coop. Econ. Insect Rept. 9, 393.  Google Scholar


F. M. Weesner 1970. Termites of the Nearctic Region. pp. 477-525 In K. Krishna and F. M. Weesner [eds.], Biology of Termites, vol. II. Academic, New York. Google Scholar


F. M. Weesner 1965. The Termites of the United States: A Handbook. Printed for the National Pest Control Association, Elizabeth, NJ. 71 pp. Google Scholar


Fig. 1.

Maximum parsimony cladogram of Reticulitermes okanaganensis and related taxa.


Fig. 2.

Distribution of Reticulitermes okanaganensis haplotypes in British Columbia, Idaho, Oregon, Nevada, and California.


Table 1.

Collection data, and haplotypes (HAP) for Reticulitermes okanaganensis


Table 2. Base pair differences between Reticulitermes okanaganensis haplotypes.

Allen L. Szalanski, James W. Austin, J. Mckern, and Matthew T. Messenger "GENETIC EVIDENCE FOR A NEW SUBTERRANEAN TERMITE SPECIES (ISOPTERA: RHINOTERMITIDAE) FROM WESTERN UNITED STATES AND CANADA," Florida Entomologist 89(3), 299-304, (1 September 2006).[299:GEFANS]2.0.CO;2
Published: 1 September 2006
16S gene
DNA sequence
Back to Top