Codon usage and codon context bias in Xanthophyllomyces dendrorhous

Comparative analysis of codon usage patterns in the chloroplast genomes of nine forage legumes

Leguminosae is one of the three largest families of angiosperms after Compositae and Orchidaceae. It is widely distributed and grows in a variety of environments, including plains, mountains, deserts, forests, grasslands, and even waters where almost all legumes can be found. It is one of the most important sources of starch, protein and oil in the food of mankind and also an important source of high-quality forage material for animals, which has important economic significance. In our study, the codon usage patterns and variation sources of the chloroplast genome of nine important forage legumes were systematically analyzed. Meanwhile, we also constructed a phylogenetic tree based on the whole chloroplast genomes and protein coding sequences of these nine forage legumes. Our results showed that the chloroplast genomes of nine forage legumes end with A/T bases, and seven identical high-frequency (HF) codons were detected among the nine forage legumes. ENC-GC3s mapping, PR2 analysis, and neutral analysis showed that the codon bias of nine forage legumes was influenced by many factors, among which natural selection was the main influencing factor. The codon usage frequency showed that the Nicotiana tabacum and Saccharomyces cerevisiae can be considered as receptors for the exogenous expression of chloroplast genes of these nine forage legumes. The phylogenetic relationships of the chloroplast genomes and protein coding genes were highly similar, and the nine forage legumes were divided into three major clades. Among the clades Melilotus officinalis was more closely related to Medicago sativa, and Galega officinalis was more closely related to Galega orientalis. This study provides a scientific basis for the molecular markers research, species identification and phylogenetic studies of forage legumes.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-024-01421-0.

Molecular Dissection of Herpes Simplex Virus Type 1 to Elucidate Molecular Mechanisms Behind Latency and Comparison of Its Codon Usage Patterns with Genes Modulated During Alzheimer's Disease as a Part of Host-Pathogen Interaction

BACKGROUND

Herpes simplex virus type 1 (HSV-1) is associated with Alzheimer's disease, which goes into a cycle of latency and reactivation. The present study was envisaged to understand the reasons for latency and specific molecular patterns present in the HSV-1.

OBJECTIVE

The objective is the molecular dissection of Herpes simplex virus type 1 to elucidate molecular mechanisms behind latency and compare its codon usage patterns with genes modulated during Alzheimer's disease as a part of host-pathogen interaction.

METHODS

In the present study, we tried to investigate the potential reasons for the latency of HSV-1 virus bioinformatically by determining the CpG patterns. Also, we investigated the codon usage pattern, the presence of rare codons, codon context, and protein properties.

RESULTS

The top 222 codon pairs graded based on their frequency in the HSV-1 genome revealed that with only one exception (CUG-UUU), all other codon pairs have codons ending with G/C. Considering it an extension of host-pathogen interaction, we compared HSV-1 codon usage with that of codon usage of genes modulated during Alzheimer's disease, and we found that CGT and TTT are only two codons that exhibited similar codon usage patterns and other codons showed statistically highly significant different codon preferences. Dinucleotide CpG tends to mutate to TpG, suggesting the presence of mutational forces and the imperative role of CpG methylation in HSV-1 latency.

CONCLUSIONS

Upon comparison of codon usage between HSV-1 and Alzheimer's disease genes, no similarities in codon usage were found as a part of host-pathogen interaction. CpG methylation plays an imperative role in latency HSV-1.

Assembly and comparative genome analysis of a Patagonian Aureobasidium pullulans isolate reveals unexpected intraspecific variation

Aureobasidium pullulans is a yeast-like fungus with remarkable phenotypic plasticity widely studied for its importance for the pharmaceutical and food industries. So far, genomic studies with strains from all over the world suggest they constitute a genetically unstructured population, with no association by habitat. However, the mechanisms by which this genome supports so many phenotypic permutations are still poorly understood. Recent works have shown the importance of sequencing yeast genomes from extreme environments to increase the repertoire of phenotypic diversity of unconventional yeasts. In this study, we present the genomic draft of A. pullulans strain from a Patagonian yeast diversity hotspot, re-evaluate its taxonomic classification based on taxogenomic approaches, and annotate its genome with high-depth transcriptomic data. Our analysis suggests this isolate could be considered a novel variant at an early stage of the speciation process. The discovery of divergent strains in a genomically homogeneous group, such as A. pullulans, can be valuable in understanding the evolution of the species. The identification and characterization of new variants will not only allow finding unique traits of biotechnological importance, but also optimize the choice of strains whose phenotypes will be characterized, providing new elements to explore questions about plasticity and adaptation.

Unraveling the Molecular Basis of Mycosporine Biosynthesis in Fungi

The Phaffia rhodozyma UCD 67-385 genome harbors a 7873 bp cluster containing DDGS, OMT, and ATPG, encoding 2-desmethy-4-deoxygadusol synthase, O-methyl transferase, and ATP-grasp ligase, respectively, of the mycosporine glutaminol (MG) biosynthesis pathway. Homozygous deletion mutants of the entire cluster, single-gene mutants, and the Δddgs^-/-;Δomt^-/- and Δomt^-/-;Δatpg^-/- double-gene mutants did not produce mycosporines. However, Δatpg^-/- accumulated the intermediate 4-deoxygadusol. Heterologous expression of the DDGS and OMT or DDGS, OMT, and ATPG cDNAs in Saccharomyces cerevisiae led to 4-deoxygadusol or MG production, respectively. Genetic integration of the complete cluster into the genome of the non-mycosporine-producing CBS 6938 wild-type strain resulted in a transgenic strain (CBS 6938_MYC) that produced MG and mycosporine glutaminol glucoside. These results indicate the function of DDGS, OMT, and ATPG in the mycosporine biosynthesis pathway. The transcription factor gene mutants Δmig1^-/-, Δcyc8^-/-, and Δopi1^-/- showed upregulation, Δrox1^-/- and Δskn7^-/- showed downregulation, and Δtup6^-/- and Δyap6^-/- showed no effect on mycosporinogenesis in glucose-containing medium. Finally, comparative analysis of the cluster sequences in several P. rhodozyma strains and the four newly described species of the genus showed the phylogenetic relationship of the P. rhodozyma strains and their differentiation from the other species of the genus Phaffia.

Insights into Synonymous Codon Usage Bias in Hepatitis C Virus and Its Adaptation to Hosts

Hepatitis C virus (HCV) is enveloped RNA virus, encoding for a polyprotein that is processed by cellular proteases. The virus is responsible for liver cirrhosis, allograft rejection, and human hepatocellular carcinoma. Based on studies including compositional analysis, odds ratio analysis, parity analysis, skew analysis, relative synonymous codon usage, codon bias, and protein properties, it was evident that codon usage bias in HCV is dependent upon the nucleotide composition. Codon context analysis revealed CTC-CTG as a preferred codon pair. While CGA and CGT codons were rare, none of the codons were rare in HCV-like viruses envisaged in the present study. Many of the preferred codon pairs were valine amino acid-initiated, which possibly infers viral infectivity; hence the role of selection forces appears to act on the HCV genome, which was further validated by neutrality analysis where selection accounted for 87.28%, while mutation accounted for 12.72% force shaping codon usage. Furthermore, codon usage was correlated with the length of the genome. HCV viruses prefer valine-initiated codon pairs, while HCV-like viruses prefer alanine-initiated codon pairs. The HCV host range is very narrow and is confined to only humans and chimpanzees. Based on indices including codon usage correlation analysis, similarity index, and relative codon deoptimization index, it is evident in the study that the chimpanzee is the primary host of the virus. The present study helped elucidate the preferred host for HCV. The information presented in the study paved the way for generating an attenuated vaccine candidate through viral recoding, with finely tuned nucleotide composition and a perfect balance of preferred and rare codons.

Comparative analysis of codon usage patterns in chloroplast genomes of ten Epimedium species

BACKGROUND

The Phenomenon of codon usage bias exists in the genomes of prokaryotes and eukaryotes. The codon usage pattern is affected by environmental factors, base mutation, gene flow and gene expression level, among which natural selection and mutation pressure are the main factors. The study of codon preference is an effective method to analyze the source of evolutionary driving forces in organisms. Epimedium species are perennial herbs with ornamental and medicinal value distributed worldwide. The chloroplast genome is self-replicating and maternally inherited which is usually used to study species evolution, gene expression and genetic transformation.

RESULTS

The results suggested that chloroplast genomes of Epimedium species preferred to use codons ending with A/U. 17 common high-frequency codons and 2-6 optimal codons were found in the chloroplast genomes of Epimedium species, respectively. According to the ENc-plot, PR2-plot and neutrality-plot, the formation of codon preference in Epimedium was affected by multiple factors, and natural selection was the dominant factor. By comparing the codon usage frequency with 4 common model organisms, it was found that Arabidopsis thaliana, Populus trichocarpa, and Saccharomyces cerevisiae were suitable exogenous expression receptors.

CONCLUSION

The evolutionary driving force in the chloroplast genomes of 10 Epimedium species probably comes from mutation pressure. Our results provide an important theoretical basis for evolutionary analysis and transgenic research of chloroplast genes.

Role of ROX1, SKN7, and YAP6 Stress Transcription Factors in the Production of Secondary Metabolites in Xanthophyllomyces dendrorhous

Xanthophyllomyces dendrorhous is a natural source of astaxanthin and mycosporines. This yeast has been isolated from high and cold mountainous regions around the world, and the production of these secondary metabolites may be a survival strategy against the stress conditions present in its environment. Biosynthesis of astaxanthin is regulated by catabolic repression through the interaction between MIG1 and corepressor CYC8–TUP1. To evaluate the role of the stress-associated transcription factors SKN7, ROX1, and YAP6, we employed an omic and phenotypic approach. Null mutants were constructed and grown in two fermentable carbon sources. The yeast proteome and transcriptome were quantified by iTRAQ and RNA-seq, respectively. The total carotenoid, sterol, and mycosporine contents were determined and compared to the wild-type strain. Each mutant strain showed significant metabolic changes compared to the wild type that were correlated to its phenotype. In a metabolic context, the principal pathways affected were glycolysis/gluconeogenesis, the pentose phosphate (PP) pathway, and the citrate (TCA) cycle. Additionally, fatty acid synthesis was affected. The absence of ROX1 generated a significant decline in carotenoid production. In contrast, a rise in mycosporine and sterol synthesis was shown in the absence of the transcription factors SKN7 and YAP6, respectively.

Response to Cold: A Comparative Transcriptomic Analysis in Eight Cold-Adapted Yeasts

Microorganisms have evolved to colonize all biospheres, including extremely cold environments, facing several stressor conditions, mainly low/freezing temperatures. In general, terms, the strategies developed by cold-adapted microorganisms include the synthesis of cryoprotectant and stress-protectant molecules, cold-active proteins, especially enzymes, and membrane fluidity regulation. The strategy could differ among microorganisms and concerns the characteristics of the cold environment of the microorganism, such as seasonal temperature changes. Microorganisms can develop strategies to grow efficiently at low temperatures or tolerate them and grow under favorable conditions. These differences can be found among the same kind of microorganisms and from the same cold habitat. In this work, eight cold-adapted yeasts isolated from King George Island, subAntarctic region, which differ in their growth properties, were studied about their response to low temperatures at the transcriptomic level. Sixteen ORFeomes were assembled and used for gene prediction and functional annotation, determination of gene expression changes, protein flexibilities of translated genes, and codon usage bias. Putative genes related to the response to all main kinds of stress were found. The total number of differentially expressed genes was related to the temperature variation that each yeast faced. The findings from multiple comparative analyses among yeasts based on gene expression changes and protein flexibility by cellular functions and codon usage bias raise significant differences in response to cold among the studied Antarctic yeasts. The way a yeast responds to temperature change appears to be more related to its optimal temperature for growth (OTG) than growth velocity. Yeasts with higher OTG prepare to downregulate their metabolism to enter the dormancy stage. In comparison, yeasts with lower OTG perform minor adjustments to make their metabolism adequate and maintain their growth at lower temperatures.

Damage response protein 1 (Dap1) functions in the synthesis of carotenoids and sterols in Xanthophyllomyces dendrorhous

Cytochrome P450s (P450s) are heme-containing proteins involved in several cellular functions, including biosynthesis of steroidal hormones, detoxification of xenobiotic compounds, among others. Damage response protein 1 (Dap1) has been described as a positive regulator of P450s through protein-protein interactions in organisms such as Schizosaccharomyces pombe. Three P450s in the carotenogenic yeast Xanthophyllomyces dendrorhous have thus far been characterized: Cyp51 and Cyp61, which are involved in ergosterol biosynthesis, and CrtS (astaxanthin synthase), which is involved in biosynthesis of the carotenoid astaxanthin. In this work, we describe the X. dendrorhous DAP1 gene, deletion of which affected yeast pigmentation by decreasing the astaxanthin fraction and increasing the β-carotene (a substrate of CrtS) fraction, which is consistent with the known role of CrtS. We found that the proportion of ergosterol was also decreased in the Δdap1 mutant. However, even though the fractions of the end products of these two pathways (the synthesis of carotenoids and sterols) were decreased in the Δdap1 mutant, the transcript levels of genes from the P450 systems involved were higher than those in the wild-type strain. We demonstrate that Dap1 coimmunoprecipitates with these three P450s, suggesting that Dap1 interacts with these three proteins. We propose that Dap1 regulates the synthesis of astaxanthin and ergosterol in X. dendrorhous, probably by regulating the P450s involved in both biosynthetic pathways at the protein level. This work suggests a new role for Dap1 in the regulation of carotenoid biosynthesis in X. dendrorhous.

Codon pair optimization (CPO): a software tool for synthetic gene design based on codon pair bias to improve the expression of recombinant proteins in Pichia pastoris

Background

Codon optimization is a common method to improve protein expression levels in Pichia pastoris and the current strategy is to replace rare codons with preferred codons to match the codon usage bias. However, codon-pair contexts have a profound effect on translation efficiency by influencing both translational elongation rates and accuracy. Until now, it remains untested whether optimized genes based on codon pair bias results in higher protein expression levels compared to codon usage bias.

Results

In this study, an algorithm based on dynamic programming was introduced to develop codon pair optimization (CPO) which is a software tool to provide simple and efficient codon pair optimization for synthetic gene design in Pichia pastoris. Two reporters (MT1-MMP E2C6 and ADAM17 A9B8 scFvs) were employed to test the effects of codon pair bias and CPO optimization on their protein expression levels. Four variants of MT1-MMP E2C6 and ADAM17 A9B8 for each were generated, one variant with the best codon-pair context, one with the worst codon-pair context, one with unbiased codon-pair context, and another optimized based on codon usage. The expression levels of variants with the worst codon-pair context were almost undetectable by Western blot and the variants with the best codon-pair context were expressed well. The expression levels on MT1-MMP E2C6 and ADAM17 A9B8 were more than five times and seven times higher in the optimized sequences based on codon-pair context compared to that based on codon usage, respectively. The results indicated that the codon-pair context-based codon optimization is more effective in enhancing expression of protein in Pichia pastoris.

Conclusions

Codon-pair context plays an important role on the protein expression in Pichia pastoris. The codon pair optimization (CPO) software developed in this study efficiently improved the protein expression levels of exogenous genes in Pichia pastoris, suggesting gene design based on codon pair bias is an alternative strategy for high expression of recombinant proteins in Pichia pastoris.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12934-021-01696-y.

Analysis of codon usage patterns of the chloroplast genome in Delphinium grandiflorum L. reveals a preference for AT-ending codons as a result of major selection constraints

Background

Codon usage bias analysis is a suitable strategy for identifying the principal evolutionary driving forces in different organisms. Delphinium grandiflorum L. is a perennial herb with high economic value and typical biological characteristics. Evolutionary analysis of D. grandiflorum can provide a rich resource of genetic information for developing hybridization resources of the genus Delphinium.

Methods

Synonymous codon usage (SCU) and related indices of 51 coding sequences from the D. grandiflorum chloroplast (cp) genome were calculated using Codon W, Cups of EMBOSS, SPSS and Microsoft Excel. Multivariate statistical analysis combined by principal component analysis (PCA), correspondence analysis (COA), PR2-plot mapping analysis and ENC plot analysis was then conducted to explore the factors affecting the usage of synonymous codons.

Results

The SCU bias of D. grandiflorum was weak and codons preferred A/T ending. A SCU imbalance between A/T and G/C at the third base position was revealed by PR2-plot mapping analysis. A total of eight codons were identified as the optimal codons. The PCA and COA results indicated that base composition (GC content, GC₃ content) and gene expression were important for SCU bias. A majority of genes were distributed below the expected curve from the ENC plot analysis and up the standard curve by neutrality plot analysis. Our results showed that with the exception of notable mutation pressure effects, the majority of genetic evolution in the D. grandiflorum cp genome might be driven by natural selection.

Discussions

Our results provide a theoretical foundation for elucidating the genetic architecture and mechanisms of D. grandiflorum, and contribute to enriching D. grandiflorum genetic resources.

Phenotypic Analysis of Mutants of Ergosterol Biosynthesis Genes (ERG3 and ERG4) in the Red Yeast Xanthophyllomyces dendrorhous

Xanthophyllomyces dendrorhous synthesizes astaxanthin, a carotenoid used in aquaculture. Astaxanthin is synthesized from metabolites of the mevalonate pathway, which are also precursors for sterols biosynthesis. The interruption of the CYP61 gene, which is involved in the synthesis of ergosterol (mutant CBS.cyp61 -), resulted in a phenotype that overproduces carotenoids due to the activation of the SREBP pathway. In this work, we constructed other mutants of ergosterol biosynthesis in this yeast to evaluate whether they have the same phenotype as mutant CBS.cyp61 -. By bioinformatic analysis, the ERG3 and ERG4 genes of X. dendrorhous were identified, and each gene was deleted in the wild-type strain. Mutants CBS.Δerg3 and CBS.Δerg4 did not produce ergosterol; CBS.Δerg3 primarily accumulated episterol, and CBS.Δerg4 primarily accumulated ergosta-5,7,22,24(28)-tetraenol. The transcription levels of the HMGS gene of the mevalonate pathway were evaluated by RT-qPCR, which showed a slight increase in CBS.Δerg4, but the transcription levels were still 10-fold lower than in strain CBS.cyp61 -. Both CBS.Δerg3 and CBS.Δerg4 did not overproduce carotenoids, even though they do not produce ergosterol. Thus, the results of this study indicate that the absence of ergosterol does not activate the SREBP pathway in X. dendrorhous, but rather it depends on other alterations in sterol composition.

Convergence between Regulation of Carbon Utilization and Catabolic Repression in Xanthophyllomyces dendrorhous

Xanthophyllomyces dendrorhous is a carotenogenic yeast with a singular metabolic capacity to produce astaxanthin, a valuable antioxidant pigment. This yeast can assimilate several carbon sources and sustain fermentation even under aerobic conditions. Since astaxanthin biosynthesis is affected by the carbon source, the study of carotenogenesis regulatory mechanisms is key for improving astaxanthin yield in X. dendrorhous This study aimed to elucidate the regulation of the metabolism of different carbon sources and the phenomenon of catabolic repression in this yeast. To this end, protein and transcript levels were quantified by iTRAQ (isobaric tags for relative and absolute quantification) and transcriptomic sequencing (RNA-seq) in the wild-type strain under conditions of glucose, maltose, or succinate treatment and in the mutant strains for genes MIG1, CYC8, and TUP1 under conditions of glucose treatment. Alternative carbon sources such as maltose and succinate affected the relative abundances of 14% of the wild-type proteins, which were mainly grouped into the carbohydrate metabolism category, with the glycolysis/gluconeogenesis and citrate cycle pathways being the most highly represented pathways. Each mutant strain showed significant proteomic profile changes, affecting approximately 2% of the total proteins identified, compared to the wild-type strain under glucose treatment conditions. Similarly to the results seen with the alternative carbon sources, the changes in the mutant strains mainly affected carbohydrate metabolism, with glycolysis/gluconeogenesis and the pentose phosphate and citrate cycle pathways being the most highly represented pathways. Our results showed convergence between carbon assimilation and catabolic repression in the strains studied. Interestingly, indications of cooperative, opposing, and overlapping processes during catabolic regulation were found. We also identified target proteins of the regulatory processes, reinforcing the likelihood of catabolic repression at the posttranscriptional level.IMPORTANCE The conditions affecting catabolic regulation in X. dendrorhous are complex and suggest the presence of an alternative mechanism of regulation. The repressors Mig1, Cyc8, and Tup1 are essential elements for the regulation of the use of glucose and other carbon sources. All play different roles but, depending on the growth conditions, can work in convergent, synergistic, and complementary ways to use carbon sources and to regulate other targets for yeast metabolism. Our results reinforced the belief that further studies in X. dendrorhous are needed to clarify a specific regulatory mechanism at the domain level of the repressors as well as its relationship with those of other metabolic repressors, i.e., the stress response, to elucidate carotenogenic regulation at the transcriptomic and proteomic levels in this yeast.

Metallopeptidase Stp1 activates the transcription factor Sre1 in the carotenogenic yeast Xanthophyllomyces dendrorhous

Xanthophyllomyces dendrorhous is a basidiomycete yeast known as a natural producer of astaxanthin, a carotenoid of commercial interest because of its antioxidant properties. Recent studies indicated that X. dendrorhous has a functional SREBP pathway involved in the regulation of isoprenoid compound biosynthesis, which includes ergosterol and carotenoids. SREBP is a major regulator of sterol metabolism and homeostasis in mammals; characterization in fungi also provides information about its role in the hypoxia adaptation response and virulence. SREBP protease processing is required to activate SREBP pathway functions in fungi. Here, we identified and described the STP1 gene, which encodes a metallopeptidase of the M50 family involved in the proteolytic activation of the transcription factor Sre1 of the SREBP pathway, in X. dendrorhous We assessed STP1 function in Δstp1 strains derived from the wild-type and a mutant of ergosterol biosynthesis that overproduces carotenoids and sterols. Bioinformatic analysis of the deduced protein predicted the presence of characteristic features identified in homologs from mammals and fungi. The Δstp1 mutation decreased yeast growth in the presence of azole drugs and reduced transcript levels of Sre1-dependent genes. This mutation also negatively affected the carotenoid- and sterol-overproducing phenotype. Western blot analysis demonstrated that Sre1 was activated in the yeast ergosterol biosynthesis mutant and that the Δstp1 mutation introduced in this strain prevented Sre1 proteolytic activation. Overall, our results demonstrate that STP1 encodes a metallopeptidase involved in proteolytic activation of Sre1 in X. dendrorhous, contributing to our understanding of fungal SREBP pathways.

Sterol Regulatory Element-Binding Protein (Sre1) Promotes the Synthesis of Carotenoids and Sterols in Xanthophyllomyces dendrorhous

Xanthophyllomyces dendrorhous is a basidiomycete yeast that synthesizes carotenoids, mainly astaxanthin, which are of great commercial interest. Currently, there are many unknown aspects related to regulatory mechanisms on the synthesis of carotenoids in this yeast. Our recent studies showed that changes in sterol levels and composition resulted in upregulation of genes in the mevalonate pathway required for the synthesis of carotenoid precursors, leading to increased production of these pigments. Sterol Regulatory Element-Binding Proteins (SREBP), called Sre1 in yeast, are conserved transcriptional regulators of sterol homeostasis and other cellular processes. Given the results linking sterols and carotenoids, we investigated the role of SREBP in sterol and carotenoid synthesis in X. dendrorhous. In this study, we present the identification and functional characterization of the X. dendrorhous SRE1 gene, which encodes the transcription factor Sre1. The deduced protein has the characteristic features of SREBP/Sre1 and binds to consensus DNA sequences in vitro. RNA-seq analysis and chromatin-immunoprecipitation experiments showed that genes of the mevalonate pathway and ergosterol biosynthesis are directly regulated by Sre1. The sre1- mutation reduced sterol and carotenoid production in X. dendrorhous, and expression of the Sre1 N-terminal domain (Sre1N) increased carotenoid production more than twofold compared to wild-type. Overall, our results indicate that in X. dendrorhous transcriptional regulation of genes in the mevalonate pathway control production of the isoprenoid derivatives, carotenoids and sterol. Our results provide new insights into the conserved regulatory functions of SREBP/Sre1 and identify pointing to the SREBP pathway as a potential target to enhance carotenoid production in X. dendrorhous.

A cross-talk on compositional dynamics and codon usage patterns of mitochondrial CYB gene in Echinodermata

Codon usage bias (CUB) refers to a phenomenon in which some synonymous codons are used in mature mRNA at a higher frequency than other members codifying the same amino acid. CUB is mainly determined by mutation pressure and natural selection. We used bioinformatic tools to analyze the protein coding sequences of mitochondrial CYB gene in different classes of Echinodermata to understand the patterns of codon usage. The ENC values of CYB gene in five different classes of Echinodermata were 41.64, 30.33, 43.63, 41.11, and 41.33, which suggested that the CUB of this gene was low. The relative synonymous codon usage (RSCU) values showed that the patterns of over-represented and under-represented codons were different among different classes. Correspondence analysis indicated that the plots of CYB gene were different across classes, suggesting that the pattern of codon usage was also different among five classes under study. Highly significant correlation (p < .01) between overall nucleotide composition and its 3rd codon position indicated that both mutational pressure and natural selection had an influence on the codon usage bias of CYB gene. Furthermore, PR-2 bias plot analysis showed that both mutation pressure and natural selection might have affected the pattern of codon usage in CYB gene of Echinodermata.

Analysis of codon usage patterns in Hirudinaria manillensis reveals a preference for GC-ending codons caused by dominant selection constraints

Background

Hirudinaria manillensis is an ephemeral, blood-sucking ectoparasite, possessing anticoagulant capacities with potential medical applications. Analysis of codon usage patterns would contribute to our understanding of the evolutionary mechanisms and genetic architecture of H. manillensis, which in turn would provide insight into the characteristics of other leeches. We analysed codon usage and related indices using 18,000 coding sequences (CDSs) retrieved from H. manillensis RNA-Seq data.

Results

We identified four highly preferred codons in H. manillensis that have G/C-endings. Points generated in an effective number of codons (ENC) plot distributed below the standard curve and the slope of a neutrality plot was less than 1. Highly expressed CDSs had lower ENC content and higher GC content than weakly expressed CDSs. Principal component analysis conducted on relative synonymous codon usage (RSCU) values divided CDSs according to GC content and divided codons according to ending bases. Moreover, by determining codon usage, we found that the majority of blood-diet related genes have undergone less adaptive evolution in H. manillensis, except for those with homologous sequences in the host species.

Conclusions

Codon usage in H. manillensis had an overall preference toward C-endings and indicated that codon usage patterns are mediated by differential expression, GC content, and biological function. Although mutation pressure effects were also notable, the majority of genetic evolution in H. manillensis was driven by natural selection.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4937-x) contains supplementary material, which is available to authorized users.

Genome-wide analysis of codon usage bias in four sequenced cotton species

Codon usage bias (CUB) is an important evolutionary feature in a genome which provides important information for studying organism evolution, gene function and exogenous gene expression. The CUB and its shaping factors in the nuclear genomes of four sequenced cotton species, G. arboreum (A₂), G. raimondii (D₅), G. hirsutum (AD₁) and G. barbadense (AD₂) were analyzed in the present study. The effective number of codons (ENC) analysis showed the CUB was weak in these four species and the four subgenomes of the two tetraploids. Codon composition analysis revealed these four species preferred to use pyrimidine-rich codons more frequently than purine-rich codons. Correlation analysis indicated that the base content at the third position of codons affect the degree of codon preference. PR2-bias plot and ENC-plot analyses revealed that the CUB patterns in these genomes and subgenomes were influenced by combined effects of translational selection, directional mutation and other factors. The translational selection (P2) analysis results, together with the non-significant correlation between GC12 and GC3, further revealed that translational selection played the dominant role over mutation pressure in the codon usage bias. Through relative synonymous codon usage (RSCU) analysis, we detected 25 high frequency codons preferred to end with T or A, and 31 low frequency codons inclined to end with C or G in these four species and four subgenomes. Finally, 19 to 26 optimal codons with 19 common ones were determined for each species and subgenomes, which preferred to end with A or T. We concluded that the codon usage bias was weak and the translation selection was the main shaping factor in nuclear genes of these four cotton genomes and four subgenomes.

Selection shapes the patterns of codon usage in three closely related species of genus Misgurnus

Neutrality plots revealed that selection probably dominates codon bias, whereas mutation plays only a minor role, in shaping the codon bias in three loaches, Misgurnus anguillicaudatus, M. mohoity, and M. bipartitus. These three species also clearly showed similar tendencies in the preferential usage of codons. Nineteen, nine, and 14 preferred codon pairs and 179, 182, and 173 avoided codon pairs were also detected in M. anguillicaudatus, M. bipartitus, and M. mohoity, respectively, and the most frequently avoided type of cP3-cA1 dinucleotide in these species was nnUAnn. The expression-linked patterns of codon usage revealed that higher expression was associated with higher GC₃, lower ENC, and a smaller proportion of amino acids with high size/complexity (S/C) scores in these three species. These results elucidate selectively driven codon bias in Misgurnus species, and reveal the potential importance of expression-mediated selection in shaping the genome evolution of fish.

Biosynthesis of Astaxanthin as a Main Carotenoid in the Heterobasidiomycetous Yeast Xanthophyllomyces dendrorhous

Carotenoids are organic lipophilic yellow to orange and reddish pigments of terpenoid nature that are usually composed of eight isoprene units. This group of secondary metabolites includes carotenes and xanthophylls, which can be naturally obtained from photosynthetic organisms, some fungi, and bacteria. One of the microorganisms able to synthesise carotenoids is the heterobasidiomycetous yeast Xanthophyllomyces dendrorhous, which represents the teleomorphic state of Phaffia rhodozyma, and is mainly used for the production of the xanthophyll astaxanthin. Upgraded knowledge on the biosynthetic pathway of the main carotenoids synthesised by X. dendrorhous, the biotechnology-based improvement of astaxanthin production, as well as the current omics approaches available in this yeast are reviewed in depth.

Characterization of the cytochrome P450 monooxygenase genes (P450ome) from the carotenogenic yeast Xanthophyllomyces dendrorhous

BACKGROUND

The cytochromes P450 (P450s) are a large superfamily of heme-containing monooxygenases involved in the oxidative metabolism of an enormous diversity of substrates. These enzymes require electrons for their activity, and the electrons are supplied by NAD(P)H through a P450 electron donor system, which is generally a cytochrome P450 reductase (CPR). The yeast Xanthophyllomyces dendrorhous has evolved an exclusive P450-CPR system that specializes in the synthesis of astaxanthin, a carotenoid with commercial potential. For this reason, the aim of this work was to identify and characterize other potential P450 genes in the genome of this yeast using a bioinformatic approach.

RESULTS

Thirteen potential P450-encoding genes were identified, and the analysis of their deduced proteins allowed them to be classified in ten different families: CYP51, CYP61, CYP5139 (with three members), CYP549A, CYP5491, CYP5492 (with two members), CYP5493, CYP53, CYP5494 and CYP5495. Structural analyses of the X. dendrorhous P450 proteins showed that all of them have a predicted transmembrane region at their N-terminus and have the conserved domains characteristic of the P450s, including the heme-binding region (FxxGxRxCxG); the PER domain, with the characteristic signature for fungi (PxRW); the ExxR motif in the K-helix region and the oxygen-binding domain (OBD) (AGxDTT); also, the characteristic secondary structure elements of all the P450 proteins were identified. The possible functions of these P450s include primary, secondary and xenobiotic metabolism reactions such as sterol biosynthesis, carotenoid synthesis and aromatic compound degradation.

CONCLUSIONS

The carotenogenic yeast X. dendrorhous has thirteen P450-encoding genes having potential functions in primary, secondary and xenobiotic metabolism reactions, including some genes of great interest for fatty acid hydroxylation and aromatic compound degradation. These findings established a basis for future studies about the role of P450s in the carotenogenic yeast X. dendrorhous and their potential biotechnological applications.

Analysis of codon usage patterns in Morus notabilis based on genome and transcriptome data

Codons play important roles in regulating gene expression levels and mRNA half-lives. However, codon usage and related studies in multicellular organisms still lag far behind those in unicellular organisms. In this study, we describe for the first time genome-wide patterns of codon bias in Morus notabilis (mulberry tree), and analyze genome-wide codon usage in 12 other species within the order Rosales. The codon usage of M. notabilis was affected by nucleotide composition, mutation pressure, nature selection, and gene expression level. Translational selection optimal codons were identified and highly expressed genes of M. notabilis tended to use the optimal codons. Genes with higher expression levels have shorter coding region and lower amino acid complexity. Housekeeping genes showed stronger translational selection, which, notably, was not caused by the large differences between the expression level of housekeeping genes and other genes.

Functional characterization of thiolase-encoding genes from Xanthophyllomyces dendrorhous and their effects on carotenoid synthesis

BACKGROUND

The basidiomycetous yeast Xanthophyllomyces dendrorhous has been described as a potential biofactory for terpenoid-derived compounds due to its ability to synthesize astaxanthin. Functional knowledge of the genes involved in terpenoid synthesis would create opportunities to enhance carotenoid production. A thiolase enzyme catalyzes the first step in terpenoid synthesis.

RESULTS

Two potential thiolase-encoding genes were found in the yeast genome; bioinformatically, one was identified as an acetyl-CoA C-acetyltransferase (ERG10), and the other was identified as a 3-ketoacyl Co-A thiolase (POT1). Heterologous complementation assays in Saccharomyces cerevisiae showed that the ERG10 gene from X. dendrorhous could complement the lack of the endogenous ERG10 gene in S. cerevisiae, thereby allowing cellular growth and sterol synthesis. X. dendrorhous heterozygous mutants for each gene were created, and a homozygous POT1 mutant was also obtained. This mutant exhibited changes in pigment composition and higher ERG10 transcript levels than the wild type strain.

CONCLUSIONS

The results support the notion that the ERG10 gene in X. dendrorhous is a functional acetyl-CoA C-acetyltransferase essential for the synthesis of mevalonate in yeast. The POT1 gene would encode a functional 3-ketoacyl Co-A thiolase that is non-essential for cell growth, but its mutation indirectly affects pigment production.

Regulation of carotenogenesis in the red yeast Xanthophyllomyces dendrorhous: the role of the transcriptional co-repressor complex Cyc8-Tup1 involved in catabolic repression

BACKGROUND

The yeast Xanthophyllomyces dendrorhous produces carotenoids of commercial interest, including astaxanthin and β-carotene. Although carotenogenesis in this yeast and the expression profiles of the genes controlling this pathway are known, the mechanisms regulating this process remain poorly understood. Several studies have demonstrated that glucose represses carotenogenesis in X. dendrorhous, suggesting that this pathway could be regulated by catabolic repression. Catabolic repression is a highly conserved regulatory mechanism in eukaryotes and has been widely studied in Saccharomyces cerevisiae. Glucose-dependent repression is mainly observed at the transcriptional level and depends on the DNA-binding regulator Mig1, which recruits the co-repressor complex Cyc8-Tup1, which then represses the expression of target genes. In this work, we studied the regulation of carotenogenesis by catabolic repression in X. dendrorhous, focusing on the role of the co-repressor complex Cyc8-Tup1.

RESULTS

The X. dendrorhous CYC8 and TUP1 genes were identified, and their functions were demonstrated by heterologous complementation in S. cerevisiae. In addition, cyc8 - and tup1 - mutant strains of X. dendrorhous were obtained, and both mutations were shown to prevent the glucose-dependent repression of carotenogenesis in X. dendrorhous, increasing the carotenoid production in both mutant strains. Furthermore, the effects of glucose on the transcript levels of genes involved in carotenogenesis differed between the mutant strains and wild-type X. dendrorhous, particularly for genes involved in the synthesis of carotenoid precursors, such as HMGR, idi and FPS. Additionally, transcriptomic analyses showed that cyc8 - and tup1 - mutations affected the expression of over 250 genes in X. dendrorhous.

CONCLUSIONS

The CYC8 and TUP1 genes are functional in X. dendrorhous, and their gene products are involved in catabolic repression and carotenogenesis regulation. This study presents the first report involving the participation of Cyc8 and Tup1 in carotenogenesis regulation in yeast.

Genome-wide comparative analysis of codon usage bias and codon context patterns among cyanobacterial genomes

With the increasing accumulation of genomic sequence information of prokaryotes, the study of codon usage bias has gained renewed attention. The purpose of this study was to examine codon selection pattern within and across cyanobacterial species belonging to diverse taxonomic orders and habitats. We performed detailed comparative analysis of cyanobacterial genomes with respect to codon bias. Our analysis reflects that in cyanobacterial genomes, A- and/or T-ending codons were used predominantly in the genes whereas G- and/or C-ending codons were largely avoided. Variation in the codon context usage of cyanobacterial genes corresponded to the clustering of cyanobacteria as per their GC content. Analysis of codon adaptation index (CAI) and synonymous codon usage order (SCUO) revealed that majority of genes are associated with low codon bias. Codon selection pattern in cyanobacterial genomes reflected compositional constraints as major influencing factor. It is also identified that although, mutational constraint may play some role in affecting codon usage bias in cyanobacteria, compositional constraint in terms of genomic GC composition coupled with environmental factors affected codon selection pattern in cyanobacterial genomes.

The Involvement of Mig1 from Xanthophyllomyces dendrorhous in Catabolic Repression: An Active Mechanism Contributing to the Regulation of Carotenoid Production

The red yeast X. dendrorhous is one of the few natural sources of astaxanthin, a carotenoid used in aquaculture for salmonid fish pigmentation and in the cosmetic and pharmaceutical industries for its antioxidant properties. Genetic control of carotenogenesis is well characterized in this yeast; however, little is known about the regulation of the carotenogenesis process. Several lines of evidence have suggested that carotenogenesis is regulated by catabolic repression, and the aim of this work was to identify and functionally characterize the X. dendrorhous MIG1 gene encoding the catabolic repressor Mig1, which mediates transcriptional glucose-dependent repression in other yeasts and fungi. The identified gene encodes a protein of 863 amino acids that demonstrates the characteristic conserved features of Mig1 proteins, and binds in vitro to DNA fragments containing Mig1 boxes. Gene functionality was demonstrated by heterologous complementation in a S. cerevisiae mig1- strain; several aspects of catabolic repression were restored by the X. dendrorhous MIG1 gene. Additionally, a X. dendrorhous mig1- mutant was constructed and demonstrated a higher carotenoid content than the wild-type strain. Most important, the mig1- mutation alleviated the glucose-mediated repression of carotenogenesis in X. dendrorhous: the addition of glucose to mig1- and wild-type cultures promoted the growth of both strains, but carotenoid synthesis was observed only in the mutant strain. Transcriptomic and RT-qPCR analyses revealed that several genes were differentially expressed between X. dendrorhous mig1- and the wild-type strain when cultured with glucose as the sole carbon source. The results obtained in this study demonstrate that catabolic repression in X. dendrorhous is an active process in which the identified MIG1 gene product plays a central role in the regulation of several biological processes, including carotenogenesis.

Analysis of codon usage bias of envelope glycoprotein genes in nuclear polyhedrosis virus (NPV) and its relation to evolution

BACKGROUND

Analysis of codon usage bias is an extremely versatile method using in furthering understanding of the genetic and evolutionary paths of species. Codon usage bias of envelope glycoprotein genes in nuclear polyhedrosis virus (NPV) has remained largely unexplored at present. Hence, the codon usage bias of NPV envelope glycoprotein was analyzed here to reveal the genetic and evolutionary relationships between different viral species in baculovirus genus.

RESULTS

A total of 9236 codons from 18 different species of NPV of the baculovirus genera were used to perform this analysis. Glycoprotein of NPV exhibits weaker codon usage bias. Neutrality plot analysis and correlation analysis of effective number of codons (ENC) values indicate that natural selection is the main factor influencing codon usage bias, and that the impact of mutation pressure is relatively smaller. Another cluster analysis shows that the kinship or evolutionary relationships of these viral species can be divided into two broad categories despite all of these 18 species are from the same baculovirus genus.

CONCLUSIONS

There are many elements that can affect codon bias, such as the composition of amino acids, mutation pressure, natural selection, gene expression level, and etc. In the meantime, cluster analysis also illustrates that codon usage bias of virus envelope glycoprotein can serve as an effective means of evolutionary classification in baculovirus genus.

Molecular characterization and heterologous expression of a Xanthophyllomyces dendrorhous α-glucosidase with potential for prebiotics production

Basidiomycetous yeast Xanthophyllomyces dendrorhous expresses an α-glucosidase with strong transglycosylation activity producing prebiotic sugars such as panose and an unusual tetrasaccharides mixture including α-(1-6) bonds as major products, which makes it of biotechnological interest. Initial analysis pointed to a homodimeric protein of 60 kDa subunit as responsible for this activity. In this study, the gene Xd-AlphaGlu was characterized. The 4131-bp-long gene is interrupted by 13 short introns and encodes a protein of 990 amino acids (Xd-AlphaGlu). The N-terminal sequence of the previously detected 60 kDa protein resides in this larger protein at residues 583-602. Functionality of the gene was proved in Saccharomyces cerevisiae, which produced a protein of about 130 kDa containing Xd-AlphaGlu sequences. All properties of the heterologously expressed protein, including thermal and pH profiles, activity on different substrates, and ability to produce prebiotic sugars were similar to that of the α-glucosidase produced in X. dendrorhous. No activity was detected in S. cerevisiae containing exclusively the 1256-bp from gene Xd-AlphaGlu that would encode synthesis of the 60 kDa protein previously detected. Data were compatible with an active monomeric α-glucosidase of 990 amino acids and an inactive hydrolysis product of 60 kDa. Protein Xd-AlphaGlu contained most of the elements characteristic of α-glucosidases included in the glycoside hydrolases family GH31 and its structural model based on the homologous human maltase-glucoamylase was obtained. Remarkably, the Xd-AlphaGlu C-terminal domain presents an unusually long 115-residue insertion that could be involved in this enzyme's activity against long-size substrates such as maltoheptaose and soluble starch.

Molecular Characterization and Functional Analysis of Cytochrome b5 Reductase (CBR) Encoding Genes from the Carotenogenic Yeast Xanthophyllomyces dendrorhous

The eukaryotic microsomal cytochrome P450 systems consist of a cytochrome P450 enzyme (P450) and a cytochrome P450 redox partner, which generally is a cytochrome P450 reductase (CPR) that supplies electrons from NADPH. However, alternative electron donors may exist such as cytochrome b5 reductase and cytochrome b5 (CBR and CYB5, respectively) via, which is NADH-dependent and are also anchored to the endoplasmic reticulum. In the carotenogenic yeast Xanthophyllomyces dendrorhous, three P450-encoding genes have been described: crtS is involved in carotenogenesis and the CYP51 and CYP61 genes are both implicated in ergosterol biosynthesis. This yeast has a single CPR (encoded by the crtR gene), and a crtR^- mutant does not produce astaxanthin. Considering that this mutant is viable, the existence of alternative cytochrome P450 electron donors like CBR and CYB5 could operate in this yeast. The aim of this work was to characterize the X. dendrorhous CBR encoding gene and to study its involvement in P450 reactions in ergosterol and carotenoid biosynthesis. Two CBRs genes were identified (CBR.1 and CBR.2), and deletion mutants were constructed. The two mutants and the wild-type strain showed similar sterol production, with ergosterol being the main sterol produced. The crtR^- mutant strain produced a lower proportion of ergosterol than did the parental strain. These results indicate that even though one of the two CBR genes could be involved in ergosterol biosynthesis, crtR complements their absence in the cbr^- mutant strains, at least for ergosterol production. The higher NADH-dependent cytochrome c reductase activity together with the higher transcript levels of CBR.1 and CYB5 in the crtR^- mutant as well as the lower NADH-dependent activity in CBS-cbr.1^- strongly suggest that CBR.1-CYB5 via participates as an alternative electron donor pathway for P450 enzymes involved in ergosterol biosynthesis in X. dendrorhous.