H12/CMS/2/3

This page describes a signal of selection found in the Cameroon An. gambiae populationusing the H12 (Garud et al. 20XX) statistic.The focus of this signal is on chromosome arm 2R between positions 40,300,000 and 40,740,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 58 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, 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), AGAP003624, AGAP029110, AGAP003626, AGAP013221, AGAP012946, AGAP003627, AGAP003629, AGAP003630, AGAP003631 (GPRGRP2 - putative gastrin/bombesin receptor 2).

The following 3 genes are within 50 kbp of the focal region: AGAP003580, AGAP003632 (U3 small nucleolar RNA-associated protein 14), AGAP003633.

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
IHS/BFS/2/3 IHS Burkina Faso An. gambiae 2R:40,320,000-40,560,000 557 99.5% nan
IHS/GNS/2/3 IHS Guinea An. gambiae 2R:40,200,000-40,420,000 450 99.0% nan
H12/BFS/2/4 H12 Burkina Faso An. gambiae 2R:40,640,000-41,060,000 384 95.9% nan
XPEHH/CMS.GWA/2/3 XPEHH Cameroon An. gambiae 2R:40,680,000-40,760,000 100 96.5% 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   = 54
    # data points      = 396
    # variables        = 4
    chi-square         = 0.084
    reduced chi-square = 0.000
    Akaike info crit   = -3341.844
    Bayesian info crit = -3325.919
[[Variables]]
    center:      0 (fixed)
    amplitude:   0.08320422 +/- 0.002716 (3.27%) (init= 0.5)
    decay:       2.99999974 +/- 0.236898 (7.90%) (init= 0.5)
    skew:       -0.93490345 +/- 0.053920 (5.77%) (init= 0)
    baseline:    0.02065198 +/- 0.001923 (9.31%) (init= 0.03)
    ceiling:     1 (fixed)
    floor:       0 (fixed)
[[Correlations]] (unreported correlations are <  0.100)
    C(decay, baseline)           = -0.778
    C(amplitude, baseline)       = -0.466
    C(skew, baseline)            = -0.329
    C(decay, skew)               =  0.307
    C(amplitude, skew)           =  0.306

Null model:

[[Model]]
    Model(constant)
[[Fit Statistics]]
    # function evals   = 11
    # data points      = 395
    # variables        = 1
    chi-square         = 0.296
    reduced chi-square = 0.001
    Akaike info crit   = -2840.165
    Bayesian info crit = -2836.186
[[Variables]]
    c:   0.05095336 +/- 0.001379 (2.71%) (init= 0.03)

Comments