RESUMEN
Nucleotide insertions/deletions are ubiquitous in eukaryotic genomes, and the resulting hemizygous (unpaired) DNA has significant, heritable effects on adjacent DNA. However, little is known about the genetic behavior of insertion DNA. Here, we describe a binary transgenic system to study the behavior of insertion DNA during meiosis. Transgenic Arabidopsis lines were generated to carry two different defective reporter genes on nonhomologous chromosomes, designated as "recipient" and "donor" lines. Double hemizygous plants (harboring unpaired DNA) were produced by crossing between the recipient and the donor, and double homozygous lines (harboring paired DNA) via self-pollination. The transfer of the donor's unmutated sequence to the recipient generated a functional ß-glucuronidase gene, which could be visualized by histochemical staining and corroborated by polymerase chain reaction amplification and sequencing. More than 673 million seedlings were screened, and the results showed that meiotic ectopic recombination in the hemizygous lines occurred at a frequency >6.49-fold higher than that in the homozygous lines. Gene conversion might have been exclusively or predominantly responsible for the gene correction events. The direct measurement of ectopic recombination events provided evidence that an insertion, in the absence of an allelic counterpart, could scan the entire genome for homologous counterparts with which to pair. Furthermore, the unpaired (hemizygous) architectures could accelerate ectopic recombination between itself and interchromosomal counterparts. We suggest that the ectopic recombination accelerated by hemizygous architectures may be a general mechanism for interchromosomal recombination through ubiquitously dispersed repeat sequences in plants, ultimately contributing to genetic renovation and eukaryotic evolution.
Asunto(s)
Arabidopsis/genética , ADN de Plantas/genética , Mutagénesis Insercional , Arabidopsis/metabolismo , Cromosomas , Intercambio Genético , Conversión Génica , Hemicigoto , Recombinación Homóloga , Homocigoto , Meiosis/genética , Plantas Modificadas Genéticamente , Recombinación Genética , PlantonesRESUMEN
Nucleotide-binding site-leucine-rich repeat (NBS-LRR) genes make up the largest plant disease resistance gene family (R genes), with hundreds of copies occurring in individual angiosperm genomes. However, the expansion history of NBS-LRR genes during angiosperm evolution is largely unknown. By identifying more than 6,000 NBS-LRR genes in 22 representative angiosperms and reconstructing their phylogenies, we present a potential framework of NBS-LRR gene evolution in the angiosperm. Three anciently diverged NBS-LRR classes (TNLs, CNLs, and RNLs) were distinguished with unique exon-intron structures and DNA motif sequences. A total of seven ancient TNL, 14 CNL, and two RNL lineages were discovered in the ancestral angiosperm, from which all current NBS-LRR gene repertoires were evolved. A pattern of gradual expansion during the first 100 million years of evolution of the angiosperm clade was observed for CNLs. TNL numbers remained stable during this period but were eventually deleted in three divergent angiosperm lineages. We inferred that an intense expansion of both TNL and CNL genes started from the Cretaceous-Paleogene boundary. Because dramatic environmental changes and an explosion in fungal diversity occurred during this period, the observed expansions of R genes probably reflect convergent adaptive responses of various angiosperm families. An ancient whole-genome duplication event that occurred in an angiosperm ancestor resulted in two RNL lineages, which were conservatively evolved and acted as scaffold proteins for defense signal transduction. Overall, the reconstructed framework of angiosperm NBS-LRR gene evolution in this study may serve as a fundamental reference for better understanding angiosperm NBS-LRR genes.
Asunto(s)
Evolución Molecular , Genes de Plantas , Variación Genética , Magnoliopsida/genética , Proteínas NLR/genética , Arabidopsis/genética , Secuencia de Bases , Sitios de Unión , Exones/genética , Intrones/genética , Motivos de Nucleótidos/genética , Filogenia , Especificidad de la EspecieRESUMEN
Very long-chain fatty acids (VLCFAs) play an important role in the survival and development of plants, and VLCFA synthesis is regulated by ß-ketoacyl-CoA synthases (KCSs), which catalyze the condensation of an acyl-CoA with malonyl-CoA. Here, we present a genome-wide survey of the genes encoding these enzymes, KCS genes, in 28 species (26 genomes and two transcriptomes), which represents a large phylogenetic scale, and also reconstruct the evolutionary history of this gene family. KCS genes were initially single-copy genes in the green plant lineage; duplication resulted in five ancestral copies in land plants, forming five fundamental monophyletic groups in the phylogenetic tree. Subsequently, KCS genes duplicated to generate 11 genes of angiosperm origin, expanding up to 20-30 members in further-diverged angiosperm species. During this process, tandem duplications had only a small contribution, whereas polyploidy events and large-scale segmental duplications appear to be the main driving force. Accompanying this expansion were variations that led to the sub- and neofunctionalization of different members, resulting in specificity that is likely determined by the 3-D protein structure. Novel functions involved in other physiological processes emerged as well, though redundancy is also observed, largely among recent duplications. Conserved sites and variable sites of KCS proteins are also identified by statistical analysis. The variable sites are likely to be involved in the emergence of product specificity and catalytic power, and conserved sites are possibly responsible for the preservation of fundamental function.
Asunto(s)
3-Oxoacil-(Proteína Transportadora de Acil) Sintasa/genética , Evolución Molecular , Familia de Multigenes/genética , Filogenia , Secuencia de Aminoácidos , Duplicación de Gen , Genoma de Planta , Magnoliopsida/genética , Poliploidía , Transcriptoma/genéticaRESUMEN
Plant genomes harbor dozens to hundreds of nucleotide-binding site-leucine-rich repeat (NBS-LRR) genes; however, the long-term evolutionary history of these resistance genes has not been fully understood. This study focuses on five Brassicaceae genomes and the Carica papaya genome to explore changes in NBS-LRR genes that have taken place in this Rosid II lineage during the past 72 million years. Various numbers of NBS-LRR genes were identified from Arabidopsis lyrata (198), A. thaliana (165), Brassica rapa (204), Capsella rubella (127), Thellungiella salsuginea (88), and C. papaya (51). In each genome, the identified NBS-LRR genes were found to be unevenly distributed among chromosomes and most of them were clustered together. Phylogenetic analysis revealed that, before and after Brassicaceae speciation events, both toll/interleukin-1 receptor-NBS-LRR (TNL) genes and non-toll/interleukin-1 receptor-NBS-LRR (nTNL) genes exhibited a pattern of first expansion and then contraction, suggesting that both subclasses of NBS-LRR genes were responding to pathogen pressures synchronically. Further, by examining the gain/loss of TNL and nTNL genes at different evolutionary nodes, this study revealed that both events often occurred more drastically in TNL genes. Finally, the phylogeny of nTNL genes suggested that this NBS-LRR subclass is composed of two separate ancient gene types: RPW8-NBS-LRR and Coiled-coil-NBS-LRR.
Asunto(s)
Brassicaceae/genética , Evolución Molecular , Genes de Plantas , Nucleótidos/metabolismo , Proteínas de Plantas/genética , Proteínas/genética , Sitios de Unión , Mapeo Cromosómico , Cromosomas de las Plantas/genética , Duplicación de Gen , Sitios Genéticos , Proteínas Repetidas Ricas en Leucina , Funciones de Verosimilitud , Familia de Multigenes , Filogenia , Especificidad de la Especie , Sintenía/genéticaRESUMEN
Proper utilization of plant disease resistance genes requires a good understanding of their short- and long-term evolution. Here we present a comprehensive study of the long-term evolutionary history of nucleotide-binding site (NBS)-leucine-rich repeat (LRR) genes within and beyond the legume family. The small group of NBS-LRR genes with an amino-terminal RESISTANCE TO POWDERY MILDEW8 (RPW8)-like domain (referred to as RNL) was first revealed as a basal clade sister to both coiled-coil-NBS-LRR (CNL) and Toll/Interleukin1 receptor-NBS-LRR (TNL) clades. Using Arabidopsis (Arabidopsis thaliana) as an outgroup, this study explicitly recovered 31 ancestral NBS lineages (two RNL, 21 CNL, and eight TNL) that had existed in the rosid common ancestor and 119 ancestral lineages (nine RNL, 55 CNL, and 55 TNL) that had diverged in the legume common ancestor. It was shown that, during their evolution in the past 54 million years, approximately 94% (112 of 119) of the ancestral legume NBS lineages experienced deletions or significant expansions, while seven original lineages were maintained in a conservative manner. The NBS gene duplication pattern was further examined. The local tandem duplications dominated NBS gene gains in the total number of genes (more than 75%), which was not surprising. However, it was interesting from our study that ectopic duplications had created many novel NBS gene loci in individual legume genomes, which occurred at a significant frequency of 8% to 20% in different legume lineages. Finally, by surveying the legume microRNAs that can potentially regulate NBS genes, we found that the microRNA-NBS gene interaction also exhibited a gain-and-loss pattern during the legume evolution.
Asunto(s)
Evolución Molecular , Fabaceae/genética , Familia de Multigenes , Eliminación de Gen , Duplicación de Gen , Genoma de Planta , MicroARNs/metabolismo , FilogeniaRESUMEN
To offer the reference and method for salt damage in the cultivation of Marsdenia tenacissima, the seeds of M. tenacissima collected from Maguan city ( Yunnan province) were taken as the test materials to study the effects of different priming materials on improving germination and growth under high-level salt stress condition. Four different treatments, which were GA3, KNO3-KH2PO4, PEG-6000, NaCl, combined with ANOVA were applied to test the performance of germination energy, germination percentage, germination index, MDA, SOD, and CAT. The results showed that the seed germination was obviously inhibited under salt stress and the soaked seeds with different priming materials could alleviate the damage of salt stress. Under these treatments, the activities of SOD, CAT the content of soluble protein significantly increased. While the content of MDA significantly decreased. The maximum index was obtained when treated with 1.20% KNO3-KH2PO4, the germination percentage increased from 52.67% to 87.33% and the activity of SOD increased from 138.01 to 219.44 respectively. Comparing with the treatment of 1.20% KNO3-KH2PO4, the germination percentage of treating with 300 mg x L(-1) GA3 increased from 52.67% to 80.67%, while the activity of SOD increased from 138.01 to 444.61.
Asunto(s)
Germinación/fisiología , Marsdenia/crecimiento & desarrollo , Cloruro de Sodio/farmacología , Germinación/efectos de los fármacos , Marsdenia/efectos de los fármacos , Nitratos/farmacología , Polietilenglicoles/farmacología , Compuestos de Potasio/farmacología , Semillas/efectos de los fármacos , Semillas/crecimiento & desarrollo , Estrés Fisiológico , Xantonas/farmacologíaRESUMEN
Dioscorea rotundata is an important food crop that is mainly cultivated in subtropical regions of the world. D. rotundata is frequently infected by various pathogens during its lifespan, which results in a substantial economic loss in terms of yield and quality. The disease resistance gene (R gene) profile of D. rotundata is largely unknown, which has greatly hampered molecular study of disease resistance in this species. Nucleotide-binding site-leucine-rich repeat (NBS-LRR) genes are the largest group of plant R genes, and they play important roles in plant defense responses to various pathogens. In this study, 167 NBS-LRR genes were identified from the D. rotundata genome. Subsequently, one gene was assigned to the resistance to powdery mildew8 (RPW8)-NBS-LRR (RNL) subclass and the other 166 genes to the coiled coil (CC)-NBS-LRR (CNL) subclass. None of the Toll/interleukin-1 receptor (TIR)-NBS-LRR (TNL) genes were detected in the genome. Among them, 124 genes are located in 25 multigene clusters and 43 genes are singletons. Tandem duplication serves as the major force for the cluster arrangement of NBS-LRR genes. Segmental duplication was detected for 18 NBS-LRR genes, although no whole-genome duplication has been documented for the species. Phylogenetic analysis revealed that D. rotundata NBS-LRR genes share 15 ancestral lineages with Arabidopsis thaliana genes. The NBS-LRR gene number increased by more than a factor of 10 during D. rotundata evolution. A conservatively evolved ancestral lineage was identified from D. rotundata, which is orthologs to the Arabidopsis RPM1 gene. Transcriptome analysis for four different tissues of D. rotundata revealed a low expression of most NBS-LRR genes, with the tuber and leaf displaying a relatively high NBS-LRR gene expression than the stem and flower. Overall, this study provides a complete set of NBS-LRR genes for D. rotundata, which may serve as a fundamental resource for mining functional NBS-LRR genes against various pathogens.
RESUMEN
OBJECTIVE: To study the level of NPK balanced fertilization on growth of Paeonia lactiflora in different growing periods. METHODS: This experiment was designed as orthogonal test of three factors and three levels. Fresh weight of root, number of bud, number of root division, length of the longest root and diameter of widest root were used as indicators. RESULTS: Nitrogenous, phoshate and potash fertilizer affected the growth of the Paeonia lactiflora from the first to the third year. And the nitrogenous was the main factor during the fourth year. CONCLUSION: NPK fertilizer, especially N and P fertilizer are needed during the first year of the Paeonia lactiflora. Only the nitrogenous fertilizer is needed in the fourth year.
Asunto(s)
Fertilizantes , Paeonia/crecimiento & desarrollo , Plantas Medicinales/crecimiento & desarrollo , Biomasa , Nitrógeno , Fósforo , Raíces de Plantas/crecimiento & desarrollo , Potasio , Factores de TiempoRESUMEN
Plant resistance conferred by nucleotide binding site (NBS)-encoding resistance genes plays a key role in the defense against various pathogens throughout the entire plant life cycle. However, comparative analyses for the systematic evaluation and determination of the evolutionary modes of NBS-encoding genes among Solanaceae species are rare. In this study, 447, 255, and 306 NBS-encoding genes were identified from the genomes of potato, tomato, and pepper, respectively. These genes usually clustered as tandem arrays on chromosomes; few existed as singletons. Phylogenetic analysis indicated that three subclasses [TNLs (TIR-NBS-LRR), CNLs (CC-NBS-LRR), and RNLs (RPW8-NBS-LRR)] each formed a monophyletic clade and were distinguished by unique exon/intron structures and amino acid motif sequences. By comparing phylogenetic and systematic relationships, we inferred that the NBS-encoding genes in the present genomes of potato, tomato, and pepper were derived from 150 CNL, 22 TNL, and 4 RNL ancestral genes, and underwent independent gene loss and duplication events after speciation. The NBS-encoding genes therefore exhibit diverse and dynamic evolutionary patterns in the three Solanaceae species, giving rise to the discrepant gene numbers observed today. Potato shows a "consistent expansion" pattern, tomato exhibits a pattern of "first expansion and then contraction," and pepper presents a "shrinking" pattern. The earlier expansion of CNLs in the common ancestor led to the dominance of this subclass in gene numbers. However, RNLs remained at low copy numbers due to their specific functions. Along the evolutionary process of NBS-encoding genes in Solanaceae, species-specific tandem duplications contributed the most to gene expansions.
Asunto(s)
Evolución Molecular , Mutación con Ganancia de Función , Modelos Genéticos , Proteínas de Plantas/genética , Solanaceae/genética , Sitios de Unión , Especiación Genética , Mutación con Pérdida de Función , Nucleótidos/metabolismo , Filogenia , Proteínas de Plantas/química , Proteínas de Plantas/metabolismo , Unión Proteica , Solanaceae/clasificaciónRESUMEN
The Toll-interleukin-1 receptor (TIR) and Nucleotide-binding site (NBS) domains are two major components of the TIR-NBS-leucine-rich repeat family plant disease resistance genes. Extensive functional and evolutionary studies have been performed on these genes; however, the characterization of a small group of genes that are composed of atypical TIR and NBS domains, namely XTNX genes, is limited. The present study investigated this specific gene family by conducting genome-wide analyses of 59 green plant genomes. A total of 143 XTNX genes were identified in 51 of the 52 land plant genomes, whereas no XTNX gene was detected in any green algae genomes, which indicated that XTNX genes originated upon emergence of land plants. Phylogenetic analysis revealed that the ancestral XTNX gene underwent two rounds of ancient duplications in land plants, which resulted in the formation of clades I/II and clades IIa/IIb successively. Although clades I and IIb have evolved conservatively in angiosperms, the motif composition difference and sequence divergence at the amino acid level suggest that functional divergence may have occurred since the separation of the two clades. In contrast, several features of the clade IIa genes, including the absence in the majority of dicots, the long branches in the tree, the frequent loss of ancestral motifs, and the loss of expression in all detected tissues of Zea mays, all suggest that the genes in this lineage might have undergone pseudogenization. This study highlights that XTNX genes are a gene family originated anciently in land plants and underwent specific conservative pattern in evolution.
RESUMEN
OBJECTIVE: To study the ITS sequence variation of Pinellia ternata from different population in China, and it correlation to geographical distribution and morpha of the plant. METHOD: The rDNA ITS regions of various P. ternata were amplified and sequenced. And they were analyzed by means of the software of CLUSTRAL and MEGA. RESULT: Complete sequence of ITS and 5.8S rDNA of 16 different P. ternata population were obtained. The sequences of ITS1, 5.8S and ITS2 are 276,162 and 246 bp, respectively. ITS1 was more conservative than ITS2. Phylogenetic tree based on ITS1 and ITS2 sequences data was conducted by Neighbor-joining method. CONCLUSION: Ribosomal DNA ITS sequence analyses can be applied to the resource research of P. ternata.
Asunto(s)
ADN de Plantas/genética , ADN Espaciador Ribosómico/genética , Pinellia/genética , Plantas Medicinales/genética , Secuencia de Bases , China , Ecosistema , Datos de Secuencia Molecular , Filogenia , Tubérculos de la Planta/genética , ARN Ribosómico 5.8S/genética , Análisis de Secuencia de ADNRESUMEN
A majority of land plants can form symbiosis with arbuscular mycorrhizal (AM) fungi. MicroRNAs (miRNAs) have been implicated to regulate this process in legumes, but their involvement in non-legume species is largely unknown. In this study, by performing deep sequencing of sRNA libraries in tomato roots and comparing with tomato genome, a total of 700 potential miRNAs were predicted, among them, 187 are known plant miRNAs that have been previously deposited in miRBase. Unlike the profiles in other plants such as rice and Arabidopsis, a large proportion of predicted tomato miRNAs was 24 nt in length. A similar pattern was observed in the potato genome but not in tobacco, indicating a Solanum genus-specific expansion of 24-nt miRNAs. About 40% identified tomato miRNAs showed significantly altered expressions upon Rhizophagus irregularis inoculation, suggesting the potential roles of these novel miRNAs in AM symbiosis. The differential expression of five known and six novel miRNAs were further validated using qPCR analysis. Interestingly, three up-regulated known tomato miRNAs belong to a known miR171 family, a member of which has been reported in Medicago truncatula to regulate AM symbiosis. Thus, the miR171 family likely regulates AM symbiosis conservatively across different plant lineages. More than 1000 genes targeted by potential AM-responsive miRNAs were provided and their roles in AM symbiosis are worth further exploring.
RESUMEN
A major soybean (Forrest cultivar) quantitative trait locus (QTL) gene, Rhg4, which controls resistance to soybean cyst nematodes (SCN), encodes the enzyme serine hydroxylmethyltransferase (SHMT). The resistant allele possesses two critical missense mutations (P130R and N358Y) compared to that of the sensitive allele, rhg4. To understand the evolutionary history of this gene, sequences of 117 SHMT family members from 18 representative plant species were used to reconstruct their phylogeny. According to this phylogeny, the plant SHMT gene family can be divided into two groups and four subgroups (Ia, Ib, IIa, and IIb). Belonging to the Subgroup Ia lineage, the rhg4 gene evolved from a recent duplication event in Glycine sp.. To further explore how the SCN-resistant allele emerged, both the rhg4 gene and its closest homolog, the rhg4h gene, were isolated from 33 cultivated and 68 wild soybean varieties. The results suggested that after gene duplication, the soybean rhg4 gene accumulated a higher number of non-synonymous mutations than rhg4h. Although a higher number of segregating sites and gene haplotypes were detected in wild soybeans than in cultivars, the SCN-resistant Rhg4 allele (represented by haplotype 4) was not found in wild varieties. Instead, a very similar allele, haplotype 3, was observed in wild soybeans at a frequency of 7.4%, although it lacked the two critical non-synonymous substitutions. Taken together, these findings support that the SCN-resistant Rhg4 allele likely emerged via artificial selection during the soybean domestication process, based on a SCN-sensitive allele inherited from wild soybeans.
RESUMEN
A genome triplication took place in the ancestor of Brassiceae species after the split of the Arabidopsis lineage. The postfragmentation and shuffling of the genome turned the ancestral hexaploid back to diploids and caused the radiation of Brassiceae species. The course of speciation was accompanied by the loss of duplicate genes and also influenced the evolution of retained genes. Of all the genes, those encoding NBS domains are typical R genes that confer resistance to invading pathogens. In this study, using the genome of Arabidopsis thaliana as a reference, we examined the loss/retention of orthologous NBS-encoding loci in the tripled Brassica rapa genome and discovered differential loss/retention frequencies. Further analysis indicated that loci of different retention ratios showed different evolutionary patterns. The loci of classesII and III (maintaining two and three syntenic loci, respectively, multi-loci) show sharper expansions by tandem duplications, have faster evolutionary rates and have more potential to be associated with novel gene functions. On the other hand, the loci that are retained at the minimal rate (keeping only one locus, class I, single locus) showed opposite patterns. Phylogenetic analysis indicated that recombination and translocation events were common among multi-loci in B. rapa, and differential evolutionary patterns between multi- and single-loci are likely the consequence of recombination. Investigations towards other gene families demonstrated different evolutionary characteristics between different gene families. The evolution of genes is more likely determined by the property of each gene family, and the whole genome triplication provided only a specific condition.
Asunto(s)
Brassica rapa/genética , Proteínas de Unión al ADN/genética , Resistencia a la Enfermedad/genética , Genes de Plantas , Proteínas de Plantas/genética , Arabidopsis/genética , Evolución Molecular , Genoma de Planta , Filogenia , Ploidias , Recombinación Genética , Análisis de Secuencia de ADN , SinteníaRESUMEN
Widely known as a severe pathogen of bean plants, the bean common mosaic virus (BCMV) has been reported to infect soybeans only sporadically and the involved strains were all found in China regions. To explore variations among soybean-infecting BCMV strains, hundreds of soybean mosaic leave samples were collected throughout China, with a total of 30 BCMV isolates detected and their genomes sequenced. These newly obtained genomes, together with 16 other BCMV genomes available in GenBank were examined from multiple aspects to characterize BCMV evolutionary processes. Phylogenetic analysis showed that both soybean-infecting BCMVs (group I) and peanut-infecting BCMVs (group II) are distantly related to other BCMVs, suggesting ancestral differentiation and host adaptation. Genetic variation analysis showed that P1, P3 and 6K2 genes and the beginning portion of CP gene showed higher levels of variation relative to other genes. Moreover, selection analyses further confirmed that a number of sites within the P1 and P3 genes have suffered positive selection. These obtained BCMV sequences also exhibit high recombination frequencies, indicating a more dynamic evolutionary history. Finally, 12 different soybean cultivars were challenged with two BCMV isolates (DXH015 and HZZB011), with most of the cultivars successfully infected. These findings suggest that BCMV is indeed a potential threat to soybean production.
Asunto(s)
Genoma Viral , Glycine max/virología , Enfermedades de las Plantas/virología , Potyvirus/aislamiento & purificación , China , Genómica , Datos de Secuencia Molecular , Filogenia , Potyvirus/clasificación , Potyvirus/genética , Proteínas Virales/genéticaRESUMEN
BACKGROUND: Dioscorea is an important plant genus in terms of food supply and pharmaceutical applications. However, its classification and identification are controversial. DNA barcoding is a recent aid to taxonomic identification and uses a short standardized DNA region to discriminate plant species. In this study, the applicability of three candidate DNA barcodes (rbcL, matK, and psbA-trnH) to identify species within Dioscorea was tested. METHODOLOGY/PRINCIPAL FINDINGS: One-hundred and forty-eight individual plant samples of Dioscorea, encompassing 38 species, seven varieties and one subspecies, representing majority species distributed in China of this genus, were collected from its main distributing areas. Samples were assessed by PCR amplification, sequence quality, extent of specific genetic divergence, DNA barcoding gap, and the ability to discriminate between species. matK successfully identified 23.26% of all species, compared with 9.30% for rbcL and 11.63% for psbA-trnH. Therefore, matK is recommended as the best DNA barcoding candidate. We found that the combination of two or three loci achieved a higher success rate of species discrimination than one locus alone. However, experimental cost would be much higher if two or three loci, rather than a single locus, were assessed. CONCLUSIONS: We conclude that matK is a strong, although not perfect, candidate as a DNA barcode for Dioscorea identification. This assessment takes into account both its ability for species discrimination and the cost of experiments.
Asunto(s)
Código de Barras del ADN Taxonómico/métodos , Dioscorea/clasificación , Dioscorea/genética , Evolución Molecular , Genes de Plantas/genética , China , Sitios Genéticos/genética , Variación Genética , Reacción en Cadena de la Polimerasa , Análisis de Secuencia de ADN , Especificidad de la Especie , Estadísticas no ParamétricasRESUMEN
A correlation method was recently adopted to identify selection-favored 'optimal' codons from 675 bacterial genomes. Surprisingly, the identities of these optimal codons were found to track the bacterial GC content, leading to a conclusion that selection would generally shape the codon usages to the same direction as the overall mutation does. Raising several concerns, here we report a thorough comparative study on 203 well-selected bacterial species, which strongly suggest that the previous conclusion is likely an illusion. Firstly, the previous study did not preclude species that are suffering weak or no selection pressures on their codon usages. For these species, as showed in this study, the optimal codon identities are prone to be incorrect and follow GC content. Secondly, the previous study only adopted the correlation method, without considering another method to test the reliability of inferred optimal codons. Actually by definition, optimal codons can also be identified by simply comparing codon usages between high- and low-expression genes. After using both methods to identify optimal codons for the selected species, we obtained highly conflicting results, suggesting at least one method is misleading. Further we found a critical problem of correlation method at the step of calculating gene bias level. Due to a failure of accurately defining the background mutation, the problem would result in wrong optimal codon identities. In other words, partial mutational effects on codon choices were mistakenly regarded as selective influences, leading to incorrect and biased optimal codon identities. Finally, considering the translational dynamics, optimal codons identified by comparison method can be well-explained by tRNA compositions, whereas optimal codons identified by correlation method can not be. For all above reasons, we conclude that real optimal codons actually do not track the genomic GC content, and correlation method is misleading in identifying optimal codons and better be avoided.