Genetic Signature Controlling Root System Architecture in Diverse Spring Wheat Germplasm.

Roots are the main sensing organ, initiating multiple signaling pathways in response to abiotic factors, including nutrients, drought, and salt stress. A focus on improving the root system architecture is a key strategy to mitigate these stresses in wheat crop. In the present study, a diversity panel comprising indigenous landraces and historical cultivars from Pakistan was characterized for the root system architecture (RSA) and important loci were identified using a genome-wide association study (GWAS). RSA of the diversity panel was characterized 30 days after sowing in brunch tubes, and root images were taken. A high-throughput root imaging analysis using Rhizovision software was performed by setting the scale to extract the eight RSA traits and four plant biomass-related traits. GWAS identified 323 association signals for 12 root and biomass traits present on all wheat chromosomes, while the most important and reliable genetic loci (based on pleotropic loci and candidate genes) were identified on chromosomes 2A, 2B, 5A, 5D, 6A, 7B, and 7D for RSA. SNP annotation and transcriptome profiling identified nine candidate genes regulating the RSA and plant biomass traits, including ROOTLESS WITH UNDETECTABLE MERISTEM1, MYB TRANSCRIPTION FACTOR4, BRASSINOSTEROID INSENSITIVE1, SLENDER RICE1, AUXIN-RESPONSIVE FACTOR25, SCARECROW, NARROW LEAF2, PIN-FORMED1 AND PHOSPHATE TRANSCRIPTION FACTOR1. This study provided pre-breeding information for deep-rooting genotypes and associated markers that will accelerate the incorporation of such traits in breeding.

Evaluation of functional kompetitive allele-specific PCR (KASP) markers for selection of drought-tolerant wheat (Triticum aestivum) genotypes.

Wheat (Triticum aestivum ) is a major crop around the globe and different techniques are being used for its productivity enhancement. Germplasm evaluation to improve crop productivity mainly depends on accurate phenotyping and selection of genotypes with a high frequency of superior alleles related to the trait of interest. Therefore, applying functional kompetitive allele-specific PCR (KASP) markers for drought-related genes is essential to characterise the genotypes for developing future climate-resilient wheat crop. In this study, eight functional KASP markers and nine morphological traits were employed to evaluate the 40 wheat genotypes for drought tolerance. Morphological traits showed significant variation (P ≤0.05) among the genotypes, except tiller count (TC), fresh root weight (FRW) and dry root weight (DRW). PCA biplot showed that 63.3% phenotypic variation was explained by the first two PCs under control treatment, while 70.8% variation was explained under drought treatment. It also indicated that root length (RL) and primary root (PR) have considerable variations among the genotypes under both treatments and are positively associated with each other. Hence, the findings of this study suggested that both these traits could be used as a selection criterion to classify the drought-tolerant wheat genotypes. KASP genotyping accompanied by morphological data revealed that genotypes Markaz, Bhakar Star, China 2, Aas and Chakwal-50 performed better under drought stress. These outperforming genotypes could be used as parents in developing drought-tolerant wheat genotypes. Hence, KASP genotyping assay for functional genes or significant haplotypes and phenotypic evaluation are prerequisites for a modern breeding program.

RNAseq of diverse spring wheat cultivars released during last 110 years.

Here, we performed RNA-seq based expression analysis of root and leaf tissues of a set of 24 historical spring wheat cultivars representing 110 years of temporal genetic variations. This huge 130 tissues RNAseq dataset was initially used to study expression pattern of 97 genes regulating root growth and development in wheat. Root system architecture (RSA) is an important target for breeding stress-resilient and high-yielding wheat cultivars under climatic fluctuations. However, root transcriptome analysis is usually obscured due to challenges in root research due to their below ground presence. We also validated the dataset by performing correlation analysis between expression of RSA related genes in roots and leaves with 25 root traits analyzed under varying moisture conditions and 10 yield-related traits. The Pearson's correlation coefficients between root phenotypes and expression of root-specific genes varied from -0.72 to 0.78, and strong correlations with genes such as DRO1, TaMOR, ARF4, PIN1 was observed. The presented datasets have multiple uses such as a) studying the change in expression pattern of genes during time, b) differential expression of genes in two very important tissues of wheat i.e., leaf and roots, and c) studying customized expression of genes associated with important phenotypes in diverse wheat cultivars. The initial findings presented here provided key insights into understanding the transcriptomic basis of phenotypic variability of RSA in wheat cultivars.

The comparative transcriptome analysis of two green super rice genotypes with varying tolerance to salt stress.

BACKGROUND: Salinity is one of the main abiotic factors that restrict plant growth, physiology, and crop productivity is salt stress. About 33% of the total irrigated land suffers from severe salinity because of intensive underground water extraction and irrigation with brackish water. Thus, it is important to understand the genetic mechanism and identify the novel genes involved in salt tolerance for the development of climate-resilient rice cultivars. METHODS AND RESULTS: In this study, two rice genotypes with varying tolerance to salt stress were used to investigate the differential expressed genes and molecular pathways to adapt under saline soil by comparative RNA sequencing at 42 days of the seedling stage. Salt-susceptible (S3) and -tolerant (S13) genotypes revealed 3982 and 3463 differentially expressed genes in S3 and S13 genotypes. The up-regulated genes in both genotypes were substantially enriched in different metabolic processes and binding activities. Biosynthesis of secondary metabolites, phenylpropanoid biosynthesis, and plant signal transduction mechanisms were highly enriched. Salt-susceptible and -tolerant genotypes shared the same salt adaptability mechanism with no significant quantitative differences at the transcriptome level. Moreover, bHLH, ERF, NAC, WRKY, and MYB transcription factors were substantially up-regulated under salt stress. 391 out of 1806 identified novel genes involved in signal transduction mechanisms. Expression profiling of six novel genes further validated the findings from RNA-seq data. CONCLUSION: These findings suggest that the differentially expressed genes and molecular mechanisms involved in salt stress adaptation are conserved in both salt-susceptible and salt-tolerant rice genotypes. Further molecular characterization of novel genes will help to understand the genetic mechanism underlying salt tolerance in rice.

Whole-genome sequencing of diverse wheat accessions uncovers genetic changes during modern breeding in China and the United States.

Breeding has dramatically changed the plant architecture of wheat (Triticum aestivum), resulting in the development of high-yielding varieties adapted to modern farming systems. However, how wheat breeding shaped the genomic architecture of this crop remains poorly understood. Here, we performed a comprehensive comparative analysis of a whole-genome resequencing panel of 355 common wheat accessions (representing diverse landraces and modern cultivars from China and the United States) at the phenotypic and genomic levels. The genetic diversity of modern wheat cultivars was clearly reduced compared to landraces. Consistent with these genetic changes, most phenotypes of cultivars from China and the United States were significantly altered. Of the 21 agronomic traits investigated, 8 showed convergent changes between the 2 countries. Moreover, of the 207 loci associated with these 21 traits, more than half overlapped with genomic regions that showed evidence of selection. The distribution of selected loci between the Chinese and American cultivars suggests that breeding for increased productivity in these 2 regions was accomplished by pyramiding both shared and region-specific variants. This work provides a framework to understand the genetic architecture of the adaptation of wheat to diverse agricultural production environments, as well as guidelines for optimizing breeding strategies to design better wheat varieties.

Evaluation of sugarcane promising clones based on the morphophysiological traits developed from fuzz.

Sugarcane is one of the critical commercial crops and principal sources of ethanol and sugar worldwide. Unfavorable conditions and poor seed setting rates hinder variety development in sugarcane. Countries like Pakistan directly import fuzz (true seed) and other propagation material from the USA, China, Brazil, etc. In this study, we imported fuzz from China, developed 29 genotypes germinating in the glasshouse, and evaluated at field conditions along with two local checks (CPF-251 and HSF-240). Morphophysiological data were recorded, including plant height (PH), cane length (CL), internodal length (IL), tiller number (TN), brix percentage (B), cane diameter (CD), chlorophyll a (Chl. a), chlorophyll b (Chl. b), and total chlorophyll (T. Chl). Results showed highly significant (p < 0.001) differences among the sugarcane accessions for all the studied traits. High broad-sense heritability (81.89% to 99.91%) was recorded for all the studied parameters. Genetic Advance (GA) ranges from 4.6% to 65.32%. The highest GA was observed for PH (65.32%), followed by CL (63.28%). Chlorophyll leaching assay was also performed at different time points (0, 50, 100, 150, and 200 min). All the genotypes showed the same leaching trend at all times, and better performing genotypes showed less leaching compared to poor performing, indicating the high amount of cutin and wax on the leaf surface. Correlation analysis showed that PH, CL, IL, and TN had significant associations. Principal components analysis (PCA) further confirms these results. Based on PCA and correlation results, PH, CL, IL, and TN can be utilized as a selection criterion for sugarcane improvement. Genotypes such as NS-4a, NS-5, NS-6, NS-8, NS-9, and NS-15 are recommended for future breeding programs related to sugarcane variety development.

Genome-wide characterization of the NLR gene family in tomato (Solanum lycopersicum) and their relatedness to disease resistance.

Nucleotide-binding leucine-rich-repeat receptors (NLR), the largest group of genes associated with plant disease resistance (R), have attracted attention due to their crucial role in protecting plants from pathogens. Genome-wide studies of NLRs have revealed conserved domains in the annotated tomato genome. The 321 NLR genes identified in the tomato genome have been randomly mapped to 12 chromosomes. Phylogenetic analysis and classification of NLRs have revealed that 211 genes share full-length domains categorized into three major clades (CNL, TNL, and RNL); the remaining 110 NLRs share partial domains and are classified in CN, TN, and N according to their motifs and gene structures. The cis-regulatory elements of NLRs exhibit the maximum number of these elements and are involved in response to biotic and abiotic stresses, pathogen recognition, and resistance. Analysis of the phylogenetic relationship between tomato NLRs and orthologs in other species has shown conservation among Solanaceae members and variation with A. thaliana. Synteny and Ka/Ks analyses of Solanum lycopersicum and Solanum tuberosum orthologs have underscored the importance of NLR conservation and diversification from ancestral species millions of years ago. RNA-seq data and qPCR analysis of early and late blight diseases in tomatoes revealed consistent NLR expression patterns, including upregulation in infected compared to control plants (with some exceptions), suggesting the role of NLRs as key regulators in early blight resistance. Moreover, the expression levels of NLRs associated with late blight resistance (Solyc04g007060 [NRC4] and Solyc10g008240 [RIB12]) suggested that they regulate S. lycopersicum resistance to P. infestans. These findings provide important fundamental knowledge for understanding NLR evolution and diversity and will empower the broader characterization of disease resistance genes for pyramiding through speed cloning to develop disease-tolerant varieties.

Genome-wide and molecular characterization of the DNA replication helicase 2 (DNA2) gene family in rice under drought and salt stress.

Rice plants experience various biotic (such as insect and pest attack) and abiotic (such as drought, salt, heat, and cold etc.) stresses during the growing season, resulting in DNA damage and the subsequent losses in rice production. DNA Replication Helicase/Nuclease2 (DNA2) is known to be involved in DNA replication and repair. In animals and yeast DNA2 are well characterized because it has the abilities of both helicase and nuclease, it plays a crucial role in DNA replication in the nucleus and mitochondrial genomes. However; they are not fully examined in plants due to less focused on plants damage repair. To fill this research gap, the current study focused on the genome-wide identification and characterization of OsDNA2 genes, along with analyses of their transcriptional expression, duplication, and phylogeny in rice. Overall, 17 OsDNA2 members were reported to be found on eight different chromosomes (2, 3, 4, 6, 7, 9, 10, and 11). Among these chromosomes (Chr), Chr4 contained a maximum of six OsDNA2 genes. Based on phylogenetic analysis, the OsDNA2 gene members were clustered into three different groups. Furthermore, the conserved domains, gene structures, and cis-regulatory elements were systematically investigated. Gene duplication analysis revealed that OsDNA2_2 had an evolutionary relationship with OsDNA2_14, OsDNA2_5 with OsDNA2_6, and OsDNA2_1 with OsDNA2_8. Moreover, results showed that the conserved domain (AAA_11 superfamily) were present in the OsDNA2 genes, which belongs to the DEAD-like helicase superfamily. In addition, to understand the post-transcriptional modification of OsDNA2 genes, miRNAs were predicted, where 653 miRNAs were reported to target 17 OsDNA2 genes. The results indicated that at the maximum, OsDNA2_1 and OsDNA2_4 were targeted by 74 miRNAs each, and OsDNA2_9 was less targeted (20 miRNAs). The three-dimensional (3D) structures of 17 OsDNA2 proteins were also predicted. Expression of OsDNA2 members was also carried out under drought and salt stresses, and conclusively their induction indicated the possible involvement of OsDNA2 in DNA repair under stress when compared with the control. Further studies are recommended to confirm where this study will offer valuable basic data on the functioning of DNA2 genes in rice and other crop plants.

Evaluation of Green Super Rice Lines for Agronomic and Physiological Traits under Salinity Stress.

Rice (Oryza sativa) is an important staple food crop worldwide, especially in east and southeast Asia. About one-third of rice cultivated area is under saline soil, either natural saline soils or irrigation with brackish water. Salinity stress is among the devastating abiotic stresses that not only affect rice growth and crop productivity but also limit its cultivation area globally. Plants adopt multiple tolerance mechanisms at the morphological, physiological, and biochemical levels to tackle salinity stress. To identify these tolerance mechanisms, this study was carried out under both a controlled glass house as well as natural saline field conditions using 22 green super rice (GSR) lines along with two local varieties ("IRRI 6 and Kissan Basmati"). Several morpho-physiological and biochemical parameters along with stress-responsive genes were used as evaluation criteria under normal and salinity stress conditions. Correlation and Principal Component Analysis (PCA) suggested that shoot-related parameters and the salt susceptible index (SSI) can be used for the identification of salt-tolerant genotypes. Based on Agglomerative Hierarchical Cluster (AHC) analysis, two saline-tolerant ("S19 and S20") and saline-susceptible ("S3 and S24") lines were selected for further molecular evaluation. Quantitative RT-PCR was performed, and results showed that expression of 1-5-phosphoribosyl -5-5-phosphoribosyl amino methylidene amino imidazole-4-carboxamide isomerase, DNA repair protein recA, and peptide transporter PTR2 related genes were upregulated in salt-tolerant genotypes, suggesting their potential role in salinity tolerance. However, additional validation using reverse genetics approaches will further confirm their specific role in salt tolerance. Identified saline-tolerant lines in this study will be useful genetic resources for future salinity breeding programs.

CRISPR/Cas9 mediated disruption of Inositol Pentakisphosphate 2-Kinase 1 (TaIPK1) reduces phytic acid and improves iron and zinc accumulation in wheat grains.

Introduction: Phytic acid (PA) is an important antinutrient agent present in cereal grains which reduces the bioavailability of iron and zinc in human body, causing malnutrition. Inositol pentakisphosphate 2- kinase 1 (IPK1) gene has been reported to be an important gene for PA biosynthesis. Objective: A recent genome editing tool CRISPR/Cas9 has been successfully applied to develop biofortified rice by disrupting IPK1 gene, however, it remained a challenge in wheat. The aim of this study was to biofortify wheat using CRISPR/Cas9. Methods: In this study, we isolated 3 TaIPK1 homeologs in wheat designated as TaIPK1.A, TaIPK1.B and TaIPK1.D and found that the expression abundance of TaIPK1.A was stronger in early stages of grain filling. Using CRISPR/Cas9, we have disrupted TaIPK1.A gene in cv. Borlaug-2016 with two guide RNAs targeting the 1st and 2nd exons. Results: We got several genome-edited lines in the T0 generation at frequencies of 12.7% and 10.8%. Sequencing analysis revealed deletion of 1-23 nucleotides and even an addition of 1 nucleotide in various lines. Analysis of the genome-edited lines revealed a significant decrease in the PA content and an increase in iron and zinc accumulation in grains compared with control plants. Conclusion: Our study demonstrates the potential application of CRISPR/Cas9 technique for the rapid generation of biofortified wheat cultivars.

Impact of Pre-Anthesis Drought Stress on Physiology, Yield-Related Traits, and Drought-Responsive Genes in Green Super Rice.

Optimum soil water availability is vital for maximum yield production in rice which is challenged by increasing spells of drought. The reproductive stage drought is among the main limiting factors leading to the drastic reduction in grain yield. The objective of this study was to investigate the molecular and morphophysiological responses of pre-anthesis stage drought stress in green super rice. The study assessed the performance of 26 rice lines under irrigated and drought conditions. Irrigated treatment was allowed to grow normally, while drought stress was imposed for 30 days at the pre-anthesis stage. Three important physiological traits including pollen fertility percentage (PFP), cell membrane stability (CMS), and normalized difference vegetative index (NDVI) were recorded at anthesis stage during the last week of drought stress. Agronomic traits of economic importance including grain yield were recorded at maturity stage. The analysis of variance demonstrated significant variation among the genotypes for most of the studied traits. Correlation and principal component analyses demonstrated highly significant associations of particular agronomic traits with grain yield, and genetic diversity among genotypes, respectively. Our study demonstrated a higher drought tolerance potential of GSR lines compared with local cultivars, mainly by higher pollen viability, plant biomass, CMS, and harvest index under drought. In addition, the molecular basis of drought tolerance in GSR lines was related to upregulation of certain drought-responsive genes including OsSADRI, OsDSM1, OsDT11, but not the DREB genes. Our study identified novel drought-responsive genes (LOC_Os11g36190, LOC_Os12g04500, LOC_Os12g26290, and LOC_Os02g11960) that could be further characterized using reverse genetics to be utilized in molecular breeding for drought tolerance.

Role of Wheat Phosphorus Starvation Tolerance 1 Genes in Phosphorus Acquisition and Root Architecture.

The wheat plant requires elevated phosphorus levels for its normal growth and yield, but continuously depleting non-renewable phosphorus reserves in the soil is one of the biggest challenges in agricultural production worldwide. The Phosphorus Starvation Tolerance 1 (PSTOL1) gene has been reported to play a key role in efficient P uptake, deeper rooting, and high yield in rice. However, the function of the PSTOL1 gene in wheat is still unclear. In this study, a total of 22 PSTOL1 orthologs were identified in the wheat genome, and found that wheat PSTOL1 orthologs are unevenly distributed on chromosomes, and these genes were under strong purifying selection. Under different phosphorus regimes, wheat PSTOL1 genes showed differential expression patterns in different tissues. These results strengthen the classification of Pakistan-13 as a P-efficient cultivar and Shafaq-06 as a P-inefficient cultivar. Phenotypic characterization demonstrated that Pakistan-13 wheat cultivar has significantly increased P uptake, root length, root volume, and root surface area compared to Shafaq-06. Some wheat PSTOL1 orthologs are co-localized with phosphorus starvation's related quantitative trait loci (QTLs), suggesting their potential role in phosphorus use efficiency. Altogether, these results highlight the role of the wheat PSTOL1 genes in wheat P uptake, root architecture, and efficient plant growth. This comprehensive study will be helpful for devising sustainable strategies for wheat crop production and adaptation to phosphorus insufficiency.

Genome-Wide Association and Genomic Prediction for Stripe Rust Resistance in Synthetic-Derived Wheats.

Stripe rust caused by Puccnina striiformis (Pst) is an economically important disease attacking wheat all over the world. Identifying and deploying new genes for Pst resistance is an economical and long-term strategy for controlling Pst. A genome-wide association study (GWAS) using single nucleotide polymorphisms (SNPs) and functional haplotypes were used to identify loci associated with stripe rust resistance in synthetic-derived (SYN-DER) wheats in four environments. In total, 92 quantitative trait nucleotides (QTNs) distributed over 65 different loci were associated with resistance to Pst at seedling and adult plant stages. Nine additional loci were discovered by the linkage disequilibrium-based haplotype-GWAS approach. The durable rust-resistant gene Lr34/Yr18 provided resistance in all four environments, and against all the five Pst races used in this study. The analysis identified several SYN-DER accessions that carried major genes: either Yr24/Yr26 or Yr32. New loci were also identified on chr2B, chr5B, and chr7D, and 14 QTNs and three haplotypes identified on the D-genome possibly carry new alleles of the known genes contributed by the Ae. tauschii founders. We also evaluated eleven different models for genomic prediction of Pst resistance, and a prediction accuracy up to 0.85 was achieved for an adult plant resistance, however, genomic prediction for seedling resistance remained very low. A meta-analysis based on a large number of existing GWAS would enhance the identification of new genes and loci for stripe rust resistance in wheat. The genetic framework elucidated here for stripe rust resistance in SYN-DER identified the novel loci for resistance to Pst assembled in adapted genetic backgrounds.

Histone Deacetylase (HDAC) Gene Family in Allotetraploid Cotton and Its Diploid Progenitors: In Silico Identification, Molecular Characterization, and Gene Expression Analysis under Multiple Abiotic Stresses, DNA Damage and Phytohormone Treatments.

Histone deacetylases (HDACs) play a significant role in a plant's development and response to various environmental stimuli by regulating the gene transcription. However, HDACs remain unidentified in cotton. In this study, a total of 29 HDACs were identified in allotetraploid Gossypium hirsutum, while 15 and 13 HDACs were identified in Gossypium arboretum and Gossypium raimondii, respectively. Gossypium HDACs were classified into three groups (reduced potassium dependency 3 (RPD3)/HDA1, HD2-like, and Sir2-like (SRT) based on their sequences, and Gossypium HDACs within each subgroup shared a similar gene structure, conserved catalytic domains and motifs. Further analysis revealed that Gossypium HDACs were under a strong purifying selection and were unevenly distributed on their chromosomes. Gene expression data revealed that G. hirsutum HDACs were differentially expressed in various vegetative and reproductive tissues, as well as at different developmental stages of cotton fiber. Furthermore, some G. hirsutum HDACs were co-localized with quantitative trait loci (QTLs) and single-nucleotide polymorphism (SNPs) of fiber-related traits, indicating their function in fiber-related traits. We also showed that G. hirsutum HDACs were differentially regulated in response to plant hormones (abscisic acid (ABA) and auxin), DNA damage agent (methyl methanesulfonate (MMS)), and abiotic stresses (cold, salt, heavy metals and drought), indicating the functional diversity and specification of HDACs in response to developmental and environmental cues. In brief, our results provide fundamental information regarding G. hirsutum HDACs and highlight their potential functions in cotton growth, fiber development and stress adaptations, which will be helpful for devising innovative strategies for the improvement of cotton fiber and stress tolerance.

Comparative Genome-wide Analysis and Expression Profiling of Histone Acetyltransferase (HAT) Gene Family in Response to Hormonal Applications, Metal and Abiotic Stresses in Cotton.

Post-translational modifications are involved in regulating diverse developmental processes. Histone acetyltransferases (HATs) play vital roles in the regulation of chromation structure and activate the gene transcription implicated in various cellular processes. However, HATs in cotton, as well as their regulation in response to developmental and environmental cues, remain unidentified. In this study, 9 HATs were identified from Gossypium raimondi and Gossypium arboretum, while 18 HATs were identified from Gossypium hirsutum. Based on their amino acid sequences, Gossypium HATs were divided into three groups: CPB, GNAT, and TAFII250. Almost all the HATs within each subgroup share similar gene structure and conserved motifs. Gossypium HATs are unevenly distributed on the chromosomes, and duplication analysis suggests that Gossypium HATs are under strong purifying selection. Gene expression analysis showed that Gossypium HATs were differentially expressed in various vegetative tissues and at different stages of fiber development. Furthermore, all the HATs were differentially regulated in response to various stresses (salt, drought, cold, heavy metal and DNA damage) and hormones (abscisic acid (ABA) and auxin (NAA)). Finally, co-localization of HAT genes with reported quantitative trait loci (QTL) of fiber development were reported. Altogether, these results highlight the functional diversification of HATs in cotton growth and fiber development, as well as in response to different environmental cues. This study enhances our understanding of function of histone acetylation in cotton growth, fiber development, and stress adaptation, which will eventually lead to the long-term improvement of stress tolerance and fiber quality in cotton.

In silico studies of C 3 metabolic pathway proteins of wheat (Triticum aestivum).

Photosynthesis is essential for plant productivity and critical for plant growth. More than 90% of plants have a C3 metabolic pathway primarily for carbon assimilation. Improving crop yields for food and fuel is a major challenge for plant biology. To enhance the production of wheat there is need to adopt the strategies that can create the change in plants at the molecular level. During the study we have employed computational bioinformatics and interactomics analysis of C3 metabolic pathway proteins in wheat. The three-dimensional protein modeling provided insight into molecular mechanism and enhanced understanding of physiological processes and biological systems. Therefore in our study, initially we constructed models for nine proteins involving C3 metabolic pathway, as these are not determined through wet lab experiment (NMR, X-ray Crystallography) and not available in RCSB Protein Data Bank and UniProt KB. On the basis of docking interaction analysis, we proposed the schematic diagram of C3 metabolic pathway. Accordingly, there also exist vice versa interactions between 3PGK and Rbcl. Future site and directed mutagenesis experiments in C3 plants could be designed on the basis of our findings to confirm the predicted protein interactions.