Estimation with Uncertainty via Conditional Generative Adversarial Networks.

Conventional predictive Artificial Neural Networks (ANNs) commonly employ deterministic weight matrices; therefore, their prediction is a point estimate. Such a deterministic nature in ANNs causes the limitations of using ANNs for medical diagnosis, law problems, and portfolio management in which not only discovering the prediction but also the uncertainty of the prediction is essentially required. In order to address such a problem, we propose a predictive probabilistic neural network model, which corresponds to a different manner of using the generator in the conditional Generative Adversarial Network (cGAN) that has been routinely used for conditional sample generation. By reversing the input and output of ordinary cGAN, the model can be successfully used as a predictive model; moreover, the model is robust against noises since adversarial training is employed. In addition, to measure the uncertainty of predictions, we introduce the entropy and relative entropy for regression problems and classification problems, respectively. The proposed framework is applied to stock market data and an image classification task. As a result, the proposed framework shows superior estimation performance, especially on noisy data; moreover, it is demonstrated that the proposed framework can properly estimate the uncertainty of predictions.

Prediction of survival risks with adjusted gene expression through risk-gene networks.

MOTIVATION: Network-based analysis of biomedical data has been extensively studied over the last decades. As a successful application, gene networks have been used to illustrate interactions among genes and explain the associated phenotypes. However, the gene network approaches have not been actively applied for survival analysis, which is one of the main interests of biomedical research. In addition, a few previous studies using gene networks for survival analysis construct networks mainly from prior knowledge, such as pathways, regulations and gene sets, while the performance considerably depends on the selection of prior knowledge. RESULTS: In this paper, we propose a data-driven construction method for survival risk-gene networks as well as a survival risk prediction method using the network structure. The proposed method constructs risk-gene networks with survival-associated genes using penalized regression. Then, gene expression indices are hierarchically adjusted through the networks to reduce the variance intrinsic in datasets. By illustrating risk-gene structure, the proposed method is expected to provide an intuition for the relationship between genes and survival risks. The risk-gene network is applied to a low grade glioma dataset, and produces a hypothesis of the relationship between genetic biomarkers of low and high grade glioma. Moreover, with multiple datasets, we demonstrate that the proposed method shows superior prediction performance compared to other conventional methods. AVAILABILITY AND IMPLEMENTATION: The R package of risk-gene networks is freely available in the web at http://cdal.korea.ac.kr/NetDA/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

mEBT: multiple-matching evidence-based translator of murine genomic responses for human immunity studies.

Summary: In this paper, we introduce multiple-matching Evidence-based Translator (mEBT) to discover genomic responses from murine expression data for human immune studies, which are significant in the given condition of mice and likely have similar responses in the corresponding condition of human. mEBT is evaluated over multiple datasets and shows improved inter-species agreement. mEBT is expected to be useful for research groups who use murine models to study human immunity. Availability and implementation: http://cdal.korea.ac.kr/mebt/. Supplementary information: Supplementary data are available at Bioinformatics online.

GAIT: Gene expression Analysis for Interval Time.

Motivation: Despite the potential usefulness, the association analysis of gene expression with interval times of two events has been hampered because the occurrence of events can be censored and the conventional survival analysis is not suitable to handle two censored events. However, the recent advances of multivariate survival analysis considering multiple censored events together provide an unprecedented chance for this problem. Based on such advances, we have developed a software tool, GAIT, for the association analysis of gene expression with interval time of two events. Results: The performance of GAIT was demonstrated by simulation studies and the real data analysis. The result indicates the usefulness of GAIT in a wide range of biomedical applications. Availability and implementation: http://cdal.korea.ac.kr/GAIT/index.html. Supplementary information: Supplementary data are available at Bioinformatics online.

Genomic responses in mouse models poorly mimic human inflammatory diseases.

A cornerstone of modern biomedical research is the use of mouse models to explore basic pathophysiological mechanisms, evaluate new therapeutic approaches, and make go or no-go decisions to carry new drug candidates forward into clinical trials. Systematic studies evaluating how well murine models mimic human inflammatory diseases are nonexistent. Here, we show that, although acute inflammatory stresses from different etiologies result in highly similar genomic responses in humans, the responses in corresponding mouse models correlate poorly with the human conditions and also, one another. Among genes changed significantly in humans, the murine orthologs are close to random in matching their human counterparts (e.g., R(2) between 0.0 and 0.1). In addition to improvements in the current animal model systems, our study supports higher priority for translational medical research to focus on the more complex human conditions rather than relying on mouse models to study human inflammatory diseases.

Density estimation on multivariate censored data with optional Pólya tree.

Analyzing the failure times of multiple events is of interest in many fields. Estimating the joint distribution of the failure times in a non-parametric way is not straightforward because some failure times are often right-censored and only known to be greater than observed follow-up times. Although it has been studied, there is no universally optimal solution for this problem. It is still challenging and important to provide alternatives that may be more suitable than existing ones in specific settings. Related problems of the existing methods are not only limited to infeasible computations, but also include the lack of optimality and possible non-monotonicity of the estimated survival function. In this paper, we proposed a non-parametric Bayesian approach for directly estimating the density function of multivariate survival times, where the prior is constructed based on the optional Pólya tree. We investigated several theoretical aspects of the procedure and derived an efficient iterative algorithm for implementing the Bayesian procedure. The empirical performance of the method was examined via extensive simulation studies. Finally, we presented a detailed analysis using the proposed method on the relationship among organ recovery times in severely injured patients. From the analysis, we suggested interesting medical information that can be further pursued in clinics.

Human transcriptome array for high-throughput clinical studies.

A 6.9 million-feature oligonucleotide array of the human transcriptome [Glue Grant human transcriptome (GG-H array)] has been developed for high-throughput and cost-effective analyses in clinical studies. This array allows comprehensive examination of gene expression and genome-wide identification of alternative splicing as well as detection of coding SNPs and noncoding transcripts. The performance of the array was examined and compared with mRNA sequencing (RNA-Seq) results over multiple independent replicates of liver and muscle samples. Compared with RNA-Seq of 46 million uniquely mappable reads per replicate, the GG-H array is highly reproducible in estimating gene and exon abundance. Although both platforms detect similar expression changes at the gene level, the GG-H array is more sensitive at the exon level. Deeper sequencing is required to adequately cover low-abundance transcripts. The array has been implemented in a multicenter clinical program and has generated high-quality, reproducible data. Considering the clinical trial requirements of cost, sample availability, and throughput, the GG-H array has a wide range of applications. An emerging approach for large-scale clinical genomic studies is to first use RNA-Seq to the sufficient depth for the discovery of transcriptome elements relevant to the disease process followed by high-throughput and reliable screening of these elements on thousands of patient samples using custom-designed arrays.

Prioritizing Disease Diagnosis in Neonatal Cohorts through Multivariate Survival Analysis: A Nonparametric Bayesian Approach.

Understanding the intricate relationships between diseases is critical for both prevention and recovery. However, there is a lack of suitable methodologies for exploring the precedence relationships within multiple censored time-to-event data, resulting in decreased analytical accuracy. This study introduces the Censored Event Precedence Analysis (CEPA), which is a nonparametric Bayesian approach suitable for understanding the precedence relationships in censored multivariate events. CEPA aims to analyze the precedence relationships between events to predict subsequent occurrences effectively. We applied CEPA to neonatal data from the National Health Insurance Service, identifying the precedence relationships among the seven most commonly diagnosed diseases categorized by the International Classification of Diseases. This analysis revealed a typical diagnostic sequence, starting with respiratory diseases, followed by skin, infectious, digestive, ear, eye, and injury-related diseases. Furthermore, simulation studies were conducted to demonstrate CEPA suitability for censored multivariate datasets compared to traditional models. The performance accuracy reached 76% for uniform distribution and 65% for exponential distribution, showing superior performance in all four tested environments. Therefore, the statistical approach based on CEPA enhances our understanding of disease interrelationships beyond competitive methodologies. By identifying disease precedence with CEPA, we can preempt subsequent disease occurrences and propose a healthcare system based on these relationships.

Automated BigSMILES conversion workflow and dataset for homopolymeric macromolecules.

The simplified molecular-input line-entry system (SMILES) has been utilized in a variety of artificial intelligence analyses owing to its capability of representing chemical structures using line notation. However, its ease of representation is limited, which has led to the proposal of BigSMILES as an alternative method suitable for the representation of macromolecules. Nevertheless, research on BigSMILES remains limited due to its preprocessing requirements. Thus, this study proposes a conversion workflow of BigSMILES, focusing on its automated generation from SMILES representations of homopolymers. BigSMILES representations for 4,927,181 records are provided, thereby enabling its immediate use for various research and development applications. Our study presents detailed descriptions on a validation process to ensure the accuracy, interchangeability, and robustness of the conversion. Additionally, a systematic overview of utilized codes and functions that emphasizes their relevance in the context of BigSMILES generation are produced. This advancement is anticipated to significantly aid researchers and facilitate further studies in BigSMILES representation, including potential applications in deep learning and further extension to complex structures such as copolymers.

Development of a genomic metric that can be rapidly used to predict clinical outcome in severely injured trauma patients.

OBJECTIVE: Many patients have complicated recoveries following severe trauma due to the development of organ injury. Physiological and anatomical prognosticators have had limited success in predicting clinical trajectories. We report on the development and retrospective validation of a simple genomic composite score that can be rapidly used to predict clinical outcomes. DESIGN: Retrospective cohort study. SETTING: Multi-institutional level 1 trauma centers. PATIENTS: Data were collected from 167 severely traumatized (injury severity score >15) adult (18-55 yr) patients. METHODS: Microarray-derived genomic data obtained from 167 severely traumatized patients over 28 days were assessed for differences in messenger RNA abundance among individuals with different clinical trajectories. Once a set of genes was identified based on differences in expression over the entire study period, messenger RNA abundance from these subjects obtained in the first 24 hours was analyzed in a blinded fashion using a rapid multiplex platform, and genomic data reduced to a single metric. RESULTS: From the existing genomic dataset, we identified 63 genes whose leukocyte expression differed between an uncomplicated and complicated clinical outcome over 28 days. Using a multiplex approach that can quantitate messenger RNA abundance in less than 12 hours, we reassessed total messenger RNA abundance from the first 24 hours after trauma and reduced the genomic data to a single composite score using the difference from reference. This composite score showed good discriminatory capacity to distinguish patients with a complicated outcome (area under a receiver-operator curve, 0.811; p <0.001). This was significantly better than the predictive power of either Acute Physiology and Chronic Health Evaluation II or new injury severity score scoring systems. CONCLUSIONS: A rapid genomic composite score obtained in the first 24 hours after trauma can retrospectively identify trauma patients who are likely to develop complicated clinical trajectories. A novel platform is described in which this genomic score can be obtained within 12 hours of blood collection, making it available for clinical decision making.

JETTA: junction and exon toolkits for transcriptome analysis.

SUMMARY: High-throughput genome-wide studies of alternatively spliced mRNA transcripts have become increasingly important in clinical research. Consequently, easy-to-use software tools are required to process data from these studies, for example, using exon and junction arrays. Here, we introduce JETTA, an integrated software package for the calculation of gene expression indices as well as the identification and visualization of alternative splicing events. We demonstrate the software using data of human liver and muscle samples hybridized on an exon-junction array. AVAILABILITY: JETTA and its demonstrations are freely available at http://igenomed.stanford.edu/~junhee/JETTA/index.html

Grid-wise simulation acceleration of the electromagnetic fields of 2D optical devices using super-resolution.

The significance of simulation has been increasing in device design due to the cost of real test. The accuracy of the simulation increases as the resolution of the simulation increases. However, the high-resolution simulation is not suited for actual device design because the amount of computing exponentially increases as the resolution increases. In this study, we introduce a model that predicts high-resolution outcomes using low-resolution calculated values which successfully achieves high simulation accuracy with low computational cost. The fast residual learning super-resolution (FRSR) convolutional network model is a model that we introduced that can simulate electromagnetic fields of optical. Our model achieved high accuracy when using the super-resolution technique on a 2D slit array under specific circumstances and achieved an approximately 18 times faster execution time than the simulator. To reduce the model training time and enhance performance, the proposed model shows the best accuracy (R2: 0.9941) by restoring high-resolution images using residual learning and a post-upsampling method to reduce computation. It has the shortest training time among the models that use super-resolution (7000 s). This model addresses the issue of temporal limitations of high-resolution simulations of device module characteristics.

CTIVA: Censored time interval variable analysis.

Traditionally, datasets with multiple censored time-to-events have not been utilized in multivariate analysis because of their high level of complexity. In this paper, we propose the Censored Time Interval Analysis (CTIVA) method to address this issue. It estimates the joint probability distribution of actual event times in the censored dataset by implementing a statistical probability density estimation technique on the dataset. Based on the acquired event time, CTIVA investigates variables correlated with the interval time of events via statistical tests. The proposed method handles both categorical and continuous variables simultaneously-thus, it is suitable for application on real-world censored time-to-event datasets, which include both categorical and continuous variables. CTIVA outperforms traditional censored time-to-event data handling methods by 5% on simulation data. The average area under the curve (AUC) of the proposed method on the simulation dataset exceeds 0.9 under various conditions. Further, CTIVA yields novel results on National Sample Cohort Demo (NSCD) and proteasome inhibitor bortezomib dataset, a real-world censored time-to-event dataset of medical history of beneficiaries provided by the National Health Insurance Sharing Service (NHISS) and National Center for Biotechnology Information (NCBI). We believe that the development of CTIVA is a milestone in the investigation of variables correlated with interval time of events in presence of censoring.

Clinical significance of urinary exosomal microRNAs in patients with IgA nephropathy.

Immunoglobulin A nephropathy (IgAN) is the most common primary glomerulonephritis worldwide. The clinical relevance of 11 urinary exosomal microRNAs (miRNAs) was evaluated in patients with IgAN. From January 2009 to November 2018, IgAN (n = 93), disease control (n = 11), and normal control (n = 19) groups were enrolled. We evaluated the expression levels of urinary exosomal miRNAs at the baseline and their relationship with clinical and pathologic features. This study aimed to discriminate statistically powerful urinary exosomal miRNAs for the prognosis of IgAN. Urinary miRNA levels of miR-16-5p, miR-29a-3p, miR-124-3p, miR-126-3p, miR-199a-3p, miR-199b-5p, and miR-335-3p showed significant correlation with both estimated glomerular filtration rate (eGFR) and urine protein-to-creatinine ratio (uPCR). In univariate regression analysis, age, body mass index, hypertension, eGFR, uPCR, Oxford classification E, and three miRNAs (miR-16-5p, miR-199a-3p, and miR-335-3p) were associated with disease progression in patients with IgAN. The area under the curve (AUC) of miR-199a-3p was high enough (0.749) without any other clinical or pathologic factors, considering that the AUC of the International IgAN Risk Prediction Tool was 0.853. Urinary exosomal miRNAs may serve as alternative prognostic biomarkers of IgAN with further research.

Corrigendum: Development and validation of urinary exosomal microRNA biomarkers for the diagnosis of acute rejection in kidney transplant recipients.

[This corrects the article DOI: 10.3389/fimmu.2023.1190576.].

Development and validation of urinary exosomal microRNA biomarkers for the diagnosis of acute rejection in kidney transplant recipients.

Introduction: Acute rejection (AR) continues to be a significant obstacle for short- and long-term graft survival in kidney transplant recipients. Herein, we aimed to examine urinary exosomal microRNAs with the objective of identifying novel biomarkers of AR. Materials and methods: Candidate microRNAs were selected using NanoString-based urinary exosomal microRNA profiling, meta-analysis of web-based, public microRNA database, and literature review. The expression levels of these selected microRNAs were measured in the urinary exosomes of 108 recipients of the discovery cohort using quantitative real-time polymerase chain reaction (qPCR). Based on the differential microRNA expressions, AR signatures were generated, and their diagnostic powers were determined by assessing the urinary exosomes of 260 recipients in an independent validation cohort. Results: We identified 29 urinary exosomal microRNAs as candidate biomarkers of AR, of which 7 microRNAs were differentially expressed in recipients with AR, as confirmed by qPCR analysis. A three-microRNA AR signature, composed of hsa-miR-21-5p, hsa-miR-31-5p, and hsa-miR-4532, could discriminate recipients with AR from those maintaining stable graft function (area under the curve [AUC] = 0.85). This signature exhibited a fair discriminative power in the identification of AR in the validation cohort (AUC = 0.77). Conclusion: We have successfully demonstrated that urinary exosomal microRNA signatures may form potential biomarkers for the diagnosis of AR in kidney transplantation recipients.

Conversion of Automated 12-Lead Electrocardiogram Interpretations to OMOP CDM Vocabulary.

BACKGROUND: A computerized 12-lead electrocardiogram (ECG) can automatically generate diagnostic statements, which are helpful for clinical purposes. Standardization is required for big data analysis when using ECG data generated by different interpretation algorithms. The common data model (CDM) is a standard schema designed to overcome heterogeneity between medical data. Diagnostic statements usually contain multiple CDM concepts and also include non-essential noise information, which should be removed during CDM conversion. Existing CDM conversion tools have several limitations, such as the requirement for manual validation, inability to extract multiple CDM concepts, and inadequate noise removal. OBJECTIVES: We aim to develop a fully automated text data conversion algorithm that overcomes limitations of existing tools and manual conversion. METHODS: We used interpretations printed by 12-lead resting ECG tests from three different vendors: GE Medical Systems, Philips Medical Systems, and Nihon Kohden. For automatic mapping, we first constructed an ontology-lexicon of ECG interpretations. After clinical coding, an optimized tool for converting ECG interpretation to CDM terminology is developed using term-based text processing. RESULTS: Using the ontology-lexicon, the cosine similarity-based algorithm and rule-based hierarchical algorithm showed comparable conversion accuracy (97.8 and 99.6%, respectively), while an integrated algorithm based on a heuristic approach, ECG2CDM, demonstrated superior performance (99.9%) for datasets from three major vendors. CONCLUSION: We developed a user-friendly software that runs the ECG2CDM algorithm that is easy to use even if the user is not familiar with CDM or medical terminology. We propose that automated algorithms can be helpful for further big data analysis with an integrated and standardized ECG dataset.

Simulator acceleration and inverse design of fin field-effect transistors using machine learning.

The simulation and design of electronic devices such as transistors is vital for the semiconductor industry. Conventionally, a device is intuitively designed and simulated using model equations, which is a time-consuming and expensive process. However, recent machine learning approaches provide an unprecedented opportunity to improve these tasks by training the underlying relationships between the device design and the specifications derived from the extensively accumulated simulation data. This study implements various machine learning approaches for the simulation acceleration and inverse-design problems of fin field-effect transistors. In comparison to traditional simulators, the proposed neural network model demonstrated almost equivalent results (R2 = 0.99) and was more than 122,000 times faster in simulation. Moreover, the proposed inverse-design model successfully generated design parameters that satisfied the desired target specifications with high accuracies (R2 = 0.96). Overall, the results demonstrated that the proposed machine learning models aided in achieving efficient solutions for the simulation and design problems pertaining to electronic devices. Thus, the proposed approach can be further extended to more complex devices and other vital processes in the semiconductor industry.

Using high-density exon arrays to profile gene expression in closely related species.

Global comparisons of gene expression profiles between species provide significant insight into gene regulation, evolutionary processes and disease mechanisms. In this work, we describe a flexible and intuitive approach for global expression profiling of closely related species, using high-density exon arrays designed for a single reference genome. The high-density probe coverage of exon arrays allows us to select identical sets of perfect-match probes to measure expression levels of orthologous genes. This eliminates a serious confounding factor in probe affinity effects of species-specific microarray probes, and enables direct comparisons of estimated expression indexes across species. Using a newly designed Affymetrix exon array, with eight probes per exon for approximately 315,000 exons in the human genome, we conducted expression profiling in corresponding tissues from humans, chimpanzees and rhesus macaques. Quantitative real-time PCR analysis of differentially expressed candidate genes is highly concordant with microarray data, yielding a validation rate of 21/22 for human versus chimpanzee differences, and 11/11 for human versus rhesus differences. This method has the potential to greatly facilitate biomedical and evolutionary studies of gene expression in nonhuman primates and can be easily extended to expression array design and comparative analysis of other animals and plants.