Next Article in Journal
Comparative Transcriptome Profiling of Skeletal Muscle from Black Muscovy Duck at Different Growth Stages Using RNA-seq
Next Article in Special Issue
Whole Genome Analysis of Environmental Pseudomonas mendocina Strains: Virulence Mechanisms and Phylogeny
Previous Article in Journal
Expression Profile of Porcine TRIM26 and Its Inhibitory Effect on Interferon-β Production and Antiviral Response
Previous Article in Special Issue
The New Klebsiella pneumoniae ST152 Variants with Hypermucoviscous Phenotype Isolated from Renal Transplant Recipients with Asymptomatic Bacteriuria—Genetic Characteristics by WGS
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Metabolic Specialization and Codon Preference of Lignocellulolytic Genes in the White Rot Basidiomycete Ceriporiopsis subvermispora

1
Laboratorio de Microbiología Ambiental y Extremófilos, Departamento de Ciencias Biológicas y Biodiversidad, Universidad de los Lagos, Osorno 5290000, Chile
2
Instituto de Ciencias Biomédicas, Facultad de Ciencias de la Salud, Universidad Autónoma de Chile, Santiago 8910132, Chile
3
Departamento de Bioquímica y Biología Molecular, Facultad de Ciencias Químicas y Farmacéuticas, Universidad de Chile, Santiago 8380492, Chile
4
Laboratorio de Metagenómica Bacteriana, Centro de Biotecnología Acuícola, Facultad de Química y Biología, Universidad de Santiago de Chile, Santiago 9170002, Chile
*
Author to whom correspondence should be addressed.
Genes 2020, 11(10), 1227; https://doi.org/10.3390/genes11101227
Submission received: 11 September 2020 / Revised: 13 October 2020 / Accepted: 14 October 2020 / Published: 20 October 2020
(This article belongs to the Special Issue Microbial Genomics and Evolution)

Abstract

:
Ceriporiopsis subvermispora is a white-rot fungus with a high specificity towards lignin mineralization when colonizing dead wood or lignocellulosic compounds. Its lignocellulose degrading system is formed by cellulose hydrolytic enzymes, manganese peroxidases, and laccases that catalyze the efficient depolymerization and mineralization of lignocellulose. To determine if this metabolic specialization has modified codon usage of the lignocellulolytic system, improving its adaptation to the fungal translational machine, we analyzed the adaptation to host codon usage (CAI), tRNA pool (tAI, and AAtAI), codon pair bias (CPB), and the number of effective codons (Nc). These indexes were correlated with gene expression of C. subvermispora, in the presence of glucose and Aspen wood. General gene expression was not correlated with the index values. However, in media containing Aspen wood, the induction of expression of lignocellulose-degrading genes, showed significantly (p < 0.001) higher values of CAI, AAtAI, CPB, tAI, and lower values of Nc than non-induced genes. Cellulose-binding proteins and manganese peroxidases presented the highest adaptation values. We also identified an expansion of genes encoding glycine and glutamic acid tRNAs. Our results suggest that the metabolic specialization to use wood as the sole carbon source has introduced a bias in the codon usage of genes involved in lignocellulose degradation. This bias reduces codon diversity and increases codon usage adaptation to the tRNA pool available in C. subvermispora. To our knowledge, this is the first study showing that codon usage is modified to improve the translation efficiency of a group of genes involved in a particular metabolic process.

1. Introduction

The main carbon source synthesized through photosynthesis which plays a central role in the carbon cycle of the planet is lignocellulose. Lignin, one of the compounds of lignocellulose, is a recalcitrant, aromatic, and amorphous polymer that protects lignocellulose from microbial attack. A small group of filamentous fungi from the basidiomycete phylum is unique in its ability to efficiently degrade lignocellulose [1]. Together they are collectively known as white rot fungi, developing an enzymatic machinery that allows degradation of the three main components of lignocellulose: lignin, cellulose, and hemicellulose. White-rot fungi mineralize lignin as a strategy to access cellulose and hemicellulose, whose sugar moieties are used as carbon and energy sources. Degradation of these aromatic and carbohydrate polymers progresses by different mechanisms. While mineralization of lignin is carried out by free radicals generated enzymatically [2,3], the degradation of cellulose and hemicellulose into its constituent sugars occurs through the combination of direct enzymatic hydrolysis and partial hydrolysis by enzymatically generated free radicals [4].
Ceriporiopsis subvermispora is a lignocellulose degrading white-rot basidiomycete that mineralizes lignin using a machinery composed of manganese peroxidase, versatile peroxidase, lignin peroxidase, and laccase, as well as an accessory system generating lipoperoxides [5,6,7,8,9]. Biochemical and genetic evidence shows that manganese peroxidases and laccases are encoded by a family of structurally related genes [6,8,10,11,12]. Hydrolysis of cellulose and hemicellulose is also achieved by a suite of cellulases and hemicellulases, which are encoded by multigene families [8]. Gene expression of enzymes involved in lignocellulose digestion and mineralization is induced by compounds present in wood, as shown by research using Northern blot, RT-PCR, and microarrays [8,13,14,15,16,17]. These studies also indicate that expression levels differ among members of a given multigene family [8,11,18]. Although selective pressure imposed by lignocellulose metabolism is reflected by an increase in the copy number of lignocellulolytic genes in C. subvermispora and other white-rot fungi [8], it is unknown whether this condition has influenced the composition of synonymous codons in these genes.
Preference for use of some synonymous codons (codon bias) is influenced by two mechanisms, mutational bias, and translational selection. Mutational bias establishes that the codon preference of a gene is determined by the GC content of the organism or the region where the gene is located in the genome [19]. On the other hand, translational selection establishes that codon bias is determined by the influence of synonymous mutations on the translational process [20,21]. Highly expressed genes, such as those coding for ribosomal proteins, tend to have a greater bias in the use of synonymous codons, preferentially using codons which are most represented in the genome [22]. In model organisms such as Escherichia coli and Saccharomyces cerevisiae, tRNAs that decode frequently used codons tend to have a larger gene copy number [23,24]. Thus, genes with frequently used codons are more adapted to the tRNA pool and consequently have more tRNAs available for their translation, which in turn improves the efficiency of protein synthesis.
In this work, we compared codon usage, adaptation to the tRNA pool, and codon bias of C. subvermispora genes involved in lignocellulose metabolism and genes related to other cell processes. Our results show that selective pressure imposed by the use of lignocellulose has specifically modified codon usage of genes involved in the utilization of lignocellulose, favoring an increase in translational efficiency with respect to genes not involved in this process.

2. Materials and Methods

2.1. Sequences

Genome sequences from C. subvermispora were downloaded from https://mycocosm.jgi.doe.gov/Cersu1/Cersu1.home.html.

2.2. Microarray Data

Microarray data were downloaded from the Gene Expression Omnibus Database, accession no. GSE34636 [8].

2.3. Determination of the Gene Copy Number Coding for tRNAs

Genes encoding tRNAs were identified using the tRNAscan-SE program, using unmasked scaffolds assembled from the genome sequence of C. subvermispora [8,25].

2.4. Determination of Codons and Codon Pair Frequencies

The frequency of codon usage and codon pairs was determined from the coding sequences of C. subvermispora. To calculate the frequency of codon usage the program JEMBOSS through the CUSP routine was used (http://emboss.sourceforge.net/Jemboss/) [26]. A written script in Python language was used for determining the frequency of codon pairs [27].

2.5. Determination of Bias in Codon-Pair Usage

Bias in codon-pair usage was determined as described by Coleman [28]. The observed frequency of amino acid pairs was deduced based on the sum of all the codon pairs that codify for the same combination of amino acids. Subsequently, we used these frequencies to calculate the CPS index (codon-pair score) using the equation described by Mueller et al. 2006 [29] (Equation (1)), where f(AB) corresponds to the frequency of the appearance of the AB codon pair expressed in parts-per-thousands with respect to the total number of codon pairs in the genome. f(A) and f(B) correspond to the abundance of each of the individual codons expressed in parts-per-thousands with respect to the total number of codons in the genome. f(X) and f(Y) correspond to the amino acids codified by the A and B codons, respectively, and f(XY) is the abundance of the amino acid pairs in all the proteins of the organism under study. The CPS values of the 3271 codon pairs in C. subvermispora can be found in the Supplementary Material (Table S1). We calculated the bias in codon usage as the arithmetic mean of the codon-pair score of a gene (CPS) as described by Coleman et al. (2006) (Equation (2)), in which i is the ith codon pair of a gene and 1 is the total length of the gene expressed in the quantity of codified amino acids.
C P S = L n   ( f ( A B ) f ( A ) × f ( B ) × f ( X ) × f ( Y ) f ( X Y ) )
C P B = i l C P S i l 1

2.6. Calculation of Relative Synonymous Codon Usage (RSCU) and Codon Adaptation Index (CAI)

The Relative Synonymous Codon Usage (RSCU) and Codon Adaptation Index (CAI) were calculated with the Emboss program [30] based on the frequency of codon usage of C. subvermispora.

2.7. Determination of Adaptation to the tRNA Pool

The adaptation to the tRNA pool of genes present in C. subvermispora was determined by calculating the values of tAI and AAtAI [31]. We calculated tAI as established by dos Reis [24] (Equations (3) and (4)), estimating the relative abundance of tRNAs from the number of copies of each tRNA gene in the C. subvermispora genome. The number of copies of each tRNA gene was determined with the tRNAscan-SE program [25] in Linux, which was used to analyze the unmasked assembled C. subvermispora genome.
W i = j = 1 n i ( 1 S i j ) t G N C i j
Wi is the relative adaptiveness of the ith codon to the tRNA pool, ni is the number of tRNA isoacceptors that recognize the ith codon, tGCNij is the number of copies of the jth tRNA gene that recognizes the ith codon, and Sij is the selective constraint in the efficiency of codon-anticodon pairing.
t A I g = ( k = 1 l g w i k g ) 1 / l g
The adaptation of a gene to the tRNA pool is calculated according to Equation (4), in which wi is defined as the ratio between Wi and Wmax (Wi/Wmax), ikg is the codon defined by the kth triplet of gene g, and lg is the length of gene g in codifying codons.

2.8. Determination of AAtAI

Briefly, AAtAI was calculated using Equation (5).
A A t A I g = ( k = 1 l g w i k g A A ) 1 / l g
AAwi is defined as the ratio between Wi and WAAmax (Wi/WAAmax), where WAAmax is the highest value of Wi among codons coding for the same amino acid. ikg is the codon defined by the kth triplet of gene g, and lg is the length of gene g in codifying codons. As AAtAI is calculated from Equation (5), which is similar to the calculation of CAI, we used the EMBOSS software [30], entering Wi data in replacement of the frequency of codon usage. Wi was calculated using the procedure described for tAI in Equation (3).

2.9. Determination of effective Number of Codons (Nc)

The number of effective codons for the C. subvermispora genes was calculated using the CodonW program (http://codonw.sourceforge.net/).

2.10. Phylogenetic Analysis

Multialignment between sequences of tRNA genes was performed using ClustalW [32]. The parameters were set up to align sequences using IUB as substitution matrix. The evolutionary history was inferred using the Neighbor-Joining method [33]. Interior branch test with a bootstrap of 1000 was used to analyze confidence of the tree [34]. The evolutionary distances were computed using the Maximum Composite Likelihood method [35]. The rate of variation among sites was modeled with a γ distribution (shape parameter = 1). Multialignment and evolutionary analyses were conducted in MEGA5 [36].

2.11. Graphs and Statistical Methods

The program SigmaPlot 11 was used for graphs and statistical tests. The significance of the differences or correlations among the data groups obtained was evaluated with the Rank Sum Test non-parametric test for comparing two groups and the non-parametric Spearman Rank Order test for correlations, using a p-value < 0.05 as a cutoff.

3. Results

3.1. C. subvermispora tRNAs

Genome analysis of C. subvermispora by tRNAScan-SE identified a total of 192 tRNAs in 32 scaffolds (Figure 1). About 72% of the tRNA genes presented introns (Table 1). The scaffold with the highest number of tRNA genes was scaffold 1, which contains 20 copies of various tRNA genes (Table S2). The tRNA with the highest number of gene copies corresponded to tRNAs charging glycine, with 17 gene copies distributed in eight scaffolds (scaffolds 1, 2, 5, 7, 9, 10, 13, and 20). The tRNAs for cysteine and tryptophan presented the lowest number of gene copies, each with three gene copies in three scaffolds (Table 1).

3.2. Phylogenetic Analysis of C. subvermispora tRNA Genes

To determine whether the high number of tRNA genes charging the same amino acid corresponds to related gene copies, a phylogenetic reconstruction and evolutionary distance calculation were performed using the tRNA sequences identified by the tRNAscan-SE program. Phylogenetic reconstruction indicated that most tRNA genes that code for the same amino acid group together, with the exception of the tRNA genes that charge arginine, valine, and alanine (Figure 1), which form two groups in each case. For tRNAs loading arginine, Group I comprises genes presenting anticodons with the WCG sequence. In contrast, group II presents a YCK anticodon consensus sequence. In tRNAs that load valine, Group I corresponds to tRNA genes without introns, while group II includes all valine tRNA genes containing introns. Genes coding for group I of alanine tRNAs exhibit anticodons with the consensus sequence YGC, whereas, the anticodon sequence is AGC for group II (Figure 1). The tRNA genes 66 and 155 also show a different pattern. tRNAScan prediction indicates that the amino acid loaded by tRNA155 should be serine, however, the sequence of this gene grouped with threonine charging tRNA genes. Additionally, the tRNA66 gene is expected to load isoleucine, though this gene does not group with isoleucine charging tRNAs (Figure 1). To identify tRNA genes that are repeated, the evolutionary distance between different tRNA genes was calculated. Genes with values of evolutionary distance equal to 0.000 were selected. This analysis identified 15 tRNA genes that are repeated between two and ten times. The group of tRNA genes that showed the greatest expansion corresponds to those tRNAs carrying glutamic acid and glycine. One glycine tRNA gene is repeated twice and the other is repeated ten times. The tRNA gene for glutamic acid is also repeated ten times (Table S3). These genes are scattered along the genome, with the exception of tRNA genes 91, 92, 94, 95, 96, and 97, which code for glycine and are located at adjacent positions.

3.3. tRNA Abundance and Codon Usage in C. subvermispora

The expansion of certain families of tRNAs in the genome of C. subvermispora could be the result of an evolutionary pressure to increase their expression. In some organisms such as E. coli and yeast, the number of copies of a tRNA gene is proportional to the abundance in the genome of the decoded codon [24]. This proportionality is explained because during translation process, aminoacyl-tRNAs that decode frequently used codons have a higher rate of consumption. To sustain adequate translation, cells must balance synthesis and consumption rates of aminoacyl-tRNAs. The increase in copy number of a gene that encodes a tRNA that recognizes a frequently used codon, allows increasing expression of this tRNA and its aminoacylated form to balance its rate of consumption.
To test if the expansion of certain tRNA families in C. subvermispora is related to the increment of the specific codons, we assessed whether there is a correlation between the frequency of codon usage and the amount of tRNA genes that decode these most highly used codons. We identify a positive correlation between these parameters (ρ = 0.406, p = 0.0016, n = 61). Dos Reis et al. show that the Relative Adaptiveness to the tRNA pool (w) which takes in account that codons can be recognized by anticodons with perfect or imperfect match (wobble codon-anticodon recognition rules) with different affinities is a better parameter to measure the adaptation of a codon to their decoding tRNAs than the absolute number of tRNA [24]. When the frequency of codon usage was correlated with the Relative Adaptiveness to the tRNA pool (w), an increased correlation was observed (ρ = 0.459, p = 2.2 × 10−4, n = 61) (Figure 2).
Synonymous codons are not used equally in an organism. As one tRNA can decode several synonymous codons with different affinities, the expansion of some tRNA families in C. subvermispora may be related to the preferential use of certain synonymous codons in coding regions of C. subvermispora. To assess this hypothesis, we correlated RSCU values with Relative Adaptiveness to the tRNA pool. Non-statistical correlation was observed, in part because w values are normalized with respect to the tRNA with the highest number of genes able to decode it (Wi/Wmax) and not with respect to the pool of tRNAs that decode the complete set of synonymous codons. When W values were normalized with respect to the total amount of tRNAs that decode a set of synonymous codons, a strong correlation with the RSCU values (ρ = 0.628, p = 0, n = 61) was observed. This suggests that among synonymous codons, those highly represented in the C. subvermispora genome tend to be decoded by those tRNAs with a high number of gene copies.

3.4. Relationship between Gene Expression Level, Codon Bias, and Translational Efficiency in C. subvermispora

The bias in codon usage and adaptation to the tRNA pool modulates translational efficiency. Thus, highly expressed genes tend to use codons that are over-represented in the genome, which in turn present greater availability of tRNAs [24]. To determine whether this relationship exists in C. subvermispora, expression levels of C. subvermispora genes were correlated with their adaptation values to the tRNA pool and to codon bias. Adaptation to the tRNA pool was evaluated using two indexes: (i) tAI, which measures adaptation of the tRNA pool compared to the relative amount of each tRNA gene, and (ii) AAtAI, which evaluates whether a gene preferentially uses the most abundant tRNA charging a particular amino acid. Codon bias was analyzed by calculation of CAI (Codon Adaptation Index), and CPB (Codon Pair-Bias), to assess if gene expression is correlated with bias in usage of codons or of codon pairs. Codon bias also was evaluated using Nc value to determine if the transcription level is associated with a decrease in the diversity of codons used. We analyzed expression levels determined by RNAseq published by Hori et al. 2014 which was obtained from C. subvermispora grown on Ball-Milled Aspen medium [17]. CAI, tAI, AAtAI, and CPB showed high degrees of correlation among them (Table S4), however, these indicators exhibit very low correlation coefficients with expression levels reported by Hori. The high correlation between the different indicators of codon bias and translational efficiency indicates that the frequency with which codons are used in C. subvermispora is related to the abundance in the genomes of the tRNAs that decode these codons.

3.5. Transcriptional Response to Growth on Ball-Milled Aspen (BMA), Codon Bias, and Translational Efficiency

Growth of C. subvermispora in natural environments is dependent on wood. In 2012, Fernandez-Fueyo et al. [8] reported microarray experiments that compared gene expression of C. subvermispora grown on glucose and on Ball-Milled Aspen (BMA) as carbon sources. Saline media with BMA has been used as a laboratory medium that mimics growth on wood to analyze expression of genes that are transcriptionally regulated by growth on wood. To analyze if genes regulated by conditions that mimic growth on wood, such as BMA have a different adaptation to the tRNA pool or codon bias, we used the microarray data published by Fernandez-Fueyo [8] and defined four groups of genes: group A corresponds to genes where expression was reduced at least 2 times with a p-value lower than 0.05. Group B includes all genes which showed increased expression of at least 2 times with a p-value lower than 0.05. Group C corresponds to all genes with non-significant differences (p > 0.05) and group D contains all genes with low changes in expression (<2 fold) that are statistically significant. (p < 0.05). Our results show that group B has lower values of Nc and higher values of CAI, tAI, AAtAI, and CPB. Groups A, C, and D show non-significant differences among them. This implies that genes induced by wood preferentially use a reduced set of codons that are better adapted to the tRNA pool present in C. subvermispora (Figure 3).
When we correlated the CAI, tAI, AAtAI, Nc, and CPB indexes with the ratio between expression in BMA and glucose culture medium, a statistically significant correlation was observed. Positive correlations were found with almost all indicators used (CAI, tAI, AAtAI, and CPB), the exception was Nc that showed a negative correlation. The higher correlation was identified in CAI and tAI indexes, and in genes that showed significant differences for expression in BMA saline medium compared to expression in glucose-supplemented saline medium (Table 2). This increase in correlation coefficients can be explained if growth on lignocellulose exerts pressure on codon usage of genes involved in the metabolization of this carbon source, thereby selecting those codons that increase the translational efficiency of these genes.
Interestingly, when genes from Group B were sorted according to their codon usage adaptation values or to the tRNA pool, we found that genes coding for ribosomal proteins presented the highest CPB, tAI, and AAtAI values (Table S5). An increase in the expression of ribosomal proteins may lead to improved ribosome biogenesis, which in turn increases the overall protein biosynthetic capacity, as observed in yeast growing in rich media [37,38]. Thus, exposure of C. subvermispora to wood or lignocellulose might lead to an increase of the overall translation rate enhancing the synthesis of proteins related to lignocellulose metabolisms.

3.6. Translational Efficiency and Codon Bias in Lignocellulolytic Genes

Genome analysis of C. subvermispora indicates that this organism presents an expansion of the number of genes directly related to the mineralization and hydrolysis of lignocellulose. The genome contains 16 annotated genes of ligninolytic peroxidases (13 Manganese peroxidase, one Versatile peroxidase, one lignin peroxidase and one generic peroxidase [9]), seven genes coding for laccases, and 14 genes coding for proteins containing a cellulose-binding domain. Moreover, it also shows an expansion of the auxiliary enzymes required for lignin degradation with four genes for cellobiose dehydrogenases (CDH), five genes for Δ-12 dehydrogenases, four genes for Δ-9 dehydrogenases, five genes for Aryl-alcohol oxidase, four genes for Methanol oxidases, two genes for Aryl alcohol dehydrogenases, three genes for copper radical oxidases, and 14 for glucose methanol choline oxidoreductase [7,8]. Genes belonging to the same family exhibit differential expression, which might reflect that they serve slightly different functions in the mineralization/hydrolysis process of lignocellulose [8]. To determine whether genes of the same family present a similar bias of codon usage and of codon adaptation or adaptation to the tRNA pool, we compared the values of Nc, tAI, AAtAI, CAI, and CPB of these genes. These values were also normalized with respect to the mean and standard deviation of the values obtained for all genes encoded in C. subvermispora (Z-values) [39]. This normalization was applied to identify if and how values for lignocellulolytic genes differ with respect to these same values in genes not directly related to the lignocellulose degradation processes (Table 3).
Genes encoding manganese peroxidase generally show above-average adaptation values to the tRNA pool and codon usage (Z-value > 0). The manganese peroxidase gene with the highest level of expression (transcript ID: 50297) proved to be the most adapted to the tRNA pool, with tAI values that are more than two standard deviations from the mean tAI values of other C. subvermispora genes (Z-tAI > 2). We also observed that the most highly induced manganese peroxidase gene in BMA medium (transcript ID:129418) also shows a high adaptation value to the tRNA pool (Z-tAI = 1.619). Genes encoding laccases showed a similar trend, as the single gene that significantly changed its expression levels after growth in BMA medium (transcript ID: 130783) showed an above average value of tAI (Z-tAI = 0.461).
Among genes encoding for proteins with cellulose binding domains, which includes cellulases and other enzyme that bind and hydrolyzes polymers related to cellulose, the gene encoding for cellulase GH7-CBM1 (transcript ID: 148588) showed the highest adaptation value to the tRNA pool (Z-tAI = 1.756) and high expression level under growth with glucose as the sole carbon source. This gene also increased its expression in BMA medium. Additionally, several genes coding for cellulases GH7-CBM1 and GH5-CBM1 showed Nc values below 40, indicating a strong bias in the use of synonymous codons (Table 3). Cellulases from the GH12 family also showed high CAI values. However, in this group of genes no association between expression levels or adaptation to the tRNA pool or any other index used in this work was found.
Genes for CDH show little bias in codon usage (Nc ~ 50). Only the CDH gene which is induced by BMA (transcript ID: 84792) shows an Nc value of 48 and slightly higher than mean (Z > 0) adaptation values to the tRNA pool (Table 3). The Δ12-dehydrogenase genes showed a weak codon bias, with the exception of gene with transcript ID: 124050; this particular gene showed a Nc value of 38 and an above average adaptation to the tRNA pool. This gene also showed strong expression in glucose-containing medium, but its expression was not modified in BMA medium. The opposite behavior was observed in the 9-dehydrogenase genes, where the gene with the greatest induction in BMA (transcript ID: 129048) showed a lower Nc value (41.46), suggesting a strong bias in codon usage (Table 3). Regarding auxiliary genes encoding for Aryl-alcohols oxidase, Methanol oxidase, Aryl-alcohol dehydrogenase and copper radical oxidases, those induced by BMA show CAI and tAI values higher than the average (Z-values > 1). An exception were the genes encoding for Methanol oxidase (transcript ID 151964) and Aryl-alcohol dehydrogenase (transcript ID: 126785) which were not induced by BMA but shown Z-values > 1.5 (Table S6).
When genes related to mineralization and digestion of lignocellulose were arranged according to their adaptation values to codon usage or the tRNA pool, genes with greater adaptation values encoded ligninolytic peroxidases, and proteins with cellulose-binding domains. Interestingly, manganese peroxidases are more adapted to the tRNA pool, while proteins with cellulose-binding domains showed higher adaptation to codon usage.
High CAI values with tAI values lower than the optimal should indicate the use of frequent codons with low availability of tRNAs to decode them. A reduction of the speed of translation has been associated with the use of rare codons or codons decoded by tRNAs with low availability. This reduction of speed allows the proper folding of the nascent protein. Proteins containing cellulose-binding domains show a complex structure, where a proper folding of at least two domains should require some reduction in the speed of translation [40,41]. Ribosome profiling experiments could help to define if these proteins require coordination between ribosome translocation and protein folding.

4. Discussion

The development of massive sequencing technologies, bioinformatics sequence analysis and synthetic biology have established that synonymous mutations, far from being silent, play an important role in the fine-tuning of protein synthesis efficiency and the role of different functional forms [42,43]. The first clues or hints of the importance of synonymous mutations for translational efficiency arose from the identification of bias in the use of synonymous codons that were detected in highly expressed ribosomal proteins from E. coli, Bacillus subtilis, and S. cerevisiae [22,44,45]. In viruses, genes coding for highly required proteins, such as those of the virus capsid or nucleoprotein, show higher adaptation values to codon usage of the host than other viral genes [39,46]. The generation of synthetic genes has revealed that the amount of synthesized protein can vary over thousand fold solely by changing the composition of synonymous codons [47]. Synthetic viruses, constructed from highly virulent viruses, show reduced replication and very low host mortality when their codon usage is changed to codon combinations present in low frequency in the host genome [48]. Further, the arrangement of codons within a gene is not random and is related to the proper folding of the nascent proteins and also to the proper recycling of ribosomes [31]. For example, the low adaptation level of the FRQ gene of Neurospora crassa to codon usage is essential to maintain the rhythm of the circadian clock [49]. Continuously increasing evidence supports the fact that bias in the use of synonymous codons plays an active role in the fine-tuning of protein production. The relationship between the use of synonymous codons and translation efficiency lies in the abundance of cognate tRNAs. Highly used codons tend to have a greater number of tRNAs that recognize them [50]. In this work, we have identified 192 genes encoding tRNAs in the white-rot basidiomycete C. subvermispora, a number that is similar to the number of tRNA genes present in other fungi such as Aspergillus fumigatus (178 genes) and Schizosaccharomyces pombe (186 genes) (http://gtrnadb.ucsc.edu/). Moreover, C. subvermispora tRNA genes present the interesting feature that each type of tRNA gene (i. e. that loads the same amino acid), groups in a different clade, indicating that the current pool of tRNAs present in the C. subvermispora monokaryotic strain B has arisen from gene duplication processes, and probably not from horizontal transfer or recombination. This result is consistent with the absence of a known C. subvermispora sexual stage [12]. Interestingly, we identified an expansion of the tRNA genes coding for amino acids glycine and glutamine. The high similarity of their sequences indicates that this expansion may have occurred recently. The expansion of this set of tRNA genes could be the result of the presence of some unidentified SINE elements, which are scored as tRNAs by the tRNAscan program [51]. Another alternative is that this expansion is guided by the need to synthesize large quantities of glycine- and glutamic acid-rich proteins, requiring an enrichment of this set of tRNA genes. This strategy is similar to that used by bacteriophages to increase their rate of protein synthesis in a new host, that involves carrying tRNA genes for those codons that are most frequent in bacteriophage genes [52]. Recently it was discovered that HIV also uses a similar strategy, specifically packaging a pool of tRNAs whose codons are poorly represented in the human genome [53]. Interestingly, the addition of glutamic acid to the culture medium increases the production of cellulases by the brown rot fungus Fomitopsis sp. RCK2010 [54], and increases the production of manganese peroxidases and laccases in some white rot fungi [55]. On the other hand, the amino acid glycine is a precursor of heme biosynthesis [56], an important cofactor present in manganese peroxidases and cytochromes, genes that are highly abundant in C. subvermispora.
We also found a strong bias in codon usage, codon pairs, and adaptation to the tRNA pool of C. subvermispora genes involved in lignocellulose degradation. Similar bias has been described in lignin peroxidase from Phaenorochaete chrysosporium [57]. Since a bioinformatic identification of these genes can only be applied to those directly involved in the mineralization/digestion of lignocellulose, such as manganese peroxidase type enzymes, laccases, and cellulases, we devised a functional definition, whereby genes involved in lignocellulose degradation are those induced in the presence of BMA, a model substrate of wood (Group B). Using this functional definition, we were able to detect biases for all evaluated parameters. In turn, the lack of a strong correlation between gene expression levels and adaptation to the tRNA pool or codon bias of C. subvermispora grown in a medium containing glucose or BMA, is consistent with this postulate. A similar correlation between gene function, levels of transcriptional induction, and adaption to the translational machinery is expected when C. subvermispora grown on wood or other media that resemble its natural substrates. Future experiments need to be addressed to test this hypothesis. Studies of phylogenetic reconstruction indicate that the lignocellulose-degrading machinery may have arisen during the Paleozoic [58]. It is reasonable to pose that efficient use of this carbon source requires a metabolic adaptation that involves the whole organism, in addition to the emergence of new types of enzymes. Therefore, it is expected that genes encoding for enzymes that are part of the general metabolism of C. subvermispora also respond to exposure of lignocellulose, as observed in microarray experiments reported by Fernandez Fueyo et al. 2012 [8] were 293 genes modify significantly its expression at least 2 fold (205 from Group B and 88 from Group A). Among them, we can find genes related to the Krebs cycle, cell division, and translation (Table S5).
Expansion of cellulase, laccase, and manganese peroxidase gene families in the C. subvermispora genome permitted that each of these genes developed differential expression; however, the precise contribution of each of these genes in the digestion and mineralization of lignocellulose is not clear. In bacteria, an increase of the copy number of a particular gene is a strategy for increasing its expression. This genomic adaptation explains the expansion of genes directly related to the digestion and mineralization of lignocellulose. However, the complex expression pattern of these genes suggests that this gene expansion process is more intricate than simply a response to the need of producing more enzymes. The high bias in codon usage, adaptation of the tRNA pool, and codon diversity, in some genes directly related to the digestion and mineralization processes of lignocellulose together with the increase in expression of ribosomal proteins in BMA, suggest that C. subvermispora also uses increased translational efficiency as an additional strategy to increase the production these specific set of proteins. Thus, the increase in copy number may be linked to the generation of a diverse array of enzymes that can process a wide range of substrates. In support of this hypothesis, it has been shown that C. subvermipora manganese peroxidases present different kinetic parameters and substrate specificities [8,9].

5. Conclusions

Our results suggest that lignocellulose degradation by C. subvermispora has modified the genome structure of this fungus, changing the bias in codon usage, the tRNA gene pool, and codon diversity in genes that are induced in the presence of wood substrates, in order to optimize the production of these proteins. This strategy may be particularly useful in slow-growing organisms, such as C. subvermispora that cannot increase the production of enzymes by increasing cell mass. To our knowledge, this study is the first example to show metabolic adaptation to a particular ecological niche by modification of the genetic structure of an organism in favor of a selective increase of the translational efficiency of genes involved in metabolizing specific substrates that determine its adaptation to a particular environment.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4425/11/10/1227/s1, Table S1: CPS values of pair-codon present in C. subvermispora, Table S2: tRNA encoded in the Draft genome of C. subvermispora, Table S3: tRNA repeated in C. subvermispora, Table S4: Correlation coefficient among indicators. Table S5: Values of CAI, tAI, AAtAI, Nc, and CPB in C. subvermipora genes expressed in Salt medium with BMA or Glucose. Table S6: Codon bias and translation efficiency of genes encoding for enzymes auxiliary in the lignocellulose degradation.

Author Contributions

Design of the in silico and bioinformatic experiments, M.T.; calculation of tAI and AAtAI values, A.G.; calculation of CAI and CPB values, G.C.; analysis of the relationship between tAI, AAtAI, CAI and CPB with the gene expression, S.L.; analysis the evolutionary implications, D.S.; writing of the article, D.S. and M.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was in part funded by USACH, grant number 021971GM_DAS to M.T., and supported by DI-ULagos to AG.

Acknowledgments

The authors acknowledge Dan Cullen for the critical review of the initial draft of this manuscript. To the authors of Fernandez-Fueyo et al. 2012, and Hori et al. 2014, for making public the expression data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Martínez, A.T.; Speranza, M.; Ruiz-Dueñas, F.J.; Ferreira, P.; Camarero, S.; Guillén, F.; Martínez, M.J.; Gutiérrez, A.; del Río, J.C. Biodegradation of lignocellulosics: Microbial, chemical, and enzymatic aspects of the fungal attack of lignin. Int. Microbiol. 2005, 8, 195–204. [Google Scholar] [PubMed]
  2. Kersten, P.; Cullen, D. Extracellular oxidative systems of the lignin-degrading Basidiomycete Phanerochaete chrysosporium. Fungal Genet. Biol. 2007, 44, 77–87. [Google Scholar] [CrossRef] [PubMed]
  3. Wan, C.; Li, Y. Fungal pretreatment of lignocellulosic biomass. Biotechnol. Adv. 2012, 30, 1447–1457. [Google Scholar] [CrossRef]
  4. Baldrian, P.; Valásková, V. Degradation of cellulose by basidiomycetous fungi. FEMS Microbiol. Rev. 2008, 32, 501–521. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Rüttimann-Johnson, C.; Salas, L.; Vicuña, R.; Kirk, T.K. Extracellular Enzyme Production and Synthetic Lignin Mineralization by Ceriporiopsis subvermispora. Appl. Environ. Microbiol. 1993, 59, 1792–1797. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Lobos, S.; Larraín, J.; Salas, L.; Cullen, D.; Vicuña, R. Isoenzymes of manganese-dependent peroxidase and laccase produced by the lignin-degrading basidiomycete Ceriporiopsis subvermispora. Microbiology 1994, 140, 2691–2698. [Google Scholar] [CrossRef] [Green Version]
  7. Enoki, M.; Watanabe, T.; Nakagame, S.; Koller, K.; Messner, K.; Honda, Y.; Kuwahara, M. Extracellular lipid peroxidation of selective white-rot fungus, Ceriporiopsis subvermispora. FEMS Microbiol. Lett. 1999, 180, 205–211. [Google Scholar] [CrossRef] [PubMed]
  8. Fernandez-Fueyo, E.; Ruiz-Dueñas, F.J.; Ferreira, P.; Floudas, D.; Hibbett, D.S.; Canessa, P.; Larrondo, L.F.; James, T.Y.; Seelenfreund, D.; Lobos, S.; et al. Comparative genomics of Ceriporiopsis subvermispora and Phanerochaete chrysosporium provide insight into selective ligninolysis. Proc. Natl. Acad. Sci. USA 2012, 109, 5458–5463. [Google Scholar] [CrossRef] [Green Version]
  9. Fernández-Fueyo, E.; Ruiz-Dueñas, F.J.; Miki, Y.; Martínez, M.J.; Hammel, K.E.; Martínez, A.T. Lignin-degrading peroxidases from genome of selective ligninolytic fungus Ceriporiopsis subvermispora. J. Biol. Chem. 2012, 287, 16903–16916. [Google Scholar] [CrossRef] [Green Version]
  10. Salas, C.; Lobos, S.; Larraín, J.; Salas, L.; Cullen, D.; Vicuña, R. Properties of laccase isoenzymes produced by the basidiomycete Ceriporiopsis subvermispora. Biotechnol. Appl. Biochem. 1995, 21, 323–333. [Google Scholar]
  11. Tello, M.; Corsini, G.; Larrondo, L.F.; Salas, L.; Lobos, S.; Vicuña, R. Characterization of three new manganese peroxidase genes from the ligninolytic basidiomycete Ceriporiopsis subvermispora. Biochim. Biophys. Acta 2000, 1490, 137–144. [Google Scholar] [CrossRef]
  12. Tello, M.; Seelenfreund, D.; Lobos, S.; Gaskell, J.; Cullen, D.; Vicuña, R. Isolation and characterization of homokaryotic strains from the ligninolytic basidiomycete Ceriporiopsis subvermispora. FEMS Microbiol. Lett. 2001, 199, 91–96. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Manubens, A.; Canessa, P.; Folch, C.; Avila, M.; Salas, L.; Vicuña, R. Manganese affects the production of laccase in the basidiomycete Ceriporiopsis subvermispora. FEMS Microbiol. Lett. 2007, 275, 139–145. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Gutiérrez, M.; Rojas, L.A.; Mancilla-Villalobos, R.; Seelenfreund, D.; Vicuña, R.; Lobos, S. Analysis of manganese-regulated gene expression in the ligninolytic basidiomycete Ceriporiopsis subvermispora. Curr. Genet. 2008, 54, 163–173. [Google Scholar] [CrossRef] [PubMed]
  15. Alvarez, J.M.; Canessa, P.; Mancilla, R.A.; Polanco, R.; Santibáñez, P.A.; Vicuña, R. Expression of genes encoding laccase and manganese-dependent peroxidase in the fungus Ceriporiopsis subvermispora is mediated by an ACE1-like copper-fist transcription factor. Fungal Genet. Biol. 2009, 46, 104–111. [Google Scholar] [CrossRef] [PubMed]
  16. Mancilla, R.A.; Canessa, P.; Manubens, A.; Vicuña, R. Effect of manganese on the secretion of manganese-peroxidase by the basidiomycete Ceriporiopsis subvermispora. Fungal Genet. Biol. 2010, 47, 656–661. [Google Scholar] [CrossRef] [PubMed]
  17. Hori, C.; Gaskell, J.; Igarashi, K.; Kersten, P.; Mozuch, M.; Samejima, M.; Cullen, D. Temporal alterations in the secretome of the selective ligninolytic fungus Ceriporiopsis subvermispora during growth on aspen wood reveal this organism’s strategy for degrading lignocellulose. Appl. Environ. Microbiol. 2014, 80, 2062–2070. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Manubens, A.; Avila, M.; Canessa, P.; Vicuña, R. Differential regulation of genes encoding manganese peroxidase (MnP) in the basidiomycete Ceriporiopsis subvermispora. Curr. Genet. 2003, 43, 433–438. [Google Scholar] [CrossRef]
  19. Sharp, P.M.; Stenico, M.; Peden, J.F.; Lloyd, A.T. Codon usage: Mutational bias, translational selection, or both? Biochem. Soc. Trans. 1993, 21, 835–841. [Google Scholar] [CrossRef] [Green Version]
  20. Supek, F.; Skunca, N.; Repar, J.; Vlahovicek, K.; Smuc, T. Translational selection is ubiquitous in prokaryotes. PLoS Genet. 2010, 6, e1001004. [Google Scholar] [CrossRef] [Green Version]
  21. Tuller, T.; Waldman, Y.Y.; Kupiec, M.; Ruppin, E. Translation efficiency is determined by both codon bias and folding energy. Proc. Natl. Acad. Sci. USA 2010, 107, 3645–3650. [Google Scholar] [CrossRef] [Green Version]
  22. Sharp, P.M.; Li, W.H. The codon Adaptation Index--a measure of directional synonymous codon usage bias, and its potential applications. Nucleic Acids Res. 1987, 15, 1281–1295. [Google Scholar] [CrossRef] [Green Version]
  23. Bulmer, M. Coevolution of codon usage and transfer RNA abundance. Nature 1987, 325, 728–730. [Google Scholar] [CrossRef] [PubMed]
  24. dos Reis, M.; Savva, R.; Wernisch, L. Solving the riddle of codon usage preferences: A test for translational selection. Nucleic Acids Res. 2004, 32, 5036–5044. [Google Scholar] [CrossRef] [Green Version]
  25. Lowe, T.M.; Eddy, S.R. tRNAscan-SE: A program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997, 25, 955–964. [Google Scholar] [CrossRef] [PubMed]
  26. Carver, T.; Bleasby, A. The design of Jemboss: A graphical user interface to EMBOSS. Bioinformatics 2003, 19, 1837–1843. [Google Scholar] [CrossRef]
  27. Tello, M.; Saavedra, J.M.; Spencer, E. Analysis of the use of codon pairs in the HE gene of the ISA virus shows a correlation between bias in HPR codon-pair use and mortality rates caused by the virus. Virol. J. 2013, 10, 180. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Coleman, J.R.; Papamichail, D.; Skiena, S.; Futcher, B.; Wimmer, E.; Mueller, S. Virus attenuation by genome-scale changes in codon pair bias. Science 2008, 320, 1784–1787. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Mueller, S.; Papamichail, D.; Coleman, J.R.; Skiena, S.; Wimmer, E. Reduction of the rate of poliovirus protein synthesis through large-scale codon deoptimization causes attenuation of viral virulence by lowering specific infectivity. J. Virol. 2006, 80, 9687–9696. [Google Scholar] [CrossRef] [Green Version]
  30. Rice, P.; Longden, I.; Bleasby, A. EMBOSS: The European Molecular Biology Open Software Suite. Trends Genet. 2000, 16, 276–277. [Google Scholar] [CrossRef]
  31. Tuller, T.; Carmi, A.; Vestsigian, K.; Navon, S.; Dorfan, Y.; Zaborske, J.; Pan, T.; Dahan, O.; Furman, I.; Pilpel, Y. An evolutionarily conserved mechanism for controlling the efficiency of protein translation. Cell 2010, 141, 344–354. [Google Scholar] [CrossRef] [Green Version]
  32. Thompson, J.D.; Higgins, D.G.; Gibson, T.J. CLUSTAL W: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22, 4673–4680. [Google Scholar] [CrossRef] [Green Version]
  33. Saitou, N.; Nei, M. The neighbor-joining method: A new method for reconstructing phylogenetic trees. Mol. Biol. Evol. 1987, 4, 406–425. [Google Scholar] [PubMed]
  34. Dopazo, J. Estimating errors and confidence intervals for branch lengths in phylogenetic trees by a bootstrap approach. J. Mol. Evol. 1994, 38, 300–304. [Google Scholar] [CrossRef]
  35. Tamura, K.; Nei, M.; Kumar, S. Prospects for inferring very large phylogenies by using the neighbor-joining method. Proc. Natl. Acad. Sci. USA 2004, 101, 11030–11035. [Google Scholar] [CrossRef] [Green Version]
  36. Tamura, K.; Peterson, D.; Peterson, N.; Stecher, G.; Nei, M.; Kumar, S. MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol. Biol. Evol. 2011, 28, 2731–2739. [Google Scholar] [CrossRef] [Green Version]
  37. Fingerman, I.; Nagaraj, V.; Norris, D.; Vershon, A.K. Sfp1 plays a key role in yeast ribosome biogenesis. Eukaryot. Cell 2003, 2, 1061–1068. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Marion, R.M.; Regev, A.; Segal, E.; Barash, Y.; Koller, D.; Friedman, N.; O’Shea, E.K. Sfp1 is a stress- and nutrient-sensitive regulator of ribosomal protein gene expression. Proc. Natl. Acad. Sci. USA 2004, 101, 14315–14322. [Google Scholar] [CrossRef] [Green Version]
  39. Tello, M.; Vergara, F.; Spencer, E. Genomic adaptation of the ISA virus to Salmo salar codon usage. Virol. J. 2013, 10, 223. [Google Scholar] [CrossRef] [Green Version]
  40. Zhao, V.; Jacobs, W.M.; Shakhnovich, E.I. Effect of protein structure on evolution of cotranslational folding. Biophys. J. 2020. [Google Scholar] [CrossRef] [PubMed]
  41. Bitran, A.; Jacobs, W.M.; Zhai, X.; Shakhnovich, E. Cotranslational folding allows misfolding-prone proteins to circumvent deep kinetic traps. Proc. Natl. Acad. Sci. USA 2020, 117, 1485–1495. [Google Scholar] [CrossRef]
  42. Aragonès, L.; Guix, S.; Ribes, E.; Bosch, A.; Pintó, R.M. Fine-tuning translation kinetics selection as the driving force of codon usage bias in the hepatitis A virus capsid. PLoS Pathog. 2010, 6, e1000797. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Plotkin, J.B.; Kudla, G. Synonymous but not the same: The causes and consequences of codon bias. Nat. Rev. Genet. 2011, 12, 32–42. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Sharp, P.M.; Tuohy, T.M.; Mosurski, K.R. Codon usage in yeast: Cluster analysis clearly differentiates highly and lowly expressed genes. Nucleic Acids Res. 1986, 14, 5125–5143. [Google Scholar] [CrossRef] [PubMed]
  45. Sharp, P.M.; Cowe, E.; Higgins, D.G.; Shields, D.C.; Wolfe, K.H.; Wright, F. Codon usage patterns in Escherichia coli, Bacillus subtilis, Saccharomyces cerevisiae, Schizosaccharomyces pombe, Drosophila melanogaster and Homo sapiens; a review of the considerable within-species diversity. Nucleic Acids Res. 1988, 16, 8207–8211. [Google Scholar] [CrossRef] [Green Version]
  46. Bahir, I.; Fromer, M.; Prat, Y.; Linial, M. Viral adaptation to host: A proteome-based analysis of codon usage and amino acid preferences. Mol. Syst. Biol. 2009, 5, 311. [Google Scholar] [CrossRef] [PubMed]
  47. Angov, E. Codon usage: Nature’s roadmap to expression and folding of proteins. Biotechnol. J. 2011, 6, 650–659. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Mueller, S.; Coleman, J.R.; Papamichail, D.; Ward, C.B.; Nimnual, A.; Futcher, B.; Skiena, S.; Wimmer, E. Live attenuated influenza virus vaccines by computer-aided rational design. Nat. Biotechnol. 2010, 28, 723–726. [Google Scholar] [CrossRef] [Green Version]
  49. Zhou, M.; Guo, J.; Cha, J.; Chae, M.; Chen, S.; Barral, J.M.; Sachs, M.S.; Liu, Y. Non-optimal codon usage affects expression, structure and function of clock protein FRQ. Nature 2013, 495, 111–115. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Novoa, E.M.; Pavon-Eternod, M.; Pan, T.; de Pouplana, L.R. A role for tRNA modifications in genome structure and codon usage. Cell 2012, 149, 202–213. [Google Scholar] [CrossRef] [Green Version]
  51. Wenke, T.; Döbel, T.; Sörensen, T.R.; Junghans, H.; Weisshaar, B.; Schmidt, T. Targeted identification of short interspersed nuclear element families shows their widespread existence and extreme heterogeneity in plant genomes. Plant Cell 2011, 23, 3117–3128. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Enav, H.; Béjà, O.; Mandel-Gutfreund, Y. Cyanophage tRNAs may have a role in cross-infectivity of oceanic Prochlorococcus and Synechococcus hosts. ISME J. 2012, 6, 619–628. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. van Weringh, A.; Ragonnet-Cronin, M.; Pranckeviciene, E.; Pavon-Eternod, M.; Kleiman, L.; Xia, X. HIV-1 modulates the tRNA pool to improve translation efficiency. Mol. Biol. Evol. 2011, 28, 1827–1834. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Deswal, D.; Khasa, Y.P.; Kuhad, R.C. Optimization of cellulase production by a brown rot fungus Fomitopsis sp. RCK2010 under solid state fermentation. Bioresour. Technol. 2011, 102, 6065–6072. [Google Scholar] [CrossRef] [PubMed]
  55. Levin, L.; Melignani, E.; Ramos, A.M. Effect of nitrogen sources and vitamins on ligninolytic enzyme production by some white-rot fungi. Dye decolorization by selected culture filtrates. Bioresour. Technol. 2010, 101, 4554–4563. [Google Scholar] [CrossRef]
  56. Heinemann, I.U.; Jahn, M.; Jahn, D. The biochemistry of heme biosynthesis. Arch. Biochem. Biophys. 2008, 474, 238–251. [Google Scholar] [CrossRef]
  57. Ritch, T.G.; Gold, M.H. Characterization of a highly expressed lignin peroxidase-encoding gene from the basidiomycete Phanerochaete chrysosporium. Gene 1992, 118, 73–80. [Google Scholar] [CrossRef]
  58. Floudas, D.; Binder, M.; Riley, R.; Barry, K.; Blanchette, R.A.; Henrissat, B.; Martínez, A.T.; Otillar, R.; Spatafora, J.W.; Yadav, J.S.; et al. The Paleozoic origin of enzymatic lignin decomposition reconstructed from 31 fungal genomes. Science 2012, 336, 1715–1719. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Phylogenetic analysis of tRNAs from Ceriporiopsis subvermispora. tRNAs were numbered according the tRNAScan-SE output. Prediction of amino acid charge, sequence of anticodon, and presence (+) or absence (−) of introns are indicated. Each tRNA type is indicated with a different symbol and color.
Figure 1. Phylogenetic analysis of tRNAs from Ceriporiopsis subvermispora. tRNAs were numbered according the tRNAScan-SE output. Prediction of amino acid charge, sequence of anticodon, and presence (+) or absence (−) of introns are indicated. Each tRNA type is indicated with a different symbol and color.
Genes 11 01227 g001
Figure 2. Correlation between codon usage and relative adaptiveness to the tRNA pool in C. subvermispora. Correlation between the relative adaptiveness to the tRNA pool and codon frequency usage (per thousand) of the 61 codons is shown.
Figure 2. Correlation between codon usage and relative adaptiveness to the tRNA pool in C. subvermispora. Correlation between the relative adaptiveness to the tRNA pool and codon frequency usage (per thousand) of the 61 codons is shown.
Genes 11 01227 g002
Figure 3. Codon bias and adaptation to tRNA pool of genes induced by Ball-Milled Aspen (BMA). The figure shows the Codon Adaptation Index (CAI), AAtAI, codon pair bias (CPB), effective Number of Codons (Nc), and tAI values of C. subvermispora genes. The genes were classified according their change in expression in BMA vs. glucose as carbon source in salt medium, in repressed (p < 0.05) at least twofold (Group A), induced (p < 0.05) at least twofold (Group B), without statistical changes (p > 0.05) (Group C), and genes with slightly (lower than two fold) but significant changes in expression (Group D). Statistical differences between groups were assessed using Mann-Whitney test. b corresponds a p < 0.001 in a comparison with Group B.
Figure 3. Codon bias and adaptation to tRNA pool of genes induced by Ball-Milled Aspen (BMA). The figure shows the Codon Adaptation Index (CAI), AAtAI, codon pair bias (CPB), effective Number of Codons (Nc), and tAI values of C. subvermispora genes. The genes were classified according their change in expression in BMA vs. glucose as carbon source in salt medium, in repressed (p < 0.05) at least twofold (Group A), induced (p < 0.05) at least twofold (Group B), without statistical changes (p > 0.05) (Group C), and genes with slightly (lower than two fold) but significant changes in expression (Group D). Statistical differences between groups were assessed using Mann-Whitney test. b corresponds a p < 0.001 in a comparison with Group B.
Genes 11 01227 g003
Table 1. Statistics of tRNAs present in the Ceriporiopsis subvermispora genome.
Table 1. Statistics of tRNAs present in the Ceriporiopsis subvermispora genome.
tRNA TypeNumber of tRNAsAnticodon and Frequency *ScaffoldsNumber of tRNAs with Introns
Ala14AGC7, CGC5, TGC23, 5, 7, 12, 14, 18, 1913
Arg16ACG8, CCG2, CCT2, TCG3, TCT12, 3, 5, 9, 14, 19, 31, 3915
Asn5GTT54, 5, 7, 13, 145
Asp10GTC101, 6, 8, 12, 16, 19, 28, 417
Cys3GCA310, 13, 201
Gln8CTG6, TTG21, 2, 6, 17, 18, 288
Glu15CTC12, TTC33, 8, 9, 19, 285
Gly17CCC2, GCC12, TCC31, 2, 5, 7, 9, 10, 13, 201
His5GTG52, 6, 7, 125
Ile9AAT8, TAT14, 5, 21, 339
Leu15AAG7, CAA2, CAG4, TAA1, TAG11, 2, 4, 5, 6, 7, 9, 13, 17, 25, 27, 3014
Lys11CTT10, TTT11, 311
Met8CAT81, 2, 3, 9, 10, 11, 12, 237
Phe5GAA51, 5, 9, 104
Pro10AGG5, CGG4, TGG14, 6, 7, 14, 18, 22, 3510
Ser13ACT1, AGA4, CGA3, GCT3, TGA21, 3, 4, 5, 6, 9, 15, 188
Thr9AGT6, CGT2, TGT11, 3, 4, 8, 9, 10, 226
Trp3CCA319, 20, 303
Tyr4GTA42, 3, 8, 304
Val11TAC1, CAC2, AAC6, CAC1, AAC15, 6, 12, 14, 15, 28, 322
* Frequency is indicated as superscript together the anticodon
Table 2. Correlation coefficient between gene expression, codon bias, and translational efficiency.
Table 2. Correlation coefficient between gene expression, codon bias, and translational efficiency.
BMA (ρ)Glu (ρ)BMA/Glu (ρ)BMA/Glu (ρ) aBMA/Glu (ρ) bBMA/Glu (ρ) c
CAI−4.95 × 10−3−8.04 × 10−2 ***3.39 × 10−1 ***4.20 × 10−1 **4.66 × 10−1 ***3.09 × 10−1 ***
CPB6.09 × 10−2 ***6.59 × 10−3 ***2.61 × 10−1 ***NA3.26 × 10−1 ***2.47 × 10−1 ***
Nc1.81 × 10−2 *8.17 × 10−2 ***−2.63 × 10−1 ***−4.51 × 10−1 ***−4.42 × 10−1 ***−2.21 × 10−1 ***
tAI3.48 × 10−2 ***−3.52 × 10−2 ***3.17 × 10−1 ***2.84 × 10−1 *4.36 × 10−1 ***2.82 × 10−1 ***
AAtAI6.54 × 10−2 ***8.34 × 10−42.67 × 10−1 ***4.75 × 10−1 ***4.17 × 10−1 ***2.27 × 10−1 ***
* p < 0.05 ** p < 0.01 *** p < 0.001. a Genes with p-value lower than 0.001, n = 52. b Genes with p-value between 0.001 and 0.05, n = 1572. c Genes with p-value higher than 0.05, n = 10,471. Glu = Glucose.
Table 3. Adaptation to codon usage, tRNA pool, and codon bias of genes from C. subvermispora involved in mineralization and digestion of lignocellulose. Bold letters indicate the five highest values.
Table 3. Adaptation to codon usage, tRNA pool, and codon bias of genes from C. subvermispora involved in mineralization and digestion of lignocellulose. Bold letters indicate the five highest values.
Ligno-Cellulolytic FunctionTranscript IDCAIAAtAItAICPBNcZ-CAIZ-AAtAIZ-CPBZ-tAIPutative FunctionMicroarray Signal (log2)
GlucoseBMABMA/Glup-Value
Fungal lignin peroxidase498630.8140.7200.3710.02840.951.2361.8480.1361.513Peroxidase, MnP9.5811.152.9700.08320
502970.8220.7350.3840.05541.461.3942.3211.5652.111Peroxidase, MnP13.1311.880.4210.19000
506860.8230.7280.3790.05141.501.4142.1001.3871.901Peroxidase, MnP9.359.421.0480.45600
1063800.8150.7160.3710.03938.671.2551.7220.7051.536Peroxidase, MnP8.938.931.0000.99900
1113640.7700.6890.3570.01550.570.3640.872−0.6110.929Peroxidase, VP-like9.189.110.9570.52800
1175210.8190.7220.3670.04240.881.3351.9110.8761.340Peroxidase, MnP9.579.490.9450.30600
1241440.6960.6530.3210.03161.00−1.101−0.2610.284−0.660Peroxidase, generic11.3011.210.9420.29300
1260180.8280.7210.3660.03040.551.5131.8800.2301.314Peroxidase, MnP 9.699.500.8740.10000
1260580.7880.6730.3440.02644.680.7210.368−0.0150.340Peroxidase, MnP 9.949.440.7070.00595
1285900.8240.7240.3670.05241.121.4341.9741.4091.352Peroxidase, MnP 9.789.520.8370.01960
1294180.7670.7090.3730.04449.060.3051.5021.0111.619Peroxidase, MnP 10.2212.926.5080.00895
1306590.7810.7210.3680.02250.530.5821.880−0.2201.396Peroxidase, LiP-like10.8810.300.6730.03900
1360580.7700.6880.3520.02950.910.3640.8410.1670.683Peroxidase, MnP 10.1610.761.5190.16700
1519470.8310.7290.3730.05240.421.5722.1321.4101.600Peroxidase, MnP 9.008.930.9550.24200
1553720.7970.7110.3690.04042.880.8991.5650.7651.443Peroxidase, MnP 8.868.840.9880.74100
1699680.8160.7180.3690.03242.961.2751.7850.3591.434Peroxidase, MnP 10.2710.160.9300.44200
Laccase841700.7410.6640.3250.03258.64−0.2100.0850.338−0.517laccase10.4510.340.9310.17000
880890.8410.6840.3360.03941.001.7700.7150.695−0.016Laccase9.389.290.9390.17800
1208340.7370.6690.3300.01157.80−0.2890.242−0.825−0.271Laccase10.7410.660.9450.50800
1270450.7520.6810.3390.02454.320.0080.620−0.1030.112Laccase10.5410.020.6940.00260
1270500.7210.6670.3340.01560.26−0.6060.179−0.581−0.098Laccase11.3511.190.8950.54300
1307830.7910.7050.3470.03743.130.7801.3760.6050.461Laccase11.0213.776.7660.00426
1496680.7750.6790.3340.01951.250.4630.557−0.397−0.120Laccase9.919.810.9310.35300
Cellulose Binding Protein597330.8120.7040.3460.02646.111.1961.344−0.0180.414GH10-CBM19.2413.6921.7830.00004
666880.8240.7160.3470.02547.021.4341.722−0.0480.488GH61-CBM19.6214.8737.9010.00002
675610.8530.7300.3660.02039.432.0082.163−0.3361.313GH10-CBM110.3013.7110.6420.00002
685690.8370.7350.3450.02443.591.6912.321−0.0950.383CE1-CBM19.8013.109.8580.00023
795570.8020.7150.3450.02242.150.9981.691−0.2310.373GH5-CBM110.2414.0413.8740.00001
875800.7830.6930.3210.01049.310.6220.998−0.884−0.692CE16-CBM110.9414.4611.5400.00004
895330.8380.7360.3530.05242.061.7112.3521.4260.719GH61-CBM110.3513.569.3060.00002
895340.8700.7310.3520.04437.582.3452.1950.9800.693GH61-CBM19.4210.131.6370.02300
1019250.8450.7230.3590.03338.221.8501.9430.3941.020GH7-CBM18.848.961.0860.09060
1067770.8040.7120.3480.03549.871.0381.5960.5110.500GH5- CBM19.5414.3528.0500.00001
1098400.8570.7410.3700.03938.162.0872.5090.7331.471GH10- CBM19.5211.855.0090.00053
1290280.8520.7330.3630.04239.911.9882.2580.8491.181GH5- CBM19.8110.791.9850.01830
1338090.8650.7130.3380.04037.082.2461.6280.7790.091GH11-CBM110.8611.421.4660.02090
1485880.8560.7420.3760.00834.452.0672.541−0.9671.756GH7-CBM111.0212.623.0150.00018
CDH847920.8030.6880.3380.02348.091.0180.841−0.1370.071CDH9.2913.7622.2410.00002
871100.7690.6790.3320.02450.710.3450.557−0.111−0.189 11.3111.040.8270.20200
1256100.7620.6650.3180.04453.340.2060.1160.981−0.814cir1 CBM110.1710.401.1700.23300
1475440.7120.6650.3180.02556.98−0.7840.116−0.060−0.798 11.2411.020.8600.19900
Delta 12 Dehidrogenase588800.7270.6700.3170.00859.40−0.4870.274−0.973−0.853Δ-12 FAD10.3610.290.9560.72900
1210740.7310.6440.3130.02255.49−0.408−0.545−0.219−1.033Δ-12 FAD10.5810.230.7830.01490
1240500.8560.7360.3680.04238.222.0672.3520.8751.416Δ-12 FAD. Cs-fad212.7412.660.9410.65300
1361010.7140.6490.3100.00756.78−0.745−0.387−1.016−1.178Δ-12 FAD11.0112.542.8950.00022
1676900.7360.6700.3250.01358.45−0.3090.274−0.706−0.484Δ-12 FAD10.6710.110.6780.00870
Delta 9 Dehidrogenase878750.8100.6940.3430.03544.181.1561.0300.4750.313Δ-9 FAD, Cs-ole18.938.941.0060.88000
1290450.7280.6790.3350.04556.92−0.4670.5571.017−0.075Δ-9 FAD, Cs-ole18.958.910.9740.52700
1290480.8480.7400.3660.05041.461.9092.4781.2911.292Δ-9 FAD, Cs-ole111.7812.351.4780.02980
1336750.7600.6690.3300.02454.620.1660.242−0.113−0.266Δ-9 FAD, Cs-ole19.649.510.9170.17900
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gonzalez, A.; Corsini, G.; Lobos, S.; Seelenfreund, D.; Tello, M. Metabolic Specialization and Codon Preference of Lignocellulolytic Genes in the White Rot Basidiomycete Ceriporiopsis subvermispora. Genes 2020, 11, 1227. https://doi.org/10.3390/genes11101227

AMA Style

Gonzalez A, Corsini G, Lobos S, Seelenfreund D, Tello M. Metabolic Specialization and Codon Preference of Lignocellulolytic Genes in the White Rot Basidiomycete Ceriporiopsis subvermispora. Genes. 2020; 11(10):1227. https://doi.org/10.3390/genes11101227

Chicago/Turabian Style

Gonzalez, Alex, Gino Corsini, Sergio Lobos, Daniela Seelenfreund, and Mario Tello. 2020. "Metabolic Specialization and Codon Preference of Lignocellulolytic Genes in the White Rot Basidiomycete Ceriporiopsis subvermispora" Genes 11, no. 10: 1227. https://doi.org/10.3390/genes11101227

APA Style

Gonzalez, A., Corsini, G., Lobos, S., Seelenfreund, D., & Tello, M. (2020). Metabolic Specialization and Codon Preference of Lignocellulolytic Genes in the White Rot Basidiomycete Ceriporiopsis subvermispora. Genes, 11(10), 1227. https://doi.org/10.3390/genes11101227

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop