Genetic Diversity of the Human Adenovirus C Isolated from Hospitalized Children in Russia (2019-2022)

The human adenovirus (HAdV) is a common pathogen in children that can cause acute respiratory virus infection (ARVI). However, the molecular epidemiological and clinical information relating to HAdV among hospitalized children with ARVI is rarely reported in Russia. A 4-year longitudinal (2019-2022) study among hospitalized children (0-17 years old) with ARVI in Novosibirsk, Russia, was conducted to evaluate the epidemiological and molecular characteristics of HAdV. Statistically significant differences in the detection rates of epidemiological and virological data of all positive viral detections of HAdV were analyzed using a two-tailed Chi-square test. The incidence of HAdV and other respiratory viruses such as human influenza A and B viruses, respiratory syncytial virus, coronavirus, parainfluenza virus, metapneumovirus, rhinovirus, bocavirus, and SARS-CoV-2 was investigated among 3190 hospitalized children using real-time polymerase chain reaction. At least one of these respiratory viruses was detected in 74.4% of hospitalized cases, among which HAdV accounted for 4%. A total of 1.3% co-infections with HAdV were also registered. We obtained full-genome sequences of 12 HAdVs, which were isolated in cell cultures. Genetic analysis revealed the circulation of adenovirus of genotypes C1, C2, C5, C89, and 108 among hospitalized children in the period from 2019-2022.

Deficiency of the ribosomal protein uS10 (RPS20) reorganizes human cells translatome according to the abundance, CDS length and GC content of mRNAs

Ribosomal protein uS10, a product of the RPS20 gene, is an essential constituent of the small (40S) subunit of the human ribosome. Disruptive mutations in its gene are associated with a predisposition to hereditary colorectal carcinoma. Here, using HEK293T cells, we show that a deficiency of this protein leads to a decrease in the level of ribosomes (ribosomal shortage). RNA sequencing of the total and polysome-associated mRNA samples reveals hundreds of genes differentially expressed in the transcriptome (t)DEGs and translatome (p)DEGs under conditions of uS10 deficiency. We demonstrate that the (t)DEG and (p)DEG sets partially overlap, determine genes with altered translational efficiency (TE) and identify cellular processes affected by uS10 deficiency-induced ribosomal shortage. We reveal that translated mRNAs of upregulated (p)DEGs and genes with altered TE in uS10-deficient cells are generally more abundant and that their GC contents are significantly lower than those of the respective downregulated sets. We also observed that upregulated (p)DEGs have longer coding sequences. Based on our findings, we propose a combinatorial model describing the process of reorganization of mRNA translation under conditions of ribosomal shortage. Our results reveal rules according to which ribosomal shortage reorganizes the transcriptome and translatome repertoires of actively proliferating cells.

A study of macroinvertebrate communities in Bolshiye Koty Bay of Lake Baikal using DNA metabarcoding

The diversity of macroinvertebrates, the structure of their communities in Bolshiye Koty Bay (Lake Baikal) was studied by a DNA metabarcoding approach using an Illumina MiSeq system. Internal primer mlCOIintF in combination with jgHCO2198 of the Folmer fragment of the COI gene were used for macroinvertebrate metabarcoding. A total of 118009 reads of the COI gene fragment (at least 313 bp in length) were obtained. The correlation of the Spearman coefficient (S = 0.6, p<0.05) with the abundance of macroinvertebrates in the samples before DNA extraction showed that the number of reads can serve as an indirect characteristic of the abundance of a species (operational taxonomic unit, OTU). 115 OTUs belonging to the higher taxa of macroinvertebrates were identified: Porifera, 1; Platyhelminthes, 3; Annelida, 38; Arthropoda, 55; Mollusca, 18. At a high level of resolution (with homology with GenBank reference sequences ≥ 95 %, coverage ≥ 90 %), 46 taxa of macroinvertebrates comprising three communities were registered: one dominated by molluscs (Choanomphalus conf. maacki) and two dominated by chironomids (Orthocladius gregarius Linev., Sergentia baicalensis Tshern.). Communities are characterized by low species diversity according to Shannon (from 0.7 to 1.2 bits), high concentration of dominance according to Simpson (from 0.5 to 0.7) and low evenness according to Pielou (from 0.3 to 0.4). Dominants and subdominants in the communities account for 91 to 96 % of COI gene fragment reads. The spatial distribution of the dominant species identified in the communities is influenced by the geomorphological features of the bottom and the composition of sediments in the area studied. The approach proposed for studying the structure of macroinvertebrate communities based on DNA metabarcoding and next generation sequencing can be recommended for express assessment of the state of aquatic ecosystems in the monitoring.

Description and Genomic Analysis of the First Facultatively Lithoautotrophic, Thermophilic Bacteria of the Genus Thermaerobacter Isolated from Low-temperature Sediments of Lake Baikal

Members of the genus Thermaerobacter belong to the phylum Firmicutes and all isolates characterised to date are strictly aerobic and thermophilic. They were isolated from a mud sample of the Challenger Deep in the Mariana Trench, hydrothermal vents, and silt compost. A novel thermophilic, facultatively lithoautotrophic bacteria of the genus Thermaerobacter, strain PB12/4term (=VKM B-3151T), with a metabolism that is uncharacteristic of the type species, was isolated from low-temperature surface sediments near the Posolsk Bank methane seep, Lake Baikal, Russia. The new strain grows with molecular hydrogen as electron donor, elemental sulfur, and thiosulfate as electron acceptors, and CO2/[Formula: see text] as carbon source. The genome of strain PB12/4term consists of one chromosome with a total length of 2.820.915 bp and the G+C content of the genomic DNA was 72.2%. The phylogenomic reconstruction based on 120 conserved bacterial single-copy proteins revealed that strain PB12/4term belongs to the genus Thermaerobacter within in the class Thermaerobacteria, phylum Firmicutes_E. The strain PB12/4term is closely related to Thermaerobacter subterraneus DSM 13965 (ANI=95.08%, AF=0.91) and Thermaerobacter marianensis DSM 12885 (ANI=84.98%, AF=0.77). Genomic and experimental data confirm the ability of the Thermaerobacter PB12/4term pure culture to facultatively lithotrophic growth, which is provided by the presence of [NiFe]hydrogenase enzymes that are absent in T. marianensis DSM 12885 and T. subterraneus DSM 13965. The data obtained on the physiological and biochemical differences of strain PB12/4term provide a deeper insight into the species diversity and functional activity of the genus Thermaerobacter.

Mutation at the Site of Hydroxylation in the Ribosomal Protein uL15 (RPL27a) Causes Specific Changes in the Repertoire of mRNAs Translated in Mammalian Cells

Ribosomal protein uL15 (RPL27a) carries a specific modification, hydroxylation, at the His39 residue, which neighbors the CCA terminus of the E-site-bound tRNA at the mammalian ribosome. Under hypoxia, the level of hydroxylation of this protein decreases. We transiently transfected HEK293T cells with constructs expressing wild-type uL15 or mutated uL15 (His39Ala) incapable of hydroxylation, and demonstrated that ribosomes containing both proteins are competent in translation. By applying RNA-seq to the total cellular and polysome-associated mRNAs, we identified differentially expressed genes (DEGs) in cells containing exogenous uL15 or its mutant form. Analyzing mRNA features of up- and down-regulated DEGs, we found an increase in the level of more abundant mRNAs and shorter CDSs in cells with uL15 mutant for both translated and total cellular mRNAs. The level of longer and rarer mRNAs, on the contrary, decreased. Our data show how ribosome heterogeneity can change the composition of the translatome and transcriptome, depending on the properties of the translated mRNAs.

Transcriptomic Analysis of CRISPR/Cas9-Mediated PARP1-Knockout Cells under the Influence of Topotecan and TDP1 Inhibitor

Topoisomerase 1 (TOP1) is an enzyme that regulates DNA topology and is essential for replication, recombination, and other processes. The normal TOP1 catalytic cycle involves the formation of a short-lived covalent complex with the 3' end of DNA (TOP1 cleavage complex, TOP1cc), which can be stabilized, resulting in cell death. This fact substantiates the effectiveness of anticancer drugs-TOP1 poisons, such as topotecan, that block the relegation of DNA and fix TOP1cc. Tyrosyl-DNA phosphodiesterase 1 (TDP1) is able to eliminate TOP1cc. Thus, TDP1 interferes with the action of topotecan. Poly(ADP-ribose) polymerase 1 (PARP1) is a key regulator of many processes in the cell, such as maintaining the integrity of the genome, regulation of the cell cycle, cell death, and others. PARP1 also controls the repair of TOP1cc. We performed a transcriptomic analysis of wild type and PARP1 knockout HEK293A cells treated with topotecan and TDP1 inhibitor OL9-119 alone and in combination. The largest number of differentially expressed genes (DEGs, about 4000 both up- and down-regulated genes) was found in knockout cells. Topotecan and OL9-119 treatment elicited significantly fewer DEGs in WT cells and negligible DEGs in PARP1-KO cells. A significant part of the changes caused by PARP1-KO affected the synthesis and processing of proteins. Differences under the action of treatment with TOP1 or TDP1 inhibitors alone were found in the signaling pathways for the development of cancer, DNA repair, and the proteasome. The drug combination resulted in DEGs in the ribosome, proteasome, spliceosome, and oxidative phosphorylation pathways.

The Cytogenetic Map of the Nile Crocodile (Crocodylus niloticus, Crocodylidae, Reptilia) with Fluorescence In Situ Localization of Major Repetitive DNAs

Tandemly arranged and dispersed repetitive DNA sequences are important structural and functional elements that make up a significant portion of vertebrate genomes. Using high throughput, low coverage whole genome sequencing followed by bioinformatics analysis, we have identified seven major tandem repetitive DNAs and two fragments of LTR retrotransposons in the genome of the Nile crocodile (Crocodylus niloticus, 2n = 32). The repeats showed great variability in structure, genomic organization, and chromosomal distribution as revealed by fluorescence in situ hybridization (FISH). We found that centromeric and pericentromeric heterochromatin of C. niloticus is composed of previously described in Crocodylus siamensis CSI-HindIII and CSI-DraI repetitive sequence families, a satellite revealed in Crocodylus porosus, and additionally contains at least three previously unannotated tandem repeats. Both LTR sequences identified here belong to the ERV1 family of endogenous retroviruses. Each pericentromeric region was characterized by a diverse set of repeats, with the exception of chromosome pair 4, in which we found only one type of satellite. Only a few repeats showed non-centromeric signals in addition to their centromeric localization. Mapping of 18S-28S ribosomal RNA genes and telomeric sequences (TTAGGG)n did not demonstrate any co-localization of these sequences with revealed centromeric and pericentromeric heterochromatic blocks.

Flow-Seq Evaluation of Translation Driven by a Set of Natural Escherichia coli 5'-UTR of Variable Length

Flow-seq is a method that combines fluorescently activated cell sorting and next-generation sequencing to deduce a large amount of data about translation efficiency from a single experiment. Here, we constructed a library of fluorescent protein-based reporters preceded by a set of 648 natural 5'-untranslated regions (5'-UTRs) of Escherichia coli genes. Usually, Flow-seq libraries are constructed using uniform-length sequence elements, in contrast to natural situations, where functional elements are of heterogenous lengths. Here, we demonstrated that a 5'-UTR library of variable length could be created and analyzed with Flow-seq. In line with previous Flow-seq experiments with randomized 5'-UTRs, we observed the influence of an RNA secondary structure and Shine-Dalgarno sequences on translation efficiency; however, the variability of these parameters for natural 5'-UTRs in our library was smaller in comparison with randomized libraries. In line with this, we only observed a 30-fold difference in translation efficiency between the best and worst bins sorted with this factor. The results correlated with those obtained with ribosome profiling.

Reorganization of the Landscape of Translated mRNAs in NSUN2-Deficient Cells and Specific Features of NSUN2 Target mRNAs

The RNA cytosine C5 methyltransferase NSUN2 has a variety of RNA substrates and plays an important role in mRNA metabolism. NSUN2 binds to specific sequences enriched in exosomal mRNAs, suggesting its possible involvement in the sorting of mRNAs into exosomes. We applied the photoactivatable.4-thiouridine-enhanced cross-linking and immunoprecipitation assay involving high-throughput RNA sequencing (RNA-seq) to HEK293T cells to determine NSUN2 mRNA targets. NSUN2 cross-linking sites were found in more than one hundred relatively abundant mRNAs with a high GC content and a pronounced secondary structure. Then, utilizing RNA-seq for the total and polysome-associated mRNA from HEK293T cells with and without the knockdown of NSUN2, we identified differentially expressed genes, as well as genes with altered translational efficiency (GATEs). It turned out that the up-regulated GATE mRNAs were much shorter on average than the down-regulated ones, and their GC content was higher; moreover, they contained motifs with C residues located in GC-rich environments. Our findings reveal the specific features of mRNAs that make them potential targets for NSUN2 and expand our understanding of the role of NSUN2 in controlling translation and, possibly, in mRNA sorting into exosomes implemented through the methylation of cytosine residues.

Deconvolution of B cell receptor repertoire in multiple sclerosis patients revealed a delay in tBreg maturation

Background

B lymphocytes play a pivotal regulatory role in the development of the immune response. It was previously shown that deficiency in B regulatory cells (Bregs) or a decrease in their anti-inflammatory activity can lead to immunological dysfunctions. However, the exact mechanisms of Bregs development and functioning are only partially resolved. For instance, only a little is known about the structure of their B cell receptor (BCR) repertoires in autoimmune disorders, including multiple sclerosis (MS), a severe neuroinflammatory disease with a yet unknown etiology. Here, we elucidate specific properties of B regulatory cells in MS.

Methods

We performed a prospective study of the transitional Breg (tBreg) subpopulations with the CD19+CD24highCD38high phenotype from MS patients and healthy donors by (i) measuring their content during two diverging courses of relapsing-remitting MS: benign multiple sclerosis (BMS) and highly active multiple sclerosis (HAMS); (ii) analyzing BCR repertoires of circulating B cells by high-throughput sequencing; and (iii) measuring the percentage of CD27+ cells in tBregs.

Results

The tBregs from HAMS patients carry the heavy chain with a lower amount of hypermutations than tBregs from healthy donors. The percentage of transitional CD24highCD38high B cells is elevated, whereas the frequency of differentiated CD27+ cells in this transitional B cell subset was decreased in the MS patients as compared with healthy donors.

Conclusions

Impaired maturation of regulatory B cells is associated with MS progression.

Deficiency of the Ribosomal Protein uL5 Leads to Significant Rearrangements of the Transcriptional and Translational Landscapes in Mammalian Cells

Protein uL5 (formerly called L11) is an integral component of the large (60S) subunit of the human ribosome, and its deficiency in cells leads to the impaired biogenesis of 60S subunits. Using RNA interference, we reduced the level of uL5 in HEK293T cells by three times, which caused an almost proportional decrease in the content of the fraction corresponding to 80S ribosomes, without a noticeable diminution in the level of polysomes. By RNA sequencing of uL5-deficient and control cell samples, which were those of total mRNA and mRNA from the polysome fraction, we identified hundreds of differentially expressed genes (DEGs) at the transcriptome and translatome levels and revealed dozens of genes with altered translational efficiency (GATEs). Transcriptionally up-regulated DEGs were mainly associated with rRNA processing, pre-mRNA splicing, translation and DNA repair, while down-regulated DEGs were genes of membrane proteins; the type of regulation depended on the GC content in the 3′ untranslated regions of DEG mRNAs. The belonging of GATEs to up-regulated and down-regulated ones was determined by the coding sequence length of their mRNAs. Our findings suggest that the effects observed in uL5-deficient cells result from an insufficiency of translationally active ribosomes caused by a deficiency of 60S subunits.

Gene Analysis, Cloning, and Heterologous Expression of Protease from a Micromycete Aspergillus ochraceus Capable of Activating Protein C of Blood Plasma

Micromycetes are known to secrete numerous enzymes of biotechnological and medical potential. Fibrinolytic protease-activator of protein C (PAPC) of blood plasma from micromycete Aspergillus ochraceus VKM-F4104D was obtained in recombinant form utilising the bacterial expression system. This enzyme, which belongs to the proteinase-K-like proteases, is similar to the proteases encoded in the genomes of Aspergillus fumigatus ATCC MYA-4609, A. oryzae ATCC 42149 and A. flavus 28. Mature PAPC-4104 is 282 amino acids long, preceded by the 101-amino acid propeptide necessary for proper folding and maturation. The recombinant protease was identical to the native enzyme from micromycete in terms of its biological properties, including an ability to hydrolyse substrates of activated protein C (pGlu-Pro-Arg-pNA) and factor Xa (Z-D-Arg-Gly-Arg-pNA) in conjugant reactions with human blood plasma. Therefore, recombinant PAPC-4104 can potentially be used in medicine, veterinary science, diagnostics, and other applications.

Knockdown of the Ribosomal Protein eL38 in HEK293 Cells Changes the Translational Efficiency of Specific Genes

The protein eL38 is one of the smallest proteins of the mammalian ribosome, which is a component of its large (60S) subunit. The haploinsufficiency of eL38 in mice leads to the Tail-short mutant phenotype characterized by defects in the development of the axial skeleton caused by the poor translation of mRNA subsets of Hox genes. Using the ribosome profiling assay applied to HEK293 cells knocked down of eL38, we examined the effects of the lack of eL38 in 60S subunits on gene expression at the level of translation. A four-fold decrease in the cell content of eL38 was shown to result in significant changes in the translational efficiencies of 150 genes. Among the genes, whose expression at the level of translation was enhanced, there were mainly those associated with basic metabolic processes; namely, translation, protein folding, chromosome organization, splicing, and others. The set of genes with reduced translation efficiencies contained those that are mostly involved in the processes related to the regulation of transcription, including the activation of Hox genes. Thus, we demonstrated that eL38 insufficiency significantly affects the expression of certain genes at the translational level. Our findings facilitate understanding the possible causes of some anomalies in eL38-deficient animals.

Extended Evaluation of Viral Diversity in Lake Baikal through Metagenomics

Lake Baikal is a unique oligotrophic freshwater lake with unusually cold conditions and amazing biological diversity. Studies of the lake's viral communities have begun recently, and their full diversity is not elucidated yet. Here, we performed DNA viral metagenomic analysis on integral samples from four different deep-water and shallow stations of the southern and central basins of the lake. There was a strict distinction of viral communities in areas with different environmental conditions. Comparative analysis with other freshwater lakes revealed the highest similarity of Baikal viromes with those of the Asian lakes Soyang and Biwa. Analysis of new data, together with previously published data allowed us to get a deeper insight into the diversity and functional potential of Baikal viruses; however, the true diversity of Baikal viruses in the lake ecosystem remains still unknown. The new metaviromic data will be useful for future studies of viral composition, distribution, and the dynamics associated with global climatic and anthropogenic impacts on this ecosystem.

Knockdown of the mRNA encoding the ribosomal protein eL38 in mammalian cells causes a substantial reorganization of genomic transcription

The ribosomal protein eL38 is a component of the mammalian translation machine. The deletion of the Rpl38 locus in mice results in the Tail-short (Ts) mutant phenotype characterized by a shortened tail and other defects in the axial skeleton development. Here, using the next-generation sequencing of total RNA from HEK293 cells knocked down of eL38 mRNA by transfection with specific siRNAs, we examined the effect of reduced eL38 content on genomic transcription. An approximately 4-fold decrease in the level of eL38 was shown to cause changes in the expression of nearly 1500 genes. Among the down-regulated genes, there were those responsible for p53 activity, Ca2+ metabolism and several signaling processes, as well as genes involved in the organization and functioning of the cytoskeleton. The genes related to rRNA processing and translation, along with many others, including those whose dysregulation is associated with developmental disorders, turned out to be up-regulated. Thus, we demonstrated that the decreased RPL38 expression leads to a significant reorganization of genomic transcription. Our findings suggest a possible link between the balance of eL38 and genes implicated in osteogenesis, thereby contributing to the elucidation of the reasons for the appearance of the above Ts mutant phenotype in animals.

In-depth characterization of ubiquitin turnover in mammalian cells by fluorescence tracking

Despite almost 40 years having passed from the initial discovery of ubiquitin (Ub), fundamental questions related to its intracellular metabolism are still enigmatic. Here we utilized fluorescent tracking for monitoring ubiquitin turnover in mammalian cells, resulting in obtaining qualitatively new data. In the present study we report (1) short Ub half-life estimated as 4 h; (2) for a median of six Ub molecules per substrate as a dynamic equilibrium between Ub ligases and deubiquitinated enzymes (DUBs); (3) loss on average of one Ub molecule per four acts of engagement of polyubiquitinated substrate by the proteasome; (4) direct correlation between incorporation of Ub into the distinct type of chains and Ub half-life; and (5) critical influence of the single lysine residue K27 on the stability of the whole Ub molecule. Concluding, our data provide a comprehensive understanding of ubiquitin-proteasome system dynamics on the previously unreachable state of the art.

Human Peripheral Blood-Derived Endothelial Colony-Forming Cells Are Highly Similar to Mature Vascular Endothelial Cells yet Demonstrate a Transitional Transcriptomic Signature

Endothelial colony-forming cells (ECFC) are currently considered as a promising cell population for the pre-endothelialization or pre-vascularization of tissue-engineered constructs, including small-diameter biodegradable vascular grafts. However, the extent of heterogeneity between ECFC and mature vascular endothelial cells (EC) is unclear. Here, we performed a transcriptome-wide study to compare gene expression profiles of ECFC, human coronary artery endothelial cells (HCAEC), and human umbilical vein endothelial cells (HUVEC). Characterization of the abovementioned cell populations was carried out by immunophenotyping, tube formation assay, and evaluation of proliferation capability while global gene expression profiling was conducted by means of RNA-seq. ECFC were similar to HUVEC in terms of immunophenotype (CD31⁺vWF⁺KDR⁺CD146⁺CD34^-CD133^-CD45^-CD90^-) and tube formation activity yet had expectedly higher proliferative potential. HCAEC and HUVEC were generally similar to ECFC with regards to their global gene expression profile; nevertheless, ECFC overexpressed specific markers of all endothelial lineages (NRP2, NOTCH4, LYVE1), in particular lymphatic EC (LYVE1), and had upregulated extracellular matrix and basement membrane genes (COL1A1, COL1A2, COL4A1, COL4A2). Proteomic profiling for endothelial lineage markers and angiogenic molecules generally confirmed RNA-seq results, indicating ECFC as an intermediate population between HCAEC and HUVEC. Therefore, gene expression profile and behavior of ECFC suggest their potential to be applied for a pre-endothelialization of bioartificial vascular grafts, whereas in terms of endothelial hierarchy they differ from HCAEC and HUVEC, having a transitional phenotype.

mRNA regions where 80S ribosomes pause during translation elongation in vivo interact with protein uS19, a component of the decoding site

In eukaryotic ribosomes, the conserved protein uS19, formerly known as S15, extends with its C-terminal tail to the decoding site. The cross-linking of uS19 to the A site codon has been detected using synthetic mRNAs bearing 4-thiouridine (s⁴U) residues. Here, we showed that the A-site tRNA prevents this cross-linking and that the P site codon does not contact uS19. Next, we focused on determining uS19-mRNA interactions in vivo by applying the photoactivatable-ribonucleoside enhancing cross-linking and immunoprecipitation method to a stable HEK293 cell line producing FLAG-tagged uS19 and grown in a medium containing s⁴U. We found that when translation was stopped by cycloheximide, uS19 was efficiently cross-linked to mRNA regions with a high frequency of Glu, Lys and, more rarely, Arg codons. The results indicate that the complexes, in which the A site codon is not involved in the formation of the mRNA-tRNA duplex, are present among the cycloheximide-arrested 80S complexes, which implies pausing of elongating ribosomes at the above mRNA regions. Thus, our findings demonstrate that the human ribosomal protein uS19 interacts with mRNAs during translation elongation and highlight the regions of mRNAs where ribosome pausing occurs, bringing new structural and functional insights into eukaryotic translation in vivo.

Degenerate consensus sequences in the 3'-untranslated regions of cellular mRNAs as specific motifs potentially involved in the YB-1-mediated packaging of these mRNAs

The multifunctional protein YB-1 has previously been shown to be the only protein of the cytoplasmic extract of HEK293 cells, which is able to specifically interact with imperfect RNA hairpins containing motifs that are often found in exosomal (e) RNAs. In addition, it has been revealed that similar hairpins formed by degenerate consensus sequences corresponding to three eRNA-specific motifs are responsible for the cooperative binding of YB-1 to RNA in vitro. Here, using the photoactivatable ribonucleoside-enhanced cross-linking and immunoprecipitation method applied to HEK293 cells producing FLAG-labeled YB-1, we identified mRNAs cross-linked to YB-1 in vivo and then carried out a search for the aforementioned sequences in the regions of the YB-1 cross-linking sites. It turned out that many of the mRNAs found cross-linked to YB-1 encode proteins associated with various regulatory processes, including responses to stress. More than half of all cross-linked mRNAs contained degenerate consensus sequences, which were preferably located in 3'-untranslated regions (UTRs), where most of the YB-1 cross-linking sites appeared, although not close to these sequences. Furthermore, YB-1 was mainly cross-linked to those mRNAs with degenerate consensus sequences, which could be classified as packaged because their translation levels were low compared to cellular levels. This suggests that the cooperative binding of YB-1 to mRNAs through the above sequences probably triggers the well-known multimerization of YB-l, leading to the packaging of these mRNAs. Thus, our findings indicate a previously unknown link between the degenerate consensus sequences present in the 3'-UTRs of many cytoplasmic mRNAs and YB-1-mediated translational silencing.

Human ribosomal protein eS1 is engaged in cellular events related to processing and functioning of U11 snRNA

Ribosomal proteins are involved in many cellular processes through interactions with various RNAs. Here, applying the photoactivatable-ribonucleoside-enhanced cross-linking and immunoprecipitation approach to HEK293 cells overproducing ribosomal protein (rp) eS1, we determined the products of RNU5A-1 and RNU11 genes encoding U5 and U11 snRNAs as the RNA partners of ribosome-unbound rp eS1. U11 pre-snRNA-associated rp eS1 was revealed in the cytoplasm and nucleus where rp eS1-bound U11/U12 di-snRNP was also found. Utilizing recombinant rp eS1 and 4-thiouridine-containing U11 snRNA transcript, we identified an N-terminal peptide contacting the U-rich sequence in the Sm site-containing RNA region. We also showed that the rp eS1 binding site on U11 snRNA is located in the cleft between stem-loops I and III and that its structure mimics the respective site on the 18S rRNA. It was found that cell depletion of rp eS1 leads to a decrease in the splicing efficiency of minor introns and to an increase in the level of U11 pre-snRNA with the unprocessed 3' terminus. Our findings demonstrate the engagement of human rp eS1 in events related to the U11 snRNA processing and to minor-class splicing. Contacts of rp eS1 with U5 snRNA in the minor pre-catalytic spliceosome are discussed.

Recognition but no repair of abasic site in single-stranded DNA by human ribosomal uS3 protein residing within intact 40S subunit

Isolated human ribosomal protein uS3 has extra-ribosomal functions including those related to base excision DNA repair, e.g. AP lyase activity that nicks double-stranded (ds) DNA 3΄ to the abasic (AP) site. However, the ability of uS3 residing within ribosome to recognize and cleave damaged DNA has never been addressed. Here, we compare interactions of single-stranded (ss) DNA and dsDNA bearing AP site with human ribosome-bound uS3 and with the isolated protein, whose interactions with ssDNA were not yet studied. The AP lyase activity of free uS3 was much higher with ssDNA than with dsDNA, whereas ribosome-bound uS3 was completely deprived of this activity. Nevertheless, an exposed peptide of ribosome-bound uS3 located far away from the putative catalytic center previously suggested for isolated uS3 cross-linked to full-length uncleaved ssDNA, but not to dsDNA. In contrast, free uS3 cross-linked mainly to the 5΄-part of the damaged DNA strand after its cleavage at the AP site. ChIP-seq analysis showed preferential uS3 binding to nucleolus-associated chromatin domains. We conclude that free and ribosome-bound uS3 proteins interact with AP sites differently, exhibiting their non-translational functions in DNA repair in and around the nucleolus and in regulation of DNA damage response in looped DNA structures, respectively.

Cytosolic YB-1 and NSUN2 are the only proteins recognizing specific motifs present in mRNAs enriched in exosomes

Exosomes, membranous vesicles secreted by various cells, are involved in intercellular communication and carry vast repertoires of RNAs and proteins. Processes mediating RNA sorting into exosomes are currently poorly understood. Using bioinformatics approaches, three structural motifs ACCAGCCU, CAGUGAGC and UAAUCCCA have been discovered as enriched in exosomal mRNAs and long noncoding RNAs. Here, utilizing short RNA hairpins, each containing one of the motifs, in a pull-down assay of cytosolic extract of human embryonic kidney 293 (HEK293) cells, we prove that multifunctional RNA-binding protein YB-1 specifically interacts with all three motifs, whereas methyltransferase NSUN2 recognizes only the motif CAGUGAGC. RNA hairpins other than those mentioned above pull out neither YB-1 nor NSUN2. Both these proteins are found in exosomes secreted by HEK293 cells. YB-1 for all that is detected as a form having a slightly higher electrophoretic mobility than that of YB-1 associated with the above RNA hairpins, assuming changes in posttranslational modifications of the protein during its transfer from cytoplasm into exosomes. Next generation sequencing of total exosomal RNA (eRNA) reveals a large representative set of RNA species, including mRNAs containing the above-mentioned motifs. The degree of enrichment in exosomes with this kind of mRNAs strongly depends on the locations of eRNA-specific motifs within the mRNA sequences. Altogether, our findings point to YB-1 and NSUN2 as possible mediators of the process of transfer of specific mRNAs into exosomes, allowing us to speculate on an involvement of these proteins in the mRNA sorting via the recognition of the above motifs.

Comparative Metagenomic Profiling of Symbiotic Bacterial Communities Associated with Ixodes persulcatus, Ixodes pavlovskyi and Dermacentor reticulatus Ticks

Ixodes persulcatus, Ixodes pavlovskyi, and Dermacentor reticulatus ticks inhabiting Western Siberia are responsible for the transmission of a number of etiological agents that cause human and animal tick-borne diseases. Because these ticks are abundant in the suburbs of large cities, agricultural areas, and popular tourist sites and frequently attack people and livestock, data regarding the microbiomes of these organisms are required. Using metagenomic 16S profiling, we evaluate bacterial communities associated with I. persulcatus, I. pavlovskyi, and D. reticulatus ticks collected from the Novosibirsk region of Russia. A total of 1214 ticks were used for this study. DNA extracted from the ticks was pooled according to tick species and sex. Sequencing of the V3-V5 domains of 16S rRNA genes was performed using the Illumina Miseq platform. The following bacterial genera were prevalent in the examined communities: Acinetobacter (all three tick species), Rickettsia (I. persulcatus and D. reticulatus) and Francisella (D. reticulatus). B. burgdorferi sensu lato and B. miyamotoi sequences were detected in I. persulcatus and I. pavlovskyi but not in D. reticulatus ticks. The pooled samples of all tick species studied contained bacteria from the Anaplasmataceae family, although their occurrence was low. DNA from A. phagocytophilum and Candidatus Neoehrlichia mikurensis was first observed in I. pavlovskyi ticks. Significant inter-species differences in the number of bacterial taxa as well as intra-species diversity related to tick sex were observed. The bacterial communities associated with the I. pavlovskyi ticks displayed a higher biodiversity compared with those of the I. persulcatus and D. reticulatus ticks. Bacterial community structure was also diverse across the studied tick species, as shown by permutational analysis of variance using the Bray-Curtis dissimilarity metric (p = 0.002). Between-sex variation was confirmed by PERMANOVA testing in I. persulcatus (p = 0.042) and I. pavlovskyi (p = 0.042) ticks. Our study indicated that 16S metagenomic profiling could be used for rapid assessment of the occurrence of medically important bacteria in tick populations inhabiting different natural biotopes and therefore the epidemic danger of studied foci.

[Mitochondrial DNA variation in Olga Bay larch (Larix olgensis A. Henry) from Primorsky Krai of Russia]

Two mitochondrial DNA fragments, nad4(3c-4r) and nad5(1-2r), were sequenced in 58 larch accessions from the range of Larix olgensis A. Henry. Combinations of the nad4 polymorphic sites formed four haplotypes, two of which (H3 and H4) were unique and two (H1, H2) were common. Haplotype H1 was found only in pure L. olgensis from the vicinity of Olga Bay and in a number of accessions from the southern part of the range. Haplotype H2 was detected in the other samples from the range of Olga Bay larch, as well as in hybrid forms. Similarly to the nad4(3c-4r) fragment, the mtDNA fragment UBC460 was able to differentiate larch populations from the range of L. olgensis examined.

Heavy-light chain interrelations of MS-associated immunoglobulins probed by deep sequencing and rational variation

The mechanisms triggering most of autoimmune diseases are still obscure. Autoreactive B cells play a crucial role in the development of such pathologies and, in particular, production of autoantibodies of different specificities. The combination of deep-sequencing technology with functional studies of antibodies selected from highly representative immunoglobulin combinatorial libraries may provide unique information on specific features in the repertoires of autoreactive B cells. Here, we have analyzed cross-combinations of the variable regions of human immunoglobulins against the myelin basic protein (MBP) previously selected from a multiple sclerosis (MS)-related scFv phage-display library. On the other hand, we have performed deep sequencing of the sublibraries of scFvs against MBP, Epstein-Barr virus (EBV) latent membrane protein 1 (LMP1), and myelin oligodendrocyte glycoprotein (MOG). Bioinformatics analysis of sequencing data and surface plasmon resonance (SPR) studies have shown that it is the variable fragments of antibody heavy chains that mainly determine both the affinity of antibodies to the parent autoantigen and their cross-reactivity. It is suggested that LMP1-cross-reactive anti-myelin autoantibodies contain heavy chains encoded by certain germline gene segments, which may be a hallmark of the EBV-specific B cell subpopulation involved in MS triggering.

[Phenotypic variation of spermatogenesis and a search for associations with genetic polymorphism in 13 inbred mouse strains]

Adult mice of the BALB/cLac, PT, CBA/Lac, DD/He, A/He, SWR, NZB, GR, DBA/2J, CC57Br, C57 B1/6J, A/Sn, and YT inbred strains were tested for the count, motility, and morphology of sperms from the caudal region of the epididymis. The protein-coding regions of the cytochrome P450 aromatase (CYP19a1), estrogen receptor 2 (ESR2), steroidogenic factor 1 (Nr5a1), and sex-determining (Sry) gene were sequenced. A substantial genetic heterogeneity for the genes was observed, as well as a phenotypic variation in spermatogenetic parameters, but the variation was rather discordant. The specifics of the interstrain variation in spermatogenetic parameters indicated that a physiological compensatory mechanism increases certain spermatogenetic parameters when other ones are low to maintain male fertility at a level sufficient for successful reproduction. For instance, a high sperm production compensated for a low sperm motility in DD/He males. In the issue of the protein-coding regions sequencing of the analyzed genes, 16 various mutations were observed. The decreases in proportion of motile sperms and in their velocity were attributed to mutations (I63T and W133L) of the Sry gene in the DD/He strain.

[Mutation analysis of the phenylalanine hydroxylase gene of phenylketonuria patients of Kemerovskaya Oblast' and Saha Republic]

Phenylketonuria (PKU) associated mutations in phenylalanine hydroxylase (PAH) gene were identified by direct DNA sequencing in 46 PKU patients and members of their families from Kemerovskaya Region and Saha Republic. Mutations found included both widespread known mutations (R158Q, R252W, R261Q, P281L, IVS10-11G>A, R408W, IVS12+1G>A) and several rare mutations (IVS2+5G>A, R155H, Y168H, W187R, E221_D222>Efs, A342T, Y386C, IVS11+1G>C). We observed the increase in diversity of PKU-associated alleles in the populations studied, probably due to their complex mixed ethnic structure.