Bioinformatic tools uncover the C-terminal strand of Rubisco's large subunit as hot-spot for specificity-enhancing mutations

Insights into adaptive evolution of plastomes in Stipa L. (Poaceae)

BACKGROUND

The study presents results of research on the evolution of plastid genomes in Stipa L. which is a large genus of the Poaceae family, comprising species diverse in terms of geographic distribution, growing under highly variated habitat conditions. Complete plastome sequences of 43 taxa from Stipeae and Ampelodesmae tribes were analyzed for the variability of the coding regions against the background of phylogenetic relationships within the genus Stipa. The research hypothesis put forward in our research was that some of coding regions are affected by a selection pressure differentiated between individual phylogenetic lines of Stipa, potentially reducing the phylogenetic informativeness of these CDS. The study aimed to answer the question, which genes evolve in Stipa most rapidly and what kind of changes in the properties of encoded amino acids this entails. Another goal of this research was to find out whether individual genes are affected by positive selection and finally, whether selective pressure is uniform within the genus or does it vary between particular evolutionary lines within the genus.

RESULTS

Results of our study proved the presence of selective pressure in 11 genes: ccsA, matK, ndhC, ndhF, ndhK, rbcL, rpoA rpoC1, rpoC2, rps8 and rps11. For the first time the effect of positive selection on the rps8, rps11, and ndhK genes was documented in grasses. The varied pace of evolution, different intensity and effects of selective pressure have been demonstrated between particular phylogenetic lines of the genus tested.

CONCLUSIONS

Positive selection in plastid genome in Stipa mostly affects photosynthetic genes. The potential strongest adaptive pressure was observed in the rbcL gene, especially in the oldest evolutionary group comprising Central Asian high-mountain species: S. basiplumosa, S. klimesii, S. penicillata and S. purpurea, where adaptive pressure probably affected the amino acids directly related to the efficiency of CO₂ assimilation.

The complete plastome of macaw palm [Acrocomia aculeata (Jacq.) Lodd. ex Mart.] and extensive molecular analyses of the evolution of plastid genes in Arecaceae

The plastome of macaw palm was sequenced allowing analyses of evolution and molecular markers. Additionally, we demonstrated that more than half of plastid protein-coding genes in Arecaceae underwent positive selection. Macaw palm is a native species from tropical and subtropical Americas. It shows high production of oil per hectare reaching up to 70% of oil content in fruits and an interesting plasticity to grow in different ecosystems. Its domestication and breeding are still in the beginning, which makes the development of molecular markers essential to assess natural populations and germplasm collections. Therefore, we sequenced and characterized in detail the plastome of macaw palm. A total of 221 SSR loci were identified in the plastome of macaw palm. Additionally, eight polymorphism hotspots were characterized at level of subfamily and tribe. Moreover, several events of gain and loss of RNA editing sites were found within the subfamily Arecoideae. Aiming to uncover evolutionary events in Arecaceae, we also analyzed extensively the evolution of plastid genes. The analyses show that highly divergent genes seem to evolve in a species-specific manner, suggesting that gene degeneration events may be occurring within Arecaceae at the level of genus or species. Unexpectedly, we found that more than half of plastid protein-coding genes are under positive selection, including genes for photosynthesis, gene expression machinery and other essential plastid functions. Furthermore, we performed a phylogenomic analysis using whole plastomes of 40 taxa, representing all subfamilies of Arecaceae, which placed the macaw palm within the tribe Cocoseae. Finally, the data showed here are important for genetic studies in macaw palm and provide new insights into the evolution of plastid genes and environmental adaptation in Arecaceae.

One-third of the plastid genes evolved under positive selection in PACMAD grasses

We demonstrate that rbcL underwent strong positive selection during the C ₃ -C ₄ photosynthetic transitions in PACMAD grasses, in particular the 3' end of the gene. In contrast, selective pressures on other plastid genes vary widely and environmental drivers remain to be identified. Plastid genomes have been widely used to infer phylogenetic relationships among plants, but the selective pressures driving their evolution have not been systematically investigated. In our study, we analyse all protein-coding plastid genes from 113 species of PACMAD grasses (Poaceae) to evaluate the selective pressures driving their evolution. Our analyses confirm that the gene encoding the large subunit of RubisCO (rbcL) evolved under strong positive selection after C₃-C₄ photosynthetic transitions. We highlight new codons in rbcL that underwent parallel changes, in particular those encoding the C-terminal part of the protein. C₃-C₄ photosynthetic shifts did not significantly affect the evolutionary dynamics of other plastid genes. Instead, while two-third of the plastid genes evolved under purifying selection or neutrality, 25 evolved under positive selection across the PACMAD clade. This set of genes encode for proteins involved in diverse functions, including self-replication of plastids and photosynthesis. Our results suggest that plastid genes widely adapt to changing ecological conditions, but factors driving this evolution largely remain to be identified.

'Dynamic Distance' Reaction Coordinate for Competing Bonds: Applications in Classical and Ab Initio Simulations

A versatile reaction coordinate, the "dynamic distance", is introduced for the study of reactions involving the rupture and formation of a series of chemical bonds or contacts. The dynamic distance is a mass-weighted mean of selected distances. When implemented as a generalized constraint, the dynamic distance is particularly suited for driving activated processes by controlled increase during a simulation. As a single constraint acting upon multiple degrees of freedom, the sequence of events along the resulting reaction pathway is determined unambiguously by the underlying energy landscape. Free energy profiles can be readily obtained from the mean constraint force. In this paper both theoretical aspects and numerical implementation are discussed, and the unique and diverse properties of this reaction coordinate are demonstrated using three examples:  In the framework of Car-Parrinello molecular dynamics, we present results for the prototypical double proton-transfer reaction in formic acid dimer and the photocycle of the guanine-cytosine DNA base pair. As a classical mechanical example, the opening of the binding pocket of the enzyme rubisco is analyzed.

Structure-function studies with the unique hexameric form II ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco) from Rhodopseudomonas palustris

The first x-ray crystal structure has been solved for an activated transition-state analog-bound form II ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco). This enzyme, from Rhodopseudomonas palustris, assembles as a unique hexamer with three pairs of catalytic large subunit homodimers around a central 3-fold symmetry axis. This oligomer arrangement is unique among all known Rubisco structures, including the form II homolog from Rhodospirillum rubrum. The presence of a transition-state analog in the active site locked the activated enzyme in a "closed" conformation and revealed the positions of critical active site residues during catalysis. Functional roles of two form II-specific residues (Ile(165) and Met(331)) near the active site were examined via site-directed mutagenesis. Substitutions at these residues affect function but not the ability of the enzyme to assemble. Random mutagenesis and suppressor selection in a Rubisco deletion strain of Rhodobacter capsulatus identified a residue in the amino terminus of one subunit (Ala(47)) that compensated for a negative change near the active site of a neighboring subunit. In addition, substitution of the native carboxyl-terminal sequence with the last few dissimilar residues from the related R. rubrum homolog increased the enzyme's kcat for carboxylation. However, replacement of a longer carboxyl-terminal sequence with termini from either a form III or a form I enzyme, which varied both in length and sequence, resulted in complete loss of function. From these studies, it is evident that a number of subtle interactions near the active site and the carboxyl terminus account for functional differences between the different forms of Rubiscos found in nature.

Ribulose-1,5-bisphosphate carboxylase/oxygenase from thermophilic cyanobacterium Thermosynechococcus elongatus

Ribulose-1,5-bisphosphate carboxylase/oxygenase (rubisco) can be divided into two branches: the "red-like type" of marine algae and the "green-like type" of cyanobacteria, green algae, and higher plants. We found that the "green-like type" rubisco from the thermophilic cyanobacterium Thermosynechococcus elongatus has an almost 2-fold higher specificity factor compared with rubiscos of mesophilic cyanobacteria, reaching the values of higher plants, and simultaneously revealing an improvement in enzyme thermostability. The difference in the activation energies at the transition stages between the oxygenase and carboxylase reactions for Thermosynechococcus elongatus rubisco is very close to that of Galdieria partita and significantly higher than that of spinach. This is the first characterization of a "green-like type" rubisco from thermophilic organism.

Ab initio study of molecular interactions in higher plant and Galdieria partita Rubiscos with the fragment molecular orbital method

Ribulose bisphosphate carboxylase/oxygenase (Rubisco) from one of the thermophilic red algae Galdieria partita with a high specificity factor shows a characteristic difference from higher plant Rubisco in structural change. We investigate such a difference by evaluating the inter-fragment interaction energy (IFIE) value with fragment molecular orbital (FMO) method in comparison to experimental structural studies. We found some important residues which determine the loop6 stability or which make difference in the structure between higher plant and G. partita Rubiscos. We found that amino acid change of LYS18 to ILE18 is important for the difference in location at which anion binding site is occupied, P1alpha or P1beta, when inorganic anions are bound to the enzyme. Occupation of P2 anion binding site makes the stabilizing interaction between LYS128 and the loop6 stronger. Amino acid change of HIS386 to GLN386 contributed to the difference in the loop6 stability, while amino acid change of MET472 to THR472 did not contribute to it. It is confirmed that the patterns of interactions among THR65, THR67, and THR462 are consistent with previous experimental discussions. However, we found a case that THR65 was not stabilized with anion at P1alpha binding site in a closed-state structure of G. partita Rubisco.