EMBO/sd-character-level-ner huggingface dataset:
Figure 5RNA-seq revealed 1067 DEGs between young HO-1-/- and HO-1+/+ LT-HSCs. RNA-seq revealed 595 DEGs between old HO-1-/- and HO-1+/+ LT-HSCs. Selected GOBP terms enriched in GSEA analysis in young HO-1-/- LT-HSCs. Selected GOBP terms enriched in GSEA analysis in old HO-1-/- LT-HSCs. Significant part of DEGs (31% total) in young HO-1-/- LT-HSCs overlaps with DEGs identified in LT-HSCs during physiological aging. PC1 in PCA analysis on selected 1148 genes separated young HO-1+/+ LT-HSCs from other groups, clustering the young HO-1-/- LT-HSCs together with old LT-HSCs. Hierarchical clustering on the 20% of most correlated 297 genes with PC1 showed 2 major clusters that includes young HO-1+/+ LT-HSCs in one and old LT-HSCs with young HO-1-/- LT-HSCs in second. Selected GOBP terms enriched in GSEA analysis based on 297 identified genes.
Figure 7Top: Confocal z-projections of dpr12∆50-81 heterozygous (A, n=15; C, n=10; E, n=10) and homozygous brains (B, n=8; D, n=18; F, n=24) expressing mCD8-GFP (CD8) driven by: (A, B) R10G03-Gal4 (PAM-DAN-γ4-Gal4); (C, D) R48H11-Gal4 (PAM-DAN-γ5-Gal4), The greyscale channels are sub-z-projections comprised of slices restricted to the γ-lobe region. White dashed line demarcates the γ-lobe. Bottom: Cartoons schematizing MB lobe structure and innervation by specific PAM-DANs or MBON. Confocal z-projections of dpr12∆50-81 heterozygous (A, n=15; C, n=10; E, n=10) and homozygous brains (B, n=8; D, n=18; F, n=24) expressing mCD8-GFP (CD8) driven by (E, F) R18H09-Gal4 (MBONγ4>γ1γ2 -Gal4). The greyscale channels are sub-z-projections comprised of slices restricted to the γ-lobe region. White dashed line demarcates the γ-lobe. Bottom: Cartoons schematizing MB lobe structure and innervation by specific PAM-DANs or MBON. Single confocal slices of WT (G, n=5/5) and dpr12∆50-81 homozygous brains (H, n=5/5) stained with anti-Brp. Dashed lines demarcate neuropil boundaries, schematic shown below. Cre, crepine, a neuropil that surrounds the medial MB lobes; AL, antenna lobe.
Figure 4A. Fatty acid analysis from untargeted MS metabolomics including free alpha-linolenic acid (ALA) and docosahexaenoic acid (DHA). B. Carnitine-bound fatty acids commonly found in serum that were detected in mass spectrometry. C. Lysophosphatidylcholine and lysophosphatidylethanolamine with C16:0 (palmitic acid) and C18:1 (oleic acid). D. Combined analysis of all found triglycerides correlating to medium-chain fatty acid lengths and long or very long chain fatty acid lengths.
Figure 3Configuration of the two-compartment model. The microenvironment of microbial swarmbot is separated from that of the chamber by its permeable hydrogel shell (orange), resulting in two separate compartments, defined as swarmbot and environment (chamber). Four major components involved in the population dynamics are explored in the model, including cell density (n), nutrient (s), carbenicillin (a) and BlaM concentration (b). Each component has a different transport rate, depending on the molecule size.Simulated microbial swarmbot dynamics under a static condition. a = 0.4, a microbial swarmbot with a low n0 cannot initiate growth (red lines) but one with a high n0 leads to growth both in the swarmbot and, after a delay, in the chamber (green lines). While (C) swarmbots with both low and high initial densities lead to eventual full growth in the chamber in the absence of the antibiotic (a = 0). Solid lines: cell density inside swarmbots; dashed lines: cell density in the chambers.(D) A microfluidic device for controlling and monitoring population dynamics of microbial swarmbot. The device consists of an input channel leading to a culturing chamber and a trapping slot. The chamber is used for culturing engineered bacterial populations both inside and outside the swarmbot, while the trapping slot is for trapping the swarmbot due to its smaller size but allowing the medium to flow through. Input flow patterns (static or pulsing) are controlled by a programmable syringe pump.(E) Experimental validation of model prediction. At 100µg/ml carbenicillin, a microbial swarmbot with a high initial density led to robust growth in the capsule and, after a delay, in the chamber (green lines). While a swarmbot with a low initial cell density did not exhibit growth either in the swarmbot or in the chamber at the same carbenicillin concentration (red lines). (F) In contrast, both the swarmbot and culturing chamber reached high cell densities in absence of carbenicillin, regardless of the initial cell densities. (G) Time-lapse images of microbial swarmbot dynamics. The population initiated growth inside the swarmbot and subsequently in the chamber at a high initial density. But no growth occurred for a swarmbot with a low initial density at the same antibiotic concentration (rightmost). The concentration of carbenicillin was 100 µg/ml and the scale bar represents 200 µm.(H) The control strain lacking BlaM circuit. The population could not survive at the same concentration of carbenicillin at the high initial density.
Figure 3a. Sorted CD4+CD25hiCD127- Tregs were stimulated in the presence of vehicle (white dots) or either SF1670 (first row), SC-79 (second row), 740 Y-P (third row) or AS1842856 (fourth row, black dots) for 3 days and IFNG and TBX21 gene expression was examined every 24 hours. Statistical analysis of n=4 independent experiments performed. b. Representative example of the frequency of IFN+Foxp3+ cells in Tregs stimulated in the presence of vehicle, SF1670, SC-79, 740 Y-P or AS1842856 for 3 days after 4 hours of stimulation with PMA and ionomycin. Numbers in brackets represent total cell numbers c. Statistical analysis of the frequency of IFN+Foxp3+ cells in Tregs stimulated in the presence of vehicle, SF1670, SC-79, 740 Y-P or AS1842856 of n=5 experiments performed. d. Proliferation of CFSE-labeled responder T cells alone (left histogram) or co-cultured with Tregs that have been pre-incubated with either vehicle, SF1670, SC-79 or 740 Y-P at day 4 after co-culture at a 1:2 Treg:Tresp ratio. e. Statistical analysis of the percentage of suppression at different Treg:Tresp ratios of n=6 experiments performed.
Figure 4A Schematic representation showing where the compounds act along the AC/PKA signalling pathway. B Lysates from HEK293T cells, transfected with 2R or 66R and treated for 24 h, immunoblotted using antibody for poly-GA expression (HA-tagged). Quantification from three biological replicates. Data are mean ± SEM. Two-tailed, unpaired t-test, ****P<0.0001, **P = 0.0027 H89 vs DMSO, ***P = 0.0004 FSK + H89 vs DMSO, *P = 0.0234 CLZ vs DMSO, *P = 0.0209 FSK + H89 vs FSK. Concentration used: 2.5 μM H89, 10 μM FSK, CLZ. Lysates from NSC34 cells, transfected with 2R or 66R and treated for 24 h, immunoblotted for poly-GP expression. (D) Quantification of PBS-insoluble poly-GP from three biological replicates. Data are mean ± SEM. 1-way ANOVA followed by uncorrected Fisher's LSD, *P = 0.0196 2R vs 66R DMSO, **P =0.0039 FSK vs DMSO, *P = 0.0418 H89 vs DMSO, **P = 0.0017 H89 vs FSK, *P = 0.0364 FSK + H89 vs FSK. Concentration used: 2.5 μM H89, 10 μM FSK (E) Quantification of poly-GP from three biological replicates. Data are mean ± SEM. 1-way ANOVA followed by uncorrected Fisher's LSD, *P = 0.0112 2R vs 66R DMSO, ****P<0.0001 FSK vs DMSO, ns (not significant) H89 vs DMSO, *P = 0.0281 FSK + H89 vs DMSO, ****P<0.0001 H89 vs FSK, ***P = 0.0004, FSK + H89 vs FSK, **P = 0.0030 H89 vs FSK +H89. Concentration used: 2.5 μM H89, 10 μM FSK F Luminescence derived from HEK293T transfected with GA-NLuc reporters WT or mutated and treated with H89 for 24 h. Mutational analysis of near-cognate CUG start codon from three biological replicates and four technical replicates, n = 12. Data are mean ± SEM. 2-way ANOVA followed by Dunnett's multiple comparisons test, *P = 0.0130 WT H89 vs WT DMSO, ****P<0.0001, CUG mutated H89 vs CUG mutated DMSO. Concentration used: 2.5 μM H89,
Cancer cells directly activate MDSCs in a TRF2-dependent manner via the TLR2/MyD88/IL-6/STAT3 pathway that inhibit NK cell degranulation and cytotoxicity.Figure 3Analysis of STAT3 phosphorylation by MSC2 cells after 18 hours' co-culture with BJcl2 overexpressing or compromised for TRF2 (A). Analysis of STAT3 phosphorylation by MSC2 cells after 18 hours' co-culture with BJcl2 overexpressing or compromised for TRF2 Histograms of the pSTAT3 mean fluorescence intensity (MFI) (B). +/- inhibitors: anti-MyD88 peptide inhibitor, anti-mouse TLR2 antibody, or lipopolysaccharide (LPS)-RS, an antagonist of TLR4 and anti-mouse IL-6 or JAK1/2 inhibitor (Ruxolitinib) +/- inhibitors: anti-MyD88 peptide inhibitor, anti-mouse TLR2 antibody, or lipopolysaccharide (LPS)-RS, an antagonist of TLR4 in (D) +/- inhibitors: anti-mouse IL-6 or JAK1/2 inhibitor (Ruxolitinib) in (E). Schematic representation of the MDSC suppression assay (F). TRF2-overexpressing or -knockdown BJcl2 cells were co-cultured with MSC2 cells for 18h +/- mouse IL6 blocking antibody or Ruxolitinib. After FACS sorting, MSC2 cells were co-cultured for 18h with primed and purified NK cells before a 4h challenge with YAK-1. NK cell degranulation and IFN-γ production were analyzed by FACS (G). Data information: Panels show all points with mean ± SEM corresponding to the percentages of each cell population among NK cells (G) (n = 3 independent experiments; *p < 0.05, **p < 0.01, and ***p < 0.001; Student's t-test). TRF2-overexpressing or -knockdown BJcl2 cells were co-cultured with MSC2 cells for 18h +/- mouse IL6 blocking antibody or Ruxolitinib. Alternatively, a 4h challenge with NIH3T3 cells was performed and the ability of NK cells to kill NIH3T3 cells were determined by colorimetric analysis (AlamarBlue) (H). Data information: Panels show all points with mean ± SEM corresponding to the percentage of specific target lysis (H) (n = 3 independent experiments; *p < 0.05, **p < 0.01, and ***p < 0.001; Student's t-test).
Effects of proteasome inhibition and ROS overproduction on APOPT1 stability A Western blot analysis of SDS-PAGE of MGM132-treated 143B APOPT1HA cells.Figure 6 A Western blot analysis of SDS-PAGE of MGM132-treated 143B APOPT1HA cells. The graph represents the densitometric quantification of the signals for the precursor and mature protein  B The upper part of the panel (Input) shows Western blot analyses of APOPT1HA in the MGM132-treated cells. Note the appearance of higher molecular weight bands upon longer exposure in the samples treated with the proteasome inhibitor. The bottom part of the panel (Purified fractions) shows the analysis of fractions from the same cells immunoprecipitated with anti-HA. Note that the higher molecular weight species are cross-reacting with both anti-HA and anti-ubiquitin  C Western blot analysis of SDS-PAGE of total lysates from 143B cells overexpressing tagged APOPT1 (as indicated) and exposed to 100 µM H2O2, as illustrated by the scheme (H2O2 treatment), for the indicated times. The graphs represent the densitometric quantification of the tagged APOPT1 signal at each time point. The graph inset shows that the increase of APOPT1 occurs in the first minutes after the exposure to H2O2  D Western blot analysis of SDS-PAGE of total lysates from 143B cells overexpressing tagged APOPT1 (as indicated) and exposed to 5 µM MitoParaquat (MitoPQ), as illustrated by the scheme (MitoPQ treatment), for the indicated times. The graphs represent the densitometric quantification of the tagged APOPT1 signal at each time point. The graph inset shows that the increase of APOPT1 occurs in the first minutes after the exposure to MitoPQ
Accumulation of autophagy adaptors NBR1 and SQSTM1 and the phagophore membrane component LC3 is induced in fibroblasts from two controls and a subject with Vici syndrome (subject 4.figf3Accumulation of autophagy adaptors NBR1 and SQSTM1 and the phagophore membrane component LC3 is induced in fibroblasts from two controls and a subject with Vici syndrome (subject 4.1) after 12 h of treatment with rapamycin or dual treatment with rapamycin and bafilomycin.
Figure 8(A) MEFs were then infected with retroviral vectors encoding the indicated B-Raf proteins. Giemsa staining revealed transformed and hyper-proliferated cells by dark purple colour. (B) Quantification of Focus formation potential of the B-Raf mutants used in (A). n = 3, mean + SEM, t-test, *** p < 0.001. (C) Following recruitment by Ras-GTP, B-Raf forms homo- or heterodimers with the assistance of 14-3-3 proteins. This loosely dimerized and membrane-bound state is imitated by the B-RafE586K and B-RafCAAX mutants, which still require an intact DIF and AL phosphorylation. Within the dimer, one protomer becomes phosphorylated at its TVKS-motif in cis (depicted by black bent arrow) leading to re-orientation of the DIF and induction of a conformational change in the other protomer (blue bent arrow) followed by phosphorylation of the AS in the other protomer in cis. B-RafV600E is locked in the highly active state and thus signals independent of an intact TVKS-motif or DIF. See text for further details.
a,b) Effects of G120A and F52A mutations, which perturbed C18-ceramide binding of LC3B, on targeting of mitochondria by LC3B-containing autophagosomes (visualized by confocal microscopy using anti-Flag and anti-Tom20) (a) or mitochondrial function (measurement of OCR) (b).figf6(a,b) Effects of G120A and F52A mutations, which perturbed C18-ceramide binding of LC3B, on targeting of mitochondria by LC3B-containing autophagosomes (visualized by confocal microscopy using anti-Flag and anti-Tom20) (a) or mitochondrial function (measurement of OCR) (b). V, vector. (c) Effects of siRNA-mediated knockdown of p62, NIX and Drp1 on OCR in the absence or presence of CerS1 and C18-ceramide induction were measured using the SeaHorse. Data shown are an average of at least three experiments ± s.d. (*P 0.05). (d) Effects of siRNA-mediated knockdown of Drp1 (lower left and right) on the localization of ceramide within mitochondrial membranes were visualized using the colocalization of anti-ceramide (red) and anti-Tom20 (green) under confocal microscopy, in the absence (upper and lower left) or presence (upper and lower right) of CerS1 and C18-ceramide induction (−/+ tet) compared to Scr siRNA-transfected controls (upper left and right). Scale bars, 10 μm in a and d. Enlarged images are shown in boxed areas in d.
Figure 4(A) Experimental design of tumor initiation and application of WBRT. (B) Representative MRI pictures of untreated and irradiated mice at d0 and d14 after treatment start. (C) Quantification of Iba1+ macrophages in IHC sections of 99LN-BrM at trial endpoint (Ctrl n=27, WBRT n=18). (D) Representative IF images of 99LN-BrM stained for the macrophage marker Iba1 (red) and the MG marker Tmem119 (white). Higher magnification images present areas in the peri-tumor region and in the tumor core. DAPI was used as nuclear counterstain. Scale bar; 100 µm. (E) Quantification of different myeloid populations in 99LN-BrM tumors with and without irradiation by flow cytometry at d14 (Ctrl n=5, WBRT n=4). (F) Representative IF images of 99LN-BrM with and without WBRT stained for DCIR2 (red) and EpCAM (green). DAPI (blue) was used as nuclear counterstain. Scale bar; 100 µm. (G) Quantification of CD11c+MHCII+ immune cells by flow cytometry in 99LN-BrM with and without WBRT. (H, I) Relative proportion of cDC1 and cDC2 in 99LN-BrM with and without WBRT. (J) cDC1:cDC2 ratio in 99LN-BrM based on flow cytometry data (H+I). (K) Representative IHC images of CD3+, FoxP3+ and CD8+ T cells in 99LN-BrM sections with and without WBRT. Scale bar; 100 µm. (L) Quantification of CD3+ cells in IHC sections of 99LN-BrM. (Ctrl n=26, WBRT n=18). (M) FoxP3:CD3 ratio in IHC sections of 99LN-BrM. (Ctrl n=25, WBRT n=18). (N) CD8:CD3 ratio in IHC sections of 99LN-BrM. (Ctrl n=25, WBRT n=18). (O-Q) Quantification of CD3+, CD4+ and CD8+ T cells in 99LN-BrM tumors with and without irradiation by flow cytometry. (R) CD4:CD8 ratio in 99LN-BrM based on data (P+Q).
Phenotypic characterisation of Δatg5.Figure 4(A) Growth curve of Leishmania promastigotes in HOMEM medium at 26°C. *, Δatg5 differed significantly from WT (p<0.05). (B) The occurrence of GFP-ATG8 puncta in promastigotes after incubation in nutrient-deprived (PBS, ND) and nutrient-rich (HOMEM medium, NR) conditions for 2 h at 26°C. Scale bar, 10 µm (C) Occurrence of GFP-ATG8 puncta in promastigotes when incubated in the conditions detailed in (B). Means ± SD from four independent experiments. * and **, occurrence of GFP-ATG8 puncta in Δatg5 were significantly different from in WT in nutrient-deprived and nutrient-rich conditions (p<0.05). (D) Western blot analysis of extracts of promastigotes expressing GFP-ATG8 at logarithmic growth under standard conditions and probed with α-GFP antibody. The faster migrating, lipidated band is labelled GFP-ATG8-II while the un-lipidated band migrating more slowly is labelled GFP-ATG8-I.
Lack of myosin VI alters the mobility and localisation of OFD1 at the centrioles.Figure 4(A) Super-resolution analysis of OFD1 localisation at the centrioles. hTERT-RPE1 cells were transfected with siRNA against myosin VI, immunostained with anti-OFD1 and anti-acetylated tubulin antibodies and visualised using structured illumination microscopy (SIM). Representative images are shown, scale bar, 500 nm. A scheme of the estimated localisation of OFD1 in the two conditions is depicted on the right side. (B) Two-colour SIM images in control and myosin VI-depleted cells illustrating the localisation and distribution of the distal appendage markers Cep164 and ODF2 at the centrioles, stained with anti-acetylated tubulin antibody. Representative images are shown; scale bar, 500 nm. (C) dSTORM super-resolution analysis of the distal appendage markers FBF1 and ODF2 in control and myosin VI-depleted cells. To emphasise the symmetry of the structures, the signals from all nine appendages were averaged (bottom). Representative images are shown; scale bar, 200 nm. FRAP analysis of centriole-associated GFP-OFD1. hTERT-RPE cells stably expressing GFP-OFD1 and centrin1-dTomato were transfected with siRNA against myosin VI. After four days, cells were treated with nocodazole (1 hour, 6 μg/ml) and subjected to live-cell imaging. (D) Left: a scheme of the localisation of GFP-OFD1 and the centriole marker centrin1-dTomato. Right: a scheme of photobleaching and recovery of GFP-OFD1 at the centrioles. FRAP analysis of centriole-associated GFP-OFD1. hTERT-RPE cells stably expressing GFP-OFD1 and centrin1-dTomato were transfected with siRNA against myosin VI. After four days, cells were treated with nocodazole (1 hour, 6 μg/ml) and subjected to live-cell imaging. (E) A representative graph of one out of three experiments. For each time point, the fraction of recovery of GFP-OFD1 is shown. Results are expressed as means with 95% confidence interval. n=12 cells (Mock), n=13 cells (MyoVI KD). (F) Quantification of the half-time of fluorescence recovery (t1/2,) and of the mobile fraction of GFP-OFD1. Results are expressed as mean ±SD. n=41 cells (Mock), n=31 cells (MyoVI KD), from three independent experiments. ns, not significant; *** P<0,0005 by Unpaired T-test.
Figure 5A: Representative AP traces measured in JLNSR594Q-CMs paced at 1 Hz, under baseline conditions (black) and after application of either AST alone (red), AST+LUF7346 3 µM (green) and AST+LUF7346 5 µM (blue). B: STV values calculated on spontaneously beating hiPSC-CM clusters over 30 consecutive beats as measured with MEA. N = 7-24. * = p<0.05 vs respective isogenic pair. C: STV changes (left) and representative examples of QT interval oscillations plotted over 30 consecutive beats (right) in the presence of AST and AST+LUF7346 20 µM. * = p<0.05 vs Baseline. ° = p<0.05 vs AST. N = 7-24.
Figure 5(A) Representative immunoblots showing NFATc1 degradation in HEK293T cells stably overexpressing NFATc1 in the presence of cycloheximide (CHX), with or without SIRT7 overexpression. (B) Quantification of NFATc1 levels in (A) from 3 independent experiments. (C) Immunoblots of Nfatc1 protein levels in Sirt7+/+ and Sirt7-/- keratinocytes. (D) Quantification of Nfatc1 levels in (C) from 3 independent experiments. (E-F) Immunoblots (E) and related quantitative data (F, n=3) showing NFATc1 protein levels in HEK293T cells overexpressing wildtype HA-SIRT7 or a SIRT7 enzyme-dead mutant (H187Y). (G-H) Immunoblots (G) and quantitative data (H, n=3) showing NFATc1 levels in HEK293T cells transfected with NFATc1-CA, with or without ectopic SIRT7, treated with MG132. (I-J) Immunoblots (I) and related quantitative data (J, n=3) showing NFATc1-WT and NFATc1-K612R protein turnover in HEK293T cells. Immunoblots showing the interaction (by co-immunoprecipitation) between ectopic NFATc1/PA28γ in HEK293T cells (K) (L) after manipulating SIRT7 levels. Immunoblots showing the interaction (by co-immunoprecipitation) between endogenous NFATc1/PA28γ in HaCaT cells (L) after manipulating SIRT7 levels. (M-N) Immunoblot (M) and quantitative data (N, n=3) analysis of NFATc1 levels upon overexpression of HA-SIRT7 with or without PA28γ oligo siRNAs treatment. (O) Co-immunoprecipitation between NFATc1 and PA28γ in HEK293T cells cotransfected with WT HA-SIRT7 or an enzyme-dead mutant (H187Y).
Figure 8A Immunoblotting of HEK293T cells co-transfected for 48 h with GFP-TLR7, plus HA-tagged SHP2 or SHP2 mutant vectors, followed by immunoprecipitation with anti-GFP. B, Immunoblotting of TLR7 in the endosome isolated from HEK293T cells that were transfected with GFP-TLR7 or GFP-TLR7 mutant vectors. C, Immunoblotting of p-p65 in lysates of HEK293T cells transfected with GFP-TLR7 and TLR7 mutant vectors and then treated with IMQ (10 µg/ml) for 30 min. D, HEK293T cells transfected with NF-κB-luciferase reporter and indicated TLR7 mutant vectors and then treated with IMQ (10 µg/ml) for 24 h were harvested for luciferase assay. Data information: Data are represented mean ± SEM. P values are determined by two-tailed unpaired Student's t-test **P<0.01, ns, not significant. E Immunoblotting of TLR7 in the endosome isolated from HEK293T cells that were transfected with GFP-TLR7 or GFP-TLR7 mutant vectors. F Immunoblotting of p-p65 in lysates of HEK293T cells transfected with GFP-TLR7 and TLR7 mutant vectors and then treated with IMQ (10 µg/ml) for 30 min. G HEK293T cells transfected with NF-κB-luciferase reporter and indicated TLR7 mutant vectors and then treated with IMQ (10 µg/ml) for 24 h were harvested for luciferase assay. Data information: Data are represented mean ± SEM. P values are determined by two-tailed unpaired Student's t-test **P<0.01, ns, not significant. H Immunofluorescence staining of HEK293T cells transfected with GFP-TLR7 and TLR7 mutant vectors, prior to treatment of R848 (10 µg/ml) for 30 min and labelled with specific antibodies. Scale bar: 50 μm. I Representative p-TLR7 Y1024 staining of skin sections from human. Scale bar: 100 μm.
SNP overview and unbiased selection according to the standard operating procedure (SOP) developed for phenotype-based genetic association study (PGAS) approaches.Figure 2Figure 2: SNP overview and unbiased selection according to the standard operating procedure (SOP) developed for phenotype-based genetic association study (PGAS) approaches. (A) Genes from the 'broader fragile X family' and their chromosomal position. (B) SNP numbers available through direct genotyping in the here used semi-custom genotyping array (Affymetrix, Santa Clara, CA, USA). (C) SNPs fulfilling some of the first round of selection criteria ('functional'=SNPs, i.e. located in promoter region, 3'UTR or coding sequence; MAF=MAF≥0.2; LD=SNPs that 'survive' after linkage disequilibrium pruning: r²<0.8). Underlined are the 13 SNPs selected for the PGAS approach using the PAUSS (selection requirements: fulfilled 2 of the above criteria or were functional). Not more than 3 SNPs per gene are selected to avoid overrepresentation of one gene. (D) SNPs with a tendency in PGAS (see Fig 3A) at single SNP basis went into the final accumulation model.