Browsing by Author "Milicevic, Ognjen (57211159715)"

Inferring transcriptional direction (orientation) of the CRISPR array is essential for many applications, including systematically investigating non-canonical CRISPR/Cas functions. The standard method, CRISPRDirection (embedded within CRISPRCasFinder), fails to predict the orientation (ND predictions) for ∼37% of the classified CRISPR arrays (>2200 loci); this goes up to >70% for the II-B subtype where non-canonical functions were first experimentally discovered. Alternatively, Potential Orientation (also embedded within CRISPRCasFinder), has a much smaller frequency of ND predictions but might have significantly lower accuracy. We propose a novel simple criterion, where the CRISPR array direction is assigned according to the direction of its associated cas genes (Cas Orientation). We systematically assess the performance of the three methods (Cas Orientation, CRISPRDirection, and Potential Orientation) across all CRISPR/Cas subtypes, by a mutual crosscheck of their predictions, and by comparing them to the experimental dataset. Interestingly, CRISPRDirection agrees much better with Cas Orientation than with Potential Orientation, despite CRISPRDirection and Potential Orientation being mutually related – Potential Orientation corresponding to one of six (heterogeneous) predictors employed by CRISPRDirection – and being unrelated to Cas Orientation. We find that Cas Orientation has much higher accuracy compared to Potential Orientation and comparable accuracy to CRISPRDirection – while accurately assigning an orientation to ∼95% of the CRISPR arrays that are non-determined by CRISPRDirection. Cas Orientation is, at the same time, simple to employ, requiring only (routine for prokaryotes) the prediction of the associated protein coding gene direction. © Copyright © 2019 Milicevic, Repac, Bozic, Djordjevic and Djordjevic.

Inferring transcriptional direction (orientation) of the CRISPR array is essential for many applications, including systematically investigating non-canonical CRISPR/Cas functions. The standard method, CRISPRDirection (embedded within CRISPRCasFinder), fails to predict the orientation (ND predictions) for ∼37% of the classified CRISPR arrays (>2200 loci); this goes up to >70% for the II-B subtype where non-canonical functions were first experimentally discovered. Alternatively, Potential Orientation (also embedded within CRISPRCasFinder), has a much smaller frequency of ND predictions but might have significantly lower accuracy. We propose a novel simple criterion, where the CRISPR array direction is assigned according to the direction of its associated cas genes (Cas Orientation). We systematically assess the performance of the three methods (Cas Orientation, CRISPRDirection, and Potential Orientation) across all CRISPR/Cas subtypes, by a mutual crosscheck of their predictions, and by comparing them to the experimental dataset. Interestingly, CRISPRDirection agrees much better with Cas Orientation than with Potential Orientation, despite CRISPRDirection and Potential Orientation being mutually related – Potential Orientation corresponding to one of six (heterogeneous) predictors employed by CRISPRDirection – and being unrelated to Cas Orientation. We find that Cas Orientation has much higher accuracy compared to Potential Orientation and comparable accuracy to CRISPRDirection – while accurately assigning an orientation to ∼95% of the CRISPR arrays that are non-determined by CRISPRDirection. Cas Orientation is, at the same time, simple to employ, requiring only (routine for prokaryotes) the prediction of the associated protein coding gene direction. © Copyright © 2019 Milicevic, Repac, Bozic, Djordjevic and Djordjevic.

Background: We aimed to assess the extent to which the buffy coat DNA methylome is representative of methylation patterns in constitutive white blood cell (WBC) types in normal pregnancy. Methods: A comparison of differential methylation of buffy coat DNA vs DNA isolated from polymorphonuclear (PMN) and lymphocytic fractions was performed for each blood sample obtained within 24 h prior to delivery from 29 normotensive pregnant women. Methylation profiles were obtained using an Illumina Human Methylation 450 BeadChip and CHaMP bioinformatics pipeline. A subset of differentially methylated probes (DMPs) showing discordant methylation were further investigated using statistical modeling and enrichment analysis. Results: The smallest number of DMPs was found between the buffy coat and the PMN fraction (2.96%). Pathway enrichment analysis of the DMPs identified biological pathways involved in the particular leukocyte lineage, consistent with perturbations during isolation. The comparisons between the buffy coat and the isolated fractions as a group using linear modeling yielded a small number of probes (∼29,000) with discordant methylation. Demethylation of probes in the buffy coat compared to derived cell lines was more common and was prevalent in shelf and open sea regions. Conclusion: Buffy coat is representative of methylation patterns in WBC types in normal pregnancy. The differential methylations are consistent with perturbations during isolation of constituent cells and likely originate in vitro due to the physical stress during cell separation and are of no physiological relevance. These findings help the interpretation of DNA methylation profiling in pregnancy and numerous other conditions. Copyright © 2022 Ghamrawi, Velickovic, Milicevic, White, Thistlethwaite, Cunningham, Milosavljevic, Milic and Garovic.

Background: We aimed to assess the extent to which the buffy coat DNA methylome is representative of methylation patterns in constitutive white blood cell (WBC) types in normal pregnancy. Methods: A comparison of differential methylation of buffy coat DNA vs DNA isolated from polymorphonuclear (PMN) and lymphocytic fractions was performed for each blood sample obtained within 24 h prior to delivery from 29 normotensive pregnant women. Methylation profiles were obtained using an Illumina Human Methylation 450 BeadChip and CHaMP bioinformatics pipeline. A subset of differentially methylated probes (DMPs) showing discordant methylation were further investigated using statistical modeling and enrichment analysis. Results: The smallest number of DMPs was found between the buffy coat and the PMN fraction (2.96%). Pathway enrichment analysis of the DMPs identified biological pathways involved in the particular leukocyte lineage, consistent with perturbations during isolation. The comparisons between the buffy coat and the isolated fractions as a group using linear modeling yielded a small number of probes (∼29,000) with discordant methylation. Demethylation of probes in the buffy coat compared to derived cell lines was more common and was prevalent in shelf and open sea regions. Conclusion: Buffy coat is representative of methylation patterns in WBC types in normal pregnancy. The differential methylations are consistent with perturbations during isolation of constituent cells and likely originate in vitro due to the physical stress during cell separation and are of no physiological relevance. These findings help the interpretation of DNA methylation profiling in pregnancy and numerous other conditions. Copyright © 2022 Ghamrawi, Velickovic, Milicevic, White, Thistlethwaite, Cunningham, Milosavljevic, Milic and Garovic.

Understanding variations in the severity of infectious diseases is essential for planning proper mitigation strategies. Determinants of COVID-19 clinical severity are commonly assessed by transverse or longitudinal studies of the fatality counts. However, the fatality counts depend both on disease clinical severity and transmissibility, as more infected also lead to more deaths. Instead, we use epidemiological modeling to propose a disease severity measure that accounts for the underlying disease dynamics. The measure corresponds to the ratio of population-averaged mortality and recovery rates (m/r), is independent of the disease transmission dynamics (i.e., the basic reproduction number), and has a direct mechanistic interpretation. We use this measure to assess demographic, medical, meteorological, and environmental factors associated with the disease severity. For this, we employ an ecological regression study design and analyze different US states during the first disease outbreak. Principal Component Analysis, followed by univariate, and multivariate analyses based on machine learning techniques, is used for selecting important predictors. The usefulness of the introduced severity measure and the validity of the approach are confirmed by the fact that, without using prior knowledge from clinical studies, we recover the main significant predictors known to influence disease severity, in particular age, chronic diseases, and racial factors. Additionally, we identify long-term pollution exposure and population density as not widely recognized (though for the pollution previously hypothesized) significant predictors. The proposed measure is applicable for inferring severity determinants not only of COVID-19 but also of other infectious diseases, and the obtained results may aid a better understanding of the present and future epidemics. Our holistic, systematic investigation of disease severity at the human-environment intersection by epidemiological dynamical modeling and machine learning ecological regressions is aligned with the One Health approach. The obtained results emphasize a syndemic nature of COVID-19 risks. © 2021 The Authors

Citations are an important, but often overlooked, part of every scientific paper. They allow the reader to trace the flow of evidence, serving as a gateway to relevant literature. Most scientists are aware of citations' errors, but few appreciate the prevalence of these problems. The purpose of the present study was to examine how often frequently cited papers in biomedical scientific literature are cited inaccurately. The study included an active participation of the first authors of included papers; to first-hand verify the citations accuracy. Findings from feasibility study, where we reviewed 1540 articles containing 2526 citations of 14 most cited articles in which the authors were affiliated with the Faculty of Medicine University of Belgrade, were further evaluated for external confirmation in an independent verification set of articles. Verification set included 4912 citations identified in 2995 articles that cited 13 most cited articles published by authors affiliated with the Mayo Clinic Division of Nephrology and Hypertension. A citation was defined as being accurate if the cited article supported or was in accordance with the statement by citing authors. At least one inaccurate citation was found in 11 and 15% of articles in the feasibility study and verification set, respectively, suggesting that inaccurate citations are common in biomedical literature. The most common problem was the citation of nonexistent findings (38.4%), followed by an incorrect interpretation of findings (15.4%). One-fifth of inaccurate citations were due to chains of inaccurate citations. Based on these findings, several actions to reduce citation inaccuracies have been proposed. © 2021 The Author(s).

Identifying the main environmental drivers of SARS-CoV-2 transmissibility in the population is crucial for understanding current and potential future outbursts of COVID-19 and other infectious diseases. To address this problem, we concentrate on the basic reproduction number R0, which is not sensitive to testing coverage and represents transmissibility in an absence of social distancing and in a completely susceptible population. While many variables may potentially influence R0, a high correlation between these variables may obscure the result interpretation. Consequently, we combine Principal Component Analysis with feature selection methods from several regression-based approaches to identify the main demographic and meteorological drivers behind R0. We robustly obtain that country's wealth/development (GDP per capita or Human Development Index) is the most important R0 predictor at the global level, probably being a good proxy for the overall contact frequency in a population. This main effect is modulated by built-up area per capita (crowdedness in indoor space), onset of infection (likely related to increased awareness of infection risks), net migration, unhealthy living lifestyle/conditions including pollution, seasonality, and possibly BCG vaccination prevalence. Also, we argue that several variables that significantly correlate with transmissibility do not directly influence R0 or affect it differently than suggested by naïve analysis. © 2021 The Authors. GeoHealth published by Wiley Periodicals LLC on behalf of American Geophysical Union.

Identifying the main environmental drivers of SARS-CoV-2 transmissibility in the population is crucial for understanding current and potential future outbursts of COVID-19 and other infectious diseases. To address this problem, we concentrate on the basic reproduction number R0, which is not sensitive to testing coverage and represents transmissibility in an absence of social distancing and in a completely susceptible population. While many variables may potentially influence R0, a high correlation between these variables may obscure the result interpretation. Consequently, we combine Principal Component Analysis with feature selection methods from several regression-based approaches to identify the main demographic and meteorological drivers behind R0. We robustly obtain that country's wealth/development (GDP per capita or Human Development Index) is the most important R0 predictor at the global level, probably being a good proxy for the overall contact frequency in a population. This main effect is modulated by built-up area per capita (crowdedness in indoor space), onset of infection (likely related to increased awareness of infection risks), net migration, unhealthy living lifestyle/conditions including pollution, seasonality, and possibly BCG vaccination prevalence. Also, we argue that several variables that significantly correlate with transmissibility do not directly influence R0 or affect it differently than suggested by naïve analysis. © 2021 The Authors. GeoHealth published by Wiley Periodicals LLC on behalf of American Geophysical Union.

The healthcare environment in Serbia has changed dramatically over a last two decades, pointing out the necessity of clinical informatics (CI) education for future MDs. Total of 77 students were enrolled and 72 (93.5%) have successfully finished this course during 4 academic years. Mean total score for all students was 83.4 ± 9.0 points, without difference between genders. We presented blended learning module as an effective way of gaining competences in CI and recommend this course to be required for future MDs. © 2020 European Federation for Medical Informatics (EFMI) and IOS Press.

The healthcare environment in Serbia has changed dramatically over a last two decades, pointing out the necessity of clinical informatics (CI) education for future MDs. Total of 77 students were enrolled and 72 (93.5%) have successfully finished this course during 4 academic years. Mean total score for all students was 83.4 ± 9.0 points, without difference between genders. We presented blended learning module as an effective way of gaining competences in CI and recommend this course to be required for future MDs. © 2020 European Federation for Medical Informatics (EFMI) and IOS Press.

Introduction: Preeclampsia (PE) is a pregnancy-associated, multi-organ, life-threatening disease that appears after the 20th week of gestation. The aim of this study was to perform a systematic review and meta-analysis to determine whether women with PE have disrupted miRNA expression compared to women who do not have PE. Methods: We conducted a systematic review and meta-analysis of studies that reported miRNAs expression levels in placenta or peripheral blood of pregnant women with vs. without PE. Studies published before October 29, 2021 were identified through PubMed, EMBASE and Web of Science. Two reviewers used predefined forms and protocols to evaluate independently the eligibility of studies based on titles and abstracts and to perform full-text screening, data abstraction and quality assessment. Standardized mean difference (SMD) was used as a measure of effect size. Results: 229 publications were included in the systematic review and 53 in the meta-analysis. The expression levels in placenta were significantly higher in women with PE compared to women without PE for miRNA-16 (SMD = 1.51,95%CI = 0.55–2.46), miRNA-20b (SMD = 0.89, 95%CI = 0.33–1.45), miRNA-23a (SMD = 2.02, 95%CI = 1.25–2.78), miRNA-29b (SMD = 1.37, 95%CI = 0.36–2.37), miRNA-155 (SMD = 2.99, 95%CI = 0.83–5.14) and miRNA-210 (SMD = 1.63, 95%CI = 0.69–2.58), and significantly lower for miRNA-376c (SMD = –4.86, 95%CI = –9.51 to –0.20). An increased level of miRNK-155 expression was found in peripheral blood of women with PE (SMD = 2.06, 95%CI = 0.35–3.76), while the expression level of miRNA-16 was significantly lower in peripheral blood of PE women (SMD = –0.47, 95%CI = –0.91 to –0.03). The functional roles of the presented miRNAs include control of trophoblast proliferation, migration, invasion, apoptosis, differentiation, cellular metabolism and angiogenesis. Conclusion: miRNAs play an important role in the pathophysiology of PE. The identification of differentially expressed miRNAs in maternal blood creates an opportunity to define an easily accessible biomarker of PE. Copyright © 2021 Cirkovic, Stanisavljevic, Milin-Lazovic, Rajovic, Pavlovic, Milicevic, Savic, Kostic Peric, Aleksic, Milic, Stanisavljevic, Mikovic, Garovic and Milic.

Introduction: Preeclampsia (PE) is a pregnancy-associated, multi-organ, life-threatening disease that appears after the 20th week of gestation. The aim of this study was to perform a systematic review and meta-analysis to determine whether women with PE have disrupted miRNA expression compared to women who do not have PE. Methods: We conducted a systematic review and meta-analysis of studies that reported miRNAs expression levels in placenta or peripheral blood of pregnant women with vs. without PE. Studies published before October 29, 2021 were identified through PubMed, EMBASE and Web of Science. Two reviewers used predefined forms and protocols to evaluate independently the eligibility of studies based on titles and abstracts and to perform full-text screening, data abstraction and quality assessment. Standardized mean difference (SMD) was used as a measure of effect size. Results: 229 publications were included in the systematic review and 53 in the meta-analysis. The expression levels in placenta were significantly higher in women with PE compared to women without PE for miRNA-16 (SMD = 1.51,95%CI = 0.55–2.46), miRNA-20b (SMD = 0.89, 95%CI = 0.33–1.45), miRNA-23a (SMD = 2.02, 95%CI = 1.25–2.78), miRNA-29b (SMD = 1.37, 95%CI = 0.36–2.37), miRNA-155 (SMD = 2.99, 95%CI = 0.83–5.14) and miRNA-210 (SMD = 1.63, 95%CI = 0.69–2.58), and significantly lower for miRNA-376c (SMD = –4.86, 95%CI = –9.51 to –0.20). An increased level of miRNK-155 expression was found in peripheral blood of women with PE (SMD = 2.06, 95%CI = 0.35–3.76), while the expression level of miRNA-16 was significantly lower in peripheral blood of PE women (SMD = –0.47, 95%CI = –0.91 to –0.03). The functional roles of the presented miRNAs include control of trophoblast proliferation, migration, invasion, apoptosis, differentiation, cellular metabolism and angiogenesis. Conclusion: miRNAs play an important role in the pathophysiology of PE. The identification of differentially expressed miRNAs in maternal blood creates an opportunity to define an easily accessible biomarker of PE. Copyright © 2021 Cirkovic, Stanisavljevic, Milin-Lazovic, Rajovic, Pavlovic, Milicevic, Savic, Kostic Peric, Aleksic, Milic, Stanisavljevic, Mikovic, Garovic and Milic.

March 6, 2020 is considered as the official date of the beginning of the COVID-19 epidemic in Serbia. In late spring and early summer 2020, Europe recorded a decline in the rate of SARS-CoV-2 infection and subsiding of the first wave. This trend lasted until the fall, when the second wave of the epidemic began to appear. Unlike the rest of Europe, Serbia was hit by the second wave of the epidemic a few months earlier. Already in June 2020, newly confirmed cases had risen exponentially. As the COVID-19 pandemic is the first pandemic in which there has been instant sharing of genomic information on isolates around the world, the aim of this study was to analyze whole SARS-CoV-2 viral genomes from Serbia, to identify circulating variants/clade/lineages, and to explore site-specific mutational patterns in the unique early second wave of the European epidemic. This analysis of Serbian isolates represents the first publication from Balkan countries, which demonstrates the importance of specificities of local transmission especially when preventive measures differ among countries. One hundred forty-eight different genome variants among 41 Serbian isolates were detected in this study. One unique and seven extremely rare mutations were identified, with locally specific continuous dominance of the 20D clade. At the same time, amino acid substitutions of newly identified variants of concern were found in our isolates from October 2020. Future research should be focused on functional characterization of novel mutations in order to understand the exact role of these variations. © Copyright © 2021 Miljanovic, Milicevic, Loncar, Abazovic, Despot and Banko.

March 6, 2020 is considered as the official date of the beginning of the COVID-19 epidemic in Serbia. In late spring and early summer 2020, Europe recorded a decline in the rate of SARS-CoV-2 infection and subsiding of the first wave. This trend lasted until the fall, when the second wave of the epidemic began to appear. Unlike the rest of Europe, Serbia was hit by the second wave of the epidemic a few months earlier. Already in June 2020, newly confirmed cases had risen exponentially. As the COVID-19 pandemic is the first pandemic in which there has been instant sharing of genomic information on isolates around the world, the aim of this study was to analyze whole SARS-CoV-2 viral genomes from Serbia, to identify circulating variants/clade/lineages, and to explore site-specific mutational patterns in the unique early second wave of the European epidemic. This analysis of Serbian isolates represents the first publication from Balkan countries, which demonstrates the importance of specificities of local transmission especially when preventive measures differ among countries. One hundred forty-eight different genome variants among 41 Serbian isolates were detected in this study. One unique and seven extremely rare mutations were identified, with locally specific continuous dominance of the 20D clade. At the same time, amino acid substitutions of newly identified variants of concern were found in our isolates from October 2020. Future research should be focused on functional characterization of novel mutations in order to understand the exact role of these variations. © Copyright © 2021 Miljanovic, Milicevic, Loncar, Abazovic, Despot and Banko.