1
|
Timofeev V, Bakhteeva I, Khlopova K, Mironova R, Titareva G, Goncharova Y, Solomentsev V, Kravchenko T, Dyatlov I, Vergnaud G. New Research on the Bacillus anthracis Genetic Diversity in Siberia. Pathogens 2023; 12:1257. [PMID: 37887773 PMCID: PMC10610006 DOI: 10.3390/pathogens12101257] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2023] [Revised: 10/13/2023] [Accepted: 10/14/2023] [Indexed: 10/28/2023] Open
Abstract
Anthrax is a particularly dangerous infection of humans and ungulates caused by the Gram-positive spore-forming bacterium Bacillus anthracis. The highly monomorphic and clonal species B. anthracis is commonly divided into three main lineages, A, B, and C, which in turn are divided into several canSNP groups. We report here a phylogenetic analysis based on the whole-genome sequence (WGS) data of fifteen strains isolated predominantly in Siberia or Central and Southern Russia. We confirm the wide distribution of the cluster of strains of the B.Br.001/002 group, endemic to the Russian Arctic, which is also present in the steppe zone of Southern Siberia. We characterize additional branches within the major A.Br.001/002 polytomy comprising the A.Br.Ames and A.Br.Sterne lineages, one of which is identified in the Arctic.
Collapse
Affiliation(s)
- Vitalii Timofeev
- State Research Center for Applied Microbiology and Biotechnology (SRCAMB), 142279 Obolensk, Russia
| | - Irina Bakhteeva
- State Research Center for Applied Microbiology and Biotechnology (SRCAMB), 142279 Obolensk, Russia
| | - Kseniya Khlopova
- State Research Center for Applied Microbiology and Biotechnology (SRCAMB), 142279 Obolensk, Russia
| | - Raisa Mironova
- State Research Center for Applied Microbiology and Biotechnology (SRCAMB), 142279 Obolensk, Russia
| | - Galina Titareva
- State Research Center for Applied Microbiology and Biotechnology (SRCAMB), 142279 Obolensk, Russia
| | - Yulia Goncharova
- State Research Center for Applied Microbiology and Biotechnology (SRCAMB), 142279 Obolensk, Russia
| | - Viktor Solomentsev
- State Research Center for Applied Microbiology and Biotechnology (SRCAMB), 142279 Obolensk, Russia
| | - Tatiana Kravchenko
- State Research Center for Applied Microbiology and Biotechnology (SRCAMB), 142279 Obolensk, Russia
| | - Ivan Dyatlov
- State Research Center for Applied Microbiology and Biotechnology (SRCAMB), 142279 Obolensk, Russia
| | - Gilles Vergnaud
- Institute for Integrative Biology of the Cell (I2BC), Université Paris-Saclay, CEA, CNRS, 91198 Gif-sur-Yvette, France
| |
Collapse
|
2
|
Pavlov V, Vakhrameeva G, Mokrievich A, Platonov ME, Titareva G, Mironova R, Kombarova T, Gapelchenkova T, Dyatlov I. Development and Properties of Francisella tularensis Subsp. holarctica 15 NIIEG Vaccine Strain without the recD Gene. Vaccines (Basel) 2022; 10:108. [PMID: 35062769 PMCID: PMC8778684 DOI: 10.3390/vaccines10010108] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2021] [Revised: 12/23/2021] [Accepted: 12/28/2021] [Indexed: 02/04/2023] Open
Abstract
The genomic analysis of all subspecies F. tularensis, as found in Gen Bank NCBI, reveals the presence of genes encoding proteins like to the multifunctional RecBCD enzyme complex in E. coli and other bacteria. To date, the role of the recD gene in F. tularensis, which encodes the alpha chain of exonuclease V, in DNA metabolism processes, has not been studied either in vitro or in vivo. F. tularensis subsp. holarctica 15 NIIEG, a vaccine strain, served as the basis to create the F. tularensis 15D strain with recD deletion. The lack of the recD gene suppresses the integration of suicide plasmids with F. tularensis genome fragments into the chromosome. The modified strain showed reduced growth in vitro and in vivo. This study shows that such deletion significantly reduces the virulence of the strain in BALB/c mice.
Collapse
Affiliation(s)
- Vitaly Pavlov
- State Research Center for Applied Microbiology and Biotechnology (SRCAMB), 142279 Obolensk, Russia; (G.V.); (A.M.); (M.E.P.); (G.T.); (R.M.); (T.K.); (T.G.); (I.D.)
| | | | | | | | | | | | | | | | | |
Collapse
|
3
|
Timofeev V, Bahtejeva I, Mironova R, Titareva G, Lev I, Christiany D, Borzilov A, Bogun A, Vergnaud G. Insights from Bacillus anthracis strains isolated from permafrost in the tundra zone of Russia. PLoS One 2019; 14:e0209140. [PMID: 31116737 DOI: 10.1101/486290] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2018] [Accepted: 05/07/2019] [Indexed: 05/28/2023] Open
Abstract
This article describes Bacillus anthracis strains isolated during an outbreak of anthrax on the Yamal Peninsula in the summer of 2016 and independently in Yakutia in 2015. A common feature of these strains is their conservation in permafrost, from which they were extracted either due to the thawing of permafrost (Yamal strains) or as the result of paleontological excavations (Yakut strains). All strains isolated on the Yamal share an identical genotype belonging to lineage B.Br.001/002, pointing to a common source of infection in a territory over 250 km in length. In contrast, during the excavations in Yakutia, three genetically different strains were recovered from a single pit. One strain belongs to B.Br.001/002, and whole genome sequence analysis showed that it is most closely related to the Yamal strains in spite of the remoteness of Yamal from Yakutia. The two other strains contribute to two different branches of A.Br.008/011, one of the remarkable polytomies described so far in the B. anthracis species. The geographic distribution of the strains belonging to A.Br.008/011 is suggesting that the polytomy emerged in the thirteenth century, in combination with the constitution of a unified Mongol empire extending from China to Eastern Europe. We propose an evolutionary model for B. anthracis recent evolution in which the B lineage spread throughout Eurasia and was subsequently replaced by the A lineage except in some geographically isolated areas.
Collapse
Affiliation(s)
- Vitalii Timofeev
- State Research Center for Applied Microbiology & Biotechnology (FBIS SRCAMB), Obolensk, Russia
| | - Irina Bahtejeva
- State Research Center for Applied Microbiology & Biotechnology (FBIS SRCAMB), Obolensk, Russia
| | - Raisa Mironova
- State Research Center for Applied Microbiology & Biotechnology (FBIS SRCAMB), Obolensk, Russia
| | - Galina Titareva
- State Research Center for Applied Microbiology & Biotechnology (FBIS SRCAMB), Obolensk, Russia
| | - Igor Lev
- State Research Center for Applied Microbiology & Biotechnology (FBIS SRCAMB), Obolensk, Russia
| | - David Christiany
- Institute for Integrative Biology of the Cell (I2BC), CEA, CNRS, Univ. Paris-Sud, Université Paris-Saclay, Gif-sur-Yvette cedex, France
| | - Alexander Borzilov
- State Research Center for Applied Microbiology & Biotechnology (FBIS SRCAMB), Obolensk, Russia
| | - Alexander Bogun
- State Research Center for Applied Microbiology & Biotechnology (FBIS SRCAMB), Obolensk, Russia
| | - Gilles Vergnaud
- Institute for Integrative Biology of the Cell (I2BC), CEA, CNRS, Univ. Paris-Sud, Université Paris-Saclay, Gif-sur-Yvette cedex, France
| |
Collapse
|
4
|
Handzhiyski Y, Mironova R, Ivanov I. Effect of Acetyl Salicyilic Acid on Glycation and Mutability ofEscherichia ColiChromosomal DNA. BIOTECHNOL BIOTEC EQ 2014. [DOI: 10.1080/13102818.2009.10817616] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022] Open
|
5
|
Ofoghi H, Mironova R, Moazami N, Domonskv N, Ivanov I. Human Calcitonin Tetrameric Gene: Comparative Expression in Yeast and Transgenic Potato Plants. BIOTECHNOL BIOTEC EQ 2014. [DOI: 10.1080/13102818.1999.10819012] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022] Open
|
6
|
Dimitrova R, Mironova R, Ivanov I. Glycation of Proteins inEscherichia Coli: Effect of Nutrient Broth Ingredients on Glycation. BIOTECHNOL BIOTEC EQ 2014. [DOI: 10.1080/13102818.2004.10817094] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022] Open
|
7
|
Atanasova A, Handzhiyski Y, Sredovska-Bozhinov A, Popova E, Odjakova M, Datsenko K, Wanner BL, Ivanov I, Mironova R. Substrate Specificity of the Escherichia ColiF RLB Amadoriase. BIOTECHNOL BIOTEC EQ 2012. [DOI: 10.5504/50yrtimb.2011.0026] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022] Open
|
8
|
Mironova R, Niwa T. Molecular heterogeneity of amyloid beta2-microglobulin and modification with advanced glycation end products. J Chromatogr B Biomed Sci Appl 2001; 758:109-15. [PMID: 11482729 DOI: 10.1016/s0378-4347(01)00139-6] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
By using liquid chromatography-electrospray ionization mass spectrometry, Western blotting and N-terminal amino acid sequence analysis, we characterized the molecular heterogeneity and advanced glycation end product (AGE) modification of beta2-microglobulin (beta2m) extracted from the amyloid tissue of a hemodialysis patient. Amyloid beta2m was composed of full-length beta2m, truncated beta2m and dimer beta2m. Truncated beta2m and dimer beta2m were modified with AGEs such as imidazolone and N(e)-(carboxymethyl)lysine, and showed fluorescence characteristic of AGE. Truncated beta2m species were formed by cleavage between amino acid residues of Pro6/Ile7, Gln/Val9 and Val9/Tyr10. Heterogeneous dimer beta2m species showed the molecular masses of 22,591 and 22 675, which resulted from cross-linking between truncated beta2m.
Collapse
Affiliation(s)
- R Mironova
- Institute of Molecular Biology, Bulgarian Academy of Sciences, Sofia
| | | |
Collapse
|
9
|
Abstract
Non-enzymatic glycosylation (glycation) is a chain of chemical reactions affecting free amino groups in proteins of long-living eukaryotes. It proceeds in several steps leading to the consecutive formation of Schiff bases, Amadori products and advanced glycation end-products (AGEs). To our knowledge, this process has not been observed in prokaryotes so far. However, the present study provides clear-cut evidence that glycation takes place in bacteria despite their short life span. We have detected AGEs in recombinant human interferon gamma (rhIFN-gamma) produced in Escherichia coli as well as in total protein of the same bacterium using three different approaches: (i) Western blotting using two monoclonal antibodies raised against AGEs; (ii) fluorescent spectroscopy; and (iii) investigation of the effect of known AGE inhibitors (such as acetyl salicylic acid and thiamine) on the glycation reaction. Our study shows that non-enzymatic glycosylation is initiated during the normal growth of E. coli and results in AGE formation even after isolation of proteins. This process seems to be tightly associated with some post-translational modifications observed in the cysteineless rhIFN-gamma, such as covalent dimerization and truncation.
Collapse
Affiliation(s)
- R Mironova
- Department of Gene Regulations, Institute of Molecular Biology, Bulgarian Academy of Sciences, 1113 Sofia, Bulgaria.
| | | | | | | | | |
Collapse
|
10
|
Abstract
Pokeweed (Phytolacca americana) antiviral protein (PAP) is a highly specific ribosome-inactivating glycosidase. The PAP gene was isolated and cloned in an expression vector containing a polylinker-derived sequence (PL) but devoid of a Shine-Dalgarno (SD) sequence. Surprisingly, E. coli cells transformed with this vector produced over twice the amount of PAP than that with the consensus SD sequence. Computer analysis of the 5' terminal region of PAP mRNA revealed a nucleotide sequence (ACCUACUCGAGUUAG) which was complementary to two domains in 16S rRNA. The heptanucleotide ACCUACU (box I) is complementary to nucleotides 1434-1440 and the GAGUUAG (box II) to nucleotides 507-513 in 16S rRNA of E. coli. To examine the role of this sequence in the translation of PAP mRNA, single or both boxes were mutated and the protein yield was measured. Mutation of box I and of box II resulted in a 2.7 and 5.3 fold decrease in protein yield respectively, indicating that the PAP gene expression was dependent on the presence of both boxes. To investigate whether PL also increases expression of other genes, human calcitonin monomeric and tetrameric genes were used as reporters. It was found that the expression level was doubled compared to that by SD. These results demonstrate that the PL is an efficient translational initiator and may be used for high level expression of certain genes in E. coli. The possible mechanisms for the high level expression are discussed.
Collapse
Affiliation(s)
- J Xu
- Department of Botany, University of Toronto, Ontario, Canada.
| | | | | | | |
Collapse
|
11
|
Golshani A, Kolev V, Mironova R, AbouHaidar MG, Ivanov IG. Enhancing activity of epsilon in Escherichia coli and Agrobacterium tumefaciens cells. Biochem Biophys Res Commun 2000; 269:508-12. [PMID: 10708584 DOI: 10.1006/bbrc.2000.2327] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
Epsilon (epsilon) sequence is a bacterial enhancer of translation found in the bacteriophage T7 gene 10. It is believed that its enhancing effect of epsilon is due to a base-pairing with the nucleotides 458-467 from the helical domain 17 of Escherichia coli 16S rRNA. To prove this we have taken advantage of the difference of this domain in Agrobacterium tumefaciens and E. coli. To evaluate the significance of nucleotide complementarity for the enhancing activity of epsilon, a series of nucleotide sequences matching either E. coli or A. tumefaciens domain 17 are cloned in a binary expression vector in front of the chloramphenicol acetyltransferase (CAT) gene. The CAT assay shows that: (i) the epsilon in combination with an SD consensus sequence increases the yield of CAT in both microorganisms over that obtained with the SD alone; (ii) the epsilon sequence complementary to the A. tumefaciens domain 17 leads to a 2.71-fold increase in the yield of CAT in homologous cells but not in E. coli cells; (iii) the yield of CAT correlates with the free energy of base-pairing with the helical domain 17 in both microorganisms.
Collapse
Affiliation(s)
- A Golshani
- Department of Botany, Virology Group, University of Toronto, 25 Willcocks Street, Toronto, Ontario, M5S 3B2, Canada
| | | | | | | | | |
Collapse
|
12
|
Vassileva-Atanassova A, Niwa T, Mironova R, Ivanov I. Mass-spectral analysis of human interferon-gamma and chloramphenicol acetyltransferase I produced in two Escherichia coli strains. J Biotechnol 2000; 78:61-7. [PMID: 10702911 DOI: 10.1016/s0168-1656(99)00236-9] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Abstract
Recombinant human interferon-gamma and chloramphenicol acetyltransferase I were isolated from two Escherichia coli strains, E. coli LE329 and E. coli XL1-blue and characterized by electrospray ionization mass spectrometry (ESI-MS). The ESI-MS analysis showed higher masses in comparison with the theoretically calculated for both proteins as well as unexpected molecular heterogeneity. The ESI-MS spectral patterns of the proteins depended on the host strain used and were more heterogenous for the proteins isolated from E. coli LE392. One of the proteins (human interferon-gamma obtained from E. coli XL1-blue) was further subjected to BrCN cleavage. The ESI-MS analysis of the polypeptide mixture revealed shift in the molecular mass for two peptides including the last 26 amino acids of the human interferon-gamma molecule.
Collapse
|
13
|
Andreeva AE, Borissova BE, Mironova R, Glykos NM, Kotsifaki D, Ivanov I, Krysteva M, Kokkinidis M. Crystallization of type I chloramphenicol acetyltransferase: an approach based on the concept of ionic strength reducers. Acta Crystallogr D Biol Crystallogr 2000; 56:101-3. [PMID: 10666642 DOI: 10.1107/s090744499901481x] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Abstract
Chloramphenicol acetyltransferase (CAT) is responsible for bacterial resistance to chloramphenicol. It catalyzes inactivation of the antibiotic by acetyl-group transfer from acetyl CoA to one or both hydroxyl groups of chloramphenicol. Type I CAT possesses some unique properties which are not observed in other CAT variants. Type I CAT overexpressed in Escherichia coli was purified and crystals with a resolution limit of 2.22 A have been obtained using a novel procedure which is based on the concept of 'ionic strength reducers'. The crystals have the symmetry of space group P1 and unit-cell parameters a = 96.46, b = 113.86, c = 114.21 A, alpha = 119.9, beta = 94.1, gamma = 98.6 degrees. These dimensions are consistent with four to six trimers per unit cell, corresponding to a solvent fraction ranging from 65 to 47%.
Collapse
Affiliation(s)
- A E Andreeva
- University of Chemical Technology and Metallurgy, Department of Biotechnology, Sofia 1756, Bulgaria
| | | | | | | | | | | | | | | |
Collapse
|
14
|
Vassileva-Atanassova A, Mironova R, Nacheva G, Ivanov I. N-terminal methionine in recombinant proteins expressed in two different Escherichia coli strains. J Biotechnol 1999; 69:63-7. [PMID: 10201115 DOI: 10.1016/s0168-1656(98)00207-7] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/01/2022]
Abstract
Two genes coding for chloramphenicol acetyltransferase and human interferon gamma, respectively, were overexpressed constitutively in two different strains of Escherichia coli (E. coli LE392 and E. coli XL1). The N-terminal amino acid analysis of the purified proteins showed that: (a) the N-terminal methionine is processed more efficiently in E. coli LE392 rather than in E. coli XL1 cells; (b) the N-terminal methionine is removed better from the heterologous human interferon gamma in comparison with the homologous chloramphenicol acetyltransferase protein: and (c) there is no strong correlation between the efficiency of N-terminal procession and the yield of recombinant protein.
Collapse
|
15
|
Golshani A, Golomehova V, Mironova R, Ivanov IG, AbouHaidar MG. Does the epsilon sequence of phage T7 function as an initiator for the translation of CAT mRNA in Escherichia coli? Biochem Biophys Res Commun 1997; 236:253-6. [PMID: 9240419 DOI: 10.1006/bbrc.1997.6842] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
Epsilon (epsilon) sequence [UUAACUUUA, complementary to nucleotides 458-466 of the 16S ribosomal RNA (rRNA)] which is naturally occurring at the 5'-untranslated leader of phage T7 gene 10 mRNA was originally described as a powerful translational enhancer in Escherichia coli. Recent studies with this sequence led to controversial conclusions about its translational initiation and enhancing capability. In this study different sequence derivatives of epsilon were constructed to evaluate its efficiency not only to enhance translation of the chloramphenicol acetyltransferase (CAT) mRNA in E. coli, but also to function as an independent initiator of translation. It was observed that the epsilon sequence in combination with the CAT natural Shine-Dalgarno (SDn) or the SD consensus sequences enhanced, as expected, the translation of CAT mRNA. The natural epsilon sequence without an SD sequence failed to initiate or enhance the translation of CAT mRNA. However, when the complementarity of epsilon to 16S rRNA was increased from 9 to 16 nucleotides, epsilon alone (without the SD sequence) became an independent translational initiator with an efficiency of about 80% that obtained with the SD consensus sequence.
Collapse
Affiliation(s)
- A Golshani
- Department of Botany, Virology Group, University of Toronto, Ontario, Canada
| | | | | | | | | |
Collapse
|
16
|
Mironova R, Philipova D, Marekov L, Venkov P, Ivanov I. Secretion of oligomeric Val8-human calcitonin by Saccharomyces cerevisiae. FEMS Microbiol Lett 1991; 67:23-8. [PMID: 1778418 DOI: 10.1016/0378-1097(91)90437-f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
Monomeric human calcitonin (hCT) gene and oligomeric hCT genes composed of two, three or four head-to-tail linked monomers were fused in-frame to the yeast alpha-factor leader coding sequence wild-type and fragile mutant Saccharomyces cerevisiae strains were transformed with the constructed plasmids and the yield of recombinant protein secreted into the culture medium was measured. The yeast cells secreted equal (molar) amounts of all of the hCT variants. The recombinant proteins remained stable in the growth medium for at least 3 days. The fragile cells secreted about 30% more hCT as compared to the wild-type yeast cells.
Collapse
Affiliation(s)
- R Mironova
- Institute of Molecular Biology, Bulgarian Academy of Sciences, Sofia
| | | | | | | | | |
Collapse
|
17
|
Abstract
Recently a convenient method for detection of specific RNA sequences in bacteria has been developed but the original protocol was inapplicable to microorganisms with a rigid cell wall. Here we report a modification of the RNA colony hybridization for use with yeast. The modified method includes the following consecutive procedures: a) treatment of the yeast colonies on the membrane filter with 10% SDS at 65% C for 30 min; b) treatment of the same filter with 3 x SSC, 10% formaldehyde at 65 degrees C for 30 min; c) hybridization with 32P-labelled oligonucleotide (or DNA) specific for the RNA sequence of interest. The intensity of the radioactive signals thus obtained is comparable with that of the E. coli colonies.
Collapse
Affiliation(s)
- I Ivanov
- Institute of Molecular Biology, Bulgarian Academy of Sciences, Sofia
| | | | | | | |
Collapse
|