Sequence analysis of mitochondrial 16S ribosomal DNA gene fragment from wasps




Ruojun Huang

Introduction

ABSTRACT

Braconidae and Ichneumonidae are two parasitic wasps with different hosts. Distinguishing these two different subfamilies of wasps by morphological features is difficult. This sequence analysis of the mitochondrial 16S rDNA gene has been used for molecular taxonomy in many insects. In this paper I analyzed a 552 bp region of the mitochondrial 16S rRNA gene in the two subfamilies of wasps. The results showed that molecular technique combined with morphological characters could be useful in identification of subfamilies.

Key words: 16S rDNA, Braconidae, Ichneumonidae


INTRODUCTION

Braconidae and ichneumonidae wasps have been extensively studied due to their importance as biological control agents. Both taxonomic groups are typically small (<15mm), dark in coloration, with two sets of wings and thin pedicel. Ichneumonidae and Braconidae are two subfamily of family Ichneumonidae. Morphological analyses suggested that two subfamilies are closely affiliated (Quicke, 1990). Precise identification of specimens is very difficult - involving very close examination of microscopic features. Molecular methods are expected to provide new and more precise means in distinguishing different subfamilies. The mitochondrial 16S ribosomal DNA gene was chosen because it has been used in many insect phylogenetic studies, at both lower and higher taxonomic levels, as well as in studies of diverse Hymenopteran groups (Simon, 1994; Dowton, 1994; Dowton, 1998).


MATERIALS AND METHODS

Wasp discrimination by morphological character
Wasps were kindly offered by Kevin Cox. They were preserved in 70% ethanol. Microscopic analysis of wing veination was used to morphologically distinguish Braconidae and Ichneumonidae.

DNA extraction
Abdomens were removed from Ichneumonidae samples and used for extraction while the abdomen and thorax were used for Braconidae samples because Ichneumonidae is bigger than Braconidae. The head region was excluded as it frequently contained PCR inhibitors. (Dowton, 1998). With each sample I, respectively: (a) removed ethanol by washing three times (15 min each) in 10 mM Tris-Hcl (pH 8) containing 100 mM mM NaCl and 1mM MgCl2, (b) homogenized tissue using a disposable minipestle, (c) added 25 ul of 50 mM DTT (Dithiothreitol), (d) added 25 ul of 0.5 M KOH, (e) spinned for 5 s at 13, 000 rpm, (f) incubated in a heat block for 5 min at 65C, (g) vortexed for 10 s, (h) incubated for 5 min at 65C, (i) spin for 10 s at 13,000 rpm, (j) add 25 ul of 0.5 M HCl, (k) add 8.4 ul of 0.5 M Tris-HCl pH. 9.0, (l) spin for 10 s at 13,000 rpm, (m) mix by hand agitation, (n) spin for 1 min at 13,000 rpm, (o) transfer 50 ul of the supernatant into a new tube (Starks, 2002).

Primers
16S rDNA primers were synthesized in an oligonucleotide synthesizer in our laboratory. The forward primer is 5’-CACCTGTTTATCAAAAACAT-3’ and reverse primer 5’-CTTATTCAACATCGAGGTC-3’ (Crozier, 1993). These primers were designed to generate a 552 bp DNA fragment.

PCR
PCR was performed in a 25 ul reaction containing 2.5 ml 10x PCR buffer, 2.0 ml 2.5mM dNTPs, 0.75 ml 50 mM MgCl2, 0.5 ml 10 pmol/ul forward primer, 0.5 ml 10 pmol/ul reverse primer, 0.175 ml Taq DNA polymerase (5 U/ ul), 87.5~175.0 ng DNA template. Tubes were placed in a Perkins-Elmer Thermal Cycler. Cycling conditions were as follows: (a) denaturation at 94C for 30 sec, (b) 50 cycles of 1min denaturation at 94C, 1 min annealing at 55C, 1 min extension at 72C, and (c) 8 min extension at 72C. (Patrick, 1999). Products were stained with ethidium bromide and visualized under UV light.

PCR purification
The PCR product was purified using Rapid PCR Purification Systems (Marligen Bioscience). For each 25 ml PCR reaction, I (a) added 400 ml of Binding solution (H1) and mix, (b) placed a cartridge in a 2ml collection tube and loaded the sample from step (a) into the cartridge, (c) centrifuged at 13000*g for 60 sec, (d) discarded flow-through and place column back into the same tube, (e) added 700 ul Wash buffer (H2) to column and centrifuge at 13000*g for 2 min, (f) discarded flow-through, placed column into tube, and centrifuge column for an additional 2 min at 13000*g speed, (g) placed column in a clean 1.5 ml microfuge tube, (h) added 50 ul dH2O to the center of column, allowed for 1 min, and then centrifuged at 13000g for 2 min.

Sequencing
Sequencing by Sanger dideoxy method was performed using AmpliCycle Sequencing kit( Perkin Elmer, Norwalk, CT). For each reaction, 4 ul of 10X buffer, 2 ul of reverse primer, 50 fmol of DNA, and 0.25 ul aP33 (dATP) was used. 6 ul of this mixture was added to each of four 200 ul Eppendorf tubes, each containing either 2 ul of ddGTP, ddATP, ddTTP, or ddCTP. Each tube was overlaid with mineral oil and placed in a Perkin-Elmer Thermal Cycler for the reaction. Cycling conditions were: 94C for 2 min, 35 cycles of 94C for 30 sec, 58C for 30 sec, 72C for 1 min. Following the sequencing reactions, 4 ul of stop solution was added to each tube and samples were boiled for 3 min prior to loading and electrophoresis on a sequencing gel.

Sequence alignment
Sequences were read into a computer from radiographs and aligned using MacVector 6.5.3.

Figures


Figure 1-Figure 1: Microscopic pictures of four specimens


Figure 2-Figure 2: PCR product of DNA extraction using 100bp ladder. Lanes 1: Ichneumonidae - 1, lane 2: Ichneumonidae –2, lane 3: Ichneumonidae –3, lane 4: Braconidae – 1, lane 5: Braconidae – 2, lane 6: Braconidae – 3, lane 7: Yellow jacket (Vespula spp.), and lane 7: negative control.


Figure 3-Fig 4: Sequence conserved percentage compariation of any two samples (%)


Figure 4


RESULTS
Six specimens were examined and based on their morphologies, and 3 wasps were classified as Ichneumonidae and 3 Braconidae (Fig 1). (Here only show the four samples with sequencing results)

Total DNA was extracted from 6 samples and one yellow jacket (Vespula spp.) as a positive control as described in Materials and Methods. PCR was performed as described in Materials Methods. Purified PCR product was then electrophosed on 1.0% agarose gel (Fig. 2). All samples showed sharp, intense bans with same size of about 500 bps.
Sequence analysis of the DNA fragments were performed as described in Material and Methods. The alignment is present in Figure 3 (Alignment of 16S rDNA sequences. ‘N’ indicated where sequence data were ambiguous. Only 5 samples were readable and comparable including Ichneumonidae - sample 1, Ichneumonidae – sample 2, Braconidae – sample 1, Braconidae – sample 2, and Yellow jacket (Vespula spp.). The result showed a sequence about 120 bp. All 5 samples aligned with honeybee mitochondrial 16S mtrDNA from which primers were designed, 43 characters (35.83%) of 120 bps were conserved. This result confirmed that PCR product of DNA fragments is 16S mitochondrial rDNA and is located at nucleotides13924-13943 in genome with 16343 bp.

Based on alignment results, a paragraph was drawn to compare the conserved percentage of any two samples of Ichneumon and Braconidae (Fig 4).

DISCUSSION

The two subfamilies can be distinguished through morphological analysis: size, ovipositors, and wing veins. Ichneumonidae is bigger and female has longer ovipositor than Braconidae. However, in our sample, it was difficult to distinguish them just by size and ovipositor as many lacked ovipositors (individuals were male) and size was highly variable. Wing veins were used to help distinguish the relatedness of the organism. The morphological classification was not completely supported by the gene sequence analysis.

Molecular data from mitochondrial genomes provide another possible method to distinguish wasps. Comparison of homology between any two of Ichneumonidae and Braconidae samples showed that 2 Braconidae samples have highest homology (75.65%), while for the 2 Ichneumonidae samples’ homology is not higher than Ichneumonidae-1 and Brconidae-2’s, Ichneumonidae-2 and Brconidae-1’s, and Ichneumonidae-1 and Brconidae-1’s. It is possible that the two Ichneumonidae samples are not same subfamilies although we thought they were by morphological features. Possible reasons to for the discrepancy include: 1) 120 of 552 bp analyzed may not be sufficient; 2) 552 bp of 16343 bp gene also may not be large enough to distinguish; 3) the analyzed region may be too variable enough to discriminate.

To increase the possibility of discrimination, several things can be considered. First, obtain more sequence information. Secondly, use different regions of gene. Last, by combine this data with other molecular markers such as nuclear (28S) rDNA, mitochondrial NADH1dehydrogenase, and mitochondrial cytochrome oxidase I (COI) gene sequences (Dowton, 1994;Dowton, 1998; Belshaw, R., 1997).

Full Paper

Acknowledgments

Sincere thanks to Jingsong Qiu and Brian Fox for their help and guidance throughout this project. Thanks to Dr. Berish Rubin for giving me this opportunity to undertake this independent piece of work. I am very grateful to Kevin Cox for providing me with wasps, to Dr. Sylvia Anderson for generating primers for me, to Dr. John. D. Wehr for providing me with minipestles.

REFERENCES

Quicke, D. L. J., and van Achterberg, C. (1990). Phylogeny of the subfamilies of the family Braconidae (Hymenoptera: Ichneumonoidea). Zool. Verh. Leiden 258:1-95
Belshaw, R., and Quicke, D.S.J., R. (1997). A molecular phylogeny of the Aphidiinae (Hymenoptera: Braconidae). Mol. Phylogenet. Evol. 7:281-293.
Crozier, R. H. , and Crozier, Y. C. (1993). The mitochondrial genome of the honebee Apis mellifera: Complete sequence and genome organization. Genetics 133: 97-117
Dowton, M., and Austin, A. D. (1994). Molecular phylogeny of the insect order Hymenoptera: Apocritan relationships. Proc. Natl. Acad. Sci. USA 91: 9911-9915
Dowton, M., and Austin, A. D. (1998). Phylogenetic relationships among the Microgastroid wasps (Hymenoptera: Braconidad): combined analysis of 16S and 28S rDNA genes and morphological data. Molecular Phylogenetics and Evolution. 10(3): 354-366.
Dowton, M., Austin, A. D., and Antolin, M. F. (1998). Evolutionary relationships among the Braconidae (Hymenoptera: Ichneumoidea ) inferred from partial 16S rDNA gene sequences. Ins. Mol. Biol., 7: 129-150.
Mardulyn, P., and Whitfield,. J. B. (1999). Phylogeneticsignal in the COI, 16S, and 28S genes for inferring relationships among genera of Microgastrinae (Hymenoptera; Braconidae): evidence of a high diversification Rate in this group of Parasitoids. Molecular Phylogenetics and Evolution.12 (3) 282-294.
Simon, C., F., Beckenbach, A., Crespi, B., Liu, H., and Flook, P. (1994). Evolution, weighting, and phylogenetic utility of mitochoncrial gene sequences and a compilation of conserved polymerase chain reaction primers. Ann. Entomol. Soc. Am. 87: 651-701.
Starks, P.T., and Peters, J.M. (2002). Semi-noondestructive genetic sampling from live eusocial wasps, Polistes dominulus and Polistes fuscatus. Insectes Soc. 49:20-22.


This document was last modified 01/31/2006.
This site is powered by the versatile Zope platform.
This is a project of the Biology Department of Fordham University
Biotechniques.org Home