@article{9910,
  abstract     = {Adult height inspired the first biometrical and quantitative genetic studies and is a test-case trait for understanding heritability. The studies of height led to formulation of the classical polygenic model, that has a profound influence on the way we view and analyse complex traits. An essential part of the classical model is an assumption of additivity of effects and normality of the distribution of the residuals. However, it may be expected that the normal approximation will become insufficient in bigger studies. Here, we demonstrate that when the height of hundreds of thousands of individuals is analysed, the model complexity needs to be increased to include non-additive interactions between sex, environment and genes. Alternatively, the use of log-normal approximation allowed us to still use the additive effects model. These findings are important for future genetic and methodologic studies that make use of adult height as an exemplar trait.},
  author       = {Slavskii, Sergei A. and Kuznetsov, Ivan A. and Shashkova, Tatiana I. and Bazykin, Georgii A. and Axenovich, Tatiana I. and Kondrashov, Fyodor and Aulchenko, Yurii S.},
  issn         = {1476-5438},
  journal      = {European Journal of Human Genetics},
  number       = {7},
  pages        = {1082--1091},
  publisher    = {Springer Nature},
  title        = {{The limits of normal approximation for adult height}},
  doi          = {10.1038/s41431-021-00836-7},
  volume       = {29},
  year         = {2021},
}

@article{9905,
  abstract     = {Vaccines are thought to be the best available solution for controlling the ongoing SARS-CoV-2 pandemic. However, the emergence of vaccine-resistant strains may come too rapidly for current vaccine developments to alleviate the health, economic and social consequences of the pandemic. To quantify and characterize the risk of such a scenario, we created a SIR-derived model with initial stochastic dynamics of the vaccine-resistant strain to study the probability of its emergence and establishment. Using parameters realistically resembling SARS-CoV-2 transmission, we model a wave-like pattern of the pandemic and consider the impact of the rate of vaccination and the strength of non-pharmaceutical intervention measures on the probability of emergence of a resistant strain. As expected, we found that a fast rate of vaccination decreases the probability of emergence of a resistant strain. Counterintuitively, when a relaxation of non-pharmaceutical interventions happened at a time when most individuals of the population have already been vaccinated the probability of emergence of a resistant strain was greatly increased. Consequently, we show that a period of transmission reduction close to the end of the vaccination campaign can substantially reduce the probability of resistant strain establishment. Our results suggest that policymakers and individuals should consider maintaining non-pharmaceutical interventions and transmission-reducing behaviours throughout the entire vaccination period.},
  author       = {Rella, Simon and Kulikova, Yuliya A. and Dermitzakis, Emmanouil T. and Kondrashov, Fyodor},
  issn         = {2045-2322},
  journal      = {Scientific Reports},
  number       = {1},
  publisher    = {Springer Nature},
  title        = {{Rates of SARS-CoV-2 transmission and vaccination impact the fate of vaccine-resistant strains}},
  doi          = {10.1038/s41598-021-95025-3},
  volume       = {11},
  year         = {2021},
}

@inproceedings{15071,
  abstract     = {A mesophilic methanogenic culture, designated JL01, was isolated from Holocene permafrost in the Russian Arctic [1]. After long-term extensive cultivation at 15°C it turned out to be a tied binary culture of archaeal (JL01) and bacterial (Sphaerochaeta associata GLS2) strains.
Strain JL01 was a strict anaerobe and grew on methanol, acetate and methylamines as energy and carbon sources. Cells were irregular coccoid, non-motile, non-spore-forming, and Gram-stainpositive. Optimum conditions for growth were 24-28 oC, pH 6.8–7.3 and 0.075-0.1 M NaCl.
Phylogenetic tree reconstructions based on 16S rRNA and concatenated alignment of broadly
conserved protein-coding genes revealed its close relation to Methanosarcina mazei S-6
T (similarity 99.5%). The comparison of whole genomic sequences (ANI) of the isolate and the type strain of M.mazei was 98.5%, which is higher than the values recommended for new species. Thus strain JL01 (=VKM B-2370=JCM 31898) represents the first M. mazei isolated from permanently subzero Arcticsediments. The long-term co-cultivation of JL01 with S. associata GLS2T showed the methane production without any additional carbon and energy sources. Genome analysis of S. associata GLS2T revealed putative genes involved in methanochondroithin catabolism.},
  author       = {Oshurkova, Viktoriia and Troshina, Olga and Trubitsyn, Vladimir and Ryzhmanova, Yana and Bochkareva, Olga and Shcherbakova, Viktoria},
  booktitle    = {Proceedings of 1st International Electronic Conference on Microbiology},
  location     = {Virtual},
  publisher    = {MDPI},
  title        = {{Characterization of methanosarcina mazei JL01 isolated from holocene arctic permafrost and study of the archaeon cooperation with bacterium Sphaerochaeta associata GLS2T}},
  doi          = {10.3390/ecm2020-07116},
  year         = {2020},
}

@article{8320,
  abstract     = {The genetic code is considered to use five nucleic bases (adenine, guanine, cytosine, thymine and uracil), which form two pairs for encoding information in DNA and two pairs for encoding information in RNA. Nevertheless, in recent years several artificial base pairs have been developed in attempts to expand the genetic code. Employment of these additional base pairs increases the information capacity and variety of DNA sequences, and provides a platform for the site-specific, enzymatic incorporation of extra functional components into DNA and RNA. As a result, of the development of such expanded systems, many artificial base pairs have been synthesized and tested under various conditions. Following many stages of enhancement, unnatural base pairs have been modified to eliminate their weak points, qualifying them for specific research needs. Moreover, the first attempts to create a semi-synthetic organism containing DNA with unnatural base pairs seem to have been successful. This further extends the possible applications of these kinds of pairs. Herein, we describe the most significant qualities of unnatural base pairs and their actual applications.},
  author       = {Mukba, S. A. and Vlasov, Petr and Kolosov, P. M. and Shuvalova, E. Y. and Egorova, T. V. and Alkalaeva, E. Z.},
  issn         = {1608-3245},
  journal      = {Molecular Biology},
  number       = {4},
  pages        = {475--484},
  publisher    = {Springer Nature},
  title        = {{Expanding the genetic code: Unnatural base pairs in biological systems}},
  doi          = {10.1134/S0026893320040111},
  volume       = {54},
  year         = {2020},
}

@article{8321,
  abstract     = {The genetic code is considered to use five nucleic bases (adenine, guanine, cytosine, thymine and uracil), which form two pairs for encoding information in DNA and two pairs for encoding information in RNA. Nevertheless, in recent years several artificial base pairs have been developed in attempts to expand the genetic code. Employment of these additional base pairs increases the information capacity and variety of DNA sequences, and provides a platform for the site-specific, enzymatic incorporation of extra functional components into DNA and RNA. As a result, of the development of such expanded systems, many artificial base pairs have been synthesized and tested under various conditions. Following many stages of enhancement, unnatural base pairs have been modified to eliminate their weak points, qualifying them for specific research needs. Moreover, the first attempts to create a semi-synthetic organism containing DNA with unnatural base pairs seem to have been successful. This further extends the possible applications of these kinds of pairs. Herein, we describe the most significant qualities of unnatural base pairs and their actual applications.},
  author       = {Mukba, S. A. and Vlasov, Petr and Kolosov, P. M. and Shuvalova, E. Y. and Egorova, T. V. and Alkalaeva, E. Z.},
  issn         = {0026-8984},
  journal      = {Molekuliarnaia biologiia},
  number       = {4},
  pages        = {531--541},
  publisher    = {Russian Academy of Sciences},
  title        = {{Expanding the genetic code: Unnatural base pairs in biological systems}},
  doi          = {10.31857/S0026898420040126},
  volume       = {54},
  year         = {2020},
}

@article{8645,
  abstract     = {Epistasis, the context-dependence of the contribution of an amino acid substitution to fitness, is common in evolution. To detect epistasis, fitness must be measured for at least four genotypes: the reference genotype, two different single mutants and a double mutant with both of the single mutations. For higher-order epistasis of the order n, fitness has to be measured for all 2n genotypes of an n-dimensional hypercube in genotype space forming a ‘combinatorially complete dataset’. So far, only a handful of such datasets have been produced by manual curation. Concurrently, random mutagenesis experiments have produced measurements of fitness and other phenotypes in a high-throughput manner, potentially containing a number of combinatorially complete datasets. We present an effective recursive algorithm for finding all hypercube structures in random mutagenesis experimental data. To test the algorithm, we applied it to the data from a recent HIS3 protein dataset and found all 199 847 053 unique combinatorially complete genotype combinations of dimensionality ranging from 2 to 12. The algorithm may be useful for researchers looking for higher-order epistasis in their high-throughput experimental data.},
  author       = {Esteban, Laura A and Lonishin, Lyubov R and Bobrovskiy, Daniil M and Leleytner, Gregory and Bogatyreva, Natalya S and Kondrashov, Fyodor and Ivankov, Dmitry N },
  issn         = {1460-2059},
  journal      = {Bioinformatics},
  number       = {6},
  pages        = {1960--1962},
  publisher    = {Oxford University Press},
  title        = {{HypercubeME: Two hundred million combinatorially complete datasets from a single experiment}},
  doi          = {10.1093/bioinformatics/btz841},
  volume       = {36},
  year         = {2020},
}

@article{8700,
  abstract     = {Translation termination is a finishing step of protein biosynthesis. The significant role in this process belongs not only to protein factors of translation termination but also to the nearest nucleotide environment of stop codons. There are numerous descriptions of stop codons readthrough, which is due to specific nucleotide sequences behind them. However, represented data are segmental and don’t explain the mechanism of the nucleotide context influence on translation termination. It is well known that stop codon UAA usage is preferential for A/T-rich genes, and UAG, UGA—for G/C-rich genes, which is related to an expression level of these genes. We investigated the connection between a frequency of nucleotides occurrence in 3' area of stop codons in the human genome and their influence on translation termination efficiency. We found that 3' context motif, which is cognate to the sequence of a stop codon, stimulates translation termination. At the same time, the nucleotide composition of 3' sequence that differs from stop codon, decreases translation termination efficiency.},
  author       = {Sokolova, E. E. and Vlasov, Petr and Egorova, T. V. and Shuvalov, A. V. and Alkalaeva, E. Z.},
  issn         = {1608-3245},
  journal      = {Molecular Biology},
  number       = {5},
  pages        = {739--748},
  publisher    = {Springer Nature},
  title        = {{The influence of A/G composition of 3' stop codon contexts on translation termination efficiency in eukaryotes}},
  doi          = {10.1134/S0026893320050088},
  volume       = {54},
  year         = {2020},
}

@article{8701,
  abstract     = {Translation termination is a finishing step of protein biosynthesis. The significant role in this process belongs not only to protein factors of translation termination but also to the nearest nucleotide environment of stop codons. There are numerous descriptions of stop codons readthrough, which is due to specific nucleotide sequences behind them. However, represented data are segmental and don’t explain the mechanism of the nucleotide context influence on translation termination. It is well known that stop codon UAA usage is preferential for A/T-rich genes, and UAG, UGA—for G/C-rich genes, which is related to an expression level of these genes. We investigated the connection between a frequency of nucleotides occurrence in 3' area of stop codons in the human genome and their influence on translation termination efficiency. We found that 3' context motif, which is cognate to the sequence of a stop codon, stimulates translation termination. At the same time, the nucleotide composition of 3' sequence that differs from stop codon, decreases translation termination efficiency.},
  author       = {Sokolova, E. E. and Vlasov, Petr and Egorova, T. V. and Shuvalov, A. V. and Alkalaeva, E. Z.},
  issn         = {0026-8984},
  journal      = {Molekuliarnaia biologiia},
  number       = {5},
  pages        = {837--848},
  publisher    = {Russian Academy of Sciences},
  title        = {{The influence of A/G composition of 3' stop codon contexts on translation termination efficiency in eukaryotes}},
  doi          = {10.31857/S0026898420050080},
  volume       = {54},
  year         = {2020},
}

@article{7622,
  abstract     = {The International Young Physicists' Tournament (IYPT) continued in 2018 in Beijing, China and 2019 in Warsaw, Poland with its 31st and 32nd editions. The IYPT is a modern scientific competition for teams of high school students, also known as the Physics World Cup. It involves long-term theoretical and experimental work focused on solving 17 publicly announced open-ended problems in teams of five. On top of that, teams have to present their solutions in front of other teams and a scientific jury, and get opposed and reviewed by their peers. Here we present a brief information about the competition with a specific focus on one of the IYPT 2018 tasks, the 'Ring Oiler'. This seemingly simple mechanical problem appeared to be of such a complexity that even the dozens of participating teams and jurying scientists were not able to solve all of its subtleties.},
  author       = {Plesch, Martin and Plesník, Samuel and Ruzickova, Natalia},
  issn         = {1361-6404},
  journal      = {European Journal of Physics},
  number       = {3},
  publisher    = {IOP Publishing},
  title        = {{The IYPT and the 'Ring Oiler' problem}},
  doi          = {10.1088/1361-6404/ab6414},
  volume       = {41},
  year         = {2020},
}

@article{7931,
  abstract     = {In the course of sample preparation for Next Generation Sequencing (NGS), DNA is fragmented by various methods. Fragmentation shows a persistent bias with regard to the cleavage rates of various dinucleotides. With the exception of CpG dinucleotides the previously described biases were consistent with results of the DNA cleavage in solution. Here we computed cleavage rates of all dinucleotides including the methylated CpG and unmethylated CpG dinucleotides using data of the Whole Genome Sequencing datasets of the 1000 Genomes project. We found that the cleavage rate of CpG is significantly higher for the methylated CpG dinucleotides. Using this information, we developed a classifier for distinguishing cancer and healthy tissues based on their CpG islands statuses of the fragmentation. A simple Support Vector Machine classifier based on this algorithm shows an accuracy of 84%. The proposed method allows the detection of epigenetic markers purely based on mechanochemical DNA fragmentation, which can be detected by a simple analysis of the NGS sequencing data.},
  author       = {Uroshlev, Leonid A. and Abdullaev, Eldar T. and Umarova, Iren R. and Il’Icheva, Irina A. and Panchenko, Larisa A. and Polozov, Robert V. and Kondrashov, Fyodor and Nechipurenko, Yury D. and Grokhovsky, Sergei L.},
  issn         = {2045-2322},
  journal      = {Scientific Reports},
  publisher    = {Springer Nature},
  title        = {{A method for identification of the methylation level of CpG islands from NGS data}},
  doi          = {10.1038/s41598-020-65406-1},
  volume       = {10},
  year         = {2020},
}

@article{7603,
  abstract     = {Plants are exposed to a variety of abiotic and biotic stresses that may result in DNA damage. Endogenous processes - such as DNA replication, DNA recombination, respiration, or photosynthesis - are also a threat to DNA integrity. It is therefore essential to understand the strategies plants have developed for DNA damage detection, signaling, and repair. Alternative splicing (AS) is a key post-transcriptional process with a role in regulation of gene expression. Recent studies demonstrate that the majority of intron-containing genes in plants are alternatively spliced, highlighting the importance of AS in plant development and stress response. Not only does AS ensure a versatile proteome and influence the abundance and availability of proteins greatly, it has also emerged as an important player in the DNA damage response (DDR) in animals. Despite extensive studies of DDR carried out in plants, its regulation at the level of AS has not been comprehensively addressed. Here, we provide some insights into the interplay between AS and DDR in plants.},
  author       = {Nimeth, Barbara Anna and Riegler, Stefan and Kalyna, Maria},
  issn         = {1664-462X},
  journal      = {Frontiers in Plant Science},
  publisher    = {Frontiers},
  title        = {{Alternative splicing and DNA damage response in plants}},
  doi          = {10.3389/fpls.2020.00091},
  volume       = {11},
  year         = {2020},
}

@article{7889,
  abstract     = {Autoluminescent plants engineered to express a bacterial bioluminescence gene cluster in plastids have not been widely adopted because of low light output. We engineered tobacco plants with a fungal bioluminescence system that converts caffeic acid (present in all plants) into luciferin and report self-sustained luminescence that is visible to the naked eye. Our findings could underpin development of a suite of imaging tools for plants.},
  author       = {Mitiouchkina, Tatiana and Mishin, Alexander S. and Gonzalez Somermeyer, Louisa and Markina, Nadezhda M. and Chepurnyh, Tatiana V. and Guglya, Elena B. and Karataeva, Tatiana A. and Palkina, Kseniia A. and Shakhova, Ekaterina S. and Fakhranurova, Liliia I. and Chekova, Sofia V. and Tsarkova, Aleksandra S. and Golubev, Yaroslav V. and Negrebetsky, Vadim V. and Dolgushin, Sergey A. and Shalaev, Pavel V. and Shlykov, Dmitry and Melnik, Olesya A. and Shipunova, Victoria O. and Deyev, Sergey M. and Bubyrev, Andrey I. and Pushin, Alexander S. and Choob, Vladimir V. and Dolgov, Sergey V. and Kondrashov, Fyodor and Yampolsky, Ilia V. and Sarkisyan, Karen S.},
  issn         = {1546-1696},
  journal      = {Nature Biotechnology},
  pages        = {944--946},
  publisher    = {Springer Nature},
  title        = {{Plants with genetically encoded autoluminescence}},
  doi          = {10.1038/s41587-020-0500-9},
  volume       = {38},
  year         = {2020},
}

@article{8707,
  abstract     = {Dynamic changes in the three-dimensional (3D) organization of chromatin are associated with central biological processes, such as transcription, replication and development. Therefore, the comprehensive identification and quantification of these changes is fundamental to understanding of evolutionary and regulatory mechanisms. Here, we present Comparison of Hi-C Experiments using Structural Similarity (CHESS), an algorithm for the comparison of chromatin contact maps and automatic differential feature extraction. We demonstrate the robustness of CHESS to experimental variability and showcase its biological applications on (1) interspecies comparisons of syntenic regions in human and mouse models; (2) intraspecies identification of conformational changes in Zelda-depleted Drosophila embryos; (3) patient-specific aberrant chromatin conformation in a diffuse large B-cell lymphoma sample; and (4) the systematic identification of chromatin contact differences in high-resolution Capture-C data. In summary, CHESS is a computationally efficient method for the comparison and classification of changes in chromatin contact data.},
  author       = { Galan, Silvia and Machnik, Nick N and Kruse, Kai and Díaz, Noelia and Marti-Renom, Marc A and Vaquerizas, Juan M},
  issn         = {1546-1718},
  journal      = {Nature Genetics},
  pages        = {1247--1255},
  publisher    = {Springer Nature},
  title        = {{CHESS enables quantitative comparison of chromatin contact data and automatic feature extraction}},
  doi          = {10.1038/s41588-020-00712-y},
  volume       = {52},
  year         = {2020},
}

@article{6506,
  abstract     = {How does environmental complexity affect the evolution of single genes? Here, we measured the effects of a set of Bacillus subtilis glutamate dehydrogenase mutants across 19 different environments—from phenotypically homogeneous single-cell populations in liquid media to heterogeneous biofilms, plant roots and soil populations. The effects of individual gene mutations on organismal fitness were highly reproducible in liquid cultures. However, 84% of the tested alleles showed opposing fitness effects under different growth conditions (sign environmental pleiotropy). In colony biofilms and soil samples, different alleles dominated in parallel replica experiments. Accordingly, we found that in these heterogeneous cell populations the fate of mutations was dictated by a combination of selection and drift. The latter relates to programmed prophage excisions that occurred during biofilm development. Overall, for each condition, a wide range of glutamate dehydrogenase mutations persisted and sometimes fixated as a result of the combined action of selection, pleiotropy and chance. However, over longer periods and in multiple environments, nearly all of this diversity would be lost—across all the environments and conditions that we tested, the wild type was the fittest allele.},
  author       = {Noda-García, Lianet and Davidi, Dan and Korenblum, Elisa and Elazar, Assaf and Putintseva, Ekaterina and Aharoni, Asaph and Tawfik, Dan S.},
  issn         = {2058-5276},
  journal      = {Nature Microbiology},
  number       = {7},
  pages        = {1221–1230},
  publisher    = {Springer Nature},
  title        = {{Chance and pleiotropy dominate genetic diversity in complex bacterial environments}},
  doi          = {10.1038/s41564-019-0412-y},
  volume       = {4},
  year         = {2019},
}

@article{7181,
  abstract     = {Multiple sequence alignments (MSAs) are used for structural1,2 and evolutionary predictions1,2, but the complexity of aligning large datasets requires the use of approximate solutions3, including the progressive algorithm4. Progressive MSA methods start by aligning the most similar sequences and subsequently incorporate the remaining sequences, from leaf-to-root, based on a guide-tree. Their accuracy declines substantially as the number of sequences is scaled up5. We introduce a regressive algorithm that enables MSA of up to 1.4 million sequences on a standard workstation and substantially improves accuracy on datasets larger than 10,000 sequences. Our regressive algorithm works the other way around to the progressive algorithm and begins by aligning the most dissimilar sequences. It uses an efficient divide-and-conquer strategy to run third-party alignment methods in linear time, regardless of their original complexity. Our approach will enable analyses of extremely large genomic datasets such as the recently announced Earth BioGenome Project, which comprises 1.5 million eukaryotic genomes6.},
  author       = {Garriga, Edgar and Di Tommaso, Paolo and Magis, Cedrik and Erb, Ionas and Mansouri, Leila and Baltzis, Athanasios and Laayouni, Hafid and Kondrashov, Fyodor and Floden, Evan and Notredame, Cedric},
  issn         = {1546-1696},
  journal      = {Nature Biotechnology},
  number       = {12},
  pages        = {1466--1470},
  publisher    = {Springer Nature},
  title        = {{Large multiple sequence alignments with a root-to-leaf regressive method}},
  doi          = {10.1038/s41587-019-0333-6},
  volume       = {37},
  year         = {2019},
}

@misc{9731,
  abstract     = {OGs with putative pseudogenes by the number of affected genomes in different chlamydial species. Frameshift and nonsense mutations located less than 60 bp upstreamof the gene end or present in a single genome from the corresponding OG were excluded. (CSV 31 kb)},
  author       = {Sigalova, Olga and Chaplin, Andrei and Bochkareva, Olga and Shelyakin, Pavel and Filaretov, Vsevolod and Akkuratov, Evgeny and Burskaia, Valentina and Gelfand, Mikhail S.},
  publisher    = {Springer Nature},
  title        = {{Additional file 11 of Chlamydia pan-genomic analysis reveals balance between host adaptation and selective pressure to genome reduction}},
  doi          = {10.6084/m9.figshare.9808772.v1},
  year         = {2019},
}

@misc{9783,
  abstract     = {Predicted frameshift and nonsense mutations in Chlamydial pan-genome. For the analysis of putative pseudogenes, events located less than 60 bp. away from gene end or present in a single genome from the corresponding OG were excluded. (CSV 600 kb)},
  author       = {Sigalova, Olga M. and Chaplin, Andrei V. and Bochkareva, Olga and Shelyakin, Pavel V. and Filaretov, Vsevolod A. and Akkuratov, Evgeny E. and Burskaia, Valentina and Gelfand, Mikhail S.},
  publisher    = {Springer Nature},
  title        = {{Additional file 10 of Chlamydia pan-genomic analysis reveals balance between host adaptation and selective pressure to genome reduction}},
  doi          = {10.6084/m9.figshare.9808760.v1},
  year         = {2019},
}

@misc{9890,
  abstract     = {Distribution of OGs with mosaic phyletic patterns across species (complete genomes only). (CSV 7 kb)},
  author       = {Sigalova, Olga M. and Chaplin, Andrei V. and Bochkareva, Olga and Shelyakin, Pavel V. and Filaretov, Vsevolod A. and Akkuratov, Evgeny E. and Burskaia, Valentina and Gelfand, Mikhail S.},
  publisher    = {Springer Nature},
  title        = {{Additional file 15 of Chlamydia pan-genomic analysis reveals balance between host adaptation and selective pressure to genome reduction}},
  doi          = {10.6084/m9.figshare.9808802.v1},
  year         = {2019},
}

@misc{9892,
  abstract     = {Distribution of OGs with mosaic phyletic patterns across species (all genomes). (CSV 10 kb)},
  author       = {Sigalova, Olga M. and Chaplin, Andrei V and Bochkareva, Olga and Shelyakin, Pavel V. and Filaretov, Vsevolod A. and Akkuratov, Evgeny E. and Burskaia, Valentina and Gelfand, Mikhail S.},
  publisher    = {Springer Nature},
  title        = {{Additional file 16 of Chlamydia pan-genomic analysis reveals balance between host adaptation and selective pressure to genome reduction}},
  doi          = {10.6084/m9.figshare.9808814.v1},
  year         = {2019},
}

@misc{9893,
  abstract     = {Summary of peripheral genesa phyletic patterns and tree concordance. (CSV 26 kb)},
  author       = {Sigalova, Olga M. and Chaplin, Andrei V. and Bochkareva, Olga and Shelyakin, Pavel V. and Filaretov, Vsevolod A. and Akkuratov, Evgeny E. and Burskaia, Valentina and Gelfand, Mikhail S.},
  publisher    = {Springer Nature},
  title        = {{Additional file 17 of Chlamydia pan-genomic analysis reveals balance between host adaptation and selective pressure to genome reduction}},
  doi          = {10.6084/m9.figshare.9808820.v1},
  year         = {2019},
}