Phylogeographic Substructuring in the Southernmost Refugium of the European Common Frog Rana temporaria.

Rana temporaria is one of the most widespread Palearctic brown frogs. We aimed to clarify distribution pattern of two main genetic clades in the understudied Balkan peninsula by using 16SrRNA and MT-CYTB sequences, already widely applied in analyses of populations from other parts of Europe, while focusing on the broad area along the Morava river (central Balkans) as a known gap in the species distribution. Additionally, we were interested in revealing the extent of haplotype diversity within the main genetic clades in the Balkans, particularly around the supposed suture zone. The results revealed a suture zone between the Western and Eastern Clades in the central part of the Balkan Peninsula. This indicated the existence of a historical barrier between the Balkan Mountain Belt and geographically close mountains surrounding the Vlasina Plateau (Rhodope/Serbian-Macedonian Massif). The overall observed haplotype diversity in populations of R. temporaria from the Balkan Peninsula seems high. Harboring both main genetic clades of R. temporaria qualifies the Balkan Peninsula as another important center of species' genetic diversity, as well as rich in unique haplotypes. This points out the necessity of applying conservation measures focused on the common European frog populations and habitats in this part of the species' distribution area.

A Cryptic Subterranean Mammal Species, the Lesser Blind Mole Rat (Nannospalax leucodon syrmiensis)-Retreated but Not Extinct.

Blind mole rats (genus Nannospalax) attract a great deal of attention because of their cancer resistance and longevity. Due to the high rate of chromosome rearrangements, 74 Nannospalax chromosomal forms have been discovered. The convergence of their external morphology complicates their taxonomy, and many cryptic species remain unrecognized. Thus, the European N. leucodon supersp. is listed in the IUCN Red List of Threatened Species with "Data Deficient" status. It is crucial for the conservation of biodiversity to clarify its taxonomy, to recognize each cryptic species, and assign to them the correct conservation status. Of the more than 20 chromosomal forms described within N. leucodon, five cryptic species occur in Serbia. The most threatened among them-N. l. syrmiensis, described and named 50 years ago in the regions of Srem, Belgrade and Macva-has been declared extinct in the literature, which may have negative consequences for the conservation of wildlife genetic diversity. Through five years of fieldwork and comparison of 16SrRNA and MT-CYTB gene segments between old, archived teeth and recently collected material, we show that N. l. syrmiensis is not extinct. However, its habitat has been fragmented and reduced, owing primarily to anthropogenic impact. Therefore, detailed surveillance, population-structure studies, risk assessment, and appropriate conservation measures are needed.

Evolutionary dynamics of Usutu virus: Worldwide dispersal patterns and transmission dynamics in Europe.

Background: Usutu virus (USUV) is an emerging mosquito-borne Flavivirus, with birds as the main zoonotic reservoir. Humans are accidental hosts and mostly develop mild or even asymptomatic infections, although severe complications such as encephalitis can also arise. Detailed characterization of the pathogen's phylogenetics may offer valuable insights into the prediction and prevention of potential epidemics; however, lack of uniformity and the number of available USUV sequences worldwide hamper comprehensive investigation. Aim: The study aimed to investigate USUV spatio-temporal dispersal inter- and intracontinentally and to estimate the dynamics of viral spread within Europe. Methods: Phylogeographic and phylodynamic analyses were done using advanced phylogenetic methods implemented in Beast 1.10.4 and Beast 2.6.4 software packages. Results: Herein, we report on a new USUV isolate from Culex pipiens collected in 2019 from Serbia. The results of this research revealed two newly described intercontinental migration events of USUV from Africa to Germany in the 1970s and from Africa to the Middle East (Israel) in the late 90s. Finally, phylodynamic analysis substantiated the ongoing active expansion of USUV in Europe. Conclusion: The data would imply a high potential for further USUV expansion in Europe. Detailed phylogenetic characterization of the pathogen may offer valuable insights into prediction and prevention of potential epidemics; however, lack of uniformity and number of available USUV sequences worldwide hampers comprehensive investigation. This study draws attention to the need for upscaling USUV surveillance.

Detection of the Xanthi Chryso-like Virus in New Geographical Area and a Novel Arthropod Carrier.

Here, we report on a serendipitous finding of a chryso-like virus associated with Culex pipiens mosquitos in the course of study aimed to detect and characterize West Nile virus (WNV) circulating in mosquitos in Serbia, Southern Europe. Upon initial detection of unexpected product in a PCR protocol for partial WNV NS5 gene amplification, further confirmation and identification was obtained through additional PCR and Sanger sequencing experiments. Bioinformatic and phylogenetic analysis identified the obtained sequences as Xanthi chryso-like virus (XCLV). The finding is particular for the fact that it associates XCLV with a new potential vector species and documents a novel geographical area of its distribution.

Cryptic Diversity of the European Blind Mole Rat Nannospalax leucodon Species Complex: Implications for Conservation.

We explored the cryptic speciation of the Nannospalax leucodon species complex, characterised by intense karyotype evolution and reduced phenotypic variability that has produced different lineages, out of which 25 are described as chromosomal forms (CFs), so many cryptic species remain unnoticed. Although some of them should be classified as threatened, they lack the official nomenclature necessary to be involved in conservation strategies. Reproductive isolation between seven CFs has previously been demonstrated. To investigate the amount and dynamics of genetic discrepancy that follows chromosomal changes, infer speciation levels, and obtain phylogenetic patterns, we analysed mitochondrial 16S rRNA and MT-CYTB nucleotide polymorphism among 17 CFs-the highest number studied so far. Phylogenetic trees delineated 11 CFs as separate clades. Evolutionary divergence values overlapped with acknowledged higher taxonomic categories, or sometimes exceeded them. The fact that CFs with higher 2n are evolutionary older corresponds to the fusion hypothesis of Nannospalax karyotype evolution. To participate in conservation strategies, N. leucodon classification should follow the biological species concept, and proposed cryptic species should be formally named, despite a lack of classical morphometric discrepancy. We draw attention towards the syrmiensis and montanosyrmiensis CFs, estimated to be endangered/critically endangered, and emphasise the need for detailed monitoring and population survey for other cryptic species.

Phylogeographic analysis of Tula hantavirus highlights a single introduction to central Europe.

Orthohantaviruses are zoonotic pathogens of humans, unique among the bunyaviruses in not being transmitted by an arthropod vector. Tula orthohantavirus (TULV) is an old-world hantavirus, of yet unclear human pathogenicity, with few reported cases of clinically relevant human infection. So far, phylogeographic studies exploring the global pathways of hantaviral migration are scarce and generally do not focus on a specific hantavirus species. The aim of the present study was to reconstruct the dispersal history of TULV lineages across Eurasia based on S segment sequences sampled from different geographic areas. Maximum-likelihood and Bayesian inference methods were used to perform the phylogenetic analysis and phylogeographic reconstructions. Sampling time and trapping localities were obtained for a total of 735 TULV S segment sequences available in public databases at the time of the study. The estimated substitution rate of the analyzed partial S segment alignment was 2.26 × 10-3 substitutions/site/year (95 per cent highest posterior density interval: 1.79 × 10-3 to 2.75 × 10-3). Continuous phylogeography of TULV S segment sequences placed the potential root and origin of TULV spread in the Black Sea region. In our study, we detect a single-lineage introduction of TULV to Europe, followed by local viral circulation further on.

Depicting the RNA Virome of Hematophagous Arthropods from Belgrade, Serbia.

Hematophagous arthropods are important vectors for zoonotic pathogens. To date, a huge number of viruses have been identified in these arthropods, with a considerable proportion of them being human pathogens. However, the viromes of hematophagous arthropods are still largely unresearched. In this study, a number of arthropods were collected from Belgrade, Serbia including mosquitoes, ticks and bedbugs. The viromes of these arthropods were identified and characterized using Illumina MiSeq sequencing. In total, 21 viruses belonging to 11 families were characterized, with 11 of them representing novel species. These results may contribute to our knowledge of RNA viruses in arthropods and the discovery of novel human pathogens.

Polymorphisms and haplotypes in VDR gene are associated with female idiopathic infertility.

Unexplained infertility refers to the absence of a definable cause of reproductive failure. Vitamin D receptor (VDR) acts as a transcription factor and regulates a number of vitamin D-responsive genes, including those involved in the immune system. Recent finding that VDR is expressed in reproductive tissues suggests a possible importance of vitamin D in pregnancy. We conducted a case-control study to examine the association of polymorphisms in VDR gene with reproductive success. DNA from 117 female patients with unexplained infertility and 130 fertile controls was isolated from peripheral blood and VDR genotypes (FokI, BsmI, ApaI and TaqI) were detected by PCR-RFLP. Haplotypes were determined using Haploview software. Our results show significant association of FokI and BsmI polymorphisms with infertility (p < 0.05). The haplotype analysis confirmed strong linkage disequilibrium between closely positioned BsmI, ApaI and TaqI polymorphisms. Two haplotypes were associated with infertility: (i) haplotype bAT was increasing the risk for secondary infertility; while (ii) haplotype BAT had a protective role against primary infertility (p < 0.05). By changing the expression and the activity of VDR gene, which leads to the change in expression of vitamin D-responsive genes, these polymorphisms and haplotypes could possibly have an effect on immune system in the female reproductive tract.

Insight into diversity of bacteria belonging to the order Rickettsiales in 9 arthropods species collected in Serbia.

Rickettsiales bacteria in arthropods play a significant role in both public health and arthropod ecology. However, the extensive genetic diversity of Rickettsiales endosymbionts of arthropods is still to be discovered. In 2016, 515 arthropods belonging to 9 species of four classes (Insecta, Chilopoda, Diplopoda and Arachnida) were collected in Serbia. The presence and genetic diversity of Rickettsiales bacteria were evaluated by characterizing the 16S rRNA (rrs), citrate synthase (gltA) and heat shock protein (groEL) genes. The presence of various Rickettsiales bacteria was identified in the majority of tested arthropod species. The results revealed co-circulation of five recognized Rickettsiales species including Rickettsia, Ehrlichia and Wolbachia, as well as four tentative novel species, including one tentative novel genus named Neowolbachia. These results suggest the remarkable genetic diversity of Rickettsiales bacteria in certain arthropod species in this region. Furthermore, the high prevalence of spotted fever group Rickettsia in Ixodes ricinus ticks highlights the potential public health risk of human Rickettsia infection.

Hemostasis-related gene polymorphisms and their epistatic relationship in women with idiopathic infertility.

: A numerous factor can cause infertility, but around one of four reproductive failure cases remain unexplained and diagnosed as idiopathic infertility. In the past few decades, analysis of gene polymorphisms takes a significant place in pathogenesis of infertility. The aim of this study was to evaluate the possible role of hemostasis-related gene polymorphisms in unexplained infertility. The study includes 117 female patients with idiopathic infertility and 130 fertile women with at least one born child. Eight polymorphisms important for hemostasis (ITGB3 1565T>C, FV 1691G>A, FII 20210G>A, MTHFR 677C>T and 1298A>C, ATIII 786G>A, PAI-14G/5G and ACE I/D) were genotyped by real-time PCR system. The frequencies of alleles and genotypes of examined polymorphisms were analyzed in SPSS statistical program, whereas gene interactions were identified using the GMDR software. Examination of etiological factors has shown that family history is a significant factor in assessing individual risk for infertility. The alleles and genotypes frequency of FV 1691G>A and FII 20210G>A polymorphisms were statistically different between control and patient group leading to a greater risk for infertility. The analysis of epistatic relationship between examined hemostasis-related gene polymorphisms identified more complex high-risk genotypes associated with infertility. Our results suggest that positive family history could be important predictive factor for fertility problems, pointing to the potential hereditary basis of this condition. Polymorphisms FVL and FII prothrombin are independent risk factors for idiopathic infertility, whereas multilocus interactions approach should be taken into consideration for the future research.

Molecular characterization of Dobrava-Belgrade hantavirus in Serbia, 2007-2011.

BACKGROUND: Hantaviruses are etiological agents of emerging zoonotic diseases worldwide, including hemorrhagic fever with renal syndrome (HFRS). A number of hantavirus species is known to be present in Europe. In Serbia, existing data on hantavirus presence and prevalence rely in serological findings. In this study, molecular analysis was performed in order to characterize HFRS causing hantaviruses in Serbia. METHODS: Sixty four serum samples of HFRS cases, previously found seropositive to anti-hantaviral antibodies, were included in the study. Partial hantaviral L and S segments were PCR amplified producing 390nt and 598nt amplicons, respectively, in parallel with human beta-actin mRNA as external reverse transcription positive control. Hantavirus specific PCR products were DNA sequenced in both direction and the obtained sequences phylogenetically confirmed and analyzed. RESULTS: PCR detection of hantavirus L and S genome segments was positive in 18/64 and 11/64 tested samples, respectively. Positive PCR results involved samples obtained from different locations, mostly from central and southern parts of Serbia. All the obtained sequences were identified as Dobrava-Belgrade virus (DOBV). In the phylogenetic analysis sequences from Serbia tended to cluster in distinctive, geographically related clusters. CONCLUSIONS: Our findings indicate DOBV as the main HFRS causing hantavirus in Serbia, the site of its initial isolation.

Y chromosome genetic data defined by 23 short tandem repeats in a Serbian population on the Balkan Peninsula.

BACKGROUND: Serbs mainly live in the territory of the recently re-established state of Serbia. However, the turbulent history in the Balkan Peninsula has led to settlement of Serbs not only within present day Serbia, but also in different parts of neighbouring countries. AIM: To define polymorphisms of 23 Y-chromosomal short tandem repeat (STR) loci in a modern Serbian population from the central part of the Balkan Peninsula. SUBJECTS AND METHODS: The reference sample consisted of 303 men declared as Serbs over three generations. Localities of the collected materials include the territories of Serbia, Bosnia and Herzegovina, Croatia and Montenegro. DNA samples were typed using the PowerPlex®Y23 amplification kit. RESULTS: The highest locus diversity was observed for DYS385 and DYS481. In this study the most abundant haplogroups were I2a, E1b1b, R1a and I1. The largest genetic distances between the Serbs and other close Southern Slavs were for the Macedonians and Slovenians. CONCLUSION: This study is the first one to define STR polymorphism of Serbian people not only from Serbia but also from other parts of the Balkan Peninsula. The presented genetic data may be useful in further examinations of the genesis and genetic structuring of the present-day Serbian gene pool.

Genetics of suspected thrombophilia in Serbian females with infertility, including three cases, homozygous for FII 20210A or FV 1691A mutations.

Reproductive failure (recurrent foetal loss, unexplained infertility and IVF implantation failure) may be, in a number of cases, explained by thrombophilia, either acquired or inherited. Several genes contribute to thrombophilia, some with major effect (Factor V, Factor II), and some with minor effect (MTHFR, PAI-1, ATIII, etc.). The aim of this study was to estimate frequency of thrombophilia-associated genotypes (FII20210G > A, FV1691G > A, MTHFR677C > T and PAI-1 -675 4G/5G) in a group of 1631 Serbian women experiencing reproductive failure, and compare it with a healthy, female control group. Our results showed marginally significant (p = 0.050) differences in allele frequencies between patients and controls for the FV1691 mutations. For the FII20210G > A, although the statistical significance was not achieved (p = 0.076), we found higher frequency of variant allele in patients compared to controls (1.87% vs. 0.38%, respectively) which may point to a possible role of this polymorphism in thrombotic events. For the MTHFR677C > T and PAI-1 -675 4G/5G, we found no difference in distributions of genotype or allele frequencies between these two groups (p > 0.05). For three subjects with very rare genotypes (two patients homozygous for FV1691G > A and one patient homozygous for FII20210G > A) we performed additional biochemical analyses for haemostasis, as well as genotyping of two polymorphisms (MTHFR1298A > C and ATIII786G > A).

Substitution rate and natural selection in parvovirus B19.

The aim of this study was to estimate substitution rate and imprints of natural selection on parvovirus B19 genotype 1. Studied datasets included 137 near complete coding B19 genomes (positions 665 to 4851) for phylogenetic and substitution rate analysis and 146 and 214 partial genomes for selection analyses in open reading frames ORF1 and ORF2, respectively, collected 1973-2012 and including 9 newly sequenced isolates from Serbia. Phylogenetic clustering assigned majority of studied isolates to G1A. Nucleotide substitution rate for total coding DNA was 1.03 (0.6-1.27) x 10-4 substitutions/site/year, with higher values for analyzed genome partitions. In spite of the highest evolutionary rate, VP2 codons were found to be under purifying selection with rare episodic positive selection, whereas codons under diversifying selection were found in the unique part of VP1, known to contain B19 immune epitopes important in persistent infection. Analyses of overlapping gene regions identified nucleotide positions under opposite selective pressure in different ORFs, suggesting complex evolutionary mechanisms of nucleotide changes in B19 viral genomes.

Exploring Supernumeraries - A New Marker for Screening of B-Chromosomes Presence in the Yellow Necked Mouse Apodemus flavicollis.

Since the density of simple sequence repeats (SSRs) may vary between different chromosomes of the same species in eukaryotic genomes, we screened SSRs of the whole genome of the yellow necked mouse, Apodemus flavicollis, in order to reveal SSR profiles specific for animals carrying B chromosomes. We found that the 2200 bp band was amplified by primer (CAG)4AC to a highly increased level in samples with B chromosomes. This quantitative difference (B-marker) between animals with (+B) and without (0B) B chromosomes was used to screen 20 populations (387 animals). The presence/absence of Bs was confirmed in 96.5% of 342 non mosaic individuals, which recommends this method for noninvasive B-presence detection. A group of 45 animals with mosaic and micro B (µB) karyotypes was considered separately and showed 55.6% of overall congruence between karyotyping and molecular screening results. Relative quantification by qPCR of two different targeted sequences from B-marker indicated that these B-specific fragments are multiplied on B chromosomes. It also confirms our assumption that different types of Bs with variable molecular composition may exist in the same individual and between individuals of this species. Our results substantiate the origin of Bs from the standard chromosomal complement. The B-marker showed 98% sequence identity with the serine/threonine protein kinase VRK1 gene, similarly to findings reported for Bs from phylogenetically highly distant mammalian species. Evolutionarily conserved protein-coding genes found in Bs, including this one in A. flavicollis, could suggest a common evolutionary pathway.

Evidence of recombination in Tula virus strains from Serbia.

Tula hantavirus (TULV) belongs to Bunyaviridae family, with negative sense RNA genome. Segmented nature of hantaviral genome allows for genetic reassortment, but the evidence of homologous recombination also exists. In this study we analyzed TULV sequences isolated in Serbia on different occasions and from different rodent hosts: 1987 strain from Microtus subterraneus and 2007 strain from Microtus arvalis. Phylogenetic analysis of both L and S segment sequences is suggestive of geographically related clustering, as previously shown for majority of hantaviruses. Reconstruction of phylogenetic tree for TULV S segment showed that both sequences from Serbia clustered together with sequences from East Slovakia, which had previously been shown to be recombinants (Kosice strain). Exploratory recombination analysis, supported by phylogenetic and amino acid pattern analysis, revealed the presence of recombination in the S segment sequences from Serbia, resulting in mosaic-like structure of TULV S segment similar to the one of Kosice strain. Although recombination is considered a rare event in molecular evolution of negative strand RNA viruses, obtained molecular data in this study support evidence of recombination in TULV, in geographically distant regions of Europe.

Molecular characterization and phylogenetic analysis of full-genome HBV subgenotype D3 sequences from Serbia.

Hepatitis B virus (HBV) is classified into 8 genotypes with distinct geographical distribution. Genotype D (HBV/D) has the widest distribution area and is comprised of 7 subgenotypes. Subgenotypes D1, D2 and D3 appear worldwide, while D4-D7 have a more restricted distribution. Within the Mediterranean area, HBV/D and subgenotype D3 are the most prevalent. The purpose of this study was to characterize the full genome of Serbian HBV/D3 isolates by comparison and phylogenetic analysis with HBV/D3 sequences (66 samples) found in GeneBank/DDBJ databases from different parts of the world. Isolates were obtained from three patients diagnosed with chronic hepatitis B (HBsAg+). All three isolates have two very rare nucleotide substitutions, A929T and T150A, which indicate the same ancestor. Phylogenetic analysis of HBV/D3 genome sequences throughout the world follows an ethno-geographical origin of isolates with rare exceptions, which could be explained by human travelling and migration. The geographically close but ethnically different Serbian and Italian isolates clustered in the same subnode, and on a common branch with strains from Northern Canada. To test the apparently close HBV phylogenetic relationship between completely separated patients from Serbia and Northern Canada we analyzed in depth a 440 bp region of the HBsAg from Canadian (n=73) and Serbian (n=70) isolates. The constructed parsimony tree revealed that strains from Serbia and Northern Canada fell along the same branch which indicates independent evolution within regions of each country. Considering that HBsAg sequence has limited variability for phylogenetic analyses, our hypothesis needs further confirmation with more HBV complete genome sequences.

Potential genotoxic effects of melted snow from an urban area revealed by the Allium cepa test.

The presence of well-known atmospheric pollutants is regularly screened for in large towns but knowledge about the effects of mixtures of different pollutants and especially their genotoxic potential is largely missing. Since falling snow collects pollutants from the air, melted snow samples could be suitable for evaluating potential genotoxicity. For this purpose the Allium cepa anaphase-telophase test was used to analyse melted snow samples from Belgrade, the capital city of Serbia. Samples of snow were taken at two sites, characterized by differences in pollution intensity, in three successive years. At the more polluted site the analyses showed a very high degree of both toxicity and genotoxicity in the first year of the study corresponding to the effects of the known mutagen used as the positive control. At the other site the situation was much better but not without warning signals. The results showed that standard analyses for the presence of certain contaminants in the air do not give an accurate picture of the possible consequences of urban air pollution because the genotoxic potential remains hidden. The A. cepa test has been demonstrated to be very convenient for evaluation of air pollution through analyses of melted snow samples.

[Role of various virus genotypes in progression of chronic hepatitis C]. / Uloga razlicitih genotipova virusa u progresiji hronicnog hepatitisa C.

The relation between HCV genotypes and the progression of chronic hepatitis is still unknown. Some studies implied more pathogenic effect of genotype 1b for the severity of liver inflammation. However, other studies did not show the association between HCV genotype 1b and the severe outcome of HCV infection. The aim of this study was to determine the most frequent genotypes in this environment and their influence on hepatitis C severity. The investigation included 34 patients with histologically confirmed chronic hepatitis C, aged 20-65 (mean 35.0 years). On the basis of pathohistological findings, applying the modern classification, the disease activity was graded as: minimal (A1) moderate (A2) and severe (A3). The extent of fibrosis was marked as: absent (F0), mild (F1), moderate (F2) and severe (F3). Genotyping was performed by nested PCR with type-specific primers and LIPA test and verified by sequencing. The most prevailing genotype in our group of patients was 1b (44.1%), followed by genotype 3a (26.4%), genotype 1a (11.7%) and 2a (2.8%). Five patients had mixed genotypes (four 3a/1b, and one 1a/1b). The severity of liver cell necrosis, measured by alanintransferase (ALT) levels in serum was not related to any of HCV genotypes. There was no statistically significant difference between histological disease activity in relation to HCV genotypes. Stage of the disease was not significantly related to the HCV genotypes. There was a strong association between the degree of fibrosis and the age of patients (p < 0.01). These results could indicate that the determination of HCV genotypes was not useful in the estimation of disease severity and that liver biopsy was the most important for the prognosis of the disease.