|Home | About | Journals | Submit | Contact Us | Français|
Conceived and designed the experiments: YM FC. Performed the experiments: YM. Analyzed the data: YM. Wrote the paper: YM FC.
In genetic epidemiology, genome-wide association studies (GWAS) are used to rapidly scan a large set of genetic variants and thus to identify associations with a particular trait or disease. The GWAS philosophy is different to that of conventional candidate-gene-based approaches, which directly test the effects of genetic variants of potentially contributory genes in an association study. One controversial question is whether GWAS provide relevant scientific outcomes by comparison with candidate-gene studies. We thus performed a bibliometric study using two citation metrics to assess whether the GWAS have contributed a capital gain in knowledge discovery by comparison with candidate-gene approaches. We selected GWAS published between 2005 and 2009 and matched them with candidate-gene studies on the same topic and published in the same period of time. We observed that the GWAS papers have received, on average, 30±55 citations more than the candidate gene papers, 1 year after their publication date, and 39±58 citations more 2 years after their publication date. The GWAS papers were, on average, 2.8±2.4 and 2.9±2.4 times more cited than expected, 1 and 2 years after their publication date; whereas the candidate gene papers were 1.5±1.2 and 1.5±1.4 times more cited than expected. While the evaluation of the contribution to scientific research through citation metrics may be challenged, it cannot be denied that GWAS are great hypothesis generators, and are a powerful complement to candidate gene studies.
The first genome-wide association study (GWAS) was published in 2005 and identified a genetic variant associated with a higher risk of age-related macular degeneration . The completion of the Human Genome Project  and the HapMap Project  yielded tools to identify common genetic variations, mainly single-nucleotide polymorphisms (SNPs), associated with many traits or diseases. In 2008, the National Human Genome Research Institute (NHGRI) published an online catalog of Genome-Wide Association Studies , listing genetic marker-trait associations published in GWAS. As of 01 November 2012, this database included 1416 publications and 7688 disease- or trait-associated genetic variants.
These “wide association” approaches, initially devoted to genetic associations, have since been extended to transcriptomic , , metabolomic  and environmental factors (chemical toxicants, pollutants or nutrients) . While examination of many common genetics variants  or environmental exposures  is relevant because of their known role in the development of diseases, one may wonder if simultaneous analysis of such an overwhelming amount of data using data-driven methods leads to a significant value in the generation of knowledge, compared with the classical hypothesis-driven candidate-gene studies.
To answer this question, we conducted a bibliometric survey using two citation metrics, the citation count and the crown index, as a surrogate marker for the contribution to the scientific research of GWAS and candidate-gene studies focusing on the same trait or disease and published during the same time period in a predefined set of journals.
We searched for GWAS published between 2005 and 2009 in Nature, Science and journals ranked among the 10 leading journals according to the 2010 Journal Impact Factors (JIFs), in the “Genetics & Heredity” or “Medicine, General & Internal” categories (see Table S1). Meta-analyses, reviews and other publication types were excluded. GWAS had to be performed in humans.
The first search was conducted in the MEDLINE database with the search terms “GWAS”, “Genome-Wide Association Study”, “Genome-Wide Association Studies”, “Genome Wide Association Study”, “Genome Wide Association Studies”, “Genomewide Association Study” and “Genomewide Association Studies”, “Whole Genome Association Study”, “Whole Genome Association Studies”, “WGAS”, “WGA study” and “WGA studies”. To ensure an exhaustive selection process, we retrieved from the NHGRI Catalog of Genome-Wide Association Studies those studies not identified by the MEDLINE search.
The search for candidate-gene studies was conducted using BioPython , a set of freely available tools for biological computation, written in Python. For each paper in the GWAS group, we searched for papers focusing on the same disease or trait and published the same year in the same set of journals. We selected the paper that allowed us to reduce both the difference between the 2010 JIFs of the journals in which the GWAS and candidate-gene study were published, and to reduce the interval (in months) between the 2 publication dates. The studied trait or disease was, if available, retrieved from the NHGRI Catalog of Genome-Wide Association Studies, or was otherwise added manually.
The “quantity” of knowledge provided by the GWAS and candidate-gene studies was assessed with 2 metrics, namely the citation count and the crown index, computed 1 and 2 years after the publication dates. The citation count was recorded in the Web of Science database. The crown index depends on the expected citation rate, which indicates how frequently cited a paper is expected to be according to its year of publication, its publication type, and the journal in which it was published. The crown index is the ratio of the citation count to the expected citation rate.
Our main objective was to explore the differences of citation metrics between the GWAS and candidate-gene studies. To help interpretation of findings, we also compared candidate-gene studies derived from GWAS with those not derived from GWAS; we investigated the association between the reported P-value and citation metrics and we evaluated the impact of residual confounding, i.e. whether residual differences in impact factor or publication dates despite matching may have influenced the results.
The statistical analysis was conducted using R 2.13.0. The alpha level for all tests was 0.05.
Wilcoxon signed-rank tests for paired samples were used to test differences of quantitative variables between GWAS and candidate gene studies. Linear mixed-effects models were used to assess the effect of the approach used, “Genome-Wide” or “candidate gene”, on the crown index. A random effect was included to account for the correlation between the citation metrics in each pair of papers. The models were fitted with and without the effect of the approach and likelihood ratio tests were performed to compare the two fits.
GWAS and candidate-gene studies statistical tests for association do not usually use the same threshold to assess the statistical significance of an association. For example, the significance threshold in GWAS is often defined between 10−8 and 10−6 while it is usually close to 10−2 in candidate-gene studies. To investigate the association of P-values with citation metrics, we therefore calculated in each study the ratio between the reported P-value and the corresponding probability threshold used to conclude statistical significance (hereafter P-ratio). We explored the association between the log-transformed P-ratio, the type of study and the crown indexes using mixed models as described above.
We finally used linear regression models to assess the impact of the difference between the 2010 JIFs and of the interval between the publication dates on the difference in the crown indexes within each pair of papers.
Note: Summary statistics are reported as mean (± sd). The estimates of the linear mixed-effects models are reported as estimate (± se).
The initial search in the MEDLINE database identified 511 articles, and 68 additional articles were identified in the NHGRI Catalog. A total of 211 GWAS were selected, of which a candidate-gene study was identified in 97 cases (62 different papers) (see Tables S2, S3 and S4 for summary statistics of the 97 GWAS/candidate-gene study pairs and Tables S5, S6, S7 and S8 for the PMID PUBMED identifiers of the articles used in the survey). The full article selection process is summarized in Figure 1. Of the 62 candidate-gene selected, 23 of them (37%) were derived from former GWAS (not necessarily from the GWAS with which the candidate-gene study was paired).
We observed that the GWAS papers were more frequently cited than the candidate-gene studies, whatever the metric used (see Table 1 and Figure 2). This was confirmed by comparing the citation counts and the crown indexes, 1 and 2 years after the publication dates. The two-sided Wilcoxon signed-rank tests were all significant (P value =1.8×10−8 and 7.8×10−11 for the citation counts; P value =2.6×10−7 and 4.7×10−12 for the crown indexes). The distribution of the crown indexes were not significantly different between the candidate-gene studies derived from GWAS, compared with those not derived from GWAS, both 1 year after the publication date (Wilcoxon test P value =0.67) and 2 years after the publication date (Wilcoxon test P value =0.46).
Assessment of the effect of the “genome-wide” or “candidate-gene” approach on the crown indexes based on linear mixed-effects models confirmed that this metric increased in GWAS papers compared with candidate-gene papers. The estimates of the “genome-wide” effect on the crown indexes were 1.31 (±0.23, likelihood ratio test P value =3.3×10−6) and 1.45 (±0.19, likelihood ratio test P value =4.3×10−7), 1 and 2 years after the publication dates.
The average log transformed P-ratio was lower in the GWAS group than in the candidate-gene group : −25.8 (±54.7) for the GWAS group, −12.6 (±24.9) for the candidate-gene group (Wilcoxon test P value =6.9×10−9). The estimates of the mixed models assessing the effect of the log transformed P-ratio on the crown index 1 and 2 years after the publication date were 3.0×10−13 (±8.6×10−8, P value =0.97) and −4.4×10−13 (±1.0×10−7, P value =0.98).
There was no evidence for residual confounding: the average difference between the 2010 JIFs of the journals in which the paired papers were published was 1.5 (±15.9, Wilcoxon test P value =0.24) and the average interval between the publication dates was 2.7 (±2.9) months (Wilcoxon test P value =0.71). The effects of these differences on the difference in the crown indexes within each pair of papers in linear models were not significant, either one year after publication (P value =0.56 for the difference in 2010 JIFs, and P value =0.37 for the interval between the publication dates) or 2 years after publication (P value =0.90 for the difference in 2010 JIFs and P value =0.36 for the interval between the publication dates).
We showed that GWAS had a higher impact in the scientific literature than candidate gene studies. We assume that our study selection process, although not systematic, produced robust results by focusing on the genetic analyses published in the papers with the most important impact factors, and that our matching process did not alter the quality of the results provided.
The use of citations metrics to assess the ability of GWAS and candidate-gene studies to generate knowledge can be questioned. Indeed it was not possible to directly link the citation metrics to the relevance or the applicability of the scientific outcomes provided by these studies, or to their contribution to understanding the disease’s etiology. An in-depth analysis of the papers citing the GWAS and the candidate-gene studies would have been required to observe how the findings provided by those papers were used. However, we can consider that the citation metrics used are good markers for “research” generation, and therefore good, albeit indirect, markers for knowledge generation. Indeed, we observed that nearly 40% of the candidate-gene studies selected in our survey were inspired by hypotheses derived from GWAS.
The investigation of the potential association of the reported P values in the GWAS and the candidate-gene studies and of the crown index did not permit to show a significant relationship. Although the statistical threshold to assess significance is stricter in GWAS because of the multiplicity of tests, their findings were not associated with a higher scientific impact, and not associated with a higher credibility for follow-up.
We have explored the impact of GWAS and candidate-gene studies through 2 citation metrics and observed that the GWAS papers were more cited than the candidate-gene papers, whatever the metrics used. If it cannot formally be claimed that the knowledge generated by GWAS outweighs that provided by classical genetic association studies, we can say that this manuscript provides evidence, at least indirect evidence, that GWAS contribute to knowledge production by allowing candidate-gene studies to focus on credible candidates.
2010 Journal Impact Factors of the Journals used for the MEDLINE database and NHGRI Catalog research.
Years of publication of the 97 pairs of papers in the “GWAS” and “candidate-gene studies” groups.
Journals of publication of the “GWAS” papers.
Journal of publication of the “candidate-gene studies” papers.
PMID PUBMED identifiers of the 511 articles identified through MEDLINE research.
PMID PUBMED identifiers of the 68 articles identified through NHGRI catalog research.
PMID PUBMED identifiers of the 211 GWAS articles conserved for candidate-gene study research.
PMID PUBMED identifiers of the 97 GWAS/candidate-gene study pairs.
We thank Rosemary M. Delabre and David Young for critical review of the manuscript.
This work was supported by the Ministry of Higher Education and Research of France. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.