De novo assembly of Dekkera bruxellensis: a multi technology approach using short and long-read sequencing and optical mapping

Relationships Among Origin, Genotype, and Oenological Traits of Brettanomyces Yeasts

Brettanomyces yeasts play a relevant role in the fermentation industry, showing controversial behavior. There is growing interest in these yeasts in the fermentation industry as beer and bioethanol production, while in winemaking, they are considered spoilage microorganisms mainly used to produce ethyl phenols. These compounds may alter wine's organoleptic characteristics, leading to significant economic loss. In this work, 45 Brettanomyces strains from seven different environments were genotyped and assayed for some oenological characters to investigate the possible relationship among sources of isolation, genotype characterization, and oenological characters. The results of biotyping showed four main clusters which were also distinguished according to the oenological characters. The oenological characters also distinguished the strains based on the isolation source, suggesting an overall relation between origin and biotypes. The negative correlation between fermentation rate and ethyl phenols production in the Brettanomyces population may indicate the adaptation to hostile environments differently from crop strains that showed the opposite behavior. The overall results contribute to clarifying some features of Brettanomyces yeasts, even if further investigations into the ability of these yeasts to colonize winemaking environments are needed.

Molecular approaches improving our understanding of Brettanomyces physiology

Brettanomyces species and particularly B. bruxellensis as the most studied representative, are strongly linked to industrial fermentation processes. This association is considered either positive or undesirable depending on the industry. While in some brewing applications and in kombucha production Brettanomyces yeasts contribute to the flavour and aroma profile of these beverages, in winemaking and bioethanol production Brettanomyces is considered a spoilage or contaminant microorganism. Nevertheless, understanding Brettanomyces biology and metabolism in detail will benefit all industries. This review discusses recent molecular biology tools including genomics, transcriptomics and genetic engineering techniques that can improve our understanding of Brettanomyces physiology and how these approaches can be used to make the industrial potential of this species a reality.

Brettanomyces bruxellensis: Overview of the genetic and phenotypic diversity of an anthropized yeast

Human-associated microorganisms are ideal models to study the impact of environmental changes on species evolution and adaptation because of their small genome, short generation time, and their colonization of contrasting and ever-changing ecological niches. The yeast Brettanomyces bruxellensis is a good example of organism facing anthropogenic-driven selective pressures. It is associated with fermentation processes in which it can be considered either as a spoiler (e.g. winemaking, bioethanol production) or as a beneficial microorganism (e.g. production of specific beers, kombucha). Besides its industrial interests, noteworthy parallels and dichotomies with Saccharomyces cerevisiae propelled B. bruxellensis as a valuable complementary yeast model. In this review, we emphasize that the broad genetic and phenotypic diversity of this species is only beginning to be uncovered. Population genomic studies have revealed the co-existence of auto- and allotriploidization events with different evolutionary outcomes. The different diploid, autotriploid and allotriploid subpopulations are associated with specific fermented processes, suggesting independent adaptation events to anthropized environments. Phenotypically, B. bruxellensis is renowned for its ability to metabolize a wide variety of carbon and nitrogen sources, which may explain its ability to colonize already fermented environments showing low-nutrient contents. Several traits of interest could be related to adaptation to human activities (e.g. nitrate metabolization in bioethanol production, resistance to sulphite treatments in winemaking). However, phenotypic traits are insufficiently studied in view of the great genomic diversity of the species. Future work will have to take into account strains of varied substrates, geographical origins as well as displaying different ploidy levels to improve our understanding of an anthropized yeast's phenotypic landscape.

Near Chromosome-Level Genome Assembly and Annotation of Rhodotorula babjevae Strains Reveals High Intraspecific Divergence

The genus Rhodotorula includes basidiomycetous oleaginous yeast species. Rhodotorula babjevae can produce compounds of biotechnological interest such as lipids, carotenoids, and biosurfactants from low value substrates such as lignocellulose hydrolysate. High-quality genome assemblies are needed to develop genetic tools and to understand fungal evolution and genetics. Here, we combined short- and long-read sequencing to resolve the genomes of two R. babjevae strains, CBS 7808 (type strain) and DBVPG 8058, at chromosomal level. Both genomes are 21 Mbp in size and have a GC content of 68.2%. Allele frequency analysis indicates that both strains are tetraploid. The genomes consist of a maximum of 21 chromosomes with a size of 0.4 to 2.4 Mbp. In both assemblies, the mitochondrial genome was recovered in a single contig, that shared 97% pairwise identity. Pairwise identity between most chromosomes ranges from 82 to 87%. We also found indications for strain-specific extrachromosomal endogenous DNA. A total of 7591 and 7481 protein-coding genes were annotated in CBS 7808 and DBVPG 8058, respectively. CBS 7808 accumulated a higher number of tandem duplications than DBVPG 8058. We identified large translocation events between putative chromosomes. Genome divergence values between the two strains indicate that they may belong to different species.

Chromosome-level genome assembly and transcriptome-based annotation of the oleaginous yeast Rhodotorula toruloides CBS 14

Rhodotorula toruloides is an oleaginous yeast with high biotechnological potential. In order to understand the molecular physiology of lipid synthesis in R. toruloides and to advance metabolic engineering, a high-resolution genome is required. We constructed a genome draft of R. toruloides CBS 14, using a hybrid assembly approach, consisting of short and long reads generated by Illumina and Nanopore sequencing, respectively. The genome draft consists of 23 contigs and 3 scaffolds, with a N50 length of 1,529,952 bp, thus largely representing chromosomal organization. The total size of the genome is 20,534,857 bp and the overall GC content is 61.83%. Transcriptomic data from different growth conditions was used to aid species-specific gene annotation. We annotated 9464 genes and identified 11,691 transcripts. Furthermore, we demonstrated the presence of a potential plasmid, an extrachromosomal circular structure of about 11 kb with a copy number about three times as high as the other chromosomes.

Draft genome assemblies using sequencing reads from Oxford Nanopore Technology and Illumina platforms for four species of North American Fundulus killifish

BACKGROUND

Whole-genome sequencing data from wild-caught individuals of closely related North American killifish species (Fundulus xenicus, Fundulus catenatus, Fundulus nottii, and Fundulus olivaceus) were obtained using long-read Oxford Nanopore Technology (ONT) PromethION and short-read Illumina platforms.

FINDINGS

Draft de novo reference genome assemblies were generated using a combination of long and short sequencing reads. For each species, the PromethION platform was used to generate 30-45× sequence coverage, and the Illumina platform was used to generate 50-160× sequence coverage. Illumina-only assemblies were fragmented with high numbers of contigs, while ONT-only assemblies were error prone with low BUSCO scores. The highest N50 values, ranging from 0.4 to 2.7 Mb, were from assemblies generated using a combination of short- and long-read data. BUSCO scores were consistently >90% complete using the Eukaryota database.

CONCLUSIONS

High-quality genomes can be obtained from a combination of using short-read Illumina data to polish assemblies generated with long-read ONT data. Draft assemblies and raw sequencing data are available for public use. We encourage use and reuse of these data for assembly benchmarking and other analyses.

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.

Assessing Population Diversity of Brettanomyces Yeast Species and Identification of Strains for Brewing Applications

Brettanomyces yeasts have gained popularity in many sectors of the biotechnological industry, specifically in the field of beer production, but also in wine and ethanol production. Their unique properties enable Brettanomyces to outcompete conventional brewer’s yeast in industrially relevant traits such as production of ethanol and pleasant flavors. Recent advances in next-generation sequencing (NGS) and high-throughput screening techniques have facilitated large population studies allowing the selection of appropriate yeast strains with improved traits. In order to get a better understanding of Brettanomyces species and its potential for beer production, we sequenced the whole genome of 84 strains, which we make available to the scientific community and carried out several in vitro assays for brewing-relevant properties. The collection includes isolates from different substrates and geographical origin. Additionally, we have included two of the oldest Carlsberg Research Laboratory isolates. In this study, we reveal the phylogenetic pattern of Brettanomyces species by comparing the predicted proteomes of each strain. Furthermore, we show that the Brettanomyces collection is well described using similarity in genomic organization, and that there is a direct correlation between genomic background and phenotypic characteristics. Particularly, genomic patterns affecting flavor production, maltose assimilation, beta-glucosidase activity, and phenolic off-flavor (POF) production are reported. This knowledge yields new insights into Brettanomyces population survival strategies, artificial selection pressure, and loss of carbon assimilation traits. On a species-specific level, we have identified for the first time a POF negative Brettanomyces anomalus strain, without the main spoilage character of Brettanomyces species. This strain (CRL-90) has lost DaPAD1, making it incapable of converting ferulic acid to 4-ethylguaiacol (4-EG) and 4-ethylphenol (4-EP). This loss of function makes CRL-90 a good candidate for the production of characteristic Brettanomyces flavors in beverages, without the contaminant increase in POF. Overall, this study displays the potential of exploring Brettanomyces yeast species biodiversity to find strains with relevant properties applicable to the brewing industry.

Brettanomyces bruxellensis wine isolates show high geographical dispersal and long persistence in cellars

Brettanomyces bruxellensis is the main wine spoiler yeast all over the world, yet the structure of the populations associated with winemaking remains elusive. In this work, we considered 1411 wine isolates from 21 countries that were genotyped using twelve microsatellite markers. We confirmed that B. bruxellensis isolates from wine environments show high genetic diversity, with 58 and 42% of putative triploid and diploid individuals respectively distributed in 5 main genetic groups. The distribution in the genetic groups varied greatly depending on the country and/or the wine-producing region. However, the two possible triploid wine groups showing sulfite resistance/tolerance were identified in almost all regions/countries. Genetically identical isolates were also identified. The analysis of these clone groups revealed that a given genotype could be isolated repeatedly in the same winery over decades, demonstrating unsuspected persistence ability. Besides cellar residency, a great geographic dispersal was also evidenced, with some genotypes isolated in wines from different continents. Finally, the study of old isolates and/or isolates from old vintages revealed that only the diploid groups were identified prior 1990 vintages. The putative triploid groups were identified in subsequent vintages, and their proportion has increased steadily these last decades, suggesting adaptation to winemaking practices such as sulfite use. A possible evolutionary scenario explaining these results is discussed.

Re-sequencing and optical mapping reveals misassemblies and real inversions on Corynebacterium pseudotuberculosis genomes

The number of draft genomes deposited in Genbank from the National Center for Biotechnology Information (NCBI) is higher than the complete ones. Draft genomes are assemblies that contain fragments of misassembled regions (gaps). Such draft genomes present a hindrance to the complete understanding of the biology and evolution of the organism since they lack genomic information. To overcome this problem, strategies to improve the assembly process are developed continuously. Also, the greatest challenge to the assembly progress is the presence of repetitive DNA regions. This article highlights the use of optical mapping, to detect and correct assembly errors in Corynebacterium pseudotuberculosis. We also demonstrate that choosing a reference genome should be done with caution to avoid assembly errors and loss of genetic information.

Assembly and Analysis of the Genome Sequence of the Yeast Brettanomyces naardenensis CBS 7540

Brettanomyces naardenensis is a spoilage yeast with potential for biotechnological applications for production of innovative beverages with low alcohol content and high attenuation degree. Here, we present the first annotated genome of B. naardenensis CBS 7540. The genome of B. naardenensis CBS 7540 was assembled into 76 contigs, totaling 11,283,072 nucleotides. In total, 5168 protein-coding sequences were annotated. The study provides functional genome annotation, phylogenetic analysis, and discusses genetic determinants behind notable stress tolerance and biotechnological potential of B. naardenensis.

Development of a genetic transformation toolkit for Brettanomyces bruxellensis

Brettanomyces bruxellensis is usually considered a spoilage microorganism, responsible for significant economic losses during the production of fermented beverages such as wine, beer and cider, though for some styles of beer its influence is essential. In recent years, the competitiveness of this yeast in bioethanol production processes has brought to attention its broader biotechnological potential. Furthermore, the species has evolved key fermentation traits in parallel with Saccharomyces cerevisiae. Attempts to better understand B. bruxellensis physiology through genomics-driven research have been hampered by a lack of functional genomics tools. Genetic transformation for B. bruxellensis has only been developed recently and with limited efficiency. Here we describe gene transformation cassettes tailored for B. bruxellensis, which provide multiple drug-resistant markers and the ability to tag B. bruxellensis with different fluorescent proteins. All marker cassettes resulted in increased transformation efficiency compared to the maximum reported in literature, with one cassette, TDH1p natMX, showing five times greater efficiency. Transformation cassettes encoding fluorescent proteins enabled discrimination between subpopulations of transformed B. bruxellensis cells by flow cytometry and fluorescent microscopy. Thus, the genetic transformation toolkit described here unlocks several molecular applications such as strain tagging, insertional mutagenesis and potentially targeted gene deletion.

Chromosomal genome assembly of the ethanol production strain CBS 11270 indicates a highly dynamic genome structure in the yeast species Brettanomyces bruxellensis

Here, we present the genome of the industrial ethanol production strain Brettanomyces bruxellensis CBS 11270. The nuclear genome was found to be diploid, containing four chromosomes with sizes of ranging from 2.2 to 4.0 Mbp. A 75 Kbp mitochondrial genome was also identified. Comparing the homologous chromosomes, we detected that 0.32% of nucleotides were polymorphic, i.e. formed single nucleotide polymorphisms (SNPs), 40.6% of them were found in coding regions (i.e. 0.13% of all nucleotides formed SNPs and were in coding regions). In addition, 8,538 indels were found. The total number of protein coding genes was 4897, of them, 4,284 were annotated on chromosomes; and the mitochondrial genome contained 18 protein coding genes. Additionally, 595 genes, which were annotated, were on contigs not associated with chromosomes. A number of genes was duplicated, most of them as tandem repeats, including a six-gene cluster located on chromosome 3. There were also examples of interchromosomal gene duplications, including a duplication of a six-gene cluster, which was found on both chromosomes 1 and 4. Gene copy number analysis suggested loss of heterozygosity for 372 genes. This may reflect adaptation to relatively harsh but constant conditions of continuous fermentation. Analysis of gene topology showed that most of these losses occurred in clusters of more than one gene, the largest cluster comprising 33 genes. Comparative analysis against the wine isolate CBS 2499 revealed 88,534 SNPs and 8,133 indels. Moreover, when the scaffolds of the CBS 2499 genome assembly were aligned against the chromosomes of CBS 11270, many of them aligned completely, some have chunks aligned to different chromosomes, and some were in fact rearranged. Our findings indicate a highly dynamic genome within the species B. bruxellensis and a tendency towards reduction of gene number in long-term continuous cultivation.

Study of the whole genome, methylome and transcriptome of Cordyceps militaris

The complete genome of Cordyceps militaris was sequenced using single-molecule real-time (SMRT) sequencing technology at a coverage over 300×. The genome size was 32.57 Mb, and 14 contigs ranging from 0.35 to 4.58 Mb with an N50 of 2.86 Mb were assembled, including 4 contigs with telomeric sequences on both ends and an additional 8 contigs with telomeric sequences on either the 5′ or 3′ end. A methylome database of the genome was constructed using SMRT and m4C and m6A methylated nucleotides, and many unknown modification types were identified. The major m6A methylation motif is GA and GGAG, and the major m4C methylation motif is GC or CG/GC. In the C. militaris genome DNA, there were four types of methylated nucleotides that we confirmed using high-resolution LCMS-IT-TOF. Using PacBio Iso-Seq, a total of 31,133 complete cDNA sequences were obtained in the fruiting body. The conserved domains of the nontranscribed regions of the genome include TATA boxes, which are the initial regions of genome replication. There were 406 structural variants between the HN and CM01 strains, and there were 1,114 structural variants between the HN and ATCC strains.

Novel genome and genome-wide SNPs reveal early fragmentation effects in an edge-tolerant songbird population across an urbanized tropical metropolis

Although edge-tolerant species are known to benefit from habitat fragmentation, less is known about the population genetic impacts fragmentation may exert on edge-tolerant species. We examined the landscape genomic structure of an edge-tolerant forest-dependent bird species, the Striped Tit-Babbler Mixornis gularis, in the heavily urbanized island of Singapore to determine if two centuries of fragmentation have led to signs of isolation and loss of population-genetic diversity in different parts of the island. We obtained a high-quality complete reference genome with 78x coverage. Using almost 4000 SNPs from double-digest RAD-Sequencing across 46 individuals, we found that the population has likely experienced a recent contraction in effective population size and presently exhibits low population genetic diversity. Using empirical and simulation-based landscape genomic analyses, we also found that the subtle population genetic structure observed in the Striped Tit-Babbler population in Singapore is likely driven by isolation by distance resulting from limited dispersal. Our results demonstrate that population genetic impoverishment and subdivision can accumulate at relatively rapid rates in edge-tolerant bird species such as the Striped Tit-Babbler as a result of fragmentation, and that subtle spatial genetic structure can be detected over fine spatial and temporal scales using relatively few multilocus genomic SNPs.

High-Quality de Novo Genome Assembly of the Dekkera bruxellensis Yeast Using Nanopore MinION Sequencing

Genetic variation in natural populations represents the raw material for phenotypic diversity. Species-wide characterization of genetic variants is crucial to have a deeper insight into the genotype-phenotype relationship. With the advent of new sequencing strategies and more recently the release of long-read sequencing platforms, it is now possible to explore the genetic diversity of any nonmodel organisms, representing a fundamental resource for biological research. In the frame of population genomic surveys, a first step is to obtain the complete sequence and high-quality assembly of a reference genome. Here, we sequenced and assembled a reference genome of the nonconventional Dekkera bruxellensis yeast. While this species is a major cause of wine spoilage, it paradoxically contributes to the specific flavor profile of some Belgium beers. In addition, an extreme karyotype variability is observed across natural isolates, highlighting that D. bruxellensis genome is very dynamic. The whole genome of the D. bruxellensis UMY321 isolate was sequenced using a combination of Nanopore long-read and Illumina short-read sequencing data. We generated the most complete and contiguous de novo assembly of D. bruxellensis to date and obtained a first glimpse into the genomic variability within this species by comparing the sequences of several isolates. This genome sequence is therefore of high value for population genomic surveys and represents a reference to study genome dynamic in this yeast species.

Genomics Perspectives on Metabolism, Survival Strategies, and Biotechnological Applications of Brettanomyces bruxellensis LAMAP2480

Wine production is an important commercial issue for the liquor industry. The global production was estimated at 275.7 million hectoliters in 2015. The loss of wine production due to Brettanomyces bruxellensis contamination is currently a problem. This yeast causes a "horse sweat" flavor in wine, which is an undesired organoleptic attribute. To date, 6 B. bruxellensis annotated genome sequences are available (LAMAP2480, AWRI1499, AWRI1608, AWRI1613, ST05.12/22, and CBS2499), and whole genome comparisons between strains are limited. In this article, we reassembled and reannotated the genome of B. bruxellensis LAMAP2480, obtaining a 27-Mb assembly with 5.5 kb of N50. In addition, the genome of B. bruxellensis LAMAP2480 was analyzed in the context of spoilage yeast and potential as a biotechnological tool. In addition, we carried out an exploratory transcriptomic analysis of this strain grown in synthetic wine. Several genes related to stress tolerance, micronutrient acquisition, ethanol production, and lignocellulose assimilation were found. In conclusion, the analysis of the genome of B. bruxellensis LAMAP2480 reaffirms the biotechnological potential of this strain. This research represents an interesting platform for the study of the spoilage yeast B. bruxellensis.