Genomic Patterns of Positive Selection at the Origin of Rust Fungi

On the effects of selection and mutation on species tree inference

The effect of selection acting on regions of the genome on the accuracy of species-level phylogenetic inference using methods that do not explicitly model selection is an open question that is relevant to most, if not all, phylogenomic studies. To address this, we derive a mathematical approximation to the Wright-Fisher model with mutation and selection in the limit as the population size becomes large. In contrast to previous approximations based on diffusion processes, our approximation can be used to study the distribution of coalescent times for an arbitrary number of lineages, allowing calculation of the probability distribution of gene genealogies under the coalescent model. We use these calculations to show that direct selection at strengths typically encountered in practice has only a small effect on the distribution of coalescent times, and hence on the distribution of gene trees. This implies that many coalescent-based methods for estimating the species tree topology will be robust to the presence of selection in a subset of the underlying genes. Selection will, however, bias the estimation of speciation times, causing them to underestimate the true speciation times. Our model captures the effects of selection on the genealogies that generate the observed sequence data, but does not model selective pressures that act only on the subsequent sequences or that negatively impact gene tree estimation.

Worldwide Population Structure of the Coffee Rust Fungus Hemileia vastatrix Is Strongly Shaped by Local Adaptation and Breeding History

The devastating disease coffee leaf rust, caused by Hemileia vastatrix, has been a major constraint to worldwide coffee production. Recently, H. vastatrix populations were shown to be structured into three divergent genetic lineages with marked host specialization (C1, C2, and C3). However, there is yet no overall understanding of the population dynamics and adaptation of the most widespread and epidemiological relevant H. vastatrix group (C3). We used restriction site-associated DNA sequencing to generate 13,804 single nucleotide polymorphisms (SNPs) across a worldwide collection of 99 H. vastatrix isolates. Phylogenetic analyses uncovered a well-supported structuring within C3, with three main subgroups (SGs; SGI, SGII, and SGIII), which seem to reflect the historical distribution, breeding, and exchange of coffee cultivars. SGI shows a ladder-like diversification pattern and occurs across all four continents sampled, SGII is mainly restricted to Africa, and SGIII is observed only in Timor, revealing a higher genetic differentiation. Outlier and association tests globally identified 112 SNPs under putative positive selection, which impacted population structure. In particular, 29 overlapping SNPs per se seemed to have an extremely strong effect on H. vastatrix population divergence. We also found exclusive and fixed alleles associated with the SGs supporting local adaptation. Functional annotation revealed that transposable elements may play a role in host adaptation. Our study provides a higher-resolution perspective on the evolutionary history of H. vastatrix on cultivated coffee, showing its strong ability to adapt and the strength of the selective force imposed by coffee hosts, which should be taken into account when designing strategies for pathogen dissemination control and selective breeding.

Exploring the therapeutic potential of modern and ancestral phenylalanine/tyrosine ammonia-lyases as supplementary treatment of hereditary tyrosinemia

Phenylalanine/tyrosine ammonia-lyases (PAL/TALs) have been approved by the FDA for treatment of phenylketonuria and may harbour potential for complementary treatment of hereditary tyrosinemia Type I. Herein, we explore ancestral sequence reconstruction as an enzyme engineering tool to enhance the therapeutic potential of PAL/TALs. We reconstructed putative ancestors from fungi and compared their catalytic activity and stability to two modern fungal PAL/TALs. Surprisingly, most putative ancestors could be expressed as functional tetramers in Escherichia coli and thus retained their ability to oligomerize. All ancestral enzymes displayed increased thermostability compared to both modern enzymes, however, the increase in thermostability was accompanied by a loss in catalytic turnover. One reconstructed ancestral enzyme in particular could be interesting for further drug development, as its ratio of specific activities is more favourable towards tyrosine and it is more thermostable than both modern enzymes. Moreover, long-term stability assessment showed that this variant retained substantially more activity after prolonged incubation at 25 °C and 37 °C, as well as an increased resistance to incubation at 60 °C. Both of these factors are indicative of an extended shelf-life of biopharmaceuticals. We believe that ancestral sequence reconstruction has potential for enhancing the properties of enzyme therapeutics, especially with respect to stability. This work further illustrates that resurrection of putative ancestral oligomeric proteins is feasible and provides insight into the extent of conservation of a functional oligomerization surface area from ancestor to modern enzyme.

The genomic determinants of adaptive evolution in a fungal pathogen

Unravelling the strength, frequency, and distribution of selective variants along the genome as well as the underlying factors shaping this distribution are fundamental goals of evolutionary biology. Antagonistic host-pathogen coevolution is thought to be a major driver of genome evolution between interacting species. While rapid evolution of pathogens has been documented in several model organisms, the genetic mechanisms of their adaptation are still poorly understood and debated, particularly the role of sexual reproduction. Here, we apply a population genomic approach to infer genome-wide patterns of selection among 13 isolates of Zymoseptoria tritici, a fungal pathogen characterized by extremely high genetic diversity, gene density, and recombination rates. We report that the genome of Z. tritici undergoes a high rate of adaptive substitutions, with 44% of nonsynonymous substitutions being adaptive on average. This fraction reaches 68% in so-called effector genes encoding determinants of pathogenicity, and the distribution of fitness effects differs in this class of genes as they undergo adaptive mutations with stronger positive fitness effects, but also more slightly deleterious mutations. Besides the globally high rate of adaptive substitutions, we report a negative relationship between pN/pS and the fine-scale recombination rate and a strong positive correlation between the rate of adaptive nonsynonymous substitutions (ωa) and recombination rate. This result suggests a pervasive role of both background selection and Hill-Robertson interference even in a species with an exceptionally high recombination rate (60 cM/Mb on average). While transposable elements (TEs) have been suggested to contribute to adaptation by creating compartments of fast-evolving genomic regions, we do not find a significant effect of TEs on the rate of adaptive mutations. Overall our study suggests that sexual recombination is a significant driver of genome evolution, even in rapidly evolving organisms subject to recurrent mutations with large positive effects.

Insights into the Genomics of Clownfish Adaptive Radiation: Genetic Basis of the Mutualism with Sea Anemones

Clownfishes are an iconic group of coral reef fishes, especially known for their mutualism with sea anemones. This mutualism is particularly interesting as it likely acted as the key innovation that triggered clownfish adaptive radiation. Indeed, after the acquisition of the mutualism, clownfishes diversified into multiple ecological niches linked with host and habitat use. However, despite the importance of this mutualism, the genetic mechanisms allowing clownfishes to interact with sea anemones are still unclear. Here, we used a comparative genomics and molecular evolutionary analyses to investigate the genetic basis of clownfish mutualism with sea anemones. We assembled and annotated the genome of nine clownfish species and one closely related outgroup. Orthologous genes inferred between these species and additional publicly available teleost genomes resulted in almost 16,000 genes that were tested for positively selected substitutions potentially involved in the adaptation of clownfishes to live in sea anemones. We identified 17 genes with a signal of positive selection at the origin of clownfish radiation. Two of them (Versican core protein and Protein O-GlcNAse) show particularly interesting functions associated with N-acetylated sugars, which are known to be involved in sea anemone discharge of toxins. This study provides the first insights into the genetic mechanisms of clownfish mutualism with sea anemones. Indeed, we identified the first candidate genes likely to be associated with clownfish protection form sea anemones, and thus the evolution of their mutualism. Additionally, the genomic resources acquired represent a valuable resource for further investigation of the genomic basis of clownfish adaptive radiation.

Whole-Genome Resequencing and Pan-Transcriptome Reconstruction Highlight the Impact of Genomic Structural Variation on Secondary Metabolite Gene Clusters in the Grapevine Esca Pathogen Phaeoacremonium minimum

The Ascomycete fungus Phaeoacremonium minimum is one of the primary causal agents of Esca, a widespread and damaging grapevine trunk disease. Variation in virulence among Pm. minimum isolates has been reported, but the underlying genetic basis of the phenotypic variability remains unknown. The goal of this study was to characterize intraspecific genetic diversity and explore its potential impact on virulence functions associated with secondary metabolism, cellular transport, and cell wall decomposition. We generated a chromosome-scale genome assembly, using single molecule real-time sequencing, and resequenced the genomes and transcriptomes of multiple isolates to identify sequence and structural polymorphisms. Numerous insertion and deletion events were found for a total of about 1 Mbp in each isolate. Structural variation in this extremely gene dense genome frequently caused presence/absence polymorphisms of multiple adjacent genes, mostly belonging to biosynthetic clusters associated with secondary metabolism. Because of the observed intraspecific diversity in gene content due to structural variation we concluded that a transcriptome reference developed from a single isolate is insufficient to represent the virulence factor repertoire of the species. We therefore compiled a pan-transcriptome reference of Pm. minimum comprising a non-redundant set of 15,245 protein-coding sequences. Using naturally infected field samples expressing Esca symptoms, we demonstrated that mapping of meta-transcriptomics data on a multi-species reference that included the Pm. minimum pan-transcriptome allows the profiling of an expanded set of virulence factors, including variable genes associated with secondary metabolism and cellular transport.

The coffee leaf rust pathogen Hemileia vastatrix: one and a half centuries around the tropics

TAXONOMY AND HISTORY

Hemileia vastatrix Berk. and Broome (Basidiomycota, Pucciniales) was described in 1869 in eastern Africa and Ceylon as the agent of coffee leaf rust and has spread to all coffee cultivation areas worldwide. Major disease outbreaks in Asia, Africa and America caused and continue to cause severe yield losses, making this the most important disease of Arabica coffee, a cash crop for many tropical and sub-tropical countries.

LIFE CYCLE AND DISEASE SYMPTOMS

Hemileia vastatrix is a hemicyclic fungus with the urediniosporic life cycle as its most important (if not only) source of inoculum. Chlorotic spots are the first macroscopic symptoms, preceding the differentiation of suprastomatal, bouquet-shaped, orange-coloured uredinia. The disease can cause yield losses of up to 35% and have a polyetic epidemiological impact on subsequent years.

DISEASE CONTROL

Although the use of fungicides is one of the preferred immediate control measures, the use of resistant cultivars is considered to be the most effective and durable disease control strategy. The discovery of 'Híbrido de Timor' provided sources of resistance that have been used in several breeding programmes and that have been proven to be effective and durable, as some have been in use for more than 30 years.

GENETIC DIVERSITY AND MOLECULAR PATHOGENICITY

Although exhibiting limited genetic polymorphism, the very large genome of H. vastatrix (c. 797 Mbp) conceals great pathological diversity, with more than 50 physiological races. Gene expression studies have revealed a very precocious activation of signalling pathways and production of putative effectors, suggesting that the plant-fungus dialogue starts as early as at the germ tube stage, and have provided clues for the identification of avr genes.

Comparative transcriptomics of Entelegyne spiders (Araneae, Entelegynae), with emphasis on molecular evolution of orphan genes

Next-generation sequencing technology is rapidly transforming the landscape of evolutionary biology, and has become a cost-effective and efficient means of collecting exome information for non-model organisms. Due to their taxonomic diversity, production of interesting venom and silk proteins, and the relative scarcity of existing genomic resources, spiders in particular are excellent targets for next-generation sequencing (NGS) methods. In this study, the transcriptomes of six entelegyne spider species from three genera (Cicurina travisae, C. vibora, Habronattus signatus, H. ustulatus, Nesticus bishopi, and N. cooperi) were sequenced and de novo assembled. Each assembly was assessed for quality and completeness and functionally annotated using gene ontology information. Approximately 100 transcripts with evidence of homology to venom proteins were discovered. After identifying more than 3,000 putatively orthologous genes across all six taxa, we used comparative analyses to identify 24 instances of positively selected genes. In addition, between ~ 550 and 1,100 unique orphan genes were found in each genus. These unique, uncharacterized genes exhibited elevated rates of amino acid substitution, potentially consistent with lineage-specific adaptive evolution. The data generated for this study represent a valuable resource for future phylogenetic and molecular evolutionary research, and our results provide new insight into the forces driving genome evolution in taxa that span the root of entelegyne spider phylogeny.

Biologically Active Secondary Metabolites from the Fungi

Many Fungi have a well-developed secondary metabolism. The diversity of fungal species and the diversification of biosynthetic gene clusters underscores a nearly limitless potential for metabolic variation and an untapped resource for drug discovery and synthetic biology. Much of the ecological success of the filamentous fungi in colonizing the planet is owed to their ability to deploy their secondary metabolites in concert with their penetrative and absorptive mode of life. Fungal secondary metabolites exhibit biological activities that have been developed into life-saving medicines and agrochemicals. Toxic metabolites, known as mycotoxins, contaminate human and livestock food and indoor environments. Secondary metabolites are determinants of fungal diseases of humans, animals, and plants. Secondary metabolites exhibit a staggering variation in chemical structures and biological activities, yet their biosynthetic pathways share a number of key characteristics. The genes encoding cooperative steps of a biosynthetic pathway tend to be located contiguously on the chromosome in coregulated gene clusters. Advances in genome sequencing, computational tools, and analytical chemistry are enabling the rapid connection of gene clusters with their metabolic products. At least three fungal drug precursors, penicillin K and V, mycophenolic acid, and pleuromutilin, have been produced by synthetic reconstruction and expression of respective gene clusters in heterologous hosts. This review summarizes general aspects of fungal secondary metabolism and recent developments in our understanding of how and why fungi make secondary metabolites, how these molecules are produced, and how their biosynthetic genes are distributed across the Fungi. The breadth of fungal secondary metabolite diversity is highlighted by recent information on the biosynthesis of important fungus-derived metabolites that have contributed to human health and agriculture and that have negatively impacted crops, food distribution, and human environments.