The European Genome-phenome Archive of human data consented for biomedical research

Inflammatory response signature score model for predicting immunotherapy response and pan-cancer prognosis

Background

Inflammatory responses influence the outcome of immunotherapy and tumorigenesis by modulating host immunity. However, systematic inflammatory response assessment models for predicting cancer immunotherapy (CIT) responses and survival across human cancers remain unexplored. Here, we investigated an inflammatory response score model to predict CIT responses and patient survival in a pan-cancer analysis.

Methods

We retrieved 12 CIT response gene expression datasets from the Gene Expression Omnibus database (GSE78220, GSE19423, GSE100797, GSE126044, GSE35640, GSE67501, GSE115821 and GSE168204), Tumor Immune Dysfunction and Exclusion database (PRJEB23709, PRJEB25780 and phs000452.v2.p1), European Genome-phenome Archive database (EGAD00001005738), and IMvigor210 cohort. The tumor samples from six cancers types: metastatic urothelial cancer, metastatic melanoma, gastric cancer, primary bladder cancer, renal cell carcinoma, and non-small cell lung cancer.We further established a binary classification model to predict CIT responses using the least absolute shrinkage and selection operator (LASSO) computational algorithm.

Findings

The model had high predictive accuracy in both the training and validation cohorts. During sub-group analysis, area under the curve (AUC) values of 0.82, 0.80, 0.71, 0.7, 0.67, and 0.64 were obtained for the non-small cell lung cancer, gastric cancer, metastatic urothelial cancer, primary bladder cancer, metastatic melanoma, and renal cell carcinoma cohorts, respectively. CIT response rates were higher in the high-scoring training cohort subjects (51%) than the low-scoring subjects (27%). The five-year survival rates in the high- and low score groups of the training cohorts were 62% and 21%, respectively, while those of the validation cohorts were 54% and 22%, respectively (P < 0·001 in all cases). Inflammatory response signature score derived from on-treatment tumor specimens are highly predictive of response to CIT in patients with metastatic melanoma. A significant correlation was observed between the inflammatory response scores and tumor purity. Regardless of the tumor purity, patients in the low score group had a significantly poorer prognosis than those in the high score group. Immune cell infiltration analysis indicated that in the high score cohort, tumor-infiltrating lymphocytes were significantly enriched, particularly effector and natural killer cells. Inflammatory response scores were positively correlated with immune checkpoint genes, suggesting that immune checkpoint inhibitors may have benefited patients with high scores. Analysis of signature scores across different cancer types from The Cancer Genome Atlas revealed that the prognostic performance of inflammatory response scores for survival in patients who have not undergone immunotherapy can be affected by tumor purity. Interleukin 21 (IL21) had the highest weight in the inflammatory response model, suggesting its vital role in the prediction mode. Since the number of metastatic melanoma patients (n = 429) was relatively large among CIT cohorts, we further performed a co-culture experiment using a melanoma cell line and CD8 + T cell populations generated from peripheral blood monocytes. The results showed that IL21 therapy combined with anti-PD1 (programmed cell death 1) antibodies (trepril monoclonal antibodies) significantly enhanced the cytotoxic activity of CD8 + T cells against the melanoma cell line.

Conclusion

In this study, we developed an inflammatory response gene signature model that predicts patient survival and immunotherapy response in multiple malignancies. We further found that the predictive performance in the non-small cell lung cancer and gastric cancer group had the highest value among the six different malignancy subgroups. When compared with existing signatures, the inflammatory response gene signature scores for on-treatment samples were more robust predictors of the response to CIT in metastatic melanoma.

CHK1 inhibitor SRA737 is active in PARP inhibitor resistant and CCNE1 amplified ovarian cancer

High-grade serous ovarian cancers (HGSOCs) with homologous recombination deficiency (HRD) are initially responsive to poly (ADP-ribose) polymerase inhibitors (PARPi), but resistance ultimately emerges. HGSOC with CCNE1 amplification (CCNE1 amp) are associated with resistance to PARPi and platinum treatments. High replication stress in HRD and CCNE1 amp HGSOC leads to increased reliance on checkpoint kinase 1 (CHK1), a key regulator of cell cycle progression and the replication stress response. Here, we investigated the anti-tumor activity of the potent, highly selective, orally bioavailable CHK1 inhibitor (CHK1i), SRA737, in both acquired PARPi-resistant BRCA1/2 mutant and CCNE1 amp HGSOC models. We demonstrated that SRA737 increased replication stress and induced subsequent cell death in vitro. SRA737 monotherapy in vivo prolonged survival in CCNE1 amp models, suggesting a potential biomarker for CHK1i therapy. Combination SRA737 and PARPi therapy increased tumor regression in both PARPi-resistant and CCNE1 amp patient-derived xenograft models, warranting further study in these HGSOC subgroups.

Low-coverage whole genome sequencing for a highly selective cohort of severe COVID-19 patients

Despite the advances in genetic marker identification associated with severe COVID-19, the full genetic characterisation of the disease remains elusive. This study explores imputation in low-coverage whole genome sequencing for a severe COVID-19 patient cohort. We generated a dataset of 79 imputed variant call format files using the GLIMPSE1 tool, each containing an average of 9.5 million single nucleotide variants. Validation revealed a high imputation accuracy (squared Pearson correlation ≍0.97) across sequencing platforms, showcasing GLIMPSE1's ability to confidently impute variants with minor allele frequencies as low as 2% in individuals with Spanish ancestry. We carried out a comprehensive analysis of the patient cohort, examining hospitalisation and intensive care utilisation, sex and age-based differences, and clinical phenotypes using a standardised set of medical terms developed to characterise severe COVID-19 symptoms. The methods and findings presented here can be leveraged for future genomic projects to gain vital insights into health challenges like COVID-19.

Tackling neurodegeneration in vitro with omics: a path towards new targets and drugs

Drug discovery is a generally inefficient and capital-intensive process. For neurodegenerative diseases (NDDs), the development of novel therapeutics is particularly urgent considering the long list of late-stage drug candidate failures. Although our knowledge on the pathogenic mechanisms driving neurodegeneration is growing, additional efforts are required to achieve a better and ultimately complete understanding of the pathophysiological underpinnings of NDDs. Beyond the etiology of NDDs being heterogeneous and multifactorial, this process is further complicated by the fact that current experimental models only partially recapitulate the major phenotypes observed in humans. In such a scenario, multi-omic approaches have the potential to accelerate the identification of new or repurposed drugs against a multitude of the underlying mechanisms driving NDDs. One major advantage for the implementation of multi-omic approaches in the drug discovery process is that these overarching tools are able to disentangle disease states and model perturbations through the comprehensive characterization of distinct molecular layers (i.e., genome, transcriptome, proteome) up to a single-cell resolution. Because of recent advances increasing their affordability and scalability, the use of omics technologies to drive drug discovery is nascent, but rapidly expanding in the neuroscience field. Combined with increasingly advanced in vitro models, which particularly benefited from the introduction of human iPSCs, multi-omics are shaping a new paradigm in drug discovery for NDDs, from disease characterization to therapeutics prediction and experimental screening. In this review, we discuss examples, main advantages and open challenges in the use of multi-omic approaches for the in vitro discovery of targets and therapies against NDDs.

Aberrant spliceosome activity via elevated intron retention and upregulation and phosphorylation of SF3B1 in chronic lymphocytic leukemia

Splicing factor 3b subunit 1 (SF3B1) is the largest subunit and core component of the spliceosome. Inhibition of SF3B1 was associated with an increase in broad intron retention (IR) on most transcripts, suggesting that IR can be used as a marker of spliceosome inhibition in chronic lymphocytic leukemia (CLL) cells. Furthermore, we separately analyzed exonic and intronic mapped reads on annotated RNA-sequencing transcripts obtained from B cells (n = 98 CLL patients) and healthy volunteers (n = 9). We measured intron/exon ratio to use that as a surrogate for alternative RNA splicing (ARS) and found that 66% of CLL-B cell transcripts had significant IR elevation compared with normal B cells (NBCs) and that correlated with mRNA downregulation and low expression levels. Transcripts with the highest IR levels belonged to biological pathways associated with gene expression and RNA splicing. A >2-fold increase of active pSF3B1 was observed in CLL-B cells compared with NBCs. Additionally, when the CLL-B cells were treated with macrolides (pladienolide-B), a significant decrease in pSF3B1, but not total SF3B1 protein, was observed. These findings suggest that IR/ARS is increased in CLL, which is associated with SF3B1 phosphorylation and susceptibility to SF3B1 inhibitors. These data provide additional support to the relevance of ARS in carcinogenesis and evidence of pSF3B1 participation in this process.

Advancing rare cancer research by MALDI mass spectrometry imaging: Applications, challenges, and future perspectives in sarcoma

MALDI mass spectrometry imaging (MALDI imaging) uniquely advances cancer research, by measuring spatial distribution of endogenous and exogenous molecules directly from tissue sections. These molecular maps provide valuable insights into basic and translational cancer research, including tumor biology, tumor microenvironment, biomarker identification, drug treatment, and patient stratification. Despite its advantages, MALDI imaging is underutilized in studying rare cancers. Sarcomas, a group of malignant mesenchymal tumors, pose unique challenges in medical research due to their complex heterogeneity and low incidence, resulting in understudied subtypes with suboptimal management and outcomes. In this review, we explore the applicability of MALDI imaging in sarcoma research, showcasing its value in understanding this highly heterogeneous and challenging rare cancer. We summarize all MALDI imaging studies in sarcoma to date, highlight their impact on key research fields, including molecular signatures, cancer heterogeneity, and drug studies. We address specific challenges encountered when employing MALDI imaging for sarcomas, and propose solutions, such as using formalin-fixed paraffin-embedded tissues, and multiplexed experiments, and considerations for multi-site studies and digital data sharing practices. Through this review, we aim to spark collaboration between MALDI imaging researchers and clinical colleagues, to deploy the unique capabilities of MALDI imaging in the context of sarcoma.

Sharing sensitive data in life sciences: an overview of centralized and federated approaches

Biomedical data are generated and collected from various sources, including medical imaging, laboratory tests and genome sequencing. Sharing these data for research can help address unmet health needs, contribute to scientific breakthroughs, accelerate the development of more effective treatments and inform public health policy. Due to the potential sensitivity of such data, however, privacy concerns have led to policies that restrict data sharing. In addition, sharing sensitive data requires a secure and robust infrastructure with appropriate storage solutions. Here, we examine and compare the centralized and federated data sharing models through the prism of five large-scale and real-world use cases of strategic significance within the European data sharing landscape: the French Health Data Hub, the BBMRI-ERIC Colorectal Cancer Cohort, the federated European Genome-phenome Archive, the Observational Medical Outcomes Partnership/OHDSI network and the EBRAINS Medical Informatics Platform. Our analysis indicates that centralized models facilitate data linkage, harmonization and interoperability, while federated models facilitate scaling up and legal compliance, as the data typically reside on the data generator's premises, allowing for better control of how data are shared. This comparative study thus offers guidance on the selection of the most appropriate sharing strategy for sensitive datasets and provides key insights for informed decision-making in data sharing efforts.

Getting your DUCs in a row - standardising the representation of Digital Use Conditions

Improving patient care and advancing scientific discovery requires responsible sharing of research data, healthcare records, biosamples, and biomedical resources that must also respect applicable use conditions. Defining a standard to structure and manage these use conditions is a complex and challenging task. This is exemplified by a near unlimited range of asset types, a high variability of applicable conditions, and differing applications at the individual or collective level. Furthermore, the specifics and granularity required are likely to vary depending on the ultimate contexts of use. All these factors confound alignment of institutional missions, funding objectives, regulatory and technical requirements to facilitate effective sharing. The presented work highlights the complexity and diversity of the problem, reviews the current state of the art, and emphasises the need for a flexible and adaptable approach. We propose Digital Use Conditions (DUC) as a framework that addresses these needs by leveraging existing standards, striking a balance between expressiveness versus ambiguity, and considering the breadth of applicable information with their context of use.

Ensemble-based classification using microRNA expression identifies a breast cancer patient subgroup with an ultralow long-term risk of metastases

BACKGROUND

Current clinical markers overestimate the recurrence risk in many lymph node negative (LNN) breast cancer (BC) patients such that a majority of these low-risk patients unnecessarily receive systemic treatments. We tested if differential microRNA expression in primary tumors allows reliable identification of indolent LNN BC patients to provide an improved classification tool for overtreatment reduction in this patient group.

METHODS

We collected freshly frozen primary tumors of 80 LNN BC patients with recurrence and 80 recurrence-free patients (mean follow-up: 20.9 years). The study comprises solely systemically untreated patients to exclude that administered treatments confound the metastasis status. Samples were pairwise matched for clinical-pathological characteristics to minimize dependence of current markers. Patients were classified into risk-subgroups according to the differential microRNA expression of their tumors via classification model building with cross-validation using seven classification methods and a voting scheme. The methodology was validated using available data of two independent cohorts (n = 123, n = 339).

RESULTS

Of the 80 indolent patients (who would all likely receive systemic treatments today) our ultralow-risk classifier correctly identified 37 while keeping a sensitivity of 100% in the recurrence group. Multivariable logistic regression analysis confirmed independence of voting results from current clinical markers. Application of the method in two validation cohorts confirmed successful classification of ultralow-risk BC patients with significantly prolonged recurrence-free survival.

CONCLUSION

Profiles of differential microRNAs expression can identify LNN BC patients who could spare systemic treatments demanded by currently applied classifications. However, further validation studies are required for clinical implementation of the applied methodology.

Applied models and molecular characteristics of small cell lung cancer

Small cell lung cancer (SCLC) is a highly aggressive type of cancer frequently diagnosed with metastatic spread, rendering it surgically unresectable for the majority of patients. Although initial responses to platinum-based therapies are often observed, SCLC invariably relapses within months, frequently developing drug-resistance ultimately contributing to short overall survival rates. Recently, SCLC research aimed to elucidate the dynamic changes in the genetic and epigenetic landscape. These have revealed distinct subtypes of SCLC, each characterized by unique molecular signatures. The recent understanding of the molecular heterogeneity of SCLC has opened up potential avenues for precision medicine, enabling the development of targeted therapeutic strategies. In this review, we delve into the applied models and computational approaches that have been instrumental in the identification of promising drug candidates. We also explore the emerging molecular diagnostic tools that hold the potential to transform clinical practice and patient care.

Update of the Minimum Information About BIobank Data Sharing (MIABIS) Core Terminology to the 3rd Version

Introduction: The Minimum Information About BIobank Data Sharing (MIABIS) is a biobank-specific terminology enabling the sharing of biobank-related data for different purposes across a wide range of database implementations. After 4 years in use and with the first version of the individual-level MIABIS component Sample, Sample donor, and Event, it was necessary to revise the terminology, especially to include biobanks that work more in the data domain than with samples. Materials & Methods: Nine use-cases representing different types of biobanks, studies, and networks participated in the development work. They represent types of data, specific sample types, or levels of organization that were not included earlier in MIABIS. To support our revision of the Biobank entity, we conducted a survey of European biobanks to chart the services they provide. An important stakeholder group for biobanks include researchers as the main users of biobanks. To be able to render MIABIS more researcher-friendly, we collected different sample/data requests to analyze the terminology adjustment needs in detail. During the update process, the Core terminology was iteratively reviewed by a large group of experts until a consensus was reached. Results: With this update, MIABIS was adjusted to encompass data-driven biobanks and to include data collections, while also describing the services and capabilities biobanks offer to their users, besides the retrospective samples. The terminology was also extended to accommodate sample and data collections of nonhuman origin. Additionally, a set of organizational attributes was compiled to describe networks. Discussion: The usability of MIABIS Core v3 was increased by extending it to cover more topics of the biobanking domain. Additionally, the focus was on a more general terminology and harmonization of attributes with the individual-level entities Sample, Sample donor, and Event to keep the overall terminology minimal. With this work, the internal semantics of the MIABIS terminology was improved.

Refined preferences of prioritizers improve intelligent diagnosis for Mendelian diseases

Phenotype-guided gene prioritizers have proved a highly efficient approach to identifying causal genes for Mendelian diseases. In our previous study, we preliminarily evaluated the performance of ten prioritizers. However, all the selected software was run based on default settings and singleton mode. With a large-scale family dataset from Deciphering Developmental Disorders (DDD) project (N = 305) and an in-house trio cohort (N = 152), the four optimal performers in our prior study including Exomiser, PhenIX, AMELIE, and LIRCIAL were further assessed through parameter optimization and/or the utilization of trio mode. The in-depth assessment revealed high diagnostic yields of the four prioritizers with refined preferences, each alone or together: (1) 83.3-91.8% of the causal genes were presented among the first ten candidates in the final ranking lists of the four tools; (2) Over 97.7% of the causal genes were successfully captured within the top 50 by either of the four software. Exomiser did best in directly hitting the target (ranking the causal gene at the very top) while LIRICAL displayed a predominant overall detection capability. Besides, cases affected by low-penetrance and high-frequency pathogenic variants were found misjudged during the automated prioritization process. The discovery of the limitations shed light on the specific directions of future enhancement for causal-gene ranking tools.

Phenotypic similarity-based approach for variant prioritization for unsolved rare disease: a preliminary methodological report

Rare diseases (RD) have a prevalence of not more than 1/2000 persons in the European population, and are characterised by the difficulty experienced in obtaining a correct and timely diagnosis. According to Orphanet, 72.5% of RD have a genetic origin although 35% of them do not yet have an identified causative gene. A significant proportion of patients suspected to have a genetic RD receive an inconclusive exome/genome sequencing. Working towards the International Rare Diseases Research Consortium (IRDiRC)'s goal for 2027 to ensure that all people living with a RD receive a diagnosis within one year of coming to medical attention, the Solve-RD project aims to identify the molecular causes underlying undiagnosed RD. As part of this strategy, we developed a phenotypic similarity-based variant prioritization methodology comparing submitted cases with other submitted cases and with known RD in Orphanet. Three complementary approaches based on phenotypic similarity calculations using the Human Phenotype Ontology (HPO), the Orphanet Rare Diseases Ontology (ORDO) and the HPO-ORDO Ontological Module (HOOM) were developed; genomic data reanalysis was performed by the RD-Connect Genome-Phenome Analysis Platform (GPAP). The methodology was tested in 4 exemplary cases discussed with experts from European Reference Networks. Variants of interest (pathogenic or likely pathogenic) were detected in 8.8% of the 725 cases clustered by similarity calculations. Diagnostic hypotheses were validated in 42.1% of them and needed further exploration in another 10.9%. Based on the promising results, we are devising an automated standardized phenotypic-based re-analysis pipeline to be applied to the entire unsolved cases cohort.

Multi-omics approaches for biomarker discovery in predicting the response of esophageal cancer to neoadjuvant therapy: A multidimensional perspective

Neoadjuvant chemoradiotherapy (NCRT) followed by surgery has been established as the standard treatment strategy for operable locally advanced esophageal cancer (EC). However, achieving pathologic complete response (pCR) or near pCR to NCRT is significantly associated with a considerable improvement in survival outcomes, while pCR patients may help organ preservation for patients by active surveillance to avoid planned surgery. Thus, there is an urgent need for improved biomarkers to predict EC chemoradiation response in research and clinical settings. Advances in multiple high-throughput technologies such as next-generation sequencing have facilitated the discovery of novel predictive biomarkers, specifically based on multi-omics data, including genomic/transcriptomic sequencings and proteomic/metabolomic mass spectra. The application of multi-omics data has shown the benefits in improving the understanding of underlying mechanisms of NCRT sensitivity/resistance in EC. Particularly, the prominent development of artificial intelligence (AI) has introduced a new direction in cancer research. The integration of multi-omics data has significantly advanced our knowledge of the disease and enabled the identification of valuable biomarkers for predicting treatment response from diverse dimension levels, especially with rapid advances in biotechnological and AI methodologies. Herein, we summarize the current status of research on the use of multi-omics technologies in predicting NCRT response for EC patients. Current limitations, challenges, and future perspectives of these multi-omics platforms will be addressed to assist in experimental designs and clinical use for further integrated analysis.

Next-generation sequencing and bioinformatics in rare movement disorders

The ability to sequence entire exomes and genomes has revolutionized molecular testing in rare movement disorders, and genomic sequencing is becoming an integral part of routine diagnostic workflows for these heterogeneous conditions. However, interpretation of the extensive genomic variant information that is being generated presents substantial challenges. In this Perspective, we outline multidimensional strategies for genetic diagnosis in patients with rare movement disorders. We examine bioinformatics tools and computational metrics that have been developed to facilitate accurate prioritization of disease-causing variants. Additionally, we highlight community-driven data-sharing and case-matchmaking platforms, which are designed to foster the discovery of new genotype-phenotype relationships. Finally, we consider how multiomic data integration might optimize diagnostic success by combining genomic, epigenetic, transcriptomic and/or proteomic profiling to enable a more holistic evaluation of variant effects. Together, the approaches that we discuss offer pathways to the improved understanding of the genetic basis of rare movement disorders.

Review of Computational Methods and Database Sources for Predicting the Effects of Coding Frameshift Small Insertion and Deletion Variations

Genetic variations (including substitutions, insertions, and deletions) exert a profound influence on DNA sequences. These variations are systematically classified as synonymous, nonsynonymous, and nonsense, each manifesting distinct effects on proteins. The implementation of high-throughput sequencing has significantly augmented our comprehension of the intricate interplay between gene variations and protein structure and function, as well as their ramifications in the context of diseases. Frameshift variations, particularly small insertions and deletions (indels), disrupt protein coding and are instrumental in disease pathogenesis. This review presents a succinct review of computational methods, databases, current challenges, and future directions in predicting the consequences of coding frameshift small indels variations. We analyzed the predictive efficacy, reliability, and utilization of computational methods and variant account, reliability, and utilization of database. Besides, we also compared the prediction methodologies on GOF/LOF pathogenic variation data. Addressing the challenges pertaining to prediction accuracy and cross-species generalizability, nascent technologies such as AI and deep learning harbor immense potential to enhance predictive capabilities. The importance of interdisciplinary research and collaboration cannot be overstated for devising effective diagnosis, treatment, and prevention strategies concerning diseases associated with coding frameshift indels variations.

Expression Atlas update: insights from sequencing data at both bulk and single cell level

Expression Atlas (www.ebi.ac.uk/gxa) and its newest counterpart the Single Cell Expression Atlas (www.ebi.ac.uk/gxa/sc) are EMBL-EBI's knowledgebases for gene and protein expression and localisation in bulk and at single cell level. These resources aim to allow users to investigate their expression in normal tissue (baseline) or in response to perturbations such as disease or changes to genotype (differential) across multiple species. Users are invited to search for genes or metadata terms across species or biological conditions in a standardised consistent interface. Alongside these data, new features in Single Cell Expression Atlas allow users to query metadata through our new cell type wheel search. At the experiment level data can be explored through two types of dimensionality reduction plots, t-distributed Stochastic Neighbor Embedding (tSNE) and Uniform Manifold Approximation and Projection (UMAP), overlaid with either clustering or metadata information to assist users' understanding. Data are also visualised as marker gene heatmaps identifying genes that help confer cluster identity. For some data, additional visualisations are available as interactive cell level anatomograms and cell type gene expression heatmaps.

FAIR data retrieval for sensitive clinical research data in Galaxy

BACKGROUND

In clinical research, data have to be accessible and reproducible, but the generated data are becoming larger and analysis complex. Here we propose a platform for Findable, Accessible, Interoperable, and Reusable (FAIR) data access and creating reproducible findings. Standardized access to a major genomic repository, the European Genome-Phenome Archive (EGA), has been achieved with API services like PyEGA3. We aim to provide a FAIR data analysis service in Galaxy by retrieving genomic data from the EGA and provide a generalized "omics" platform for FAIR data analysis.

RESULTS

To demonstrate this, we implemented an end-to-end Galaxy workflow to replicate the findings from an RD-Connect synthetic dataset Beyond the 1 Million Genomes (synB1MG) available from the EGA. We developed the PyEGA3 connector within Galaxy to easily download multiple datasets from the EGA. We added the gene.iobio tool, a diagnostic environment for precision genomics, to Galaxy and demonstrate that it provides a more dynamic and interpretable view for trio analysis results. We developed a Galaxy trio analysis workflow to determine the pathogenic variants from the synB1MG trios using the GEMINI and gene.iobio tool. The complete workflow is available at WorkflowHub, and an associated tutorial was created in the Galaxy Training Network, which helps researchers unfamiliar with Galaxy to run the workflow.

CONCLUSIONS

We showed the feasibility of reusing data from the EGA in Galaxy via PyEGA3 and validated the workflow by rediscovering spiked-in variants in synthetic data. Finally, we improved existing tools in Galaxy and created a workflow for trio analysis to demonstrate the value of FAIR genomics analysis in Galaxy.

Ontologies for increasing the FAIRness of plant research data

The importance of improving the FAIRness (findability, accessibility, interoperability, reusability) of research data is undeniable, especially in the face of large, complex datasets currently being produced by omics technologies. Facilitating the integration of a dataset with other types of data increases the likelihood of reuse, and the potential of answering novel research questions. Ontologies are a useful tool for semantically tagging datasets as adding relevant metadata increases the understanding of how data was produced and increases its interoperability. Ontologies provide concepts for a particular domain as well as the relationships between concepts. By tagging data with ontology terms, data becomes both human- and machine- interpretable, allowing for increased reuse and interoperability. However, the task of identifying ontologies relevant to a particular research domain or technology is challenging, especially within the diverse realm of fundamental plant research. In this review, we outline the ontologies most relevant to the fundamental plant sciences and how they can be used to annotate data related to plant-specific experiments within metadata frameworks, such as Investigation-Study-Assay (ISA). We also outline repositories and platforms most useful for identifying applicable ontologies or finding ontology terms.

Loss of p53-DREAM-mediated repression of cell cycle genes as a driver of lymph node metastasis in head and neck cancer

BACKGROUND

The prognosis for patients with head and neck cancer (HNC) is poor and has improved little in recent decades, partially due to lack of therapeutic options. To identify effective therapeutic targets, we sought to identify molecular pathways that drive metastasis and HNC progression, through large-scale systematic analyses of transcriptomic data.

METHODS

We performed meta-analysis across 29 gene expression studies including 2074 primary HNC biopsies to identify genes and transcriptional pathways associated with survival and lymph node metastasis (LNM). To understand the biological roles of these genes in HNC, we identified their associated cancer pathways, as well as the cell types that express them within HNC tumor microenvironments, by integrating single-cell RNA-seq and bulk RNA-seq from sorted cell populations.

RESULTS

Patient survival-associated genes were heterogenous and included drivers of diverse tumor biological processes: these included tumor-intrinsic processes such as epithelial dedifferentiation and epithelial to mesenchymal transition, as well as tumor microenvironmental factors such as T cell-mediated immunity and cancer-associated fibroblast activity. Unexpectedly, LNM-associated genes were almost universally associated with epithelial dedifferentiation within malignant cells. Genes negatively associated with LNM consisted of regulators of squamous epithelial differentiation that are expressed within well-differentiated malignant cells, while those positively associated with LNM represented cell cycle regulators that are normally repressed by the p53-DREAM pathway. These pro-LNM genes are overexpressed in proliferating malignant cells of TP53 mutated and HPV + ve HNCs and are strongly associated with stemness, suggesting that they represent markers of pre-metastatic cancer stem-like cells. LNM-associated genes are deregulated in high-grade oral precancerous lesions, and deregulated further in primary HNCs with advancing tumor grade and deregulated further still in lymph node metastases.

CONCLUSIONS

In HNC, patient survival is affected by multiple biological processes and is strongly influenced by the tumor immune and stromal microenvironments. In contrast, LNM appears to be driven primarily by malignant cell plasticity, characterized by epithelial dedifferentiation coupled with EMT-independent proliferation and stemness. Our findings postulate that LNM is initially caused by loss of p53-DREAM-mediated repression of cell cycle genes during early tumorigenesis.

MiR-181a targets STING to drive PARP inhibitor resistance in BRCA- mutated triple-negative breast cancer and ovarian cancer

BACKGROUND

Poly (ADP-ribose) polymerase inhibitors (PARPi) are approved for the treatment of BRCA-mutated breast cancer (BC), including triple-negative BC (TNBC) and ovarian cancer (OvCa). A key challenge is to identify the factors associated with PARPi resistance; although, previous studies suggest that platinum-based agents and PARPi share similar resistance mechanisms.

METHODS

Olaparib-resistant (OlaR) cell lines were analyzed using HTG EdgeSeq miRNA Whole Transcriptomic Analysis (WTA). Functional assays were performed in three BRCA-mutated TNBC cell lines. In-silico analysis were performed using multiple databases including The Cancer Genome Atlas, the Genotype-Tissue Expression, The Cancer Cell Line Encyclopedia, Genomics of Drug Sensitivity in Cancer, and Gene Omnibus Expression.

RESULTS

High miR-181a levels were identified in OlaR TNBC cell lines (p = 0.001) as well as in tumor tissues from TNBC patients (p = 0.001). We hypothesized that miR-181a downregulates the stimulator of interferon genes (STING) and the downstream proinflammatory cytokines to mediate PARPi resistance. BRCA1 mutated TNBC cell lines with miR-181a-overexpression were more resistant to olaparib and showed downregulation in STING and the downstream genes controlled by STING. Extracellular vesicles derived from PARPi-resistant TNBC cell lines horizontally transferred miR-181a to parental cells which conferred PARPi-resistance and targeted STING. In clinical settings, STING levels were positively correlated with interferon gamma (IFNG) response scores (p = 0.01). In addition, low IFNG response scores were associated with worse response to neoadjuvant treatment including PARPi for high-risk HER2 negative BC patients (p = 0.001). OlaR TNBC cell lines showed resistance to platinum-based drugs. OvCa cell lines resistant to platinum showed resistance to olaparib. Knockout of miR-181a significantly improved olaparib sensitivity in OvCa cell lines (p = 0.001).

CONCLUSION

miR-181a is a key factor controlling the STING pathway and driving PARPi and platinum-based drug resistance in TNBC and OvCa. The miR-181a-STING axis can be used as a potential marker for predicting PARPi responses in TNBC and OvCa tumors.

Genomic Biomarker Discovery in Disease Progression and Therapy Response in Bladder Cancer Utilizing Machine Learning

Cancer in all its forms of expression is a major cause of death. To identify the genomic reason behind cancer, discovery of biomarkers is needed. In this paper, genomic data of bladder cancer are examined for the purpose of biomarker discovery. Genomic biomarkers are indicators stemming from the study of the genome, either at a very low level based on the genome sequence itself, or more abstractly such as measuring the level of gene expression for different disease groups. The latter method is pivotal for this work, since the available datasets consist of RNA sequencing data, transformed to gene expression levels, as well as data on a multitude of clinical indicators. Based on this, various methods are utilized such as statistical modeling via logistic regression and regularization techniques (elastic-net), clustering, survival analysis through Kaplan-Meier curves, and heatmaps for the experiments leading to biomarker discovery. The experiments have led to the discovery of two gene signatures capable of predicting therapy response and disease progression with considerable accuracy for bladder cancer patients which correlates well with clinical indicators such as Therapy Response and T-Stage at surgery with Disease Progression in a time-to-event manner.

Whole-genome sequencing across 449 samples spanning 47 ethnolinguistic groups provides insights into genetic diversity in Nigeria

African populations have been drastically underrepresented in genomics research, and failure to capture the genetic diversity across the numerous ethnolinguistic groups (ELGs) found on the continent has hindered the equity of precision medicine initiatives globally. Here, we describe the whole-genome sequencing of 449 Nigerian individuals across 47 unique self-reported ELGs. Population structure analysis reveals genetic differentiation among our ELGs, consistent with previous findings. From the 36 million SNPs and insertions or deletions (indels) discovered in our dataset, we provide a high-level catalog of both novel and medically relevant variation present across the ELGs. These results emphasize the value of this resource for genomics research, with added granularity by representing multiple ELGs from Nigeria. Our results also underscore the potential of using these cohorts with larger sample sizes to improve our understanding of human ancestry and health in Africa.

Identifying New Contributors to Brain Metastasis in Lung Adenocarcinoma: A Transcriptomic Meta-Analysis

Lung tumors frequently metastasize to the brain. Brain metastasis (BM) is common in advanced cases, and a major cause of patient morbidity and mortality. The precise molecular mechanisms governing BM are still unclear, in part attributed to the rarity of BM specimens. In this work, we compile a unique transcriptomic dataset encompassing RNA-seq, microarray, and single-cell analyses from BM samples obtained from patients with lung adenocarcinoma (LUAD). By integrating this comprehensive dataset, we aimed to enhance understanding of the molecular landscape of BM, thereby facilitating the identification of novel and efficient treatment strategies. We identified 102 genes with significantly deregulated expression levels in BM tissues, and discovered transcriptional alterations affecting the key driver 'hub' genes CD69 (a type II C-lectin receptor) and GZMA (Granzyme A), indicating an important role of the immune system in the development of BM from primary LUAD. Our study demonstrated a BM-specific gene expression pattern and revealed the presence of dendritic cells and neutrophils in BM, suggesting an immunosuppressive tumor microenvironment. These findings highlight key drivers of LUAD-BM that may yield therapeutic targets to improve patient outcomes.

Applications of multi-omics analysis in human diseases

Multi-omics usually refers to the crossover application of multiple high-throughput screening technologies represented by genomics, transcriptomics, single-cell transcriptomics, proteomics and metabolomics, spatial transcriptomics, and so on, which play a great role in promoting the study of human diseases. Most of the current reviews focus on describing the development of multi-omics technologies, data integration, and application to a particular disease; however, few of them provide a comprehensive and systematic introduction of multi-omics. This review outlines the existing technical categories of multi-omics, cautions for experimental design, focuses on the integrated analysis methods of multi-omics, especially the approach of machine learning and deep learning in multi-omics data integration and the corresponding tools, and the application of multi-omics in medical researches (e.g., cancer, neurodegenerative diseases, aging, and drug target discovery) as well as the corresponding open-source analysis tools and databases, and finally, discusses the challenges and future directions of multi-omics integration and application in precision medicine. With the development of high-throughput technologies and data integration algorithms, as important directions of multi-omics for future disease research, single-cell multi-omics and spatial multi-omics also provided a detailed introduction. This review will provide important guidance for researchers, especially who are just entering into multi-omics medical research.

Integrating Omics Data in Genome-Scale Metabolic Modeling: A Methodological Perspective for Precision Medicine

Recent advancements in omics technologies have generated a wealth of biological data. Integrating these data within mathematical models is essential to fully leverage their potential. Genome-scale metabolic models (GEMs) provide a robust framework for studying complex biological systems. GEMs have significantly contributed to our understanding of human metabolism, including the intrinsic relationship between the gut microbiome and the host metabolism. In this review, we highlight the contributions of GEMs and discuss the critical challenges that must be overcome to ensure their reproducibility and enhance their prediction accuracy, particularly in the context of precision medicine. We also explore the role of machine learning in addressing these challenges within GEMs. The integration of omics data with GEMs has the potential to lead to new insights, and to advance our understanding of molecular mechanisms in human health and disease.

Assessing transcriptomic reidentification risks using discriminative sequence models

Gene expression data provide molecular insights into the functional impact of genetic variation, for example, through expression quantitative trait loci (eQTLs). With an improving understanding of the association between genotypes and gene expression comes a greater concern that gene expression profiles could be matched to genotype profiles of the same individuals in another data set, known as a linking attack. Prior works show such a risk could analyze only a fraction of eQTLs that is independent owing to restrictive model assumptions, leaving the full extent of this risk incompletely understood. To address this challenge, we introduce the discriminative sequence model (DSM), a novel probabilistic framework for predicting a sequence of genotypes based on gene expression data. By modeling the joint distribution over all known eQTLs in a genomic region, DSM improves the power of linking attacks with necessary calibration for linkage disequilibrium and redundant predictive signals. We show greater linking accuracy of DSM compared with existing approaches across a range of attack scenarios and data sets including up to 22,288 individuals, suggesting that DSM helps uncover a substantial additional risk overlooked by previous studies. Our work provides a unified framework for assessing the privacy risks of sharing diverse omics data sets beyond transcriptomics.

The Human Extracellular Matrix Diseasome Reveals Genotype-Phenotype Associations with Clinical Implications for Age-Related Diseases

The extracellular matrix (ECM) is earning an increasingly relevant role in many disease states and aging. The analysis of these disease states is possible with the GWAS and PheWAS methodologies, and through our analysis, we aimed to explore the relationships between polymorphisms in the compendium of ECM genes (i.e., matrisome genes) in various disease states. A significant contribution on the part of ECM polymorphisms is evident in various types of disease, particularly those in the core-matrisome genes. Our results confirm previous links to connective-tissue disorders but also unearth new and underexplored relationships with neurological, psychiatric, and age-related disease states. Through our analysis of the drug indications for gene-disease relationships, we identify numerous targets that may be repurposed for age-related pathologies. The identification of ECM polymorphisms and their contributions to disease will play an integral role in future therapeutic developments, drug repurposing, precision medicine, and personalized care.

Consent Codes: Maintaining Consent in an Ever-expanding Open Science Ecosystem

We previously proposed a structure for recording consent-based data use 'categories' and 'requirements' - Consent Codes - with a view to supporting maximum use and integration of genomic research datasets, and reducing uncertainty about permissible re-use of shared data. Here we discuss clarifications and subsequent updates to the Consent Codes (v4) based on new areas of application (e.g., the neurosciences, biobanking, H3Africa), policy developments (e.g., return of research results), and further practical considerations, including developments in automated approaches to consent management.

Evaluation of Epigenetic Age Acceleration Scores and Their Associations with CVD-Related Phenotypes in a Population Cohort

We evaluated associations between nine epigenetic age acceleration (EAA) scores and 18 cardiometabolic phenotypes using an Eastern European ageing population cohort richly annotated for a diverse set of phenotypes (subsample, n = 306; aged 45-69 years). This was implemented by splitting the data into groups with positive and negative EAAs. We observed strong association between all EAA scores and sex, suggesting that any analysis of EAAs should be adjusted by sex. We found that some sex-adjusted EAA scores were significantly associated with several phenotypes such as blood levels of gamma-glutamyl transferase and low-density lipoprotein, smoking status, annual alcohol consumption, multiple carotid plaques, and incident coronary heart disease status (not necessarily the same phenotypes for different EAAs). We demonstrated that even after adjusting EAAs for sex, EAA-phenotype associations remain sex-specific, which should be taken into account in any downstream analysis involving EAAs. The obtained results suggest that in some EAA-phenotype associations, negative EAA scores (i.e., epigenetic age below chronological age) indicated more harmful phenotype values, which is counterintuitive. Among all considered epigenetic clocks, GrimAge was significantly associated with more phenotypes than any other EA scores in this Russian sample.

TogoVar: A comprehensive Japanese genetic variation database

TogoVar ( https://togovar.org ) is a database that integrates allele frequencies derived from Japanese populations and provides annotations for variant interpretation. First, a scheme to reanalyze individual-level genome sequence data deposited in the Japanese Genotype-phenotype Archive (JGA), a controlled-access database, was established to make allele frequencies publicly available. As more Japanese individual-level genome sequence data are deposited in JGA, the sample size employed in TogoVar is expected to increase, contributing to genetic study as reference data for Japanese populations. Second, public datasets of Japanese and non-Japanese populations were integrated into TogoVar to easily compare allele frequencies in Japanese and other populations. Each variant detected in Japanese populations was assigned a TogoVar ID as a permanent identifier. Third, these variants were annotated with molecular consequence, pathogenicity, and literature information for interpreting and prioritizing variants. Here, we introduce the newly developed TogoVar database that compares allele frequencies among Japanese and non-Japanese populations and describes the integrated annotations.

DNA Data Bank of Japan (DDBJ) update report 2022

The Bioinformation and DNA Data Bank of Japan (DDBJ) Center (https://www.ddbj.nig.ac.jp) maintains database archives that cover a wide range of fields in life sciences. As a founding member of the International Nucleotide Sequence Database Collaboration (INSDC), our primary mission is to collect and distribute nucleotide sequence data, as well as their study and sample information, in collaboration with the National Center for Biotechnology Information in the United States and the European Bioinformatics Institute. In addition to INSDC resources, the Center operates databases for functional genomics (GEA: Genomic Expression Archive), metabolomics (MetaboBank), and human genetic and phenotypic data (JGA: Japanese Genotype-Phenotype Archive). These databases are built on the supercomputer of the National Institute of Genetics, whose remaining computational capacity is actively utilized by domestic researchers for large-scale biological data analyses. Here, we report our recent updates and the activities of our services.

Missed colorectal cancers in a fecal immunochemical test-based screening program: Molecular profiling of interval carcinomas

BACKGROUND

For optimizing fecal immunochemical test (FIT)-based screening programs, reducing the rate of missed colorectal cancers (CRCs) by FIT (FIT-interval CRCs) is an important aspect. Knowledge of the molecular make-up of these missed lesions could facilitate more accurate detection of all (precursor) lesions.

AIM

To compare the molecular make-up of FIT-interval CRCs to lesions that are detected by FIT [screen-detected CRCs (SD-CRCs)].

METHODS

FIT-interval CRCs observed in a Dutch pilot-program of FIT-based screening were compared to a control group of SD-CRCs in a 1:2 ratio, resulting in 27 FIT-interval CRC and 54 SD-CRCs. Molecular analyses included microsatellite instability (MSI), CpG island methylator phenotype (CIMP), DNA sequence mutations and copy number alterations (CNAs).

RESULTS

Although no significant differences were reached, FIT-interval CRCs were more often CIMP positive and MSI positive (33% CIMP in FIT-interval CRCs vs 21% in SD-CRCs (P = 0.274); 19% MSI in FIT-interval CRCs vs 12% in SD-CRCs (P = 0.469)), and showed more often serrated pathway associated features such as BRAF (30% vs 12%, P = 0.090) and PTEN (15% vs 2.4%, P = 0.063) mutations. APC mutations, a classic feature of the adenoma-carcinoma-sequence, were more abundant in SD-CRCs (68% vs 40% in FIT-interval CRCs P = 0.035). Regarding CNAs differences between the two groups; FIT-interval CRCs less often showed gains at the regions 8p11.22-q24.3 (P = 0.009), and more often gains at 20p13-p12.1 (P = 0.039).

CONCLUSION

Serrated pathway associated molecular features seem to be more common in FIT-interval CRCs, while classic adenoma carcinoma pathway associated molecular features seem to be more common in SD-CRCs. This indicates that proximal serrated lesions may be overrepresented among FIT-interval CRCs.

Sensitive and reproducible cell-free methylome quantification with synthetic spike-in controls

Cell-free methylated DNA immunoprecipitation sequencing (cfMeDIP-seq) identifies genomic regions with DNA methylation, using a protocol adapted to work with low-input DNA samples and with cell-free DNA (cfDNA). We developed a set of synthetic spike-in DNA controls for cfMeDIP-seq to provide a simple and inexpensive reference for quantitative normalization. We designed 54 DNA fragments with combinations of methylation status (methylated and unmethylated), fragment length (80 bp, 160 bp, 320 bp), G + C content (35%, 50%, 65%), and fraction of CpG dinucleotides within the fragment (1/80 bp, 1/40 bp, 1/20 bp). Using 0.01 ng of spike-in controls enables training a generalized linear model that absolutely quantifies methylated cfDNA in MeDIP-seq experiments. It mitigates batch effects and corrects for biases in enrichment due to known biophysical properties of DNA fragments and other technical biases.

Investigation of Anti-Liver Cancer Activity of the Herbal Drug FDY003 Using Network Pharmacology

Globally, liver cancer (LC) is the sixth-most frequently occurring and the second-most fatal malignancy, responsible for 0.83 million deaths annually. Although the application of herbal drugs in cancer therapies has increased, their anti-LC activity and relevant mechanisms have not been fully studied from a systems perspective. To address these issues, we conducted a system-perspective network pharmacological investigation into the activity and mechanisms underlying the action of the herbal drug. FDY003 reduced the viability of human LC treatment. FDY003 reduced the viability of human LC cells and elevated their chemosensitivity. There were a total of 16 potential bioactive chemical components in FDY003 and they had 91 corresponding targets responsible for the pathological processes in LC. These FDY003 targets were functionally involved in regulating the survival, proliferation, apoptosis, and cell cycle of LC cells. Additionally, we found that FDY003 may target key signaling cascades connected to diverse LC pathological mechanisms, namely, PI3K-Akt, focal adhesion, IL-17, FoxO, MAPK, and TNF pathways. Overall, this study contributed to integrative mechanistic insights into the anti-LC potential of FDY003.

Teaching genomics to life science undergraduates using cloud computing platforms with open datasets

The final year of a biochemistry degree is usually a time to experience research. However, laboratory-based research projects were not possible during COVID-19. Instead, we used open datasets to provide computational research projects in metagenomics to biochemistry undergraduates (80 students with limited computing experience). We aimed to give the students a chance to explore any dataset, rather than use a small number of artificial datasets (~60 published datasets were used). To achieve this, we utilized Google Colaboratory (Colab), a virtual computing environment. Colab was used as a framework to retrieve raw sequencing data (analyzed with QIIME2) and generate visualizations. Setting up the environment requires no prior experience; all students have the same drive structure and notebooks can be shared (for synchronous sessions). We also used the platform to combine multiple datasets, perform a meta-analysis, and allowed the students to analyze large datasets with 1000s of subjects and factors. Projects that required increased computational resources were integrated with Google Cloud Compute. In future, all research projects can include some aspects of reanalyzing public data, providing students with data science experience. Colab is also an excellent environment in which to develop data skills in multiple languages (e.g., Perl, Python, Julia).

Genome-wide mutational signatures in low-coverage whole genome sequencing of cell-free DNA

Mutational signatures accumulate in somatic cells as an admixture of endogenous and exogenous processes that occur during an individual’s lifetime. Since dividing cells release cell-free DNA (cfDNA) fragments into the circulation, we hypothesize that plasma cfDNA might reflect mutational signatures. Point mutations in plasma whole genome sequencing (WGS) are challenging to identify through conventional mutation calling due to low sequencing coverage and low mutant allele fractions. In this proof of concept study of plasma WGS at 0.3–1.5x coverage from 215 patients and 227 healthy individuals, we show that both pathological and physiological mutational signatures may be identified in plasma. By applying machine learning to mutation profiles, patients with stage I-IV cancer can be distinguished from healthy individuals with an Area Under the Curve of 0.96. Interrogating mutational processes in plasma may enable earlier cancer detection, and might enable the assessment of cancer risk and etiology.

Detection of mutational signatures in cell-free DNA (cfDNA) is challenging due to low sequence coverage and low mutant allele fractions. Here, the authors identify mutational signatures in plasma whole genome sequencing of cancer patients and use machine learning to distinguish them from healthy individuals.

An in-library ligation strategy and its application in CRISPR/Cas9 screening of high-order gRNA combinations

Simultaneous targeting multiple genes is a big advantage of CRISPR (clustered regularly interspaced short palindromic repeats) genome editing but challenging to achieve in CRISPR screening. The crosstalk among genes or gene products is a common and fundamental mechanism to ensure cellular stability and functional diversity. However, the screening approach to map high-order gene combinations to the interesting phenotype is still lacking. Here, we developed a universal in-library ligation strategy and applied it to generate multiplexed CRISPR library, which could perturb four pre-designed targets in a cell. We conducted in vivo CRISPR screening for potential guide RNA (gRNA) combinations inducing anti-tumor immune responses. Simultaneously disturbing a combination of three checkpoints in CD8+ T cells was demonstrated to be more effective than disturbing Pdcd1 only for T cell activation in the tumor environment. This study developed a novel in-library ligation strategy to facilitate the multiplexed CRISPR screening, which could extend our ability to explore the combinatorial outcomes from coordinated gene behaviors.

Opportunities and challenges for the use of common controls in sequencing studies

Genome-wide association studies using large-scale genome and exome sequencing data have become increasingly valuable in identifying associations between genetic variants and disease, transforming basic research and translational medicine. However, this progress has not been equally shared across all people and conditions, in part due to limited resources. Leveraging publicly available sequencing data as external common controls, rather than sequencing new controls for every study, can better allocate resources by augmenting control sample sizes or providing controls where none existed. However, common control studies must be carefully planned and executed as even small differences in sample ascertainment and processing can result in substantial bias. Here, we discuss challenges and opportunities for the robust use of common controls in high-throughput sequencing studies, including study design, quality control and statistical approaches. Thoughtful generation and use of large and valuable genetic sequencing data sets will enable investigation of a broader and more representative set of conditions, environments and genetic ancestries than otherwise possible.

Genetic Profiling of Colorectal Carcinomas of Patients with Primary Sclerosing Cholangitis and Inflammatory Bowel Disease

BACKGROUND

Patients with primary sclerosing cholangitis (PSC) and inflammatory bowel disease (IBD) run a 10-fold increased risk of developing colorectal cancer (CRC) compared to patients with IBD only. The aim of this study was to perform an extensive screen of known carcinogenic genomic alterations in patients with PSC-IBD, and to investigate whether such changes occur already in nondysplastic mucosa.

METHODS

Archival cancer tissue and nondysplastic mucosa from resection specimens of 19 patients with PSC-IBD-CRC were characterized, determining DNA copy-number variations, microsatellite instability (MSI), mutations on 48 cancer genes, and CpG island methylator phenotype (CIMP). Genetic profiles were compared with 2 published cohorts of IBD-associated CRC (IBD-CRC; n = 11) and sporadic CRC (s-CRC; n = 100).

RESULTS

Patterns of chromosomal aberrations in PSC-IBD-CRC were similar to those observed in IBD-CRC and s-CRC, MSI occurred only once. Mutation frequencies were comparable between the groups, except for mutations in KRAS, which were less frequent in PSC-IBD-CRC (5%) versus IBD-CRC (38%) and s-CRC (31%; P = .034), and in APC, which were less frequent in PSC-IBD-CRC (5%) and IBD-CRC (0%) versus s-CRC (50%; P < .001). Cases of PSC-IBD-CRC were frequently CIMP positive (44%), at similar levels to cases of s-CRC (34%; P = .574) but less frequent than in cases with IBD-CRC (90%; P = .037). Similar copy number aberrations and mutations were present in matched cancers and adjacent mucosa in 5/15 and 7/11 patients, respectively.

CONCLUSIONS

The excess risk of CRC in patients with PSC-IBD was not explained by copy number aberrations, mutations, MSI, nor CIMP status, in cancer tissue, nor in adjacent mucosa. These findings set the stage for further exome-wide and epigenetic studies.

Exploration of the System-Level Mechanisms of the Herbal Drug FDY003 for Pancreatic Cancer Treatment: A Network Pharmacological Investigation

Pancreatic cancer (PC) is the most lethal cancer with the lowest survival rate globally. Although the prescription of herbal drugs against PC is gaining increasing attention, their polypharmacological therapeutic mechanisms are yet to be fully understood. Based on network pharmacology, we explored the anti-PC properties and system-level mechanisms of the herbal drug FDY003. FDY003 decreased the viability of human PC cells and strengthened their chemosensitivity. Network pharmacological analysis of FDY003 indicated the presence of 16 active phytochemical components and 123 PC-related pharmacological targets. Functional enrichment analysis revealed that the PC-related targets of FDY003 participate in the regulation of cell growth and proliferation, cell cycle process, cell survival, and cell death. In addition, FDY003 was shown to target diverse key pathways associated with PC pathophysiology, namely, the PIK3-Akt, MAPK, FoxO, focal adhesion, TNF, p53, HIF-1, and Ras pathways. Our network pharmacological findings advance the mechanistic understanding of the anti-PC properties of FDY003 from a system perspective.

A quality control portal for sequencing data deposited at the European genome-phenome archive

Since its launch in 2008, the European Genome-Phenome Archive (EGA) has been leading the archiving and distribution of human identifiable genomic data. In this regard, one of the community concerns is the potential usability of the stored data, as of now, data submitters are not mandated to perform any quality control (QC) before uploading their data and associated metadata information. Here, we present a new File QC Portal developed at EGA, along with QC reports performed and created for 1 694 442 files [Fastq, sequence alignment map (SAM)/binary alignment map (BAM)/CRAM and variant call format (VCF)] submitted at EGA. QC reports allow anonymous EGA users to view summary-level information regarding the files within a specific dataset, such as quality of reads, alignment quality, number and type of variants and other features. Researchers benefit from being able to assess the quality of data prior to the data access decision and thereby, increasing the reusability of data (https://ega-archive.org/blog/data-upcycling-powered-by-ega/).

FAIR Genomes metadata schema promoting Next Generation Sequencing data reuse in Dutch healthcare and research

The genomes of thousands of individuals are profiled within Dutch healthcare and research each year. However, this valuable genomic data, associated clinical data and consent are captured in different ways and stored across many systems and organizations. This makes it difficult to discover rare disease patients, reuse data for personalized medicine and establish research cohorts based on specific parameters. FAIR Genomes aims to enable NGS data reuse by developing metadata standards for the data descriptions needed to FAIRify genomic data while also addressing ELSI issues. We developed a semantic schema of essential data elements harmonized with international FAIR initiatives. The FAIR Genomes schema v1.1 contains 110 elements in 9 modules. It reuses common ontologies such as NCIT, DUO and EDAM, only introducing new terms when necessary. The schema is represented by a YAML file that can be transformed into templates for data entry software (EDC) and programmatic interfaces (JSON, RDF) to ease genomic data sharing in research and healthcare. The schema, documentation and MOLGENIS reference implementation are available at https://fairgenomes.org .

Whole-genome sequencing of 1,171 elderly admixed individuals from São Paulo, Brazil

As whole-genome sequencing (WGS) becomes the gold standard tool for studying population genomics and medical applications, data on diverse non-European and admixed individuals are still scarce. Here, we present a high-coverage WGS dataset of 1,171 highly admixed elderly Brazilians from a census-based cohort, providing over 76 million variants, of which ~2 million are absent from large public databases. WGS enables identification of ~2,000 previously undescribed mobile element insertions without previous description, nearly 5 Mb of genomic segments absent from the human genome reference, and over 140 alleles from HLA genes absent from public resources. We reclassify and curate pathogenicity assertions for nearly four hundred variants in genes associated with dominantly-inherited Mendelian disorders and calculate the incidence for selected recessive disorders, demonstrating the clinical usefulness of the present study. Finally, we observe that whole-genome and HLA imputation could be significantly improved compared to available datasets since rare variation represents the largest proportion of input from WGS. These results demonstrate that even smaller sample sizes of underrepresented populations bring relevant data for genomic studies, especially when exploring analyses allowed only by WGS.

Whole genome sequencing (WGS) data on non-European and admixed individuals remains scarce. Here, the authors analyse WGS data from 1,171 admixed elderly Brazilians from a census cohort, characterising population-specific genetic variation and exploring the clinical utility of this expanded dataset.

Integrating whole genome sequencing, methylation, gene expression, topological associated domain information in regulatory mutation prediction: a study of follicular lymphoma

A major challenge in human genetics is of the analysis of the interplay between genetic and epigenetic factors in a multifactorial disease like cancer. Here, a novel methodology is proposed to investigate genome-wide regulatory mechanisms in cancer, as studied with the example of follicular Lymphoma (FL). In a first phase, a new machine-learning method is designed to identify Differentially Methylated Regions (DMRs) by computing six attributes. In a second phase, an integrative data analysis method is developed to study regulatory mutations in FL, by considering differential methylation information together with DNA sequence variation, differential gene expression, 3D organization of genome (e.g., topologically associated domains), and enriched biological pathways. Resulting mutation block-gene pairs are further ranked to find out the significant ones. By this approach, BCL2 and BCL6 were identified as top-ranking FL-related genes with several mutation blocks and DMRs acting on their regulatory regions. Two additional genes, CDCA4 and CTSO, were also found in top rank with significant DNA sequence variation and differential methylation in neighboring areas, pointing towards their potential use as biomarkers for FL. This work combines both genomic and epigenomic information to investigate genome-wide gene regulatory mechanisms in cancer and contribute to devising novel treatment strategies.

Evaluation of phenotype-driven gene prioritization methods for Mendelian diseases

It’s challenging work to identify disease-causing genes from the next-generation sequencing (NGS) data of patients with Mendelian disorders. To improve this situation, researchers have developed many phenotype-driven gene prioritization methods using a patient’s genotype and phenotype information, or phenotype information only as input to rank the candidate’s pathogenic genes. Evaluations of these ranking methods provide practitioners with convenience for choosing an appropriate tool for their workflows, but retrospective benchmarks are underpowered to provide statistically significant results in their attempt to differentiate. In this research, the performance of ten recognized causal-gene prioritization methods was benchmarked using 305 cases from the Deciphering Developmental Disorders (DDD) project and 209 in-house cases via a relatively unbiased methodology. The evaluation results show that methods using Human Phenotype Ontology (HPO) terms and Variant Call Format (VCF) files as input achieved better overall performance than those using phenotypic data alone. Besides, LIRICAL and AMELIE, two of the best methods in our benchmark experiments, complement each other in cases with the causal genes ranked highly, suggesting a possible integrative approach to further enhance the diagnostic efficiency. Our benchmarking provides valuable reference information to the computer-assisted rapid diagnosis in Mendelian diseases and sheds some light on the potential direction of future improvement on disease-causing gene prioritization methods.

Integrative OMICS Data-Driven Procedure Using a Derivatized Meta-Analysis Approach

The wealth of high-throughput data has opened up new opportunities to analyze and describe biological processes at higher resolution, ultimately leading to a significant acceleration of scientific output using high-throughput data from the different omics layers and the generation of databases to store and report raw datasets. The great variability among the techniques and the heterogeneous methodologies used to produce this data have placed meta-analysis methods as one of the approaches of choice to correlate the resultant large-scale datasets from different research groups. Through multi-study meta-analyses, it is possible to generate results with greater statistical power compared to individual analyses. Gene signatures, biomarkers and pathways that provide new insights of a phenotype of interest have been identified by the analysis of large-scale datasets in several fields of science. However, despite all the efforts, a standardized regulation to report large-scale data and to identify the molecular targets and signaling networks is still lacking. Integrative analyses have also been introduced as complementation and augmentation for meta-analysis methodologies to generate novel hypotheses. Currently, there is no universal method established and the different methods available follow different purposes. Herein we describe a new unifying, scalable and straightforward methodology to meta-analyze different omics outputs, but also to integrate the significant outcomes into novel pathways describing biological processes of interest. The significance of using proper molecular identifiers is highlighted as well as the potential to further correlate molecules from different regulatory levels. To show the methodology’s potential, a set of transcriptomic datasets are meta-analyzed as an example.

CIBERER: Spanish National Network for Research on Rare Diseases: a highly productive collaborative initiative

CIBER (Center for Biomedical Network Research; Centro de Investigación Biomédica En Red) is a public national consortium created in 2006 under the umbrella of the Spanish National Institute of Health Carlos III (ISCIII). This innovative research structure comprises 11 different specific areas dedicated to the main public health priorities in the National Health System. CIBERER, the thematic area of CIBER focused on rare diseases (RDs) currently consists of 75 research groups belonging to universities, research centers, and hospitals of the entire country. CIBERER's mission is to be a center prioritizing and favoring collaboration and cooperation between biomedical and clinical research groups, with special emphasis on the aspects of genetic, molecular, biochemical, and cellular research of RDs. This research is the basis for providing new tools for the diagnosis and therapy of low‐prevalence diseases, in line with the International Rare Diseases Research Consortium (IRDiRC) objectives, thus favoring translational research between the scientific environment of the laboratory and the clinical setting of health centers. In this article, we intend to review CIBERER's 15‐year journey and summarize the main results obtained in terms of internationalization, scientific production, contributions toward the discovery of new therapies and novel genes associated to diseases, cooperation with patients' associations and many other topics related to RD research.

Genome and transcriptome profiling of spontaneous preterm birth phenotypes

Preterm birth (PTB) occurs before 37 weeks of gestation. Risk factors include genetics and infection/inflammation. Different mechanisms have been reported for spontaneous preterm birth (SPTB) and preterm birth following preterm premature rupture of membranes (PPROM). This study aimed to identify early pregnancy biomarkers of SPTB and PPROM from the maternal genome and transcriptome. Pregnant women were recruited at the Liverpool Women’s Hospital. Pregnancy outcomes were categorised as SPTB, PPROM (≤ 34 weeks gestation, n = 53), high-risk term (HTERM, ≥ 37 weeks, n = 126) or low-risk (no history of SPTB/PPROM) term (LTERM, ≥ 39 weeks, n = 188). Blood samples were collected at 16 and 20 weeks gestation from which, genome (UK Biobank Axiom array) and transcriptome (Clariom D Human assay) data were acquired. PLINK and R were used to perform genetic association and differential expression analyses and expression quantitative trait loci (eQTL) mapping. Several significant molecular signatures were identified across the analyses in preterm cases. Genome-wide significant SNP rs14675645 (ASTN1) was associated with SPTB whereas microRNA-142 transcript and PPARG1-FOXP3 gene set were associated with PPROM at week 20 of gestation and is related to inflammation and immune response. This study has determined genomic and transcriptomic candidate biomarkers of SPTB and PPROM that require validation in diverse populations.

The cell cycle-related genes RHAMM, AURKA, TPX2, PLK1, and PLK4 are associated with the poor prognosis of breast cancer patients

Breast cancer is the third most common type of cancer diagnosed. Cell cycle is a complex but highly organized and controlled process, in which normal cells sense mitogenic growth signals that instruct them to enter and progress through their cell cycle. This process culminates in cell division generating two daughter cells with identical amounts of genetic material. Uncontrolled proliferation is one of the hallmarks of cancer. In this study, we analyzed the expression of the cell cycle-related genes receptor for hyaluronan (HA)-mediated motility (RHAMM), AURKA, TPX2, PLK1, and PLK4 and correlated them with the prognosis in a collective of 3952 breast cancer patients. A high messenger RNA expression of all studied genes correlated with a poor prognosis. Stratifying the patients according to the expression of hormonal receptors, we found that in patients with estrogen and progesterone receptor-positive and human epithelial growth factor receptor 2-negative tumors, and Luminal A and Luminal B tumors, the expression of the five analyzed genes correlates with worse survival. qPCR analysis of a panel of breast cancer cell lines representative of major molecular subtypes indicated a predominant expression in the luminal subtype. In vitro experiments showed that radiation influences the expression of the five analyzed genes both in luminal and triple-negative model cell lines. Functional analysis of MDA-MB-231 cells showed that small interfering RNA knockdown of PLK4 and TPX2 and pharmacological inhibition of PLK1 had an impact on the cell cycle and colony formation. Looking for a potential upstream regulation by microRNAs, we observed a differential expression of RHAMM, AURKA, TPX2, PLK1, and PLK4 after transfecting the MDA-MB-231 cells with three different microRNAs. Survival analysis of miR-34c-5p, miR-375, and miR-142-3p showed a different impact on the prognosis of breast cancer patients. Our study suggests that RHAMM, AURKA, TPX2, PLK1, and PLK4 can be used as potential targets for treatment or as a prognostic value in breast cancer patients.