G3

Fungal diseases cause devastating crop infestations and millions of deaths annually. To improve disease treatment development, we must understand the biological processes these treatments target and the fungal cellular responses they elicit. Current methods to identify genetic responses of fungi to antifungal drugs and fungicides typically involve laborious molecular genetics experiments or transcriptomics. We utilized high-throughput screens in the basidiomycete yeast Rhodotorula (Rhodosporidium) toruloides to rapidly identify genes associated with the fungal response to a copper-zinc fungicide developed by VM Agritech and compared identified genes to those transcriptionally activated in response to the fungicide in R. toruloides and the filamentous ascomycete Neurospora crassa. These 'omics technologies combined with molecular genetics and biochemical analyses showed copper and zinc synergize to inhibit fungal growth due, at least in part, to an increase in cellular zinc concentrations when copper is present. Intriguingly, although we identified genes involved in the copper-zinc fungicide response using both high-throughput screens and transcriptomics, there was little overlap in genes identified through the two methods. Our findings demonstrate the power of combining functional genomics and transcriptomics approaches in identifying a more comprehensive map of the genetic mechanisms necessary for responding to fungicides.

The Candida Genome Database (CGD; www.candidagenome.org) is both a model organism database and a fungal pathogen database. As a model organism database, CGD stores data for Candida albicans, which serves as a model species both for other Candida spp. and for non-Candida fungi that form biofilms and undergo routine morphogenic switching. As a fungal pathogen database, CGD now hosts locus pages for six species of the best-studied pathogenic fungi in the Candida group. Pathogenic Candida species have become increasingly drug resistant and there is thus a pressing need for research into basic Candida biology, epidemiology, phylogeny, and potential new antifungals, as well as a single location where all of the available data are collected, curated, and made easily searchable. CGD curates the gene-based Candida experimental literature in real time, extracting, organizing and standardizing gene annotations. CGD also links clinical data on disease to relevant Literature Topics to improve searchability for clinical researchers. Because CGD curates the literature for multiple species and most research focuses on aspects related to pathogenicity, we focus our curation efforts on assigning Literature Topic tags, collecting detailed mutant phenotype data, and assigning controlled Gene Ontology terms with accompanying evidence codes. Our Summary pages for each locus include the primary name and all aliases for that locus, a description of the gene and/or gene product, detailed ortholog information with links, a synteny view, a JBrowse window with a visual view of the gene on its chromosome, links to Phenotype, Gene Ontology, Interactions, and Expression pages, as well as sequence information, references cited on the summary page itself, and any locus notes. The database also serves as a community hub, where we link to various types of reference material of relevance to Candida researchers, including colleague information, news, and notice of upcoming meetings. We routinely survey the community to learn how the field is evolving and how needs may have changed. Here we describe CGDs new modern web interface and multiple new tools that have been added in the last 6 months, allowing, among other things, users to better understand the available expression data for a locus and seamlessly switch between species for a given locus.

Biofilm formation by the fungal pathogen Candida albicans is a major source of infection. The biofilm regulator Ume6 forms complexes with partner transcription factors Efg1, Ndt80, and Upc2 to drive expression of biofilm-related genes. Here we present chromatin immunoprecipitation with sequencing (ChIP-seq) data for Ume6-chromatin association in both wild-type and efg1{Delta}/{Delta} ndt80{Delta}/{Delta} upc2{Delta}/{Delta} backgrounds. Ume6 associates with two promoter region classes. For one class, Ume6 association is significantly decreased in the efg1{Delta}/{Delta} ndt80{Delta}/{Delta} upc2{Delta}/{Delta} background. This class includes 577 genes, many with well-established roles in biofilm formation or the related process of filamentation. For the second class, Ume6 association is unperturbed in the triple mutant background. This class includes 534 genes, some with roles in iron homeostasis (e.g., SFU1) and ergosterol synthesis (e.g., ERG11), which may contribute to the impact of Ume6 on ndt80{Delta}/{Delta} mutant azole drug sensitivity. It also includes FLO9, a putative adhesin gene that is required for the emergent biofilm state produced by efg1{Delta}/{Delta} ndt80{Delta}/{Delta} double mutants. The data show that Ume6 binding to ~half of its targets requires known partners, and Ume6 binding to ~half of its targets does not. Ume6 may bind to the latter set of promoter regions independently of any partner, or perhaps with additional partners that have yet to be discovered.

Environmental bacteria have abundant signaling systems wired into complex gene regulatory networks to adapt to fluctuating conditions. In Myxococcus xanthus, starvation triggers a developmental program (specialized biofilm) that produces spore-filled multicellular fruiting bodies surrounded by a distinct quiescent state termed peripheral rods. Fruiting body structure as well as the proportion of cells following each fate can be tuned by a large repertoire of signaling proteins, including numerous orphan histidine kinases. Here, we focus on the histidine kinase TodK which was previously demonstrated to influence developmental progression. We find that loss of TodK produces distinct developmental phenotypes that vary with environmental conditions. To quantify these effects, we developed an image-analysis pipeline that measures aggregation and fruiting body patterning during development on nutrient-limited agar. These analyses revealed the todK mutant precociously aggregates particularly at the peripheries of the colony. Under submerged-culture conditions, initial production of aggregates was not accelerated but aggregates exhibited accelerated progression to mature fruiting bodies. Overexpression of active TodK completely blocked fruiting body formation. Molecular analyses demonstrated that TodK overproduction suppressed expression of core developmental regulators including FruA and CsgA (C-signal). Interestingly, protein accumulation of MrpC, necessary for expression of both FruA and the C-signal was not significantly perturbed suggesting TodK silences MrpC transcriptional activity. Together, these findings establish TodK as a modulator of developmental progression and demonstrate how quantitative phenotyping approaches can reveal biologically meaningful functions for orphan histidine kinases whose mutant phenotypes might otherwise appear subtle. Summary StatementQuantitative analysis of multicellular development reveals previously hidden functions of an orphan histidine kinase, highlighting the importance of robust phenotyping approaches for understanding bacterial signaling networks.

Recent work has established septation as a critical process for virulence in A. fumigatus, as strains lacking septa are avirulent. While our work has focused on the well-studied Septation Initiation Network (SIN), understanding of septation machinery upstream and downstream of SIN remains limited in filamentous fungi. Proximity labeling techniques are powerful tools to study pathway interactions, especially those that are transient. Recent advances have seen TurboID used to study pathways and processes important for fungal pathogenesis. In this study, we adapt TurboID for use in Aspergillus fumigatus and use overlapping datasets of two SIN components, the terminal NDR kinase SidB and its binding partner MobA, to identify putative septation effectors. Many known septation effectors, including septins and cell wall synthesis proteins, were enriched in our datasets, validating the approach. Phenotypic characterization of hits identified by LC-MS/MS revealed several previously uncharacterized gene products involved in septum formation and cell wall stress tolerance. These included the PP2A regulatory subunit PabA and IQGAP SepG, which have been shown to regulate septation in Aspergillus nidulans. Thus, TurboID is a useful tool to study fungal signaling and physiology in A. fumigatus and may be used to improve understanding of pathogenic processes.

During meiosis, accurate chromosome segregation requires significant condensation and compaction. These processes are mediated by condensins, cohesins, and histone tail modifications. We identified that MET-2, a histone methyltransferase that catalyzes the dimethylation of histone H3 lysine 9 (H3K9me2), differentially impacts chromosome size in the male vs. female C. elegans germline. In met-2 null worms, autosomes during spermatogenesis are significantly larger than wild type, while chromosome size during oogenesis is unaffected. X-univalent size in males is also unaffected by loss of MET-2, indicating MET-2 differentially regulates autosomal and X-chromosome compaction in male spermatogenesis. Autosome size is not changed when males harbor a catalytically deficient MET-2 (met-2CD) or have mutations preventing germline histone H3K9 methylation (H3K9R). In addition, met-2 males, in contrast to met-2CD or H3K9R males, have more active RNA pol II in later stages of meiosis. These data suggest MET-2 plays a noncatalytic role in mediating chromosome structure and transcription. In met-2 male germ lines, genes on the X chromosome, which is typically enriched in H3K9me2, are significantly more likely to be upregulated than genes on autosomes, even though X-univalent size is unchanged. These results suggest that MET-2 plays a sex-specific role that is not limited to its enzymatic activity.

Anthropogenic copper (Cu) contamination from agriculture, mining, and industrial runoff creates environmental gradients affecting physiology and behavior in wild populations. While Cu toxicity in Drosophila melanogaster is well characterized, it remains unclear whether Cu resistance is one integrated trait or several independently evolving components. Using a subset of recombinant inbred lines (RILs) from the Drosophila Synthetic Population Resource (DSPR), we measured three components of Cu response: feeding avoidance, oviposition avoidance, and physiological tolerance (median lethal time, LT50) under sustained Cu exposure. All three traits showed substantial phenotypic variation among RILs. Feeding and oviposition avoidance were both highly heritable (H 2 ~ 0.88), and RIL identity accounted for 49.5% of the variance in LT50. However, the three traits showed no significant correlation across RILs, indicating distinct genetic architecture. We identified a single male specific quantitative trait locus (QTL) on chromosome 2R that explained 17.7% of the variation in feeding preference; the interval included candidate detoxification genes Jheh1, Jheh2, Jheh3 and sano, the latter of which is associated with olfactory behavior. No significant QTL were detected for oviposition preference, suggesting a highly polygenic structure that may difficult to detect with our limited panel size. Together, these results indicate that Cu resistance in D. melanogaster is genetically modular. Behavioral avoidance during feeding, oviposition, and physiological tolerance are heritable but architecturally distinct components, each with potential to respond to selection independently.

Introduction. Histoplasma capsulatum is an important fungal pathogen of humans, especially in immunocompromised patients, and the World Health Organisation has listed Histoplasma as a High Priority fungal pathogen. The taxonomy of the genus Histoplasma is still under development, but its only species H. capsulatum has at least of seven lineages, and a range of sublineages. Gap Statement. Allelic multilocus sequence typing (MLST) has been very successfully applied for typing and epidemiology of pathogenic bacteria, but its adoption for fungal pathogens is more limited. Genome sequence-based MLST approaches promise much a higher resolution for comparative phylogenetics, but there are only a small number of genome assemblies available for H. capsulatum. Aim. To develop and validate genome-based MLST typing for H. capsulatum. Methodology. A total of 406 genome assemblies were generated from publicly available Illumina short-read sequencing samples. Selected genome assemblies were subsequently used to generate MLST schemes with chewBBACA, based on whole genome sequences (wgMLST), NCBI-obtained annotations (annMLST) and annotations obtained with Funannotate (fanMLST), and used for lineage identification, and comparative analyses of population structure. Results. Genome assemblies were successfully generated for 406 H. capsulatum samples, with BUSCO completeness scores of over 98%. A subset of 60 genome assemblies was used to generate and validate MLST schemes based on the 60 genome assemblies, 5 genome annotations obtained from Genbank, and 13 genome annotations generated using Funannotate, resulting in schemes with > 8,000 markers. Each of the MLST schemes allowed reproducible separation of all 7 lineages and 6 H. capsulatum Suramericanum sublineages, without the need for high-performance clusters or powerful computing equipment. Comparison of the different schemes using tanglegrams and scatter plots showed good reproducibility (r2 values of 0.8-1.0). Comparison of total sizes of genome assemblies showed significant lineage-dependent size differences, with the Suramericanum, LAmB, India, Capsulatum and Mississippiense lineages averaging a total genome size of 30-33 million bp, and the Africa and Ohiense lineages averaging at 39 million bp. Conclusion. Genome sequence- and especially genome annotation-based MLST analyses are a promising approach that can be adapted as a rapid, high-throughput method for genomic clustering and epidemiological tracing of pathogenic fungi.

During pathogen infection, the C. elegans transcription factor SKN-1 is activated through the p38 MAPK cascade to protect against oxidative damage and promote host survival. SKN-1, the functional ortholog of the mammalian Nrf family of transcription factors, participates in various biological processes and is subject to complex regulation. In this study, we identify a previously unrecognized role for LIN-23 in regulating SKN-1 during adult stress conditions. LIN-23 is an F-box protein that functions as the substrate recognition component of the Skp-Cullin-F-box (SCF) E3 ubiquitin ligase complex and has been implicated in diverse cellular processes, including cell cycle regulation, neurite outgrowth, and centrosome duplication. Although LIN-23 has previously been reported to negatively regulate SKN-1 in other contexts, our findings demonstrate that during pathogen exposure, LIN-23 acts as a positive regulator of SKN-1 activity. Specifically, loss of LIN-23 reduced SKN-1 activity and decreased survival of infected adult animals. Having established this novel relationship between LIN-23 and SKN-1, we investigated the mechanism by which LIN-23 regulates SKN-1 activity. Because SKN-1 activation occurs via the p38 MAPK signaling pathway, followed by nuclear localization, we examined whether LIN-23 influences these events, but observed no decrease in p38 MAPK phosphorylation or SKN-1 nuclear localization. Instead, we provide evidence that LIN-23 function is dependent on WDR-23, a well-established negative regulator of SKN-1. A model is proposed in which LIN-23 promotes SKN-1 activity by targeting nuclear WDR-23 for degradation. SUMMARYSKN-1 is a C. elegans transcription factor and the ortholog of mammalian Nrf proteins. Under oxidative stress conditions, including those induced by infection, SKN-1 plays a protective role. Since SKN-1 regulation is complex, understanding the mechanisms that modulate its activity is important for defining stress response pathways. The authors demonstrate that the F-box protein LIN-23 functions as a positive regulator of SKN-1. Their genetic analyses indicate that LIN-23 does not influence the canonical SKN-1 activation cascade. Instead, LIN-23 appears to regulate SKN-1 activity by modulating a negative regulator, WDR-23.

Filamentation, a central virulence trait of the fungal pathogen Candida albicans, is required for host cell damage, tissue invasion, and biofilm formation. The transcription factor Bcr1 was first found to be required for biofilm formation but not filamentation. Subsequent studies revealed a negative role for Bcr1 in filamentation in opaque cells, a cell type that is required for mating, under hypoxic conditions, and in a wor2del/del mutant background. Here we characterize a new context in which Bcr1 negatively regulates filamentation and present its associated gene expression impact. We compared the wild-type and bcr1del/del mutant cells in the SC5314 reference strain background at 30 degrees C in glucose media, conditions that do not induce filamentation. The bcr1del/del mutant displayed increased invasive growth in a solid medium, and some filamentation ability in a liquid medium. The bcr1del/del effect on filamentation was augmented by overexpression of hyphal cyclin gene HGC1, which caused pseudohyphal growth in wild-type, efg1del/del, or brg1del/del strains, yet caused hyphal growth in the bcr1del/del strain. RNA-sequencing (RNA-seq) analysis of 30 degrees C cells shows that the bcr1del/del mutant has elevated expression of two genes that drive filamentation: HGC1 and UME6. In published data these same two genes have decreased expression in the bcr1del/del mutant at 37 degrees C, the temperature at which biofilm formation is typically assayed. A growing cadre of biofilm/hyphal regulators can have both positive and negative effects on hypha-associated gene expression, including Bcr1, Efg1, Ndt80, and Nrg1.

A complete reference genome for species with chromosomally-determined separate sexes should contain scaffolds for all sex chromosome homologs. However, sex chromosomes present distinct computational challenges compared to autosomes. Here we present a k-mer based analysis that utilizes whole-genome sequencing of a few sex-identified isolates: Cytogenetics-By-Sequencing (CBS) tools. Unlike other approaches that typically address one aspect of the sex chromosomes, CBS-tools strives to guide users from the discovery of the heterogametic sex through identifying the sex-determination region (SDR). The core of CBS-tools is automated quantification of sex-specific k-mers in order to predict the heterogametic sex. Using publicly-available datasets, CBS-tools correctly identified the known sex-system of the 31 species tested. Additionally, we used these k-mers to verify and correct phasing of sex chromosomes between haplotypes in species representing different sex-systems. Finally, we used these k-mers to delimit the SDR boundary using an interactive web platform. CBS-tools was developed with previously unexplored sex chromosome systems in mind, but is also suitable for well-examined sex chromosome pairs.

To face the fungal infections that affect millions of people every year, a better understanding of resistance-conferring mutations is necessary. One obstacle is the considerable time and effort necessary to construct point mutations in pathogenic species, which is required to validate their role in resistance. To overcome this problem, we have constructed a series of heterologous expression plasmids that allow expression of the ERG3, ERG6 and ERG11 orthologs from Candida albicans, Candidozyma auris, Nakaseomyces glabratus, Candida parapsilosis, Candida tropicalis, and the CYP51A and CYP51B genes of Aspergillus fumigatus. All orthologs are codon-optimized and synthesized for expression in the budding yeast Saccharomyces cerevisiae. We also constructed host strains in which the native ERG genes are either deleted or under the control of a repressible promoter. We showed that all orthologs complement the function of the native genes when expressed from the S. cerevisiae native promoters. We recreated previously reported resistance-conferring mutations, and showed that expression in S. cerevisiae replicates known resistance phenotypes. We then used the system to test corresponding mutations across orthologs, and showed that most resistance phenotypes are conserved. The heterologous expression system allows rapid construction of mutants, by taking advantage of the abundant genetic tools available for S. cerevisiae. Screening of mutants can be done within a few experimental steps, greatly increasing the throughput compared to the introduction of mutations in pathogenic fungi. As a proof-of-concept, we used the system to discover new resistance-conferring mutations, using random mutagenesis. All plasmids are made publicly available.

Long breeding cycles constrain genetic gain in Scots pine, while mature progeny trials can provide reference populations for genomic selection. We combined phenotypic records from two German trials established in 1990 with dense SNP data. We complemented these data with an offspring-level parentage audit, duplicate filtering, and simulations of breeding strategies. Of 1,986 phenotype-matched genotyped trees, 1,668 formed the parentage-green set and 1,651 remained after exclusion of 17 near-duplicate samples. The final population comprised 1,521 supported controlled-cross offspring and 130 mother-known, father-unknown offspring from 87 progeny labels. Across 14 trait-by-age measurements, GBLUP gave the clearest improvements for diameter and volume from age 20 onwards, whereas height was more mixed. At age 35, ABLUP versus GBLUP heritability was 0.155 versus 0.235 for diameter, 0.265 versus 0.263 for height, and 0.172 versus 0.236 for volume. Leave-progeny-out predictive ability increased from 0.205 to 0.238, 0.249 to 0.265, and 0.203 to 0.239, respectively. SNPscan_breeder simulations compared phenotypic selection, progeny testing, cross-generation genomic selection, and genomic selection with phenotypic thinning under five diversity variants. Progeny testing produced the greatest cumulative gain, but its 45-year cycle reduced annual response. Under the base assumptions, genomic selection with phenotypic thinning gave the highest annual gains for height and fungal resistance, whereas pure genomic selection gave the highest annual diameter gain. Genomic strategies accumulated more kinship than conventional strategies, although a {lambda} = 0.10 kinship penalty improved founder retention with little loss of gain; no mitigation option simultaneously maximised gain, prediction accuracy, and diversity. These results support genomic shortlisting within a field-tested programme with explicit reference-population updating and diversity management.

Pseudogymnoascus destructans is a psychrophilic ascomycete fungus that causes White-Nose Syndrome in North American hibernating bats. P. destructans invades the tissue present in bats' noses and wing membranes, causing the characteristic white nose and cupping erosion in the wings. P. destructans exhibits characteristics of a hemibiotrophic fungal pathogen, in which P. destructans can invade the host cells without inducing damage. During this stage, P. destructans hyphae assume a morphology similar to fungal plant pathogens during the biotrophic stage of infection. However, the mechanism by which P. destructans regulates its morphology and virulence is not fully understood. The phosphatase calcineurin is a key regulator of fungal stress response and virulence and could play a role in P. destructans virulence. Calcineurin can be inhibited using FK506 or Cyclosporin A (CSA). FK506 interacts with FKBP12, while CSA interacts with Cyclosporin A, and the protein-drug complex then inhibits calcineurin activity. Here, we show that P. destructans is intrinsically resistant to one calcineurin inhibitor (FK506), possibly because P. destructans FKBP12 has a lysine instead of the conserved arginine at position 89, which has been shown to reduce the affinity of the FKBP12-FK506 complex for calcineurin. Using cyclosporin A, we found that calcineurin is important for fungal growth. Additionally, calcineurin inhibition increases voriconazole inhibition. In contrast, calcineurin inhibition was antagonistic to the anti-cell wall drug micafungin. Taken together, calcineurin's role in fungal growth and response to voriconazole is conserved in P. destructans but may play a unique role in P. destructans response to micafungin.

The cultivated dahlia has a complex evolutionary history that is likely the driver for its diverse floral forms. Previous studies have proposed cultivated dahlia as an allotetraploid, autotetraploid, and in one study an auto-allooctoploid. Here, we used PacBio HiFi long-reads and Dovetail Omni-C reads to assemble and annotate a genome for a tetraploid Dahlia variabilis cultivar Edna C. We identified Edna C as a segmental tetraploid with signature of auto- and allo-polyploidy events based on subgenome-specific repeat partitioning and single-copy gene tree construction within homoeologous groups. Additionally, we identified a cultivated dahlia-specific centromeric monomer shared between Dahlia variabilis Edna C and the Dahlia pinnata Kelvin Floodlight. Based on finding three major chromosome clusters, we propose that Dahlia variabilis Edna C is a segmental tetraploid with evidence of auto- and allopolyploidy.

Leptosphaeria maculans is a plant-pathogenic fungus that infects Brassica species, including Brassica napus (oilseed rape). Breeding oilseed rape varieties with genetic resistance is an efficient way to control the disease; however, L. maculans can adapt and overcome these resistances. Understanding the mechanisms that enable L. maculans to adapt is crucial for managing the emergence of better-adapted isolates. Brassica carinata, the Ethiopian Mustard, although closely related to B. napus, is considered a nonhost species of L. maculans because this fungus cannot infect it. Despite the extreme resistance of B. carinata, one natural L. maculans isolate has been identified as unable to infect B. napus, causing moderate and atypical symptoms on this species. We performed a cross between this isolate and an isolate adapted to B. napus, followed by a QTL analysis, which identified seven QTL, each encompassing candidate genes involved in L. maculans adaptation to B. carinata or B. napus. Additionally, we observed transgression in the progeny, wherein a few strains caused significantly more or less aggressive symptoms on both species of Brassica. We found that epistasis within the L. maculans genome contributes to the observed transgression. These initial findings provide further opportunities to study the adaptive capacities of L. maculans, as well as data to initiate analysis of the extreme resistance of B. carinata to L. maculans. HighlightsO_LISeven pathogenicity QTL identified, carrying several interesting candidate genes C_LIO_LITransgression of some progeny isolates on B. napus and B. carinata was reported C_LIO_LIEpistasis plays a significant role in the adaptation of L. maculans toward host and nonhost Brassica species C_LI

Genetic rearrangement in the absence of selection is poorly understood. Long-read, high-depth sequencing has enabled analysis of the genetic landscape of spontaneous structural variants. We have applied PacBio Revio sequencing to identify the major classes of genetic alterations during long-term vegetative growth in Aspergillus flavus. These classes include a) precise deletions frequently associated with the flanking region of homology, b) clusters of short deletions that are heavily enriched within centromeric and subtelomeric regions that are typically 94% A+T, and c) deletions and expansions within a subset of homopolymer tracts. We find 70 homopolymer tracts longer than 44 bp, dispersed throughout the genome in noncoding regions and unassociated with transposable elements. A broader analysis of third-generation sequenced datasets from fungi found within the Phylum Ascomycota showed a high but variable prevalence of homopolymeric tracts across species. The adventitious identification of a putative flbA mutation, responsible for an altered conidial phenotype, by the Deep Variant AI system also demonstrates the power of this approach. These studies show that high-depth, long-read sequencing can probe SNPs, indels, homopolymeric tracts, and large-scale rearrangements in an unbiased fashion, and demonstrate that short deletions and homopolymer tract expansions arise spontaneously in A. flavus and other filamentous fungi.

Domestication of molds for cheese production has repeatedly shaped Penicillium fungi, most notably giving rise to the emblematic P. camemberti lineage, derived from P. biforme. Here, we identified a new P. biforme lineage, named cheesy, likely selected from P. biforme for food fermentation, including cheese and sausage production. This lineage exhibits evidence of a severe bottleneck, with little nucleotide polymorphism and a single mating type. The cheesy lineage has evolved advantageous traits for cheesemaking: compared to other P. biforme strains and its wild relative P. fuscoglaucum, it displays faster growth on cheese, galactose and lactose media, higher sporulation and germination rates on cheese, elevated lipolytic activity, enhanced inhibition capacities, and produced specific volatile organic compounds. Additionally, P. camemberti and P. biforme cheesy differ in their content of Starship mobile elements, acquired through horizontal transfers. These elements carry cargo genes potentially relevant for adaptation to cheese. Notably, the cheesy lineage has acquired a 20 kb Starship element (Rattus), nested within a much larger 160 kb Starship (Bilge), and carrying cargo genes with predicted functions involved in antagonistic interactions among micro-organisms. SignificanceO_LIWe identified a new Penicillium biforme lineage, named cheesy, which was likely selected from P. biforme for food fermentation (cheese and sausage), genetically and phenotypically different from the emblematic P. camemberti lineage. This lineage suffered from a severe bottleneck, displaying very little nucleotide polymorphism and a single mating type. C_LIO_LIThis newly identified lineage has evolved advantageous traits for cheesemaking; compared to other P. biforme strains and its wild relative P. fuscoglaucum, it displayed faster growth on cheese, galactose and lactose media, higher sporulation and germination rates on cheese, higher lipolytic activity, better inhibition capacities and specific volatile organic compounds. C_LIO_LIPenicillium camemberti and P. biforme cheesy also differed by their content in Starship mobile elements, acquired through horizontal transfers, and carrying cargo genes potentially relevant for adaptation to cheese. C_LIO_LINotably, the cheesy lineage has acquired a specific 20 kb Starship nested within a much larger 160 kb Starship, and likely involved in antagonistic interactions. C_LI

Transposable elements (TEs) are major contributors to genomic and regulatory variation and can influence the expression of nearby genes in response to environmental stress. The invasive fruit fly Drosophila suzukii, whose genome is composed of nearly 50% TE-derived sequences, provides an excellent model to investigate the regulatory impact of TEs on stress responses. Although previous studies found no global association between TEs and oxidative stress transcriptomic response in this species, the effects of individual TE insertions remain largely unexplored. Here, we investigated the role of two polymorphic intronic TE insertions, Mrp4-TE and inaE-TE, located within genes implicated in the oxidative stress response in D. melanogaster. Using wild-type strains of D. suzukii, we assessed the effect of these insertions on gene expression following oxidative stress induced by paraquat exposure. We found that their corresponding nearby genes Mrp4 and inaE were upregulated after paraquat treatment and that this response was significantly stronger in flies carrying the corresponding TE insertions. Allele-specific expression analyses further showed that the alleles containing the TE insertions exhibited a tendency towards higher expression levels than alleles lacking the insertions, both under control and oxidative stress conditions, supporting a cis-regulatory effect. In addition, Mrp4-TE was associated with increased sensitivity to oxidative stress. Together, our results provide the first evidence that TE insertions can modulate oxidative stress-responsive gene expression in D. suzukii. These findings highlight the importance of investigating individual TE insertions to understand the genetic and regulatory mechanism underlying stress responses, adaptation, and invasion success in this species.

Poly(ADP-ribose) polymerases (PARPs) catalyze ADP-ribosylation, a conserved post- translational modification involved in DNA repair, transcriptional regulation, and chromatin remodeling. Although extensively studied in animals, the evolution and diversification of PARPs across the fungal kingdom remain largely unexplored. Here, we present the first kingdom-wide comparative genomic analysis of PARP proteins across 534 fungal species spanning eight phyla. We identified two primary conserved fungal PARP protein types corresponding to human PARP1 and PARP6. Both exhibit highly dynamic evolutionary histories characterized by frequent independent gain and loss events. Ancestral state reconstruction supports the presence of PARP1 in the last common ancestor of fungi, whereas the PARP6-like family has undergone repeated lineage-specific gains and losses. Fungal PARP6-like proteins retain a compact PARP catalytic domain fused to a C-terminal E2 ubiquitin-conjugating domain, whereas the PARP1 family displays extensive structural diversification through domain shuffling and lineage-specific fusions associated with DNA metabolism, chromatin remodeling, and signal transduction. Reconstruction of ancestral catalytic motifs across fungi and other eukaryotes revealed convergent evolution of a non-canonical H-Y-Y catalytic triad, with multiple motif variants co-occurring within individual proteins, suggesting functional diversification. In the Fusarium oxysporum species complex, we identified a lineage-specific expansion of the PARP family, driven exclusively by accessory chromosomes. Genomes with expanded PARP1 repertoires exhibited elevated basal PARylation, increased resistance to DNA-damaging agents that induce single strand breaks, and DNA damage- induced expression of accessory Parp genes. These findings reveal fungal PARPs as evolutionarily dynamic proteins that likely contribute to genome stability, adaptation, and pathogenicity.