Search

We present avidity sequencing, a sequencing chemistry that separately optimizes the processes of stepping along a DNA template and that of identifying each nucleotide within the template. Nucleotide identification uses multivalent nucleotide ligands on dye-labeled cores to form polymerase–polymer–nucleotide complexes bound to clonal copies of DNA targets. These polymer–nucleotide substrates, termed avidites, decrease the required concentration of reporting nucleotides from micromolar to nanomolar and yield negligible dissociation rates. Avidity sequencing achieves high accuracy, with 96.2% and 85.4% of base calls having an average of one error per 1,000 and 10,000 base pairs, respectively. We show that the average error rate of avidity sequencing remained stable following a long homopolymer.

We present a novel sequencing chemistry implemented as part of the AVITI system. Relying on the proximal DNA binding sites created through DNA amplification on a solid support, avidity sequencing uses multivalent nucleotide ligands on dye-labeled cores to simultaneously form polymerase-polymer nucleotide complexes bound to clonal copies of DNA targets. These polymer-nucleotide substrates, termed avidites, decrease the required concentration of reporting nucleotides by 100x and yield a negligible dissociation rate. We demonstrate the use of avidites within a novel sequencing technology that surpasses Q40 accuracy and enables a diversity of applications that include single cell RNA-seq and whole human genome sequencing.

Motivation Data normalization is essential to ensure accurate inference and comparability of gene expression measures across samples or conditions. Ideally, gene expression data should be rescaled based on consistently expressed reference genes. However, to normalize biologically diverse samples, the most commonly used reference genes exhibit striking expression variability and size-factor or distribution-based normalization methods can be problematic when the amount of asymmetry in differential expression is significant. Results We report an efficient and accurate data-driven method—Cosine score-based iterative normalization (Cosbin)—to normalize biologically diverse samples. Based on the Cosine scores of cross-condition expression patterns, the Cosbin pipeline iteratively eliminates asymmetric differentially expressed genes, identifies consistently expressed genes, and calculates sample-wise normalization factors. We demonstrate the superior performance and enhanced utility of Cosbin compared with six representative peer methods using both simulation and real multi-omics expression datasets. Implemented in open-source R scripts and specifically designed to address normalization bias due to significant asymmetry in differential expression across multiple conditions, the Cosbin tool complements rather than replaces the existing methods and will allow biologists to more accurately detect true molecular signals among diverse phenotypic groups. Availability and implementation The R scripts of Cosbin pipeline are freely available at https://github.com/MinjieSh/Cosbin. Supplementary information Supplementary data are available at Bioinformatics Advances online.

Motivation Ideally, a molecularly distinct subtype would be composed of molecular features that are expressed uniquely in the subtype of interest but in no others—so-called marker genes (MGs). MG plays a critical role in the characterization, classification or deconvolution of tissue or cell subtypes. We and others have recognized that the test statistics used by most methods do not exactly satisfy the MG definition and often identify inaccurate MG. Results We report an efficient and accurate data-driven method, formulated as a Cosine-based One-sample Test (COT) in scatter space, to detect MG among many subtypes using subtype expression profiles. Fundamentally different from existing approaches, the test statistic in COT precisely matches the mathematical definition of an ideal MG. We demonstrate the performance and utility of COT on both simulated and real gene expression and proteomics data. The open source Python/R tool will allow biologists to efficiently detect MG and perform a more comprehensive and unbiased molecular characterization of tissue or cell subtypes in many biomedical contexts. Nevertheless, COT complements not replaces existing methods. Availability and implementation The Python COT software with a detailed user’s manual and a vignette are freely available at https://github.com/MintaYLu/COT. Supplementary information Supplementary data are available at Bioinformatics Advances online.

Missing values are a major issue in quantitative proteomics analysis. While many methods have been developed for imputing missing values in high-throughput proteomics data, a comparative assessment of imputation accuracy remains inconclusive, mainly because mechanisms contributing to true missing values are complex and existing evaluation methodologies are imperfect. Moreover, few studies have provided an outlook of future methodological development. We first re-evaluate the performance of eight representative methods targeting three typical missing mechanisms. These methods are compared on both simulated and masked missing values embedded within real proteomics datasets, and performance is evaluated using three quantitative measures. We then introduce fused regularization matrix factorization, a low-rank global matrix factorization framework, capable of integrating local similarity derived from additional data types. We also explore a biologically-inspired latent variable modeling strategy - convex analysis of mixtures - for missing value imputation and present preliminary experimental results. While some winners emerged from our comparative assessment, the evaluation is intrinsically imperfect because performance is evaluated indirectly on artificial missing or masked values not authentic missing values. Nevertheless, we show that our fused regularization matrix factorization provides a novel incorporation of external and local information, and the exploratory implementation of convex analysis of mixtures presents a biologically plausible new approach.

Complex tissues are composite ecological systems whose components interact with each other to create a unique physiological or pathophysiological state distinct from that found in other tissue microenvironments. To explore this ground yet dynamic state, molecular profiling of bulk tissues and mathematical deconvolution can be jointly used to characterize heterogeneity as an aggregate of molecularly distinct tissue or cell subtypes. We first introduce an efficient and fully unsupervised deconvolution method, namely the Convex Analysis of Mixtures – CAM3.0, that may aid biologists to confirm existing or generate novel scientific hypotheses about complex tissues in many biomedical contexts. We then evaluate the CAM3.0 functional pipelines using both simulations and benchmark data. We also report diverse case studies on bulk tissues with unknown number, proportion and expression patterns of the molecular archetypes. Importantly, these preliminary results support the concept that expression patterns of molecular archetypes often reflect the interactive not individual contributions of many known or novel cell types, and unsupervised deconvolution would be more powerful in uncovering novel multicellular or subcellular archetypes.

Background: Missing values are a major issue in quantitative proteomics data analysis. While many methods have been developed for imputing missing values in high-throughput proteomics data, comparative assessment on the accuracy of existing methods remains inconclusive, mainly because the true missing mechanisms are complex and the existing evaluation methodologies are imperfect. Moreover, few studies have provided an outlook of current and future development.Results: We first report an assessment of eight representative methods collectively targeting three typical missing mechanisms. The selected methods are compared on both realistic simulation and real proteomics datasets, and the performance is evaluated using three quantitative measures. We then discuss fused regularization matrix factorization, a popular low-rank matrix factorization framework with similarity and/or biological regularization, which is extendable to integrating multi-omics data such as gene expressions or clinical variables. We further explore the potential application of convex analysis of mixtures, a biologically-inspired latent variable modeling strategy, to missing value imputation. The preliminary results on proteomics data are provided together with an outlook into future development directions.Conclusion: While a few winners emerged from our comparative assessment, data-driven evaluation of imputation methods is imperfect because performance is evaluated indirectly on artificial missing or masked values not authentic missing values. Imputation accuracy may vary with signal intensity. Fused regularization matrix factorization provides a possibility of incorporating external information. Convex analysis of mixtures presents a biologically plausible new approach.

MotivationComplex biological tissues are often a heterogeneous mixture of several molecularly distinct cell or tissue subtypes. Both subtype compositions and expressions in individual samples can vary across different biological states or conditions. Computational deconvolution aims to dissect patterns of bulk gene expression data into subtype compositions and subtype-specific expressions. Typically, existing deconvolution methods can only estimate averaged subtype-specific expressions in a population, while detecting differential expressions or co-expression networks in particular subtypes requires unique subtype expression estimates in individual samples. Different from population-level deconvolution, however, individual-level deconvolution is mathematically an underdetermined problem because there are more variables than observations. ResultsWe report a sample-wise Convex Analysis of Mixtures (swCAM) method that can estimate subtype proportions and subtype-specific expressions in individual samples from bulk tissue transcriptomes. We extend our previous CAM framework to include a new term accounting for between-sample variations and formulate swCAM as a nuclear-norm and ℓ2,1-norm regularized matrix factorization problem. We determine hyperparameter values using a cross-validation scheme with random entry exclusion and obtain a swCAM solution using an efficient alternating direction method of multipliers. The swCAM is implemented in open-source R scripts. Experimental results on realistic simulation data show that swCAM can accurately estimate subtype-specific expressions in individual samples and successfully extract co-expression networks in particular subtypes that are otherwise unobtainable using bulk expression data. Application of swCAM to bulk-tissue data of 320 samples from bipolar disorder patients and controls identified changes in cell proportions, expression and coexpression modules in patient neurons. Mitochondria related genes showed significant changes suggesting an important role of energy dysregulation in bipolar disorder.Availability and implementationThe R Scripts of swCAM is freely available at https://github.com/Lululuella/swCAM. A user’s guide and a vignette are provided.Contactyuewang@vt.eduSupplementary informationSupplementary data are available at Bioinformatics online.

We develop an accurate and efficient method to detect marker genes among many subtypes using subtype-enriched expression profiles. We implement a Cosine based One-sample Test (COT) Python software that is easy to use and applicable to multi-omics data. We demonstrate the performance and utility of COT on gene expression and proteomics data acquired from tissue or cell subtypes. Formulated as a one-sample test with Cosine similarity test statistic in scatter space, the detected de novo marker genes will allow biologists to perform a more comprehensive and unbiased molecular characterization, deconvolution and classification of complex tissue or cell subtypes.

Among multiple subtypes of tissue or cell, subtype-specific differentially-expressed genes (SDEGs) are defined as being most-upregulated in only one subtype but not in any other. Detecting SDEGs plays a critical role in the molecular characterization and deconvolution of multicellular complex tissues. Classic differential analysis assumes a null hypothesis whose test statistic is not subtype-specific, thus can produce a high false positive rate and/or lower detection power. Here we first introduce a One-Versus-Everyone Fold Change (OVE-FC) test for detecting SDEGs. We then propose a scaled test statistic (OVE-sFC) for assessing the statistical significance of SDEGs that applies a mixture null distribution model and a tailored permutation test. The OVE-FC/sFC test was validated on both type 1 error rate and detection power using extensive simulation data sets generated from real gene expression profiles of purified subtype samples. The OVE-FC/sFC test was then applied to two benchmark gene expression data sets of purified subtype samples and detected many known or previously unknown SDEGs. Subsequent supervised deconvolution results on synthesized bulk expression data, obtained using the SDEGs detected from the independent purified expression data by the OVE-FC/sFC test, showed superior performance in deconvolution accuracy when compared with popular peer methods. National Institutes of HealthUnited States Department of Health & Human ServicesNational Institutes of Health (NIH) - USA [HL111362-05A1, HL133932, W81XWH-18-1-0723, BC171885P1, U01NS115658-01] Published version This work was funded in part by the National Institutes of Health under Grants HL111362-05A1, HL133932, W81XWH-18-1-0723 (BC171885P1), and U01NS115658-01.

Summary We develop a fully unsupervised deconvolution method to dissect complex tissues into molecularly distinctive tissue or cell subtypes based on bulk expression profiles. We implement an R package, deconvolution by Convex Analysis of Mixtures (debCAM) that can automatically detect tissue/cell-specific markers, determine the number of constituent subtypes, calculate subtype proportions in individual samples and estimate tissue/cell-specific expression profiles. We demonstrate the performance and biomedical utility of debCAM on gene expression, methylation, proteomics and imaging data. With enhanced data preprocessing and prior knowledge incorporation, debCAM software tool will allow biologists to perform a more comprehensive and unbiased characterization of tissue remodeling in many biomedical contexts. Availability and implementation http://bioconductor.org/packages/debCAM. Supplementary information Supplementary data are available at Bioinformatics online.

Motivation Synapses are essential to neural signal transmission. Therefore, quantification of synapses and related neurites from images is vital to gain insights into the underlying pathways of brain functionality and diseases. Despite the wide availability of synaptic punctum imaging data, several issues are impeding satisfactory quantification of these structures by current tools. First, the antibodies used for labeling synapses are not perfectly specific to synapses. These antibodies may exist in neurites or other cell compartments. Second, the brightness of different neurites and synaptic puncta is heterogeneous due to the variation of antibody concentration and synapse-intrinsic differences. Third, images often have low signal to noise ratio due to constraints of experiment facilities and availability of sensitive antibodies. These issues make the detection of synapses challenging and necessitates developing a new tool to easily and accurately quantify synapses. Results We present an automatic probability-principled synapse detection algorithm and integrate it into our synapse quantification tool SynQuant. Derived from the theory of order statistics, our method controls the false discovery rate and improves the power of detecting synapses. SynQuant is unsupervised, works for both 2D and 3D data, and can handle multiple staining channels. Through extensive experiments on one synthetic and three real datasets with ground truth annotation or manually labeling, SynQuant was demonstrated to outperform peer specialized unsupervised synapse detection tools as well as generic spot detection methods. Availability and implementation Java source code, Fiji plug-in, and test data are available at https://github.com/yu-lab-vt/SynQuant. Supplementary information Supplementary data are available at Bioinformatics online.

MotivationLiquid chromatography - mass spectrometry (LC-MS) is a standard method for proteomics and metabolomics analysis of biological samples. Unfortunately, it suffers from various changes in the retention times (RT) of the same compound in different samples, and these must be subsequently corrected (aligned) during data processing. Classic alignment methods such as in the popular XCMS package often assume a single time-warping function for each sample. Thus, the potentially varying RT drift for compounds with different masses in a sample is neglected in these methods. Moreover, the systematic change in RT drift across run order is often not considered by alignment algorithms. Therefore, these methods cannot completely correct misalignments. For a large-scale experiment involving many samples, the existence of misalignment becomes inevitable and concerning.ResultsHere we describe an integrated reference-free profile alignment method, neighbor-wise compound-specific Graphical Time Warping (ncGTW), that can detect misaligned features and align profiles by leveraging expected RT drift structures and compound-specific warping functions. Specifically, ncGTW uses individualized warping functions for different compounds and assigns constraint edges on warping functions of neighboring samples. Validated with both realistic synthetic data and internal quality control samples, ncGTW applied to two large-scale metabolomics LC-MS datasets identifies many misaligned features and successfully realigns them. These features would otherwise be discarded or uncorrected using existing methods. The ncGTW software tool is developed currently as a plug-in to the XCMS package.Availability and ImplementationAn R package of ncGTW is freely available at https://github.com/ChiungTingWu/ncGTW. A detailed user’s manual and a vignette are provided within the package.Contactyuewang@vt.edu, yug@vt.eduSupplementary informationSupplementary data are available at Bioinformatics online.

MotivationIdentification of biological pathways plays a central role in understanding both human health and diseases. Although much work has previously been done to explore the biological pathways by using single omics data, little effort has been reported using multi-omics data integration, mainly due to methodological and technological limitations. Compared to single omics data, multi-omics data will help identifying disease specific functional pathways with both higher sensitivity and specificity, thus gaining more comprehensive insights into the molecular architecture of disease processes.ResultsIn this paper, we propose two computational approaches that integrate multi-omics data and identify disease-specific biological pathways with high sensitivity and specificity. Applying our methods to an experimental multi-omics data dataset on muscular dystrophy subtypes, we identified disease-specific pathways of high biological plausibility. The developed methodology will likely have a broad impact on improving the molecular characterization of many common diseases.Contactyuewang@vt.eduSupplementary informationSupplementary information attached.

Most genetic or environmental factors work together in determining complex disease risk. Detecting gene-environment interactions may allow us to elucidate novel and targetable molecular mechanisms on how environmental exposures modify genetic effects. Unfortunately, standard logistic regression (LR) assumes a convenient mathematical structure for the null hypothesis that however results in both poor detection power and type 1 error, and is also susceptible to missing factor, imperfect surrogate, and disease heterogeneity confounding effects. Here we describe a new baseline framework, the asymmetric independence model (AIM) in case-control studies, and provide mathematical proofs and simulation studies verifying its validity across a wide range of conditions. We show that AIM mathematically preserves the asymmetric nature of maintaining health versus acquiring a disease, unlike LR, and thus is more powerful and robust to detect synergistic interactions. We present examples from four clinically discrete domains where AIM identified interactions that were previously either inconsistent or recognized with less statistical certainty. National Institutes of Health [HL111362, HL133932, BC171885P1, U24CA160036-05S1, MH110504] This work was supported by the National Institutes of Health under Grants HL111362, HL133932, BC171885P1, U24CA160036-05S1, and MH110504.

We demonstrate a novel method for reducing saturation artifacts in spectral-domain optical coherence tomography (SD-OCT) systems. This method is based on a two-level SD-OCT system with a dual-line charge-coupled device (CCD) camera. We compensate the saturated signal detected by the first line using the unsaturated signal detected by the second line. The Fourier transform of the compensated spectrum shows effective suppression of saturation artifacts. This method was also successfully performed on phantom material and skin on a human finger. Our method causes neither back-scattering power loss nor signal-to-noise ratio (SNR) degradation. The only difference between the traditional system and our two-level system is our utilization of the dual-line CCD camera; no additional devices or complex designs are needed.

A procedure for computer analyzing an optical coherence tomography (OCT) image of normal and precancerous oral mucosae is demonstrated to reasonably plot the boundary between epithelium (EP) and lamina propria (LP) layers, determine the EP thickness, and estimate the range of dysplastic cell distribution based on standard deviation (SD) mapping. In this study, 54 normal oral mucosa, 39 oral mild dysplasia, and 44 oral moderate dysplasia OCT images are processed for evaluating the diagnosis statistics. Based on SD mapping in an OCT image, it is found that the laterally average range percentages of 70% SD maximum level in the EP layer is a reasonably good threshold for differentiating moderate dysplasia from mild dysplasia oral lesion based on the OCT image analysis. The sensitivity and specificity in diagnosis statistics can reach 82 and 90%, respectively.

Saturation artifact is a significant image noise due to the limited detectable maximum light intensity of the CCD camera in the spectrometer module of spectral-domain optical coherence tomography (SD-OCT). This effect makes the OCT image appear with saturation noise in certain areas along with even more obvious vertical lines of noise. The unwanted artifacts cover the real sample structure information and can render the image unclear. This artifact problem is significant when using highly focused objective lenses for improving image transverse resolution. In this study, we developed a parabolic-based reconstruction method to estimate the saturated parts of the measured spectra in order to obtain the corresponding pseudo spectra. The saturation artifact is clearly reduced when the OCT images are processed from the pseudo spectra.

The significantly less stringent operation of a two-reference swept-source optical coherence tomography (OCT) system for suppressing the mirror image is demonstrated based on the spatially localized image processing method. With this method, the phase difference between the two reference signals is not limited to 90 degrees. Based on the current experimental operation, the mirror image can be effectively suppressed as long as the phase difference is larger than 20 degrees. In other words, the adjustment of the beam splitter orientation for controlling the phase difference becomes much more flexible. Also, based on a phantom experiment, the combination the spatially localized mirror image suppression method with the two-reference OCT operation leads to the implementation of full-range optical Doppler tomography.

An improved image processing procedure for suppressing the phase noise due to a motion artifact acquired during optical coherence tomography scanning and effectively illustrating the blood vessel distribution in a living tissue is demonstrated. This new processing procedure and the widely used procedure for micro-angiography application are based on the selection of high-frequency components in the spatial-frequency spectrum of B-mode scanning (x-space), which are contributed from the image portions of moving objects. However, by switching the processing order between the x-space and k-space, the new processing procedure shows the superior function of effectively suppressing the phase noise due to a motion artifact. After the blood vessel positions are precisely acquired based on the new processing procedure, the projected blood flow speed can be more accurately calibrated based on a previously reported method. The demonstrated new procedure is useful for clinical micro-angiography application, in which a stepping motor of generating motion artifacts is usually used in the scanning probe.

A method, novel to our knowledge, for effective mirror image suppression in Fourier-domain optical coherence tomography based on a phase shift between neighboring A-mode scans is demonstrated. By realizing that the phase shifts of the real and mirror images are mutually reversed and assuming that the real image intensities of the two successive A-mode scans are the same, we can solve a set of two coupled equations to obtain the real image signals. The images based on the scanning of a high-resolution spectral-domain optical coherence tomography system are processed to show effective mirror image suppression results. Compared with a similar method of broad application, our approach has the advantages of shorter process time and higher flexibility in selecting the concerned image portions for processing.

The theory and experimental results of a computation time-saving mirror image suppression method in Fourier-domain optical coherence tomography, which utilizes the property of reversed system phase shift between the real and mirror images, for differentiating one from the other are demonstrated. By solving a set of two equations based on a reasonable approximation, the real image signal can be obtained. The theoretical backgrounds and the improved real image quality of the average and iteration procedures in this method are particularly illustrated. Also, the mirror image suppression ratios under various process conditions, including different process iteration numbers and different system phase shifts between two neighboring A-mode scans, are evaluated. Meanwhile, the mirror image suppression results based on our method are compared with those obtained from the widely used BM-scan technique. It is found that when a process procedure of two iterations is used, the mirror image suppression quality based on our method can be higher than that obtained from the BM-scan technique. The computation time of our method is significantly shorter than that of the BM-scan technique.