IHS/BFS/2/3

This page describes a signal of selection found in the Burkina Faso An. gambiae populationusing the IHS (Cite et al. 20XX) statistic.The focus of this signal is on chromosome arm 2R between positions 40,320,000 and 40,560,000. The evidence supporting this signal is strong (\Delta_{i} >= 100 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 43 genes overlap the focal region: AGAP0035811 (D-xylulose reductase A), AGAP0035821 (D-xylulose reductase A), AGAP0035831 (L-iditol 2-dehydrogenase), AGAP0035841 (L-iditol 2-dehydrogenase), AGAP003585, AGAP003586 (Phosphate carrier, mitochondrial), AGAP013333 (Phosphate carrier, mitochondrial), AGAP003587 (solute carrier family 26 (sodium-independent sulfate anion transporter)), AGAP003588 (VATH - V-type H -transporting ATPase subunit H), AGAP003589 (prolyl-tRNA synthetase), AGAP003590, AGAP003591 (LRR-repeat protein 1), AGAP003592 (RpLP0 - 60S ribosomal protein LP0), AGAP003593 (Protein arginine methyltransferase RmtB), AGAP003594, AGAP003595 (rRNA biogenesis protein RRP5), AGAP003596, AGAP003597 (histone-lysine N-methyltransferase SUV39H), AGAP003598 (template-activating factor I), AGAP003599 (F-type H -transporting ATPase oligomycin sensitivity conferral protein), AGAP003600 (Elongation of very long chain fatty acids protein), AGAP003601 (alpha-1,2-mannosyltransferase), AGAP003602 (methyltransferase-like protein 9), AGAP003603, AGAP003604, AGAP003605 (cullin 2), AGAP003606 (Fatty acyl-CoA reductase), AGAP003607 (Cad96Cb, isoform B), AGAP0036081 (CYP4AA1 - cytochrome P450), AGAP003610, AGAP003611 (Fatty acyl-CoA reductase), AGAP003612, AGAP003613 (palmitoyltransferase ZDHHC3), AGAP003614, AGAP013196, AGAP013432, AGAP029070, AGAP029099, AGAP029064, AGAP029085, AGAP029083, AGAP003615 (toll-interacting protein), AGAP003616.

The following 7 genes are within 50 kbp of the focal region: AGAP003580, AGAP003618 (serine/threonine-protein kinase LATS1/2), AGAP003619 (transcription elongation factor SPT6), AGAP003620, AGAP003621 (coiled-coil domain-containing protein 25), AGAP003622 (Queuine tRNA-ribosyltransferase catalytic subunit 1), AGAP003623 (long-chain acyl-CoA synthetase).

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
IHS/CMS/2/2 IHS Cameroon An. gambiae 2R:40,060,000-40,780,000 643 100.0% nan
H12/CMS/2/3 H12 Cameroon An. gambiae 2R:40,300,000-40,740,000 501 98.7% nan
IHS/GNS/2/3 IHS Guinea An. gambiae 2R:40,200,000-40,420,000 450 99.0% 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_exponential_peak)
[[Fit Statistics]]
    # function evals   = 61
    # data points      = 543
    # variables        = 4
    chi-square         = 183.521
    reduced chi-square = 0.340
    Akaike info crit   = -581.035
    Bayesian info crit = -563.846
[[Variables]]
    center:      0 (fixed)
    amplitude:   2.73190578 +/- 0.088779 (3.25%) (init= 3)
    decay:       2.99999999 +/- 0.001870 (0.06%) (init= 0.5)
    skew:       -0.90857473 +/- 0.051831 (5.70%) (init= 0)
    baseline:    2.05968896 +/- 0.064085 (3.11%) (init= 1)
    ceiling:     100 (fixed)
    floor:       0 (fixed)
[[Correlations]] (unreported correlations are <  0.100)
    C(decay, baseline)           =  0.769
    C(amplitude, baseline)       = -0.509
    C(skew, baseline)            = -0.276
    C(amplitude, skew)           =  0.239
    C(decay, skew)               = -0.227

Null model:

[[Model]]
    Model(constant)
[[Fit Statistics]]
    # function evals   = 11
    # data points      = 542
    # variables        = 1
    chi-square         = 516.663
    reduced chi-square = 0.955
    Akaike info crit   = -23.948
    Bayesian info crit = -19.653
[[Variables]]
    c:   3.06005740 +/- 0.041974 (1.37%) (init= 1)

Comments