Search

Abstract Background As genomes of many eukaryotic species, especially plants, are large and complex, their de novo sequencing and assembly is still a difficult task despite progress in sequencing technologies. An alternative to genome assembly is the assembly of transcriptome, the set of RNA products of the expressed genes. While a bunch of de novo transcriptome assemblers exists, the challenges of transcriptomes (the existence of isoforms, the uneven expression levels across genes) complicates the generation of high-quality assemblies suitable for downstream analyses. Results We developed Trans2express – a web-based tool and a pipeline of de novo hybrid transcriptome assembly and postprocessing based on rnaSPAdes with a set of subsequent filtrations. The pipeline was tested on Arabidopsis thaliana cDNA sequencing data obtained using Illumina and Oxford Nanopore Technologies platforms and three non-model plant species. The comparison of structural characteristics of the transcriptome assembly with reference Arabidopsis genome revealed the high quality of assembled transcriptome with 86.1% of Arabidopsis expressed genes assembled as a single contig. We tested the applicability of the transcriptome assembly for gene expression analysis. For both Arabidopsis and non-model species the results showed high congruence of gene expression levels and sets of differentially expressed genes between analyses based on genome and based on the transcriptome assembly. Conclusions We present Trans2express – a protocol for de novo hybrid transcriptome assembly aimed at recovering of a single transcript per gene. We expect this protocol to promote the characterization of transcriptomes and gene expression analysis in non-model plants and web-based tool to be of use to a wide range of plant biologists.

Abstract Motivation Alignment of reads to a reference genome sequence is one of the key steps in the analysis of human whole-genome sequencing data obtained through Next-generation sequencing (NGS) technologies. The quality of the subsequent steps of the analysis, such as the results of clinical interpretation of genetic variants or the results of a genome-wide association study, depends on the correct identification of the position of the read as a result of its alignment. The amount of human NGS whole-genome sequencing data is constantly growing. There are a number of human genome sequencing projects worldwide that have resulted in the creation of large-scale databases of genetic variants of sequenced human genomes. Such information about known genetic variants can be used to improve the quality of alignment at the read alignment stage when analysing sequencing data obtained for a new individual, for example, by creating a genomic graph. While existing methods for aligning reads to a linear reference genome have high alignment speed, methods for aligning reads to a genomic graph have greater accuracy in variable regions of the genome. The development of a read alignment method that takes into account known genetic variants in the linear reference sequence index allows combining the advantages of both sets of methods. Results In this paper, we present the minimap2_index_modifier tool, which enables the construction of a modified index of a reference genome using known single nucleotide variants and insertions/deletions (indels) specific to a given human population. The use of the modified minimap2 index improves variant calling quality without modifying the bioinformatics pipeline and without significant additional computational overhead. Using the PrecisionFDA Truth Challenge V2 benchmark data (for HG002 short-read data aligned to the GRCh38 linear reference (GCA_000001405.15) with parameters k = 27 and w = 14) it was demonstrated that the number of false negative genetic variants decreased by more than 9500, and the number of false positives decreased by more than 7000 when modifying the index with genetic variants from the Human Pangenome Reference Consortium.

Abstract Background Capsella bursa-pastoris, a cosmopolitan weed of hybrid origin, is an emerging model object for the study of early consequences of polyploidy, being a fast growing annual and a close relative of Arabidopsis thaliana. The development of this model is hampered by the absence of a reference genome sequence. Results We present here a subgenome-resolved chromosome-scale assembly and a genetic map of the genome of Capsella bursa-pastoris. It shows that the subgenomes are mostly colinear, with no massive deletions, insertions, or rearrangements in any of them. A subgenome-aware annotation reveals the lack of genome dominance—both subgenomes carry similar number of genes. While most chromosomes can be unambiguously recognized as derived from either paternal or maternal parent, we also found homeologous exchange between two chromosomes. It led to an emergence of two hybrid chromosomes; this event is shared between distant populations of C. bursa-pastoris. The whole-genome analysis of 119 samples belonging to C. bursa-pastoris and its parental species C. grandiflora/rubella and C. orientalis reveals introgression from C. orientalis but not from C. grandiflora/rubella. Conclusions C. bursa-pastoris does not show genome dominance. In the earliest stages of evolution of this species, a homeologous exchange occurred; its presence in all present-day populations of C. bursa-pastoris indicates on a single origin of this species. The evidence coming from whole-genome analysis challenges the current view that C. grandiflora/rubella was a direct progenitor of C. bursa-pastoris; we hypothesize that it was an extinct (or undiscovered) species sister to C. grandiflora/rubella.

The flyways of many different wild waterfowl pass through the Caspian Sea region. The western coast of the middle Caspian Sea is an area with many wetlands, where wintering grounds with large concentrations of birds are located. It is known that wild waterfowl are a natural reservoir of the influenza A virus. In the mid-2000s, in the north of this region, the mass deaths of swans, gulls, and pelicans from high pathogenicity avian influenza virus (HPAIV) were noted. At present, there is still little known about the presence of avian influenza virus (AIVs) and different avian paramyxoviruses (APMVs) in the region’s waterfowl bird populations. Here, we report the results of monitoring these viruses in the wild waterfowl of the western coast of the middle Caspian Sea from 2017 to 2020. Samples from 1438 individuals of 26 bird species of 7 orders were collected, from which 21 strains of AIV were isolated, amounting to a 1.46% isolation rate of the total number of samples analyzed (none of these birds exhibited external signs of disease). The following subtypes were determined and whole-genome nucleotide sequences of the isolated strains were obtained: H1N1 (n = 2), H3N8 (n = 8), H4N6 (n = 2), H7N3 (n = 2), H8N4 (n = 1), H10N5 (n = 1), and H12N5 (n = 1). No high pathogenicity influenza virus H5 subtype was detected. Phylogenetic analysis of AIV genomes did not reveal any specific pattern for viruses in the Caspian Sea region, showing that all segments belong to the Eurasian clades of classic avian-like influenza viruses. We also did not find the amino acid substitutions in the polymerase complex (PA, PB1, and PB2) that are critical for the increase in virulence or adaptation to mammals. In total, 23 hemagglutinating viruses not related to influenza A virus were also isolated, of which 15 belonged to avian paramyxoviruses. We were able to sequence 12 avian paramyxoviruses of three species, as follows: Newcastle disease virus (n = 4); Avian paramyxovirus 4 (n = 5); and Avian paramyxovirus 6 (n = 3). In the Russian Federation, the Newcastle disease virus of the VII.1.1 sub-genotype was first isolated from a wild bird (common pheasant) in the Caspian Sea region. The five avian paramyxovirus 4 isolates obtained belonged to the common clade in Genotype I, whereas phylogenetic analysis of three isolates of Avian paramyxovirus 6 showed that two isolates, isolated in 2017, belonged to Genotype I and that an isolate identified in 2020 belonged to Genotype II. The continued regular monitoring of AIVs and APMVs, the obtaining of data on the biological properties of isolated strains, and the accumulation of information on virus host species will allow for the adequate planning of epidemiological measures, suggest the most likely routes of spread of the virus, and assist in the prediction of the introduction of the viruses in the western coastal region of the middle Caspian Sea.

Neural crest cells are highly multipotent stem cells, but it remains unclear how their fate restriction to specific fates occurs. Here, the authors show in zebrafish that broad multipotency is retained even after migration, suggesting that fate restriction occurs directly, but dynamically.

The wetlands of southwestern Siberia (SWS) are a crossroads of bird migration routes, bringing avian influenza (AIV) strains that were previously isolated in different regions of the continent to Siberia. It is known that Anseriformes that breed in SWS migrate for the winter to central Hindustan or further west, while their migration routes to southeast Asia (SEA) remain unconfirmed. Here, we mapped the molting sites of the migrating Common Teals (Anas crecca) via analyzing stable hydrogen isotope content in feathers of hunters’ prey and supplemented the analysis with the genetic structure of viruses isolated from teals in the same region. Post-breeding molt of autumn teals most likely occurred within the study region, whereas probable pre-breeding molting grounds of spring teals were in the south of Hindustan. This link was supported by viral phylogenetic analysis, which showed a close relationship between SWS isolates and viruses from south and southeast Asia. Most viral segments have the highest genetic similarity and the closest phylogenetic relationships with viruses from teal wintering areas in southeast Asian countries, including India and Korea. We assume that the winter molt of SWS breeding teals on the Hindustan coast suggests contacts with the local avifauna, including species migrating along the coast to SEA. Perhaps this is one of the vectors of AIV transmission within Eurasia.

The article deals with decommunization in Ukraine from 2014 to 2020, examining the motives and actions of the main actors of this policy and societal responses and recounting its most significant developments and outcomes. The decommunization in Ukraine in 2014–2020 was a successful, though belated, effort initiated at a time when Communist symbols had lost their original ideological sense as well as social, cultural, and political significance. The promoters of decommunization reanimated it and managed to contextualize it under actual threat from Russia to justify their policy, which they used to promote a nationalist narrative of memory. [ABSTRACT FROM AUTHOR]

Grass pollen is one of the leading causes of pollinosis, affecting 10–30% of the world’s population. The allergenicity of pollen from different Poaceae species is not the same and is estimated from moderate to high. Aerobiological monitoring is a standard method that allows one to track and predict the dynamics of allergen concentration in the air. Poaceae is a stenopalynous family, and thus grass pollen can usually be identified only at the family level with optical microscopy. Molecular methods, in particular the DNA barcoding technique, can be used to conduct a more accurate analysis of aerobiological samples containing the DNA of various plant species. This study aimed to test the possibility of using the ITS1 and ITS2 nuclear loci for determining the presence of grass pollen from air samples via metabarcoding and to compare the analysis results with the results of phenological observations. Based on the high-throughput sequencing data, we analyzed the changes in the composition of aerobiological samples taken in the Moscow and Ryazan regions for three years during the period of active flowering of grasses. Ten genera of the Poaceae family were detected in airborne pollen samples. The representation for most of them for ITS1 and ITS2 barcodes was similar. At the same time, in some samples, the presence of specific genera was characterized by only one sequence: either ITS1 or ITS2. Based on the analysis of the abundance of both barcode reads in the samples, the following order could describe the change with time in the dominant species in the air: Poa, Alopecurus, and Arrhenatherum in early mid-June, Lolium, Bromus, Dactylis, and Briza in mid-late June, Phleum, Elymus in late June to early July, and Calamagrostis in early mid-July. In most samples, the number of taxa found via metabarcoding analysis was higher compared to that in the phenological observations. The semi-quantitative analysis of high-throughput sequencing data well reflects the abundance of only major grass species at the flowering stage.

Abstract Background Salivary cell secretion (SCS) plays a critical role in blood feeding by medicinal leeches, making them of use for certain medical purposes even today. Results We annotated the Hirudo medicinalis genome and performed RNA-seq on salivary cells isolated from three closely related leech species, H. medicinalis, Hirudo orientalis, and Hirudo verbana. Differential expression analysis verified by proteomics identified salivary cell-specific gene expression, many of which encode previously unknown salivary components. However, the genes encoding known anticoagulants have been found to be expressed not only in salivary cells. The function-related analysis of the unique salivary cell genes enabled an update of the concept of interactions between salivary proteins and components of haemostasis. Conclusions Here we report a genome draft of Hirudo medicinalis and describe identification of novel salivary proteins and new homologs of genes encoding known anticoagulants in transcriptomes of three medicinal leech species. Our data provide new insights in genetics of blood-feeding lifestyle in leeches.

The approach to managing patients after femoropopliteal bypass surgery among vascular surgeons, even of one country, is sometimes not the same. The main difficulty is the problem of the late postoperative period patency of femoropopliteal bypasses (FPB). Aim. Analysis of the risk factors for occlusion of the FPB and the effects of postoperative diagnostic ultrasound on the primary, primary-assisted patency and secondary patency of such bypasses. Materials and methods. The Pubmed database was analyzed and the studies were selected from 2007 till 2017. The effect of preoperative risk factors on the FPB’s primary patency was evaluated in patients of the department of vascular surgery at the Orel regional hospital in comparison with the of foreign researchers’ data. The second stage was the evaluation of the effect of ultrasound diagnostic on the primary, primary-assisted and secondary patency of such bypasses. Results. The greatest impact on the primary patency of FPB was exerted by the presence of CLI, arterial hypertension and the male gender, and a slightly fewer by smoking and angina pectoris (χ2 cf. = -0.983-0.985*). The effect of diabetes was 2 times weaker (χ2 = -0.547*). Ultrasound control, in turn, had the greatest impact on the frequency of primary-assisted patency compared with the group without ultrasound control (81.3% vs. 76.1%). The frequencies of primary patency (67.4% vs. 64.9%) and secondary patency (83.2% vs. 82.8%) turned out to be almost comparable between the two groups of patients. Conclusions. Thus, the male sex, arterial hypertension and CLI are the most significant pre-operative risk factors affecting the primary patency of FPB. Regular diagnostic ultrasound in the postoperative period, especially in the period of 4-6 weeks, 3 months, 6 months and 1 year, allows to improve long-term results of such bypasses’ patency.

The vast diversity of Orchidaceae together with sophisticated adaptations to pollinators and other unique features make this family an attractive model for evolutionary and functional studies. The sequenced genome of Phalaenopsis equestris facilitates Orchidaceae research. Here, we present an RNA-seq-based transcriptome map of P. equestris that covers 19 organs of the plant, including leaves, roots, floral organs and the shoot apical meristem. We demonstrated the high quality of the data and showed the similarity of the P. equestris transcriptome map with the gene expression atlases of other plants. The transcriptome map can be easily accessed through our database Transcriptome Variation Analysis (TraVA) for visualizing gene expression profiles. As an example of the application, we analyzed the expression of Phalaenopsis “orphan” genes–those that do not have recognizable similarity with the genes of other plants. We found that approximately half of these genes were not expressed; the ones that were expressed were predominantly expressed in reproductive structures.

Naturally occurring mutants whose phenotype recapitulates the changes that distinguish closely related species are of special interest from the evolutionary point of view. They can give a key about the genetic control of the changes that led to speciation. In this study, we described lepidium-like (lel), a naturally occurring variety of an allotetraploid species Capsella bursa-pastoris that is characterized by the typical loss of all four petals. In some cases, one or two basal flowers in the raceme had one or two small petals. The number and structure of other floral organs are not affected. Our study of flower development in the mutant showed that once initiated, petals either cease further development and cannot be traced in anthetic flowers or sometimes develop to various degrees. lel plants showed an earlier beginning of floral organ initiation and delayed petal initiation compared to the wild-type plants. lel phenotype has a wide geographical distribution, being found at the northern extremity of the species range as well as in the central part. The genetic analysis of inheritance demonstrated that lel phenotype is controlled by two independent loci. While the flower in the family Cruciferae generally has a very stable structure (i.e., four sepals, four petals, six stamens, and two carpels), several deviations from this ground plan are known, in particular in the genus Lepidium, C. bursa-pastoris is an emerging model for the study of polyploidy (which is also very widespread in Cruciferae); the identification and characterization of the apetalous mutant lays a foundation for further research of morphological evolution in polyploids.

Abstract Background All living organisms experience physiological changes regulated by endogenous circadian rhythms. The main factor controlling the circadian clock is the duration of daylight. The aim of this research was to identify the impact of various lighting conditions on physiological parameters and gut microbiota composition in rats. 3 groups of outbred rats were subjected to normal light-dark cycles, darkness and constant lighting. Results After 1 and 3 months we studied urinary catecholamine levels in rats; indicators of lipid peroxidation and antioxidant activity in the blood; protein levels of BMAL1, CLOCK and THRA in the hypothalamus; composition and functional activity of the gut microbiota. Subjecting the rats to conditions promoting desynchronosis for 3 months caused disruptions in homeostasis. Conclusions Changing the lighting conditions led to changes in almost all the physiological parameters that we studied. Catecholamines can be regarded as a synchronization super system of split-level circadian oscillators. We established a correlation between hypothalamic levels of Bmal1 and urinary catecholamine concentrations. The magnitude of changes in the GM taxonomic composition was different for LL/LD and DD/LD but the direction of these changes was similar. As for the predicted functional properties of the GM which characterize its metabolic activity, they didn’t change as dramatically as the taxonomic composition. All differences may be viewed as a compensatory reaction to new environmental conditions and the organism has adapted to those conditions.

Abstract Background Transcriptome map is a powerful tool for a variety of biological studies; transcriptome maps that include different organs, tissues, cells and stages of development are currently available for at least 30 plants. Some of them include samples treated by environmental or biotic stresses. However, most studies explore only limited set of organs and developmental stages (leaves or seedlings). In order to provide broader view of organ-specific strategies of cold stress response we studied expression changes that follow exposure to cold (+ 4 °C) in different aerial parts of plant: cotyledons, hypocotyl, leaves, young flowers, mature flowers and seeds using RNA-seq. Results The results on differential expression in leaves are congruent with current knowledge on stress response pathways, in particular, the role of CBF genes. In other organs, both essence and dynamics of gene expression changes are different. We show the involvement of genes that are confined to narrow expression patterns in non-stress conditions into stress response. In particular, the genes that control cell wall modification in pollen, are activated in leaves. In seeds, predominant pattern is the change of lipid metabolism. Conclusions Stress response is highly organ-specific; different pathways are involved in this process in each type of organs. The results were integrated with previously published transcriptome map of Arabidopsis thaliana and used for an update of a public database TraVa: http://travadb.org/browse/Species=AthStress.

Common buckwheat (Fagopyrum esculentum) is an important non-cereal grain crop and a prospective component of functional food. Despite this, the genomic resources for this species and for the whole family Polygonaceae, to which it belongs, are scarce. Here, we report the assembly of the buckwheat genome using long-read technology and a high-resolution expression atlas including 46 organs and developmental stages. We found that the buckwheat genome has an extremely high content of transposable elements, including several classes of recently (0.5–1 Mya) multiplied TEs (“transposon burst”) and gradually accumulated TEs. The difference in TE content is a major factor contributing to the three-fold increase in the genome size of F. esculentum compared with its sister species F. tataricum. Moreover, we detected the differences in TE content between the wild ancestral subspecies F. esculentum ssp. ancestrale and buckwheat cultivars, suggesting that TE activity accompanied buckwheat domestication. Expression profiling allowed us to test a hypothesis about the genetic control of petaloidy of tepals in buckwheat. We showed that it is not mediated by B-class gene activity, in contrast to the prediction from the ABC model. Based on a survey of expression profiles and phylogenetic analysis, we identified the MYB family transcription factor gene tr_18111 as a potential candidate for the determination of conical cells in buckwheat petaloid tepals. The information on expression patterns has been integrated into the publicly available database TraVA: http://travadb.org/browse/Species=Fesc/. The improved genome assembly and transcriptomic resources will enable research on buckwheat, including practical applications.

As permanent residents of the normal gut microbiota, bifidobacteria have evolved to adapt to the host’s immune response whose priority is to eliminate pathogenic agents. The mechanisms that ensure the survival of commensals during inflammation and maintain the stability of the core component of the normal gut microbiota in such conditions remain poorly understood. We propose a new in vitro approach to study the mechanisms of resistance to immune response factors based on high-throughput sequencing followed by transcriptome analysis. This approach allowed us to detect differentially expressed genes associated with inflammation. In this study, we demonstrated that the presence of the pro-inflammatory cytokines IL-6 and TNFα to the growth medium of the B. longum subsp. longum GT15 strain changes the latter’s growth rate insignificantly while affecting the expression of certain genes. We identified these genes and performed a COG and a KEGG pathway enrichment analysis. Using phylogenetic profiling we predicted the operons of genes whose expression was triggered by the cytokines TNFα and IL-6 in vitro. By mapping the transcription start points, we experimentally validated the predicted operons. Thus, in this study, we predicted the genes involved in a putative signaling pathway underlying the mechanisms of resistance to inflammatory factors in bifidobacteria. Since bifidobacteria are a major component of the human intestinal microbiota exhibiting pronounced anti-inflammatory properties, this study is of great practical and scientific relevance.

Grass pollen is one of the major causes of allergy. Aerobiological monitoring is a necessary element of the complex of anti-allergic measures, but the similar pollen morphology of Poaceae species makes it challenging to discriminate species in airborne pollen mixes, which impairs the quality of aerobiological monitoring. One of the solutions to this problem is the metabarcoding approach employing DNA barcodes for taxonomical identification of species in a mix by high-throughput sequencing of the pollen DNA. A diverse set of 14 grass species of different genera were selected to create a local reference database of nuclear ITS1, ITS2, 5′-ETS, and plastome trnL-F DNA barcodes. Sequences for the database were Sanger sequenced from live field and herbarium specimens and collected from GenBank. New Poaceae-specific primers for 5′-ETS were designed and tested to obtain a 5′-ETS region less than 600 bp long, suitable for high-throughput sequencing. The DNA extraction method for single-species pollen samples and mixes was optimized to increase the yield for amplification and sequencing of pollen DNA. Barcode sequences were analyzed and compared by the barcoding gap and intra- and interspecific distances. Their capability to correctly identify grass pollen was tested on artificial pollen mixes of various complexity. Metabarcoding analysis of the artificial pollen mixes showed that nuclear DNA barcodes ITS1, ITS2, and 5′-ETS proved to be more efficient than the plastome barcode in both amplification from pollen DNA and identification of grass species. Although the metabarcoding results were qualitatively congruent with the actual composition of the pollen mixes in most cases, the quantitative results based on read-counts did not match the actual ratio of pollen grains in the mixes.

Abstract 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.