Search

In this white paper, we describe the founding of a new ELIXIR Community - the Systems Biology Community - and its proposed future contributions to both ELIXIR and the broader community of systems biologists in Europe and worldwide. The Community believes that the infrastructure aspects of systems biology - databases, (modelling) tools and standards development, as well as training and access to cloud infrastructure - are not only appropriate components of the ELIXIR infrastructure, but will prove key components of ELIXIR’s future support of advanced biological applications and personalised medicine. By way of a series of meetings, the Community identified seven key areas for its future activities, reflecting both future needs and previous and current activities within ELIXIR Platforms and Communities. These are: overcoming barriers to the wider uptake of systems biology; linking new and existing data to systems biology models; interoperability of systems biology resources; further development and embedding of systems medicine; provisioning of modelling as a service; building and coordinating capacity building and training resources; and supporting industrial embedding of systems biology. A set of objectives for the Community has been identified under four main headline areas: Standardisation and Interoperability, Technology, Capacity Building and Training, and Industrial Embedding. These are grouped into short-term (3-year), mid-term (6-year) and long-term (10-year) objectives.

Fungi of the genus Aspergillus are widespread in the environment, where they produce large quantities of airborne conidia. Inhalation of Aspergillus spp. conidia in immunocompromised individuals can cause a wide spectrum of diseases, ranging from hypersensitivity responses to lethal invasive infections. Upon deposition in the lung epithelial surface, conidia encounter and interact with complex microbial communities that constitute the lung microbiota. The lung microbiota has been suggested to influence the establishment and growth of Aspergillus spp. in the human airways. However, the mechanisms underlying this interaction have not yet been sufficiently investigated. In this study, we aimed to enrich and isolate bacterial strains capable of inhibiting the germination and growth of A. fumigatus conidia from bronchoalveolar lavage fluid (BALF) samples of lung transplant recipients using a novel enrichment method. This method is based on a soft agar overlay plate assay in which bacteria are directly in contact with conidia, allowing inhibition to be readily observed during enrichment. We isolated a total of five clonal bacterial strains with identical genotypic fingerprints, as shown by random amplified polymorphic DNA PCR (RAPD–PCR). All strains were identified as Pseudomonas aeruginosa (strains b1–b5). The strains were able to inhibit the germination and growth of Aspergillus fumigatus in a soft agar confrontation assay, as well as in a germination multiplate assay. Moreover, when compared with ten P. aeruginosa strains isolated from expectoration through standard methods, no significant differences in inhibitory potential were observed. Additionally, we showed inhibition of A. fumigatus growth on Calu-3 cell culture monolayers. However, the isolated P. aeruginosa strains were shown to cause significant damage to the cell monolayers. Overall, although P. aeruginosa is a known opportunistic lung pathogen and antagonist of A. fumigatus, we validated this novel one-step enrichment approach for the isolation of bacterial strains antagonistic to A. fumigatus from BALF samples as a proof-of-concept. This opens up a new venue for the targeted enrichment of antagonistic bacterial strains against specific fungal pathogens.

Abstract Surface antigenic variation is crucial for major pathogens that infect humans. To escape the immune system, they exploit various mechanisms. Understanding these mechanisms is important to better prevent and fight the deadly diseases caused. Those used by the fungus Pneumocystis jirovecii that causes life-threatening pneumonia in immunocompromised individuals remain poorly understood. Here, though this fungus is currently not cultivable, our detailed analysis of the subtelomeric sequence motifs and genes encoding surface proteins suggests that the system involves the reassortment of the repertoire of ca. 80 non-expressed genes present in each strain, from which single genes are retrieved for mutually exclusive expression. Dispersion of the new repertoires, supposedly by healthy carrier individuals, appears very efficient because identical alleles are observed in patients from different countries. Our observations reveal a unique strategy of antigenic variation. They also highlight the possible role in genome rearrangements of small imperfect mirror sequences forming DNA triplexes.

In this white paper, we describe the founding of a new ELIXIR Community - the Systems Biology Community - and its proposed future contributions to both ELIXIR and the broader community of systems biologists in Europe and worldwide. The Community believes that the infrastructure aspects of systems biology - databases, (modelling) tools and standards development, as well as training and access to cloud infrastructure - are not only appropriate components of the ELIXIR infrastructure, but will prove key components of ELIXIR’s future support of advanced biological applications and personalised medicine. By way of a series of meetings, the Community identified seven key areas for its future activities, reflecting both future needs and previous and current activities within ELIXIR Platforms and Communities. These are: overcoming barriers to the wider uptake of systems biology; linking new and existing data to systems biology models; interoperability of systems biology resources; further development and embedding of systems medicine; provisioning of modelling as a service; building and coordinating capacity building and training resources; and supporting industrial embedding of systems biology. A set of objectives for the Community has been identified under four main headline areas: Standardisation and Interoperability, Technology, Capacity Building and Training, and Industrial Embedding. These are grouped into short-term (3-year), mid-term (6-year) and long-term (10-year) objectives.

In this white paper, we describe the founding of a new ELIXIR Community - the Systems Biology Community - and its proposed future contributions to both ELIXIR and the broader community of systems biologists in Europe and worldwide. The Community believes that the infrastructure aspects of systems biology - databases, (modelling) tools and standards development, as well as training and access to cloud infrastructure - are not only appropriate components of the ELIXIR infrastructure, but will prove key components of ELIXIR’s future support of advanced biological applications and personalised medicine. By way of a series of meetings, the Community identified seven key areas for its future activities, reflecting both future needs and previous and current activities within ELIXIR Platforms and Communities. These are: overcoming barriers to the wider uptake of systems biology; linking new and existing data to systems biology models; interoperability of systems biology resources; further development and embedding of systems medicine; provisioning of modelling as a service; building and coordinating capacity building and training resources; and supporting industrial embedding of systems biology. A set of objectives for the Community has been identified under four main headline areas: Standardisation and Interoperability, Technology, Capacity Building and Training, and Industrial Embedding. These are grouped into short-term (3-year), mid-term (6-year) and long-term (10-year) objectives.

In this white paper, we describe the founding of a new ELIXIR Community - the Systems Biology Community - and its proposed future contributions to both ELIXIR and the broader community of systems biologists in Europe and worldwide. The Community believes that the infrastructure aspects of systems biology - databases, (modelling) tools and standards development, as well as training and access to cloud infrastructure - are not only appropriate components of the ELIXIR infrastructure, but will prove key components of ELIXIR’s future support of advanced biological applications and personalised medicine. By way of a series of meetings, the Community identified seven key areas for its future activities, reflecting both future needs and previous and current activities within ELIXIR Platforms and Communities. These are: overcoming barriers to the wider uptake of systems biology; linking new and existing data to systems biology models; interoperability of systems biology resources; further development and embedding of systems medicine; provisioning of modelling as a service; building and coordinating capacity building and training resources; and supporting industrial embedding of systems biology. A set of objectives for the Community has been identified under four main headline areas: Standardisation and Interoperability, Technology, Capacity Building and Training, and Industrial Embedding. These are grouped into short-term (3-year), mid-term (6-year) and long-term (10-year) objectives.

Abstract New genomic tools open doors to study ecology, evolution, and population genomics of wild animals. For the Barn owl species complex, a cosmopolitan nocturnal raptor, a very fragmented draft genome was assembled for the American species (Tyto furcata pratincola) (Jarvis et al. 2014). To improve the genome, we assembled de novo Illumina and Pacific Biosciences (PacBio) long reads sequences of its European counterpart (Tyto alba alba). This genome assembly of 1.219 Gbp comprises 21,509 scaffolds and results in a N50 of 4,615,526 bp. BUSCO (Universal Single‐Copy Orthologs) analysis revealed an assembly completeness of 94.8% with only 1.8% of the genes missing out of 4,915 avian orthologs searched, a proportion similar to that found in the genomes of the zebra finch (Taeniopygia guttata) or the collared flycatcher (Ficedula albicollis). By mapping the reads of the female American barn owl to the male European barn owl reads, we detected several structural variants and identified 70 Mbp of the Z chromosome. The barn owl scaffolds were further mapped to the chromosomes of the zebra finch. In addition, the completeness of the European barn owl genome is demonstrated with 94 of 128 proteins missing in the chicken genome retrieved in the European barn owl transcripts. This improved genome will help future barn owl population genomic investigations.

As part of the human microbiota, the fungus Candida albicans colonizes the oral cavity and other mucosal surfaces of the human body. Commensalism is tightly controlled by complex interactions of the fungus and the host to preclude fungal elimination but also fungal overgrowth and invasion, which can result in disease. As such, defects in antifungal T cell immunity render individuals susceptible to oral thrush due to interrupted immunosurveillance of the oral mucosa. The factors that promote commensalism and ensure persistence of C. albicans in a fully immunocompetent host remain less clear. Using an experimental model of C. albicans oral colonization in mice we explored fungal determinants of commensalism in the oral cavity. Transcript profiling of the oral isolate 101 in the murine tongue tissue revealed a characteristic metabolic profile tailored to the nutrient poor conditions in the stratum corneum of the epithelium where the fungus resides. Metabolic adaptation of isolate 101 was also reflected in enhanced nutrient acquisition when grown on oral mucosa substrates. Persistent colonization of the oral mucosa by C. albicans also correlated inversely with the capacity of the fungus to induce epithelial cell damage and to elicit an inflammatory response. Here we show that these immune evasive properties of isolate 101 are explained by a strong attenuation of a number of virulence genes, including those linked to filamentation. De-repression of the hyphal program by deletion or conditional repression of NRG1 abolished the commensal behaviour of isolate 101, thereby establishing a central role of this factor in the commensal lifestyle of C. albicans in the oral niche of the host.

Dictyostelium discoideum amoebae feed by ingesting bacteria, then killing them in phagosomes. Ingestion and killing of different bacteria have been shown to rely on largely different molecular mechanisms. One would thus expect that D. discoideum adapts its ingestion and killing machinery when encountering different bacteria. In this study, we investigated by RNA sequencing if and how D. discoideum amoebae respond to the presence of different bacteria by modifying their gene expression patterns. Each bacterial species analyzed induced a specific modification of the transcriptome. Bacteria such as Bacillus subtilis, Klebsiella pneumoniae, or Mycobacterium marinum induced a specific and different transcriptional response, while Micrococcus luteus did not trigger a significant gene regulation. Although folate has been proposed to be one of the key molecules secreted by bacteria and recognized by hunting amoebae, it elicited a very specific and restricted transcriptional signature, distinct from that triggered by any bacteria analyzed here. Our results indicate that D. discoideum amoebae respond in a highly specific, almost non-overlapping manner to different species of bacteria. We additionally identify specific sets of genes that can be used as reporters of the response of D. discoideum to different bacteria.

Tubercular Mycobacteria and Legionella pneumophila are the causative agents of potentially fatal respiratory diseases due to their intrinsic pathogenesis but also due to the emergence of antibiotic resistance that limits treatment options. The aim of our study was to explore the antimicrobial activity of a small ligand-based chemical library of 1255 structurally diverse compounds. These compounds were screened in a combination of three assays, two monitoring the intracellular growth of the pathogenic bacteria, Mycobacterium marinum and L. pneumophila, and one assessing virulence of M. marinum. We set up these assays using two amoeba strains, the genetically tractable social amoeba Dictyostelium discoideum and the free-living amoeba Acanthamoeba castellanii. In summary, 64 (5.1%) compounds showed anti-infective/anti-virulence activity in at least one of the three assays. The intracellular assays hit rate varied between 1.7% (n = 22) for M. marinum and 2.8% (n = 35) for L. pneumophila with seven compounds in common for both pathogens. In parallel, 1.2% (n = 15) of the tested compounds were able to restore D. discoideum growth in the presence of M. marinum spiked in a lawn of food bacteria. We also validated the generality of the hits identified in the A. castellanii–M. marinum anti-infective screen using the D. discoideum–M. marinum host–pathogen model. The characterization of anti-infective and antibacterial hits in the latter infection model revealed compounds able to reduce intracellular growth more than 50% at 30 μM. Moreover, the chemical space and physico-chemical properties of the anti-M. marinum hits were compared to standard and candidate Mycobacterium tuberculosis (Mtb) drugs using ChemGPS-NP. A principle component analysis identified separate clusters for anti-M. marinum and anti-L. pneumophila hits unveiling the potentially new physico-chemical properties of these hits compared to standard and candidate M. tuberculosis drugs. Our studies underscore the relevance of using a combination of low-cost and low-complexity assays with full 3R compliance in concert with a rationalized focused library of compounds to identify new chemical scaffolds and to dissect some of their properties prior to taking further steps toward compound development.

Controlled immune activation in response to commensal microbes is critical for the maintenance of stable colonization and prevention of microbial overgrowth on epithelial surfaces. Our understanding of the host mechanisms that regulate bacterial commensalism has increased substantially, however, much less data exist regarding host responses to members of the fungal microbiota on colonized surfaces. Using a murine model of oropharyngeal candidiasis, we have recently shown that differences in immune activation in response to diverse natural isolates of Candida albicans are associated with different outcomes of the host-fungal interaction. Here we applied a genome-wide transcriptomic approach to show that rapid induction of a strong inflammatory response characterized by neutrophil-associated genes upon C. albicans colonization inversely correlated with the ability of the fungus to persist in the oral mucosa. Surprisingly, persistent fungal isolates showed no signs of a compensatory regulatory immune response. By combining RNA-seq data, genetic mouse models, and co-infection experiments, we show that attenuation of the inflammatory response at the onset of infection with a persistent isolate is not a consequence of enhanced immunosuppression. Importantly, depletion of regulatory T cells or deletion of the immunoregulatory cytokine IL-10 did not alter host-protective type 17 immunity nor did it impair fungal survival in the oral mucosa, indicating that persistence of C. albicans in the oral mucosa is not a consequence of suppressed antifungal immunity.

Arbuscular mycorrhizal fungi (AMF) are important symbionts of plants. Recently, studies of the AMF Rhizophagus irregularis recorded within-isolate genetic variation that does not completely match the proposed homokaryon or heterokaryon state (where heterokaryons comprise a population of two distinct nucleus genotypes). We re-analysed published data showing that bi-allelic sites (and their frequencies), detected in proposed homo- and heterokaryote R. irregularis isolates, were similar across independent studies using different techniques. This indicated that observed within-fungus genetic variation was not an artefact of sequencing and that such within- fungus genetic variation possibly exists. We then looked to see if bi-allelic transcripts from three R. irregularis isolates matched those observed in the genome as this would give a strong indication of whether bi-allelic sites recorded in the genome were reliable variants. In putative homokaryon isolates, very few bi-allelic transcripts matched those in the genome. In a putative heterokaryon, a large number of bi-allelic transcripts matched those in the genome. Bi-allelic transcripts also occurred in the same frequency in the putative heterokaryon as predicted from allele frequency in the genome. Our results indicate that while within-fungus genome variation in putative homokaryon and heterokaryon AMF was highly similar in 2 independent studies, there was little support that this variation is transcribed in homokaryons. In contrast, within-fungus variation thought to be segregated among two nucleus genotypes in a heterokaryon isolate was indeed transcribed in a way that is proportional to that seen in the genome.

Sleep is essential for optimal brain functioning and health, but the biological substrates through which sleep delivers these beneficial effects remain largely unknown. We used a systems genetics approach in the BXD genetic reference population (GRP) of mice and assembled a comprehensive experimental knowledge base comprising a deep "sleep-wake" phenome, central and peripheral transcriptomes, and plasma metabolome data, collected under undisturbed baseline conditions and after sleep deprivation (SD). We present analytical tools to interactively interrogate the database, visualize the molecular networks altered by sleep loss, and prioritize candidate genes. We found that a one-time, short disruption of sleep already extensively reshaped the systems genetics landscape by altering 60%-78% of the transcriptomes and the metabolome, with numerous genetic loci affecting the magnitude and direction of change. Systems genetics integrative analyses drawing on all levels of organization imply α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) receptor trafficking and fatty acid turnover as substrates of the negative effects of insufficient sleep. Our analyses demonstrate that genetic heterogeneity and the effects of insufficient sleep itself on the transcriptome and metabolome are far more widespread than previously reported.

ABSTRACT Pneumocystis species are opportunistic mammalian pathogens that cause severe pneumonia in immunocompromised individuals. These fungi are highly host specific and uncultivable in vitro. Human Pneumocystis infections present major challenges because of a limited therapeutic arsenal and the rise of drug resistance. To investigate the diversity and demographic history of natural populations of Pneumocystis infecting humans, rats, and mice, we performed whole-genome and large-scale multilocus sequencing of infected tissues collected in various geographic locations. Here, we detected reduced levels of recombination and variations in historical demography, which shape the global population structures. We report estimates of evolutionary rates, levels of genetic diversity, and population sizes. Molecular clock estimates indicate that Pneumocystis species diverged before their hosts, while the asynchronous timing of population declines suggests host shifts. Our results have uncovered complex patterns of genetic variation influenced by multiple factors that shaped the adaptation of Pneumocystis populations during their spread across mammals. IMPORTANCE Understanding how natural pathogen populations evolve and identifying the determinants of genetic variation are central issues in evolutionary biology. Pneumocystis, a fungal pathogen which infects mammals exclusively, provides opportunities to explore these issues. In humans, Pneumocystis can cause a life-threatening pneumonia in immunosuppressed individuals. In analysis of different Pneumocystis species infecting humans, rats, and mice, we found that there are high infection rates and that natural populations maintain a high level of genetic variation despite low levels of recombination. We found no evidence of population structuring by geography. Our comparisons of the times of divergence of these species to their respective hosts suggest that Pneumocystis may have undergone recent host shifts. The results demonstrate that Pneumocystis strains are widely disseminated geographically and provide a new understanding of the evolution of these pathogens.

ABSTRACT Microbial pathogens commonly escape the human immune system by varying surface proteins. We investigated the mechanisms used for that purpose by Pneumocystis jirovecii. This uncultivable fungus is an obligate pulmonary pathogen that in immunocompromised individuals causes pneumonia, a major life-threatening infection. Long-read PacBio sequencing was used to assemble a core of subtelomeres of a single P. jirovecii strain from a bronchoalveolar lavage fluid specimen from a single patient. A total of 113 genes encoding surface proteins were identified, including 28 pseudogenes. These genes formed a subtelomeric gene superfamily, which included five families encoding adhesive glycosylphosphatidylinositol (GPI)-anchored glycoproteins and one family encoding excreted glycoproteins. Numerical analyses suggested that diversification of the glycoproteins relies on mosaic genes created by ectopic recombination and occurs only within each family. DNA motifs suggested that all genes are expressed independently, except those of the family encoding the most abundant surface glycoproteins, which are subject to mutually exclusive expression. PCR analyses showed that exchange of the expressed gene of the latter family occurs frequently, possibly favored by the location of the genes proximal to the telomere because this allows concomitant telomere exchange. Our observations suggest that (i) the P. jirovecii cell surface is made of a complex mixture of different surface proteins, with a majority of a single isoform of the most abundant glycoprotein, (ii) genetic mosaicism within each family ensures variation of the glycoproteins, and (iii) the strategy of the fungus consists of the continuous production of new subpopulations composed of cells that are antigenically different. IMPORTANCE Pneumocystis jirovecii is a fungus causing severe pneumonia in immunocompromised individuals. It is the second most frequent life-threatening invasive fungal infection. We have studied the mechanisms of antigenic variation used by this pathogen to escape the human immune system, a strategy commonly used by pathogenic microorganisms. Using a new DNA sequencing technology generating long reads, we could characterize the highly repetitive gene families encoding the proteins that are present on the cellular surface of this pest. These gene families are localized in the regions close to the ends of all chromosomes, the subtelomeres. Such chromosomal localization was found to favor genetic recombinations between members of each gene family and to allow diversification of these proteins continuously over time. This pathogen seems to use a strategy of antigenic variation consisting of the continuous production of new subpopulations composed of cells that are antigenically different. Such a strategy is unique among human pathogens.

The stratigraphic succession of the Thakkhola region (central Nepal) begins in the Cambro-Ordovician with thick shallow-water carbonates, followed in the Siluro-Devonian by pelitic units deposited in deeper-water environments and at much higher accumulation rates with respect to the northwestern Himalaya sequence. During Carboniferous and Permian times, shallow-water limestones and coarse quartzose arenites were deposited, with major unconformities related to the initial opening of Neotethys. The following Triassic to Jurassic sedimentary units display comparable features all along the Tethys Himalaya. Owing to strong thermo-tectonic subsidence, the newly-formed Indian margin was rapidly drowned in the Early Triassic, and covered by thick marly limestones in Anisian to Carnian times. The Norian was characterized by very fine- grained quartzo-feldspathic detritus derived from the rejuvenated Indian foreland, while exclusive quartzose detritus suggests latitudinal drift towards more humid climates or more subdued relief in the Rhaetian. Shallow subtidal Early Jurassic limestones were covered in the Bajocian- Bathonian by bioclastic hybrid arenites, capped by a widespread condensed ironstone in the Early Callovian. In the Late Jurassic, black shales were deposited in outer shelf-slope conditions. The Cretaceous sedimentary evolution was related to the final fragmentation of Gondwanaland and initial opening of the Indian Ocean. Renewed quartzose detritus in the early? Neocomian was followed by deposition of thick deltaic to shallow-marine sandstones invariably characterized by abundant volcanic rock fragments. The northward-prograding deltaic complexes were overlain by a glauconitic condensed section, followed by late Albian pelagic marls. After the onset of collision berween India and Asia, the sedimentary units were involved in fold-thrust deformation at pre-anchimetamorphic to medium metamorphic grade, followed by uplift, formation of the Thakkhola graben and deposition of thick fluviolacustrine sediment in the Plio-Pleistocene.

The Carnian sedimentary succession of the Lierna area, more than 4OO m thick and belonging to the Coltignone tectonic unit (Grigne Group), has considerable paleogeographic significance, for it documents the western termination of the deltalagoon-carbonate platform depositional system of the Lombard Prealps. The Carnian succession, directly overlying basinal (Perledo-Varenna Formation) to shallow-marine (Lierna Formation) carbonates of Ladinian age, is characterized by deltaic terrigenous deposits (Val Sabbia Sandstone), progressively replaced by lagoonal (Gorno Formation) to peritidal (Breno Formation) limestones. In the studied stratigraphic sections these three lithofacies alternate vertically at decametric scale, and are each characterized by different types of cyclothems at metric scale. High-frequency "punctuated aggradational cycles" (PAC) are best recognized at the top of the Val Sabbia clastic body, where the final transition to the Gorno lagoon is marked by arenite layers containing ferruginous ooids. Grains of this type have nor been previously reported from Southalpine arenites. The petrographic composition of the Val Sabbia Sandstone points to provenance from an active volcanic chain, but abundance of quartzose detritus and occurrence of perthitic K-feldspar and sedimentary rock fragments suggest extensive erosion of the Hercynian basement and Upper Carboniferous to Ladinian sedimentary cover. These peculiar characteristics confirm the existence of a small "Lierna lobe" separated from the deltaic systems of central Lombardy, and fed by a basement- high located southwest of Como lake.

ABSTRACT Dermatophytes are the most common agents of superficial mycoses in humans and animals. The aim of the present investigation was to systematically identify the extracellular, possibly secreted, proteins that are putative virulence factors and antigenic molecules of dermatophytes. A complete gene expression profile of Arthroderma benhamiae was obtained during infection of its natural host (guinea pig) using RNA sequencing (RNA-seq) technology. This profile was completed with those of the fungus cultivated in vitro in two media containing either keratin or soy meal protein as the sole source of nitrogen and in Sabouraud medium. More than 60% of transcripts deduced from RNA-seq data differ from those previously deposited for A. benhamiae. Using these RNA-seq data along with an automatic gene annotation procedure, followed by manual curation, we produced a new annotation of the A. benhamiae genome. This annotation comprised 7,405 coding sequences (CDSs), among which only 2,662 were identical to the currently available annotation, 383 were newly identified, and 15 secreted proteins were manually corrected. The expression profile of genes encoding proteins with a signal peptide in infected guinea pigs was found to be very different from that during in vitro growth when using keratin as the substrate. Especially, the sets of the 12 most highly expressed genes encoding proteases with a signal sequence had only the putative vacuolar aspartic protease gene PEP2 in common, during infection and in keratin medium. The most upregulated gene encoding a secreted protease during infection was that encoding subtilisin SUB6, which is a known major allergen in the related dermatophyte Trichophyton rubrum. IMPORTANCE Dermatophytoses (ringworm, jock itch, athlete’s foot, and nail infections) are the most common fungal infections, but their virulence mechanisms are poorly understood. Combining transcriptomic data obtained from growth under various culture conditions with data obtained during infection led to a significantly improved genome annotation. About 65% of the protein-encoding genes predicted with our protocol did not match the existing annotation for A. benhamiae. Comparing gene expression during infection on guinea pigs with keratin degradation in vitro, which is supposed to mimic the host environment, revealed the critical importance of using real in vivo conditions for investigating virulence mechanisms. The analysis of genes expressed in vivo, encoding cell surface and secreted proteins, particularly proteases, led to the identification of new allergen and virulence factor candidates.

ABSTRACT In vivo transcriptional analyses of microbial pathogens are often hampered by low proportions of pathogen biomass in host organs, hindering the coverage of full pathogen transcriptome. We aimed to address the transcriptome profiles of Candida albicans, the most prevalent fungal pathogen in systemically infected immunocompromised patients, during systemic infection in different hosts. We developed a strategy for high-resolution quantitative analysis of the C. albicans transcriptome directly from early and late stages of systemic infection in two different host models, mouse and the insect Galleria mellonella. Our results show that transcriptome sequencing (RNA-seq) libraries were enriched for fungal transcripts up to 1,600-fold using biotinylated bait probes to capture C. albicans sequences. This enrichment biased the read counts of only ~3% of the genes, which can be identified and removed based on a priori criteria. This allowed an unprecedented resolution of C. albicans transcriptome in vivo, with detection of over 86% of its genes. The transcriptional response of the fungus was surprisingly similar during infection of the two hosts and at the two time points, although some host- and time point-specific genes could be identified. Genes that were highly induced during infection were involved, for instance, in stress response, adhesion, iron acquisition, and biofilm formation. Of the in vivo-regulated genes, 10% are still of unknown function, and their future study will be of great interest. The fungal RNA enrichment procedure used here will help a better characterization of the C. albicans response in infected hosts and may be applied to other microbial pathogens. IMPORTANCE Understanding the mechanisms utilized by pathogens to infect and cause disease in their hosts is crucial for rational drug development. Transcriptomic studies may help investigations of these mechanisms by determining which genes are expressed specifically during infection. This task has been difficult so far, since the proportion of microbial biomass in infected tissues is often extremely low, thus limiting the depth of sequencing and comprehensive transcriptome analysis. Here, we adapted a technology to capture and enrich C. albicans RNA, which was next used for deep RNA sequencing directly from infected tissues from two different host organisms. The high-resolution transcriptome revealed a large number of genes that were so far unknown to participate in infection, which will likely constitute a focus of study in the future. More importantly, this method may be adapted to perform transcript profiling of any other microbes during host infection or colonization.