Quantitative trait loci in <it>Anopheles gambiae </it>controlling the encapsulation response against <it>Plasmodium cynomolgi </it>Ceylon

1471-2156-4-16 1471-2156 Research article Quantitative trait loci in <it>Anopheles gambiae </it>controlling the encapsulation response against <it>Plasmodium cynomolgi </it>Ceylon Zheng Liangbiao liangbiao.zheng@yale.edu Wang Shuang shuang.wang@yale.edu Romans Patricia victorie@zoo.utoronto.ca Zhao Hongyu hongyu.zhao@yale.edu Luna Coralia coralia.luna@yale.edu Benedict Q Mark mqb0@cdc.gov

Yale University School of Medicine, Department of Epidemiology and Public Health, 60 College Street, New Haven, CT 06520, USA

Department of Zoology, University of Toronto, Toronto, ON, Canada M5S 3G5

Centers for Disease Control and Prevention, Mailstop F22, Chamblee, GA 30334, USA

BMC Genetics 1471-2156 2003 4 1 16 http://www.biomedcentral.com/1471-2156/4/16 10.1186/1471-2156-4-16 14577840 10 7 2003 24 10 2003 24 10 2003 2003 Zheng et al; licensee BioMed Central Ltd. This is an Open Access article: verbatim copying and redistribution of this article are permitted in all media for any purpose, provided this notice is preserved along with the article's original URL. Anopheles gambiae Plasmodium cynomolgi Ceylon malaria refractoriness melanotic encapsulation QTL analysis.

Abstract

Background

Anopheles gambiae females are the world's most successful vectors of human malaria. However, a fraction of these mosquitoes is refractory to Plasmodium development. L3-5, a laboratory selected refractory strain, encapsulates transforming ookinetes/early oocysts of a wide variety of Plasmodium species. Previous studies on these mosquitoes showed that one major (Pen1) and two minor (Pen2, Pen3) autosomal dominant quantitative trait loci (QTLs) control the melanotic encapsulation response against P. cynomolgi B, a simian malaria originating in Malaysia.

Results

We have investigated the response of L3-5 to infection with P. cynomolgi Ceylon, a different but related parasite species, in crosses with the susceptible strain 4Arr. Refractoriness to this parasite is incompletely recessive. Infection and genotyping of F2 intercross females at genome-spanning microsatellite loci revealed that 3 autosomal QTLs control encapsulation of this species. Two loci map to the regions containing Pen2 and Pen3. The novel QTL maps to chromosome 3R, probably to polytene division 32 or 33. Thus the relative contribution of any QTL to oocyst encapsulation varies with the species of parasite. Further, different QTLs were most readily identified in different F2 families. This, like the F1 data, suggests that L3-5 is not genetically homogeneous and that somewhat different pathways may be used to achieve an encapsulation response.

Conclusion

We have shown here that different QTLs are involved in responses against different Plasmodium parasites.

Background

Melanotic encapsulation of invading parasites and parasitoids is a widespread insect defense mechanism 1234. Depending on the insect and the parasite involved, this can be either a cellular or a humoral process. The humoral form of melanotic encapsulation occurs primarily, if not exclusively in the order Diptera, and is correlated with low numbers of circulating hemocytes. Melanotic encapsulation of malaria parasites was first observed within a year of the discovery that mosquitoes are the vectors of malaria 5. Since then, melanotic encapsulation of plasmodia has been observed in a wide variety of mosquito species. Depending on the mosquito and the parasite species, encapsulation is directed against different stages of development 6.

Anopheles gambiae is the principal vector for human malaria in Africa. A laboratory-selected strain derived from the G3 colony of A. gambiae was found to deposit melanin capsules around transforming ookinetes or very early oocysts of a variety of Plasmodium species 78910. Since hemocytes do not appear to participate in formation of these capsules, the mechanism underlying encapsulation is likely to be humoral and/or to involve the midgut epithelium or its basal lamina. One locus, Pif-B, was postulated to control encapsulation of P. cynomolgi B. Comparisons between mosquitoes from this original refractory strain and a counter-selected susceptible strain showed that refractoriness was highly correlated with the 2L+^achromosomal arrangement and esterase A. Susceptibility was correlated with inversion 2La and esterase C. A similar correlation was found for the parental G3 colony. Thus it was suggested that Pif-B and the esterase A/C polymorphism are located within or very close to the 2La region (divisions 23–26 of ovarian nurse cell polytene chromosome 2L) 111213. Genetic mapping experiments with the current refractory (L3-5) and susceptible strains (4Arr) revealed three QTLs controlling the encapsulation of P. cynomolgi B. These loci were named Pen1-3 for

Plasmodium

encapsulation 1, 2 and 3. Pen1 is the major QTL whereas Pen2 and Pen3 are minor QTLs. The refractory alleles at these loci are largely dominant and the contributions of Pen1 and Pen2 to the encapsulation response appeared additive 14. The Pen1 QTL is also the major QTL for encapsulation of P. berghei and for melanin deposition on Sephadex-CM C-25 beads 15. Pen1 has been mapped to chromosome 2R, division 8C/D 16, while Pen2 and Pen3 locations remain less precisely defined, on distal chromosome 3L in or near division 43, and on proximal chromosome 2R, possibly in division 14, respectively 1718. None of these loci affects the intensity of Plasmodium infection as measured by number of oocysts on the midgut. Recently, the two loci on chromosome 2L were found to control the P. falciparum infection intensity in natural populations of A. gambiae 19.

In contrast, refractoriness to P. cynomolgi Ceylon in the original refractory strain seemed to be largely but not completely recessive 712. P. cynomolgi Ceylon was originally described as a Ceylonese subspecies of P. cynomolgi 20, whereas P. cynomolgi B originated in Malaysia. The two are currently regarded as being different species (W.E. Collins, personal communication) and molecular data support this contention 2122. Genetic mapping suggested that Pif-C, the locus controlling encapsulation of this parasite, assorted independently from the 2La inversion polymorphism, and thus was not linked to the left arm of chromosome 2. However, it was not possible to map this locus further because of the lack of polymorphic markers. More recently, a locus near Diphenol oxidase-A2 (3R; 33B) on the right arm of the chromosome 3 was also suggested to be important in the melanization of P. cynomolgi Ceylon 23. The refractory allele at this Dox-A2 linked locus seemed to be most highly expressive when the mosquitoes were heavily infected (geometric mean 178.3, average 212.1 ± SD 104.9 oocysts per midgut). It is not clear whether this locus is the same as Pif-C. It is definitely not Pen2. Further, it was suggested that a locus or loci unlinked to Dox-A2, but also likely to be on chromosome 3, was involved in refractoriness at lower infection intensity, (geometric mean 28.3, average 40.6 ± SD32.7 oocysts per midgut), and may have contributed to parasite encapsulation at high infection intensity as well.

In this study we have investigated the genetics of the encapsulation response against P. cynomolgi Ceylon in crosses between the extant refractory L3-5 and susceptible 4Arr strains using microsatellite markers spanning the genome and QTL analysis.

Results

Refractoriness to P. cynomolgi Ceylon is incompletely recessive

Infection of parental line mosquitoes showed that L3-5 refractory strain mosquitoes encapsulated P. cynomolgi Ceylon, whereas 4Arr did not (Fig. 1A). In three separate infection experiments on a total of 54 individual L3-5 mosquitoes, the overall infection prevalence was 81.5%. The number of oocysts in infected individuals ranged between 1 and 147. All oocysts were encapsulated. The infection intensities were generally low and quite variable. The average (± standard deviation) of the intensity of infection varied from a low of 4.93 ± 2.62, to a high of 47.80 ± 40.95 between infections. The overall intensity of infection was 26.36 ± 34.09.

Figure 1

Melanotic encapsulation of P. cynomolgi Ceylon by A. gambiae

Melanotic encapsulation of P. cynomolgi Ceylon by A. gambiae. (A) Responses of individual parental L3-5 refractory strain (squares and circles) and 4Arr susceptible strain females (diamonds and triangles) to P. cynomolgi Ceylon infection. Exp.1 and exp. 2 indicate two separate infection experiments. (B) Responses of F1 female progeny from one cross between L3-5 females and 4Arr males (closed diamonds) and three crosses between 4Arr females and L3-5 males (open squares, circles and triangles).

In parallel infection experiments on a total of 60 4Arr susceptible strain mosquitoes blood-fed on the same monkeys used to feed the refractory ones, the infection prevalence was 70.0%. Despite the apparent difference in infection prevalence, the overall intensity of infection was very similar to that of the L3-5 females, 21.50 ± 19.66. Further, the lowest and highest infection intensities, 7.00 ± 7.51 and 32.94 ± 23.12 in 4Arr corresponded to the lowest and highest infection intensities for L3-5 that fed on the same monkeys on the same days. Oocyst numbers ranged between 1 and 79. One encapsulated oocyst was found in each of two 4Arr mosquitoes. These data suggest that if there are differences in P. cynomolgi Ceylon infection intensities, they are more likely to be due to the infection parameters of the vertebrate host than to differences intrinsic to the parental mosquito lines.

Reciprocal crosses between L3-5 and 4Arr yielded F1 progeny that were mostly susceptible by the criterion they had no or a low fraction of encapsulated oocysts (Fig. 1B). When L3-5 females were mated with 4Arr males, 33 tested F1 progeny had an average of 27% encapsulated oocysts (infection intensity of 12.46 ± 7.31). The reciprocal cross yielded 101 F1 progeny with an average of 39% encapsulated oocysts (infection intensity of 18.33 ± 17.10). Thus refractoriness to P. cynomolgi Ceylon in L3-5 mosquitoes appeared to be largely recessive, but less so than observed in the originally selected strain 24. Further, these data give no evidence for any maternal effects on the encapsulation response.

Phenotypic responses of F2 females against P. cynomolgi Ceylon

The encapsulation response of 167 females in six F2 intercross families was tested in three separate feeding experiments (Table 2). Family 8–6 females fed on one infected monkey while females in the other five families fed on a different monkey on two successive days of ascending parasitemia. The infection intensity for the combined F2 progeny was 22.96 ± 36.14. The majority of the F2 progeny were either highly refractory or highly susceptible (Table 2; Figure 2), with 58 (35%) and 47 (28%) individuals exhibiting ≥ 80% or ≤ 20% oocysts encapsulated. 62 (37%) of the females had intermediate fractions of encapsulated oocysts. These data suggest that alleles controlling melanotic encapsulation are segregating at more than one locus.

Table 1

New microsatellite markers and primers developed for this study.

Marker name

Genetic (Cytogenetic) location

Forward primer (5' – 3')

Repeat type

Reverse primer (5' – 3')

Sequenced allele size (bp)

Chromosome 2

H682

(2R, 7B)

CGGGCGGTCGTTATTAGGAG

(AC)7

CTGAGTTGACACGAGTTGCC

3F11-C1

(2R,9B)

AGATGGTGTCCTTCGGTGTC

(GT)17

GTGGAGGAGAAGCCACACAT

156

3B20-B1

(2L, 24B)

GAGACCAAAACATGTCAACCG

(CA)17

GTATGTCTGCGCATGAACGAT

122

R007

71.1 (2L, 26D)

GCACCAGTCAAATACGCACT

(CA)3+9

GTCCGTGTTGTTCGAAGATG

184

Chromosome 3

H249

11.9 (3R, 30A)

GGAAGCAGCTGCCTTTTATG

(TG)15

ACAATGTAACACAGCGCGAG

161

3A13-12

(3R, 30E)

GAACAAAACGTGCCACAGC

(GT)9

TCCCCGACTTACGTGAACC

168

3A7-37

(3R, 31B)

CGATTGGAACCAACAAGTGG

(CA)2+9

AAGCCATCTCTTGACGTTCG

151

Ddc*

(3R, 33B)

CGTTGAGAACCGAAAACAGG

(GT)10

ATCGGGACACTCAGCAACAG

121

3A6-17

(3R, 33C)

TGTGAAGAATGTTTTCGCTGAC

(AC)17

CTCCCCCGTGTAATTTGAAC

120

9D17-42

(3R, 33C)

TAAACCTCACCGTCCACATC

(TG)23+5

AAATTGAAACCCGTCCAAAG

150

6B10-1

(3R, 33D)

CATGAGTGCCGCTGTAAATG

(TG)2+16

ATCGCTCTCCTCCTCTAGCC

120

3A23-100

(3R, 34A)

AATGGACCATCGTTAGCACAG

(TG)12

TCCAAAGTTTTCCCTTCCAC

136

3L09-C1

(3L, 38C)

GAAACTTCCCTGAATTTCCGC

(CA)10+8

ATGGAGGCGACTAAACGTTGC

169

3A2-29

(3L, 41A)

GGCTGCTCTAGGAGGATGTG

(CT)16

AATCGAGGGTGGTGTGAATG

189

3O10-C2

(3L, 41D)

TTAATCACCCACCTTGATCGC

(CA)18

CCTAGAGCTTGGAATCGAATC

101

3A14-11

(3L, 43D)

GTAATAGCGAGCGAGCGAAC

(TG)36

GAGCGATCCGAATCGTTAGT

180

3A14-47

(3L, 43D)

GGAGGCACTCTACCACAACC

(GT)10

CCGAACGTCATCGCTAATG

104

3A19-32

(3L, 43D)

CGTGCTGGCATTAAAAAGTG

(CA)15

CTACGCATGCATCTTTCTCC

138

3A5-25

(3L, 45A)

CCTTGGGCGCTATAAAACAC

(TG)10

TATTAGCAGCAGGGGAAACC

131

19A18-54

(3L, 46A)

TAATGCAGTGGTGCGGTTC

(GT)5+6

AATCACAGCCCAGTCAGTCC

121

19C20-B1

(3L, 46B)

ACCAAGTGCTCTTTGGGTTG

(TG)15

ACGAATGCTGGACAATGTGA

141

The cytogenetic location of H682 is A. gambiae scaffold predicted by it's location in AAAB01008987, contig 566 (GenBank accession AAAB01008987). New primer pairs were designed for markers R007 44 and H249 14. Sequenced allele size is expressed in base pairs (bp). * This Dopa decarboxylase microsatellite is contained within the large first intron of predicted Ddc isoform 1. DNA sequence is from G3 strain (GenBank # AF063021; P. Romans, unpublished).

Table 2

Prevalence and intensity of oocyst infections in the F2 progeny

Family

Total # of F2 female progeny

Prevalence*

Infection range#

Average ± SD

Parasite/μl †

8-3

0.67

1 – 152

35.25 ± 39.73

34,750

8-4

0.75

3 – 76

23.8 ± 24.28

34,750

8-6

0.55

2 – 30

8.50 ± 7.22

19,080

10-1

0.97

3 – 149

44.89 ± 34.42

34,750

12-1

0.65

9 – 180

65.15 ± 44.38

34,750

13-4

0.76

4 – 104

41.11 ± 26.78

34,750

Total

167

0.74 (123/167)

1 – 180

39.04 ± 36.14

* Fraction of F2 females that became infected. # Number of oocysts in the infected midguts. † Number of parasites per microliter of blood at the time of feeding.

Figure 2

P. cynomolgi Ceylon infection phenotypes of F2 female progeny from crosses of six F1 females with their sibling males

P. cynomolgi Ceylon infection phenotypes of F2 female progeny from crosses of six F1 females with their sibling males. (A) Intensity of infection was plotted against the total number of oocysts, both normal and encapsulated, per female midgut. N = 167. (B) Encapsulation responses of 123 infected F2 females. The number of F2 females was plotted against the fraction of oocysts encapsulated per midgut.

Genotyping and genotype frequencies

Genotyping was performed on these 167 female progeny and their six female F1 parents at 52 microsatellite loci spanning the current genetic map of A. gambiae. These included 6, 22, and 24 loci on the X, second and third chromosomes, respectively (see Materials and Methods). Six markers were monomorphic in all families and thus uninformative (H678 on the X; H682 and H117 on the second chromosome; H776, H525, and H544 on the third chromosome). Marker H197 on the second chromosome gave uninterpretable data.

Only two alleles were found for each informative marker on the X chromosome. However, among the six F1 female parents of the F2 families, more than two alleles were found at 11 autosomal loci, indicating that one or both of the parental strains was polymorphic at these loci. The H750 locus, which maps to 65.7 centiMorgan (cM) on the third chromosome, had four alleles. Furthermore, the two parental strains also share some common alleles. This is exhibited by the observation that one of the six F1 females (12-1) was homozygous at 10 autosomal loci. In total, 5, 16 and 14 loci were informative on the X, second and the third chromosomes, respectively. Genotypes of the F2 female progeny were tested for deviations from the expected Mendelian ratios at each locus within each family. Although significant deviations were found for some loci, these loci were not clustered in any chromosomal region (data not shown).

Quantitative trait loci for encapsulation of P. cynomolgi Ceylon

No QTL was identified for infection intensity in the previous study 14. Nor did our infection data suggest any genetic variation for infection intensity between L3-5 and 4Arr. We thus focus our attention on the melanotic encapsulation phenotype, which we have expressed as fraction of encapsulated oocysts in a single mosquito. Linear regression analysis was performed for each locus within each family and for all families combined. Considering the combined data, no significant association was found between the oocyst encapsulation phenotype and any marker on the X chromosome (Table 3). This is similar to previous results obtained for P. cynomolgi B and all other parasites tested.

Table 3

Probability of a QTL for P. cynomolgi Ceylon encapsulation on the X chromosome.

Locus

Genetic (Cytogenetic) Location

Maximum F statistic*

R ^2†

Average of Family Maximum F statistics ^‡

P ^§

P' ^§

Average of Family R²

H145

1.9 (4C)

0.64

0.0052

0.31

0.956

1.0

0.0194

H503

11.8 (4B)

2.71

0.0217

0.96

0.612

1.0

0.0365

H036

13.1 (4B)

0.37

0.0032

0.94

0.605

1.0

0.0557

H053

25.0 (3D)

4.16

0.0335

2.67

0.041

0.205

0.1040

H711

40.8 (2A)

5.94

0.0495

1.57

0.273

1.0

0.0677

* Maximum F statistic calculated using all informative data. †Calculated from all informative data. ‡ Average of the maximum F statistics for all informative families. ^§Probability values P and P' were calculated based on 1000 permutation tests before and after Bonferroni correction respectively.

Surprisingly, no major QTL was identified on chromosome 2 near H290 or H175 (Table 4), the markers that most precisely defined Pen1 14. Thus it is unlikely that a Pen1 allele has a significant role in the encapsulation of P. cynomolgi Ceylon. 3F11-C1 (2R, 9B) was the only marker near the Pen1 cytogenetic region that might be associated with an encapsulation QTL. 3F11-C1 had a very high average maximum F statistic of 200.88 and gave an F statistic calculated from data on all informative individuals of 11.74, though the P' value associated with it, 0.16, is not significant. More proximally, H135 and the closely linked marker H095 appear to define a QTL, which is termed Pcen2R (for

P. cynomolgi

Ceylon

encapsulation, chromosome arm

2R). These two markers are the same markers most closely associated with Pen3, a minor QTL for P. cynomolgi B encapsulation. H095 exhibited the highest average maximum F statistic of this data set, 287.99, with an overall maximum of 12.70, the highest of any marker on chromosome 2. However these overall maximum F statistics were lower than similar statistics associated with many chromosome 3 loci (Table 5).

Table 4

Probability of a QTL for P. cynomolgi Ceylon encapsulation on chromosome 2.

Locus

Genetic (Cytogenetic) Location

Maximum F statistic*

R ^2†

Average of Family Maximum F statistics ^‡

P ^§

P' ^§

Average of Family R²

H417

0.0 (2R, 7B)

4.19

0.0398

3.48

0.069

1.0

0.1602

H784b

2.1(2R, 7A – 10D)

2.23

0.0304

2.55

0.37

1.0

0.1639

H175

3.1 (2R, 8C)

6.93

0.0901

5.76

0.039

0.624

0.2347

H290

4.3 (2R, 8C)

0.32

0.0043

3.21

0.167

1.0

0.1733

H799

18.7 (2R, 9B)

3.30

0.0839

3.30

0.074

1.0

0.0839

3F11-C1

(2R, 9B)

11.74

0.1240

200.88

0.01

0.16

0.4640

H187

24.7 (2R, 11A)

6.56

0.0588

5.15

0.024

0.348

0.1755

H095

42.3 (2R,13B)

12.70

0.1313

287.99

0.001

0.016

0.4184

H135

47.3 (2R, 14D)

9.48

0.0721

258.33

0.000

0.3389

H770

50.1 (2R, 15D – 16A)

5.85

0.0542

3.06

0.071

1.0

0.1399

3B20-B1

(2L, 24B)

0.13

0.0011

3.14

0.171

1.0

0.1297

H772

52.9 (2L, 25D – 28D)

1.19

0.0097

2.06

0.219

1.0

0.0901

H603

54.7 (2L, 26D)

0.04

0.0004

3.35

0.132

1.0

0.1381

H819

55.1(2L, 25D – 28D)

2.61

0.0305

3.28

0.143

1.0

0.1572

H802

57.0 (2L, 26D)

1.52

0.0127

4.69

0.027

0.432

0.1606

R007

71.1 (2L, 26D)

1.26

0.0105

3.67

0.109

1.0

0.1578

Conventions are as for Table 3. Underlined markers

H175

and

H135

indicate maximum LOD scores for Pen1 and Pen3 respectively 14. Cytogenetic map positions have been used to order loci on 2L where possible. Cytogenetic and genetic map positions are derived from 14444546 and unpublished. Cytogenetic locations were further confirmed by comparison to the current scaffolds of the A. gambiae genome (GenBank accessions AAAB01000001–AAAB01008987). Where the precise location of a marker is not known, its cytogenetic location is shown as the range covered by the corresponding scaffold.

Table 5

Probability of a QTL for P. cynomolgi Ceylon encapsulation on chromosome 3.

Locus

Genetic (Cytogenetic) Location

Maximum F statistic*

R ^2†

Average of Family Maximum F statistics ^‡

P ^§

P' ^§

Average of Family R²

H249

11.9 (3R, 30A)

0.40

0.0451

6.03

0.019

0.285

0.2217

3A13-12

(3R, 30E)

2.27

0.0451

1.96

0.18

1.0

0.1087

H158

46.8 (3R, 32D)

5.43

0.1576

8.73

0.005

0.075

0.2690

H555

37.1 (3R, 33B)

5.97

0.0786

7.03

0.002

0.03

0.2343

3A6-17

(3R, 33C)

32.59

0.2208

9.77

0.024

0.36

0.3208

6B10-1

(3R, 33D)

6.36

0.0696

5.40

0.054

0.81

0.2120

H750

68.0 (3L, 39A)

26.59

0.1789

5.73

0.006

0.09

0.2063

H127

68.2 (3L, 39A)

30.20

0.1984

7.14

0.003

0.045

0.2539

3A2-29

(3L, 41A)

31.13

0.2658

30.4096

0.004

0.06

0.4906

H577

71.9 (3L, 42A)

7.74

0.1770

7.74

0.015

0.225

0.1770

H758

81.0 (3L, 43B)

18.24

0.1401

9.2693

0.006

0.09

0.3094

3A19-32

(3L, 43D)

16.22

0.1209

8.8599

0.01

0.15

0.3033

H753

88.7 (3L, 43D – 46D)

9.24

0.1036

7.5450

0.014

0.21

0.2060

3A14-47

(3L, 43D)

7.86

0.1792

7.8612

0.02

0.30

0.1792

Conventions are as for Table 3. Underlined marker

H758

indicates maximum LOD score for Pen2 14

There are at least two encapsulation QTLs on chromosome 3 (Table 5). Markers H555/H158 define a novel QTL in the central region of 3R, in division 32 or 33 that we have named Pcen3R. Marker 3A2-29 (3L, 41A) defines another QTL on 3L, Pcen3L, which covers the region spanned by divisions 39 – 43. This QTL may explain the relatively high F statistics exhibited for markers H750/H127 in division 39 to H758 in division 43B. It should be noted that H758 is the marker most closely associated with Pen2 and that marker 3A2-29 was not available at the time of the P. cynomolgi B encapsulation studies.

When the complete set of progeny was analyzed as a whole, the QTLs on 2R, 3R and 3L explained 13%, 16% and 26% of the refractoriness to P. cynomolgi Ceylon, respectively. Similar results were obtained when only all the families except 8-6 (which fed on a different monkey with different parasite level) were analyzed (data not shown).

The refractory and susceptible lines were not highly inbred except for malaria refractoriness or susceptibility and In(2L)+^a(see Methods). This posed a significant challenge to determining the parental origin of alleles at some loci, in some families. In these cases, the phase information of each allele was predicted using a maximum likelihood method. When a Hidden Markov Model was used to determine the number of QTLs based on the derived phase information, it yielded similar results with two QTLs on the third chromosome and one on the second chromosome. These QTLs were linked with the same marker loci as those identified using single locus linear regression analysis (data not shown, S.W and H. Z, unpublished).

F statistics were also calculated for each marker within each individual F2 family. The results for all significant QTLs and markers associated with high F statistics in any family are presented in Table 6. Data for H175, the marker closest to Pen1, are included for comparison. There appeared to be considerable F2 family specificity as to which QTLs contributed to refractoriness. Much of the H095 F statistic was due to family 8-4.

Table 6

F statistics for markers associated with significant QTL loci within each individual family.

Family

8-3

8-4

8-6

10-1

12-1

13-4

Locus

H175

9.56

0.0

1.15

6.56

3F11-C1

4.16*

757.55

42.00

3.30

1.17

0.69

H095

0.89

1408.89

18.75

7.05*

3.16

0.90

H135

8.27

1522.87

8.09

7.05

3.16

0.57

H158

0.003

5.95

17.45

1.23

11.02

1.27

H555

0.05

3.98*

17.45

11.11

9.35

1.28

3A6-17

8.59

3.98

5.80

26.23

8.02

5.99

H750

8.59

5.95

4.69

14.05

1.09

0.04

H127

8.59

5.95

4.69

14.15

9.35

0.23

3A2-29

8.59*

8.03

4.69*

37.62

69.99

5.99

H758

8.27

2.94

8.40

10.19

23.84

1.98

3A19-32

8.27

2.94

9.39

8.78

20.66

3.14

*The marker was not informative in this family. The F statistic shown was from the cytogenetically closest neighboring informative locus in the same family. Underlined markers

H135

and

H758

defined P. cynomolgi B encapsulation QTLs Pen3 and Pen 2 respectively 14

Pcen2R explained the majority of the encapsulation phenotype in family 8-4, less of the phenotype in family 8-6, very little in family 13-4 and contributed marginally to the phenotype in the other three families. Pcen3R appeared to be most important in determining refractoriness in family 8-6 and less important in families 10-1 and 12-1. Pcen3L on the other hand, appeared to be the major determinant of refractoriness in families 10-1 and 12-1. It probably played a minor role in the other families.

Discussion

Melanotic encapsulation of P. cynomolgi Ceylon parasites is controlled by incompletely recessive alleles in A. gambiae L3-5. Similar results were obtained in previous analysis of the original R and S strains with this parasite 1223. Refractoriness of L3-5 to P. cynomolgi Ceylon is shown here to be controlled by at least three quantitative trait loci that we have named Pcen2R, 3R and 3L. Two of the QTLs (Pcen2R and Pcen3L) are located near Pen3 and Pen2 respectively, previously identified as minor QTLs involved in the encapsulation response against P. cynomolgi B. These loci may actually be Pen2 and Pen3. A third QTL, Pcen3R, is most closely associated with microsatellite locus H158. This QTL, therefore, is located in the central region of chromosome 3R, in or near division 32 or 33. This is a novel QTL. A nearby locus, H555, was associated with a local peak in LOD score in the P. cynomolgi B studies 14. However, in that analysis, this peak in LOD score did not appear as a QTL separate from Pen2. The previously identified major QTL for encapsulation of P. cynomolgi B, Pen1, contributed little to the encapsulation of P. cynomolgi Ceylon. Thus, it appears that different genetic loci may be involved in responses to infections by different malaria parasites or, perhaps more precisely, may make malaria species-dependent contributions to the genetic variance in encapsulation.

Initial studies of A. gambiae refractoriness to P. cynomolgi Ceylon used non-homosequential mosquito strains in which a major difference between the refractory and susceptible strains was the absence or presence of inversion 2La. This inversion encompasses polytene chromosome divisions 23 through 26. It was absent in the refractory strain and fixed in the counter-selected susceptible one. Those studies suggested that melanotic encapsulation of this parasite was controlled by a largely recessive allele at one locus, Pif-C 12. This locus assorted independently from Pif-B, a major locus for encapsulation of P. cynomolgi B that appeared to be tightly associated the 2L+^achromosomal arrangement 111213. Vernick et al. 12 also found that a refractory allele at Pif-B acted as a dominant enhancer of expressivity of the refractoriness to P. cynomolgi Ceylon. A later study using the same R strain and a backcross strategy ensuring that each offspring received at least one refractory allele of Pif-B showed that the refractory phenotype is controlled by at least two, largely recessive loci on chromosome 3 23. One of these loci appeared linked to Dox-A2 (3R; 33B). Additional analysis of those data suggests that this locus maps approximately 32 cM from Dox-A2. (P. Romans, unpublished). Therefore this locus is expected to be located extremely distally on chromosome 3R, or else to be tightly linked to the centromere on 3R or 3L. The other locus or loci assorted independently of Dox-A2 but did not assort independently of the Dox-A2 linked refractoriness locus. It was not mapped further because of lack of genetic markers. This study provided strong evidence for variable penetrance and expressivity of refractory alleles at the Dox-A2 linked locus. As indicated by Romans et al. 23, this locus was more completely expressive when the geometric mean of the parasite load was 178 oocysts per midgut, rather than when it was 28. Since the infection intensities in the present study were similar (geometric mean 22.96, average 39.04 ± SD 36.14) to this lower number, it is possible that the Dox-A2 linked locus could have escaped detection even though most of the third chromosome was subjected to QTL analysis. Nevertheless, the 3A2-29 QTL in proximal 3L maps at the appropriate distance from Dox-A2 to be the Dox-A2-linked gene. Further, it is linked to H758 in division 43, a marker that would assort independently of Dox-A2. Therefore, if H758 is not part of the QTL associated with marker 3A2-29, but is a separate QTL, markers 3A2-29 and H758 could define the loci inferred in the previous study.

Our observation that different genes may be responsible for encapsulation of different parasite species raises the possibility that different QTL(s) might be required for melanotic encapsulation of human malaria parasites, such Plasmodium falciparum. The previous selected refractory strain encapsulated P. falciparum isolates from Africa less effectively 7. Whether refractoriness to P. falciparum is dominant or recessive has not been determined for the previous selected or the L3-5 refractory strain of A. gambiae.

Our observation that different genes may be responsible for encapsulation of different parasite species is not unique. In Drosophila melanogaster, a strain that encapsulates eggs of the parasitic wasp Leptopilina boulardi L104, a Congolese strain, is susceptible to parasitism by L. boulardi G317 from Tunisia and by L. heterotoma 25262728. Further, genetic mapping experiments have shown that the encapsulation responses to L. boulardi and Asobara tabida are controlled by single dominant, but different loci. The locus for L. boulardi encapsulation, Rlb, maps to chromosome 2 at approximately 86.7 cM, in the region 55E3 – E6;F3 on the polytene map 272930. Rat, the locus for A. tabida encapsulation, also maps to chromosome 2R, but at 51.3 cM 2630. The identities of these genes have not yet been determined. Benassi et al., 26 have suggested that the products of these loci should be involved in parasitoid recognition.

In the G3 parental A. gambiae strain from which the previous and current inbred refractory strains were derived, >50% of females encapsulate P. cynomolgi B, though only 1 – 3 % of females encapsulate P. cynomolgi Ceylon. Recently, field-caught A. gambiae sensu lato mosquitoes from Tanzania were tested for their ability to melanize P. falciparum and CM Sephadex C-25 beads 31. While fewer than 1% of these mosquitoes were able to encapsulate P. falciparum oocysts, 90% melanized the beads. Even though these results with a natural population of mosquitoes could be confounded by species, age, nutritional and reproductive status 32, they are intriguingly similar to results obtained with G3. It is possible that genetic polymorphism in the encapsulation response may be maintained by selection in natural populations.

The identities of the genes contributing to malaria refractoriness in A. gambiae are currently unknown. Both susceptible and refractory strains deposit melanin in egg chorions and in the cuticles of larvae and adults. Thus it is likely that refractoriness loci encode variant molecules that somehow expose parasites as foreign, or directly affect the signaling pathways involved in activation of the phenoloxidase cascade. Encapsulation occurs just as ookinetes penetrate through the basal labyrinth of the midgut, contact the basal lamina and become exposed to the contents of the hemocoel 8910. It is also possible that defects in the basal lamina could be responsible for exposing parasites to soluble receptors and components of the phenoloxidase cascade.

It is quite possible that soluble receptors are involved in parasite recognition and initiation of encapsulation. Paskewitz and Riehl 33 showed that if Sephadex CM C-25 beads had been exposed first to hemolymph from susceptible A. gambiae 4Arr, they could not subsequently be encapsulated by refractory L3-5 mosquitoes. This result could indicate that the receptors are soluble, although other interpretations are possible.

L3-5 refractory mosquitoes appear somewhat darker than 4Arr susceptible ones at all stages beyond L2. Mortality associated with molting from L3 to L4, from L4 to pupa and from pupa to adult is characterized by extreme melanization in the L3-5 strain (P. Romans and C. Brotherston, unpublished). It is unlikely that this difference is due to phenoloxidase enzymes themselves: Brey et al. 34 showed that the original refractory line of A. gambiae and its G3 parent did not differ significantly in phenoloxidase activity. Thus it is also possible that some refractoriness loci encode proteins involved in regulating activation of the phenoloxidase cascade, putting it on a sensitive trigger downstream of any receptor, or that they encode enzymes involved in production of phenoloxidase substrates.

Methods

Mosquitoes and parasites

Mosquitoes and parasites were maintained at the Centers for Diseases Control and Prevention, Atlanta GA as described previously 3536. The derivation of the A. gambiae malaria refractory L3-5 and malaria susceptible 4Arr strains used in these crosses is given by 1437. Reciprocal crosses between L3-5 and 4Arr were carried out en mass. F1 offspring of each reciprocal cross were sib-mated and females were then allowed to oviposit individually in order to establish F2 families of mosquitoes. F2 females were then allowed to feed on rhesus monkeys infected with P. cynomolgi Ceylon. These infections were carried out under the approved protocol 676-COL-MON-CH. Midguts were dissected from F2 female mosquitoes 8 – 10 days post-blood meal and stained with 2% mercurochrome (Sigma, St Louis, USA). The numbers of normal and encapsulated oocysts on each midgut were determined under a microscope 7.

Microsatellite Analysis

Genomic DNA extraction and genotyping of A. gambiae were carried out as described 1438. Most of the microsatellite marker alleles were scored using an ALFexpress DNA sequencer (Amersham Biosciences, Piscataway, NJ) with Cy-5 labeled primers from Sigma Genosys (The Woodlands, TX), Genset (La Jolla, CA), or Integrated DNA Technologies (Coralville, IA). Some markers were scored using ³²P-end-labeled primers as described 39. The previously described markers we employed were H145, H503, H036, H053, H711 and H678 for the X chromosome; H682, H784b, H417, H290, H175, H197, H187, H799, H095, H135, H770, H772, H603, H819, H802, H117, and R007 for chromosome 2;H776, H249, H525, H555, H158, H750, H127, H577, H758, H753, and H544 for chromosome 33840. New primer pairs designed for R007 and H249 are listed in Table 1. Additional new microsatellite markers and primer pairs designed for this study are also described in Table 1. With one exception these markers were identified among A. gambiae PEST strain bacterial artificial chromosome clones (BACs) from the Notre Dame-1 library (Collins et al., unpublished), referred to as ND-1 by 41, either from end sequences available at Genoscope http://www.cns.fr/externe/English/Projets/Resultats/rapport.html or by screening HindIII or EcoRI subclones of BACs (that have been mapped to the chromosome 3 by in situ hybridization) with ³²P-end labeled d(GT)₁₅and d(CT)₁₅. The Ddc microsatellite was deduced from G3 sequence, GenBank # AF063021. Genotypes of the F1 female parents and their F2 female progeny were determined.

Statistical Analyses

Both parental strains were polymorphic at many microsatellite loci, except for those on the X chromosome. This led to different numbers of alleles being found at different loci. Interval mapping 42 was not used for linkage mapping because the statistical analysis requires knowledge of the cis- or trans-relationships between alleles at different loci, and this could not always be determined. Encapsulation phenotypes were tested for statistical association with genotypes of the marker loci within families and for the combined informative data. A single locus association test was performed by fitting a linear regression between the trait value (fraction of oocysts encapsulated) and the genotype at each locus within a family. The significance of each F statistic was established by calculating the proportion of 1000 permutation samples exhibiting an F statistic as large as or larger than the F statistic determined from the observed data both before and after Bonferroni correction. The contribution of a given QTL to the genetic variance was estimated by R²43.

Authors' Contributions

MQB performed genetic crosses and infection studies. Additional genetic markers were developed by PR, CL and LZ. LZ and CL performed genotyping experiments. Statistical analyses were performed by SW, HZ. Writing was carried out primarily by PR and LZ.

Acknowledgements

We thank Frank Collins and Anton Cornel for provision of cytogenetically mapped BAC clones from the Notre Dame-1 library. We also thank JoAnn Sullivan for excellent technical support and Helen Rodd for additional statistical analysis. This work was supported by NIH grant AI43053 and an award from the Burroughs Wellcome Fund to L.Z. C.L. was supported in part by a grant (Project 990478) from UNPD/World Bank/WHO Special Programme for Research and Training in Tropical Diseases (TDR). Additional support was provided by an operating grant to P.R. from the Natural Sciences and Engineering Research Council of Canada.

Humoral encapsulation in insects Götz P Vey A Hemocytic and Humoral Immunity in Arthropods New York, Wiley and Sons AP Gupta 1986 407 429 Insect immunity Götz P Boman HG Comprehensive Insect Physiology, Biochemistry and Pharmacology Oxford, Pergamon Press GA Kerkut and LI Gilbert 1985 3 Mechanisms of encapsulation in dipteran hosts. Götz P Immune Mechanisms in Invertebrate Vectors Oxford, Oxford University Press AM Lackie 1986 pp.1 20 Immune mechanisms in insects Lackie AM Parasitol Today 1988 4 98 105 10.1016/0169-4758(88)90035-X On some peculiar pigmented cells found in two mosquitoes fed on malarial blood Ross R The British Medical Journal 1897 2 1786 1788 Genetics of Encapsulation Response of Anopheles gambiae against Plasmodium Parasites Zheng L Benedict MQ Romans P Paskewitz SM Collins FH Kafatos FC IX International Congress of Parasitology Monduzzi Editore, Bologna, Italy I Tada, S Kojima and M Tsuji 1998 103 109 Genetic selection of a Plasmodium-refractory strain of the malaria vector Anopheles gambiae Collins FH Sakai RK Vernick KD Paskewitz S Seeley DC Miller LH Collins WE Campbell CC Gwadz RW Science 1986 234 607 610 3532325 Ultrastructure of the encapsulation of Plasmodium cynomolgi (B strain) on the midgut of a refractory strain of Anopheles gambiae Paskewitz SM Brown MR Lea AO Collins FH J. Parasitol. 1988 74 432 439 3379524 Ultrastructural localizatioin of phenoloxidase in the midgut of refractory Anopheles gambiae and association of the enzyme with encapsulated Plasmodium cynomolgi Paskewitz SM Brown MR Collins FH Lea AO J. Parasitol. 1989 75 594 600 2503600 Immune responses of vectors Paskewitz SM Christensen BM The biology of disease vectors Niwot, Colorado, University Press of Colorado B J Beaty and W C Marquardt 1996 Association of a Plasmodium-refractory phenotype with an esterase locus in Anopheles gambiae Vernick KD Collins FH Am. J. Trop. Med. Hyg. 1989 40 593 597 2742037 A general system of resistance to malaria infection in Anopheles gambiae controlled by two main genetic loci Vernick KD Collins FH Gwadz RW Am. J. Trop. Med. Hyg. 1989 40 585 592 2742036 Association of two esterase genes, a chromosomal inversion, and susceptibility to Plasmodium cynomolgi in the African malaria vector Anopheles gambiae Crews-Oyen AE Kumar VK Collins FH Am. J. Trop. Med. Hyg. 1993 49 341 347 8372956 Quantitative trait loci for refractoriness of Anopheles gambiae to Plasmodium cynomolgi B Zheng L Cornel AJ Wang R Erfle H Voss H Ansorge W Kafatos FC Collins FH Science 1997 276 425 428 10.1126/science.276.5311.425 9103203 Mapping a quantitative trait locus involved in melanotic encapsulation of foreign bodies in the malaria vector, Anopheles gambiae Gorman MJ Severson DW Cornel AJ Collins FH Paskewitz SM Genetics 1997 146 965 971 9215900 Comparative genomic analysis in the region of a major Plasmodium-refractoriness locus of Anopheles gambiae Thomasova D Ton LQ Copley RR Zdobnov EM Wang X Hong YS Sim C Bork P Kafatos FC Collins FH Proc Natl Acad Sci U S A 2002 99 8179 8184 123041 12060762 10.1073/pnas.082235599 Progress in map-based cloning of Anopheles gambiae genes responsible for encapsulation of malaria parasites Collins FH Zheng L Paskewitz SM Kafatos FC Annals Trop. Med. Parasitol. 1997 91 517 521 10.1080/00034989760888 Genetics in the study of mosquito susceptibility to Plasmodium Collins FH Saunders RD Kafatos FC Roth C Ke Z Wang X Dymbrowski K Ton L Hogan J Parassitologia 1999 41 163 168 10697850 Genetic loci affecting resistance to human malaria parasites in a West African mosquito vector population Niare O Markianos K Volz J Oduol F Toure A Bagayoko M Sangare D Traore SF Wang R Blass C Dolo G Bouare M Kafatos FC Kruglyak L Toure YT Vernick KD Science 2002 298 213 216 10.1126/science.1073420 12364806 Two new malaria parasites, Plasmodium cynomolgi ceylonensis subsp. nov. and Plasmodium fragile sp. nov., from monkeys in Ceylon Dissanaike AS Nelson P Garnham PCC Ceylon J Med Sci (D) 1965 14 The circumsporozoite gene of the Plasmodium cynomolgi complex Galinski MR Arnot DE Cochrane AH Barnwell JW Nussenzweig RS Enea V Cell 1987 48 311 319 3802196 Evolutionary origin of human and primate malarias: evidence from the circumsporozoite protein gene Escalante AA Barrio E Ayala FJ Mol Biol & Evo 1995 12 616 626 A gene effecting malaria refractoriness is linked to Diphenol oxidase A2 on Chromosome 3 of Anopheles gambiae Romans P IV W. C. Black Sakai RK Gwadz RW Am. J. Trop. Med. Hyg. 1999 60 22 29 9988317 The genetics and expression of an esterase locus in Anopheles gambiae Vernick KD Collins FH Seeley DC Gwadz RW Miller LH Biochem. Genet. 1988 26 367 379 3214416 Leptopilina heterotoma and L. boulardi: strategies to avoid cellular defense responses of Drosophila melanogaster Rizki TM Rizki RM Carton Y Experimental Parasitology 1990 70 466 475 10.1016/0014-4894(90)90131-U 2108875 A new specific gene for wasp cellular immune resistance in Drosophila Benassi V Frey F Carton Y Heredity 1998 80 347 352 10.1038/sj.hdy.6883030 9569639 Drosophila Cellular immunity against parasitoids Carton Y Nappi AJ Parasitol. Today 1997 13 218 227 10.1016/S0169-4758(97)01058-2 Genetic basis of encapsulation response in Anopheles gambiae Zheng L Parassitologia 1999 41 181 184 10697853 Genetic localization of a Drosophila melanogaster resistance gene to a parasitoid wasp and physical mapping of the region Hita MT Poirie M Leblanc N Lemeunier F Lutcher F Frey F Periquet G Carton Y Genome Res 1999 9 471 481 10330127 Drosophila resistance genes to parasitoids: chromosomal location and linkage analysis Poirie M Frey F Hita M Huguet E Lemeunier F Periquet G Carton Y Proc R Soc Lond B Biol Sci 2000 267 1417 1421 10.1098/rspb.2000.1158 10983825 Melanization of plasmodium falciparum and C-25 sephadex beads by field-caught Anopheles gambiae (Diptera: Culicidae) from southern Tanzania Schwartz A Koella JC J Med Entomol 2002 39 84 88 11931276 Effect of mosquito age and reproductive status on melanization of sephadex beads in Plasmodium-refractory and -susceptible strains of Anopheles gambiae Chun J Riehle M Paskewitz SM J Invertebr Pathol 1995 66 11 17 10.1006/jipa.1995.1054 7544819 Response of Plasmodium refractory and susceptible strains of Anopheles gambiae to inoculated Sephadex beads Paskewitz SM Riehle MA Dev. Comp. Immunol. 1994 18 369 375 10.1016/0145-305X(94)90002-7 7535253 Tyrosinase-type prophenoloxidase distribution in the alimentary canal of strains of Anopheles gambiae refractory and susceptible to Plasmodium infection Brey PT Ahmed A Lee W-J Ashida M Lehane MJ Exp. Parasitol. 1995 80 654 664 10.1006/expr.1995.1081 7758546 Care and maintenance of anopheline mosquito colonies Benedict MQ The Molecular Biology of Insect Disease Vectors London, Chapman and Hall J M Crampton C B Beard C Louis 1997 3 12 Infection of mosquitoes with primate malaria Collins WE The Molecular Biology of Insect Disease Vectors London, Chapman and Hall J M Crampton C B Beard C Louis 1997 92 100 A genetic study of a melanization response to Sephadex beads in Plasmodium-refractory and -susceptible strains of Anopheles gambiae Gorman MJ Paskewitz SM Am J Trop Med Hyg 1997 56 Microsatellite markers and genotyping procedures for Anopheles gambiae, a malaria vector in Africa Wang R Kafatos FC Zheng L Parasitol. Today 1998 15 33 37 10.1016/S0169-4758(98)01360-X A detailed genetic map for the X chromosome of the malaria vector, Anopheles gambiae Zheng L Collins FH Kumar V Kafatos FC Science 1993 261 605 608 8342025 An integrated genetic map of the African human malaria vector mosquito, Anopheles gambiae Zheng L Benedict MQ Cornel AJ Collins FH Kafatos FC Genetics 1996 143 941 952 8725240 Construction of a BAC library and generation of BAC end sequence-tagged connectors for genome sequencing of the African malaria mosquito Anopheles gambiae Hong YS Hogan JR Wang X Sarkar A Sim C Loftus BJ Ren C Huff ER Carlile JL Black K Zhang HB Gardner MJ Collins FH Mol Genet Genomics 2003 268 720 728 12655398 Mapmaker: an interactive computer package for constructing primary genetic linkage maps of experimental and natural populations Lander E Green P Abrahamson J Barlow A Daley M Lincoln S Newburg L Genomics 1987 1 174 181 3692487 Mathematical statistics and data analysis Rice JA Belmont, CA, Duxbury Press (Wadsworth Publishing Company) 2nd 1995 Integrated genetic map of Anopheles gambiae: use of RAPD polymorphisms for genetic, cytogenetic and STS landmarks Dimopoulos G Zheng L Kumar V dellaTorre A Kafatos FC Louis C Genetics 1996 143 953 960 8725241 Molecular approaches to mosquito-parasite interactions Zheng L Arch. Insect Biochem. Physiol. 1996 34 1 18 10.1002/(SICI)1520-6327(1997)34:1<1::AID-ARCH1>3.3.CO;2-R When genetic distance matters: measuring genetic differentiation at microsatellite loci in whole-genome scans of recent and incipient mosquito species Wang R Zheng L Toure YT Dandekar T Kafatos FC Proc Natl Acad Sci U S A 2001 98 10769 10774 58550 11553812 10.1073/pnas.191003598