Stepwise recombination suppression around the mating-type locus in an ascomycete fungus with self-fertile spores.

Recombination is often suppressed at sex-determining loci in plants and animals, and at self-incompatibility or mating-type loci in plants and fungi. In fungal ascomycetes, recombination suppression around the mating-type locus is associated with pseudo-homothallism, i.e. the production of self-fertile dikaryotic sexual spores carrying the two opposite mating types. This has been well studied in two species complexes from different families of Sordariales: Podospora anserina and Neurospora tetrasperma. However, it is unclear whether this intriguing association holds in other species. We show here that Schizothecium tetrasporum, a fungus from a third family in the order Sordariales, also produces mostly self-fertile dikaryotic spores carrying the two opposite mating types. This was due to a high frequency of second meiotic division segregation at the mating-type locus, indicating the occurrence of a single and systematic crossing-over event between the mating-type locus and the centromere, as in P. anserina. The mating-type locus has the typical Sordariales organization, plus a MAT1-1-1 pseudogene in the MAT1-2 haplotype. High-quality genome assemblies of opposite mating types and segregation analyses revealed a suppression of recombination in a region of 1.47 Mb around the mating-type locus. We detected three evolutionary strata, indicating a stepwise extension of recombination suppression. The three strata displayed no rearrangement or transposable element accumulation but gene losses and gene disruptions were present, and precisely at the strata margins. Our findings indicate a convergent evolution of self-fertile dikaryotic sexual spores across multiple ascomycete fungi. The particular pattern of meiotic segregation at the mating-type locus was associated with recombination suppression around this locus, that had extended stepwise. This association between pseudo-homothallism and recombination suppression across lineages and the presence of gene disruption at the strata limits are consistent with a recently proposed mechanism of sheltering deleterious alleles to explain stepwise recombination suppression.

Unraveling coevolutionary dynamics using ecological genomics.

Coevolutionary interactions, from the delicate co-dependency in mutualistic interactions to the antagonistic relationship of hosts and parasites, are a ubiquitous driver of adaptation. Surprisingly, little is known about the genomic processes underlying coevolution in an ecological context. However, species comprise genetically differentiated populations that interact with temporally variable abiotic and biotic environments. We discuss the recent advances in coevolutionary theory and genomics as well as shortcomings, to identify coevolving genes that take into account this spatial and temporal variability of coevolution, and propose a practical guide to understand the dynamic of coevolution using an ecological genomics lens.

Tempo of Degeneration Across Independently Evolved Nonrecombining Regions.

Recombination is beneficial over the long term, allowing more effective selection. Despite long-term advantages of recombination, local recombination suppression can evolve and lead to genomic degeneration, in particular on sex chromosomes. Here, we investigated the tempo of degeneration in nonrecombining regions, that is, the function curve for the accumulation of deleterious mutations over time, leveraging on 22 independent events of recombination suppression identified on mating-type chromosomes of anther-smut fungi, including newly identified ones. Using previously available and newly generated high-quality genome assemblies of alternative mating types of 13 Microbotryum species, we estimated degeneration levels in terms of accumulation of nonoptimal codons and nonsynonymous substitutions in nonrecombining regions. We found a reduced frequency of optimal codons in the nonrecombining regions compared with autosomes, that was not due to less frequent GC-biased gene conversion or lower ancestral expression levels compared with recombining regions. The frequency of optimal codons rapidly decreased following recombination suppression and reached an asymptote after ca. 3 Ma. The strength of purifying selection remained virtually constant at dN/dS = 0.55, that is, at an intermediate level between purifying selection and neutral evolution. Accordingly, nonsynonymous differences between mating-type chromosomes increased linearly with stratum age, at a rate of 0.015 per My. We thus develop a method for disentangling effects of reduced selection efficacy from GC-biased gene conversion in the evolution of codon usage and we quantify the tempo of degeneration in nonrecombining regions, which is important for our knowledge on genomic evolution and on the maintenance of regions without recombination.

Convergent recombination cessation between mating-type genes and centromeres in selfing anther-smut fungi.

The degree of selfing has major impacts on adaptability and is often controlled by molecular mechanisms determining mating compatibility. Changes in compatibility systems are therefore important evolutionary events, but their underlying genomic mechanisms are often poorly understood. Fungi display frequent shifts in compatibility systems, and their small genomes facilitate elucidation of the mechanisms involved. In particular, linkage between the pre- and postmating compatibility loci has evolved repeatedly, increasing the odds of gamete compatibility under selfing. Here, we studied the mating-type chromosomes of two anther-smut fungi with unlinked mating-type loci despite a self-fertilization mating system. Segregation analyses and comparisons of high-quality genome assemblies revealed that these two species displayed linkage between mating-type loci and their respective centromeres. This arrangement renders the same improved odds of gamete compatibility as direct linkage of the two mating-type loci under the automictic mating (intratetrad selfing) of anther-smut fungi. Recombination cessation was found associated with a large inversion in only one of the four linkage events. The lack of trans-specific polymorphism at genes located in nonrecombining regions and linkage date estimates indicated that the events of recombination cessation occurred independently in the two sister species. Our study shows that natural selection can repeatedly lead to similar genomic patterns and phenotypes, and that different evolutionary paths can lead to distinct yet equally beneficial responses to selection. Our study further highlights that automixis and gene linkage to centromeres have important genetic and evolutionary consequences, while being poorly recognized despite being present in a broad range of taxa.

Onset and stepwise extensions of recombination suppression are common in mating-type chromosomes of Microbotryum anther-smut fungi.

Sex chromosomes and mating-type chromosomes can display large genomic regions without recombination. Recombination suppression often extended stepwise with time away from the sex- or mating-type-determining genes, generating evolutionary strata of differentiation between alternative sex or mating-type chromosomes. In anther-smut fungi of the Microbotryum genus, recombination suppression evolved repeatedly, linking the two mating-type loci and extended multiple times in regions distal to the mating-type genes. Here, we obtained high-quality genome assemblies of alternative mating types for four Microbotryum fungi. We found an additional event of independent chromosomal rearrangements bringing the two mating-type loci on the same chromosome followed by recombination suppression linking them. We also found, in a new clade analysed here, that recombination suppression between the two mating-type loci occurred in several steps, with first an ancestral recombination suppression between one of the mating-type locus and its centromere; later, completion of recombination suppression up to the second mating-type locus occurred independently in three species. The estimated dates of recombination suppression between the mating-type loci ranged from 0.15 to 3.58 million years ago. In total, this makes at least nine independent events of linkage between the mating-type loci across the Microbotryum genus. Several mating-type locus linkage events occurred through the same types of chromosomal rearrangements, where similar chromosome fissions at centromeres represent convergence in the genomic changes leading to the phenotypic convergence. These findings further highlight Microbotryum fungi as excellent models to study the evolution of recombination suppression.

Higher Gene Flow in Sex-Related Chromosomes than in Autosomes during Fungal Divergence.

Nonrecombining sex chromosomes are widely found to be more differentiated than autosomes among closely related species, due to smaller effective population size and/or to a disproportionally large-X effect in reproductive isolation. Although fungal mating-type chromosomes can also display large nonrecombining regions, their levels of differentiation compared with autosomes have been little studied. Anther-smut fungi from the Microbotryum genus are castrating pathogens of Caryophyllaceae plants with largely nonrecombining mating-type chromosomes. Using whole genome sequences of 40 fungal strains, we quantified genetic differentiation among strains isolated from the geographically overlapping North American species and subspecies of Silene virginica and S. caroliniana. We inferred that gene flow likely occurred at the early stages of divergence and then completely stopped. We identified large autosomal genomic regions with chromosomal inversions, with higher genetic divergence than the rest of the genomes and highly enriched in selective sweeps, supporting a role of rearrangements in preventing gene flow in genomic regions involved in ecological divergence. Unexpectedly, the nonrecombining mating-type chromosomes showed lower divergence than autosomes due to higher gene flow, which may be promoted by adaptive introgressions of less degenerated mating-type chromosomes. The fact that both mating-type chromosomes are always heterozygous and nonrecombining may explain such patterns that oppose to those found for XY or ZW sex chromosomes. The specific features of mating-type chromosomes may also apply to the UV sex chromosomes determining sexes at the haploid stage in algae and bryophytes and may help test general hypotheses on the evolutionary specificities of sex-related chromosomes.

Independent domestication events in the blue-cheese fungus Penicillium roqueforti.

Domestication provides an excellent framework for studying adaptive divergence. Using population genomics and phenotypic assays, we reconstructed the domestication history of the blue cheese mould Penicillium roqueforti. We showed that this fungus was domesticated twice independently. The population used in Roquefort originated from an old domestication event associated with weak bottlenecks and exhibited traits beneficial for pre-industrial cheese production (slower growth in cheese and greater spore production on bread, the traditional multiplication medium). The other cheese population originated more recently from the selection of a single clonal lineage, was associated with all types of blue cheese worldwide except Roquefort, and displayed phenotypes more suited for industrial cheese production (high lipolytic activity, efficient cheese cavity colonization ability and salt tolerance). We detected genomic regions affected by recent positive selection and putative horizontal gene transfers. This study sheds light on the processes of rapid adaptation and raises questions about genetic resource conservation.

Evolutionary strata on young mating-type chromosomes despite the lack of sexual antagonism.

Sex chromosomes can display successive steps of recombination suppression known as "evolutionary strata," which are thought to result from the successive linkage of sexually antagonistic genes to sex-determining genes. However, there is little evidence to support this explanation. Here we investigate whether evolutionary strata can evolve without sexual antagonism using fungi that display suppressed recombination extending beyond loci determining mating compatibility despite lack of male/female roles associated with their mating types. By comparing full-length chromosome assemblies from five anther-smut fungi with or without recombination suppression in their mating-type chromosomes, we inferred the ancestral gene order and derived chromosomal arrangements in this group. This approach shed light on the chromosomal fusion underlying the linkage of mating-type loci in fungi and provided evidence for multiple clearly resolved evolutionary strata over a range of ages (0.9-2.1 million years) in mating-type chromosomes. Several evolutionary strata did not include genes involved in mating-type determination. The existence of strata devoid of mating-type genes, despite the lack of sexual antagonism, calls for a unified theory of sex-related chromosome evolution, incorporating, for example, the influence of partially linked deleterious mutations and the maintenance of neutral rearrangement polymorphism due to balancing selection on sexes and mating types.

A genome scan of diversifying selection in Ophiocordyceps zombie-ant fungi suggests a role for enterotoxins in co-evolution and host specificity.

Identification of the genes underlying adaptation sheds light on the biological functions targeted by natural selection. Searches for footprints of positive selection, in the form of rapid amino acid substitutions, and the identification of species-specific genes have proved to be powerful approaches to identifying the genes involved in host specialization in plant-pathogenic fungi. We used an evolutionary comparative genomic approach to identify genes underlying host adaptation in the ant-infecting genus Ophiocordyceps, which manipulates ant behaviour. A comparison of the predicted genes in the genomes of species from three species complexes-O. unilateralis, O. australis and O. subramanianii-revealed an enrichment in pathogenesis-associated functions, including heat-labile enterotoxins, among species-specific genes. Furthermore, these genes were overrepresented among those displaying significant footprints of positive selection. Other categories of genes suspected to be important for virulence and pathogenicity in entomopathogenic fungi (e.g., chitinases, lipases, proteases, core secondary metabolism genes) were much less represented, although a few candidate genes were found to evolve under positive selection. An analysis including orthologs from other entomopathogenic fungi in a broader context showed that positive selection on enterotoxins was specific to the ant-infecting genus Ophiocordyceps. Together with previous studies reporting the overexpression of an enterotoxin during behavioural manipulation in diseased ants, our findings suggest that heat-labile enterotoxins are important effectors in host adaptation and co-evolution in the Ophiocordyceps entomopathogenic fungi.

Fungal evolutionary genomics provides insight into the mechanisms of adaptive divergence in eukaryotes.

Fungi are ideal model organisms for dissecting the genomic bases of adaptive divergence in eukaryotes. They have simple morphologies and small genomes, occupy contrasting, well-identified ecological niches and tend to have short generation times, and many are amenable to experimental approaches. Fungi also display diverse lifestyles, from saprotrophs to pathogens or mutualists, and they play extremely important roles in both ecosystems and human activities, as wood decayers, mycorrhizal fungi, lichens, endophytes, plant and animal pathogens, and in fermentation or drug production. We review here recent insights into the patterns and mechanisms of adaptive divergence in fungi, including sources of divergence, genomic variation and, ultimately, speciation. We outline the various ecological sources of divergent selection and genomic changes, showing that gene loss and changes in gene expression and in genomic architecture are important adaptation processes, in addition to the more widely recognized processes of amino acid substitution and gene duplication. We also review recent findings regarding the interspecific acquisition of genomic variation and suggesting an important role for introgression, hybridization and horizontal gene transfers (HGTs). We show that transposable elements can mediate several of these genomic changes, thus constituting important factors for adaptation. Finally, we review the consequences of divergent selection in terms of speciation, arguing that genetic incompatibilities may not be as widespread as generally thought and that pleiotropy between adaptation and reproductive isolation is an important route of speciation in fungal pathogens.

Repeated loss of function at HD mating-type genes and of recombination suppression without mating-type locus linkage in anther-smut fungi.

A wide diversity of mating systems occur in nature, with frequent evolutionary transitions in mating-compatibility mechanisms. Basidiomycete fungi typically have two mating-type loci controlling mating compatibility, HD and PR, usually residing on different chromosomes. In Microbotryum anther-smut fungi, there have been repeated events of linkage between the two mating-type loci through chromosome fusions, leading to large non-recombining regions. By generating high-quality genome assemblies, we found that two sister Microbotryum species parasitizing Dianthus plants, M. superbum and M. shykoffianum, as well as the distantly related M. scorzonarae, have their HD and PR mating-type loci on different chromosomes, but with the PR mating-type chromosome fused with part of the ancestral HD chromosome. Furthermore, progressive extensions of recombination suppression have generated evolutionary strata. In all three species, rearrangements suggest the existence of a transient stage of HD-PR linkage by whole chromosome fusion, and, unexpectedly, the HD genes lost their function. In M. superbum, multiple natural diploid strains were homozygous, and the disrupted HD2 gene was hardly expressed. Mating tests confirmed that a single genetic factor controlled mating compatibility (i.e. PR) and that haploid strains with identical HD alleles could mate and produce infectious hyphae. The HD genes have therefore lost their function in the control of mating compatibility in these Microbotryum species. While the loss of function of PR genes in mating compatibility has been reported in a few basidiomycete fungi, these are the first documented cases for the loss of mating-type determination by HD genes in heterothallic fungi. The control of mating compatibility by a single genetic factor is beneficial under selfing and can thus be achieved repeatedly, through evolutionary convergence in distant lineages, involving different genomic or similar pathways.

New tricks of an old pattern: structural versatility of scorpion toxins with common cysteine spacing.

Scorpion venoms are a rich source of K(+) channel-blocking peptides. For the most part, they are structurally related small disulfide-rich proteins containing a conserved pattern of six cysteines that is assumed to dictate their common three-dimensional folding. In the conventional pattern, two disulfide bridges connect an α-helical segment to the C-terminal strand of a double- or triple-stranded ß-sheet, conforming a cystine-stabilized α/ß scaffold (CSα/ß). Here we show that two K(+) channel-blocking peptides from Tityus scorpions conserve the cysteine spacing of common scorpion venom peptides but display an unconventional disulfide pattern, accompanied by a complete rearrangement of the secondary structure topology into a CS helix-loop-helix fold. Sequence and structural comparisons of the peptides adopting this novel fold suggest that it would be a new elaboration of the widespread CSα/ß scaffold, thus revealing an unexpected structural versatility of these small disulfide-rich proteins. Acknowledgment of such versatility is important to understand how venom structural complexity emerged on a limited number of molecular scaffolds.

Dynamics of transposable element accumulation in the non-recombining regions of mating-type chromosomes in anther-smut fungi.

In the absence of recombination, the number of transposable elements (TEs) increases due to less efficient selection, but the dynamics of such TE accumulations are not well characterized. Leveraging a dataset of 21 independent events of recombination cessation of different ages in mating-type chromosomes of Microbotryum fungi, we show that TEs rapidly accumulated in regions lacking recombination, but that TE content reached a plateau at ca. 50% of occupied base pairs by 1.5 million years following recombination suppression. The same TE superfamilies have expanded in independently evolved non-recombining regions, in particular rolling-circle replication elements (Helitrons). Long-terminal repeat (LTR) retrotransposons of the Copia and Ty3 superfamilies also expanded, through transposition bursts (distinguished from gene conversion based on LTR divergence), with both non-recombining regions and autosomes affected, suggesting that non-recombining regions constitute TE reservoirs. This study improves our knowledge of genome evolution by showing that TEs can accumulate through bursts, following non-linear decelerating dynamics.

Domestication in dry-cured meat Penicillium fungi: Convergent specific phenotypes and horizontal gene transfers without strong genetic subdivision.

Some fungi have been domesticated for food production, with genetic differentiation between populations from food and wild environments, and food populations often acquiring beneficial traits through horizontal gene transfers (HGTs). Studying their adaptation to human-made substrates is of fundamental and applied importance for understanding adaptation processes and for further strain improvement. We studied here the population structures and phenotypes of two distantly related Penicillium species used for dry-cured meat production, P. nalgiovense, the most common species in the dry-cured meat food industry, and P. salamii, used locally by farms. Both species displayed low genetic diversity, lacking differentiation between strains isolated from dry-cured meat and those from other environments. Nevertheless, the strains collected from dry-cured meat within each species displayed slower proteolysis and lipolysis than their wild conspecifics, and those of P. nalgiovense were whiter. Phenotypically, the non-dry-cured meat strains were more similar to their sister species than to their conspecific dry-cured meat strains, indicating an evolution of specific phenotypes in dry-cured meat strains. A comparison of available Penicillium genomes from various environments revealed HGTs, particularly between P. nalgiovense and P. salamii (representing almost 1.5 Mb of cumulative length). HGTs additionally involved P. biforme, also found in dry-cured meat products. We further detected positive selection based on amino acid changes. Our findings suggest that selection by humans has shaped the P. salamii and P. nalgiovense populations used for dry-cured meat production, which constitutes domestication. Several genetic and phenotypic changes were similar in P. salamii, P. nalgiovense and P. biforme, indicating convergent adaptation to the same human-made environment. Our findings have implications for fundamental knowledge on adaptation and for the food industry: the discovery of different phenotypes and of two mating types paves the way for strain improvement by conventional breeding, to elucidate the genomic bases of beneficial phenotypes and to generate diversity.

Structure, function, and chemical synthesis of Vaejovis mexicanus peptide 24: a novel potent blocker of Kv1.3 potassium channels of human T lymphocytes.

Animal venoms are rich sources of ligands for studying ion channels and other pharmacological targets. Proteomic analyses of the soluble venom from the Mexican scorpion Vaejovis mexicanus smithi showed that it contains more than 200 different components. Among them, a 36-residue peptide with a molecular mass of 3864 Da (named Vm24) was shown to be a potent blocker of Kv1.3 of human lymphocytes (K(d) ∼ 3 pM). The three-dimensional solution structure of Vm24 was determined by nuclear magnetic resonance, showing the peptide folds into a distorted cystine-stabilized α/ß motif consisting of a single-turn α-helix and a three-stranded antiparallel ß-sheet, stabilized by four disulfide bridges. The disulfide pairs are formed between Cys6 and Cys26, Cys12 and Cys31, Cys16 and Cys33, and Cys21 and Cys36. Sequence analyses identified Vm24 as the first example of a new subfamily of α-type K(+) channel blockers (systematic number α-KTx 23.1). Comparison with other Kv1.3 blockers isolated from scorpions suggests a number of structural features that could explain the remarkable affinity and specificity of Vm24 toward Kv1.3 channels of lymphocytes.

Vm24, a natural immunosuppressive peptide, potently and selectively blocks Kv1.3 potassium channels of human T cells.

Blockade of Kv1.3 K(+) channels in T cells is a promising therapeutic approach for the treatment of autoimmune diseases such as multiple sclerosis and type 1 diabetes mellitus. Vm24 (α-KTx 23.1) is a novel 36-residue Kv1.3-specific peptide isolated from the venom of the scorpion Vaejovis mexicanus smithi. Vm24 inhibits Kv1.3 channels of human lymphocytes with high affinity (K(d) = 2.9 pM) and exhibits >1500-fold selectivity over other ion channels assayed. It inhibits the proliferation and Ca(2+) signaling of human T cells in vitro and reduces delayed-type hypersensitivity reactions in rats in vivo. Our results indicate that Vm24 has exceptional pharmacological properties that make it an excellent candidate for treatment of certain autoimmune diseases.

Molecular evidence for an Asian origin of monitor lizards followed by Tertiary dispersals to Africa and Australasia.

Monitor lizards are emblematic reptiles that are widely distributed in the Old World. Although relatively well studied in vertebrate research, their biogeographic history is still controversial. We constructed a molecular dataset for 54 anguimorph species, including representatives of all families with detailed sampling of the Varanidae (38 species). Our results are consistent with an Asian origin of the Varanidae followed by a dispersal to Africa 41 (49-33) Ma, possibly via an Iranian route. Another major event was the dispersal of monitors to Australia in the Late Eocene-Oligocene 32 (39-26) Ma. This divergence estimate adds to the suggestion that Australia was colonized by several squamate lineages prior to the collision of the Australian plate with the Asian plate starting 25 Ma.

Host diversification is concurrent with linear motif evolution in a Mastadenovirus hub protein.

Over one hundred Mastadenovirus types infect seven orders of mammals. Virus-host coevolution may involve cospeciation, duplication, host switch and partial extinction events. We reconstruct Mastadenovirus diversification, finding that while cospeciation is dominant, the other three events are also common in Mastadenovirus evolution. Linear motifs are fast-evolving protein functional elements and key mediators of virus-host interactions, thus likely to partake in adaptive viral evolution. We study the evolution of eleven linear motifs in the Mastadenovirus E1A protein, a hub of virus-host protein-protein interactions, in the context of host diversification. The reconstruction of linear motif gain and loss events shows fast linear motif turnover, corresponding a virus-host protein-protein interaction turnover orders of magnitude faster than in model host proteomes. Evolution of E1A linear motifs is coupled, indicating functional coordination at the protein scale, yet presents motif-specific patterns suggestive of convergent evolution. We report a pervasive association between Mastadenovirus host diversification events and the evolution of E1A linear motifs. Eight of 17 host switches associate with the gain of one linear motif and the loss of four different linear motifs, while five of nine partial extinctions associate with the loss of one linear motif. The specific changes in E1A linear motifs during a host switch or a partial extinction suggest that changes in the host molecular environment lead to modulation of the interactions with the retinoblastoma protein and host transcriptional regulators. Altogether, changes in the linear motif repertoire of a viral hub protein are associated with adaptive evolution events during Mastadenovirus evolution.