Natural history and evolutionary principles of gene duplication in fungi

Exploring whole-genome duplicate gene retention with complex genetic interaction analysis

Whole-genome duplication has played a central role in the genome evolution of many organisms, including the human genome. Most duplicated genes are eliminated, and factors that influence the retention of persisting duplicates remain poorly understood. We describe a systematic complex genetic interaction analysis with yeast paralogs derived from the whole-genome duplication event. Mapping of digenic interactions for a deletion mutant of each paralog, and of trigenic interactions for the double mutant, provides insight into their roles and a quantitative measure of their functional redundancy. Trigenic interaction analysis distinguishes two classes of paralogs: a more functionally divergent subset and another that retained more functional overlap. Gene feature analysis and modeling suggest that evolutionary trajectories of duplicated genes are dictated by combined functional and structural entanglement factors.

Supplementation of Aspergillus glaucus with gfdB gene encoding a glycerol 3-phosphate dehydrogenase in Aspergillus nidulans

In Aspergillus nidulans, there are two putative glycerol 3-phosphate dehydrogenases encoded by the genes gfdA and gfdB, while the genome of the osmophilic Aspergillus glaucus harbors only the ortholog of the A. nidulans gfdA gene. Our aim was to insert the gfdB gene into the genome of A. glaucus, and we reached this goal with the adaptation of the Agrobacterium tumefaciens-mediated transformation method. We tested the growth of the gfdB-complemented A. glaucus strains on a medium containing 2 mol l^-1 sorbitol in the presence of oxidative stress generating agents such as tert-butyl hydroperoxide, H₂ O₂ , menadione sodium bisulfite, as well as the cell wall integrity stress-inducing agent Congo Red and the heavy metal stress eliciting CdCl₂ . The growth of the complemented strains was significantly higher than that of the wild-type strain on media supplemented with these stress generating agents. The A. nidulans ΔgfdB mutant was also examined under the same conditions and resulted in a considerably lower growth than that of the control strain in all stress exposure experiments. Our results shed light on the fact that the gfdB gene from A. nidulans was also involved in the stress responses of the complemented A. glaucus strains supporting our hypothesis on the antioxidant function of GfdB in the Aspergilli. Nevertheless, the osmotolerant nature of A. glaucus could not be explained by the lack of the gfdB gene in A. glaucus, as we hypothesized earlier.

The lasting after-effects of an ancient polyploidy on the genomes of teleosts

The ancestor of most teleost fishes underwent a whole-genome duplication event three hundred million years ago. Despite its antiquity, the effects of this event are evident both in the structure of teleost genomes and in how the surviving duplicated genes still operate to drive form and function. I inferred a set of shared syntenic regions that survive from the teleost genome duplication (TGD) using eight teleost genomes and the outgroup gar genome (which lacks the TGD). I then phylogenetically modeled the TGD's resolution via shared and independent gene losses and applied a new simulation-based statistical test for the presence of bias toward the preservation of genes from one parental subgenome. On the basis of that test, I argue that the TGD was likely an allopolyploidy. I find that duplicate genes surviving from this duplication in zebrafish are less likely to function in early embryo development than are genes that have returned to single copy at some point in this species' history. The tissues these ohnologs are expressed in, as well as their biological functions, lend support to recent suggestions that the TGD was the source of a morphological innovation in the structure of the teleost retina. Surviving duplicates also appear less likely to be essential than singletons, despite the fact that their single-copy orthologs in mouse are no less essential than other genes.

LysM Proteins Regulate Fungal Development and Contribute to Hyphal Protection and Biocontrol Traits in Clonostachys rosea

Lysin motif (LysM) modules are approximately 50 amino acids long and bind to peptidoglycan, chitin and its derivatives. Certain LysM proteins in plant pathogenic and entomopathogenic fungi are shown to scavenge chitin oligosaccharides and thereby dampen host defense reactions. Other LysM proteins can protect the fungal cell wall against hydrolytic enzymes. In this study, we investigated the biological function of LysM proteins in the mycoparasitic fungus Clonostachys rosea. The C. rosea genome contained three genes coding for LysM-containing proteins and gene expression analysis revealed that lysm1 and lysm2 were induced during mycoparasitic interaction with Fusarium graminearum and during colonization of wheat roots. Lysm1 was suppressed in germinating conidia, while lysm2 was induced during growth in chitin or peptidoglycan-containing medium. Deletion of lysm1 and lysm2 resulted in mutants with increased levels of conidiation and conidial germination, but reduced ability to control plant diseases caused by F. graminearum and Botrytis cinerea. The Δlysm2 strain showed a distinct, accelerated mycelial disintegration phenotype accompanied by reduced biomass production and hyphal protection against hydrolytic enzymes including chitinases, suggesting a role of LYSM2 in hyphal protection against chitinases. The Δlysm2 and Δlysm1Δlysm2 strains displayed reduced ability to colonize wheat roots, while only Δlysm1Δlysm2 failed to suppress expression of the wheat defense response genes PR1 and PR4. Based on our data, we propose a role of LYSM1 as a regulator of fungal development and of LYSM2 in cell wall protection against endogenous hydrolytic enzymes, while both are required to suppress plant defense responses. Our findings expand the understanding of the role of LysM proteins in fungal-fungal interactions and biocontrol.

A systematic pipeline for classifying bacterial operons reveals the evolutionary landscape of biofilm machineries

In bacteria functionally related genes comprising metabolic pathways and protein complexes are frequently encoded in operons and are widely conserved across phylogenetically diverse species. The evolution of these operon-encoded processes is affected by diverse mechanisms such as gene duplication, loss, rearrangement, and horizontal transfer. These mechanisms can result in functional diversification, increasing the potential evolution of novel biological pathways, and enabling pre-existing pathways to adapt to the requirements of particular environments. Despite the fundamental importance that these mechanisms play in bacterial environmental adaptation, a systematic approach for studying the evolution of operon organization is lacking. Herein, we present a novel method to study the evolution of operons based on phylogenetic clustering of operon-encoded protein families and genomic-proximity network visualizations of operon architectures. We applied this approach to study the evolution of the synthase dependent exopolysaccharide (EPS) biosynthetic systems: cellulose, acetylated cellulose, poly-β-1,6-N-acetyl-D-glucosamine (PNAG), Pel, and alginate. These polymers have important roles in biofilm formation, antibiotic tolerance, and as virulence factors in opportunistic pathogens. Our approach revealed the complex evolutionary landscape of EPS machineries, and enabled operons to be classified into evolutionarily distinct lineages. Cellulose operons show phyla-specific operon lineages resulting from gene loss, rearrangement, and the acquisition of accessory loci, and the occurrence of whole-operon duplications arising through horizonal gene transfer. Our evolution-based classification also distinguishes between PNAG production from Gram-negative and Gram-positive bacteria on the basis of structural and functional evolution of the acetylation modification domains shared by PgaB and IcaB loci, respectively. We also predict several pel-like operon lineages in Gram-positive bacteria and demonstrate in our companion paper (Whitfield et al PLoS Pathogens, in press) that Bacillus cereus produces a Pel-dependent biofilm that is regulated by cyclic-3',5'-dimeric guanosine monophosphate (c-di-GMP).

The Evolution of Gene Duplicates in Angiosperms and the Impact of Protein-Protein Interactions and the Mechanism of Duplication

Gene duplicates, generated through either whole genome duplication (WGD) or small-scale duplication (SSD), are prominent in angiosperms and are believed to play an important role in adaptation and in generating evolutionary novelty. Previous studies reported contrasting evolutionary and functional dynamics of duplicate genes depending on the mechanism of origin, a behavior that is hypothesized to stem from constraints to maintain the relative dosage balance between the genes concerned and their interaction context. However, the mechanisms ultimately influencing loss and retention of gene duplicates over evolutionary time are not yet fully elucidated. Here, by using a robust classification of gene duplicates in Arabidopsis thaliana, Solanum lycopersicum, and Zea mays, large RNAseq expression compendia and an extensive protein-protein interaction (PPI) network from Arabidopsis, we investigated the impact of PPIs on the differential evolutionary and functional fate of WGD and SSD duplicates. In all three species, retained WGD duplicates show stronger constraints to diverge at the sequence and expression level than SSD ones, a pattern that is also observed for shared PPI partners between Arabidopsis duplicates. PPIs are preferentially distributed among WGD duplicates and specific functional categories. Furthermore, duplicates with PPIs tend to be under stronger constraints to evolve than their counterparts without PPIs regardless of their mechanism of origin. Our results support dosage balance constraint as a specific property of genes involved in biological interactions, including physical PPIs, and suggest that additional factors may be differently influencing the evolution of genes following duplication, depending on the species, time, and mechanism of origin.

Evolutionary Histories of Type III Polyketide Synthases in Fungi

Type III polyketide synthases (PKSs) produce secondary metabolites with diverse biological activities, including antimicrobials. While they have been extensively studied in plants and bacteria, only a handful of type III PKSs from fungi has been characterized in the last 15 years. The exploitation of fungal type III PKSs to produce novel bioactive compounds requires understanding the diversity of these enzymes, as well as of their biosynthetic pathways. Here, phylogenetic and reconciliation analyses of 522 type III PKSs from 1,193 fungal genomes revealed complex evolutionary histories with massive gene duplications and losses, explaining their discontinuous distribution in the fungal tree of life. In addition, horizontal gene transfer events from bacteria to fungi and, to a lower extent, between fungi, could be inferred. Ancestral gene duplication events have resulted in the divergence of eight phylogenetic clades. Especially, two clades show ancestral linkage and functional co-evolution between a type III PKS and a reducing PKS genes. Investigation of the occurrence of protein domains in fungal type III PKS predicted gene clusters highlighted the diversity of biosynthetic pathways, likely reflecting a large chemical landscape. Type III PKS genes are most often located next to genes encoding cytochrome P450s, MFS transporters and transcription factors, defining ancestral core gene clusters. This analysis also allowed predicting gene clusters for the characterized fungal type III PKSs and provides working hypotheses for the elucidation of the full biosynthetic pathways. Altogether, our analyses provide the fundamental knowledge to motivate further characterization and exploitation of fungal type III PKS biosynthetic pathways.

A Novel cis Element Achieves the Same Solution as an Ancestral cis Element During Thiamine Starvation in Candida glabrata

Regulatory networks often converge on very similar cis sequences to drive transcriptional programs due to constraints on what transcription factors are present. To determine the role of constraint loss on cis element evolution, we examined the recent appearance of a thiamine starvation regulated promoter in Candida glabrata. This species lacks the ancestral transcription factor Thi2, but still has the transcription factor Pdc2, which regulates thiamine starvation genes, allowing us to determine the effect of constraint change on a new promoter. We identified two different cis elements in C. glabrata - one present in the evolutionarily recent gene called CgPMU3, and the other element present in the other thiamine (THI) regulated genes. Reciprocal swaps of the cis elements and incorporation of the S. cerevisiaeThi2 transcription factor-binding site into these promoters demonstrate that the two elements are functionally different from one another. Thus, this loss of an imposed constraint on promoter function has generated a novel cis sequence, suggesting that loss of trans constraints can generate a non-convergent pathway with the same output.

Location, location! cellular relocalization primes specialized metabolic diversification

Specialized metabolites are structurally diverse and cell- or tissue-specific molecules produced in restricted plant lineages. In contrast, primary metabolic pathways are highly conserved in plants and produce metabolites essential for all of life, such as amino acids and nucleotides. Substrate promiscuity - the capacity to accept non-native substrates - is a common characteristic of enzymes, and its impact is especially apparent in generating specialized metabolite variation. However, promiscuity only leads to metabolic diversity when alternative substrates are available; thus, enzyme cellular and subcellular localization directly influence chemical phenotypes. We review a variety of mechanisms that modulate substrate availability for promiscuous plant enzymes. We focus on examples where evolution led to modification of the 'cellular context' through changes in cell-type expression, subcellular relocalization, pathway sequestration, and cellular mixing via tissue damage. These varied mechanisms contributed to the emergence of structurally diverse plant specialized metabolites and inform future metabolic engineering approaches.

Interspecific hybridization facilitates niche adaptation in beer yeast

Hybridization between species often leads to non-viable or infertile offspring, yet examples of evolutionarily successful interspecific hybrids have been reported in all kingdoms of life. However, many questions on the ecological circumstances and evolutionary aftermath of interspecific hybridization remain unanswered. In this study, we sequenced and phenotyped a large set of interspecific yeast hybrids isolated from brewing environments to uncover the influence of interspecific hybridization in yeast adaptation and domestication. Our analyses demonstrate that several hybrids between Saccharomyces species originated and diversified in industrial environments by combining key traits of each parental species. Furthermore, posthybridization evolution within each hybrid lineage reflects subspecialization and adaptation to specific beer styles, a process that was accompanied by extensive chimerization between subgenomes. Our results reveal how interspecific hybridization provides an important evolutionary route that allows swift adaptation to novel environments.

Inferring biosynthetic and gene regulatory networks from Artemisia annua RNA sequencing data on a credit card-sized ARM computer

Prediction of gene function and gene regulatory networks is one of the most active topics in bioinformatics. The accumulation of publicly available gene expression data for hundreds of plant species, together with advances in bioinformatical methods and affordable computing, sets ingenuity as one of the major bottlenecks in understanding gene function and regulation. Here, we show how a credit card-sized computer retailing for <50 USD can be used to rapidly predict gene function and infer regulatory networks from RNA sequencing data. To achieve this, we constructed a bioinformatical pipeline that downloads and allows quality-control of RNA sequencing data; and generates a gene co-expression network that can reveal enzymes and transcription factors participating and controlling a given biosynthetic pathway. We exemplify this by first identifying genes and transcription factors involved in the biosynthesis of secondary cell wall in the plant Artemisia annua, the main natural source of the anti-malarial drug artemisinin. Networks were then used to dissect the artemisinin biosynthesis pathway, which suggest potential transcription factors regulating artemisinin biosynthesis. We provide the source code of our pipeline (https://github.com/mutwil/LSTrAP-Lite) and envision that the ubiquity of affordable computing, availability of biological data and increased bioinformatical training of biologists will transform the field of bioinformatics. This article is part of a Special Issue entitled: Transcriptional Profiles and Regulatory Gene Networks edited by Dr. Dr. Federico Manuel Giorgi and Dr. Shaun Mahony.

Gene Repositioning Is Under Constraints After Evolutionary Conserved Gene Neighborhood Separate

Genes are not randomly distributed on eukaryotic chromosomes. Some neighboring genes show order conservation among species, while some neighboring genes separate during evolution even though their neighborhoods are conserved in some species. Here, I investigated whether after-separation gene repositioning is under natural selection for evolutionary conserved gene neighborhoods compared with nonconserved neighborhoods. After separation, genes with conserved neighborhoods show low-expression divergence between the after-separation species and the before-separation species. After genes separate from their conserved gene neighbors, their after-separation gene neighbors tend to show coexpression and coprotein complex with their before-separation gene neighbors. These results indicate evolutionary constraints on the selection of neighboring genes after evolutionary conserved gene neighborhoods separate.

Selective deployment of transcription factor paralogs with submaximal strength facilitates gene regulation in the immune system

In multicellular organisms, duplicated genes can diverge through tissue-specific gene expression patterns, as exemplified by highly regulated expression of Runx transcription factor paralogs with apparent functional redundancy. Here we asked what cell type-specific biologies might be supported by the selective expression of Runx paralogs during Langerhans cell and inducible regulatory T cell differentiation. We uncovered functional non-equivalence between Runx paralogs. Selective expression of native paralogs allowed integration of transcription factor activity with extrinsic signals, while non-native paralogs enforced differentiation even in the absence of exogenous inducers. DNA-binding affinity was controlled by divergent amino acids within the otherwise highly conserved RUNT domain, and evolutionary reconstruction suggested convergence of RUNT domain residues towards sub-maximal strength. Hence, the selective expression of gene duplicates in specialized cell types can synergize with the acquisition of functional differences to enable appropriate gene expression, lineage choice and differentiation in the mammalian immune system.

Genome-Wide Identification and Expression Analyses of the bZIP Transcription Factor Genes in moso bamboo (Phyllostachys edulis)

The basic leucine zipper (bZIP) transcription factor (TF) family is one of the largest gene families, and play crucial roles in many processes, including stress responses, hormone effects. The TF family also participates in plant growth and development. However, limited information is available for these genes in moso bamboo (Phyllostachys edulis), one of the most important non-timber forest products in the world. In the present study, 154 putative PhebZIP genes were identified in the moso bamboo genome. The phylogenetic analyses indicate that the PhebZIP gene proteins classify into 9 subfamilies and the gene structures and conserved motifs that analyses identified among all PhebZIP proteins suggested a high group-specificity. Microsynteny and evolutionary patterns analyses of the non-synonymous (Ka) and synonymous (Ks) substitution rates and their ratios indicated that paralogous pairs of PhebZIP genes in moso bamboo underwent a large-scale genome duplication event that occurred 7-15 million years ago (MYA). According to promoter sequence analysis, we further selected 18 genes which contain the higher number of cis-regulatory elements for expression analysis. The result showed that these genes are extensively involved in GA-, ABA- and MeJA-responses, with possibly different mechanisms. The tissue-specific expression profiles of PhebZIP genes in five plant tissues/organs/developmental stages suggested that these genes are involved in moso bamboo organ development, especially seed development. Subcellular localization and transactivation activity analysis showed that PhebZIP47 and PhebZIP126 were localized in the nucleus and PhebZIP47 with no transcriptional activation in yeast. Our research provides a comprehensive understanding of PhebZIP genes and may aid in the selection of appropriate candidate genes for further cloning and functional analysis in moso bamboo growth and development, and improve their resistance to stress during their life.

Invasive invertebrates associated with highly duplicated gene content

Invasion of alien species has led to serious problems, including the destruction of native ecosystems. In general, invasive species adapt to new environments rapidly, suggesting that they have high genetic diversity that can directly influence environmental adaptability. However, it is not known how genomic architecture causes genetic diversity that leads to invasiveness. Recent studies have showed that the proportion of duplicated genes (P_D ) in whole animal genomes correlate with environmental variability within a habitat. Here, we show that P_D and propagule size significantly explain the differences in species categories (invasive species, noninvasive species, and parasites). P_D correlated negatively with the propagule size. The residual values of regression of P_D on propagule size revealed that the invasive species had higher P_D values and larger propagule size than those of the noninvasive species, whereas the parasites had lower P_D values and smaller propagule size than those of others. There were no correlations between the invasive species and other genomic factors including the genome size, number of genes, and certain gene families. Our results suggest that the P_D values of a genome might be a potential genomic source causing genetic variations for adaptation to diverse environments. The results also showed that the invasiveness status of a species would be predicted by the residual values of regression of P_D on propagule size. Our innovative approach provides a measure to estimate the environmental adaptability of organisms based on genomic data.

Evidence for Highly Variable, Region-Specific Patterns of T-Cell Epitope Mutations Accumulating in Mycobacterium tuberculosis Strains

Vaccines that confer protection through induction of adaptive T-cell immunity rely on understanding T-cell epitope (TCE) evolution induced by immune escape. This is poorly understood in tuberculosis (TB), an ancient, chronic disease, where CD4 T-cell immunity is of recognized importance. We probed 905 functionally validated, curated human CD4 T cell epitopes in 79 Mycobacterium tuberculosis (Mtb) whole genomes from India. This screen resulted in identifying 64 mutated epitopes in these strains initially using a computational pipeline and subsequently verified by single nucleotide polymorphism (SNP) analysis. SNP based phylogeny revealed the 79 Mtb strains to cluster to East African Indian (EAI), Central Asian Strain (CAS), and Beijing (BEI) lineages. Eighty-nine percent of the mutated T-cell epitopes (mTCEs) identified in the 79 Mtb strains from India has not previously been reported. These mTCEs were encoded by genes with high nucleotide diversity scores including seven mTCEs encoded by six antigens in the top 10% of rapidly divergent Mtb genes encoded by these strains. Using a T cell functional assay readout, we demonstrate 62% of mTCEs tested to significantly alter CD4 T-cell IFNγ and/or IL2 secretion with associated changes in predicted HLA-DR binding affinity: the gain of function mutations displayed higher predicted HLA-DR binding affinity and conversely mutations resulting in loss of function displayed lower predicted HLA-DR binding affinity. Most mutated antigens belonged to the cell wall/cell processes, and, intermediary metabolism and respiration families though all known Mtb proteins encoded mutations. Analysis of the mTCEs in an SNP database of 5,310 global Mtb strains identified 82% mTCEs to be significantly more prevalent in Mtb strains isolated from India, including 36 mTCEs identified exclusively in strains from India. These epitopes had a significantly higher predicted binding affinity to HLA-DR alleles that were highly prevalent in India compared to HLA-DR alleles rare in India, highlighting HLA-DR maybe an important driver of these mutations. This first evidence of region-specific TCE mutations potentially employed by Mtb to escape host immunity has important implications for TB vaccine design.

Genome Sequencing of Cladobotryum protrusum Provides Insights into the Evolution and Pathogenic Mechanisms of the Cobweb Disease Pathogen on Cultivated Mushroom

Cladobotryum protrusum is one of the mycoparasites that cause cobweb disease on cultivated edible mushrooms. However, the molecular mechanisms of evolution and pathogenesis of C. protrusum on mushrooms are largely unknown. Here, we report a high-quality genome sequence of C. protrusum using the single-molecule, real-time sequencing platform of PacBio and perform a comparative analysis with closely related fungi in the family Hypocreaceae. The C. protrusum genome, the first complete genome to be sequenced in the genus Cladobotryum, is 39.09 Mb long, with an N50 of 4.97 Mb, encoding 11,003 proteins. The phylogenomic analysis confirmed its inclusion in Hypocreaceae, with its evolutionary divergence time estimated to be ~170.1 million years ago. The genome encodes a large and diverse set of genes involved in secreted peptidases, carbohydrate-active enzymes, cytochrome P450 enzymes, pathogen⁻host interactions, mycotoxins, and pigments. Moreover, C. protrusum harbors arrays of genes with the potential to produce bioactive secondary metabolites and stress response-related proteins that are significant for adaptation to hostile environments. Knowledge of the genome will foster a better understanding of the biology of C. protrusum and mycoparasitism in general, as well as help with the development of effective disease control strategies to minimize economic losses from cobweb disease in cultivated edible mushrooms.

Functional genetic discovery of enzymes using full-scan mass spectrometry metabolomics

Our understanding of metabolic networks is incomplete, and new enzymatic activities await discovery in well-studied organisms. Mass spectrometric measurement of cellular metabolites reveals compounds inside cells that are unexplained by current maps of metabolic reactions, and existing computational models are unable to account for all activities observed within cells. Additional large-scale genetic and biochemical approaches are required to elucidate metabolic gene function. We have used full-scan mass spectrometry metabolomics of polar small molecules to examine deletion mutants of candidate enzymes in the model yeast Saccharomyces cerevisiae. We report the identification of 25 genes whose deletion results in focal metabolic changes consistent with loss of enzymatic activity and describe the informatic approaches used to enrich for candidate enzymes from uncharacterized open reading frames. Triumphs and pitfalls of metabolic phenotyping screens are discussed, including estimates of the frequency of uncharacterized eukaryotic genes that affect metabolism and key issues to consider when searching for new enzymatic functions in other organisms.

Genetic Basis of Variation in Heat and Ethanol Tolerance in Saccharomyces cerevisiae

Saccharomyces cerevisiae has the capability of fermenting sugar to produce concentrations of ethanol that are toxic to most organisms. Other Saccharomyces species also have a strong fermentative capacity, but some are specialized to low temperatures, whereas S. cerevisiae is the most thermotolerant. Although S. cerevisiae has been extensively used to study the genetic basis of ethanol tolerance, much less is known about temperature dependent ethanol tolerance. In this study, we examined the genetic basis of ethanol tolerance at high temperature among strains of S. cerevisiae. We identified two amino acid polymorphisms in SEC24 that cause strong sensitivity to ethanol at high temperature and more limited sensitivity to temperature in the absence of ethanol. We also identified a single amino acid polymorphism in PSD1 that causes sensitivity to high temperature in a strain dependent fashion. The genes we identified provide further insight into genetic variation in ethanol and temperature tolerance and the interdependent nature of these two traits in S. cerevisiae.

Diversification of DNA binding specificities enabled SREBP transcription regulators to expand the repertoire of cellular functions that they govern in fungi

The Sterol Regulatory Element Binding Proteins (SREBPs) are basic-helix-loop-helix transcription regulators that control the expression of sterol biosynthesis genes in higher eukaryotes and some fungi. Surprisingly, SREBPs do not regulate sterol biosynthesis in the ascomycete yeasts (Saccharomycotina) as this role was handed off to an unrelated transcription regulator in this clade. The SREBPs, nonetheless, expanded in fungi such as the ascomycete yeasts Candida spp., raising questions about their role and evolution in these organisms. Here we report that the fungal SREBPs diversified their DNA binding preferences concomitantly with an expansion in function. We establish that several branches of fungal SREBPs preferentially bind non-palindromic DNA sequences, in contrast to the palindromic DNA motifs recognized by most basic-helix-loop-helix proteins (including SREBPs) in higher eukaryotes. Reconstruction and biochemical characterization of the likely ancestor protein suggest that an intrinsic DNA binding promiscuity in the family was resolved by alternative mechanisms in different branches of fungal SREBPs. Furthermore, we show that two SREBPs in the human commensal yeast Candida albicans drive a transcriptional cascade that inhibits a morphological switch under anaerobic conditions. Preventing this morphological transition enhances C. albicans colonization of the mammalian intestine, the fungus' natural niche. Thus, our results illustrate how diversification in DNA binding preferences enabled the functional expansion of a family of eukaryotic transcription regulators.

Comparative evolutionary histories of fungal proteases reveal gene gains in the mycoparasitic and nematode-parasitic fungus Clonostachys rosea

Background

The ascomycete fungus Clonostachys rosea (order Hypocreales) can control several important plant diseases caused by plant pathogenic fungi and nematodes. Subtilisin-like serine proteases are considered to play an important role in pathogenesis in entomopathogenic, mycoparasitic, and nematophagous fungi used for biological control. In this study, we analysed the evolutionary histories of protease gene families, and investigated sequence divergence and regulation of serine protease genes in C. rosea.

Results

Proteases of selected hypocrealean fungal species were classified into families based on the MEROPS peptidase database. The highest number of protease genes (590) was found in Fusarium solani, followed by C. rosea with 576 genes. Analysis of gene family evolution identified non-random changes in gene copy numbers in the five serine protease gene families S1A, S8A, S9X, S12 and S33. Four families, S1A, S8A, S9X, and S33, displayed gene gains in C. rosea. A gene-tree / species-tree reconciliation analysis of the S8A family revealed that the gene copy number increase in C. rosea was primarily associated with the S08.054 (proteinase K) subgroup. In addition, regulatory and predicted structural differences, including twelve sites evolving under positive selection, among eighteen C. rosea S8A serine protease paralog genes were also observed. The C. rosea S8A serine protease gene prs6 was induced during interaction with the plant pathogenic species F. graminearum.

Conclusions

Non-random increases in S8A, S9X and S33 serine protease gene numbers in the mycoparasitic species C. rosea, Trichoderma atroviride and T. virens suggests an involvement in fungal-fungal interactions. Regulatory and predicted structural differences between C. rosea S8A paralogs indicate that functional diversification is driving the observed increase in gene copy numbers. The induction of prs6 expression in C. rosea during confrontation with F. graminearum suggests an involvement of the corresponding protease in fungal-fungal interactions. The results pinpoint the importance of serine proteases for ecological niche adaptation in C. rosea, including a potential role in the mycoparasitic attack on fungal prey.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1291-1) contains supplementary material, which is available to authorized users.

SynerClust: a highly scalable, synteny-aware orthologue clustering tool

Accurate orthologue identification is a vital component of bacterial comparative genomic studies, but many popular sequence-similarity-based approaches do not scale well to the large numbers of genomes that are now generated routinely. Furthermore, most approaches do not take gene synteny into account, which is useful information for disentangling paralogues. Here, we present SynerClust, a user-friendly synteny-aware tool based on synergy that can process thousands of genomes. SynerClust was designed to analyse genomes with high levels of local synteny, particularly prokaryotes, which have operon structure. SynerClust’s run-time is optimized by selecting cluster representatives at each node in the phylogeny; thus, avoiding the need for exhaustive pairwise similarity searches. In benchmarking against Roary, Hieranoid2, PanX and Reciprocal Best Hit, SynerClust was able to more completely identify sets of core genes for datasets that included diverse strains, while using substantially less memory, and with scalability comparable to the fastest tools. Due to its scalability, ease of installation and use, and suitability for a variety of computing environments, orthogroup clustering using SynerClust will enable many large-scale prokaryotic comparative genomics efforts.

Systematic gene overexpression in Candida albicans identifies a regulator of early adaptation to the mammalian gut

Candida albicans is part of the human gastrointestinal (GI) microbiota. To better understand how C. albicans efficiently establishes GI colonisation, we competitively challenged growth of 572 signature-tagged strains (~10% genome coverage), each conditionally overexpressing a single gene, in the murine gut. We identified CRZ2, a transcription factor whose overexpression and deletion respectively increased and decreased early GI colonisation. Using clues from genome-wide expression and gene-set enrichment analyses, we found that the optimal activity of Crz2p occurs under hypoxia at 37°C, as evidenced by both phenotypic and transcriptomic analyses following CRZ2 genetic perturbation. Consistent with early colonisation of the GI tract, we show that CRZ2 overexpression confers resistance to acidic pH and bile salts, suggesting an adaptation to the upper sections of the gut. Genome-wide location analyses revealed that Crz2p directly modulates the expression of many mannosyltransferase- and cell-wall protein-encoding genes, suggesting a link with cell-wall function. We show that CRZ2 overexpression alters cell-wall phosphomannan abundance and increases sensitivity to tunicamycin, suggesting a role in protein glycosylation. Our study reflects the powerful use of gene overexpression as a complementary approach to gene deletion to identify relevant biological pathways involved in C. albicans interaction with the host environment.

Patterns of Population Variation in Two Paleopolyploid Eudicot Lineages Suggest That Dosage-Based Selection on Homeologs Is Long-Lived

Genes that are inherently subject to strong selective constraints tend to be overretained in duplicate after polyploidy. They also continue to experience similar, but somewhat relaxed, constraints after that polyploidy event. We sought to assess for how long the influence of polyploidy is felt on these genes’ selective pressures. We analyzed two nested polyploidy events in Brassicaceae: the At-α genome duplication that is the most recent polyploidy in the model plant Arabidopsis thaliana and a more recent hexaploidy shared by the genus Brassica and its relatives. By comparing the strength and direction of the natural selection acting at the population and at the species level, we find evidence for continued intensified purifying selection acting on retained duplicates from both polyploidies even down to the present. The constraint observed in preferentially retained genes is not a result of the polyploidy event: the orthologs of such genes experience even stronger constraint in nonpolyploid outgroup genomes. In both the Arabidopsis and Brassica lineages, we further find evidence for segregating mildly deleterious variants, confirming that the population-level data uncover patterns not visible with between-species comparisons. Using the A. thaliana metabolic network, we also explored whether network position was correlated with the measured selective constraint. At both the population and species level, nodes/genes tended to show similar constraints to their neighbors. Our results paint a picture of the long-lived effects of polyploidy on plant genomes, suggesting that even yesterday’s polyploids still have distinct evolutionary trajectories.

Divergent Expression Patterns and Function of Two cxcr4 Paralogs in Hermaphroditic Epinephelus coioides

Chemokine receptor Cxcr4 evolved two paralogs in the teleost lineage. However, cxcr4a and cxcr4b have been characterized only in a few species. In this study, we identified two cxcr4 paralogs from the orange-spotted grouper, Epinephelus coioides. The phylogenetic relationship and gene structure and synteny suggest that the duplicated cxcr4a/b should result from the teleost-specific genome duplication (Ts3R). The teleost cxcr4 gene clusters in two paralogous chromosomes exhibit a complementary gene loss/retention pattern. Ec_cxcr4a and Ec_cxcr4b show differential and biased expression patterns in grouper adult tissue, gonads, and embryos at different stages. During embryogenesis, Ec_cxcr4a/b are abundantly transcribed from the neurula stage and mainly expressed in the neural plate and sensory organs, indicating their roles in neurogenesis. Ec_Cxcr4a and Ec_Cxcr4b possess different chemotactic migratory abilities from the human SDF-1α, Ec_Cxcl12a, and Ec_Cxcl12b. Moreover, we uncovered the N-terminus and TM5 domain as the key elements for specific ligand⁻receptor recognition of Ec_Cxcr4a-Ec_Cxcl12b and Ec_Cxcr4b-Ec_Cxcl12a. Based on the biased and divergent expression patterns of Eccxcr4a/b, and specific ligand⁻receptor recognition of Ec_Cxcl12a/b⁻Ec_Cxcr4b/a, the current study provides a paradigm of sub-functionalization of two teleost paralogs after Ts3R.

A prior-based integrative framework for functional transcriptional regulatory network inference

Transcriptional regulatory networks specify regulatory proteins controlling the context-specific expression levels of genes. Inference of genome-wide regulatory networks is central to understanding gene regulation, but remains an open challenge. Expression-based network inference is among the most popular methods to infer regulatory networks, however, networks inferred from such methods have low overlap with experimentally derived (e.g. ChIP-chip and transcription factor (TF) knockouts) networks. Currently we have a limited understanding of this discrepancy. To address this gap, we first develop a regulatory network inference algorithm, based on probabilistic graphical models, to integrate expression with auxiliary datasets supporting a regulatory edge. Second, we comprehensively analyze our and other state-of-the-art methods on different expression perturbation datasets. Networks inferred by integrating sequence-specific motifs with expression have substantially greater agreement with experimentally derived networks, while remaining more predictive of expression than motif-based networks. Our analysis suggests natural genetic variation as the most informative perturbation for network inference, and, identifies core TFs whose targets are predictable from expression. Multiple reasons make the identification of targets of other TFs difficult, including network architecture and insufficient variation of TF mRNA level. Finally, we demonstrate the utility of our inference algorithm to infer stress-specific regulatory networks and for regulator prioritization.

Intrinsic cooperativity potentiates parallel cis-regulatory evolution

Convergent evolutionary events in independent lineages provide an opportunity to understand why evolution favors certain outcomes over others. We studied such a case where a large set of genes-those coding for the ribosomal proteins-gained cis-regulatory sequences for a particular transcription regulator (Mcm1) in independent fungal lineages. We present evidence that these gains occurred because Mcm1 shares a mechanism of transcriptional activation with an ancestral regulator of the ribosomal protein genes, Rap1. Specifically, we show that Mcm1 and Rap1 have the inherent ability to cooperatively activate transcription through contacts with the general transcription factor TFIID. Because the two regulatory proteins share a common interaction partner, the presence of one ancestral cis-regulatory sequence can 'channel' random mutations into functional sites for the second regulator. At a genomic scale, this type of intrinsic cooperativity can account for a pattern of parallel evolution involving the fixation of hundreds of substitutions.

Duplications and losses of genes encoding known elements of the stress defence system of the Aspergilli contribute to the evolution of these filamentous fungi but do not directly influence their environmental stress tolerance

The contribution of stress protein duplication and deletion events to the evolution of the Aspergilli was studied. We performed a large-scale homology analysis of stress proteins and generated and analysed three stress defence system models based on Saccharomyces cerevisiae, Schizosaccharomyces pombe and Aspergillus nidulans. Although both yeast-based and A. nidulans-based models were suitable to trace evolutionary changes, the A. nidulans-based model performed better in mapping stress protein radiations. The strong Mantel correlation found between the positions of species in the phylogenetic tree on the one hand and either in the A. nidulans-based or S. cerevisiae-based models on the other hand demonstrated that stress protein expansions and reductions contributed significantly to the evolution of the Aspergilli. Interestingly, stress tolerance attributes correlated well with the number of orthologs only for a few stress proteins. Notable examples are Ftr1 iron permease and Fet3 ferro-O2-oxidoreductase, elements of the reductive iron assimilation pathway, in the S. cerevisiae-based model, as well as MpkC, a HogA-like mitogen activated protein kinase in the A. nidulans-based model. In the case of the iron assimilation proteins, the number of orthologs showed a positive correlation with H2O2-induced stress tolerance while the number of MpkC orthologs correlated positively with Congo Red induced cell wall stress, sorbitol induced osmotic stress and H2O2 induced oxidative stress tolerances. For most stress proteins, changes in the number of orthologs did not correlate well with any stress tolerance attributes. As a consequence, stress tolerance patterns of the studied Aspergilli did not correlate with either the sets of stress response proteins in general or with the phylogeny of the species studied. These observations suggest that stress protein duplication and deletion events significantly contributed to the evolution of stress tolerance attributes of Aspergilli. In contrast, there are other processes, which may counterbalance the effects of stress gene duplications or deletions including (i) alterations in the structures of stress proteins leading to changes in their biological activities, (ii) varying biosynthesis of stress proteins, (iii) rewiring stress response regulatory networks or even (iv) acquiring new stress response genes by horizontal gene transfer. All these multilevel changes are indispensable for the successful adaptation of filamentous fungi to altering environmental conditions, especially when these organisms are entering new ecological niches.

Single-cell RNA sequencing technologies and bioinformatics pipelines

Rapid progress in the development of next-generation sequencing (NGS) technologies in recent years has provided many valuable insights into complex biological systems, ranging from cancer genomics to diverse microbial communities. NGS-based technologies for genomics, transcriptomics, and epigenomics are now increasingly focused on the characterization of individual cells. These single-cell analyses will allow researchers to uncover new and potentially unexpected biological discoveries relative to traditional profiling methods that assess bulk populations. Single-cell RNA sequencing (scRNA-seq), for example, can reveal complex and rare cell populations, uncover regulatory relationships between genes, and track the trajectories of distinct cell lineages in development. In this review, we will focus on technical challenges in single-cell isolation and library preparation and on computational analysis pipelines available for analyzing scRNA-seq data. Further technical improvements at the level of molecular and cell biology and in available bioinformatics tools will greatly facilitate both the basic science and medical applications of these sequencing technologies.

Showing which genes are expressed, or switched on, in individual cells may help to reveal the first signs of disease. Each cell in an organism contains the same genetic information, but cell type and behavior depend on which genes are expressed. Previously, researchers could only sequence cells in batches, averaging the results, but technological improvements now allow sequencing of the genes expressed in an individual cell, known as single-cell RNA sequencing (scRNA-seq). Ji Hyun Lee (Kyung Hee University, Seoul) and Duhee Bang and Byungjin Hwang (Yonsei University, Seoul) have reviewed the available scRNA-seq technologies and the strategies available to analyze the large quantities of data produced. They conclude that scRNA-seq will impact both basic and medical science, from illuminating drug resistance in cancer to revealing the complex pathways of cell differentiation during development.

Functions of the COPII gene paralogs SEC23A and SEC23B are interchangeable in vivo

Approximately one-third of the mammalian proteome is transported from the endoplasmic reticulum-to-Golgi via COPII-coated vesicles. SEC23, a core component of coat protein-complex II (COPII), is encoded by two paralogous genes in vertebrates (Sec23a and Sec23b). In humans, SEC23B deficiency results in congenital dyserythropoietic anemia type-II (CDAII), while SEC23A deficiency results in a skeletal phenotype (with normal red blood cells). These distinct clinical disorders, together with previous biochemical studies, suggest unique functions for SEC23A and SEC23B. Here we show indistinguishable intracellular protein interactomes for human SEC23A and SEC23B, complementation of yeast Sec23 by both human and murine SEC23A/B, and rescue of the lethality of sec23b deficiency in zebrafish by a sec23a-expressing transgene. We next demonstrate that a Sec23a coding sequence inserted into the murine Sec23b locus completely rescues the lethal SEC23B-deficient pancreatic phenotype. We show that SEC23B is the predominantly expressed paralog in human bone marrow, but not in the mouse, with the reciprocal pattern observed in the pancreas. Taken together, these data demonstrate an equivalent function for SEC23A/B, with evolutionary shifts in the transcription program likely accounting for the distinct phenotypes of SEC23A/B deficiency within and across species, a paradigm potentially applicable to other sets of paralogous genes. These findings also suggest that enhanced erythroid expression of the normal SEC23A gene could offer an effective therapeutic approach for CDAII patients.

Systematic Study of Nucleosome-Displacing Factors in Budding Yeast

Nucleosomes present a barrier for the binding of most transcription factors (TFs). However, special TFs known as nucleosome-displacing factors (NDFs) can access embedded sites and cause the depletion of the local nucleosomes as well as repositioning of the neighboring nucleosomes. Here, we developed a novel high-throughput method in yeast to identify NDFs among 104 TFs and systematically characterized the impact of orientation, affinity, location, and copy number of their binding motifs on the nucleosome occupancy. Using this assay, we identified 29 NDF motifs and divided the nuclear TFs into three groups with strong, weak, and no nucleosome-displacing activities. Further studies revealed that tight DNA binding is the key property that underlies NDF activity, and the NDFs may partially rely on the DNA replication to compete with nucleosome. Overall, our study presents a framework to functionally characterize NDFs and elucidate the mechanism of nucleosome invasion.

Dynamic Changes in Yeast Phosphatase Families Allow for Specialization in Phosphate and Thiamine Starvation

Convergent evolution is often due to selective pressures generating a similar phenotype. We observe relatively recent duplications in a spectrum of Saccharomycetaceae yeast species resulting in multiple phosphatases that are regulated by different nutrient conditions – thiamine and phosphate starvation. This specialization is both transcriptional and at the level of phosphatase substrate specificity. In Candida glabrata, loss of the ancestral phosphatase family was compensated by the co-option of a different histidine phosphatase family with three paralogs. Using RNA-seq and functional assays, we identify one of these paralogs, CgPMU3, as a thiamine phosphatase. We further determine that the 81% identical paralog CgPMU2 does not encode thiamine phosphatase activity; however, both are capable of cleaving the phosphatase substrate, 1-napthyl-phosphate. We functionally demonstrate that members of this family evolved novel enzymatic functions for phosphate and thiamine starvation, and are regulated transcriptionally by either nutrient condition, and observe similar trends in other yeast species. This independent, parallel evolution involving two different families of histidine phosphatases suggests that there were likely similar selective pressures on multiple yeast species to recycle thiamine and phosphate. In this work, we focused on duplication and specialization, but there is also repeated loss of phosphatases, indicating that the expansion and contraction of the phosphatase family is dynamic in many Ascomycetes. The dynamic evolution of the phosphatase gene families is perhaps just one example of how gene duplication, co-option, and transcriptional and functional specialization together allow species to adapt to their environment with existing genetic resources.

The mycoparasitic fungus Clonostachys rosea responds with both common and specific gene expression during interspecific interactions with fungal prey

Clonostachys rosea is a necrotrophic mycoparasitic fungus, used for biological control of plant pathogenic fungi. A better understanding of the underlying mechanisms resulting in successful biocontrol is important for knowledge-based improvements of the application and use of biocontrol in agricultural production systems. Transcriptomic analyses revealed that C. rosea responded with both common and specific gene expression during interactions with the fungal prey species Botrytis cinerea and Fusarium graminearum. Genes predicted to encode proteins involved in membrane transport, biosynthesis of secondary metabolites and carbohydrate-active enzymes were induced during the mycoparasitic attack. Predicted major facilitator superfamily (MFS) transporters constituted 54% of the induced genes, and detailed phylogenetic and evolutionary analyses showed that a majority of these genes belonged to MFS gene families evolving under selection for increased paralog numbers, with predicted functions in drug resistance and transport of carbohydrates and small organic compounds. Sequence analysis of MFS transporters from family 2.A.1.3.65 identified rapidly evolving loop regions forming the entry to the transport tunnel, indicating changes in substrate specificity as a target for selection. Deletion of the MFS transporter gene mfs464 resulted in mutants with increased growth inhibitory activity against F. graminearum, providing evidence for a function in interspecific fungal interactions. In summary, we show that the mycoparasite C. rosea can distinguish between fungal prey species and modulate its transcriptomic responses accordingly. Gene expression data emphasize the importance of secondary metabolites in mycoparasitic interactions.

Genomic Structure and Tissue Expression of the NK-Lysin Gene Family in Bison

Antimicrobial peptides (AMPs) are a class of natural peptides with varying numbers of amino acids. They are principal components of innate immunity in vertebrates, encoding natural antibiotics and providing a protective response against a broad range of microbes including those responsible for tuberculosis, an important disease in bison. NK-lysins are AMPs that have been described in various organisms and are coded by a single gene in several mammalian species, including human. Recently, we described a family of 4 NK-lysin genes in cattle. Here, we examined NK-lysin genes in bison and identified 4 bison paralogs (NK1, NK2A, NK2B, and NK2C), although the current bison genome assembly annotates only 2 (NK1 and NK2). Sequence and phylogenetic analysis support the triplication of NK2 prior to the most recent common ancestor of bison and cattle. Comparative mapping of bison and cattle paralogs indicates that the NK-lysin family is located on bison chromosome 11 with well-conserved synteny of flanking genes relative to cattle. The 3 bison NK-lysin2 genes share high sequence similarity with each other. RNA-seq analysis demonstrates that NK2A, NK2B, and NK2C are expressed primarily in the lung, whereas NK1 is expressed at low levels in all tissues studied. This tissue expression pattern differs from that previously reported for cattle, suggesting some divergence in function since the evolutionary separation of the 2 species.

The evolution of a series of behavioral traits is associated with autism-risk genes in cavefish

Background

An essential question in evolutionary biology is whether shifts in a set of polygenic behaviors share a genetic basis across species. Such a behavioral shift is seen in the cave-dwelling Mexican tetra, Astyanax mexicanus. Relative to surface-dwelling conspecifics, cavefish do not school (asocial), are hyperactive and sleepless, adhere to a particular vibration stimulus (imbalanced attention), behave repetitively, and show elevated stress hormone levels. Interestingly, these traits largely overlap with the core symptoms of human autism spectrum disorder (ASD), raising the possibility that these behavioral traits are underpinned by a similar set of genes (i.e. a repeatedly used suite of genes).

Result

Here, we explored whether modification of ASD-risk genes underlies cavefish evolution. Transcriptomic analyses revealed that > 58.5% of 3152 cavefish orthologs to ASD-risk genes are significantly up- or down-regulated in the same direction as genes in postmortem brains from ASD patients. Enrichment tests suggest that ASD-risk gene orthologs in A. mexicanus have experienced more positive selection than other genes across the genome. Notably, these positively selected cavefish ASD-risk genes are enriched for pathways involved in gut function, inflammatory diseases, and lipid/energy metabolism, similar to symptoms that frequently coexist in ASD patients. Lastly, ASD drugs mitigated cavefish’s ASD-like behaviors, implying shared aspects of neural processing.

Conclusion

Overall, our study indicates that ASD-risk genes and associated pathways (especially digestive, immune and metabolic pathways) may be repeatedly used for shifts in polygenic behaviors across evolutionary time.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1199-9) contains supplementary material, which is available to authorized users.

Evaluation of functionality for serine and threonine phosphorylation with different evolutionary ages in human and mouse

BACKGROUND

Rapid evolution of phosphorylation sites could provide raw materials of natural selection to fit the environment by rewiring the regulation of signal pathways. However, a large part of phosphorylation sites was suggested to be non-functional. Although the new-arising phosphorylation sites with little functional implications prevailed in fungi, the evolutionary performance of vertebrate phosphorylation sites remained elusive.

RESULTS

In this study, we evaluated the functionality of human and mouse phosphorylation sites by dividing them into old, median and young age groups based on the phylogeny of vertebrates. We found the sites in the old group were more likely to be functional and involved in signaling pathways than those in the young group. A smaller proportion of sites in the young group originated from aspartate/glutamate, which could restore the ancestral functions. In addition, both the phosphorylation level and breadth was increased with the evolutionary age. Similar to cases in fungi, these results implied that the newly emerged phosphorylation sites in vertebrates were also more likely to be non-functional, especially for serine and threonine phosphorylation in disordered regions.

CONCLUSIONS

This study provided not only insights into the dynamics of phosphorylation evolution in vertebrates, but also new clues to identify the functional phosphorylation sites from massive noisy data.

Formal description of sequence-based voucherless Fungi: promises and pitfalls, and how to resolve them

There is urgent need for a formal nomenclature of sequence-based, voucherless Fungi, given that environmental sequencing has accumulated more than one billion fungal ITS reads in the Sequence Read Archive, about 1,000 times as many as fungal ITS sequences in GenBank. These unnamed Fungi could help to bridge the gap between 115,000 to 140,000 currently accepted and 2.2 to 3.8 million predicted species, a gap that cannot realistically be filled using specimen or culture-based inventories. The Code never aimed at placing restrictions on the nature of characters chosen for taxonomy, and the requirement for physical types is now becoming a constraint on the advancement of science. We elaborate on the promises and pitfalls of sequence-based nomenclature and provide potential solutions to major concerns of the mycological community. Types of sequence-based taxa, which by default lack a physical specimen or culture, could be designated in four alternative ways: (1) the underlying sample ('bag' type), (2) the DNA extract, (3) fluorescent in situ hybridization (FISH), or (4) the type sequence itself. Only (4) would require changes to the Code and the latter would be the most straightforward approach, complying with three of the five principal functions of types better than physical specimens. A fifth way, representation of the sequence in an illustration, has been ruled as unacceptable in the Code. Potential flaws in sequence data are analogous to flaws in physical types, and artifacts are manageable if a stringent analytical approach is applied. Conceptual errors such as homoplasy, intragenomic variation, gene duplication, hybridization, and horizontal gene transfer, apply to all molecular approaches and cannot be used as a specific argument against sequence-based nomenclature. The potential impact of these phenomena is manageable, as phylogenetic species delimitation has worked satisfactorily in Fungi. The most serious shortcoming of sequence-based nomenclature is the likelihood of parallel classifications, either by describing taxa that already have names based on physical types, or by using different markers to delimit species within the same lineage. The probability of inadvertently establishing sequence-based species that have names available is between 20.4 % and 1.5 % depending on the number of globally predicted fungal species. This compares favourably to a historical error rate of about 30 % based on physical types, and this rate could be reduced to practically zero by adding specific provisions to this approach in the Code. To avoid parallel classifications based on different markers, sequence-based nomenclature should be limited to a single marker, preferably the fungal ITS barcoding marker; this is possible since sequence-based nomenclature does not aim at accurate species delimitation but at naming lineages to generate a reference database, independent of whether these lineages represent species, closely related species complexes, or infraspecies. We argue that clustering methods are inappropriate for sequence-based nomenclature; this approach must instead use phylogenetic methods based on multiple alignments, combined with quantitative species recognition methods. We outline strategies to obtain higher-level phylogenies for ITS-based, voucherless species, including phylogenetic binning, 'hijacking' species delimitation methods, and temporal banding. We conclude that voucherless, sequence-based nomenclature is not a threat to specimen and culture-based fungal taxonomy, but a complementary approach capable of substantially closing the gap between known and predicted fungal diversity, an approach that requires careful work and high skill levels.

Sources of Fungal Genetic Variation and Associating It with Phenotypic Diversity

The first eukaryotic genome to be sequenced was fungal, and there continue to be more sequenced genomes in the kingdom Fungi than in any other eukaryotic kingdom. Comparison of these genomes reveals many sources of genetic variation, from single nucleotide polymorphisms to horizontal gene transfer and on to changes in the arrangement and number of chromosomes, not to mention endofungal bacteria and viruses. Population genomics shows that all sources generate variation all the time and implicate natural selection as the force maintaining genome stability. Variation in wild populations is a rich resource for associating genetic variation with phenotypic variation, whether through quantitative trait locus mapping, genome-wide association studies, or reverse ecology. Subjects of studies associating genetic and phenotypic variation include model fungi, e.g., Saccharomyces and Neurospora, but pioneering studies have also been made with fungi pathogenic to plants, e.g., Pyricularia (= Magnaporthe), Zymoseptoria, and Fusarium, and to humans, e.g., Coccidioides, Cryptococcus, and Candida.