IHS/BFM/X/1

This page describes a signal of selection found in the Burkina Faso An. coluzzii populationusing the IHS (Cite et al. 20XX) statistic.The focus of this signal is on chromosome arm X between positions 14,760,000 and 15,620,000. The evidence supporting this signal is strong (\Delta_{i} >= 100 on both flanks).

This signal overlaps the Cyp9k1 locus, a genome region with prior evidence of an association with insecticide resistance and/or recent positive selection in Anopheles mosquitoes.

Signal location. Blue markers show the values of the selection statistic. The dashed black line shows the fitted peak model. The shaded red area shows the focus of the selection signal. The shaded blue area shows the genomic region in linkage with the selection event. Use the mouse wheel or the controls at the top right of the plot to zoom in, and hover over genes to see gene names and descriptions.

Genes

The following 41 genes overlap the focal region: AGAP000806 (Angiopoietin-like 1), AGAP000807 (helix-loop-helix transcription factor), AGAP000808 (DNA damage-regulated autophagy modulator protein 2), AGAP013022, AGAP000809 (Proteasome 26S non-ATPase subunit 10), AGAP000810, AGAP000812 (calcium binding protein), AGAP0008131 (Frataxin homolog, mitochondrial), AGAP000814, AGAP000815 (INTB - integrin beta subunit), AGAP000817, AGAP000816, AGAP013474, AGAP013173, AGAP013424, AGAP012997, AGAP0008181 (CYP9K1 - cytochrome P450), AGAP000819 (nuclear receptor subfamily 2 group E member (Tailless)), AGAP0008204 (CPR125 - cuticular protein RR-2 family 125), AGAP000821, AGAP000822, AGAP000823 (CD81 antigen), AGAP000824 (bone morphogenetic protein 5), AGAP000825, AGAP000826 (cap-specific mRNA (nucleoside-2’-O-)-methyltransferase 1), AGAP000829 (calpain-15), AGAP000830 (CASPS7 - short caspase 7), AGAP000831 (DnaJ homolog subfamily C member 25), AGAP000832 (Derlin-2/3), AGAP000833 (MIP - myoinhibitory-like peptide), AGAP000834, AGAP000835, AGAP028655, AGAP013506, AGAP000840 (amiloride-sensitive sodium channel, other), AGAP000841 (Ras-related protein Rab-40), AGAP013101, AGAP000842 (NADH dehydrogenase (ubiquinone) 1 alpha subcomplex, assembly factor 1), AGAP000843 (cardiolipin synthase), AGAP000844 (Progestin and adipoQ receptor family member 4), AGAP000847 (CDC-like kinase).

The following 4 genes are within 50 kbp of the focal region: AGAP0008012 (GLURIIb - ionotropic receptor GLURIIb), AGAP0008032 (GLURIIa - ionotropic receptor GLURIIa), AGAP0008041 (GPXH2 - glutathione peroxidase 2), AGAP000805 (BTB/POZ domain-containing protein KCTD16).

Key to insecticide resistance candidate gene types: 1 metabolic; 2 target-site; 3 behavioural; 4 cuticular.

Overlapping selection signals

The following selection signals have a focus which overlaps with the focus of this signal.

Signal Statistic Population Focus Peak model \Delta_{i} Max. percentile Known locus
H12/BFM/X/1 H12 Burkina Faso An. coluzzii X:15,100,000-15,300,000 1,267 100.0% Cyp9k1
XPEHH/BFS.UGS/X/1 XPEHH Burkina Faso An. gambiae X:14,960,000-15,140,000 961 99.9% nan
XPEHH/BFS.GWA/X/1 XPEHH Burkina Faso An. gambiae X:15,060,000-15,220,000 866 99.6% Cyp9k1
XPEHH/BFM.GWA/X/1 XPEHH Burkina Faso An. coluzzii X:15,060,000-15,260,000 852 99.4% Cyp9k1
H12/BFS/X/1 H12 Burkina Faso An. gambiae X:15,100,000-15,240,000 815 98.5% Cyp9k1
XPEHH/BFM.BFS/X/1 XPEHH Burkina Faso An. coluzzii X:15,580,000-15,800,000 718 99.9% nan
IHS/GNS/X/1 IHS Guinea An. gambiae X:15,240,000-15,780,000 553 99.3% Cyp9k1
IHS/UGS/X/1 IHS Uganda An. gambiae X:14,640,000-15,360,000 517 96.7% Cyp9k1
XPEHH/BFM.AOM/X/1 XPEHH Burkina Faso An. coluzzii X:15,040,000-15,180,000 511 100.0% nan
IHS/BFS/X/2 IHS Burkina Faso An. gambiae X:14,660,000-14,980,000 477 99.9% nan
H12/GNS/X/1 H12 Guinea An. gambiae X:14,960,000-15,160,000 419 97.8% nan
H12/UGS/X/1 H12 Uganda An. gambiae X:15,160,000-15,480,000 348 94.9% Cyp9k1
XPEHH/GAS.GWA/X/1 XPEHH Gabon An. gambiae X:14,460,000-14,860,000 301 98.5% nan
XPEHH/UGS.CMS/X/1 XPEHH Uganda An. gambiae X:15,240,000-15,460,000 284 99.9% Cyp9k1
XPEHH/AOM.GWA/X/2 XPEHH Angola An. coluzzii X:14,500,000-15,180,000 228 98.5% nan
XPEHH/UGS.BFS/X/1 XPEHH Uganda An. gambiae X:14,960,000-15,320,000 204 89.8% Cyp9k1
IHS/CMS/X/1 IHS Cameroon An. gambiae X:15,520,000-15,560,000 162 96.6% nan

Diagnostics

The information below provides some diagnostics from the Peak modelling algorithm.

Selection signal in context. @@TODO

Peak targetting. @@TODO

Peak fitting diagnostics. @@TODO

Model fit reports

Peak model:

[[Model]]
    Model(skewed_gaussian)
[[Fit Statistics]]
    # function evals   = 43
    # data points      = 314
    # variables        = 4
    chi-square         = 153.412
    reduced chi-square = 0.495
    Akaike info crit   = -216.908
    Bayesian info crit = -201.910
[[Variables]]
    center:      0 (fixed)
    amplitude:   4.81018313 +/- 0.149632 (3.11%) (init= 3)
    sigma:       2.04839171 +/- 0.071299 (3.48%) (init= 0.5)
    skew:       -0.22832452 +/- 0.026888 (11.78%) (init= 0)
    baseline:    2.58879421 +/- 0.067752 (2.62%) (init= 1)
    ceiling:     100 (fixed)
    floor:       0 (fixed)
[[Correlations]] (unreported correlations are <  0.100)
    C(sigma, baseline)           = -0.615
    C(amplitude, sigma)          = -0.450
    C(amplitude, baseline)       = -0.193

Null model:

[[Model]]
    Model(constant)
[[Fit Statistics]]
    # function evals   = 12
    # data points      = 313
    # variables        = 1
    chi-square         = 770.857
    reduced chi-square = 2.471
    Akaike info crit   = 284.107
    Bayesian info crit = 287.853
[[Variables]]
    c:   3.85278548 +/- 0.088845 (2.31%) (init= 1)

Comments