Highlights d The ability for fast acquisition of drug resistance is widespread in Candida glabrata d Resistance-conferring mutations are very diverse but affect a small number of genes d Cross-resistance to fluconazole is common in strains adapted to anidulafungin d ERG3 mutations often drive fluconazole resistance in anidulafungin-adapted strains

Autoregulatory feedback loops occur in the regulation of molecules ranging from ATP to MAP kinases to zinc. Negative feedback loops can increase a system's robustness, while positive feedback loops can mediate transitions between cell states. Recent genome-wide experimental and computational studies predict hundreds of novel feedback loops. However, not all physical interactions are regulatory, and many experimental methods cannot detect self-interactions. Our understanding of regulatory feedback loops is therefore hampered by the lack of high-throughput methods to experimentally quantify the presence, strength and temporal dynamics of autoregulatory feedback loops. Here we present a mathematical and experimental framework for high-throughput quantification of feedback regulation and apply it to RNA binding proteins (RBPs) in yeast. Our method is able to determine the existence of both direct and indirect positive and negative feedback loops, and to quantify the strength of these loops. We experimentally validate our model using two RBPs which lack native feedback loops and by the introduction of synthetic feedback loops. We find that RBP Puf3 does not natively participate in any direct or indirect feedback regulation, but that replacing the native 3'UTR with that of COX17 generates an auto-regulatory negative feedback loop which reduces gene expression noise. Likewise, RBP Pub1 does not natively participate in any feedback loops, but a synthetic positive feedback loop involving Pub1 results in increased expression noise. Our results demonstrate a synthetic experimental system for quantifying the existence and strength of feedback loops using a combination of high-throughput experiments and mathematical modeling. This system will be of great use in measuring auto-regulatory feedback by RNA binding proteins, a regulatory motif that is difficult to quantify using existing high-throughput methods.

Fungal pathogens pose an increasingly worrying threat to human health, food security and ecosystem diversity. To tackle fungal infections and improve current diagnostic and therapeutic tools it is necessary to understand virulence and antifungal drug resistance mechanisms in diverse species. Recent advances in genomics approaches have provided a suitable framework to understand these phenotypes, which ultimately depend on genetically encoded determinants. In this work, we review how the study of genome sequences has been key to ascertain the bases of virulence and drug resistance traits. We focus on the contribution of comparative genomics, population genomics and directed evolution studies. In addition, we discuss how different types of genomic mutations (small or structural variants) contribute to intraspecific differences in virulence or drug resistance. Finally, we review current challenges in the field and anticipate future directions to solve them. In summary, this work provides a short overview of how genomics can be used to understand virulence and drug resistance in fungal pathogens.

Structural variants (SVs) like translocations, deletions, and other rearrangements underlie genetic and phenotypic variation. SVs are often overlooked due to difficult detection from short-read sequencing. Most algorithms yield low recall on humans, but the performance in other organisms is unclear. Similarly, despite remarkable differences across species’ genomes, most approaches use parameters optimized for humans. To overcome this and enable species-tailored approaches, we developed perSVade (personalized Structural Variation Detection), a pipeline that identifies SVs in a way that is optimized for any input sample. Starting from short reads, perSVade uses simulations on the reference genome to choose the best SV calling parameters. The output includes the optimally-called SVs and the accuracy, useful to assess the confidence in the results. In addition, perSVade can call small variants and copy-number variations. In summary, perSVade automatically identifies several types of genomic variation from short reads using sample-optimized parameters. We validated that perSVade increases the SV calling accuracy on simulated variants for six diverse eukaryotes, and on datasets of validated human variants. Importantly, we found no universal set of “optimal” parameters, which underscores the need for species-specific parameter optimization. PerSVade will improve our understanding about the role of SVs in non-human organisms.

Oxidative phosphorylation is among the most conserved mitochondrial pathways. However, one of the cornerstones of this pathway, the multi-protein complex NADH : ubiquinone oxidoreductase (complex I) has been lost multiple independent times in diverse eukaryotic lineages. The causes and consequences of these convergent losses remain poorly understood. Here, we used a comparative genomics approach to reconstruct evolutionary paths leading to complex I loss and infer possible evolutionary scenarios. By mining available mitochondrial and nuclear genomes, we identified eight independent events of mitochondrial complex I loss across eukaryotes, of which six occurred in fungal lineages. We focused on three recent loss events that affect closely related fungal species, and inferred genomic changes convergently associated with complex I loss. Based on these results, we predict novel complex I functional partners and relate the loss of complex I with the presence of increased mitochondrial antioxidants, higher fermentative capabilities, duplications of alternative dehydrogenases, loss of alternative oxidases and adaptation to antifungal compounds. To explain these findings, we hypothesize that a combination of previously acquired compensatory mechanisms and exposure to environmental triggers of oxidative stress (such as hypoxia and/or toxic chemicals) induced complex I loss in fungi.

Organisms regulate gene expression through changes in the activity of transcription factors (TFs). In yeast, the response of genes to changes in TF activity is generally assumed to be encoded in the promoter. To directly test this assumption, we chose 42 genes and, for each, replaced the promoter with a synthetic inducible promoter and measured how protein expression changes as a function of TF activity. Most genes exhibited gene-specific TF dose-response curves not due to differences in mRNA stability, translation, or protein stability. Instead, most genes have an intrinsic ability to buffer the effects of promoter activity. This can be encoded in the open reading frame and the 3' end of genes and can be implemented by both autoregulatory feedback and by titration of limiting trans regulators. We show experimentally and computationally that, when misexpression of a gene is deleterious, this buffering insulates cells from fitness defects due to misregulation.

Structural variants (SVs) underlie genomic variation but are often overlooked due to difficult detection from short reads. Most algorithms have been tested on humans, and it remains unclear how applicable they are in other organisms. To solve this, we develop perSVade (personalized structural variation detection), a sample-tailored pipeline that provides optimally called SVs and their inferred accuracy, as well as small and copy number variants. PerSVade increases SV calling accuracy on a benchmark of six eukaryotes. We find no universal set of optimal parameters, underscoring the need for sample-specific parameter optimization. PerSVade will facilitate SV detection and study across diverse organisms.

These authors contributed equally to this work.Auto regulatory feedback loops occur in the regulation of molecules ranging from ATP to MAP kinases to zinc. Negative feedback loops can increase a system's robustness, while positive feedback loops can mediate transitions between cell states. Recent genome-wide experimental and computational studies predict hundreds of novel feedback loops. However, not all physical interactions are regulatory, and many experimental methods cannot detect self-interactions. Our understanding of regulatory feedback loops is therefore hampered by the lack of high-throughput methods to experimentally quantify the presence, strength, and temporal dynamics of auto regulatory feedback loops. Here we present a mathematical and experimental framework for high-throughput quantification of feedback regulation, and apply it to RNA binding proteins (RBPs) in yeast. Our method is able to determine the existence of both direct and indirect positive and negative feedback loops, and to quantify the strength of these loops. We experimentally validate our model using two RBPs which lack native feedback loops, and by the introduction of synthetic feedback loops. We find that the the RBP Puf3 does not natively participate in any direct or indirect feedback regulation, but that replacing the native 3'UTR with that of COX17 generates an auto-regulatory negative feedback loop which reduces gene expression noise. Likewise, the RBP Pub1 does not natively participate in any feedback loops, but a synthetic positive feedback loop involving Pub1 results in increased expression noise. Our results demonstrate a synthetic experimental system for quantifying the existence and strength of feedback loops using a combination of high-throughput experiments and mathematical modeling. This system will be of great use in measuring auto-regulatory feedback by RNA binding proteins, a regulatory motif that is difficult to quantify using existing high-throughput methods.