Genetic variability of Roma population in Serbia: The perspective from autosomal STR markers.

Genetic variability of Roma population was shaped by the strong influence of genetic drift and gene flow during the migrations from their ancestral homeland in Indian subcontinent towards Europe. In addition, social stigmatization in many European countries, as a consequence of different cultural heritage and social practices, induced further genetic differentiation and sub structuring within the population. Although many populations genetic studies on European Roma were carried out, the genetic structure of the Serbian Roma has not been described yet, since only the modest number of individuals from this territory was analyzed. The main aim of this study was the characterization of genetic variability of the Roma and the assessment of intrapopulation genetic differentiation based on the analysis of 21 autosomal STR loci of 259 self-identified unrelated individuals from Serbia. Intrapopulation analysis revealed divergence of Roma groups illustrating the effect of the historical events after their arrival on Balkan Peninsula and emphasizing significance of the religious affiliation on admixture with autochthonous population. Genetic distance analysis showed the greatest similarity of the studied population with the Middle Eastern populations, while South Asian and European population were more distant. Our results demonstrate that Roma groups in this region of Balkan Peninsula do not represent completely isolated, but rather admixed populations with different proportion of gene flow with other Roma and non-Roma groups.

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

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

Distribution of Y-chromosome haplogroups in Serbian population groups originating from historically and geographically significant distinct parts of the Balkan Peninsula.

Our study enrolled 1200 Serbian males originating from three geographical regions in the Balkan Peninsula inhabited by Serbs: present-day Serbia, regions of Old Herzegovina and Kosovo and Metohija. These samples were genotyped using the combination of 23 Y-chromosomal short tandem repeats (Y-STRs) loci and 17 Ychromosomal single nucleotide polymorphisms (Y-SNPs) loci for the haplotype and haplogroup analysis in order to characterize in detail Y chromosome flow in the recent history. Serbia's borders have changed through history, forcing Serbs constantly to migrate to different regions of Balkan Peninsula. The most significant migration waves in the recent history towards present-day Serbia occurred from the regions of Old- Herzegovina and Kosovo and Metohija that lie in the south-west/south. High haplotype diversity and discrimination capacity were observed in all three datasets, with the highest number of unique haplotypes (381) and discrimination capacity (0.97) detected in the samples originating from the present-day Serbia. Haplogroup composition didn't differ significantly among datasets, with three dominant haplogroups (I-M170, E-P170 and R-M198), and haplogroup I-M170 being the most frequent in all three datasets. Haplogroup E-P170 was the second most dominant in the dataset originating from geographical region of Kosovo and Metohija, whereas haplogroup R-M198 was the second most prevalent in the dataset from historical region of Old Herzegovina. Based on the phylogenetic three for haplogroup I constructed within this study, haplogroup I2a1-P37.2 was the most dominant within all three datasets, especially in the dataset from historical region of Old Herzegovina, where 182 out of 400 samples were derived for SNP P37.2. Genetic distances between three groups of samples, evaluated by the Fst and Rst statistical values, and further visualized through multidimensional scaling plot, showed great genetic similarity between datasets from Old Herzegovina and present-day Serbia. Genetic difference in the haplogroup distribution and frequency between datasets from historical region of Old Herzegovina and from geographical region of Kosovo and Metohija was confirmed with highest Fst and Rst vaules. In this study we have distinguished genetic structure, diversity and haplogroup frequencies within 1200 Serbian males from three datasets, relationships among them as well as with other Balkan and European populations, which is useful for studying recent demographic history.

A comprehensive mutation study in wide deep-rooted R1b Serbian pedigree: mutation rates and male relative differentiation capacity of 36 Y-STR markers.

Haplotyping of Y-chromosomal short tandem repeats (Y-STRs) reflects the paternal lineage, although, the father-son pair profiles may differ due to the germline mutations. In order to discriminate between closely related males in criminal cases, as well as for the correct application of Y-STRs in the paternity/kinship analysis and determination of the most recent common ancestor in the familial searching or genealogy research, the assessment of mutation rates of routinely used Y-STRs is of a great importance. We genotyped 120 males belonging to one wide deep-rooted pedigree separated by 1-20 meiosis. The haplotypes of analyzed males distributed over 12 different families (according to their surnames), with 113 originating from one ancestor, and the remaining 7 from the second, closely related to the previous one, belong to the R1b haplogroup. The analysis was performed using Powerplex® Y23 kit, Yfiler™ plus kit and 13 rapidly mutating (RM13) Y-STRs. In 20,855 allele transmissions, 175 mutations (61% repeat losses and 39% gains) and one gene conversion event were found at 25 out of 36 markers. The medians of locus-specific mutation rates estimated using the Bayesian approach ranged from 1.42 × 10-3 (95% credible interval (CI): 0.05 × 10-3 - 7.56 × 10-3) for loci with no observed mutations to 130.91 × 10-3 (95% CI: 102.91 × 10-3 - 162.78 × 10-3) for DYF399S1, with a median rate across all 36 markers of 10.06 × 10-3 (95% CI: 8.65 × 10-3 - 11.61 × 10-3). In 6349 male relative pairs, the 36 Y-STR set distinguished 98.4% relative pairs by at least one mutation, compared to 95.9%, 65.5% and 57.4% for RM13, Yfiler™ plus, and Powerplex® Y23 set, respectively. The extra-pair paternity rate was estimated at 11.9 × 10-3 (95% CI: 4.4 × 10-3 - 25.8 × 10-3) fitting within the range reported for some European populations. A significant positive correlation was observed between fathers' ages at the time of the Y chromosome transmission and mutability rates (R2 = 0.9495, p = 0.0256), with more significant results when analyzing RM markers (R2 = 0.9827, p = 0.0087).

Identification of a broad spectrum of mammalian and avian species using the short fragment of the mitochondrially encoded cytochrome b gene.

Mitochondrial DNA (mtDNA), especially the gene for cytochrome b (MT-CYB), has been found to be highly informative for species identification. In this study, we present the results of the analysis of a 127 bp long fragment of MT-CYB, amplified using universal primers, variable enough to be used for species identification and discrimination, even in highly degraded animal samples. The total number of analyzed species in this study was 30, including 17 mammalian and 13 bird species. Using a newly created primer pair, we successfully amplified and sequenced the target sequence in almost all tested species. The amplification was incomplete in just two species, and as a result, partial, but still variable sequences, were obtained. Using the target fragment we successfully identified all tested samples. Initial results suggested that the intraspecies genetic diversity of the target region, in all tested species, was low - from 0 to 4.72%. The interspecies genetic diversity of the target region, crucial for successful discrimination, showed relatively high diversity, ranging from 8.36% to 42.52%. Given its short length, the target region should be used for species determination, particularly in samples that are degraded or are low in DNA quantity.

Comparison of different methods of DNA recovery and PCR amplification in STR profiling of casings-a retrospective study.

Casings represent common evidence in a forensic laboratory, due to high frequency of firearms usage during perpetration of criminal offenses. Possible DNA evidence from casings is compromised by degradation, inhibition, and initial low-quantity deposition of biological material. For that reason, in the last 15 years, scientists have been trying to optimize procedures for recovery and amplification of DNA possibly present on its surface. In this study, we share our 12-year experience done on a total of 698 casework casings, comparing two DNA recovery methods commonly used-soaking and swabbing, as well as efficacy of two commercially available DNA amplification kits (AmpFLSTR® Identifiler® and AmpFLSTR® Identifiler® Plus kits). Of all analyzed casings, 30 were excluded as 28 (4%) matched the victims' DNA profiles and 2 (0.3%) samples were proved to be contaminated by technicians. Overall success in obtaining interpretable DNA profiles was 15.6% (104/668) (13.8% (55/399) for AmpFLSTR® Identifiler® Plus combined with soaking, 22% (33/150) for AmpFLSTR® Identifiler® Plus combined with swabbing, and 13.4% (16/119) using AmpFLSTR® Identifiler® kit and swabbing recovery method). Our data suggest the importance of both DNA recovery methods and amplification kits used, and point out swabbing of casings combined with AmpFLSTR® Identifiler® Plus kit as methods of choice. Nonetheless, our results are based on real casework and are prone to uncontrolled variables.

Differentiation of Cannabis subspecies by THCA synthase gene analysis using RFLP.

Cannabis sativa subspecies, known as industrial hemp (C. sativa sativa) and marijuana (C. sativa indica) show no evident morphological distinctions, but they contain different levels of psychoactive Δ-9-tetrahidrocanabinol (THC), with considerably higher concentration in marijuana than in hemp. C. sativa subspecies differ in sequence of tetrahydrocannabinolic acid (THCA) synthase gene, responsible for THC production, and only one active copy of the gene, distinctive for marijuana, is capable of producing THC in concentration more then 0,3% in dried plants, usually punishable by the law. Twenty different samples of marijuana that contain THC in concentration more then 0,3% and three varieties of industrial hemp were analyzed for presence of an active copy of THCA synthase gene using in-house developed restriction fragment length polymorphism (RFLP) method All twenty samples of marijuana were positive for the active copy of THCA synthase gene, 16 of them heterozygous. All three varieties of industrial hemp were homozygous for inactive copy. An algorithm for the fast and accurate forensic analysis of samples suspected to be marijuana was constructed, answering the question if an analyzed sample is capable of producing THC in concentrations higher than 0.3%.

Genetics of Lafora progressive myoclonic epilepsy: current perspectives.

Lafora disease (LD) is a fatal neurodegenerative disorder caused by loss-of-function mutations in either laforin glycogen phosphatase gene (EPM2A) or malin E3 ubiquitin ligase gene (NHLRC1). LD is associated with gradual accumulation of Lafora bodies (LBs). LBs are aggregates of polyglucosan, a long, linear, poorly branched, hyperphosphorylated, insoluble form of glycogen. Loss-of-function mutations either in the EPM2A or in the NHLRC1 gene lead to polyglucosan formation. One hypothesis on LB formation is based on findings that laforin-malin complex downregulates glycogen synthase (GS) through malin-mediated ubiquitination, and the other one is based on findings that laforin dephosphorylates glycogen. According to the first hypothesis, polyglucosan formation is a result of increased GS activity, and according to the second, an increased glycogen phosphate leads to glycogen conformational change, unfolding, precipitation, and conversion to polyglucosan, while GS remains bound to the precipitating glycogen. In this review, we summarize all the recent findings that have important implications for the treatment of LD, all of them showing that partial inhibition of GS activity may be sufficient to prevent the progression of the disease. The current perspective in LD is high-throughput screening for small molecules that act on the disease pathway, that is, partial inhibitors of GS, which opens a therapeutic window for potential treatment of this fatal disease.

Analysis of PMP22 duplication and deletion using a panel of six dinucleotide tandem repeats.

BACKGROUND: Charcot-Marie-Tooth type 1A (CMT1A) is the most common type of hereditary motor and sensory neuropathies (HMSN), caused by the duplication of the 17p11.2 region that includes the PMP22 gene. Reciprocal deletion of the same region is the main cause of hereditary neuropathy with liability to pressure palsies (HNPP). CMT1A accounts for approximately 50% of HMSN patients. Diagnostics of CMT1A and HNPP are based on quantitative analysis of the affected region or RFLP detection of breakage points. The aim of this study was to improve the sensitivity and efficiency of CMT1A and HNPP genetic diagnostics by introducing analysis of six STR markers (D17S261-D17S122-D17S839-D17S1358-D17S955-D17S921) spanning the duplicated region. METHODS: Forty-six CMT1A and seven HNPP patients, all genetically diagnosed by RFLP analysis, were tested for duplication or deletion using six STR markers. RESULTS: In all CMT1A and HNPP patients, microsatellite analysis comprising six STR markers confirmed the existence of a duplication or deletion. In 89% (41/46) CMT1A patients the confirmation was based on detecting three alleles on at least one locus. In the remaining 11% (5) CMT1A patients, duplication was also confirmed based on two peaks with clear dosage difference for at least two different markers. All HNPP patients (7/7) displayed only one allele for each analyzed locus. CONCLUSIONS: Microsatellite analysis using six selected STR loci showed a high level of sensitivity and specificity for genetic diagnostics of CMT1A and HNPP. The results here strongly suggest STR marker analysis as a method of choice in PMP22 duplication/deletion testing.

A shared haplotype indicates a founder event in Unverricht-Lundborg disease patients from Serbia.

Unverricht-Lundborg disease (ULD) is an autosomal recessive disorder caused by dodecamer repeat expansion in the promoter region of the cystatin B (CSTB) gene in approximately 90% of the disease alleles worldwide. This study presents results of genetic findings in four Serbian unrelated patients with clinical and molecular diagnosis of ULD. Using newly established PCR protocol with betaine, we detected a homozygous expansion of dodecamer repeats in the CSTB gene in four patients with clinical diagnosis of ULD. Our results are in agreement with previous studies showing that dodecamer repeats expansion is the most common mutation associated with ULD. Haplotype analysis of eight unrelated ULD chromosomes was performed using seven markers flanking CSTB gene and one intragenic variant. We demonstrated the existence of a founder effect, strongly supported by LD calculations. Size of the minimal common haplotype implies that the most recent common ancestor of the Serbian ULD patients lived about 110 generations ago. We showed that Serbian ULD patients share the same common ancestor with patients from Baltic countries and North Africa. In the light of our data, we proposed extended minimal common haplotype, which could be considered as initial haplotype of the founder event common for Serbian, Baltic, and North African ULD patients.

Lafora disease: severe phenotype associated with homozygous deletion of the NHLRC1 gene.

Lafora disease (LD) is a severe, autosomal recessive, latechildhood- to teenage-onset, progressive myoclonic epilepsy. It is due to either EPM2A or NHLRC1 mutations. We describe a patient with homozygous deletion encompassing the entire NHLRC1 gene, not previously reported, and with clinical course more progressive than in the most patients with NHLRC1 mutations. The diagnosis of LD in our patient was based on the typical clinic, neurophysiological presentation, as well as skin biopsy followed by molecular genetics findings. She developed normally until the age of 15, when she had her first occipital and generalized seizures. Four years after the first seizure the patient became bedridden, demented and presented with severe clinical condition. She died of pneumonia at age 20. This report is the first case of homozygosity for NHLRC1 deletion and thus adds to mutational heterogeneity of LD. Besides, it widens the spectrum of LD patients with severe phenotype and NHLRC1 mutations.

An algorithm for genetic testing of Serbian patients with demyelinating Charcot-Marie-Tooth.

Charcot-Marie Tooth (CMT) is a clinically and genetically heterogeneous group of diseases with rough genotype-phenotype correlation, so the final diagnosis requires extensive clinical and electrophysiological examination, family data, and gene mutation analysis. Although there is a common pattern of genetic basis of CMT, there could be some population differences that should be taken into account to facilitate analyses. Here we present the algorithm for genetic testing in Serbian patients with demyelinating CMT, based on their genetic specificities: in cases of no PMP22 duplication, and if -X-linked CMT (CMTX) is not contraindicated by pattern of inheritance (male-to-male transmission), one should test for c.94A>G GJB founder mutation, first. Also, when a patient is of Romani ethnicity, or if there is an autosomal recessive inheritance in a family and unclear ethnicity, c.442C>T mutation in NDRG1 should be tested.

A novel P66S mutation in exon 3 of the SOD1 gene with early onset and rapid progression.

Amyotrophic lateral sclerosis (ALS) is a fatal motor neuron disease in adults of unknown origin in most cases. Here we report a novel P66S mutation in exon 3 of the SOD1 gene in an apparently sporadic ALS patient with unusual early onset and rapid disease progression. Our data widen the spectrum of SOD1 mutations and clinical presentations of ALS.

Polymorphisms of the prion protein gene (PRNP) in a Serbian population.

Prion diseases are a group of etiologically heterogeneous neurodegenerative disorders. We have analyzed the coding region of PRNP gene in 121 healthy citizens of Serbia to determine whether the frequencies of M129V, E219K, and octapeptide repeat number polymorphism. For Serbian population, polymorphism of PRNP gene at codon 129 does not differ from healthy European populations. Also codon 219 is monomorphic for the Glu allele both in Serbian population and other European populations. On the contrary, in Serbian population we did not detect any deletions or insertions in octapeptide repeat region, whereas deletions were detected in other European populations.

Schizophrenia and apolipoprotein E gene polymorphism in Serbian population.

Apolipoprotein E (APOE) gene variants are associated with alterations in brain function and increased risk of Alzheimer's disease (AD) and conflicting results have been reported in schizophrenia. Our results showed no significant differences in APOE allele or genotype frequencies between the Serbian schizophrenic patients and control individuals. However, we observed a possible association between particular subtypes of schizophrenia and APOE epsilon3/epsilon3 genotype (p = .01221) and epsilon4 allele showed a tendency toward positive association with responding to typical neuroleptics. APOE genotypes have no major influence on risk of schizophrenia, treatment and response to conventional antipsychotics, and age of onset in schizophrenia.

Mutational analysis of GJB1, MPZ, PMP22, EGR2, and LITAF/SIMPLE in Serbian Charcot-Marie-Tooth patients.

We report the results of mutational analysis in the following genes: GJB1, MPZ, PMP22, EGR2, and LITAF/SIMPLE in 57 Charcot-Marie-Tooth (CMT) patients of Serbian origin without the PMP22 duplication. We found 10 different mutations in 14 CMT patients: 6 mutations in GJB1, 3 in MPZ, and 1 in PMP22. Five of six GJB1 mutations are reported for the first time, and the most frequent one appears to be a founder mutation in the Serbian population. No mutations were found in EGR2 or LITAF. Thus, GJB1 mutation analysis should be done in patients without the PMP22 duplication and male-to-male transmission of CMT.

Coexistence of Unverricht-Lundborg disease and congenital deafness: molecular resolution of a complex comorbidity.

PURPOSE: We report on genetic analysis of a complex condition in a Serbian family of four siblings, wherein two had progressive myoclonic epilepsy (PME) and congenital deafness (CD), one had isolated congenital deafness (ICD), and one was healthy. METHODS AND RESULTS: Molecular diagnosis performed by Southern blotting confirmed Unverricht-Lundborg disease in the available sibling with PME/CD. In the sibling with ICD (heterozygote for expansion mutation in CSTB) we demonstrated recombination event between the D21S2040 marker and the CSTB gene and identified c.207delC (p.T70Xfs) mutation in the fourth exon of the transmembrane protease, serine-3 (TMPRSS3) gene (maps in close proximity to CSTB), responsible for nonsyndromic deafness in the sibling with PME/CD as well. DISCUSSION: To the best of our knowledge this is the first genetic confirmation of the coexistence of these two mutations.

Human Y-specific STR haplotypes in population of Serbia and Montenegro.

Nine Y chromosome short tandem repeat (STR) loci (DYS19, DYS385, DYS389I, DYS389II, DYS390, DYS391, DYS392 and DYS393) were analyzed in group of 237 unrelated healthy males from population of Serbia and Montenegro in order to assess the frequencies of Y haplotypes. We observed 174 different haplotypes of which 146 (61.6%) were seen only once. Y-STR allelic frequencies in Serbia and Montenegro, in general, correspond to other European populations, except for the DYS19, DYS385 and DYS389II loci.

Clinical case report atypical myopathy in a young girl with 91 CTG repeats in DM1 locus and a positive DM1 family history.

Myotonic dystrophy type 1 (DM1) is an autosomal dominant inheritable disease associated with an expansion of CTG repeats in the 3' UTR of the DMPK gene. The subject is an 11-year-old girl with atypical myopathy. Because the proband's family has a positive DM1 history, a molecular-genetic analysis for DM1 was performed. This study showed that proband had a small DMPK expansion (91 CTG repeats) although the observed myopathy would not normally be associated with DM1. These results show how the phenotypic manifestation of DM1 can have unusual symptoms with a completely unexpected relationship to genotype.