Genome Sequence Analysis of the Oleaginous Yeast, Rhodotorula diobovata, and Comparison of the Carotenogenic and Oleaginous Pathway Genes and Gene Products with Other Oleaginous Yeasts

Extracellular oil production by Rhodotorula paludigena BS15 for biorefinery without complex downstream processes

To realize biomass refinery without complex downstream processes, we extensively screened for microbial strains that efficiently produce extracellular oil from sugars. Rhodotorula paludigena (formerly Rhodosporidium paludigenum) BS15 was found to efficiently produce polyol esters of fatty acids (PEFAs), which mainly comprised of 3-acetoxypalmitic acid and partially acetylated mannitol/arabinitol. To evaluate the performance of this strain, fed-batch fermentation was demonstrated on a flask scale, and 110 g/L PEFA and 103 g/L dry cells were produced in 12 days. To the best of our knowledge, the strain BS15 exhibited the highest PEFA titer (g/L) ever to be reported so far. Because the PEFA precipitated at the bottom of the culture broth, it could be easily recovered by simply discarding the upper phase. Various carbon sources can be utilized for cell growth and/or PEFA production, which signifies the potential for converting diverse biomass sources. Two different types of next-generation sequencers, Illumina HiSeq and Oxford Nanopore PromethION, were used to analyze the whole-genome sequence of the strain BS15. The integrative data analysis generated a high-quality and reliable reference genome for PEFA-producing R. paludigena. The 22.5-M base genome sequence and the estimated genes were registered in Genbank (accession numbers BQKY01000001-BQKY01000019). KEY POINTS: • R. paludigena BS15 was isolated after an extensive screening of extracellular oil producers from natural sources. • Fed-batch fermentation of R. paludigena BS15 yielded 110 g/L of PEFA, which is the highest titer ever reported to date. • Combined analysis using Illumina and Oxford Nanopore sequencers produced the near-complete genome sequence.

Diversity, astaxanthin production, and genomic analysis of Rhodotorula paludigena SP9-15

Astaxanthin is a carotenoid known for its powerful antioxidant properties. This study focused on isolating yeast strains capable of producing astaxanthin from flower and fruit samples collected in Thailand. Out of 115 isolates, 11 strains were identified that produced astaxanthin. Molecular identification techniques revealed that these isolates belonged to two species: Rhodotorula paludigena (5 isolates) and Rhodosporidiobolus ruineniae (6 isolates). Whole-genome analysis of one representative strain, R. paludigena SP9-15, identified putative candidate astaxanthin synthesis-associated genes, such as CrtE, CrtYB, CrtI, CrtS, CrtR, CrtW, CrtO, and CrtZ. High-performance liquid chromatography (HPLC) and liquid chromatography-mass spectrometry (LC-MS) confirmed astaxanthin production. Further optimization of astaxanthin production was carried out by investigating the effects of various factors on the growth rate and astaxanthin production. The optimal conditions were 40 g/L glucose as a carbon source, pH 7.5, and cultivation at 25 °C with 200 rpm for 3 days. Under these conditions, R. paludigena SP9-15 synthesized biomass of 11.771 ± 0.003 g/L, resulting in astaxanthin with a content of 0.558 ± 0.018 mg/g DCW (dry cell weight), an astaxanthin yield of 6.565 ± 0.238 mg/L, and astaxanthin productivity of 2.188 ± 0.069 g/L/day. These findings provide insights into astaxanthin production using red yeast strains from Thailand and highlight the potential of R. paludigena SP9-15 for further application.

Integrative analysis of genomic and metabolomic data reveals key metabolic pathways involved in lipid and carotenoid biosynthesis in oleaginous red yeast Rhodosporidiobolus odoratus XQR

Rhodosporidiobolus odoratus, one of the oleaginous red yeasts, is gaining biotechnological importance for their ability to produce microbial lipids and carotenoids. However, to date, the genomic resource underling lipid and carotenoid biosynthesis in R. odoratus has not been reported. Here, we reported the first genome assembly of R. odoratus using a combination of PacBio and Illumina techniques. The final genome assembly is 23.74 Mb in size, containing 52 scaffolds with a N50 length of 1200,460 bp and a GC content of 56.90%. Benchmarking Universal Single-Copy Orthologs (BUSCO) assessment showed that our assembly contains 94.23% of Basidiomycota universal single-copy orthologs. The genome was predicted to contain 4986 protein-coding genes, 4967 of which were functionally annotated. Metabolomic profiling identified 574 lipids, 3 carotenoids, and 208 volatile organic compounds synthesized by R. odoratus. Integrative analysis of genomics and metabolomics provides insights into the biosynthesis of lipid, carotenoid, and other bioactive compounds in R. odoratus. Collectively, the results presented herein greatly enhance our understanding of R. odoratus in lipids and carotenoids biosynthesis, and thus further accelerate its fundamental molecular investigations and biotechnological applications.

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.