Pharmacogenetics (PGx) aims to identify the genetic factors that determine inter-individual differences in response to drug treatment maximizing efficacy while decreasing the risk of adverse events. Estimating the prevalence of PGx variants involved in drug response, is a critical preparatory step for large-scale implementation of a personalized medicine program in a target population. Here, we profiled pharmacogenetic variation in fourteen clinically relevant genes in a representative sample set of 1,577 unrelated sequenced Sardinians -- an ancient island population that well accounts for genetic variation in Europe as a whole and, at the same time are enriched in genetic variants that are very rare elsewhere. To this end, we used PGxPOP, a PGx allele caller based on the guidelines created by the Clinical Pharmacogenetics Implementation Consortium (CPIC), to identify the main phenotypes associated with the PGx alleles most represented in Sardinians. We estimate that 99.43% of Sardinian individuals may potentially respond atypically to at least one drug, that on average each individual is expected to have an abnormal response to about 17 drugs, and that for 27 drugs the fraction of the population at risk of atypical responses to therapy is more than 40%. Finally, we highlighted 174 pharmacogenetic variants for which the minor allele frequency is at least 10% higher among Sardinians as compared to other European populations, a fact that may contribute to substantial interpopulation variability in drug response phenotypes. This study provides baseline information for further large-scale pharmacogenomic investigations in the Sardinian population and underlines the importance of the PGx characterization of diverse than European populations as Sardinians.

Rheumatoid arthritis (RA); is a chronic systemic autoimmune disease which causes loss of joint function and significantly reduces quality of life. Plasma metabolite concentrations of disease-modifying anti-rheumatic drugs (DMARDs) can influence treatment efficacy and toxicity. This study explored the relationship between DMARD metabolising gene variants and plasma metabolite levels in RA patients. DMARD metabolite concentrations were determined by tandem mass-spectrometry in plasma samples from 100 RA patients with actively flaring disease, collected at two intervals. Taqman-probes were used to discriminate single nucleotide polymorphism (SNP) genotypes in cohort genomic DNA: rs246240 (ABCC1), rs1476413 (MTHFR), rs2231142 (ABCG2), rs3740065 (ABCC2), rs4149081 (SLCO1B1), rs4846051 (MTHFR), rs10280623 (ABCB1), rs16853826 (ATIC), rs17421511 (MTHFR) and rs717620 (ABCC2). Mean plasma concentrations of methotrexate (MTX) and MTX-7-OH metabolites were higher (p<0.05, C.I. 95%) at baseline in rs4149081 GA genotype patients. Patients with rs1476413 SNP TT or CT alleles in the have significantly higher (p<0.001, 95% C.I) plasma poly-glutamate metabolites at both study time points and correspondingly elevated disease activity scores. Patients with the rs17421511 SNP AA allele reported significantly lower pain scores (p<0.05, 95% C.I.) at both study intervals. Genotyping strategies could help prioritise treatments to RA patients most likely to gain clinical benefit, whilst minimizing toxicity.

Personalized medicine is an aspect of the P4 medicine (predictive, preventive, personalized and participatory) based precisely on the customization of all medical characters of each subject. In personalized medicine, the development of medical treatments and drugs is tailored to the individual characteristics and needs of each subject, according to the study of diseases at different scales from genotype to phenotype scale. To make concrete the goal of personalized medicine, it is necessary to employ high-throughput methodologies such as Next Generation Sequencing (NGS), Genome-Wide Association Studies (GWAS), Mass Spectrometry or Microarrays, that are able to investigate a single disease from a broader perspective. For example, by using genotyping microarrays (e.g. collections of Single Nucleotide Polymorphism - SNP) it is possible to uncover the reasons (i.e. mutation in genes) because a treatment works properly in some patients (for example absence of mutated genes), but it does not work (presence of mutated genes) in others. A side effect of high-throughput methodologies is the massive amount of data produced for each single experiment, that poses several challenges (e.g. high execution time and required memory) to bioinformatic software. Thus a main requirement of modern bioinformatic software is the use of good software engineering methods and efficient programming techniques, able to face those challenges, that include the use of parallel programming and efficient and compact data structures. Thus, to exploit all the potential of this massive amount of data in the short possible time (before that data becomes obsolete), the necessity to develop parallel software tools for efficient data collection and analysis arise. Moreover, due to the heterogeneity of the data produced by the different kinds of experimental platforms, it is necessary to automatize in a comprehensive software pipeline, the various steps that compose a bioinformatic analysis, such as: the preprocessing of raw data to remove noise or corrupted data; the annotation of data with external knowledge (e.g. Gene Ontology), and the integration of molecular data with clinical data. It should be noted that such steps are necessary to make statistical or data mining analysis more effective. This paper presents the design and the experimentation of a comprehensive software pipeline, named microPipe, for the preprocessing, annotation and analysis of microarray-based SNP genotyping data. A case study in pharmacogenomics is presented. The main advantages of using microPipe are: the reduction of errors that may happen when trying to make data compatible among different tools; the possibility to analyze in parallel huge datasets; the easy annotation and integration of data.

Genome sequencing is enabling precision medicine—tailoring treatment to the unique constellation of variants in an individual’s genome. The impact of recurrent pathogenic variants is often understood, leaving a long tail of rare genetic variants that are uncharacterized. The problem of uncharacterized rare variation is especially acute when it occurs in genes of known clinical importance with functionally consequent frequent variants and associated mechanisms. Variants of unknown significance (VUS) in these genes are discovered at a rate that outpaces current ability to classify them using databases of previous cases, experimental evaluation, and computational predictors. Clinicians are thus left without guidance about the significance of variants that may have actionable consequences. Computational prediction of the impact of rare genetic variation is increasingly becoming an important capability. In this paper, we review the technical and ethical challenges of interpreting the function of rare variants in two settings: inborn errors of metabolism in newborns, and pharmacogenomics. We propose a framework for a genomic learning healthcare system with an initial focus on early-onset treatable disease in newborns and actionable pharmacogenomics. We argue that (1) a genomic learning healthcare system must allow for continuous collection and assessment of rare variants, (2) emerging machine learning methods will enable algorithms to predict the clinical impact of rare variants on protein function, and (3) ethical considerations must inform the construction and deployment of all rare-variation triage strategies, particularly with respect to health disparities arising from unbalanced ancestry representation.

Background: Accurate prediction of patients’ response to therapy is clinically indispensable, howbeit challenging. With increased understanding of the human genome and malignancies, there is the renaissance of in silico pharmacogenomics with renewed interest in drug response predictability based on gene-drug interaction. Objective: Evidence-based transcript-proteome profiling is essential for synthesizing clinically applicable algorithms for predicting response to anticancer therapy, including immune checkpoint blockade (ICBT); thus, saving physicians’ time, reducing polypharmacy, and curtailing unnecessary treatment expense. In this study, we tested and validated the hypothesis that a selected proteomic signature in ICBT-naïve patients is sufficient for the prediction of response to ICBT. Methods: Using a multimodal approach consisting of computational pharmacogenomics, transcript-proteome analytics, mathematical modeling, and machine learning systems; we delineated therapy-sensitivity and stratified patients into graduated response groups based on their proteomic profile. Protein expression levels in our cohort tissue specimens were evaluated based on T cell- and non-T cell- inflamed phenotypes by immunohistochemistry. Results: We established β-catenin, PDL1, CD3 and CD8 expression-based ICBT response model. Statistical regression models validated the predictive association between our predefined algorithms and therapeutic outcome. Interestingly, our 4-gene prediction classifier was constitutively independent of tumor tissue origin, correctly stratified patients into high-, low-, and non- responders pre-treatment, with high prediction accuracy, and exhibited good association with patients’ performance status and prognosis (p < 0.01). Conclusion: Our findings demonstrate the possibility of accurate proteomics based ICBT response prediction and provide a putative basis for drug response prediction based on selective proteome profile in untreated cancer patients.

Fibromyalgia (FM) and Myalgic Encephalomyelitis/Chronic Fatigue Syndrome (ME/CFS) are diseases of unknown etiology presenting complex and often overlapping symptomatology. Despite promising advances on the study of miRNomes of these diseases, no validated molecular diagnostic biomarker yet exists. Since FM and ME/CFS patient treatments commonly include polypharmacy it is of concern that biomarker miRNAs are masked by drug interactions. Aiming at discriminating between drug-effects and true disease-associated differential miRNA expression, we evaluated the potential impact of commonly prescribed drugs on disease miRNomes, as reported by the literature. By using the web search tools SM2miR, Pharmaco-miR and repoDB, we found a list of commonly prescribed drugs that impact on FM and ME/CFS miRNomes and therefore could be interfering in the process of biomarker discovery. On another end, disease-associated miRNomes may incline patient´s response to treatment and toxicity. Here, we explored treatments for diseases in general that could be affected by FM and ME/CFS miRNomes finding a long list of them, including treatments for lymphoma, a type of cancer affecting ME/CFS patients at a higher rate than healthy population. We conclude that FM and ME/CFS miRNomes could help refine pharmacogenomic/pharmacoepigenomic analysis to elevate future personalized medicine and precision medicine programs in the clinic.