Tuberculosis (TB) poses a worldwide threat due to advancing multidrug-resistant strains and deadly co-infections with Human immunodeficiency virus. Today large amounts of Mycobacterium tuberculosis whole genome sequencing data are being assessed broadly and yet there exists no comprehensive online resource that connects M. tuberculosis genome variants with geographic origin, with drug resistance or with clinical outcome. Here we describe a broadly inclusive unifying Genome-wide Mycobacterium tuberculosis Variation (GMTV) database, ( http://mtb.dobzhanskycenter.org ) that catalogues genome variations of M. tuberculosis strains collected across Russia. GMTV contains a broad spectrum of data derived from different sources and related to M. tuberculosis molecular biology, epidemiology, TB clinical outcome, year and place of isolation, drug resistance profiles and displays the variants across the genome using a dedicated genome browser. GMTV database, which includes 1084 genomes and over 69,000 SNP or Indel variants, can be queried about M. tuberculosis genome variation and putative associations with drug resistance, geographical origin, and clinical stages and outcomes. Implementation of GMTV tracks the pattern of changes of M. tuberculosis strains in different geographical areas, facilitates disease gene discoveries associated with drug resistance or different clinical sequelae, and automates comparative genomic analyses among M. tuberculosis strains.

B chromosomes are dispensable and variable karyotypic elements found in some species of animals, plants and fungi. They often originate from duplications and translocations of host genomic regions or result from hybridization. In most species, little is known about their DNA content. Here we perform high-throughput sequencing and analysis of B chromosomes of roe deer and brocket deer, the only representatives of Cetartiodactyla known to have B chromosomes. In this study we developed an approach to identify genomic regions present on chromosomes by high-throughput sequencing of DNA generated from flow-sorted chromosomes using degenerate-oligonucleotide-primed PCR. Application of this method on small cattle autosomes revealed a previously described KIT gene region translocation associated with colour sidedness. Implementing this approach to B chromosomes from two cervid species, Siberian roe deer (Capreolus pygargus) and grey brocket deer (Mazama gouazoubira), revealed dramatically different genetic content: roe deer B chromosomes consisted of two duplicated genomic regions (a total of 1.42-1.98 Mbp) involving three genes, while grey brocket deer B chromosomes contained 26 duplicated regions (a total of 8.28-9.31 Mbp) with 34 complete and 21 partial genes, including KIT and RET protooncogenes, previously found on supernumerary chromosomes in canids. Sequence variation analysis of roe deer B chromosomes revealed a high frequency of mutations and increased heterozygosity due to either amplification within B chromosomes or divergence between different Bs. In contrast, grey brocket deer B chromosomes were found to be more homogeneous and resembled autosomes in patterns of sequence variation. Similar tendencies were observed in repetitive DNA composition. Our data demonstrate independent origins of B chromosomes in the grey brocket and roe deer. We hypothesize that the B chromosomes of these two cervid species represent different stages of B chromosome sequences evolution: probably nascent and similar to autosomal copies in brocket deer, highly derived in roe deer. Based on the presence of the same orthologous protooncogenes in canids and brocket deer Bs we argue that genomic regions involved in B chromosome formation are not random. In addition, our approach is also applicable to the characterization of other evolutionary and clinical rearrangements.

Russian Republic of Karelia is located at the Russian-Finnish border. It contains most of the historical Karelia land inhabited with autochthonous Karels and more recently migrated Russians. Although tuberculosis (TB) incidence in Karelia is decreasing, it remains high (45.8/100 000 in 2014) with the rate of multi-drug resistance (MDR) among newly diagnosed TB patients reaching 46.5 %. The study aimed to genetically characterize Mycobacterium tuberculosis isolates obtained at different time points from TB patients from Karelia to gain insight into the phylogeographic specificity of the circulating genotypes and to assess trends in evolution of drug resistant subpopulations. The sample included 150 M. tuberculosis isolates: 78 isolated in 2013–2014 (“new” collection) and 72 isolated in 2006 (“old” collection). Drug susceptibility testing was done by the method of absolute concentrations. DNA was subjected to spoligotyping and analysis of genotype-specific markers of the Latin-American-Mediterranean (LAM) family and its sublineages and Beijing B0/W148-cluster. The largest spoligotypes were SIT1 (Beijing family, n = 42) and SIT40 (T family, n = 5). Beijing family was the largest (n = 43) followed by T (n = 11), Ural (n = 10) and LAM (n = 8). Successful Russian clone, Beijing В0/W148, was identified in 15 (34.9 %) of 43 Beijing isolates; all В0/W148 isolates were drug-resistant. Seven of 8 LAM isolates belonged to the RD115/LAM-RUS branch, 1 - to the LAM RD174/RD-Rio sublineage. MDR was found in Beijing (32/43), Ural (3/10), and LAM (3/8). In contrast, all T isolates were pansusceptible. Comparison of drug resistant subgroups of the new and old collections showed an increasing prevalence of the B0/W148 clonal cluster, from 18.0 % (mainly polyresistant) in 2006 to 32.6 % in 2014 (mainly MDR and pre-XDR). The West–east increasing gradient is observed for the Ural genotype that may be defined a ‘Russian’ strain. In contrast, the spoligotype SIT40 of the T family appears to be a historical Karelian strain. Circulation of the MDR M. tuberculosis isolates of the Beijing genotype and its B0/W148 cluster continues to critically influence the current situation with the MDR-TB control in northwestern Russia including the Republic of Karelia. Revealed phylogeographic patterns of some genotypes reflect a complex demographic history of Karelia within the course of the 20th century.

In addition to the obligatory pathogenic species of the Mycobacterium tuberculosis complex and Mycobacterium leprae , the genus Mycobacterium also includes conditionally pathogenic species that in rare cases can lead to the development of nontuberculous mycobacterial diseases. Because tuberculosis and mycobacteriosis have similar clinical signs, the accurate identification of the causative agent in a clinical microbiology laboratory is important for diagnostic verification and appropriate treatment. This report describes a low-density hydrogel-based microarray containing oligonucleotide probes based on the species-specific sequences of the gyrB gene fragment for mycobacterial species identification. The procedure included the amplification of a 352-nucleotide fragment of the gene and its hybridization on a microarray. The triple-species-specific probe design and the algorithm for hybridization profile recognition based on the calculation of Pearson correlation coefficients, followed by the construction of a profile database, allowed for the reliable and accurate identification of mycobacterial species, including mixed-DNA samples. The assay was used to evaluate 543 clinical isolates from two regions of Russia, demonstrating its ability to detect 35 mycobacterial species, with 99.8% sensitivity and 100% specificity when using gyrB , 16S, and internal transcribed spacer (ITS) fragment sequencing as the standard. The testing of clinical samples showed that the sensitivity of the assay was 89% to 95% for smear-positive samples and 36% for smear-negative samples. The large number of identified species, the high level of sensitivity, the ability to detect mycobacteria in clinical samples, and the up-to-date profile database make the assay suitable for use in routine laboratory practice.

B chromosomes (Bs) represent a variable addition to the main karyotype in some lineages of animals and plants. Bs accumulate through non-Mendelian inheritance and become widespread in populations. Despite the presence of multiple genes, most Bs lack specific phenotypic effects, although their influence on host genome epigenetic status and gene expression are recorded. Previously, using sequencing of isolated Bs of ruminants and rodents, we demonstrated that Bs originate as segmental duplications of specific genomic regions, and subsequently experience pseudogenization and repeat accumulation. Here, we used a similar approach to characterize Bs of the red fox (Vulpes vulpes L.) and the Chinese raccoon dog (Nyctereutes procyonoides procyonoides Gray). We confirm the previous findings of the KIT gene on Bs of both species, but demostrate an independent origin of Bs in these species, with two reused regions. Comparison of gene ensembles in Bs of canids, ruminants, and rodents once again indicates enrichment with cell-cycle genes, development-related genes, and genes functioning in the neuron synapse. The presence of B-chromosomal copies of genes involved in cell-cycle regulation and tissue differentiation may indicate importance of these genes for B chromosome establishment.