Structure of Kluyveromyces lactis subtelomeres: duplications and gene content

Mechanisms of adaptation and evolution in Toxoplasma gondii

Toxoplasma has high host flexibility, infecting all nucleated cells of mammals and birds. This implies that during its infective process the parasite must constantly adapt to different environmental situations, which in turn leads to modifications in its metabolism, regulation of gene transcription, translation of mRNAs and stage specific factors. There are conserved pathways that support these adaptations, which we aim to elucidate in this review. We begin by exploring the widespread epigenetic mechanisms and transcription regulators, continue with the supportive role of Heat Shock Proteins (Hsp), the translation regulation, stress granules, and finish with the emergence of contingency genes in highly variable genomic domains, such as subtelomeres. Within epigenetics, the discovery of a new histone variant of the H2B family (H2B.Z), contributing to T. gondii virulence and differentiation, but also gene expression regulation and its association with the metabolic state of the parasite, is highlighted. Associated with the regulation of gene expression are transcription factors (TFs). An overview of the main findings on TF and development is presented. We also emphasize the role of Hsp90 and Tgj1 in T. gondii metabolic fitness and the regulation of protein translation. Translation regulation is also highlighted as a mechanism for adaptation to conditions encountered by the parasite as well as stress granules containing mRNA and proteins generated in the extracellular tachyzoite. Another important aspect in evolution and adaptability are the subtelomeres because of their high variability and gene duplication rate. Toxoplasma possess multigene families of membrane proteins and contingency genes that are associated with different metabolic stresses. Among them parasite differentiation and environmental stresses stand out, including those that lead tachyzoite to bradyzoite conversion. Finally, we are interested in positioning protozoa as valuable evolution models, focusing on research related to the Extended Evolutionary Synthesis, based on models recently generated, such as extracellular adaptation and ex vivo cyst recrudescence.

Dynamics of the Streptomyces chromosome: chance and necessity

Streptomyces are prolific producers of specialized metabolites with applications in medicine and agriculture. Remarkably, these bacteria possess a large linear chromosome that is genetically compartmentalized: core genes are grouped in the central part, while the ends are populated by poorly conserved genes including antibiotic biosynthetic gene clusters. The genome is highly unstable and exhibits distinct evolutionary rates along the chromosome. Recent chromosome conformation capture (3C) and comparative genomics studies have shed new light on the interplay between genome dynamics in space and time. Here, we review insights that illustrate how the balance between chance (random genome variations) and necessity (structural and functional constraints) may have led to the emergence of spatial structuring of the Streptomyces chromosome.

Evolution and functional diversification of yeast sugar transporters

While simple sugars such as monosaccharides and disaccharide are the typical carbon source for most yeasts, whether a species can grow on a particular sugar is generally a consequence of presence or absence of a suitable transporter to enable its uptake. The most common transporters that mediate sugar import in yeasts belong to the major facilitator superfamily (MFS). Some of these, for example the Saccharomyces cerevisiae Hxt proteins have been extensively studied, but detailed information on many others is sparce. In part, this is because there are many lineages of MFS transporters that are either absent from, or poorly represented in, the model S. cerevisiae, which actually has quite a restricted substrate range. It is important to address this knowledge gap to gain better understanding of the evolution of yeasts and to take advantage of sugar transporters to exploit or engineer yeasts for biotechnological applications. This article examines the full repertoire of MFS proteins in representative budding yeasts (Saccharomycotina). A comprehensive analysis of 139 putative sugar transporters retrieved from 10 complete genomes sheds new light on the diversity and evolution of this family. Using the phylogenetic lens, it is apparent that proteins have often been misassigned putative functions and this can now be corrected. It is also often seen that patterns of expansion of particular genes reflects the differential importance of transport of specific sugars (and related molecules) in different yeasts, and this knowledge also provides an improved resource for the selection or design of tailored transporters.

Subtelomeres as Specialized Chromatin Domains

Specificities associated with chromosomal linearity are not restricted to telomeres. Here, recent results obtained on fission and budding yeast are summarized and an attempt is made to define subtelomeres using chromatin features extending beyond the heterochromatin emanating from telomeres. Subtelomeres, the chromosome domains adjacent to telomeres, differ from the rest of the genome by their gene content, rapid evolution, and chromatin features that together contribute to organism adaptation. However, current definitions of subtelomeres are generally based on synteny and are largely gene-centered. Taking into consideration both the peculiar gene content and dynamics as well as the chromatin properties of those domains, it is discussed how chromatin features can contribute to subtelomeric properties and functions, and play a pivotal role in the emergence of subtelomeres.

Origin of Lactose Fermentation in Kluyveromyces lactis by Interspecies Transfer of a Neo-functionalized Gene Cluster during Domestication

Humans have used yeasts to make cheese and kefir for millennia, but the ability to ferment the milk sugar lactose is found in only a few yeast species, of which the foremost is Kluyveromyces lactis [1]. Two genes, LAC12 (lactose permease) and LAC4 (lactase), are sufficient for lactose uptake and hydrolysis to glucose and galactose [2]. Here, we show that these genes have a complex evolutionary history in the genus Kluyveromyces that is likely the result of human activity during domestication. We show that the ancestral Lac12 was bifunctional, able to import both lactose and cellobiose into the cell. These disaccharides were then hydrolyzed by Lac4 in the case of lactose or Cel2 in the case of cellobiose. A second cellobiose transporter, Cel1, was also present ancestrally. In the K. lactis lineage, the ancestral LAC12 and LAC4 were lost and a separate upheaval in the sister species K. marxianus resulted in loss of CEL1 and quadruplication of LAC12. One of these LAC12 genes became neofunctionalized to encode an efficient lactose transporter capable of supporting fermentation, specifically in dairy strains of K. marxianus, where it formed a LAC4-LAC12-CEL2 gene cluster, although another remained a cellobiose transporter. Then, the ability to ferment lactose was acquired very recently by K. lactis var. lactis by introgression of LAC12 and LAC4 on a 15-kb subtelomeric region from a dairy strain of K. marxianus. The genomic history of the LAC genes shows that strong selective pressures were imposed on yeasts by early dairy farmers.

Genome Mining of Non-Conventional Yeasts: Search and Analysis of MAL Clusters and Proteins

Genomic clustering of functionally related genes is rare in yeasts and other eukaryotes with only few examples available. Here, we summarize our data on a nontelomeric MAL cluster of a non-conventional methylotrophic yeast Ogataea (Hansenula) polymorpha containing genes for α-glucosidase MAL1, α-glucoside permease MAL2 and two hypothetical transcriptional activators. Using genome mining, we detected MAL clusters of varied number, position and composition in many other maltose-assimilating non-conventional yeasts from different phylogenetic groups. The highest number of MAL clusters was detected in Lipomyces starkeyi while no MAL clusters were found in Schizosaccharomyces pombe and Blastobotrys adeninivorans. Phylograms of α-glucosidases and α-glucoside transporters of yeasts agreed with phylogenesis of the respective yeast species. Substrate specificity of unstudied α-glucosidases was predicted from protein sequence analysis. Specific activities of Scheffersomycesstipitis α-glucosidases MAL7, MAL8, and MAL9 heterologously expressed in Escherichia coli confirmed the correctness of the prediction-these proteins were verified promiscuous maltase-isomaltases. α-Glucosidases of earlier diverged yeasts L. starkeyi, B. adeninivorans and S. pombe showed sequence relatedness with α-glucosidases of filamentous fungi and bacilli.

Polymorphisms in the LAC12 gene explain lactose utilisation variability in Kluyveromyces marxianus strains

Kluyveromyces marxianus is a safe yeast used in the food and biotechnology sectors. One of the important traits that sets it apart from the familiar yeasts, Saccharomyces cerevisiae, is its capacity to grow using lactose as a carbon source. Like in its close relative, Kluyveromyces lactis, this requires lactose transport via a permease and intracellular hydrolysis of the disaccharide. Given the importance of the trait, it was intriguing that most, but not all, strains of K. marxianus are reported to consume lactose efficiently. In this study, primarily through heterologous expression in S. cerevisiae and K. marxianus, it was established that a single gene, LAC12, is responsible for lactose uptake in K. marxianus. Strains that failed to transport lactose showed variation in 13 amino acids in the Lac12p protein, rendering the protein non-functional for lactose transport. Genome analysis showed that the LAC12 gene is present in four copies in the subtelomeric regions of three different chromosomes but only the ancestral LAC12 gene encodes a functional lactose transporter. Other copies of LAC12 may be non-functional or have alternative substrates. The analysis raises some interesting questions regarding the evolution of sugar transporters in K. marxianus.

Mild Telomere Dysfunction as a Force for Altering the Adaptive Potential of Subtelomeric Genes

Subtelomeric regions have several unusual characteristics, including complex repetitive structures, increased rates of evolution, and enrichment for genes involved in niche adaptation. The adaptive telomere failure hypothesis suggests that certain environmental stresses can induce a low level of telomere failure, potentially leading to elevated subtelomeric recombination that could result in adaptive mutational changes within subtelomeric genes. Here, we tested a key prediction of the adaptive telomere failure hypothesis-that telomere dysfunction mild enough to have little or no overall effect on cell fitness could still lead to substantial increases in the mutation rates of subtelomeric genes. Our results show that a mutant of Kluyveromyces lactis with stably short telomeres produced a large increase in the frequency of mutations affecting the native subtelomeric β-galactosidase (LAC4) gene. All lac4 mutants examined from strains with severe telomere dysfunction underwent terminal deletion/duplication events consistent with being due to break-induced replication. In contrast, although cells with mild telomere dysfunction also exhibited similar terminal deletion and duplication events, up to 50% of lac4 mutants from this background unexpectedly contained base changes within the LAC4 coding region. This mutational bias for producing base changes demonstrates that mild telomere dysfunction can be well suited as a force for altering the adaptive potential of subtelomeric genes.

Complete Genome Sequence of Kluyveromyces lactis Strain GG799, a Common Yeast Host for Heterologous Protein Expression

We report the genome sequence of the dairy yeast Kluyveromyces lactis strain GG799 obtained using the Pacific Biosciences RS II platform. K. lactis strain GG799 is a common host for the expression of proteins at both laboratory and industrial scales.

Genome Diversity and Evolution in the Budding Yeasts (Saccharomycotina)

Considerable progress in our understanding of yeast genomes and their evolution has been made over the last decade with the sequencing, analysis, and comparisons of numerous species, strains, or isolates of diverse origins. The role played by yeasts in natural environments as well as in artificial manufactures, combined with the importance of some species as model experimental systems sustained this effort. At the same time, their enormous evolutionary diversity (there are yeast species in every subphylum of Dikarya) sparked curiosity but necessitated further efforts to obtain appropriate reference genomes. Today, yeast genomes have been very informative about basic mechanisms of evolution, speciation, hybridization, domestication, as well as about the molecular machineries underlying them. They are also irreplaceable to investigate in detail the complex relationship between genotypes and phenotypes with both theoretical and practical implications. This review examines these questions at two distinct levels offered by the broad evolutionary range of yeasts: inside the best-studied Saccharomyces species complex, and across the entire and diversified subphylum of Saccharomycotina. While obviously revealing evolutionary histories at different scales, data converge to a remarkably coherent picture in which one can estimate the relative importance of intrinsic genome dynamics, including gene birth and loss, vs. horizontal genetic accidents in the making of populations. The facility with which novel yeast genomes can now be studied, combined with the already numerous available reference genomes, offer privileged perspectives to further examine these fundamental biological questions using yeasts both as eukaryotic models and as fungi of practical importance.

Contrasting evolutionary genome dynamics between domesticated and wild yeasts

Structural rearrangements have long been recognized as an important source of genetic variation, with implications in phenotypic diversity and disease, yet their detailed evolutionary dynamics remain elusive. Here we use long-read sequencing to generate end-to-end genome assemblies for 12 strains representing major subpopulations of the partially domesticated yeast Saccharomyces cerevisiae and its wild relative Saccharomyces paradoxus. These population-level high-quality genomes with comprehensive annotation enable precise definition of chromosomal boundaries between cores and subtelomeres and a high-resolution view of evolutionary genome dynamics. In chromosomal cores, S. paradoxus shows faster accumulation of balanced rearrangements (inversions, reciprocal translocations and transpositions), whereas S. cerevisiae accumulates unbalanced rearrangements (novel insertions, deletions and duplications) more rapidly. In subtelomeres, both species show extensive interchromosomal reshuffling, with a higher tempo in S. cerevisiae. Such striking contrasts between wild and domesticated yeasts are likely to reflect the influence of human activities on structural genome evolution.

Production of the antibiotic secondary metabolite solanapyrone A by the fungal plant pathogen Ascochyta rabiei during fruiting body formation in saprobic growth

Fungi are noted producers of a diverse array of secondary metabolites, many of which are of pharmacological importance. However, the biological roles of the vast majority of these molecules during the fungal life cycle in nature remain elusive. Solanapyrones are polyketide-derived secondary metabolites produced by diverse fungal species including the plant pathogen Ascochyta rabiei. This molecule was originally thought to function as a phytotoxin facilitating pathogenesis of A. rabiei. Chemical profiling and gene expression studies showed that solanapyrone A was specifically produced during saprobic, but not parasitic growth of A. rabiei. Expression of the gene encoding the final enzymatic step in solanapyrone biosynthesis was specifically associated with development of the asexual fruiting bodies of the fungus on certain substrates. In confrontation assays with saprobic fungi that were commonly found in chickpea debris in fields, A. rabiei effectively suppressed the growth of all competing fungi, such as Alternaria, Epicoccum and Ulocladium species. Solanapyrone A was directly detected in the inhibitory zone using a MALDI-imaging mass spectrometry, and the purified compound showed significant antifungal activities against the potential saprobic competitors. These results suggest that solanapyrone A plays an important role for competition and presumably the survival of the fungus.

A Novel Type Pathway-Specific Regulator and Dynamic Genome Environments of a Solanapyrone Biosynthesis Gene Cluster in the Fungus Ascochyta rabiei

Secondary metabolite genes are often clustered together and situated in particular genomic regions, like the subtelomere, that can facilitate niche adaptation in fungi. Solanapyrones are toxic secondary metabolites produced by fungi occupying different ecological niches. Full-genome sequencing of the ascomycete Ascochyta rabiei revealed a solanapyrone biosynthesis gene cluster embedded in an AT-rich region proximal to a telomere end and surrounded by Tc1/Mariner-type transposable elements. The highly AT-rich environment of the solanapyrone cluster is likely the product of repeat-induced point mutations. Several secondary metabolism-related genes were found in the flanking regions of the solanapyrone cluster. Although the solanapyrone cluster appears to be resistant to repeat-induced point mutations, a P450 monooxygenase gene adjacent to the cluster has been degraded by such mutations. Among the six solanapyrone cluster genes (sol1 to sol6), sol4 encodes a novel type of Zn(II)2Cys6 zinc cluster transcription factor. Deletion of sol4 resulted in the complete loss of solanapyrone production but did not compromise growth, sporulation, or virulence. Gene expression studies with the sol4 deletion and sol4-overexpressing mutants delimited the boundaries of the solanapyrone gene cluster and revealed that sol4 is likely a specific regulator of solanapyrone biosynthesis and appears to be necessary and sufficient for induction of the solanapyrone cluster genes. Despite the dynamic surrounding genomic regions, the solanapyrone gene cluster has maintained its integrity, suggesting important roles of solanapyrones in fungal biology.

Complete DNA sequence of Kuraishia capsulata illustrates novel genomic features among budding yeasts (Saccharomycotina)

The numerous yeast genome sequences presently available provide a rich source of information for functional as well as evolutionary genomics but unequally cover the large phylogenetic diversity of extant yeasts. We present here the complete sequence of the nuclear genome of the haploid-type strain of Kuraishia capsulata (CBS1993^T), a nitrate-assimilating Saccharomycetales of uncertain taxonomy, isolated from tunnels of insect larvae underneath coniferous barks and characterized by its copious production of extracellular polysaccharides. The sequence is composed of seven scaffolds, one per chromosome, totaling 11.4 Mb and containing 6,029 protein-coding genes, ∼13.5% of which being interrupted by introns. This GC-rich yeast genome (45.7%) appears phylogenetically related with the few other nitrate-assimilating yeasts sequenced so far, Ogataea polymorpha, O. parapolymorpha, and Dekkera bruxellensis, with which it shares a very reduced number of tRNA genes, a novel tRNA sparing strategy, and a common nitrate assimilation cluster, three specific features to this group of yeasts. Centromeres were recognized in GC-poor troughs of each scaffold. The strain bears MAT alpha genes at a single MAT locus and presents a significant degree of conservation with Saccharomyces cerevisiae genes, suggesting that it can perform sexual cycles in nature, although genes involved in meiosis were not all recognized. The complete absence of conservation of synteny between K. capsulata and any other yeast genome described so far, including the three other nitrate-assimilating species, validates the interest of this species for long-range evolutionary genomic studies among Saccharomycotina yeasts.

MFS transporters required for multidrug/multixenobiotic (MD/MX) resistance in the model yeast: understanding their physiological function through post-genomic approaches

Multidrug/Multixenobiotic resistance (MDR/MXR) is a widespread phenomenon with clinical, agricultural and biotechnological implications, where MDR/MXR transporters that are presumably able to catalyze the efflux of multiple cytotoxic compounds play a key role in the acquisition of resistance. However, although these proteins have been traditionally considered drug exporters, the physiological function of MDR/MXR transporters and the exact mechanism of their involvement in resistance to cytotoxic compounds are still open to debate. In fact, the wide range of structurally and functionally unrelated substrates that these transporters are presumably able to export has puzzled researchers for years. The discussion has now shifted toward the possibility of at least some MDR/MXR transporters exerting their effect as the result of a natural physiological role in the cell, rather than through the direct export of cytotoxic compounds, while the hypothesis that MDR/MXR transporters may have evolved in nature for other purposes than conferring chemoprotection has been gaining momentum in recent years. This review focuses on the drug transporters of the Major Facilitator Superfamily (MFS; drug:H⁺ antiporters) in the model yeast Saccharomyces cerevisiae. New insights into the natural roles of these transporters are described and discussed, focusing on the knowledge obtained or suggested by post-genomic research. The new information reviewed here provides clues into the unexpectedly complex roles of these transporters, including a proposed indirect regulation of the stress response machinery and control of membrane potential and/or internal pH, with a special emphasis on a genome-wide view of the regulation and evolution of MDR/MXR-MFS transporters.

Discovering the secondary metabolite potential encoded within entomopathogenic fungi

This highlight discusses the secondary metabolite potential of the insect pathogensMetarhiziumandBeauveria, including a bioinformatics analysis of secondary metabolite genes for which no products are yet identified. (Top picture is a mole cricket infected withBeauveria bassianaand the bottom picture is a wasp infected withBeauveria bassiana.)

The drug:H⁺ antiporters of family 2 (DHA2), siderophore transporters (ARN) and glutathione:H⁺ antiporters (GEX) have a common evolutionary origin in hemiascomycete yeasts

BACKGROUND

The Saccharomyces cerevisiae 14-spanner Drug:H+ Antiporter family 2 (DHA2) are transporters of the Major Facilitator Superfamily (MFS) involved in multidrug resistance (MDR). Although poorly characterized, DHA2 family members were found to participate in the export of structurally and functionally unrelated compounds or in the uptake of amino acids into the vacuole or the cell. In S. cerevisiae, the four ARN/SIT family members encode siderophore transporters and the two GEX family members encode glutathione extrusion pumps. The evolutionary history of DHA2, ARN and GEX genes, encoding 14-spanner MFS transporters, is reconstructed in this study.

RESULTS

The translated ORFs of 31 strains from 25 hemiascomycetous species, including 10 pathogenic Candida species, were compared using a local sequence similarity algorithm. The constraining and traversing of a network representing the pairwise similarity data gathered 355 full size proteins and retrieved ARN and GEX family members together with DHA2 transporters, suggesting the existence of a close phylogenetic relationship among these 14-spanner major facilitators. Gene neighbourhood analysis was combined with tree construction methodologies to reconstruct their evolutionary history and 7 DHA2 gene lineages, 5 ARN gene lineages, and 1 GEX gene lineage, were identified. The S. cerevisiae DHA2 proteins Sge1, Azr1, Vba3 and Vba5 co-clustered in a large phylogenetic branch, the ATR1 and YMR279C genes were proposed to be paralogs formed during the Whole Genome Duplication (WGD) whereas the closely related ORF YOR378W resides in its own lineage. Homologs of S. cerevisiae DHA2 vacuolar proteins Vba1, Vba2 and Vba4 occur widespread in the Hemiascomycetes. Arn1/Arn2 homologs were only found in species belonging to the Saccharomyces complex and are more abundant in the pre-WGD species. Arn4 homologs were only found in sub-telomeric regions of species belonging to the Sacharomyces sensu strictu group (SSSG). Arn3 type siderophore transporters are abundant in the Hemiascomycetes and form an ancient gene lineage extending to the filamentous fungi.

CONCLUSIONS

The evolutionary history of DHA2, ARN and GEX genes was reconstructed and a common evolutionary root shared by the encoded proteins is hypothesized. A new protein family, denominated DAG, is proposed to span these three phylogenetic subfamilies of 14-spanner MFS transporters.

Extensive intra-kingdom horizontal gene transfer converging on a fungal fructose transporter gene

Comparative genomics revealed in the last decade a scenario of rampant horizontal gene transfer (HGT) among prokaryotes, but for fungi a clearly dominant pattern of vertical inheritance still stands, punctuated however by an increasing number of exceptions. In the present work, we studied the phylogenetic distribution and pattern of inheritance of a fungal gene encoding a fructose transporter (FSY1) with unique substrate selectivity. 109 FSY1 homologues were identified in two sub-phyla of the Ascomycota, in a survey that included 241 available fungal genomes. At least 10 independent inter-species instances of horizontal gene transfer (HGT) involving FSY1 were identified, supported by strong phylogenetic evidence and synteny analyses. The acquisition of FSY1 through HGT was sometimes suggestive of xenolog gene displacement, but several cases of pseudoparalogy were also uncovered. Moreover, evidence was found for successive HGT events, possibly including those responsible for transmission of the gene among yeast lineages. These occurrences do not seem to be driven by functional diversification of the Fsy1 proteins because Fsy1 homologues from widely distant lineages, including at least one acquired by HGT, appear to have similar biochemical properties. In summary, retracing the evolutionary path of the FSY1 gene brought to light an unparalleled number of independent HGT events involving a single fungal gene. We propose that the turbulent evolutionary history of the gene may be linked to the unique biochemical properties of the encoded transporter, whose predictable effect on fitness may be highly variable. In general, our results support the most recent views suggesting that inter-species HGT may have contributed much more substantially to shape fungal genomes than heretofore assumed.

Detection and characterization of megasatellites in orthologous and nonorthologous genes of 21 fungal genomes

Megasatellites are large DNA tandem repeats, originally described in Candida glabrata, in protein-coding genes. Most of the genes in which megasatellites are found are of unknown function. In this work, we extended the search for megasatellites to 20 additional completely sequenced fungal genomes and extracted 216 megasatellites in 203 out of 142,121 genes, corresponding to the most exhaustive description of such genetic elements available today. We show that half of the megasatellites detected encode threonine-rich peptides predicted to be intrinsically disordered, suggesting that they may interact with several partners or serve as flexible linkers. Megasatellite motifs were clustered into several families. Their distribution in fungal genes shows that different motifs are found in orthologous genes and similar motifs are found in unrelated genes, suggesting that megasatellite formation or spreading does not necessarily track the evolution of their host genes. Altogether, these results suggest that megasatellites are created and lost during evolution of fungal genomes, probably sharing similar functions, although their primary sequences are not necessarily conserved.

Functional divergence of gene duplicates through ectopic recombination

This report reveals that duplicated genes undergo ectopic recombination, which leads to new chimaeric alleles. Mimicking these intergenic recombination events creates chimaera with phenotypes that differ from those of their parental genes.

Gene duplication stimulates evolutionary innovation as the resulting paralogs acquire mutations that lead to sub- or neofunctionalization. A comprehensive in silico analysis of paralogs in Saccharomyces cerevisiae reveals that duplicates of cell-surface and subtelomeric genes also undergo ectopic recombination, which leads to new chimaeric alleles. Mimicking such intergenic recombination events in the FLO (flocculation) family of cell-surface genes shows that chimaeric FLO alleles confer different adhesion phenotypes than the parental genes. Our results indicate that intergenic recombination between paralogs can generate a large set of new alleles, thereby providing the raw material for evolutionary adaptation and innovation.

Maintenance of very long telomeres by recombination in the Kluyveromyces lactis stn1-M1 mutant involves extreme telomeric turnover, telomeric circles, and concerted telomeric amplification

Some cancers utilize the recombination-dependent process of alternative lengthening of telomeres (ALT) to maintain long heterogeneous telomeres. Here, we studied the recombinational telomere elongation (RTE) of the Kluyveromyces lactis stn1-M1 mutant. We found that the total amount of the abundant telomeric DNA in stn1-M1 cells is subject to rapid variation and that it is likely to be primarily extrachromosomal. Rad50 and Rad51, known to be required for different RTE pathways in Saccharomyces cerevisiae, were not essential for the production of either long telomeres or telomeric circles in stn1-M1 cells. Circles of DNA containing telomeric repeats (t-circles) either present at the point of establishment of long telomeres or introduced later into stn1-M1 cells each led to the formation of long tandem arrays of the t-circle's sequence, which were incorporated at multiple telomeres. These tandem arrays were extraordinarily unstable and showed evidence of repeated rounds of concerted amplification. Our results suggest that the maintenance of telomeres in the stn1-M1 mutant involves extreme turnover of telomeric sequences from processes including both large deletions and the copying of t-circles.

Intermittency as a universal characteristic of the complete chromosome DNA sequences of eukaryotes: from protozoa to human genomes

Large-scale dynamical properties of complete chromosome DNA sequences of eukaryotes are considered. Using the proposed deterministic models with intermittency and symbolic dynamics we describe a wide spectrum of large-scale patterns inherent in these sequences, such as segmental duplications, tandem repeats, and other complex sequence structures. It is shown that the recently discovered gene number balance on the strands is not of a random nature, and certain subsystems of a complete chromosome DNA sequence exhibit the properties of deterministic chaos.

The duplicated deacetylases Sir2 and Hst1 subfunctionalized by acquiring complementary inactivating mutations

Protein families are generated by successive rounds of gene duplication and subsequent diversification. However, the paths by which duplicated genes acquire distinct functions are not well characterized. We focused on a pair of duplicated deacetylases from Saccharomyces cerevisiae, Sir2 and Hst1, that subfunctionalized after duplication. As a proxy for the ancestral, nonduplicated deacetylase, we studied Sir2 from another yeast, Kluyveromyces lactis. We compared the interaction domains of these deacetylases for the Sir transcriptional silencing complex, which acts with ScSir2, and the Sum1 repressor, which acts with ScHst1, and found that these interaction domains have been retained over the course of evolution and can be disrupted by simple amino acid substitutions. Therefore, Sir2 and Hst1 subfunctionalized by acquiring complementary inactivating mutations in these interaction domains.

The hidden duplication past of the plant pathogen Phytophthora and its consequences for infection

BACKGROUND

Oomycetes of the genus Phytophthora are pathogens that infect a wide range of plant species. For dicot hosts such as tomato, potato and soybean, Phytophthora is even the most important pathogen. Previous analyses of Phytophthora genomes uncovered many genes, large gene families and large genome sizes that can partially be explained by significant repeat expansion patterns.

RESULTS

Analysis of the complete genomes of three different Phytophthora species, using a newly developed approach, unveiled a large number of small duplicated blocks, mainly consisting of two or three consecutive genes. Further analysis of these duplicated genes and comparison with the known gene and genome duplication history of ten other eukaryotes including parasites, algae, plants, fungi, vertebrates and invertebrates, suggests that the ancestor of P. infestans, P. sojae and P. ramorum most likely underwent a whole genome duplication (WGD). Genes that have survived in duplicate are mainly genes that are known to be preferentially retained following WGDs, but also genes important for pathogenicity and infection of the different hosts seem to have been retained in excess. As a result, the WGD might have contributed to the evolutionary and pathogenic success of Phytophthora.

CONCLUSIONS

The fact that we find many small blocks of duplicated genes indicates that the genomes of Phytophthora species have been heavily rearranged following the WGD. Most likely, the high repeat content in these genomes have played an important role in this rearrangement process. As a consequence, the paucity of retained larger duplicated blocks has greatly complicated previous attempts to detect remnants of a large-scale duplication event in Phytophthora. However, as we show here, our newly developed strategy to identify very small duplicated blocks might be a useful approach to uncover ancient polyploidy events, in particular for heavily rearranged genomes.

Rapid expansion and functional divergence of subtelomeric gene families in yeasts

BACKGROUND

Subtelomeres, regions proximal to telomeres, exhibit characteristics unique to eukaryotic genomes. Genes residing in these loci are subject to epigenetic regulation and elevated rates of both meiotic and mitotic recombination. However, most genome sequences do not contain assembled subtelomeric sequences, and, as a result, subtelomeres are often overlooked in comparative genomics.

RESULTS

We studied the evolution and functional divergence of subtelomeric gene families in the yeast lineage. Our computational results show that subtelomeric families are evolving and expanding much faster than families that do not contain subtelomeric genes. Focusing on three related subtelomeric MAL gene families involved in disaccharide metabolism that show typical patterns of rapid expansion and evolution, we show experimentally how frequent duplication events followed by functional divergence yield novel alleles that allow the metabolism of different carbohydrates.

CONCLUSIONS

Taken together, our computational and experimental analyses show that the extraordinary instability of eukaryotic subtelomeres supports rapid adaptation to novel niches by promoting gene recombination and duplication followed by functional divergence of the alleles.

Origin and fate of pseudogenes in Hemiascomycetes: a comparative analysis

Background

Pseudogenes are ubiquitous genetic elements that derive from functional genes after mutational inactivation. Characterization of pseudogenes is important to understand genome dynamics and evolution, and its significance increases when several genomes of related organisms can be compared. Among yeasts, only the genome of the S. cerevisiae reference strain has been analyzed so far for pseudogenes.

Results

We present here the first comparative analysis of pseudogenes within the fully sequenced and annotated genomes of eight yeast species, spanning the entire phylogenetic range of Hemiascomycetes. A total of 871 pseudogenes were found, out of which mutational degradation patterns and consequences on the genetic repertoire of each species could be identified. We found that most pseudogenes in yeasts originate from mutational degradation of gene copies formed after species-specific duplications but duplications of pseudogenes themselves are also encountered. In all yeasts, except in Y. lipolytica, pseudogenes tend to cluster in subtelomeric regions where they can outnumber the number of functional genes from 3 to 16 times. Pseudogenes are generally not conserved between the yeast species studied (except in two cases), consistent with their large evolutionary distances, but tend to be conserved among S. cerevisiae strains. Reiterated pseudogenization of some genes is often observed in different lineages and may affect functions essential in S. cerevisiae, which are, therefore, lost in other species. Although a variety of functions are affected by pseudogenization, there is a bias towards functions involved in the adaptation of the yeasts to their environment, and towards genes of unknown functions.

Conclusions

Our work illustrates for the first time the formation of pseudogenes in different branches of hemiascomycetous yeasts, showing their limited conservation and how they testify for the adaptation of the yeasts functional repertoires.

Dynamic evolution of megasatellites in yeasts

Megasatellites are a new family of long tandem repeats, recently discovered in the yeast Candida glabrata. Compared to shorter tandem repeats, such as minisatellites, megasatellite motifs range in size from 135 to more than 300 bp, and allow calculation of evolutionary distances between individual motifs. Using divergence based on nucleotide substitutions among similar motifs, we determined the smallest distance between two motifs, allowing their subsequent clustering. Motifs belonging to the same cluster are recurrently found in different megasatellites located on different chromosomes, showing transfer of genetic information between megasatellites. In comparison, evolution of the few similar tandem repeats in Saccharomyces cerevisiae FLO genes mainly involves subtelomeric homologous recombination. We estimated selective constraints acting on megasatellite motifs and their host genes, and found that motifs are under strong purifying selection. Surprisingly, motifs inserted within pseudogenes are also under purifying selection, whereas the pseudogenes themselves evolve neutrally. We propose that megasatellite motifs propagate by a combination of three different molecular mechanisms: (i) gene duplication, (ii) ectopic homologous recombination and (iii) transfer of motifs from one megasatellite to another one. These mechanisms actively cooperate to create new megasatellites, that may play an important role in the adaptation of Candida glabrata to its human host.

Characterization and expression analysis of a maltose-utilizing (MAL) cluster in Aspergillus oryzae

Starch and maltooligosaccharides such as maltose and maltotriose induce the production of amylolytic enzymes including alpha-amylase in Aspergillus oryzae. A transcriptional activator gene amyR, required for maltose induction of amylolytic enzymes, has been cloned and characterized. The amyR gene deletion mutant showed significantly poor growth on starch medium but normal growth on maltose medium. This indicated the existence of another maltose-utilizing system, whose expression might not be controlled by amyR. We have identified a gene cluster homologous to the MAL cluster of Saccharomyces cerevisiae in the A. oryzae genome. The cluster consists of a MAL61 homolog (designated malP), a MAL62 homolog (designated malT), and a MAL63 homolog (designated malR). Overexpression of malT in A. oryzae resulted in a significant increase in intracellular alpha-glucosidase activity, and that of malP allowed S. cerevisiaemal61Delta to grow on maltose. The expression of both malP and malT genes was highly up-regulated in the presence of maltose, but malR expressed constitutively irrespective of carbon sources. Disruption of malR resulted in the loss of malP and malT expression and thus in restricted growth on maltose medium. In addition, a malP disruptant showed a significantly reduced expression of malT and malR and exhibited a growth defect on maltose similar to the malR disruptant. These results suggest that the MAL cluster of A. oryzae is responsible for the assimilation of maltose in A. oryzae.

Combined phylogeny and neighborhood analysis of the evolution of the ABC transporters conferring multiple drug resistance in hemiascomycete yeasts

Background

Pleiotropic Drug Resistant transporters (PDR) are members of the ATP-Binding Cassette (ABC) subfamily which export antifungals and other xenobiotics in fungi and plants. This subfamily of transmembrane transporters has nine known members in Saccharomyces cerevisiae.

We have analyzed the complex evolution of the pleiotropic drug resistance proteins (Pdrp) subfamily where gene duplications and deletions occur independently in individual genomes.

This study was carried out on 62 Pdrp from nine hemiascomycetous species, seven of which span 6 of the 14 clades of the Saccharomyces complex while the two others species, Debaryomyces hansenii and Yarrowia lipolytica, are further apart from an evolutive point of view.

Results

Combined phylogenetic and neighborhood analyses enabled us to identify five Pdrp clusters in the Saccharomyces complex. Three of them comprise orthologs of the Pdrp sensu stricto, Pdr5p, Pdr10p, Pdr12p, Pdr15p, Snq2p and YNR070wp. The evolutive pathway of the orthologs of Snq2 and YNR070w is particularly complex due to a tandem gene array in Eremothecium gossypii, Kluyveromyces lactis and Saccharomyces (Lachancea) kluyveri. This pathway and different cases of duplications and deletions were clarified by using a neighborhood analysis based on synteny.

For the two distant species, Yarrowia lipolytica and Debaryomyces hansenii, no neighborhood evidence is available for these clusters and many homologs of Pdr5 and Pdr15 are phylogenetically assigned to species-based clusters.

Two other clusters comprise the orthologs of the sensu lato Pdrp, Aus1p/Pdr11p and YOL075cp respectively. The evolutionary pathway of these clusters is simpler. Nevertheless, orthologs of these genes are missing in some species.

Conclusion

Numerous duplications were traced among the Hemiascomycetous Pdrp studied. The role of the Whole Genome Duplication (WGD) is sorted out and our analyses confirm the common ancestrality of Pdr5p and Pdr15p. A tandem gene array is observed in Eremothecium gossypii. One of the copies is the ortholog of Snq2 while the other one is lost in the post-WGD species.

The neighborhood analysis provides an efficient method to trace the history of genes and disentangle the orthology and paralogy relationships.

Hidden chromosome symmetry: in silico transformation reveals symmetry in 2D DNA walk trajectories of 671 chromosomes

Maps of 2D DNA walk of 671 examined chromosomes show composition complexity change from symmetrical half-turn in bacteria to pseudo-random trajectories in archaea, fungi and humans. In silico transformation of gene order and strand position returns most of the analyzed chromosomes to a symmetrical bacterial-like state with one transition point. The transformed chromosomal sequences also reveal remarkable segmental compositional symmetry between regions from different strands located equidistantly from the transition point. Despite extensive chromosome rearrangement the relation of gene numbers on opposite strands for chromosomes of different taxa varies in narrow limits around unity with Pearson coefficient r = 0.98. Similar relation is observed for total genes' length (r = 0.86) and cumulative GC (r = 0.95) and AT (r = 0.97) skews. This is also true for human coding sequences (CDS), which comprise only several percent of the entire chromosome length. We found that frequency distributions of the length of gene clusters, continuously located on the same strand, have close values for both strands. Eukaryotic gene distribution is believed to be non-random. Contribution of different subsystems to the noted symmetries and distributions, and evolutionary aspects of symmetry are discussed.

Comparative genomics of protoploid Saccharomycetaceae

Our knowledge of yeast genomes remains largely dominated by the extensive studies on Saccharomyces cerevisiae and the consequences of its ancestral duplication, leaving the evolution of the entire class of hemiascomycetes only partly explored. We concentrate here on five species of Saccharomycetaceae, a large subdivision of hemiascomycetes, that we call "protoploid" because they diverged from the S. cerevisiae lineage prior to its genome duplication. We determined the complete genome sequences of three of these species: Kluyveromyces (Lachancea) thermotolerans and Saccharomyces (Lachancea) kluyveri (two members of the newly described Lachancea clade), and Zygosaccharomyces rouxii. We included in our comparisons the previously available sequences of Kluyveromyces lactis and Ashbya (Eremothecium) gossypii. Despite their broad evolutionary range and significant individual variations in each lineage, the five protoploid Saccharomycetaceae share a core repertoire of approximately 3300 protein families and a high degree of conserved synteny. Synteny blocks were used to define gene orthology and to infer ancestors. Far from representing minimal genomes without redundancy, the five protoploid yeasts contain numerous copies of paralogous genes, either dispersed or in tandem arrays, that, altogether, constitute a third of each genome. Ancient, conserved paralogs as well as novel, lineage-specific paralogs were identified.

Additions, losses, and rearrangements on the evolutionary route from a reconstructed ancestor to the modern Saccharomyces cerevisiae genome

Comparative genomics can be used to infer the history of genomic rearrangements that occurred during the evolution of a species. We used the principle of parsimony, applied to aligned synteny blocks from 11 yeast species, to infer the gene content and gene order that existed in the genome of an extinct ancestral yeast about 100 Mya, immediately before it underwent whole-genome duplication (WGD). The reconstructed ancestral genome contains 4,703 ordered loci on eight chromosomes. The reconstruction is complete except for the subtelomeric regions. We then inferred the series of rearrangement steps that led from this ancestor to the current Saccharomyces cerevisiae genome; relative to the ancestral genome we observe 73 inversions, 66 reciprocal translocations, and five translocations involving telomeres. Some fragile chromosomal sites were reused as evolutionary breakpoints multiple times. We identified 124 genes that have been gained by S. cerevisiae in the time since the WGD, including one that is derived from a hAT family transposon, and 88 ancestral loci at which S. cerevisiae did not retain either of the gene copies that were formed by WGD. Sites of gene gain and evolutionary breakpoints both tend to be associated with tRNA genes and, to a lesser extent, with origins of replication. Many of the gained genes in S. cerevisiae have functions associated with ethanol production, growth in hypoxic environments, or the uptake of alternative nutrient sources.

Regulation of secondary metabolite production in filamentous ascomycetes

Fungi are renowned for their ability to produce bioactive small molecules otherwise known as secondary metabolites. These molecules have attracted much attention due to both detrimental (e.g. toxins) and beneficial (e.g. pharmaceuticals) effects on human endeavors. Once the topic only of chemical and biochemical studies, secondary metabolism research has reached a sophisticated level in the realm of genetic regulation. This review covers the latest insights into the processes regulating secondary metabolite production in filamentous fungi.

Histone deacetylase activity regulates chemical diversity in Aspergillus

Bioactive small molecules are critical in Aspergillus species during their development and interaction with other organisms. Genes dedicated to their production are encoded in clusters that can be located throughout the genome. We show that deletion of hdaA, encoding an Aspergillus nidulans histone deacetylase (HDAC), causes transcriptional activation of two telomere-proximal gene clusters--and subsequent increased levels of the corresponding molecules (toxin and antibiotic)--but not of a telomere-distal cluster. Introduction of two additional HDAC mutant alleles in a DeltahdaA background had minimal effects on expression of the two HdaA-regulated clusters. Treatment of other fungal genera with HDAC inhibitors resulted in overproduction of several metabolites, suggesting a conserved mechanism of HDAC repression of some secondary-metabolite gene clusters. Chromatin regulation of small-molecule gene clusters may enable filamentous fungi to successfully exploit environmental resources by modifying chemical diversity.

Analysis of the VSG gene silent archive in Trypanosoma brucei reveals that mosaic gene expression is prominent in antigenic variation and is favored by archive substructure

Trypanosoma brucei evades host acquired immunity through differential activation of its large archive of silent variant surface glycoprotein (VSG) genes, most of which are pseudogenes in subtelomeric arrays. We have analyzed 940 VSGs, representing one half to two thirds of the arrays. Sequence types A and B of the VSG N-terminal domains were confirmed, while type C was found to be a constituent of type A. Two new C-terminal domain types were found. Nearly all combinations of domain types occurred, with some bias to particular combinations. One-third of encoded N-terminal domains, but only 13% of C-terminal domains, are intact, indicating a particular need for silent VSGs to gain a functional C-terminal domain to be expressed. About 60% of VSGs are unique, the rest occurring in subfamilies of two to four close homologs (>50%-52% peptide identity). We found a subset of VSG-related genes, differing from VSGs in genomic environment and expression patterns, and predict they have distinct function. Almost all (92%) full-length array VSGs have the partially conserved flanks associated with the duplication mechanism that activates silent genes, and these sequences have also contributed to archive evolution, mediating most of the conversions of segments, containing >/=1 VSG, within and between arrays. During infection, intact array genes became activated by duplication after two weeks, and mosaic VSGs assembled from pseudogenes became expressed by week three and predominated by week four. The small subfamily structure of the archive appears to be fundamental in providing the interacting donors for mosaic formation.