RNA-seq Libraries Research Articles

Transcriptomic and intervention evidence reveals domestic dogs as a promising model for anti-inflammatory investigation.

Domestic dogs have great potential to expand our understanding of the determinants of aging. To understand the aging pattern of domestic dogs and evaluate whether they can be used as an aging model, we performed RNA sequencing on white blood cells from domestic dogs aged 1-9 years and treated aged dogs with classical antiaging approaches. We obtained 30 RNA sequencing libraries and identified 61 age-associated genes with dynamic changes, the majority of which were related to metabolism and immune function, which may be predominant biomarkers for aging in dogs. We next treated aged dogs with canine mesenchymal stem cells (cMSCs), nicotinamide mononucleotide, and rapamycin to determine whether and how they responded to the antiaging interventions. The results showed that these treatments can significantly reduce the level of inflammatory factors (IL-6 and TNF-α). MSCs effectively improved the heart functions of aged dogs. Three key potential age-related genes (PYCR1, CCRL2, and TOX) were reversed by MSC treatment, two of which (CCRL2 and TOX) are implicated in immunity. Overall, we profiled the transcriptomic pattern of domestic dogs and revealed that they may be a good model of aging, especially in anti-inflammatory investigations.

Host adaptive radiation is associated with rapid virus diversification and cross-species transmission in African cichlid fishes

Adaptive radiations are generated through a complex interplay of biotic and abiotic factors. Although adaptive radiations have been widely studied in the context of animal and plant evolution, little is known about how they impact the evolution of the viruses that infect these hosts, which in turn may provide insights into the drivers of cross-species transmission and hence disease emergence. We examined how the rapid adaptive radiation of the cichlid fishes of African Lake Tanganyika over the last 10 million years has shaped the diversity and evolution of the viruses they carry. Through metatranscriptomic analysis of 2,242 RNA sequencing libraries, we identified 121 vertebrate-associated viruses among various tissue types that fell into 13 RNA and4 DNA virus groups. Host-switching was commonplace, particularly within the Astroviridae, Metahepadnavirus, Nackednavirus, Picornaviridae, and Hepacivirus groups, occurring more frequently than in other fish communities. A time-calibrated phylogeny revealed that hepacivirus diversification was not constant throughout the cichlid radiation but accelerated 2-3 million years ago, coinciding with a period of rapid cichlid diversification and niche packing in Lake Tanganyika, thereby providing more closely related hosts for viral infection. These data depict a dynamic virus ecosystem within the cichlids of Lake Tanganyika, characterized by rapid virus diversification and frequent host jumping, and likely reflecting their close phylogenetic relationships that lower the barriers to cross-species virus transmission.

The underlying molecular mechanism of ciliated epithelium dysfunction and TGF-β signaling in children with congenital pulmonary airway malformations

The aim of this study was to investigate the variation in gene expression in the complete transcripts of Congenitalpulmonary airwaymalformation (CPAM) of the lung using Next Generation Sequencing (NGS) technology. There were 20 cases involving children with CPAM were used for selection of study sample. NGS was used to establish RNA-Seq libraries for the two groups of samples separately, and both groups were conducted to differential expression analysis and Gene Ontology (GO) functional enrichment analysis. The pathways of the differential genes were analyzed to find the enriched target pathways. A total of 592 genes were expressed with significant differences (CPAM vs. normal tissue, P < 0.05). GO functional analysis of DEGs indicated that abnormal ciliary function played a role in the development of CPAM. Subsequently, analysis of these genes pathways showed the TGF-β signaling pathway was significantly enriched. Finally, the results of immunohistochemical analysis of some DEGs showed that a significant reduction in the expression of SMAD6, a gene related to the TGF-β signaling pathway, led to abnormal activation of the pathway. TGF-β signaling pathway involved in the evolution of the disease obtained by DEGs enrichment pathway analysis. SMAD6, a gene involved in this pathway, might be a potential biomarker for the diagnosis and treatment of CPAM.

Common mitochondrial deletions in RNA-Seq: evaluation of bulk, single-cell, and spatial transcriptomic datasets

Common mitochondrial DNA (mtDNA) deletions are large structural variants in the mitochondrial genome that accumulate in metabolically active tissues with age and have been investigated in various diseases. We applied the Splice-Break2 pipeline (designed for high-throughput quantification of mtDNA deletions) to human RNA-Seq datasets and describe the methodological considerations for evaluating common deletions in bulk, single-cell, and spatial transcriptomics datasets. A robust evaluation of 1570 samples from 14 RNA-Seq studies showed: (i) the abundance of some common deletions detected in PCR-amplified mtDNA correlates with levels observed in RNA-Seq data; (ii) RNA-Seq library preparation method has a strong effect on deletion detection; (iii) deletions had a significant, positive correlation with age in brain and muscle; (iv) deletions were enriched in cortical grey matter, specifically in layers 3 and 5; and (v) brain regions with dopaminergic neurons (i.e., substantia nigra, ventral tegmental area, and caudate nucleus) had remarkable enrichment of common mtDNA deletions.

First Report of Lily Virus X Infecting Paris polyphylla var. yunnanensis in Yunnan, China

Lily virus X (LVX) is a positive-sense ssRNA virus belonging to the genus Potexvirus in the family Alphaflexiviridae. LVX is known to infect plants of the genera Lilium and Tricyrtis in the family Liliacea. LVX was first reported in an asymptomatic lily (Lilium formosanum) from England (Stone, 1980), but has been shown to infect plants in the Netherlands (Chen et al. 2005), the United States (Jordan et al. 2008) and Japan (Nijo et al. 2018). To date, the complete genomes of two LVX isolates from the Netherlands and Japan have been reported. Paris polyphylla var. yunnanensis, known as Dianchonglou in China, is a perennial plant of the family Melanthiaceae (formerly belonging to the family Trillium). In China, its rhizome is commonly used as an antispasmodic agent for stroke and cancer treatment (Chang et al. 2017). From 2019 to 2022, leaf mottle and shrinkage which are typical symptoms of viral infections were observed on the leaves of P. polyphylla var. yunnanensis plants in Dianchonglou fields in Qujing, Yunnan. Disease incidence ranged from 19% to 45% across 5 fields (90 plants per field) in Qujing. To identify the possible viral pathogen(s) associated with the disease, the mirVanaTM miRNA isolation Kit was used to extract total RNA was from a mixed sample pool of 5 symptomatic leaf samples collected from the 5 fields. RNA sequencing library was constructed using TruSeqTM RNA sample preparation kit. Sequencing on the Illumina HiSeqTM 2500 platform (Illumina, USA) with 125-bp paired-end reads yielded 23,077,786 raw reads. 22,534,100 clean reads were obtained by removing reads of low quality and poly-N using Trimmomatic software (Bolger et al. 2014). By utilizing the paired-end splicing method in Trinity software (Grabherr et al. 2011) the the raw reads were De novo assembled into 184,596 contigs, of which 303 were related to viruses, including Paris mosaic necrosis virus (PMNV), Pear alphapartitivirus (PAPV), Dahlia mosaic virus (DMV), and Lily virus X (LVX). BLASTn analysis revealed that 12 contigs (lengths ranging from 344 nt to 5,981 nt, query cover 6% to 99%) were most similar (57.32% to 91.67% nt identities) to the genome sequences of LVX, suggesting a possible infection of LVX in the plants. To confirm the result, a full-length genomic sequence of LVX was obtained by reverse transcription polymerase chain reaction (RT-PCR) using specific primers designed based on the sequence of the assembled contigs. The PCR products were cloned into pGEM-T vector (Promega Corporation, USA) and sequenced using the Sanger method (Sangon Biotech, Shanghai, China). The obtained full-length genomic sequence of the LVX isolate (LVX-PP, accession number OM100017) was 5,981 nt in length. BLASTp analysis demonstrated that the putative Rep and CP of LVX-PP shared 76.27% to 81.05% and 80.81% to 81.82% aa sequence similarities with that of other LVX isolates, respectively. Maximum-likelihood phylogenetic trees inferred from the Rep and CP aa sequences showed that LVX-PP clustered closely with LVX isolates. The leaf samples were further analyzed using a lily virus X (LVX) ELISA kit (DEIAPV181, Creative Diagnostics, U.S.A.). Healthy P. polyphylla var. yunnanensis leaves were taken as a negative control and buffer solution as a blank control. The results showed a positive reaction for all five symptomatic plants (OD = 1.259 ± 0.007) relative to the negative (OD = 0.099) and blank (OD = 0.073) controls. These results indicate that LVX can infect P. polyphylla var. yunnanensis. To our knowledge, this is the first report that LVX has been detected in P. polyphylla var. yunnannensis. This study will serve as an important reference for the study of the host range of LVX. Further studies will be required to determine how LVX spreads between P. polyphylla var. yunnannensis and other host plants.

Abstract 147: Novel Inflammatory and Pro-Thrombotic mRNA Pathways Associated With Worse 3-month Outcome After Intracerebral Hemorrhage

Background: Intracerebral hemorrhage (ICH) has the highest disability among stroke survivors. The Recovery and Outcomes from Stroke (ROSE) study (U01NS100417) utilized RNA sequencing (RNAseq) to identify biomarkers associated with poor outcome at 3 months. Methods: Supratentorial, non-moribund, ≥ 18 years of age, spontaneous ICH patients were prospectively recruited in a multi-center study. Whole blood was obtained at baseline (&lt;15 days; RNA PAXgene tubes, Qiagen). 3-month outcomes were determined by modified Rankin Scale (mRS) 0-3 (good) vs 4-6 (poor). RNAseq libraries were constructed using the Illumina TruSeq Stranded Total RNA kit with Ribo-Zero, sequenced on an llumina NovaSeq 6000 platform using 1 x100 bp single end reads to generate 25-30 million reads at &gt;90% &gt;Q30 reads. Differential expression was tested using DESeq2 (negative binomial), adjusting for sequencing batch, age, sex, race, ICH volume, intraventricular hemorrhage, presenting Glasgow Coma Scale, and lobar vs non-lobar location. Pathway analyses and protein-protein interactions were computed on a set of genes (P &lt; 0.0001) using STRING, Cytoscape, MCODE, and Reactome. Results: We analyzed baseline blood samples from 278 patients (N = 156 vs 122; good vs poor outcomes at 3 months). Multiple pathways were associated with poor outcome, including those involved in clot formation and cell-cell/cell-matrix interactions (not shown). The statistically most significant pathway was observed for transcripts involved in neutrophil activation and degranulation. Conclusion: Variation in mRNA expression pathways associated with inflammation, coagulation and cellular interactions were associated with poor outcome after controlling for severity measures.

First Report of Three Bunya-Like Viruses, Apple Luteovirus 1, and Apple Hammerhead Viroid in Apples from Hakkari, Türkiye

Türkiye is a major apple fruit producer in the crossroads of Europe and the Middle East. Several reports have described the presence of multiple viruses affecting apple production in Türkiye, including apple stem grooving virus (ASGV), apple stem pitting virus (ASPV), apple chlorotic leafspot virus (ACLSV), and apple mosaic virus (ApMV) (Kurçman 1977; Fidan 1994; Çağlayan et al. 2003). However, there are no reports of the presence of the recently discovered bunya-like viruses citrus concave gum-associated virus (CCGaV), and apple rubbery wood viruses 1 and 2 (ARWV1 and 2), as well as apple luteovirus 1 (ALV-1), and apple hammerhead viroid (AHVd) in Türkiye, all of which have been previously reported in other apple-producing countries (Wright et al. 2018; Liu et al. 2018; Zhang et al. 2014). Leaves from one Gala, two Granny Smith, and one Golden Delicious apple trees showing mild symptoms of curling, chlorosis, and yellowing were collected from four different orchards in the province of Hakkari, southeast Türkiye during June 2022 and sent to USDA APHIS Plant Germplasm Quarantine Program (under permit) for virus and viroid HTS-based diagnostics. Total RNA was isolated using the RNeasy Plant Mini Kit (Qiagen) following the manufacturer’s guidelines to prepare RNAseq libraries using the TruSeq Stranded Total RNA Library Plant Kit (Illumina, Inc) as described in Malapi-Wight et al. (2021). Libraries were sequenced on the NextSeq500 sequencer (PE 2x75), and approximately 45 million reads were obtained per each sample on average. Bioinformatic analysis was performed as described in Costa et al. (2022) using Phytopipe, where unclassified pathogen-derived reads were de novo assembled and contigs were compared to the NCBI viral nucleotide and protein databases by BlastN and BlastX respectively using a 10-4 e-value cutoff. Nearly complete genome contigs were obtained for ACLSV (OR640150) and ASPV (OR640151) in all four samples and for ASGV (OR640152) in 3 of the 4 samples. The average BlastN identity to sequences in GenBank was 92.3% for ACLSV, ranging from 89-94 %. BlastN identity for ASPV was 86%, ranging from 81-92 % while the ASGV average BlastN identity was 98.2%. Nearly complete genomes with average genome coverage of 92.4% and 95.6% for RNA1 and RNA2 of CCGaV (OR640153 and OR640154), were found in two of the four samples with BlastN identity of 94.7% and 94.8% to GenBank sequences. Additionally, nearly complete genome of the large (L), medium (M), and small (S) segments for ARWV1 were found in two samples with average genome coverage of 99.9%, 99.4%, and 100% respectively and BlastN identity of 98.8%, 95.2%, and 98.4% (OR640155, OR640156, OR640157). ARWV2 contigs were also found in 1 sample where M and S segments had a coverage of 99.8% and BlastN identity of 95.4% (OR640158 and OR640159). The nearly complete genome of ALV-1 was also found in two of four samples with genome coverage of 94.1% and an average BlastN identity of 93.4% (OR640160). AHVd was found in one of the Granny Smith trees with 19,260 mapped reads to the reference GenBank MH049335.1 and identity of 98.3% (OR640149). The HTS findings of CCGaV, ARWV1, ARWV2, and ALV-1, from Türkiye were later confirmed by Sanger sequencing using custom-designed primers targeting the coat protein, the RNA-dependent RNA polymerase, or ~390bp for the AHVd genome (Supplementary Table 1). To further learn about the incidence of these agents, we tested 12 other apple samples from six different neighboring orchards and found them at 18.8% rate for CCGaV, 12.5% for both ARWV1 and ARWV2, 25% for ALV-1, and 37.5% for AHVd respectively. To our knowledge, this is the first report of the apple viruses CCGaV, ARWV1, ARWV2, and ALV-1, and the AHVd viroid in Türkiye. Further studies of the impact of these agents on orchard’s health are necessary, including their prevalence in high apple production regions of Türkiye.

Transfer RNA-derived small RNA tRF-Glu-CTC attenuates neointimal formation via inhibition of fibromodulin

Neointimal hyperplasia is a pathological vascular remodeling caused by abnormal proliferation and migration of subintimal vascular smooth muscle cells (VSMCs) following intimal injury. There is increasing evidence that tRNA-derived small RNA (tsRNA) plays an important role in vascular remodeling. The purpose of this study is to search for tsRNAs signature of neointima formation and to explore their potential functions. The balloon injury model of rat common carotid artery was replicated to induce intimal hyperplasia, and the differentially expressed tsRNAs (DE-tsRNAs) in arteries with intimal hyperplasia were screened by small RNA sequencing and tsRNA library. A total of 24 DE-tsRNAs were found in the vessels with intimal hyperplasia by small RNA sequencing. In vitro, tRF-Glu-CTC inhibited the expression of fibromodulin (FMOD) in VSMCs, which is a negative modulator of TGF-β1 activity. tRF-Glu-CTC also increased VSMC proliferation and migration. In vivo experiments showed that inhibition of tRF-Glu-CTC expression after balloon injury of rat carotid artery can reduce the neointimal area. In conclusion, tRF-Glu-CTC expression is increased after vascular injury and inhibits FMOD expression in VSMCs, which influences neointima formation. On the other hand, reducing the expression of tRF-Glu-CTC after vascular injury may be a potential approach to prevent vascular stenosis.

Precision run-on sequencing (PRO-seq) for microbiome transcriptomics.

Bacteria respond to environmental stimuli through precise regulation of transcription initiation and elongation. Bulk RNA sequencing primarily characterizes mature transcripts, so to identify actively transcribed loci we need to capture RNA polymerase (RNAP) complexed with nascent RNA. However, such capture methods have only previously been applied to culturable, genetically tractable organisms such as E. coli and B. subtilis. Here we apply precision run-on sequencing (PRO-seq) to profile nascent transcription in cultured E. coli and diverse uncultured bacteria. We demonstrate that PRO-seq can characterize the transcription of small, structured, or post-transcriptionally modified RNAs, which are often absent from bulk RNA-seq libraries. Applying PRO-seq to the human microbiome highlights taxon-specific RNAP pause motifs and pause-site distributions across non-coding RNA loci that reflect structure-coincident pausing. We also uncover concurrent transcription and cleavage of CRISPR guide RNAs and transfer RNAs. We demonstrate the utility of PRO-seq for exploring transcriptional dynamics in diverse microbial communities.

Exploring large-scale gene coexpression networks in peach (Prunus persica L.): a new tool for predicting gene function.

Peach is a model for Prunus genetics and genomics, however, identifying and validating genes associated to peach breeding traits is a complex task. A gene coexpression network (GCN) capable of capturing stable gene-gene relationships would help researchers overcome the intrinsic limitations of peach genetics and genomics approaches and outline future research opportunities. In this study, we created four GCNs from 604 Illumina RNA-Seq libraries. We evaluated the performance of every GCN in predicting functional annotations using an algorithm based on the 'guilty-by-association' principle. The GCN with the best performance was COO300, encompassing 21 956 genes. To validate its performance predicting gene function, we performed two case studies. In case study 1, we used two genes involved in fruit flesh softening: the endopolygalacturonases PpPG21 and PpPG22. Genes coexpressing with both genes were extracted and referred to as melting flesh (MF) network. Finally, we performed an enrichment analysis of MF network and compared the results with the current knowledge regarding peach fruit softening. The MF network mostly included genes involved in cell wall expansion and remodeling, and with expressions triggered by ripening-related phytohormones, such as ethylene, auxin, and methyl jasmonate. In case study 2, we explored potential targets of the anthocyanin regulator PpMYB10.1 by comparing its gene-centered coexpression network with that of its grapevine orthologues, identifying a common regulatory network. These results validated COO300 as a powerful tool for peach and Prunus research. This network, renamed as PeachGCN v1.0, and the scripts required to perform a function prediction analysis are available at https://github.com/felipecobos/PeachGCN.

Purification of Circular RNAs Using Poly(A) Tailing Followed by RNase R Digestion.

Thousands of eukaryotic protein-coding genes can be alternatively spliced to yield linear mRNAs and circular RNAs (circRNAs). Some circRNAs accumulate to higher levels than their cognate linear mRNAs, but the vast majority are expressed at low levels. Hence, for most circRNAs, only a handful of sequencing reads, if any, that span the backsplicing junction are observed in standard RNA-seq libraries. It thus has become common to use the 3'-5' exonuclease ribonuclease R (RNase R) to selectively degrade linear RNAs when aiming to prove transcript circularity or biochemically enrich circRNAs. However, RNase R fails to degrade linear RNAs with structured 3' ends or internal G-quadruplex structures. To overcome these shortcomings, we describe an improved protocol for circRNA purification from total RNA that employs a poly(A) tailing step prior to RNase R digestion, which is performed in a Li+ containing buffer (rather than K+) to destabilize G-quadruplexes. This biochemical method enables higher enrichment (two- to threefold) of circRNAs to be obtained compared to standard RNase R protocols due to more efficient removal of linear RNAs. By then performing quantitative RT-PCR (RT-qPCR) or generating RNA-seq libraries, the expression of individual circRNAs can be examined or the entire set of expressed circRNAs defined using established annotation algorithms. We describe step-by-step methods for annotating circRNAs using the CIRI2 and CIRCexplorer2 algorithms. In total, this overall approach can be used to enrich for circRNAs from any total RNA sample, thereby enabling one to quickly identify and validate circRNAs of interest for functional studies.

Single Cell Micro RNA Sequencing Library Preparation.

Micro RNAs represent important post-transcriptional regulators in health and are involved in the onset of many diseases. Therefore, the further characterization of physiological miRNA functions is an important basic research question, and miRNAs even have high potential as biomarkers both for prognosis and diagnosis. In order to exploit this potential, it is mandatory to accurately quantify the miRNA expression not only in bulk but also on the single-cell level. Here, we describe a protocol, which facilitates miRNA sequencing library preparation of very low input samples, single cells, and even clinical samples such as circulating tumor cells. The protocol can be combined with different single-cell isolation methods (e.g., micromanipulation and FACS sorting). After cell lysis, sequencing adapters are ligated to the miRNAs, other ncRNA species, and adapter dimers are reduced by exonuclease digest, the miRNA library is reverse transcribed, amplified, and purified. Furthermore, quality controls are described to select only high-quality samples for sequencing.

RNA Isolation from Ctenophores.

RNA-seq or transcriptome analysis of individual cells and small cell populations is essential for virtually any biomedical field. Here, we examine and discuss the different methods of RNA isolation specific to ctenophores. We present a convenient, inexpensive, and reproducible protocol for RNA-seq libraries that are designed forlow quantities of samples. We demonstrated these methods on early (one, two, four, eight cells) embryonic and developmental stages, tissues, and even a single aboral organ from the ctenophore Pleurobrachia bachei and other ctenophore species (e.g., Mnemiopsis, Bolinopsis, and Beroe).

Blockade of histone deacetylase activity affects transcription and splicing of neuronal and glial genes

The study of the molecular mechanisms underlying plastic processes in the nervous system is of great interest in modern neuroscience. It is important to understand that epigenetic modifications, which are crucial for the development and cellular differentiation, can also be involved in plastic processes in the adult nervous system. In our early work, we provided evidence that the expression of important memory-related genes, such as Prkcz and Prkci, can be regulated epigenetically [1]. In the current study, we extended the previous work to the systemic level by applying the RNA sequencing approach to evaluate global changes in the expression patterns of various genes during the induction of epigenetic rearrangements. Rat cortical neuron cultures were incubated with one of the nonselective histone deacetylase inhibitors (trichostatin A, TSA; sodium butyrate, NaB) to change the level of epigenetic regulation. Next, the total RNA was extracted and subjected to RNA-Seq libraries preparation and subsequent NGS-sequencing. Bioinformatics analysis of transcriptomic data revealed substantial overlapping of differentially expressed genes (DEGs) in NaB-treated and TSA-treated groups, indicating that different histone deacetylase (HDAC) inhibitors induce transcriptional changes in primary neuron cultures through common regulatory pathways irrespective of chemical structure of applied inhibitor. We found that histone deacetylase blockade is accompanied by a transition from proliferative processes to cellular differentiation. Gene Ontology (GO) analysis of DEGs datasets revealed that the upregulated genes were engaged in cell differentiation and specialization, tissue and embryonic morphogenesis, and the development of various peripheral tissues and organs. On the contrary, genes that reduce the expression under induction of epigenetic rearrangements were involved in biological processes associated with cell proliferation and, most interestingly, the specialization of various brain cells (neurons, astrocytes, oligodendrocytes). It was shown that the expression of a number of glial markers typical for astrocytes and oligodendrocytes was significantly reduced after application of HDAC inhibitors, which was also confirmed by quantitative PCR using specific primer pairs on selected target genes. During the data analysis, we also found a significant decrease in the expression of various neuronal markers associated with the cytoskeleton, the organization of pre- and postsynaptic endings, synaptic transmission. It is known that fine-tuning of various processes in the central nervous system is due to the production of different isoforms of proteins from the same gene due to the process of alternative splicing of the resulting mRNA. According to the literature, epigenetic rearrangements create a certain environment for the regulation of alternative splicing of different genes [2, 3]. It is shown that production of alternative isoforms can play an important role in various plastic processes [4, 5]. Therefore, using IsoformSwitchAnalyzeR and DEXSeq packages we analyzed the possibility of alternative splicing during the induction of epigenetic rearrangements in rat cortical neuron cultures by evaluating the abundance of various transcripts based on exon usage. We found that some glial genes and a large number of neuronal genes, especially those associated with postsynaptic organization and cell communication, were alternatively spliced after application of histone deacetylase inhibitors. Inhibition of HDAC activity in cortical neuron cultures mainly affected the choice of alternative transcription starts (ATSS) and terminators (ATTS), and to a lesser extent alternative splicing of exons. Obtained results were selectively confirmed by the quantitative PCR using specific pairs of primers for individual exons of different transcript isoforms. Thus, within this study, it was found that histone deacetylases play an important role in the specialization of various brain cells, and the suppression of their activity affects the expression and alternative splicing of various glial and neuronal marker genes. We do not exclude that global transcriptome changes caused by alternative splicing will lead to qualitative rearrangements of the neuron network, and this is the direction of future research.

Discovery and Identification of Viruses Infecting Oomycetes.

This chapter describes protocols suitable for the detection and identification of RNA viruses infecting oomycetes (so-called water molds of Kingdom Heterokonta, Stramenopila), focusing on species of Phytophthora and exemplified by P. fragariae. The protocol includes laboratory procedures for oomycete cultivation and total RNA extraction from harvested mycelia, followed by instructions on suitable parameters given for sequencing companies on ribosomal RNA depletion, cDNA library preparation, and total RNA-sequencing (RNA-Seq). We also describe the bioinformatics steps needed for de novo assembly of raw reads into contigs, removal of host-associated contigs, and virus identification by database searches, as well as host validation by RT-PCR. All steps are described using an exemplar RNA-Seq library containing a yet undescribed fusagravirus hosted by a P. fragariae isolate.

Distinct cultured microglia characteristics distinguish cognitive SuperAgers from their cognitively normal peers

AbstractBackgroundSuperAgers are individuals ≥80 years with exceptional episodic memory equal to or better than individuals 20 to 30 years their junior. SuperAgers display resistance to age‐related neurofibrillary degeneration. Qualitative observations suggested that SuperAger cultured microglia display higher rates of proliferation compared to age‐matched controls. The purpose of this study was a comprehensive and quantitative investigation of potential differences between SuperAger and control microglia.MethodApproximately 6.6×103 microglia from the frontal cortex of 5 SuperAgers and 5 controls were seeded and allowed to proliferate to 80% confluence prior to passage (P). Number of cells at confluence, number of days to confluence, and highest passage reached were recorded. Gene expression profiles were investigated in P4 cultures. RNASeq libraries were prepared using the Lexogen QuantSeq 3’ RNA Library kit and assessed using the Agilent High Sensitivity DNA chip. Bioinformatic analysis was performed using the BlueBee integrated data analysis pipeline. Microglia function was assayed through uptake of fluorescently tagged oligomeric and fibrillar amyloid‐β (Aβ42) and production of reactive oxygen species (ROS) in response to the two Aβ conformations.ResultAt each passage, SuperAger cultures contained greater numbers of microglia when compared with controls (by 150% at P5, p&lt;0.01) and reached 80% confluent faster (by 108%). On average, control microglia ceased proliferating at P9, while SuperAger microglia proliferated up to P12. RNAseq analysis identified 434 genes differentially expressed in the two groups (p&lt;0.05). Upregulated genes in SuperAgers included genes involved in DNA repair and in membrane trafficking. SuperAgers exhibited downregulation of several transcripts involved in age‐related disease pathogenesis, including Adap2, MAP2, RBM20, CNTN3, and genes that regulate mitochondrial oxidative burst. Preliminary functional analysis revealed more efficient uptake and response to oligomeric when compared with fibrillar Aβ. SuperAger microglia displayed greater ROS production when compared with controls. No differences were observed in Aβ uptake.ConclusionThese results point to distinct characteristics of microglia that may contribute to the biological bases of the SuperAging phenotype. Upregulation of genes involved in DNA repair is likely to slow microglia senescence, allowing increased proliferation in SuperAgers. Additional studies are needed to validate and extend these findings.

Heparin-binding epidermal growth factor and fibroblast growth factor 2 rescue Müller glia-derived progenitor cell formation in microglia- and macrophage-ablated chick retinas.

Recent studies have demonstrated the impact of pro-inflammatory signaling and reactive microglia/macrophages on the formation of Müller glial-derived progenitor cells (MGPCs) in the retina. In chick retina, ablation of microglia/macrophages prevents the formation of MGPCs. Analyses of single-cell RNA-sequencing chick retinal libraries revealed that quiescent and activated microglia/macrophages have a significant impact upon the transcriptomic profile of Müller glia (MG). In damaged monocyte-depleted retinas, MG fail to upregulate genes related to different cell signaling pathways, including those related to Wnt, heparin-binding epidermal growth factor (HBEGF), fibroblast growth factor (FGF) and retinoic acid receptors. Inhibition of GSK3β, to simulate Wnt signaling, failed to rescue the deficit in MGPC formation, whereas application of HBEGF or FGF2 completely rescued the formation of MGPCs in monocyte-depleted retinas. Inhibition of Smad3 or activation of retinoic acid receptors partially rescued the formation of MGPCs in monocyte-depleted retinas. We conclude that signals produced by reactive microglia/macrophages in damaged retinas stimulate MG to upregulate cell signaling through HBEGF, FGF and retinoic acid, and downregulate signaling through TGFβ/Smad3 to promote the reprogramming of MG into proliferating MGPCs.

Study of microRNAs and their targets in exosome from Alzheimer’s and Frontotemporal Dementia patients

AbstractBackgroundA central role of RNA metabolism and non‐coding RNAs biogenesis has emerged in Alzheimer’s disease (AD) and Frontotemporal Dementia (FTD). We investigated miRNA cargo of Small Extracellular Vesicles (SEVs) derived from plasma of FTD and AD patients, and Healthy Controls (HC). The purpose was to evaluate deregulated miRNAs in patients to identify new peripheral biomarkers and to determine mRNA targets involved in FTD and AD pathogenesis.MethodSEVs were isolated from plasma of 13 FTD patients, 20 AD and 20 HC by differential centrifugation and characterized by Nanosight and Western Blot. MicroRNA libraries were generated using Small RNA Seq Library Prep Kit (Lexogen) and sequenced on a NextSeq 500 (Illumina). Interaction prediction was carried out on TarBase v.8 database.ResultWe found a total of 339 Differentially Expressed (DE) microRNAs in FTD, and 291 DE microRNAs in AD compared to controls. Interestingly, 74 miRNAs were commonly deregulated: 38 had the same trend and 36 had an opposite trend. In order to find a miRNA biomarker able to discern FTD and AD, we selected 8 DE miRNAs with opposite trend based on fold changes difference and targets. In fact these miRNAs were found to have validated mRNAs target that may be of interest for molecular study of dementia, such as APP and PTEN.ConclusionOur data highlight the importance of miRNAs cargo examination in EVs of FTD and AD patients. In fact, their potential is exploitable both for biomarkers discovery and for study of gene expression alteration in dementia pathogenesis.

Computational Analysis Tutorial for Chimeric Small Noncoding RNA: Target RNA Sequencing Libraries.

An understanding of the in vivo gene regulatory interactions of small noncoding RNAs (sncRNAs), such as microRNAs (miRNAs), with their target RNAs has been advanced in recent years by biochemical approaches which use cross-linking followed by ligation to capture sncRNA:target RNA interactions through the formation of chimeric RNAs and subsequent sequencing libraries. While datasets from chimeric RNA sequencing provide genome-wide and substantially less ambiguous input than miRNA prediction software, distilling this data into meaningful and actionable information requires additional analyses and may dissuade investigators lacking a computational background. This report provides a tutorial to support entry-level computational biologists in installing and applying a recent open-source software tool: Small Chimeric RNA Analysis Pipeline (SCRAP). Platform requirements, updates, and an explanation of pipeline steps and manipulation of key user-input variables is provided. Reducing a barrier for biologists to gain insights from chimeric RNA sequencing approaches has the potential to springboard discovery-based investigations of regulatory sncRNA:target RNA interactions in multiple biological contexts.

Molecular Mechanisms Underlying Abnormal Activation of T Cells in Idiopathic Aplastic Anemia and Hypoplastic Myelodysplastic Neoplasms

Idiopathic aplastic anemia (AA) and hypoplastic myelodysplastic neoplasms (MDS-h) are severe hematopoietic disorders characterized by pancytopenia and hypoplastic bone marrow. There is compelling evidence that these distinct clinical entities share a common pathophysiology based on the damage of hematopoietic stem and progenitor cells (HSPCs) by cytotoxic T cells. Expanded T cells overproduce proinflammatory cytokines, resulting in decreased proliferation and increased apoptosis of HSPCs. To uncover the molecular mechanisms underlying this abnormal immune response, we used RNA-Seq for the transcriptome analysis of T cells from patients with idiopathic AA and MDS-h. CD3+ cells were isolated from peripheral blood of 15 patients with AA or MDS-h at diagnosis and 6 healthy blood donors. The RNA-Seq library was constructed by NEB Next Ultra II Directional RNA Library Prep Kit and the data were analysed by R software 4.0.2. Functional annotation of differently expressed genes was performed using DAVID database and Gene Set Enrichment Analysis (GSEA). Principal component analysis-based clustering of RNA-Seq data defined the three major components for AA and MDS-h patients, and healthy individuals (Fig. 1A). As expected, the control samples formed a distinct group, meanwhile patient samples partially clustered together. The same results were obtained using unsupervised hierarchical clustering. The differential expression analysis (DEA) identified 256 significantly upregulated and 96 downregulated genes (ІlogFCІ&amp;gt;1, FDR&amp;lt;0.05) in patients compared to controls. Notably, the DEA detected many long noncoding RNAs (lncRNAs) including upregulated lncRNAs associated with immunological disturbances (e.g. CDC42-IT1 and NEAT1), and oncogenesis (e.g. ACAP2-IT1, and FAM238A). Particularly, NEAT1 positively regulates differentiation of CD4+ T cells into Th17 cells through STAT3 protein. Autoreactive T cells may be further stimulated by increased expression of S100A8/A9 detected in the patient cells. Gene Ontology (GO) analysis annotated the upregulated genes in AA/MDS-h to biological processes associated with oxygen transport, B cell receptor signaling pathway, innate immune response, positive regulation of inflammatory response, and apoptotic process, etc. (adjusted p&amp;lt;0.05) (Fig. 1B). The downregulated genes were significantly enriched in processes related to adaptive immune response, cell surface receptor signaling pathway, mitochondrial ATP synthesis coupled proton transport, etc. (adjusted p&amp;lt;0.05) (Fig.1B). B cell receptor signaling pathway was enriched by upregulated genes such as BANK1, MEF2C, and SYK whose deregulation is associated with autoimmune disorders. Herein, the activation of apoptosis is likely driven via AP-1 positive modulators that showed increased expressions in the patients. These modulators are essential for T cell activation and regulatory T cell differentiation. Pathway enrichment analysis determined the following signaling pathways (adjusted p&amp;lt;0.05): B cell receptor signaling pathway, osteoclast differentiation, oxidative phosphorylation, chemical carcinogenesis - reactive oxygen species, and cytokine-cytokine receptor interaction. The GSEA using the Hallmark gene sets showed that DNA Repair and Oxidative Phosphorylation were significantly negatively correlated with patient phenotype (Fig. 1C). Both GSEA and pathway analysis revealed that oxidative phosphorylation was supressed in patient-derived T cells. This is consistent with previous metabolic findings demonstrating that a transition from oxidative phosphorylation to glycolysis is a sign of T cell activation and it is critical for support of rapid cell growth. In conclusion, our study suggests the possible mechanisms of aberrant T cellular-immunity in idiopathic AA and MDS-h. Consistent with abnormal biological features of T cells, a large number of deregulated genes in the patients was involved in immune and inflammatory responses. Furthermore, our analyses revealed novel pathways, such as oxidative phosphorylation, as well as candidates ( TLR2, and CDC42-IT1) for functional studies. The deregulation of multiple lncRNAs in the patients indicates their implication in the molecular pathogenesis and thus their roles in bone marrow failure need to be explored. Supported by AZV (NU21-03-00565), and MH CZ-DRO (UHKT, 00023736).