Next Article in Journal
NO and GSH Alleviate the Inhibition of Low-Temperature Stress on Cowpea Seedlings
Next Article in Special Issue
The Shift in Key Functional Traits Caused by Precipitation under Nitrogen and Phosphorus Deposition Drives Biomass Change in Leymus chinensis
Previous Article in Journal
Evaluation of Amino Acid Profiles of Rice Genotypes under Different Salt Stress Conditions
Previous Article in Special Issue
Genome-Wide Identification and Evolutionary Analysis of Gossypium YTH Domain-Containing RNA-Binding Protein Family and the Role of GhYTH8 in Response to Drought Stress
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Plants
Volume 12
Issue 6
10.3390/plants12061309
plants-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Transcriptomic Analysis Reveals Panicle Heterosis in an Elite Hybrid Rice ZZY10 and Its Parental Lines

by
Zhengzheng Zhong
1,†,
Yawen Wu
1,†,
Peng Zhang
1,
Guocheng Hu
1,
Dong Fu
2,
Guoping Yu
3 and
Hanhua Tong
1,*
1
State Key Laboratory of Rice Biology and Breeding, China National Rice Research Institute, Hangzhou 311400, China
2
State Key Laboratory of Biocatalysis and Enzyme Engineering, School of Life Sciences, Hubei University, Wuhan 430062, China
3
National Nanfan Research Institute, Chinese Academy of Agricultural Sciences, Sanya 572025, China
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Plants 2023, 12(6), 1309; https://doi.org/10.3390/plants12061309
Submission received: 4 December 2022 / Revised: 10 March 2023 / Accepted: 12 March 2023 / Published: 14 March 2023
(This article belongs to the Special Issue The Trade-Offs between Growth and Development and Stress in Plants)
Browse Figures
Versions Notes

Abstract

:
Heterosis is the phenomenon in which some hybrid traits are superior to those of their parents. Most studies have analyzed the heterosis of agronomic traits of crops; however, heterosis of the panicles can improve yield and is important for crop breeding. Therefore, a systematic study of panicle heterosis is needed, especially during the reproductive stage. RNA sequencing (RNA Seq) and transcriptome analysis are suitable for further study of heterosis. Using the Illumina Nova Seq platform, the transcriptome of ZhongZheYou 10 (ZZY10), an elite rice hybrid, the maintainer line ZhongZhe B (ZZB), and the restorer line Z7-10 were analyzed at the heading date in Hangzhou, 2022. 581 million high-quality short reads were obtained by sequencing and were aligned against the Nipponbare reference genome. A total of 9000 differential expression genes were found between the hybrids and their parents (DGHP). Of the DGHP, 60.71% were up-regulated and 39.29% were down-regulated in the hybrid. Comparative transcriptome analysis revealed that 5235 and 3765 DGHP were between ZZY10 and ZhongZhe B and between ZZY10 and Z7-10, respectively. This result is consistent with the transcriptome profile of ZZY10 and was similar to Z7-10. The expression patterns of DGHP mainly exhibited over-dominance, under-dominance, and additivity. Among the DGHP-involved GO terms, pathways such as photosynthesis, DNA integration, cell wall modification, thylakoid, and photosystem were significant. 21 DGHP, which were involved in photosynthesis, and 17 random DGHP were selected for qRT-PCR validation. The up-regulated PsbQ and down-regulated subunits of PSI and PSII and photosynthetic electron transport in the photosynthesis pathway were observed in our study. Extensive transcriptome data were obtained by RNA-Seq, providing a comprehensive overview of panicle transcriptomes at the heading stage in a heterotic hybrid.

1. Introduction

Rice (Oryza sativa L.) is a model cereal crop and one of the main global food crops. Heterosis is a phenomenon in which hybrids show superior performance to one or both parental phenotypes. In the past fifty years, heterosis has significantly increased the yield of global rice crops [1,2]. Hybrid rice varieties, such as Shanyou63 and Liangyoupei9, are widely cultivated in China. One of the best approaches for achieving super-high yields of rice is through inter-subspecific hybrids [1]. The yield of hybrid varieties of indica and japonica increased by approximately 25% [3]. Zhongzheyou 10 (ZZY10), which shows very strong yield heterosis in some agronomic traits and is widely adapted, is an elite rice hybrid released in 2016. Hybrid ZZY10 was developed using Zhongzhe A (ZZA) and Z7-10 as its two parental lines. The Cytoplasmic Male Sterile (CMS) line ZZA, which is widely used as the parent, has created several excellent hybrids cultivated in China. The most attractive feature of ZZY10 is strong yield heterosis, including panicle length (PL), 1000 grain weight (1000-GW), number of secondary branches (SBN), tiller number (TN), plant height (PH), and grain width (GW), which is termed optimal heterosis and is extremely required in breeding. This hybrid ZZY10 can be systematically used for studying the molecular mechanism of heterosis.
Improving photosynthesis is one of the most effective ways of increasing crop yield [4,5]. However, most previous studies analyzed improvements in leaf photosynthesis, but few have assessed non-foliar green tissue [6,7]. Panicles of rice and ears of wheat have strong photosynthetic capacities. Photosystem I (PSI) and photosystem II (PSII) protein complexes play a crucial role during photosynthesis [8,9]. All protein factors involved in the kinetic reactions of PSI and PSII complexes must be characterized in crops to improve photosynthetic efficiency and productivity [10,11]. It is speculated that improving ears photosynthesis in wheat will help increase the yield by 12–60% [12,13]. Therefore, improving photosynthetic efficiency in the ear is crucial for yield [14].
RNA sequencing-based transcriptomic analysis is a valuable strategy for exploring different pathways and regulatory networks and can help better understand biological processes [15,16,17]. Using transcriptome data as an informative platform, we can investigate the molecular mechanism of heterosis and the complex gene expression patterns associated with heterosis in crops (like rice, maize, and wheat) [18,19,20,21,22,23]. The classical genetic pattern of heterosis results in dominance, over-dominance, additivity, and epistasis [23,24,25,26].
In the present study, we concentrated on transcriptomic data from the Maintainer Line Zhongzhe B (ZZB), Z7-10, and ZZY10 to analyze the molecular aspects of heterosis. The analysis provides some explanations for the optimal heterosis of increased panicle length (PL) in hybrid ZZY10. The comparative transcriptomes of ZZB, Z7-10, and the hybrid ZZY10 were comprehensively analyzed to provide insights into genetic hypotheses, expression divergence, and the regulation pathway of genes between the parents and the hybrid.

2. Results

2.1. Characterization of ZZY10 and Its Two Parental Lines

ZZY10 is a rice hybrid produced by crossing the maternal line ZZB and the paternal line Z7-10. As suggested by Pickett [27], many yield-related traits of ZZY10 could show heterosis (Table 1). Interestingly, the panicle length of ZZY10 exceeded the panicle length of both parents (Table 2). We investigated the gene expression and phenotypes of ZZY10, ZZB, and Z7-10 using RNA-seq. The yield-related traits of ZZY10 showed heterosis (Table 1), as suggested by Pickett [27]. ZZY10 also has a longer panicle than its parents (Table 2). According to the phenotype of F1 hybrid ZZY10, mid-parent heterosis (MPH) and high-parent heterosis (HPH) were calculated to measure heterosis. The values of MPH and HPH are shown in Table 2. At the heading date, the MPH and HPH of the panicle length were 10.79% and −0.77%, respectively. In comparison, the MPH and HPH of the panicle length were 1.44% and −3.05% at the ripening stage, respectively (Table 2). The difference in panicle length’s MPH was significant at the heading stage.

2.2. Mapping Short Reads to the Rice Genome

To perform the transcriptomes of the panicle tissue in hybrid ZZY10 and the parents at the heading date, cDNA libraries were constructed, and RNA-Seq was performed on the tissue of rice panicles using the Illumina Nova-Seq platform. We obtained more than 581 million short reads for rice varieties with the three biological replicates, which aligned with gene expression levels at 0.81 < R2 < 0.98 (Figure S1). The raw data have been uploaded into the NCBI database SRA (PRJNA934656). We aligned the short reads with the Nipponbare reference genome (NCBI build 4.0), and our results showed that exonic regions covered 89.54–91.76% of reads, intronic regions covered 2.87–3.34%, and intergenic regions covered 5.30–7.17% (Table S1). In our RNA-Seq analysis, 49,044 annotated transcripts were identified in the Nipponbare reference genome (Table S2). From these annotated transcripts, we finally confirmed 9000 differentially expressed genes from the heading stage panicles (Table S3). We identified a total of 18,503 novel transcript isoforms in ZZB, Z7, and ZZY10 (Table S4).
GOSlim term analysis of the novel transcripts revealed that novel genes involved in biological processes (≥6%), followed by metabolic processes (≥4%), organic substance metabolic processes (≥3%), cellular processes (≥3%), primary metabolic processes (≥3%), cellular metabolic processes (≥3%), macromolecule metabolic processes (≥2%), cellular macromolecule metabolic processes (≥2%), macromolecule metabolic processes (≥2%), and single-organism processes (≥2%) represented approximately 30% of the novel transcripts (Figure S2). In the molecular function category, more than 40% of the novel transcripts were classified as molecular function (≥11%), binding (≥7%), catalytic activity (≥6%), heterocyclic compound binding (≥4%), and organic cyclic compound binding (≥4%) (Figure S2). A total of 40% of novel transcripts belonged to cellular components (≥12%), membrane (≥7%), cell (≥6%), cell parts (≥6%), intracellular (≥5%), and intracellular parts (≥5%) of the cellular component GOSlim categories (Figure S2). The number of novel transcripts that were categorized in biological processes was the highest compared to the others.

2.3. Differentially Expressed Genes (DEGs) Identified by RNA-Seq

FPKM (Fragments Per Kilo bases per Millionreads) involves the influence of sequencing depth and gene length on read count. Since three biological replicates of the experiments were performed, we used DESeq2 software to analyze the differentially expressed genes and used a p value less than or equal to 0.05 as the standard of a significant difference. Based on these criteria, 9000 reliable DEGs were identified among the three genotypes. DEGs between the hybrid ZZY10 and its two parents, ZZB and Z7-10, were designated as DGHP. Considering that different expression between the hybrid and its parents is the underlying basis of their different phenotypes, DGHP could be relevant to the phenotypes of heterosis [23]. Comparative transcriptome analysis revealed that a total of 5235 transcripts, including 3438 up-regulated and 1797 down-regulated transcripts, were between ZZY10 and ZZB (Figure 1A). A total of 3765 transcripts (2026 and 1739 transcripts commonly up-regulated and down-regulated, respectively) were between ZZY10 and Z7-10 and were differentially expressed (Figure 1A). In addition, 1468 transcripts between ZZY10 vs. Z7-10 and ZZY10 vs. ZZB were commonly regulated (Figure 1B). More transcripts showed differences between ZZY10 and ZZB. These results showed that the gene expression profile of ZZY10 is more similar to that of Z7-10 than that of ZZB. Principal component analysis demonstrated that ZZY10 is much closer to Z7-10 (Figure S3). Moreover, the transcriptome profile of ZZY10 was similar to Z7-10 from the hierarchical cluster analysis (Figure S4). These results consistently demonstrated that ZZY10 is closer to Z7-10 than ZZB.
When comparing the expression levels of all the DEGs in ZZB, Z7-10, and ZZY10, we found a composite expression pattern, including five sub-patterns: over-dominance, additivity, low-parent dominance, high-parent dominance, and under-dominance, which suggested multiple modes of action (Table 3). The percentages of up-regulated, additivity, and down-regulated genes in DGHP are 35.6%, 38.8%, and 25.6%, respectively. As shown in Figure 2 and Table 3, the distributions of DGHP are over-dominance (28.9%), high-parent dominance (6.7%), additivity (38.8%), low-parent dominance (2.1%), and under-dominance (23.5%). First, the majority of DGHP exhibits additivity. Second, the percentages of over-dominant DGHP and under-dominant DGHP are 28.9% and 23.5%, respectively (Figure 2 and Table 3). Finally, the percentages of additivity, over-dominance, and under-dominance account for 91.2% of all DGHP (Figure 2 and Table 3).

2.4. Gene Ontology (GO) Analysis

GO enrichment analysis was performed for the functional classification of DGHP to better understand their roles in different biological activities. In total, 2845 of the 9000 DGHP between the hybrid and its parents were assigned to at least one term among the categories of cellular component (CC), biological process (BP), and molecular function (MP). Transcripts were further divided into 36 subcategories to provide an overview of ontology content (Figure 3). In the MP category, the most represented groups were binding and catalytic activity. In the CC category, the membrane and cell were prominently represented, while metabolic processes and cellular processes dominated the BP category.
We further identified that DGHP were significantly enriched in GO terms (FDR-corrected p-values < 0.05) (Table 4). GO terms such as photosynthesis, DNA integration, cell wall modification, thylakoid, photosystem, and photosynthetic membrane were significant, suggesting that these significant pathways may be related to panicle length.

2.5. Kyoto Encyclopedia of Genes and Genomes (KEGG) Enrichment Analysis

To identify metabolic pathways in which DEGs between the hybrid and its parents were enriched, we performed a KEGG analysis. Our results demonstrated that 584 of the 9000 DGHP were classified into eight functional categories and were mainly involved in the pathways of brite hierarchy, metabolism, genetic information processing, cellular processes, environmental information processing, and organismal systems (Figure 4). These eight categories were further categorized into 97 subcategories corresponding to their functions, and the significant DGHP-involved KEGG Orthology (KO) terms were identified, including photosynthesis, plant-pathogen interaction, alpha-linolenic acid metabolism, photosynthesis-antenna proteins, and protein processing in the endoplasmic reticulum. These annotations may supply a useful resource for investigating specific processes, pathways, and functions related to heterosis of panicle length.

2.6. Quantitative Real-Time PCR (qRT-PCR)

To prove the reliability of the transcriptome results, we selected a subset of 38 DGHP that was used for qRT-PCR validation. We randomly selected 17 DGHP, and the other DGHP were pointed at photosynthesis. Primer sequences are displayed in Table S5. The results of the transcriptome and quantitative PCR were analyzed. Our comparison results showed that the expression of transcripts obtained from qRT-PCR and RNA-Seq had consistent trends (R2 = 0.8579) (Figure 5).

3. Discussion

Heterosis, which has been widely applied in crop breeding, plays an important role in crop production, though its functional mechanisms must still be clarified. DGHP that causes different phenotypes could be related to heterosis [23,28,29,30,31]. In this study, our RNA-Seq analysis demonstrated that 581 million high-quality short reads were generated from panicle tissues of ZZY10 and its parents at the heading stage, and 49,044 annotated transcripts were aligned in the Nipponbare reference genome (Tables S1 and S2). Transcriptional profiling of F1 hybrid rice ZZY10 and its parents by serial analysis of gene expression revealed 9000 DGHP, among which DGHP were found significantly enriched in pathways such as photosynthesis, DNA integration, cell wall modification, thylakoid, and photosystem. Two of the key genes involved in the photosynthesis pathway exhibited up-regulated expression in F1 hybrid rice. Additionally, 38 DGHP of RNA-Seq were verified by qRT-PCR.

3.1. Comparative Analysis of the Annotated DGHP

A subset of DGHP was observed in our comparative transcriptome analysis at the heading stage. Among the significant DGHP-involved GO terms, the top five significant GO terms were photosynthesis, DNA integration, cell wall modification, thylakoid, and photosystem, respectively. In the two significant pathways, 21 DGHP were detected in the photosynthesis pathway and 16 DGHP in the photosystem pathway. 2 DGHP both in the photosynthesis pathway and in the photosystem pathway were up-regulated in the F1 hybrid. In this study, the panicle length of ZZY10 showed mid-parent heterosis at the heading stage. Interestingly, transcriptional expression levels of PsbP (Os12g0564400) in ZZY10 are higher than Z7-10 and lower than ZZB, which is consistent with the phenotype of panicle length in ZZY10. The PsbP proteins are a key component of photosystem II in eukaryotic photosynthetic organisms [32]. The functional characterization of PsbP homologs revealed that Arabidopsis mutants play important roles in maintaining photosynthetic electron transfer [33]. It can be assumed that DGHP involved in the photosynthesis pathway may regulate the heterosis of panicle length.

3.2. Photosynthesis Regulation in Heterosis

Photosynthetic electron transport complex contains two light-energy-driven photosystems (PSI and PSII), the ATP synthase, and the cytochrome (Cyt) b6f [32]. By comparative analysis of transcriptomes, significantly different photosynthesis was observed in the panicle tissues of the hybrid and the parents (Figure 6). Many DGHP-encoded proteins involved in photosynthesis have complex organizations (Figure 6). Oxygenic photosynthesis generates molecular oxygen via the water-oxidizing reaction under the catalysis of PSII, which consists of intrinsic and extrinsic subunits of membrane proteins [34,35]. In this regard, transcripts encoding PsbO, PsbP, PsbR, PsbS, PsbW, PsbY, and Psb27, which are subunits of PSII (Os01g0501800, PsbO; Os07g0141400, PsbP; Os07g0147500 and Os08g0200300, PsbR; Os01g0869800 and Os04g0690800, PsbS; Os01g0773700, PsbW; Os08g0119800, PsbY and Os03g0333400, Psb27), were down-regulated in ZZY10 compared with the parental lines. However, Os02g0578400 and Os07g0105600, which encode PsbQ, were more up-regulated in ZZY10 than the parental lines. The high level of PsbQ expression in ZZY10 may play a role in photosynthesis and promote panicle length.
In the linear electron transport (LET), the water-extracted electrons in PSII are delivered to PSI via Cyt b6f and eventually generate NADPH. Transcriptional expression levels of PsaD (Os08g0560900), PsaG (Os09g0481200), PsaH (Os05g0560000), PsaK (Os07g0148900), PsaL (Os12g0420400), PsaN (Os12g0189400), and PsaO (Os04g0414700), which are subunits involved in photosystem I (PSI), were down-regulated in ZZY10. Two DGHP (Os06g0101600 and Os03g0835900) that encode PetE and PetF in the photosynthetic electron transport showed down-regulated expression in ZZY10. The up-regulated PsbQ, down-regulated subunits, and photosynthetic electron transport in the photosynthesis pathway observed in our study may co-regulate the panicle length of ZZY10. It can be assumed that photosynthesis is involved in the heterosis of panicle length, but a more accurate characterization will be done in the future to prove it.

3.3. Implications for the Mechanism of Heterosis

The detection of multiple expression patterns of DEGs, such as under-dominance, additivity, and over-dominance (Figure 2), suggested that heterosis is a complex issue. It is difficult for a single hypothesis to verify its molecular basis since additivity hypotheses, dominance, or over-dominance could all contribute to rice heterosis [36,37]. Among the 9000 DGHP, 3536 are down-regulated and 5464 are up-regulated. In total, 3492 (38.8%), 2601 (28.9%), and 2115 (23.5%) DGHP showed additivity, over-dominance, and under-dominance, respectively. In this study, we found that most of the DGHP displayed additivity. The number of over-dominant DGHP was almost equal to the number of under-dominant DGHP. Since ZZB and Z7-10 have some excellent minor and dominant genes, respectively, ZZY10 is a large panicle hybrid rice variety that has excellent minor and dominant genes. More than half of the DEGs differ between ZZY10 and ZZB (58.2%), while less than half show differences between ZZY10 and Z7-10 (41.8%). This result could explain why the panicle phenotype of ZZY10 resembles Z7-10 more than ZZB. Therefore, the hybrid vigor of ZZY10 could be produced by beneficial allele combinations of its parents, which could complement the relatively weaker alleles, or it could also come from the accumulation of multiple minor genes between two parents.

4. Materials and Methods

4.1. Plants and Growing Conditions

ZZY10, ZZB, and Z7-10, which are hybrid rice, a maintainer line, and a restorer line, respectively, obtained in our breeding used in this study. Before sowing, rice seeds of ZZY10, ZZB, and Z7-10 were placed in the electrothermal dryer at 35 °C for 5 days and soaked in distilled water (37 °C, 2 d) to germinate before being finally sown in the fields of Fuyang, Hangzhou, and Zhejiang on May 20, 2022 (30°08 N, 119°92 E, and an altitude of 10 m). Thirty days later, 40 day seedlings of ZZY10, ZZB, and Z7-10 were transplanted into the paddy field. Panicles were collected with nine replicates to evaluate heterosis. Finally, three panicles of each genotype at the heading date were stored in −80 °C refrigerators until RNA-Seq analysis.

4.2. Agronomic Trait Evaluation

In this study, several agronomic traits, including panicle length, primary branch number, plant height, secondary branch number, 1000-grain weight, tiller number, grain length, and grain width, were measured according to the protocols previously described [38]. HPH and MPH were calculated as follows:
HPH = (F1 − HP)/HP × 100%
and
MPH = (F1 − MP)/MP × 100%
F1, HP, and MP stand for the performance of the first hybrid generation, the best parent’s performance, and the average performance of the parents, respectively.

4.3. Transcriptome Analysis

Transcriptome data was analyzed through reference genome mapping of clean reads by Bowtie2 and Tophat2 tools [39]. DESeq2 software was used to analyze DEGs. DEGs were identified using the difference significance criteria of FDR ≤ 0.05 and an estimated absolute log2 (FC) ≥ 1. GO function enrichment analysis and KEGG pathway enrichment analysis for DEGs were performed using the clusterProfiler tool. Cluster analysis was performed by Cluster 3.0 software. Annotations for differentially expressed genes were obtained from the Rice Genome Annotation Project on 8 October 2022 (http://rice.plantbiology.msu.edu/ accessed on 1 January 2020).

4.4. qRT-PCR

Total RNA was extracted from the panicle tissues of ZZY10, ZZB, and Z7-10. RNA of the three samples was reverse transcribed into double-stranded cDNA with a PrimeScript™ RT reagent kit with a gDNA Eraser kit (TaKaRa Bio, Ohtsu, Japan). SYBR-based qRT-PCR (Takara, Japan) was carried out with the SYBR® Premix Ex Taq™ II kit (Takara Bio, Ohtsu, Japan) and a LightCycler 480 system. The reaction procedure was as follows: 95 °C for 30 s, 95 °C for 30 s, and 60 °C for 30 s (40 cycles). The results of the relative quantification experiments were analyzed by the ∆Ct method. The rice gene ACTIN1 was used as an internal control.

5. Conclusions

We systematically investigated the global transcriptomes of the heading-stage panicles from the hybrid rice ZZY10 and its parents using RNA-seq. The results of GO and KEGG analysis and the heterosis functional mechanism exploration create a valuable resource for rice breeding and reveal putative transcripts that could contribute to panicle length. Expression changes of candidate genes could provide helpful information for future studies assessing possible mechanisms of panicle heterosis.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/plants12061309/s1, Figure S1: Pearson correlation of gene expression from biological replicates in ZZB, Z7-10, and ZZY10; Figure S2: GOSlim term analysis in BP, MF, and CC; Figure S3: PCA of the hybrid ZZY10 and its parents; Figure S4: The cluster analysis of all differentially expressed transcripts in ZZB, Z7-10, and ZZY10; Table S1: Summary of filtered, mapped, and reference-based assembled read data; Table S2: RPKM of annotated transcripts; Table S3: List of differentially expressed genes; Table s4: List of novel transcripts; Table S5: Primers for qPCR.

Author Contributions

Z.Z. and H.T.: Conceptualization, methodology, writing and funding acquisition. Y.W.: Data curation, software. P.Z. and G.Y.: Review and editing draft. G.H. and D.F.: Investigation. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Natural Science Foundation of Zhejiang province (grant no. LY20C130010), the National Natural Science Foundation of China (grant no. 31872862 and 32172077), the Nanfan special project, CAAS (grant no. YBXM06), and the Open Program of the State Key Laboratory of Rice Biology (grant no. 20210206).

Data Availability Statement

Not applicable.

Acknowledgments

We are grateful to Weixun Wu (China National Rice Research Institute) for editing the final manuscript and Yue Feng (China National Rice Research Institute) for the calculation method of MPH and HPH.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cheng, S.H.; Zhuang, J.Y.; Fan, Y.Y.; Du, J.H.; Cao, L.Y. Progress in research and development on hybrid rice: A super-domesticate in China. Ann. Bot. 2007, 100, 959–966. [Google Scholar] [CrossRef] [PubMed]
  2. Ma, G.H.; Yuan, L.P. Hybrid rice achievements, development and prospect in China. J. Integr. Agric. 2015, 14, 197–205. [Google Scholar] [CrossRef]
  3. Khush, G.S. Breaking the yield frontier of rice. Geojournal 1995, 35, 329–332. [Google Scholar] [CrossRef]
  4. Lawson, T.; Flexas, J. Fuelling life: Recent advances in photosynthesis research. Plant J. 2020, 101, 753–755. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Long, S.P.; Marshall-Colon, A.; Zhu, X.G. Meeting the Global Food Demand of the Future by Engineering Crop Photosynthesis and Yield Potential. Cell 2015, 161, 56–66. [Google Scholar] [CrossRef] [Green Version]
  6. Aschan, G.; Pfanz, H. Non-foliar photosynthesis and its contribution to the overall carbon balance of plants. Ann. Des Sci. 2003, 25, 3–23. [Google Scholar] [CrossRef] [Green Version]
  7. Simkin, A.J.; Faralli, M.; Ramamoorthy, S.; Lawson, T. Photosynthesis in non-foliar tissues: Implications for yield. Plant J. 2020, 101, 1001–1015. [Google Scholar] [CrossRef]
  8. Nelson, N.; Ben-Shem, A. The complex architecture of oxygenic photosynthesis. Nat. Rev. Mol. Cell Biol. 2004, 5, 971–982. [Google Scholar] [CrossRef]
  9. Wu, G.X.; Ma, L.; Yuan, C.; Dai, J.H.; Luo, L.; Poudyal, R.S.; Sayre, R.T.; Lee, C.H. Formation of light-harvesting complex II aggregates from LHCII–PSI–LHCI complexes in rice plants under high light. J. Exp. Bot. 2021, 72, 4938–4948. [Google Scholar] [CrossRef]
  10. Gao, P.; Xia, H.Q.; Li, Q.; Li, Z.Z.; Zhai, C.; Weng, L.; Mi, H.L.; Yan, S.; Datla, R.; Wang, H.; et al. PALE-GREEN LEAF 1, a rice cpSRP54 protein, is essential for the assembly of the PSI-LHCI supercomplex. Plant Direct 2022, 6, e451. [Google Scholar] [CrossRef]
  11. Reddy, M.R.; Mangrauthia, S.K.; Reddy, S.V.; Manimaran, P.; Yugandhar, P.; Naresh Babu, P.; Vishnukiran, T.; Subrahmanyam, D.; Sundaram, R.M.; Balachandran, S.M. PAP90, a novel rice protein plays a critical role in regulation of D1 protein stability of PSII. J. Adv. Res. 2021, 30, 197–211. [Google Scholar] [CrossRef] [PubMed]
  12. Maydup, M.L.; Antonietta, M.; Guiamet, J.J.; Graciano, C.; Lopez, J.R.; Tambussi, E.A. The contribution of ear photosynthesis to grain filling in bread wheat (Triticum aestivum L.). Field Crops Res. 2010, 119, 48–58. [Google Scholar] [CrossRef]
  13. Sanchez-Bragado, R.; Molero, G.; Reynolds, M.P.; Araus, J.L. Photosynthetic contribution of the ear to grain filling in wheat: A comparison of different methodologies for evaluation. J. Exp. Bot. 2016, 67, 2787–2798. [Google Scholar] [CrossRef]
  14. Maydup, M.L.; Antonietta, M.; Graciano, C.; Guiamet, J.J.; Tambussi, E.A. The contribution of the awns of bread wheat (Triticum aestivum L.) to grain filling: Responses to water deficit and the effects of awns on ear temperature and hydraulic conductance. Field Crops Res. 2014, 167, 102–111. [Google Scholar] [CrossRef]
  15. Coolen, S.; Proietti, S.; Hickman, R.; Davila Olivas, N.H.; Huang, P.P.; Van Verk, M.C.; Van Pelt, J.A.; Wittenberg, A.H.; De Vos, M.; Prins, M.; et al. Transcriptome dynamics of Arabidopsis during sequential biotic and abiotic stresses. Plant J. 2016, 86, 249–267. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Digel, B.; Pankin, A.; Korff, M.V. Global Transcriptome Profiling of Developing Leaf and Shoot Apices Reveals Distinct Genetic and Environmental Control of Floral Transition and Inflorescence Development in Barley. Plant Cell 2015, 27, 2318–2334. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Vogel, J.T.; Zarka, D.G.; Van Buskirk, H.A.; Fowler, S.G.; Thomashow, M.F. Roles of the CBF2 and ZAT12 transcription factors in configuring the low temperature transcriptome of Arabidopsis. Plant J. 2005, 41, 195–211. [Google Scholar] [CrossRef]
  18. Ding, H.P.; Cheng, Q.; Luo, X.R.; Li, L.J.; Chen, Z.; Liu, H.J.; Gao, J.; Lin, H.J.; Shen, Y.O.; Zhao, M.J.; et al. Heterosis in Early Maize Ear Inflorescence Development: A Genome-Wide Transcription Analysis for Two Maize Inbred Lines and Their Hybrid. Int. J. Mol. Sci. 2014, 15, 13892–13915. [Google Scholar] [CrossRef] [Green Version]
  19. Li, A.; Liu, D.C.; Wu, J.; Zhao, X.B.; Hao, M.; Geng, S.F.; Yan, J.; Jiang, X.X.; Zhang, L.Q.; Wu, J.Y.; et al. mRNA and Small RNA Transcriptomes Reveal Insights into Dynamic Homoeolog Regulation of Allopolyploid Heterosis in Nascent Hexaploid Wheat. Plant Cell 2014, 26, 1878–1900. [Google Scholar] [CrossRef] [Green Version]
  20. Shankar, R.; Bhattacharjee, A.; Jain, M. Transcriptome analysis in different rice cultivars provides novel insights into desiccation and salinity stress responses. Sci. Rep. 2016, 6, 23719. [Google Scholar] [CrossRef] [Green Version]
  21. Song, G.Y.; Guo, Z.B.; Liu, Z.W.; Qu, X.F.; Jiang, D.M.; Wang, W.; Zhu, Y.G.; Yang, D.C. The phenotypic predisposition of the parent in F1hybrid is correlated with transcriptome preference of the positive general combining ability parent. BMC Genom. 2014, 15, 297. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Swanson-Wagner, R.A.; Jia, Y.; DeCook, R.; Borsuk, L.A.; Nettleton, D.; Schnable, P.S. All possible modes of gene action are observed in a global comparison of gene expression in a maize F-1 hybrid and its inbred parents. Proc. Natl. Acad. Sci. USA 2006, 103, 6805–6810. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Wei, G.; Tao, Y.; Liu, G.Z.; Chen, C.; Luo, R.Y.; Xia, H.; Gan, Q.; Zeng, H.P.; Lu, Z.K.; Han, Y.N.; et al. A transcriptomic analysis of superhybrid rice LYP9 and its parents. Proc. Natl. Acad. Sci. USA 2009, 106, 7695–7701. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Li, Z.K.; Luo, L.J.; Mei, H.W.; Wang, D.L.; Shu, Q.Y.; Tabien, R.; Zhong, D.B.; Ying, C.S.; Stansel, J.W.; Khush, G.S.; et al. Overdominant epistatic loci are the primary genetic basis of inbreeding depression and heterosis in rice. I. Biomass and grain yield. Genetics 2001, 158, 1737–1753. [Google Scholar] [CrossRef] [PubMed]
  25. Shull, G.H. The Composition of a Field of Maize. J. Hered. 1908, 4, 296–301. [Google Scholar] [CrossRef] [Green Version]
  26. Yu, R.; Xin, C. Gene location, chromosomal marker and inheritance of coat colour in pigs. Acta Vet. Zootech. Sin. 1997, 28, 97–102. [Google Scholar]
  27. Pickett, A.A. Hybrid wheat: Results and problems. J. Agric. Sci. 1993, 121, 1–259. [Google Scholar] [CrossRef]
  28. Ge, X.M.; Chen, W.H.; Song, S.H.; Wang, W.W.; Hu, S.N.; Yu, J. Transcriptomic profiling of mature embryo from an elite super-hybrid rice LYP9 and its parental lines. BMC Plant Biol. 2008, 8, 114. [Google Scholar] [CrossRef] [Green Version]
  29. He, G.M.; Zhu, X.P.; Elling, A.A.; Chen, L.B.; Wang, X.F.; Guo, L.; Liang, M.Z.; He, H.; Zhang, H.Y.; Chen, F.F.; et al. Global Epigenetic and Transcriptional Trends among Two Rice Subspecies and Their Reciprocal Hybrids. Plant Cell 2010, 22, 17–33. [Google Scholar] [CrossRef] [Green Version]
  30. Song, G.S.; Zhai, H.L.; Peng, Y.G.; Zhang, L.; Wei, G.; Chen, X.Y.; Xiao, Y.G.; Wang, L.L.; Chen, Y.J.; Wu, B.; et al. Comparative Transcriptional Profiling and Preliminary Study on Heterosis Mechanism of Super-Hybrid Rice. Mol. Plant 2010, 3, 1012–1025. [Google Scholar] [CrossRef]
  31. Zhang, H.Y.; He, H.; Chen, L.B.; Li, L.; Liang, M.Z.; Wang, X.F.; Liu, X.G.; He, G.M.; Chen, R.S.; Ma, L.G.; et al. A genome-wide transcription analysis reveals a close correlation of promoter INDEL polymorphism and heterotic gene expression in rice hybrids. Mol. Plant 2008, 1, 720–731. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Ifuku, K. The PsbP and PsbQ family proteins in the photosynthetic machinery of chloroplasts. Plant Physiol. Biochem. 2014, 81, 108–114. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Ifuku, K.; Ishihara, S.; Shimamoto, R.; Ido, K.; Fumihiko, S. Structure, function, and evolution of the PsbP protein family in higher plants. Photosynth. Res. 2008, 98, 427–437. [Google Scholar] [CrossRef] [PubMed]
  34. Pagliano, C.; Saracco, G.; Barber, J. Structural, functional and auxiliary proteins of photosystem II. Photosynth. Res. 2013, 116, 167–188. [Google Scholar] [CrossRef]
  35. Shen, J.R. The Structure of Photosystem II and the Mechanism of Water Oxidation in Photosynthesis. Annu. Rev. Plant Biol. 2015, 66, 23–48. [Google Scholar] [CrossRef] [Green Version]
  36. Bruce, A.B. The mendelian theory of heredity and the augmentation of vigor. Science 1910, 32, 627–628. [Google Scholar] [CrossRef]
  37. Yu, S.B.; Li, J.X.; Xu, C.G.; Tan, Y.F.; Gao, Y.J.; Li, X.H.; Zhang, Q.; Maroof, M.A.S. Importance of epistasis as the genetic basis of heterosis in an elite rice hybrid. Proc. Natl. Acad. Sci. USA 1997, 94, 9226–9231. [Google Scholar] [CrossRef] [Green Version]
  38. Guo, H.; Liu, X. The Research on Autotetraploid Rice; Guangzhou South China University of Technology Press: Guangzhou, China, 2014; ISBN 9787562341208. [Google Scholar]
  39. Kim, D.; Pertea, G.; Trapnell, C.; Pimentel, H.; Kelley, R.; Salzberg, S.L. TopHat2: Accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013, 14, R36. [Google Scholar] [CrossRef] [Green Version]
Plants 12 01309 g001 550
Figure 1. Venn diagram of DGHP. (A) Number of up-regulated and down-regulated genes expressed in two cross combinations; (B) Number of genes expressed in two cross combinations.
Figure 1. Venn diagram of DGHP. (A) Number of up-regulated and down-regulated genes expressed in two cross combinations; (B) Number of genes expressed in two cross combinations.
Plants 12 01309 g001
Plants 12 01309 g002 550
Figure 2. Breakdown of the DGHP according to the expression pattern.
Figure 2. Breakdown of the DGHP according to the expression pattern.
Plants 12 01309 g002
Plants 12 01309 g003 550
Figure 3. Comparison of GO term classifications of DGHP at the heading stage.
Figure 3. Comparison of GO term classifications of DGHP at the heading stage.
Plants 12 01309 g003
Plants 12 01309 g004 550
Figure 4. (A) KEGG pathway assignments into eight functional categories. (B) The significant enrichment list of KEGG pathways in subcategories. The significantly enriched pathway (p < 0.05) and the number of enriched transcripts in each pathway are shown. The x-axis in Figure B indicates the number of differentially expressed genes of these terms that are significantly enriched.
Figure 4. (A) KEGG pathway assignments into eight functional categories. (B) The significant enrichment list of KEGG pathways in subcategories. The significantly enriched pathway (p < 0.05) and the number of enriched transcripts in each pathway are shown. The x-axis in Figure B indicates the number of differentially expressed genes of these terms that are significantly enriched.
Plants 12 01309 g004
Plants 12 01309 g005 550
Figure 5. Association analysis on 38 selected transcripts between RNA-Seq and qRT-PCR.
Figure 5. Association analysis on 38 selected transcripts between RNA-Seq and qRT-PCR.
Plants 12 01309 g005
Plants 12 01309 g006 550
Figure 6. KEGG enrichment pathway of differential expression genes in the photosynthetic system. Black and yellow boxes indicate up-regulated and down-regulated genes, respectively.
Figure 6. KEGG enrichment pathway of differential expression genes in the photosynthetic system. Black and yellow boxes indicate up-regulated and down-regulated genes, respectively.
Plants 12 01309 g006
Table
Table 1. MPH and HPH of panicle traits at the heading date.
Table 1. MPH and HPH of panicle traits at the heading date.
ZZBZ7-10ZZY10HPH (%)MPH (%)
PL (cm)26.95 ± 0.8229.83 ± 0.9431.45 ± 1.61−0.7710.79 **
1000-GW (g)1.54 ± 0.031.69 ± 0.051.94 ± 0.0415.03 **20.34 **
PBN14.89 ± 0.9313.87 ± 0.7814.36 ± 1.01−2.241.16
SBN53.22 ± 5.7866.56 ± 5.5559.22 ± 8.53−11.02 *−1.11
TN25 ± 3.1611 ± 2.8313 ± 1.79−47.20 **−26.67 **
PH (cm)123.36 ± 1.55123.02 ± 3.43146.07 ± 1.1516.81 **18.66 **
GL (cm)9.62 ± 0.279.24 ± 0.209.55 ± 0.09−0.800.80
GW (cm)2.32 ± 0.062.62 ± 0.092.42 ± 0.07−7.49 **−1.88
** Significant difference with p < 0.01; * Significant difference with p < 0.05. MPH and HPH stand for mid-parent heterosis and high-parent heterosis, respectively. PL, 1000-GW, PBN, SBN, TN, PH, GL, and GW refer to panicle length, 1000-grain weight, primary branch number, number of secondary branches, tiller number, plant height, grain length, and grain width, respectively. ± means standard deviation (SD).
Table
Table 2. MPH and HPH heterosis of panicle length during male gametogenesis.
Table 2. MPH and HPH heterosis of panicle length during male gametogenesis.
StageZZB (cm)Z7-10 (cm)ZZY10 (cm)HPH (%)MPH (%)
meiosis stage11.05 ± 0.4612.03 ± 0.4113.84 ± 0.5815.0719.98
uninucleate stage13.65 ± 0.6114.94 ± 0.3815.70 ± 0.485.119.85
bicellular stage20.29 ± 0.8621.08 ± 0.5124.03 ± 0.6914.0316.19
tricellular stage23.94 ± 0.6726.58 ± 0.5028.33 ± 0.726.5612.13
Heading stage26.95 ± 0.8229.83 ± 0.9431.45 ± 1.61−0.7710.79 **
Ripening stage29.84 ± 1.5029.76 ± 0.9730.23 ± 1.14−3.05 *1.44
** Significant difference with p < 0.01, * Significant difference with p < 0.05, MPH and HPH stand for mid-parent heterosis and high-parent heterosis, respectively. ± means standard deviation (SD).
Table
Table 3. DEGs between the hybrid ZZY10 and its parents.
Table 3. DEGs between the hybrid ZZY10 and its parents.
Genetic ClassSubclassExpression Patternp-Value (Total Number of DEGs)
0.05 (9000)
Number%
Up-regulatedOver-dominanceZZY10 > Z7 > ZZB126114.1
ZZY10 > ZZB > Z7132914.8
ZZY10 > ZZB = Z700
High-parent dominanceZZY10 = ZZB > Z73654.1
ZZY10 = Z7 > ZZB2362.6
AdditivityAdditivityZZB < ZZY10 < Z7258728.7
Z7 < ZZY10 < ZZB91010.1
Down-regulatedLow-parent dominance ZZY10 = ZZB < Z7580.6
ZZY10 = Z7 < ZZB1391.5
Under-dominanceZZY10 < ZZB < Z77318.1
ZZY10 < Z7 < ZZB138415.4
ZZY10 < ZZB = Z700
Table
Table 4. Top ten significant GO terms of DGHP.
Table 4. Top ten significant GO terms of DGHP.
CategoryGO_TermGO_Term_Annotationp Value
BPGO:0015979photosynthesis2.98 × 10−6
BPGO:0015074DNA integration5.93 × 10−5
BPGO:0042545cell wall modification0.000219
CCGO:0009579thylakoid6.08 × 10−6
CCGO:0044436thylakoid part6.08 × 10−6
CCGO:0009521photosystem1.31 × 10−5
CCGO:0034357photosynthetic membrane2.04 × 10−5
CCGO:0009522photosystem I0.000582
CCGO:0009523photosystem II0.000645
CCGO:0005618cell wall0.001305
BP and CC stand for biological process and cellular component, respectively.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhong, Z.; Wu, Y.; Zhang, P.; Hu, G.; Fu, D.; Yu, G.; Tong, H. Transcriptomic Analysis Reveals Panicle Heterosis in an Elite Hybrid Rice ZZY10 and Its Parental Lines. Plants 2023, 12, 1309. https://doi.org/10.3390/plants12061309

AMA Style

Zhong Z, Wu Y, Zhang P, Hu G, Fu D, Yu G, Tong H. Transcriptomic Analysis Reveals Panicle Heterosis in an Elite Hybrid Rice ZZY10 and Its Parental Lines. Plants. 2023; 12(6):1309. https://doi.org/10.3390/plants12061309

Chicago/Turabian Style

Zhong, Zhengzheng, Yawen Wu, Peng Zhang, Guocheng Hu, Dong Fu, Guoping Yu, and Hanhua Tong. 2023. "Transcriptomic Analysis Reveals Panicle Heterosis in an Elite Hybrid Rice ZZY10 and Its Parental Lines" Plants 12, no. 6: 1309. https://doi.org/10.3390/plants12061309

APA Style

Zhong, Z., Wu, Y., Zhang, P., Hu, G., Fu, D., Yu, G., & Tong, H. (2023). Transcriptomic Analysis Reveals Panicle Heterosis in an Elite Hybrid Rice ZZY10 and Its Parental Lines. Plants, 12(6), 1309. https://doi.org/10.3390/plants12061309

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

Article Metrics

Back to TopTop