A total of 59 contigs with A. darlingi transcriptome SSRs were identified (34 dinucleotides, 23 trinucleotides, 1 tetranucleotide and 1 hexanucleotide). Thirty-seven pairs of primers were designed. The type of repetition of most of the SSRs was perfect simple and the most frequent repeat motif was dinucleotide (AC/AG).
The characterization of 40 A. darlingi in the 37 pairs of primers generated the amplification of 28 SSR loci of which 15 were polymorphic and amplified at 56 ºC to 65 ºC with size ranged from 115 to 368 bp (Table 1).
In all the 15 polymorphic SSR loci the dinucleotide repeat motif (AC/AG) was the most frequent. This motif was also found in A. darlingi from Coari, Amazonas state (Brazil) [12], in A. gambiae from Africa [21], in A. culicifacies from India [22], and in A. triannulatus from Manaus, Amazonas state, Brazil) [8]. However, in A. darlingi from Capanema, Pará state (Brazil) [10], the repeat motif that stood out was the (TC/GA). This indicates that dinucleotide is the most frequent type of SSR, and this is irrespective of the species.
The 15 SSR loci had a total of 76 alleles, ranging from 2 to 9, with a mean of 5.0 alleles per locus; variation of observed heterozygosity was 0.026–0.769 with a mean HO of 0.402 and an expected HE of 0.025–0.776, with a mean HE of 0.509 (Table 2).
Considering the genetic variability, the allelic frequency found in A. marajoara (11 to 52) [11], in A. albitarsis (2 to 10) [7], in A. triannulatus (3 to 10) [8], in A. stephensi (3 to 16) [4], in A. darlingi (7 to 44) [10], (4 to 11) [12] and from 2 to 9 (Table 2) in this study, indicated that SSRs are efficient for the investigation of the genetic structure of populations due to their high polymorphism.
The variations of observed (HO = 0.026–0.769) and expected heterozygosity (HE = 0.025–0.776) found in this study were similar (Table 2) when compared to A. culicifacies (HO = 0.045–0.607 and HE = 0.170–0.841) [22], A. gambiae (HO = 0.216–0.894) [21], and A. nuneztovari (HO = 0.354–0.866 and HE = 0.613–0.932) [9]. Comparing the heterozygosity of A. darlingi in this study with those of other populations of the same species, such as those from Capanema, in the Pará state (HO = 0.368–0.769 and HE = 0.785–0.987) [10] and Coari, Amazonas state (HO = 0.037–0.833 and HE = 0.177–0.871) [12], higher values are noted.
Based on the values of the PIC scale [23], eight highly informative loci (Adar4. Adar11, Adar13, Adar16, Adar24, Adar30, Adar34, Adar36) and five moderately informative loci (Adar1, Adar3, Adar6, Adar19, Adar32) were observed. Thus, new loci with this performance are made available for future studies of the species.
Eight loci presented significant HWE after Bonferroni correction (P: (5%) ≤ 0.0005) [24]. The HWE deviation present in the other loci may have been caused by the presence of null alleles, the bottleneck effect, or due to mating between related individuals. No linkage disequilibrium was found among the loci (Table 2).
HWE was found in all SSR loci of A. moucheti with no null alleles in 3 loci of A. maculipennis, which presented some null alleles due to the deficiency in heterozygotes [13], and in 2 loci of A. triannulatus, with the presence of null alleles due to excess homozygotes [8]. On the other hand, HWE was found in 11 loci of A. albitarsis [7], in 6 loci of A. nuneztovari [9], in 10 loci of A. darlingi [12] and in 7 short tandem repeat (STR) loci of A. darlingi [5], which is similar to what was found in this study. This is due to the presence of null alleles and/or a bottleneck effect.
In addition, ESTs such as cuticular protein 23 RR-1 family, troponin C, ribosomal genes, and myosin light chain 2 from a cDNA library of A. darlingi (Rafael, personal communication; [13]), were submitted to the basic local alignment search tool (BLAST) against A. gambiae, the primary malaria vector in Africa, via an in silico hybridization method. These gene sequences showed regions of similarity in the 2L, 2R, 3L and X chromosome arms between both species, as well as with ESTs from other insects such as Aedes aegypti, Drosophila melanogaster and Culex spp (L. Bridi, personal communication). Furthermore, troponin (L. Bridi, personal communication) and actin [25] were hybridized as probes via a fluorescent in situ hybridization method in the polytene chromosomes of A. darlingi, which showed markers located in the arm 2R, section 14B, and in the arm 2L, section 23, respectively. These results provide a suitable basis for establishing chromosomal and evolutionary genomics homologies among major malaria vectors.
The present work provides 15 new microsatellite markers for A. darlingi, totaling 64 SSRs, when considering the studies that developed and characterized SSR loci, which are the first EST-SSR markers.
Considering the great capacity of A. darlingi for the transmission of malaria and its resistance to synthetic insecticide, the EST-SSR markers characterized are potential tools in future populational studies of this mosquito and actions for its control, especially in the Amazon.