Methotrexate (MTX) is one of the staples of pediatric acute lymphoblastic leukemia (ALL) treatment. MTX targets the folate metabolic pathway (FMP). Abnormal function of the enzymes in FMP, due to genetic aberrations, leads to adverse drug reactions. The aim of this study was to investigate variants in pharmacogenes involved in FMP and their association with MTX pharmacokinetics (MTX elimination profile) and toxicity in the consolidation therapy phase of pediatric ALL patients. Eleven variants in the thymidylate synthetase (TYMS), methylenetetrahydrofolate reductase (MTHFR), dihydrofolate reductase (DHFR), SLC19A1 and SLCO1B genes were analyzed in 148 patients, using PCR- and sequencing-based methodology. For the Serbian and European control groups, data on allele frequency distribution were extracted from in-house and public databases. Our results show that the A allele of SLC19A1 c.80 variant contributes to slow MTX elimination. Additionally, the AA genotype of the same variant is a predictor of MTX-related hepatotoxicity. Patients homozygous for TYMS 6bp deletion were more likely to experience gastrointestinal toxicity. No allele frequency dissimilarity was found for the analyzed variants between Serbian and European populations. Statistical modelling did not show a joint effect of analyzed variants. Our results indicate that SLC19A1 c.80 variant and TYMS 6bp deletion are the most promising pharmacogenomic markers of MTX response in pediatric ALL patients.

Host genetics, an important contributor to the COVID-19 clinical susceptibility and severity, currently is the focus of multiple genome-wide association studies (GWAS) in populations affected by the pandemic. This is the first study from Serbia that performed a GWAS of COVID-19 outcomes to identify genetic risk markers of disease severity. A group of 128 hospitalized COVID-19 patients from the Serbian population was enrolled in the study. We conducted a GWAS comparing (1) patients with pneumonia (n = 80) against patients without pneumonia (n = 48), and (2) severe (n = 34) against mild disease (n = 48) patients, using a genotyping array followed by imputation of missing genotypes. We have detected a significant signal associated with COVID-19 related pneumonia at locus 13q21.33, with a peak residing upstream of the gene KLHL1 (p = 1.91 × 10⁻⁸). Our study also replicated a previously reported COVID-19 risk locus at 3p21.31, identifying lead variants in SACM1L and LZTFL1 genes suggestively associated with pneumonia (p = 7.54 × 10⁻⁶) and severe COVID-19 (p = 6.88 × 10⁻⁷), respectively. Suggestive association with COVID-19 pneumonia has also been observed at chromosomes 5p15.33 (IRX, NDUFS6, MRPL36, p = 2.81 × 10⁻⁶), 5q11.2 (ESM1, p = 6.59 × 10⁻⁶), and 9p23 (TYRP1, LURAP1L, p = 8.69 × 10⁻⁶). The genes located in or near the risk loci are expressed in neural or lung tissues, and have been previously associated with respiratory diseases such as asthma and COVID-19 or reported as differentially expressed in COVID-19 gene expression profiling studies. Our results revealed novel risk loci for pneumonia and severe COVID-19 disease which could contribute to a better understanding of the COVID-19 host genetics in different populations.

Systemic sclerosis (SSc) is a rare connective tissue disorder with highest morbidity and mortality among rheumatologic diseases. Disease progression is highly heterogeneous between patients, implying a strong need for individualization of therapy. Four pharmacogenetic variants, namely TPMT rs1800460, TPMT rs1142345, MTHFR rs1801133 and SLCO1B1 rs4149056 were tested for association with severe disease outcomes in 102 patients with SSc from Serbia treated either with immunosuppressants azathioprine (AZA) and methotrexate (MTX) or with other types of medications. Genotyping was performed using PCR-RFLP and direct Sanger sequencing. R software was used for statistical analysis and development of polygenic risk score (PRS) model. Association was found between MTHFR rs1801133 and higher risk for elevated systolic pressure in all patients except those prescribed with MTX, and higher risk for kidney insufficiency in patients prescribed with other types of drugs. In patients treated with MTX, variant SLCO1B1 rs4149056 was protective against kidney insufficiency. For patients receiving MTX a trend was shown for having a higher PRS rank and elevated systolic pressure. Our results open a door wide for more extensive research on pharmacogenomics markers in patients with SSc. Altogether, pharmacogenomics markers could predict the outcome of patients with SSc and help in prevention of adverse drug reactions.

Background: Since there are no certified therapeutics to treat COVID-19 patients, drug repurposing became important. With lack of time to test individual pharmacogenomics markers, population pharmacogenomics could be helpful in predicting a higher risk of developing adverse reactions and treatment failure in COVID-19 patients. Aim of our study was to identify pharmacogenes and pharmacogenomics markers associated with drugs recommended for COVID-19 treatment, chloroquine/hydroxychloroquine, azithromycin, lopinavir and ritonavir, in population of Serbia and other world populations. Methods: Genotype information of 143 individuals of Serbian origin was extracted from database previously obtained using TruSight One Gene Panel (Illumina). Genotype data of individuals from different world populations were extracted from the 1000 Genome Project. Fisher’s exact test was used for comparison of allele frequencies. Results: We have identified 11 potential pharmacogenomics markers in 7 pharmacogenes relevant for COVID-19 treatment. Based on high alterative allele frequencies in population and the functional effect of the variants, ABCB1 rs1045642 and rs2032582 could be relevant for reduced clearance of azithromycin, lopinavir and ritonavir drugs and UGT1A7 rs17868323 for hyperbilirubinemia in ritonavir treated COVID-19 patients in Serbian population. SLCO1B1 rs4149056 is a potential marker of lopinavir response, especially in Italian population. Our results confirmed that pharmacogenomics profile of African population is different from the rest of the world. Conclusions: Considering population specific pharmacogenomics landscape, preemptive testing for pharmacogenes relevant for drugs used in COVID-19 treatment could contribute to better understanding of the inconsistency in therapy response and could be applied to improve the outcome of the COVID-19 patients.

Background/Objectives: Vedolizumab (VDZ), a monoclonal antibody targeting α4β7 integrin, is used in Crohn’s disease (CD) management, yet patients’ responses vary, underscoring the need for pharmacogenomic (PGx) markers. This study aimed to identify PGx pathways associated with suboptimal VDZ response using a rare-variant analytical framework. Methods: DNA from 63 CD patients treated with VDZ as first-line advanced therapy underwent whole-exome sequencing. Clinical response at week 14 classified patients as optimal responders (ORs) or suboptimal responders (SRs). Sequencing data were processed using GATK Best Practices, annotated with variant effect predictors, and filtered for rare damaging variants (damaging missense and high-confidence loss-of-function; minor allele frequency < 0.05). Variants were mapped to genes specific for SRs and ORs, and analyzed for pathway enrichment using the Reactome database. Rare-variant burden and composition differences were assessed with Fisher’s exact test and SKAT-O gene-set association analysis. Results: Suboptimal VDZ response was associated with pathways related to membrane transport (ABC-family proteins, ion channels), L1–ankyrin interactions, and bile acid recycling, while optimal response was associated with pathways involving MET signaling. SKAT-O identified lipid metabolism-related pathways as significantly different—SRs harbored variants in pro-inflammatory lipid signaling and immune cell trafficking genes (e.g., PIK3CG, CYP4F2, PLA2R1), whereas ORs carried variants in fatty acid oxidation and detoxification genes (e.g., ACADM, CYP1A1, ALDH3A2, DECR1, MMUT). Conclusions: This study underscores the potential of exome-based rare-variant analysis to stratify CD patients and guide precision medicine approaches. The identified genes and pathways are potential PGx markers for CD patients treated with VDZ.