Pan-transcriptome reveals a large accessory genome contribution to gene expression variation in yeast.

Gene expression is an essential step in the translation of genotypes into phenotypes. However, little is known about the transcriptome architecture and the underlying genetic effects at the species level. Here we generated and analyzed the pan-transcriptome of ~1,000 yeast natural isolates across 4,977 core and 1,468 accessory genes. We found that the accessory genome is an underappreciated driver of transcriptome divergence. Global gene expression patterns combined with population structure showed that variation in heritable expression mainly lies within subpopulation-specific signatures, for which accessory genes are overrepresented. Genome-wide association analyses consistently highlighted that accessory genes are associated with proportionally more variants with larger effect sizes, illustrating the critical role of the accessory genome on the transcriptional landscape within and between populations.

Telomere-to-telomere assemblies of 142 strains characterize the genome structural landscape in Saccharomyces cerevisiae.

Pangenomes provide access to an accurate representation of the genetic diversity of species, both in terms of sequence polymorphisms and structural variants (SVs). Here we generated the Saccharomyces cerevisiae Reference Assembly Panel (ScRAP) comprising reference-quality genomes for 142 strains representing the species' phylogenetic and ecological diversity. The ScRAP includes phased haplotype assemblies for several heterozygous diploid and polyploid isolates. We identified circa (ca.) 4,800 nonredundant SVs that provide a broad view of the genomic diversity, including the dynamics of telomere length and transposable elements. We uncovered frequent cases of complex aneuploidies where large chromosomes underwent large deletions and translocations. We found that SVs can impact gene expression near the breakpoints and substantially contribute to gene repertoire evolution. We also discovered that horizontally acquired regions insert at chromosome ends and can generate new telomeres. Overall, the ScRAP demonstrates the benefit of a pangenome in understanding genome evolution at population scale.

Contrasting Genomic Evolution Between Domesticated and Wild Kluyveromyces lactis Yeast Populations.

The process of domestication has variable consequences on genome evolution leading to different phenotypic signatures. Access to the complete genome sequences of a large number of individuals makes it possible to explore the different facets of this domestication process. Here, we sought to explore the genome evolution of Kluyveromyces lactis, a yeast species well known for its involvement in dairy processes and also present in natural environments. Using a combination of short- and long-read sequencing strategies, we investigated the genomic variability of 41 K. lactis isolates and found that the overall genetic diversity of this species is very high (θw = 3.3 × 10-2) compared with other species such as Saccharomyces cerevisiae (θw = 1.6 × 10-2). However, the domesticated dairy population shows a reduced level of diversity (θw = 1 × 10-3), probably due to a domestication bottleneck. In addition, this entire population is characterized by the introgression of the LAC4 and LAC12 genes, responsible for lactose fermentation and coming from the closely related species, Kluyveromyces marxianus, as previously described. Our results highlighted that the LAC4/LAC12 gene cluster was acquired through multiple and independent introgression events. Finally, we also identified several genes that could play a role in adaptation to dairy environments through copy number variation. These genes are involved in sugar consumption, flocculation, and drug resistance, and may play a role in dairy processes. Overall, our study illustrates contrasting genomic evolution and sheds new light on the impact of domestication processes on it.

Dissection of quantitative trait loci in the Lachancea waltii yeast species highlights major hotspots.

Dissecting the genetic basis of complex trait remains a real challenge. The budding yeast Saccharomyces cerevisiae has become a model organism for studying quantitative traits, successfully increasing our knowledge in many aspects. However, the exploration of the genotype-phenotype relationship in non-model yeast species could provide a deeper insight into the genetic basis of complex traits. Here, we have studied this relationship in the Lachancea waltii species which diverged from the S. cerevisiae lineage prior to the whole-genome duplication. By performing linkage mapping analyses in this species, we identified 86 quantitative trait loci (QTL) impacting the growth in a large number of conditions. The distribution of these loci across the genome has revealed two major QTL hotspots. A first hotspot corresponds to a general growth QTL, impacting a wide range of conditions. By contrast, the second hotspot highlighted a trade-off with a disadvantageous allele for drug-free conditions which proved to be advantageous in the presence of several drugs. Finally, a comparison of the detected QTL in L. waltii with those which had been previously identified for the same trait in a closely related species, Lachancea kluyveri was performed. This analysis clearly showed the absence of shared QTL across these species. Altogether, our results represent a first step toward the exploration of the genetic architecture of quantitative trait across different yeast species.

Pervasive Phenotypic Impact of a Large Nonrecombining Introgressed Region in Yeast.

To explore the origin of the diversity observed in natural populations, many studies have investigated the relationship between genotype and phenotype. In yeast species, especially in Saccharomyces cerevisiae, these studies are mainly conducted using recombinant offspring derived from two genetically diverse isolates, allowing to define the phenotypic effect of genetic variants. However, large genomic variants such as interspecies introgressions are usually overlooked even if they are known to modify the genotype-phenotype relationship. To have a better insight into the overall phenotypic impact of introgressions, we took advantage of the presence of a 1-Mb introgressed region, which lacks recombination and contains the mating-type determinant in the Lachancea kluyveri budding yeast. By performing linkage mapping analyses in this species, we identified a total of 89 loci affecting growth fitness in a large number of conditions and 2,187 loci affecting gene expression mostly grouped into two major hotspots, one being the introgressed region carrying the mating-type locus. Because of the absence of recombination, our results highlight the presence of a sexual dimorphism in a budding yeast for the first time. Overall, by describing the phenotype-genotype relationship in the Lachancea kluyveri species, we expanded our knowledge on how genetic characteristics of large introgression events can affect the phenotypic landscape.

Variation of the meiotic recombination landscape and properties over a broad evolutionary distance in yeasts.

Meiotic recombination is a major factor of genome evolution, deeply characterized in only a few model species, notably the yeast Saccharomyces cerevisiae. Consequently, little is known about variations of its properties across species. In this respect, we explored the recombination landscape of Lachancea kluyveri, a protoploid yeast species that diverged from the Saccharomyces genus more than 100 million years ago and we found striking differences with S. cerevisiae. These variations include a lower recombination rate, a higher frequency of chromosomes segregating without any crossover and the absence of recombination on the chromosome arm containing the sex locus. In addition, although well conserved within the Saccharomyces clade, the S. cerevisiae recombination hotspots are not conserved over a broader evolutionary distance. Finally and strikingly, we found evidence of frequent reversal of commitment to meiosis, resulting in return to mitotic growth after allele shuffling. Identification of this major but underestimated evolutionary phenomenon illustrates the relevance of exploring non-model species.

de novo assembly and population genomic survey of natural yeast isolates with the Oxford Nanopore MinION sequencer.

BACKGROUND: Oxford Nanopore Technologies Ltd (Oxford, UK) have recently commercialized MinION, a small single-molecule nanopore sequencer, that offers the possibility of sequencing long DNA fragments from small genomes in a matter of seconds. The Oxford Nanopore technology is truly disruptive; it has the potential to revolutionize genomic applications due to its portability, low cost, and ease of use compared with existing long reads sequencing technologies. The MinION sequencer enables the rapid sequencing of small eukaryotic genomes, such as the yeast genome. Combined with existing assembler algorithms, near complete genome assemblies can be generated and comprehensive population genomic analyses can be performed. RESULTS: Here, we resequenced the genome of the Saccharomyces cerevisiae S288C strain to evaluate the performance of nanopore-only assemblers. Then we de novo sequenced and assembled the genomes of 21 isolates representative of the S. cerevisiae genetic diversity using the MinION platform. The contiguity of our assemblies was 14 times higher than the Illumina-only assemblies and we obtained one or two long contigs for 65 % of the chromosomes. This high contiguity allowed us to accurately detect large structural variations across the 21 studied genomes. CONCLUSION: Because of the high completeness of the nanopore assemblies, we were able to produce a complete cartography of transposable elements insertions and inspect structural variants that are generally missed using a short-read sequencing strategy. Our analyses show that the Oxford Nanopore technology is already usable for de novo sequencing and assembly; however, non-random errors in homopolymers require polishing the consensus using an alternate sequencing technology.

Dissection of quantitative traits by bulk segregant mapping in a protoploid yeast species.

Since more than a decade ago, Saccharomyces cerevisiae has been used as a model to dissect complex traits, revealing the genetic basis of a large number of traits in fine detail. However, to have a more global view of the genetic architecture of traits across species, the examination of the molecular basis of phenotypes within non-conventional species would undoubtedly be valuable. In this respect, the Saccharomycotina yeasts represent ideal and potential non-model organisms. Here we sought to assess the feasibility of genetic mapping by bulk segregant analysis in the protoploid Lachancea kluyveri (formerly S. kluyveri) yeast species, a distantly related species to S. cerevisiae For this purpose, we designed a fluorescent mating-type marker, compatible with any mating-competent strains representative of this species, to rapidly create a large population of haploid segregants (>10(5) cells). Quantitative trait loci can be mapped by selecting and sequencing an enriched pool of progeny with extreme phenotypic values. As a test bed, we applied this strategy and mapped the causal loci underlying halotolerance phenotypes in L. kluyveri Overall, this study demonstrates that bulk segregant mapping is a powerful way for investigating the genetic basis of natural variations in non-model yeast organisms and more precisely in L. kluyveri.

Genetic background determines metabolic phenotypes in the mouse.

To evaluate the contribution of genetic background to phenotypic variation, we compared a large range of biochemical and metabolic parameters at different ages of four inbred mice strains, C57BL/6J, 129SvPas, C3HeB/FeJ, and Balb/cByJ. Our results demonstrate that important metabolic, hematologic, and biochemical differences exist between these different inbred strains. Most of these differences are gender independent and are maintained or accentuated throughout life. It is therefore imperative that the genetic background is carefully defined in phenotypic studies. Our results also argue that certain backgrounds are more suited to study a given physiologic phenomenon, as distinct mouse strains have a different propensity to develop particular biochemical, hematologic, and metabolic abnormalities. These genetic differences can furthermore be exploited to identify new genes/proteins that contribute to phenotypic abnormalities. The choice of the genetic background in which to generate and analyze genetically engineered mutant mice is important as it is, together with environmental factors, one of the most important contributors to the variability of phenotypic results.

Dissociation of analgesic and hormonal responses to forced swim stress using opioid receptor knockout mice.

Exposure to stress triggers hormonal and behavioral responses. It has been shown that the endogenous opioid system plays a role in some physiological reactions to stress. The opioid system was described to mediate analgesia induced by mild stressors and to modulate the activation of the hypothalamic-pituitary-adrenal axis. Our study assessed the contribution of opioid receptors in stress-induced analgesia and adrenocorticotropic hormone (ACTH) and corticosterone release by a genetic approach. We performed a parallel analysis of mice deficient in mu, delta, or kappa opioid receptors, as well as of triple opioid receptor knockout mice, following exposure to a mild stress (3-min swim at 32 degrees C). In wild-type mice, stress elicited an increase in jumping latency on the hot plate, which was influenced by gender and genetic background. This analgesic response was reversed both by naloxone and by the triple mutation, and decreased in mu and delta opioid receptor knockout females. In wild-type females, stress also delayed front- and hindpaw behaviors in the hot plate test and increased tail-flick latency in the tail immersion test. Opioid receptor deletion however did not affect these stress responses. In addition, stress produced an increase in ACTH and corticosterone plasma levels. This endocrine response remained unchanged in all mutant strains. Therefore our data indicate that, under our stress conditions, the endogenous opioid system is recruited to produce some analgesia whereas it does not influence hypothalamic-pituitary-adrenal axis activity. This implies that brain circuits mediating analgesic and hormonal responses to stress can be dissociated.

Mouse functional genomics requires standardization of mouse handling and housing conditions.

The study of mouse models is crucial for the functional annotation of the human genome. The recent improvements in mouse genetics now moved the bottleneck in mouse functional genomics from the generation of mutant mice lines to the phenotypic analysis of these mice lines. Simple, validated, and reproducible phenotyping tests are a prerequisite to improving this phenotyping bottleneck. We analyzed here the impact of simple variations in animal handling and housing procedures, such as cage density, diet, gender, length of fasting, as well as site (retro-orbital vs. tail), timing, and anesthesia used during venipuncture, on biochemical, hematological, and metabolic/endocrine parameters in adult C57BL/6J mice. Our results, which show that minor changes in procedures can profoundly affect biological variables, underscore the importance of establishing uniform and validated animal procedures to improve reproducibility of mouse phenotypic data.