Logo

Publikacije (45085)

Nazad
V. Šuljagić, M. Bajčetić, V. Mioljević, G. Dragovac, B. Mijovic, I. Janicijevic, Zorana Đorđevic, G. Krtinić et al.

As the only non-European Union (EU) country, Serbia participated in a second point prevalence survey of healthcare-associated infections (HAIs) and antimicrobial use (AMU) organized by the European Centre for Disease Prevention and Control (ECDC) in the EU countries. Here, we aimed to estimate the prevalence of HAI and AMU in patients who had recently undergone a surgery and to compare risk profile, HAI rates, and AMU among surgical patients and non-surgical patients. A national PPS was performed in 65 Serbian acute-care hospitals, in November 2017. In this paper, the data of 61 hospitals for adult acute-care were analyzed. To ensure the comparability of study design we used the Serbian translation of ECDC case definitions and ECDC PPS protocol. The trained infection control staff, led by a hospital coordinator, reviewed medical records to identify HAI active at the time of the survey and AMU. Only inpatients admitted to the ward before 8 a.m. on the day if the survey were included. A total of 12,478 patients from 61 hospitals for adult acute-care were eligible for inclusion in this study. Significantly higher proportions of surgical patients were female, belonged to the 60-to-79 age group, and were less severely ill. Also, extrinsic factors (invasive devices, hospitalization at the ICU, and prior antibiotics therapy) were more frequent in surgical patients. Prevalence of HAIs was higher among surgical patients (261/3626; 7.2%) than among non-surgical patients (258/8852; 2.9%) (p < 0.0001). The highest prevalence of all HAIs was noted in patients who had kidney transplantation (4/11; 36.4%), while SSIs were the most prevalent among patients who had peripheral vascular bypass surgery (3/15; 20.0%). Non-surgical patients received treatment for community-acquired infections in significantly higher proportion (2664/8852; 64.3) (p < 0.001). Surgical prophylaxis for more than 1 day was applied in 71.4% of surgical patients. We have provided an insight into the burden of HAIs and AMU among Serbia acute-care hospitals, and highlighted several priority areas and targets for quality improvement.

Dženan Pleho, Amra Mačak Hadžiomerović, Kenan Pleho, Jasmina Pleho, Dinko Remić, Davor Arslanagić, Milojko Lazić, Aldina Alibegović

Introduction: Work-related musculoskeletal disorders (WRMSDs) are the most common work-related diseases and describe a wide range of degenerative and inflammatory conditions affecting blood vessels, peripheral nerves, joints, ligaments, tendons, and muscles. WRMSDs are becoming an increasing problem in modern society. They are the second biggest cause of short-term or temporary incapacity for work just following a cold. At the workplace, health professionals represent a very vulnerable category in terms of musculoskeletal disorders (MSDs) development. It is estimated that almost one-third of all cases of absence from work among health care professionals are related to MSDs. Studies also show that a large number of health professionals report the occurrence of MSDs in one or more regions of the body, with a problem with the lower back being one of the most common.Methods: This article presents a non-experimental (qualitative) research, or a scientific review of the published literaturewhere the databases were reviewed in which the keywords for the review were: MSDs, work, health, intervention program,and ergonomics. Various databases were used in the preparation of this article, including PubMed, Google Scholar, Medline, Hrcak, Dabar, Science Direct, and Science Citation.Results: The results include a review and analysis of eighteen published scientific articles in the period 2001-2020. The studies published in these articles has been conducted in the United States, Australia, Switzerland, Portugal, Slovenia, Turkey, China, Nigeria, Israel, Tunisia, Iran, Croatia, Taiwan, Bangladesh, Serbia, Macedonia, and Greece.Conclusion: This article should indicate the magnitude of the problems of WRMSDs in health professionals and that this topic is an inexhaustible and very interesting basis for further studies by current and future researchers to create strategies for the prevention and treatment of this disorder and to eliminate its causes.

N. Stefanovic, Dejan Jerkan

U radu je razvijen matematički model trofaznog, dvonamotajnog energetskog transformatora sprege Dyn5, podesan za proračune ustaljenih složeno­periodičnih režima. Model obezbeđuje uvažavanje postojanja nenultog sprežnog broja transformatora, omogućavajući proračune ustaljenih složenoperiodičnih radnih režima uz istovremeno prisustvo nesimetrije i/ili neuravnoteženosti. Na ilustrativnim primerima je pokazano kakav je uticaj energetskih transformatora na uspostavljanje složenoperiodičnih režima usled harmo­nijskog izobličenja napona napajanja na njegovim primarnim priključcima.

Geoflly L. Adonias, H. Šiljak, Michael Taynnan Barros, S. Balasubramaniam

The COVID-19 pandemic has shaken the world unprecedentedly, where it has affected the vast global population both socially and economically. The pandemic has also opened our eyes to the many threats that novel virus infections can pose for humanity. While numerous unknowns are being investigated in terms of the distributed damage that the virus can do to the human body, recent studies have also shown that the infection can lead to lifelong sequelae that could affect other parts of the body, and one example is the brain. As part of this work, we investigate how viral infection can affect the brain by modelling and simulating a neuron's behaviour under demyelination that is affected by the cytokine storm. We quantify the effects of cytokine-induced demyelination on the propagation of action potential signals within a neuron. We used information and communication theory analysis on the signal propagated through the axonal pathway under different intensity levels of demyelination to analyse these effects. Our simulations demonstrate that virus-induced degeneration can play a role in the signal power and spiking rate and the probability of releasing neurotransmitters and compromising the propagation and processing of information between the neurons. We also propose a transfer function that models these attenuation effects that degenerates the action potential, where this model has the potential to be used as a framework for the analysis of virus-induced neurodegeneration that can pave the way to improved understanding of virus-induced demyelination.

M. Kovačič, Shpetim Salihu, G. Gantar, U. Župerl

In this paper, machinability influences from the start to end of final product production in a steel plant were analyzed, including chemical composition, deoxidizing agents and casting parameters, which drastically influence the macrostructure and segregation (i.e., chemical nonhomogeneity) of continuously cast and subsequently rolled material. The data (seven parameters from secondary metallurgy, four parameters from the casting process and the content of ten chemical elements) from the serial production of calcium-treated steel grades (254 batches of 25 different steel grades from January 2018 to March 2020) were used for predicting machinability. Machinability was determined based on ISO 3685:1993, where the machinability of each individual batch is represented as the cutting speed and the tool is worn out within fifteen minutes. For the prediction of these cutting speeds, linear regression and genetic programming were used. Out of 25 analyzed steel grades, 20MnV6 steel grade was the most problematic and also the most often produced. Out of 57 produced batches of 20MnVS6 steel, 23 batches had nonconforming machinability. Based on the modeling results, the steelmaking process was optimized. Consequently, 40 additional batches of 20MnV6 (from March 2020 to July 2020) were subsequently produced based on an optimized steelmaking process. In all 40 cases, the required machinability was achieved without changing other properties required by the customers.

Muamer Dervisevic, M. Alba, Li Yan, M. Şenel, T. Gengenbach, B. Prieto‐Simón, N. Voelcker

Current technology for blood glucose level monitoring is mainly based on the invasive finger‐prick extraction of a small drop of blood using a lancet and measured via a handheld glucometer, which is not conducive to continuous measurements. Interstitial fluid (ISF) is gaining attention as an alternative biofluid. Its biochemical composition is very similar to that of blood and it can be monitored in a continuous manner via minimally invasive methods that cause no pain and minimize any risk of infection. Herein, a microneedle array (MNA) based transdermal sensing system for the pain free monitoring of ISF glucose is presented. High‐density silicon microneedles (≈9500 microneedles cm−2) are used to prepare a three‐electrode patch for the electrochemical monitoring of glucose. The MNA glucose patch shows very good selectivity when tested in artificial ISF, with a sensitivity of 0.1622 µA mm−1 cm−2 and a detection limit of 0.66 mm. In vivo application of the microneedle array in mice shows that the ISF glucose concentrations obtained with the MNA sensor gave very good correlation with the blood glucose levels determined with a commercial glucometer. This microneedle‐based sensing system hence provides an alternative transdermal diagnostic tool to the invasive existing techniques.

Race walking (RW) is a cyclic athletic discipline, of aerobic character, which is characterized by high intensity and high energy consumption. The result is correlated with the anthropological dimensions of the walker as well as good technical performance. The study includes 25 male students, third year of study, Faculty of Physical Education and Sports in East Sarajevo. For the purposes of the research, the results of RW at 2km (criterion variable) were measured, on the basis of which average values of walking speed (WS) were calculated. The aim of the research was to determine the correlation between the achieved result depending on the walking speed. By applying the Pearson product of the correlation moment, the obtained results confirmed a significant and high correlation between the result and the walking speed (r = -0.988; p<0.05) with an inverse function that confirms the high dependence of the result on the WS. Article visualizations:

J. Kamberović, M. Gligora Udovič, M. Kahlert, K. Tapolczai, Zoran Lukić, A. Ahmić, Anita Dedić, Antonija Kulaš et al.

The process of travertine formation and carbonate deposition in the rivers is unique, delicate, and depends on the activity of algae and mosses. Although diatoms have been used extensively in hydrobiological studies, the comparative analysis data on diatom communities of the travertine barriers in karstic rivers are still scarce. ‡ § | ¶ ‡ ‡ # § ¤,« »,˄ ˄,» »,˄ © Kamberović J et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. The study aimed to detect the diatom composition on travertine barriers in the Una River, the large karstic river in Bosnia and Herzegovina. An integrated classical morphological identification approach with metabarcoding was applied on eight samples across the river length profile. Morphological analyses were performed using both light and scanning electron microscopes. Subsequent DNA metabarcoding of the chloroplastic gene 312bp rbcL was done. The DADA2 pipeline was used for the bioinformatic treatment of the demultiplexed MiSeq reads to infer Amplicon Sequence Variants (ASVs). ASVs were taxonomically assigned using the Diat.barcode v7 reference database. A total of 126 species were identified using the morphological approach, while 133 ASVs were taxonomically assigned to 58 unique taxa with the molecular approach. Diatom community structures in terms of molecular and morphological approaches were congruent with 49 shared species. Species from genera Gomphonema, Navicula and Encyonema were less assigned in molecular analysis. The most abundant taxa in the Una River are alkaliphilous, belonging to the genera Gomphonema, Nitzshia and Navicula. Although specific for their extremely good chemical status, the travertine barriers of the Una River are largely inhabited with meso-eutraphentic taxa.

Dalila Destanović, Lejla Ušanović, J. Hanjalić, Lejla Lasić, B. Kalamujić Stroil

Bosnia and Herzegovina has valuable natural resources with a high percentage of endemic and autochthonous species (Kučinić et al. 2008, Đug and Drešković 2012). The freshwater fauna of Trichoptera in this area is under-investigated, with a lack of morphological description of different life stages and DNA barcode data. Public data show 58,993 barcode entries for Trichoptera in the Barcode of Life Data Systems (BOLD) submitted from 92 countries, and none from Bosnia and Herzegovina (B&H) (BOLD 2021). Previous research in Bosnia and Herzegovina has provided the first DNA barcode for the endemic species Rhyacophila bosnica, stored in GeneBank, under accession number MK211322 by a domestic institution (Kalamujić Stroil et al. 2018). A few DNA barcodes of adult individuals of Trichoptera from Bosnia and Herzegovina were found in BOLD. However, these specimens were collected on B&H territory, but analyzed, processed, and stored by foreign institutions. To change the current state of DNA barcoding of Trichoptera in Bosnia and Herzegovina, we aimed to employ this approach in investigating caddisflies in selected habitats in the Sarajevo Canton. Our fieldwork was done in all five protected areas (spring of the Bosna River, Bijambare, Trebević, Skakavac, and Bentbaša) in which larvae samples were collected according to ‡,§ ‡,§ ‡,§ ‡,§ ‡ © Destanović D et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. the AQEM sampling methodology. The standard animal DNA barcode was successfully obtained using degenerated primers LCO1490JJ and HCO2198-JJ (Astrin and Stüben 2008). Out of 684 collected individuals (313 Trebević, 130 spring of the Bosna River, 117 Bijambare, 71 Bentbaša, 53 Skakavac), a subset of specimens were sequenced. We uncovered 14 different taxa, 11 genera and six families (Limnephilidae, Glossosomatidae, Rhyacophilidae, Goeridae, Hydropsychidae, Polycentropodidae). The preliminary data of Trichoptera composition in the Sarajevo Canton indicated species richness. Based on our sequential data, a new subspecies was discovered in two investigated areas (Valladolid et al. 2020), proving that Trichoptera species diversity in our country is far from entirely uncovered. The benefit and power of the DNA barcoding approach are that it can pinpoint the areas of vast and unknown species diversity more economically, both financially and temporarily, than the morphological approach. Therefore, we believe that it is critical to support the development of DNA barcoding for the bioassessment of freshwater ecosystems in Bosnia and Herzegovina. Several problems prevented us from exploiting sequential data to the fullest. Despite a general notion among scientists that European Trichoptera species are well covered in the BOLD database, most of the sequences we obtained were absent from the database. Secondly, we recognized that morphological data about the larval developmental stage of B&H Trichoptera species are largely missing. The unified, updated, and complete data on this order of insects is urgently needed. However, insufficient financial support by governmental institutions and lack of systematic approach to barcoding the wildlife of Bosnia and Herzegovina hampers this process. Further attempts to collaborate with the stakeholders can be crucial with profound and substantial implications for biomonitoring of aquatic macroinvertebrates in general. New approaches, such as novel DNA barcoding-based methodology can fill an important gap in our knowledge of Balkan caddisflies haplotypes, lineages, and their diversification and distribution patterns.

Dalila Destanović, Lejla Ušanović, Lejla Lasić, J. Hanjalić, B. Kalamujić Stroil

Chaetopteryx villosa (Fabricius, 1798) is a caddisfly species distributed throughout Europe, except in the Balkan and Apennine Peninsula. However, phylogenetically close species belonging to the C. villosa group are widespread throughout entire Europe. Species of this group (C. villosa, C. gessneri, C. fusca, C. sahlbergi, C. atlantica, C. bosniaca, C. vulture, and C. trinacriae) have distinct distributions with some overlaps. Adult forms of these species are morphologically similar, whereas larval morphology is only known for some species. There are also indications of species hybridization (e.g., C. villosa x fusca). Presumably, the molecular approach for the species determination of this group would be highly beneficial. In the BOLD database, there are 154 specimens with COI-5P barcodes of C. villosa species. Out of the remaining species, C. sahlbergi has 27 specimens with a barcode, C. fusca 20, C. gessneri 5, C. bosniaca 5, and C. atlantica 1, whereas sequences from the species C. vulture and C. trinacriae are missing. Therefore, we tested the power of discrimination of the COI-5P marker in the C. villosa group, as the most common barcoding markers for species identification in animals. Only sequences from public records originating from experienced research groups or taxonomists and containing a specimen photograph were taken as input. A total of 75 ‡,§ ‡,§ ‡,§ ‡,§ ‡,§ © Destanović D et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. sequences from the BOLD database were obtained. Out of these sequences, 11 belonged to C. fusca, 5 to C. gessneri, 52 to C. villosa, 5 to C. bosniaca, and 2 to C. sahlbergi. For the generation of overview trees, COI-5P barcodes of Rhyacophila fasciata and Rh. nubila were used as outgroups. All sequences were trimmed at 5’ and 3’ ends, resulting in a final alignment length of 516 base pairs. Multiple sequence alignments and editing were done in the MEGA-X software. Analysis of nucleotide polymorphism was done in DNASP6 software. MEGA-X was used to calculate the pairwise distance and overall mean pdistance, and to construct the overview trees. Analysis of DNA polymorphism revealed 14 haplotypes of C. villosa, 3 haplotypes of C. fusca, 2 haplotypes of C. gessneri, and one for species C. bosniaca and C. sahlbergi. There were no significant interspecific and intraspecific differences among haplotypes based on pairwise distances. The p-distance between one of the haplotypes of C. fusca and C. villosa was 0.000, whereas the p-distance among haplotypes of C. villosa varied from 0.001 to about 0.055. The mean overall p-distance among haplotypes of all species equaled 0.03. No species-specific clusters were observed when phylogenetic trees were constructed except for C. gessneri, regardless of the method used (i.e., NJ, UPGMA, ML, ME, or MP). To minimize the possibility of species misidentification, we used only records submitted by NTNU-Norwegian University of Science and Technology (Norway), SNSB-Zoologische Staatssammlung Muenchen (Germany), Zoologisches Forschungsmuseum Alexander Koenig (Germany), University of Oulu, Zoological Museum (Finland), prof Hans Malicky and prof Mladen Kučinić. No records identified as hybrids were included in the analyses. With the exception of C. gessneri, COI-5P marker failed to separate the species of the C. villosa group. However, it is highly unlikely that poor species determination was the basis for such a result. To enable the comprehensive and unbiased evaluation of the relationships within this group, data coverage in BOLD database for most of the studied species should be enhanced, encompassing different geographical distribution of samples. Further studies are needed to detect the array of molecular markers suitable for the species delineation in a complex group such as C. villosa.

Darija Bilandžija, Marija Galić, Ž. Zgorelec

<p>In order to mitigate climate change and reduce the anthropogenic greenhouse gas (GHG) emissions, the Kyoto protocol has been adopted in 1997 and the Paris Agreement entered into force in 2016. The Paris Agreement have ratified 190 out of 197 Parties of the United Nations Framework Convention on Climate Change (UNFCCC) and Croatia is one of them as well. Each Party has obliged regularly to submit the national inventory report (NIR) providing the information on the national anthropogenic GHG emissions by sources and removals by sinks to the UNFCCC. Reporting under the NIR is divided into six categories / sectors, and one of them is land use, land use change and forestry (LULUCF) sector, where an issue of uncertainty estimates on carbon emissions and removals occurs. As soil respiration represents the second-largest terrestrial carbon flux, the national studies on soil respiration can reduce the uncertainty and improve the estimation of country-level carbon fluxes. Due to the omission of national data, the members of the University of Zagreb Faculty of Agriculture, Department of General Agronomy have started to study soil respiration rates in 2012, and since then many different studies on soil respiration under different agricultural land uses (i.e. annual crops, energy crop and vineyard), management practices (i.e. tillage and fertilization) and climate conditions (i.e. continental and mediterranean) in Croatia have been conducted. The obtained site specific results on field measurements of soil carbon dioxide concentrations by <em>in situ</em> closed static chamber method will be presented in this paper.</p>

S. Štrbac, M. Kašanin-Grubin, G. Veselinović, G. Gajica, S. Stojadinovic, A. Šajnović, D. Dimovic

<p>Human activities have changed ecosystems and today &#8776; 60% of the world&#8217;s ecosystems are already degraded. These changes have caused growing environmental costs, including biodiversity loss and land degradation, which in turn has resulted in many economic, social and cultural losses. Protected areas (PAs) are the key tool in biodiversity conservation, moreover they may help to maintain water supplies and food security, strengthen climate resilience and improve human health and well-being. International Union for Conservation of Nature (IUCN) defined PA as &#8222;a clearly defined geographical space, recognized, dedicated, and managed, through legal or other effective means, to achieve the long-term conservation of nature with associated ecosystem services (ES) and cultural values&#8221;. Such areas represent Earth systems in which influence of human interactions with preserved ecosystems are readily evident. The coverage of PA is a widely used indicator of sustainable development, because the loss of biodiversity is recognized as one of the most serious global environmental threats. The &#8220;Big Five&#8221; threats to global biodiversity are fragmentation, habitat loss, overexploitation of natural resources, pollution, and the spread of invasive alien species. New interventions for governing nature are captured by the umbrella of nature-based solutions (NBS) in the European Union (EU) policy context. NBS can offer accessible, sustainable, and feasible benefits via a range of areas affecting public health and social well-being. According to IUCN NBS are defined as &#8220;actions to protect, sustainably manage, and restore natural or modified ecosystems, that address societal challenges effectively and adaptively, simultaneously providing human well-being and biodiversity benefits&#8221;. NBS address these societal challenges over the delivery of ES. The main objective of this study is to use the effect of NBS to enhance the sustainability of management of the PAs that would have environmental, social and economic benefits. The methodology includes determination of heavy metals in soils and needles of Picea alba, and quantification and qualification of PAs benefits based on Protected Areas Benefits Assessment Tool + (PA-BAT+) in six sites: Zlatibor, Golija, Tara, &#272;erdap, Stara planina, and Fru&#353;ka gora. Zlatibor, Golija, and Stara planina are protected as a Nature Park &#8211; protected areas of international, national, i.e., exceptional importance Category I (first) in accordance with the Law on Nature Protection ("Off. Gazette of RS", No. 36/2009, 88/2010 , 91/2010 and 14/2016). By the decision of the UNESCO commission within the MAB program in 2001, Golija was declared as Biosphere Reserve &#8221;Golija - Studenica&#8221;. Tara, &#272;erdap, and Fru&#353;ka gora are protected as National Parks &#8211; protected area of international, national, i.e., exceptional importance Category I (first) in accordance with the Law on National Parks ("Off. Gazette of RS", No. 39/1993, 44/1993-correction, 53/1993, 67/1993, 48/1994, 101/2005 and 36/2009). According to categorization of the IUCN Zlatibor, Golija, and Stara planina are classified in Category V, while Tara, &#272;erdap, and Fru&#353;ka gora are classified in Category II. Based on heavy metals content in soils and needles, different interventions in managed ecosystems are proposed.</p>

M. Drilleau, R. Garcia, H. Samuel, A. Rivoldini, M. Wieczorek, P. Lognonné, M. Panning, C. Perrin et al.

<p>The InSight (Interior Exploration using Seismic Investigations, Geodesy and Heat Transport) lander successfully delivered a geophysical instrument package to the Martian surface on November 26, 2018, including a broadband seismometer called SEIS (Seismic Experiment for Interior Structure). After two years of recording, seismic body waves phases of a small number of high-quality marsquakes have been clearly identified. In this work, we will present how we estimate the body waves arrival times, and how we handle them to constrain the locations of the marsquakes and the interior structure. The inverse problem relies on a Bayesian approach, to investigate a large range of possible locations and interior models. Due to the small number of data, the advantage of using such a method is to provide a quantitative measure of the uncertainties and the non-uniqueness. In order to take into account the strong variations of the crustal thickness due to the crustal dichotomy, and thus consider the seismic lateral variations, which could cause significant misinterpretations, arrival times corrections are added using crustal thickness maps obtained from gravity and topography data.</p><p>&#160;</p>

N. Compaire, L. Margerin, R. Garcia, M. Calvet, B. Pinot, G. Orhand-Mainsant, Doyeon Kim, V. Lekić et al.

<p>Since early February 2019, the SEIS seismometer deployed at the surface of Mars in the framework of the NASA-InSight mission has been continuously recording the ground motion at Elysium Planitia. In this work, we take advantage of this exceptional dataset to put constraints on the crustal properties of Mars using seismic interferometry (SI). This method use the seismic waves, either from background vibrations of the planet or from quakes, that are scattered in the medium in order to recover the ground response between two seismic sensors. Applying the principles of SI to the single-station configuration of SEIS, we compute, for each Sol (martian day) and each local hour, all the components of the time-domain autocorrelation tensor of random ambient vibrations in various frequency bands. A similar computation is performed on the diffuse waveforms generated by more than a hundred Marsquakes. For imaging application a careful signal-to-noise ratio analysis and an inter-comparison between the two datasets are applied. These analyses suggest that the reconstructed ground responses are most reliable in a relatively narrow frequency band around 2.4Hz, where an amplification of both ambient vibrations and seismic events is observed. The average Auto-Correlation Functions (ACFs) from both ambient vibrations and seismic events contain well identifiable seismic arrivals, that are very consistent between the two datasets. We interpret the vertical and horizontal ACFs as the ground reflection response below InSight for the compressional waves and the shear waves respectively. We propose a simple stratified velocity model of the crust, which is most compatible with the arrival times of the detected phases, as well as with previous seismological studies of the SEIS record. The hourly computation of the ACFs over one martian year also allows us to study the diurnal and seasonal variations of the reconstructed ground response with a technique call Passive Image Interferometry (PII). In this study we present measurements of the relative stretching coefficient between consecutive ACF waveforms and discuss the potential origins of the observed temporal variations.</p>

Nema pronađenih rezultata, molimo da izmjenite uslove pretrage i pokušate ponovo!

Pretplatite se na novosti o BH Akademskom Imeniku

Ova stranica koristi kolačiće da bi vam pružila najbolje iskustvo

Saznaj više