Sourdough yeast-bacteria interactions can change ferulic acid metabolism during fermentation.

The metabolism of ferulic acid (FA) was studied during fermentation with different species and strains of lactic acid bacteria (LAB) and yeasts, in synthetic sourdough medium. Yeast strains of Kazachstania humilis, Kazachstania bulderi, and Saccharomyces cerevisiae, as well as lactic acid bacteria strains of Fructilactobacillus sanfranciscensis, Lactiplantibacillus plantarum, Lactiplantibacillus xiangfangensis, Levilactobacillus hammesii, Latilactobacillus curvatus and Latilactobacillus sakei were selected from French natural sourdoughs. Fermentation in presence or absence of FA was carried out in LAB and yeasts monocultures, as well as in LAB/yeast co-cultures. Our results indicated that FA was mainly metabolized into 4-vinylguaiacol (4-VG) by S. cerevisiae strains, and into dihydroferulic acid (DHFA) and 4-VG in the case of LAB. Interactions of LAB and yeasts led to the modification of FA metabolism, with a major formation of DHFA, even by the strains that do not produce it in monoculture. Interestingly, FA was almost completely consumed by the F. sanfranciscensis bFs17 and K. humilis yKh17 pair and converted into DHFA in 89.5 ± 19.6% yield, while neither bFs17, nor yKh17 strains assimilated FA in monoculture.

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.

Characterization of hexose transporters in Yarrowia lipolytica reveals new groups of Sugar Porters involved in yeast growth.

Sugar assimilation has been intensively studied in the model yeast S. cerevisiae, and for two decades, it has been clear that the homologous HXT genes, which encode a set of hexose transporters, play a central role in this process. However, in the yeast Yarrowia lipolytica, which is well-known for its biotechnological applications, sugar assimilation is only poorly understood, even though this yeast exhibits peculiar intra-strain differences in fructose uptake: some strains (e.g., W29) are known to be slow-growing in fructose while others (e.g., H222) grow rapidly under the same conditions. Here, we retrieved 24 proteins of the Sugar Porter family from these two strains, and determined that at least six of these proteins can function as hexose transporters in the heterologous host Saccharomyces cerevisiae EBY.VW4000. Transcriptional studies and deletion analysis in Y. lipolytica indicated that two genes, YHT1 and YHT4, are probably the main players in both strains, with a similar role in the uptake of glucose, fructose, and mannose at various concentrations. The other four genes appear to constitute a set of 'reservoir' hexose transporters with an as-yet unclear physiological role. Furthermore, through examining Sugar Porters of the entire Yarrowia clade, we show that they constitute a dynamic family, within which hexose transport genes have been duplicated and lost several times. Our phylogenetic analyses support the existence of at least three distinct evolutionary groups of transporters which allow yeasts to grow on hexoses. In addition to the well-known and widespread Hxt-type transporters (which are not essential in Y. lipolytica), we highlight a second group of transporters, represented by Yht1, which are phylogenetically related to sensors that play a regulatory role in S. cerevisiae, and a third group, represented by Yht4, previously thought to contain only high-affinity glucose transporters related to Hgt1of Kluyveromyces lactis.

Truncation of Gal4p explains the inactivation of the GAL/MEL regulon in both Saccharomyces bayanus and some Saccharomyces cerevisiae wine strains.

In the past, the galactose-negative (Gal(-)) phenotype was a key physiological character used to distinguish Saccharomyces bayanus from S. cerevisiae In this work, we investigated the inactivation of GAL gene networks in S. bayanus, which is an S. uvarum/S. eubayanus hybrid, and in S. cerevisiae wine strains erroneously labelled 'S. bayanus'. We made an inventory of their GAL genes using genomes that were either available publicly, re-sequenced by us, or assembled from public data and completed with targeted sequencing. In the S. eubayanus/S. uvarum CBS 380(T) hybrid, the GAL/MEL network is composed of genes from both parents: from S. uvarum, an otherwise complete set that lacks GAL4, and from S. eubayanus, a truncated version of GAL4 and an additional copy of GAL3 and GAL80 Similarly, two different truncated GAL4 alleles were found in S. cerevisiae wine strains EC1118 and LalvinQA23. The lack of GAL4 activity in these strains was corrected by introducing a full-length copy of S. cerevisiae GAL4 on a CEN4/ARS plasmid. Transformation with this plasmid restored galactose utilisation in Gal(-) strains, and melibiose fermentation in strain CBS 380(T) The melibiose fermentation phenotype, formerly regarded as characteristic of S. uvarum, turned out to be widespread among Saccharomyces species.

Comparative genomics of emerging pathogens in the Candida glabrata clade.

BACKGROUND: Candida glabrata follows C. albicans as the second or third most prevalent cause of candidemia worldwide. These two pathogenic yeasts are distantly related, C. glabrata being part of the Nakaseomyces, a group more closely related to Saccharomyces cerevisiae. Although C. glabrata was thought to be the only pathogenic Nakaseomyces, two new pathogens have recently been described within this group: C. nivariensis and C. bracarensis. To gain insight into the genomic changes underlying the emergence of virulence, we sequenced the genomes of these two, and three other non-pathogenic Nakaseomyces, and compared them to other sequenced yeasts. RESULTS: Our results indicate that the two new pathogens are more closely related to the non-pathogenic N. delphensis than to C. glabrata. We uncover duplications and accelerated evolution that specifically affected genes in the lineage preceding the group containing N. delphensis and the three pathogens, which may provide clues to the higher propensity of this group to infect humans. Finally, the number of Epa-like adhesins is specifically enriched in the pathogens, particularly in C. glabrata. CONCLUSIONS: Remarkably, some features thought to be the result of adaptation of C. glabrata to a pathogenic lifestyle, are present throughout the Nakaseomyces, indicating these are rather ancient adaptations to other environments. Phylogeny suggests that human pathogenesis evolved several times, independently within the clade. The expansion of the EPA gene family in pathogens establishes an evolutionary link between adhesion and virulence phenotypes. Our analyses thus shed light onto the relationships between virulence and the recent genomic changes that occurred within the Nakaseomyces. SEQUENCE ACCESSION NUMBERS: Nakaseomyces delphensis: CAPT01000001 to CAPT01000179Candida bracarensis: CAPU01000001 to CAPU01000251Candida nivariensis: CAPV01000001 to CAPV01000123Candida castellii: CAPW01000001 to CAPW01000101Nakaseomyces bacillisporus: CAPX01000001 to CAPX01000186.

Nonlinear SAR Modelling of Mosquito Repellents for Skin Application.

Finding new marketable mosquito repellents is a complex and time-consuming process that can be optimized via modelling. In this context, a SAR (Structure-Activity Relationship) model was designed from a set of 2171 molecules whose actual repellent activity against Aedes aegypti was available. Information-rich descriptors were used as input neurons of a three-layer perceptron (TLP) to compute the models. The most interesting classification model was a 20/6/2 TLP showing 94% and 89% accuracy on the training set and test set, respectively. A total of 57 other artificial neural network models based on the same architecture were also computed. This allowed us to consider all chemicals both as training and test set members in order to better interpret the results obtained with the selected model. Most of the wrong predictions were explainable. The 20/6/2 TLP model was then used for predicting the potential repellent activity of new molecules. Among them, two were successfully evaluated in vivo.

Insights into the Genomic and Phenotypic Landscape of the Oleaginous Yeast Yarrowia lipolytica.

Although Yarrowia lipolytica is a model yeast for the study of lipid metabolism, its diversity is poorly known, as studies generally consider only a few standard laboratory strains. To extend our knowledge of this biotechnological workhorse, we investigated the genomic and phenotypic diversity of 56 natural isolates. Y. lipolytica is classified into five clades with no correlation between clade membership and geographic or ecological origin. A low genetic diversity (π = 0.0017) and a pan-genome (6528 genes) barely different from the core genome (6315 genes) suggest Y. lipolytica is a recently evolving species. Large segmental duplications were detected, totaling 892 genes. With three new LTR-retrotransposons of the Gypsy family (Tyl4, Tyl9, and Tyl10), the transposable element content of genomes appeared diversified but still low (from 0.36% to 3.62%). We quantified 34 traits with substantial phenotypic diversity, but genome-wide association studies failed to evidence any associations. Instead, we investigated known genes and found four mutational events leading to XPR2 protease inactivation. Regarding lipid metabolism, most high-impact mutations were found in family-belonging genes, such as ALK or LIP, and therefore had a low phenotypic impact, suggesting that the huge diversity of lipid synthesis and accumulation is multifactorial or due to complex regulations.

Draft Genome Sequence of Candida railenensis Strain CLIB 1423, Isolated from Papaya Fruit in French Guiana.

Here, we report the draft genome sequence and annotation of the yeast Candida railenensis strain CLIB 1423. The assembly consists of 57 nuclear scaffolds and 1 complete mitochondrial chromosome, for a total of 13.8 Mb (N50, 0.54 Mb; L50, 9). The annotation contains 6,013 coding DNA sequences (CDSs) (BUSCO completeness, 99.6%).

Whole-Genome Sequences of Two Kazachstania barnettii Strains Isolated from Anthropic Environments.

Recent studies have suggested that species of the Kazachstania genus may be interesting models of yeast domestication. Among these, Kazachstania barnettii has been isolated from various microbially transformed foodstuffs such as sourdough bread and kefir. In the present work, we sequence, assemble, and annotate the complete genomes of two K. barnettii strains: CLIB 433, being one of the two reference strains for K. barnettii that was isolated as a spoilage organism in soft drink, and CLIB 1767, recently isolated from artisan bread-making sourdough. Both assemblies are of high quality with N50 statistics greater than 1.3 Mb and BUSCO score greater than 99%. An extensive comparison of the two obtained genomes revealed very few differences between the two K. barnettii strains, considering both genome structure and gene content. The proposed genome assemblies will constitute valuable references for future comparative genomic, population genomic, or transcriptomic studies of the K. barnettii species.

Lethal and Sublethal Effects of Pyriproxyfen on Apis and Non-Apis Bees.

Pyriproxyfen is a juvenile hormone mimic used extensively worldwide to fight pests in agriculture and horticulture. It also has numerous applications as larvicide in vector control. The molecule disrupts metamorphosis and adult emergence in the target insects. The same types of adverse effects are expected on non-target insects. In this context, the objective of this study was to evaluate the existing information on the toxicity of pyriproxyfen on the honey bee (Apis mellifera) and non-Apis bees (bumble bees, solitary bees, and stingless bees). The goal was also to identify the gaps necessary to fill. Thus, whereas the acute and sublethal toxicity of pyriproxyfen against A. mellifera is well-documented, the information is almost lacking for the non-Apis bees. The direct and indirect routes of exposure of the non-Apis bees to pyriproxyfen also need to be identified and quantified. More generally, the impacts of pyriproxyfen on the reproductive success of the different bee species have to be evaluated as well as the potential adverse effects of its metabolites.

High Complexity and Degree of Genetic Variation in Brettanomyces bruxellensis Population.

Genome-wide characterization of genetic variants of a large population of individuals within the same species is essential to have a deeper insight into its evolutionary history as well as the genotype-phenotype relationship. Population genomic surveys have been performed in multiple yeast species, including the two model organisms, Saccharomyces cerevisiae and Schizosaccharomyces pombe. In this context, we sought to characterize at the population level the Brettanomyces bruxellensis yeast species, which is a major cause of wine spoilage and can contribute to the specific flavor profile of some Belgium beers. We have completely sequenced the genome of 53 B. bruxellensis strains isolated worldwide. The annotation of the reference genome allowed us to define the gene content of this species. As previously suggested, our genomic data clearly highlighted that genetic diversity variation is related to ploidy level, which is variable in the B. bruxellensis species. Genomes are punctuated by multiple loss-of-heterozygosity regions, whereas aneuploidies as well as segmental duplications are uncommon. Interestingly, triploid genomes are more prone to gene copy number variation than diploids. Finally, the pangenome of the species was reconstructed and was found to be small with few accessory genes compared with S. cerevisiae. The pangenome is composed of 5,409 ORFs (open reading frames) among which 5,106 core ORFs and 303 ORFs that are variable within the population. All these results highlight the different trajectories of species evolution and consequently the interest of establishing population genomic surveys in more populations.

Genome Sequence of the Oleaginous Yeast Yarrowia lipolytica H222.

Here, we report the genome sequence of the oleaginous yeast Yarrowia lipolytica H222. De novo genome assembly shows three main chromosomal rearrangements compared to that of strain E150/CLIB122. This genomic resource will help integrate intraspecies diversity into synthetic biology projects that utilize Yarrowia as a biotechnological chassis for value-added chemical productions.

Exon junction complex components Y14 and Mago still play a role in budding yeast.

Since their divergence from Pezizomycotina, the mRNA metabolism of budding yeasts have undergone regressive evolution. With the dramatic loss of introns, a number of quality control mechanisms have been simplified or lost during evolution, such as the exon junction complex (EJC). We report the identification of the core EJC components, Mago, Y14, and eIF4A3, in at least seven Saccharomycotina species, including Yarrowia lipolytica. Peripheral factors that join EJC, either to mediate its assembly (Ibp160 or Cwc22), or trigger downstream processes, are present in the same species, forming an evolutionary package. Co-immunoprecipitation studies in Y. lipolytica showed that Mago and Y14 have retained the capacity to form heterodimers, which successively bind to the peripheral factors Upf3, Aly/REF, and Pym. Phenotypes and RNA-Seq analysis of EJC mutants showed evidence of Y14 and Mago involvement in mRNA metabolism. Differences in unspliced mRNA levels suggest that Y14 binding either interferes with pre-mRNA splicing or retains mRNA in the nucleus before their export and translation. These findings indicate that yeast could be a relevant model for understanding EJC function.

Investigation of Genetic Relationships Between Hanseniaspora Species Found in Grape Musts Revealed Interspecific Hybrids With Dynamic Genome Structures.

Hanseniaspora, a predominant yeast genus of grape musts, includes sister species recently reported as fast evolving. The aim of this study was to investigate the genetic relationships between the four most closely related species, at the population level. A multi-locus sequence typing strategy based on five markers was applied on 107 strains, confirming the clear delineation of species H. uvarum, H. opuntiae, H. guilliermondii, and H. pseudoguilliermondii. Huge variations were observed in the level of intraspecific nucleotide diversity, and differences in heterozygosity between species indicate different life styles. No clear population structure was detected based on geographical or substrate origins. Instead, H. guilliermondii strains clustered into two distinct groups, which may reflect a recent step toward speciation. Interspecific hybrids were detected between H. opuntiae and H. pseudoguilliermondii. Their characterization using flow cytometry, karyotypes and genome sequencing showed different genome structures in different ploidy contexts: allodiploids, allotriploids, and allotetraploids. Subculturing of an allotriploid strain revealed chromosome loss equivalent to one chromosome set, followed by an auto-diploidization event, whereas another auto-diploidized tetraploid showed a segmental duplication. Altogether, these results suggest that Hanseniaspora genomes are not only fast evolving but also highly dynamic.

Identification of telomerase RNAs in species of the Yarrowia clade provides insights into the co-evolution of telomerase, telomeric repeats and telomere-binding proteins.

Telomeric repeats in fungi of the subphylum Saccharomycotina exhibit great inter- and intra-species variability in length and sequence. Such variations challenged telomeric DNA-binding proteins that co-evolved to maintain their functions at telomeres. Here, we compare the extent of co-variations in telomeric repeats, encoded in the telomerase RNAs (TERs), and the repeat-binding proteins from 13 species belonging to the Yarrowia clade. We identified putative TER loci, analyzed their sequence and secondary structure conservation, and predicted functional elements. Moreover, in vivo complementation assays with mutant TERs showed the functional importance of four novel TER substructures. The TER-derived telomeric repeat unit of all species, except for one, is 10 bp long and can be represented as 5'-TTNNNNAGGG-3', with repeat sequence variations occuring primarily outside the vertebrate telomeric motif 5'-TTAGGG-3'. All species possess a homologue of the Yarrowia lipolytica Tay1 protein, YlTay1p. In vitro, YlTay1p displays comparable DNA-binding affinity to all repeat variants, suggesting a conserved role among these species. Taken together, these results add significant insights into the co-evolution of TERs, telomeric repeats and telomere-binding proteins in yeasts.

An agent-based model for predicting the prevalence of Trypanosoma cruzi I and II in their host and vector populations.

Trypanosoma cruzi is the causative agent of Chagas disease, an endemic human parasitosis in Latin America. This protozoan is transmitted to human and other mammals by blood-feeding bugs belonging to the Triatominae subfamily. There are two strains (T. cruzi I and T. cruzi II) presenting different biological and ecological characteristics. An original agent-based model (ABM) was designed for predicting the prevalence (i.e., proportion of infected individuals in the total population at a given time) of T. cruzi I and II during single and mixed infections. The ABM was calibrated from experimental data retrieved from literature. It was shown that inclusion of reservoir hosts as supplementary type of agent in the model was necessary for obtaining realistic simulation results of the prevalence of the two strains. This is totally in agreement with experimental and field observations on the importance of reservoirs in the parasite transmission cycle. Proposals were made for refining the model. More generally, the advantages and limitations of the ABM in parasitology modeling have been discussed.

Genome Sequence of Torulaspora microellipsoides CLIB 830T.

We report here the genome sequence of the ascomycetous yeast Torulaspora microellipsoides CLIB 830T A reference genome for this species, which has been found as a donor of genetic material in wine strains of Saccharomyces cerevisiae, will undoubtedly give clues to our understanding of horizontal transfer mechanisms between species in the wine environment.

EUF1 - a newly identified gene involved in erythritol utilization in Yarrowia lipolytica.

The gene YALI0F01562g was identified as an important factor involved in erythritol catabolism of the unconventional yeast Yarrowia lipolytica. Its putative role was identified for the first time by comparative analysis of four Y. lipolytica strains: A-101.1.31, Wratislavia K1, MK1 and AMM. The presence of a mutation that seriously damaged the gene corresponded to inability of the strain Wratislavia K1 to utilize erythritol. RT-PCR analysis of the strain MK1 demonstrated a significant increase in YALI0F01562g expression during growth on erythritol. Further studies involving deletion and overexpression of the selected gene showed that it is indeed essential for efficient erythritol assimilation. The deletion strain Y. lipolytica AMM∆euf1 was almost unable to grow on erythritol as the sole carbon source. When the strain was applied in the process of erythritol production from glycerol, the amount of erythritol remained constant after reaching the maximal concentration. Analysis of the YALI0F01562g gene sequence revealed the presence of domains characteristic for transcription factors. Therefore we suggest naming the studied gene Erythritol Utilization Factor - EUF1.

Enhancing Structural Annotation of Yeast Genomes with RNA-Seq Data.

The number of fully sequenced genomes of yeasts is dramatically increasing but both structural and functional annotation quality are usually neglected, as most frequently based on automatic annotation transfer tools from reference genomes. RNA sequencing technologies offer the possibility to better characterize yeast transcriptomes and to correct or improve the prediction of mRNA, ncRNA, or miscellaneous RNA. We describe a computational approach to enhance structural annotation of yeast genomes based on RNA-Seq data exploitation. The proposed pipeline is primarily based on read mapping with TopHat2. Mapping outputs are then used for various applications such as: (1) validation of exon-exon junctions of predicted transcripts, (2) definition of new transcribed features, (3) prediction of 3' UTR, and (4) identification of extra features absent from the genome assembly. We strongly encourage curators to proceed to a manual validation and editing of the reference genome. Releasing genomes with high-quality annotation is an important issue, as they will be considered as references for further predictions.

Whole-Genome Sequencing and Intraspecific Analysis of the Yeast Species Lachancea quebecensis.

The gold standard in yeast population genomics has been the model organism Saccharomyces cerevisiae. However, the exploration of yeast species outside the Saccharomyces genus is essential to broaden the understanding of genome evolution. Here, we report the analyses of whole-genome sequences of nineisolates from the recently described yeast species Lachancea quebecensis. The genome of one isolate was assembled and annotated, and the intraspecific variability within L. quebecensis was surveyed by comparing the sequences from the eight other isolates to this reference sequence. Our study revealed that these strains harbor genomes with an average nucleotide diversity of π = 2 × 10(-3) which is slightly lower, although on the same order of magnitude, as that previously determined for S. cerevisiae (π = 4 × 10(-3)). Our results show that even though these isolates were all obtained from a relatively isolated geographic location, the same ecological source, and represent a smaller sample size than is available for S. cerevisiae, the levels of divergence are similar to those observed in this model species. This divergence is essentially linked to the presence of two distinct clusters delineated according to geographic location. However, even with relatively similar ranges of genome divergence, L. quebecensis has an extremely low global phenotypic variance of 0.062 compared with 0.59 previously determined in S. cerevisiae.