Genomic analysis reveals cryptic diversity in aphelids and sheds light on the emergence of Fungi.

Over the past decade, molecular phylogenetics has reshaped our understanding of the fungal tree of life by unraveling a hitherto elusive diversity of the protistan relatives of Fungi. Aphelida constitutes one of these novel deep branches that precede the emergence of osmotrophic fungal lifestyle and hold particular significance as the pathogens of algae. Here, we obtain and analyze the genomes of aphelid species Amoeboaphelidium protococcarum and Amoeboaphelidium occidentale. Genomic data unmask the vast divergence between these species, hidden behind their morphological similarity, and reveal hybrid genomes with a complex evolutionary history in two strains of A. protococcarum. We confirm the proposed sister relationship between Aphelida and Fungi using phylogenomic analysis and chart the reduction of characteristic proteins involved in phagocytic activity in the evolution of Holomycota. Annotation of aphelid genomes demonstrates the retention of actin nucleation-promoting complexes associated with phagocytosis and amoeboid motility and also reveals a conspicuous expansion of receptor-like protein kinases, uncharacteristic of fungal lineages. We find that aphelids possess multiple carbohydrate-processing enzymes that are involved in fungal cell wall synthesis but do not display rich complements of algal cell-wall-processing enzymes, suggesting an independent origin of fungal plant-degrading capabilities. Aphelid genomes show that the emergence of Fungi from phagotrophic ancestors relied on a common cell wall synthetic machinery but required a different set of proteins for digestion and interaction with the environment.

Comparative analysis of alignment-free genome clustering and whole genome alignment-based phylogenomic relationship of coronaviruses.

The SARS-CoV-2 is the third coronavirus in addition to SARS-CoV and MERS-CoV that causes severe respiratory syndrome in humans. All of them likely crossed the interspecific barrier between animals and humans and are of zoonotic origin, respectively. The origin and evolution of viruses and their phylogenetic relationships are of great importance for study of their pathogenicity and development of antiviral drugs and vaccines. The main objective of the presented study was to compare two methods for identifying relationships between coronavirus genomes: phylogenetic one based on the whole genome alignment followed by molecular phylogenetic tree inference and alignment-free clustering of triplet frequencies, respectively, using 69 coronavirus genomes selected from two public databases. Both approaches resulted in well-resolved robust classifications. In general, the clusters identified by the first approach were in good agreement with the classes identified by the second using K-means and the elastic map method, but not always, which still needs to be explained. Both approaches demonstrated also a significant divergence of genomes on a taxonomic level, but there was less correspondence between genomes regarding the types of diseases they caused, which may be due to the individual characteristics of the host. This research showed that alignment-free methods are efficient in combination with alignment-based methods. They have a significant advantage in computational complexity and provide valuable additional alternative information on the genomes relationships.

The Mitochondrial Genome of a Freshwater Pelagic Amphipod Macrohectopus branickii Is among the Longest in Metazoa.

There are more than 350 species of amphipods (Crustacea) in Lake Baikal, which have emerged predominantly through the course of endemic radiation. This group represents a remarkable model for studying various aspects of evolution, one of which is the evolution of mitochondrial (mt) genome architectures. We sequenced and assembled the mt genome of a pelagic Baikalian amphipod species Macrohectopus branickii. The mt genome is revealed to have an extraordinary length (42,256 bp), deviating significantly from the genomes of other amphipod species and the majority of animals. The mt genome of M. branickii has a unique gene order within amphipods, duplications of the four tRNA genes and Cox2, and a long non-coding region, that makes up about two thirds of the genome's size. The extension of the mt genome was most likely caused by multiple duplications and inversions of regions harboring ribosomal RNA genes. In this study, we analyzed the patterns of mt genome length changes in amphipods and other animal phyla. Through a statistical analysis, we demonstrated that the variability in the mt genome length may be a characteristic of certain phyla and is primarily conferred by expansions of non-coding regions.

The diversity of the Baikal lineage of Hydra oligactis Pallas, 1766: molecular and morphological evidence.

In this paper, molecular analyses of Baikal hydras from the 'oligactis group', based on COI and ITS1-5.8S-ITS2, and morphological analysis of their holotrichous isorhizas, were performed. Low genetic diversity and shared haplotypes were found between Hydra oligactis Pallas, 1766 and Hydra baikalensis Swarczewsky, 1923 specimens, which is evidence of the mixing of these lineages. Genetic distances among all Baikal hydras (0.006) were less than the interspecific distances of other hydras. The size of hydras and proportions of their holotrichous isorhizas varied depending on microhabitat and environmental conditions. Our combined molecular and morphological approach proves that H. baikalensis is synonymous with H. oligactis.

Bioinformatic tools for tRNA gene analyses in mitochondrial DNA sequence data.

The data presented here are related to the research article entitled "Hidden cases of tRNA genes duplication and remolding in mitochondrial genomes of amphipods" (Romanova et al., 2020) [1]. Correct tRNA gene sequence annotation in mitochondrial (mt) and nuclear genomes sometimes can be a challenging task because of the differential performances of tRNA annotation/prediction programmes. These programmes may cause false positive or false negative predictions. Moreover, additional difficulties with annotation may be caused by the presence of duplicated tRNA genes and those coding tRNAs with altered identities occurring as due to a mutation in their anticodon sequence (tRNA gene remolding/recruitment). We developed an R script automating the diagnosis of ancestor tRNA gene coding specificity regardless of anticodon sequence based on genetic distance comparison. Some of the predicted tRNA genes from the mt genomes of amphipods are presented. We also developed an R script for estimation of the best mode of sequence alignment, which was applied to determine the best alignment of tRNA genes in [1], but is also suitable for testing of any nucleotide alignment sets used in phylogenetic inferences.

Hidden cases of tRNA gene duplication and remolding in mitochondrial genomes of amphipods.

The evolution of tRNA genes in mitochondrial (mt) genomes is a complex process that includes duplications, degenerations, and transpositions, as well as a specific process of identity change through mutations in the anticodon (tRNA gene remolding or tRNA gene recruitment). Using amphipod-specific tRNA models for annotation, we show that tRNA duplications are more common in the mt genomes of amphipods than what was revealed by previous annotations. Seventeen cases of tRNA gene duplications were detected in the mt genomes of amphipods, and ten of them were tRNA genes that underwent remolding. The additional tRNA gene findings were verified using phylogenetic analysis and genetic distance analysis. The majority of remolded tRNA genes (seven out of ten cases) were found in the mt genomes of endemic amphipod species from Lake Baikal. All additional mt tRNA genes arose independently in the Baikalian amphipods, indicating the unusual plasticity of tRNA gene evolution in these species assemblages. The possible reasons for the unusual abundance of additional tRNA genes in the mt genomes of Baikalian amphipods are discussed. The amphipod-specific tRNA models developed for MiTFi refine existing predictions of tRNA genes in amphipods and reveal additional cases of duplicated tRNA genes overlooked by using less specific Metazoa-wide models. The application of these models for mt tRNA gene prediction will be useful for the correct annotation of mt genomes of amphipods and probably other crustaceans.

"886-84-like" tick-borne encephalitis virus strains: Intraspecific status elucidated by comparative genomics.

Tick-borne encephalitis virus (TBEV) can cause severe meningitis, encephalitis, and meningoencephalitis. TBEV represents a pathogen of high zoonotic potential and an emerging global threat. There are three known subtypes of TBEV: Far-Eastern, Siberian and European. Since 2001 there have been suggestions that two new subtypes may be distinguished: "178-79" and "886-84". These assumptions are based on the results of the envelope gene fragment sequencing (Zlobin et al., 2001; Kovalev and Mukhacheva, 2017) and genotype-specific probes molecular hybridization (Demina et al., 2010). There is only one full-genome sequence of "178-79" strain and two identical ones of "886-84" strain can be found in GenBank. For clarification of the intraspecific position of the "886-84-like" strains group we completely sequenced six previously unknown "886-84-like" strains isolated in Eastern Siberia. As a result of applying different bioinformatics approaches, we can confirm that "886-84-like" strains group is a distinct subtype of TBEV.

Population genetic structure and phylogeography of sterlet (Acipenser ruthenus, Acipenseridae) in the Ob and Yenisei river basins.

The sterlet (Acipenser ruthenus Linnaeus, 1758) is a relatively small sturgeon widely distributed in Eurasian rivers from the Danube to the Yenisei. During the twentieth century, all wild sterlet populations have declined due to anthropogenic factors including: overfishing, poaching, construction of dams, and pollution. Despite the necessity of characterization both wild and captive stocks, few studies of population genetics have been performed thus far. Here we studied the genetic diversity and geographic structure of sterlet populations across the eastern range - Ob-Irtysh and Yenisei basins - by sequencing a 628-bp fragment of mitochondrial DNA control region. We identified 98 new haplotypes, delineated 12 haplogroups and estimated the time of basal haplogroup divergence within the species as over 8 million years ago. Our data suggest that Ob-Irtysh and Yenisei populations are isolated from each other and much lower genetic diversity is present in the Yenisei population than in the Ob-Irtysh population. Our data imply that sterlet populations in Siberian rivers underwent bottleneck or fragmentation, followed by subsequent population expansion. The data obtained here are important for sterlet population monitoring and restocking management.

A new, genetically divergent species of Pseudobaikalia Lindholm, 1909 (Caenogastropoda, Baicaliidae).

A new gastropod species, Pseudobaikalia michelae Sitnikoiva & Kovalenkova, sp. n., (family Baicaliidae) is described from Lake Baikal. This is the first new species from the Baicaliidae for forty years. The new species is distinguished from its sister taxa by means of comparative morphology as well as analyses of DNA sequences (mtDNA and an intron of alpha-subunit gene of ATP-synthase). It was found in the southern and central-eastern parts of the lake where it occurs sympatrically with three other baicaliid species. Characters of the female reproductive system (i.e., a long oviduct loop with 2-3 narrow tube-like evaginations) and the aperture (i.e., oval shape with a simple outer lip) place these snails in the genus Pseudobaikalia Lindholm, 1909. The new species is most similar in its shell morphology to the northern Baikal species Pseudobaikalia jentteriana (smooth elongated shape) but differs by a more oval aperture that is slightly angled to the columella. Combined mitochondrial and nuclear sequences in a Bayesian analysis showed that all specimens of Pseudobaikalia michelae sp. n. form a well-supported clade.

The complete mitochondrial genome of a deep-water Baikalian amphipoda Brachyuropus grewingkii (Dybowsky, 1874).

In this study, we present a complete mitochondrial genome of a deep-water amphipoda Brachyuropus grewingkii (Dybowsky, 1874) from Lake Baikal. A circular mitochondrial DNA has 17,118 bp in length and contains 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes, a putative control region, and five intergenic spacers. An extended control region and altered positions of some tRNA genes distinguish mitochondrial genome of B. grewingkii from the mitochondrial genomes described for other Baikalian amphipoda species.

The complete mitochondrial genome of Baikalian amphipoda Eulimnogammarus vittatus Dybowsky, 1874.

A complete mitochondrial genome sequence of amphipoda Eulimnogammarus vittatus Dybowsky, 1874 from Lake Baikal was obtained using next-generation sequencing approach. Mitochondrial DNA with the length of 15,534 bp contains 13 protein-coding genes, 2 ribosomal RNA, 23 transfer RNA and non-coding sequences: a putative control region and 7 intergenic spacers. A brief comparative analysis of mitochondrial genomes of E. vittatus and its sister species Eulimnogammarus verrucosus was performed.

Evolution of mitochondrial genomes in Baikalian amphipods.

BACKGROUND: Amphipods (Crustacea) of Lake Baikal are a very numerous and diverse group of invertebrates generally believed to have originated by adaptive radiation. The evolutionary history and phylogenetic relationships in Baikalian amphipods still remain poorly understood. Sequencing of mitochondrial genomes is a relatively feasible way for obtaining a set of gene sequences suitable for robust phylogenetic inferences. The architecture of mitochondrial genomes also may provide additional information on the mechanisms of evolution of amphipods in Lake Baikal. RESULTS: Three complete and four nearly complete mitochondrial genomes of Baikalian amphipods were obtained by high-throughput sequencing using the Illumina platform. A phylogenetic inference based on the nucleotide sequences of all mitochondrial protein coding genes revealed the Baikalian species to be a monophyletic group relative to the nearest non-Baikalian species with a completely sequenced mitochondrial genome - Gammarus duebeni. The phylogeny of Baikalian amphipods also suggests that the shallow-water species Eulimnogammarus has likely evolved from a deep-water ancestor, however many other species have to be added to the analysis to test this hypothesis. The gene order in all mitochondrial genomes of studied Baikalian amphipods differs from the pancrustacean ground pattern. Mitochondrial genomes of four species possess 23 tRNA genes, and in three genomes the extra tRNA gene copies have likely undergone remolding. Widely varying lengths of putative control regions and other intergenic spacers are typical for the mitochondrial genomes of Baikalian amphipods. CONCLUSIONS: The mitochondrial genomes of Baikalian amphipods display varying organization suggesting an intense rearrangement process during their evolution. Comparison of complete mitochondrial genomes is a potent approach for studying the amphipod evolution in Lake Baikal.

The complete mitochondrial genome of Baikalospongia intermedia (Lubomirskiidae): description and phylogenetic analysis.

The complete mitochondrial genome of the Lake Baikal sponge Baikalospongia intermedia was sequenced. The circular mitochondrial genome is 28,327 bp in length and includes 14 protein-coding genes, 2 ribosomal RNA genes and 25 transfer RNA genes. Bayesian comparative analysis of molecular evolution rates was found no acceleration of the mtDNA evolution of B. intermedia. This species clustered with other species of the genus Baikalospongia on the Bayesian tree.

Genotyping of Capreolus pygargus fossil DNA from Denisova cave reveals phylogenetic relationships between ancient and modern populations.

BACKGROUND: The extant roe deer (Capreolus Gray, 1821) includes two species: the European roe deer (C. capreolus) and the Siberian roe deer (C. pygargus) that are distinguished by morphological and karyotypical differences. The Siberian roe deer occupies a vast area of Asia and is considerably less studied than the European roe deer. Modern systematics of the Siberian roe deer remain controversial with 4 morphological subspecies. Roe deer fossilized bones are quite abundant in Denisova cave (Altai Mountains, South Siberia), where dozens of both extant and extinct mammalian species from modern Holocene to Middle Pleistocene have been retrieved. METHODOLOGY/PRINCIPAL FINDINGS: We analyzed a 629 bp fragment of the mitochondrial control region from ancient bones of 10 Holocene and four Pleistocene Siberian roe deer from Denisova cave as well as 37 modern specimen belonging to populations from Altai, Tian Shan (Kyrgyzstan), Yakutia, Novosibirsk region and the Russian Far East. Genealogical reconstructions indicated that most Holocene haplotypes were probably ancestral for modern roe deer populations of Western Siberia and Tian Shan. One of the Pleistocene haplotypes was possibly ancestral for modern Yakutian populations, and two extinct Pleistocene haplotypes were close to modern roe deer from Tian Shan and Yakutia. Most modern geographical populations (except for West Siberian Plains) are heterogeneous and there is some tentative evidence for structure. However, we did not find any distinct phylogenetic signal characterizing particular subspecies in either modern or ancient samples. CONCLUSION/SIGNIFICANCE: Analysis of mitochondrial DNA from both ancient and modern samples of Siberian roe deer shed new light on understanding the evolutionary history of roe deer. Our data indicate that during the last 50,000 years multiple replacements of populations of the Siberian roe deer took place in the Altai Mountains correlating with climatic changes. The Siberian roe deer represent a complex and heterogeneous species with high migration rates and without evident subspecies structure. Low genetic diversity of the West Siberian Plain population indicates a recent bottleneck or founder effect.

When environmental changes do not cause geographic separation of fauna: differential responses of Baikalian invertebrates.

BACKGROUND: While the impact of climate fluctuations on the demographic histories of species caused by changes in habitat availability is well studied, populations of species from systems without geographic isolation have received comparatively little attention. Using CO1 mitochondrial sequences, we analysed phylogeographic patterns and demographic histories of populations of five species (four gastropod and one amphipod species) co-occurring in the southwestern shore of Lake Baikal, an area where environmental oscillations have not resulted in geographical isolation of habitats. RESULTS: Species with stronger habitat preferences (gastropods B. turriformis, B. carinata and B. carinatocostata) exhibit rather stable population sizes through their evolutionary history, and their phylogeographic pattern indicates moderate habitat fragmentation. Conversely, species without strong habitat preference (gastropod M. herderiana and amphipod G. fasciatus) exhibit haplotype networks with a very abundant and widespread central haplotype and a big number of singleton haplotypes, while their reconstructed demographic histories show a population expansion starting about 25-50 thousand years ago, a period marked by climate warming and increase in diatom abundance as inferred from bottom-lake sedimentary cores. CONCLUSIONS: In agreement with previous studies, we found that species reacted differently to the same environmental changes. Our results highlight the important role of dispersal ability and degree of ecological specialization in defining a species' response to environmental changes.

Genetic flows in a structured one-dimensional population: simulation and real data on Baikalian polychaetes M. godlewskii.

In this work, we apply population dynamic approach enhanced with simulation of the fate of neutrally evolving "DNA sequences" included into each individual in the computer experiment to the case of the Baikalian endemic polychaetes Manavunkia. These animals inhabit a narrow littoral zone around whole the lake perimeter and are of very limited mobility. Accordingly, the general model was modified by addition of a "geographic barrier" of different isolating power and length of existence. Using this model, we simulated the process of genetic differentiation of groups in this organism taking into account isolation by distance and geographical barriers. Wright's F(st) test was used in order to estimate of genetic flow. Relevant sample sizes were estimated, which would produce most important population parameters of the precision required to describe micro-evolutionary processes in Manayunkia.