Analysis Pipeline Research Articles

Computational analysis and NGS pipeline for Human GRK2 in complex with G-beta-gamma in cardiovascular disease with Biopython

Computational analysis and NGS pipeline for Human GRK2 in complex with G-beta-gamma in cardiovascular disease with Biopython

An analytical approach for identifying trend-seasonal components and detecting unexpected behaviour in psychological time-series.

The recent advances in technological capabilities have led to a massive production of time-series data and remarkable progress in longitudinal designs and analyses within psychological research. However, implementing time-series analysis can be challenging due to the various characteristics and complexities involved, as well as the need for statistical expertise. This paper introduces a statistical pipeline on time-series analysis for studying the changes in a single process over time at either a population or individual level, both retrospectively and prospectively. This is achieved through systemization and extension of existing modelling and inference techniques. This analytical approach enables practitioners not only to track but also to model and evaluate emerging trends and apparent seasonality. It also allows for the detection of unexpected events, quantifying their deviations from baseline and forecasting future values. Given that other discernible population- and individual-level changes in psychological and behavioural processes have not yet emerged, continued surveillance is warranted. A near real-time monitoring tool of time-series data could guide community psychological responses across multiple ecological levels, making it a valuable resource for field practitioners and psychologists. An empirical study is conducted to illustrate the implementation of the introduced analytical pipeline in practice and to demonstrate its capabilities.

Design and implementation of a metagenomic analytical pipeline for respiratory pathogen detection

ObjectiveWe developed an in-house bioinformatics pipeline to improve the detection of respiratory pathogens in metagenomic sequencing data. This pipeline addresses the need for short-time analysis, high accuracy, scalability, and reproducibility in a high-performance computing environment.ResultsWe evaluated our pipeline using ninety synthetic metagenomes designed to simulate nasopharyngeal swab samples. The pipeline successfully identified 177 out of 204 respiratory pathogens present in the compositions, with an average processing time of approximately 4 min per sample (processing 1 million paired-end reads of 150 base pairs). For the estimation of all the 470 taxa included in the compositions, the pipeline demonstrated high accuracy, identifying 420 and achieving a correlation of 0.9 between their actual and predicted relative abundances. Among the identified taxa, 27 were significantly underestimated or overestimated, including only three clinically relevant pathogens. We also validated the pipeline by applying it to a clinical dataset from a study on metagenomic pathogen characterization in patients with acute respiratory infections and successfully identified all pathogens responsible for the diagnosed infections. These findings underscore the pipeline’s effectiveness in pathogen detection and highlight its potential utility in respiratory pathogen surveillance.

Multiplexed chromatin immunoprecipitation sequencing for quantitative study of histone modifications and chromatin factors.

ChIP-seq is a widely used technique for studying histone post-translational modifications and DNA-binding proteins. DNA fragments associated with a specific protein or histone modification epitope are captured by using antibodies, sequenced and mapped to a reference genome. Albeit versatile and popular, performing many parallel ChIP-seq experiments to compare different conditions, replicates and epitopes is laborious, is prone to experimental variation and does not allow quantitative comparisons unless adequate spike-in chromatin is included. We present a detailed protocol for performing and analyzing a multiplexed quantitative chromatin immunoprecipitation-sequencing experiment (MINUTE-ChIP), in which multiple samples are profiled against multiple epitopes in a single workflow. Multiplexing not only dramatically increases the throughput of ChIP-seq experiments (e.g., profiling 12 samples against multiple histone modifications or DNA-binding proteins in a single experiment), but also enables accurate quantitative comparisons. The protocol consists of four parts: sample preparation (i.e., lysis, chromatin fragmentation and barcoding of native or formaldehyde-fixed material), pooling and splitting of the barcoded chromatin into parallel immunoprecipitation reactions, preparation of next-generation sequencing libraries from input and immunoprecipitated DNA and data analysis using our dedicated analysis pipeline. This pipeline autonomously generates quantitatively scaled ChIP-seq tracks for downstream analysis and visualization, alongside necessary quality control indicators. The entire workflow requires basic knowledge in molecular biology and bioinformatics and can be completed in 1 week. MINUTE-ChIP empowers biologists to perform every ChIP-seq experiment with an appropriate number of replicates and control conditions, delivering more statistically robust, exquisitely quantitative and biologically meaningful results.

B-146 Unlocking the potential of large-cohort proteomics studies with the Orbitrap Astral mass platform

Abstract Background Large-cohort proteomics analysis using mass spectrometry is a powerful approach to discover and validate new biomarkers. In combination with clinical data and computational analysis, large-cohort proteomics study brings opportunities in improving early diagnosis, refining patient stratification, and predicting/monitoring treatment response. Yet, to achieve meaningful biological insights in large-cohort studies, robust, reproducible, and comprehensive proteome profiling in a high-throughput manner still remains challenging. Here, we use a high-resolution accurate mass (HRAM) Oribtrap Astral platform to enable high-quality and robust protein quantification across thousands of LC-MS/MS analyses. With a throughput of 100 samples/day, we can reproducibly profile ∼9000 proteins from human cell line and ∼800 proteins from undepleted plasma across multiple instruments and more than 10 consecutive days in a 24/7 operation mode. Methods Experiments were performed on multiple Orbitrap Astral mass spectrometers with and without FAIMS Pro interface. Chromatographic separations were performed using a trap/elute method on Vanquish Neo system at a 11-minute gradient and a flow rate of 1ul/min, resulting in a throughput of 100 samples/day. Undepleted plasma digest was analyzed in a 24/7 mode for > 1000 injections on each LC-MS/MS setup. In addition, HeLa digest as quality control (QC) was analyzed every 12hours. The Oribtrap Astral platform was operated in data-independent acquisition (DIA) mode using a full scan with m/z 380-980, and DIA MS/MS scans with an isolation window of 2Th. Resulting DIA raw files were automatically transferred and processed using Chimerys in a beta version of Proteome Discoverer software. Results Multiple LC-MS/MS systems were operated in DIA mode either with or without FAIMS Pro device in a 24/7 operation mode at a throughput of 100SPD. Undepleted plasma digest was analyzed with >1000 injections on each LC-MS/MS setup. To monitor the baseline performance, HeLa digest served as QC were analyzed periodically every 12 hours with 3 technical replicates each time. To effectively manage and analyze these thousands of data files generated, we developed an automated data transfer and analysis pipeline. Automatically, the resulting DIA raw data files were immediately transferred to a server, then processed by state-of-the-art intelligent search algorithm Chimerys. Benefiting from the ultra-high scan speed of up to 200Hz on the Orbitrap Astral mass spectrometer, a much narrower isolation window width of 2Th was applied in the DIA method, comparing to 10-20Th on the classic DIA method. As a result, ∼9000 proteins from HeLa digest and ∼800 proteins from undepleted plasma digest were identified within only a 11-minutes gradient, respectively. More than 80% of the plasma proteins were reproducibly identified and quantified from all the runs on each LC-MS/MS setup, indicating a great reproducibility from run-to-run longitudinally. Stable and robust peptide quantitation was observed by extracting peptides with high, medium, and low abundant across the runs. Importantly, QC showed no performance degradation throughout the entire study, indicating high robustness of the entire LC-MS/MS setup. Conclusions Orbitrap Astral mass platform can comprehensively analyze the proteome of >1000s of sample robustly and reproducibly in a high-throughput manner, addressing the needs in large-cohort studies.

Quantification of Total Ketone Bodies in Biological Matrices Using Stable Isotope-Labeled Internal Standards and RP-UHPLC-MS/MS.

Acetoacetate (AcAc) and D-beta-hydroxybutyrate (D-βOHB), the two major ketone bodies found in circulation, are linked to multiple physiological and pathophysiological states. Therefore, analytical methodologies surrounding the quantification of total ketone body (TKB) concentrations in biological matrices are paramount. Traditional methods to quantify TKBs relied on indirect spectrophotometric assays with narrow dynamic ranges, which have been significantly improved upon by modern mass spectrometry (MS)-based approaches. However, the lack of stable isotope-labeled internal standards (ISs) for AcAc and the need to distinguish D-βOHB from its closely related structural and enantiomeric isomers pose significant obstacles. Here, we provide a protocol to synthesize and quantify a [13C] stable isotope-labeled IS for AcAc, which, in conjunction with a commercially available [2H] stable isotope-labeled IS for βOHB, allows TKBs to be measured across multiple biological matrices. This rapid (7min) analysis employs reverse phase ultra-high performance liquid chromatography (RP-UHPLC) coupled to tandem MS (MS/MS) to distinguish βOHB from three structural isomers using parallel reaction monitoring (PRM), providing excellent specificity and selectivity. Finally, a method is provided that distinguishes D-βOHB from L-βOHB using a simple one-step derivatization to produce the corresponding diastereomers, which can be chromatographically resolved using the same rapid RP-UHPLC separation with new PRM transitions. In summary, this method provides a rigorous analytical pipeline for the analysis of TKBs in biological matrices via leveraging two authentic stable isotope-labeled ISs and RP-UHPLC-MS/MS.

Distinguishing Artifactual Fatty Acid Dimers from Fatty Acid Esters of Hydroxy Fatty Acids in Untargeted LC-MS Pipelines.

Untargeted metabolomics is a powerful profiling tool for the discovery of possible biomarkers of disease onset and progression. Analytical pipelines applying liquid chromatography (LC) and mass spectrometry (MS)-based methods are widely used to survey a broad range of metabolites within various metabolic pathways, including organic acids, amino acids, nucleosides, and lipids. Accurate and complete identification of putative metabolites is an ongoing challenge in untargeted metabolomics studies. Highly sensitive instrumentation can result in the detection of adduct and fragment ions that form reproducibly and contain identifiable ions that are difficult to distinguish from metabolic pathway intermediates, which may result in false-positive identification. At concentrations as low as 10μM, free fatty acids have been found to form homo- and heterodimers in untargeted metabolomics pipelines that resemble the lipid class fatty acid esters of hydroxy fatty acids (FAHFAs), resulting in misidentification. This chapter details a protocol for LC-MS-based untargeted metabolomics using hydrophilic interaction chromatography (HILIC) that specifically aids in distinguishing artifactual fatty acid dimers from endogenous FAHFAs.

Sensorimotor brain–computer interface performance depends on signal-to-noise ratio but not connectivity of the mu rhythm in a multiverse analysis of longitudinal data

Objective.Serving as a channel for communication with locked-in patients or control of prostheses, sensorimotor brain-computer interfaces (BCIs) decode imaginary movements from the recorded activity of the user's brain. However, many individuals remain unable to control the BCI, and the underlying mechanisms are unclear. The user's BCI performance was previously shown to correlate with the resting-state signal-to-noise ratio (SNR) of the mu rhythm and the phase synchronization (PS) of the mu rhythm between sensorimotor areas. Yet, these predictors of performance were primarily evaluated in a single BCI session, while the longitudinal aspect remains rather uninvestigated. In addition, different analysis pipelines were used to estimate PS in source space, potentially hindering the reproducibility of the results.Approach.To systematically address these issues, we performed an extensive validation of the relationship between pre-stimulus SNR, PS, and session-wise BCI performance using a publicly available dataset of 62 human participants performing up to 11 sessions of BCI training. We performed the analysis in sensor space using the surface Laplacian and in source space by combining 24 processing pipelines in a multiverse analysis. This way, we could investigate how robust the observed effects were to the selection of the pipeline.Main results.Our results show that SNR had both between- and within-subject effects on BCI performance for the majority of the pipelines. In contrast, the effect of PS on BCI performance was less robust to the selection of the pipeline and became non-significant after controlling for SNR.Significance.Taken together, our results demonstrate that changes in neuronal connectivity within the sensorimotor system are not critical for learning to control a BCI, and interventions that increase the SNR of the mu rhythm might lead to improvements in the user's BCI performance.

Tracking Small Animals in Complex Landscapes: AComparison of Localisation Workflows for Automated Radio Telemetry Systems.

Automated radio telemetry systems (ARTS) have the potential to revolutionise our understanding of animal movement by providing a near-continuous record of individual locations in the wild. However, localisation errors in ARTS data can be very high, especially in natural landscapes with complex vegetation structure and topography. This curtails the research questions that may be addressed with this technology. We set up an ARTS grid in a valley with heterogeneous vegetation cover in the Colombian high Andes and applied an analytical pipeline to test the effectiveness of localisation methods. We performed calibration trials to simulate animal movement in high- or low-flight, or walking on the ground, and compared workflows with varying decisions related to signal cleaning, selection, smoothing, and interpretation, along with four multilateration approaches. We also quantified the influence of spatial features on the system's accuracy. Results showed large variation in localisation error, ranging between 0.4-43.4 m and 474-1929 m, depending on the localisation method used. We found that the selection of higher radio signal strengths and data smoothing based on the temporal autocorrelation are useful tools to improve accuracy. Moreover, terrain ruggedness, height of movement, vegetation type, and the location of animals inside or outside the grid area influence localisation error. In the case of our study system, thousands of location points were successfully estimated for two high-altitude hummingbird species that previously lacked movement data. Our case study on hummingbirds suggests ARTS grids can be used to estimate small animals' home ranges, associations with vegetation types, and seasonality in occurrence. We present a comparative localisation pipeline, highlighting the variety of possible decisions while processing radio signal data. Overall, this study provides guidance to improve the resolution of location estimates, broadening the application of this tracking technology in the study of the spatial ecology of wild populations.

Small RNAs in plasma extracellular vesicles define biomarkers of premanifest changes in Huntington's disease.

Despite the advances in the understanding of Huntington's disease (HD), there is a need for molecular biomarkers to categorize mutation carriers during the preclinical stage of the disease preceding functional decline. Small RNAs (sRNAs) are a promising source of biomarkers since their expression levels are highly sensitive to pathobiological processes. Here, using an optimized method for plasma extracellular vesicles (EVs) purification and an exhaustive analysis pipeline of sRNA sequencing data, we show that EV-sRNAs are downregulated early in mutation carriers and that this deregulation is associated with premanifest cognitive performance. Seven candidate sRNAs (tRF-Glu-CTC, tRF-Gly-GCC, miR-451a, miR-21-5p, miR-26a-5p, miR-27a-3p and let7a-5p) were validated in additional subjects, showing a significant diagnostic accuracy at premanifest stages. Of these, miR-21-5p was significantly decreased over time in a longitudinal study; and miR-21-5p and miR-26a-5p levels correlated with cognitive changes in the premanifest cohort. In summary, the present results suggest that deregulated plasma EV-sRNAs define an early biosignature in mutation carriers with specific species highlighting the progression and cognitive changes occurring at the premanifest stage.

Highly-Multiplexed Immunofluorescence PhenoCycler Panel for Murine FFPE Yields Insight into Tumor Microenvironment Immunoengineering

Spatial proteomics profiling is an emerging set of technologies that has the potential to elucidate the cell types, interactions, and molecular signatures that make up complex tissue microenvironments, with applications in the study of cancer, immunity, and much more. An emerging technique in the field is Co-Detection-by-indEXing (CODEX), recently renamed as the PhenoCycler system. This is a highly-multiplexed immunofluorescence imaging technology that relies on oligonucleotide-barcoded antibodies and cyclic immunofluorescence to visualize many antibody markers in a single specimen while preserving tissue architecture. Existing PhenoCycler panels are primarily designed for fresh-frozen tissues. Formalin-fixed paraffin-embedded (FFPE) blocks offer several advantages in preclinical research, but few antibody clones have been identified in this setting for PhenoCycler imaging. Here, we present a novel PhenoCycler panel of 28 validated antibodies for murine FFPE tissues. We describe our workflow for selecting and validating clones, barcoding antibodies, designing our panel, and performing multiplex imaging. We further detail our analysis pipeline for comparing marker expressions, clustering and phenotyping single-cell proteomics data, and quantifying spatial relationships. We then apply our panel and analysis protocol to profile the effects of three gene-delivery nanoparticle formulations, in combination with systemic anti-PD1, on the murine melanoma tumor immune microenvironment. Intralesional delivery of genes expressing the costimulatory molecule 4-1BBL and the cytokine IL-12 led to a shift towards intratumoral M1 macrophage polarization and promoted closer associations between intratumoral CD8 T cells and macrophages. Delivery of IFNγ, in addition to 4-1BBL and IL-12, further increased markers of antigen presentation on tumor cells and intratumoral antigen-presenting cells but also promoted greater expression of checkpoint marker PD-L1 and closer associations between intratumoral CD8 T cells and PD-L1-expressing tumor cells. These findings help to explain the benefits of 4-1BBL and IL-12 delivery while offering additional mechanistic insights into the limitations of IFNγ therapeutic efficacy.

Presenting the framework of the whole slide image file Babel fish: An OCR-based file labeling tool

IntroductionMetadata extraction from digitized slides or whole slide image files is a frequent, laborious, and tedious task. In this work, we present a tool to automatically extract all relevant slide information, such as case number, year, slide number, block number, and staining from the macro-images of the scanned slide.We named the tool Babel fish as it helps translate relevant information printed on the slide. It is written to contain certain basic assumptions regarding, for example, the location of certain information. This can be adapted to the respective location. The extracted metadata can then be used to sort digital slides into databases or to link them with associated case IDs from laboratory information systems. Material and methodsThe tool is based on optical character recognition (OCR). For most information, the easyOCR tool is used. For the block number and cases with insufficient results in the first OCR round, a second OCR with pytesseract is applied.Two datasets are used: one for tool development has 342 slides; and another for one for testing has 110 slides. ResultsFor the testing set, the overall accuracy for retrieving all relevant information per slide is 0.982. Of note, the accuracy for most information parts is 1.000, whereas the accuracy for the block number detection is 0.982. ConclusionThe Babel fish tool can be used to rename vast amounts of whole slide image files in an image analysis pipeline. Furthermore, it could be an essential part of DICOM conversion pipelines, as it extracts relevant metadata like case number, year, block ID, and staining.

Possible Carbon Dioxide above the Thick Aerosols of GJ 1214 b

Sub-Neptune planets with radii smaller than Neptune (3.9 R ⊕) are the most common type of planet known to exist in the Milky Way, even though they are absent in the solar system. These planets can potentially have a large diversity of compositions as a result of different mixtures of rocky material, icy material, and gas accreted from a protoplanetary disk. However, the bulk density of a sub-Neptune, informed by its mass and radius alone, cannot uniquely constrain its composition; atmospheric spectroscopy is necessary. GJ 1214 b, which hosts an atmosphere that is potentially the most favorable for spectroscopic detection of any sub-Neptune, is instead enshrouded in aerosols (thus showing no spectroscopic features), hiding its composition from view at previously observed wavelengths in its terminator. Here, we present a JWST NIRSpec transmission spectrum from 2.8 to 5.1 μm that shows signatures of CO2 and CH4, expected at high metallicity. A model containing both these molecules is preferred by 3.3σ and 3.6σ as compared to a featureless spectrum for two different data analysis pipelines, respectively. Given the low signal-to-noise of the features compared to the continuum, however, more observations are needed to confirm the CO2 and CH4 signatures and better constrain other diagnostic features in the near-infrared. Further modeling of the planet’s atmosphere, interior structure and origins will provide valuable insights about how sub-Neptunes like GJ 1214 b form and evolve.

A population-scale analysis of 36 gut microbiome studies reveals universal species signatures for common diseases

The gut microbiome has been implicated in various human diseases, though findings across studies have shown considerable variability. In this study, we reanalyzed 6314 publicly available fecal metagenomes from 36 case-control studies on different diseases to investigate microbial diversity and disease-shared signatures. Using a unified analysis pipeline, we observed reduced microbial diversity in many diseases, while some exhibited increased diversity. Significant alterations in microbial communities were detected across most diseases. A meta-analysis identified 277 disease-associated gut species, including numerous opportunistic pathogens enriched in patients and a depletion of beneficial microbes. A random forest classifier based on these signatures achieved high accuracy in distinguishing diseased individuals from controls (AUC = 0.776) and high-risk patients from controls (AUC = 0.825), and it also performed well in external cohorts. These results offer insights into the gut microbiome’s role in common diseases in the Chinese population and will guide personalized disease management strategies.

Processing pipelines and analytical methods for single-cell DNA methylation sequencing data.

Single-cell DNA methylation sequencing technology has seen rapid advancements in recent years, playing a crucial role in uncovering cellular heterogeneity and the mechanisms of epigenetic regulation. As sequencing technologies have progressed, the quality and quantity of single-cell methylation data have also increased, making standardized preprocessing workflows and appropriate analysis methods essential for ensuring data comparability and result reliability. However, a comprehensive data analysis pipeline to guide researchers in mining existing data has yet to be established. This review systematically summarizes the preprocessing steps and analysis methods for single-cell methylation data, introduces relevant algorithms and tools, and explores the application prospects of single-cell methylation technology in neuroscience, hematopoietic differentiation, and cancer research. The aim is to provide guidance for researchers in data analysis and to promote the development and application of single-cell methylation sequencing technology.

ScRNAbox: empowering single-cell RNA sequencing on high performance computing systems

BackgroundSingle-cell RNA sequencing (scRNAseq) offers powerful insights, but the surge in sample sizes demands more computational power than local workstations can provide. Consequently, high-performance computing (HPC) systems have become imperative. Existing web apps designed to analyze scRNAseq data lack scalability and integration capabilities, while analysis packages demand coding expertise, hindering accessibility.ResultsIn response, we introduce scRNAbox, an innovative scRNAseq analysis pipeline meticulously crafted for HPC systems. This end-to-end solution, executed via the SLURM workload manager, efficiently processes raw data from standard and Hashtag samples. It incorporates quality control filtering, sample integration, clustering, cluster annotation tools, and facilitates cell type-specific differential gene expression analysis between two groups. We demonstrate the application of scRNAbox by analyzing two publicly available datasets.ConclusionScRNAbox is a comprehensive end-to-end pipeline designed to streamline the processing and analysis of scRNAseq data. By responding to the pressing demand for a user-friendly, HPC solution, scRNAbox bridges the gap between the growing computational demands of scRNAseq analysis and the coding expertise required to meet them.

AAV vector-derived elements integrate into Cas9-generated double-strand breaks and disrupt gene transcription

We previously developed an adeno-associated virus (AAV) Cas9 gene therapy for Angelman syndrome that integrated into the genome and prematurely terminated Ube3a-ATS. Here, we assessed the performance of 3 additional AAV vectors containing S.aureus Cas9 invitro and invivo, and 25 vectors containing N. meningitidis Cas9 invitro, all targeting single sites within Ube3a-ATS. We found that none of these single-target gRNA vectors were as effective as multi-target gRNA vectors at reducing Ube3a-ATS expression in neurons. We also developed an anchored multiplex PCR sequencing method and analysis pipeline to quantify the relative frequency of all possible editing events at target sites, including AAV integration and unresolved double-strand breaks. We found that integration of AAV was the most frequent editing event (67%-89% of all edits) at three different single target sites, surpassing insertions and deletions (indels). None of the most frequently observed indels were capable of blocking transcription when incorporated into a Ube3a-ATS minigene reporter, whereas two vector derived elements-the poly(A) and reverse promoter-reduced downstream transcription by up to 50%. Our findings suggest that the probability that a gene trapping AAV integration event occurs is influenced by which vector-derived element(s) are integrated and by the number of target sites.

ThermoNet: advanced deep neural network-based thermogram processing pipeline for automatic time series analysis of specific skin areas in moving legs

Infrared thermography is an emerging technique in biomedical research, potentially providing diagnostic insights into psychological stress, physical strain, muscle fatigue, inflammation, tissue damage, and diseases with thermogenic effects. However, manual analysis strategies are frequently applied causing incomparable, non-reproducible results and hampering standardization. Moreover, widely applied manual analysis cannot recognize blood vessel-related thermal radiation patterns during physical exercise. Therefore, an enhanced processing pipeline, “ThermoNet”, has been developed to automatically process thermograms captured during running. For acquisition, an automatic temperature calibration technique has been introduced to obtain reliable pixel-temperature mapping. The thermograms are semantically segmented in the processing pipeline to extract the anatomical regions of interest (ROIs) by a state-of-the-art deep neural network rather than considering both legs as a single area. A second neural network further examines the ROIs to identify different venous and arterial (perforator) patterns. Within the segments, advanced statistical features are computed to provide time series data. Separate analysis of venous and perforator vessel patterns is carried out on individual connected components, resulting in the extraction of 276 features for each thermogram. The enhanced ROI extraction achieved a high accuracy for the left and right calf on the manually annotated test set. Each step of the ThermoNet pipeline represents a significant improvement over previous analysis methods. Finally, ThermoNet is a transferable pipeline for automatic, reproducible, and objective analysis of ROIs in thermal image sequences of moving test individuals.

The GENESIS database and tools: A decade of discovery in Mendelian genomics

In the past decade, human genetics research saw an acceleration of disease gene discovery and further dissection of the genetic architectures of many disorders. Much of this progress was enabled via data aggregation projects, collaborative data sharing among researchers, and the adoption of sophisticated and standardized bioinformatics analyses pipelines. In 2012, we launched the GENESIS platform, formerly known as GEM.app, with the aims to 1) empower clinical and basic researchers without bioinformatics expertise to analyze and explore genome level data and 2) facilitate the detection of novel pathogenic variation and novel disease genes by leveraging data aggregation and genetic matchmaking. The GENESIS database has grown to over 20,000 datasets from rare disease patients, which were provided by multiple academic research consortia and many individual investigators. Some of the largest global collections of genome-level data are available for Charcot-Marie-Tooth disease, hereditary spastic paraplegia, and cerebellar ataxia. A number of rare disease consortia and networks are archiving their data in this database. Over the past decade, more than 1500 scientists have registered and used this resource and published over 200 papers on gene and variant identifications, which garnered >6000 citations. GENESIS has supported >100 gene discoveries and contributed to approximately half of all gene identifications in the fields of inherited peripheral neuropathies and spastic paraplegia in this time frame. Many diagnostic odysseys of rare disease patients have been resolved. The concept of genomes-to-therapy has borne out for a number of such discoveries that let to rapid clinical trials and expedited natural history studies. This marks GENESIS as one of the most impactful data aggregation initiatives in rare monogenic diseases.

A Robust [18F]-PSMA-1007 Radiomics Ensemble Model for Prostate Cancer Risk Stratification.

The aim of this study is to investigate the role of [18F]-PSMA-1007 PET in differentiating high- and low-risk prostate cancer (PCa) through a robust radiomics ensemble model. This retrospective study included 143 PCa patients who underwent [18F]-PSMA-1007 PET/CT imaging. PCa areas were manually contoured on PET images and 1781 image biomarker standardization initiative (IBSI)-compliant radiomics features were extracted. A 30 times iterated preliminary analysis pipeline, comprising of the least absolute shrinkage and selection operator (LASSO) for feature selection and fivefold cross-validation for model optimization, was adopted to identify the most robust features to dataset variations, select candidate models for ensemble modelling, and optimize hyperparameters. Thirteen subsets of selected features, 11 generated from the preliminary analysis plus two additional subsets, the first based on the combination of robust and fine-tuning features, and the second only on fine-tuning features were used to train the model ensemble. Accuracy, area under curve (AUC), sensitivity, specificity, precision, and f-score values were calculated to provide models' performance. Friedman test, followed by post hoc tests corrected with Dunn-Sidak correction for multiple comparisons, was used to verify if statistically significant differences were found in the different ensemble models over the 30 iterations. The model ensemble trained with the combination of robust and fine-tuning features obtained the highest average accuracy (79.52%), AUC (85.75%), specificity (84.29%), precision (82.85%), and f-score (78.26%). Statistically significant differences (p < 0.05) were found for some performance metrics. These findings support the role of [18F]-PSMA-1007 PET radiomics in improving risk stratification for PCa, by reducing dependence on biopsies.