Metagenomic profiling of viral and microbial communities from the pox lesions of lumpy skin disease virus and sheeppox virus-infected hosts

Introduction

It has been recognized that capripoxvirus infections have a strong cutaneous tropism with the manifestation of skin lesions in the form of nodules and scabs in the respective hosts, followed by necrosis and sloughing off. Considering that the skin microbiota is a complex community of commensal bacteria, fungi and viruses that are influenced by infections leading to pathological states, there is no evidence on how the skin microbiome is affected during capripoxvirus pathogenesis.

Methods

In this study, shotgun metagenomic sequencing was used to investigate the microbiome in pox lesions from hosts infected with lumpy skin disease virus and sheep pox virus.

Results

The analysis revealed a high degree of variability in bacterial community structures across affected skin samples, indicating the importance of specific commensal microorganisms colonizing individual hosts. The most common and abundant bacteria found in scab samples were Fusobacterium necrophorum, Streptococcus dysgalactiae, Helcococcus ovis and Trueperella pyogenes, irrespective of host. Bacterial reads belonging to the genera Moraxella, Mannheimia, Corynebacterium, Staphylococcus and Micrococcus were identified.

Discussion

This study is the first to investigate capripox virus-associated changes in the skin microbiome using whole-genome metagenomic profiling. The findings will provide a basis for further investigation into capripoxvirus pathogenesis. In addition, this study highlights the challenge of selecting an optimal bioinformatics approach for the analysis of metagenomic data in clinical and veterinary practice. For example, direct classification of reads using a kmer-based algorithm resulted in a significant number of systematic false positives, which may be attributed to the peculiarities of the algorithm and database selection. On the contrary, the process of de novo assembly requires a large number of target reads from the symbiotic microbial community. In this work, the obtained sequencing data were processed by three different approaches, including direct classification of reads based on k-mers, mapping of reads to a marker gene database, and de novo assembly and binning of metagenomic contigs. The advantages and disadvantages of these techniques and their practicality in veterinary settings are discussed in relation to the results obtained.

Koban culture genome-wide and archeological data open the bridge between Bronze and Iron Ages in the North Caucasus

The North Caucasus played a key role during the ancient colonization of Eurasia and the formation of its cultural and genetic ancestry. Previous archeogenetic studies described a relative genetic and cultural continuity of ancient Caucasus societies, since the Eneolithic period. The Koban culture, which formed in the Late Bronze Age on the North Caucasian highlands, is considered as a cultural "bridge" between the ancient and modern autochthonous peoples of the Caucasus. Here, we discuss the place of this archeological culture and its representatives in the genetic orbit of Caucasian cultures using genome-wide SNP data from five individuals of the Koban culture and one individual of the early Alanic culture as well as previously published genomic data of ancient and modern North Caucasus individuals. Ancient DNA analysis shows that an ancient individual from Klin-Yar III, who was previously described as male, was in fact a female. Additional studies on well-preserved ancient human specimens are necessary to determine the level of local mobility and kinship between individuals in ancient societies of North Caucasus. Further studies with a larger sample size will allow us gain a deeper understanding of this topic.

Search for differentially methylated regions in ancient and modern genomes

Currently, active research is focused on investigating the mechanisms that regulate the development of various pathologies and their evolutionary dynamics. Epigenetic mechanisms, such as DNA methylation, play a significant role in evolutionary processes, as their changes have a faster impact on the phenotype compared to mutagenesis. In this study, we attempted to develop an algorithm for identifying differentially methylated regions associated with metabolic syndrome, which have undergone methylation changes in humans during the transition from a hunter-gatherer to a sedentary lifestyle. The application of existing whole-genome bisulfite sequencing methods is limited for ancient samples due to their low quality and fragmentation, and the approach to obtaining DNA methylation profiles differs significantly between ancient hunter-gatherer samples and modern tissues. In this study, we validated DamMet, an algorithm for reconstructing ancient methylomes. Application of DamMet to Neanderthal and Denisovan genomes showed a moderate level of correlation with previously published methylation profiles and demonstrated an underestimation of methylation levels in the reconstructed profiles by an average of 15-20 %. Additionally, we developed a new Python-based algorithm that allows for the comparison of methylomes in ancient and modern samples, despite the absence of methylation profiles in modern bone tissue within the context of obesity. This analysis involves a two-step data processing approach, where the first step involves the identification and filtration of tissue-specific methylation regions, and the second step focuses on the direct search for differentially methylated regions in specific areas associated with the researcher's target condition. By applying this algorithm to test data, we identified 38 differentially methylated regions associated with obesity, the majority of which were located in promoter regions. The pipeline demonstrated sufficient efficiency in detecting these regions. These results confirm the feasibility of reconstructing DNA methylation profiles in ancient samples and comparing them with modern methylomes. Furthermore, possibilities for further methodological development and the implementation of a new step for studying differentially methylated positions associated with evolutionary processes are discussed.

Mid-holocene Brown Bear (Ursus arctos) from the Bolshoy Lyakhovsky Island (New Siberian Islands)

A morphological description is provided for a unique find of a frozen mummified subfossil brown bear (Ursus arctos L., 1758), found for the first time ever. The find is a well-preserved bear carcass of approximately 3500 years in age. Results of computed tomography and DNA testing are discussed.

The Rurikids: The First Experience of Reconstructing the Genetic Portrait of the Ruling Family of Medieval Rus' Based on Paleogenomic Data

The Rurikids were the reigning house of Rus', its principalities and, ultimately the Tsardom of Russia, for seven centuries: from the IX to the end of the XVI century. According to the Primary Chronicle (the Tale of Bygone Years), the main chronicle of Rus', the Rurik dynasty was founded by the Varangian prince Rurik, invited to reign in Novgorod in 862, but still there is no direct genetic evidence of the origin of the early Rurikids. This research, for the first time, provides a genome-wide paleogenetic analysis of bone remains belonging to one of the Rurikids, Prince Dmitry Alexandrovich (?-1294), the son of the Grand Prince of Vladimir Alexander Yaroslavich Nevsky (1221-1263). It has been established that his Y chromosome belongs to the N1a haplogroup. Most of the modern Rurikids, according to their genealogies, belonging to the N1a haplogroup, have the most similar variants of Y chromosomes to each other, as well as to the Y chromosome of Prince Dmitry Alexandrovich. Genome-wide data of the medieval and modern Rurikids unequivocally indicates that they belong to the N1a haplogroup of the Y chromosome, starting at least from the XI century (since the time of Prince Yaroslav the Wise). All the other alleged Rurikids, both ancient and modern, being carriers of other haplogroups (R1a, I2a), possess high heterogeneity of the sequence of Y chromosomes, meaning that we cannot confirm their common ancestry. The most probable ancestors of Prince Dmitry Alexandrovich in the male line were the men who left the burial ground Bolshoy Oleny Island on the coast of the Kola Peninsula about 3,600 years ago. The reconstruction of the genome of Prince Dmitry Alexandrovich indicates the contribution of three ancestral components to his origin: (1) the early medieval population of the east of Scandinavia from the island of Oland, (2) representatives of the steppe nomadic peoples of the Eurasian steppes of the Iron Age or the early medieval population of central Europe (steppe nomads from the territory of Hungary), and (3) the ancient East-Eurasian component. Reliable statistics were also obtained when the Scandinavians were replaced with the Medieval Russian Slavic populations of the XI century. Thus, for the first time, we have shown the complex nature of interethnic interactions in the formation of the nobility of medieval Rus' on the example of the ancient Rurikid.

SNP-Based Analysis Reveals Authenticity and Genetic Similarity of Russian Indigenous V. vinifera Grape Cultivars

9 Russian Vitis vinifera grape varieties and the European variety Muscat Hamburg were sequenced and genotyped using 527 SNPs (single nucleotide polymorphisms) with high minor allele frequency for the first time. The data were coupled with previously identified genotypes of 783 varieties and subjected to parentage and population analysis. As a result, contrary to the historical and ampelographic data published in many sources from 1800 to 2012, only two of the nine Russian varieties (Pukhlyakovskiy Belyi and Sibirkovyi) were related to foreign ones and were obviously imported from Europe to the Russian Empire. The remaining seven varieties, led by Krasnostop Zolotovskiy, are not directly related either in the Caucasus or in Europe, they form separate clusters on the genetic distance-based dendrogram and the world parentage network of V. vinifera. The resulting pedigree of Muscat Hamburg and its descendants is in accordance with SSR-based (simple sequence repeats) studies and the described pedigree of this variety which confirms the use of the reduced SNP set for further studies.

Caucasian treasure: Genomics sheds light on the evolution of half-extinct Sevan trout, Salmo ischchan, species flock

Five ecologically and phenotypically divergent ecomorphs of the genus Salmo are known from a landlocked alpine lake in the Caucasus, Lake Sevan. It is an example of sympatric diversification within a species-rich lineage with predominate mode of speciation being allopatric. The diversification of Sevan trouts was accompanied by spawning resource partitioning. Four lacustrine ecomorphs with different temporal-spatial spawning strategies and divergent morphology and coloration evolved along with a fifth ecomorph, brook trout, inhabiting the tributaries. Unfortunately, the Sevan trout diversity was almost destroyed by human activity, with two ecomorphs becoming extinct in the 1980s. We performed reconstruction of the evolutionary history of Sevan trouts based on high-throughput sequencing of both contemporary and historical DNA (∼ 50 y.o.) of all Sevan trout ecomorphs. Our study of complete mitogenomes along with genome-wide SNP data revealed the monophyly of four lacustrine ecomorphs and local brook trout, all derived from the anadromous form Caspian salmon, S. caspius. The species tree suggests a scenario of stepwise evolution from riverine to lacustrine spawning. Three genomic clusters were revealed, of which two refer to the riverine and lacustrine spawners within the flock of Sevan trouts (with F_ST value = 0.069). A few SNP outliers under selection were discovered that could be responsible for assortative mating based on visual recognition. The Holocene climatic oscillations and the desiccation of tributaries could have played an important role in the origin of lacustrine spawning. The relationships between lacustrine ecomorphs were not yet fully resolved. This radiation warrants further investigation.

Steller's sea cow genome suggests this species began going extinct before the arrival of Paleolithic humans

Anthropogenic activity is the top factor directly related to the extinction of several animal species. The last Steller's sea cow (Hydrodamalis gigas) population on the Commander Islands (Russia) was wiped out in the second half of the 18^th century due to sailors and fur traders hunting it for the meat and fat. However, new data suggests that the extinction process of this species began much earlier. Here, we present a nuclear de novo assembled genome of H. gigas with a 25.4× depth coverage. Our results demonstrate that the heterozygosity of the last population of this animal is low and comparable to the last woolly mammoth population that inhabited Wrangel Island 4000 years ago. Besides, as a matter of consideration, our findings also demonstrate that the extinction of this marine mammal starts along the North Pacific coastal line much earlier than the first Paleolithic humans arrived in the Bering sea region.

Molecular phylogeny of one extinct and two critically endangered Central Asian sturgeon species (genus Pseudoscaphirhynchus) based on their mitochondrial genomes

The enigmatic and poorly studied sturgeon genus Pseudoscaphirhynchus (Scaphirhynchinae: Acipenseridae) comprises three species: the Amu Darya shovelnose sturgeon (Pseudoscaphirhynchus kaufmanni (Bogdanow)), dwarf Amu Darya shovelnose sturgeon P. hermanni (Kessler), and Syr Darya shovelnose sturgeon (P. fedtschenkoi (Bogdanow). Two species - P. hermanni and P. kaufmanni - are critically endangered due to the Aral Sea area ecological disaster, caused by massive water use for irrigation to support cotton agriculture, subsequent pesticide pollution and habitat degradation. For another species - P. fedtschenkoi - no sightings have been reported since 1960-s and it is believed to be extinct, both in nature and in captivity. In this study, complete mitochondrial (mt) genomes of these three species of Pseudoscaphirhynchus were characterized using Illumina and Sanger sequencing platforms. Phylogenetic analyses showed the significant divergence between Amu Darya and Syr Darya freshwater sturgeons and supported the monophyletic origin of the Pseudoscaphirhynchus species. We confirmed that two sympatric Amu Darya species P. kaufmanni and P. hermanni form a single genetic cluster, which may require further morphological and genetic study to assess possible hybridization, intraspecific variation and taxonomic status and to develop conservation measures to protect these unique fishes.

A partial genome assembly of the miniature parasitoid wasp, Megaphragma amalphitanum

Body size reduction, also known as miniaturization, is an important evolutionary process that affects a number of physiological and phenotypic traits and helps animals conquer new ecological niches. However, this process is poorly understood at the molecular level. Here, we report genomic and transcriptomic features of arguably the smallest known insect-the parasitoid wasp, Megaphragma amalphitanum (Hymenoptera: Trichogrammatidae). In contrast to expectations, we find that the genome and transcriptome sizes of this parasitoid wasp are comparable to other members of the Chalcidoidea superfamily. Moreover, compared to other chalcid wasps the gene content of M. amalphitanum is remarkably conserved. Intriguingly, we observed significant changes in M. amalphitanum transposable element dynamics over time, in which an initial burst was followed by suppression of activity, possibly due to a recent reinforcement of the genome defense machinery. Overall, while the M. amalphitanum genomic data reveal certain features that may be linked to the unusual biological properties of this organism, miniaturization is not associated with a large decrease in genome complexity.

The complete mitochondrial genome of the extinct Pleistocene horse (Equus cf. lenensis) from Kotelny Island (New Siberian Islands, Russia) and its phylogenetic assessment

The complete mitochondrial genome from the Pleistocene stallion horse (Equus cf. lenensis) which complete skull was found in 1901 on Kotelny Island (New Siberian Archipelago, Sakha Republic, Russia) is published in this paper. The mitochondrial DNA (mtDNA) is 16,584 base pairs (bp) in length and contained 13 protein-coding genes, 2 rRNA genes, 22 tRNA genes. The overall base composition of the genome in descending order was 32.3% - A, 28.5% - C, 13.4% - G, 25.8% - T without a significant AT bias of 58.2%.

Phylogenetic position of the presumably extinct slender-billed curlew, Numenius tenuirostris

The high-capacity DNA analysis of museum samples opens new opportunities, associated with the investigation of extinct species evolution. Here, the complete mitochondrial genome of the presumably extinct bird species, the slender-billed curlew Numenius tenuirostris (Charadriiformes: Scolopacidae) is presented. Our results showed that mitochondrial DNA (mtDNA) is 16,705 base pairs (bp) in length and contain 13 protein-coding genes, two rRNA genes, and 22 tRNA genes. The overall base composition of the genome is 30.8% - A, 29.8% - C, 25.4% - T, 14.0% - G, and without a significant GC bias of 43.7%. Phylogenetic analyses based on the cytochrome B (cytB) gene and the whole mtDNA sequences revealed that N. tenuirostris had a close genetic relationship to Eurasian curlew (N. arquata), Far Eastern curlew (N. madagascariensis), and long-billed curlew - N. americanus. Besides, it reveals that Numenius genus is genetically distant from other Scolopacidae taxons. Together, these results provide a clear genetic perspective into the speciation process among the curlew genus members and points to a clear taxonomic position of N. tenuirostris.

Sequencing of two mitochondrial genomes of endangered form of the Sevan trout Salmo ischchan aestivalis

The two complete mitochondrial genomes of endangered form of the Sevan trout Salmo ischchan aestivalis are published in this paper. The mitochondrial DNA (mtDNA) is 16,677 base pairs (bp) in length and contained 13 protein-coding genes, 2 rRNA genes, and 22 tRNA genes. The overall base composition of the genome in descending order was 29.4% - C, 27.9% - A, 26.0% - T, 16.7% - G, without a significant AT bias of 53.9%.

Two complete mitochondrial genomes of extinct form of the Sevan trout Salmo ischchan danilewskii

The mitochondrial genomes from two individuals of the extinct subspecies of the Sevan trout Salmo ischchan danilewskii are published in this paper. The mitochondrial DNA (mtDNA) is 16,665 base pairs (bp) in length and contained 13 protein-coding genes, 2 rRNA genes, and 22 tRNA genes. The overall base composition of the genome in descending order was 27.9% of A, 29.4% of C, 16.7% of G, and 26.0% of T without a significant AT bias of 53.9%.

[Hsp70 Genes of the Megaphragma amalphitanum (Hymenoptera: Trichogrammatidae) Parasitic Wasp]

Miniaturization is an evolutionary process that is widely represented in both invertebrates and vertebrates. Miniaturization frequently affects not only the size of the organism and its constituent cells, but also changes the genome structure and functioning. The structure of the main heat shock genes (hsp70 and hsp83) was studied in one of the smallest insects, the Megaphragma amalphitanum (Hymenoptera: Trichogrammatidae) parasitic wasp, which is comparable in size with unicellular organisms. An analysis of the sequenced genome has detected six genes that relate to the hsp70 family, some of which are apparently induced upon heat shock. Both induced and constitutively expressed hsp70 genes contain a large number of introns, which is not typical for the genes of this family. Moreover, none of the found genes form clusters, and they are all very heterogeneous (individual copies are only 75-85% identical), which indicates the absence of gene conversion, which provides the identity of genes of this family in Drosophila and other organisms. Two hsp83 genes, one of which contains an intron, have also been found in the M. amalphitanum genome.

Metagenomic analysis of microbial community of a parasitoid wasp Megaphragma amalphitanum

The vast majority of multicellular organisms coexist with bacterial symbionts that may play various roles during their life cycle. Parasitoid wasp Megaphragma amalphitanum (Hymenoptera: Trichogrammatidae) belongs to the smallest known insects whose size is comparable with some bacteria. Using 16S rRNA gene sequencing and Whole Genome Sequencing (WGS), we described microbiota diversity for this arthropod and its potential impact on their lifecycle. Metagenomic sequences were deposited to SRA database which is available at NCBI with accession number SRX2363723 and SRX2363724. We found that small body size and limited lifespan do not lead to a significant reduction of bacterial symbionts diversity. At the same time, we show here a specific feature of microbiota composition in M. amalphitanum – the absence of the Rickettsiaceae family representatives that are known to cause sex-ratio distortion in arthropods and well represented in other populations of parasitoid wasps.

Identification of novel microRNA genes in freshwater and marine ecotypes of the three-spined stickleback (Gasterosteus aculeatus)

The three-spined stickleback (Gasterosteus aculeatus L.) is an important model organism for studying the molecular mechanisms of speciation and adaptation to salinity. Despite increased interest to microRNA discovery and recent publication on microRNA prediction in the three-spined stickleback using bioinformatics approaches, there is still a lack of experimental support for these data. In this paper, high-throughput sequencing technology was applied to identify microRNA genes in gills of the three-spined stickleback. In total, 595 miRNA genes were discovered; half of them were predicted in previous computational studies and were confirmed here as microRNAs expressed in gill tissue. Moreover, 298 novel microRNA genes were identified. The presence of miRNA genes in selected 'divergence islands' was analysed and 10 miRNA genes were identified as not randomly located in 'divergence islands'. Regulatory regions of miRNA genes were found enriched with selective SNPs that may play a role in freshwater adaptation.

Mitochondrial genome of Megaphragma amalphitanum (Hymenoptera: Trichogrammatidae)

The mitochondrial genome of the parasitic wasp Megaphragma amalphitanum is published in this paper. The mitochondrial DNA (mtDNA) is 15 041 base pairs (bp) in length and contained 13 protein-coding genes, 2 rRNA genes and 22 tRNA genes. The overall base composition of the genome in descending order was 45.67% A, 8.71% C, 39.63% T and 5.99% G, with a significant AT bias of 85.30%.