|
Mem Inst Oswaldo Cruz, Rio de Janeiro, Vol. 96(2) Feb. 2001, pp. 189-196 Genetic Variability among Populations of Lutzomyia (Psathyromyia) shannoni (Dyar 1929) (Diptera: Psychodidae: Phlebotominae) in Colombia Estrella Cárdenas/+, Leonard E Munstermann*, Orlando Martínez**, Darío Corredor**, Cristina Ferro Laboratorio de Entomología, Instituto Nacional de Salud, Avenida Eldorado, Carrera 50, Zona Postal 6, Apartado Aéreo 80080, Bogotá DC, Colombia *Department of Epidemiology and Public Health, School of Medicine, Yale University, New Haven, CT, USA **Facultad de Agronomía, Universidad Nacional de Colombia, Bogotá DC, Colombia +Corresponding author. Fax: +57-1-315.7341.
E-mail: ecardenas95@latinmail.com
Received 27 January 2000 Code Number: oc01027 Polyacrylamide gel electrophoresis was used to elucidate genetic variation at 13 isozyme loci among forest populations of Lutzomyia shannoni from three widely separated locations in Colombia: Palambí (Nariño Department), Cimitarra (Santander Department) and Chinácota (Norte de Santander Department). These samples were compared with a laboratory colony originating from the Magdalena Valley in Central Colombia. The mean heterozygosity ranged from 16 to 22%, with 2.1 to 2.6 alleles detected per locus. Nei's genetic distances among populations were low, ranging from 0.011 to 0.049. The estimated number of migrants (Nm=3.8) based on Wright's F-Statistic, FST, indicated low levels of gene flow among Lu. shannoni forest populations. This low level of migration indicates that the spread of stomatitis virus occurs via infected host, not by infected insect. In the colony sample of 79 individuals, the Gpi locus was homozygotic (0.62/0.62) in all females and heterozygotic (0.62/0.72) in all males. Although this phenomenon is probably a consequence of colonization, it indicates that Gpi is linked to a sex determining locus. Key words: sand flies - Lutzomyia shannoni - isozyme electrophoresis - genetic variability - dispersal gene flow - Colombia Of the more than 350 species of New World phlebotomine sand flies, only the visceral leishmaniasis vector, Lutzomyia longipalpis, has been extensively examined in terms of its population genetic structure among local populations (Morrison et al. 1995, Munstermann et al. 1998, Mutebi et al. 1998), and over a broader geographic range (Mukhopadhyay et al. 1998, Lanzaro et al. 1998, Mutebi et al. 1998, 1999). The forest sand fly Lu. shannoni, is a zoophilic species that has an even greater geographical distribution in the Western Hemisphere, extending from the southeastern United States to northern Argentina. In Colombia, this species is widely distributed from sea level up to 1,300 m (Young 1979). It is a vector of vesicular stomatitis virus (VSV) (Comer et al. 1990, 1991) and can support the development of at least three species of Leishmania: Le. mexicana (Lawyer et al. 1987), Le. chagasi (Endris et al. 1982) and Le. panamensis (Ferro et al. unpublished). The biology of Lu. shannoni has been studied by Comer et al. (1994) and Memmott (1991, 1992) in field conditions and under laboratory conditions by Ferro et al. (1998) and Cárdenas et al. (1999). However, no information on population genetic structure for this phlebotomine is available. In the present study, baseline parameters of genetic variation and interpopulation migration rates were estimated and genetic diversity was compared at 13 isozyme loci among forest populations of Lu. shannoni from three widely separated locations in Colombia and compared with a laboratory colony originating from the Magdalena Valley in Central Colombia. MATERIALS AND METHODS Sampling information - Sand flies were collected during day time (10 h-14 h) by aspirating them from tree trunks in a traditional coffee plantation at 1,235 m elevation near Chinácota, Department of Norte de Santander (7o34'N, 72o37'W); in a primary forest in Palambí at 68 m elevation (1o45'N, 78o28'W) near the city of Tumaco, Department of Nariño; and in a secondary forest called San Miguel adjacent to a livestock farm at 113 m elevation (6o23'N, 74o22'W) near Cimitarra, Department of Santander. The three sites are separated by eastern and western branches of the Andean mountain ranges, the Cordillera Oriental and Cordillera Occidental. Chinácota lies on the eastern slope, Cimitarra in the Magdalena Valley in Central Colombia, and Palambí on the Pacific Coast on the border with Ecuador (Fig. 1). After collection, the sand flies were separated by sex for visual identification, then the specimens were placed in 2 ml cryovials, and preserved in liquid nitrogen. In addition, flies from a laboratory colony were analyzed for comparison. The colony was derived from flies collected in 1992 from secondary forests, either in Carmen de Bolívar, Bolívar Department (9o43'N, 75o7'W) or Magdalena Valley Puerto Boyacá, Boyacá Department (5o58'N, 74o22'W) (see Results and Discussion). Taxonomic determination - Field specimens of Lu. shannoni were identified to the species level by morphological characteristics of the male terminalia and female spermathecae (Young & Duncan 1994), both obtained by removing the terminal three segments of the abdomen. The remainder of the specimen was preserved in a microcen-trifuge vial and frozen at -80ºC until used for isozyme analysis. Specimens identified as other species were preserved as voucher collections. Electrophoretic method - Samples were analyzed for electrophoretic variation at 13 isozyme loci. Individuals of Lu. shannoni were prepared for electrophoresis by homogenization in 10 or 20 µl of grinding solution (10% sucrose, 0.1% Triton X-100, 0.02M tris-citrate, and bromophenol blue). One to 3 µl of the homogenate were applied to vertical 6% polyacrylamide gels with the volume depending on the enzyme to be assayed. Gel preparation was based on standard procedures (Munstermann 1994a). Two buffer systems (0.025 M tris-citrate, pH 7.0 and 0.1 M tris-borate-EDTA, pH 8.9) maximized electrophoretic separation of isozymes (Table I). Enzyme-specific histochemical staining procedures were employed to visualize isozyme phenotypes on gels (Munstermann 1994a). Samples placed on gels migrated from the cathode to the anode. The loci elucidated on a single gel were numbered sequentially from slowest to most rapidly migrating. In each locus, the allozymes were designated according to their relative mobilities (Rf) (May 1992). The Rf values were estimated as a ratio of distance between the origin of the gel and the leading edge of an enzyme activity zone divided by the location of the activity zone of a commercially available enzyme standard (Table I). Data analysis - Gene and genotype data were analyzed with the BIOSYS-1 computer program (Swofford & Selander 1981) to determine gene frequencies, genotype frequencies in equilibrium Hardy-Weinberg, average heterozygosities and genetic distances. The number of migrants among populations was based on Wright's FST value and estimated by Nm=1-FST/4FST (Tabachnick & Black 1996), where Nm = number of migrants every two generations and FST inbreeding coefficient among populations. RESULTS AND DISCUSSION Genetic variability - Eleven of the 13 isozyme loci were polymorphic in all populations. The most variable locus was Est-2 with seven alleles, followed by Fum and Gpi with four alleles per locus. The following alleles showed gene frequencies <5%: Ak7.65 and Est-22.00 in the laboratory colony; Gpi0.50 in Palambí; Idh2.30 in Chinácota; and Ak5.45, Est-21.58, Gpi0.85, Hk0.83, and Mdh-11.25 in Cimitarra. The Ark and Est-1 loci were monomorphic in all samples. The remaining 11 loci were classified into three groups (Munstermann et al. 1998) based on mean heterozygosity (H): low (H= 0.01-0.1) Idh, Hk, Mdh-1; medium (H=0.1-0.2) Ak, Mdh-2; and high (H> 0.2) Aat, Me, Gpi, Gpd, Fum, Est-2. The laboratory colony contained seven monomorphic loci, whereas fewer were monomorphic in field samples: four in Palambi, five in Chinácota, and three in Cimitarra (Table II). In the laboratory colony, the locus Gpi was homozygotic in all females, whereas all males were heterozygotic (Fig. 2A). This phenomenon was never observed in field populations, although in the Palambí, Chinácota and Cimitarra samples, many more males than females were heterozygotes, (Fig. 2B). According to Bello et al. (1997), Lu. shannoni has four pairs of chromosomes with one pair heteromorphic (Jiménez et al. 2001). During the colonization of Lu. shannoni, the chromosome bearing the male sex locus apparently has trapped the 0.72 allele of Gpi, and the chromosomes bearing the female locus became homozygous for the 0.62 allele. In 79 colony individuals examined (39 males, 40 females), no female heterozygotes were observed, indicating no crossovers and a very tight linkage between the sex locus and Gpi. One consequence of the tight linkage is a forced heterozygosity at the Gpi locus, artificially increasing the "average heterozygosity" for the colonized strain. In other diptera, entrapment of variant alleles is not uncommon in the course of laboratory colonization and inbreeding (Munstermann 1994b). The average number of alleles per locus over all populations was 2.3 (range, 2.1-2.6). The percentage of polymorphic loci for population ranged from 38.5 to 69%; the mean heterozygosity ranged from 16 to 21.7%. The mean heterozygosity of the laboratory colony was very similar to the observed at Chinácota and Cimitarra, but the percentage of polymorphic loci and average number of alleles per locus were lower than in any of the wild populations (Table III). Conformation to Hardy-Weinberg Equilibrium - Among the 52 strain-locus combinations, 19 were monomorphic. Of the remaining 33 variable loci, 16 did not conform to frequencies predicted by the Hardy-Weinberg equilibrium. Explanations for these deviations are numerous, and some or all may apply to these 16 loci. Two assumptions in applying the H-W equilibrium are that (1) the samples were random and representative, and (2) that populations are in equilibrium. The first assumption is difficult to test (although samples were taken in from the field randomly as they were encountered) and the second is probably seldom true in a field setting. Furthermore, several of the 16 loci had common elements. The Est-2 and Fum loci were in disequilibrium in all of the four populations and therefore account for half of the 16 combinations (Table IV). These loci were also the most polymorphic and most prone to generation of large Chi-squares. This was probably a function of sample size and is a frequent result in population genetic studies with highly variable loci. Attempts to resolve these cases by linkage disequilibrium analyses have been largely unsuccessful (Tabachnick & Black 1996). Ak locus was in disequilibrium in Cimitarra and Palambí (monomorphic in Chinácota), and the remainder of the significantly deviant loci were distributed among populations and among loci without apparent pattern. However, all disequilibria are a consequence of heterozygote deficiency, with the exception of the Gpi deviation in the colony (Table IV). Given the direction of the disequilibrium, sampling error is an unlikely explanation for these occurrences, and consequently the observed heterozygote deficiencies may be due to localized inbreeding or possibly selection against heterozygous genotypes (Murphy et al. 1996). Finally, by the p>0.05 rejection of the null hypothesis alone, 1-2 may in fact not be in disequilibrium. Genetic distance - Nei's genetic distances among forest populations (Palambí, Chinácota and Cimitarra) ranged between 0.023 and 0.049, suggesting that the Colombian populations are genetically isolated, but still well with the range of divergence expected within a single species (see for example Tabachnick et al. 1979 for comparable distances among other haematophagus diptera). These relationships are illustrated in Fig. 3. The genetic distances are highly correlated with the geographical distances among populations: Palambí and Chinácota-0.049/950 km; Palambí and Cimitarra-0.038/680 km; Cimitarra and Chinácota-0.023/250 km. The values obtained are similar to those obtained by Lanzaro et al. (1998) for Colombia populations of Lu. longipalpis separated by the Andean ranges. Among local populations of Lu. longipalpis, Nm of 17 or greater were obtained (Munstermann et al. 1998); for more distant populations separated by mountain ranges, Nm was 2 (Lanzaro et al. 1998). The latter results were similar to those for the Lu. shannoni populations (Nm=3.8) described here, separated by 250 to 950 km. Although "in some migration models the genetic distance is linearly related to geographical distance or area" (Nei 1972), these models are not applicable for populations separated by highly heterogeneous landscapes. Despite a correlation between the genetic distances and geographic distances measured for Colombian Lu. shannoni, the routes of migration cannot be linear between sites _ particularly in view of the weak flight capabilities of sand flies (Chaniotis et al. 1974, Morrison et al. 1993) and separation by mountain ranges exceeding 3,000 m. F-statistics and gene flow - To better understand the breeding structure of local populations, the Wright statistical derivations were obtained. The laboratory colony was not included in this analysis. The inbreeding coefficient within subpopulations, FIS, was 0.300. An FIS > 0 indicates that the subpopulations are not in Hardy-Weinberg equilibrium, and a low value suggests reduction of heterozygosity due to inbreeding within the subpopulations (Tabachnick & Black 1996). The inbreeding coefficient among populations FST , was 0.061. The positive value suggested that gene flow was restricted among field populations. With the chi-square test, the null hypothesis was rejected at the probability level of 0.05 for the FST. This indicated the occurrence of gene flow among field populations, and therefore, a relatively low genetic divergence among them. The total population inbreeding coefficient, FIT, was 0.343. This value indicated a degree of inbreeding within subpopulations, producing a Wahlund effect among subpopulations (Tabachnick & Black 1996). Relevance to disease epidemiology - The migration rate estimates the dispersion potential of a pathogen transported from host to host by an insect vector. The low Nm=3.8 (migrants every two generations) for Lu. shannoni predicts that during an epidemic outbreak of VSV, the virus is more likely to be dispersed by the vertebrate host than by the insect vector. A high potential exists for implicating Lu. shannoni in transmission of Colombian arboviruses. Already, VSV virus (Alagoas serotype) has been isolated from Lu. shannoni collected in the secondary forests of Arboledas and Durania, Norte de Santander Department (Tesh et al. 1987). Note that Lu. shannoni is associated almost invariably with primary and secondary forests-many of which are situated adjacent to livestock regions (Entomology Laboratory, INS). Furthermore, the Instituto Colombiano Agropecuario (ICA) has recorded large numbers of cattle infected with the New Jersey and Indiana serotypes of VSV in nearly 100 Colombian municipalities (ICA 1999). The definitive vector has not been identified, but Lu. shannoni remains a primary suspect. Colony origin - Although initially the colony was recorded as originating from Carmen de Bolívar, Bolívar Department, the isozyme comparisons indicated a closer relationship to the Cimitarra population. When the archival collecting records were reviewed, we found that in March 1992 sand flies were collected in two forests - one near Puerto Boyacá, Boyacá Department, and the other near Carmen de Bolívar, Bolívar Department. Since San Miguel forest in Cimitarra is closer to Los Corales forest in Puerto Boyacá (70 km) than Carmen the Bolívar (480 km), the laboratory colony may have been derived from Los Corales, Puerto Boyacá instead of Carmen de Bolívar. Nei's genetic distances closely group the laboratory colony with Cimitarra (Fig. 3). Since this colony is the only colonized Lu. shannoni in the world and serves as source of insects for studies in behavior and vector competence, correct placement of its geographic origins may be relevant and may contrast significantly with populations of this species in other areas. ACKNOWLEDGMENTS To Marco Fidel Suarez, to personnel of the Secretaria de Salud of Norte de Santander Department, and to personnel of CIDEIN of Valle Department for their collaboration in field collections of Lutzomyia shannoni specimens. REFERENCES
This research was supported by the Colombian Instituto Nacional de Salud (Project No. 3100302 to CF of the Entomology Laboratory) and the United States National Institutes of Health (AI-34521 to LEM). Copyright 2001 Fundacao Oswaldo Cruz Fiocruz The following images related to this document are available:Photo images[oc01027t4.jpg] [oc01027t1.jpg] [oc01027f3.jpg] [oc01027t3.jpg] [oc01027t2.jpg] [oc01027f1.jpg] [oc01027f2.jpg] |
|