Tuberculosis in Ukrainian War Refugees and Migrants in the Czech Republic and Slovakia: A Molecular Epidemiological Study

BACKGROUND

The war in Ukraine has led to significant migration to neighboring countries, raising public health concerns. Notable tuberculosis (TB) incidence rates in Ukraine emphasize the immediate requirement to prioritize approaches that interrupt the spread and prevent new infections.

METHODS

We conducted a prospective genomic surveillance study to assess migration's impact on TB epidemiology in the Czech Republic and Slovakia. Mycobacterium tuberculosis isolates from Ukrainian war refugees and migrants, collected from September 2021 to December 2022 were analyzed alongside 1574 isolates obtained from Ukraine, the Czech Republic, and Slovakia.

RESULTS

Our study revealed alarming results, with historically the highest number of Ukrainian tuberculosis patients detected in the host countries. The increasing number of cases of multidrug-resistant TB, significantly linked with Beijing lineage 2.2.1 (p < 0.0001), also presents substantial obstacles to control endeavors. The genomic analysis identified the three highly related genomic clusters, indicating the recent TB transmission among migrant populations. The largest clusters comprised war refugees diagnosed in the Czech Republic, TB patients from various regions of Ukraine, and incarcerated individuals diagnosed with pulmonary TB specialized facility in the Kharkiv region, Ukraine, pointing to a national transmission sequence that has persisted for over 14 years.

CONCLUSIONS

The data showed that most infections were likely the result of reactivation of latent disease or exposure to TB before migration rather than recent transmission occurring within the host country. However, close monitoring, appropriate treatment, careful surveillance, and social support are crucial in mitigating future risks, though there is currently no evidence of local transmission in EU countries.

Analysis of Cdx2 VDR gene polymorphism rs11568820 in association with multiple sclerosis in Slovaks

Vitamin D deficiency is involved in the pathogenesis of multiple sclerosis (MS), a severe autoimmune demyelinating disease of the central nervous system. The gene polymorphism Cdx-2 (rs11568820, G/A) seriously influences the trancriptional activity of the vitamin D receptor (VDR) that binds the vitamin D responsive elements of target genes including HLA-DRB1*15. The aim of the present study in Slovaks was to analyse the association of Cdx-2 variants with the risk of MS and disability progression, and to assess the DRB1*15:01 allele as a possible confounding factor. In total, 493 MS patients and 417 healthy controls were involved in this study. The genotyping of Cdx-2 was performed using restriction analysis; DRB1*15:01 positivity was determined by a high-resolution melting analysis of its surrogate marker rs3135388 (G/A). Our results did not prove any allelic association between Cdx-2 and a risk of MS (minor allele A - 0.181 in patients vs. 0.161 in controls, OR = 1.15, .95 CI = 0.90-1.47, p = 0.289). The logistic regression analysis, adjusted for sex and age, showed no differences in Cdx-2 genotype counts when using an additive, dominant or recessive genetic model (p = 0.351, 0.150, 0.240 respectively). The Cdx-2 variants were also not associated with disease disability progression, evaluated using the Multiple Sclerosis Severity Score. The HLA-DRB1*15:01 allele was found to strongly increase the risk of MS in our study (0.300 in patients vs. 0.101 in controls, OR = 3.83, .95 CI = 2.94-4.99, p = 1.016 × 10-26, dominant genetic model OR = 4.62, .95 CI = 3.40-6.26, p = 9.1 × 10-23). In summary, we found the Cdx-2 as a single genetic marker not to be associated with MS development or progression in Slovaks, independently of HLA-DRB1*15:01 status.

Resistance patterns and transmission of mono- and polyresistant TB: clinical impact of WGS

Objectives

Rapidly diagnosing drug-resistant TB is crucial for improving treatment and transmission control. WGS is becoming increasingly accessible and has added value to the diagnosis and treatment of TB. The aim of the study was to perform WGS to determine the rate of false-positive results of phenotypic drug susceptibility testing (pDST) and characterize the molecular mechanisms of resistance and transmission of mono- and polyresistant Mycobacterium (M.) tuberculosis.

Methods

WGS was performed on 53 monoresistant and 25 polyresistant M. tuberculosis isolates characterized by pDST. Sequencing data were bioinformatically processed to infer mutations encoding resistance and determine the origin of resistance and phylogenetic relationship between isolates studied.

Results

The data showed the variable sensitivity and specificity of WGS in comparison with pDST as the gold standard: isoniazid 92.7% and 92.3%; streptomycin 41.9% and 100.0%; pyrazinamide 15% and 94.8%; and ethambutol 75.0% and 98.6%, respectively. We found novel mutations encoding resistance to streptomycin (in gidB) and pyrazinamide (in kefB). Most isolates belonged to lineage 4 (80.1%) and the overall clustering rate was 11.5%. We observed lineage-specific gene variations encoding resistance to streptomycin and pyrazinamide.

Conclusions

This study highlights the clinical potential of WGS in ruling out false-positive drug resistance following phenotypic or genetic drug testing, and recommend this technology together with the WHO catalogue in designing an optimal individualized treatment regimen and preventing the development of MDR TB. Our results suggest that resistance is primarily developed through spontaneous mutations or selective pressure.

Small-scale variants and large deletions in BRCA1/2 genes in Slovak high-grade serous ovarian cancer

The role of PARP inhibitors is to prevent the polymerase from repairing the single-strand break that occurred due to tumor growth and thus induce cell apoptosis when the homologous recombination deficiency (HRD) system is disabled. The eliminated system can be monitored especially in patients with serous ovarian epithelial tumors. Current studies still show the highest progression-free survival (PFS) in the examined groups with BRCA mutant status, even though they are also effective in the case of a disrupted HRD system, apart from BRCA genes. The study cohort consists of women diagnosed with high-grade serous ovarian cancer (HGSOC), after at least two lines of chemotherapy and after relapse of the disease, as determined by ESMO standards and guidelines. The commercially available tool SOPHIA DDM™ (SophiaGenetics, Switzerland) was used to classify the variants after sequencing. The most common variants (pathogenic or likely pathogenic) were in BRCA1 c.1067 A>G (rs1799950) and c.5266dupC (rs80357906) and in BRCA2 c.9976 A>T (rs11571833). Large deletions were detected in one and three cases in the BRCA1 and BRCA2 genes, respectively. A mutation in the BRCA1/2 genes was confirmed in 50% of the examined patients. In the study, we focused on the identification of mutated BRCA genes by a commercially available Sophia DDM™ system to identify a pathogenic or probable pathogenic variant in a cohort of patients with HGSOC in the Slovak population, which could result in better management and stratification of the individual.

Polymorphisms in the gene encoding CYP1A2 influence prostate cancer risk and progression

The role of the cytochrome P450 1A2 (CYP1A2) rs2472299, rs2470890 and rs11072508 polymorphisms in prostate cancer risk, disease progression and tumour development remains unclear. The potential associations of these three CYP1A2 polymorphisms and haplotypes with prostate cancer susceptibility and its clinicopathological characteristics were therefore investigated. The present case-control study consisted of 522 patients with prostate cancer and 554 healthy controls. High-resolution melting analysis was used to determine the CYP1A2 polymorphisms. No significant association in prostate cancer risk was seen for CYP1A2 rs2472299 and rs11072508. However, a significantly decreased risk of prostate cancer was found for CYP1A2 rs2470890 [odds ratio (OR), 0.67; P=0.02] in the recessive model. After analysis of the associations of clinical status and these three CYP1A2 polymorphisms, the CYP1A2 rs2470890 and rs11072508 polymorphisms showed a positive association with a higher Gleason score (rs2470890 OR, 1.36, P=0.04 in the allelic model; rs11072508 OR, 1.37, P=0.04 in the allelic model and OR, 1.60, P=0.03 in the dominant model). All three polymorphisms showed a significant positive association with pathological T stage in the additive, allelic and dominant genetic models (P<0.05). Haplotype analysis revealed that the most common haplotypes 'GTT' and 'ACC' were significantly associated with pathological T stages 3 and 4 (OR, 0.62; P=0.02 and OR, 1.54; P=0.03, respectively). A significant association was found between the 'GTT' haplotype and the Gleason score (OR, 0.71; P=0.03). In conclusion, these CYP1A2 polymorphisms and haplotypes have the potential to predict prostate cancer disease progression.

Systematic Genomic Surveillance of SARS-CoV-2 Virus on Illumina Sequencing Platforms in the Slovak Republic-One Year Experience

To explore a genomic pool of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) during the pandemic, the Ministry of Health of the Slovak Republic formed a genomics surveillance workgroup, and the Public Health Authority of the Slovak Republic launched a systematic national epidemiological surveillance using whole-genome sequencing (WGS). Six out of seven genomic centers implementing Illumina sequencing technology were involved in the national SARS-CoV-2 virus sequencing program. Here we analyze a total of 33,024 SARS-CoV-2 isolates collected from the Slovak population from 1 March 2021, to 31 March 2022, that were sequenced and analyzed in a consistent manner. Overall, 28,005 out of 30,793 successfully sequenced samples met the criteria to be deposited in the global GISAID database. During this period, we identified four variants of concern (VOC)-Alpha (B.1.1.7), Beta (B.1.351), Delta (B.1.617.2) and Omicron (B.1.1.529). In detail, we observed 165 lineages in our dataset, with dominating Alpha, Delta and Omicron in three major consecutive incidence waves. This study aims to describe the results of a routine but high-level SARS-CoV-2 genomic surveillance program. Our study of SARS-CoV-2 genomes in collaboration with the Public Health Authority of the Slovak Republic also helped to inform the public about the epidemiological situation during the pandemic.

Whole genome sequencing of multidrug-resistant Mycobacterium tuberculosis isolates collected in the Czech Republic, 2005-2020

The emergence and spread of resistant tuberculosis (TB) pose a threat to public health, so it is necessary to diagnose the drug-resistant forms in a clinically short time frame and closely monitor their transmission. In this study, we carried out a first whole genome sequencing (WGS)-based analysis of multidrug resistant (MDR) M. tuberculosis strains to explore the phylogenetic lineages diversity, drug resistance mechanisms, and ongoing transmission chains within the country. In total, 65 isolates phenotypically resistant to at least rifampicin and isoniazid collected in the Czech Republic in 2005-2020 were enrolled for further analysis. The agreement of the results obtained by WGS with phenotypic drug susceptibility testing (pDST) in the determination of resistance to isoniazid, rifampicin, pyrazinamide, streptomycin, second-line injectables and fluoroquinolones was more than 80%. Phylogenetic analysis of WGS data revealed that the majority of MDR M. tuberculosis isolates were the Beijing lineage 2.2.1 (n = 46/65; 70.8%), while the remaining strains belonged to Euro-American lineage. Cluster analysis with a predefined cut-off distance of less than 12 single nucleotide polymorphisms between isolates showed 19 isolates in 6 clusters (clustering rate 29.2%), located mainly in the region of the capital city of Prague. This study highlights the utility of WGS as a high-resolution approach in the diagnosis, characterization of resistance patterns, and molecular-epidemiological analysis of resistant TB in the country.

Anti-tuberculosis drug resistance in Slovakia, 2018-2019: The first whole-genome epidemiological study

•

The Euro-American lineage 4.7 was the most represented among the XDR-TB/MDR-TB resistant strains.

•

Seven of 12 isolates (58%) belonged to the same recent transmission chain.

•

Strains originally from Eastern European countries were diagnosed in patients without any epidemiological links with these countries.

•

Incidence of resistant XDR-TB/MDR-TB in Slovakia, while low, still should be taken to monitor.

Objective

The resistance of Mycobacterium (M.) tuberculosis to antituberculosis drugs poses a major threat to global public health. Whole genome sequencing (WGS) is an increasingly preferred method in the diagnostics and monitoring of the transmission dynamics of resistant forms of tuberculosis (TB). The aim of the study was to, for the first time, use the sequencing-based analysis to study the transmission and resistance patterns of a systematic and recent collection of extensively drug resistant (XDR) and multidrug resistant tuberculosis (MDR-TB) isolates and to expand our knowledge about drug resistant (DR) TB epidemiological dynamics in Slovakia.

Design

A total of 495 patients with pulmonary TB, who were referred to National Reference Laboratory for Mycobacteriology (Vyšné Hágy, Slovakia) in the years 2018–2019, were studied. Out of the total of 495 patients, 4 XDR-TB (0.8%) and 8 (1.6%) MDR-TB isolates were identified by conventional drug susceptibility testing on Löwenstein-Jensen solid medium and subjected to whole genome sequencing. Sequencing data were evaluated for molecular-epidemiological analysis and identification of resistance patterns.

Results

Phylogenetic and cluster analysis showed extensive recent transmission events and the predominance of Euro-American lineage 4.7 in Slovakia. However, phylogenetic analysis revealed the circulation of several lineages that originally occurred in Eastern European countries. Resistance patterns for first- and second-line antituberculosis drugs characterized by whole genome sequencing were in high concordance with the results of phenotypic drug susceptibility testing.

Conclusion

Forty percent of at least MDR-TB isolates were not genetically linked, indicating that appropriate measures should be taken to monitor and prevent the spread of drug-resistant tuberculosis within the country as well as in other regions.

Common gene haplotypes of gelatinases and their tissue inhibitors in abdominal aortic aneurysm

Abdominal aortic aneurysm (AAA) involves complex dynamic remodeling processes in the aortic wall. Gelatinases (MMP2 and MMP9) and their respective tissue inhibitors (TIMP1 and TIMP2) play a crucial role during extracellular matrix (ECM) turnover in aortic tissue. In this study we characterized associations between the haplotypes of genes encoding gelatinase/inhibitor pairs and pathways involved in AAA, a total of 100 AAA patients and 192 controls were enrolled. For males, a significant decrease in the distribution of the minor G allele of the TIMP2 rs8082025 was observed in AAA patients (p = 0.01, 23.1% controls vs. 13.1% AAA). In addition, in males, the major TIMP2 GA haplotype was associated with AAA (86.9% AAA vs. 76.9% control; p = 0.009, OR = 1.997), whereas the TIMP2 GG haplotype (7.7% AAA vs. 13.9% control) was associated with protection against AAA (p = 0.046, OR = 0.518). The minor GAGC MMP9 haplotype was related to AAA for all study subjects as well as the males only subset (p = 0.011, OR = 2.202 and p = 0.025, OR = 2.156, respectively). Small differences in the distribution of gene haplotypes could be associated with different levels of gene expression and in turn influence gelatinases activity in AAA.

DPYD genotype and haplotype analysis and colorectal cancer susceptibility in a case-control study from Slovakia

Dihydropyrimidine dehydrogenase (DPD) acts as the first-step enzyme catabolizing pyrimidines in vivo. DPYD gene mutations interfere with the breakdown of uracil and thymine. Genetic variations of DPYD can cause an enzyme deficiency state, which results in severe toxicity or other adverse side effects such as DNA damage or RNA damage caused by imbalance of the nucleotide pool. Our case-control study investigates the possible association between seven DPYD gene polymophisms (rs1801267, rs72547602, rs1801160, rs3918290, rs1801159, rs1801158, rs1801265) and risk of colorectal cancer (CRC). The association analysis for DPD was performed on 273 CRC patients and 187 healthy controls. There is significant allele association of SNP rs1801160 with colorectal cancer (p = 0.003, OR = 3.264, 95% CI = 1.425-7.475) in present analysis. Haplotype analysis of four DPYD polymorphisms showed significant difference in the distribution "IISt" haplotype between cases and controls. In comparison to the most common haplotype (VISt), the "IISt" haplotype was associated with increased risk for CRC (p = 0.038, OR = 2.733, 95% CI = 1.019-7.326). The present study suggests that the SNP rs1801160 and the "IISt" haplotype in the DPYD gene may also have a role in colorectal cancer risk.

Association of EGF, IGFBP-3 and TP53 Gene Polymorphisms with Major Depressive Disorder in Slovak Population

BACKGROUND

Major depressive disorder (MDD) is a main public health concern worldwide. Despite extensive investigations, the exact mechanisms responsible for MDD have not been identified. Epidermal growth factor (EGF) and insulin growth factor binding protein-3 (IGFBP-3) are involved in brain function. Tumour suppressor protein p53 is widely involved in neuronal death in response to different forms of acute insults and neurological disorders. The present study focuses on the possible associations of the single-nucleotide polymorphisms (SNP) of EGF A61G (rs4444903), IGFBP-3 C32G (rs2854746) and TP53 G72C (rs1042522) genes with MDD risk in the Slovak population.

METHODS

The present case-control association study was carried out in 111 confirmed MDD patients and 207 healthy subjects. Genotyping was performed by polymerase chain reaction-restriction fragment length polymorphism methods.

RESULTS

Logistic regression analysis showed no association between SNPs of selected genes and MDD risk in the Slovak population. However, the stratification of individuals by gender revealed that males carrying IGFBP-3 G alleles (G32G or GG) had marginally increased risk for developing MDD as compared to CC homozygous males (p=0.09). In women, inverse association was observed between SNP rs1042522 and MDD risk (p=0.04 for recessive model).

CONCLUSION

Our results suggest the protective effect of minor allele 72C of TP53 gene towards MDD. The disruption of mechanisms involved in cell survival and death regulation may be involved in pathophysiology of MDD.

GP6 Haplotype of Missense Variants is Associated with Sticky Platelet Syndrome Manifested by Fetal Loss

Disequilibrium of hemostasis is central to the pathogenesis of all thromboses, and platelets are essential for primary hemostasis. The platelet membrane glycoprotein receptor is involved in the clot formation in blood; therefore, the changes in related genes could impair platelet aggregation in patients with sticky platelet syndrome (SPS). Patients with SPS who experienced fetal loss were shown to harbor a risk haplotype at GP6 locus. The aim of the study was to examine the genetic linkage of this selected risk haplotype with single nucleotide variations (SNVs) in the coding sequence of the GP6 gene in order to identify possible functional SNVs in association with SPS and fetal loss. A total of 37 patients with SPS manifested fetal loss, and 42 healthy controls were enrolled in the study. The SPS was diagnosed with platelet aggregometry. The SNVs were determined by dideoxy sequencing and high-resolution melting analysis. The missense variations were detected in patients with risk haplotype only. The association analysis showed association of the minor alleles with the SPS manifested by fetal loss as follows-rs1671152 (odds ratio [OR]: 4.667, 95% confidence interval [CI]: 1.462-14.89, P = .006), rs2304167 (OR: 5.085, 95% CI: 1.605-16.10, P = .003), and rs1654416 (OR: 5.085, 95% CI: 1.605-16.10, P = .003). Using the Expectation-Maximization (EM) algorithm, the estimated minor haplotype with predicted protein residue PEAN was significantly associated with the given phenotype (OR: 4.746, 95% CI: 1.486-15.15, P = .005). We have shown that haplotype PEAN associated with SPS and manifested by fetal loss and suggest that the mechanism involved in the action of GPVI has significant effect on GPVI-mediated signal transduction through Syk-phosphorylation.

ApaI, BsmI and TaqI VDR gene polymorphisms in association with multiple sclerosis in Slovaks

OBJECTIVE

Vitamin D acts through vitamin D receptor (VDR) and has promising beneficial effects in the development and progression of multiple sclerosis (MS). The purpose of our study was to investigate the possible association of the VDR gene polymorphisms ApaI, BsmI and TaqI with the MS susceptibility and with the rate of disease disability progression in the Central European Slovak population.

METHODS

The allele and genotype variants of ApaI, TaqI and BsmI VDR gene polymorphisms were analysed in 270 clinically diagnosed MS patients and 303 healthy controls. Genotyping was performed by polymerase chain reaction and restriction analysis. Patients were stratified by the rate of disease disability progression using MSSS scores.

RESULTS

By logistic regression analysis, we revealed that genotype BB (AA) of BsmI VDR gene polymorphism is decreasing the risk of MS (BB (AA) vs Bb (AG) + bb (GG); OR = 0.59, 95% CI = 0.39-0.90, plog = 0.014). We did not identify any association of ApaI and TaqI polymorphisms neither with MS development nor with the disease progression.

DISCUSSION

We showed for the first time that BsmI genotype BB (AA) is associated with the decreased susceptibility to MS in Slovak population. We propose the BsmI gene polymorphism to be one of the important genetic markers in evaluation of the risk of MS. However, our data suggest that VDR gene polymorphisms ApaI, BsmI and TaqI are not useful in the prediction of disease disability progression rate in MS in Slovaks.