Search

The lack of samples for generating standardized DNA datasets for setting up a sequencing pipeline or benchmarking the performance of different algorithms limits the implementation and uptake of cancer genomics. Here, we describe reference call sets obtained from paired tumor–normal genomic DNA (gDNA) samples derived from a breast cancer cell line—which is highly heterogeneous, with an aneuploid genome, and enriched in somatic alterations—and a matched lymphoblastoid cell line. We partially validated both somatic mutations and germline variants in these call sets via whole-exome sequencing (WES) with different sequencing platforms and targeted sequencing with >2,000-fold coverage, spanning 82% of genomic regions with high confidence. Although the gDNA reference samples are not representative of primary cancer cells from a clinical sample, when setting up a sequencing pipeline, they not only minimize potential biases from technologies, assays and informatics but also provide a unique resource for benchmarking ‘tumor-only’ or ‘matched tumor–normal’ analyses.

Clinical applications of precision oncology require accurate tests that can distinguish true cancer-specific mutations from errors introduced at each step of next-generation sequencing (NGS). To date, no bulk sequencing study has addressed the effects of cross-site reproducibility, nor the biological, technical and computational factors that influence variant identification. Here we report a systematic interrogation of somatic mutations in paired tumor–normal cell lines to identify factors affecting detection reproducibility and accuracy at six different centers. Using whole-genome sequencing (WGS) and whole-exome sequencing (WES), we evaluated the reproducibility of different sample types with varying input amount and tumor purity, and multiple library construction protocols, followed by processing with nine bioinformatics pipelines. We found that read coverage and callers affected both WGS and WES reproducibility, but WES performance was influenced by insert fragment size, genomic copy content and the global imbalance score (GIV; G > T/C > A). Finally, taking into account library preparation protocol, tumor content, read coverage and bioinformatics processes concomitantly, we recommend actionable practices to improve the reproducibility and accuracy of NGS experiments for cancer mutation detection.

As an indispensable tool for transcriptome-wide analysis of differential gene expression, RNA sequencing (RNAseq) has demonstrated great potential in clinical applications such as companion diagnostics and prognostics. However, there is a lack of certified RNA reference materials and the corresponding reference datasets of differential expression for assessing the reliability of RNAseq for its intended use in detecting intrinsically small biological differences in clinical settings such as molecular subtyping of diseases. As part of the Quartet Project for quality control and data integration of multiomics profiling, we established four RNA reference materials derived from immortalized B-lymphoblastoid cell lines from four members of a monozygotic twin family. Additionally, we constructed ratio-based transcriptome-wide reference datasets between two reference samples, providing "ground truth" for cross-platform and cross-lab proficiency test. Moreover, Quartet-sample-based quality metrics were developed for assessing the reliability of RNAseq technology in terms of intra-batch measurement and cross-batch data integration. The small intrinsic biological differences among the Quartet samples enabled sensitive assessment of performance of transcriptomic measurements and their cross-batch integration at the ratio level. The Quartet RNA reference materials combined with the ratio-based reference datasets can serve as unique resources for assessing data quality and improving reliability of transcriptomic profiling.

Multiomics profiling is a powerful tool to characterize the same samples with complementary features orchestrating the genome, epigenome, transcriptome, proteome, and metabolome. However, the lack of ground truth hampers the objective assessment of and subsequent choice from a plethora of measurement and computational methods aiming to integrate diverse and often enigmatically incomparable omics datasets. Here we establish and characterize the first suites of publicly available multiomics reference materials of matched DNA, RNA, proteins, and metabolites derived from immortalized cell lines from a family quartet of parents and monozygotic twin daughters, providing built-in truth defined by family relationship and the central dogma. We demonstrate that the "ratio"-based omics profiling data, i.e., by scaling the absolute feature values of a study sample relative to those of a concurrently measured universal reference sample, were inherently much more reproducible and comparable across batches, labs, platforms, and omics types, thus empower the horizontal (within-omics) and vertical (cross-omics) data integration in multiomics studies. Our study identifies "absolute" feature quantitation as the root cause of irreproducibility in multiomics measurement and data integration, and urges a paradigm shift from "absolute" to "ratio"-based multiomics profiling with universal reference materials.

Batch effects are notorious technical variations that are common in multiomic data and may result in misleading outcomes. With the era of big data, tackling batch effects in multiomic integration is urgently needed. As part of the Quartet Project for quality control and data integration of multiomic profiling, we comprehensively assess the performances of seven batch-effect correction algorithms (BECAs) for mitigating the negative impact of batch effects in multiomic datasets, including transcriptomics, proteomics, and metabolomics. Performances are evaluated based on accuracy of identifying differentially expressed features, robustness of predictive models, and the ability of accurately clustering cross-batch samples into their biological sample groups. Ratio-based method is more effective and widely applicable than others, especially in cases when batch effects are highly confounded with biological factors of interests. We further provide practical guidelines for the implementation of ratio-based method using universal reference materials profiled with study samples. Our findings show the promise for eliminating batch effects and enhancing data integration in increasingly large-scale, cross-batch multiomic studies.

The implementation of quality control for multiomic data requires the widespread use of well-characterized reference materials, reference datasets, and related resources. The Quartet Data Portal was built to facilitate community access to such rich resources established in the Quartet Project. A convenient platform is provided for users to request the DNA, RNA, protein, and metabolite reference materials, as well as multi-level datasets generated across omics, platforms, labs, protocols, and batches. Interactive visualization tools are offered to assist users to gain a quick understanding of the reference datasets. Crucially, the Quartet Data Portal continuously collects, evaluates, and integrates the community-generated data of the distributed Quartet multiomic reference materials. In addition, the portal provides analysis pipelines to assess the quality of user-submitted multiomic data. Furthermore, the reference datasets, performance metrics, and analysis pipelines will be improved through periodic review and integration of multiomic data submitted by the community. Effective integration of the evolving technologies via active interactions with the community will help ensure the reliability of multiomics-based biological discoveries. The Quartet Data Portal is accessible athttps://chinese-quartet.org.Graphical Abstract

Coronavirus disease 2019 (COVID-19) is an ongoing pandemic and there is an urgent need for safe and effective drugs for COVID-19 treatment. Since developing a new drug is time consuming, many approved or investigational drugs have been repurposed for COVID-19 treatment in clinical trials. Therefore, selection of safe drugs for COVID-19 patients is vital for combating this pandemic. Our goal was to evaluate the safety concerns of drugs by analyzing adverse events reported in post-market surveillance. We collected 296 drugs that have been evaluated in clinical trials for COVID-19 and identified 28,597,464 associated adverse events at the system organ classes (SOCs) level in the FDA adverse events report systems (FAERS). We calculated Z-scores of SOCs that statistically quantify the relative frequency of adverse events of drugs in FAERS to quantitatively measure safety concerns for the drugs. Analyzing the Z-scores revealed that these drugs are associated with different significantly frequent adverse events. Our results suggest that this safety concern metric may serve as a tool to inform selection of drugs with favorable safety profiles for COVID-19 patients in clinical practices. Caution is advised when administering drugs with high Z-scores to patients who are vulnerable to associated adverse events.

Since growing evidence has manifested that bisphenol A (BPA) may adversely affect human health, numerous BPA replacement compounds have been gradually introduced into the industry. Although BPA rep...

The rat is an important model organism in biomedical research for studying human disease mechanisms and treatments, but its annotated transcriptome is far from complete. We constructed a Rat Transcriptome Re-annotation named RTR using RNA-seq data from 320 samples in 11 different organs generated by the SEQC consortium. Totally, there are 52 807 genes and 114 152 transcripts in RTR. Transcribed regions and exons in RTR account for ∼42% and ∼6.5% of the genome, respectively. Of all 73 074 newly annotated transcripts in RTR, 34 213 were annotated as high confident coding transcripts and 24 728 as high confident long noncoding transcripts. Different tissues rather than different stages have a significant influence on the expression patterns of transcripts. We also found that 11 715 genes and 15 852 transcripts were expressed in all 11 tissues and that 849 house-keeping genes expressed different isoforms among tissues. This comprehensive transcriptome is freely available at http://www.unimd.org/rtr/. Our new rat transcriptome provides essential reference for genetics and gene expression studies in rat disease and toxicity models.

Liver toxicity is a major adverse drug reaction that accounts for drug failure in clinical trials and withdrawal from the market. Therefore, predicting potential liver toxicity at an early stage in drug discovery is crucial to reduce costs and the potential for drug failure. However, current in vivo animal toxicity testing is very expensive and time consuming. As an alternative approach, various machine learning models have been developed to predict potential liver toxicity in humans. This chapter reviews current advances in the development and application of machine learning models for prediction of potential liver toxicity in humans and discusses possible improvements to liver toxicity prediction.

The wide application of nanomaterials in consumer and medical products has raised concerns about their potential adverse effects on human health. Thus, more and more biological assessments regarding the toxicity of nanomaterials have been performed. However, the different ways the evaluations were performed, such as the utilized assays, cell lines, and the differences of the produced nanoparticles, make it difficult for scientists to analyze and effectively compare toxicities of nanomaterials. Fortunately, machine learning has emerged as a powerful tool for the prediction of nanotoxicity based on the available data. Among different types of toxicity assessments, nanomaterial cytotoxicity was the focus here because of the high sensitivity of cytotoxicity assessment to different treatments without the need for complicated and time-consuming procedures. In this review, we summarized recent studies that focused on the development of machine learning models for prediction of cytotoxicity of nanomaterials. The goal was to provide insight into predicting potential nanomaterial toxicity and promoting the development of safe nanomaterials.

Endocrine-disrupting chemicals (EDCs) can inadvertently interact with 12 classic nuclear receptors (NRs) that disrupt the endocrine system and cause adverse effects. There is no widely accepted understanding about what structural features make thousands of EDCs able to activate different NRs as well as how these structural features exert their functions and induce different outcomes at the cellular level. This paper applies the hierarchical characteristic fragment methodology and high-throughput screening molecular docking to comprehensively explore the structural and functional features of EDCs for the 12 NRs based on more than 7000 chemicals from curated datasets. EDCs share three levels of key fragments. The primary and secondary fragments are associated with the binding of EDCs to four groups of receptors: steroidal nuclear receptors (SNRs, including androgen, estrogen, glucocorticoid, mineralocorticoid, and progesterone), retinoic acid receptors, thyroid hormone receptors, and vitamin D receptors. The tertiary fragments determine the activity type by interacting with two key locations in the ligand-binding domains of NRs (N-H5-H3-C and N-H7-H11-C for SNRs and N-H5-H5'-H2'-H3-C and N-H6'-H11-C for non-SNRs). The resulting compiled structural fragments of EDCs together with elucidated compound NR binding modes provide a framework for understanding the interactions between EDCs and NRs, facilitating faster and more accurate screening of EDCs for multiple NRs in the future.

Since November 2021, Omicron has been the dominant severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variant that causes the coronavirus disease 2019 (COVID-19) and has continuously impacted human health. Omicron sublineages are still increasing and cause increased transmission and infection rates. The additional 15 mutations on the receptor binding domain (RBD) of Omicron spike proteins change the protein conformation, enabling the Omicron variant to evade neutralizing antibodies. For this reason, many efforts have been made to design new antigenic variants to induce effective antibodies in SARS-CoV-2 vaccine development. However, understanding the different states of Omicron spike proteins with and without external molecules has not yet been addressed. In this review, we analyze the structures of the spike protein in the presence and absence of angiotensin-converting enzyme 2 (ACE2) and antibodies. Compared to previously determined structures for the wildtype spike protein and other variants such as alpha, beta, delta, and gamma, the Omicron spike protein adopts a partially open form. The open-form spike protein with one RBD up is dominant, followed by the open-form spike protein with two RBD up, and the closed-form spike protein with the RBD down. It is suggested that the competition between antibodies and ACE2 induces interactions between adjacent RBDs of the spike protein, which lead to a partially open form of the Omicron spike protein. The comprehensive structural information of Omicron spike proteins could be helpful for the efficient design of vaccines against the Omicron variant.

BackgroundReproducible detection of inherited variants with whole genome sequencing (WGS) is vital for the implementation of precision medicine and is a complicated process in which each step affects variant call quality. Systematically assessing reproducibility of inherited variants with WGS and impact of each step in the process is needed for understanding and improving quality of inherited variants from WGS.ResultsTo dissect the impact of factors involved in detection of inherited variants with WGS, we sequence triplicates of eight DNA samples representing two populations on three short-read sequencing platforms using three library kits in six labs and call variants with 56 combinations of aligners and callers. We find that bioinformatics pipelines (callers and aligners) have a larger impact on variant reproducibility than WGS platform or library preparation. Single-nucleotide variants (SNVs), particularly outside difficult-to-map regions, are more reproducible than small insertions and deletions (indels), which are least reproducible when > 5 bp. Increasing sequencing coverage improves indel reproducibility but has limited impact on SNVs above 30×.ConclusionsOur findings highlight sources of variability in variant detection and the need for improvement of bioinformatics pipelines in the era of precision medicine with WGS.

With the rapid advancement of sequencing technologies, next generation sequencing (NGS) analysis has been widely applied in cancer genomics research. More recently, NGS has been adopted in clinical oncology to advance personalized medicine. Clinical applications of precision oncology require accurate tests that can distinguish tumor-specific mutations from artifacts introduced during NGS processes or data analysis. Therefore, there is an urgent need to develop best practices in cancer mutation detection using NGS and the need for standard reference data sets for systematically measuring accuracy and reproducibility across platforms and methods. Within the SEQC2 consortium context, we established paired tumor-normal reference samples and generated whole-genome (WGS) and whole-exome sequencing (WES) data using sixteen library protocols, seven sequencing platforms at six different centers. We systematically interrogated somatic mutations in the reference samples to identify factors affecting detection reproducibility and accuracy in cancer genomes. These large cross-platform/site WGS and WES datasets using well-characterized reference samples will represent a powerful resource for benchmarking NGS technologies, bioinformatics pipelines, and for the cancer genomics studies., Measurement(s)Somatic Mutation AnalysisTechnology Type(s)whole genome sequencing • Whole Exome SequencingFactor Type(s)sequencing platform • sample prepration • library preparation • bioinformatics methodSample Characteristic - OrganismHomo sapiens Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.16713655

Peroxisome proliferator activator receptor gamma (PPARγ) agonist activity of chemicals is an indicator of concerned health conditions such as fatty liver and obesity. In silico screening PPARγ agonists based on quantitative structure-activity relationship (QSAR) models could serve as an efficient and pragmatic strategy. Owing to the broad research interests in discovery of PPARγ-targeted drugs, a large amount of PPARγ agonist activity data has been produced in the field of medicinal chemistry, facilitating development of robust QSAR models. In this study, random forest classifiers were developed based on the binary-category data transformed from the heterogeneous PPARγ agonist activity data of drug-like compounds. Coupling with applicability domains, capability of the established classifiers for predicting environmental chemicals was evaluated using two external data sets. Our results demonstrated that applicability domains could enhance application of the developed classifiers to predict environmental PPARγ agonists.

Estrogen receptor alpha (ERα) is a ligand-dependent transcriptional factor in the nuclear receptor superfamily. Many structures of ERα bound with agonists and antagonists have been determined. However, the dynamic binding patterns of agonists and antagonists in the binding site of ERα remains unclear. Therefore, we performed molecular docking, molecular dynamics (MD) simulations, and quantum mechanical calculations to elucidate agonist and antagonist dynamic binding patterns in ERα. 17β-estradiol (E2) and 4-hydroxytamoxifen (OHT) were docked in the ligand binding pockets of the agonist and antagonist bound ERα. The best complex conformations from molecular docking were subjected to 100 nanosecond MD simulations. Hierarchical clustering was conducted to group the structures in the trajectory from MD simulations. The representative structure from each cluster was selected to calculate the binding interaction energy value for elucidation of the dynamic binding patterns of agonists and antagonists in the binding site of ERα. The binding interaction energy analysis revealed that OHT binds ERα more tightly in the antagonist conformer, while E2 prefers the agonist conformer. The results may help identify ERα antagonists as drug candidates and facilitate risk assessment of chemicals through ER-mediated responses.

Opioid receptors (OPRs) are the main targets for the treatment of pain and related disorders. The opiate compounds that activate these receptors are effective analgesics but their use leads to adverse effects and they often are highly addictive drugs of abuse. There is an urgent need for alternative chemicals that are analgesic and reducing/avoiding the unwanted effects in order to relieve the public health crisis of opioid addiction. Here, we aim to develop computational models to predict the OPR activity of small molecule compounds based on chemical structures and apply these models to identify novel OPR active compounds. We used four different machine learning algorithms to build models based on quantitative high throughput screening (qHTS) datasets of three OPRs in both agonist and antagonist modes. The best performing models were applied to virtually screen a large collection of compounds. The model predicted active compounds which were experimentally validated using the same qHTS assays that generated the training data. Random forest was the best classifier with the highest performance metrics and the mu OPR (OPRM)-agonist model achieved the best performance with AUC-ROC (0.88) and MCC (0.7) values. The model predicted actives resulted in hit rates ranging from 2.3% (delta OPR-agonist) to 15.8% (OPRM-agonist) after experimental confirmation. Comparing to the original assay hit rate, all models enriched hit rate by ≥ 2-fold. Our approach produced robust OPR prediction models that can be applied to prioritize compounds identified from large libraries for further experimental validation. The models identified several novel potent compounds as activators/inhibitors of OPRs that were confirmed experimentally. The potent hits were further investigated using molecular docking to find the interactions of the novel ligands in the active site of the corresponding OPR.

Chemicals may cause adverse effects on human health through binding to peroxisome proliferator-activated receptor γ (PPARγ). Hence, binding affinity is useful for evaluating chemicals with potential endocrine-disrupting effects. Quantitative structure-activity relationship (QSAR) regression models with defined applicability domains (ADs) are important to enable efficient screening of chemicals with PPARγ binding activity. However, lack of large data sets hindered the development of QSAR models. In this study, based on PPARγ binding affinity data sets curated from various sources, 30 QSAR models were developed using molecular fingerprints, two-dimensional descriptors, and five machine learning algorithms. Structure-activity landscapes (SALs) of the training compounds were described by network-like similarity graphs (NSGs). Based on the NSGs, local discontinuity scores were calculated and found to be positively correlated with the cross-validation absolute prediction errors of the models using the different training sets, descriptors, and algorithms. Moreover, innovative ADs were defined based on pairwise similarities between compounds and were found to outperform some conventional ADs. The curated data sets and developed regression models could be useful for evaluating PPARγ-involved adverse effects of chemicals. The SAL analysis and the innovative ADs could facilitate understanding of prediction results from QSAR models.

Severe acute respiratory syndrome coronavirus 2 (SARS−CoV−2) has caused the ongoing global COVID-19 pandemic that began in late December 2019. The rapid spread of SARS−CoV−2 is primarily due to person-to-person transmission. To understand the epidemiological traits of SARS−CoV−2 transmission, we conducted phylogenetic analysis on genome sequences from &gt, 54K SARS−CoV−2 cases obtained from two public databases. Hierarchical clustering analysis on geographic patterns in the resulting phylogenetic trees revealed a co-expansion tendency of the virus among neighboring countries with diverse sources and transmission routes for SARS−CoV−2. Pairwise sequence similarity analysis demonstrated that SARS−CoV−2 is transmitted locally and evolves during transmission. However, no significant differences were seen among SARS−CoV−2 genomes grouped by host age or sex. Here, our identified epidemiological traits provide information to better prevent transmission of SARS−CoV−2 and to facilitate the development of effective vaccines and therapeutics against the virus.

Background Prescription opioids are powerful pain-reducing medications, but they may cause a variety of adverse effects. Long-term prescription opioid use (POU) is contributing to an opioid-related epidemic of addiction and death, and the scope of the opioid crisis continues to expand. As such, there is a need to identify the adverse effects associated with prescription opioid use (POU). Thousands of articles that focus on POU and its associated medical disorders have been published. However, it is time-consuming and labor-intensive to extract and understand the information of all POU-related published articles.Methods In this study, we applied the well-adapted topic modeling method, Latent Dirichlet Allocation (LDA), to perform text mining on POU-related literature. We compiled six large academic abstract datasets by searching PubMed using the Medical Subject Headings (MeSH): prescription opioid, codeine, morphine, hydrocodone, oxycodone, and methadone. We then applied topic modeling to identify topics and analyze topic similarities/differences in these six datasets. Word clouds and histograms were used to depict the distribution of vocabularies over each topic in which the most prevalent words conveyed a topic’s substance.Results The LDA topics recaptured the search keywords in PubMed, and further revealed relevant themes, such as patients, drugs, side effects, and association links between different POU and risk factors, such as gender and age. Moreover, based on the topic modeling results, TreeMap was used to fingerprint abstracts, which revealed the possibility of constructing a visualized literature index by combining topic modeling and visualization tools such as TreeMap. Meanwhile, while performing trend analysis to explore the prevalent topic dynamics in the POU-related literature, we found that an increasing trend in opioid prescription and its associated health risks are assessed as the most central issues.Conclusion The topic modeling results presented in this study not only convey an understandable and thematic structure of the POU literature, but also provide a means to discover which documents contain information about medical disorders associated with POU, thus, reducing the time and effort needed to review the literature for relevant articles. These results can be used as a preliminary study to systematically understand the risk factors related to increased POU-associated medical disorders.

Indeterminate acute liver failure (IND-ALF) is a rare clinical syndrome with a high mortality rate. Lacking a known etiology makes rapid evaluation and treatment difficult, with liver transplantation often considered as the only therapeutic option. Our aim was to identify genetic variants from whole exome sequencing data that might be associated with IND-ALF clinical outcomes.Bioinformatics analysis was performed on whole exome sequencing data for 22 patients with IND-ALF. A 2-tier approach was used to identify significant single-nucleotide polymorphisms (SNPs) associated with IND-ALF clinical outcomes. Tier 1 identified the SNPs with a higher relative risk in the IND-ALF population compared with those identified in control populations. Tier 2 determined the SNPs connected to transplant-free survival and associated with model for end-stage liver disease serum sodium and Acute Liver Failure Study Group prognostic scores.Thirty-one SNPs were found associated with a higher relative risk in the IND-ALF population compared with those in controls, of which 11 belong to the human leukocyte antigen (HLA) class II genes but none for the class I. Further analysis showed that 5 SNPs: rs796202376, rs139189937, and rs113473719 of HLA-DRB5; rs9272712 of HLA-DQA1; and rs747397929 of IDO1 were associated with a higher probability of IND-ALF transplant-free survival. Using 3 selected SNPs, a model for the polygenic risk score was developed to predict IND-ALF prognoses, which are comparable with those by model for end-stage liver disease serum sodium and Acute Liver Failure Study Group prognostic scores.Certain gene variants in HLA-DRB5, HLA-DQA1, and IDO1 were found associated with IND-ALF transplant-free survival. Once validated, these identified SNPs may help elucidate the mechanism of IND-ALF and assist in its diagnosis and management.

Serious harmful effects have been reported for thiophenols, which are widely used industrial materials. To date, little information is available on whether such chemicals can elicit endocrine-related detrimental effects. Herein the potential binding affinity and underlying mechanism of action between human transthyretin (hTTR) and seven halogenated-thiophenols were examined experimentally and computationally. Experimental results indicated that the halogenated-thiophenols, except for pentafluorothiophenol, were powerful hTTR binders. The differentiated hTTR binding affinity of halogenated-thiophenols and halogenated-phenols were observed. The hTTR binding affinity of mono- and di-halo-thiophenols was higher than that of corresponding phenols; while the opposite relationship was observed for tri- and penta-halo-thiophenols and phenols. Our results also confirmed that the binding interactions were influenced by the degree of ligand dissociation. Molecular modeling results implied that the dominant noncovalent interactions in the molecular recognition processes between hTTR and halogenated-thiophenols were ionic pair, hydrogen bonds and hydrophobic interactions. Finally, a model with acceptable predictive ability was developed, which can be used to computationally predict the potential hTTR binding affinity of other halogenated-thiophenols and phenols. Taken together, our results highlighted that more research is needed to determine their potential endocrine-related harmful effects and appropriate management actions should be taken to promote their sustainable use.

Recent studies have revealed that circulating microRNAs are promising biomarkers for detecting toxicity or disease. Quantitative real-time polymerase chain reaction (qPCR) is often used to measure the levels of microRNAs. Besides complete and certain data, investigators inevitably have observed technically incomplete or uncertain qPCR data. Investigators usually set incomplete observations equal to the maximum quality number of qPCR cycles, apply the complete-observation method, or choose not to analyze targets with incomplete observations. Using biostatistical knowledge and published studies, we show that three commonly applied methods tend to cause biased inference and decrease reproducibility in biomarker detection. More efforts are needed to address the challenges to identify and detect reliable, novel circulating biomarkers in liquid biopsies.