The high clinical mortality and economic burden posed by invasive fungal infections (IFIs), along with significant agricultural crop loss caused by various fungal species, has resulted in the widespread use of antifungal agents. Selective drug pressure, fungal attributes, and host- and drug-related factors have counteracted the efficacy of the limited systemic antifungal drugs and changed the epidemiological landscape of IFIs. Species belonging to Candida, Aspergillus, Cryptococcus, and Pneumocystis are among the fungal pathogens showing notable rates of antifungal resistance. Drug-resistant fungi from the environment are increasingly identified in clinical settings. Furthermore, we have a limited understanding of drug class-specific resistance mechanisms in emerging Candida species. The establishment of antifungal stewardship programs in both clinical and agricultural fields and the inclusion of species identification, antifungal susceptibility testing, and therapeutic drug monitoring practices in the clinic can minimize the emergence of drug-resistant fungi. New antifungal drugs featuring promising therapeutic profiles have great promise to treat drug-resistant fungi in the clinical setting. Mitigating antifungal tolerance, a prelude to the emergence of resistance, also requires the development of effective and fungal-specific adjuvants to be used in combination with systemic antifungals.

Fungi are a large group of eukaryotes found in nearly all ecosystems. More than 250 fungal genomes have already been sequenced, greatly improving our understanding of fungal evolution, physiology, and development. However, for the Pezizomycetes, an early-diverging lineage of filamentous ascomycetes, there is so far only one genome available, namely that of the black truffle, Tuber melanosporum, a mycorrhizal species with unusual subterranean fruiting bodies. To help close the sequence gap among basal filamentous ascomycetes, and to allow conclusions about the evolution of fungal development, we sequenced the genome and assayed transcriptomes during development of Pyronema confluens, a saprobic Pezizomycete with a typical apothecium as fruiting body. With a size of 50 Mb and ~13,400 protein-coding genes, the genome is more characteristic of higher filamentous ascomycetes than the large, repeat-rich truffle genome; however, some typical features are different in the P. confluens lineage, e.g. the genomic environment of the mating type genes that is conserved in higher filamentous ascomycetes, but only partly conserved in P. confluens. On the other hand, P. confluens has a full complement of fungal photoreceptors, and expression studies indicate that light perception might be similar to distantly related ascomycetes and, thus, represent a basic feature of filamentous ascomycetes. Analysis of spliced RNA-seq sequence reads allowed the detection of natural antisense transcripts for 281 genes. The P. confluens genome contains an unusually high number of predicted orphan genes, many of which are upregulated during sexual development, consistent with the idea of rapid evolution of sex-associated genes. Comparative transcriptomics identified the transcription factor gene pro44 that is upregulated during development in P. confluens and the Sordariomycete Sordaria macrospora. The P. confluens pro44 gene (PCON_06721) was used to complement the S. macrospora pro44 deletion mutant, showing functional conservation of this developmental regulator.

Many filamentous ascomycetes develop three-dimensional fruiting bodies for production and dispersal of sexual spores. Fruiting bodies are among the most complex structures differentiated by ascomycetes; however, the molecular mechanisms underlying this process are insufficiently understood. Previous comparative transcriptomics analyses of fruiting body development in different ascomycetes suggested that there might be a core set of genes that are transcriptionally regulated in a similar manner across species. Conserved patterns of gene expression can be indicative of functional relevance, and therefore such a set of genes might constitute promising candidates for functional analyses. In this study, we have sequenced the genome of the Pezizomycete Ascodesmis nigricans, and performed comparative transcriptomics of developing fruiting bodies of this fungus, the Pezizomycete Pyronema confluens, and the Sordariomycete Sordaria macrospora With only 27 Mb, the A. nigricans genome is the smallest Pezizomycete genome sequenced to date. Comparative transcriptomics indicated that gene expression patterns in developing fruiting bodies of the three species are more similar to each other than to nonsexual hyphae of the same species. An analysis of 83 genes that are upregulated only during fruiting body development in all three species revealed 23 genes encoding proteins with predicted roles in vesicle transport, the endomembrane system, or transport across membranes, and 13 genes encoding proteins with predicted roles in chromatin organization or the regulation of gene expression. Among four genes chosen for functional analysis by deletion in S. macrospora, three were shown to be involved in fruiting body formation, including two predicted chromatin modifier genes.

Candida parapsilosis is one of the most commen causes of life-threatening candidaemia, particularly in premature neonates, individuals with cancer of the haematopoietic system, and recipients of organ transplants. Historically, drug-susceptible strains have been linked to clonal outbreaks. However, worldwide studies started since 2018 have reported severe outbreaks among adults caused by fluconazole-resistant strains. Outbreaks caused by fluconazole-resistant strains are associated with high mortality rates and can persist despite strict infection control strategies. The emergence of resistance threatens the efficacy of azoles, which is the most widely used class of antifungals and the only available oral treatment option for candidaemia. The fact that most patients infected with fluconazole-resistant strains are azole-naive underscores the high potential adaptability of fluconazole-resistant strains to diverse hosts, environmental niches, and reservoirs. Another concern is the multidrug-resistant and echinocandin-tolerant C parapsilosis isolates, which emerged in 2020. Raising awareness, establishing effective clinical interventions, and understanding the biology and pathogenesis of fluconazole-resistant C parapsilosis are urgently needed to improve treatment strategies and outcomes.

Many filamentous ascomycetes develop three-dimensional fruiting bodies for production and dispersal of sexual spores. Fruiting bodies are among the most complex structures differentiated by ascomycetes; however, the molecular mechanisms underlying this process are insufficiently understood. Previous comparative transcriptomics analyses of fruiting body development in different ascomycetes suggested that there might be a core set of genes that are transcriptionally regulated in a similar manner across species. Conserved patterns of gene expression can be indicative of functional relevance, and therefore such a set of genes might constitute promising candidates for functional analyses. In this study, we have sequenced the genome of the Pezizomycete Ascodesmis nigricans, and performed comparative transcriptomics of developing fruiting bodies of this fungus, the Pezizomycete Pyronema confluens, and the Sordariomycete Sordaria macrospora With only 27 Mb, the A. nigricans genome is the smallest Pezizomycete genome sequenced to date. Comparative transcriptomics indicated that gene expression patterns in developing fruiting bodies of the three species are more similar to each other than to nonsexual hyphae of the same species. An analysis of 83 genes that are upregulated only during fruiting body development in all three species revealed 23 genes encoding proteins with predicted roles in vesicle transport, the endomembrane system, or transport across membranes, and 13 genes encoding proteins with predicted roles in chromatin organization or the regulation of gene expression. Among four genes chosen for functional analysis by deletion in S. macrospora, three were shown to be involved in fruiting body formation, including two predicted chromatin modifier genes.

Background

Modern civilization depends on only a few plant species for its nourishment. These crops were derived via several thousands of years of human selection that transformed wild ancestors into high-yielding domesticated descendants. Among cultivated plants, common bean (Phaseolus vulgaris L.) is the most important grain legume. Yet, our understanding of the origins and concurrent shaping of the genome of this crop plant is limited.

Results

We sequenced the genomes of 29 accessions representing 12 Phaseolus species. Single nucleotide polymorphism-based phylogenomic analyses, using both the nuclear and chloroplast genomes, allowed us to detect a speciation event, a finding further supported by metabolite profiling. In addition, we identified ~1200 protein coding genes (PCGs) and ~100 long non-coding RNAs with domestication-associated haplotypes. Finally, we describe asymmetric introgression events occurring among common bean subpopulations in Mesoamerica and across hemispheres.

Conclusions

We uncover an unpredicted speciation event in the tropical Andes that gave rise to a sibling species, formerly considered the "wild ancestor" of P. vulgaris, which diverged before the split of the Mesoamerican and Andean P. vulgaris gene pools. Further, we identify haplotypes strongly associated with genes underlying the emergence of domestication traits. Our findings also reveal the capacity of a predominantly autogamous plant to outcross and fix loci from different populations, even from distant species, which led to the acquisition by domesticated beans of adaptive traits from wild relatives. The occurrence of such adaptive introgressions should be exploited to accelerate breeding programs in the near future.

Abstract Background Recent genomic analyses of arthropod defense mechanisms suggest conservation of key elements underlying responses to pathogens, parasites and stresses. At the center of pathogen-induced immune responses are signaling pathways triggered by the recognition of fungal, bacterial and viral signatures. These pathways result in the production of response molecules, such as antimicrobial peptides and lysozymes, which degrade or destroy invaders. Using the recently sequenced genome of the pea aphid (Acyrthosiphon pisum), we conducted the first extensive annotation of the immune and stress gene repertoire of a hemipterous insect, which is phylogenetically distantly related to previously characterized insects models. Results Strikingly, pea aphids appear to be missing genes present in insect genomes characterized to date and thought critical for recognition, signaling and killing of microbes. In line with results of gene annotation, experimental analyses designed to characterize immune response through the isolation of RNA transcripts and proteins from immune-challenged pea aphids uncovered few immune-related products. Gene expression studies, however, indicated some expression of immune and stress-related genes. Conclusions The absence of genes suspected to be essential for the insect immune response suggests that the traditional view of insect immunity may not be as broadly applicable as once thought. The limitations of the aphid immune system may be representative of a broad range of insects, or may be aphid specific. We suggest that several aspects of the aphid life style, such as their association with microbial symbionts, could facilitate survival without strong immune protection.

DNA N6-adenine methylation (6mA) has recently gained importance as an epigenetic modification in eukaryotes. Its function in lineages with high levels, such as early-diverging fungi (EDF), is of particular interest. Here, we investigated the biological significance and evolutionary implications of 6mA in EDF, which exhibit divergent evolutionary patterns in 6mA usage. The analysis of two Mucorales species displaying extreme 6mA usage reveals that species with high 6mA levels show symmetric methylation enriched in highly expressed genes. In contrast, species with low 6mA levels show mostly asymmetric 6mA. Interestingly, transcriptomic regulation throughout development and in response to environmental cues is associated with changes in the 6mA landscape. Furthermore, we identify an EDF-specific methyltransferase, likely originated from endosymbiotic bacteria, as responsible for asymmetric methylation, while an MTA-70 methylation complex performs symmetric methylation. The distinct phenotypes observed in the corresponding mutants reinforced the critical role of both types of 6mA in EDF.

The great cats of the genus Panthera comprise a recent radiation whose evolutionary history is poorly understood. Their rapid diversification poses challenges to resolving their phylogeny while offering opportunities to investigate the historical dynamics of adaptive divergence. We report the sequence, de novo assembly, and annotation of the jaguar (Panthera onca) genome, a novel genome sequence for the leopard (Panthera pardus), and comparative analyses encompassing all living Panthera species. Demographic reconstructions indicated that all of these species have experienced variable episodes of population decline during the Pleistocene, ultimately leading to small effective sizes in present-day genomes. We observed pervasive genealogical discordance across Panthera genomes, caused by both incomplete lineage sorting and complex patterns of historical interspecific hybridization. We identified multiple signatures of species-specific positive selection, affecting genes involved in craniofacial and limb development, protein metabolism, hypoxia, reproduction, pigmentation, and sensory perception. There was remarkable concordance in pathways enriched in genomic segments implicated in interspecies introgression and in positive selection, suggesting that these processes were connected. We tested this hypothesis by developing exome capture probes targeting ~19,000 Panthera genes and applying them to 30 wild-caught jaguars. We found at least two genes (DOCK3 and COL4A5, both related to optic nerve development) bearing significant signatures of interspecies introgression and within-species positive selection. These findings indicate that post-speciation admixture has contributed genetic material that facilitated the adaptive evolution of big cat lineages.

Skates are cartilaginous fish whose novel body plan features remarkably enlarged wing-like pectoral fins that allow them to thrive in benthic environments. The molecular underpinnings of this unique trait, however, remain elusive. Here we investigate the origin of this phenotypic innovation by developing the little skate Leucoraja erinacea as a genomically enabled model. Analysis of a high-quality chromosome-scale genome sequence for the little skate shows that it preserves many ancestral jawed vertebrate features compared with other sequenced genomes, including numerous ancient microchromosomes. Combining genome comparisons with extensive regulatory datasets in developing fins (gene expression, chromatin occupancy and three-dimensional (3D) conformation) we find skate-specific genomic rearrangements that alter the 3D regulatory landscape of genes involved in the planar cell polarity (PCP) pathway. Functional inhibition of PCP signaling resulted in marked reduction of anterior fin size, confirming this pathway as a major contributor of batoid fin morphology. We also identified a fin-specific enhancer that interacts with 3' HOX genes, consistent with the redeployment of Hox gene expression in anterior pectoral fins, and confirmed the potential of this element to activate transcription in the anterior fin using zebrafish reporter assays. Our findings underscore the central role of genome reorganizations and regulatory variation in the evolution of phenotypes, shedding light on the molecular origin of an enigmatic trait.