Sequencing the transcriptome of single cells: a path to discovering cellular biodiversity / Sekwencjonowanie transkryptomu pojedynczych komórek: droga do odkrywania bioróznorodnosci komórkowej.

Single-cell transcriptomics (scRNA-Seq) is a breakthrough technology that has opened the way to characterizing gene expression with unprecedented resolution. It has enabled the discovery of the cellular diversity of organisms and tracing their developmental processes. A range of technological solutions have been developed to allow analysis of tens of thousands to even a million cells in a single experiment, as well as an extensive set of tools for bioinformatics analysis of the generated data. The wealth of information provided by scRNA-Seq and the possibility of using this method to study cells, organoids, tissues and even entire organisms determine its wide range of applications. In this paper, we present the experimental and computational parts of the scRNA-Seq procedure, as well as the most important applications of this technology in biomedicine, developmental biology and plant biology.

Small RNA fragments derived from multiple RNA classes - the missing element of multi-omics characteristics of the hepatitis C virus cell culture model.

BACKGROUND: A pool of small RNA fragments (RFs) derived from diverse cellular RNAs has recently emerged as a rich source of functionally relevant molecules. Although their formation and accumulation has been connected to various stress conditions, the knowledge on RFs produced upon viral infections is very limited. Here, we applied the next generation sequencing (NGS) to characterize RFs generated in the hepatitis C virus (HCV) cell culture model (HCV-permissive Huh-7.5 cell line). RESULTS: We found that both infected and non-infected cells contained a wide spectrum of RFs derived from virtually all RNA classes. A significant fraction of identified RFs accumulated to similar levels as miRNAs. Our analysis, focused on RFs originating from constitutively expressed non-coding RNAs, revealed three major patterns of parental RNA cleavage. We found that HCV infection induced significant changes in the accumulation of low copy number RFs, while subtly altered the levels of high copy number ones. Finally, the candidate RFs potentially relevant for host-virus interactions were identified. CONCLUSIONS: Our results indicate that RFs should be considered an important component of the Huh-7.5 transcriptome and suggest that the main factors influencing the RF biogenesis are the RNA structure and RNA protection by interacting proteins. The data presented here significantly complement the existing transcriptomic, miRnomic, proteomic and metabolomic characteristics of the HCV cell culture model.

Alterations in N-glycosylation of HCV E2 Protein in Children Patients with IFN-RBV Therapy Failure.

The glycosylation of viral envelope proteins plays an important role in virus biology and the immune response of the host to infection. Hepatitis C virus (HCV) envelope proteins E1 and E2, key players in virus entry and spread, are highly N-glycosylated and possess 4 (5 in certain genotypes) to 11 conserved glycosylation sites, respectively. Many published results based on recombinant proteins indicate that the glycan shield can mask the epitopes targeted by neutralizing antibodies. Glycan shifting within the conserved linear E2 region (412-423) could be one of the escape strategies used by HCV. In the present report, we isolated E2 genes from samples (collected before the IFN-RBV therapy) originating from pediatric patients infected with HCV gt 1a. We analyzed the biochemical properties of cloned E2 glycoprotein variants and investigated their glycosylation status. The sequencing of E2 genes isolated from patients who did not respond to therapy revealed mutations at N-glycosylation sites, thus leading to a lower molecular weight and a low affinity to both linear and conformational neutralizing antibodies. The loss of the glycosylation site within the conserved epitope (amino acid 417) impaired the binding with AP33, an antibody that potently neutralizes all genotypes of HCV. Our findings, based on clinical samples, confirm the influence of N-glycosylation aberrations on the antigenic and conformational properties of HCV E1/E2, which may possibly correlate with the outcome of therapy in patients.

Targeting sgRNA N secondary structure as a way of inhibiting SARS-CoV-2 replication.

SARS-CoV-2 is a betacoronavirus that causes COVID-19, a global pandemic that has resulted in many infections, deaths, and socio-economic challenges. The virus has a large positive-sense, single-stranded RNA genome of ∼30 kb, which produces subgenomic RNAs (sgRNAs) through discontinuous transcription. The most abundant sgRNA is sgRNA N, which encodes the nucleocapsid (N) protein. In this study, we probed the secondary structure of sgRNA N and a shorter model without a 3' UTR in vitro, using the SHAPE (selective 2'-hydroxyl acylation analyzed by a primer extension) method and chemical mapping with dimethyl sulfate and 1-cyclohexyl-(2-morpholinoethyl) carbodiimide metho-p-toluene sulfonate. We revealed the secondary structure of sgRNA N and its shorter variant for the first time and compared them with the genomic RNA N structure. Based on the structural information, we designed gapmers, siRNAs and antisense oligonucleotides (ASOs) to target the N protein coding region of sgRNA N. We also generated eukaryotic expression vectors containing the complete sequence of sgRNA N and used them to screen for new SARS-CoV-2 gene N expression inhibitors. Our study provides novel insights into the structure and function of sgRNA N and potential therapeutic tools against SARS-CoV-2.

Antibody-assisted selective isolation of Purkinje cell nuclei from mouse cerebellar tissue.

We developed a method that utilizes fluorescent labeling of nuclear envelopes alongside cytometry sorting for the selective isolation of Purkinje cell (PC) nuclei. Beginning with SUN1 reporter mice, we GFP-tagged envelopes to confirm that PC nuclei could be accurately separated from other cell types. We then developed an antibody-based protocol to make PC nuclear isolation more robust and adaptable to cerebellar tissues of any genotypic background. Immunofluorescent labeling of the nuclear membrane protein RanBP2 enabled the isolation of PC nuclei from C57BL/6 cerebellum. By analyzing the expression of PC markers, nuclear size, and nucleoli number, we confirmed that our method delivers a pure fraction of PC nuclei. To demonstrate its applicability, we isolated PC nuclei from spinocerebellar ataxia type 7 (SCA7) mice and identified transcriptional changes in known and new disease-associated genes. Access to pure PC nuclei offers insights into PC biology and pathology, including the nature of selective neuronal vulnerability.

Identification of stable, high copy number, medium-sized RNA degradation intermediates that accumulate in plants under non-stress conditions.

It is becoming increasingly evident that the RNA degradome is a crucial component of the total cellular RNA pool. Here, we present an analysis of the medium-sized RNAs (midi RNAs) that form in Arabidopsis thaliana. Our analyses revealed that the midi RNA fraction contained mostly 20-70-nt-long fragments derived from various RNA species, including tRNA, rRNA, mRNA and snRNA. The majority of these fragments could be classified as stable RNA degradation intermediates (RNA degradants). Using two dimensional polyacrylamide gel electrophoresis, we demonstrated that high copy number RNA (hcn RNA) degradants appear in plant cells not only during stress, as it was earlier suggested. They are continuously produced also under physiological conditions. The data collected indicated that the accumulation pattern of the hcn RNA degradants is organ-specific and can be affected by various endogenous and exogenous factors. In addition, we demonstrated that selected degradants efficiently inhibit translation in vitro. Thus, the results of our studies suggest that hcn RNA degradants are likely to be involved in the regulation of gene expression in plants.

RNA degradome--its biogenesis and functions.

RNA degradation is among the most fundamental processes that occur in living cells. The continuous decay of RNA molecules is associated not only with nucleotide turnover, but also with transcript maturation and quality control. The efficiency of RNA decay is ensured by a broad spectrum of both specific and non-specific ribonucleases. Some of these ribonucleases participate mainly in processing primary transcripts and in RNA quality control. Others preferentially digest mature, functional RNAs to yield a variety of molecules that together constitute the RNA degradome. Recently, it has become increasingly clear that the composition of the cellular RNA degradome can be modulated by numerous endogenous and exogenous factors (e.g. by stress). In addition, instead of being hydrolyzed to single nucleotides, some intermediates of RNA degradation can accumulate and function as signalling molecules or participate in mechanisms that control gene expression. Thus, RNA degradation appears to be not only a process that contributes to the maintenance of cellular homeostasis but also an underestimated source of regulatory molecules.

Engineered deaminases as a key component of DNA and RNA editing tools.

Over recent years, zinc-dependent deaminases have attracted increasing interest as key components of nucleic acid editing tools that can generate point mutations at specific sites in either DNA or RNA by combining a targeting module (such as a catalytically impaired CRISPR-Cas component) and an effector module (most often a deaminase). Deaminase-based molecular tools are already being utilized in a wide spectrum of therapeutic and research applications; however, their medical and biotechnological potential seems to be much greater. Recent reports indicate that the further development of nucleic acid editing systems depends largely on our ability to engineer the substrate specificity and catalytic activity of the editors themselves. In this review, we summarize the current trends and achievements in deaminase engineering. The presented data indicate that the potential of these enzymes has not yet been fully revealed or understood. Several examples show that even relatively minor changes in the structure of deaminases can give them completely new and unique properties.

How to explore what is hidden? A review of techniques for vascular tissue expression profile analysis.

The evolution of plants to efficiently transport water and assimilates over long distances is a major evolutionary success that facilitated their growth and colonization of land. Vascular tissues, namely xylem and phloem, are characterized by high specialization, cell heterogeneity, and diverse cell components. During differentiation and maturation, these tissues undergo an irreversible sequence of events, leading to complete protoplast degradation in xylem or partial degradation in phloem, enabling their undisturbed conductive function. Due to the unique nature of vascular tissue, and the poorly understood processes involved in xylem and phloem development, studying the molecular basis of tissue differentiation is challenging. In this review, we focus on methods crucial for gene expression research in conductive tissues, emphasizing the importance of initial anatomical analysis and appropriate material selection. We trace the expansion of molecular techniques in vascular gene expression studies and discuss the application of single-cell RNA sequencing, a high-throughput technique that has revolutionized transcriptomic analysis. We explore how single-cell RNA sequencing will enhance our knowledge of gene expression in conductive tissues.

2D-PAGE as an effective method of RNA degradome analysis.

The continuously growing interest in small regulatory RNA exploration is one of the important factors that have inspired the recent development of new high throughput techniques such as DNA microarrays or next generation sequencing. Each of these methods offers some significant advantages but at the same time each of them is expensive, laborious and challenging especially in terms of data analysis. Therefore, there is still a need to develop new analytical methods enabling the fast, simple and cost-effective examination of the complex RNA mixtures. Recently, increasing attention has been focused on the RNA degradome as a potential source of riboregulators. Accordingly, we attempted to employ a two-dimensional gel electrophoresis as a quick and uncomplicated method of profiling RNA degradome in plant or human cells. This technique has been successfully used in proteome analysis. However, its application in nucleic acids studies has been very limited. Here we demonstrate that two dimensional electrophoresis is a technique which allows one to quickly and cost-effectively identify and compare the profiles of 10-90 nucleotide long RNA accumulation in various cells and organs.

Dimensionality reduction by UMAP for visualizing and aiding in classification of imaging flow cytometry data.

Recent advances in imaging flow cytometry (IFC) have revolutionized high-throughput multiparameter analyses at single-cell resolution. Although enabling the discovery of population heterogeneities and the detection of rare events, IFC generates hyperdimensional datasets that demand innovative analytical approaches. Current methods work in a supervised manner, utilize only limited information content, or require large annotated reference datasets. Dimensionality reduction algorithms, including uniform manifold approximation and projection (UMAP), have been successfully applied to analyze the large number of parameters generated in various high-throughput techniques. Here, we apply a workflow incorporating UMAP to analyze different IFC datasets. We demonstrate that it out-competes other popular dimensionality reduction methods in speed and accuracy. Moreover, it enables fast visualization, clustering, and tagging of unannotated objects in large-scale experiments. We anticipate that our workflow will be a robust method to address complex IFC datasets, either alone or as an upstream addition to the deep learning approaches.

Molecular Composition of Serum Exosomes Could Discriminate Rectal Cancer Patients with Different Responses to Neoadjuvant Radiotherapy.

Identification of biomarkers that could be used for the prediction of the response to neoadjuvant radiotherapy (neo-RT) in locally advanced rectal cancer remains a challenge addressed by different experimental approaches. Exosomes and other classes of extracellular vesicles circulating in patients' blood represent a novel type of liquid biopsy and a source of cancer biomarkers. Here, we used a combined proteomic and metabolomic approach based on mass spectrometry techniques for studying the molecular components of exosomes isolated from the serum of rectal cancer patients with different responses to neo-RT. This allowed revealing several proteins and metabolites associated with common pathways relevant for the response of rectal cancer patients to neo-RT, including immune system response, complement activation cascade, platelet functions, metabolism of lipids, metabolism of glucose, and cancer-related signaling pathways. Moreover, the composition of serum-derived exosomes and a whole serum was analyzed in parallel to compare the biomarker potential of both specimens. Among proteins that the most properly discriminated good and poor responders were GPLD1 (AUC = 0.85, accuracy of 74%) identified in plasma as well as C8G (AUC = 0.91, accuracy 81%), SERPINF2 (AUC = 0.91, accuracy 79%) and CFHR3 (AUC = 0.90, accuracy 81%) identified in exosomes. We found that the proteome component of serum-derived exosomes has the highest capacity to discriminate samples of patients with different responses to neo-RT when compared to the whole plasma proteome and metabolome. We concluded that the molecular components of exosomes are associated with the response of rectal cancer patients to neo-RT and could be used for the prediction of such response.

At-C-RNA database, a one-stop source for information on circRNAs in Arabidopsis thaliana in a unified format.

Circular RNAs (circRNAs) are a large class of noncoding RNAs with functions that, in most cases, remain unknown. Recent genome-wide analysis of circRNAs using RNA-Seq has revealed that circRNAs are abundant and some of them conserved in plants. Furthermore, it has been shown that the expression of circRNAs in plants is regulated in a tissue-specific manner. Arabidopsis thaliana circular RNA database is a new resource designed to integrate and standardize the data available for circRNAs in a model plant A. thaliana, which is currently the best-characterized plant in terms of circRNAs. The resource integrates all applicable publicly available RNA-seq datasets. These datasets were subjected to extensive reanalysis and curation, yielding results in a unified format. Moreover, all data were normalized according to our optimized approach developed for circRNA identification in plants. As a result, the database accommodates circRNAs identified across organs and seedlings of wild-type A. thaliana and its single-gene knockout mutants for genes related to splicing. The database provides free access to unified data and search functionalities, thus enabling comparative analyses of A. thaliana circRNAs between organs, variants and studies for the first time. Database URLhttps://plantcircrna.ibch.poznan.pl/.

Hepatitis C virus quasispecies in chronically infected children subjected to interferon-ribavirin therapy.

Accumulating evidence suggests that certain features of hepatitis C virus (HCV), especially its high genetic variability, might be responsible for the low efficiency of anti-HCV treatment. Here, we present a bioinformatic analysis of HCV-1a populations isolated from 23 children with chronic hepatitis C (CHC) subjected to interferon-ribavirin therapy. The structures of the viral quasispecies were established based on a 132-amino-acid sequence derived from E1/E2 protein, including hypervariable region 1 (HVR1). Two types of HCV populations were identified. The first type, found in non-responders, contained a small number of closely related variants. The second type, characteristic for sustained responders, was composed of a large number of distantly associated equal-rank variants. Comparison of 445 HVR1 sequences showed that a significant number of variants present in non-responding patients are closely related, suggesting that certain, still unidentified properties of the pathogen may be key factors determining the result of CHC treatment.

Quasispecies nature of Pepino mosaic virus and its evolutionary dynamics.

Genetic variability is an essential feature of RNA viruses. It allows them to adapt to the ever-changing environmental conditions. Important biological properties of the viruses, their infectivity, adaptability, and host range, may also depend on the level of quasispecies diversity. Here, we present the analysis of the genetic polymorphism of Pepino mosaic virus (PepMV). The examined populations were isolated from the naturally infected tomato plants (Solanum lycopersicum). In order to determine the complexity of the PepMV populations, the number of different viral variants and their genetic diversity was established. Moreover, phylogenetic trees were created to depict relations between the identified variants. For the first time we have shown that the PepMV exists as a quasispecies. The observed level of genetic variability allows PepMV for a quick and flexible adaptation to different hosts. Our results suggest that the level of PepMV variability possibly influences the course of infection.

Expression Landscape of circRNAs in Arabidopsis thaliana Seedlings and Adult Tissues.

RNA-seq is currently the only method that can provide a comprehensive landscape of circular RNA (circRNAs) in the whole organism and its particular organs. Recent years have brought an increasing number of RNA-seq-based reports on plant circRNAs. Notably, the picture they revealed is questionable and depends on the applied circRNA identification and quantification techniques. In consequence, little is known about the biogenesis and functions of circRNAs in plants. In this work, we tested two experimental and six bioinformatics procedures of circRNA analysis to determine the optimal approach for studying the profiles of circRNAs in Arabidopsis thaliana. Then using the optimized strategy, we determined the accumulation of circular and corresponding linear transcripts in plant seedlings and organs. We observed that only a small fraction of circRNAs was reproducibly generated. Among them, two groups of circRNAs were discovered: ubiquitous and organ-specific. The highest number of circRNAs with significantly increased accumulation in comparison to other organs/seedlings was found in roots. The circRNAs in seedlings, leaves and flowers originated mainly from genes involved in photosynthesis and the response to stimulus. The levels of circular and linear transcripts were not correlated. Although RNase R treatment enriches the analyzed RNA samples in circular transcripts, it may also have a negative impact on the stability of some of the circRNAs. We also showed that the normalization of NGS data by the library size is not proper for circRNAs quantification. Alternatively, we proposed four other normalization types whose accuracy was confirmed by ddPCR. Moreover, we provided a comprehensive characterization of circRNAs in A. thaliana organs and in seedlings. Our analyses revealed that plant circRNAs are formed in both stochastic and controlled processes. The latter are less frequent and likely engage circRNA-specific mechanisms. Only a few circRNAs were organ-specific. The lack of correlation between the accumulation of linear and circular transcripts indicated that their biogenesis depends on different mechanisms.

Arabidopsis thaliana cbp80, c2h2, and flk Knockout Mutants Accumulate Increased Amounts of Circular RNAs.

Circular RNAs (circRNAs) are the products of the non-canonical splicing of pre-mRNAs. In contrast to humans and animals, our knowledge of the biogenesis and function of circRNAs in plants is very scarce. To identify proteins involved in plant circRNA generation, we characterized the transcriptomes of 18 Arabidopsis thaliana knockout mutants for genes related to splicing. The vast majority (>90%) of circRNAs were formed in more than one variant; only a small fraction of circRNAs was mutant-specific. Five times more circRNA types were identified in cbp80 and three times more in c2h2 mutants than in the wild-type. We also discovered that in cbp80, c2h2 and flk mutants, the accumulation of circRNAs was significantly increased. The increased accumulation of circular transcripts was not accompanied by corresponding changes in the accumulation of linear transcripts. Our results indicate that one of the roles of CBP80, C2H2 and FLK in splicing is to ensure the proper order of the exons. In the absence of one of the above-mentioned factors, the process might be altered, leading to the production of circular transcripts. This suggests that the transition toward circRNA production can be triggered by factors sequestering these proteins. Consequently, the expression of linear transcripts might be regulated through circRNA production.

Human- and virus-encoded microRNAs as potential targets of antiviral therapy.

It has recently been demonstrated that short RNA molecules, called microRNAs, are one of the major factors regulating the expression of human genes. There are several lines of evidence that microRNAs also play a key role in host-virus interactions. It is believed that both human- and virus-encoded miRNA will, in the nearest future, become very attractive targets of antiviral therapy.

Functional characterization of RNA fragments using high-throughput interactome screening.

Populations of small eukaryotic RNAs, in addition to relatively well recognized molecules such as miRNAs or siRNAs, also contain fragments derived from all classes of constitutively expressed non-coding RNAs. It has been recently demonstrated that the formation and accumulation of RNA fragments (RFs) is cell-/tissue-specific and depends on internal and external stimuli. Unfortunately, the mechanisms underlying RF biogenesis and function remain unclear. To better understand them, we employed RNA pull-down and mass spectrometry methods to characterize the interactions of seven RFs originating from tRNA, snoRNA and snRNA. By integrating our results with publicly available data on physical protein-protein interactions, we constructed an RF interactome network. We determined that the RF interactome comprises proteins generally different from those that interact with their parental full length RNAs. Proteins captured by the RFs were involved in mRNA splicing, tRNA processing, DNA recombination/replication, protein biosynthesis and carboxylic acid metabolism. Our data suggest that RFs can be endogenous aptamer-like molecules and potential players in recently revealed RNA-protein regulatory networks. SIGNIFICANCE: In the recent decade it has become evident that RNAs with well-known functions (for example tRNA, snoRNA or rRNA) can be cleaved to yield short fragments, whose role in cells remains only partially characterized. At the same time, unconventional interactions between mRNA and proteins without RNA-binding domains have been demonstrated, revealing novel layers of possible RNA-mediated regulation. Considering the above, we hypothesized that RNA fragments (RFs) can be endogenous aptamer-like molecules that unconventionally interact with proteins. In this study we identified protein partners of seven selected RFs. We found that RFs bind different set of proteins than their parental full length RNAs and identified proteins differentially bound by the particular RFs. These observations suggest biological relevance of the discovered interactions. Our data provide a novel perspective on the significance of RFs and point to this pool of molecules as to a rich collection of potential components of the recently discovered RNA-protein regulatory networks.

Mutations in human AID differentially affect its ability to deaminate cytidine and 5-methylcytidine in ssDNA substrates in vitro.

Activation-induced cytidine deaminase (AID) is known for its established role in antibody production. AID induces the diversification of antibodies by deaminating deoxycytidine (C) within immunoglobulin genes. The capacity of AID to deaminate 5-methyldeoxycytidine (5 mC) and/or 5-hydroxymethyldeoxycytidine (5 hmC), and consequently AID involvement in active DNA demethylation, is not fully resolved. For instance, structural determinants of AID activity on different substrates remain to be identified. To better understand the latter issue, we tested how mutations in human AID (hAID) influence its ability to deaminate C, 5 mC, and 5 hmC in vitro. We showed that each of the selected mutations differentially affects hAID's ability to deaminate C and 5 mC. At the same time, we did not observe hAID activity on 5 hmC. Surprisingly, we found that the N51A hAID mutant, with no detectable activity on C, efficiently deaminated 5 mC, which may suggest different requirements for C and 5 mC deamination. Homology modeling and molecular dynamics simulations revealed that the pattern of enzyme-substrate recognition is one of the important factors determining enzyme activity on C and 5 mC. Consequently, we have proposed mechanisms that explain why wild type hAID more efficiently deaminates C than 5 mC in vitro and why 5 hmC is not deaminated.