Cis-acting mutation and duplication: History of molecular evolution in a P450 haplotype responsible for insecticide resistance in Culex quinquefasciatus

A cytochrome P450 gene, Cyp9m10, is more than 200-fold overexpressed in a pyrethroid resistant strain of Culex quinquefasciatus, JPal-per. The haplotype of this strain contains two copies of Cyp9m10 resulted from recent tandem duplication. In this study, we discovered and isolated a Cyp9m10 haplotype closely related to this duplicated Cyp9m10 haplotype from JHB, a strain used for the recent genome project for this mosquito species. The isolated haplotype (JHB-NIID-B haplotype) shared the same insertion of a transposable element upstream of the coding region with JPal-per strain but not duplicated. The JHB-NIID-B haplotype was considered to have diverged from the JPal-per lineage just before the duplication event. Cyp9m10 was moderately overexpressed in larvae with the JHB-NIID-B haplotype. The overexpressions in JHB-NIID-B and JPal-per haplotypes were developmentally regulated in similar pattern indicating both haplotypes share a common cis-acting mutation responsible for the overexpressions. The isolated moderately overexpressed haplotype conferred resistance, however, its efficacy was relatively small. We hypothesized that the first cis-acting mutation modified the consequence of the subsequent duplication in JPal-per lineage to confer stronger phenotypic effect than that if it occurred before the first cis-acting mutation.     

Allelic Variation of Cytochrome P450s Drives Resistance to Bednet Insecticides in a Major Malaria Vector

Scale up of Long Lasting Insecticide Nets (LLINs) has massively contributed to reduce malaria mortality across Africa. However, resistance to pyrethroid insecticides in malaria vectors threatens its continued effectiveness. Deciphering the detailed molecular basis of such resistance and designing diagnostic tools is critical to implement suitable resistance management strategies. Here, we demonstrated that allelic variation in two cytochrome P450 genes is the most important driver of pyrethroid resistance in the major African malaria vector Anopheles funestus and detected key mutations controlling this resistance. An Africa-wide polymorphism analysis of the duplicated genes CYP6P9a and CYP6P9b revealed that both genes are directionally selected with alleles segregating according to resistance phenotypes. Modelling and docking simulations predicted that resistant alleles were better metabolizers of pyrethroids than susceptible alleles. Metabolism assays performed with recombinant enzymes of various alleles confirmed that alleles from resistant mosquitoes had significantly higher activities toward pyrethroids. Additionally, transgenic expression in Drosophila showed that flies expressing resistant alleles of both genes were significantly more resistant to pyrethroids compared with those expressing the susceptible alleles, indicating that allelic variation is the key resistance mechanism. Furthermore, site-directed mutagenesis and functional analyses demonstrated that three amino acid changes (Val¹⁰⁹Ile, Asp³³⁵Glu and Asn³⁸⁴Ser) from the resistant allele of CYP6P9b were key pyrethroid resistance mutations inducing high metabolic efficiency. The detection of these first DNA markers of metabolic resistance to pyrethroids allows the design of DNA-based diagnostic tools to detect and track resistance associated with bednets scale up, which will improve the design of evidence-based resistance management strategies.     

Mitochondrial DNA analyses reveal low genetic diversity in Culex quinquefasciatus from residential areas in Malaysia

The present study explored the intraspecific genetic diversity, dispersal patterns and phylogeographic relationships of Culex quinquefasciatus Say (Diptera: Culicidae) in Malaysia using reference data available in GenBank in order to reveal this species' phylogenetic relationships. A statistical parsimony network of 70 taxa aligned as 624 characters of the cytochrome c oxidase subunit I (COI) gene and 685 characters of the cytochrome c oxidase subunit II (COII) gene revealed three haplotypes (A1–A3) and four haplotypes (B1–B4), respectively. The concatenated sequences of both COI and COII genes with a total of 1309 characters revealed seven haplotypes (AB1–AB7). Analysis using tcs indicated that haplotype AB1 was the common ancestor and the most widespread haplotype in Malaysia. The genetic distance based on concatenated sequences of both COI and COII genes ranged from 0.00076 to 0.00229. Sequence alignment of Cx. quinquefasciatus from Malaysia and other countries revealed four haplotypes (AA1–AA4) by the COI gene and nine haplotypes (BB1–BB9) by the COII gene. Phylogenetic analyses demonstrated that Malaysian Cx. quinquefasciatus share the same genetic lineage as East African and Asian Cx. quinquefasciatus. This study has inferred the genetic lineages, dispersal patterns and hypothetical ancestral genotypes of Cx. quinquefasciatus.     

Allele frequency distribution of CYP2C9*2 and CYP2C9*3 polymorphisms in six Mexican populations

Allele frequency differences of functional CYP2C9 polymorphisms are responsible for some of the variation in drug response observed in human populations. The most relevant CYP2C9 functional variants are CYP2C9*2 (rs1799853) and CYP2C9*3 (rs1057910). These polymorphisms show variation in allele frequencies among different population groups. The present study aimed to analyze these polymorphisms in 947 Mexican-Mestizo from Mexico City and 483 individuals from five indigenous Mexican populations: Nahua, Teenek, Tarahumara, Purepecha and Huichol. The CYP2C9*2 allele frequencies in the Mestizo, Nahua and Teenek populations were 0.051, 0.007 and 0.005, respectively. As for CYP2C9*3, the allelic frequencies in the Mestizo, Nahua and Teenek populations were 0.04, 0.005 and 0.005, respectively. The CYP2C9*2 and CYP2C9*3 alleles were not observed in the Tarahumara, Purepecha and Huichol populations. These findings are in agreement with previous studies reporting very low allele frequencies for these polymorphisms in American Indigenous populations.     

Permethrin resistance variation and susceptible reference line isolation in a field population of the mosquito,Culex quinquefasciatus (Diptera: Culicidae)

This study examines the genetic variations and mechanisms involved in the development of permethrin resistance in individual mosquitoes from a field population of Culex quinquefasciatus, HAmCq^G0, and characterizes susceptible reference lines of mosquitoes with a similar genetic background to the field HAmCq^G0 strain. Six upregulated cytochrome P450 genes, CYP9M10, CYP9J34, CYP6P14, CYP9J40, CYP6AA7, and CYP4C52v1, previously identified as being upregulated in the larvae of resistant HAmCq^G8 mosquitoes were examined in the larvae of 3 strains (susceptible S‐Lab, parental HAmCq^G0 and permethrin‐selected highly resistant HAmCq^G8) and 8 HAmCq^G0 single‐egg raft colonies, covering a range of levels of susceptibility/resistance to permethrin and exhibiting different variations in the expression of A and/or T alleles at the L‐to‐F kdr locus of the sodium channel. The 2 lines with the lowest tolerance to permethrin and bearing solely the susceptible A allele at the L‐to‐F kdr locus of the sodium channels, from colonies Cx_SERC5 and Cx_SERC8, showed lower or similar levels of all 6 of the P450 genes tested compared with the S‐Lab strain, suggesting that these 2 lines could be used as the reference mosquitoes in future studies characterizing insecticide resistance in HAmCq mosquitoes. This study also provides a detailed investigation of the mechanisms involved in insecticide resistance in individuals within a population: individuals with elevated levels of resistance to permethrin all displayed one or more potential resistance mechanisms–either elevated levels of P450 gene expression, or L‐to‐F mutations in the sodium channel, or both.     

Expression variation of the porcine ADRB2 has a complex genetic background

Porcine adrenergic receptor beta 2 (ADRB2) gene exhibits differential allelic expression in skeletal muscle, and its genetic variation has been associated with muscle pH. Exploring the molecular–genetic background of expression variation for porcine ADRB2 will provide insight into the mechanisms driving its regulatory divergence and may also contribute to unraveling the genetic basis of muscle-related traits in pigs. In the present study, we therefore examined haplotype effects on the expression of porcine ADRB2 in four tissues: longissimus dorsi muscle, liver, subcutaneous fat, and spleen. The diversity and structure of haplotypes of the proximal gene region segregating in German commercial breeds were characterized. Seven haplotypes falling into three clades were identified. Two clades including five haplotypes most likely originated from introgression of Asian genetics during formation of modern breeds. Expression analyses revealed that the Asian-derived haplotypes increase expression of the porcine ADRB2 compared to the major, wild-type haplotype independently of tissue type. In addition, several tissue-specific differences in the expression of the Asian-derived haplotypes were found. Inspection of haplotype sequences showed that differentially expressed haplotypes exhibit polymorphisms in a polyguanine tract located in the core promoter region. These findings demonstrate that expression variation of the porcine ADRB2 has a complex genetic basis and suggest that the promoter polyguanine tract is causally involved. This study highlights the challenges of finding causal genetic variants underlying complex traits.     

Development of permethrin resistance in Culex quinquefasciatus Say in Kuala Lumpur,Malaysia

The resistance status towards permethrin among the laboratory strain, the permethrin-selected strain and four field strains of Culex quinquefasciatus collected in Kuala Lumpur, Malaysia was determined using three standard laboratory methods: WHO larval bioassay, WHO adult bioassay and biochemical microplate assay. Cx. quinquefasciatus permethrin-selected strain larvae were the least susceptible to permethrin with a resistance ratio of 47.28-folds, whereas all field strain larvae of the same species were tolerant to permethrin with resistance ratios of more than 3-folds. In contrast, in adult stage, the permethrin exposed permethrin-selected strain (resistance ratio = 1.27) was found to be more susceptible to permethrin than all permethrin-exposed field strains (resistance ratios = 2.23–2.48). Complete mortalities for all strains of Cx. quinquefasciatus adults proved the effectiveness of the synergist; piperonyl butoxide (PBO). For the biochemical microplate assay, the reduction of the mean optical density of elevated oxidase activity of three field strains upon exposure to PBO confirmed the association between oxidase activity and permethrin tolerance. On the other hand, irregular patterns of the mean optical density of elevated oxidase activity in the laboratory strain, permethrin-selected strain and Jalan Fletcher strain illustrated the gene variation within these mosquito colonies as well as the involvement of other enzyme activities in the permethrin resistance occurred.     

Single nucleotide deletion of cqm1 gene results in the development of resistance to Bacillus sphaericus in Culex quinquefasciatus

The entomopathogen Bacillus sphaericus is one of the most effective biolarvicides used to control the Culex species of mosquito. The appearance of resistance in mosquitoes to this bacterium, however, remains a threat to its continuous use in integrated mosquito control programs. Previous work showed that the resistance to B. sphaericus in Culex colonies was associated with the absence of the 60-kDa binary toxin receptor (Cpm1/Cqm1), an alpha-glucosidase present in the larval midgut microvilli. In this work, we studied the molecular basis of the resistance developed by Culex quinquefasciatus to B. sphaericus C3-41. The cqm1 genes were cloned from susceptible (CqSL) and resistant (CqRL/C3-41) colonies, respectively. The sequence of the cDNA and genomic DNA derived from CqRL/C3-41 colony differed from that of CqSL one by a one-nucleotide deletion which resulted in a premature stop codon, leading to production of a truncated protein. Recombinant Cqm1S from the CqSL colony expressed in Escherichia coli specifically bound to the Bin toxin and had α-glucosidase activity, whereas the Cqm1R from the CqRL/C3-41 colony, with a deletion of three quarters of the receptor’s C-terminal lost its α-glucosidase activity and could not bind to the binary toxin. Immunoblotting experiments showed that Cqm1 was undetectable in CqRL/C3-41 larvae, although the gene was correctly transcribed. Thus, the cqm1R represents a new allele in C. quinquefasciatus that confers resistance to B. sphaericus.     

Very low mitochondrial variability in a stingless bee endemic to cerrado

Partamona mulata is a stingless bee species endemic to cerrado, a severely threatened phytogeographical domain. Clearing for pasture without proper soil treatment in the cerrado facilitates the proliferation of termite ground nests, which are the nesting sites for P. mulata. The genetic consequences of these changes in the cerrado environment for bee populations are still understudied. In this work, we analyzed the genetic diversity of 48 colonies of P. mulata collected throughout the species’ distribution range by sequencing two mitochondrial genes, cytochrome oxidase I and cytochrome B. A very low polymorphism rate was observed when compared to another Partamona species from the Atlantic forest. Exclusive haplotypes were observed in two of the five areas sampled. The sharing of two haplotypes between collection sites separated by a distance greater than the flight range of queens indicates an ancient distribution for these haplotypes. The low haplotype and nucleotide diversity observed here suggests that P. mulata is either a young species or one that has been through population bottlenecks. Locally predominant and exclusive haplotypes (H2 and H4) may have been derived from local remnants through cerrado deforestation and the expansion of a few colonies with abundant nesting sites.     

TP53 PIN3 and PEX4 polymorphisms and infertility associated with endometriosis or with post-in vitro fertilization implantation failure

p53 has a crucial role in human fertility by regulating the expression of leukemia inhibitory factor (LIF), a secreted cytokine critical for blastocyst implantation. To examine whether TP53 polymorphisms may be involved with in vitro fertilization (IVF) failure and endometriosis (END), we have assessed the associations between TP53 polymorphism in intron 2 (PIN2; G/C, intron 2), PIN3 (one (N, non-duplicated) or two (D, duplicated) repeats of a 16-bp motif, intron 3) and polymorphism in exon 4 (PEX4; C/G, p.P72R, exon 4) in 98 women with END and 115 women with post-IVF failure. In addition, 134 fertile women and 300 women unselected with respect to fertility-related features were assessed. TP53 polymorphisms and haplotypes were identified by amplification refractory mutation system polymerase chain reaction. TP53 PIN3 and PEX4 were associated with both END (P=0.042 and P=0.007, respectively) and IVF (P=0.004 and P=0.009, respectively) when compared with women both selected and unselected for fertility-related features. Haplotypes D-C and N-C were related to higher risk for END (P=0.002, P=0.001, respectively) and failure of IVF (P=0.018 and P=0.002, respectively) when compared with the Fertile group. These results support that specific TP53 haplotypes are associated with an increased risk of END-associated infertility and with post-IVF failure.     

Identification of the amino acids in the Major Histocompatibility Complex class II region of Scottish Blackface sheep that are associated with resistance to nematode infection

Lambs with the Major Histocompatibility Complex DRB1*1101 allele have been shown to produce fewer nematode eggs following natural and deliberate infection. These sheep also possess fewer adult Teladorsagia circumcincta than sheep with alternative alleles at the DRB1 locus. However, it is unclear if this allele is responsible for the reduced egg counts or merely acts as a marker for a linked gene. This study defined the MHC haplotypes in a population of naturally infected Scottish Blackface sheep by PCR amplification and sequencing, and examined the associations between MHC haplotypes and faecal egg counts by generalised linear mixed modelling. The DRB1*1101 allele occurred predominately on one haplotype and a comparison of haplotypes indicated that the causal mutation or mutations occurred in or around this locus. Additional comparisons with another resistant haplotype indicated that mutations in or around the DQB2*GU191460 allele were also responsible for resistance to nematode infections. Further analyses identified six amino acid substitutions in the antigen binding site of DRB1*1101 that were significantly associated with reductions in the numbers of adult T. circumcincta.     

Circadian control of permethrin-resistance in the mosquito Aedes aegypti

Daily fluctuation of permethrin-resistance was found in adult mosquito Aedes aegypti, the major vector of dengue viruses in Taiwan. We hypothesized there is a relationship between resistance and the circadian clock. To test our hypothesis we correlated changes in the knock-down time (KT₅₀) response to permethrin with the expression of the pyrethroid-resistant gene CYP9M9 and the clock gene period (per) during a 12:12 h photoperiodic cycle. Rhythmic expression of per peaked at early scotophase of the light-dark cycle and at early subjective night in constant darkness. The values of KT₅₀ and the expression of CYP9M9 also exhibited circadian rhythms in both susceptible and permethrin-resistant mosquito strains, from which we inferred a link to the circadian clock. The KT₅₀ was significantly longer in the light than in the dark phase, and the level of CYP9M9 mRNA was maximal in early scotophase, dropped to a minimum in the midnight and then slowly increased through the photophase. Existence of a clock control over mosquito sensitivity to permethrin was further indicated by reduced expression of CYP9M9 and reduced mosquito resistance to permethrin after temporal silencing of the per gene. These data provide the first evidence on the circadian control of insect resistance to permethrin.     

Population structure of North African honey bees is influenced by both biological and anthropogenic factors

Honey bee diversity is under threat due to anthropogenic factors as the use of pesticides and the replacement of local colonies to recover from colony losses. To assess the effect of these activities on the genetic diversity and structure of North African honey bee colonies, we studied colonies from the north (Tellian) and the south (Saharan) regions in Algeria, by determining their mitochondrial haplotype and the variation of ten microsatellite loci. Particular haplotypes have been found with a high frequency in each region that may constitute subspecies-specific markers for Apis mellifera intermissa (haplotype A8 at the north) and A. m. sahariensis (haplotype A9 at the south). Moreover, the presence of the haplotype A8 in some Saharan colonies may reflect recent introductions of A. m. intermissa. Structure analysis suggests that a natural differentiation between honey bee populations from Saharan and Tellian regions still exists despite increased colony movements (migration, queen purchases, etc.) during the last decades. One apiary established for the conservation of A. m. sahariensis showed no indication of maternal introgression since all the colonies bear the same haplotype A9. Furthermore, Hardy–Weinberg equilibrium observed in this population indicates that this apiary is appropriate for conservation programs of A. m. sahariensis.     

KIR diversity in Māori and Polynesians: populations in which HLA-B is not a significant KIR ligand

HLA class I molecules and killer cell immunoglobulin-like receptors (KIR) form a diverse system of ligands and receptors that individualize human immune systems in ways that improve the survival of individuals and populations. Human settlement of Oceania by island-hopping East and Southeast Asian migrants started ~3,500 years ago. Subsequently, New Zealand was reached ~750 years ago by ancestral Māori. To examine how this history impacted KIR and HLA diversity, and their functional interaction, we defined at high resolution the allelic and haplotype diversity of the 13 expressed KIR genes in 49 Māori and 34 Polynesians. Eighty KIR variants, including four ‘new’ alleles, were defined, as were 35 centromeric and 22 telomeric KIR region haplotypes, which combine to give >50 full-length KIR haplotypes. Two new and divergent variant KIR form part of a telomeric KIR haplotype, which appears derived from Papua New Guinea and was probably obtained by the Asian migrants en route to Polynesia. Māori and Polynesian KIR are very similar, but differ significantly from African, European, Japanese, and Amerindian KIR. Māori and Polynesians have high KIR haplotype diversity with corresponding allotype diversity being maintained throughout the KIR locus. Within the population, each individual has a unique combination of HLA class I and KIR. Characterizing Māori and Polynesians is a paucity of HLA-B allotypes recognized by KIR. Compensating for this deficiency are high frequencies (>50 %) of HLA-A allotypes recognized by KIR. These HLA-A allotypes are ones that modern humans likely acquired from archaic humans at a much earlier time.     

Association between MT-CO3 haplotypes and high-altitude adaptation in Tibetan chicken

Genetic mutation in cytochrome c oxidase subunit III gene (MT-CO3) could influence the kinetics of cytochrome c oxidase (COX), which catalyzes oxygen transport capacity in oxidative phosphorylation. However, the potential relationship between MT-CO3 variants and high-altitude adaptation remains poorly understood in Tibetan chicken. Here, we sequenced MT-CO3 gene of 125 Tibetan chickens and 144 Chinese domestic chickens in areas at a low elevation (below 1000 m). Eight single nucleotide polymorphisms (SNPs) were detected; and five of them (m.10081A>G, m.10115G>A, m.10270G>A, m.10336A>G and m.10447C>T) shared by Tibetan chicken and lowland chicken with the significant difference in their respective allele frequencies. Nine haplotypes (H1–H9) were finally defined. Among them, haplotype H4 was positively associated with high-altitude adaptation whereas haplotypes H6, H7 and H8 had negative association with high-altitude adaptation. The Median-joining profile suggested that haplotype H5 had the ancestral position to the other haplotypes but had no significant relationship with high-altitude adaptation. However, there was only m.10081A>G mutation differed from haplotype H4 and H5. Results also suggested that chickens with A allele at m.10081A>G, had over 2.6 times than those with G allele in the probability of the ability to adapt hypoxia. It suggests that the synonymous mutation m.10081A>G may be a prerequisite for shaping high-altitude adaptation-specific haplotypes.     

CYP6 P450 Enzymes and ACE-1 Duplication Produce Extreme and Multiple Insecticide Resistance in the Malaria Mosquito Anopheles gambiae

Malaria control relies heavily on pyrethroid insecticides, to which susceptibility is declining in Anopheles mosquitoes. To combat pyrethroid resistance, application of alternative insecticides is advocated for indoor residual spraying (IRS), and carbamates are increasingly important. Emergence of a very strong carbamate resistance phenotype in Anopheles gambiae from Tiassalé, Côte d''Ivoire, West Africa, is therefore a potentially major operational challenge, particularly because these malaria vectors now exhibit resistance to multiple insecticide classes. We investigated the genetic basis of resistance to the most commonly-applied carbamate, bendiocarb, in An. gambiae from Tiassalé. Geographically-replicated whole genome microarray experiments identified elevated P450 enzyme expression as associated with bendiocarb resistance, most notably genes from the CYP6 subfamily. P450s were further implicated in resistance phenotypes by induction of significantly elevated mortality to bendiocarb by the synergist piperonyl butoxide (PBO), which also enhanced the action of pyrethroids and an organophosphate. CYP6P3 and especially CYP6M2 produced bendiocarb resistance via transgenic expression in Drosophila in addition to pyrethroid resistance for both genes, and DDT resistance for CYP6M2 expression. CYP6M2 can thus cause resistance to three distinct classes of insecticide although the biochemical mechanism for carbamates is unclear because, in contrast to CYP6P3, recombinant CYP6M2 did not metabolise bendiocarb in vitro. Strongly bendiocarb resistant mosquitoes also displayed elevated expression of the acetylcholinesterase ACE-1 gene, arising at least in part from gene duplication, which confers a survival advantage to carriers of additional copies of resistant ACE-1 G119S alleles. Our results are alarming for vector-based malaria control. Extreme carbamate resistance in Tiassalé An. gambiae results from coupling of over-expressed target site allelic variants with heightened CYP6 P450 expression, which also provides resistance across contrasting insecticides. Mosquito populations displaying such a diverse basis of extreme and cross-resistance are likely to be unresponsive to standard insecticide resistance management practices.