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.

Crosstalk between Glycogen-Selective Autophagy, Autophagy and Apoptosis as a Road towards Modifier Gene Discovery and New Therapeutic Strategies for Glycogen Storage Diseases.

Glycogen storage diseases (GSDs) are rare metabolic monogenic disorders characterized by an excessive accumulation of glycogen in the cell. However, monogenic disorders are not simple regarding genotype-phenotype correlation. Genes outside the major disease-causing locus could have modulatory effect on GSDs, and thus explain the genotype-phenotype inconsistencies observed in these patients. Nowadays, when the sequencing of all clinically relevant genes, whole human exomes, and even whole human genomes is fast, easily available and affordable, we have a scientific obligation to holistically analyze data and draw smarter connections between genotype and phenotype. Recently, the importance of glycogen-selective autophagy for the pathophysiology of disorders of glycogen metabolism have been described. Therefore, in this manuscript, we review the potential role of genes involved in glycogen-selective autophagy as modifiers of GSDs. Given the small number of genes associated with glycogen-selective autophagy, we also include genes, transcription factors, and non-coding RNAs involved in autophagy. A cross-link with apoptosis is addressed. All these genes could be analyzed in GSD patients with unusual discrepancies between genotype and phenotype in order to discover genetic variants potentially modifying their phenotype. The discovery of modifier genes related to glycogen-selective autophagy and autophagy will start a new chapter in understanding of GSDs and enable the usage of autophagy-inducing drugs for the treatment of this group of rare-disease patients.

Identification and Classification of Novel Genetic Variants: En Route to the Diagnosis of Primary Ciliary Dyskinesia.

Primary ciliary dyskinesia (PCD) is a disease caused by impaired function of motile cilia. PCD mainly affects the lungs and reproductive organs. Inheritance is autosomal recessive and X-linked. PCD patients have diverse clinical manifestations, thus making the establishment of proper diagnosis challenging. The utility of next-generation sequencing (NGS) technology for diagnostic purposes allows for better understanding of the PCD genetic background. However, identification of specific disease-causing variants is difficult. The main aim of this study was to create a unique guideline that will enable the standardization of the assessment of novel genetic variants within PCD-associated genes. The designed pipeline consists of three main steps: (1) sequencing, detection, and identification of genes/variants; (2) classification of variants according to their effect; and (3) variant characterization using in silico structural and functional analysis. The pipeline was validated through the analysis of the variants detected in a well-known PCD disease-causing gene (DNAI1) and the novel candidate gene (SPAG16). The application of this pipeline resulted in identification of potential disease-causing variants, as well as validation of the variants pathogenicity, through their analysis on transcriptional, translational, and posttranslational levels. The application of this pipeline leads to the confirmation of PCD diagnosis and enables a shift from candidate to PCD disease-causing gene.

Association of Vitamin D, Zinc and Selenium Related Genetic Variants With COVID-19 Disease Severity.

Background: COVID-19 pandemic has proved to be an unrelenting health threat for more than a year now. The emerging amount of data indicates that vitamin D, zinc and selenium could be important for clinical presentation of COVID-19. Here, we investigated association of genetic variants related to the altered level and bioavailability of vitamin D, zinc and selenium with clinical severity of COVID-19. Methods: We analyzed variants in genes significant for the status of vitamin D (DHCR7/NADSYN1 rs12785878, GC rs2282679, CYP2R1 rs10741657, and VDR rs2228570), zinc (PPCDC rs2120019) and selenium (DMGDH rs17823744) in 120 Serbian adult and pediatric COVID-19 patients using allelic discrimination. Furthermore, we carried out comparative population genetic analysis among European and other worldwide populations to investigate variation in allelic frequencies of selected variants. Results: Study showed that DHCR7/NADSYN rs12785878 and CYP2R1 rs10741657 variants were associated with severe COVID-19 in adults (p = 0.03, p = 0.017, respectively); carriers of DHCR7/NADSYN TG+GG and CYP2R1 GG genotypes had 0.21 and 5.9 the odds for developing severe disease, OR 0.21 (0.05-0.9) and OR 5.9 (1.4-25.2), respectively. There were no associations between selected genetic variants and disease severity in pediatric patients. Comparative population genetic analysis revealed that Serbian population had the lowest frequency of CYP2R1 rs10741657 G allele compared to other non-Finish Europeans (0.58 compared to 0.69 and 0.66 in Spanish and Italian population, respectively), suggesting that other populations should also investigate the relationship of CYP2R1 variant and the COVID-19 disease course. Conclusion: The results of the study indicated that vitamin D related genetic variants were implicated in severe COVID-19 in adults. This could direct prevention strategies based on population specific nutrigenetic profiles.

A novel 9 bp deletion (c.1271_1279delGTGCCCGCG) in exon 10 of CYP21A2 gene causing severe congenital adrenal hyperplasia.

BACKGROUND: Congenital adrenal hyperplasia (CAH) is an autosomal recessive disorder of adrenal steroidogenesis with a broad spectrum of clinical presentations, ranging from the severe classical salt-wasting (SW) and simple-virilizing (SV) form, to the mild nonclassical form. A large variety of CYP21A2 genotypes in correlation with phenotype have been described. MATERIALS AND METHODS: DNA samples from a 14-day-old male newborn with clinical and laboratory signs of SW CAH and family members were subjected for molecular analysis of the nine most common point CYP21A2 mutations by ACRS/PCR method. Direct DNA sequencing of the whole CYP21A2 gene was performed to detect the second mutant allele in the patient. The in silico predicting analysis and the crystal structure analysis of the mutated CYP21A2 protein have been performed. RESULTS: Molecular analysis confirmed that the patient was compound heterozygote carrying p.Q318X mutation inherited from the mother and a novel c.1271_1279delGTGCCCGCG (p.G424_R426del) variant in exon 10 inherited from the father. The in silico predicting software tools classified the novel mutation as pathogenic. Crystal structure analysis showed that the three residues affected by the novel in-frame deletion form several hydrogen bonds that could lead to impaired stability and function of the CYP21A2 protein. These findings were concordant with the patient's phenotype. The need of several molecular methods to elucidate the genotype in this patient has also been discussed. CONCLUSIONS: A novel 9 bp deletion in CYP21A2 gene with predicted pathogenic effect on the enzyme activity was detected in neonatal patient causing severe SW CAH.

Genetic variants in TNFA, LTA, TLR2 and TLR4 genes and risk of sepsis in patients with severe trauma: nested case-control study in a level-1 trauma centre in SERBIA.

INTRODUCTION: Single nucleotide variants (SNVs) represent important genetic risk factors for susceptibility to posttraumatic sepsis and a potential target for immunotherapy. We aimed to evaluate the association between 8 different SNVs within tumor necrosis factor alpha (TNFA), lymphotoxin alpha (LTA) and Toll-like receptor (TLR2 and TLR4) genes and the risk of posttraumatic sepsis. METHODS: Nested case-control study was conducted in the emergency department of the Clinical Centre of Serbia including 228 traumatized patients (44 with sepsis and 184 without sepsis). To compare the results of trauma subjects with the data from the general population, a control group of 101 healthy persons was included in the study. Genotyping of TNFA (rs1800629 and rs361525), LTA (rs909253), TLR2 (rs3804099, rs4696480 and rs3804100), and TLR4 (rs4986790 and rs4986791) was performed for all patients within all three groups using the real-time PCR method. MutationTaster database and in silico software SIFT were used to predict the variant pathogenic effect. RESULTS: Carriage of the G allele of the TNFA rs1800629 gene variant (OR 2.1, 95%CI 1.06-4.16) and T allele-carriage of the TLR4 rs4986791 genetic variant (OR 3.02, 95%CI 1.31-6.57) were associated with significantly higher risk of sepsis in trauma patients when compared to the general population prone to sepsis and traumatized patients without developing a sepsis, respectively. Of these two variants, only variant in TLR4 gene (rs4986791) has been labeled as disease causing by both the MutationTaster database and the in-silico software SIFT, which further supports the role of this variant in various pathologies including sepsis. For the remaining six variants no significant association with the susceptibility to sepsis was detected. CONCLUSIONS: Carriage of the G allele of the TNFA rs1800629 gene variant and T allele-carriage of the TLR4 rs4986791 genetic variant confer significant risk of posttraumatic sepsis. TLR4 gene variants (rs4986790 and rs4986791) has been labelled as disease causing.

Impact of genotype on neutropenia in a large cohort of Serbian patients with glycogen storage disease type Ib.

BACKGROUND: Glycogen storage disease type Ib (GSD-Ib) is an inherited metabolic disorder caused by autosomal recessive mutations in SLC37A4 coding for the glucose-6-phosphate transporter. Neutropenia represents major feature of GSD-Ib along with metabolic disturbances. Previous research in GSD-Ib patients did not reveal significant genotype-phenotype correlation. Our objective was to explore the frequency and severity of neutropenia and it's complications in relation to genotype of GSD-Ib patients. METHODS: We estimated cumulative incidence of neutropenia and severe neutropenia, relation of genotype to absolute neutrophil count (ANC), and dynamics of ANC during serious bacterial infections (SBI) in a cohort of Serbian GSD Ib patients. Impact of genotype on GSD Ib-related inflammatory bowel disease (IBD) was also assessed. RESULTS: Absolute neutrophil count (ANC) < 1500/mm3 was present in all 33 patients, with severe neutropenia (ANC<500/mm3) occurring in 60.6% of patients. The median age at neutropenia onset was 24 months, while severe neutropenia developed at median of 4.5 years. The ANC was elevated during 90.5% episodes of SBI. Genotypes c.81T>A/c.785G>A and c.81T>A/c.1042_1043delCT are associated with earlier onset of neutropenia. Patients carrying c.785G>A mutation express a higher capacity for ANC increase during SBI. Inflammatory bowel disease was diagnosed in 8 patients (24.2% of total) with median age of onset at 7 years. Risk for IBD occurrence was not significantly affected by gender, genotype and severity of neutropenia. CONCLUSIONS: We may conclude that certain mutations in SLC37A4 influence the risk for severe neutropenia occurrence but also affect the capacity to increase ANC during SBI.

CRISPR/Cas9 genome editing of SLC37A4 gene elucidates the role of molecular markers of endoplasmic reticulum stress and apoptosis in renal involvement in glycogen storage disease type Ib.

Glycogen storage disease type Ib (GSD Ib) is an autosomal recessive disorder, caused by a deficiency of ubiquitously expressed SLC37A4 protein. Deficiency of SLC37A4 leads to abnormal storage of glycogen in the liver and kidneys, resulting in long-term complications of renal disease and hepatocellular adenomas, whose mechanisms are poorly understood. Molecular markers of the adaptive responses to the metabolic stress caused by a deficiency of SLC37A4, such as markers related to the endoplasmic reticulum (ER) stress and unfolded protein response (UPR), have not been extensively studied. The aim of this study was to investigate the expression of molecular markers of the UPR response and apoptosis related to a deficiency of SLC37A4 in kidney cells. For that purpose, we intended to establish a human kidney cell model system for GSD Ib. The novel variant c.248G>A, found in GSD Ib patients, was introduced into the Flp-In™T-REx™-293 cell line using CRISPR/Cas9-mediated precise gene editing method, resulting in significant decrease of SLC37A4 gene expression. In this model system we used RT-qPCR analysis to investigate the expression of molecular markers of the UPR response (ATF4, DDIT3, HSPA5, and XBP1s) and apoptosis (BCL2, BAX). We demonstrated that under chronic metabolic stress conditions caused by SLC37A4 deficiency, the ER stress-induced UPR was triggered, resulting in suppression of the UPR molecular markers and cell survival promotion (decreased expression levels of ATF4, DDIT3, HSPA5, with the exception of XBP1s). However, persistent metabolic stress overrides an adaptation and induces apoptosis through increased expression of pro-apoptotic markers (decreased ratio of BCL2/BAX genes). We established a cellular model system characterized by a deficiency of SLC37A4, which presents pathological manifestations of GSD Ib in the kidney. Expression analysis in a novel model system supports the hypothesis that renal dysfunction in the GSD Ib is partly due to the ER stress and increased apoptosis.

Genomic profiling supports the diagnosis of primary ciliary dyskinesia and reveals novel candidate genes and genetic variants.

Primary ciliary dyskinesia (PCD) is a rare inherited autosomal recessive or X-linked disorder that mainly affects lungs. Dysfunction of respiratory cilia causes symptoms such as chronic rhinosinusitis, coughing, rhinitis, conductive hearing loss and recurrent lung infections with bronchiectasis. It is now well known that pathogenic genetic changes lead to ciliary dysfunction. Here we report usage of clinical-exome based NGS approach in order to reveal underlying genetic causes in cohort of 21 patient with diagnosis of PCD. By detecting 18 (12 novel) potentially pathogenic genetic variants, we established the genetic cause of 11 (9 unrelated) patients. Genetic variants were detected in six PCD disease-causing genes, as well as in SPAG16 and SPAG17 genes, that were not detected in PCD patients so far, but were related to some symptoms of PCD. The most frequently mutated gene in our cohort was DNAH5 (27.77%). Identified variants were in homozygous, compound heterozygous and trans-heterozygous state. For detailed characterization of one novel homozygous genetic variant in DNAI1 gene (c. 947_948insG, p. Thr318TyrfsTer11), RT-qPCR and Western Blot analysis were performed. Molecular diagnostic approach applied in this study enables analysis of 29 PCD disease-causing and related genes. It resulted in mutation detection rate of 50% and enabled discovery of twelve novel mutations and pointed two possible novel PCD candidate genes.

Functional Characterization of Novel Phenylalanine Hydroxylase p.Gln226Lys Mutation Revealed Its Non-responsiveness to Tetrahydrobiopterin Treatment in Hepatoma Cellular Model.

Treatment with tetrahydrobiopterin (BH4) is the latest therapeutic option approved for patients with phenylketonuria (PKU)-one of the most frequent inborn metabolic diseases. PKU or phenylalanine hydroxylase (PAH) deficiency is caused by mutations in the PAH gene. Given that some PAH mutations are responsive to BH4 treatment while others are non-responsive, for every novel mutation that is discovered it is essential to confirm its pathogenic effect and to assess its responsiveness to a BH4 treatment in vitro, before the drug is administered to patients. We found a c.676C>A (p.Gln226Lys) mutation in the PAH gene in two unrelated patients with PKU. The corresponding aberrant protein has never been functionally characterized in vitro and its response to BH4 treatment is unknown. Computational analyses proposed that glutamine at position 226 is an important, evolutionary conserved amino acid while the substitution with lysine probably disturbs tertiary protein structure and impacts posttranslational PAH modifications. Using hepatoma cellular model, we demonstrated that the amount of mutant p.Gln226Lys PAH detected by Western blot was only 1.2% in comparison to wild-type PAH. The addition of sepiapterin, intracellular precursor of BH4, did not increase PAH protein yield thus marking p.Gln226Lys as BH4-non-responsive mutation. Therefore, computational, experimental, and clinical data were all in accordance showing that p.Gln226Lys is a severe pathogenic PAH mutation. Its non-responsiveness to BH4 treatment in hepatoma cellular model should be considered when deciding treatment options for PKU patients carrying this mutation. Consequently, our study will facilitate clinical genetic practice, particularly genotype-based stratification of PKU treatment.