This page describes a signal of selection found in the
Burkina Faso An. gambiae population using the
H12 (Garud et al. 20XX) statistic.The focus of this signal is on chromosome arm
2R between positions 28,440,001 and
28,520,000.
The evidence supporting this signal is
strong ( >= 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.
The following 12 genes overlap the focal region: AGAP002859 (solute carrier family 8 (sodium/calcium exchanger)), AGAP0028621 (CYP6AA1 - cytochrome P450), AGAP0131281 (CYP6AA2 - cytochrome P450), AGAP0028631 (COEAE6O - carboxylesterase alpha esterase), AGAP0028641 (CYP6P15P - cytochrome P450), AGAP0028651 (CYP6P3 - cytochrome P450), AGAP0028661 (CYP6P5 - cytochrome P450), AGAP0028671 (CYP6P4 - cytochrome P450), AGAP0028681 (CYP6P1 - cytochrome P450), AGAP0028691 (CYP6P2 - cytochrome P450), AGAP0028701 (CYP6AD1 - cytochrome P450), AGAP013202.
The following 6 genes are within 50 kbp of the focal region: AGAP000586, AGAP002872 (lipase), AGAP002873, AGAP013069, AGAP002874, AGAP002875 (protein HEXIM1/2).
Key to insecticide resistance candidate gene types: 1 metabolic; 2 target-site; 3 behavioural; 4 cuticular.
The following selection signals have a focus which overlaps with the focus of this signal.
Signal | Focus | Score (Left | Right) |
---|---|---|
H12 / Uganda An. gambiae / Chromosome 2 / #1 | 2R:28,420,001-28,520,000 | 1230 (582 | 647) |
H12 / Guinea An. gambiae / Chromosome 2 / #2 | 2R:28,420,001-28,540,000 | 848 (493 | 354) |
H12 / Burkina Faso An. coluzzii / Chromosome 2 / #3 | 2R:28,420,001-28,560,000 | 557 (369 | 188) |
The information below provides some diagnostics from the Peak modelling algorithm.
Selection signal in context. @@TODO
Peak targetting. @@TODO
Peak fitting diagnostics. @@TODO
Left flank, peak model:
[[Model]]
Model(exponential)
[[Fit Statistics]]
# function evals = 35
# data points = 151
# variables = 3
chi-square = 0.031
reduced chi-square = 0.000
Akaike info crit = -1274.992
Bayesian info crit = -1265.940
[[Variables]]
amplitude: 0.38063638 +/- 0.010058 (2.64%) (init= 0.5)
decay: 0.24559296 +/- 0.010400 (4.23%) (init= 0.5)
c: 0.02019622 +/- 0.001292 (6.40%) (init= 0.03)
cap: 1 (fixed)
[[Correlations]] (unreported correlations are < 0.100)
C(amplitude, decay) = -0.607
C(decay, c) = -0.298
Right flank, peak model:
[[Model]]
Model(exponential)
[[Fit Statistics]]
# function evals = 39
# data points = 149
# variables = 3
chi-square = 0.028
reduced chi-square = 0.000
Akaike info crit = -1273.906
Bayesian info crit = -1264.894
[[Variables]]
amplitude: 0.34898913 +/- 0.010487 (3.01%) (init= 0.5)
decay: 0.35025840 +/- 0.015079 (4.31%) (init= 0.5)
c: 0.02480874 +/- 0.001290 (5.20%) (init= 0.03)
cap: 1 (fixed)
[[Correlations]] (unreported correlations are < 0.100)
C(amplitude, decay) = -0.706
C(decay, c) = -0.365
Left flank, null model:
[[Model]]
Model(constant)
[[Fit Statistics]]
# function evals = 6
# data points = 150
# variables = 1
chi-square = 0.358
reduced chi-square = 0.002
Akaike info crit = -903.690
Bayesian info crit = -900.679
[[Variables]]
c: 0.03437788 +/- 0.004002 (11.64%) (init= 0.03)
Right flank, null model:
[[Model]]
Model(constant)
[[Fit Statistics]]
# function evals = 6
# data points = 148
# variables = 1
chi-square = 0.353
reduced chi-square = 0.002
Akaike info crit = -891.818
Bayesian info crit = -888.821
[[Variables]]
c: 0.04210692 +/- 0.004026 (9.56%) (init= 0.03)