Mouse chromocenters DNA content: sequencing and in silico analysis.

BACKGROUND: Chromocenters are defined as a punctate condensed blocks of chromatin in the interphase cell nuclei of certain cell types with unknown biological significance. In recent years a progress in revealing of chromocenters protein content has been made although the details of DNA content within constitutive heterochromatin still remain unclear. It is known that these regions are enriched in tandem repeats (TR) and transposable elements. Quick improvement of genome sequencing does not help to assemble the heterochromatic regions due to lack of appropriate bioinformatics techniques. RESULTS: Chromocenters DNA have been isolated by a biochemical approach from mouse liver cells nuclei and sequenced on the Illumina MiSeq resulting in ChrmC dataset. Analysis of ChrmC dataset by the bioinformatics tools available revealed that the major component of chromocenter DNA are TRs: ~ 66% MaSat and ~ 4% MiSat. Other previously classified TR families constitute ~ 1% of ChrmC dataset. About 6% of chromocenters DNA are mostly unannotated sequences. In the contigs assembled with IDBA_UD there are many fragments of heterochromatic Y-chromosome, rDNA and other pseudo-genes and non-coding DNA. A protein coding sfi1 homolog gene fragment was also found in contigs. The Sfi1 homolog gene is located on the chromosome 11 in the reference genome very close to the Golden Pass Gap (a ~ 3 Mb empty region reserved to the pericentromeric region) and proves the purity of chromocenters isolation. The second major fraction are non-LTR retroposons (SINE and LINE) with overwhelming majority of LINE - ~ 11% of ChrmC. Most of the LINE fragments are from the ~ 2 kb region at the end of the 2nd ORF and its' flanking region. The precise LINEs' segment of ~ 2 kb is the necessary mouse constitutive heterohromatin component together with TR. The third most abundant fraction are ERVs. The ERV distribution in chromocenters differs from the whole genome: IAP (ERV2 class) is the most numerous in ChrmC while MaLR (ERV3 class) prevails in the reference genome. IAP and its LTR also prevail in TR containing contigs extracted from the WGS dataset. In silico prediction of IAP and LINE fragments in chromocenters was confirmed by direct fluorescent in situ hybridization (FISH). CONCLUSION: Our data of chromocenters' DNA (ChrmC) sequencing demonstrate that IAP with LTR and a precise ~ 2 kb fragment of LINE represent a substantial fraction of mouse chromocenters (constitutive heteroсhromatin) along with TRs.

Contrasting origin of B chromosomes in two cervids (Siberian roe deer and grey brocket deer) unravelled by chromosome-specific DNA sequencing.

BACKGROUND: 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. RESULTS: 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. CONCLUSIONS: 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.

Evaluation of a low-density hydrogel microarray technique for mycobacterial species identification.

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.

Genome-wide Mycobacterium tuberculosis variation (GMTV) database: a new tool for integrating sequence variations and epidemiology.

BACKGROUND: 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. DESCRIPTION: 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. CONCLUSIONS: 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.

SARS-CoV-2 Recombination and Coinfection Events Identified in Clinical Samples in Russia.

Recombination is one of the mechanisms of SARS-CoV-2 evolution along with the occurrence of point mutations, insertions, and deletions. Recently, recombinant variants of SARS-CoV-2 have been registered in different countries, and some of them have become circulating forms. In this work, we performed screening of SARS-CoV-2 genomic sequences to identify recombination events and co-infections with various strains of the SARS-CoV-2 virus detected in Russia from February 2020 to March 2022. The study included 9336 genomes of the COVID-19 pathogen obtained as a result of high-throughput sequencing on the Illumina platform. For data analysis, we used an algorithm developed by our group that can identify viral recombination variants and cases of co-infections by estimating the frequencies of characteristic substitutions in raw read alignment files and VCF files. The detected cases of recombination were confirmed by alternative sequencing methods, principal component analysis, and phylogenetic analysis. The suggested approach allowed for the identification of recombinant variants of strains BA.1 and BA.2, among which a new recombinant variant was identified, as well as a previously discovered one. The results obtained are the first evidence of the spread of recombinant variants of SARS-CoV-2 in Russia. In addition to cases of recombination we identified cases of coinfection: eight of them contained the genome of the Omicron line as one of the variants, six of them the genome of the Delta line, and two with the genome of the Alpha line.

Genomic Variability of SARS-CoV-2 Omicron Variant Circulating in the Russian Federation during Early December 2021 and Late January 2022.

Analysis of genomic variability of pathogens associated with heightened public health concerns is an opportunity to track transmission routes of the disease and helps to develop more effective vaccines and specific diagnostic tests. We present the findings of a detailed genomic analysis of the genomic variability of the SARS-CoV-2 Omicron variant that spread in Russia between 8 December 2021 and 30 January 2022. We performed phylogenetic analysis of Omicron viral isolates collected in Moscow (n = 589) and downloaded from GISAID (n = 397), and identified that the BA.1 lineage was predominant in Russia during this period. The BA.2 lineage was also identified early in December 2021. We identified three cases of BA.1/BA.2 coinfections and one case of Delta/Omicron coinfection. A comparative genomic analysis of SARS-CoV-2 viral variants that spread in other countries allowed us to identify possible cases of transmission. We also found that some mutations that are quite rare in the Global Omicron dataset have a higher incidence rate, and identified genetic markers that could be associated with ways of Omicron transmission in Russia. We give the genomic variability of single nucleotide variations across the genome and give a characteristic of haplotype variability of Omicron strains in both Russia and around the world, and we also identify them.

Sequencing of Supernumerary Chromosomes of Red Fox and Raccoon Dog Confirms a Non-Random Gene Acquisition by B Chromosomes.

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.