Abstract
Motivation We are witnessing an enormous growth in the amount of molecular profiling (-omics) data. The integration of multi-omics data is challenging. Moreover, human multi-omics data may be privacy-sensitive and misused to de-anonymize and (re-)identify individuals. Hence, most data is kept in secure and protected silos. Therefore, it remains a challenge to reuse these data without infringing the privacy of the individuals from which the data were derived. Federated analysis of FAIR data is a privacy-preserving solution to make optimal use of these multi-omics data and transform them into actionable knowledge.
Results The Netherlands X-omics Initiative is a National Roadmap Large-Scale Research Infrastructure aiming for efficient integration of data generated within X-omics and external datasets. To facilitate this, we developed the FAIR Data Cube (FDCube), which adopts and applies the FAIR principles and helps researchers to create FAIR data and metadata, facilitate reuse of their data, and make their data analysis workflows transparent. The FDCube also meets security-by-design and privacy-by-design principles.
Availability https://github.com/Xomics/FAIRDataCube
Contact Xiaofeng.Liao{at}radboudumc.nl, Peter-Bram.tHoen{at}radboudumc.nl
Supplementary information Supplementary data are available at Bioinformatics online.
Introduction
It is now widely acknowledged that understanding the mechanisms underlying health and disease requires the concerted study of different molecular levels (DNA, RNA, proteins, metabolites). Moreover, a transition from static simplified views to dynamic comprehensive views on molecular pathways (encompassing (e.g. genomics, proteomics and metabolomics) is needed. Currently, this is not simple nor scalable. There is an increasing need to combine -omics data from different sources, but the data and their associated metadata are not always findable, accessible, interoperable, and reusable (FAIR) [31]. For that reason, the Netherlands X-omics Initiative has developed a multi-omics data infrastructure that facilitates FAIR-compliant multi-omics data storage and analysis. The proposed data infrastructure provides an analysis environment for (federated) data handling and analysis meeting the security-by-design and privacy-by-design principles.
This paper introduces our solution of integrated analysis on FAIR multi-omics data in decentralized databases. In the remainder of this paper, section 2 investigates existing work in this research direction. Section 3 presents the design and implementation of the FAIR Data Cube (FDCube) and section 4 showcases the use of FDCube in the Trusted World of Corona project[11]. Finally, section 5 discusses further developments.
Related work
There are several tools that aid researchers in managing research metadata in a FAIR manner, for instance the FAIR Data Station[23], the FAIR-in-a-box[5] approach, and the DataFAIRifier[1]. Most of these tools focus on the production of FAIR data, including ingestion, generation, and publication. For a more comprehensive coverage of FAIR processes including data management, data security, data exchange, and federated analysis, additional tools are required. For example, MOLGENIS is an open-source web-application covering the typical flow of human genomics data including data collection, management, analysis, visualization, and sharing, as well as offering support to make data FAIR[30]. MOLGENIS can be hosted on-site and stores the data locally in a PostgreSQL database. This offers all the advantages of a true database including a local access control system (in light of the European General Data Protection Regulation) with detailed data management.
The Personal Health Train (PHT)[14] concept is underlying a number of approaches for decentralised analysis of health-related data. The essence of the PHT approach is the analogy of a station representing the data source and a train representing the research question (or a computational request) visiting the data stations. Stations range from very large databases to small personal lockers containing the data of one person. Each station has its own set of house rules describing what a visiting ‘train’ is allowed to do with its data[14]. By moving trains towards stations rather than moving data, copying of data is avoided, data remains under complete control of the person or institute generating the data, and privacy concerns around data sharing are alleviated.
DataSHIELD[18]implements the idea of bringing algorithms to the data to ensure data privacy and security. DataSHIELD facilitates (co-)analysis of (harmonised) biomedical, healthcare and social-science data stored at one or multiple locations. The analysis requests are sent from a central analysis machine to several data-holding machines, which store the harmonised data to be co-analysed. The datasets are then analysed simultaneously, but in parallel. MOLGENIS developed a DataSHIELD implementation called Armadillo in its MOLGENIS suite.
Vantage6[22, 28] is a different implementation of the PHT concept. Vantage6 enables collaboration between multiple parties to participate in one or multiple studies across multiple data stations.
In terms of programming language, DataSHIELD restricts itself to a single language (R)[24] and to a pre-defined library of functions and algorithms. By contrast, using Vantage6, the researcher can pose a request to use their preferred programming language, as long as the language is supported by the targeted data station.
To advance and further build upon the currently available federated, FAIR solutions for the scientific community, we here present the FAIR Data Cube (FDCube) for public use under an open MIT license. In contrast to the more generic MOLGENIS Armadillo approach, FDCube contains specialised services for the analysis of multi-omics data. The FDCube is developed based on the principle that data should be “as open as possible and as closed as necessary” [17]. By incorporating a FAIR Data Point (FDP) component, the metadata can be as open as possible and made FAIR-at-the-source. By integrating a Vantage6 component, the data security/privacy can be ensured by collaborated federated analysis.
Result
The FDCube is a technological framework for the storage, analysis and integration of multi-omics data. The FDcube reuses and extends existing open software components/modules and initiatives. This includes the FAIR Data Point[16] and Vantage6[22]. Further elements of the FDCube are the Investigation-Study-Assay (ISA) metadata framework[27, 20] for capturing general study metadata, sample (including basic sample characteristics), and assay metadata, and the Phenopackets[21] standards for capturing phenotypic description of a patient/sample. The concept of the FDCube is illustrated in Fig 1 and detailed below from the perspective of a dataset owner and a researcher as a user of that dataset, respectively. The complete and detailed documentation on the FDCube can also be found at https://github.com/Xomics/FAIRDataCube/wiki.
The concept of the FDCube. Dataset owners and researchers as a dataset consumer can both benefit from FDCube on various aspects, including creation of multi-omics data and metadata, querying multi-omics studies and analyzing access-protected data via federated analysis.
Dataset owner
A dataset owner registers their dataset by publishing the metadata on a FAIR Data Point (FDP). The FDP is a metadata repository that provides public access to metadata in accordance with the FAIR principles[16]. The FDP helps dataset owners to publish the metadata of their dataset, and facilitates researchers (dataset users) to find and access information (metadata) about the registered datasets, including pointers to that data (irrespective of data access restrictions and licenses, which is typically arranged at the location of the data store/source).
Considering the various metadata formats adopted by the different X-omics communities, it is reasonable to adopt a standard metadata format as a template for submitting the metadata. To this purpose, we employed the Investigation-Study-Assay (ISA) metadata framework[27, 20] as our basic framework to capture and standardize study (design) information from the different -omics metadata schemes. The ISA metadata schema is commonly adopted by the research community for submission of metabolomics data, for example by EMBL-EBI’s MetaboLights repository[8].
In biomedical studies, clinical characteristics and phenotype information of the study subjects may be collected in addition to (-omics or other) measurements data. This information is essential for making interpretations from research experimental data. Thus, phenotype data need to be standardized as well, so that researchers and clinicians can more easily link phenotypes to experimental data. To achieve this, the Phenopackets framework[21] developed by Global Alliance for Genomics and Health (GA4GH) was adopted. This framework comprises a comprehensive data structure (model), using common ontology terms, to categorise and connect different types of phenotype data.
Researcher
The researcher can be both a data set owner and a data set consumer. As a dataset consumer, the research can search a FDP, which is part of a FDCube, for any dataset of interest. Since all metadata is represented in a linked data format, the researcher can conduct semantic searches on datasets and their corresponding study information by using the SPARQL Protocol and RDF Query Language (SPARQL) query interface. The information that can be queried is the ontologized description of, for instance: samples and their (biological) source; sample preparation; methods and techniques applied; (-omics) measurement and (data) analysis strategies, workflows and reports, including the detected (molecular) data features; research group affiliations. Example questions that may be asked are:
Find all studies which use mass spectrometry-based metabolomics and study a specific metabolic disorder;
Find datasets with more than two -omics types and more than a 100 individuals;
Find measurements for proteins and metabolites that belong to a particular metabolic pathway.
To explore more complex research questions, the researcher could raise a computational request to the dataset owner. This is achieved by the Vantage6 component of the FDCube.
Demonstration of FDCube in TWOC
We adopted the Trusted World of Corona (TWOC) project to demonstrate how to utilize the FDCube for integrated multi-omics federated analysis. The TWOC project aims to contribute to a more sustainable, innovative high-quality and person-oriented healthcare system. To this end, they created a platform in which humans and machines can meet based on FAIR data, protocols and algorithms.
In Fig 2, we provide an example of the creation and application of the FDCube based on a public dataset on COVID-19 featuring multi-omics patient data by Su et al., 2020[29], which was FAIRified as part of the TWOC project.
Use of FDCube in the TWOC demonstrator. The FDCube workflows covers various functions including creating and publishing metadata, browsing and querying the metadata on FDP, and creating and running federated data analysis.
Below is an overview of the workflows for creating, filling, and using the FDCube.
Storage of raw and processed omics data
A multi-modal dataset from COVID-19 patients[29] was prepared, harmonized and FAIRified as part of the TWOC project. The dataset consists of paired omics data layers describing transcriptomics, proteomics and metabolomics of blood samples, and includes comprehensive phenotype information. The dataset is publicly accessible at TWOC’s GitHub repository[12].
To allow interactive and joint querying of data and metadata, we store the processed -omics data along with their feature annotation files. These are both stored in a flat-text tabular .csv format, with features as rows and samples as columns.
Creation of metadata
In the TWOC project, both the ISA metadata schema and Phenopackets schema are adopted. The ISA metadata schema is used as a standard metadata schema to capture metadata about (-omics) experiments, and serializes in an ISA-json file using ISA tools[26, 20]. The ISA tools API provides additional functionalities to convert the ISA objects into linked data.
The FAIR Data Point adopts the W3C’s Data Catalog Vocabulary (DCAT)[15] as its basic metadata schema to capture generic information of the registered datasets and their distributions. To host the experimental metadata in the ISA schema, we defined extra resources and Shapes Constraint Language (SHACL) shapes for the investigation and study. Fig 3 shows all the FDP resources, including the additional investigation and study files. A detailed SHACL shape of the new resources can be found on the FDCube GitHub repository[10].
Overview of all FDP resources. The Investigation and Study resources are defined in addition to the FDP’s default resources.
Example scripts[13] are provided to assist researchers in using these frameworks to capture study and experimental (meta)data as well as phenotype information and to share it on a FDP server.
The Investigation and Study part of ISA is made DCAT-compatible and is used to create an input form to publish metadata on the FDP. Given the flexibility in the Assay part, a potential solution is to first import the ISA metadata and the accompanying phenotype metadata into the triplestore behind the FDP, like GraphDB or Blazegraph. After that, a subset of the metadata can be selected and publically displayed for browsing. The selection of a triplestore is an option that can be selected by the FDCube user. A containerized environment to utilize the ISA-API[7], coupled with the ISA cookbook [6], was created for researchers to FAIRify experimental metadata that is used as input for the FDCube. Moreover, we developed a containerized workflow for the automatic submission of the TWOC clinical metadata to FDP[3].
For phenotype data, a python script[13] was developed based on the phenopackets data schema to automatically convert unFAIRified phenotype information into csv format. We then wrote a YARRRML[19] template that embedded the phenopackets RDF-schema [9], making use of the transformation service in FAIR-in-a-box[5]. This converts the csv file into linked data.
Querying of metadata
The FAIR Data Point can display complete/partial metadata in a human-readable portal for browsing, searching and querying. The FAIRified metadata of the TWOC dataset was published on a FDP portal [4] as shown in Fig 4. A SPARQL query can be run against the metadata via the a query portal to gain deeper knowledge of a dataset, as illustrated in Fig 5. The FDP portal provides a user interface where users can design SPARQL queries. After finding an interesting dataset via browsing or by SPARQL, the researcher could further run follow-up analyses on the target dataset by raising a computation request to the Vantage6 server and retrieve the returning results from the data station via Vantage6.
FAIRified metadata of TWOC dataset published on FAIR Data Point portal
The SPARQL query portal.
Running a data analysis script
Vantage6 delivers the user’s computational request to a data station. A computation request consists of:
A reference to a Docker image, which contains the code (computation) that the researcher would like to run on the target dataset;
A list describing the dataset of interest and its purpose-of-use.
The Vantage6 server handles authentication, keeps track of all computation requests, assigns them to nodes for computation, and stores the returning results of the analyses. The Vantage6 server could also host a private Docker registry.
A Vantage6 node is typically installed at a dataset station. For security reason, the dataset station could stay in an access-protected environment, for example, in a Digital Research Environment (DRE)[2], which is a cloud based, globally available research environment.
Fig 6 shows the Vantage6 user interface at which a researcher can create a task.
Create a task in the Vantage6 user interface
In this example, we used an averaging algorithm hosted on Docker Hub1. This algorithm expects an argument column name to be defined, and will compute the average over that column. We specified in the kwargs fields the parameter ‘column name’ with value ‘age’. The averaging algorithm is dispatched to run on a Vantage6 node, where the dataset is stored. In this example, the dataset is a .csv file prepared from TWOC, which contains a column titled age. The Database field in Fig 6 is labeled default, which is configurable in the Vantage6 node configuration file. For simplicity, this task is created for a collaboration with only one organization (in our example Radboudumc).
Fig 7 shows the result of running the averaging algorithm on the patients’ age in the TWOC dataset, which specifically calculates the average value in the column labelled age. This result can be passed back as the response to the computation request.
Vantage 6 task running result
Conclusion
We have created the FAIR Data Cube, a software and programmatic infrastructure to make -omics data FAIR, and to facilitate the management, reuse, integration and analysis of biomedical (-omics) data, while ensuring data sovereignty, by utilizing Vantage6’ capability of ‘bringing research questions to data’ rather than ‘sending data to research questions’. Vantage6’s management capability covers comprehensive aspects (including organization, collaboration, users, roles, nodes and tasks), and makes FDCube a useful platform to carry out cross-organization federated analysis on decentralized datasets.
We used the FDCube in the TWOC project to demonstrate its capability and usage on, creating and publishing ISA and phenotype meta data, browsing and querying the metadata on FDP, and creating and running federated data analysis on a real dataset.
There are several ways to improve and extend the design and implementation of the current FDCube. For example, a Beacon[25] component can be integrated into FDCube. The reason for this integration is that a FDP (by design) only exposes metadata of datasets. In contrast, Beacon allows for more insights about the presence/absence of specific genomic mutation in a set of data[25]. The combined information from both metadata (via FDP) and real data (via Beacon query), would help a researcher to get more insights into possibly available datasets before designing a data analysis request as dictated by the researcher’s study questions.
Another potential work would be to integrate DataSHIELD and Vantage6 to grant users of Vantage6 access to rich analysis algorithms in DataSHIELD.
Data Availability
All data produced in the present study are available upon reasonable request to the authors
Competing interests
No competing interest is declared.
Author contributions statement
P.A.C.H., A.J.G, M.A.S conceived the project. J.H. worked on phenotype data modelling. A.N., C.V worked on ISA metadata. T.E. managed connection to the TWOC project and FAIRification of the presented dataset. P.K worked on lipidomics metadata. C.D. promoted FDCube and provided scientific feedback. M.B supported the hosting environment. K.J.V provided insights from MOLGENIS perspective. A.N. presented the high level concept diagram. X.L. implemented and set up the architecture with help from all team members. X.L. wrote the manuscript with critical input and revisions from A.N., C.D., C.V., J.H., T.E., P.A.C.H, P.K., K.J.V., A.J.G. All authors reviewed the manuscript.
Acknowledgments
This work was funded by a Dutch Research Council (NWO) grant to The Netherlands X-omics Initiative (project 184.034.019), a Horizon2020 grant to the European Joint Programme on Rare Diseases (grant agreement Number 82557), a Horizon2020 grant to the EATRIS-Plus project (grant agreement Number 871096), and a LSH HealthHolland grant to the Trusted World of Corona (TWOC) consortium.
Footnotes
↵1 harbor2.vantage6.ai/demo/average