Genetic determinants controlling maize rubisco activase gene expression and a comparison with rice counterparts
BMC Plant Biology volume 19, Article number: 351 (2019)
Rubisco activase (RCA) regulates the activity of Rubisco and is a key enzyme of photosynthesis. RCA expression was widely reported to affect plant photosynthesis and crop yield, but the molecular basis of natural variation in RCA expression in a wide range of maize materials has not been fully elucidated.
In this study, correlation analysis in approximately 200 maize inbred lines revealed a significantly positive correlation between the expression of maize RCA gene ZmRCAβ and grain yield. A genome-wide association study revealed both cis-expression quantitative trait loci (cis-eQTLs) and trans-eQTLs underlying the expression of ZmRCAβ, with the latter playing a more important role. Further allele mining and genetic transformation analysis showed that a 2-bp insertion and a 14-bp insertion in the promoter of ZmRCAβ conferred increased gene expression. Because rice is reported to have higher RCA gene expression than does maize, we subsequently compared the genetic factors underlying RCA gene expression between maize and rice. The promoter activity of the rice RCA gene was shown to be stronger than that of the maize RCA gene, suggesting that replacing the maize RCA gene promoter with that of the rice RCA gene would improve the expression of RCA in maize.
Our results revealed two DNA polymorphisms regulating maize RCA gene ZmRCAβ expression, and the RCA gene promoter activity of rice was stronger than that of maize. This work increased understanding of the genetic mechanism that underlies RCA gene expression and identify new targets for both genetic engineering and selection for maize yield improvement.
Ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco) catalyzes the first step of photosynthetic carbon assimilation and is a rate-limiting enzyme of photosynthesis. In higher plants, Rubisco activity is regulated by a second enzyme named Rubisco activase (RCA) [1,2,3,4,5,6]. RCA removes inhibitory sugar phosphates from the active sites of Rubisco in an ATP-dependent way and thus activates Rubisco . The endogenous levels of RCA have been widely reported to affect plant photosynthesis and crop yield. In rice, RCA plays an important role in the regulation of non-steady-state photosynthesis . Overexpression of the RCA gene improves rice grain yield, whereas anti-sense expression shows the opposite effect [9, 10]. In wheat, RCA expression showed a significantly positive linear correlation with plant productivity . In soybean, the expression levels of two RCA genes showed significant correlations with photosynthesis and grain yield [12, 13]. In maize, high-yielding population showed higher RCA expression [14, 15], and the RCA expression displayed significant and positive correlations with grain yield in a wide range of maize inbred lines .
Characterization of the genetic elements regulating RCA gene expression could facilitate the modulation of the RCA gene for improving crop yield. For this purpose, genetic determinants that control the expression of RCA genes have been identified in various crops. Functional analysis of the rice RCA gene promoter revealed several cis-elements that regulate its expression . Expression quantitative trait locus (eQTL) mapping identified several eQTLs that control the expression of soybean RCA genes [12, 13]. A screen of a leaf cDNA library using a yeast one-hybrid (Y1H) system identified two bZIP transcription factors that interact with the promoter of the soybean RCA gene . Maize is one of the most important C4 crops worldwide. Although C4 monocot plants including maize show higher photosynthetic efficiency due to an extra biochemical CO2-concentrating mechanism (CCM) yielding higher CO2 partial pressure (pCO2) at Rubisco catalytic sites, their Rubisco activation status is generally lower than that of C3 monocot plants . Therefore, there should be much room to increase maize Rubisco activation state and subsequently photosynthesis. Because of the primary role of RCA to activate Rubisco in vivo, mapping the genetic factors underlying RCA expression in maize is of paramount importance.
Recently, based on linkage analysis using a maize recombinant inbred line (RIL) population derived from two parental inbred lines, we identified several eQTLs underlying RCA gene expression . Although these eQTLs can help to understand the regulation of the RCA gene in maize, the molecular basis of natural variation in RCA expression in a wide range of maize materials has not been fully elucidated. Genome-wide association analysis (GWAS) based on linkage disequilibrium (LD) studies in a natural population allows for evaluation of multiple alleles underlying the natural variation of a target trait. Additionally, GWAS can increase resolution and facilitate gene discovery .
Rice, one of the most important C3 crops, shows a higher level of RCA gene expression than that of maize . The difference in RCA gene expression between maize and rice most likely required the evolution of a different regulatory mechanism. Therefore, the higher level of RCA gene expression in rice than that in maize suggests that rice might have stronger RCA gene regulatory factors, in the form of either cis- or trans-regulatory factors, than maize. To date, no information is available on the similarities or differences in the RCA gene regulatory elements between maize and rice. If rice is confirmed to have a stronger RCA regulator than maize, future high-yield breeding of maize can focus on replacing its weaker RCA gene regulator with the stronger one from rice.
In maize, two RCA genes, ZmRCAα and ZmRCAβ, have been identified and cloned. The ZmRCAβ expression level was approximately 10-fold that of ZmRCAα. Additionally, compared with ZmRCAα, ZmRCAβ showed a much higher correlation with grain yield . In most studies, only one RCA gene, OsRCA, has been identified in the rice genome . Alternative splicing of this gene produces two RCA isoforms, the short β form and the long α form . In this study, GWAS was employed to map eQTLs for the maize RCA gene, ZmRCAβ, and the rice RCA gene, OsRCA. Subsequently, sequencing of target DNA segments and candidate gene association were employed to screen for elite alleles of the RCA gene promoter in these two species. Finally, genetic transformation experiments were conducted to examine the differences in RCA gene promoter activity between maize and rice. Our findings showed that a 2-bp insertion and a 14-bp insertion in the promoter of ZmRCAβ conferred increased gene expression and that the rice RCA gene promoter had higher activity than that of the maize RCA gene promoter.
Variation of ZmRCAβ expression and its correlation with grain yield
Because alternative splicing of the ZmRCAβ precursor RNA creates two transcripts , we measured the transcript expression of this gene by designing primers that detect the two transcripts simultaneously (Additional file 1: Table S1). The transcript expression of ZmRCAβ and grain yield exhibited wide phenotypic variations in the maize association populations (Table 1 and Fig. 1). The means, standard deviations, ranges, skewness, broad-sense heritability values, and analysis of variance (ANOVA) for ZmRCAβ transcript expression and grain yield are presented in Table 1. The average relative value of ZmRCAβ transcript expression were 1.46 and 1.38 in 2015 and 2016 growing season, respectively. The variation range of ZmRCAβ expression was larger, 95% of lines fell in the range of 0.21–3.13 in 2015 and 0.17–2.93 in 2016. The average of grain yield were 45.60 and 45.11 g per plant (g/plant) in 2015 and 2016 growing season, respectively. The variation range of grain yield was larger, 95% of lines fell in the range of 19.83–71.37 g/plant in 2015 and 17.32–72.90 g/plant in 2016. ANOVA indicated that ZmRCAβ transcript expression and grain yield were significantly affected by genotype (P < 0.01). The heritability (H2) estimates for ZmRCAβ transcript expression ranged from 75 to 80%, and for grain yield from 76 to 86% (Table 1). Overall, the transcript expression of ZmRCAβ and grain yield clearly exhibited wide natural variation in the maize association populations.
In our previous study, we determined the transcript and protein expressions of ZmRCAβ and a newly discovered RCA gene ZmRCAα in 123 maize inbred lines. The transcript expression of ZmRCAβ was observed to correlate well with its protein expression; both the transcript and protein expressions of ZmRCAβ were significantly correlated with grain yield, and more importantly, the correlations were much higher than those between ZmRCAα and grain yield . In the current study, we focused on ZmRCAβ transcript expression, and observed that during the two growing seasons of 2015 and 2016, the ZmRCAβ transcript expression also displayed a significant correlation with the corresponding grain yield (P < 0.01). Pearson phenotypic correlation coefficient between ZmRCAβ transcript expression and grain yield was 0.38 in the 2015 growing season and 0.36 in the 2016 growing season (Fig. 2). The above results indicate that the ZmRCAβ transcript expression could affect maize grain yield.
Genome-wide association study for ZmRCAβ expression
To identify the genetic loci associated with ZmRCAβ expression, GWAS was performed by using 558,629 genomic single nucleotide polymorphisms (SNP) and the ZmRCAβ transcript expression data generated in an association population comprising 222 and 182 maize lines in 2015 and 2016, respectively (see Plant materials and methods). To avoid spurious associations, such as the historical relationships and selection patterns of the maize lines in the association population, multiple statistical models were evaluated in the GWAS. Visual observation of quantile-quantile (Q-Q) plots provides an indication of the accuracy of the model used to analyze the data. As shown in the Q-Q plots (Fig. 3A, C), the P-values of the mixed linear model (MLM) were close to the expected values, indicating this model was suitable for reducing the effects of population structure on ZmRCAβ transcript expression. Thus, we conducted GWAS for ZmRCAβ transcript expression with the MLM model.
Three and 12 SNPs were identified with significant marker-trait associations with ZmRCAβ transcript expression at the Bonferroni-adjusted significance threshold in 2015 and 2016 growing seasons, respectively (Table 2 and Fig. 3B, D). These SNPs were distributed on chromosomes 1, 4, 6, 8, and 9, and the phenotypic variation explained by each SNP ranged from 13.65 to 29.71%. Interestingly, GWAS for grain yield showed that the significant marker chr6.S-28129864 for ZmRCAβ expression (Table 2) was also significantly associated with grain yield (Additional file 9: Table S6). Of the significant marker-trait associations for ZmRCAβ expression, the marker PZE-104005992 was repeatedly detected in 2015 and 2016 growing seasons. A survey of the maize reference genome sequence (https://www.maizegdb.org/) revealed that this marker was located at 3.6 Mb upstream of ZmRCAβ on chromosome 4.
Candidate gene association analysis of the polymorphisms in the promoter region of ZmRCAβ
The 1700-bp region of the ZmRCAβ promoter, which was previously shown to be sufficient for the proper regulation of expression in plants [25, 26], was cloned and sequenced in the maize association population. The promoter sequences of 186 inbred lines of different geographic origins were successfully cloned and sequenced. Multiple sequence alignment of the promoter region revealed 79 polymorphic sites in the 186 maize inbred lines (Additional file 2: Table S2). To determine whether these sequence polymorphisms were related to the observed variation in ZmRCAβ transcript expression, candidate-gene association analysis was performed. Two polymorphisms, a 2-bp indel and 14-bp indel, were repeatedly detected to be significantly associated with the expression of ZmRCAβ in both 2015 and 2016 growing seasons at P < 0.01 (Table 3). The 2-bp indel located 144-bp upstream of the start codon of ZmRCAβ explained 4.24 and 6.15% of the phenotypic variation in 2015 and 2016, respectively. The 14-bp indel located 210-bp upstream of the start codon of ZmRCAβ explained 3.80 and 5.98% of the phenotypic variation in 2015 and 2016, respectively.
Haplotype analysis for ZmRCAβ promoter
Based on the polymorphisms significantly related with RCA expression, the promoter sequence of ZmRCAβ was divided into 3 haplotypes, ZmHap1, ZmHap2 and ZmHap3 (Fig. 4A). The inbred lines harboring each of these haplotypes are shown in Additional file 3: Table S3. The inbred lines carrying the ZmHap1 showed the highest ZmRCAβ expression and the highest grain yield, and those carrying the ZmHap3 showed the lowest ZmRCAβ expression and the lowest grain yield (Fig. 4B, C).
A dual-luciferase reporter system was used to verify the ZmHap1 promoter of ZmRCAβ had higher activity than that of the ZmHap3 promoter. We cloned ZmHap1 and ZmHap3 promoters from inbred lines CIMBL57 and CIMBL88, respectively. These two inbred lines showed the highest and lowest levels of ZmRCAβ expression, respectively. The ZmHap1 or ZmHap3 promoter-driven LUC (ZmHap1 or ZmHap3 pro-LUC) and CaMV35S promoter-driven REN (35S-REN; as an internal control) were constructed in the same plasmid (Fig. 5A), and expressed in maize protoplasts, rice protoplasts, and tobacco (Nicotiana benthamiana) leaves. The dual luciferase expression system with no promoter to drive LUC was used as a negative control (Fig. 5A). The LUC/REN ratio, which reflects in vivo promoter activity, was monitored. The LUC/REN value of ZmHap1 promoter was significantly higher than that of the ZmHap3 promoter in maize protoplasts, rice protoplasts, and the leaves of tobacco (Fig. 5B, C). We also compared the activity of ZmHap1 and ZmHap3 promoters using rice callus genetic transformation experiment, and observed a similar result (Additional file 4: Figure S1). These results confirmed that the ZmHap1 promoter had higher activity than that of the ZmHap3 promoter. Therefore, ZmHap1 can be regarded as an elite allele for the ZmRCAβ promoter.
Comparison of RCA gene regulatory elements between maize and rice
To compare the regulating strength of RCA gene expression determinants between maize and rice, the transcript expression of rice RCA gene OsRCA was determined in an association population containing 114 cultivars with extensive genetic diversity, and the genetic determinants underlying the transcript expression of OsRCA were identified using the same method employed in maize as described above. Similar to the measurement of ZmRCAβ expression, because alternative splicing of the OsRCA precursor RNA produces two transcripts , we also measured the transcript expression of this gene by designing primers that detect the two transcripts simultaneously (Additional file 1: Table S1). The transcript expression of OsRCA exhibited wide natural variations in the rice association population (Table 1 and Fig. 1). A 12.75-fold difference occurred among the cultivars in the rice association population (Table 1). GWAS using the rice association population revealed 5 SNPs significantly associated with OsRCA transcript expression (Additional file 5: Table S4 and Additional file 6: Figure S2).
The promoter sequence of 1500-bp was cloned and sequenced in 107 rice cultivars. Candidate-gene association analysis revealed 12 polymorphisms significantly associated with OsRCA transcript expression in its promoter region (Additional file 7: Table S5). Based on these polymorphisms, the OsRCA promoter was divided into 5 haplotypes (Additional file 3: Table S3). The cultivars carrying OsHap1 had the highest OsRCA expression and those carrying OsHap5 had the lowest OsRCA expression (Additional file 8: Figure S3B).
The aforementioned experiments of protoplast transient expression, leaf transient expression, and callus genetic transformation were conducted to compare the promoter activity of maize ZmRCAβ and rice OsRCA. We cloned OsHap1 and OsHap5 promoters from rice cultivars Jing 185–7 and Saber, respectively. These two cultivars showed the highest and lowest levels of OsRCA expression, respectively, in the rice association population. Each of the OsHap1 and OsHap5 promoter sequences was fused with the LUC or GUS reporter gene in plasmid (Fig. 5A and Additional file 4: Figure S1A). In the protoplast and leaf transient expression experiments, the LUC/REN value of OsHap1 was significantly higher than that of OsHap5, confirming that the OsHap1 promoter had higher activity than that of the OsHap5 promoter (Fig. 5B, C). The LUC/REN value of rice stronger OsHap1 promoter was significantly higher than maize stronger ZmHap1 promoter. Notably, the LUC/REN value of the rice weaker OsHap5 promoter was even higher than the maize stronger ZmHap1 promoter although the difference observed in maize protoplast experiment was not significant (Fig. 5B, C). The similar results were also observed in the callus genetic transformation experiment (Additional file 4: Figure S1B). These results suggested that the activity of the rice RCA promoter was higher than that of the maize RCA promoter.
Mapping of favorable alleles underlying ZmRCAβ expression for maize high-yield breeding
Strong evidence exists that two traits are causally related when the two traits are significantly correlated with one another and have coincident QTLs that are consistent with the correlation . In the present study, the expression level of ZmRCAβ was positively correlated with grain yield in an association population consisting of approximately 200 maize inbred lines. Furthermore, GWAS for grain yield showed that one QTL for grain yield coincided with the eQTL for ZmRCAβ expression (Additional file 9: Table S6). We also observed similar results in our previous study using 123 maize inbred lines . These data suggest that the RCA gene might be a potential factor affecting plant productivity. Therefore, identifying the key genomic regions that regulate the expression level of ZmRCAβ could aid in maize high-yield breeding via the marker-assisted-selection (MAS) method.
The success of MAS breeding depends on the gene mapping resolution and germplasm diversity. In our previous study, we detected five eQTLs for ZmRCAβ in a RIL population containing 242 maize lines . However, these eQTLs were mapped to relatively large mapping intervals, and no causal explicit genetic elements could be inferred for regulating ZmRCAβ expression. Additionally, these eQTLs do not represent the elite alleles underlying natural variation of ZmRCAβ expression because they were mapped in a bi-parental population. GWAS combined with target region association analysis in a diverse association population can provide the opportunity to discover elite alleles of interest and delimit them into a relatively small genomic region . This method takes advantage of current and historical recombination events for gene fine mapping, particularly in maize, a species that shows high diversity and rapid LD decay. Therefore, in the present study, GWAS combined with target region association analysis was employed to discover favorable alleles underlying ZmRCAβ expression in a maize association population. The results show that the favorable alleles underlying ZmRCAβ expression were finely mapped.
In the present study, GWAS initially detected one major eQTL, PZE-104005992, on chromosome 4 in approximately 200 inbred lines with extensive genetic diversity. This eQTL was repeatedly detected across two field-growth conditions and explained a large phenotypic variation (Table 2). We also detected this eQTL in our previous study . Therefore, PZE-104005992 can be considered as a primary breeding target for modulation of ZmRCAβ expression. Subsequently, we focused on refined mapping of this eQTL. Considering that the expression of one gene is usually affected by the sequence variation in its promoter and that the genomic location of the ZmRCAβ promoter coincided with the eQTL PZE-104005992, we employed candidate-gene association analysis to the ZmRCAβ promoter. The results show that two polymorphic sites, a 2-bp indel and a 14-bp indel, in the promoter were significantly associated with gene expression (Table 3). The inbred lines carrying the 2-bp and 14-bp insertion showed higher ZmRCAβ expression and higher yields than those of other inbred lines (Fig. 4). Finally, genetic transformation experiments showed the promoter with the 2-bp and 14-bp insertion had higher activity than that of other promoters (Fig. 5 and Additional file 4: Figure S1). Based on these results, we conclude that the promoter carrying the 14-bp and 2-bp insertion might be a favorable allele for maize high-yield breeding. Functional markers are derived from polymorphic sites within genes that are causally involved in phenotypic trait variation . The 2-bp and 14-bp indels identified in the present study can be used for development of functional markers.
Genetic mechanism underlying RCA gene expression
When an eQTL is within 5 Mb of its target gene, a cis-eQTL is defined; otherwise, the eQTL is a trans-eQTL [30,31,32,33]. The identification of a cis-eQTL suggests that the genetic determinants for gene expression are adjacent to or coincident with the genomic region of the target gene. A trans-eQTL indicates that genetic determinants far from the target gene control its expression . In the present study, GWAS revealed 14 and 5 SNPs that were significantly associated with the expression of ZmRCAβ and OsRCA, respectively (Table 2 and Additional file 5: Table S4). Of the 14 significant SNPs for ZmRCAβ, the marker PZE-104005992 was 3.6 Mb from ZmRCAβ in the maize reference genome (https://www.maizegdb.org/) and was called a cis-eQTL, whereas the other SNPs were trans-eQTLs. Of the 5 significant SNPs for OsRCA, the marker id11011480 was 0.7 Mb from OsRCA in the rice reference genome (http://www.ricedata.cn/gene/) and was a cis-eQTL, whereas the other SNPs were trans-eQTL. Notably, the cis-eQTL showed a much stronger association with gene expression and explained a much larger proportion of phenotypic variation than any individual trans-eQTL in both rice and maize (Table 2 and Additional file 5: Table S4). We also obtained similar results in our previous study: in maize, we identified 1 cis-eQTL and 4 trans-eQTLs for the expression of ZmRCAβ in a RIL population under two field-growth environments and observed that the cis-eQTL explained the most phenotypic variation and was stably expressed across the different environments . In soybean, we identified 13 SNPs significantly associated with the expression of the RCA gene GmRCAβ and found that the cis-eQTL had a stronger association with gene expression than that of the trans-eQTLs . These results indicate that the expression of the RCA gene might be controlled by a combination of both cis-acting and trans-acting eQTLs and that the cis-eQTL played a dominant role in regulating RCA gene expression.
Promoters can affect the action of both cis- and trans-eQTLs on the expression of a target gene via alterations in its cis-acting regulatory elements or transcription factor binding sites . Therefore, DNA sequence variation in promoters is a potential factor explaining the diversity of ZmRCAβ and OsRCA expression. In the present study, candidate-gene association analysis revealed 2 polymorphic sites in the ZmRCAβ promoter and 12 polymorphic sites in the OsRCA promoter that were significantly associated with gene expression (Table 3 and Additional file 7: Table S5). The RCA gene expression was significantly different in haplotypes of maize inbred lines or rice cultivars divided according to the detected significant association polymorphic sites (Fig. 4 and Additional file 8: Figure S3). Furthermore, genetic transformation experiments confirmed that the two types of promoters of ZmRCAβ and OsRCA had different activity (Fig. 5 and Additional file 4: Figure S1). These results suggest that promoter sequence variation can influence RCA gene expression. A similar result was also observed in soybean in which sequence variation in the GmRCAβ promoter contributed to diversity in GmRCAβ expression .
Comparison of the promoter activity between maize and rice provides a potential approach for maize high-photosynthetic-efficiency breeding
In the present study, we used qRT-PCR to determine maize and rice RCA gene expression. Although gene expression can be measured quickly and accurately, this method provides only the relative expression data of the targeted gene compared with a reference gene. Because the absolute values of reference gene expression are not available in maize and rice, we could not compare the RCA gene expression magnitude between these two species. Therefore, we focused on comparing the promoter activity of maize ZmRCAβ and rice OsRCA using genetic transformation experiments.
Different results have been reported on the cell type-specific expression pattern of maize RCA [35,36,37]. A high density of immune-gold particles for RCA was observed to be located in chloroplast stroma of both bundle sheath (BS) and mesophyll (M) cells in maize leaf . Similarly, although the authors did not point out, an immunoblot analysis of BS and M cells total proteins in maize revealed that RCA accumulated to a similar level in the two types of cells (according to Fig. 7B in the paper by Wostrikoff et al. ). In contrast, maize RCA was shown to express at a much higher level in BS cells than in M cells as revealed by large-scale quantitative proteomics  and RNA-seq analysis . If the promoter of ZmRCAβ and OsRCA has cell type-specific expression pattern between BS and M cells, it should show a different expression magnitude between maize and rice leaf tissues that contain a different ratio of M and BS cell number. In the present study, neither the promoter activity of a maize haplotype, either ZmHap1 or ZmHap3, nor the promoter activity of a rice haplotype, either OsHap1 or OsHap5, showed significantly different activity between maize and rice protoplasts (Fig. 5B). These results suggest that the ZmRCAβ and OsRCA promoter did not yield significant cell type-specific activity in BS and M cells.
Genetic transformation experiments revealed that the activity of the rice OsRCA promoter was significantly higher than that of the maize ZmRCAβ promoter (Fig. 5 and Additional file 4: Figure S1). This result is consistent with the finding that the RCA gene in maize showed a lower expression level than that in rice . Compared with maize, rice Rubisco operates in a much lower CO2 environment and at a higher activation state . To maintain the high activation state of Rubisco, the rice RCA gene promoter likely evolved a high regulating strength to keep the RCA gene expression at a high level. Considering C4 plants have an extra CCM and consequently much higher level of CO2 concentration at their Rubisco catalytic sites than C3 plants, we attempted to search for the difference in the RCA gene promoter between C3 and C4 species. A survey of the PlantCARE database  showed that, compared with C4 species including maize, the promoter from C3 species including rice have a significantly higher number of cis-acting elements associated with light regulation (Additional file 10: Table S7). Because the RCA gene is highly regulated by light , the increase in number of light regulation-related cis-elements might lead to an RCA gene promoter with increased activity. Future studies are required to compare the regulating strength of the RCA promoter between C3 and C4 plants using multiple C3 and C4 species.
As it was shown that the activity of the rice OsRCA promoter was higher than that of the maize ZmRCAβ promoter (Fig. 5 and Additional file 4: Figure S1), we suppose that maize photosynthesis could be improved by replacing the maize ZmRCAβ promoter with the rice OsRCA promoter. This replacement can combine the high CO2 environment at the Rubisco catalytic sites with a high RCA gene expression in maize plants. A concern for this replacement is that it might affect the action of some heat-stress related regulatory elements that control maize RCA gene expression. RCA is known to be particularly susceptible to heat stress [41, 42], and it was reported that heat stress induced a new RCA isoform with higher molecular weight in maize . Another concern is that the replacement might lead to a decrease of maize Rubisco content. Recently, it was reported that RCA expression negatively affected Rubisco content in transgenic rice . Further studies are required to investigate the effect of the replacement on maize RCA gene expression, Rubisco content, photosynthesis, and grain yield under both common and heat-stressed conditions.
This study of eQTL mapping and genetic transformations of promoter regions were performed to characterize the genetic factors underlying the expression of the ZmRCAβ gene in maize. Furthermore, we compared the regulating strength of the identified factor between the C4-plant maize and the C3-plant rice. Our study showed that the RCA gene might play a role in maize grain yield. A 2-bp insertion and a 14-bp insertion in the promoter conferred higher expression of maize RCA gene. The promoter activity of the rice RCA gene was stronger than that of the maize RCA gene. The two identified indels associated with maize RCA gene can be used in marker-assisted selection breeding for high photosynthetic efficiency. More importantly, people usually transfer C4-specific genes from C4-crop maize to C3-crop rice in traditional high-photosynthetic-efficiency breeding. Our study suggested a possible and novel method in high-photosynthetic-efficiency breeding, which is, in the other way around, transferring gene from C3-crop rice to C4-crop maize.
Plant materials and plant growth conditions
A panel of 437 maize lines described previously  was used. All 437 lines were genotyped for SNP markers using high-throughput genotyping platforms, and approximately 558,629 polymorphisms with minor allele frequencies (MAF) ≥0.05 were available for a GWAS . The 437 lines were planted in two growing seasons in 2015 and 2016 on the experimental farm of the Agricultural College of Yangzhou University. In each growing season, the 437 lines were planted in a randomized block design with two replicates. Each line was planted in one row per block. The length of each row was 2.5 m, the spacing between plants in each row was 0.25 m, and the spacing between rows was 0.55 m. Leaf tissues for gene expression measurement were prepared according to our previous study . At 32 d after anthesis, the leaves located nearest to the ear were harvested from three randomly selected plants of each line per block in the morning (9:00–11:30 a.m.) on a sunny day, frozen immediately in liquid nitrogen, and stored at − 80 °C before measuring ZmRCAβ expression. To minimize the impacts of differences in plant growth stage on gene expression measurements, we chose 222 and 182 lines with similar anthesis time in 2015 and 2016, respectively. The origin of the selected lines is shown in Additional file 3: Table S3.
A panel of 218 rice cultivars was used to evaluate OsRCA gene expression. All 218 cultivars were genotyped for SNP markers using high-throughput genotyping platforms, and approximately 36,900 polymorphisms with MAF ≥ 0.05 were available for GWAS . All cultivars were planted on the experimental farm of the Agricultural College of Yangzhou University in 2016. The cultivars were planted in a randomized block design with two replicates. Each cultivar was planted in one row per block. Row length was 1.5 m, row distance was 0.2 m, and plant distance in each row was 0.132 m. At the late filling stage of development, the rice flag leaves were collected individually from three randomly selected plants of each cultivar per block in the morning (9:00–11:30 a.m.) on a sunny day, frozen immediately in liquid nitrogen, and stored at − 80 °C until further use. To minimize the impacts of differences in growth stage on the measurements of gene expression, we chose 114 cultivars with similar anthesis time. The origin of the selected cultivars is shown in Additional file 3: Table S3.
Trait measurement and phenotypic data collection
Total RNA was extracted using TRIzol reagent (Vazyme, China). We used HiScript reverse transcriptase (Vazyme, China) to prepare cDNA. Gene expression levels were determined by quantitative real-time PCR (qRT-PCR) using an SYBR Premix Ex Taq kit (Vazyme, China) on a Bio-Rad CFX96 system (Applied Biosystems). The expression of endogenous reference genes ZmActin (GenBank accession number J01238) and OsActin (GenBank accession number AB047313) was used to normalize the transcript levels in each maize and rice sample, respectively. Gene specific primers (Additional file 1: Table S1) were used for qRT-PCR. The following formula was used to obtain the normalized expression of target gene for each sample: ΔΔCT = (CT, Target – CT, Actin) genotype – (CT, Target – CT, Actin) calibrator. Each sample was prepared by mixing equal amounts of the collected leaves from the three plants described above for each line (cultivar) per block. Therefore, in each growing season, each maize line or rice cultivar contained two samples from two blocks.
Maize grain yield for each line per block was estimated using the average yield of five plants in the middle of each row. At maturity, the ears of the corresponding plants were hand harvested, dried to a constant mass, and threshed; and the mean grain yield per plant was recorded.
Statistical analyses of phenotypes
Descriptive statistical analysis and correlation analyses were performed using the statistical software package SPSS Statistics 17.0 for Windows (SPSS, Inc., Chicago, IL, USA). Two-tailed ANOVA was used for comparisons between blocks and genotypes. The broad-sense heritability (H2) for RCA gene expression or grain yield in each growing season was implemented by SAS 9.2 (SAS institute, Cray, USA) and was estimated by the following equation H2 = Vg/(Vg + Ve/r), where Vg and Ve represent genetic variance and error variance, respectively, and were calculated from the ANOVA as described above; r is the block replicate number in field experiment. Frequency distribution analysis was performed using GraphPad Prism (version 7, GraphPad Software, San Diego, CA).
Genome-wide association study
TASSEL software 5.0 was used to perform GWAS . To account for the effects of population structure on the mapping panel and genetic relatedness among panel members, the kinship matrix (K) of the association panel was calculated as previously described . According to the Q-Q plots from the TASSEL 5.0 output, the MLM mode with K was chosen for the association mapping of ZmRCAβ and OsRCA expression. Markers were identified as significantly associated with expression of ZmRCAβ or OsRCA by comparison with the Bonferroni threshold (P < 1/n, −logP > 5.75 for ZmRCAβ, −logP > 4.57 for OsRCA).
Promoter-based association analysis
CTAB method was employed to extract genomic DNA (gDNA) of maize and rice leaves . The promoters of ZmRCAβ and OsRCA were amplified using specific primers (Additional file 1: Table S1) and the PCR procedure as described in our previous study . The sample PCR products were then purified and sequenced by BGI (Shanghai, China). Tested sequences were aligned using CLUSTALX version 1.83 . The identification of DNA variations among these genotypes and the candidate-gene association analysis between polymorphisms and phenotypes were performed using Tassel 5.0 . In the candidate-gene association analysis for the RCA gene promoter, all polymorphisms with a MAF ≥0.05 were considered .
Promoter activity analysis in rice callus
Maize and rice RCA promoters, which are 1700 bp and 1500 bp in length, respectively, were amplified from genomic DNA with specific primers (Additional file 1: Table S1) added with the cutting sites of HindIII and Sa1I enzymes. The amplified fragment was separated on 1% agarose gels and purified with a gel extraction kit (Axygen) according to the manufacturer’s protocol. The purified fragment was fused with the GUS reporter gene in the binary vector pCAMBIA1381Z digested with HindIII and Sa1I to construct recombinant GUS vector. Promoter-less vector pCAMBIA1381Z and CaMV 35S promoter-driving GUS vector pCAMBIA1301 were used as negative and positive controls, respectively. The recombinant GUS vector, negative control, and positive control were individually transformed into rice callus via Agrobacterium-mediated transformation as described by Hiei et al. .
GUS activity was determined in the rice positive callus according to the method as described previously  with some modifications. Plant materials were immersed in 100 mM sodium phosphate buffer (pH 7.0) containing 1% Triton X-100, 1% DMSO, 10 mM EDTA and 0.5 mg/ml 5-bromo-4-chloro-3-indolyl glucuronide and incubated overnight at 37 °C. The staining solution was then removed, and the samples were dehydrated using 75% ethanol. GUS staining was observed under an Olympus SZX12 stereomicroscope and photographed with a digital camera (CoolSNAP, RS photometrics).
Promoter activity analysis using dual luciferase assay
A dual luciferase assay vector pGreenII0800-LUC was used to analyze the activity of the target promoter. This vector contains firefly luciferase (LUC) reporter gene that can be driven by the target promoter and Renilla luciferase (REN) reporter gene driven by 35S. The purified DNA fragment of the target promoter was fused with LUC reporter gene in the vector digested with HindIII and Sa1I enzymes to construct the recombinant vector. The vector pGreenII0800-LUC without promoter insertion before LUC reporter gene was used as negative control. The recombinant and negative control vectors were individually transformed into maize and rice protoplasts or tobacco (Nicotiana benthamiana) leaves. In the protoplast transient expression experiments, the isolation of protoplasts from maize or rice green leaves, PEG-calcium transfection of plasmid DNA, and protoplast culture were performed according to standard protocols . In the leaf transient expression experiment, tobacco plants were grown until at least six leaves were available. The transfection of plasmid DNA into tobacco leaves by infiltration with Agrobacterium and transient expression assay followed standard protocols .
The ratio of LUC and REN activity (LUC/REN) was used to reflect the activity of the target promoter. The LUC/REN value was determined using the dual luciferase reporter assay system (Promega). Briefly, the transformed protoplasts or homogenate prepared from transformed tobacco leaves were centrifuged at 12000×g for 15 s at room temperature, and the supernatant was removed. Next, 100 μl of passive lysis buffer was added for further homogenization. Twenty microliters of lysate was mixed with 100 μl of LAR II, and then the LUC activity was measured using a GloMax 20/20 luminometer (Promega). Finally, 100 μl of Top & Glo reagent was added to the reaction, and the REN activity was measured.
Availability of data and materials
All data generated or analysed during this study are included in this published article (and its supplementary information files).
analysis of variance
Expression quantitative trait locus
Genome-wide association analysis
- H 2 :
the broad-sense heritability
minor allele frequencies
the mixed linear model
- pCO2 :
CO2 partial pressure
recombinant inbred line
single nucleotide polymorphisms
Salvucci ME, Ogren WL. The mechanism of rubisco activase: insights from studies of the properties and structure of the enzyme. Photosynth Res. 1996;47(1):1–11.
Streusand VJ, Portis AR. Rubisco Activase mediates ATP-dependent activation of ribulose bisphosphate carboxylase. Plant Physiol. 1987;85(1):152–4.
Portis J. Regulation of ribulose 1,5-bisphosphate carboxylase/oxygenase activity. Annu Rev Plant Biol. 1992;43(1):415–37.
Somerville CR, Portis AR, Ogren WL. A mutant of Arabidopsis thaliana which lacks activation of RuBP carboxylase in vivo. Plant Physiol. 1982;70(2):381.
Bracher A, Whitney SM, Hartl FU, Hayer-Hartl M. Biogenesis and metabolic maintenance of rubisco. Annu Rev Plant Biol. 2017;68:29–60.
Fukayama H, Mizumoto A, Ueguchi C, Katsunuma J, Morita R, Sasayama D, et al. Expression level of rubisco activase negatively correlates with rubisco content in transgenic rice. Photosynth Res. 2018;137(3):465–74.
Portis AR. Rubisco activase - Rubisco's catalytic chaperone. Photosynth Res. 2003;75(1):11–27.
Yamori W, Masumoto C, Fukayama H, Makino A. Rubisco activase is a key regulator of non-steady-state photosynthesis at any leaf temperature and, to a lesser extent, of steady-state photosynthesis at high temperature. Plant J. 2012;71(6):871–80.
Wu HR, Li LB, Jing YX, Kuang TY. Over-and anti-sense expressions of the large isoform of ribulose-1,5-bisphosphate carboxylase/oxygenase activase gene in Oryza sativa affect the photosynthetic capacity. Photosynthetica. 2007;45(2):194–201.
Fukayama H, Ueguchi C, Nishikawa K, Katoh N, Ishikawa C, Masumoto C, et al. Overexpression of rubisco activase decreases the photosynthetic CO2 assimilation rate by reducing rubisco content in rice leaves. Plant Cell Physiol. 2012;53(6):976–86.
Ristic Z, Momcilovic I, Bukovnik U, Prasad PV, Fu J, Deridder BP, et al. Rubisco activase and wheat productivity under heat-stress conditions. J Exp Bot. 2009;60(14):4003–14.
Chao M, Yin Z, Hao D, Zhang J, Song H, Ning A, et al. Variation in rubisco activase (RCAbeta) gene promoters and expression in soybean [Glycine max (L.) Merr]. J Exp Bot. 2014;65(1):47–59.
Yin Z, Meng F, Song H, Wang X, Xu X, Yu D. Expression quantitative trait loci analysis of two genes encoding rubisco activase in soybean. Plant Physiol. 2010;152(3):1625–37.
Martin M, Casano LM, Zapata JM, Guera A, del Campo EM, Schmitz-Linneweber C, et al. Role of thylakoid Ndh complex and peroxidase in the protection against photo-oxidative stress: fluorescence and enzyme activities in wild-type and ndhF-deficient tobacco. Physiol Plant. 2004;122(4):443–52.
Morales A, Ortega-Delgado ML, Molina-Galán J, Jiménez ESD. Importance of rubisco activase in maize productivity based on mass selection procedure. J Exp Bot. 1999;50(335):823–9.
Yin Z, Zhang Z, Deng D, Chao M, Gao Q, Wang Y, et al. Characterization of rubisco activase genes in maize: an alpha-isoform gene functions alongside a beta-isoform gene. Plant Physiol. 2014;164(4):2096–106.
Yang Z, Lu Q, Wen X, Chen F, Lu C. Functional analysis of the rice rubisco activase promoter in transgenic Arabidopsis. Biochem Bioph Res Co. 2012;418(3):565–70.
Zhang J, Du H, Chao M, Yin Z, Yang H, Li Y, et al. Identification of two bZIP transcription factors interacting with the promoter of soybean rubisco Activase gene (GmRCAalpha). Front Plant Sci. 2016;7:628.
Sharwood RE, Ghannoum O, Kapralov MV, Gunn LH, Whitney SM. Temperature responses of rubisco from Paniceae grasses provide opportunities for improving C3 photosynthesis. Nat Plants. 2016;2(12):16186.
Sun Q, Zhang Y, Chen B, Jia B, Zhang ZL, Cui M, et al. Expression quantitative trait loci analysis of the rubisco activase gene in maize. Photosynthetica. 2017;55(2):329–37.
Zhang Y, Cui M, Zhang J, Zhang L, Li C, Kan X, et al. Confirmation and fine mapping of a major QTL for aflatoxin resistance in maize using a combination of linkage and association mapping. Toxins. 2016;8(9):E258.
Xu J, Brautigam A, Weber AP, Zhu XG. Systems analysis of cis-regulatory motifs in C4 photosynthesis genes using maize and rice leaf transcriptomic data during a process of de-etiolation. J Exp Bot. 2016;67(17):5105–17.
To KY, Suen DF, Chen SC. Molecular characterization of ribulose-1,5-bisphosphate carboxylase/oxygenase activase in rice leaves. Planta. 1999;209(1):66–76.
Zhang ZL, Komatsu S. Molecular cloning and characterization of cDNAs encoding two isoforms of ribulose-1,5-bisphosphate carboxylase/oxygenase activase in rice (Oryza sativa L.). J Biochem. 2000;128(3):383–9.
Liu Z, Taub CC, McClung CR. Identification of an Arabidopsis thaliana ribulose-1,5-bisphosphate carboxylase/oxygenase activase (RCA) minimal promoter regulated by light and the circadian clock. Plant Physiol. 1996;112(1):43–51.
Orozco BM, Ogren WL. Localization of light-inducible and tissue-specific regions of the spinach ribulose bisphosphate carboxylase/oxygenase (rubisco) activase promoter in transgenic tobacco plants. Plant Mol Biol. 1993;23(6):1129–38.
Thumma BR, Naidu BP, Chandra A, Cameron DF, Bahnisch LM, Liu C. Identification of causal relationships among traits related to drought resistance in Stylosanthes scabra using QTL analysis. J Exp Bot. 2001;52(355):203–14.
Yan J, Warburton M, Crouch J. Association mapping for enhancing maize ( L.) genetic improvement. Crop Sci. 2011;51(2):433.
Andersen JR, Lübberstedt T. Functional markers in plants. Trends Plant Sci Trends Plant Sci. 2003;8(11):554–60.
Morley M, Molony CM, Weber TM, Devlin JL, Ewens KG, Spielman RS, et al. Genetic analysis of genome-wide variation in human gene expression. Nature. 2004;430(7001):743–7.
Ghazalpour A, Doss S, Zhang B, Wang S, Plaisier C, Castellanos R, et al. Integrating genetic and network analysis to characterize genes related to mouse weight. PLoS Genet. 2006;2(8):e130.
Gatti DM, Lu L, Williams RW, Sun W, Wright FA, Threadgill DW, et al. MicroRNA expression in the livers of inbred mice. Mutati Res-fund Mol M. 2011;714(1–2):126–33.
Grundberg E, Kwan T, Pastinen TM. Analysis of the impact of genetic variation on human gene expression. Methods Mol Biol. 2010;628:321–39.
Vuylsteke M, Van Den Daele H, Vercauteren A, Zabeau M, Kuiper M. Genetic dissection of transcriptional regulation by cDNA-AFLP. Plant J. 2006;45(3):439–46.
Hong J, Wang W, Jiang D, Hu D. The immunogold localization of rubisco and its activase in chloroplasts of barley and maize leaves. Acta Photophysiologica Sinica. 2004;30(5):561.
Wostrikoff K, Clark A, Sato S, Clemente T, Stern D. Ectopic expression of rubisco subunits in maize mesophyll cells does not overcome barriers to cell type-specific accumulation. Plant Physiol. 2012;160(1):419–32.
Friso G, Majeran W, Huang M, Sun Q, van Wijk KJ. Reconstruction of metabolic pathways, protein expression, and homeostasis machineries across maize bundle sheath and mesophyll chloroplasts: large-scale quantitative proteomics using the first maize genome assembly. Plant Physiol. 2010;152(3):1219–50.
Li P, Ponnala L, Gandotra N, Wang L, Si Y, Tausta SL, et al. The developmental dynamics of the maize leaf transcriptome. Nat Genet. 2010;42(12):1060–7.
Lescot M, Dehais P, Thijs G, Marchal K, Moreau Y, Van de Peer Y, et al. PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res. 2002;30(1):325–7.
Qu D, Song Y, Li WM, Pei XW, Wang ZX, Jia SR, et al. Isolation and characterization of the organ-specific and light-inducible promoter of the gene encoding rubisco activase in potato (Solanum tuberosum). Genet Mol Res. 2011;10(2):621–31.
Salvucci ME, CraftsBrandner SJ. Relationship between the heat tolerance of photosynthesis and the thermal stability of rubisco Activase in plants from contrasting thermal environments. Plant Physiol. 2004;134(4):1460–70.
Crafts-Brandner SJ, Salvucci ME. Rubisco activase constrains the photosynthetic potential of leaves at high temperature and CO2. Proc Natl Acad Sci U S A. 2000;97(24):13430–5.
Sánchez dJE, Medrano L, Martínezbarajas E. Rubisco activase, a possible new member of the molecular chaperone family. Biochemistry. 1995;34(9):2826–31.
Li H, Peng Z, Yang X, Wang W, Fu J, Wang J, Han Y, Chai Y, Guo T, Yang N. Genome-wide association study dissects the genetic architecture of oil biosynthesis in maize kernels. Nat Genet. 2013;45(1):43–50.
Zhao K, Tung CW, Eizenga GC, Wright MH, Ali ML, Price AH, et al. Genome-wide association mapping reveals a rich genetic architecture of complex traits in Oryza sativa. Nat Commun. 2011;2:467.
Bradbury PJ, Zhang Z, Kroon DE, Casstevens TM, Ramdoss Y, Buckler ES, Bishop PM. Genetics and population analysis TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics. 2008;23(19):2633–5.
Kumar GR, Sakthivel K, Sundaram RM, Neeraja CN, Balachandran SM, Rani NS, et al. Allele mining in crops: prospects and potentials. Biotechnol Adv. 2010;28(4):451–61.
Murray MG, Thompson WF. Rapid isolation of high molecular weight plant DNA. Nucleic Acids Res. 1980;8(19):4321–5.
Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG. The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997;25(25):4876–82.
Hiei Y, Ohta S, Komari T, Kumashiro T. Efficient transformation of rice (Oryza sativa L.) mediated by agrobacterium and sequence analysis of the boundaries of the T-DNA. Plant J. 1994;6(2):271–82.
Jefferson RA, Kavanagh TA, Bevan MW. GUS fusions: beta-glucuronidase as a sensitive and versatile gene fusion marker in higher plants. EMBO J. 1987;6(13):3901.
Yoo SD, Cho YH, Sheen J. Arabidopsis mesophyll protoplasts: a versatile cell system for transient gene expression analysis. Nat Protoc. 2007;2(7):1565.
Zahur M, Maqbool A, Irfan M, Barozai MYK, Qaiser U, Rashid B, et al. Functional analysis of cotton small heat shock protein promoter region in response to abiotic stresses in tobacco using agrobacterium-mediated transient assay. Mol Biol Rep. 2009;36(7):1915–21.
This work received grant support in part from the National Natural Science Foundation of China (31571669), the National Transgenic Research Project of China (2018ZX08009-04B), the Jiangsu Key Research and Development Program (BE2017365), Jiangsu Agriculture Science and Technology Innovation Fund (CX (18)1001), the Jiangsu Major Variety Breeding Project (PZCZ201710), the Yangzhou University Outstanding Young Teacher Training Project, and the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD). The funding agencies provided only the experimental cost and publication fee for this work. However, the experimental design and data collection and analysis were managed by the contributing authors.
Ethics approval and consent to participate
Consent for publication
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Primer pairs used in this study. (DOCX 108 kb)
Table S2. Polymorphic sites in the promoter region of ZmRCAβ and OsRCA. (XLSX 9 kb)
Table S3. Origin and haplotype of plant materials of maize and rice association populations in this study. (XLSX 11 kb)
Figure S1. Activity of different types of RCA promoters. (A) Schematic diagrams of the GUS vector. GUS, β-glucuronidase; ter, terminator. (B) Histochemical staining of rice callus transformed with recombinant plasmids containing different types of ZmRCAβ and OsRCA promoters, positive control, and negative control. The result showed that GUS expression driven by the ZmHap1 promoter was stronger than that by ZmHap3 promoter, GUS expression driven by the OsHap1 promoter was stronger than that by the OsHap5 promoter, and the GUS expression driven by rice promoters was stronger than that by the maize promoters. (DOCX 14 kb)
Table S4. SNP markers associated with the expression of OsRCA in association population. (XLSX 9 kb)
Figure S2. GWAS for the expression of OsRCA. (A) Quantile-quantile plot for the GWAS under a mixed linear model (MLM). (B) Manhattan plot for the GWAS. The red line indicates a significant association signal (−logP > 4.57). (DOCX 248 kb)
Table S5. Polymorphisms in the OsRCA promoter region significantly associated with gene expression. (XLSX 10 kb)
Figure S3. OsRCA promoter haplotypes and the average OsRCA expression for each haplotype. (A) Haplotype analysis of the OsRCA promoter regions. The dark gray-shaded cells represent the favorable alleles. (B) The average expression of OsRCA for different haplotypes. Error bars represent the standard error. The number of data points used for calculating the standard error is the product of number of rice cultivars belonging to each haplotype and biological replicates of each cultivar. Multiple comparison of phenotypic data was performed using LSD test. Different English alphabet means significant difference at P < 0.05. (DOCX 878 kb)
Table S6. SNP markers associated with the maize grain yield in the association population. (XLSX 36 kb)
Table S7. Cis-acting elements in the RCA promoter region of C3 and C4 species. (XLSX 12 kb)
About this article
Cite this article
Zhang, Y., Zhou, Y., Sun, Q. et al. Genetic determinants controlling maize rubisco activase gene expression and a comparison with rice counterparts. BMC Plant Biol 19, 351 (2019). https://doi.org/10.1186/s12870-019-1965-x
- Rubisco activase (RCA)
- Quantitative trait locus (eQTL)