IHS/AOM/2/3

This page describes a signal of selection found in the Angola An. coluzzii populationusing the IHS (Cite et al. 20XX) statistic.The focus of this signal is on chromosome arm 2R between positions 33,940,000 and 34,440,000. The evidence supporting this signal is moderate (\Delta_{i} >= 50 on both flanks).

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 32 genes overlap the focal region: AGAP003212 (Adrenodoxin), AGAP003213 (nijmegen breakage syndrome protein 1), AGAP003214, AGAP003215 (26S proteasome regulatory subunit T2), AGAP003216 (26S proteasome regulatory subunit T2), AGAP003217 (mRpS6 - 28S ribosomal protein S6, mitochondrial), AGAP003218, AGAP003219, AGAP003220, AGAP003221 (ABCC3 - ATP-binding cassette transporter (ABC transporter) family C member 3), AGAP003222 (phosphatidylinositol glycan, class N), AGAP003224, AGAP003225, AGAP003226 (dynein heavy chain 2), AGAP003227 (nucleoside-diphosphate kinase), AGAP003228 (histone acetyltransferase type B subunit 2), AGAP003229, AGAP003230 (poly (ADP-ribose) polymerase), AGAP0032311 (NADH dehydrogenase (ubiquinone) iron-sulfur protein 7, mitochondrial), AGAP003232 (splicing factor 4), AGAP0133041, AGAP0032331, AGAP003235 (lachesin), AGAP003236, AGAP003237 (tubulin, alpha 1), AGAP003238 (N-myc downstream regulated), AGAP003239 (meiotic chromosome segregation protein), AGAP003240 (Protein jagunal), AGAP003241, AGAP003242 (RNA polymerase-associated protein LEO1), AGAP003243, AGAP003244 (GPRGHP2 - putative growth hormone releasing hormone receptor 2).

The following 15 genes are within 50 kbp of the focal region: AGAP003208 (solute carrier family 25, member 44), AGAP0032091 (C-4 methylsterol oxidase), AGAP003211, AGAP003245 (CLIPA19 - CLIP-domain serine protease), AGAP003246 (CLIPB2 - CLIP-domain serine protease), AGAP003247 (CLIPB19 - CLIP-domain serine protease), AGAP003248, AGAP013487, AGAP003249 (CLIPB3 - CLIP-domain serine protease), AGAP003250 (CLIPB4 - CLIP-domain serine protease), AGAP003251 (CLIPB1 - CLIP-domain serine protease), AGAP003252 (CLIPB6 - CLIP-domain serine protease), AGAP013184 (CLIPB36 - CLIP-domain serine protease), AGAP003253 (Gr15 - gustatory receptor 15), AGAP003254 (Gr16 - gustatory receptor 16).

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/AOM/2/7 H12 Angola An. coluzzii 2R:34,040,000-34,220,000 222 96.1% nan
H12/BFS/2/6 H12 Burkina Faso An. gambiae 2R:34,420,000-34,460,000 114 92.9% nan
XPEHH/AOM.BFM/2/7 XPEHH Angola An. coluzzii 2R:34,320,000-34,440,000 108 86.9% nan
XPEHH/BFS.UGS/2/5 XPEHH Burkina Faso An. gambiae 2R:34,420,000-34,480,000 95 96.7% 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   = 78
    # data points      = 373
    # variables        = 4
    chi-square         = 80.937
    reduced chi-square = 0.219
    Akaike info crit   = -561.909
    Bayesian info crit = -546.222
[[Variables]]
    center:      0 (fixed)
    amplitude:   1.48777673 +/- 0.090539 (6.09%) (init= 3)
    sigma:       1.02675162 +/- 0.082145 (8.00%) (init= 0.5)
    skew:       -0.84473510 +/- 0.091140 (10.79%) (init= 0)
    baseline:    2.37067481 +/- 0.031211 (1.32%) (init= 1)
    ceiling:     100 (fixed)
    floor:       0 (fixed)
[[Correlations]] (unreported correlations are <  0.100)
    C(sigma, skew)               =  0.481
    C(sigma, baseline)           = -0.411
    C(amplitude, sigma)          = -0.351
    C(amplitude, baseline)       = -0.229

Null model:

[[Model]]
    Model(constant)
[[Fit Statistics]]
    # function evals   = 11
    # data points      = 372
    # variables        = 1
    chi-square         = 145.014
    reduced chi-square = 0.391
    Akaike info crit   = -348.448
    Bayesian info crit = -344.529
[[Variables]]
    c:   2.61272388 +/- 0.032414 (1.24%) (init= 1)

Comments