Machine-learned analysis of the association of next-generation sequencing-based human TRPV1 and TRPA1 genotypes with the sensitivity to heat stimuli and topically applied capsaicin

Artificial intelligence and database for NGS-based diagnosis in rare disease

Rare diseases (RDs) are rare complex genetic diseases affecting a conservative estimate of 300 million people worldwide. Recent Next-Generation Sequencing (NGS) studies are unraveling the underlying genetic heterogeneity of this group of diseases. NGS-based methods used in RDs studies have improved the diagnosis and management of RDs. Concomitantly, a suite of bioinformatics tools has been developed to sort through big data generated by NGS to understand RDs better. However, there are concerns regarding the lack of consistency among different methods, primarily linked to factors such as the lack of uniformity in input and output formats, the absence of a standardized measure for predictive accuracy, and the regularity of updates to the annotation database. Today, artificial intelligence (AI), particularly deep learning, is widely used in a variety of biological contexts, changing the healthcare system. AI has demonstrated promising capabilities in boosting variant calling precision, refining variant prediction, and enhancing the user-friendliness of electronic health record (EHR) systems in NGS-based diagnostics. This paper reviews the state of the art of AI in NGS-based genetics, and its future directions and challenges. It also compare several rare disease databases.

Pain modalities in the body and brain: Current knowledge and future perspectives

Development and validation of pain biomarkers has become a major issue in pain research. Recent advances in multimodal data acquisition have allowed researchers to gather multivariate and multilevel whole-body measurements in patients with pain conditions, and data analysis techniques such as machine learning have led to novel findings in neural biomarkers for pain. Most studies have focused on the development of a biomarker to predict the severity of pain with high precision and high specificity, however, a similar approach to discriminate different modalities of pain is lacking. Identification of more accurate and specific pain biomarkers will require an in-depth understanding of the modality specificity of pain. In this review, we summarize early and recent findings on the modality specificity of pain in the brain, with a focus on distinct neural activity patterns between chronic clinical and acute experimental pain, direct, social, and vicarious pain, and somatic and visceral pain. We also suggest future directions to improve our current strategy of pain management using our knowledge of modality-specific aspects of pain.

Unearthing of Key Genes Driving the Pathogenesis of Alzheimer's Disease via Bioinformatics

Alzheimer’s disease (AD) is a neurodegenerative disease with unelucidated molecular pathogenesis. Herein, we aimed to identify potential hub genes governing the pathogenesis of AD. The AD datasets of GSE118553 and GSE131617 were collected from the NCBI GEO database. The weighted gene coexpression network analysis (WGCNA), differential gene expression analysis, and functional enrichment analysis were performed to reveal the hub genes and verify their role in AD. Hub genes were validated by machine learning algorithms. We identified modules and their corresponding hub genes from the temporal cortex (TC), frontal cortex (FC), entorhinal cortex (EC), and cerebellum (CE). We obtained 33, 42, 42, and 41 hub genes in modules associated with AD in TC, FC, EC, and CE tissues, respectively. Significant differences were recorded in the expression levels of hub genes between AD and the control group in the TC and EC tissues (P < 0.05). The differences in the expressions of FCGRT, SLC1A3, PTN, PTPRZ1, and PON2 in the FC and CE tissues among the AD and control groups were significant (P < 0.05). The expression levels of PLXNB1, GRAMD3, and GJA1 were statistically significant between the Braak NFT stages of AD. Overall, our study uncovered genes that may be involved in AD pathogenesis and revealed their potential for the development of AD biomarkers and appropriate AD therapeutics targets.

Computational Functional Genomics-Based AmpliSeq™ Panel for Next-Generation Sequencing of Key Genes of Pain

The genetic background of pain is becoming increasingly well understood, which opens up possibilities for predicting the individual risk of persistent pain and the use of tailored therapies adapted to the variant pattern of the patient's pain-relevant genes. The individual variant pattern of pain-relevant genes is accessible via next-generation sequencing, although the analysis of all "pain genes" would be expensive. Here, we report on the development of a cost-effective next generation sequencing-based pain-genotyping assay comprising the development of a customized AmpliSeq™ panel and bioinformatics approaches that condensate the genetic information of pain by identifying the most representative genes. The panel includes 29 key genes that have been shown to cover 70% of the biological functions exerted by a list of 540 so-called "pain genes" derived from transgenic mice experiments. These were supplemented by 43 additional genes that had been independently proposed as relevant for persistent pain. The functional genomics covered by the resulting 72 genes is particularly represented by mitogen-activated protein kinase of extracellular signal-regulated kinase and cytokine production and secretion. The present genotyping assay was established in 61 subjects of Caucasian ethnicity and investigates the functional role of the selected genes in the context of the known genetic architecture of pain without seeking functional associations for pain. The assay identified a total of 691 genetic variants, of which many have reports for a clinical relevance for pain or in another context. The assay is applicable for small to large-scale experimental setups at contemporary genotyping costs.

Epidermal expression of human TRPM8, but not of TRPA1 ion channels, is associated with sensory responses to local skin cooling

Human cold perception and nociception play an important role in persisting pain. However, species differences in the target temperature of thermosensitive ion channels expressed in peripheral nerve endings have fueled discussions about the mechanism of cold nociception in humans. Most frequently implicated thermosensors are members of the transient receptor potential (TRP) ion channel family TRPM8 and TRPA1. Regularly observed, distinct cold pain phenotype groups suggested the existence of interindividually differing molecular bases. In 28 subjects displaying either high or medium sensitivity to local cooling of the skin, the density at epidermal nerve fibers of TRPM8, but not that of TRPA1 expression, correlated significantly with the cold pain threshold. Moreover, reproducible grouping of the subjects, based on high or medium sensitivity to cooling, was reflected in an analogous grouping based on high or low TRPM8 expression at epidermal nerve fibers. The distribution of TRPM8 expression in epidermal nerve fibers provided an explanation for the previously observed (bi)modal distribution of human cold pain thresholds which was reproduced in this study. In the light of current controversies on the role of human TRPA1 ion channels in cold pain perception, the present observations demonstrating a lack of association of TRPA1 channel expression with cold sensitivity-related measures reinforce doubts about involvement of this channel in cold pain in humans. Since TRP inhibitors targeting TRPM8 and TRPA1 are currently entering clinical phases of drug development, the existence of known species differences, in particular in the function of TRPA1, emphasizes the increasing importance of new methods to directly approach the roles of TRPs in humans.

Machine-learned analysis of the association of next-generation sequencing-based genotypes with persistent pain after breast cancer surgery

Cancer and its surgical treatment are among the most important triggering events for persistent pain, but additional factors need to be present for the clinical manifestation, such as variants in pain-relevant genes. In a cohort of 140 women undergoing breast cancer surgery, assigned based on a 3-year follow-up to either a persistent or nonpersistent pain phenotype, next-generation sequencing was performed for 77 genes selected for known functional involvement in persistent pain. Applying machine-learning and item categorization techniques, 21 variants in 13 different genes were found to be relevant to the assignment of a patient to either the persistent pain or the nonpersistent pain phenotype group. In descending order of importance for correct group assignment, the relevant genes comprised DRD1, FAAH, GCH1, GPR132, OPRM1, DRD3, RELN, GABRA5, NF1, COMT, TRPA1, ABHD6, and DRD4, of which one in the DRD4 gene was a novel discovery. Particularly relevant variants were found in the DRD1 and GPR132 genes, or in a cis-eCTL position of the OPRM1 gene. Supervised machine-learning-based classifiers, trained with 2/3 of the data, identified the correct pain phenotype group in the remaining 1/3 of the patients at accuracies and areas under the receiver operator characteristic curves of 65% to 72%. When using conservative classical statistical approaches, none of the variants passed α-corrected testing. The present data analysis approach, using machine learning and training artificial intelligences, provided biologically plausible results and outperformed classical approaches to genotype-phenotype association.

Machine-Learned Association of Next-Generation Sequencing-Derived Variants in Thermosensitive Ion Channels Genes with Human Thermal Pain Sensitivity Phenotypes

Genetic association studies have shown their usefulness in assessing the role of ion channels in human thermal pain perception. We used machine learning to construct a complex phenotype from pain thresholds to thermal stimuli and associate it with the genetic information derived from the next-generation sequencing (NGS) of 15 ion channel genes which are involved in thermal perception, including ASIC1, ASIC2, ASIC3, ASIC4, TRPA1, TRPC1, TRPM2, TRPM3, TRPM4, TRPM5, TRPM8, TRPV1, TRPV2, TRPV3, and TRPV4. Phenotypic information was complete in 82 subjects and NGS genotypes were available in 67 subjects. A network of artificial neurons, implemented as emergent self-organizing maps, discovered two clusters characterized by high or low pain thresholds for heat and cold pain. A total of 1071 variants were discovered in the 15 ion channel genes. After feature selection, 80 genetic variants were retained for an association analysis based on machine learning. The measured performance of machine learning-mediated phenotype assignment based on this genetic information resulted in an area under the receiver operating characteristic curve of 77.2%, justifying a phenotype classification based on the genetic information. A further item categorization finally resulted in 38 genetic variants that contributed most to the phenotype assignment. Most of them (10) belonged to the TRPV3 gene, followed by TRPM3 (6). Therefore, the analysis successfully identified the particular importance of TRPV3 and TRPM3 for an average pain phenotype defined by the sensitivity to moderate thermal stimuli.

Computational functional genomics-based reduction of disease-related gene sets to their key components

MOTIVATION

The genetic architecture of diseases becomes increasingly known. This raises difficulties in picking suitable targets for further research among an increasing number of candidates. Although expression based methods of gene set reduction are applied to laboratory-derived genetic data, the analysis of topical sets of genes gathered from knowledge bases requires a modified approach as no quantitative information about gene expression is available.

RESULTS

We propose a computational functional genomics-based approach at reducing sets of genes to the most relevant items based on the importance of the gene within the polyhierarchy of biological processes characterizing the disease. Knowledge bases about the biological roles of genes can provide a valid description of traits or diseases represented as a directed acyclic graph (DAG) picturing the polyhierarchy of disease relevant biological processes. The proposed method uses a gene importance score derived from the location of the gene-related biological processes in the DAG. It attempts to recreate the DAG and thereby, the roles of the original gene set, with the least number of genes in descending order of importance. This obtained precision and recall of over 70% to recreate the components of the DAG charactering the biological functions of n=540 genes relevant to pain with a subset of only the k=29 best-scoring genes.

CONCLUSIONS

A new method for reduction of gene sets is shown that is able to reproduce the biological processes in which the full gene set is involved by over 70%; however, by using only ∼5% of the original genes.

AVAILABILITY AND IMPLEMENTATION

The necessary numerical parameters for the calculation of gene importance are implemented in the R package dbtORA at https://github.com/IME-TMP-FFM/dbtORA.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Machine-learned analysis of global and glial/opioid intersection-related DNA methylation in patients with persistent pain after breast cancer surgery

BACKGROUND

Glial cells in the central nervous system play a key role in neuroinflammation and subsequent central sensitization to pain. They are therefore involved in the development of persistent pain. One of the main sites of interaction of the immune system with persistent pain has been identified as neuro-immune crosstalk at the glial-opioid interface. The present study examined a potential association between the DNA methylation of two key players of glial/opioid intersection and persistent postoperative pain.

METHODS

In a cohort of 140 women who had undergone breast cancer surgery, and were assigned based on a 3-year follow-up to either a persistent or non-persistent pain phenotype, the role of epigenetic regulation of key players in the glial-opioid interface was assessed. The methylation of genes coding for the Toll-like receptor 4 (TLR4) as a major mediator of glial contributions to persistent pain or for the μ-opioid receptor (OPRM1) was analyzed and its association with the pain phenotype was compared with that conferred by global genome-wide DNA methylation assessed via quantification of the methylation in the retrotransposon LINE1.

RESULTS

Training of machine learning algorithms indicated that the global DNA methylation provided a similar diagnostic accuracy for persistent pain as previously established non-genetic predictors. However, the diagnosis can be based on a single DNA based marker. By contrast, the methylation of TLR4 or OPRM1 genes could not contribute further to the allocation of the patients to the pain-related phenotype groups.

CONCLUSIONS

While clearly supporting a predictive utility of epigenetic testing, the present analysis cannot provide support for specific epigenetic modulation of persistent postoperative pain via methylation of two key genes of the glial-opioid interface.

[Generating knowledge from complex data sets in human experimental pain research]

Pain has a complex pathophysiology that is expressed in multifaceted and heterogeneous clinical phenotypes. This makes research on pain and its treatment a potentially data-rich field as large amounts of complex data are generated. Typical sources of such data are investigations with functional magnetic resonance imaging, complex quantitative sensory testing, next-generation DNA sequencing and functional genomic research approaches, such as those aimed at analgesic drug discovery or repositioning of drugs known from other indications as new analgesics. Extracting information from these big data requires complex data scientific-based methods belonging more to computer science than to statistics. A particular interest is currently focused on machine learning, the methods of which are used for the detection of interesting and biologically meaningful structures in high-dimensional data. Subsequently, classifiers can be created that predict clinical phenotypes from, e.g. clinical or genetic features acquired from subjects. In addition, knowledge discovery in big data accessible in electronic knowledge bases, can be used to generate hypotheses and to exploit the accumulated knowledge about pain for the discovery of new analgesic drugs. This enables so-called data-information-knowledge-wisdom (DIKW) approaches to be followed in pain research. This article highlights current examples from pain research to provide an overview about contemporary data scientific methods used in this field of research.

Genetic study in patients operated dentally and anesthetized with articaine-epinephrine

Aims: In this study we wanted to figure out if there was a correlation between OPRM1 N40D, TRPV1 I316M, TRPV1 I585V, NOS3 −786T>C and IL6 −174C>G polymorphisms and the response to locally applied articaine-epinephrine anesthetic.

Methods: In this observational study, 114 oral cell samples of patients anesthetized with articaine-epinephrine (54 from men 60 from women), were collected from dental centers in Madrid (Spain). High molecular weight DNA was obtained from oral mucosa cells. The analysis of OPRM1 N40D (rs1799971), TRPV1 I316M (rs222747), TRPV1 I585V (rs8065080) and IL6 −174C>G polymorphism was performed through real-time PCR allelic discrimination using TaqMan probes. Polymorphism NOS3 −786T> C (rs2070744) was analyzed using RFLP-PCR.

Results: The studied polymorphisms are involved neither in the response to the anesthetic, nor in the intensity of perceived dental pain. However, in a subset of female patients we found that TRPV1 I316M was associated with a delayed onset of anesthesia.

Conclusions: There is no association among these polymorphisms and the time elapsed between the application of the anesthetic and the onset of its effect.

Machine Learning in Human Olfactory Research

The complexity of the human sense of smell is increasingly reflected in complex and high-dimensional data, which opens opportunities for data-driven approaches that complement hypothesis-driven research. Contemporary developments in computational and data science, with its currently most popular implementation as machine learning, facilitate complex data-driven research approaches. The use of machine learning in human olfactory research included major approaches comprising 1) the study of the physiology of pattern-based odor detection and recognition processes, 2) pattern recognition in olfactory phenotypes, 3) the development of complex disease biomarkers including olfactory features, 4) odor prediction from physico-chemical properties of volatile molecules, and 5) knowledge discovery in publicly available big databases. A limited set of unsupervised and supervised machine-learned methods has been used in these projects, however, the increasing use of contemporary methods of computational science is reflected in a growing number of reports employing machine learning for human olfactory research. This review provides key concepts of machine learning and summarizes current applications on human olfactory data.

Human pharmacological approaches to TRP-ion-channel-based analgesic drug development

The discovery of novel analgesic drug targets is an active research topic owing to insufficient treatment options for persisting pain. Modulators of temperature-sensing transient receptor potential ion channels (thermoTRPs), in particular TRPV1, TRPV2, TRPM8 and TRPA1, have reached clinical development. This requires access for TRP channels and the effects of specific modulators in humans. This is currently possible via (i) the study of TRP channel function in human-derived cell lines, (ii) immunohistochemical visualization of TRP channel expression in human tissues, (iii) human experimental pain models employing sensitization by means of topical application of TRP channel activators including capsaicin (TRPV1), menthol (TRPM8), mustard oil and cinnamaldehyde (TRPA1), and (iv) the study of phenotypic consequences of human TRP gene variants.