Pseudomonas aeruginosa Activates Quorum Sensing, Antioxidant Enzymes and Type VI Secretion in Response to Oxidative Stress to Initiate Biofilm Formation and Wound Chronicity.

Pseudomonas aeruginosa (PA) is an opportunistic pathogen frequently isolated from cutaneous chronic wounds. How PA, in the presence of oxidative stress (OS), colonizes chronic wounds and forms a biofilm is still unknown. The purpose of this study is to investigate the changes in gene expression seen when PA is challenged with the high levels of OS present in chronic wounds. We used a biofilm-forming PA strain isolated from the chronic wounds of our murine model (RPA) and performed a qPCR to obtain gene expression patterns as RPA developed a biofilm in vitro in the presence of high levels of OS, and then compared the findings in vivo, in our mouse model of chronic wounds. We found that the planktonic bacteria under OS conditions overexpressed quorum sensing genes that are important for the bacteria to communicate with each other, antioxidant stress genes important to reduce OS in the microenvironment for survival, biofilm formation genes and virulence genes. Additionally, we performed RNAseq in vivo and identified the activation of novel genes/pathways of the Type VI Secretion System (T6SS) involved in RPA pathogenicity. In conclusion, RPA appears to survive the high OS microenvironment in chronic wounds and colonizes these wounds by turning on virulence, biofilm-forming and survival genes. These findings reveal pathways that may be promising targets for new therapies aimed at disrupting PA-containing biofilms immediately after debridement to facilitate the treatment of chronic human wounds.

spatialHeatmap: visualizing spatial bulk and single-cell assays in anatomical images.

Visualizing spatial assay data in anatomical images is vital for understanding biological processes in cell, tissue, and organ organizations. Technologies requiring this functionality include traditional one-at-a-time assays, and bulk and single-cell omics experiments, including RNA-seq and proteomics. The spatialHeatmap software provides a series of powerful new methods for these needs, and allows users to work with adequately formatted anatomical images from public collections or custom images. It colors the spatial features (e.g. tissues) annotated in the images according to the measured or predicted abundance levels of biomolecules (e.g. mRNAs) using a color key. This core functionality of the package is called a spatial heatmap plot. Single-cell data can be co-visualized in composite plots that combine spatial heatmaps with embedding plots of high-dimensional data. The resulting spatial context information is essential for gaining insights into the tissue-level organization of single-cell data, or vice versa. Additional core functionalities include the automated identification of biomolecules with spatially selective abundance patterns and clusters of biomolecules sharing similar abundance profiles. To appeal to both non-expert and computational users, spatialHeatmap provides a graphical and a command-line interface, respectively. It is distributed as a free, open-source Bioconductor package (https://bioconductor.org/packages/spatialHeatmap) that users can install on personal computers, shared servers, or cloud systems.

MYC acetylated lysine residues drive oncogenic cell transformation and regulate select genetic programs for cell adhesion-independent growth and survival.

The MYC oncogenic transcription factor is acetylated by the p300 and GCN5 histone acetyltransferases. The significance of MYC acetylation and the functions of specific acetylated lysine (AcK) residues have remained unclear. Here, we show that the major p300-acetylated K148(149) and K157(158) sites in human (or mouse) MYC and the main GCN5-acetylated K323 residue are reversibly acetylated in various malignant and nonmalignant cells. Oncogenic overexpression of MYC enhances its acetylation and alters the regulation of site-specific acetylation by proteasome and deacetylase inhibitors. Acetylation of MYC at different K residues differentially affects its stability in a cell type-dependent manner. Lysine-to-arginine substitutions indicate that although none of the AcK residues is required for MYC stimulation of adherent cell proliferation, individual AcK sites have gene-specific functions controlling select MYC-regulated processes in cell adhesion, contact inhibition, apoptosis, and/or metabolism and are required for the malignant cell transformation activity of MYC. Each AcK site is required for anchorage-independent growth of MYC-overexpressing cells in vitro, and both the AcK148(149) and AcK157(158) residues are also important for the tumorigenic activity of MYC transformed cells in vivo. The MYC AcK site-specific signaling pathways identified may offer new avenues for selective therapeutic targeting of MYC oncogenic activities.

Genome-wide identification of core components of ABA signaling and transcriptome analysis reveals gene circuits involved in castor bean (Ricinus communis L.) response to drought.

Castor bean (Ricinus communis L.) can withstand long periods of water deficit and high temperatures, and therefore has been recognized as a drought-resistant plant species, allowing the study of gene networks involved in drought response and tolerance. The identification of genes networks related to drought response in this plant may yield important information in the characterization of molecular mechanisms correlating changes in the gene expression with the physiological adaptation processes. In this context, gene families related to abscisic acid (ABA) signaling play a crucial role in developmental and environmental adaptation processes of plants to drought stress. However, the families that function as the core components of ABA signaling, as well as genes networks related to drought response, are not well understood in castor bean. In this study 7 RcPYL, 63 RcPP2C, and 6 RcSnRK2 genes were identified in castor bean genome, which was further supported by chromosomal distribution, gene structure, evolutionary relationships, and conserved motif analyses. The castor bean general expression profile was investigated by RNAseq in root and leaf tissues in response to drought stress. These analyses allowed the identification of genes differentially expressed, including genes from the ABA signaling core, genes related to photosynthesis, cell wall, energy transduction, antioxidant response, and transcription factors. These analyses provide new insights into the core components of ABA signaling in castor bean, allow the identification of several molecular responses associated with the high physiological adaptation of castor bean to drought stress, and contribute to the identification of candidate genes for genetic improvement.

Transcriptomic Evidence That Switching from Tobacco to Electronic Cigarettes Does Not Reverse Damage to the Respiratory Epithelium.

The health benefits of switching from tobacco to electronic cigarettes (ECs) are neither confirmed nor well characterized. To address this problem, we used RNA-seq analysis to compare the nasal epithelium transcriptome from the following groups (n = 3 for each group): (1) former smokers who completely switched to second generation ECs for at least 6 months, (2) current tobacco cigarette smokers (CS), and (3) non-smokers (NS). Group three included one former cigarette smoker. The nasal epithelial biopsies from the EC users vs. NS had a higher number of differentially expressed genes (DEGs) than biopsies from the CS vs. NS and CS vs. EC sets (1817 DEGs total for the EC vs. NS, 407 DEGs for the CS vs. NS, and 116 DEGs for the CS vs. EC comparison). In the EC vs. NS comparison, enriched gene ontology terms for the downregulated DEGs included cilium assembly and organization, whereas gene ontologies for upregulated DEGs included immune response, keratinization, and NADPH oxidase. Similarly, ontologies for cilium movement were enriched in the downregulated DEGs for the CS vs. NS group. Reactome pathway analysis gave similar results and also identified keratinization and cornified envelope in the upregulated DEGs in the EC vs. NS comparison. In the CS vs. NS comparison, the enriched Reactome pathways for upregulated DEGs included biological oxidations and several metabolic processes. Regulator effects identified for the EC vs. NS comparison were inflammatory response, cell movement of phagocytes and degranulation of phagocytes. Disease Ontology Sematic Enrichment analysis identified lung disease, mouth disease, periodontal disease and pulmonary fibrosis in the EC vs. NS comparison. Squamous metaplasia associated markers, keratin 10, keratin 13 and involucrin, were increased in the EC vs. NS comparison. Our transcriptomic analysis showed that gene expression profiles associated with EC use are not equivalent to those from non-smokers. EC use may interfere with airway epithelium recovery by promoting increased oxidative stress, inhibition of ciliogenesis, and maintaining an inflammatory response. These transcriptomic alterations may contribute to the progression of diseases with chronic EC use.

Gene regulatory networks shape developmental plasticity of root cell types under water extremes in rice.

Understanding how roots modulate development under varied irrigation or rainfall is crucial for development of climate-resilient crops. We established a toolbox of tagged rice lines to profile translating mRNAs and chromatin accessibility within specific cell populations. We used these to study roots in a range of environments: plates in the lab, controlled greenhouse stress and recovery conditions, and outdoors in a paddy. Integration of chromatin and mRNA data resolves regulatory networks of the following: cycle genes in proliferating cells that attenuate DNA synthesis under submergence; genes involved in auxin signaling, the circadian clock, and small RNA regulation in ground tissue; and suberin biosynthesis, iron transporters, and nitrogen assimilation in endodermal/exodermal cells modulated with water availability. By applying a systems approach, we identify known and candidate driver transcription factors of water-deficit responses and xylem development plasticity. Collectively, this resource will facilitate genetic improvements in root systems for optimal climate resilience.

High-throughput small molecule screening reveals Nrf2-dependent and -independent pathways of cellular stress resistance.

Aging is the dominant risk factor for most chronic diseases. Development of antiaging interventions offers the promise of preventing many such illnesses simultaneously. Cellular stress resistance is an evolutionarily conserved feature of longevity. Here, we identify compounds that induced resistance to the superoxide generator paraquat (PQ), the heavy metal cadmium (Cd), and the DNA alkylator methyl methanesulfonate (MMS). Some rescue compounds conferred resistance to a single stressor, while others provoked multiplex resistance. Induction of stress resistance in fibroblasts was predictive of longevity extension in a published large-scale longevity screen in Caenorhabditis elegans, although not in testing performed in worms and flies with a more restricted set of compounds. Transcriptomic analysis and genetic studies implicated Nrf2/SKN-1 signaling in stress resistance provided by two protective compounds, cardamonin and AEG 3482. Small molecules identified in this work may represent attractive tools to elucidate mechanisms of stress resistance in mammalian cells.

signatureSearch: environment for gene expression signature searching and functional interpretation.

signatureSearch is an R/Bioconductor package that integrates a suite of existing and novel algorithms into an analysis environment for gene expression signature (GES) searching combined with functional enrichment analysis (FEA) and visualization methods to facilitate the interpretation of the search results. In a typical GES search (GESS), a query GES is searched against a database of GESs obtained from large numbers of measurements, such as different genetic backgrounds, disease states and drug perturbations. Database matches sharing correlated signatures with the query indicate related cellular responses frequently governed by connected mechanisms, such as drugs mimicking the expression responses of a disease. To identify which processes are predominantly modulated in the GESS results, we developed specialized FEA methods combined with drug-target network visualization tools. The provided analysis tools are useful for studying the effects of genetic, chemical and environmental perturbations on biological systems, as well as searching single cell GES databases to identify novel network connections or cell types. The signatureSearch software is unique in that it provides access to an integrated environment for GESS/FEA routines that includes several novel search and enrichment methods, efficient data structures, and access to pre-built GES databases, and allowing users to work with custom databases.

HOTAIRM1 lncRNA is downregulated in clear cell renal cell carcinoma and inhibits the hypoxia pathway.

HOXA Transcript Antisense RNA, Myeloid-Specific 1 (HOTAIRM1) is a conserved long non-coding RNA (lncRNA) involved in myeloid and neural differentiation that is deregulated in acute myeloid leukemia and other cancers. Previous studies focused on the nuclear unspliced HOTAIRM1 transcript, however cytoplasmic splice variants exist whose roles have remained unknown. Here, we report novel functions of HOTAIRM1 in the kidney. HOTAIRM1 transcripts are induced during renal lineage differentiation of embryonic stem cells and required for expression of specific renal differentiation genes. We show that the major HOTAIRM1 transcript in differentiated cells is the spliced cytoplasmic HM1-3 isoform and that HM1-3 is downregulated in >90% of clear cell renal cell carcinomas (ccRCCs). Knockdown of HM1-3 in renal cells deregulates hypoxia-responsive and angiogenic genes, including ANGPTL4. Furthermore, HOTAIRM1 transcripts are downregulated by hypoxia-mimetic stress and knockdown of the cytoplasmic HM1-3 isoform in normoxic cells post-transcriptionally induces Hypoxia-Inducible Factor 1α (HIF1α) protein, a key activator of ANGPTL4. Our results demonstrate the pervasive downregulation of the specific HOTAIRM1 cytoplasmic isoform HM1-3 in ccRCC and suggest possible roles of HOTAIRM1 in kidney differentiation and suppression of HIF1-dependent angiogenic pathways.

Genetic Support for Longevity-Enhancing Drug Targets: Issues, Preliminary Data, and Future Directions.

Interventions meant to promote longevity and healthy aging have often been designed or observed to modulate very specific gene or protein targets. If there are naturally occurring genetic variants in such a target that affect longevity as well as the molecular function of that target (eg, the variants influence the expression of the target, acting as "expression quantitative trait loci" or "eQTLs"), this could support a causal relationship between the pharmacologic modulation of the target and longevity and thereby validate the target at some level. We considered the gene targets of many pharmacologic interventions hypothesized to enhance human longevity and explored how many variants there are in those targets that affect gene function (eg, as expression quantitative trait loci). We also determined whether variants in genes associated with longevity-related phenotypes affect gene function or are in linkage disequilibrium with variants that do, and whether pharmacologic studies point to compounds exhibiting activity against those genes. Our results are somewhat ambiguous, suggesting that integrating genetic association study results with functional genomic and pharmacologic studies is necessary to shed light on genetically mediated targets for longevity-enhancing drugs. Such integration will require more sophisticated data sets, phenotypic definitions, and bioinformatics approaches to be useful.

Experimental Acute Exposure to Thirdhand Smoke and Changes in the Human Nasal Epithelial Transcriptome: A Randomized Clinical Trial.

Importance: No previous studies have shown that acute inhalation of thirdhand smoke (THS) activates stress and survival pathways in the human nasal epithelium. Objective: To evaluate gene expression in the nasal epithelium of nonsmoking women following acute inhalation of clean air and THS. Design, Setting, and Participants: Nasal epithelium samples were obtained from participants in a randomized clinical trial (2011-2015) on the health effects of inhaled THS. In a crossover design, participants were exposed, head only, to THS and to conditioned, filtered air in a laboratory setting. The order of exposures was randomized and exposures were separated by at least 21 days. Ribonucleic acid was obtained from a subset of 4 healthy, nonsmoking women. Exposures: By chance, women in the subset were randomized to receive clean air exposure first and THS exposure second. Exposures lasted 3 hours. Main Outcomes and Measures: Differentially expressed genes were identified using RNA sequencing with a false-discovery rate less than 0.1. Results: Participants were 4 healthy, nonsmoking women aged 27 to 49 years (mean [SD] age, 42 [10.2] years) with no chronic diseases. A total of 389 differentially expressed genes were identified in nasal epithelium exposed to THS, while only 2 genes, which were not studied further, were affected by clean air. Enriched gene ontology terms associated with stress-induced mitochondrial hyperfusion were identified, such as respiratory electron transport chain (q = 2.84 × 10-3) and mitochondrial inner membrane (q = 7.21 × 10-6). Reactome pathway analysis identified terms associated with upregulation of DNA repair mechanisms, such as nucleotide excision repair (q = 1.05 × 10-2). Enrichment analyses using ingenuity pathway analysis identified canonical pathways related to stress-induced mitochondrial hyperfusion (eg, increased oxidative phosphorylation) (P = .001), oxidative stress (eg, glutathione depletion phase II reactions) (P = .04), and cell survival (z score = 5.026). Conclusions and Relevance: This study found that acute inhalation of THS caused cell stress that led to the activation of survival pathways. Some responses were consistent with stress-induced mitochondrial hyperfusion and similar to those demonstrated previously in vitro. These data may be valuable to physicians treating patients exposed to THS and may aid in formulating regulations for the remediation of THS-contaminated environments.

Global isoform-specific transcript alterations and deregulated networks in clear cell renal cell carcinoma.

Extensive genome-wide analyses of deregulated gene expression have now been performed for many types of cancer. However, most studies have focused on deregulation at the gene-level, which may overlook the alterations of specific transcripts for a given gene. Clear cell renal cell carcinoma (ccRCC) is one of the best-characterized and most pervasive renal cancers, and ccRCCs are well-documented to have aberrant RNA processing. In the present study, we examine the extent of aberrant isoform-specific RNA expression by reporting a comprehensive transcript-level analysis, using the new kallisto-sleuth-RATs pipeline, investigating coding and non-coding differential transcript expression in ccRCC. We analyzed 50 ccRCC tumors and their matched normal samples from The Cancer Genome Altas datasets. We identified 7,339 differentially expressed transcripts and 94 genes exhibiting differential transcript isoform usage in ccRCC. Additionally, transcript-level coexpression network analyses identified vasculature development and the tricarboxylic acid cycle as the most significantly deregulated networks correlating with ccRCC progression. These analyses uncovered several uncharacterized transcripts, including lncRNAs FGD5-AS1 and AL035661.1, as potential regulators of the tricarboxylic acid cycle associated with ccRCC progression. As ccRCC still presents treatment challenges, our results provide a new resource of potential therapeutics targets and highlight the importance of exploring alternative methodologies in transcriptome-wide studies.

Sequence analysis of the potato aphid Macrosiphum euphorbiae transcriptome identified two new viruses.

The potato aphid, Macrosiphum euphorbiae, is an important agricultural pest that causes economic losses to potato and tomato production. To establish the transcriptome for this aphid, RNA-Seq libraries constructed from aphids maintained on tomato plants were used in Illumina sequencing generating 52.6 million 75-105 bp paired-end reads. The reads were assembled using Velvet/Oases software with SEED preprocessing resulting in 22,137 contigs with an N50 value of 2,003bp. After removal of contigs from tomato host origin, 20,254 contigs were annotated using BLASTx searches against the non-redundant protein database from the National Center for Biotechnology Information (NCBI) as well as IntereProScan. This identified matches for 74% of the potato aphid contigs. The highest ranking hits for over 12,700 contigs were against the related pea aphid, Acyrthosiphon pisum. Gene Ontology (GO) was used to classify the identified M. euphorbiae contigs into biological process, cellular component and molecular function. Among the contigs, sequences of microbial origin were identified. Sixty five contigs were from the aphid bacterial obligate endosymbiont Buchnera aphidicola origin and two contigs had amino acid similarities to viruses. The latter two were named Macrosiphum euphorbiae virus 2 (MeV-2) and Macrosiphum euphorbiae virus 3 (MeV-3). The highest sequence identity to MeV-2 had the Dysaphis plantaginea densovirus, while to MeV-3 is the Hubei sobemo-like virus 49. Characterization of MeV-2 and MeV-3 indicated that both are transmitted vertically from adult aphids to nymphs. MeV-2 peptides were detected in the aphid saliva and only MeV-2 and not MeV-3 nucleic acids were detected inside tomato leaves exposed to virus-infected aphids. However, MeV-2 nucleic acids did not persist in tomato leaf tissues, after clearing the plants from aphids, indicating that MeV-2 is likely an aphid virus.

Transcriptome-wide microRNA and target dynamics in the fat body during the gonadotrophic cycle of Aedes aegypti.

The mosquito Aedes aegypti is a major vector of numerous viral diseases, because it requires a blood meal to facilitate egg development. The fat body, a counterpart of mammalian liver and adipose tissues, is the metabolic center, playing a key role in reproduction. Therefore, understanding of regulatory networks controlling its functions is critical, and the role of microRNAs (miRNAs) in the process is largely unknown. We aimed to explore miRNA expression and potential targets in the female fat body of Ae. aegypti, as well as their changes postblood meal (PBM). Small RNA library analysis revealed five unique miRNA patterns sequentially expressed at five sampled time points, likely responding to, and affecting, waves of upstream hormonal signals and gene expression in the same period. To link miRNA identities with downstream targets, transcriptome-wide mRNA 3' UTR interaction sites were experimentally determined at 72 h posteclosion and 24 h PBM through Argonaute 1 cross-linking and immunoprecipitation followed by high-throughput sequencing. Several target sites were validated by means of in vitro luciferase assays with wild-type and mutated 3' UTRs for six miRNA families. With established transgenic lines, consistent results were observed with spatiotemporal knockdown of miR-8 and luciferase assays. We further investigated miRNAs potentially regulating various physiological processes based on Clusters of Orthologous Groups functional categories. Hence, the present work comprehensively elucidated miRNA expression and target dynamics in the female mosquito fat body, providing a solid foundation for future functional studies of miRNA regulation during the gonadotrophic cycle.

Large-scale bioactivity analysis of the small-molecule assayed proteome.

This study presents an analysis of the small molecule bioactivity profiles across large quantities of diverse protein families represented in PubChem BioAssay. We compared the bioactivity profiles of FDA approved drugs to non-FDA approved compounds, and report several distinct patterns characteristic of the approved drugs. We found that a large fraction of the previously reported higher target promiscuity among FDA approved compounds, compared to non-FDA approved bioactives, was frequently due to cross-reactivity within rather than across protein families. We identified 804 potentially novel protein target candidates for FDA approved drugs, as well as 901 potentially novel target candidates with active non-FDA approved compounds, but no FDA approved drugs with activity against these targets. We also identified 486348 potentially novel compounds active against the same targets as FDA approved drugs, as well as 153402 potentially novel compounds active against targets without active FDA approved drugs. By quantifying the agreement among replicated screens, we estimated that more than half of these novel outcomes are reproducible. Using biclustering, we identified many dense clusters of FDA approved drugs with enriched activity against a common set of protein targets. We also report the distribution of compound promiscuity using a Bayesian statistical model, and report the sensitivity and specificity of two common methods for identifying promiscuous compounds. Aggregator assays exhibited greater accuracy in identifying highly promiscuous compounds, while PAINS substructures were able to identify a much larger set of "middle range" promiscuous compounds. Additionally, we report a large number of promiscuous compounds not identified as aggregators or PAINS. In summary, the results of this study represent a rich reference for selecting novel drug and target protein candidates, as well as for eliminating candidate compounds with unselective activities.

Genomic and functional analysis of the host response to acute simian varicella infection in the lung.

Varicella Zoster Virus (VZV) is the causative agent of varicella and herpes zoster. Although it is well established that VZV is transmitted via the respiratory route, the host-pathogen interactions during acute VZV infection in the lungs remain poorly understood due to limited access to clinical samples. To address these gaps in our knowledge, we leveraged a nonhuman primate model of VZV infection where rhesus macaques are intrabronchially challenged with the closely related Simian Varicella Virus (SVV). Acute infection is characterized by immune infiltration of the lung airways, a significant up-regulation of genes involved in antiviral-immunity, and a down-regulation of genes involved in lung development. This is followed by a decrease in viral loads and increased expression of genes associated with cell cycle and tissue repair. These data provide the first characterization of the host response required to control varicella virus replication in the lung and provide insight into mechanisms by which VZV infection can cause lung injury in an immune competent host.

Acute Simian Varicella Virus Infection Causes Robust and Sustained Changes in Gene Expression in the Sensory Ganglia.

Primary infection with varicella-zoster virus (VZV), a neurotropic alphaherpesvirus, results in varicella. VZV establishes latency in the sensory ganglia and can reactivate later in life to cause herpes zoster. The relationship between VZV and its host during acute infection in the sensory ganglia is not well understood due to limited access to clinical specimens. Intrabronchial inoculation of rhesus macaques with simian varicella virus (SVV) recapitulates the hallmarks of VZV infection in humans. We leveraged this animal model to characterize the host-pathogen interactions in the ganglia during both acute and latent infection by measuring both viral and host transcriptomes on days postinfection (dpi) 3, 7, 10, 14, and 100. SVV DNA and transcripts were detected in sensory ganglia 3 dpi, before the appearance of rash. CD4 and CD8 T cells were also detected in the sensory ganglia 3 dpi. Moreover, lung-resident T cells isolated from the same animals 3 dpi also harbored SVV DNA and transcripts, suggesting that T cells may be responsible for trafficking SVV to the ganglia. Transcriptome sequencing (RNA-Seq) analysis showed that cessation of viral transcription 7 dpi coincides with a robust antiviral innate immune response in the ganglia. Interestingly, a significant number of genes that play a critical role in nervous system development and function remained downregulated into latency. These studies provide novel insights into host-pathogen interactions in the sensory ganglia during acute varicella and demonstrate that SVV infection results in profound and sustained changes in neuronal gene expression. IMPORTANCE: Many aspects of VZV infection of sensory ganglia remain poorly understood, due to limited access to human specimens and the fact that VZV is strictly a human virus. Infection of rhesus macaques with simian varicella virus (SVV), a homolog of VZV, provides a robust model of the human disease. Using this model, we show that SVV reaches the ganglia early after infection, most likely by T cells, and that the induction of a robust innate immune response correlates with cessation of virus transcription. We also report significant changes in the expression of genes that play an important role in neuronal function. Importantly, these changes persist long after viral replication ceases. Given the homology between SVV and VZV, and the genetic and physiological similarities between rhesus macaques and humans, our results provide novel insight into the interactions between VZV and its human host and explain some of the neurological consequences of VZV infection.

systemPipeR: NGS workflow and report generation environment.

BACKGROUND: Next-generation sequencing (NGS) has revolutionized how research is carried out in many areas of biology and medicine. However, the analysis of NGS data remains a major obstacle to the efficient utilization of the technology, as it requires complex multi-step processing of big data demanding considerable computational expertise from users. While substantial effort has been invested on the development of software dedicated to the individual analysis steps of NGS experiments, insufficient resources are currently available for integrating the individual software components within the widely used R/Bioconductor environment into automated workflows capable of running the analysis of most types of NGS applications from start-to-finish in a time-efficient and reproducible manner. RESULTS: To address this need, we have developed the R/Bioconductor package systemPipeR. It is an extensible environment for both building and running end-to-end analysis workflows with automated report generation for a wide range of NGS applications. Its unique features include a uniform workflow interface across different NGS applications, automated report generation, and support for running both R and command-line software on local computers and computer clusters. A flexible sample annotation infrastructure efficiently handles complex sample sets and experimental designs. To simplify the analysis of widely used NGS applications, the package provides pre-configured workflows and reporting templates for RNA-Seq, ChIP-Seq, VAR-Seq and Ribo-Seq. Additional workflow templates will be provided in the future. CONCLUSIONS: systemPipeR accelerates the extraction of reproducible analysis results from NGS experiments. By combining the capabilities of many R/Bioconductor and command-line tools, it makes efficient use of existing software resources without limiting the user to a set of predefined methods or environments. systemPipeR is freely available for all common operating systems from Bioconductor ( http://bioconductor.org/packages/devel/systemPipeR ).

bioassayR: Cross-Target Analysis of Small Molecule Bioactivity.

Despite a large and rapidly growing body of small molecule bioactivity screens available in the public domain, systematic leverage of the data to assess target druggability and compound selectivity has been confounded by a lack of suitable cross-target analysis software. We have developed bioassayR, a computational tool that enables simultaneous analysis of thousands of bioassay experiments performed over a diverse set of compounds and biological targets. Unique features include support for large-scale cross-target analyses of both public and custom bioassays, generation of high throughput screening fingerprints (HTSFPs), and an optional preloaded database that provides access to a substantial portion of publicly available bioactivity data. bioassayR is implemented as an open-source R/Bioconductor package available from https://bioconductor.org/packages/bioassayR/ .

Transcriptomes of Eight Arabidopsis thaliana Accessions Reveal Core Conserved, Genotype- and Organ-Specific Responses to Flooding Stress.

Climate change has increased the frequency and severity of flooding events, with significant negative impact on agricultural productivity. These events often submerge plant aerial organs and roots, limiting growth and survival due to a severe reduction in light reactions and gas exchange necessary for photosynthesis and respiration, respectively. To distinguish molecular responses to the compound stress imposed by submergence, we investigated transcriptomic adjustments to darkness in air and under submerged conditions using eight Arabidopsis (Arabidopsis thaliana) accessions differing significantly in sensitivity to submergence. Evaluation of root and rosette transcriptomes revealed an early transcriptional and posttranscriptional response signature that was conserved primarily across genotypes, although flooding susceptibility-associated and genotype-specific responses also were uncovered. Posttranscriptional regulation encompassed darkness- and submergence-induced alternative splicing of transcripts from pathways involved in the alternative mobilization of energy reserves. The organ-specific transcriptome adjustments reflected the distinct physiological status of roots and shoots. Root-specific transcriptome changes included marked up-regulation of chloroplast-encoded photosynthesis and redox-related genes, whereas those of the rosette were related to the regulation of development and growth processes. We identified a novel set of tolerance genes, recognized mainly by quantitative differences. These included a transcriptome signature of more pronounced gluconeogenesis in tolerant accessions, a response that included stress-induced alternative splicing. This study provides organ-specific molecular resolution of genetic variation in submergence responses involving interactions between darkness and low-oxygen constraints of flooding stress and demonstrates that early transcriptome plasticity, including alternative splicing, is associated with the ability to cope with a compound environmental stress.