An interconnected data infrastructure to support large-scale rare disease research.

The Solve-RD project brings together clinicians, scientists, and patient representatives from 51 institutes spanning 15 countries to collaborate on genetically diagnosing ("solving") rare diseases (RDs). The project aims to significantly increase the diagnostic success rate by co-analyzing data from thousands of RD cases, including phenotypes, pedigrees, exome/genome sequencing, and multiomics data. Here we report on the data infrastructure devised and created to support this co-analysis. This infrastructure enables users to store, find, connect, and analyze data and metadata in a collaborative manner. Pseudonymized phenotypic and raw experimental data are submitted to the RD-Connect Genome-Phenome Analysis Platform and processed through standardized pipelines. Resulting files and novel produced omics data are sent to the European Genome-Phenome Archive, which adds unique file identifiers and provides long-term storage and controlled access services. MOLGENIS "RD3" and Café Variome "Discovery Nexus" connect data and metadata and offer discovery services, and secure cloud-based "Sandboxes" support multiparty data analysis. This successfully deployed and useful infrastructure design provides a blueprint for other projects that need to analyze large amounts of heterogeneous data.

Common conditions of use elements. Atomic concepts for consistent and effective information governance.

Myriad policy, ethical and legal considerations underpin the sharing of biological resources, implying the need for standardised and yet flexible ways to digitally represent diverse 'use conditions'. We report a core lexicon of terms that are atomic, non-directional 'concepts of use', called Common Conditions of use Elements. This work engaged biobanks and registries relevant to the European Joint Programme for Rare Diseases and aimed to produce a lexicon that would have generalised utility. Seventy-six concepts were initially identified from diverse real-world settings, and via iterative rounds of deliberation and user-testing these were optimised and condensed down to 20 items. To validate utility, support software and training information was provided to biobanks and registries who were asked to create Sharing Policy Profiles. This succeeded and involved adding standardised directionality and scope annotations to the employed terms. The addition of free-text parameters was also explored. The approach is now being adopted by several real-world projects, enabling this standard to evolve progressively into a universal basis for representing and managing conditions of use.

Getting your DUCs in a row - standardising the representation of Digital Use Conditions.

Improving patient care and advancing scientific discovery requires responsible sharing of research data, healthcare records, biosamples, and biomedical resources that must also respect applicable use conditions. Defining a standard to structure and manage these use conditions is a complex and challenging task. This is exemplified by a near unlimited range of asset types, a high variability of applicable conditions, and differing applications at the individual or collective level. Furthermore, the specifics and granularity required are likely to vary depending on the ultimate contexts of use. All these factors confound alignment of institutional missions, funding objectives, regulatory and technical requirements to facilitate effective sharing. The presented work highlights the complexity and diversity of the problem, reviews the current state of the art, and emphasises the need for a flexible and adaptable approach. We propose Digital Use Conditions (DUC) as a framework that addresses these needs by leveraging existing standards, striking a balance between expressiveness versus ambiguity, and considering the breadth of applicable information with their context of use.

Genome-wide association study of treatment-related toxicity two years following radiotherapy for breast cancer.

BACKGROUND AND PURPOSE: Up to a quarter of breast cancer patients treated by surgery and radiotherapy experience clinically significant toxicity. If patients at high risk of adverse effects could be identified at diagnosis, their treatment could be tailored accordingly. This study was designed to identify common single nucleotide polymorphisms (SNPs) associated with toxicity two years following whole breast radiotherapy. MATERIALS AND METHODS: A genome-wide association study (GWAS) was performed in 1,640 breast cancer patients with complete SNP, clinical, treatment and toxicity data, recruited across 18 European and US centres into the prospective REQUITE cohort study. Toxicity data (CTCAE v4.0) were collected at baseline, end of radiotherapy, and annual follow-up. A total of 7,097,340 SNPs were tested for association with the residuals of toxicity endpoints, adjusted for clinical, treatment co-variates and population substructure. RESULTS: Quantile-quantile plots showed more associations with toxicity above the p < 5 × 10-5 level than expected by chance. Eight SNPs reached genome-wide significance. Nipple retraction grade ≥ 2 was associated with the rs188287402 variant (p = 2.80 × 10-8), breast oedema grade ≥ 2 with rs12657177 (p = 1.12 × 10-10), rs75912034 (p = 1.12 × 10-10), rs145328458 (p = 1.06 × 10-9) and rs61966612 (p = 1.23 × 10-9), induration grade ≥ 2 with rs77311050 (p = 2.54 × 10-8) and rs34063419 (p = 1.21 × 10-8), and arm lymphoedema grade ≥ 1 with rs643644 (p = 3.54 × 10-8). Heritability estimates across significant endpoints ranged from 25% to 39%. Our study did not replicate previously reported SNPs associated with breast radiation toxicity at the pre-specified significance level. CONCLUSIONS: This GWAS for long-term breast radiation toxicity provides further evidence for significant association of common SNPs with distinct toxicity endpoints.

Beacon v2 and Beacon networks: A "lingua franca" for federated data discovery in biomedical genomics, and beyond.

Beacon is a basic data discovery protocol issued by the Global Alliance for Genomics and Health (GA4GH). The main goal addressed by version 1 of the Beacon protocol was to test the feasibility of broadly sharing human genomic data, through providing simple "yes" or "no" responses to queries about the presence of a given variant in datasets hosted by Beacon providers. The popularity of this concept has fostered the design of a version 2, that better serves real-world requirements and addresses the needs of clinical genomics research and healthcare, as assessed by several contributing projects and organizations. Particularly, rare disease genetics and cancer research will benefit from new case level and genomic variant level requests and the enabling of richer phenotype and clinical queries as well as support for fuzzy searches. Beacon is designed as a "lingua franca" to bridge data collections hosted in software solutions with different and rich interfaces. Beacon version 2 works alongside popular standards like Phenopackets, OMOP, or FHIR, allowing implementing consortia to return matches in beacon responses and provide a handover to their preferred data exchange format. The protocol is being explored by other research domains and is being tested in several international projects.

Variation in human herpesvirus 6B telomeric integration, excision, and transmission between tissues and individuals.

Human herpesviruses 6A and 6B (HHV-6A/6B) are ubiquitous pathogens that persist lifelong in latent form and can cause severe conditions upon reactivation. They are spread by community-acquired infection of free virus (acqHHV6A/6B) and by germline transmission of inherited chromosomally integrated HHV-6A/6B (iciHHV-6A/6B) in telomeres. We exploited a hypervariable region of the HHV-6B genome to investigate the relationship between acquired and inherited virus and revealed predominantly maternal transmission of acqHHV-6B in families. Remarkably, we demonstrate that some copies of acqHHV-6B in saliva from healthy adults gained a telomere, indicative of integration and latency, and that the frequency of viral genome excision from telomeres in iciHHV-6B carriers is surprisingly high and varies between tissues. In addition, newly formed short telomeres generated by partial viral genome release are frequently lengthened, particularly in telomerase-expressing pluripotent cells. Consequently, iciHHV-6B carriers are mosaic for different iciHHV-6B structures, including circular extra-chromosomal forms that have the potential to reactivate. Finally, we show transmission of an HHV-6B strain from an iciHHV-6B mother to her non-iciHHV-6B son. Altogether, we demonstrate that iciHHV-6B can readily transition between telomere-integrated and free virus forms.

Array-based dynamic allele specific hybridization (Array-DASH): Optimization-free microarray processing for multiple simultaneous genomic assays.

We report proof-of-principle experiments regarding a dynamic microarray protocol enabling accurate and semi-quantitative DNA analysis for re-sequencing, fingerprinting and genotyping. Single-stranded target molecules hybridise to surface-bound probes during initial gradual cooling with high-fidelity. Real-time tracking of target denaturation (via fluorescence) during a 'dynamic' gradual heating phase permits 'melt-curve' analysis. The probe most closely matching the target sequence is identified based on the highest melting temperature. We demonstrated a >99% re-sequencing accuracy and a potential detection rate of 1% for SNPs. Experiments employing Hypericum ribosomal ITS regions and HIV genomes illustrated a reliable detection level of 5% plus simultaneous re-sequencing and genotyping. Such performance suggests a range of potential real-world applications involving rapid sequence interrogation, for example, in the Covid-19 pandemic. Guidance is offered towards the development of a commercial platform and dedicated software required to bring this technique into mainstream science.

Prescreening for European Prevention of Alzheimer Dementia (EPAD) trial-ready cohort: impact of AD risk factors and recruitment settings.

BACKGROUND: Recruitment is often a bottleneck in secondary prevention trials in Alzheimer disease (AD). Furthermore, screen-failure rates in these trials are typically high due to relatively low prevalence of AD pathology in individuals without dementia, especially among cognitively unimpaired. Prescreening on AD risk factors may facilitate recruitment, but the efficiency will depend on how these factors link to participation rates and AD pathology. We investigated whether common AD-related factors predict trial-ready cohort participation and amyloid status across different prescreen settings. METHODS: We monitored the prescreening in four cohorts linked to the European Prevention of Alzheimer Dementia (EPAD) Registry (n = 16,877; mean ± SD age = 64 ± 8 years). These included a clinical cohort, a research in-person cohort, a research online cohort, and a population-based cohort. Individuals were asked to participate in the EPAD longitudinal cohort study (EPAD-LCS), which serves as a trial-ready cohort for secondary prevention trials. Amyloid positivity was measured in cerebrospinal fluid as part of the EPAD-LCS assessment. We calculated participation rates and numbers needed to prescreen (NNPS) per participant that was amyloid-positive. We tested if age, sex, education level, APOE status, family history for dementia, memory complaints or memory scores, previously collected in these cohorts, could predict participation and amyloid status. RESULTS: A total of 2595 participants were contacted for participation in the EPAD-LCS. Participation rates varied by setting between 3 and 59%. The NNPS were 6.9 (clinical cohort), 7.5 (research in-person cohort), 8.4 (research online cohort), and 88.5 (population-based cohort). Participation in the EPAD-LCS (n = 413 (16%)) was associated with lower age (odds ratio (OR) age = 0.97 [0.95-0.99]), high education (OR = 1.64 [1.23-2.17]), male sex (OR = 1.56 [1.19-2.04]), and positive family history of dementia (OR = 1.66 [1.19-2.31]). Among participants in the EPAD-LCS, amyloid positivity (33%) was associated with higher age (OR = 1.06 [1.02-1.10]) and APOE ɛ4 allele carriership (OR = 2.99 [1.81-4.94]). These results were similar across prescreen settings. CONCLUSIONS: Numbers needed to prescreen varied greatly between settings. Understanding how common AD risk factors link to study participation and amyloid positivity is informative for recruitment strategy of studies on secondary prevention of AD.

Do infants with transposition of the great arteries born outside a specialist centre have different outcomes?

BACKGROUND: Infants born with undiagnosed transposition of the great arteries continue to be born in district general hospitals despite the improvements made in antenatal scanning. Evidence indicates improved outcomes with early definitive treatment after birth, hence the recommendation of delivery in a tertiary centre. The role of specialist paediatric and neonatal transport teams, to advise, stabilise, and transport the infants to a tertiary centre in a timely manner, is critical for those infants born in a district general hospital. This pilot study aims to compare outcomes between infants born in district general hospitals and those who were born in a tertiary maternity unit in South West England and South Wales. METHODS: This was a secondary data analysis of data collected from the local Paediatric Intensive Care Audit Network and the local transport database. Infants born with a confirmed diagnosis of transposition of the great arteries, that required an arterial switch operation as the definitive procedure between April, 2012 and March 2018 were included. RESULTS: Forty-five infants with a confirmed diagnosis of transposition of the great arteries were included. Statistical analysis demonstrated there were no significant differences in the time to balloon atrial septostomy (p = 0.095), time to arterial switch operation (p = 0.461), length of paediatric ICU stay (p = 0.353), and hospital stay (p = 0.095) or mortality between the two groups. CONCLUSIONS: We found no significant differences in outcomes between infants delivered outside the specialist centre, who were transferred in by a specialist team.

Responsible sharing of biomedical data and biospecimens via the "Automatable Discovery and Access Matrix" (ADA-M).

Given the data-rich nature of modern biomedical research, there is a pressing need for a systematic, structured, computer-readable way to capture, communicate, and manage sharing rules that apply to biomedical resources. This is essential for responsible recording, versioning, communication, querying, and actioning of resource sharing plans. However, lack of a common "information model" for rules and conditions that govern the sharing of materials, methods, software, data, and knowledge creates a fundamental barrier. Without this, it can be virtually impossible for Research Ethics Committees (RECs), Institutional Review Boards (IRBs), Data Access Committees (DACs), biobanks, and end users to confidently track, manage, and interpret applicable legal and ethical requirements. This raises costs and burdens of data stewardship and decreases efficient and responsible access to data, biospecimens, and other resources. To address this, the GA4GH and IRDiRC organizations sponsored the creation of the Automatable Discovery and Access Matrix (ADA-M, read simply as "Adam"). ADA-M is a comprehensive information model that provides the basis for producing structured metadata "Profiles" of regulatory conditions, thereby enabling efficient application of those conditions across regulatory spheres. Widespread use of ADA-M will aid researchers in globally searching and prescreening potential data and/or biospecimen resources for compatibility with their research plans in a responsible and efficient manner, increasing likelihood of timely DAC approvals while also significantly reducing time and effort DACs, RECs, and IRBs spend evaluating resource requests and research proposals. Extensive online documentation, software support, video guides, and an Application Programming Interface (API) for ADA-M have been made available.

European Prevention of Alzheimer's Dementia Registry: Recruitment and prescreening approach for a longitudinal cohort and prevention trials.

INTRODUCTION: It is a challenge to find participants for Alzheimer's disease (AD) prevention trials within a short period of time. The European Prevention of Alzheimer's Dementia Registry (EPAD) aims to facilitate recruitment by preselecting subjects from ongoing cohort studies. This article introduces this novel approach. METHODS: A virtual registry, with access to risk factors and biomarkers for AD through minimal data sets of ongoing cohort studies, was set up. RESULTS: To date, ten cohorts have been included in the EPAD. Around 2500 participants have been selected, using variables associated with the risk for AD. Of these, 15% were already recruited in the EPAD longitudinal cohort study, which serves as a trial readiness cohort. DISCUSSION: This study demonstrates that a virtual registry can be used for the preselection of participants for AD studies.

Inherited Chromosomally Integrated Human Herpesvirus 6 Genomes Are Ancient, Intact, and Potentially Able To Reactivate from Telomeres.

The genomes of human herpesvirus 6A (HHV-6A) and HHV-6B have the capacity to integrate into telomeres, the essential capping structures of chromosomes that play roles in cancer and ageing. About 1% of people worldwide are carriers of chromosomally integrated HHV-6 (ciHHV-6), which is inherited as a genetic trait. Understanding the consequences of integration for the evolution of the viral genome, for the telomere, and for the risk of disease associated with carrier status is hampered by a lack of knowledge about ciHHV-6 genomes. Here, we report an analysis of 28 ciHHV-6 genomes and show that they are significantly divergent from the few modern nonintegrated HHV-6 strains for which complete sequences are currently available. In addition, ciHHV-6B genomes in Europeans are more closely related to each other than to ciHHV-6B genomes from China and Pakistan, suggesting regional variation of the trait. Remarkably, at least one group of European ciHHV-6B carriers has inherited the same ciHHV-6B genome, integrated in the same telomere allele, from a common ancestor estimated to have existed 24,500 ± 10,600 years ago. Despite the antiquity of some, and possibly most, germ line HHV-6 integrations, the majority of ciHHV-6B (95%) and ciHHV-6A (72%) genomes contain a full set of intact viral genes and therefore appear to have the capacity for viral gene expression and full reactivation.IMPORTANCE Inheritance of HHV-6A or HHV-6B integrated into a telomere occurs at a low frequency in most populations studied to date, but its characteristics are poorly understood. However, stratification of ciHHV-6 carriers in modern populations due to common ancestry is an important consideration for genome-wide association studies that aim to identify disease risks for these people. Here, we present full sequence analysis of 28 ciHHV-6 genomes and show that ciHHV-6B in many carriers with European ancestry most likely originated from ancient integration events in a small number of ancestors. We propose that ancient ancestral origins for ciHHV-6A and ciHHV-6B are also likely in other populations. Moreover, despite their antiquity, all of the ciHHV-6 genomes appear to retain the capacity to express viral genes, and most are predicted to be capable of full viral reactivation. These discoveries represent potentially important considerations in immunocompromised patients, in particular in organ transplantation and in stem cell therapy.

Copy number variation of scavenger-receptor cysteine-rich domains within DMBT1 and Crohn's disease.

Previous work has shown that the gene DMBT1, which encodes a large secreted epithelial glycoprotein known as salivary agglutinin, gp340, hensin or muclin, is an innate immune defence protein that binds bacteria. A deletion variant of DMBT1 has been previously associated with Crohn's disease, and a DMBT1(-/-) knockout mouse has increased levels of colitis induced by dextran sulphate. DMBT1 has a complex copy number variable structure, with two, independent, rapidly mutating copy number variable regions, called CNV1 and CNV2. Because the copy number variable regions are predicted to affect the number of bacteria-binding domains, different alleles may alter host-microbe interactions in the gut. Our aim was to investigate the role of this complex variation in susceptibility to Crohn's disease by assessing the previously reported association. We analysed the association of both copy number variable regions with presence of Crohn's disease, and its severity, on three case-control cohorts. We also reanalysed array comparative genomic hybridisation data (aCGH) from a large case-control cohort study for both copy number variable regions. We found no association with a linear increase in copy number, nor when the CNV1 is regarded as presence or absence of a deletion allele. Taken together, we show that the DMBT1 CNV does not affect susceptibility to Crohn's disease, at least in Northern Europeans.

New technologies for DNA analysis--a review of the READNA Project.

The REvolutionary Approaches and Devices for Nucleic Acid analysis (READNA) project received funding from the European Commission for 41/2 years. The objectives of the project revolved around technological developments in nucleic acid analysis. The project partners have discovered, created and developed a huge body of insights into nucleic acid analysis, ranging from improvements and implementation of current technologies to the most promising sequencing technologies that constitute a 3(rd) and 4(th) generation of sequencing methods with nanopores and in situ sequencing, respectively.

Cafe Variome: general-purpose software for making genotype-phenotype data discoverable in restricted or open access contexts.

Biomedical data sharing is desirable, but problematic. Data "discovery" approaches-which establish the existence rather than the substance of data-precisely connect data owners with data seekers, and thereby promote data sharing. Cafe Variome (http://www.cafevariome.org) was therefore designed to provide a general-purpose, Web-based, data discovery tool that can be quickly installed by any genotype-phenotype data owner, or network of data owners, to make safe or sensitive content appropriately discoverable. Data fields or content of any type can be accommodated, from simple ID and label fields through to extensive genotype and phenotype details based on ontologies. The system provides a "shop window" in front of data, with main interfaces being a simple search box and a powerful "query-builder" that enable very elaborate queries to be formulated. After a successful search, counts of records are reported grouped by "openAccess" (data may be directly accessed), "linkedAccess" (a source link is provided), and "restrictedAccess" (facilitated data requests and subsequent provision of approved records). An administrator interface provides a wide range of options for system configuration, enabling highly customized single-site or federated networks to be established. Current uses include rare disease data discovery, patient matchmaking, and a Beacon Web service.

A 129-kb deletion on chromosome 12 confers substantial protection against rheumatoid arthritis, implicating the gene SLC2A3.

We describe a copy-number variant (CNV) for which deletion alleles confer a protective affect against rheumatoid arthritis (RA). This CNV reflects net unit deletions and expansions to a normal two-unit tandem duplication located on human chr12p13.31, a region with conserved synteny to the rat RA susceptibility quantitative trait loci Oia2. Genotyping, using the paralogue ratio test and SNP intensity data, in Swedish samples (2,403 cases, 1,269 controls) showed that the frequency of deletion variants is significantly lower in cases (P = 0.0012, OR = 0.442 [95%CI 0.258-0.755]). Reduced frequencies of deletion variants were also seen in replication materials comprising 9,201 UK samples (1,846 cases, 7,355 controls) and 2,963 US samples (906 controls, 1,967 cases) (Mantel-Haenszel P = 0.036, OR = 0.559 [95%CI 0.323-0.966]). Combining the three datasets produces a Mantel-Haenszel OR of 0.497 (P < 0.0002). The deletion variant lacks 129-kb of DNA containing SLC2A3, NANOGP1, and SLC2A14. SLC2A3 encodes a high-affinity glucose transporter important in the immune response and chondrocyte metabolism, both key aspects of RA pathogenesis. The large effect size of this association, its potential relevance to other diseases in which SLC2A3 is implicated, and the possibility of targeting drugs to inhibit SLC2A3, argue for further examination of the genetics and the biology of this CNV.

Automated design of paralogue ratio test assays for the accurate and rapid typing of copy number variation.

MOTIVATION: Genomic copy number variation (CNV) can influence susceptibility to common diseases. High-throughput measurement of gene copy number on large numbers of samples is a challenging, yet critical, stage in confirming observations from sequencing or array Comparative Genome Hybridization (CGH). The paralogue ratio test (PRT) is a simple, cost-effective method of accurately determining copy number by quantifying the amplification ratio between a target and reference amplicon. PRT has been successfully applied to several studies analyzing common CNV. However, its use has not been widespread because of difficulties in assay design. RESULTS: We present PRTPrimer (www.prtprimer.org) software for automated PRT assay design. In addition to stand-alone software, the web site includes a database of pre-designed assays for the human genome at an average spacing of 6 kb and a web interface for custom assay design. Other reference genomes can also be analyzed through local installation of the software. The usefulness of PRTPrimer was tested within known CNV, and showed reproducible quantification. This software and database provide assays that can rapidly genotype CNV, cost-effectively, on a large number of samples and will enable the widespread adoption of PRT. AVAILABILITY: PRTPrimer is available in two forms: a Perl script (version 5.14 and higher) that can be run from the command line on Linux systems and as a service on the PRTPrimer web site (www.prtprimer.org).

A mechanistic basis for amplification differences between samples and between genome regions.

BACKGROUND: For many analytical methods the efficiency of DNA amplification varies across the genome and between samples. The most affected genome regions tend to correlate with high C + G content, however this relationship is complex and does not explain why the direction and magnitude of effects varies considerably between samples. RESULTS: Here, we provide evidence that sequence elements that are particularly high in C + G content can remain annealed even when aggressive melting conditions are applied. In turn, this behavior creates broader 'Thermodynamically Ultra-Fastened' (TUF) regions characterized by incomplete denaturation of the two DNA strands, so reducing amplification efficiency throughout these domains. CONCLUSIONS: This model provides a mechanistic explanation for why some genome regions are particularly difficult to amplify and assay in many procedures, and importantly it also explains inter-sample variability of this behavior. That is, DNA samples of varying quality will carry more or fewer nicks and breaks, and hence their intact TUF regions will have different lengths and so be differentially affected by this amplification suppression mechanism - with 'higher' quality DNAs being the most vulnerable. A major practical consequence of this is that inter-region and inter-sample variability can be largely overcome by employing routine fragmentation methods (e.g. sonication or restriction enzyme digestion) prior to sample amplification.

Mutations in FRMD7, a newly identified member of the FERM family, cause X-linked idiopathic congenital nystagmus.

Idiopathic congenital nystagmus is characterized by involuntary, periodic, predominantly horizontal oscillations of both eyes. We identified 22 mutations in FRMD7 in 26 families with X-linked idiopathic congenital nystagmus. Screening of 42 singleton cases of idiopathic congenital nystagmus (28 male, 14 females) yielded three mutations (7%). We found restricted expression of FRMD7 in human embryonic brain and developing neural retina, suggesting a specific role in the control of eye movement and gaze stability.

Genetic association of the serotonin transporter in pulmonary arterial hypertension.

RATIONALE: The bone morphogenetic receptor type II gene is the major genetic determinant for the inherited form of pulmonary arterial hypertension. However, deleterious mutations of this gene are not observed in the majority of subjects who develop the condition spontaneously and familial disease displays age- and sex-dependent penetrance, indicating the requirement for additional environmental and/or genetic modifiers for disease development. METHODS: We investigated polymorphic variation of the serotonin transporter gene, a biological candidate for predisposition to this vascular disorder. RESULTS: No significant evidence of association between alleles of the serotonin transporter gene and pulmonary hypertension was detected, nor did we observe a relationship with age of onset in familial and idiopathic disease. CONCLUSIONS: Variation of the serotonin transporter gene appears unlikely to confer significant susceptibility to pulmonary arterial hypertension. This study emphasizes the need for adequately powered cohorts for association analyses to identify not only genetic determinants of disease susceptibility but also inherited modifiers for disease development.