Reconstruction of ancestral chromosome architecture and gene repertoire reveals principles of genome evolution in a model yeast genus.

Reconstructing genome history is complex but necessary to reveal quantitative principles governing genome evolution. Such reconstruction requires recapitulating into a single evolutionary framework the evolution of genome architecture and gene repertoire. Here, we reconstructed the genome history of the genus Lachancea that appeared to cover a continuous evolutionary range from closely related to more diverged yeast species. Our approach integrated the generation of a high-quality genome data set; the development of AnChro, a new algorithm for reconstructing ancestral genome architecture; and a comprehensive analysis of gene repertoire evolution. We found that the ancestral genome of the genus Lachancea contained eight chromosomes and about 5173 protein-coding genes. Moreover, we characterized 24 horizontal gene transfers and 159 putative gene creation events that punctuated species diversification. We retraced all chromosomal rearrangements, including gene losses, gene duplications, chromosomal inversions and translocations at single gene resolution. Gene duplications outnumbered losses and balanced rearrangements with 1503, 929, and 423 events, respectively. Gene content variations between extant species are mainly driven by differential gene losses, while gene duplications remained globally constant in all lineages. Remarkably, we discovered that balanced chromosomal rearrangements could be responsible for up to 14% of all gene losses by disrupting genes at their breakpoints. Finally, we found that nonsynonymous substitutions reached fixation at a coordinated pace with chromosomal inversions, translocations, and duplications, but not deletions. Overall, we provide a granular view of genome evolution within an entire eukaryotic genus, linking gene content, chromosome rearrangements, and protein divergence into a single evolutionary framework.

Population genomics reveals chromosome-scale heterogeneous evolution in a protoploid yeast.

Yeast species represent an ideal model system for population genomic studies but large-scale polymorphism surveys have only been reported for species of the Saccharomyces genus so far. Hence, little is known about intraspecific diversity and evolution in yeast. To obtain a new insight into the evolutionary forces shaping natural populations, we sequenced the genomes of an expansive worldwide collection of isolates from a species distantly related to Saccharomyces cerevisiae: Lachancea kluyveri (formerly S. kluyveri). We identified 6.5 million single nucleotide polymorphisms and showed that a large introgression event of 1 Mb of GC-rich sequence in the chromosomal arm probably occurred in the last common ancestor of all L. kluyveri strains. Our population genomic data clearly revealed that this 1-Mb region underwent a molecular evolution pattern very different from the rest of the genome. It is characterized by a higher recombination rate, with a dramatically elevated A:T → G:C substitution rate, which is the signature of an increased GC-biased gene conversion. In addition, the predicted base composition at equilibrium demonstrates that the chromosome-scale compositional heterogeneity will persist after the genome has reached mutational equilibrium. Altogether, the data presented herein clearly show that distinct recombination and substitution regimes can coexist and lead to different evolutionary patterns within a single genome.

Genetic Basis of Ammonium Toxicity Resistance in a Sake Strain of Yeast: A Mendelian Case.

High concentrations of ammonium at physiological concentrations of potassium are toxic for the standard laboratory strain of Saccharomyces cerevisiae In the original description of this metabolic phenotype, we focused on the standard laboratory strains of Saccharomyces In this study, we screened a large collection of S. cerevisiae natural isolates and identified one strain that is resistant to high concentrations of ammonium. This strain, K12, was isolated in sake breweries. When the K12 strain was crossed to the standard laboratory strain (FY4), the resulting tetrads displayed 2:2 segregation of the resistance phenotype, suggesting a single gene trait. Using a bulk segregant analysis strategy, we mapped this trait to a 150-kb region on chromosome X containing the TRK1 gene. This gene encodes a transporter required for high-affinity potassium transport in S. cerevisiae Data from reciprocal hemizygosity experiments with TRK1 deletion strains in K12 and BY backgrounds, as well as analysis of the deletion of this gene in the K12 strain, demonstrate that the K12 allele of TRK1 is responsible for ammonium toxicity resistance. Furthermore, we determined the minimal amount of potassium required for both the K12 and laboratory strain needed for growth. These results demonstrate that the gene encoded by the K12 allele of TRK1 has a greater affinity for potassium than the standard allele of TRK1 found in Saccharomyces strains. We hypothesize that this greater-affinity allele of the potassium transporter reduces the flux of ammonium into the yeast cells under conditions of ammonium toxicity. These findings further refine our understanding of ammonium toxicity in yeast and provide an example of using natural variation to understand cellular processes.

Mitochondrial genome evolution in a single protoploid yeast species.

Mitochondria are organelles, which play a key role in some essential functions, including respiration, metabolite biosynthesis, ion homeostasis, and apoptosis. The vast numbers of mitochondrial DNA (mtDNA) sequences of various yeast species, which have recently been published, have also helped to elucidate the structural diversity of these genomes. Although a large corpus of data are now available on the diversity of yeast species, little is known so far about the mtDNA diversity in single yeast species. To study the genetic variations occurring in the mtDNA of wild yeast isolates, we performed a genome-wide polymorphism survey on the mtDNA of 18 Lachancea kluyveri (formerly Saccharomyces kluyveri) strains. We determined the complete mt genome sequences of strains isolated from various geographical locations (in North America, Asia, and Europe) and ecological niches (Drosophila, tree exudates, soil). The mt genome of the NCYC 543 reference strain is 51,525 bp long. It contains the same core of genes as Lachancea thermotolerans, the nearest relative to L. kluyveri. To explore the mt genome variations in a single yeast species, we compared the mtDNAs of the 18 isolates. The phylogeny and population structure of L. kluyveri provide clear-cut evidence for the existence of well-defined geographically isolated lineages. Although these genomes are completely syntenic, their size and the intron content were found to vary among the isolates studied. These genomes are highly polymorphic, showing an average diversity of 28.5 SNPs/kb and 6.6 indels/kb. Analysis of the SNP and indel patterns showed the existence of a particularly high overall level of polymorphism in the intergenic regions. The dN/dS ratios obtained are consistent with purifying selection in all these genes, with the noteworthy exception of the VAR1 gene, which gave a very high ratio. These data suggest that the intergenic regions have evolved very fast in yeast mitochondrial genomes.