Abstract
Functional connectivity magnetic resonance imaging (fcMRI) is a well-established technique for studying brain networks in both healthy and diseased individuals. However, no fcMRI-based biomarker has yet achieved clinical relevance. To establish better understanding of the state of the art in quantifying abnormal connectivity in comparison to a reference distribution, for potential use in individual patients, we have conducted a scoping review over 5672 entries from the last 10 years. We have located five publications proposing methods of abnormal connectivity quantification, reported these methods and formalized them. We also illustrated the emerging trends and technical innovations in fcMRI research that may facilitate development of individualized fcMRI-based abnormal connectivity metrics.
1 Introduction
Functional connectivity magnetic resonance imaging (fcMRI), first used for connectivity analysis in humans by Biswal et al. [1] and based on the blood oxygen level dependent (BOLD) signal [2, 3, 4], is widely regarded as a valuable imaging method for the inquiry into connectivity in human [5, 6] and non-human [7] brain research alike. With the scientific community increasingly reconceptualizing neurodegenerative [8, 9], psychiatric [10] and neuro-oncological [11, 12, 13] disorders as “network disorders”, fcMRI-based biomarkers that quantify abnormal connectivity in relation to the distribution in a healthy reference sample may pave a way for a connectivity metrics suited for validation and application in clinical diagnostics.
To date, no fcMRI biomarker has achieved clinical relevance. This can be linked to two major challenges: (1) limited interpretability of the acquired signal in consequence of intra-subject variability and device- and procedure-related confounds [14, 15, 16] and (2) a lack of well-established and readily accessible reference values for functional connectivity in individuals despite available datasets (e.g. Human Connectome Project [52] and 1000connectomes [17]). Alleviating these issues through systematic use of reference samples and normative modeling may permit consistent data interpretation and pave the way for an fcMRI biomarker accessible enough for potential incorporation into diagnostic practice.
In light of the potential benefits of establishing such a normative model for fcMRI, and considering the successful biomarker normalization attempts in other brain imaging modalities [18, 19, 20], two assertions can be made.
Firstly, there exists an apparent unmet medical need for validated and clinically implemented fcMRI-based abnormality metrics that satisfy the criteria of relationality and countability. Herein, a relational metric may be defined as a metric that relies on a control cohort sufficiently representative of the target individual, allowing to establish a normative model of connectivity that compares a given individual to a distribution of controls and quantifies the discrepancy, while a countable metric may be defined as an interval or rational metric that can be used as grounds for grading or comparison.
Secondly, there is minimal study coverage pertaining to the introduction and validation of such metrics, which limits current insight into individualized abnormality detection in functional connectivity.
An initial step toward addressing the question of normative modelling in fcMRI consists in a scoping review of fcMRI-based metrics of connectivity abnormality, the results of which we present here. Within the scope of this paper, we review and analyze the fcMRI abnormality metrics yielded by our search, explore the degree of their refinement, and determine their readiness for clinical validation. Moreover, we discuss the need of moving beyond group comparison and towards quantitative fcMRI anomaly metrics for application in individual patients. We also elucidate emerging trends and technical innovations in fcMRI research that may facilitate development of relational fcMRI-based abnormality metrics.
2. Methods
2.1 Overall Protocol
We have conducted our review in adherence to the general framework of scoping reviews proposed by Arksey and O’Malley [21] and refined by Levac et al. [22]. We reported our results in compliance with the Preferred Reporting Items for Systematic reviews and Meta-Analyses (PRISMA) extension for scoping reviews (PRISMA-ScR) [23]. The PRISMA-ScR compliance checklist can be accessed in the Supplementary Materials.
2.2 Review Objectives
Within the scope of this review, we intended to determine (1) whether there exist metrics to quantify the deviation of functional connectivity in an individual patient from a reference population, (2) whether they are validated to guarantee sufficient technology readiness and clinical utility and (3) whether they satisfy the criteria of relationality and countability outlined in the introduction.
In pursuit of this objective, we have reviewed the state-of-the-art (SOTA) in fcMRI connectivity abnormality detection, analyzed the results, formalized them, and reported our findings.
2.3 Information Sources, Search Strategy, Data Acquisition and Handling
We have leveraged the Google Scholar database for our seearch. We set the query year range at the years 2014-2024 and employed Publish or Perish 8.10.4612.8838 [24] to automate our query. We searched in 1-year batches to yield the most entries and circumvent the internal limit of 1000 entries per query. We input the following search request: “fcMRI connectivity connectome abnormality detection anomaly map deviation individual reference metric.”
All data was aggregated using pandas 2.1.1 [25] and NumPy 1.23.5 [26], exported as comma-separated values, and uploaded for subsequent group review on a secure team space in Notion [27]. Using Notion’s integrated tools and functions, we removed damaged or empty entries. The remaining entries were subjected to screening and eligibility assessment (see below).
2.4 Study Screening and Selection
We employed a 2-phase screening and eligibility selection strategy. During the screening phase, we excluded sources that (1) did not report research based on fMRI or did not use BOLD signal, (2) reported experiments on participants under 18 years of age, (3) did not have a healthy reference cohort against which the patients would be gauged, (4) were reviews, (5) were preprints, (6) were book chapters, (7) did not report research on resting-state fcMRI, (8) were not accessible for full text, (9) reported research on data acquired with a field strength under 3.0 T, (10) were theses or dissertations, (11) were meta-analyses, (12) reported research conducted on non-human data, (13) were citation records, (14) were abstract almanacs or miscellaneous publications, (15) were conference papers, (16) were study protocols or (17) were not in English.
Eligibility assessment phase consisted in elimination of articles that did not||||||||||||||| report metrics that satisfy the criteria of relationality and countability outlined in the introduction. Eligibility assessment relied on an in-depth inspection of the “Methods” section and a deeper examination of other paper sections in cases where it was necessary. Edge cases were resolved by reviewer consensus.
2.5 Study Analysis
The sources which passed screening and selection were fully studied. Subsequently, we extracted the metric computation methods reported by the respective authors, described them, and formalized them. To explore the degree of their refinement, state of validation, and level of applicability in a clinical setting, we chose to follow the citations of the articles in question (for better narration consistency and text legibility, these searches will be reported within the results section). Subsequently, we integrated these findings to yield our statements. We additionally assigned to every metric a Technology Readiness Level (TRL) as specified by ISO 16290:2013 [28] in the edition of EU Commission Decision C(2017)7124 [29], elucidated for fcMRI-based abnormality detection applications as per Table 1.
3 Results
3.1 Query Results
Our query cumulatively returned 5696 entries, 5672 of them valid (non-empty, not damaged or fragmentary) entries. After screening, 4964 sources were excluded (Fig. 1), while 708 sources were deemed eligible for selection. Only 5 passed selection and were subjected to a full-depth analysis. A PRISMA flow diagram is available in Fig. 2.
3.2 State of the Art and its Aspects
3.2.1 The Nenning Index
Nenning et al. [30] introduced a voxel-level connectivity abnormality metric in their 2020 glioblastoma paper. Briefly, it is computed as follows: (1) voxel-wise connectivity matrices for both patients and controls (80 control subjects) are built using z-scored Pearson correlations; (2) element-wise average of control population connectivity matrices is computed to yield a group average “normal” connectivity matrix; (3) a vector of voxel-wise differences is computed between the patients and group average as row-wise cosine similarity; (4) for every voxel in controls’ connectivity matrices and the group average matrix, cosine similarities are computed to yield voxel-wise distribution; from that distribution, the median and mean absolute deviation (MAD) are computed (the “voxel mean” and “voxel MAD” respectively); (5) for every patient and for every patient voxel’s cosine similarity, an abnormality score is computed as the difference of cosine similarity and voxel mean, subsequently divided by the voxel MAD.
Analytically, this can be summarized as follows: where is the connectivity vector of voxel v for patient is the connectivity vector of voxel v for control subject ci, with i = 1, 2, …, N and N = 80 being the number of control subjects, is the average connectivity vector of voxel v across all control subjects, ∥· ∥ denotes the Euclidean norm of a vector, · represents the dot product between two vectors, median(·) computes the median of a set of values and MAD(·) computes the median absolute deviation of a set of values.
It is important to mention that Nenning’s team focused on reporting abnormality in non-infiltrated regions but pointed out that the inclusion of tumor infiltrated regions did not significantly alter the overall connectivity signature. Additionally, they demonstrate that in glioblastoma, functional proximity to the tumor tends to be reflected stronger than structural proximity in coefficients derived from fcMRI signal, while visual, somatomotor, and limbic networks tend to exhibit anomaly coefficients more evenly informed by both spatial and functional distance alike. Finally, Nenning’s team demonstrate precedence of network anomalies before tumor recurrence, highlighting a potential prognostic capacity for abnormality index computation.
PubMed citation check revealed no further studies employing this index in their computations; however, the longitudinal character of the study in focus supports the assignment to this index of a TRL 5 out of 9.
3.2.2 The Dysconnectivity Index
Stoecklein and Liu [31] present another voxel-level connectivity abnormality metric in their publication on gliomas. It is computed as follows: (1) voxel-wise connectivity matrices are built for both patients and controls (1000 control subjects) using Pearson correlations; (2) for every control subject connectivity matrix, every voxel position in the matrix, and every element in the voxel, a distribution of connectivity coefficients is built; (3) the distribution’s mean and standard deviation are computed to yield respective elements of the mean and standard deviation vectors; (4) for every patient connectivity matrix, every voxel position in the matrix, and every element in the voxel, a z-score is computed for using the elements of the mean and standard deviation vectors computed before (i.e., for i-th element in the patient’s voxel, respective i-th element of the mean and standard deviation vector is used) to yield a vector of z-scores; (5) a sum of z-scores higher than a specific threshold is computed to yield the voxel-level “abnormality coefficient.”
Analytically, for the voxel at the position i this can be summarized as follows: where P ij is the connectivity coefficient at voxel position i, j for the patient, is the connectivity coefficient at voxel position i, j for control subject c, N is the number of control subjects, T is the specific threshold, and I(·) is the indicator function, which evaluates to 1 if the condition is true and 0 otherwise. The authors have conducted computations for the entire brain (without tumor mask exclusion) and demonstrated not only that tumor sites can be captured by their index, but that abnormality can be detected far beyond the lesion itself, even in the contralateral hemisphere, particularly in high grade gliomas. They have also shown that, in glioma, their abnormality index correlates with neurocognitive performance, WHO grade, PET metabolic data, and IDH mutation status. Additionally, the authors hypothesized that abnormal connectivity may not only originate from tumor functional or structural proximity but also indicate sub-clinical tumor cell infiltration and speculated that functional disruption also indicates possible tumor cell infiltration.
PubMed citation check revealed two studies based on this index. In the first publication [32], the authors demonstrated that their abnormality index (in more recent sources referred to as DCI - the “dysconnectivity index”) can be employed to assess immune effector cell-associated neurotoxicity syndrome (ICANS) in patients under CAR-T therapy and hypothesized that it may be used to objectify damage to functional networks in encephalopathies; furthermore, the authors stated that their index may provide an imaging correlate to trace and possibly predict neurotoxic side-effects of oncologic treatment. In the second publication [33], the authors show a direct association between the DCI and the perifocal edema volume in meningiomas, as well as neurocognitive performance (i.e., higher DCI implies larger edema and more degraded cognition). The sizable body of knowledge amassed in relation to this index, as well as validation for different diseases of the human brain and their sequelae, allows us to assign to this index a TRL of 6 out of 9.
3.2.3 The Doucet Normative Person-Based Similarity Index
In their publication, Doucet et al. [34] report the normative person-based similarity index (nPBSI). Computed from both functional connectivity and cortical morphometry per aspect, their index explicitly seeks to make a patient’s condition relative to a set control population (93 control subjects). Doucet’s group presents four indices for which clinical, genetic, demographic, and environmental correlates have been described - normative cortical thickness PBSI (nPBSI-CT), normative subcortical volume PBSI (nPBSI-SV), normative module cohesion PBSI (nPBSI-MC) and normative module integrations (nPBSI-MI).
Within the scope of this review, our attention was focused on the fcMRI-based module cohesion and module integration metrics, computed as follows: (1) the patient’s brain is parcellated into default mode, central executive, salience, sensorimotor, and visual networks; (2) within-module connectivity is represented as the average value of a voxel wise z-transformed Pearson correlation coefficient between all of the module’s voxel pairs and used to build a patient’s module cohesion profile, encoded as a module cohesion feature vector; (3) between-module connectivity is represented as z-transformed Pearson correlation coefficients of the modules’ averaged time series and used to build a patient’s module integrations profile, encoded as a module integrations feature vector, and finally, (4) the nPBSI-MC or nPBSI-MI are computed as averaged Spearman correlations between the patient and the healthy controls’ respective (module cohesion or module integrations) feature vectors.
Analytically, for the patient p this can be summarized as follows: where N represents the number of brain modules (default mode, central executive, salience, sensorimotor, and visual networks), Mi is the set of voxels in module i, Ki is the number of voxel pairs in module i, is the Pearson correlation coefficient between voxels vp|and vq for the patient p, is the Pearson correlation coefficient between voxels vp and vq for a healthy control h, is the Pearson correlation coefficient between the average time series of modules i and j for the patient p, is the same for a healthy control h, is the Fisher z-transformation, ρ denotes the Spearman correlation coefficient, H is the set of healthy controls and |H| is the number of healthy controls.
PubMed citation check revealed no studies employing the normative index from this publication in their computations of functional connectivity metrics. The closest possible match [35] relied on computing both the within- and between-network connectivity but did not compute the nPBSI itself. Modest validation for bipolar disorder and lack of nPBSI validation for other disorders justifies the assignment to this metric of a TRL 4 out of 9.
3.2.4 The Network Topography Spatial Similarity Index
Silvestri and Corbetta present a spatial similarity index (SSI) for network to-pographies derived from independent component analysis (ICA) in their 2022 publication on gliomas [36]. Briefly, it is computed as follows: (1) rs-fcMRI data of the control population (308 individuals) are subjected to a group ICA (G-ICA) to yield group-level template independent component (IC) maps for ten functional networks (specifically, visual, sensorimotor, auditory, cingulo-opercular, dorsal attention, fronto-parietal, default mode, cognitive control, frontal and language networks); (2) the group-level template IC maps are used as spatial constraints for group information-guided ICA (GIG-ICA) of both controls and patients (24 individuals) to produce individual-specific, single-subject level IC maps; (3) for each IC in subject, a cosine similarity is computed between a single-subject IC map and a template IC map thresholded at a value of 1 and is yielded as the network topography spatial similarity index.
Analytically, this can be expressed as follows: where SSIIC is the spatial similarity index for a given independent component, represents the rs-fMRI data of the control population, are the group-level template IC maps for the ten functional networks obtained from group ICA, Ds is the rs-fMRI data of subject s, GIG-ICA(Ds; produces the single-subject IC map Ss,j for subject s and component j using the group-level templates as spatial constraints, Threshold1(Tj) denotes the template IC map Tj thresholded at a value of 1, the numerator (·) represents the dot product between the two vectors and the denominator (∥ · ∥) represents the Euclidean norm (magnitude) of the vectors.
The team around Silvestri and Corbetta reported testing structural and functional proximity of their index to the tumor sites, describing partial overlap of index abnormalities and glioma-infiltrated areas and highlighting index abnormalities in non-infiltrated areas. They also analyzed changes in network topography scores and neuropsychological performance and were able to capture a statistically relevant relationship between the SSI and the attention domain. PubMed citation check revealed no studies employing this normative index in their computations of functional connectivity metrics. Modest validation for gliomas and lack of validation for other disorders justifies the assignment to this metric of a TRL 4 out of 9.
3.2.5 The Morgan Network Topology Method
Morgan et al. present various metrics and indices in their publication on the role of anterior hippocampus in mesial temporal lobe epilepsy (mTLE) [37]. Their computations rely on multi-modal data and operate within four topologies: the streamline length (TLEN), structural connectivity (TSC), functional connectivity (TF C) and resting-state network topology (TRSN). Within the scope of our review, we will focus on the functional connectivity topology and its respective distance index, as no similar index has been reported for the resting-state network topology.
Briefly, it is computed as follows: (1) functional connectivity maps are built for controls (70 individuals) and patients (40 individuals, of them 29 with right mTLE and 11 with left mTLE) from z-transformed functional connectivity matrices through age regression and subsequent averaging of signal over 109 anatomical ROIs; (2) a topology is built from the functional connectivity maps by selecting 55 ROIs of a single hemisphere for patients and controls; (3) a seed vector is used to slice anterior hippocampal connectivity from the topology into a connectivity vector for both patients and controls; (4) the connectivity vector is stratified along connectivity intensity into “bins” to yield their respective connectivity vectors of k elements for both patients and controls; (5) for patient and bin, the Mahalanobis distance between the patient’s bin connectivity vector and the mean of controls’ bin connectivity vectors is computed and yielded as connectivity deviation metric.
Analytically, this can be summarized as follows: with patient’s connectivity vector in bin, controls’ mean vector in bin and controls’ covariance matrix in bin as, respectively, and where Mi is a functional connectivity matrix (size 109 × 109) for individual i, S(Mi) denotes a selection operator extracting a 55 × 55 hemisphere-specific submatrix from Mi, R is a seed vector (size 1 × 55) with 1 at the anterior hippocampus position and 0 elsewhere, Bb(·) symbolizes the binning function that selects elements belonging to bin b based on connectivity intensity, Nc = 70 is the number of control individuals, µb is the mean vector of controls’ connectivity vectors in bin b and Sb is the covariance matrix of controls’ connectivity vectors in bin b.
PubMed citation check revealed two studies which reported intriguing use of the logic behind this computational approach. The first publication of interest by Morgan et al. [38] reports use of similar connectivity profiling techniques and the Mahalanobis distance for outcome prediction in mTLE patients by means of distance computation between a patient’s connectivity profile and a normative population of individuals who achieved seizure-free status after mesial temporal lobe surgery. Notably, the team around Morgan reported sensitivity of 100% and specificity of 90% for their prediction approach.
The second publication by Guerrero-Gonzalez et al. [39] does not pertain to functional MRI, but describes use of the comparable logic of normative modeling and Mahalanobis distance computing to quantify abnormality in tractography of traumatic brain injury patients.
The epilepsy-specific focus of Morgan’s distance-based approach limits the scope of potential use of this metric; however, success of similar computational approaches in other modalities and remarkable performance of the Mahalanobis distance-based index in the surgical outcome prediction task support the assignment to this metric of a TRL 5 out of 9.
4 Discussion
4.1 Group Comparison Currently Prevails in Studies of Abnormal Connectivity
In this scoping review, we have been able to show that, despite the strong knowledge base to support the concept of neurodegenerative [8, 9], psychiatric [10] and neuro-oncological [11, 12, 13] as “network disorders”, a metric capable of evaluating and quantifying large-scale functional brain network disruptions in individual patients is yet to be developed, validated and made accessible enough for potential incorporation into diagnostic practice.
We also demonstrated that, despite the significant benefits of relational metrics as integral elements of normative modeling [40], we could only retrieve five such metrics of functional connectivity deviation that have been proposed within the last ten years. Of note, in many studies that we evaluated for this review, the findings and the hypotheses that lead to these findings were built around the aspiration to illustrate binary differences between patients and healthy controls, which resulted in reports of metrics being increased or decreased in patients without a clearly specified relation between the increment of metric and increment of pathological state. The development of patient-centric fcMRI markers requires moving beyond group comparison and toward relational metrics based on normative populations that span variability in demographic and procedural factors.
4.2 Artificial Intelligence and Big Data Emerge as Methods in fcMRI Research
The advent of big data and artificial intelligence-based methods in fcMRI research may boost the development of relational connectivity metrics by enhancing the current computational approaches and data accessibility.
The drastic progress in computing technology [41] has made possible the widespread use of industrial-grade hardware acceleration of previously strictly linear computing through parallel computing with the help of much more readily accessible graphical processing units (GPUs) [42, 43]. Improved hardware-software synergy now permits optimization of both speed and efficiency of data engineering and machine learning, allowing for faster simultaneous read-/write operations and deeper insight into highly complex multidimensional data. This is well-manifested by the packages for accelerated Python computing (e.g. CuPy[44] or Dask [45]), optimized tensor storage solutions (e.g. Zarr [46] or Xarray [47]), new Neuroimaging Informatics Technology Initiative (NifTI) image manipulation modules (e.g. Xibabel [48]) or the advancements in the field of machine learning (ML) frameworks [49, 50, 51].
Simultaneously, high-quality data can be accessed freely by virtue of recognized cohorts (e.g. Human Connectome Project, Alzheimer’s Disease Neuroimaging Initiative or Brain Genomics Superstruct Project [52, 53, 54]) and open-access data repositories (e.g. OpenNeuro [55]), which permits compilation of harmonized, statistically powerful reference datasets, capturing variability across demographics and technical parameters. The utility of accounting for these factors is well-substantiated by evidence of variables such as age [56, 57, 58], sex [59, 60] and scan parameters [61, 62] having significant influence on fcMRI metrics. Therefore, creation of large-scale reference datasets augmented by technical and demographic parameters may help pave the way for normative modelling in fcMRI.
Moreover, the current rise of deep learning models for operations on fcMRI data can help streamline previously time-consuming elements of data preprocessing and enrichment, potentially accelerating research on relational fcMRI-based metrics manyfold. This is prominently exemplified by ML breakthroughs in the area of structural image preprocessing with algorithms such as FastSurfer [63], a deep learning pipeline for brain segmentation, cortical surface reconstruction, cortical label mapping and thickness analyses. Similar advancements have also been reported for affine registration with tools such as SynthMorph [64], a model that resolves a tensor-to-tensor mapping problem for an image pair, yielding a compatible spatial transform. Lastly, experimental ML-boosted integrated pipelines for fcMRI image preprocessing (e.g. DeepPrep [65]) have also been proposed.
In summary, the current circumstances create a uniquely favorable setting for more practical progress on relational fcMRI-based metrics of abnormal con-||| nectivity.
4.3 Limitations
Our search only comprises sources released before mid-May 2024. Additionally, our search terms might not include all relevant publications. In particular, preprints, theses and dissertations have been excluded as reports that have not undergone a peer review process. Additionally, not all publications could be accessed for full text. Furthermore, due to considerably less generalizable dynamics of neurobiological development in pediatric and adolescent individuals, a decision was made not to consider publications that concerned persons under 18 years of age. Finally, if a publication matched more than one exclusion criterion during screening, its exclusion was attributed to a single most prominently matching criterion in an effort to prevent redundant statistical entries.
5 Summary
Patients suffering from neuro-oncological, psychiatric and neurodegenerative disorders can benefit from individualized detection and quantification of abnormal functional connectivity. However, no fcMRI-derived biomarkers have yet seen widespread adoption in clinical research or practice. Within the scope of this scoping review, we have asserted both the necessity and the current absence of a well-established relational and countable metric for abnormal functional connectivity in individuals. We have subsequently leveraged the Google Scholar database to retrieve sources that matched our search criteria and subjected them to PRISMA-compliant screening and selection to yield items for subsequent in-depth analysis. We have yielded and demonstrated five currently reported methods/metrics for relational, normative quantification of abnormal connectivity and formalized their computation methods. Building upon our results, we have discussed the need of moving beyond group comparison and toward quantitative fcMRI anomaly metrics for application in individual patients and briefly elucidated the emerging trends and technical innovations in fcMRI research that may facilitate development of relational metrics of functional connectivity.
Funding Information
S.S. received support through the LMU Investment Fund (LMU Excellence AOST: 865105-7). Funding sources had no role in the design, implementation, analysis, interpretation, or reporting of this research.
Conflict of Interest
The authors have no relevant conflict of interest to declare.
Data Availability Statement
Data sharing is not applicable in the context of this publication, as no datasets were generated or analyzed during this scoping review. The tabular reports of the included and excluded articles are available from the corresponding authors upon reasonable request.
Code Availability Statement
No novel code was generated during the current study. Minimal scripting was done to support data aggregation.
Inclusion and Ethics Statement
This scoping review concerns peer-reviewed publications and therefore does not require ethical approval.
Author Contributions
A. T. - Conceptualization, Methodology Selection & Implementation, Data Collection, Entry Screening, Source Eligibility Selection, Source Analysis, Formalization & Integration of Findings, Original Draft Preparation, Visualization, Review and Editing, Project Administration.
D. V. - Data Collection, Entry Screening, Source Eligibility Selection, Source Analysis, Formalization & Integration of Findings, Original Draft Preparation, Visualization, Review and Editing, Project Administration.
H. K. - Data Collection, Entry Screening, Source Eligibility Selection, Source Analysis, Formalization & Integration of Findings, Original Draft Preparation, Visualization, Review and Editing, Project Administration.
R. L. - Entry Screening, Source Eligibility Selection, Source Analysis, Formalization & Integration of Findings.
P. M. - Entry Screening, Source Eligibility Selection.
D. R. - Entry Screening, Source Eligibility Selection.
A. D. - Entry Screening.
A. V. - Entry Screening.
V. P. - Entry Screening.
S. W. - Source Eligibility Selection, Source Analysis.
S. S. - Conceptualization, Source Analysis, Formalization & Integration of Findings, Review and Editing, Supervision, Funding Acquisition, Project Administration, Resources, Oversight and Approvals.
D.V. and H.K. contributed equally to this publication.
D.R. and P.M. contributed equally to this publication.
A.D. and A.V. contributed equally to this publication.
Acknowledgements
We would like to thank Julia Ruat for her invaluable support in the management of this project.
References
- [1].↵
- [2].↵
- [3].↵
- [4].↵
- [5].↵
- [6].↵
- [7].↵
- [8].↵
- [9].↵
- [10].↵
- [11].↵
- [12].↵
- [13].↵
- [14].↵
- [15].↵
- [16].↵
- [17].↵
- [18].↵
- [19].↵
- [20].↵
- [21].↵
- [22].↵
- [23].↵
- [24].↵
- [25].↵
- [26].↵
- [27].↵
- [28].↵
- [29].↵
- [30].↵
- [31].↵
- [32].↵
- [33].↵
- [34].↵
- [35].↵
- [36].↵
- [37].↵
- [38].↵
- [39].↵
- [40].↵
- [41].↵
- [42].↵
- [43].↵
- [44].↵
- [45].↵
- [46].↵
- [47].↵
- [48].↵
- [49].↵
- [50].↵
- [51].↵
- [52].↵
- [53].↵
- [54].↵
- [55].↵
- [56].↵
- [57].↵
- [58].↵
- [59].↵
- [60].↵
- [61].↵
- [62].↵
- [63].↵
- [64].↵
- [65].↵