Leveraging European infrastructures to access 1 million human genomes by 2022.

Human genomics is undergoing a step change from being a predominantly research-driven activity to one driven through health care as many countries in Europe now have nascent precision medicine programmes. To maximize the value of the genomic data generated, these data will need to be shared between institutions and across countries. In recognition of this challenge, 21 European countries recently signed a declaration to transnationally share data on at least 1 million human genomes by 2022. In this Roadmap, we identify the challenges of data sharing across borders and demonstrate that European research infrastructures are well-positioned to support the rapid implementation of widespread genomic data access.

Author Correction: Leveraging European infrastructures to access 1 million human genomes by 2022.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

GWAS of self-reported mosquito bite size, itch intensity and attractiveness to mosquitoes implicates immune-related predisposition loci.

Understanding the interaction between humans and mosquitoes is a critical area of study due to the phenomenal burdens on public health from mosquito-transmitted diseases. In this study, we conducted the first genome-wide association studies (GWAS) of self-reported mosquito bite reaction size (n = 84,724), itchiness caused by bites (n = 69,057), and perceived attractiveness to mosquitoes (n = 16,576). In total, 15 independent significant (P < 5×10-8) associations were identified. These loci were enriched for immunity-related genes that are involved in multiple cytokine signalling pathways. We also detected suggestive enrichment of these loci in enhancer regions that are active in stimulated T-cells, as well as within loci previously identified as controlling central memory T-cell levels. Egger regression analysis between the traits suggests that perception of itchiness and attractiveness to mosquitoes is driven, at least in part, by the genetic determinants of bite reaction size.Our findings illustrate the complex genetic and immunological landscapes underpinning human interactions with mosquitoes.

Functional confirmation that the R1488* variant in SCN9A results in complete loss-of-function of Nav1.7.

BACKGROUND: Individuals with an extremely rare inherited condition, termed Congenital Insensitivity to Pain (CIP), do not feel pain in response to noxious stimuli. Variants in SCN9A, encoding the transmembrane voltage-gated sodium channel Nav1.7, have previously been reported in subjects with CIP accompanied by anosmia, which are typically transmitted in a recessive pattern. Functional characterisations of some of these SCN9A mutations show that they result in complete loss-of-function of Nav1.7. METHODS: In a consanguineous family we performed whole exome sequencing of three members who have a diagnosis of CIP and one unaffected family member. The functional effects of the segregating variant in SCN9A were determined using patch clamp electrophysiology in human embryonic kidney (HEK) 293 cells transfected with the variant. RESULTS: We found that each CIP subject was homozygous for a putatively nonsense variant, R1488*, in SCN9A. This variant was reported elsewhere in a subject with CIP, though the functional effect was not determined. Using electrophysiology, we confirm that this variant results in a complete loss-of-function of Nav1.7. CONCLUSIONS: We confirm through electrophysiological analysis that this R1488* variant in SCN9A results in complete loss-of-function of Nav1.7, which is consistent with reports on other variants in this gene in subjects with CIP.

Genetic and Environmental Risk for Chronic Pain and the Contribution of Risk Variants for Major Depressive Disorder: A Family-Based Mixed-Model Analysis.

BACKGROUND: Chronic pain is highly prevalent and a significant source of disability, yet its genetic and environmental risk factors are poorly understood. Its relationship with major depressive disorder (MDD) is of particular importance. We sought to test the contribution of genetic factors and shared and unique environment to risk of chronic pain and its correlation with MDD in Generation Scotland: Scottish Family Health Study (GS:SFHS). We then sought to replicate any significant findings in the United Kingdom Biobank study. METHODS AND FINDINGS: Using family-based mixed-model analyses, we examined the contribution of genetics and shared family environment to chronic pain by spouse, sibling, and household relationships. These analyses were conducted in GS:SFHS (n = 23,960), a family- and population-based study of individuals recruited from the Scottish population through their general practitioners. We then examined and partitioned the correlation between chronic pain and MDD and estimated the contribution of genetic factors and shared environment in GS:SFHS. Finally, we used data from two independent genome-wide association studies to test whether chronic pain has a polygenic architecture and examine whether genomic risk of psychiatric disorder predicted chronic pain and whether genomic risk of chronic pain predicted MDD. These analyses were conducted in GS:SFHS and repeated in UK Biobank, a study of 500,000 from the UK population, of whom 112,151 had genotyping and phenotypic data. Chronic pain is a moderately heritable trait (heritability = 38.4%, 95% CI 33.6% to 43.9%) that is significantly concordant in spouses (variance explained 18.7%, 95% CI 9.5% to 25.1%). Chronic pain is positively correlated with depression (ρ = 0.13, 95% CI 0.11 to 0.15, p = 2.72x10-68) and shows a tendency to cluster within families for genetic reasons (genetic correlation = 0.51, 95%CI 0.40 to 0.62, p = 8.24x10-19). Polygenic risk profiles for pain, generated using independent GWAS data, were associated with chronic pain in both GS:SFHS (maximum ß = 6.18x10-2, 95% CI 2.84 x10-2 to 9.35 x10-2, p = 4.3x10-4) and UK Biobank (maximum ß = 5.68 x 10-2, 95% CI 4.70x10-2 to 6.65x10-2, p < 3x10-4). Genomic risk of MDD is also significantly associated with chronic pain in both GS:SFHS (maximum ß = 6.62x10-2, 95% CI 2.82 x10-2 to 9.76 x10-2, p = 4.3x10-4) and UK Biobank (maximum ß = 2.56x10-2, 95% CI 1.62x10-2 to 3.63x10-2, p < 3x10-4). Limitations of the current study include the possibility that spouse effects may be due to assortative mating and the relatively small polygenic risk score effect sizes. CONCLUSIONS: Genetic factors, as well as chronic pain in a partner or spouse, contribute substantially to the risk of chronic pain for an individual. Chronic pain is genetically correlated with MDD, has a polygenic architecture, and is associated with polygenic risk of MDD.

A common coding variant in CASP8 is associated with breast cancer risk.

The Breast Cancer Association Consortium (BCAC) has been established to conduct combined case-control analyses with augmented statistical power to try to confirm putative genetic associations with breast cancer. We genotyped nine SNPs for which there was some prior evidence of an association with breast cancer: CASP8 D302H (rs1045485), IGFBP3 -202 C --> A (rs2854744), SOD2 V16A (rs1799725), TGFB1 L10P (rs1982073), ATM S49C (rs1800054), ADH1B 3' UTR A --> G (rs1042026), CDKN1A S31R (rs1801270), ICAM5 V301I (rs1056538) and NUMA1 A794G (rs3750913). We included data from 9-15 studies, comprising 11,391-18,290 cases and 14,753-22,670 controls. We found evidence of an association with breast cancer for CASP8 D302H (with odds ratios (OR) of 0.89 (95% confidence interval (c.i.): 0.85-0.94) and 0.74 (95% c.i.: 0.62-0.87) for heterozygotes and rare homozygotes, respectively, compared with common homozygotes; P(trend) = 1.1 x 10(-7)) and weaker evidence for TGFB1 L10P (OR = 1.07 (95% c.i.: 1.02-1.13) and 1.16 (95% c.i.: 1.08-1.25), respectively; P(trend) = 2.8 x 10(-5)). These results demonstrate that common breast cancer susceptibility alleles with small effects on risk can be identified, given sufficiently powerful studies.

A genome-wide copy number association study of osteoporotic fractures points to the 6p25.1 locus.

BACKGROUND: Osteoporosis is a systemic skeletal disease characterised by reduced bone mineral density and increased susceptibility to fracture; these traits are highly heritable. Both common and rare copy number variants (CNVs) potentially affect the function of genes and may influence disease risk. AIM: To identify CNVs associated with osteoporotic bone fracture risk. METHOD: We performed a genome-wide CNV association study in 5178 individuals from a prospective cohort in the Netherlands, including 809 osteoporotic fracture cases, and performed in silico lookups and de novo genotyping to replicate in several independent studies. RESULTS: A rare (population prevalence 0.14%, 95% CI 0.03% to 0.24%) 210 kb deletion located on chromosome 6p25.1 was associated with the risk of fracture (OR 32.58, 95% CI 3.95 to 1488.89; p = 8.69 × 10(-5)). We performed an in silico meta-analysis in four studies with CNV microarray data and the association with fracture risk was replicated (OR 3.11, 95% CI 1.01 to 8.22; p = 0.02). The prevalence of this deletion showed geographic diversity, being absent in additional samples from Australia, Canada, Poland, Iceland, Denmark, and Sweden, but present in the Netherlands (0.34%), Spain (0.33%), USA (0.23%), England (0.15%), Scotland (0.10%), and Ireland (0.06%), with insufficient evidence for association with fracture risk. CONCLUSIONS: These results suggest that deletions in the 6p25.1 locus may predispose to higher risk of fracture in a subset of populations of European origin; larger and geographically restricted studies will be needed to confirm this regional association. This is a first step towards the evaluation of the role of rare CNVs in osteoporosis.

Genes contributing to pain sensitivity in the normal population: an exome sequencing study.

Sensitivity to pain varies considerably between individuals and is known to be heritable. Increased sensitivity to experimental pain is a risk factor for developing chronic pain, a common and debilitating but poorly understood symptom. To understand mechanisms underlying pain sensitivity and to search for rare gene variants (MAF<5%) influencing pain sensitivity, we explored the genetic variation in individuals' responses to experimental pain. Quantitative sensory testing to heat pain was performed in 2,500 volunteers from TwinsUK (TUK): exome sequencing to a depth of 70× was carried out on DNA from singletons at the high and low ends of the heat pain sensitivity distribution in two separate subsamples. Thus in TUK1, 101 pain-sensitive and 102 pain-insensitive were examined, while in TUK2 there were 114 and 96 individuals respectively. A combination of methods was used to test the association between rare variants and pain sensitivity, and the function of the genes identified was explored using network analysis. Using causal reasoning analysis on the genes with different patterns of SNVs by pain sensitivity status, we observed a significant enrichment of variants in genes of the angiotensin pathway (Bonferroni corrected p = 3.8×10(-4)). This pathway is already implicated in animal models and human studies of pain, supporting the notion that it may provide fruitful new targets in pain management. The approach of sequencing extreme exome variation in normal individuals has provided important insights into gene networks mediating pain sensitivity in humans and will be applicable to other common complex traits.

The FAIR Cookbook - the essential resource for and by FAIR doers.

The notion that data should be Findable, Accessible, Interoperable and Reusable, according to the FAIR Principles, has become a global norm for good data stewardship and a prerequisite for reproducibility. Nowadays, FAIR guides data policy actions and professional practices in the public and private sectors. Despite such global endorsements, however, the FAIR Principles are aspirational, remaining elusive at best, and intimidating at worst. To address the lack of practical guidance, and help with capability gaps, we developed the FAIR Cookbook, an open, online resource of hands-on recipes for "FAIR doers" in the Life Sciences. Created by researchers and data managers professionals in academia, (bio)pharmaceutical companies and information service industries, the FAIR Cookbook covers the key steps in a FAIRification journey, the levels and indicators of FAIRness, the maturity model, the technologies, the tools and the standards available, as well as the skills required, and the challenges to achieve and improve data FAIRness. Part of the ELIXIR ecosystem, and recommended by funders, the FAIR Cookbook is open to contributions of new recipes.

Tackling the translational challenges of multi-omics research in the realm of European personalised medicine: A workshop report.

Personalised medicine (PM) presents a great opportunity to improve the future of individualised healthcare. Recent advances in -omics technologies have led to unprecedented efforts characterising the biology and molecular mechanisms that underlie the development and progression of a wide array of complex human diseases, supporting further development of PM. This article reflects the outcome of the 2021 EATRIS-Plus Multi-omics Stakeholder Group workshop organised to 1) outline a global overview of common promises and challenges that key European stakeholders are facing in the field of multi-omics research, 2) assess the potential of new technologies, such as artificial intelligence (AI), and 3) establish an initial dialogue between key initiatives in this space. Our focus is on the alignment of agendas of European initiatives in multi-omics research and the centrality of patients in designing solutions that have the potential to advance PM in long-term healthcare strategies.

Association of ESR1 gene tagging SNPs with breast cancer risk.

We have conducted a three-stage, comprehensive single nucleotide polymorphism (SNP)-tagging association study of ESR1 gene variants (SNPs) in more than 55,000 breast cancer cases and controls from studies within the Breast Cancer Association Consortium (BCAC). No large risks or highly significant associations were revealed. SNP rs3020314, tagging a region of ESR1 intron 4, is associated with an increase in breast cancer susceptibility with a dominant mode of action in European populations. Carriers of the c-allele have an odds ratio (OR) of 1.05 [95% Confidence Intervals (CI) 1.02-1.09] relative to t-allele homozygotes, P = 0.004. There is significant heterogeneity between studies, P = 0.002. The increased risk appears largely confined to oestrogen receptor-positive tumour risk. The region tagged by SNP rs3020314 contains sequence that is more highly conserved across mammalian species than the rest of intron 4, and it may subtly alter the ratio of two mRNA splice forms.

Seq4SNPs: new software for retrieval of multiple, accurately annotated DNA sequences, ready formatted for SNP assay design.

BACKGROUND: In moderate-throughput SNP genotyping there was a gap in the workflow, between choosing a set of SNPs and submitting their sequences to proprietary assay design software, which was not met by existing software. Retrieval and formatting of sequences flanking each SNP, prior to assay design, becomes rate-limiting for more than about ten SNPs, especially if annotated for repetitive regions and adjacent variations. We routinely process up to 50 SNPs at once. IMPLEMENTATION: We created Seq4SNPs, a web-based, walk-away software that can process one to several hundred SNPs given rs numbers as input. It outputs a file of fully annotated sequences formatted for one of three proprietary design softwares: TaqMan's Primer-By-Design FileBuilder, Sequenom's iPLEX or SNPstream's Autoprimer, as well as unannotated fasta sequences. We found genotyping assays to be inhibited by repetitive sequences or the presence of additional variations flanking the SNP under test, and in multiplexes, repetitive sequence flanking one SNP adversely affects multiple assays. Assay design software programs avoid such regions if the input sequences are appropriately annotated, so we used Seq4SNPs to provide suitably annotated input sequences, and improved our genotyping success rate. Adjacent SNPs can also be avoided, by annotating sequences used as input for primer design. CONCLUSION: The accuracy of annotation by Seq4SNPs is significantly better than manual annotation (P < 1e-5).Using Seq4SNPs to incorporate all annotation for additional SNPs and repetitive elements into sequences, for genotyping assay designer software, minimizes assay failure at the design stage, reducing the cost of genotyping. Seq4SNPs provides a rapid route for replacement of poor test SNP sequences. We routinely use this software for assay sequence preparation. Seq4SNPs is available as a service at (http://moya.srl.cam.ac.uk/oncology/bio/s4shome.html) and (http://moya.srl.cam.ac.uk/cgi-bin/oncology/srl/ncbi/seq4snp1.pl), currently for human SNPs, but easily extended to include any species in dbSNP.

Publisher Correction: Federated discovery and sharing of genomic data using Beacons.

In the version of this article initially published, Lena Dolman's second affiliation was given as Wellcome Trust Sanger Institute, Wellcome Genome Campus, Hinxton, Cambridge, UK. The correct second affiliation is Ontario Institute for Cancer Research, Toronto, Ontario, Canada. The error has been corrected in the HTML and PDF versions of the article.

Large-scale analysis of association between polymorphisms in the transforming growth factor beta 1 gene (TGFB1) and osteoporosis: the GENOMOS study.

INTRODUCTION: The TGFB1 gene which encodes transforming growth factor beta 1, is a strong candidate for susceptibility to osteoporosis and several studies have reported associations between bone mineral density (BMD), osteoporotic fractures and polymorphisms of TGFB1, although these studies have yielded conflicting results. METHODS: We investigated associations between TGFB1 polymorphisms and BMD and fracture in the GENOMOS study: a prospective multicenter study involving 10 European research studies including a total of 28,924 participants. Genotyping was conducted for known TGFB1 polymorphisms at the following sites: G-1639-A (G-800-A, rs1800468), C-1348-T (C-509-T, rs1800469), T29-C (Leu10Pro, rs1982073), G74-C (Arg25Pro, rs1800471) and C788-T (Thr263Ile, rs1800472). These polymorphisms were genotyped prospectively and methodology was standardized across research centers. Genotypes and haplotypes were related to BMD at the lumbar sine and femoral neck and fractures. RESULTS: There were no significant differences in either women or men at either skeletal site for any of the examined polymorphisms with the possible exception of a weak association with reduced BMD (-12 mg/cm2) in men with the T-1348 allele (p<0.05). None of the haplotypes was associated with BMD and none of the polymorphisms or haplotypes significantly affected overall risk of fractures, however, the odds ratio for incident vertebral fracture in carriers of the rare T788 allele was 1.64 (95% CI: 1.09-2.64), p<0.05. CONCLUSIONS: This study indicates that polymorphic variation in the TGFB1 gene does not play a major role in regulating BMD or susceptibility to fractures. The weak associations we observed between the C-1348-T and lumbar spine BMD in men and between C788-T and risk of incident vertebral fractures are of interest but could be chance findings and will need replication in future studies.

Large-scale analysis of association between LRP5 and LRP6 variants and osteoporosis.

CONTEXT: Mutations in the low-density lipoprotein receptor-related protein 5 (LRP5) gene cause rare syndromes characterized by altered bone mineral density (BMD). More common LRP5 variants may affect osteoporosis risk in the general population. OBJECTIVE: To generate large-scale evidence on whether 2 common variants of LRP5 (Val667Met, Ala1330Val) and 1 variant of LRP6 (Ile1062Val) are associated with BMD and fracture risk. DESIGN AND SETTING: Prospective, multicenter, collaborative study of individual-level data on 37,534 individuals from 18 participating teams in Europe and North America. Data were collected between September 2004 and January 2007; analysis of the collected data was performed between February and May 2007. Bone mineral density was assessed by dual-energy x-ray absorptiometry. Fractures were identified via questionnaire, medical records, or radiographic documentation; incident fracture data were available for some cohorts, ascertained via routine surveillance methods, including radiographic examination for vertebral fractures. MAIN OUTCOME MEASURES: Bone mineral density of the lumbar spine and femoral neck; prevalence of all fractures and vertebral fractures. RESULTS: The Met667 allele of LRP5 was associated with reduced lumbar spine BMD (n = 25,052 [number of participants with available data]; 20-mg/cm2 lower BMD per Met667 allele copy; P = 3.3 x 10(-8)), as was the Val1330 allele (n = 24,812; 14-mg/cm2 lower BMD per Val1330 copy; P = 2.6 x 10(-9)). Similar effects were observed for femoral neck BMD, with a decrease of 11 mg/cm2 (P = 3.8 x 10(-5)) and 8 mg/cm2 (P = 5.0 x 10(-6)) for the Met667 and Val1330 alleles, respectively (n = 25 193). Findings were consistent across studies for both LRP5 alleles. Both alleles were associated with vertebral fractures (odds ratio [OR], 1.26; 95% confidence interval [CI], 1.08-1.47 for Met667 [2001 fractures among 20 488 individuals] and OR, 1.12; 95% CI, 1.01-1.24 for Val1330 [1988 fractures among 20,096 individuals]). Risk of all fractures was also increased with Met667 (OR, 1.14; 95% CI, 1.05-1.24 per allele [7876 fractures among 31,435 individuals)]) and Val1330 (OR, 1.06; 95% CI, 1.01-1.12 per allele [7802 fractures among 31 199 individuals]). Effects were similar when adjustments were made for age, weight, height, menopausal status, and use of hormone therapy. Fracture risks were partly attenuated by adjustment for BMD. Haplotype analysis indicated that Met667 and Val1330 variants both independently affected BMD. The LRP6 Ile1062Val polymorphism was not associated with any osteoporosis phenotype. All aforementioned associations except that between Val1330 and all fractures and vertebral fractures remained significant after multiple-comparison adjustments. CONCLUSIONS: Common LRP5 variants are consistently associated with BMD and fracture risk across different white populations. The magnitude of the effect is modest. LRP5 may be the first gene to reach a genome-wide significance level (a conservative level of significance [herein, unadjusted P < 10(-7)] that accounts for the many possible comparisons in the human genome) for a phenotype related to osteoporosis.

Registered access: authorizing data access.

The Global Alliance for Genomics and Health (GA4GH) proposes a data access policy model-"registered access"-to increase and improve access to data requiring an agreement to basic terms and conditions, such as the use of DNA sequence and health data in research. A registered access policy would enable a range of categories of users to gain access, starting with researchers and clinical care professionals. It would also facilitate general use and reuse of data but within the bounds of consent restrictions and other ethical obligations. In piloting registered access with the Scientific Demonstration data sharing projects of GA4GH, we provide additional ethics, policy and technical guidance to facilitate the implementation of this access model in an international setting.

Common variants in the ATM, BRCA1, BRCA2, CHEK2 and TP53 cancer susceptibility genes are unlikely to increase breast cancer risk.

INTRODUCTION: Certain rare, familial mutations in the ATM, BRCA1, BRCA2, CHEK2 or TP53 genes increase susceptibility to breast cancer but it has not, until now, been clear whether common polymorphic variants in the same genes also increase risk. METHODS: We have attempted a comprehensive, single nucleotide polymorphism (SNP)- and haplotype-tagging association study on each of these five genes in up to 4,474 breast cancer cases from the British, East Anglian SEARCH study and 4,560 controls from the EPIC-Norfolk study, using a two-stage study design. Nine tag SNPs were genotyped in ATM, together with five in BRCA1, sixteen in BRCA2, ten in CHEK2 and five in TP53, with the aim of tagging all other known, common variants. SNPs generating the common amino acid substitutions were specifically forced into the tagging set for each gene. RESULTS: No significant breast cancer associations were detected with any individual or combination of tag SNPs. CONCLUSION: It is unlikely that there are any other common variants in these genes conferring measurably increased risks of breast cancer in our study population.