Genome sequence diversity of SARS-CoV-2 in Serbia: insights gained from a 3-year pandemic study.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), responsible for the COVID-19 pandemic, has been evolving rapidly causing emergence of new variants and health uncertainties. Monitoring the evolution of the virus was of the utmost importance for public health interventions and the development of national and global mitigation strategies. Here, we report national data on the emergence of new variants, their distribution, and dynamics in a 3-year study conducted from March 2020 to the end of January 2023 in the Republic of Serbia. Nasopharyngeal and oropharyngeal swabs from 2,398 COVID-19-positive patients were collected and sequenced using three different next generation technologies: Oxford Nanopore, Ion Torrent, and DNBSeq. In the subset of 2,107 SARS-CoV-2 sequences which met the quality requirements, detection of mutations, assignment to SARS-CoV-2 lineages, and phylogenetic analysis were performed. During the 3-year period, we detected three variants of concern, namely, Alpha (5.6%), Delta (7.4%), and Omicron (70.3%) and one variant of interest-Omicron recombinant "Kraken" (XBB1.5) (<1%), whereas 16.8% of the samples belonged to other SARS-CoV-2 (sub)lineages. The detected SARS-CoV-2 (sub)lineages resulted in eight COVID-19 pandemic waves in Serbia, which correspond to the pandemic waves reported in Europe and the United States. Wave dynamics in Serbia showed the most resemblance with the profile of pandemic waves in southern Europe, consistent with the southeastern European location of Serbia. The samples were assigned to sixteen SARS-CoV-2 Nextstrain clades: 20A, 20B, 20C, 20D, 20E, 20G, 20I, 21J, 21K, 21L, 22A, 22B, 22C, 22D, 22E, and 22F and six different Omicron recombinants (XZ, XAZ, XAS, XBB, XBF, and XBK). The 10 most common mutations detected in the coding and untranslated regions of the SARS-CoV-2 genomes included four mutations affecting the spike protein (S:D614G, S:T478K, S:P681H, and S:S477N) and one mutation at each of the following positions: 5'-untranslated region (5'UTR:241); N protein (N:RG203KR); NSP3 protein (NSP3:F106F); NSP4 protein (NSP4:T492I); NSP6 protein (NSP6: S106/G107/F108 - triple deletion), and NSP12b protein (NSP12b:P314L). This national-level study is the most comprehensive in terms of sequencing and genomic surveillance of SARS-CoV-2 during the pandemic in Serbia, highlighting the importance of establishing and maintaining good national practice for monitoring SARS-CoV-2 and other viruses circulating worldwide.

Can Pharmacogenetic Variants in TPMT, MTHFR and SLCO1B1 Genes Be Used as Potential Markers of Outcome Prediction in Systemic Sclerosis Patients?

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.

Genome-Wide Association Study of COVID-19 Outcomes Reveals Novel Host Genetic Risk Loci in the Serbian Population.

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-8). 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-6) and severe COVID-19 (p = 6.88 × 10-7), respectively. Suggestive association with COVID-19 pneumonia has also been observed at chromosomes 5p15.33 (IRX, NDUFS6, MRPL36, p = 2.81 × 10-6), 5q11.2 (ESM1, p = 6.59 × 10-6), and 9p23 (TYRP1, LURAP1L, p = 8.69 × 10-6). 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.

Pharmacogenomics landscape of COVID-19 therapy response in Serbian population and comparison with worldwide populations.

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 alternative 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.

Analysis of the promoter regions of disease-causing genes in maturity-onset diabetes of the young patients.

Maturity-onset diabetes of the young (MODY) is a form of monogenic diabetes caused by the variants in MODY-related genes. In addition to coding variants, variants in the promoter region of MODY-related genes can cause the disease as well. In this study, we screened the promoter regions of the most common MODY-related genes GCK, HNF1A, HNF4A and HNF1B in our cohort of 29 MODY patients. We identified one genetic variant in the HNF1A gene, a 7 bp insertion c.-154-160insTGGGGGT, and three variants in the GCK gene, -282C>T; -194A>G; 402C>G appearing as set. Chloramphenicol acetyltransferase (CAT) assay was performed to test the effect of the 7 bp insertion and the variant set on the activity of the reporter gene in HepG2 and RIN-5F cell, respectively, where a decreasing trend was observed for both variants. In silico analysis and electrophoretic mobility shift assay showed that the 7 bp insertion did not create the binding site for new transcriptional factors, but gave rise to additional binding sites for the existing ones. Results from our study indicated that the 7 bp insertion in the HNF1A gene could be associated with the patient's diabetes. As for the GCK variant set, it is probably not associated with diabetes in patients, but it may modify the fasting glucose level by causing small elevation in variant set carriers. We have presented two promoter variants in MODY-related genes. Variant in the HNF1A gene is presumed to be disease-causing and the GCK promoter variant set could be a phenotype modifier.

Pharmacogenomic Markers of Methotrexate Response in the Consolidation Phase of Pediatric Acute Lymphoblastic Leukemia Treatment.

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.