Kai Pi
Ying Huang1,2
Renxiang Liu- 1College of Tobacco, Guizhou University, Guiyang, China
- 2Key Laboratory of Tobacco Quality in Guizhou Province, Guiyang, China
- 3College of Agriculture, Guizhou University, Guiyang, China
Heterosis has greatly improved the yield and quality of crops. However, previous studies often focused on improving the yield and quality of the shoot system, while research on the root system was neglected. We determined the root numbers of 12 F1 hybrids, all of which showed strong heterosis, indicating that tobacco F1 hybrids have general heterosis. To understand its molecular mechanism, we selected two hybrids with strong heterosis, GJ (G70 × Jiucaiping No.2) and KJ (K326 × Jiucaiping No.2), and their parents for transcriptome analysis. There were 84.22% and 90.25% of the differentially expressed genes were overdominantly expressed. The enrichment analysis of these overdominantly expressed genes showed that “Plant hormone signal transduction”, “Phenylpropanoid biosynthesis”, “MAPK signaling pathway - plant”, and “Starch and sucrose metabolism” pathways were associated with root development. We focused on the analysis of the biosynthetic pathways of auxin(AUX), cytokinins(CTK), abscisic acid(ABA), ethylene(ET), and salicylic acid(SA), suggesting that overdominant expression of these hormone signaling pathway genes may enhance root development in hybrids. In addition, Nitab4.5_0011528g0020、Nitab4.5_0003282g0020、Nitab4.5_0004384g0070 may be the genes involved in root growth. Genome-wide comparative transcriptome analysis enhanced our understanding of the regulatory network of tobacco root development and provided new ideas for studying the molecular mechanisms of tobacco root development.
Introduction
The growth period of tobacco is about 120 to 130 days. The plant is tall and needs a lot of water and nutrients for growth and development, and these substances are transported to the shoot via the root system (Zhaojun, 2021). The root system is an indispensable organ for tobacco growth and consists of three parts: main root, lateral root, and adventitious root. Many lateral roots can produce secondary lateral roots and tertiary lateral roots, which become the main part of the root system (Yang et al., 2002). Tobacco roots, like most plant roots, have the functions of absorbing and transporting soil moisture and inorganic salts, fixing and supporting plants, and synthesizing amino acids and proteins (Meister et al., 2014), which have a great impact on the growth and development of tobacco plants and the yield and quality of tobacco leaves. In addition, the iconic tobacco product, nicotine, is synthesized in the root system and then transported to the leaves (Solt, 1959) and is extremely important for the smoking quality of tobacco (Tong et al., 2020). Heterosis is a common natural phenomenon, which refers to when heterozygotes generated by hybridization are superior to the two parents in terms of growth potential, biomass, stress resistance, and adaptability (Shull, 1908). In most previous studies, heterosis is usually demonstrated in terms of biomass, yield, and stress resistance (Fujimoto et al., 2018; Li et al., 2021). In recent years, however, some studies have shown that the level of heterosis in root traits may be higher than that in the overground part (Ju et al., 2018; Dafna et al., 2021), and some efforts have been made to reveal the molecular mechanism of heterosis in root development (Paschold et al., 2010; Shalby et al., 2021), which means that the root system may be the perfect organ to study the genetic basis of tobacco heterosis. Although there are several genetic models used to describe heterosis, including dominant, overdominant, and epistasis, but its molecular basis remains poorly understood (Chen, 2013). To better understand heterosis, Thiemann et al. (2014) proposed two models of gene expression correlation: additive and non-additive expression patterns. In different crops, gene expression shows different additive and non-additive expression patterns (Tian et al., 2018; Shahzad et al., 2020; Shalby et al., 2021). In recent years, a large number of transcriptome research advances have provided new insights into the molecular basis of heterosis in species such as Arabidopsis (Fujimoto et al., 2012; Wang et al., 2019; Liu et al., 2020), rice (Shao et al., 2019; Ren et al., 2020), maize (Ko et al., 2016), and rape (Xiong et al., 2022).
Root development is the result of the synergistic action of various plant hormones. Hormones play an important role in root development through regulation of gene expression (Jia et al., 2022). Auxin (AUX), one of the most important plant hormones, regulates the growth of primary and lateral roots by promoting cell division and elongation (Rivas et al., 2022). Furthermore, cytokinins (CTKs) have been reported to have important control roles in the regulation of root structure and nutrients (Muraro et al., 2011). Abscisic acid (ABA) has an important control role in regulating root adaptive responses (Belda-Palazon et al., 2018; Ma et al., 2018). Ethylene (ETH) plays an important control role in root cell elongation (Song et al., 2019). Salicylic acid (SA) controls root growth by regulating the activity of root apical meristem (Bagautdinova et al., 2022), and studies have shown that Arabidopsis hybrids with superior SA content exhibit root growth heterosis (Zhang et al., 2016). The research shows that the number of root cells, cell differentiation and other processes will affect the root development (Cuadrado et al., 1987). Hormone content can also affect root development by affecting the number of root cells, cell differentiation and other processes (McCarthy-Suárez, 2021).
In tobacco, heterosis can be used to obtain tobacco hybrids with higher nicotine and potassium contents than the parents, and our previous studies have shown that nicotine heterosis is due to the higher efficiency of nicotine synthesis in hybrids (Tian et al., 2018; Mo et al., 2021). Since nicotine is synthesized in the root system (Solt, 1959), tobacco roots may be the basis for the heterosis of nicotine and potassium contents, which directly affects the development and quality of tobacco leaves. Therefore, this study preliminarily revealed the expression of early root heterosis in tobacco and performed transcriptome sequencing of two hybrids with strong heterosis to analyze the differentially expressed genes (DEGs), gene expression patterns, and biological processes that mediate heterosis in root growth. These new findings may help to reveal the biological mechanism of tobacco root heterosis and discover related candidate genes. These findings may also be significant in guiding the selection of new tobacco germplasm with high yield and high quality.
Materials and methods
Plant materials, growth conditions, and sample preparation
Based on our previous research (Tian et al., 2018; Pi et al., 2022), we selected the hybrids and their parents with different nicotine and potassium contents, and 12 hybrids were made from 7 parents according to the incomplete diallel crossing (NCII) method(Supplementary File 1, Table 1). All materials were provided by Guizhou Key Laboratory of Tobacco Quality Research. From 2021-2022, in the tobacco research base of Guizhou University, Yangwu Township, Anshun City, Guizhou Province, China, Seeds were sown in greenhouses using the floating seedling method, and the seedlings with six real leaves were sampled. The methods were as follows: Firstly, 3 seedlings with similar growth vigor were selected for root phenotype determination. After that, fresh biological samples were frozen with liquid nitrogen and stored in an ultra-low temperature refrigerator at - 80 ℃. These fresh samples were used for transcriptome sequencing analysis and RT-qPCR.
Determination of root phenotype and calculation of heterosis
The cleaned complete root system was scanned with the root scanner LA2400 to obtain a clear image of the root structure. Use WinRhizo Pro 2021a software (Regent, Canada) to count tobacco root tips.
RNA isolation and sequencing
Total RNA was extracted from the root tissue using TRIzol® Reagent (Plant RNA Purification Reagent for plant tissue) according to the manufacturer’s instructions (Invitrogen), and genomic DNA was removed using DNase I (TaKara). Then RNA quality was determined by 2100 Bioanalyser (Agilent) and quantified using the ND-2000 (NanoDrop Technologies). Only high-quality RNA sample (OD260/280 = 1.8~2.2, OD260/230≥2.0, RIN≥6.5, 28S:18S≥1.0, >1μg) was used to construct sequencing library. Paired-end RNA-seq sequencing library was sequenced with the NovaSeq 6000 sequencer (2 × 150bp read length).
Transcriptomics data processing and analysis
By preprocessing the raw reads, short sequences with a length <25 nt and low quality sequences were removed. After preprocessing, the obtained reads were mapped to the Nicotiana tabacum sequenced cultivar K326 genome (Edwards et al., 2017) using the splice-aware mapping tool, Tophat2 (Kim et al., 2019). RSEM was used to quantify gene abundances (Li and Dewey, 2011). Essentially, differential expression analysis was performed using the DESeq2 (Love et al., 2014), DEGs with |log2FC| ≥ 1 and p ≤ 0.05 were considered to be significantly different expressed genes. Gene ontology (GO) functional-enrichment analysis were carried out by Goatools (Klopfenstein et al., 2018). In addition, we use PlantTFDB (http://planttfdb.gao-lab.org/) Transcription factor analysis of genes.
Classification of gene expression patterns
In order to divide the differentially expressed genes into 12 expression modes, according to the definition of Rapp et al. (2009). we used STEM software (Ernst and Bar-Joseph, 2006) to perform additional analysis on the differentially expressed genes. To allow clustering on a reasonable number of possible model profiles, the parameter for “STEM clustering method,” model profiles “was set to 50 and 2 was selected as the” maximum unit change between time points “, and” Minimum Absolute Expression Change “was set to 0 to classify the expression patterns of all differentially expressed genes.
Identification of core genes in tobacco root development
The typical genes of Arabidopsis root development were used to identify core genes in tobacco. BLASTP identified 6 core genes with E-values truncated at 1E-10. To further classify root development genes in tobacco, phylogenetic analyses were performed using homologous protein sequences from Arabidopsis thaliana and Triticum aestivum L. The phylogenetic tree was constructed by MEGA 11.0(https://www.megasoftware.net/), The ultrafast bootstrap with 1000 replicates was conducted to obtain the supporting values for each bootstrap of the tree. Finally, gene trees were visualized and colored using iTOL(https://itol.embl.de/).
RT-qPCR validation
For the validation of the transcriptome data, we randomly selected 6 DEGs for RT-qPCR (PCR quantitative real-time) analysis. Supplementary File 1: Table 16 lists the genes and corresponding primers used in RT-qPCR. The total RNA used in the sequencing was reversely transcribed to obtain cDNA, which was used as a template to amplify the target genes, and the RT-qPCR experiment was conducted. The RNA was reversely transcribed into cDNA using the FastKing reverse transcription kit (Tiangen, China) as per the manual’s instructions. qPCR was performed using the BIO-RAD CFX96 Real-Time qPCR system. The relative expression of each gene was calculated using 2−△△Ct (Livak and Schmittgen, 2001).
Statistical analysis
SPSS 25.0 software was used for statistical analyses. The variance analysis of root phenotype was carried out using the Duncan’s new multiple range method (P < 0.05). Based on the root phenotype, the values of mid-parent heterosis (MPH) were calculated according to the following method, , where F1represents the value of the first generation of hybrid and MP represents the average value of parents.
Results
Heterosis in the root number of F1 hybrids
Based on the screening results of tobacco root number heterosis (Supplementary File 1, Table 1), the root numbers of the 12 F1 hybrids all showed strong heterosis values, indicating that tobacco F1 hybrids have general heterosis. Two highly dominant F1 hybrids, GJ (G70 × Jiucaiping No.2) and KJ (K326 × Jiucaiping No.2), were selected from these F1 hybrids for further analysis (their median parent values were not significantly different). At the tobacco seedling stage, we observed that the root development of the two hybrids were significantly better than that in their parents (Figure 1A). The number of root systems in GJ was 623, whereas G and J had 515 and 343, respectively; the number of roots in KJ was 878, whereas K and J had 711 and 343, respectively. The results showed that the average number of roots of both hybrids were significantly higher than that of their parents (Figure 1B), and the heterosis was 45.66% (GJ) and 67.12% (KJ) (Figure 1C).
Figure 1 Phenotypic parameters and Heterosis of two tobacco F1 hybrids and their parents. (A) Photographs of the two hybrids and their parents used in this study. (B) Root number of F1 hybrid and its parents(G = G70, J = Jiucaiping No.2, GJ = G70 × Jiucaiping No.2, K = K326, KJ =K326 × Jiucaiping No.2). (C) Heterosis of two tobacco F1 hybrids (GJ and KJ). Error bars indicate standard error. Different small letters showed significant differences (P<0.05).
Global transcriptome analysis of heterosis in root number
Using two hybrids, GJ and KJ, with strong heterosis and their three parents as sequencing materials, 15 cDNA libraries were constructed and sequenced on the Illumina Novaseq 6000 platform. Each sample generated 76-98 million clean reads in RNA-Seq, which were then used for further analysis. Clean reads were mapped to the K326 genome, and about 95% of clean reads could be mapped to the genome (Supplementary File 1, Table 2). A total of 39,345 genes were expressed in 15 analyzed samples (Supplementary File 1, Table 3). Among them, 8.30% of the genes were very highly expressed (FPKM ≥ 50), 64.19% and 27.51% of the genes were moderately (1 ≤ FPKM ≤ 10) and highly expressed (10 ≤ FPKM ≤ 50), respectively. These results reflect that there is no significant difference in expression coverage between hybrids and parents during development. However, according to hierarchical clustering analysis(Figure 2A) and principal component analysis of parents and hybrids (Figure 2B), both parents are closely clustered, and hybrids can be clearly divided into two different expression situations. These results show that there are great differences in gene expression between parents and hybrids during root development.
Figure 2 Global gene expression profiling during root development. (A) Pearson correlation analysis of the two hybrids and their parents, and the color code on the right represents Pearson correlation coefficient. (B) Principal component analysis (PCA) of RNA-Seq data. Each point shows the average of three repeats.
Transcriptome differences between roots of hybrids and their parents
In the case of significance level P ≤ 0.05 and | log2 fold change | ≥ 1, we identified 1430 genes up-regulated and 1697 genes down-regulated between G (female parent) and GJ (hybrid), 2015 genes up-regulated and 1809 genes down-regulated between J (male parent) and GJ (hybrid) (Figures 3A, B; Supplementary File 1, Table 4), 1613 genes up-regulated and 1922 genes down-regulated between the middle parent and GJ hybrid, and 302 genes were up-regulated and 291 genes were down-regulated between G and J (Supplementary File 1, Table 5). We identified 3082 genes up-regulated and 3007 genes down-regulated between K (female parent) and KJ (hybrid), 4118 genes up-regulated and 3810 genes down-regulated between J (male parent) and KJ (hybrid) (Figures 3C, D), 3683 up-regulated genes and 3874 down-regulated genes between middle parent and KJ hybrid, and 283 genes up-regulated and 342 genes down-regulated between K and J (Supplementary File 1, Table 6). This demonstrates that there are relatively few differential genes between the parents of the two strong heterosis hybrids, but there are a large number of differential genes between the parents and the hybrids, which may be the reason for root number heterosis.
Figure 3 Distribution of differentially expressed genes (DEGs) between hybrids and parents. MP represents the average value of parents. (A) The number of different expressed genes between GJ and its parents during root development. (B) Venn shows the distribution of GJ and parental DEGs. (C) The number of different expressed genes between KJ and its parents during root development. (D) Venn shows the distribution of KJ and parental DEGs.
F1 hybrid showed a overdominant gene expression pattern
To further analyze the DEGs of hybrids and parents, the genes were divided into 12 expression patterns (P1-P12, Figure 4A) according to the definition of Rapp et al. (2009). The genes in the P1 and P2 modes are additive, P3-P6 modes are dominant, and P7-P12 modes show overdominant expression, in which the genes in P7-P9 are down-regulated overdominant and P10-P12 are up-regulated overdominant. Among the overdominant genes, 238, 376, and 1564 genes of GJ hybrid showed up-regulated overdominant expression patterns(Supplementary File 1, Table 7), and 516, 857, and 733 genes showed down-regulated overdominant expression patterns; there were 373, 1160, and 2750 genes in KJ hybrid that showed up-regulated overdominant expression patterns(Supplementary File 1, Table 8), and 1375, 1995, and 772 genes that showed down-regulated overdominant expression patterns (Figure 4B). Among these non-additively expressed genes (P3-P12), GJ hybrid had the highest proportion (84.22%) in the expression pattern of overdominant (P7-P12), and KJ hybrid had the highest proportion (90.25%) in the expression pattern of overdominant (P7-P12) (Figure 4C). Therefore, our results show that the overdominant expression advantage is the main reason for the formation of heterosis of tobacco roots.
Figure 4 12 gene expression patterns of DEGs between hybrids and parents. (A) Classification of 12 gene expression patterns (♂: paternal; H: hybrid; ♀: maternal). (B) The number of genes of GJ and KJ hybrids in each model. (C) The number and proportion of the five total expression patterns after classification in DEGs.
Enrichment analysis of overdominant genes
Gene Ontology (GO) enrichment analysis was performed on the overdominant expression gene sets in GJ and KJ hybrids. We found that most overdominant genes were involved in “biological process” (55.3 – 56.1%). In GJ hybrids, most of them were significantly enriched in the jasmonic acid biosynthetic process, riboflavin metabolic process, riboflavin biosynthetic process, aminoglycan metabolic process, and jasmonic acid metabolic process (Supplementary File 1, Table 9; P < 0.05). In KJ hybrids, the overdominant genes were significantly enriched in the glucosamine containing compound catabolic process, chitin metallic process, amino sugar catabolic process, chitin catabolic process, and aminoglycan catabolic process (Supplementary File 1, Table 10; P < 0.05). According to the GO enrichment analysis of GJ and KJ co overdominant expression genes, most of the genes involved in cell differentiation, proliferation, cell growth and other processes were significantly enriched. It indicated that the overdominant expression of genes in the process of cell differentiation and proliferation was related to the formation of root heterosis.
According to the analysis of Kyoto Encyclopedia of Genes and Genomes (KEGG) metabolic pathways, 4284 overdominant genes in GJ were annotated to 117 pathways (Supplementary File 1, Table 11; P < 0.05), whereas 8425 dominant genes in KJ were annotated to 123 pathways (Supplementary File 12, Table 12; P < 0.05). As shown in Figures 5A, B, the overdominant gene enrichment pathways of the two hybrids are mainly “phylpropanoid biosynthesis” (9.72 and 5.27%, respectively), “plant hormone signal transduction” (4.97 and 4.93%, respectively), “MAPK signaling pathway - plant” (4.91 and 4.06%, respectively), and “start and cross metadata” (2.48 and 2.56%, respectively). In addition, 84 and 158 overdominant genes of the two hybrids were related to “plant hormone signal transduction”, and 42 and 82 genes were related to “starch and sucrose metabolism”. Based on KEGG analysis, we further studied specific functional genes related to root development.
Figure 5 Enrichment analysis of KEGG pathway of GJ and KJ overdominant expression genes. (A) KEGG analysis of overdominant up-regulated and down-regulated genes in GJ hybrids. (B) KEGG analysis of overdominant up-regulated and down-regulated genes in KJ hybrids.
Overdominant expression of genes related to phytohormone signaling pathway during root development
Combined with functional analysis, in order to strengthen our understanding of hormone signaling during heterosis development in tobacco root growth, we focused on the overdominant genes of hormone signal transduction related to the biosynthesis of AUX, CTK, ABA, ET, and SA (Figure 6; Supplementary File 1, Table 13). In the AUX signal transduction pathway, a total of 11 genes were overexpressed in F1 hybrids (GJ and KJ), including 3 that were annotated to the AUX1 gene, 5 that were annotated to the GH3 gene, and 3 that were annotated to the AUX response gene SAUR. In the CTK signal transduction pathway, one receptor histidine kinase (CRE1) gene was negatively overexpressed in F1 hybrids (GJ and KJ). One gene was annotated to A-type response regulator (A-ARR), which showed up-regulated overdominant expression in hybrids. In the ABA signal transduction pathway, 9 genes were found to be overexpressed in F1 hybrids, including 5 that were annotated to protein phosphatase (PP2C) and 4 encoding Subcrosse non-fermenting 1-related protein kinase 2 (SnRK2). Interestingly, the genes related to PP2C showed up-regulated overdominant expression in both hybrids, of which 2 genes were highly expressed (10 ≤ FPKM ≤ 50). In the ET signal transduction pathway, 4 genes were annotated to ET receptor sensor (ETR2), one was annotated to EIN3, three were annotated to bind F-box protein 1/2 (EBF1/2), and four were annotated to ET responsive transcription factor 1/2 (ERF1/2). These ET signal transduction pathway results show that the genes related to ETR2, EIN3 and EBF1/2 are up-regulated and overdominantly expressed in both hybrids. For the SA signal transduction pathway, a total of 11 genes were overdominantly expressed in F1 hybrids (GJ and KJ). Among them, 1 gene was annotated to TGA transcriptional regulator, which showed down-regulated overdominant expression in the roots of hybrids, and 7 genes were annotated to protein PR-1. It is worth noting that these 7 PR-1 related genes not only showed an overdominant expression pattern in the hybrid, but also had high differential expression multiples and highly expressed (FPKM ≥ 50).
Figure 6 Overdominant expression genes of hormone signal transduction related to AUX, CTK, ABA, ET and SA biosynthesis. These genes showed overdominant expression in both hybrids.
Overdominant expression of genes related to cell development during root development
Among the common overdominant expression genes of the two hybrids, various GO categories such as cell population promotion, cell differentiation, multiple cell growth and primary meristem tissue development are significantly enriched. In the process of cell differentiation, most genes show up-regulated overdominant expression (Figure 7), in which Nitab4.5_ 0002714g0010、Nitab4.5_ 0008353g0010、Nitab4.5_ 0002918g0030 showed high expression in hybrids. In the process of cell population promotion and cell population promotion, 3 and 5 genes were up-regulated overdominant expression in the hybrid, indicating that cell population promotion and cell population promotion were important processes of root heterosis. In the process of primary meristem issue development, the two genes are both expressed down-regulated in the hybrid. In general, a lot of changes have taken place in cell population promotion, cell differentiation and multiple cell growth of hybrids, which may be the reason for the formation of the heterosis of root.
Figure 7 Overdominant expression genes related to cell population promotion, cell differentiation, multiple cell growth and primary commodity issue development. These genes were overdominant expression in both hybrids.
Identification of gene co-expression network modules
A weighted gene co-expression network analysis (WGCNA) was performed on all genes to thoroughly characterize the expression of genes related to root number. WGCNA identified 10 vivid modules (including 19338 genes) (Figure 8A; Supplementary File 1, Table 14). The module analysis revealed that there was a significant correlation between the brown module and the heterosis of the number of hybrid roots (Figure 8B). The genes in the “brown” module are mainly related to cysteine biological process, protein serine/threonine kinase activity, ethylene binding, etc. Further analysis of 6 genes with degree>700 in the brown module, Nitab4.5_ 0000119g0080 gene is Lysine tRNA ligase, Nitab4.5_ 0004384g0070、Nitab4.5_ 0011528g0020、Nitab4.5_ 0003282g0020 gene is Serine/threonine -/dual specificity protein kinase, Nitab4.5_ 0002818g0040 gene is NAD dependent emergency/dehydratase, Nitab4.5_ 0009232g0030 gene is Mitochondrial carrier protein (Figure 8C). Interestingly, these genes are up-regulated overdominant expression patterns. We believe that these genes may be related to the formation of root heterosis.
Figure 8 Weighted gene co-expression network analysis (WGCNA). (A) Hierarchical cluster tree shows 10 modules of co-expressed genes. The lower heatmap shows modules in designated colors. The module grey is for unassigned genes. (B) Heat maps showing the correlation of module-trait, the Pearson correlation coefficient and P values of significant modules are given. (C) Connectivity between 6 genes.
RT-qPCR results
RT-qPCR was used to analyze 6 overdominant genes that have been shown to be involved in root development in other crops. In brief, protein sequences of well-characterized genes of each family from Arabidopsis thaliana were used as BLAST query sequences (Supplementary File 1, Table 15), and the following sequences in tobacco were filtered by stringent threshold (E-value <1E-25, % identity>30 and % query coverage>30). In tobacco, protein sequences of root development exhibited a clear topological relationship, with high bootstrap values (Arabidopsis thaliana, Triticum aestivum L.)(Figure 9A). Actin gene was used as a reference gene to standardize the expression levels of each gene. Supplementary File 1: Table 16 lists the genes and corresponding primers used in the RT-qPCR. The results showed that (Figure 9B), the genes of root development also showed a trend of overdominant expression. Although the expression multiples of RT-qPCR data and RNA-seq sequencing results are slightly different, they basically show a consistent expression trend, so RNA-seq data has certain reliability.
Figure 9 (A) Phylogenetic trees of root development genes in tobacco. (B) Comparison of RNA-seq and RT-qPCR gene expression levels of root related genes in two hybrids and their parents.
Discussion
Heterosis refers to the superiority of hybrid F1 to either parent for a specific trait. The utilization of heterosis has improved crop yield and quality worldwide (Meyer et al., 2004). Many studies have shown that heterosis can be observed in the seedling stage of plants, such as Arabidopsis (Zhu et al., 2016; Liu et al., 2020), rice (Ma et al., 2011; Zhu et al., 2016; Liu et al., 2018), corn (Paschold et al., 2012), and wheat (Liu et al., 2018). The application of transcriptome analysis in heterosis research has effectively revealed the molecular basis of heterosis in Arabidopsis (Groszmann et al., 2015; Liu et al., 2021), rice (Katara et al., 2020), corn (Song et al., 2013; Guttikonda et al., 2020), cotton (Shahzad et al., 2020), and soybean (Zhang et al., 2017), but the molecular basis of heterosis in tobacco root growth has not been reported in the literature. In the early stage of tobacco growth, improved root development plays an important role in enhancing the absorption of nutrients and water by plants, promoting their growth and increasing their yield (Zhang et al., 2020; Shalby et al., 2021). From this point of view, through comparative transcriptome analysis of roots, the genetic basis of root development heterosis of two strong heterosis hybrids and their parents was studied.
Plant hormones are signal molecules produced by plants, which are used to regulate various aspects of plant growth and development (Santner et al., 2009). Plant hormone signals are sensed and transmitted to the nucleus through a series of signal transduction elements which regulate gene expression and trigger a series of physiological processes. Through the KEGG pathway enrichment analysis of the overdominantly expressed genes, we found that the overdominant genes encoding hormone signal transduction between F1 hybrids and parents were significantly enriched. In order to reveal that hormone signal transduction is involved in the growth and development of tobacco roots, we focused on AUX, CTK, ABA, ET, and SA.
Auxin, including biosynthesis and transport of AUX and its signaling, plays a crucial role in regulating root growth and development (Saini et al., 2013). A total of 11 overdominant expression genes (ODEs) encoding AUX signal components in the two hybrids were identified in this study, including 3 AUX1 genes, 5 GH3 genes, and 3 SAUR genes. GH3s and SAURs are AUX early response genes (Yang et al., 2015; Stortenbeker and Bemer, 2019). With 8 of the 11 ODEs being AUX early response genes, this shows that AUX plays a core role in regulating the growth and development of tobacco roots.
CTK regulates many developmental processes of plants, including chloroplast formation, root growth, and nutrient absorption (Zubo and Schaller, 2020). A-ARRs have been reported as negative regulators of the CTK signaling pathway (To et al., 2004; Kieber and Schaller, 2014), inhibiting CTK signaling, whereas CTK receptor CRE1 plays an important role in regulating lateral root development (Gonzalez-Rizzo et al., 2006). In our study, CRE1 was shown to be down-regulated overdominant and type A-ARRs were up-regulated overdominant, which indicated that CTKs accumulated less in the roots of the hybrid F1 generation.
Abscisic acid signaling is also important during root development (Brookbank et al., 2021). Some studies have shown that ABA is a negative regulator of lateral root appearance, but genetic evidence also suggests that ABA and AUX have a regulatory interaction in lateral root formation (Wang et al., 2019). SnRK2 has been reported to organize root development under non-stressful conditions and is necessary in all root tissues (Kawa et al., 2020). Overexpression of NtSnrk2.2 increased soluble sugar accumulation, increased lateral roots, and improved root development in tobacco plants (Liu et al., 2020). Group A protein phosphorase 2Cs (PP2Cs) are ABA co-receptors that negatively regulate the ABA signaling pathway by inhibiting the downstream SnRK2 protein kinase. 9 ODEs encoding for ABA signaling were found in the two hybrids.
Ethylene is another important plant hormone. It has been reported that ET stimulates the root growth of many plant species (Qin et al., 2019). Qin and Huang, (2018) reported that ET regulates the root development by changing the synthesis of AUX in Arabidopsis thaliana; jasmonic acid and ET jointly regulate AUX signaling and promote root development (Xu et al., 2020). The current results show that 12 ODEs in the two hybrids jointly encode AUX signal transduction, including 4 ETR genes (3 were AUX1 genes), 3 EIN3 genes, 3 EBF1/2 genes, and 4 ERF1/2 genes. Among the 12 ODEs, 4 ODEs encoding ERF1/2 were found to be highly expressed, indicating that ERF1/2 may play a role in regulating root growth and development.
Salicylic acid was identified as the sixth plant hormone in 1992 (Raskin, 1992). In recent years, a large number of studies have confirmed the important role of SA in plant root genesis (Bagautdinova et al., 2022). Arabidopsis hybrid varieties with sub optimal and super optimal SA content show root growth heterosis (Zhang et al., 2016), and chromatin remodeling reduces DNA methylation 1 (DDM1), linking heterosis with endogenous SA levels. In our study, the genes for salicylic acid synthesis showed overdominant expression. In the future, we can further study the effect of salicylic acid on the formation of root heterosis.
A weighted gene co-expression network analysis (WGCNA), including three genes Nitab4.5_ 0011528g0020、Nitab4.5_ 0003282g0020、Nitab4.5_ 0004384g0070 has a serial threonine/tyrosine protein kinase catalytic domain. The results showed that this domain affected root development of Arabidopsis thaliana through microtubules in primary root zone cells (Sheremet et al., 2010). OsESG1 gene in rice regulates the initiation and development of crown and root by controlling auxin response and distribution. The gene also has a Serine threonine/tyrosine protein kinase catalytic domain (Pan et al., 2020). It is worth noting that Nitab4.5_ 0003282g0020 gene has a cysteine rich plant receptor like kinase domain. Research shows that CRK28 in Arabidopsis plays an important role in root growth and epidermal cell differentiation (Pelagio-Flores et al., 2019). Therefore, we believe Nitab4.5_ 0011528g0020, Nitab4.5_ 0003282g0020 and Nitab4.5_ 0004384g0070 may be the key gene regulating the root growth of hybrid.
Conclusions
In this study, universal heterosis in the root system of hybrids was proved by measuring the number of roots. Comparative transcriptome analysis of hybrids and their parents showed that heterosis was related to the overall gene expression pattern. GJ and KJ hybrids had 5.08% and 10.87% DEGs compared with their parents, respectively; 84.22% and 90.25% of these differential genes showed an overdominant expression pattern, respectively. Through the GO function enrichment analysis of overdominantly expressed genes, it was found that the “plant hormone signal transduction”, “phylpropanoid biosynthesis”, “MAPK signaling pathway plant”, and “starch and supra metadata” pathways were related to root development. We focused on analyzing the biosynthetic pathways of AUX, CTK, ABA, ET, and SA and showed that the overdominantly expressed genes of these hormone signal transduction pathways may enhance root development in hybrids. Cell population promotion, cell differentiation and multiple cell growth in hybrids may be the reason of root heterosis. In addition, Nitab4.5_0011528g0020、Nitab4.5_0003282g0020、Nitab4.5_0004384g0070 may be genes involved in root growth. To sum up, this study provides a new understanding of the relevant mechanisms of tobacco root heterosis formation. However, further gene function studies are needed to elucidate the development of tobacco root heterosis.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE218712.
Author contributions
KP planned and designed the research, analyzed the data, and wrote the manuscript. WL and SZ performed most of the field work; ZM and LD carried out the molecular biology studies. RL and YH conceived the study and participated in the design and coordination. All authors reviewed the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This study was funded by the National Science Foundation of China (32060510), the Department of Science and Technology of Guizhou Province ([2016]5663 and [2019]1405), and the Guizhou Tobacco Company (2022XM02). The funding institution Guizhou tobacco company was not involved in the study design, collection, analysis, interpretation of data, the writing of this article, or the decision to submit it for publication.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2023.1107550/full#supplementary-material
References
Bagautdinova, Z. Z., Omelyanchuk, N., Tyapkin, A. V., Kovrizhnykh, V. V., Lavrekha, V. V., Zemlyanskaya, E.V. (2022). Salicylic acid in root growth and development. Int. J. Mol. Sci. 23, 2228. doi: 10.3390/ijms23042228
Belda-Palazon, B., Gonzalez-Garcia, M. P., Lozano-Juste, J., Coego, A., Antoni, R., Julian, J., et al. (2018). Pyl8 mediates aba perception in the root through non-cell-autonomous and ligand-stabilization-based mechanisms. Proc. Natl. Acad. Sci. U S 115, E11857–E11863. doi: 10.1073/pnas.1815410115
Brookbank, B. P., Patel, J., Gazzarrini, S., Nambara, E. (2021). Role of basal aba in plant growth and development. Genes (Basel) 12, 1936. doi: 10.3390/genes12121936
Chen, Z. J. (2013). Genomic and epigenetic insights into the molecular bases of heterosis. Nat. Rev. Genet. 14, 471–482. doi: 10.1038/nrg3503
Cuadrado, A., Cánovas, J., Navarrete, M. H. (1987). Influence of cell size on differentiation of root meristem cells. Environ. Exp. Bot. 27, 273–277. doi: 10.1016/0098-8472(87)90036-0
Dafna, A., Halperin, I., Oren, E., Isaacson, T., Tzuri, G., Meir, A., et al. (2021). Underground heterosis for yield improvement in melon. J. Exp. Bot. 72, 6205–6218. doi: 10.1093/jxb/erab219
Edwards, K. D., Fernandez-Pozo, N., Drake-Stowe, K., Humphry, M., Evans, A. D., Bombarely, A., et al. (2017). A reference genome for nicotiana tabacum enables map-based cloning of homeologous loci implicated in nitrogen utilization efficiency. BMC Genomics 18, 448. doi: 10.1186/s12864-017-3791-6
Ernst, J., Bar-Joseph, Z. (2006). Stem: a tool for the analysis of short time series gene expression data. BMC Bioinf. 7, 191. doi: 10.1186/1471-2105-7-191
Fujimoto, R., Taylor, J. M., Shirasawa, S., Peacock, W. J., Dennis, E. S. (2012). Heterosis of arabidopsis hybrids between c24 and col is associated with increased photosynthesis capacity. Proc. Natl. Acad. Sci. U S 109, 7109–7114. doi: 10.1073/pnas.1204464109
Fujimoto, R., Uezono, K., Ishikura, S., Osabe, K., Peacock, W. J., Dennis, E. S. (2018). Recent research on the mechanism of heterosis is important for crop and vegetable breeding systems. Breed Sci. 68, 145–158. doi: 10.1270/jsbbs.17155
Gonzalez-Rizzo, S., Crespi, M., Frugier, F. (2006). The medicago truncatula cre1 cytokinin receptor regulates lateral root development and early symbiotic interaction with sinorhizobium meliloti. Plant Cell 18, 2680–2693. doi: 10.1105/tpc.106.043778
Groszmann, M., Gonzalez-Bayon, R., Lyons, R. L., Greaves, I. K., Kazan, K., Peacock, W. J., et al. (2015). Hormone-regulated defense and stress response networks contribute to heterosis in arabidopsis f1 hybrids. Proc. Natl. Acad. Sci. U S 112, E6397–E6406. doi: 10.1073/pnas.1519926112
Guttikonda, H., Thummala, S. R., Agarwal, S., Mangrauthia, S. K., Ramanan, R., Neelamraju, S., et al. (2020). Genome-wide transcriptome profile of rice hybrids with and without oryza rufipogon introgression reveals candidate genes for yield. Sci. Rep. 10, 4873. doi: 10.1038/s41598-020-60922-6
Jia, Z., Giehl, R., von Wiren, N. (2022). Nutrient-hormone relations: driving root plasticity in plants. Mol. Plant 15, 86–103. doi: 10.1016/j.molp.2021.12.004
Ju, C., Zhang, W., Liu, Y., Gao, Y., Wang, X., Yan, J., et al. (2018). Genetic analysis of seedling root traits reveals the association of root trait with other agronomic traits in maize. BMC Plant Biol. 18, 171. doi: 10.1186/s12870-018-1383-5
Katara, J. L., Verma, R. L., Parida, M., Ngangkham, U., Molla, K. A., Barbadikar, K. M., et al. (2020). Differential expression of genes at panicle initiation and grain filling stages implied in heterosis of rice hybrids. Int. J. Mol. Sci. 21, 1080. doi: 10.3390/ijms21031080
Kawa, D., Meyer, A. J., Dekker, H. L., Abd-El-Haliem, A. M., Gevaert, K., Van De Slijke, E., et al. (2020). Snrk2 protein kinases and mrna decapping machinery control root development and response to salt. Plant Physiol. 182, 361–377. doi: 10.1104/pp.19.00818
Kim, D., Paggi, J. M., Park, C., Bennett, C., Salzberg, S. L. (2019). Graph-based genome alignment and genotyping with hisat2 and hisat-genotype. Nat. Biotechnol. 37, 907–915. doi: 10.1038/s41587-019-0201-4
Klopfenstein, D. V., Zhang, L., Pedersen, B. S., Ramírez, F., Warwick, V. A., Naldi, A., et al. (2018). Goatools: a python library for gene ontology analyses. Sci. Rep. 8, 10872. doi: 10.1038/s41598-018-28948-z
Ko, D. K., Rohozinski, D., Song, Q., Taylor, S. H., Juenger, T. E., Harmon, F. G., et al. (2016). Temporal shift of circadian-mediated gene expression and carbon fixation contributes to biomass heterosis in maize hybrids. PloS Genet. 12, e1006197. doi: 10.1371/journal.pgen.1006197
Li, B., Dewey, C. N. (2011). Rsem: accurate transcript quantification from rna-seq data with or without a reference genome. BMC Bioinf. 12, 323. doi: 10.1186/1471-2105-12-323
Li, P., Su, T., Zhang, D., Wang, W., Xin, X., Yu, Y., et al. (2021). Genome-wide analysis of changes in mirna and target gene expression reveals key roles in heterosis for chinese cabbage biomass. Hortic. Res. 8, 39. doi: 10.1038/s41438-021-00474-6
Liu, Y. J., Gao, S. Q., Tang, Y. M., Gong, J., Zhang, X., Wang, Y. B., et al. (2018). Transcriptome analysis of wheat seedling and spike tissues in the hybrid jingmai 8 uncovered genes involved in heterosis. Planta 247, 1307–1321. doi: 10.1007/s00425-018-2848-3
Liu, W., He, G., Deng, X. W. (2021). Biological pathway expression complementation contributes to biomass heterosis in arabidopsis. Proc. Natl. Acad. Sci. U S 118, 16. doi: 10.1073/pnas.2023278118
Liu, P. C., Peacock, W. J., Wang, L., Furbank, R., Larkum, A., Dennis, E. S.. (2020). Leaf growth in early development is key to biomass heterosis in arabidopsis. J. Exp. Bot. 71, 2439–2450. doi: 10.1093/jxb/eraa006
Liu, M., Wang, J., Gou, J., Wang, X., Li, Z., Yang, X., et al. (2020). Overexpression of ntsnrk2.2 enhances salt tolerance in nicotiana tabacum by regulating carbohydrate metabolism and lateral root development. Funct. Plant Biol. 47, 537–543. doi: 10.1071/FP19299
Livak, K. J., Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative pcr and the 2–δδct method. Methods 25, 402–408. doi: 10.1006/meth.2001.1262
Love, M. I., Huber, W., Anders, S. (2014). Moderated estimation of fold change and dispersion for rna-seq data with deseq2. Genome Biol. 15, 550. doi: 10.1186/s13059-014-0550-8
Ma, Q., Hedden, P., Zhang, Q. (2011). Heterosis in rice seedlings: Its relationship to gibberellin content and expression of gibberellin metabolism and signaling genes. Plant Physiol. 156, 1905–1920. doi: 10.1104/pp.111.178046
Ma, H., Liu, C., Li, Z., Ran, Q., Xie, G., Wang, B., et al. (2018). Zmbzip4 contributes to stress resistance in maize by regulating aba synthesis and root development. Plant Physiol. 178, 753–770. doi: 10.1104/pp.18.00436
McCarthy-Suárez, I. (2021). Supra-physiological levels of gibberellins/dellas modify the root cell size/number and the root architecture in root tips of a. thaliana seedlings. connections to the root hair patterning and abundance. Biorxiv, 2021–2027. doi: 10.1101/2021.07.25.453699
Meister, R., Rajani, M. S., Ruzicka, D., Schachtman, D. P. (2014). Challenges of modifying root traits in crops for agriculture. Trends Plant Sci. 19, 779–788. doi: 10.1016/j.tplants.2014.08.005
Meyer, R. C., Torjek, O., Becher, M., Altmann, T. (2004). Heterosis of biomass production in arabidopsis. establishment during early development. Plant Physiol. 134, 1813–1823. doi: 10.1104/pp.103.033001
Mo, Z., Pu, Y., Zhou, J., Tian, Z., Teng, J., Chen, Q., et al. (2021). Effect of the over-dominant expression of proteins on nicotine heterosis via proteomic analysis. Sci. Rep. 11, 21063. doi: 10.1038/s41598-021-00614-x
Muraro, D., Wilson, M., Bennett, M. J. (2011). Root development: cytokinin transport matters, too! Curr. Biol. 21, R423–R425. doi: 10.1016/j.cub.2011.04.048
Pan, J., Li, Z., Wang, Q., Yang, L., Yao, F., Liu, W., et al. (2020). An s-domain receptor-like kinase, osesg1, regulates early crown root development and drought resistance in rice. Plant Sci. 290, 110318. doi: 10.1016/j.plantsci.2019.110318
Paschold, A., Jia, Y., Marcon, C., Lund, S., Larson, N. B., Yeh, C. T., et al. (2012). Complementation contributes to transcriptome complexity in maize (zea mays l.) hybrids relative to their inbred parents. Genome Res. 22, 2445–2454. doi: 10.1101/gr.138461.112
Paschold, A., Marcon, C., Hoecker, N., Hochholdinger, F. (2010). Molecular dissection of heterosis manifestation during early maize root development. Theor. Appl. Genet. 120, 383–388. doi: 10.1007/s00122-009-1082-6
Pelagio-Flores, R., Munoz-Parra, E., Barrera-Ortiz, S., Ortiz-Castro, R., Saenz-Mata, J., Ortega-Amaro, M. A., et al. (2019). The cysteine-rich receptor-like protein kinase crk28 modulates arabidopsis growth and development and influences abscisic acid responses. Planta 251, 2. doi: 10.1007/s00425-019-03296-y
Pi, K., Luo, W., Mo, Z., Duan, L., Ke, Y., Wang, P., et al. (2022). Overdominant expression of related genes of ion homeostasis improves k(+) content advantage in hybrid tobacco leaves. BMC Plant Biol. 22, 335. doi: 10.1186/s12870-022-03719-1
Qin, H., Huang, R. (2018). Auxin controlled by ethylene steers root development. Int. J. Mol. Sci. 19, 3656. doi: 10.3390/ijms19113656
Qin, H., Wang, J., Chen, X., Wang, F., Peng, P., Zhou, Y., et al. (2019). Rice osdof15 contributes to ethylene-inhibited primary root elongation under salt stress. New Phytol. 223, 798–813. doi: 10.1111/nph.15824
Rapp, R. A., Udall, J. A., Wendel, J. F. (2009). Genomic expression dominance in allopolyploids. BMC Biol. 7, 18. doi: 10.1186/1741-7007-7-18
Raskin, I. (1992). Salicylate, a new plant hormone. Plant Physiol. 99, 799–803. doi: 10.1104/pp.99.3.799
Ren, J., Zhang, F., Gao, F., Zeng, L., Lu, X., Zhao, X., et al. (2020). Transcriptome and genome sequencing elucidates the molecular basis for the high yield and good quality of the hybrid rice variety chuanyou6203. Sci. Rep. 10, 19935. doi: 10.1038/s41598-020-76762-3
Rivas, M. A., Friero, I., Alarcon, M. V., Salguero, J. (2022). Auxin-cytokinin balance shapes maize root architecture by controlling primary root elongation and lateral root development. Front. Plant Sci. 13. doi: 10.3389/fpls.2022.836592
Saini, S., Sharma, I., Kaur, N., Pati, P. K. (2013). Auxin: a master regulator in plant root development. Plant Cell Rep. 32, 741–757. doi: 10.1007/s00299-013-1430-5
Santner, A., Calderon-Villalobos, L. I., Estelle, M. (2009). Plant hormones are versatile chemical regulators of plant growth. Nat. Chem. Biol. 5, 301–307. doi: 10.1038/nchembio.165
Shahzad, K., Zhang, X., Guo, L., Qi, T., Bao, L., Zhang, M., et al. (2020). Comparative transcriptome analysis between inbred and hybrids reveals molecular insights into yield heterosis of upland cotton. BMC Plant Biol. 20, 239. doi: 10.1186/s12870-020-02442-z
Shalby, N., Mohamed, I., Xiong, J., Hu, K., Yang, Y., Nishawy, E., et al. (2021). Overdominance at the gene expression level plays a critical role in the hybrid root growth of brassica napus. Int. J. Mol. Sci. 22, 9246. doi: 10.3390/ijms22179246
Shao, L., Xing, F., Xu, C., Zhang, Q., Che, J., Wang, X., et al. (2019). Patterns of genome-wide allele-specific expression in hybrid rice and the implications on the genetic basis of heterosis. Proc. Natl. Acad. Sci. U S 116, 5653–5658. doi: 10.1073/pnas.1820513116
Sheremet, I., Emets, A. I., Vissenberg, K., Verbelen, J., Blium, I. (2010). [the effect of inhibitors of serinethreonine protein kinases on arabidopsis thaliana root morphology and microtubules organization in its cells]. Tsitologiia 52, 389–398.
Shull, G. H. (1908). The composition of a field of maize. J. Heredity 4, 296–301. doi: 10.1093/jhered/os-4.1.296
Solt, D. (1959). Estimated contributions of root and shoot to the nicotine content of the tobacco plant. Plant Physiol. 34, 656–661.
Song, G., Guo, Z., Liu, Z., Cheng, Q., Qu, X., Chen, R., et al. (2013). Global rna sequencing reveals that genotype-dependent allele-specific expression contributes to differential expression in rice f1 hybrids. BMC Plant Biol. 13, 221. doi: 10.1186/1471-2229-13-221
Song, W., Wang, F., Chen, L., Ma, R., Zuo, X., Cao, A., et al. (2019). Ghvtc1, the key gene for ascorbate biosynthesis in gossypium hirsutum, involves in cell elongation under control of ethylene. Cells 8, 1039. doi: 10.3390/cells8091039
Stortenbeker, N., Bemer, M. (2019). The saur gene family: the plant's toolbox for adaptation of growth and development. J. Exp. Bot. 70, 17–27. doi: 10.1093/jxb/ery332
Thiemann, A., Fu, J., Seifert, F., Grant-Downton, R. T., Schrag, T. A., Pospisil, H., et al. (2014). Genome-wide meta-analysis of maize heterosis reveals the potential role of additive gene expression at pericentromeric loci. BMC Plant Biol. 14, 88. doi: 10.1186/1471-2229-14-88
Tian, M., Nie, Q., Li, Z., Zhang, J., Liu, Y., Long, Y., et al. (2018). Transcriptomic analysis reveals overdominance playing a critical role in nicotine heterosis in nicotiana tabacum l. BMC Plant Biol. 18, 48. doi: 10.1186/s12870-018-1257-x
To, J. P., Haberer, G., Ferreira, F. J., Deruere, J., Mason, M. G., Schaller, G. E., et al. (2004). Type-a arabidopsis response regulators are partially redundant negative regulators of cytokinin signaling. Plant Cell 16, 658–671. doi: 10.1105/tpc.018978
Tong, Z., Fang, D., Chen, X., Jiao, F., Zhang, Y., Li, Y., et al. (2020). Genome-wide association study of leaf chemistry traits in tobacco. Breed Sci. 70, 253–264. doi: 10.1270/jsbbs.19067
Wang, X., Guo, C., Peng, J., Li, C., Wan, F., Zhang, S., et al. (2019). Abre-binding factors play a role in the feedback regulation of aba signaling by mediating rapid aba induction of aba co-receptor genes. New Phytol. 221, 341–355. doi: 10.1111/nph.15345
Wang, L., Liu, P. C., Wu, L. M., Tan, J., Peacock, W. J., Dennis, E. S. (2019). Cotyledons contribute to plant growth and hybrid vigor in arabidopsis. Planta 249, 1107–1118. doi: 10.1007/s00425-018-3068-6
Xiong, J., Hu, K., Shalby, N., Zhuo, C., Wen, J., Yi, B., et al. (2022). Comparative transcriptomic analysis reveals the molecular mechanism underlying seedling biomass heterosis in brassica napus. BMC Plant Biol. 22, 283. doi: 10.1186/s12870-022-03671-0
Xu, P., Zhao, P. X., Cai, X. T., Mao, J. L., Miao, Z. Q., Xiang, C. B., et al. (2020). Integration of jasmonic acid and ethylene into auxin signaling in root development. Front. Plant Sci. 11, 271. doi: 10.3389/fpls.2020.00271
Yang, G., Chen, S., Wang, S., Liu, G., Li, H., Huang, H., et al. (2015). Bpgh3.5, an early auxin-response gene, regulates root elongation in betula platyphylla × betula pendula. Plant Cell Tissue Organ Culture (Pctoc) 120, 239–250. doi: 10.1007/s11240-014-0599-9
Yang, L., Shao, H., Zhang, X., Yu, J., Ma, X. (2002). Research progress of tobacco root system. Tobacco Technol. 10, 45–48.
Zhang, C., Lin, C., Fu, F., Zhong, X., Peng, B., Yan, H., et al. (2017). Comparative transcriptome analysis of flower heterosis in two soybean f1 hybrids by rna-seq. PloS One 12, e181061. doi: 10.1371/journal.pone.0181061
Zhang, C., Li, X., Wang, Z., Zhang, Z., Wu, Z. (2020). Identifying key regulatory genes of maize root growth and development by rna sequencing. Genomics 112, 5157–5169. doi: 10.1016/j.ygeno.2020.09.030
Zhang, Q., Li, Y., Xu, T., Srivastava, A. K., Wang, D., Zeng, L., et al. (2016). The chromatin remodeler ddm1 promotes hybrid vigor by regulating salicylic acid metabolism. Cell Discovery 2, 16027. doi: 10.1038/celldisc.2016.27
Zhaojun, D. (2021). The current and future studies on plant root development and root microbiota. Chin. Sci.: Life Sci. 51, 1447–1456.
Zhu, A., Greaves, I. K., Liu, P. C., Wu, L., Dennis, E. S., Peacock, W. J. (2016). Early changes of gene activity in developing seedlings of arabidopsis hybrids relative to parents may contribute to hybrid vigour. Plant J. 88, 597–607. doi: 10.1111/tpj.13285
Zhu, D., Zhou, G., Xu, C., Zhang, Q. (2016). Genetic components of heterosis for seedling traits in an elite rice hybrid analyzed using an immortalized f2 population. J. Genet. Genomics 43, 87–97. doi: 10.1016/j.jgg.2016.01.002
Keywords: heterosis, overdominant expression, root system, transcriptomics, WGCNA
Citation: Pi K, Huang Y, Luo W, Zeng S, Mo Z, Duan L and Liu R (2023) Overdominant expression of genes plays a key role in root growth of tobacco hybrids. Front. Plant Sci. 14:1107550. doi: 10.3389/fpls.2023.1107550
Received: 25 November 2022; Accepted: 17 January 2023;
Published: 31 January 2023.
Edited by:
Helena Fernández, University of Oviedo, SpainReviewed by:
Mengnan An, Shenyang Agricultural University, ChinaIldiko Matusikova, University of St. Cyril and Methodius, Slovakia
Copyright © 2023 Pi, Huang, Luo, Zeng, Mo, Duan and Liu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Renxiang Liu, rxliu@gzu.edu.cn