Innate immune responses vary by pathogen and host genetics. We analyze quantitative trait loci (eQTLs) and transcriptomes of monocytes from 215 individuals stimulated by fungal, Gram-negative or Gram-positive bacterial pathogens. We identify conserved monocyte responses to bacterial pathogens and a distinct antifungal response. These include 745 response eQTLs (reQTLs) and corresponding genes with pathogen-specific effects, which we find first in samples of male donors and subsequently confirm for selected reQTLs in females. reQTLs affect predominantly upregulated genes that regulate immune response via e.g., NOD-like, C-type lectin, Toll-like and complement receptor-signaling pathways. Hence, reQTLs provide a functional explanation for individual differences in innate response patterns. Our identified reQTLs are also associated with cancer, autoimmunity, inflammatory and infectious diseases as shown by external genome-wide association studies. Thus, reQTLs help to explain interindividual variation in immune response to infection and provide candidate genes for variants associated with a range of diseases.

Aspergillus fumigatus is a major fungal pathogen that causes severe problems due to its increasing resistance to many therapeutic agents. Fludioxonil is a compound that triggers a lethal activation of the fungal-specific High Osmolarity Glycerol pathway. Its pronounced antifungal activity against A. fumigatus and other pathogenic molds renders this agent an attractive lead substance for the development of new therapeutics. The group III hydride histidine kinase TcsC and its downstream target Skn7 are key elements of the multistep phosphorelay that represents the initial section of the High Osmolarity Glycerol pathway. Loss of tcsC results in resistance to fludioxonil, whereas a Δskn7 mutant is partially, but not completely resistant. In this study, we compared the fludioxonil-induced transcriptional responses in the ΔtcsC and Δskn7 mutant and their parental A. fumigatus strain. The number of differentially expressed genes correlates well with the susceptibility level of the individual strains. The wild type and, to a lesser extend also the Δskn7 mutant, showed a multi-faceted stress response involving genes linked to ribosomal and peroxisomal function, iron homeostasis and oxidative stress. A marked difference between the sensitive wild type and the largely resistant Δskn7 mutant was evident for many cell wall-related genes and in particular those involved in the biosynthesis of chitin. Biochemical data corroborate this differential gene expression that does not occur in response to hyperosmotic stress. Our data reveal that fludioxonil induces a strong and TcsC-dependent stress that affects many aspects of the cellular machinery. The data also demonstrate a link between Skn7 and the cell wall reorganizations that foster the characteristic ballooning and the subsequent lysis of fludioxonil-treated cells.

Candida species are a major cause of invasive fungal infections. While Candida albicans, C. glabrata, C. parapsilosis, and C. tropicalis are the most dominant species causing life-threatening candidiasis, C. auris recently emerged as a new species causing invasive infections with high rates of clinical treatment failures. To mimic initial phases of systemic Candida infections with dissemination via the bloodstream and to elucidate the pathogenic potential of C. auris, we used an ex vivo whole blood infection model. Similar to other clinically relevant Candida spp., C. auris is efficiently killed in human blood, but showed characteristic patterns of immune cell association, survival rates, and cytokine induction. Dual-species transcriptional profiling of C. auris-infected blood revealed a unique C. auris gene expression program during infection, while the host response proofed similar and conserved compared to other Candida species. C. auris-specific responses included adaptation and survival strategies, such as counteracting oxidative burst of immune cells, but also expression of potential virulence factors, (drug) transporters, and cell surface-associated genes. Despite comparable pathogenicity to other Candida species in our model, C. auris-specific transcriptional adaptations as well as its increased stress resistance and long-term environmental survival, likely contribute to the high risk of contamination and distribution in a nosocomial setting. Moreover, infections of neutrophils with pre-starved C. auris cells suggest that environmental preconditioning can have modulatory effects on the early host interaction. In summary, we present novel insights into C. auris pathogenicity, revealing adaptations to human blood and environmental niches distinctive from other Candida species.

The human pathogenic fungus Aspergillus fumigatus is readily eradicated by the innate immunity of immunocompetent human hosts, but can cause severe infections, such as invasive aspergillosis (IA), in immunocompromised individuals. During infection, the fungal redox homeostasis can be challenged by reactive oxygen species (ROS), either derived from the oxidative burst of innate immune cells or the action of antifungal drugs. The peroxiredoxin Asp f3 was found to be essential to cause IA in mice, but how Asp f3 integrates with fungal redox homeostasis remains unknown. Here, we show that in vivo, Asp f3 acts as a sensor for ROS. While global transcription in fungal hyphae under minimal growth conditions was fully independent of Asp f3, a robust induction of the oxidative stress response required the presence of the peroxiredoxin. Hyphae devoid of Asp f3 failed to activate several redox active genes, like members of the gliotoxin biosynthesis gene cluster and integral members of the Afyap1 regulon, the central activator of the ROS defense machinery in fungi. Upon deletion of the asp f3 gene Afyap1 displayed significantly reduced nuclear localization during ROS exposure, indicating that Asp f3 can act as an intracellular redox sensor for several target proteins.

Rhinovirus (RV) and influenza virus are the most frequently detected respiratory viruses among adult patients with community acquired pneumonia. Previous clinical studies have identified major differences in the clinical presentations and inflammatory or immune response during these infections. A systematic transcriptomic analysis directly comparing influenza and RV is lacking. Here, we sought to compare the transcriptomic response to these viral infections. Human airway epithelial Calu-3 cells were infected with contemporary clinical isolates of RV, influenza A virus (IAV), or influenza B virus (IBV). Host gene expression was determined using RNA-seq. Differentially expressed genes (DEGs) with respect to mock-infected cells were identified using the overlapping gene-set of four different statistical models. Transcriptomic analysis showed that RV-infected cells have a more blunted host response with fewer DEGs than IAV or IBV-infected cells. IFNL1 and CXCL10 were among the most upregulated DEGs during RV, IAV, and IBV infection. Other DEGs that were highly expressed for all 3 viruses were mainly genes related to type I or type III interferons (RSAD2, IDO1) and chemokines (CXCL11). Notably, ICAM5, a known receptor for enterovirus D68, was highly expressed during RV infection only. Gene Set Enrichment Analysis (GSEA) confirmed that pathways associated with interferon response, innate immunity, or regulation of inflammatory response, were most perturbed for all three viruses. Network analysis showed that steroid-related pathways were enriched. Taken together, our data using contemporary virus strains suggests that genes related to interferon and chemokine predominated the host response associated with RV, IAV, and IBV infection. Several highly expressed genes, especially ICAM5 which is preferentially-induced during RV infection, deserve further investigation.

High-throughput RNA sequencing (RNA-seq) is routinely applied to study diverse biological processes; however, when performed separately on interacting organisms, systemic noise intrinsic to RNA extraction, library preparation, and sequencing hampers the identification of cross-species interaction nodes. Here, we develop triple RNA-seq to simultaneously detect transcriptomes of monocyte-derived dendritic cells (moDCs) infected with the frequently co-occurring pulmonary pathogens Aspergillus fumigatus and human cytomegalovirus (CMV). Comparing expression patterns after co-infection with those after single infections, our data reveal synergistic effects and mutual interferences between host responses to the two pathogens. For example, CMV attenuates the fungus-mediated activation of pro-inflammatory cytokines through NF-κB (nuclear factor κB) and NFAT (nuclear factor of activated T cells) cascades, while A. fumigatus impairs viral clearance by counteracting viral nucleic acid-induced activation of type I interferon signaling. Together, the analytical power of triple RNA-seq proposes molecular hubs in the differential moDC response to fungal/viral single infection or co-infection that contribute to our understanding of the etiology and, potentially, clearance of post-transplant infections.

Only four species, Candida albicans, C. glabrata, C. parapsilosis, and C. tropicalis, together account for about 90% of all Candida bloodstream infections and are among the most common causes of invasive fungal infections of humans. However, virulence potential varies among these species, and the phylogenetic tree reveals that their pathogenicity may have emerged several times independently during evolution. We therefore tested these four species in a human whole-blood infection model to determine, via comprehensive dual-species RNA-sequencing analyses, which fungal infection strategies are conserved and which are recent evolutionary developments. The ex vivo infection progressed from initial immune cell interactions to nearly complete killing of all fungal cells. During the course of infection, we characterized important parameters of pathogen-host interactions, such as fungal survival, types of interacting immune cells, and cytokine release. On the transcriptional level, we obtained a predominantly uniform and species-independent human response governed by a strong upregulation of proinflammatory processes, which was downregulated at later time points after most of the fungal cells were killed. In stark contrast, we observed that the different fungal species pursued predominantly individual strategies and showed significantly different global transcriptome patterns. Among other findings, our functional analyses revealed that the fungal species relied on different metabolic pathways and virulence factors to survive the host-imposed stress. These data show that adaptation of Candida species as a response to the host is not a phylogenetic trait, but rather has likely evolved independently as a prerequisite to cause human infections.IMPORTANCE To ensure their survival, pathogens have to adapt immediately to new environments in their hosts, for example, during the transition from the gut to the bloodstream. Here, we investigated the basis of this adaptation in a group of fungal species which are among the most common causes of hospital-acquired infections, the Candida species. On the basis of a human whole-blood infection model, we studied which genes and processes are active over the course of an infection in both the host and four different Candida pathogens. Remarkably, we found that, while the human host response during the early phase of infection is predominantly uniform, the pathogens pursue largely individual strategies and each one regulates genes involved in largely disparate processes in the blood. Our results reveal that C. albicans, C. glabrata, C. parapsilosis, and C. tropicalis all have developed individual strategies for survival in the host. This indicates that their pathogenicity in humans has evolved several times independently and that genes which are central for survival in the host for one species may be irrelevant in another.

Predatory interactions among microbes are considered to be a major evolutionary driving force for biodiversity and the defense against phagocytic killing. The fungivorous amoeba Protostelium aurantium has a wide fungal food spectrum but strongly discriminates among major pathogenic members of the Saccharomycotina. While C. albicans is not recognized, C. glabrata is rapidly internalized, but remains undigested. Phagocytic killing and feeding by P. aurantium is highly effective for the third major fungal pathogen, C. parapsilosis. Here we show that the different prey patterns of the three yeasts were reflected by distinct transcriptional responses, indicating fungal copper and redox homeostasis as primary targets during intracellular killing of C. parapsilosis. Gene deletions in this fungus for the highly expressed copper exporter Crp1 and the peroxiredoxin Prx1 confirmed their role in copper and redox homeostasis, respectively and identified methionine biosynthesis as a ROS sensitive metabolic target during predation. Both, intact Cu export and redox homeostasis contributed to the survival of C. parapsilosis not only when encountering P. aurantium, but also in the presence of human macrophages. As both genes were found to be widely conserved within the entire Candida clade, our results suggest that they could be part of a basic tool-kit to survive phagocytic attacks by environmental predators.

ABSTRACTIn transcriptomics, the study of the total set of RNAs transcribed by the cell, RNA sequencing (RNA-seq) has become the standard tool for analysing gene expression. The primary goal is the detection of genes whose expression changes significantly between two or more conditions, either for a single species or for two or more interacting species at the same time (dual RNA-seq, triple RNA-seq and so forth). The analysis of RNA-seq can be simplified as many steps of the data pre-processing can be standardised in a pipeline.In this publication we present the “GEO2RNAseq” pipeline for complete, quick and concurrent pre-processing of single, dual, and triple RNA-seq data. It covers all pre-processing steps starting from raw sequencing data to the analysis of differentially expressed genes, including various tables and figures to report intermediate and final results. Raw data may be provided in FASTQ format or can be downloaded automatically from the Gene Expression Omnibus repository. GEO2RNAseq strongly incorporates experimental as well as computational metadata. GEO2RNAseq is implemented in R, lightweight, easy to install via Conda and easy to use, but still very flexible through using modular programming and offering many extensions and alternative workflows.GEO2RNAseq is publicly available at https://anaconda.org/xentrics/r-geo2rnaseq and https://bitbucket.org/thomas_wolf/geo2rnaseq/overview, including source code, installation instruction, and comprehensive package documentation.

SummaryPredatory interactions among microbes are considered to be a major evolutionary driving force for biodiversity and the defense against phagocytic killing. The fungivorous amoebaProtostelium aurantiumhas a wide fungal food spectrum but strongly discriminates among major pathogenic members of theSaccharomycotina. WhileC. albicansis not recognized,C. glabratais rapidly internalized, but remains undigested. Phagocytic killing and feeding byP. aurantiumis highly effective for the third major fungal pathogen,C. parapsilosis.Here we show that the different prey patterns of the three yeasts were reflected by distinct transcriptional responses, indicating fungal copper and redox homeostasis as primary targets during intracellular killing ofC. parapsilosis. Gene deletions in this fungus for the highly expressed copper exporter Crp1 and the peroxiredoxin Prx1 confirmed their role in copper and redox homeostasis, respectively and identified methionine biosynthesis as a ROS sensitive metabolic target during predation. Both, intact Cu export and redox homeostasis contributed to the survival ofC. parapsilosisnot only when encounteringP. aurantium, but also in the presence of human macrophages. As both genes were found to be widely conserved within the entireCandidaclade, our results suggest that they could be part of a basic tool-kit to survive phagocytic attacks by environmental predators.

Organisms do not exist isolated from each other, but constantly interact. Cells can sense the presence of interaction partners by a range of receptors and, via complex regulatory networks, specifically react by changing the expression of many of their genes. Technological advances in next-generation sequencing over the recent years now allow us to apply RNA sequencing to two species at the same time (dual RNA-seq), and thus to directly study the gene expression of two interacting species without the need to physically separate cells or RNA. In this review, we give an overview over the latest studies in interspecies interactions made possible by dual RNA-seq, ranging from pathogenic to symbiotic relationships. We summarize state-of-the-art experimental techniques, bioinformatic data analysis and data interpretation, while also highlighting potential problems and pitfalls starting from the selection of meaningful time points and number of reads to matters of rRNA depletion. A short outlook on new trends in the field of dual RNA-seq concludes this review, looking at sequencing of non-coding RNAs during host-pathogen interactions and the prediction of molecular interspecies interactions networks.

Many antibiotics, chemotherapeutics, crop protection agents and food preservatives originate from molecules produced by bacteria, fungi or plants. In recent years, genome mining methodologies have been widely adopted to identify and characterize the biosynthetic gene clusters encoding the production of such compounds. Since 2011, the 'antibiotics and secondary metabolite analysis shell-antiSMASH' has assisted researchers in efficiently performing this, both as a web server and a standalone tool. Here, we present the thoroughly updated antiSMASH version 4, which adds several novel features, including prediction of gene cluster boundaries using the ClusterFinder method or the newly integrated CAS-SIS algorithm, improved substrate specificity prediction for non-ribosomal peptide synthetase adenylation domains based on the new SANDPUMA algorithm, improved predictions for terpene and ribosomally synthesized and post-translationally modified peptides cluster products, reporting of sequence similarity to proteins encoded in experimentally characterized gene clusters on a per-protein basis and a domain-level alignment tool for comparative analysis of trans-AT polyketide synthase assembly line architectures. Additionally, several usability features have been updated and improved. Together, these improvements make antiSMASH up-to-date with the latest developments in natural product research and will further facilitate computational genome mining for the discovery of novel bioactive molecules.

Background: The Rhynchosporium species complex consists of hemibiotrophic fungal pathogens specialized to different sweet grass species including the cereal crops barley and rye. A sexual stage has not been described, but several lines of evidence suggest the occurrence of sexual reproduction. Therefore, a comparative genomics approach was carried out to disclose the evolutionary relationship of the species and to identify genes demonstrating the potential for a sexual cycle. Furthermore, due to the evolutionary very young age of the five species currently known, this genus appears to be well-suited to address the question at the molecular level of how pathogenic fungi adapt to their hosts.Results: The genomes of the different Rhynchosporium species were sequenced, assembled and annotated using ab initio gene predictors trained on several fungal genomes as well as on Rhynchosporium expressed sequence tags. Structures of the rDNA regions and genome-wide single nucleotide polymorphisms provided a hypothesis for intra-genus evolution. Homology screening detected core meiotic genes along with most genes crucial for sexual recombination in ascomycete fungi. In addition, a large number of cell wall-degrading enzymes that is characteristic for hemibiotrophic and necrotrophic fungi infecting monocotyledonous hosts were found. Furthermore, the Rhynchosporium genomes carry a repertoire of genes coding for polyketide synthases and non-ribosomal peptide synthetases. Several of these genes are missing from the genome of the closest sequenced relative, the poplar pathogen Marssonina brunnea, and are possibly involved in adaptation to the grass hosts. Most importantly, six species-specific genes coding for protein effectors were identified in R. commune. Their deletion yielded mutants that grew more vigorously in planta than the wild type.Conclusion: Both cryptic sexuality and secondary metabolites may have contributed to host adaptation. Most importantly, however, the growth-retarding activity of the species-specific effectors suggests that host adaptation of R. commune aims at extending the biotrophic stage at the expense of the necrotrophic stage of pathogenesis. Like other apoplastic fungi Rhynchosporium colonizes the intercellular matrix of host leaves relatively slowly without causing symptoms, reminiscent of the development of endophytic fungi. Rhynchosporium may therefore become an object for studying the mutualism-parasitism transition.

ABSTRACT The assembly and annotation of the draft genome sequence of Shewanella sp. strain P1-14-1 are reported here to investigate the genes responsible for interkingdom interactions, secondary metabolite production, and microbial electrogenesis.

Basal lineages of the kingdom of fungi comprise terrestrial and aquatic fungi, which are traditionally summarized as Chytridiomycota and Zygomycota in a colloquial sense. The zygomycetes build up the most basal lineage among the terrestrial fungi. As the earliest fungi to conquer terrestrial habitats, they facilitated early coevolution with other land-dwelling microorganisms with implications on the primary and secondary metabolism which is adaptive to the life in soil. Advances in next-generation genome sequencing allow us to gain a new layer of knowledge with respect to their metabolism and its regulation. Recent revisions of the metabolic potential of zygomycetes using genome mining suggested a hidden potential for the production of secondary metabolites. First insights suggest a revision of our understanding of the importance and mode of action of metabolic pathways revealing a fascinating potential to produce novel natural compounds. The aim of this review is to draw attention to the zygomycetes as a unique fungal group and to inspire the scientific community to explore its potential as participants in complex ecosystems for innovative metabolic regulatory frameworks and novel compounds.

Motivation: Secondary metabolites (SM) are structurally diverse natural products of high pharmaceutical importance. Genes involved in their biosynthesis are often organized in clusters, i.e., are co-localized and co-expressed. In silico cluster prediction in eukaryotic genomes remains problematic mainly due to the high variability of the clusters’ content and lack of other distinguishing sequence features. Results: We present Cluster Assignment by Islands of Sites (CASSIS), a method for SM cluster prediction in eukaryotic genomes, and Secondary Metabolites by InterProScan (SMIPS), a tool for genome-wide detection of SM key enzymes (‘anchor’ genes): polyketide synthases, non-ribosomal peptide synthetases and dimethylallyl tryptophan synthases. Unlike other tools based on protein similarity, CASSIS exploits the idea of co-regulation of the cluster genes, which assumes the existence of common regulatory patterns in the cluster promoters. The method searches for ‘islands’ of enriched cluster-specific motifs in the vicinity of anchor genes. It was validated in a series of cross-validation experiments and showed high sensitivity and specificity. Availability and implementation: CASSIS and SMIPS are freely available at https://sbi.hki-jena.de/cassis. Contact: thomas.wolf@leibniz-hki.de or ekaterina.shelest@leibniz-hki.de Supplementary information: Supplementary data are available at Bioinformatics online.

To gain a broader understanding of the importance of a surface-associated lifestyle and morphogenic capability, we have assembled and annotated the genome sequences of Pseudoalteromonas strains P1-7a, P1-9, P1-13-1a, P1-16-1b, P1-25, and P1-26, isolated from Hydractinia echinata. These genomes will allow detailed studies on bacterial factors mediating interkingdom communication.

ABSTRACT The genomes of three Pseudoalteromonas strains (P1-8, P1-11, and P1-30) were sequenced and assembled. These genomes will inform future study of the genes responsible for the production of biologically active compounds responsible for these strains’ antimicrobial, biofouling, and algicidal activities.

Following antifungal treatment, Candida albicans, and other human pathogenic fungi can undergo microevolution, which leads to the emergence of drug resistance. However, the capacity for microevolutionary adaptation of fungi goes beyond the development of resistance against antifungals. Here we used an experimental microevolution approach to show that one of the central pathogenicity mechanisms of C. albicans, the yeast-to-hyphae transition, can be subject to experimental evolution. The C. albicans cph1Δ/efg1Δ mutant is nonfilamentous, as central signaling pathways linking environmental cues to hyphal formation are disrupted. We subjected this mutant to constant selection pressure in the hostile environment of the macrophage phagosome. In a comparatively short time-frame, the mutant evolved the ability to escape macrophages by filamentation. In addition, the evolved mutant exhibited hyper-virulence in a murine infection model and an altered cell wall composition compared to the cph1Δ/efg1Δ strain. Moreover, the transcriptional regulation of hyphae-associated, and other pathogenicity-related genes became re-responsive to environmental cues in the evolved strain. We went on to identify the causative missense mutation via whole genome- and transcriptome-sequencing: a single nucleotide exchange took place within SSN3 that encodes a component of the Cdk8 module of the Mediator complex, which links transcription factors with the general transcription machinery. This mutation was responsible for the reconnection of the hyphal growth program with environmental signals in the evolved strain and was sufficient to bypass Efg1/Cph1-dependent filamentation. These data demonstrate that even central transcriptional networks can be remodeled very quickly under appropriate selection pressure.