Evolution of the pathogenic mold Aspergillus fumigatus on high copper levels identifies novel resistance genes.

Aspergillus fumigatus is the leading cause of severe mold infections in immunocompromised patients. This common fungus possesses innate attributes that allow it to evade the immune system, including its ability to survive the high copper (Cu) levels in phagosomes. Our previous work has revealed that under high Cu levels, the A. fumigatus transcription factor AceA is activated, inducing the expression of the copper exporter CrpA to expel excess Cu. To identify additional elements in Cu resistance, we evolved A. fumigatus wild-type and mutant ΔaceA or ΔcrpA strains under increasing Cu concentrations. Sequencing of the resultant resistant strains identified both shared and unique evolutionary pathways to resistance. Reintroduction of three of the most common mutations in genes encoding Pma1 (plasma membrane H+-ATPase), Gcs1 (glutamate cysteine-ligase), and Cpa1 (carbamoyl-phosphate synthetase), alone and in combination, into wild-type A. fumigatus confirmed their additive role in conferring Cu resistance. Detailed analysis indicated that the pma1 mutation L424I preserves Pma1 H+-ATPase activity under high Cu concentrations and that the cpa1 mutation A37V confers a survival advantage to conidia in the presence of Cu. Interestingly, simultaneous mutations of all three genes did not alter virulence in infected mice. Our work has identified novel Cu-resistance pathways and provides an evolutionary approach for dissecting the molecular basis of A. fumigatus adaptation to diverse environmental challenges.IMPORTANCEAspergillus fumigatus is the most common mold infecting patients with weakened immunity. Infection is caused by the inhalation of mold spores into the lungs and is often fatal. In healthy individuals, spores are engulfed by lung immune cells and destroyed by a combination of enzymes, oxidants, and high levels of copper. However, the mold can protect itself by pumping out excess copper with specific transporters. Here, we evolved A. fumigatus under high copper levels and identified new genetic mutations that help it resist the toxic effects of copper. We studied how these mutations affect the mold's ability to resist copper and how they impact its ability to cause disease. This is the first such study in a pathogenic mold, and it gives us a better understanding of how it manages to bypass our body's defenses during an infection.

Two ARGONAUTE proteins loaded with transposon-derived small RNAs are associated with the reproductive cell lineage in Arabidopsis.

In sexually propagating organisms, genetic, and epigenetic mutations are evolutionarily relevant only if they occur in the germline and are hence transmitted to the next generation. In contrast to most animals, plants are considered to lack an early segregating germline, implying that somatic cells can contribute genetic information to progeny. Here we demonstrate that 2 ARGONAUTE proteins, AGO5 and AGO9, mark cells associated with sexual reproduction in Arabidopsis (Arabidopsis thaliana) throughout development. Both AGOs are loaded with dynamically changing small RNA populations derived from highly methylated, pericentromeric, long transposons. Sequencing of single stem cell nuclei revealed that many of these transposons are co-expressed within an AGO5/9 expression domain in the shoot apical meristem (SAM). Co-occurrence of transposon expression and specific ARGONAUTE (AGO) expression in the SAM is reminiscent of germline features in animals and supports the existence of an early segregating germline in plants. Our results open the path to investigating transposon biology and epigenome dynamics at cellular resolution in the SAM stem cell niche.

Identification and Analysis of Fungal-Specific Regions in the Aspergillus fumigatus Cu Exporter CrpA That Are Essential for Cu Resistance but Not for Virulence.

The opportunistic fungus Aspergillus fumigatus is the primary invasive mold pathogen in humans, and is responsible for an estimated 200,000 yearly deaths worldwide. Most fatalities occur in immunocompromised patients who lack the cellular and humoral defenses necessary to halt the pathogen's advance, primarily in the lungs. One of the cellular responses used by macrophages to counteract fungal infection is the accumulation of high phagolysosomal Cu levels to destroy ingested pathogens. A. fumigatus responds by activating high expression levels of crpA, which encodes a Cu+ P-type ATPase that actively transports excess Cu from the cytoplasm to the extracellular environment. In this study, we used a bioinformatics approach to identify two fungal-unique regions in CrpA that we studied by deletion/replacement, subcellular localization, Cu sensitivity in vitro, killing by mouse alveolar macrophages, and virulence in a mouse model of invasive pulmonary aspergillosis. Deletion of CrpA fungal-unique amino acids 1-211 containing two N-terminal Cu-binding sites, moderately increased Cu-sensitivity but did not affect expression or localization to the endoplasmic reticulum (ER) and cell surface. Replacement of CrpA fungal-unique amino acids 542-556 consisting of an intracellular loop between the second and third transmembrane helices resulted in ER retention of the protein and strongly increased Cu-sensitivity. Deleting CrpA N-terminal amino acids 1-211 or replacing amino acids 542-556 also increased sensitivity to killing by mouse alveolar macrophages. Surprisingly, the two mutations did not affect virulence in a mouse model of infection, suggesting that even weak Cu-efflux activity by mutated CrpA preserves fungal virulence.

Transcriptional response of Aspergillus fumigatus to copper and the role of the Cu chaperones.

Aspergillus fumigatus is the leading cause of life-threatening invasive mold infections in immunocompromised individuals. This ubiquitous saprophyte possesses several natural attributes allowing it to evade the immune system, including the ability to withstand high toxic Cu concentrations within the phagosomes of macrophages and neutrophils. We previously established that at high levels, Cu binds and activates the A. fumigatus transcription factor AceA, which upregulates the expression of the Cu exporter CrpA to expel excess Cu. Deletion of aceA or crpA result in extreme Cu sensitivity and attenuated virulence.To identify other elements participating in resistance to Cu, we performed a genome-wide analysis of the transcriptome by RNAseq to analyze the AceA-dependent response of A. fumigatus to excess Cu. We deleted key genes whose transcription was strongly upregulated by high Cu, including those encoding homologs of the three Cu chaperones cox17, atx1 and ccs1. Detailed analysis of these genes indicates that in A. fumigatus, cox17 is an essential gene with a possible role in respiration, the atxA gene product participates in reductive iron uptake and ccsA encodes the Cu chaperone activating A. fumigatus Sod1. Interestingly, although the ccsA-null strain was extremely sensitive to high Cu and oxidative stress, it was not attenuated in virulence in a mouse model of invasive pulmonary aspergillosis.Our work provides (i) a detailed view of the genome-wide transcriptional response of A. fumigatus to excess Cu, (ii) identification of the AceA-dependent transcriptome and (iii) analysis of the roles of the three Cu chaperones cox17, atxA and ccsA.

Point Mutation or Overexpression of Aspergillus fumigatus cyp51B, Encoding Lanosterol 14α-Sterol Demethylase, Leads to Triazole Resistance.

Aspergillus fumigatus is the most common cause of invasive fungal mold infections in immunocompromised individuals. Current antifungal treatment relies heavily on the triazole antifungals which inhibit fungal Erg11/Cyp51 activity and subsequent ergosterol biosynthesis. However, resistance, due primarily to cyp51 mutation, is rapidly increasing. A. fumigatus contains two Cyp51 isoenzymes, Cyp51A and Cyp51B. Overexpression and mutation of Cyp51A is a major cause of triazole resistance in A. fumigatus. The role of Cyp51B in generating resistance is unclear. Here, we show that overexpression or mutation of cyp51B results in triazole resistance. We demonstrate that introduction of a G457S Cyp51B mutation identified in a resistant clinical isolate results in voriconazole resistance in a naive recipient strain. Our results indicate that mutations in cyp51B resulting in clinical resistance do exist and should be monitored.

Dissection of floral transition by single-meristem transcriptomes at high temporal resolution.

The vegetative-to-floral transition is a dramatic developmental change of the shoot apical meristem, promoted by the systemic florigen signal. However, poor molecular temporal resolution of this dynamic process has precluded characterization of how meristems respond to florigen induction. Here, we develop a technology that allows sensitive transcriptional profiling of individual shoot apical meristems. Computational ordering of hundreds of tomato samples reconstructed the floral transition process at fine temporal resolution and uncovered novel short-lived gene expression programs that are activated before flowering. These programs are annulled only when both florigen and a parallel signalling pathway are eliminated. Functional screening identified genes acting at the onset of pre-flowering programs that are involved in the regulation of meristem morphogenetic changes but dispensable for the timing of floral transition. Induced expression of these short-lived transition-state genes allowed us to determine their genetic hierarchies and to bypass the need for the main flowering pathways. Our findings illuminate how systemic and autonomous pathways are integrated to control a critical developmental switch.

IRS1 phosphorylation underlies the non-stochastic probability of cancer cells to persist during EGFR inhibition therapy.

Stochastic transition of cancer cells between drug-sensitive and drug-tolerant persister phenotypes has been proposed to play a key role in non-genetic resistance to therapy. Yet, we show here that cancer cells actually possess a highly stable inherited chance to persist (CTP) during therapy. This CTP is non-stochastic, determined pre-treatment and has a unimodal distribution ranging from 0 to almost 100%. Notably, CTP is drug specific. We found that differential serine/threonine phosphorylation of the insulin receptor substrate 1 (IRS1) protein determines the CTP of lung and of head and neck cancer cells under epidermal growth factor receptor inhibition, both in vitro and in vivo. Indeed, the first-in-class IRS1 inhibitor NT219 was highly synergistic with anti-epidermal growth factor receptor therapy across multiple in vitro and in vivo models. Elucidation of drug-specific mechanisms that determine the degree and stability of cellular CTP may establish a framework for the elimination of cancer persisters, using new rationally designed drug combinations.

Single-cell analysis of clonal maintenance of transcriptional and epigenetic states in cancer cells.

Propagation of clonal regulatory programs contributes to cancer development. It is poorly understood how epigenetic mechanisms interact with genetic drivers to shape this process. Here, we combine single-cell analysis of transcription and DNA methylation with a Luria-Delbrück experimental design to demonstrate the existence of clonally stable epigenetic memory in multiple types of cancer cells. Longitudinal transcriptional and genetic analysis of clonal colon cancer cell populations reveals a slowly drifting spectrum of epithelial-to-mesenchymal transcriptional identities that is seemingly independent of genetic variation. DNA methylation landscapes correlate with these identities but also reflect an independent clock-like methylation loss process. Methylation variation can be explained as an effect of global trans-acting factors in most cases. However, for a specific class of promoters-in particular, cancer-testis antigens-de-repression is correlated with and probably driven by loss of methylation in cis. This study indicates how genetic sub-clonal structure in cancer cells can be diversified by epigenetic memory.

MetaCell: analysis of single-cell RNA-seq data using K-nn graph partitions.

scRNA-seq profiles each represent a highly partial sample of mRNA molecules from a unique cell that can never be resampled, and robust analysis must separate the sampling effect from biological variance. We describe a methodology for partitioning scRNA-seq datasets into metacells: disjoint and homogenous groups of profiles that could have been resampled from the same cell. Unlike clustering analysis, our algorithm specializes at obtaining granular as opposed to maximal groups. We show how to use metacells as building blocks for complex quantitative transcriptional maps while avoiding data smoothing. Our algorithms are implemented in the MetaCell R/C++ software package.

Multiple Auxin-Response Regulators Enable Stability and Variability in Leaf Development.

Auxin-signal transduction is mediated by the antagonistic activity of transcriptional activators and repressors. Both activators and repressors belong to gene families, but the biological importance of this complexity is not clear. Here, we addressed this question using tomato leaf development as a model by generating and analyzing mutants in multiple auxin-response components. In developing compound tomato leaves, auxin promotes leaflet formation and blade growth, and in the intercalary regions between leaflets, auxin response is inhibited by the Aux/IAA protein ENTIRE (E). e mutants form simple leaves due to ectopic blade growth in the intercalary domain. Using this unique loss-of-function phenotype and genome editing of auxin-response factor (ARF) genes, encoding auxin-response activators, we identified the contribution of specific ARFs to the e phenotype. Mutations in the related ARFs SlMP, SlARF19A, and SlARF19B, but not SlARF7, reduced the leaf blade and suppressed the e phenotype in a dosage-dependent manner that correlated with their relative expression, leading to a continuum of shapes. While single e and slmp mutants affected blade growth in an opposite manner, leaves of e slmp double mutants were similar to those of the wild type. However, the leaf shape of e slmp was more variable than that of the wild type, and it showed increased sensitivity to auxin. Our findings demonstrate that the existence of multiple auxin-response repressors and activators stabilizes the developmental output of auxin and that tuning their activity enables shape variability. The increased complexity of the auxin response therefore balances stability and flexibility in leaf patterning.

Riboflavin and pantothenic acid biosynthesis are crucial for iron homeostasis and virulence in the pathogenic mold Aspergillus fumigatus.

BACKGROUND: Aspergillus fumigatus is the most prevalent airborne fungal pathogen, causing invasive fungal infections mainly in immunosuppressed individuals. Death rates from invasive aspergillosis remain high because of limited treatment options and increasing antifungal resistance. The aim of this study was to identify key fungal-specific genes participating in vitamin B biosynthesis in A. fumigatus. Because these genes are absent in humans they can serve as possible novel targets for antifungal drug development. METHODS: By sequence homology we identified, deleted and analysed four key A. fumigatus genes (riboB, panA, pyroA, thiB) involved respectively in the biosynthesis of riboflavin (vitamin B2), pantothenic acid (vitamin B5), pyridoxine (vitamin B6) and thiamine (vitamin B1). RESULTS: Deletion of riboB, panA, pyroA or thiB resulted in respective vitamin auxotrophy. Lack of riboflavin and pantothenic acid biosynthesis perturbed many cellular processes including iron homeostasis. Virulence in murine pulmonary and systemic models of infection was severely attenuated following deletion of riboB and panA, strongly reduced after pyroA deletion and weakly attenuated after thiB deletion. CONCLUSIONS: This study reveals the biosynthetic pathways of the vitamins riboflavin and pantothenic acid as attractive targets for novel antifungal therapy. Moreover, the virulence studies with auxotrophic mutants serve to identify the availability of nutrients to pathogens in host niches. ABBREVIATIONS: BPS: bathophenanthrolinedisulfonate; BSA: bovine serum albumin; CFU: colony forming unit; -Fe: iron starvation; +Fe: iron sufficiency; hFe: high iron; NRPSs: nonribosomal peptide synthetases; PKSs: polyketide synthaseses; wt: wild type.

Vitamin Biosynthesis as an Antifungal Target.

The large increase in the population of immunosuppressed patients, coupled with the limited efficacy of existing antifungals and rising resistance toward them, have dramatically highlighted the need to develop novel drugs for the treatment of invasive fungal infections. An attractive possibility is the identification of possible drug targets within essential fungal metabolic pathways not shared with humans. Here, we review the vitamin biosynthetic pathways (vitamins A⁻E, K) as candidates for the development of antifungals. We present a set of ranking criteria that identify the vitamin B2 (riboflavin), B5 (pantothenic acid), and B9 (folate) biosynthesis pathways as being particularly rich in new antifungal targets. We propose that recent scientific advances in the fields of drug design and fungal genomics have developed sufficiently to merit a renewed look at these pathways as promising sources for the development of novel classes of antifungals.

A subset of conserved mammalian long non-coding RNAs are fossils of ancestral protein-coding genes.

BACKGROUND: Only a small portion of human long non-coding RNAs (lncRNAs) appear to be conserved outside of mammals, but the events underlying the birth of new lncRNAs in mammals remain largely unknown. One potential source is remnants of protein-coding genes that transitioned into lncRNAs. RESULTS: We systematically compare lncRNA and protein-coding loci across vertebrates, and estimate that up to 5% of conserved mammalian lncRNAs are derived from lost protein-coding genes. These lncRNAs have specific characteristics, such as broader expression domains, that set them apart from other lncRNAs. Fourteen lncRNAs have sequence similarity with the loci of the contemporary homologs of the lost protein-coding genes. We propose that selection acting on enhancer sequences is mostly responsible for retention of these regions. As an example of an RNA element from a protein-coding ancestor that was retained in the lncRNA, we describe in detail a short translated ORF in the JPX lncRNA that was derived from an upstream ORF in a protein-coding gene and retains some of its functionality. CONCLUSIONS: We estimate that ~ 55 annotated conserved human lncRNAs are derived from parts of ancestral protein-coding genes, and loss of coding potential is thus a non-negligible source of new lncRNAs. Some lncRNAs inherited regulatory elements influencing transcription and translation from their protein-coding ancestors and those elements can influence the expression breadth and functionality of these lncRNAs.

Systematic Identification of Anti-Fungal Drug Targets by a Metabolic Network Approach.

New antimycotic drugs are challenging to find, as potential target proteins may have close human orthologs. We here focus on identifying metabolic targets that are critical for fungal growth and have minimal similarity to targets among human proteins. We compare and combine here: (I) direct metabolic network modeling using elementary mode analysis and flux estimates approximations using expression data, (II) targeting metabolic genes by transcriptome analysis of condition-specific highly expressed enzymes, and (III) analysis of enzyme structure, enzyme interconnectedness ("hubs"), and identification of pathogen-specific enzymes using orthology relations. We have identified 64 targets including metabolic enzymes involved in vitamin synthesis, lipid, and amino acid biosynthesis including 18 targets validated from the literature, two validated and five currently examined in own genetic experiments, and 38 further promising novel target proteins which are non-orthologous to human proteins, involved in metabolism and are highly ranked drug targets from these pipelines.