Abstract
Motivation Genetic risk scores (GRS) summarise genetic data into a single number and allow for discrimination between cases and controls. Many applications of GRSs would benefit from comparisons with multiple datasets to assess quality of the GRS across different groups. However, genetic data is often unavailable. If summary statistics of the genetic data could be used to simulate GRSs more comparisons could be made, potentially leading to improved research.
Results We present a methodology that utilises only summary statistics of genetic data to simulate GRSs with an example of a type 1 diabetes (T1D) GRS. An example on European populations of the mean T1D GRS for real and simulated data are 10.31 (10.12-10.48) and 10.38 (10.24-10.53) respectively. An example of a case-control set for T1D has a area under the receiver operating characteristic curve of 0.917 (0.903-0.93) for real data and 0.914 (0.898-0.929) for simulated data.
Availability The code is available at https://github.com/stevensquires/simulating_genetic_risk_scores.
Contact s.squires{at}exeter.ac.uk
Competing Interest Statement
The authors have declared no competing interest.
Funding Statement
This study was supported by a grant from the Randox Corporation. This study was supported by the National Institute for Health and Care Research Exeter Biomedical Research Centre. The views expressed are those of the author(s) and not necessarily those of the NIHR or the Department of Health and Social Care.
Author Declarations
I confirm all relevant ethical guidelines have been followed, and any necessary IRB and/or ethics committee approvals have been obtained.
Yes
The details of the IRB/oversight body that provided approval or exemption for the research described are given below:
The study used only openly available human data. The three datasets used were: 1000G available at https://www.internationalgenome.org/; UK Biobank available at https://www.ukbiobank.ac.uk/; T1DGC available at https://repository.niddk.nih.gov/studies/t1dgc/.
I confirm that all necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived, and that any patient/participant/sample identifiers included were not known to anyone (e.g., hospital staff, patients or participants themselves) outside the research group so cannot be used to identify individuals.
Yes
I understand that all clinical trials and any other prospective interventional studies must be registered with an ICMJE-approved registry, such as ClinicalTrials.gov. I confirm that any such study reported in the manuscript has been registered and the trial registration ID is provided (note: if posting a prospective study registered retrospectively, please provide a statement in the trial ID field explaining why the study was not registered in advance).
Yes
I have followed all appropriate research reporting guidelines, such as any relevant EQUATOR Network research reporting checklist(s) and other pertinent material, if applicable.
Yes
Data Availability
The study used only openly available human data. The three datasets used were: 1000G available at https://www.internationalgenome.org/; UK Biobank available at https://www.ukbiobank.ac.uk/; T1DGC available at https://repository.niddk.nih.gov/studies/t1dgc/.
