Transcriptomic Response of Mussel Gills After a Vibrio splendidus Infection Demonstrates Their Role in the Immune Response

Mussels (Mytilus galloprovincialis) are filter feeder bivalves that are constantly in contact with a wide range of microorganisms, some of which are potentially pathogenic. How mussels recognize and respond to pathogens has not been fully elucidated to date; therefore, we investigated the immune mechanisms that these animals employ in response to a bacterial bath infection from the surrounding water, mimicking the response that mussels mount under natural conditions. After the bath infection, mussels were able to remove the bacteria from their bodies and from the water tank. Accordingly, antibacterial activity was detected in gill extracts, demonstrating that this tissue plays a central role in removing and clearing potential pathogens. A transcriptomic study performed after a bath infection with Vibrio splendidus identified a total of 1,156 differentially expressed genes. The expression levels of genes contributing to a number of biological processes, such as immune response activation pathways and their regulation with cytokines, cell recognition, adhesion and apoptosis, were significantly modulated after infection, suggesting that the gills play important roles in pathogen recognition, as well as being activators and regulators of the mussel innate immune response. In addition to RNA-seq analysis, long non-coding RNAs and their neighboring genes were also analyzed and exhibited modulation after the bacterial challenge. The response of gills against bath infection was compared with the findings of a previous transcriptomic study on hemocytes responding to systemic infection, demonstrating the different and specific functions of gills. The results of this study indicate that recognition processes occur in the gill, thereby activating the effector agents of the immune response to overcome bacterial infection.

NDRindex: a method for the quality assessment of single-cell RNA-Seq preprocessing data

BACKGROUND

Single-cell RNA sequencing can be used to fairly determine cell types, which is beneficial to the medical field, especially the many recent studies on COVID-19. Generally, single-cell RNA data analysis pipelines include data normalization, size reduction, and unsupervised clustering. However, different normalization and size reduction methods will significantly affect the results of clustering and cell type enrichment analysis. Choices of preprocessing paths is crucial in scRNA-Seq data mining, because a proper preprocessing path can extract more important information from complex raw data and lead to more accurate clustering results.

RESULTS

We proposed a method called NDRindex (Normalization and Dimensionality Reduction index) to evaluate data quality of outcomes of normalization and dimensionality reduction methods. The method includes a function to calculate the degree of data aggregation, which is the key to measuring data quality before clustering. For the five single-cell RNA sequence datasets we tested, the results proved the efficacy and accuracy of our index.

CONCLUSIONS

This method we introduce focuses on filling the blanks in the selection of preprocessing paths, and the result proves its effectiveness and accuracy. Our research provides useful indicators for the evaluation of RNA-Seq data.

Optimized workflow for single-cell transcriptomics on infectious diseases including COVID-19

In December 2019, a new coronavirus, SARS-CoV-2, which causes the respiratory illness that led to the COVID-19 pandemic, was reported. In the face of such a new pathogen, special precautions must be taken to examine potentially infectious materials due to the lack of knowledge on disease transmissibility, infectivity, and molecular pathogenicity. Here, we present a complete and safe workflow for performing scRNA-seq experiments on blood samples of infected patients from cell isolation to data analysis using the micro-well based BD Rhapsody platform. For complete information on the use and execution of this protocol, please refer to Schulte-Schrepping et al. (2020).

RNA-Seq transcriptome data of human cells infected with influenza A/Puerto Rico/8/1934 (H1N1) virus

Human influenza remains a serious public health problem. This data article reports the transcriptome analysis data of human cell lines infected with influenza A/Puerto Rico/8/1934 (H1N1) virus. Mock-infected cells were included as controls. Human embryonic fibroblasts (MRC-5) and immortalized cell lines (A549, HEK293FT, WI-38 VA-13) were selected for RNA sequencing using Illumina NextSeq500 platform. Raw data were applied to the bioinformatic pipeline, which includes quality control with FastQC and MultiQC, adapter and quality trimming with Cutadapt, filtering to the genome of influenza A with STAR, transcript quantification with Salmon tool (GRCh38_RefSeq_Transcripts). Differential expressed genes were identified using R package DESeq2 with FDR-adjusted p-value < 0.001 and absolute value of log2(FC) > 1. Lists of differentially expressed genes is provided. The raw and processed RNA-seq data presented in this article were deposited to the European Nucleotide Archive via the ArrayExpress partner repository with the dataset accession number E-MTAB-9511 .

Nothobranchius annual killifishes

Annual fishes of the genus Nothobranchius inhabit ephemeral habitats in Eastern and Southeastern Africa. Their life cycle is characterized by very rapid maturation, a posthatch lifespan of a few weeks to months and embryonic diapause to survive the dry season. The species N. furzeri holds the record of the fastest-maturing vertebrate and of the vertebrate with the shortest captive lifespan and is emerging as model organism in biomedical research, evolutionary biology, and developmental biology. Extensive characterization of age-related phenotypes in the laboratory and of ecology, distribution, and demography in the wild are available. Species/populations from habitats differing in precipitation intensity show parallel evolution of lifespan and age-related traits that conform to the classical theories on aging. Genome sequencing and the establishment of CRISPR/Cas9 techniques made this species particularly attractive to investigate the effects genetic and non-genetic intervention on lifespan and aging-related phenotypes. At the same time, annual fishes are a very interesting subject for comparative approaches, including genomics, transcriptomics, and proteomics. The N. furzeri community is highly diverse and rapidly expanding and organizes a biannual meeting.

Unexpected Transcriptional Programs Contribute to Hippocampal Memory Deficits and Neuronal Stunting after Early-Life Adversity

Early-life adversity (ELA) is associated with lifelong memory deficits, yet the responsible mechanisms remain unclear. We impose ELA by rearing rat pups in simulated poverty, assess hippocampal memory, and probe changes in gene expression, their transcriptional regulation, and the consequent changes in hippocampal neuronal structure. ELA rats have poor hippocampal memory and stunted hippocampal pyramidal neurons associated with ~140 differentially expressed genes. Upstream regulators of the altered genes include glucocorticoid receptor and, unexpectedly, the transcription factor neuron-restrictive silencer factor (NRSF/REST). NRSF contributes critically to the memory deficits because blocking its function transiently following ELA rescues spatial memory and restores the dendritic arborization of hippocampal pyramidal neurons in ELA rats. Blocking NRSF function in vitro augments dendritic complexity of developing hippocampal neurons, suggesting that NRSF represses genes involved in neuronal maturation. These findings establish important, surprising contributions of NRSF to ELA-induced transcriptional programming that disrupts hippocampal maturation and memory function.

Expression analysis of microRNAs and their target mRNAs of testes with high and low sperm motility in domestic pigeons (Columba livia)

Sperm motility is one of the most important indicators to evaluate poultry fertility. In order to explore key molecular regulation roles related to sperm motility, we employed testicular RNA sequencing of pigeon. A total of 705 known and 385 novel microRNAs were identified. Compared with the low sperm motility group, four upregulated and two downregulated miRNAs in the high sperm motility group were identified. A total of 3567 target mRNAs were predicted and four target mRNAs were selected to validate by qPCR. The miRNA-mRNA interaction network analysis, indicated that mmu-miR-183-5p /FOXO1 and PC-3p-244994_31/CHDH pairs might affect sperm motility. GO and KEGG annotation analysis showed that target genes of differentially expressed miRNAs were related to serine/threonine kinase activity, ATP binding, Wnt and MAPK signaling pathway. The study provided a global miRNAs transcriptome of pigeon and a novel insight into the expression of the miRNAs in testes that associated with sperm motility.

Interleukin-8 Receptor 2 (IL-8R2)-Deficient Mice Are More Resistant to Pulmonary Coccidioidomycosis than Control Mice

The pathology of human coccidioidomycosis is granulomatous inflammation with many neutrophils surrounding ruptured spherules, but the chemotactic pathways that draw neutrophils into the infected tissues are not known. We previously showed that formalin-killed spherules (FKS) stimulate mouse macrophages to secret macrophage inflammatory protein 2 (MIP-2), which suggested that CXC ELR+ chemokines might be involved in neutrophil recruitment in vivo To test that hypothesis, we intranasally infected interleukin-8R2 (IL-8R2) (Cxcr2)-deficient mice on a BALB/c background with Coccidioides immitis RS. IL-8R2-deficient mice had fewer neutrophils in infected lungs than controls, but unexpectedly the IL-8R2-deficient mice had fewer organisms in their lungs than the control mice. Infected IL-8R2-deficient mouse lungs had higher expression of genes associated with lymphocyte activation, including the Th1 and Th17-related cytokines Ifnγ and Il17a and the transcription factors Stat1 and Rorc Additionally, bronchial alveolar lavage fluid from infected IL-8R2-deficient mice contained more IL-17A and interferon-γ (IFN-γ). We postulate that neutrophils in the lung directly or indirectly interfere with the development of a protective Th1/Th17 immune response to C. immitis at the site of infection.

Transcriptome analysis reveals modulation of the STAT family in PEDV-infected IPEC-J2 cells

BACKGROUND

Porcine epidemic diarrhea virus (PEDV) is a causative agent of serious viral enteric disease in suckling pigs. Such diseases cause considerable economic losses in the global swine industry. Enhancing our knowledge of PEDV-induced transcriptomic responses in host cells is imperative to understanding the molecular mechanisms involved in the immune response. Here, we analyzed the transcriptomic profile of intestinal porcine epithelial cell line J2 (IPEC-J2) after infection with a classical strain of PEDV to explore the host response.

RESULTS

In total, 854 genes were significantly differentially expressed after PEDV infection, including 716 upregulated and 138 downregulated genes. Functional annotation analysis revealed that the differentially expressed genes were mainly enriched in the influenza A, TNF signaling, inflammatory response, cytokine receptor interaction, and other immune-related pathways. Next, the putative promoter regions of the 854 differentially expressed genes were examined for the presence of transcription factor binding sites using the MEME tool. As a result, 504 sequences (59.02%) were identified as possessing at least one binding site of signal transducer and activator of transcription (STAT), and five STAT transcription factors were significantly induced by PEDV infection. Furthermore, we revealed the regulatory network induced by STAT members in the process of PEDV infection.

CONCLUSION

Our transcriptomic analysis described the host genetic response to PEDV infection in detail in IPEC-J2 cells, and suggested that STAT transcription factors may serve as key regulators in the response to PEDV infection. These results further our understanding of the pathogenesis of PEDV.

Transcriptional profiling of the L. monocytogenes PrfA regulon identifies six novel putative PrfA-regulated genes

The transcriptional activator Positive Regulatory Factor A (PrfA) regulates expression of genes essential for virulence in Listeria monocytogenes. To define the PrfA regulon, the 10403S wildtype (WT) strain, a constitutively active prfA* mutant, and an isogenic ∆prfA mutant were grown under PrfA-inducing conditions in a medium containing glucose-1-phosphate and pre-treated with 0.2% activated charcoal. RNA-seq-generated transcript levels were compared as follows: (i) prfA* and WT; (ii) WT and ∆prfA and (iii) prfA* and ∆prfA. Significantly higher transcript levels in the induced WT or constitutively active PrfA* were identified for 18 genes and 2 ncRNAs in at least one of the three comparisons. These genes included: (i) 10/12 of the genes previously identified as directly PrfA-regulated; (ii) 2 genes previously identified as PrfA-regulated, albeit likely indirectly; and (iii) 6 genes newly identified as PrfA-regulated, including one (LMRG_0 2046) with a σA-dependent promoter and PrfA box located within an upstream open reading frame. LMRG_0 2046, which encodes a putative cyanate permease, is reported to be downregulated by a σB-dependent anti-sense RNA. This newly identified overlap between the σB and PrfA regulons highlights the complexity of regulatory networks important for fine-tuning bacterial gene expression in response to the rapidly changing environmental conditions associated with infection.

annoFuse: an R Package to annotate, prioritize, and interactively explore putative oncogenic RNA fusions

BACKGROUND

Gene fusion events are significant sources of somatic variation across adult and pediatric cancers and are some of the most clinically-effective therapeutic targets, yet low consensus of RNA-Seq fusion prediction algorithms makes therapeutic prioritization difficult. In addition, events such as polymerase read-throughs, mis-mapping due to gene homology, and fusions occurring in healthy normal tissue require informed filtering, making it difficult for researchers and clinicians to rapidly discern gene fusions that might be true underlying oncogenic drivers of a tumor and in some cases, appropriate targets for therapy.

RESULTS

We developed annoFuse, an R package, and shinyFuse, a companion web application, to annotate, prioritize, and explore biologically-relevant expressed gene fusions, downstream of fusion calling. We validated annoFuse using a random cohort of TCGA RNA-Seq samples (N = 160) and achieved a 96% sensitivity for retention of high-confidence fusions (N = 603). annoFuse uses FusionAnnotator annotations to filter non-oncogenic and/or artifactual fusions. Then, fusions are prioritized if previously reported in TCGA and/or fusions containing gene partners that are known oncogenes, tumor suppressor genes, COSMIC genes, and/or transcription factors. We applied annoFuse to fusion calls from pediatric brain tumor RNA-Seq samples (N = 1028) provided as part of the Open Pediatric Brain Tumor Atlas (OpenPBTA) Project to determine recurrent fusions and recurrently-fused genes within different brain tumor histologies. annoFuse annotates protein domains using the PFAM database, assesses reciprocality, and annotates gene partners for kinase domain retention. As a standard function, reportFuse enables generation of a reproducible R Markdown report to summarize filtered fusions, visualize breakpoints and protein domains by transcript, and plot recurrent fusions within cohorts. Finally, we created shinyFuse for algorithm-agnostic interactive exploration and plotting of gene fusions.

CONCLUSIONS

annoFuse provides standardized filtering and annotation for gene fusion calls from STAR-Fusion and Arriba by merging, filtering, and prioritizing putative oncogenic fusions across large cancer datasets, as demonstrated here with data from the OpenPBTA project. We are expanding the package to be widely-applicable to other fusion algorithms and expect annoFuse to provide researchers a method for rapidly evaluating, prioritizing, and translating fusion findings in patient tumors.

Genetic Rescue of X-Linked Retinoschisis Mouse (Rs1-/y) Retina Induces Quiescence of the Retinal Microglial Inflammatory State Following AAV8-RS1 Gene Transfer and Identifies Gene Networks Underlying Retinal Recovery

To understand RS1 gene interaction networks in the X-linked retinoschisis (XLRS) mouse retina (Rs1^-/y), we analyzed the transcriptome by RNA sequencing before and after in vivo expression of exogenous retinoschisin (RS1) gene delivered by AAV8. RS1 is a secreted cell adhesion protein that is critical for maintaining structural lamination and synaptic integrity of the neural retina. RS1 loss-of-function mutations cause XLRS disease in young boys and men, with splitting ("schisis") of retinal layers and synaptic dysfunction that cause progressive vision loss with age. Analysis of differential gene expression profiles and pathway enrichment analysis of Rs1-KO (Rs1^-/y) retina identified cell surface receptor signaling and positive regulation of cell adhesion as potential RS1 gene interaction networks. Most importantly, it also showed massive dysregulation of immune response genes at early age, with characteristics of a microglia-driven proinflammatory state. Delivery of AAV8-RS1 primed the Rs1-KO retina toward structural and functional recovery. The disease transcriptome transitioned toward a recovery phase with upregulation of genes implicated in wound healing, anatomical structure (camera type eye) development, metabolic pathways, and collagen IV networks that provide mechanical stability to basement membrane. AAV8-RS1 expression also attenuated the microglia gene signatures to low levels toward immune quiescence. This study is among the first to identify RS1 gene interaction networks that underlie retinal structural and functional recovery after RS1 gene therapy. Significantly, it also shows that providing wild-type RS1 gene function caused the retina immune status to transition from a degenerative inflammatory phenotype toward immune quiescence, even though the transgene is not directly linked to microglia function. This study indicates that inhibition of microglial proinflammatory responses is an integral part of therapeutic rescue in XLRS gene therapy, and gene therapy might realize its full potential if delivered before microglia activation and photoreceptor cell death. Clinical Trials. gov Identifier NTC 02317887.

Comparison of transcriptional change of B. melitensis M5-90 after macrophage infection highlights the role of ribosome gene L31 in virulence

Brucella, a facultative intracellular bacterium, can survive and replicate in various cell types such as epithelial cell, fibroblasts and macrophage. Macrophage is the most important sites for the survival of Brucella in vivo. The mechanisms of pathogenesis are difficult to address, since the unknown virulence genes are still exist. RNA-seq is available to study transcriptional changes that occur during disease as a way to identify important virulence-related genes. Here we described and analyzed the transcriptional change of avirulent strain Brucella melitensis M5-90 (B. melitensis M5-90) during macrophage infection using RNA-seq technology. We detected 601 significant changed genes of which 428 were upregulated after infection. The upregulated gene L31 which involved in ribosome KEGG pathway was selected to illustrate its effect on virulence in a vaccine strain B. melitensis M5-90 and a virulent strain B. melitensis M28. Deletion of L31 significant attenuates the spleen colonization in model of M5-90 or M28 infection mouse at 7, 21 and 35 days post-infection (P < 0.05). We further examine the role of L31 in a macrophage cell infection model, and the result showed a significant reduction of intracellular M28ΔL31 cells at 48 h post-infection (P < 0.001). In total, our study provided a view of transcriptional landscape of B. melitensis M5-90 intracellular, and found L31 gene is required for the full virulence of B. melitensis.

Transcriptome Analyses Reveal Differential Transcriptional Profiles in Early- and Late-Dividing Porcine Somatic Cell Nuclear Transfer Embryos

Somatic cell nuclear transfer (SCNT) is not only a valuable tool for understanding nuclear reprogramming, but it also facilitates the generation of genetically modified animals. However, the development of SCNT embryos has remained an uncontrollable process. It was reported that the SCNT embryos that complete the first cell division sooner are more likely to develop to the blastocyst stage, suggesting their better developmental competence. Therefore, to better understand the underlying molecular mechanisms, RNA-seq of pig SCNT embryos that were early-dividing (24 h postactivation) and late-dividing (36 h postactivation) was performed. Our analysis revealed that early- and late-dividing embryos have distinct RNA profiles, and, in all, 3077 genes were differentially expressed. Gene ontology (GO)and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses revealed that early-dividing embryos exhibited higher expression in genes that participated in the meiotic cell cycle, while enrichment of RNA processing- and translation-related genes was found in late-dividing embryos. There are also fewer somatic memory genes such as FLRT2, ADAMTS1, and FOXR1, which are abnormally activated or suppressed in early-dividing cloned embryos. These results show that early-dividing SCNT embryos have different transcriptional profiles than late-dividing embryos. Early division of SCNT embryos may be associated with their better reprogramming capacity, and somatic memory genes may act as a reprogramming barrier in pig SCNT reprogramming.

Weighted Gene Co-Expression Network Coupled with a Critical-Time-Point Analysis during Pathogenesis for Predicting the Molecular Mechanism Underlying Blast Resistance in Rice

BACKGROUND

Rice blast, caused by the ascomycete fungus M. oryzae, is one of the most important diseases of rice. Although many blast resistance (R) genes have been identified and deployed in rice varieties, the molecular mechanisms responsible for the R gene-mediated defense responses are yet not fully understood.

RESULTS

In this study, we used comparative transcriptomic analysis to explore the molecular mechanism involved in Piz-t-mediated resistance in a transgenic line containing Piz-t (NPB-Piz-t) compared to Nipponbare (NPB). Clustering and principal component analysis (PCA) revealed that the time-point at 24-h post inoculation (hpi) was the most important factor distinguishing the four time-points, which consisted of four genes of mitogen-activated protein kinases (MAPKs) signaling pathway, one gene related to WRKY DNA-binding domain containing protein, five pathogenesis-related protein (OsPR1s) genes, and three genes of R proteins involving in the most significant protein-protein interaction (PPI) pathway. Using weighted gene co-expression network analysis (WGCNA) to investigate RNA-seq data across 0, 24, 48, and 72 hpi, nine modules with similar patterns expression pattern (SEP) and three modules with differential expression pattern (DEP) between NPB-Piz-t and NPB across 0, 24, 48, and 72 hpi with KJ201 (referred to as Piz-t-KJ201 and NPB-KJ201) were identified. Among these the most representative SEP green-yellow module is associated with photosynthesis, and DEP pink module comprised of two specific expressed nucleotide-binding domain and leucine-rich repeat (NLR) genes of LOC_Os06g17900 and LOC_Os06g17920 of Pi2/9 homologous, three NLR genes of LOC_Os11g11810, LOC_Os11g11770, and LOC_Os11g11920 which are putatively associated with important agronomic traits, and a B3 DNA binding domain containing protein related genes (LOC_Os10g39190). Knockout of LOC_Os10g39190 via CRISPR-Cas9 resulted in plant death at the seedling stage.

CONCLUSIONS

The research suggested that Piz-t and multiple NLR network might play important roles in the regulation of the resistance response in the Piz-t-KJ201 interaction system. The identified genes provide an NLR repository to study the rice-M. oryzae interaction system and facilitate the breeding of blast-resistant cultivars in the future.

Identification and Validation of Potential Pathogenic Genes and Prognostic Markers in ESCC by Integrated Bioinformatics Analysis

Esophageal squamous cell carcinoma (ESCC) is one of the most fatal malignancies of the digestive tract, but its underlying molecular mechanisms are not known. We aim to identify the genes involved in ESCC carcinogenesis and discover potential prognostic markers using integrated bioinformatics analysis. Three pairs of ESCC tissues and paired normal tissues were sequenced by high-throughput RNA sequencing (RNA-seq). Integrated bioinformatics analysis was used to identify differentially expressed coding genes (DECGs) and differentially expressed long non-coding RNA (lncRNA) genes (DELGs). A protein-protein interaction (PPI) network of DECGs was established using the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) website and visualized with Cytoscape. Survival analysis was conducted by log-rank tests to identify "hub" genes with potential prognostic value, and real-time reverse transcription-quantitative polymerase chain reaction (RT-qPCR) was conducted to assess expression of these genes in ESCC tissues. Transwell^TM assays were employed to examine the migration ability of cells after knockdown of LINC01614 expression, followed by investigation of epithelial-mesenchymal transition (EMT) by western blotting (WB). A total of 106 upregulated genes and 42 downregulated genes were screened out from the ESCC data sets. Survival analysis showed two hub protein-coding genes with higher expression in module 1 of the PPI network (SPP1 and BGN) and another three upregulated lncRNAs (LINC01614, LINC01415, NKILA) that were associated with a poor prognosis. High expression of SPP1, BGN, LINC01614, and LINC01415 in tumor samples was validated further by RT-qPCR. In vitro experiments show that knockdown of LINC01614 expression could significantly inhibit the migration of ESCC cells by regulating EMT, which was confirmed by WB. These results indicate that BGN, SPP1, LINC01614, and LINC01415 might be critical genes in ESCC and potential prognostic biomarkers.

Brain α-Tocopherol Concentration and Stereoisomer Profile Alter Hippocampal Gene Expression in Weanling Mice

BACKGROUND

Alpha-tocopherol (αT), the bioactive constituent of vitamin E, is essential for fertility and neurological development. Synthetic αT (8 stereoisomers; all rac-αT) is added to infant formula at higher concentrations than natural αT (RRR-αT only) to adjust for bio-potency differences, but its effects on brain development are poorly understood.

OBJECTIVES

The objective was to determine the impact of bio-potency-adjusted dietary all rac-αT versus RRR-αT, fed to dams, on the hippocampal gene expression in weanling mice.

METHODS

Male/female pairs of C57BL/6J mice were fed AIN 93-G containing RRR-αT (NAT) or all rac-αT (SYN) at 37.5 or 75 IU/kg (n = 10/group) throughout gestation and lactation. Male pups were euthanized at 21 days. Half the brain was evaluated for the αT concentration and stereoisomer distribution. The hippocampus was dissected from the other half, and RNA was extracted and sequenced. Milk αT was analyzed in separate dams.

RESULTS

A total of 797 differentially expressed genes (DEGs) were identified in the hippocampi across the 4 dietary groups, at a false discovery rate of 10%. Comparing the NAT-37.5 group to the NAT-75 group or the SYN-37.5 group to the SYN-75 group, small differences in brain αT concentrations (10%; P < 0.05) led to subtle changes (<10%) in gene expression of 600 (NAT) or 487 genes (SYN), which were statistically significant. Marked differences in brain αT stereoisomer profiles (P < 0.0001) had a small effect on fewer genes (NAT-37.5 vs. SYN-37.5, 179; NAT-75 vs. SYN-75, 182). Most of the DEGs were involved in transcription regulation and synapse formation. A network analysis constructed around known vitamin E interacting proteins (VIPs) revealed a group of 32 DEGs between NAT-37.5 vs. SYN-37.5, explained by expression of the gene for the VIP, protein kinase C zeta (Pkcz).

CONCLUSIONS

In weanling mouse hippocampi, a network of genes involved in transcription regulation and synapse formation was differentially affected by dam diet αT concentration and source: all rac-αT or RRR-αT.

Deep sequencing of the transcriptome from murine lung infected with H5N8 subtype avian influenza virus with combined substitutions I283M and K526R in PB2 gene

H5N8 subtype highly pathogenic avian influenza viruses (HPAIVs) pose a huge threat to poultry industry and general public health. Our previous study demonstrated that synergistic effect of 283M and 526R in PB2 gene was a critical factor for viral high pathogenicity in mammals. However, the potential pathogenic mechanism of the mutant virus is still unclear. Here, RNA-seq method was used to analyze the global host response of murine lungs after infecting with parental r-JY virus and JY-PB2-I283M-K526R mutant virus. We found that both amounts and the expression levels of host differentially expressed genes (DEGs) were higher in mutant virus-infected mice compared with the group of parental virus. Furthermore, the DEGs mainly related with innate immune response by GO and KEGG analysis. Especially, PB2-I283M-K526R mutation strongly induced a sharp expression of cytokine storm-related genes, including MX1, CXCL10, and IFN-γ, performed by qRT-PCR. We also found that PB2-I283M-K526R mutation accelerated the level of cell apoptosis by heat map analysis of apoptosis-related DEGs in lungs and apoptosis assay in vitro. Taken together, our data demonstrated that PB2-I283M-K526R of H5N8 subtype HPAIV exacerbated the innate immune response and the level of cell apoptosis, which might be a key pathogenic mechanism for the enhanced pathogenicity of mutants in mammals.

RNA-Seq-Based Transcriptomic Analysis of Saccharopolyspora spinosa Revealed the Critical Function of PEP Phosphonomutase in the Replenishment Pathway

Spinosyns, the secondary metabolites produced by Saccharopolyspora spinosa, are the active ingredients in a family of novel biological insecticides. Although the complete genome sequence of S. spinosa has been published, the transcriptome of S. spinosa remains poorly characterized. In this study, high-throughput RNA sequencing (RNA-seq) technology was applied to dissect the transcriptome of S. spinosa. Through transcriptomic analysis of different periods of S. spinosa growth, we found large numbers of differentially expressed genes and classified them according to their different functions. Based on the RNA-seq data, the CRISPR-Cas9 method was used to knock out the PEP phosphonomutase gene (orf 06952-4171). The yield of spinosyns A and D in S. spinosa-ΔPEP was 178.91 mg/L and 42.72 mg/L, which was 2.14-fold and 1.76-fold higher than that in the wild type (83.51 and 24.34 mg/L), respectively. The analysis of the mutant strains also verified the validity of the transcriptome data. The deletion of the PEP phosphonomutase gene leads to an increase in pyruvate content and affects the biosynthesis of spinosad. The replenishment of phosphoenol pyruvate in S. spinosa provides the substrate for the production of spinosad. We envision that these transcriptomic analysis results will contribute to the further study of secondary metabolites in actinomycetes.

Disruption of Plasmodium falciparum histidine-rich protein 2 may affect haem metabolism in the blood stage

Background

Haem is a key metabolic factor in the life cycle of the malaria parasite. In the blood stage, the parasite acquires host haemoglobin to generate amino acids for protein synthesis and the by-product haem for metabolic use. The malaria parasite can also synthesize haem de novo on its own. Plasmodium falciparum-specific histidine-rich protein 2 (PfHRP2) has a haem-binding site to mediate the formation of haemozoin, a biocrystallized form of haem aggregates. Notably, the gene regulates the mechanism of haemoglobin-derived haem metabolism and the de novo haem biosynthetic pathway in the Pfhrp2-disrupted parasite line during the intraerythrocytic stages.

Methods

The CRISPR/Cas9 system was used to disrupt the gene locus of Pfhrp2. DNA was extracted from the transgenic parasite, and PCR, Southern blotting and Western blotting were used to confirm the establishment of transgenic parasites. RNA-sequencing and comparative transcriptome analysis were performed to identify differences in gene expression between 3D7 and Pfhrp2^--3D7 parasites.

Results

Pfhrp2^- transgenic parasites were successfully established by the CRISPR/Cas9 system. A total of 964, 1261, 3138, 1064, 2512 and 1778 differentially expressed genes (DEGs) were identified in the six comparison groups, respectively, with 373, 520, 1499, 353, 1253 and 742 of these DEGs upregulated and 591, 741, 1639, 711, 1259 and 1036 of them downregulated, respectively. Five DEGs related to haem metabolism and synthesis were identified in the comparison groups at six time points (0, 8, 16, 24, 32, and 40 h after merozoite invasion). The genes encoding delta-aminolevulinic acid synthetase and ferrochelatase, both related to haem biosynthesis, were found to be significantly upregulated in the comparison groups, and those encoding haem oxygenase, stromal-processing peptidase and porphobilinogen deaminase were found to be significantly downregulated. No GO terms were significantly enriched in haem-related processes (Q value = 1).

Conclusion

Our data revealed changes in the transcriptome expression profile of the Pfhrp2^--3D7 parasite during the intraerythrocytic stages. The findings provide insight at the gene transcript level that will facilitate further research on and development of anti-malaria drugs.

Prediction of novel non-coding RNAs relevant for the growth of Pseudomonas putida in a bioreactor

Pseudomonas putida is a micro-organism with great potential for industry due to its stress-endurance traits and easy manipulation of the metabolism. However, optimization is still required to improve production yields. In the last years, manipulation of bacterial small non-coding RNAs (ncRNAs) has been recognized as an effective tool to improve the production of industrial compounds. So far, very few ncRNAs are annotated in P. putida beyond the generally conserved. In the present study, P. putida was cultivated in a two-compartment scale-down bioreactor that simulates large-scale industrial bioreactors. We performed RNA-Seq of samples collected at distinct locations and time-points to predict novel and potentially important ncRNAs for the adaptation of P. putida to bioreactor stress conditions. Instead of using a purely genomic approach, we have rather identified regions of putative ncRNAs with high expression levels using two different programs (Artemis and sRNA detect). Only the regions identified with both approaches were considered for further analysis and, in total, 725 novel ncRNAs were predicted. We also found that their expression was not constant throughout the bioreactor, showing different patterns of expression with time and position. This is the first work focusing on the ncRNAs whose expression is triggered in a bioreactor environment. This information is of great importance for industry, since it provides possible targets to engineer more effective P. putida strains for large-scale production.

Polo like kinase 1 expression in cervical cancer tissues generated from multiple detection methods

Background

Existing studies of PLK1 in cervical cancer had several flaws. The methods adopted by those studies of detecting PLK1 expression in cervical cancer were single and there lacks comprehensive evaluation of the clinico-pathological significance of PLK1 in cervical cancer.

Methods

A total of 303 cervical tissue samples were collected for in-house tissue microarrays. Immunohistochemistry was performed for evaluating PLK1 expression between cervical cancer (including cervical squamous cell carcinoma (CESC) and cervical adenocarcinoma) and non-cancer samples. The Expression Atlas database was searched for querying PLK1 expression in different cervical cancer cell lines and different tissues in the context of pan-cancer. Standard mean difference (SMD) was calculated and the summarized receiver's operating characteristics (SROC) curves were plotted for integrated tissue microarrays, exterior high-throughput microarrays and RNA sequencing data as further verification. The effect of PLK1 expression on the overall survival, disease-free survival and event-free survival of cervical cancer patients was analyzed through Kaplan Meier survival curves for cervical cancer patients from RNA-seq and GSE44001 datasets. The gene mutation and alteration status of PLK1 in cervical cancer was inspected in COSMIC and cBioPortal databases. Functional enrichment analysis was performed for genes correlated with PLK1 from aggregated RNA-seq and microarrays.

Results

A total of 963 cervical cancer samples and 178 non-cancer samples were collected from in-house tissue microarrays and exterior microarrays and RNA-seq datasets. The combined expression analysis supported overexpression of PLK1 in CESC, cervical adenocarcinoma and all types of cervical cancer (SMD = 1.59, 95%CI [0.56-2.63]; SMD = 2.99, 95%CI [0.75-5.24]; SMD = 1.57, 95% CI [0.85-2.29]) and the significant power of PLK1 expression in distinguishing CESC or all types of cervical cancer samples from non-cancer samples (AUC = 0.94, AUC = 0.92). Kaplan-Meier survival curves showed that the event-free survival rate of cervical cancer patients with higher expression of PLK1 was shorter than that of patients with lower PLK1 (HR = 2.020, P = 0.0197). Genetic alteration of PLK1 including missense mutation and mRNA low occurred in 6% of cervical cancer samples profiled in mRNA expression. Genes positively or negatively correlated with PLK1 were mainly assembled in pathways such as DNA replication, cell cycle, mismatch repair, Ras signaling pathway, melanoma, EGFR tyrosine kinase inhibitor resistance and homologous recombination (P < 0.05).

Conclusions

Here, we provided sufficient evidence of PLK1 overexpression in cervical cancer. The overexpression of PLK1 in cervical cancer and the contributory effect of it on clinical progression indicated the hopeful prospect of PLK1 as a biomarker for cervical cancer.

Transcriptome Profiling Provides Insights Into Potential Antagonistic Mechanisms Involved in Chaetomium globosum Against Bipolaris sorokiniana

Chaetomium globosum Kunze is recognized as a potential biocontrol fungus against spot blotch of wheat caused by Bipolaris sorokiniana. Its molecular mechanism of biocontrol activity and the biosynthetic pathways involved have not been yet elucidated. Here, global transcriptome profiling of C. globosum strain Cg2 during interaction with B. sorokiniana isolate BS112 using RNA-seq was performed in order to gain insights into the potential mechanisms of antagonism. The Illumina HiSeq platform (2 × 150 bp) yielded an average of 20-22 million reads with 50-58% GC. De novo assembly generated 45,582 transcripts with 27,957 unigenes. Transcriptome analysis displayed distinct expression profiles in the interaction (Cg2-BS112), out of which 6,109 unique differentially expressed genes were present. The predominant transcripts classified as genes involved in "catalytic activity" constituted 45.06%, of which 10.02% were associated with "hydrolytic activity" (GO:0008152), and similarly, in the biological process, 29.18% of transcripts were involved in "metabolic activity" (GO:0004096 and GO:0006979). Heat map and cluster categorization suggested an increase in the expression levels of genes encoding secondary metabolites like polyketide synthase (GO:0009058), S-hydroxymethyl glutathione dehydrogenase (GO:0006069), terpene cyclase (EC 4.2.3.-), aminotran_1_2 domain-containing protein (GO:0009058), and other hydrolytic CAZYmes such as the glycosyl hydrolase (GH) family (GH 13, GH 2, GH 31, and GH 81; GO:0005975), cellulase domain-containing protein, chitinases, β-1, 3-glucanases (GO:0004565), glucan endo-1,3-beta-glucanase (GO:0052861), and proteases (GO:0004177). The obtained RNA-seq data were validated by RT-qPCR using 20 randomly chosen genes, showing consistency with the RNA-seq results. The present work is worldwide the first effort to unravel the biocontrol mechanism of C. globosum against B. sorokiniana. It generated a novel dataset for further studies and facilitated improvement of the gene annotation models in the C. globosum draft genome.

Characterization and Transcriptome Analysis of Exosomal and Nonexosomal RNAs in Bovine Adipocytes

Exosomes are endosome-derived extracellular vesicles that allow intercellular communication. However, the biological significance of adipocyte exosomal RNAs remains unclear. To determine the role of RNAs from bovine adipocytes and exosomes in bovine adipogenesis, exosomal and nonexosomal RNAs were extracted from three bovine primary white adipocyte samples and then profiles were generated using DNBSEQ/BGISEQ-500 technology. The RNAome of adipocytes consisted of 12,082 mRNAs, 8589 lncRNAs, and 378 miRNAs for a higher complexity that that detected in exosomes, with 1083 mRNAs, 105 lncRNAs, and 48 miRNAs. Exosomal miRNA-mRNA and lncRNA–miRNA–mRNA networks were constructed and enrichment analysis was performed to predict functional roles and regulatory mechanisms. Our study provides the first characterization of RNAs from bovine adipocyte and exosomes. The findings reveal that some RNAs are specifically packaged in adipocyte-derived exosomes, potentially enabling crosstalk between adipocytes and/or other cells that is mediated by exosomes. Our results greatly expand our understanding of exosomal RNAs from bovine adipocytes, and provide a reference for future functional investigations of adipocyte exosomal RNAs under normal physiological conditions.

Comparison of lncRNA and mRNA expression in mouse brains infected by a wild-type and a lab-attenuated Rabies lyssavirus

Rabies is a lethal disease caused by Rabies lyssavirus, commonly known as rabies virus (RABV), and results in nearly 100 % death once clinical symptoms occur in human and animals. Long non-coding RNAs (lncRNAs) have been reported to be associated with viral infection. But the role of lncRNAs involved in RABV infection is still elusive. In this study, we performed global transcriptome analysis of both of lncRNA and mRNA expression profiles in wild-type (WT) and lab-attenuated RABV-infected mouse brains by using next-generation sequencing. The differentially expressed lncRNAs and mRNAs were analysed by using the edgeR package. We identified 1422 differentially expressed lncRNAs and 4475 differentially expressed mRNAs by comparing WT and lab-attenuated RABV-infected brains. Then we predicted the enriched biological pathways by the Gene Ontology (GO) and Kyoto Encyclopaedia of Genes and Genomes (KEGG) database based on the differentially expressed lncRNAs and mRNAs. Our analysis revealed the relationships between lncRNAs and RABV-infection-associated immune response and ion transport-related pathways, which provide a fresh insight into the potential role of lncRNA in immune evasion and neuron injury induced by WT RABV.