Transcriptomic and Proteomic Changes Driving Pulmonary Fibrosis Resolution in Young and Old Mice.

Bleomycin-induced pulmonary fibrosis in mice mimics major hallmarks of idiopathic pulmonary fibrosis. Yet in this model, it spontaneously resolves over time. We studied molecular mechanisms of fibrosis resolution and lung repair, focusing on transcriptional and proteomic signatures and the effect of aging. Old mice showed incomplete and delayed lung function recovery 8 weeks after bleomycin instillation. This shift in structural and functional repair in old bleomycin-treated mice was reflected in a temporal shift in gene and protein expression. We reveal gene signatures and signaling pathways that underpin the lung repair process. Importantly, the downregulation of WNT, BMP, and TGFß antagonists Frzb, Sfrp1, Dkk2, Grem1, Fst, Fstl1, and Inhba correlated with lung function improvement. Those genes constitute a network with functions in stem cell pathways, wound, and pulmonary healing. We suggest that insufficient and delayed downregulation of those antagonists during fibrosis resolution in old mice explains the impaired regenerative outcome. Together, we identified signaling pathway molecules with relevance to lung regeneration that should be tested in-depth experimentally as potential therapeutic targets for pulmonary fibrosis.

Induced Systemic Resistance in the Bacillus spp.-Capsicum chinense Jacq.-PepGMV Interaction, Elicited by Defense-Related Gene Expression.

Induced systemic resistance (ISR) is a mechanism involved in the plant defense response against pathogens. Certain members of the Bacillus genus are able to promote the ISR by maintaining a healthy photosynthetic apparatus, which prepares the plant for future stress situations. The goal of the present study was to analyze the effect of the inoculation of Bacillus on the expression of genes involved in plant responses to pathogens, as a part of the ISR, during the interaction of Capsicum chinense infected with PepGMV. The effects of the inoculation of the Bacillus strains in pepper plants infected with PepGMV were evaluated by observing the accumulation of viral DNA and the visible symptoms of pepper plants during a time-course experiment in greenhouse and in in vitro experiments. The relative expression of the defense genes CcNPR1, CcPR10, and CcCOI1 were also evaluated. The results showed that the plants inoculated with Bacillus subtilis K47, Bacillus cereus K46, and Bacillus sp. M9 had a reduction in the PepGMV viral titer, and the symptoms in these plants were less severe compared to the plants infected with PepGMV and non-inoculated with Bacillus. Additionally, an increase in the transcript levels of CcNPR1, CcPR10, and CcCOI1 was observed in plants inoculated with Bacillus strains. Our results suggest that the inoculation of Bacillus strains interferes with the viral replication, through the increase in the transcription of pathogenesis-related genes, which is reflected in a lowered plant symptomatology and an improved yield in the greenhouse, regardless of PepGMV infection status.

Antifibrotic Drug Nintedanib Inhibits CSF1R to Promote IL-4-associated Tissue Repair Macrophages.

Profibrotic and prohomeostatic macrophage phenotypes remain ill-defined, both in vivo and in vitro, impeding the successful development of drugs that reprogram macrophages as an attractive therapeutic approach to manage fibrotic disease. The goal of this study was to reveal profibrotic and prohomeostatic macrophage phenotypes that could guide the design of new therapeutic approaches targeting macrophages to treat fibrotic disease. This study used nintedanib, a broad kinase inhibitor approved for idiopathic pulmonary fibrosis, to dissect lung macrophage phenotypes during fibrosis-linked inflammation by combining in vivo and in vitro bulk and single-cell RNA-sequencing approaches. In the bleomycin model, nintedanib drove the expression of IL-4/IL-13-associated genes important for tissue regeneration and repair at early and late time points in lung macrophages. These findings were replicated in vitro in mouse primary bone marrow-derived macrophages exposed to IL-4/IL-13 and nintedanib. In addition, nintedanib promoted the expression of IL-4/IL-13 pathway genes in human macrophages in vitro. The molecular mechanism was connected to inhibition of the colony stimulating factor 1 (CSF1) receptor in both human and mouse macrophages. Moreover, nintedanib counterbalanced the effects of TNF on IL-4/IL-13 in macrophages to promote expression of IL-4/IL-13-regulated tissue repair genes in fibrotic contexts in vivo and in vitro. This study demonstrates that one of nintedanib's antifibrotic mechanisms is to increase IL-4 signaling in macrophages through inhibition of the CSF1 receptor, resulting in the promotion of tissue repair phenotypes.

miR449 Protects Airway Regeneration by Controlling AURKA/HDAC6-Mediated Ciliary Disassembly.

Airway mucociliary regeneration and function are key players for airway defense and are impaired in chronic obstructive pulmonary disease (COPD). Using transcriptome analysis in COPD-derived bronchial biopsies, we observed a positive correlation between cilia-related genes and microRNA-449 (miR449). In vitro, miR449 was strongly increased during airway epithelial mucociliary differentiation. In vivo, miR449 was upregulated during recovery from chemical or infective insults. miR0449-/- mice (both alleles are deleted) showed impaired ciliated epithelial regeneration after naphthalene and Haemophilus influenzae exposure, accompanied by more intense inflammation and emphysematous manifestations of COPD. The latter occurred spontaneously in aged miR449-/- mice. We identified Aurora kinase A and its effector target HDAC6 as key mediators in miR449-regulated ciliary homeostasis and epithelial regeneration. Aurora kinase A is downregulated upon miR449 overexpression in vitro and upregulated in miR449-/- mouse lungs. Accordingly, imaging studies showed profoundly altered cilia length and morphology accompanied by reduced mucociliary clearance. Pharmacological inhibition of HDAC6 rescued cilia length and coverage in miR449-/- cells, consistent with its tubulin-deacetylating function. Altogether, our study establishes a link between miR449, ciliary dysfunction, and COPD pathogenesis.

CD206+ tumor-associated macrophages cross-present tumor antigen and drive antitumor immunity.

In many solid cancers, tumor-associated macrophages (TAM) represent the predominant myeloid cell population. Antigen (Ag) cross-presentation leading to tumor Ag-directed cytotoxic CD8+ T cell responses is crucial for antitumor immunity. However, the role of recruited monocyte-derived macrophages, including TAM, as potential cross-presenting cells is not well understood. Here, we show that primary human as well as mouse CD206+ macrophages are effective in functional cross-presentation of soluble self-Ag and non-self-Ag, including tumor-associated Ag (TAA), as well as viral Ag. To confirm the presence of cross-presenting TAM in vivo, we performed phenotypic and functional analysis of TAM from B16-F10 and CT26 syngeneic tumor models and have identified CD11b+F4/80hiCD206+ TAM to effectively cross-present TAA. We show that CD11b+CD206+ TAM represent the dominant tumor-infiltrating myeloid cell population, expressing a unique cell surface repertoire, promoting Ag cross-presentation and Ag-specific CD8+ T cell activation comparable with cross-presenting CLEC9A+ DCs (cDC1). The presence of cross-presenting CD206+ TAM is associated with reduced tumor burden in mouse syngeneic tumor models and with improved overall survival in cutaneous melanoma patients. Therefore, the demonstration of effective Ag cross-presentation capabilities of CD206+ TAM, including their clinical relevance, expands our understanding of TAM phenotypic diversity and functional versatility.

Automation enables high-throughput and reproducible single-cell transcriptomics library preparation.

Next-generation sequencing (NGS) has revolutionized genomics, decreasing sequencing costs and allowing researchers to draw correlations between diseases and DNA or RNA changes. Technical advances have enabled the analysis of RNA expression changes between single cells within a heterogeneous population, known as single-cell RNA-seq (scRNA-seq). Despite resolving transcriptomes of cellular subpopulations, scRNA-seq has not replaced RNA-seq, due to higher costs and longer hands-on time. Here, we developed an automated workflow to increase throughput (up to 48 reactions) and to reduce by 75% the hands-on time of scRNA-seq library preparation, using the 10X Genomics Single Cell 3' kit. After gel bead-in-emulsion (GEM) generation on the 10X Genomics Chromium Controller, cDNA amplification was performed, and the product was normalized and subjected to either the manual, standard library preparation method or a fully automated, walk-away method using a Biomek i7 Hybrid liquid handler. Control metrics showed that both quantity and quality of the single-cell gene expression libraries generated were equivalent in size and yield. Key scRNA-seq downstream quality metrics, such as unique molecular identifiers count, mitochondrial RNA content, and cell and gene counts, further showed high correlations between automated and manual workflows. Using the UMAP dimensionality reduction technique to visualize all cells, we were able to further correlate the results observed between the manual and automated methods (R=0.971). The method developed here allows for the fast, error-free, and reproducible multiplex generation of high-quality single-cell gene expression libraries.

pyGenomeTracks: reproducible plots for multivariate genomic datasets.

MOTIVATION: Generating publication ready plots to display multiple genomic tracks can pose a serious challenge. Making desirable and accurate figures requires considerable effort. This is usually done by hand or using a vector graphic software. RESULTS: pyGenomeTracks (PGT) is a modular plotting tool that easily combines multiple tracks. It enables a reproducible and standardized generation of highly customizable and publication ready images. AVAILABILITY AND IMPLEMENTATION: PGT is available through a graphical interface on https://usegalaxy.eu and through the command line. It is provided on conda via the bioconda channel, on pip and it is openly developed on github: https://github.com/deeptools/pyGenomeTracks. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A single-cell RNA-sequencing training and analysis suite using the Galaxy framework.

BACKGROUND: The vast ecosystem of single-cell RNA-sequencing tools has until recently been plagued by an excess of diverging analysis strategies, inconsistent file formats, and compatibility issues between different software suites. The uptake of 10x Genomics datasets has begun to calm this diversity, and the bioinformatics community leans once more towards the large computing requirements and the statistically driven methods needed to process and understand these ever-growing datasets. RESULTS: Here we outline several Galaxy workflows and learning resources for single-cell RNA-sequencing, with the aim of providing a comprehensive analysis environment paired with a thorough user learning experience that bridges the knowledge gap between the computational methods and the underlying cell biology. The Galaxy reproducible bioinformatics framework provides tools, workflows, and trainings that not only enable users to perform 1-click 10x preprocessing but also empower them to demultiplex raw sequencing from custom tagged and full-length sequencing protocols. The downstream analysis supports a range of high-quality interoperable suites separated into common stages of analysis: inspection, filtering, normalization, confounder removal, and clustering. The teaching resources cover concepts from computer science to cell biology. Access to all resources is provided at the singlecell.usegalaxy.eu portal. CONCLUSIONS: The reproducible and training-oriented Galaxy framework provides a sustainable high-performance computing environment for users to run flexible analyses on both 10x and alternative platforms. The tutorials from the Galaxy Training Network along with the frequent training workshops hosted by the Galaxy community provide a means for users to learn, publish, and teach single-cell RNA-sequencing analysis.

Dataset on application of electrochemical and photochemical processes for sulfacetamide antibiotic elimination in water.

Sulfonamide-class antibiotics are recognized as water pollutants, which have negative environmental impacts. A strategy to deal with sulfonamides is throughout the application of oxidation processes. This work presents the treatment of the sulfacetamide (SAM) antibiotic by electrochemical oxidation, UV-C/H2O2 and photo-Fenton process. It was established the main degradation routes during each process action. A DFT computational analysis for SAM structure was done and mass spectra of primary transformation products were determined. Chemical oxygen demand (COD), total organic carbon (TOC) and biochemical oxygen demand (BOD5) were also followed. Additionally, SAM treatment in simulated seawater and hospital wastewater was measured. These data can be useful for comparative purposes about degradation of sulfonamide-class antibiotics by electrochemical and advanced oxidation processes.

Hi-C guided assemblies reveal conserved regulatory topologies on X and autosomes despite extensive genome shuffling.

Genome rearrangements that occur during evolution impose major challenges on regulatory mechanisms that rely on three-dimensional genome architecture. Here, we developed a scaffolding algorithm and generated chromosome-length assemblies from Hi-C data for studying genome topology in three distantly related Drosophila species. We observe extensive genome shuffling between these species with one synteny breakpoint after approximately every six genes. A/B compartments, a set of large gene-dense topologically associating domains (TADs), and spatial contacts between high-affinity sites (HAS) located on the X chromosome are maintained over 40 million years, indicating architectural conservation at various hierarchies. Evolutionary conserved genes cluster in the vicinity of HAS, while HAS locations appear evolutionarily flexible, thus uncoupling functional requirement of dosage compensation from individual positions on the linear X chromosome. Therefore, 3D architecture is preserved even in scenarios of thousands of rearrangements highlighting its relevance for essential processes such as dosage compensation of the X chromosome.

CAPRI enables comparison of evolutionarily conserved RNA interacting regions.

RNA-protein complexes play essential regulatory roles at nearly all levels of gene expression. Using in vivo crosslinking and RNA capture, we report a comprehensive RNA-protein interactome in a metazoan at four levels of resolution: single amino acids, domains, proteins and multisubunit complexes. We devise CAPRI, a method to map RNA-binding domains (RBDs) by simultaneous identification of RNA interacting crosslinked peptides and peptides adjacent to such crosslinked sites. CAPRI identifies more than 3000 RNA proximal peptides in Drosophila and human proteins with more than 45% of them forming new interaction interfaces. The comparison of orthologous proteins enables the identification of evolutionary conserved RBDs in globular domains and intrinsically disordered regions (IDRs). By comparing the sequences of IDRs through evolution, we classify them based on the type of motif, accumulation of tandem repeats, conservation of amino acid composition and high sequence divergence.

Monoclonal Antibodies Against the Capsular Polysaccharides A, C, Y, W, and X of Neisseria meningitidis: A Platform for the Quality Control of Meningococcal Vaccines.

Vaccination has reduced morbidity and mortality of many diseases that previously caused devastating epidemics and deaths globally. Vaccines as a biological product may contain microorganisms or their derivatives. This aspect together with the fact that they are administered to healthy individuals (mainly children) means that approximately 70% of vaccines development time is dedicated to quality control. Monoclonal antibodies (MAbs) have become essential analytical tools for application in ELISAs, Western and Dot blotting, immunoprecipitation, and flow cytometric assays that ensure the quality control of vaccines. The aim of this work is to present a review of the methods used to obtain a platform of MAbs against Neisseria meningitidis polysaccharide antigens to use as an analytical tool for quality control of anti-meningococcal polysaccharide (Ps) vaccines. The MAbs obtained are used in five sandwich ELISAs developed for Ps quantification. The assays showed good reproducibility and repeatability, with quantitation and detection limits below 1 ng/mL. Dot Blot, as the Identity test of the Ps vaccine, was carried out to positively identify licensed and experimental vaccines. All assays described are suitable for the screening of multiple vaccine samples and could be useful for monitoring lot-to-lot consistency and stability.

RELACS nuclei barcoding enables high-throughput ChIP-seq.

Chromatin immunoprecipitation followed by deep sequencing (ChIP-seq) is an invaluable tool for mapping chromatin-associated proteins. Current barcoding strategies aim to improve assay throughput and scalability but intense sample handling and lack of standardization over cell types, cell numbers and epitopes hinder wide-spread use in the field. Here, we present a barcoding method to enable high-throughput ChIP-seq using common molecular biology techniques. The method, called RELACS (restriction enzyme-based labeling of chromatin in situ) relies on standardized nuclei extraction from any source and employs chromatin cutting and barcoding within intact nuclei. Barcoded nuclei are pooled and processed within the same ChIP reaction, for maximal comparability and workload reduction. The innovative barcoding concept is particularly user-friendly and suitable for implementation to standardized large-scale clinical studies and scarce samples. Aiming to maximize universality and scalability, RELACS can generate ChIP-seq libraries for transcription factors and histone modifications from hundreds of samples within three days.

A comprehensive analysis of 195 DNA methylomes reveals shared and cell-specific features of partially methylated domains.

BACKGROUND: Partially methylated domains are extended regions in the genome exhibiting a reduced average DNA methylation level. They cover gene-poor and transcriptionally inactive regions and tend to be heterochromatic. We present a comprehensive comparative analysis of partially methylated domains in human and mouse cells, to identify structural and functional features associated with them. RESULTS: Partially methylated domains are present in up to 75% of the genome in human and mouse cells irrespective of their tissue or cell origin. Each cell type has a distinct set of partially methylated domains, and genes expressed in such domains show a strong cell type effect. The methylation level varies between cell types with a more pronounced effect in differentiating and replicating cells. The lowest level of methylation is observed in highly proliferating and immortal cancer cell lines. A decrease of DNA methylation within partially methylated domains tends to be linked to an increase in heterochromatic histone marks and a decrease of gene expression. Characteristic combinations of heterochromatic signatures in partially methylated domains are linked to domains of early and middle S-phase and late S-G2 phases of DNA replication. CONCLUSIONS: Partially methylated domains are prominent signatures of long-range epigenomic organization. Integrative analysis identifies them as important general, lineage- and cell type-specific topological features. Changes in partially methylated domains are hallmarks of cell differentiation, with decreased methylation levels and increased heterochromatic marks being linked to enhanced cell proliferation. In combination with broad histone marks, partially methylated domains demarcate distinct domains of late DNA replication.

GOATOOLS: A Python library for Gene Ontology analyses.

The biological interpretation of gene lists with interesting shared properties, such as up- or down-regulation in a particular experiment, is typically accomplished using gene ontology enrichment analysis tools. Given a list of genes, a gene ontology (GO) enrichment analysis may return hundreds of statistically significant GO results in a "flat" list, which can be challenging to summarize. It can also be difficult to keep pace with rapidly expanding biological knowledge, which often results in daily changes to any of the over 47,000 gene ontologies that describe biological knowledge. GOATOOLS, a Python-based library, makes it more efficient to stay current with the latest ontologies and annotations, perform gene ontology enrichment analyses to determine over- and under-represented terms, and organize results for greater clarity and easier interpretation using a novel GOATOOLS GO grouping method. We performed functional analyses on both stochastic simulation data and real data from a published RNA-seq study to compare the enrichment results from GOATOOLS to two other popular tools: DAVID and GOstats. GOATOOLS is freely available through GitHub: https://github.com/tanghaibao/goatools .

Community-Driven Data Analysis Training for Biology.

The primary problem with the explosion of biomedical datasets is not the data, not computational resources, and not the required storage space, but the general lack of trained and skilled researchers to manipulate and analyze these data. Eliminating this problem requires development of comprehensive educational resources. Here we present a community-driven framework that enables modern, interactive teaching of data analytics in life sciences and facilitates the development of training materials. The key feature of our system is that it is not a static but a continuously improved collection of tutorials. By coupling tutorials with a web-based analysis framework, biomedical researchers can learn by performing computation themselves through a web browser without the need to install software or search for example datasets. Our ultimate goal is to expand the breadth of training materials to include fundamental statistical and data science topics and to precipitate a complete re-engineering of undergraduate and graduate curricula in life sciences. This project is accessible at https://training.galaxyproject.org.

Galaxy HiCExplorer: a web server for reproducible Hi-C data analysis, quality control and visualization.

Galaxy HiCExplorer is a web server that facilitates the study of the 3D conformation of chromatin by allowing Hi-C data processing, analysis and visualization. With the Galaxy HiCExplorer web server, users with little bioinformatic background can perform every step of the analysis in one workflow: mapping of the raw sequence data, creation of Hi-C contact matrices, quality assessment, correction of contact matrices and identification of topological associated domains (TADs) and A/B compartments. Users can create publication ready plots of the contact matrix, A/B compartments, and TADs on a selected genomic locus, along with additional information like gene tracks or ChIP-seq signals. Galaxy HiCExplorer is freely usable at: https://hicexplorer.usegalaxy.eu and is available as a Docker container: https://github.com/deeptools/docker-galaxy-hicexplorer.

High-resolution TADs reveal DNA sequences underlying genome organization in flies.

Despite an abundance of new studies about topologically associating domains (TADs), the role of genetic information in TAD formation is still not fully understood. Here we use our software, HiCExplorer (hicexplorer.readthedocs.io) to annotate >2800 high-resolution (570 bp) TAD boundaries in Drosophila melanogaster. We identify eight DNA motifs enriched at boundaries, including a motif bound by the M1BP protein, and two new boundary motifs. In contrast to mammals, the CTCF motif is only enriched on a small fraction of boundaries flanking inactive chromatin while most active boundaries contain the motifs bound by the M1BP or Beaf-32 proteins. We demonstrate that boundaries can be accurately predicted using only the motif sequences at open chromatin sites. We propose that DNA sequence guides the genome architecture by allocation of boundary proteins in the genome. Finally, we present an interactive online database to access and explore the spatial organization of fly, mouse and human genomes, available at http://chorogenome.ie-freiburg.mpg.de .

Functional interplay between MSL1 and CDK7 controls RNA polymerase II Ser5 phosphorylation.

Proper gene expression requires coordinated interplay among transcriptional coactivators, transcription factors and the general transcription machinery. We report here that MSL1, a central component of the dosage compensation complex in Drosophila melanogaster and Drosophila virilis, displays evolutionarily conserved sex-independent binding to promoters. Genetic and biochemical analyses reveal a functional interaction of MSL1 with CDK7, a subunit of the Cdk-activating kinase (CAK) complex of the general transcription factor TFIIH. Importantly, MSL1 depletion leads to decreased phosphorylation of Ser5 of RNA polymerase II. In addition, we demonstrate that MSL1 is a phosphoprotein, and transgenic flies expressing MSL1 phosphomutants show mislocalization of the histone acetyltransferase MOF and histone H4 K16 acetylation, thus ultimately causing male lethality due to a failure of dosage compensation. We propose that, by virtue of its interaction with components of the general transcription machinery, MSL1 exists in different phosphorylation states, thereby modulating transcription in flies.

deepTools2: a next generation web server for deep-sequencing data analysis.

We present an update to our Galaxy-based web server for processing and visualizing deeply sequenced data. Its core tool set, deepTools, allows users to perform complete bioinformatic workflows ranging from quality controls and normalizations of aligned reads to integrative analyses, including clustering and visualization approaches. Since we first described our deepTools Galaxy server in 2014, we have implemented new solutions for many requests from the community and our users. Here, we introduce significant enhancements and new tools to further improve data visualization and interpretation. deepTools continue to be open to all users and freely available as a web service at deeptools.ie-freiburg.mpg.de The new deepTools2 suite can be easily deployed within any Galaxy framework via the toolshed repository, and we also provide source code for command line usage under Linux and Mac OS X. A public and documented API for access to deepTools functionality is also available.