De novo sequencing, assembly and functional annotation of Armillaria borealis genome

BACKGROUND

Massive forest decline has been observed almost everywhere as a result of negative anthropogenic and climatic effects, which can interact with pests, fungi and other phytopathogens and aggravate their effects. Climatic changes can weaken trees and make fungi, such as Armillaria more destructive. Armillaria borealis (Marxm. & Korhonen) is a fungus from the Physalacriaceae family (Basidiomycota) widely distributed in Eurasia, including Siberia and the Far East. Species from this genus cause the root white rot disease that weakens and often kills woody plants. However, little is known about ecological behavior and genetics of A. borealis. According to field research data, A. borealis is less pathogenic than A. ostoyae, and its aggressive behavior is quite rare. Mainly A. borealis behaves as a secondary pathogen killing trees already weakened by other factors. However, changing environment might cause unpredictable effects in fungus behavior.

RESULTS

The de novo genome assembly and annotation were performed for the A. borealis species for the first time and presented in this study. The A. borealis genome assembly contained ~ 68 Mbp and was comparable with ~ 60 and ~ 79.5 Mbp for the A. ostoyae and A. mellea genomes, respectively. The N50 for contigs equaled 50,544 bp. Functional annotation analysis revealed 21,969 protein coding genes and provided data for further comparative analysis. Repetitive sequences were also identified. The main focus for further study and comparative analysis will be on the enzymes and regulatory factors associated with pathogenicity.

CONCLUSIONS

Pathogenic fungi such as Armillaria are currently one of the main problems in forest conservation. A comprehensive study of these species and their pathogenicity is of great importance and needs good genomic resources. The assembled genome of A. borealis presented in this study is of sufficiently good quality for further detailed comparative study on the composition of enzymes in other Armillaria species. There is also a fundamental problem with the identification and classification of species of the Armillaria genus, where the study of repetitive sequences in the genomes of basidiomycetes and their comparative analysis will help us identify more accurately taxonomy of these species and reveal their evolutionary relationships.

Antibiotic resistance in Vibrio cholerae El Tor strains isolated during cholera complications in Siberia and the Far East of Russia

Currently, the spread of antimicrobial resistance (AMR) is a global trend and poses a severe threat to public health. The causative agent of cholera, a severe infectious disease with pandemic expansion, becomes more and more resistant to a wider range of drugs with every coming year. The Vibrio cholerae genome is highly flexible and adaptive; the acquisition of the SXT mobile element with a cluster of antibiotic resistance genes on it has marked a new stage in the adaptive evolution of the pathogen. The territory of Siberia and the Russian Far East is free of cholera; however, in the 1970s and 1990s a number of infection importation cases and acute outbreaks associated with the cholera importation were reported. The aim of this study was to describe the phenotypic characteristics and genetic determinants of AMR in V. cholerae strains isolated during epidemic complications in Siberia and the Far East of Russia, as well as to clarify the origin of the strains. The present research comprises analysis of nine V. cholerae El Tor strains isolated from patients and water sources during epidemic complications in Siberia and the Russian Far East in the 1990s. Here, we compared the phenotypic manifestations of antibiotic resistance among strains, harbored the resistance patterns in genomes; we also determined the structure, the type of SXT elements, and the mobilome profile based on the accepted classification. We identified that strains that caused outbreaks in Vladivostok and Yuzhno-Sakhalinsk in 1999 had ICEVchCHN4210 type SXT element with deletion of some loci. The research shows that the integration of the genome, SNP and the mobilome, associated with antibiotic resistance, analyses is necessary to understand the cholera epidemiology, it also helps to establish the origin of strains. The study of resistance determinants features allowed to make a conclusion about the heterogeneity of V. cholerae strains that were isolated during outbreaks in Vladivostok and Yuzhno-Sakhalinsk in 1999.

Composition and Structure of Gut Microbiome in Adolescents with Obesity and Different Breastfeeding Duration

Gut microbiome of adolescents with obesity and different duration of breastfeeding was analyzed by metagenomic analysis of V3-V4 variable domains of the 16S rRNA gene. In subgroup with breastfeeding duration <3 months, intrapopulation structure of gut microbiome by alpha diversity indices was similar in adolescents with obesity and normal body weight. The decrease in phylotype abundance in the structure of communities was associated only with obesity, while dysbiotic state persisted in both lean and overweight adolescents, which confirmed the effect of breastfeeding duration on stability of gut microbiome.

Stepwise large genome assembly approach: a case of Siberian larch (Larix sibirica Ledeb)

Background

De novo assembling of large genomes, such as in conifers (~ 12–30 Gbp), which also consist of ~ 80% of repetitive DNA, is a very complex and computationally intense endeavor. One of the main problems in assembling such genomes lays in computing limitations of nucleotide sequence assembly programs (DNA assemblers). As a rule, modern assemblers are usually designed to assemble genomes with a length not exceeding the length of the human genome (3.24 Gbp). Most assemblers cannot handle the amount of input sequence data required to provide sufficient coverage needed for a high-quality assembly.

Results

An original stepwise method of de novo assembly by parts (sets), which allows to bypass the limitations of modern assemblers associated with a huge amount of data being processed, is presented in this paper. The results of numerical assembling experiments conducted using the model plant Arabidopsis thaliana, Prunus persica (peach) and four most popular assemblers, ABySS, SOAPdenovo, SPAdes, and CLC Assembly Cell, showed the validity and effectiveness of the proposed stepwise assembling method.

Conclusion

Using the new stepwise de novo assembling method presented in the paper, the genome of Siberian larch, Larix sibirica Ledeb. (12.34 Gbp) was completely assembled de novo by the CLC Assembly Cell assembler. It is the first genome assembly for larch species in addition to only five other conifer genomes sequenced and assembled for Picea abies, Picea glauca, Pinus taeda, Pinus lambertiana, and Pseudotsuga menziesii var. menziesii.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2570-y) contains supplementary material, which is available to authorized users.

Comparative genomics of Vibrio cholerae El Tor strains isolated at epidemic complications in Siberia and at the Far East

The territory of Siberia and the Far East of Russia is classified as epidemically safe for cholera; however, in the 1970s and 1990s a number of infection importation cases and acute outbreaks associated with the cholera importation were reported. Here, we analyze genomes of four Vibrio cholerae El Tor strains isolated from humans during epidemic complications (imported cases, an outbreak) in the 1990s. The analyzed strains harbor the classical allele of the cholera toxin subunit B gene (ctxB1); thus, belong to genetically altered variants of the El Tor biotype. Analysis of the genomes revealed their high homology with the V. cholerae N16961 reference strain: 85-93 SNPs were identified in the core genome as compared to the reference. The determined features of SNPs in the CTX prophage made it possible to propose the presence of a new subtype - CTX-2a in two strains; the other two strains carried the prophage of CTX-3 type. Results of phylogenetic analysis based on SNP-typing demonstrated that two strains belonged to the second wave, and two - to the early third wave of cholera dissemination in the world. Phylogenetic reconstruction in combination with epidemiological data permitted to trace the origin of the strains and the way of their importation to the Russian Federation directly or through temporary cholera foci.

De novo assembling and primary analysis of genome and transcriptome of gray whale Eschrichtius robustus

Background

Gray whale, Eschrichtius robustus (E. robustus), is a single member of the family Eschrichtiidae, which is considered to be the most primitive in the class Cetacea. Gray whale is often described as a “living fossil”. It is adapted to extreme marine conditions and has a high life expectancy (77 years). The assembly of a gray whale genome and transcriptome will allow to carry out further studies of whale evolution, longevity, and resistance to extreme environment.

Results

In this work, we report the first de novo assembly and primary analysis of the E. robustus genome and transcriptome based on kidney and liver samples. The presented draft genome assembly is complete by 55% in terms of a total genome length, but only by 24% in terms of the BUSCO complete gene groups, although 10,895 genes were identified. Transcriptome annotation and comparison with other whale species revealed robust expression of DNA repair and hypoxia-response genes, which is expected for whales.

Conclusions

This preliminary study of the gray whale genome and transcriptome provides new data to better understand the whale evolution and the mechanisms of their adaptation to the hypoxic conditions.

Electronic supplementary material

The online version of this article (doi: 10.1186/s12862-017-1103-z) contains supplementary material, which is available to authorized users.

NS2B/NS3 protease: allosteric effect of mutations associated with the pathogenicity of tick-borne encephalitis virus

The sequences of the protease domain of the tick-borne encephalitis (TBE) virus NS3 protein have two amino acid substitutions, 16 R→K and 45 S→F, in the highly pathogenic and poorly pathogenic strains of the virus, respectively. Two models of the NS2B-NS3 protease complex for the highly pathogenic and poorly pathogenic strains of the virus were constructed by homology modeling using the crystal structure of West Nile virus NS2B-NS3 protease as a template; 20 ns molecular dynamic simulations were performed for both models, the trajectories of the dynamic simulations were compared, and the averaged distance between the two models was calculated for each residue. Conformational differences between two models were revealed in the identified pocket. The different conformations of the pocket resulted in different orientations of the NS2B segment located near the catalytic triad. In the model of the highly pathogenic TBE virus the identified pocket had a more open conformation compared to the poorly pathogenic model. We propose that conformational changes in the active protease center, caused by two amino acid substitutions, can influence enzyme functioning and the virulence of the virus.

[GeneBee-NET: An Internet based server for biopolymer structure analysis]

A network server providing biopolymer structure databank retrieval as well as some other biocomputing procedures for Internet users is described. Its basic procedures consist in looking for sequence and 3D homologies (similarities). Found homologies are used for constructing multiple alignment, for predicting RNA and protein secondary structures as well as for constructing phylogenetic trees. Alongside traditional methods of sequence homology retrieval, a "matrix-free" (correlation) method is proposed. A similar procedure is used to locate protein 3D similarities. For novel procedures algorithm ideas and their possible applications are discussed. The service ideology is based on the interaction of server and client programs. The client program (GeneBee for IBM PC) can be used to form queries to the server as well as to manipulate a treatment result. In the absence of the client program the interaction with the server can be in the text mode. The E-mail and WWW addresses for the server are as follows: SERVE/INDY.GENEBEE.MSU.SU and WWW.GENEBEE.MSU.SU.

A novel method of multiple alignment of biopolymer sequences

A novel algorithm for multiple alignment of biological sequences is suggested. At the first step the DotHelix procedure is employed for construction of motifs, i.e. continuous fragments of local similarity of various "thickness" and strength, and then these motifs are concatenated into chains consistent with the order of letters in the sequences. The algorithm is implemented in the MA-Tools program of the GeneBee package. An example illustrating the effectivity of the algorithm is presented.