A Boolean algebra for genetic variants.

MOTIVATION: Beyond identifying genetic variants, we introduce a set of Boolean relations, which allows for a comprehensive classification of the relations of every pair of variants by taking all minimal alignments into account. We present an efficient algorithm to compute these relations, including a novel way of efficiently computing all minimal alignments within the best theoretical complexity bounds. RESULTS: We show that these relations are common, and many non-trivial, for variants of the CFTR gene in dbSNP. Ultimately, we present an approach for the storing and indexing of variants in the context of a database that enables efficient querying for all these relations. AVAILABILITY AND IMPLEMENTATION: A Python implementation is available at https://github.com/mutalyzer/algebra/tree/v0.2.0 as well as an interface at https://mutalyzer.nl/algebra.

Three-tiered EGFr domain risk stratification for individualized NOTCH3-small vessel disease prediction.

Cysteine-altering missense variants (NOTCH3cys) in one of the 34 epidermal growth-factor-like repeat (EGFr) domains of the NOTCH3 protein are the cause of NOTCH3-associated small vessel disease (NOTCH3-SVD). NOTCH3-SVD is highly variable, ranging from cerebral autosomal dominant arteriopathy with subcortical infarcts and leukoencephalopathy (CADASIL) at the severe end of the spectrum to non-penetrance. The strongest known NOTCH3-SVD modifier is NOTCH3cys variant position: NOTCH3cys variants located in EGFr domains 1-6 are associated with a more severe phenotype than NOTCH3cys variants located in EGFr domains 7-34. The objective of this study was to further improve NOTCH3-SVD genotype-based risk prediction by using relative differences in NOTCH3cys variant frequencies between large CADASIL and population cohorts as a starting point. Scientific CADASIL literature, cohorts and population databases were queried for NOTCH3cys variants. For each EGFr domain, the relative difference in NOTCH3cys variant frequency (NVFOR) was calculated using genotypes of 2574 CADASIL patients and 1647 individuals from population databases. Based on NVFOR cut-off values, EGFr domains were classified as either low (LR-EGFr), medium (MR-EGFr) or high risk (HR-EGFr). The clinical relevance of this new three-tiered EGFr risk classification was cross-sectionally validated by comparing SVD imaging markers and clinical outcomes between EGFr risk categories using a genotype-phenotype data set of 434 CADASIL patients and 1003 NOTCH3cys positive community-dwelling individuals. CADASIL patients and community-dwelling individuals harboured 379 unique NOTCH3cys variants. Nine EGFr domains were classified as an HR-EGFr, which included EGFr domains 1-6, but additionally also EGFr domains 8, 11 and 26. Ten EGFr domains were classified as MR-EGFr and 11 as LR-EGFr. In the population genotype-phenotype data set, HR-EGFr individuals had the highest risk of stroke [odds ratio (OR) = 10.81, 95% confidence interval (CI): 5.46-21.37], followed by MR-EGFr individuals (OR = 1.81, 95% CI: 0.84-3.88) and LR-EGFr individuals (OR = 1 [reference]). MR-EGFr individuals had a significantly higher normalized white matter hyperintensity volume (nWMHv; P = 0.005) and peak width of skeletonized mean diffusivity (PSMD; P = 0.035) than LR-EGFr individuals. In the CADASIL genotype-phenotype data set, HR-EGFr domains 8, 11 and 26 patients had a significantly higher risk of stroke (P = 0.002), disability (P = 0.041), nWMHv (P = 1.8 × 10-8), PSMD (P = 2.6 × 10-8) and lacune volume (P = 0.006) than MR-EGFr patients. SVD imaging marker load and clinical outcomes were similar between HR-EGFr 1-6 patients and HR-EGFr 8, 11 and 26 patients. NVFOR was significantly associated with vascular NOTCH3 aggregation load (P = 0.006), but not with NOTCH3 signalling activity (P = 0.88). In conclusion, we identified three clinically distinct NOTCH3-SVD EGFr risk categories based on NFVOR cut-off values, and identified three additional HR-EGFr domains located outside of EGFr domains 1-6. This EGFr risk classification will provide an important key to individualized NOTCH3-SVD disease prediction.

Legal constraints on genetic data processing in European grids.

European laws on privacy and data security are not explicit about the storage and processing of genetic data. Especially whole-genome data is identifying and contains a lot of personal information. Is processing of such data allowed in computing grids? To find out, we looked at legal precedents in related fields, current literature, and interviews with legal experts. We found that processing of genetic data is only allowed on distributed systems with specific security measures, both technical and organizational. Informed consent, although important, offers no substitute for such requirements.

Provenance for distributed biomedical workflow execution.

Scientific research has become very data and compute intensive because of the progress in data acquisition and measurement devices, which is particularly true in Life Sciences. To cope with this deluge of data, scientists use distributed computing and storage infrastructures. The use of such infrastructures introduces by itself new challenges to the scientists in terms of proper and efficient use. Scientific workflow management systems play an important role in facilitating the use of the infrastructure by hiding some of its complexity. Although most scientific workflow management systems are provenance-aware, not all of them come with provenance functionality out of the box. In this paper we describe the improvement and integration of a provenance system into an e-infrastructure for biomedical research based on the MOTEUR workflow management system. The main contributions of the paper are: presenting an OPM implementation using relational database backend for the provenance store, providing an e-infrastructure with a comprehensive provenance system, defining a generic approach to provenance implementation, potentially suitable for other workflow systems and application domains and demonstrating the value of this system based on use cases presenting the provenance data through a user-friendly web interface.