Employing a Systematic Approach to Biobanking and Analyzing Clinical and Genetic Data for Advancing COVID-19 Research

Study design

The purpose of the GEN-COVID Multicenter Study is to make the best use of the widest possible sets of patient data and genetic material in order to identify potential links between patient genetic variation and clinical variability, patient presentation and disease severity. By exposing the potential links between genetic variability and disease variability, the study believes it can contribute to improved patient-level diagnostics, prognosis, and personalized treatment of COVID-19. To achieve this overall aim, the following specific objectives are being pursued: i) to perform sequencing (WES) on 2,000 COVID-19 patient samples [performed by the University of Siena (UNISI)]; ii) to perform genotyping (GWAS) on 2000 COVID-19 patients [performed by the Institute for Molecular Medicine of Finland (FIMM)]; iii) to associate the host genetic data obtained on 2,000 COVID-19 patients with severity and prognosis; iv) to share phenotypic data and samples across the GEN-COVID consortium platform as well as in cooperation with research institutions and national platforms through the GEN-COVID Disease Registry and Biobank; v) to share genetic data through the Network of Italian Genome (NIG: http://www.nig.cineca.it/, NIG database: http://nigdb.cineca.it) at CINECA, the largest Italian computing center.

Planned key deliverables of the project are i) to develop a state-of-the-art Patient Registry and Biobank for COVID-19 clinical research with access for academic and industry partners; ii) to understand the genetic and molecular basis of susceptibility to SARS-CoV-2 infection and (susceptibility to a potentially more severe clinical outcome [prognosis] within 12 months); and iii) to understand the genetic profile of patients, contributing to the rapid identification of medicines to be repurposed for personalized therapeutic approaches that demonstrate greater efficacy against the COVID-19 virus. As the initial starting point of this process, the ACE2 gene has already been extensively investigated in the Italian population [5].

The GEN-COVID Multicenter Study includes a network of 22 Italian hospitals, 13 of which from Northern Italy, 5 from Central Italy, and 4 from Southern Italy. It also includes local healthcare units and departments of preventative medicine (https://sites.google.com/dbm.unisi.it/gen-covid). The network continues to grow as more hospitals and healthcare centers express an interest in contributing samples and data. It started its activity on March 16, 2020, following approval by the Ethical Review Board of the Promoter Center, University of Siena (Protocol n. 16929, approval dated March 16, 2020). Written informed consent was obtained from all individuals who contributed samples and data. Detailed clinical and laboratory characteristics (data), specifically related to COVID-19, were collected for all subjects.

Study participants and recruitment

In order to ensure a collection that could be, as much as possible, comprehensive and representative of the Italian population, hospitals from across Italy, local healthcare units, and departments of preventative medicine have been involved in collecting samples and associated patient-level data for the GEN-COVID Multicenter Study. The inclusion criteria for the study are PCR-positive SARS-CoV-2 infection, age ≥ 18 years, and appropriately given informed consent. In addition to the samples collection, an extensive questionnaire is used to assess disease severity and collect basic demographic information from each patient (Supplementary Table 1).

As of July 16th, 2020, we have collected samples and data from 1,033 individuals (1021 without family ties and 12 with family relations). All were positively diagnosed with SARS-CoV-2 and they represent a wide range of disease severity, ranging from hospitalized patients with severe COVID-19 disease to asymptomatic individuals. Infection status was confirmed by SARS-CoV-2 viral RNA polymerase-chain-reaction (PCR) test collected mainly from nasopharyngeal swabs. Recruitment remains ongoing with the goal of including samples and data from 2000 individuals by the end of September 2020. So far recruiting an averages of 200 patients per week.

Data collection and storage

The GEN-COVID registry was designed in order to guarantee data accuracy and, at the same time, to ensure ease of data entry in order to facilitate compliance and save clinicians time. The highest data integrity and data privacy standards, with reference to the EU General Data Protection Regulation (GDPR) [6], were also built into the training for personnel. Samples and data were collected and systematized in order to meet the FAIR Data Principles requirements.

The socio-demographic information included sex, age, and ethnicity. Information about family history, (pre-existing) chronic conditions, and SARS-CoV-2 related symptoms were collected through a detailed core clinical questionnaire as previously reported [7]. This clinical data was continually updated accordingly as new information appeared regarding COVID-19 (Supplementary Table 1). More than 150 clinical items have been collected and synthesized in a binary mode for each involved organ/system: heart, liver, pancreas, kidney, and olfactory/gustatory and lymphoid systems. The collection and organizing methodologies allowed for rapid statistical analysis. Data were handled and stored in accordance with the EU GDPR [6].

Peripheral blood samples in ethylenediamine tetraacetic acid (EDTA)-containing tubes were collected for all subjects. Genomic DNA was centrally isolated from peripheral blood samples using the MagCore®Genomic DNA Whole Blood Kit (Diatech Pharmacogenetics, Jesi, Italy) according to the manufacturer’s protocol at the Promoter Center. For all subjects, aliquots of plasma and serum are also available. Whenever possible, leukocytes were isolated from whole blood by density gradient centrifugation and stored in dimethyl sulfoxide (DMSO) solution and frozen using liquid nitrogen. For the majority of cohort, swab specimens are also available and stored at the reference hospitals.

Genetic data from GWAS and WES were generated for all patients. The generation of such a massive amount of sequencing data required sufficient computing resources able to store and analyse large quantities of data. For this purpose, GEN-COVID took advantage of University of Siena’s participation in the Network for Italian Genomes (NIG, http://www.nig.cineca.it/, NIG database, http://nigdb.cineca.it/), which collects genome sequencing data from the Italian population. NIG has a specific agreement with CINECA, the largest computing centre in Italy and one of the largest in Europe, for the use of the CINECA facility for the storage and analysis of data. Data upload followed quality and regulatory requirements already in place to ensure adequate uniformity and homogeneity levels. Data were formatted to meet the requirements of the FAIR Data Principles and thus made interoperable with other FAIR omics data and reference databases.

Collected laboratory and instrumental data

A continuous quantitative respiratory score, the PaO₂/FiO₂ [Partial pressure of oxygen/Fraction of inspired oxygen ratio (P/F)] was assigned to each patient as an indicator of the respiratory involvement. Taking the normal value >300 as the threshold, we defined four grades of severity score for the PaO₂/FiO₂ ratio: P/F less than or equal to 100, between 101 and 200, between 201 and 300, and greater than 300. A P/F value is not available for the non-hospitalized subjects because the test is only performed in hospitalized patients when needed. Heart involvement was considered on the basis of one or more of the following abnormal data: a cardiac Troponin T (cTnT) value higher than the reference range (<15 ng/L) (indicative of ischemic disorder), an increase in the N-terminal (NT)-pro hormone BNP (NT-proBNP) value (reference value <88 pg/ml for males and <153 pg/ml for females) (indicative of heart failure), and the presence of arrhythmias (indicative of electric disorder). Hepatic involvement was defined on the basis of a clear liver enzymes elevation as alanine transaminase (ALT) and aspartate transaminase (AST) higher than the gender specific reference value (for ALT <41 UI/L in males and <31 UI/L in females; for AST <37 UI/L in males and <31 UI/L in females). Pancreatic involvement was considered on the basis of pancreatic enzymes as pancreatic amylase (PA) and lipase (PL) higher or lower than their specific reference range (13-53 UI/l for PA and 13-60UI/l per PL). Kidney involvement was defined in the presence of a creatinine value higher than the gender specific reference value (0,7-1,20 mg/dl in males and 0,5-1,10 mg/dl in females). Lymphoid system involvement was designated as Natural killer (NK) cells and/or peripheral CD4⁺ T cells below reference value (NK cells>90 cell/ul (mm ^{^}3); CD4+T cells>400 cell/ul (mm^{^}3)). For each patient a numerical grading for the olfactory and gustatory dysfunction was defined through a clinical questionnaire administered by ENT specialists. D-Dimer values of >10X, with or without low Fibrinogen level, were used to interpret the involvement of the blood clotting system. Interleukin 6 (IL6), lactate dehydrogenase (LDH), and c-reactive protein (CRP) values above the reference range (<0,5 mg/dl for CRP and 135-225 UI/l in males and 135-214 UI/l in females for LDH) were used to determine proinflammatory cytokines system involvement.

Whole Exome sequencing

Whole Exome Sequencing with at least 97% coverage at 20x was performed using the Illumina NovaSeq6000 System (Illumina, San Diego, CA, USA). Library preparation was performed using the Illumina Exome Panel (Illumina) according to the manufacturer’s protocol. Library enrichment was tested by qPCR and the size distribution and concentration were determined using Agilent Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA, USA). The Novaseq6000 System (Illumina) was used for DNA sequencing through 150 bp paired-end reads.

Genotyping

Genotyping data on 700,000 genetic markers were obtained on genomic DNA using the Illumina Global Screening Array (Illumina) according to the manufacturer’s protocol. Homo sapiens (human) Genome Reference Consortium Human Build 38 (GRCh38) was used. Quality checks (SNP calling quality, cluster separation, and Mendelian and replication error) were done using GenomeStudio analysis software (Illumina). The computer package Plink v1.90 [8] was used to process 700k SNP-genotyping data and to calculate SNP genotype statistics.

Statistical analysis

Descriptive statistics were calculated to determine the distribution of clinical features by sex, age, and ethnicity. Chi-square tests were used to evaluate the statistical association between the clinical severity of the disease (from no hospitalization to intubation) and the categorical clinical variables: gender, ethnicity, blood group, respiratory severity, taste/smell involvement, heart involvement, liver involvement, pancreas involvement, kidney involvement, lymphoid involvement, cytokines trigger, D-dimer, and number of comorbidities. A linear regression model was used to test the statistical association between COVID-19 severity and age.

The variability within clinical features and their relative relationships have been summarised and described by principal component analysis (PCA). Only numerical variables with a missing rate lower than 50% were selected; these included hyposmia, neutrophils, CRP, fibrinogen, LDH, D-dimer, and number of comorbidities. Missing data were imputed using KNN (k-nearest neighbour) imputation [9], based on Gower distances [10]. After imputation, variables were centered and scaled prior to PCA. Descriptive statistics, chi-square tests, linear regression, and PCA were performed with the R environment for statistical computing [11].

A descriptive analysis of the phenotypes by using a hierarchically-clustered heatmap was performed. In particular, both patients and phenotypes are clusterized with the agglomerative hierarchical clustering methodology, where the chosen metric is the hamming distance and the linkage criterion is the “average” one (unweighted pair group method with arithmetic mean, UPGMA). The corresponding dendrograms of the clusterization are reported in the upper and in the left part of the heat plot. Then the information for the grading of severity of the patients is added a posteriori on the left strip. The resulting plot is obtained with the Python Seaborn package.

The GEN-COVID Biobank (GCB)

The GEN-COVID Biobank (GCB), a collection of bio-specimens from patients affected by COVID-19 and the associated GEN-COVID Patient Registry (GCPR) were established and maintained at the University of Siena using the infrastructure of an already well-established biobank (est. 1998) (http://www.biobank.unisi.it/ScegliArchivio.asp).

The Biobank is closely linked to national and international biobanking efforts aimed at collecting high quality samples and patient data in a uniform manner and ensuring their FAIR (Findable, Accessible, Interoperable and Reusable) management. It is part of the BBMRI-IT [12], EuroBioBank (EBB; [13]), Telethon Network of Genetic Biobanks (TNGB; [14]), and RD-Connect [15]. The biobank and registry are ISO-certified (certificate 199556-2016-AQ-ITA-ACCREDIA) and accredited according to SIGU (the Italian Society of Human Genetics) requirements (Certificate 204107-2016-AQ-ITA-DNV).

Collected biological samples include peripheral blood, plasma, serum, primary leukocytes, and DNA samples. Samples were stored in a dedicated biobank section while associated clinical data were entered in the related registry. The biobank and registry were organized according to the highest scientific standards, preserving patients’ and citizens’ privacy, while providing services to the healthcare and scientific community to develop better treatments, test diagnostic tools, and advance COVID-19 and coronavirus research. Biobank personnel are responsible for sample pseudonymization, storage, and insertion in the online biobank catalogue.

Geographical coverage

The GEN-COVID Multicenter Study reached a large number of subjects throughout Italy. Tuscany, which is the region in which the study is carried out, contributes presently 22.8% of enrolled patients. The Northern Italian regions, particularly Lombardy and Venetia, currently contribute 52.3% of enrolled patients (Figure 2). This distribution reflects closely the incidence of SARS-CoV-2 infection per 100,000 inhabitants for each Italian region, as updated to 4 July 2020 [3].

The GEN-COVID Patient Registry (GCPR)

From April 7, 2020 to July 16, 2020, the GEN-COVID Patient Registry (GCPR) collected clinical data from a total of 1033 Italian SARS-Cov-2 PCR-positive individuals. For each individual, we collected clinical information using standardized clinical schedules (Supplementary Table 1). The study protocol also provides access to patients’ medical records and continual clinical data updating in order to secure continuity for patient follow-up.

The mean age of the entire cohort is presently 58.7 years (range 18-99). The cohort is presently predominantly male (57.1%) with a mean age of 59.5 years (range 18-99); the mean age of the females is 57.6 years (range 19-98) (Table 1). About 40.3% of the cohort has no chronic conditions. The overall case-fatality rate (CFR) is 3.6% (37 deaths among 1,033 cases with a mean age of 75.2 years [range 62-91]. Regarding the ethnicity, the cohort is composed of 998 White (96.61%), 21 Hispanic (2.03%), 4 Black (0.38%), and 10 Asian (0.96%) patients (Table 1).

Subjects have been divided into five qualitative severity clinical categories depending on the need for hospitalization, the respiratory impairment and, consequently, the type of ventilation required: i) hospitalized and intubated (9.5%); ii) hospitalized and CPAP-BiPAP and high-flows oxygen treated (18.4%); iii) hospitalized and treated with conventional oxygen support only (31.55%); iv) hospitalized without respiratory support (14.8%); v) not hospitalized pauci/asymptomatic individuals (25.75%) (Group 4 to 0 in Table 1).

Gender distribution was statistically significantly different among the 5 groups (p- value=7.81×10’⁶). In the group with high care intensity (Group 4), 72.4% of subjects were male, while in the group with the milder phenotype (Group 0) 59.8% of subjects were female (Table 2). Hyposmia and/or hypogeusia were present in 13.9% of cases in Group 4, 25.3% in Group 3, 31.6% in Group 2, 19.3% in Group 1, and in 57.1% of Group 0. A slight statistically significant difference among the 5 groups was found regarding the presence of comorbidities (p-value=0.012). No statistically significant difference was present for ethnicity and blood group distribution (Table 2).

Figure 3 shows the relationships between continually updated laboratory variables from PCA. The first two principal components explain 42.4% of the variability in the data (PC1: 23.8%; PC2: 18.6%). Neutrophils, LDH, and D-dimer appear to be positively correlated, while fibrinogen and CRP, and hyposmia and the number of comorbidities have been found, pairwise, to be negatively correlated. The largest contributors to PC1 were LDH (24.4%), neutrophils (24%), D-dimer (23.8%), and hyposmia (15.3%); the largest contributors to PC2 were hyposmia (25.8%), fibrinogen (22.1%), CRP (19.8%), and the number of comorbidities (15.4%) (Figure 3).

The continually updated laboratory values used in Figure 3 can be further mined through clinical reasoning and represented as a binary clinical classification for organ/system damage (Table 3).

Table 4 shows the prevalence of different organ/systems damage in the 5 different clinical categories based on respiratory failure (Table 4). Heart involvement was detected in 55% of subjects in Group 4, 39% of subjects in Group 3, 34.1% in Group 2, and 21.6% in Group 1. Liver involvement was present in 72.4% of cases in Group 4, 59.3% in Group 3, 46% in Group 2, and 33.7% in Group 1. Statistically significant difference among the 5 groups was found for all organs/systems, except for the lymphoid system.

Finally, Figure 4 shows by dendrogram COVID-19 phenotype can be clustered using the above reported clinical data representations. Hierarchical Clustering analysis identified five main clinical categories and several subcategories: A) severe multisystemic, with either thromboembolic (A1) or pancreatic variant (A2); B) cytokine storm, either moderate (B1) or severe with liver involvement (B2); C) mild, either with (C1) or without hyposmia (C2); D) moderate, either without (D1) or with (D2) liver damage; E) heart type, either with (E1) or without (E2) liver damage (Figure 4).

GEN-COVID Genetic Data Repository (GCGDR)

WES and Genotype (GWAS) data were generated within the GEN-COVID Genetic Data Repository (GCGDR). In order to be able to store and analyse the massive amount of genomic data (mainly WES with coverage > 97% at 20x, and prospetically including also WGS) generated with the analysis of the entire cohort of samples populating the biobank, we relied on the NIG. External users can upload and analyse data using the NIG pipeline by registering and creating a specific project. A section dedicated to COVID-19 samples has been created within the NIG database (http://nigdb.cineca.it/) that provides variant frequencies as a free tool for both clinicians and researchers.

The data from WES are available both in Variant Call Format (VCF) file or as binarized file, according to different classes of variants: i) rare variants (minor allele frequency (MAF)<1%); ii) low frequency variants (MAF<5%); iii) common polymorphisms (MAF>5%) in either homozygosity or supposed compound heterozygosity, with rare or low frequency variants. The distribution of these 3 classes of variants according to mutated genes in our cohort is shown in Supplementary Figure 1.

From WES 580,688 variants have been called: of these, 543,138 are SNP and 37,550 are MNP (multi-nucleotide polymorphisms). Exonic SNPs were distributed over the 22 autosomes of the human genome, plus the sex chromosomes. The average missing rate was 0.01, with per-sample maximum value of 0.017. 15,285 SNP loci had a missing rate greater than 5%. The average MAF was 0.032 (std. dev. 0.091), with a right-skewed distribution (median MAF = 0.0007). Only 1,041 SNPs were monomorphic (0.2%), but 437,246 (80.5%) had a frequency < 0.01. From the genotype perspective, the average observed heterozygosity was 0.047.

The data from high-density (700k) SNP genotyping are also generated on the same cohort and shared with international collaborations, including the COVID-19 Host Genetics Initiative (https://covid-19genehostinitiative.net/) and with GoFAIR VODAN [16]. From this analysis, SNP genotypes at 730,059 loci distributed over the entire human genome have been obtained. The average missing-rate was 0.015, with per-sample maximum value of 0.042. 11,163 SNP loci had a missing rate greater than 5%. The average MAF was 0.113 (std. dev. 0.145), with a right-skewed distribution (median MAF = 0.035). In total, 147,579 SNPs were monomorphic (20.2%). From the genotype perspective, the average observed heterozygosity was 0.155.

Clinical data representation and interpretation

Clinical data may be represented and consequently interpreted in different ways. The simplest way of representation is using the raw data of laboratory/instrument values. In this case, reasoning about which value has to be considered and/or at which time of clinical evolution the value needs to be measured is necessary in order to have consistency within the cohort. PCA analysis using the WORSEN score at the time of admission has shown the expected variability with hyposmia to be juxtaposed to the number of comorbidities and thus representing a marker of less severity. The fibrinogen value is juxtaposed to inflammatories markers, such as CRP (and D-Dimer and LDH) because it is consumed during the prothrombotic state. We can conclude that such raw laboratory values are fairly good for representing the clinical variability of the cohort in classical PCA analysis.

A more elaborate way of representing clinical data is to filter the raw laboratory/instrument values by clinical reasoning, which often requires a face-to-face meeting with organ reference specialists and direct access to the patients’ medical records. The proposed mediation of such a clinical methodology for COVID-19 is represented in Table 3 and its distribution against lung dysfunction synthesised in Table 4.

Involvement of relevant organs or systems is represented in binary and is then used for representing COVID-19 as a systemic disorder (Figure 4). We propose this representation as one of the best, being closer to the real complexity of the disease. It should be considered for use in further data mining and correlation with genetic data. The emerging clinical categories from Hierarchical Cluster Analysis point to specific types and subtypes that are more likely to have common genetic factors.

As unmasked by our dendrogram (group A), there is indeed a growing body of evidence suggesting that, in addition to the common respiratory symptoms (fever, cough and dyspnea), COVID-19 severely-ill patients can often have symptoms of a multisystemic disorder [21]. Multiple organ failure due to diffuse microvascular damage is an important cause of death in COVID-19 severely affected patients [22]. In line with our definition of an A1 subgroup, a retrospective study on 21 deaths after SARS-Co-V2 infection recently reported that 71% of the patients who died had disseminated intravascular coagulation (DIC), while the incidence of DIC in surviving patients was 0.6% [23]. These data suggest that DIC is an important risk factor for increased in-hospital mortality and special attention should be paid to its early diagnosis and treatment.

While a debate still exists about the significance of pancreatic enzyme elevations during COVID-19 infection and the capability of SARS-CoV-2 virus to induce pancreatic injury due to cytotoxic effects [24, 25], it is worth noting that among patients with a multisystemic involvement we observe a subclass of individuals (group A2) with pancreatic damage, likely suggesting a secondary effect of SARS-CoV-2 infection on a subgroup of genetically predisposed individuals. Inflammatory cytokine “storm,” has been reported as playing a key role in the severe immune injury to the lungs caused by T-cell over-activation (group B) [26]. While some investigators have suggested a potential mechanism of myocardial injury due to COVID-19-induced cytokine storm that is mediated by a mixed T helper cell response in combination to hypoxia [27], our findings indicate rather a distinct class of patients (group E) presenting with heart involvement in the absence of an inflammatory cascade. This would tend to support the hypothesis that SARS-CoV-2 may directly damage myocardial tissue and induce a major cardiovascular event. Thus, as currently recommended, our research reinforces the need to monitor plasma cTnT and NT-proBNP levels in COVID-19 patients. In line with current evidence [28; 29], although liver injury seems to occur more frequently among critically ill patients with COVID-19 (group B), it can also be present in non-critically ill patients (groups D and E) and, as suggested, it could be mostly related to prolonged hospitalization and viral shedding duration. This allows defining, for each group, a clinical subclass according to this organ involvement.

A recent extensive review determined the prevalence of chemosensory deficits based on pooling together forty-two studies reporting on 23,353 patients [30]. Estimated random prevalence was 38.5% for olfactory dysfunction, 30.4% for taste dysfunction, and 50.2% for overall chemosensory dysfunction. No correlation with age was detected, but anosmia/hypogeusia decreased with disease severity and ethnicity turned out to play a significant role with Caucasians having a 3 to 6 times higher prevalence of chemosensory deficits than East Asians. In accordance with evidence found in the literature, hyposmia was mostly represented among patients in group C with mild clinical symptoms [31].

Genetic data representation and interpretation

Similar to the clinical data, large aggregates of genetic data derived from WES may be represented, and consequently interpreted, in different ways. After variant calling, it is possible to use data as such, or variants can be prioritized and filtered according to standard bioinformatics procedures [32], such as damaging effect predictions, healthy population allele frequency, and gene constraints to variation.

Alternatively, it is also possible to represent data in a binary mode as follows: i) select missense, splicing, and loss of function variants below 1% (rare variants); ii) select missense, splicing, and loss of function variants between 1% and 5% (low frequency variants); iii) select missense, splicing, and loss of function variants above 5% (common polymorphisms) in either homozygosity or supposed compound heterozygosity with rare or low frequency variants. The majority of patients showed about 3% of mutated genes in the above class i), 5% in class ii) and 28% in class iii) variants (Supplementary Figure 1A). No patients showed variants in more than 8,000 genes (Supplementary Figure 1B).

Protein interaction network and pathway analysis have been widely used to uncover and describe genetic relationships in complex diseases, such as cancer [33;34]. For example, over-representation analysis of the biological processes and pathways significantly affected by mutations will be instrumental to empower the statistical detection of genetic signatures associated to specific COVID-19 phenotypes and to reduce the number of parameters to consider (e.g. dimensionality reduction) with the purpose of developing robust algorithms for the prediction of genetic susceptibility to COVID-19 infection and response. Variants, genes, or biological processes will be employed as features to train interpretable, supervised machine learning classifiers (e.g. gradient boosting decision trees [35;36]), which will ease the identification of the genetic factors associated with clinical phenotypes.

While data collection is being consolidated and brought to completion according to the study design, we have started to work on a relatively new methodology based on Topological Data Analysis to provide a detailed multidimensional and multiscale exploration of the whole exome data that can drive an AI selection of genes that provide higher predictive power in a machine learning model. The method will be presented, together with the results, in a forthcoming paper.

Post-Mendelian model of complex diseases

Previous attempts to interpret the genetic bases of complex disorders have failed with very few exceptions, even in those disorders in which (like COVID-19) twin studies demonstrated a very high rate of heritability, such as in psychiatric disorders. The reason for this story containing such a lack of scientific success resides in several weak points in the overall genetic approach to complex diseases: i) the method used to represent the complexity of the phenotype; ii) the procedure employed to represent the huge amount of different genetic data; and iii) the absence of a robust mathematical model able to interpret genetic data in non-Mendelian (non-rare) disorders. This paper provides a contribution to the first 2 points, likely paving the way for a solution to the third.

Frequently the phenotype of common (complex) disorders is oversimplified, thus attenuating reliable correlation with genetic data. Limiting the representation to differences of single parameters, such as respiratory assistance (intubation, CPAP-BiPAP, oxygen supplementation, etc.), is a possible trap for studies on complex disorders, as may be the case with COVID-19. Similarly, genetic data are often too large to be mined and fragmented in different non-communicating methods, betting on either the power of common polymorphisms (GWAS) or the power of variant accumulation (burden gene test for WES). The binary representation we are proposing here, together with network propagation for feature reduction, and followed by machine learning approaches, may help in this task. A rare disorder called TAR (OMIM # 274000) is teaching us that combinatorial rules of rare variant(s) with more common polymorphism(s) is what we are looking for [37].

The GEN-COVID Multicenter Study with its Registry (GCPR), Biobank (GCB), and Genetic Data Repository (GCGDR) is structured to continually link with leading European and international research organizations, public and private, as well as with regulatory and public health authorities for developing COVID-19 and SARS-related medicines research and treatment protocols. The success of the developing research and understanding of COVID-19 and the underlying SARS-CoV-2 virus will rely in large part on human biological materials and patient-level data that is comprehensively collected and systematically organized with careful attention to sample and data integrity as well as the FAIR Data Principles. Improving diagnostics, developing existing or new therapeutics, improving treatment protocols, and even developing public health policies relies upon a foundation of evidence that requires the comprehensive, patient, and systematic collection and organizing of COVID-19 patient biological samples and data of high integrity, confidence, and interoperability. The GEN-COVID Multicenter Study’s GCPR, GCB, and GCGDR present a model that can be further explored as a systematic approach to sample and data collection while also being immediately deployable in our collective fight against COVID-19.