Transcriptome analysis of rice root heterosis by RNA-Seq

Transcriptome Analysis Revealed the Paternal Importance to Vegetative Growth Heterosis in Populus

Parental selection is important for heterosis formation during crossbreeding of Populus. However, in poplar hybrids, the effect of parents on vegetative growth heterosis is not well understood. In this study, one female parent (P. simonii XY4) and two male parents (P. nigra OH and P. deltoides × P. nigra BJLY3#) were used to produce two progenies (Hyb1 and Hyb2). Vegetative growth investigation showed that both Hyb1 and Hyb2 performed heterosis in plant growth and ground diameter. The vegetative growth of hybrids was strongly correlated with the male parents but not with the female parents. The gene expression levels in the hybrids were more biased toward the male parents. In Hyb1 and Hyb2, 51.93% and 45.03% of the expressed genes showed the non-additive effect, respectively, and over 65% of the non-additively expressed genes showed the dominant effect. It is noteworthy that genes of paternal expression dominant effect (ELD_♂) account for the majority of dominantly expressed genes, suggesting the paternal contribution to heterosis. KEGG enrichment analysis indicated that a large number of non-additively expressed genes were enriched in the plant hormone signal transduction pathway. WGCNA analysis showed that MEcyan was significantly correlated with the traits of hybrids, and 12 plant hormone signal transduction pathway genes were enriched in this module. Transcription factors (TFs) MYB88, LHY, and TCP4 may be involved in the regulation of these pathway genes. This finding supported that the male parents play an important role in the formation of vegetative growth heterosis of Populus. In addition, the non-additively expressed genes of the signal transduction pathway and the regulation of TFs related to these pathway genes may be one of the reasons for the generation of heterosis.

hybridexpress: an R/Bioconductor package for comparative transcriptomic analyses of hybrids and their progenitors

Hybridization, the process of crossing individuals from diverse genetic backgrounds, plays a pivotal role in evolution, biological invasiveness, and crop breeding. At the transcriptional level, hybridization often leads to complex nonadditive effects, presenting challenges for understanding its consequences. Although standard transcriptomic analyses exist to compare hybrids to their progenitors, such analyses have not been implemented in a software package, hindering reproducibility. We introduce hybridexpress, an R/Bioconductor package designed to facilitate the analysis, visualization, and comparison of gene expression patterns in hybrid triplets (hybrids and their progenitors). hybridexpress provides users with a user-friendly and comprehensive workflow that includes all standard comparative analyses steps, including data normalization, calculation of midparent expression values, sample clustering, expression-based gene classification into categories and classes, and overrepresentation analysis for functional terms. We illustrate the utility of hybridexpress through comparative transcriptomic analyses of cotton allopolyploidization and rice root trait heterosis. hybridexpress is designed to streamline comparative transcriptomic studies of hybrid triplets, advancing our understanding of evolutionary dynamics in allopolyploids, and enhancing plant breeding strategies. hybridexpress is freely accessible from Bioconductor (https://bioconductor.org/packages/HybridExpress) and its source code is available on GitHub (https://github.com/almeidasilvaf/HybridExpress).

Comparative transcriptome analysis provides molecular insights into heterosis of waterlogging tolerance in Chrysanthemum indicum

BACKGROUND

Heterosis breeding is one of the most important breeding methods for chrysanthemum. To date, the genetic mechanisms of heterosis for waterlogging tolerance in chrysanthemum are still unclear. This study aims to analyze the expression profiles and potential heterosis-related genes of two hybrid lines and their parents with extreme differences in waterlogging tolerance under control and waterlogging stress conditions by RNA-seq.

RESULTS

A population of 140 F1 progeny derived from Chrysanthemum indicum (Nanchang) (waterlogging-tolerant) and Chrysanthemum indicum (Nanjing) (waterlogging-sensitive) was used to characterize the extent of genetic variation in terms of seven waterlogging tolerance-related traits across two years. Lines 98 and 95, respectively displaying positive and negative overdominance heterosis for the waterlogging tolerance traits together with their parents under control and waterlogging stress conditions, were used for RNA-seq. In consequence, the maximal number of differentially expressed genes (DEGs) occurred in line 98. Gene ontology (GO) enrichment analysis revealed multiple stress-related biological processes for the common up-regulated genes. Line 98 had a significant increase in non-additive genes under waterlogging stress, with transgressive up-regulation and paternal-expression dominant patterns being the major gene expression profiles. Further, GO analysis identified 55 and 95 transgressive up-regulation genes that overlapped with the up-regulated genes shared by two parents in terms of responses to stress and stimulus, respectively. 6,640 genes in total displaying maternal-expression dominance patterns were observed in line 95. In addition, 16 key candidate genes, including SAP12, DOX1, and ERF017 which might be of significant importance for the formation of waterlogging tolerance heterosis in line 98, were highlighted.

CONCLUSION

The current study provides a comprehensive overview of the root transcriptomes among F1 hybrids and their parents under waterlogging stress. These findings lay the foundation for further studies on molecular mechanisms underlying chrysanthemum heterosis on waterlogging tolerance.

RNA Sequencing Analysis and Verification of Paeonia ostii 'Fengdan' CuZn Superoxide Dismutase (PoSOD) Genes in Root Development

Tree peony (Paeonia suffruticosa) is a significant medicinal plant. However, the low rooting number is a bottleneck problem in the micropropagation protocols of P. ostii 'Fengdan'. The activity of superoxide dismutase (SOD) is closely related to root development. But research on the SOD gene's impact on rooting is still lacking. In this study, RNA sequencing (RNA-seq) was used to analyze the four crucial stages of root development in P. ostii 'Fengdan' seedlings, including the early root primordium formation stage (Gmfq), root primordium formation stage (Gmf), root protrusion stage (Gtq), and root outgrowth stage (Gzc). A total of 141.77 GB of data were obtained; 71,718, 29,804, and 24,712 differentially expressed genes (DEGs) were identified in the comparison groups of Gmfq vs. Gmf, Gmf vs. Gtq, and Gtq vs. Gzc, respectively. Among the 20 most highly expressed DEGs in the three comparison groups, only the CuZnSOD gene (SUB13202229, PoSOD) was found to be significantly expressed in Gtq vs. Gzc. The overexpression of PoSOD increased the number of adventitious roots and promoted the activities of peroxidase (POD) and SOD in P. ostii 'Fengdan'. The gene ADVENTITIOUS ROOTING RELATED OXYGENASE1 (PoARRO-1), which is closely associated with the development of adventitious roots, was also significantly upregulated in overexpressing PoSOD plants. Furthermore, PoSOD interacted with PoARRO-1 in yeast two-hybrid (Y2H) and biomolecular luminescence complementation (BiFC) assays. In conclusion, PoSOD could interact with PoARRO-1 and enhance the root development of tube plantlets in P. ostii 'Fengdan'. This study will help us to preliminarily understand the molecular mechanism of adventitious root formation and improve the root quality of tree peony and other medicinal plants.

Gene Expression Analysis of Different Organs and Identification of AP2 Transcription Factors in Flax (Linum usitatissimum L.)

Flax (Linum usitatissimum L.) is an important oilseed crop widely cultivated for its oil and fiber. This study conducted transcriptome analysis to analyze the gene expression profiles of roots, leaves, stamens, pistils, and fruits in the flax cultivar Longya10. A total of 43,471 genes were detected in the RNA-seq data, with 34,497 genes showing differential expression levels between different organs. Gene expression patterns varied across different organs, with differences observed in expression-regulating genes within specific organs. However, 23,448 genes were found to be commonly expressed across all organs. Further analysis revealed organ-specific gene expressions, with 236, 690, 544, 909, and 1212 genes identified in pistils, fruits, leaves, roots, and stamens, respectively. Gene Ontology (GO) enrichment analysis was performed on these organ-specific genes, and significant enrichment was observed in various biological processes, cellular components, and molecular functions, providing new insights for the specific growth patterns of flax organs. Furthermore, we investigated the expression differences of AP2 transcription factors in various tissues and organs of Longya10. We identified 96 AP2 genes that were differentially expressed in different organs and annotated them into various biological pathways. Our results suggest that AP2 transcription factors may play important roles in regulating the growth and development of flax organs including stress response. In summary, our study provides a comprehensive analysis of gene expression patterns in different organs and tissues of flax plant and identifies potential critical regulators of flax organ growth and development. These findings contribute to a better understanding of the molecular mechanisms underlying flax organ development and may have important implications for the genetic improvement of flax crops.

Transcriptomics and metabolomics reveal the changes induced by arbuscular mycorrhizal fungi in Panax quinquefolius L

BACKGROUND

Panax quinquefolius L. is one of the most important foods and herbs because of its high nutritional value and medicinal potential. In our previous study we found that the ginsenoside content in P. quinquefolius was improved by arbuscular mycorrhizal fungi (AMFs). However, little research has been conducted on the molecular mechanisms in P. quinquefolius roots induced by AMFs colonization. To identify the metabolomic and transcriptomic mechanisms of P. quinquefolius induced by AMFs, non-mycorrhized (control) and mycorrhized (AMF) P. quinquefolius were used as experimental materials for comparative analysis of the transcriptome and metabolome.

RESULTS

Compared with the control, 182 metabolites and 545 genes were significantly changed at the metabolic and transcriptional levels in AMFs treatment. The metabolic pattern of AMFs was changed, and the contents of ginsenosides (Rb1, Rg2), threonine, and glutaric acid were significantly increased. There were significant differences in the expression of genes involved in plant hormone signal transduction, glutathione metabolism, and the plant-pathogen interaction pathway. In addition, several transcription factors from the NAC, WRKY, and basic helix-loop-helix families were identified in AMFs versus the control. Furthermore, the combined analysis of 'transcriptomics-metabolomics' analysis showed that 'Plant hormone signal transduction', 'Amino sugar and nucleotide sugar metabolism' and 'Glutathione metabolism' pathways were the important enriched pathways in response to AMFs colonization.

CONCLUSION

Overall, these results provide new insights into P. quinquefolius response to AMFs, which improve our understanding of the molecular mechanisms of P. quinquefolius induced by AMFs. © 2023 Society of Chemical Industry.

Transcriptomic and metabolomic analyses unveil the growth advantage mechanism conferred by heterosis of Michelia 'Zhongshanhanxiao'

Michelia compressa (Maxim.) Sarg. is one of the important timber trees in Taiwan province, P. R. China. Michelia 'Zhongshanhanxiao' is a group of variants found among the progeny of M. compressa that exhibit higher growth rates compared with normal individuals, with a significantly increased stem diameter and height, as well as enlarged leaves and flowers. However, the molecular mechanisms fostering the growth advantage and morphological variations are unknown and deserve further study. Through analysing the transcriptome, metabolome and physiological processes of leaves, we identified remarkable differences in gene expression and metabolic profiles between Michelia 'Zhongshanhanxiao' and both the maternal M. compressa and its normal progeny. These differences were widely associated with a plant-pathogen interaction, phenylpropanoid biosynthesis, cyanoamino acid metabolism, carbon fixation in photosynthetic organisms and plant hormone signal transduction. Additionally, physiological measurements showed that Michelia 'Zhongshanhanxiao' possesses stronger photosynthetic capacity and higher plant hormone content. These results suggest that the heterosis of Michelia 'Zhongshanhanxiao' is regulated by candidates related to cell division, resistance to pathogens and the accumulation of organic compounds. The findings of this study provide crucial information on the molecular mechanisms underlying the growth advantages conferred by heterosis in trees.

Molecular dissection of heterosis in cereal roots and their rhizosphere

Heterosis is already manifested early in root development. Consistent with the dominance model of heterosis, gene expression complementation is a general mechanism that contributes to phenotypic heterosis in maize hybrids. Highly heterozygous F1-hybrids outperform their parental inbred lines, a phenomenon known as heterosis. Utilization of heterosis is of paramount agricultural importance and has been widely applied to increase yield in many crop cultivars. Plant roots display heterosis for many traits and are an important target for further crop improvement. To explain the molecular basis of heterosis, several genetic hypotheses have been proposed. In recent years, high-throughput gene expression profiling techniques have been applied to investigate hybrid vigor. Consistent with the classical genetic dominance model, gene expression complementation has been demonstrated to be a general mechanism to contribute to phenotypic heterosis in diverse maize hybrids. Functional classification of these genes supported the notion that gene expression complementation can dynamically promote hybrid vigor under fluctuating environmental conditions. Hybrids tend to respond differently to available nutrients in the soil. It was hypothesized that hybrid vigor is promoted through a higher nutrient use efficiency which is linked to an improved root system performance of hybrids in comparison to their inbred parents. Recently, the interaction between soil microbes and their plant host was added as further dimension to disentangle heterosis in the belowground part of plants. Soil microbes influenced the performance of maize hybrids as illustrated in comparisons of sterile soil and soil inhabited by beneficial microorganisms.

CsBPC2 is a key regulator of root growth and development

BASIC PENTACYSTEINE (BPCs) transcription factors are important regulators of plant growth and development. However, the regulatory mechanism of BPC2 in roots remains unclear. In our previous study, we created Csbpc2 cucumber mutants by the CRISPR/Cas9 system, and our studies on the phenotype of Csbpc2 mutants showed that the root growth was inhibited compared with wide-type (WT). Moreover, the surface area, volume and number of roots decreased significantly, with root system architecture changing from dichotomous branching to herringbone branching. Compared with WT, the leaf growth of the Csbpc2 mutants was not affected. However, the palisade and spongy tissue were significantly thinner, which was not beneficial for photosynthesis. The metabolome of root exudates showed that compared with WT, amino acids and their derivatives were significantly decreased, and the enriched pathways were mainly regulated by amino acids and their derivatives, indicating that knockout of CsBPC2 mainly affected the amino acid content in root exudates. Importantly, transcriptome analysis showed that knockout of CsBPC2 mainly affected root gene expression. Knockout of CsBPC2 significantly reduced the gene expression of gibberellins synthesis. However, the expression of genes related to amino acid synthesis, nitrogen fixation and PSII-related photosynthesis increased significantly, which may be due to the effect of knocking out CsBPC2 on gibberellins synthesis, resulting in the inhibition of seedling growth, thus forming negative feedback regulation. Generally, we showed for the first time that BPC2 is a key regulator gene of root growth and development, laying the foundation for future mechanisms of BPC2 regulation in roots.

Genetic Analysis for the Flag Leaf Heterosis of a Super-Hybrid Rice WFYT025 Combination Using RNA-Seq

The photosynthetic capacity of flag leaf plays a key role in grain yield in rice. Nevertheless, there are few studies on the heterosis of the rice flag leaf. Therefore, this study focuses on investigating the genetic basis of heterosis for flag leaf in the indica super hybrid rice combination WFYT025 in China using a high-throughput next-generation RNA-seq strategy. We analyzed the gene expression of flag leaf in different environments and different time periods between WFYT025 and its female parent. After obtaining the gene expression profile of the flag leaf, we further investigated the gene regulatory network. Weighted gene expression network analysis (WGCNA) was used to identify the co-expressed gene sets, and a total of 5000 highly expressed genes were divided into 24 co-expression groups. In CHT025, we found 13 WRKY family transcription factors in SDG_hps under the environment of early rice and 16 WRKY family genes in SDG_hps of under the environment of middle rice. We found nine identical transcription factors in the two stages. Except for five reported TFs, the other four TFs might play an important role in heterosis for grain number and photosynthesis. Transcription factors such as WRKY3, WRKY68, and WRKY77 were found in both environments. To eliminate the influence of the environment, we examined the metabolic pathway with the same SDG_hp (SSDG_hp) in two environments. There were 312 SSDG_hps in total. These SSDG_hps mainly focused on the phosphorus metallic process, phosphorylation, plasma membrane, etc. These results provide resources for studying heterosis during super hybrid rice flag leaf development.

Dissecting the genetic basis of heterosis in elite super-hybrid rice

Y900 is one of the top hybrid rice (Oryza sativa) varieties, with its yield exceeding 15 t·hm-2. To dissect the mechanism of heterosis, we sequenced the male parent line R900 and female parent line Y58S using long-read and Hi-C technology. High-quality reference genomes of 396.41 Mb and 398.24 Mb were obtained for R900 and Y58S, respectively. Genome-wide variations between the parents were systematically identified, including 1,367,758 single-nucleotide polymorphisms, 299,149 insertions/deletions, and 4,757 structural variations. The level of variation between Y58S and R900 was the lowest among the comparisons of Y58S with other rice genomes. More than 75% of genes exhibited variation between the two parents. Compared with other two-line hybrids sharing the same female parent, the portion of Geng/japonica (GJ)-type genetic components from different male parents increased with yield increasing in their corresponding hybrids. Transcriptome analysis revealed that the partial dominance effect was the main genetic effect that constituted the heterosis of Y900. In the hybrid, both alleles from the two parents were expressed, and their expression patterns were dynamically regulated in different tissues. The cis-regulation was dominant for young panicle tissues, while trans-regulation was more common in leaf tissues. Overdominance was surprisingly prevalent in stems and more likely regulated by the trans-regulation mechanism. Additionally, R900 contained many excellent GJ haplotypes, such as NARROW LEAF1, Oryza sativa SQUAMOSA PROMOTER BINDING PROTEIN-LIKE13, and Grain number, plant height, and heading date8, making it a good complement to Y58S. The fine-tuned mechanism of heterosis involves genome-wide variation, GJ introgression, key functional genes, and dynamic gene/allele expression and regulation pattern changes in different tissues and growth stages.

Transcriptome Analysis of Reciprocal Hybrids Between Crassostrea gigas and C. angulata Reveals the Potential Mechanisms Underlying Thermo-Resistant Heterosis

Heterosis, also known as hybrid vigor, is widely used in aquaculture, but the molecular causes for this phenomenon remain obscure. Here, we conducted a transcriptome analysis to unveil the gene expression patterns and molecular bases underlying thermo-resistant heterosis in Crassostrea gigas ♀ × Crassostrea angulata ♂ (GA) and C. angulata ♀ × C. gigas ♂ (AG). About 505 million clean reads were obtained, and 38,210 genes were identified, of which 3779 genes were differentially expressed between the reciprocal hybrids and purebreds. The global gene expression levels were toward the C. gigas genome in the reciprocal hybrids. In GA and AG, 95.69% and 92.00% of the differentially expressed genes (DEGs) exhibited a non-additive expression pattern, respectively. We observed all gene expression modes, including additive, partial dominance, high and low dominance, and under- and over-dominance. Of these, 77.52% and 50.00% of the DEGs exhibited under- or over-dominance in GA and AG, respectively. The over-dominance DEGs common to reciprocal hybrids were significantly enriched in protein folding, protein refolding, and intrinsic apoptotic signaling pathway, while the under-dominance DEGs were significantly enriched in cell cycle. As possible candidate genes for thermo-resistant heterosis, GRP78, major egg antigen, BAG, Hsp70, and Hsp27 were over-dominantly expressed, while MCM6 and ANAPC4 were under-dominantly expressed. This study extends our understanding of the thermo-resistant heterosis in oysters.

Transcriptome Analysis Reveals a Major Gene Expression Pattern and Important Metabolic Pathways in the Control of Heterosis in Chinese Cabbage

Although heterosis is commonly used in Chinese cabbage, its molecular basis is poorly understood. In this study, 16Chinese cabbage hybrids were utilized as test subjects to explore the potential molecular mechanism of heterosis. RNA sequencing revealed 5815-10,252 differentially expressed genes (DEGs) (female parent vs. male parent), 1796-5990 DEGs (female parent-vs-hybrid), and 2244-7063 DEGs (male parent vs. hybrid) in 16 cross combinations at the middle stage of heading. Among of them, 72.83-84.20% DEGs conformed to the dominant expression pattern, which is the predominant expression pattern in hybrids. There were 13 pathways in which DEGs were significantly enriched in most cross combinations. Among them, the plant-pathogen interaction (ko04626) and circadian rhythm-plant (ko04712)were significantly enriched by DEGs in strong heterosis hybrids. WGCNA also proved that the two pathways were significantly related to heterosis in Chinese cabbage.

Transcriptomic profiling reveals candidate allelopathic genes in rice responsible for interactions with barnyardgrass

Echinochloa crus-galli (barnyardgrass) is one of the most damaging weeds in rice fields worldwide. Allelopathy has been considered a possible application for weed management. Thus understanding its molecular mechanisms is important for rice production. This study generated transcriptomes from rice under mono- and co-culture with barnyardgrass at two-time points to identify the candidate genes controlling allelopathic interactions between rice and barnyardgrass. A total of 5,684 differentially expressed genes (DEGs) were detected, amongst which 388 genes were transcription factors. These DEGs include genes associated with momilactone and phenolic acid biosynthesis, which play critical roles in allelopathy. Additionally, we found significantly more DEGs at 3 hours than at 3 days, suggesting a quick allelopathic response in rice. Up-regulated DEGs involve diverse biological processes, such as response to stimulus and pathways related to phenylpropanoid and secondary metabolites biosynthesis. Down-regulated DEGs were involved in developmental processes, indicating a balance between growth and stress response to allelopathy from barnyardgrass. Comparison of DEGs between rice and barnyardgrass shows few common genes, suggesting different mechanisms underlying allelopathic interaction in these two species. Our results offer an important basis for identifying of candidate genes responsible for rice and barnyardgrass interactions and contribute valuable resources for revealing its molecular mechanisms.

Comparative Transcriptome Analysis of the Heterosis of Salt Tolerance in Inter-Subspecific Hybrid Rice

Soil salinity is one of the major abiotic stresses limiting rice growth. Hybrids outperform their parents in salt tolerance in rice, while its mechanism is not completely understood. In this study, a higher seedling survival was observed after salt treatment in an inter-subspecific hybrid rice, Zhegengyou1578 (ZGY1578), compared with its maternal japonica Zhegeng7A (ZG7A) and paternal indica Zhehui1578 (ZH1578). A total of 2584 and 3061 differentially expressed genes (DEGs) with at least twofold changes were identified between ZGY1578 and ZG7A and between ZGY1578 and ZH1578, respectively, in roots under salt stress using the RNA sequencing (RNA-Seq) approach. The expressions of a larger number of DEGs in hybrid were lower or higher than those of both parents. The DEGs associated with transcription factors, hormones, and reactive oxygen species (ROS)-related genes might be involved in the heterosis of salt tolerance. The expressions of the majority of transcription factors and ethylene-, auxin-, and gibberellin-related genes, as well as peroxidase genes, were significantly higher in the hybrid ZGY1578 compared with those of both parents. The identified genes provide valuable clues to elucidate the heterosis of salt tolerance in inter-subspecific hybrid rice.

Comparative physiological and root transcriptome analysis of two annual ryegrass cultivars under drought stress

Annual ryegrass is a widely cultivated forage grass with rapid growth and high productivity. However, drought is one of the abiotic stresses affecting ryegrass growth and quality. In this study, we compared the physiological and transcriptome responses of Chuansi No.1 (drought-tolerant, DT) and Double Barrel (drought-sensitive, DS) under drought stress simulated by PEG-6000 for 7 days. The results showed that Chuansi No. 1 had stronger physiological and biochemical parameters such as root properties, water content, osmotic adjustment ability and antioxidant ability. In addition, RNA-seq was used to elucidate the molecular mechanism of root drought resistance. We identified 8588 differentially expressed genes related to drought tolerance in root, which were mainly enriched in oxidation-reduction process, carbohydrate metabolic process, apoplast, arginine and proline metabolism, and phenylpropanoid biosynthesis pathways. The expression levels of DEGs were consistent with physiological changes of ryegrass under drought stress. We found that genes related to sucrose and starch synthesis, root development, osmotic adjustment, ABA signal regulation and specifically up-regulated transcription factors such as WRKY41, WRKY51, ERF7, ERF109, ERF110, NAC43, NAC68, bHLH162 and bHLH148 in Chuansi No. 1 may be the reason for its higher drought tolerance. This study revealed the underlying physiological and molecular mechanisms of root response to drought stress in ryegrass and provided some new candidate genes for breeding rye drought tolerant varieties.

Transcriptome profiling of two super hybrid rice provides insights into the genetic basis of heterosis

BACKGROUND

Heterosis is a phenomenon that hybrids show superior performance over their parents. The successful utilization of heterosis has greatly improved rice productivity, but the molecular basis of heterosis remains largely unclear.

RESULTS

Here, the transcriptomes of young panicles and leaves of the two widely grown two-line super hybrid rice varieties (Jing-Liang-You-Hua-Zhan (JLYHZ) and Long-Liang-You-Hua-Zhan (LLYHZ)) and their parents were analyzed by RNA-seq. Transcriptome profiling of the hybrids revealed 1,778 ~ 9,404 differentially expressed genes (DEGs) in two tissues, which were identified by comparing with their parents. GO, and KEGG enrichment analysis showed that the pathways significantly enriched in both tissues of two hybrids were all related to yield and resistance, like circadian rhythm (GO:0,007,623), response to water deprivation (GO:0,009,414), and photosynthetic genes (osa00196). Allele-specific expression genes (ASEGs) were also identified in hybrids. The ASEGs were most significantly enriched in ionotropic glutamate receptor signaling pathway, which was hypothesized to be potential amino acid sensors in plants. Moreover, the ASEGs were also differentially expressed between parents. The number of variations in ASEGs is higher than expected, especially for large effect variations. The DEGs and ASEGs are the potential reasons for the formation of heterosis in the two elite super hybrid rice.

CONCLUSIONS

Our results provide a comprehensive understanding of the heterosis of two-line super hybrid rice and facilitate the exploitation of heterosis in hybrid rice breeding with high yield heterosis.

Twenty years of rice genomics research: From sequencing and functional genomics to quantitative genomics

Since the completion of the rice genome sequencing project in 2005, we have entered the era of rice genomics, which is still in its ascendancy. Rice genomics studies can be classified into three stages: structural genomics, functional genomics, and quantitative genomics. Structural genomics refers primarily to genome sequencing for the construction of a complete map of rice genome sequence. This is fundamental for rice genetics and molecular biology research. Functional genomics aims to decode the functions of rice genes. Quantitative genomics is large-scale sequence- and statistics-based research to define the quantitative traits and genetic features of rice populations. Rice genomics has been a transformative influence on rice biological research and contributes significantly to rice breeding, making rice a good model plant for studying crop sciences.

From hybrid genomes to heterotic trait output: Challenges and opportunities

Heterosis (or hybrid vigor) has been widely used in crop seed breeding to improve many key economic traits. Nevertheless, the genetic and molecular basis of this important phenomenon has long remained elusive, constraining its flexible and effective exploitation. Advanced genomic approaches are efficient in characterizing the mechanism of heterosis. Here, we review how the omics approaches, including genomic, transcriptomic, and population genetics methods such as genome-wide association studies, can reveal how hybrid genomes outperform parental genomes in plants. This information opens up opportunities for genomic exploration and manipulation of heterosis in crop breeding.

Full-length transcriptome reconstruction reveals genetic differences in hybrids of Oryza sativa and Oryza punctata with different ploidy and genome compositions

BACKGROUND

Allopolyploid breeding is an efficient technique for improving the low seed setting rate of autotetraploids in plant breeding and one of the most promising breeding methods. However, there have been few comprehensive studies of the posttranscriptional mechanism in allopolyploids.

RESULTS

By crossing cultivated rice (Oryza sativa, genome AA) with wild rice (Oryza punctata, genome BB), we created hybrid rice lines with different ploidy and genome compositions [diploid hybrid F01 (AB), allotetraploid hybrid F02 (AABB) and F03 (AAAB)]. The genetic differences of the hybrids and the mechanism of allopolyploid breeding dominance were revealed through morphological and cytological observations and single molecule real-time sequencing techniques. The tissues and organs of allotetraploid hybrid F02 exhibited "gigantism" and the highest levels of fertility. The numbers of non-redundant transcripts, gene loci and new isoforms in the polyploid rice lines were higher and the isoform lengths greater than those of the diploid line. Moreover, alternative splicing (AS) events occurred twice as often in the polyploid rice lines than the diploid line. During these events, intron retention dominated. Furthermore, a large number of new genes and isoforms specific to the lines of different ploidy were discovered.

CONCLUSIONS

The results indicated that alternative polyadenylation (APA) and AS events contributed to the complexity and superiority of polyploids in the activity of translation regulators, nucleic acid binding transcription factor activities and the regulation of molecular function. Therefore, these APA and AS events in allopolyploid rice were found to play a role in regulation. Our study provides new germplasm for polyploid rice breeding and reveals complex regulatory mechanisms that may be related to heterosis and fertility.

Comprehensive Identification and Profiling of miRNAs Involved in Terpenoid Synthesis of Gleditsia sinensis Lam

Gleditsia sinensis Lam. is a tree with worldwide distribution and important economic and medicinal values; its pods contain terpenoids including gleditsioside, thiamine, and brassinosteroids. However, thus far, there are few studies on the terpenoid regulation of G. sinensis at the molecular level. microRNA (miRNA) is a class of small RNAs with conserved and crucial roles in the regulation of diverse biological processes during plant growth and development. To identify the miRNAs of G. sinensis and evaluate their involvement in terpenoid synthesis, this investigation quantified the content changes in saponins in pods at three developmental stages: May (pod-setting stage), July (elongation stage), and September (browning stage), and then we performed genome-wide miRNA profiles during the three development stages of the G. sinensis pods. A total of 351 conserved miRNAs belonging to 216 families were identified, among which 36 conserved miRNAs exist specifically in legumes. Through target analysis, 708 unigenes were predicted to be candidate targets of 37 differentially expressed miRNAs. The targets of miR838-3p and miR2093-5p were involved in the derived branches of monoterpenes and gleditsioside, in brassinosteroid biosynthesis (BRB), and in indole alkaloid biosynthesis (IAB). Intriguingly, the targets of miR829-3p.1 were predicted to take part in thiamine biosynthesis, and the targets of miR4414b and miR5037a were involved in the main process of cytokinin synthesis. The corresponding targets participated in BRB, IAB, and terpenoid backbone biosynthesis, which were enriched significantly, suggesting that miR2093-5p, miR4414b, miR5037a, miR829-3p.1, and miR838-3p play indispensable roles in the regulation of triterpenoid saponin and monoterpenoid biosynthesis. To date, this is the first report of miRNA identification in G. sinensis and miRNA expression profiles at different developmental stages of G. sinensis pods, which provides a basis for further uncovering the molecular regulation of terpenoid synthesis in G. sinensis and new insights into the role of miRNAs in legumes.

Combined full-length transcriptomic and metabolomic analysis reveals the regulatory mechanisms of adaptation to salt stress in asparagus

Soil salinity is a very serious abiotic stressor that affects plant growth and threatens crop yield. Thus, it is important to explore the mechanisms of salt tolerance of plant and then to stabilize and improve crop yield. Asparagus is an important cash crop, but its salt tolerance mechanisms are largely unknown. Full-length transcriptomic and metabolomic analyses were performed on two asparagus genotypes: 'jx1502' (a salt-tolerant genotype) and 'gold crown' (a salt-sensitive genotype). Compared with the distilled water treatment (control), 877 and 1610 differentially expressed genes (DEGs) were identified in 'jx1502' and 'gold crown' under salt stress treatment, respectively, and 135 and 73 differentially accumulated metabolites (DAMs) were identified in 'jx1502' and 'gold crown' under salt stress treatment, respectively. DEGs related to ion transport, plant hormone response, and cell division and growth presented differential expression profiles between 'jx1502' and 'gold crown.' In 'jx1502,' 11 ion transport-related DEGs, 8 plant hormone response-related DEGs, and 12 cell division and growth-related DEGs were upregulated, while 7 ion transport-related DEGs, 4 plant hormone response-related DEGs, and 2 cell division and growth-related DEGs were downregulated. Interestingly, in 'gold crown,' 14 ion transport-related DEGs, 2 plant hormone response-related DEGs, and 6 cell division and growth-related DEGs were upregulated, while 45 ion transport-related DEGs, 13 plant hormone response-related DEGs, and 16 cell division and growth-related DEGs were downregulated. Genotype 'jx1502' can modulate K⁺/Na⁺ and water homeostasis and maintain a more constant transport system for nutrient uptake and distribution than 'gold crown' under salt stress. Genotype 'jx1502' strengthened the response to auxin (IAA), as well as cell division and growth for root remodeling and thus salt tolerance. Therefore, the integration analysis of transcriptomic and metabolomic indicated that 'jx1502' enhanced sugar and amino acid metabolism for energy supply and osmotic regulatory substance accumulation to meet the demands of protective mechanisms against salt stress. This work contributed to reveal the underlying salt tolerance mechanism of asparagus at transcription and metabolism level and proposed new directions for asparagus variety improvement.

Transcriptome and DNA methylome analyses provide insight into the heterosis in flag leaf of inter-subspecific hybrid rice

Flag leaf heterosis of inter-subspecific hybrid rice is suggested to be related to leaf area, gene expression pattern and allele-specific expression, putatively related to DNA methylation differences between the hybrid and its parents. Inter-subspecific hybrid rice combinations of indica × japonica have great potential to broaden genetic diversity and enhance the heterosis. The genetic and epigenetic molecular mechanism of its heterosis is not completely understood. Here, the dissection of gene expression and epigenetic regulation of an elite inter-subspecific hybrid rice were reported. In the hybrid, plant height, flag leaf area and Pn showed significant heterosis at the heading stage. Chloroplast-related differentially expressed genes (DEGs) and 530 allele-specific expression genes in hybrid were identified. Analysis of the genome-wide distribution of DNA methylation (5-methylcytosine, 5mC) and its association with transcription showed that there were variant DNA methylation maps and that the regulation of gene expression levels was negatively regulated by DNA methylation in the inter-subspecific hybrid rice. Differentially methylated DEGs were significantly enriched in photosynthetic functions. Moreover, distinct 5mC sequence contexts and distinct functional elements (promoter/gene body) may have different influences on heterosis related genes. The data identified heterosis related molecular mechanisms in inter-subspecific hybrid rice and suggested that epigenetic changes could extensively influence the flag leaf gene expression and heterosis.

Transcriptomic analysis of salt tolerance-associated genes and diversity analysis using indel markers in yardlong bean (Vigna unguiculata ssp. sesquipedialis)

BACKGROUND

High salinity is a devastating abiotic stresses for crops. To understand the molecular basis of salinity stress in yardlong bean (Vigna unguiculata ssp. sesquipedalis), and to develop robust markers for improving this trait in germplasm, whole transcriptome RNA sequencing (RNA-seq) was conducted to compare the salt-tolerant variety Suzi 41 and salt-sensitive variety Sujiang 1419 under normal and salt stress conditions.

RESULTS

Compared with controls, 417 differentially expressed genes (DEGs) were identified under exposure to high salinity, including 42 up- and 11 down-regulated DEGs in salt-tolerant Suzi 41 and 186 up- and 197 down-regulated genes in salt-sensitive Sujiang 1419, validated by qRT-PCR. DEGs were enriched in "Glycolysis/Gluconeogenesis" (ko00010), "Cutin, suberine and wax biosynthesis" (ko00073), and "phenylpropanoid biosynthesis" (ko00940) in Sujiang 1419, although "cysteine/methionine metabolism" (ko00270) was the only pathway significantly enriched in salt-tolerant Suzi 41. Notably, AP2/ERF, LR48, WRKY, and bHLH family transcription factors (TFs) were up-regulated under high salt conditions. Genetic diversity analysis of 84 yardlong bean accessions using 26 InDel markers developed here could distinguish salt-tolerant and salt-sensitive varieties.

CONCLUSIONS

These findings show a limited set of DEGs, primarily TFs, respond to salinity stress in V. unguiculata, and that these InDels associated with salt-inducible loci are reliable for diversity analysis.

Integrated Analysis of mRNA and Non-coding RNA Transcriptome in Pepper (Capsicum chinense) Hybrid at Seedling and Flowering Stages

Pepper is an important vegetable in the world. In this work, mRNA and ncRNA transcriptome profiles were applied to understand the heterosis effect on the alteration in the gene expression at the seedling and flowering stages between the hybrid and its parents in Capsicum chinense. Our phenotypic data indicated that the hybrid has dominance in leaf area, plant scope, plant height, and fruit-related traits. Kyoto Encyclopedia of Genes and Genomes analysis showed that nine members of the plant hormone signal transduction pathway were upregulated in the seedling and flowering stages of the hybrid, which was supported by weighted gene coexpression network analysis and that BC332_23046 (auxin response factor 8), BC332_18317 (auxin-responsive protein IAA20), BC332_13398 (ethylene-responsive transcription factor), and BC332_27606 (ethylene-responsive transcription factor WIN1) were candidate hub genes, suggesting the important potential role of the plant hormone signal transduction in pepper heterosis. Furthermore, some transcription factor families, including bHLH, MYB, and HSF were greatly over-dominant. We also identified 2,525 long ncRNAs (lncRNAs), 47 micro RNAs (miRNAs), and 71 circle RNAs (circRNAs) in the hybrid. In particular, downregulation of miR156, miR169, and miR369 in the hybrid suggested their relationship with pepper growth vigor. Moreover, we constructed some lncRNA–miRNA–mRNA regulatory networks that showed a multi-dimension to understand the ncRNA relationship with heterosis. These results will provide guidance for a better understanding of the molecular mechanism involved in pepper heterosis.

Overdominance at the Gene Expression Level Plays a Critical Role in the Hybrid Root Growth of Brassica napus

Despite heterosis contributing to genetic improvements in crops, root growth heterosis in rapeseed plants is poorly understood at the molecular level. The current study was performed to discover key differentially expressed genes (DEGs) related to heterosis in two hybrids with contrasting root growth performance (FO; high hybrid and FV; low hybrid) based on analysis of the root heterosis effect. Based on comparative transcriptomic analysis, we believe that the overdominance at the gene expression level plays a critical role in hybrid roots’ early biomass heterosis. Our findings imply that a considerable increase in up-regulation of gene expression underpins heterosis. In the FO hybrid, high expression of DEGs overdominant in the starch/sucrose and galactose metabolic pathways revealed a link between hybrid vigor and root growth. DEGs linked to auxin, cytokinin, brassinosteroids, ethylene, and abscisic acid were also specified, showing that these hormones may enhance mechanisms of root growth and the development in the FO hybrid. Moreover, transcription factors such as MYB, ERF, bHLH, NAC, bZIP, and WRKY are thought to control downstream genes involved in root growth. Overall, this is the first study to provide a better understanding related to the regulation of the molecular mechanism of heterosis, which assists in rapeseed growth and yield improvement.

Rapid Generation and Analysis of a Barley Doubled Haploid Line with Higher Nitrogen Use Efficiency Than Parental Lines by F1 Microspore Embryogenesis

Creating varieties with high nitrogen use efficiency (NUE) is crucial for sustainable agriculture development. In this study, a superior barley doubled haploid line (named DH45) with improved NUE was produced via F₁ microspore embryogenesis with three rounds of screening in different nitrogen levels by hydroponic and field experiments. The molecular mechanisms responsible for the NUE of DH45 surpassing that of its parents were investigated by RNA-seq analysis. A total of 1027 differentially expressed genes (DEGs) were identified that were up- or down-regulated in DH45 under low nitrogen conditions but showed no significant differences in the parents. GO analysis indicated that genes involved in nitrogen compound metabolic processes were significantly enriched in DH45 compared with the parents. KEGG analysis showed the MAPK signaling pathway plant to be highly enriched in DH45 relative to its parents, as well as genes involved in alanine, aspartate and glutamate metabolism, and arginine biosynthesis. In conclusion, our study revealed the potential to fix trait superiority in a line by combining crossing with F₁ microspore culture technologies in future crop breeding and also identified several candidate genes that are expressed in shoots and may enable barley to cope with low-nitrogen stress.

Genome-wide identification and characterization of the OFP gene family in Chinese cabbage (Brassica rapa L. ssp. pekinensis)

Ovate family proteins (OFPs) are a class of proteins with a conserved OVATE domain that contains approximately 70 amino acid residues. OFP proteins are plant-specific transcription factors that participate in regulating plant growth and development and are widely distributed in many plants. Little is known about OFPs in Brassica rapa to date. We identified 29 OFP genes in Brassica rapa and found that they were unevenly distributed on 10 chromosomes. Intron gain events may have occurred during the structural evolution of BraOFP paralogues. Syntenic analysis verified Brassica genome triplication, and whole genome duplication likely contributed to the expansion of the OFP gene family. All BraOFP genes had light responsive- and phytohormone-related cis-acting elements. Expression analysis from RNA-Seq data indicated that there were obvious changes in the expression levels of six OFP genes in the Brassica rapa hybrid, which may contribute to the formation of heterosis. Finally, we found that the paralogous genes had different expression patterns among the hybrid and its parents. These results provide the theoretical basis for the further analysis of the biological functions of OFP genes across the Brassica species.

Beyond gene ontology (GO): using biocuration approach to improve the gene nomenclature and functional annotation of rice S-domain kinase subfamily

The S-domain subfamily of receptor-like kinases (SDRLKs) in plants is poorly characterized. Most members of this subfamily are currently assigned gene function based on the S-locus Receptor Kinase from Brassica that acts as the female determinant of self-incompatibility (SI). However, Brassica like SI mechanisms does not exist in most plants. Thus, automated Gene Ontology (GO) pipelines are not sufficient for functional annotation of SDRLK subfamily members and lead to erroneous association with the GO biological process of SI. Here, we show that manual bio-curation can help to correct and improve the gene annotations and association with relevant biological processes. Using publicly available genomic and transcriptome datasets, we conducted a detailed analysis of the expansion of the rice (Oryza sativa) SDRLK subfamily, the structure of individual genes and proteins, and their expression.The 144-member SDRLK family in rice consists of 82 receptor-like kinases (RLKs) (67 full-length, 15 truncated),12 receptor-like proteins, 14 SD kinases, 26 kinase-like and 10 GnK2 domain-containing kinases and RLKs. Except for nine genes, all other SDRLK family members are transcribed in rice, but they vary in their tissue-specific and stress-response expression profiles. Furthermore, 98 genes show differential expression under biotic stress and 98 genes show differential expression under abiotic stress conditions, but share 81 genes in common.Our analysis led to the identification of candidate genes likely to play important roles in plant development, pathogen resistance, and abiotic stress tolerance. We propose a nomenclature for 144 SDRLK gene family members based on gene/protein conserved structural features, gene expression profiles, and literature review. Our biocuration approach, rooted in the principles of findability, accessibility, interoperability and reusability, sets forth an example of how manual annotation of large-gene families can fill in the knowledge gap that exists due to the implementation of automated GO projections, thereby helping to improve the quality and contents of public databases.

Comparative Transcriptomic Analysis of Gene Expression Inheritance Patterns Associated with Cabbage Head Heterosis

The molecular mechanism of heterosis or hybrid vigor, where F1 hybrids of genetically diverse parents show superior traits compared to their parents, is not well understood. Here, we studied the molecular regulation of heterosis in four F1 cabbage hybrids that showed heterosis for several horticultural traits, including head size and weight. To examine the molecular mechanisms, we performed a global transcriptome profiling in the hybrids and their parents by RNA sequencing. The proportion of genetic variations detected as single nucleotide polymorphisms and small insertion–deletions as well as the numbers of differentially expressed genes indicated a larger role of the female parent than the male parent in the genetic divergence of the hybrids. More than 86% of hybrid gene expressions were non-additive. More than 81% of the genes showing divergent expressions showed dominant inheritance, and more than 56% of these exhibited maternal expression dominance. Gene expression regulation by cis-regulatory mechanisms appears to mediate most of the gene expression divergence in the hybrids; however, trans-regulatory factors appear to have a higher effect compared to cis-regulatory factors on parental expression divergence. These observations bring new insights into the molecular mechanisms of heterosis during the cabbage head development.

The influence of soil drought stress on the leaf transcriptome of faba bean (Vicia faba L.) in the Qinghai-Tibet Plateau

Water deficit has a significant impact on growth, development and yield of fava bean (Vicia fava L.) in arid and semi-arid climates. The aim of this study was to identify differentially expressed genes in the Qinghai 13 genotype under soil drought through leaf transcriptome analysis. A total of 256.95 M clean reads were obtained and assembled into 176334 unigenes, with an average length of 766 bp. A total of 9126 (4439 upregulated and 4687 downregulated) differentially expressed genes (DEGs) were identified in faba bean leaves under soil drought. In total, 324 putative transcription factors were identified and classified as belonging to different transcription factor families. According to GO and KEGG analysis, the soil drought stress-inducible DEGs encoded proteins mainly involved in regulating photosynthesis, osmotic adjustment, detoxification, autophagy and other functions. In addition, a large portion of DEGs appeared to be novel because they could not be annotated in any functional databases, therefore, suggesting a specific response to soil drought in faba bean. Finally, RNA-seq analysis was validated by quantitative reverse-transcription PCR analysis. This work provides comprehensive and valuable information for understanding the molecular mechanisms which faba bean uses to respond to soil drought.

The interplay between microRNA and alternative splicing of linear and circular RNAs in eleven plant species

MOTIVATION

MicroRNA (miRNA) and alternative splicing (AS)-mediated post-transcriptional regulation has been extensively studied in most eukaryotes. However, the interplay between AS and miRNAs has not been explored in plants. To our knowledge, the overall profile of miRNA target sites in circular RNAs (circRNA) generated by alternative back splicing has never been reported previously. To address the challenge, we identified miRNA target sites located in alternatively spliced regions of the linear and circular splice isoforms using the up-to-date single-molecule real-time (SMRT) isoform sequencing (Iso-Seq) and Illumina sequencing data in eleven plant species.

RESULTS

In total, we identified 399 401 and 114 574 AS events from linear and circular RNAs, respectively. Among them, there were 64 781 and 41 146 miRNA target sites located in linear and circular AS region, respectively. In addition, we found 38 913 circRNAs to be overlapping with 45 648 AS events of its own parent isoforms, suggesting circRNA regulation of AS of linear RNAs by forming R-loop with the genomic locus. Here, we present a comprehensive database of miRNA targets in alternatively spliced linear and circRNAs (ASmiR) and a web server for deposition and identification of miRNA target sites located in the alternatively spliced region of linear and circular RNAs. This database is accompanied by an easy-to-use web query interface for meaningful downstream analysis. Plant research community can submit user-defined datasets to the web service to search AS regions harboring small RNA target sites. In conclusion, this study provides an unprecedented resource to understand regulatory relationships between miRNAs and AS in both gymnosperms and angiosperms.

AVAILABILITY AND IMPLEMENTATION

The readily accessible database and web-based tools are available at http://forestry.fafu.edu.cn/bioinfor/db/ASmiR.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Transcriptome Analysis by RNA-Seq Reveals Genes Related to Plant Height in Two Sets of Parent-hybrid Combinations in Easter lily (Lilium longiflorum)

In this study, two different hybrids of Easter lily (Lilium longiflorum), obtained from two cross combinations, along with their four parents were sequenced by high–throughput RNA–sequencing (RNA–Seq) to find out differentially expressed gene in parent-hybrid combinations. The leaf mRNA profiles of two hybrids and their four parents were RNA–sequenced with a view to identify the potential candidate genes related to plant height heterosis. In both cross combinations, based to morphological traits mid–parent heterosis (MPH) was higher than high–parent heterosis (HPH) for plant height, leaf length, and number of flowers whereas HPH was higher than MPH for flowering time. A total of 4,327 differentially expressed genes (DEGs) were identified through RNA–Seq between the hybrids and their parents based on fold changes (FC) ≥ 2 for up– and ≤ –2 for down–regulation. Venn diagram analysis revealed that there were 703 common DEGs in two hybrid combinations, those were either up– or down–regulated. Most of the commonly expressed DEGs exhibited higher non–additive effects especially overdominance (75.9%) rather than additive (19.4%) and dominance (4.76%) effects. Among the 384 functionally annotated DEGs identified through Blast2GO tool, 12 DEGs were up–regulated and 16 of them were down–regulated in a similar fashion in both hybrids as revealed by heat map analysis. These 28 universally expressed DEGs were found to encode different types of proteins and enzymes those might regulate heterosis by modulating growth, development and stress–related functions in lily. In addition, gene ontology (GO) analysis of 260 annotated DEGs revealed that biological process might play dominant role in heterotic expression. In this first report of transcriptome sequencing in Easter lily, the notable universally up-regulated DEGs annotated ABC transporter A family member–like, B3 domain–containing, disease resistance RPP13/1, auxin–responsive SAUR68–like, and vicilin–like antimicrobial peptides 2–2 proteins those were perhaps associated with plant height heterosis. The genes expressed universally due to their overdominace function perhaps influenced MPH for greater plant height― largely by modulating biological processes involved therein. The genes identified in this study might be exploited in heterosis breeding for plant height of L. longiflorum.

The Role of Transcriptional Regulation in Hybrid Vigor

The genetic basis of hybrid vigor in plants remains largely unsolved but strong evidence suggests that variation in transcriptional regulation can explain many aspects of this phenomenon. Natural variation in transcriptional regulation is highly abundant in virtually all species and thus a potential source of heterotic variability. Allele Specific Expression (ASE), which is tightly linked to parent of origin effects and modulated by complex interactions in cis and in trans, is generally considered to play a key role in explaining the differences between hybrids and parental lines. Here we discuss the recent developments in elucidating the role of transcriptional variation in a number of aspects of hybrid vigor, thereby bridging old paradigms and hypotheses with contemporary research in various species.

Heterosis analysis and underlying molecular regulatory mechanism in a wide-compatible neo-tetraploid rice line with long panicles

BACKGROUND

Neo-tetraploid rice, which is a new germplasm developed from autotetraploid rice, has a powerful biological and yield potential and could be used for commercial utilization. The length of panicle, as a part of rice panicle architecture, contributes greatly to high yield. However, little information about long panicle associated with heterosis or hybrid vigor is available in neo-tetraploid rice.

RESULTS

In the present study, we developed a neo-tetraploid rice line, Huaduo 8 (H8), with long panicles and harboring wide-compatibility genes for pollen and embryo sac fertility. All the hybrids generated by H8 produced significant high-parent yield heterosis and displayed long panicles similar to H8. RNA-seq analysis detected a total of 4013, 7050, 6787 and 6195 differentially expressed genes uniquely belonging to F₁ and specifically (DEGFu-sp) associated with leaf, sheath, main panicle axis and spikelet in the two hybrids, respectively. Of these DEGFu-sp, 279 and 89 genes were involved in kinase and synthase, and 714 cloned genes, such as GW8, OsGA20ox1, Ghd8, GW6a, and LP1, were identified and validated by qRT-PCR. A total of 2925 known QTLs intervals, with an average of 1~100 genes per interval, were detected in both hybrids. Of these, 109 yield-related QTLs were associated with seven important traits in rice. Moreover, 1393 non-additive DEGs, including 766 up-regulated and 627 down-regulated, were detected in both hybrids. Importantly, eight up-regulated genes associated with panicle were detected in young panicles of the two hybrids compared to their parents by qRT-PCR. Re-sequencing analysis depicted that LP (a gene controlling long panicle) sequence of H8 was different from many other neo-tetraploid rice and most of the diploid and autotetraploid lines. The qRT-PCR results showed that LP was up-regulated in the hybrid compared to its parents at very young stage of panicle development.

CONCLUSIONS

These results suggested that H8 could overcome the intersubspecific autotetraploid hybrid rice sterility caused by embryo sac and pollen sterility loci. Notably, long panicles of H8 showed dominance phenomenon and played an important role in yield heterosis, which is a complex molecular mechanism. The neo-tetraploid rice is a useful germplasm to attain high yield of polyploid rice.

Comparative transcriptome analysis of inbred lines and contrasting hybrids reveals overdominance mediate early biomass vigor in hybrid cotton

Background

Heterosis breeding is the most useful method for yield increase around the globe. Heterosis is an intriguing process to develop superior offspring to either parent in the desired character. The biomass vigor produced during seedling emergence stage has a direct influence on yield heterosis in plants. Unfortunately, the genetic basis of early biomass vigor in cotton is poorly understood.

Results

Three stable performing F₁ hybrids varying in yield heterosis named as high, medium and low hybrids with their inbred parents were used in this study. Phenotypically, these hybrids established noticeable biomass heterosis during the early stage of seedling growth in the field. Transcriptome analysis of root and leaf revealed that hybrids showed many differentially expressed genes (DEGs) relative to their parents, while the comparison of inbred parents showed limited number of DEGs indicating similarity in their genetic constitution. Further analysis indicated expression patterns of most DEGs were overdominant in both tissues of hybrids. According to GO results, functions of overdominance genes in leaf were enriched for chloroplast, membrane, and protein binding, whereas functions of overdominance genes in root were enriched for plasma membrane, extracellular region, and responses to stress. We found several genes of circadian rhythm pathway related to LATE ELONGATED HYPOCOTYL (LHY) showed downregulated overdominant expressions in both tissues of hybrids. In addition to circadian rhythm, several leaf genes related to Aux/IAA regulation, and many root genes involved in peroxidase activity also showed overdominant expressions in hybrids. Twelve genes involved in circadian rhythm plant were selected to perform qRT-PCR analysis to confirm the accuracy of RNA-seq results.

Conclusions

Through genome-wide comparative transcriptome analysis, we strongly predict that overdominance at gene expression level plays a pivotal role in early biomass vigor of hybrids. The combinational contribution of circadian rhythm and other metabolic process may control vigorous growth in hybrids. Our result provides an important foundation for dissecting molecular mechanisms of biomass vigor in hybrid cotton.

Differential Expression of Genes at Panicle Initiation and Grain Filling Stages Implied in Heterosis of Rice Hybrids

RNA-Seq technology was used to analyze the transcriptome of two rice hybrids, Ajay (based on wild-abortive (WA)-cytoplasm) and Rajalaxmi (based on Kalinga-cytoplasm), and their respective parents at the panicle initiation (PI) and grain filling (GF) stages. Around 293 and 302 million high quality paired-end reads of Ajay and Rajalaxmi, respectively, were generated and aligned against the Nipponbare reference genome. Transcriptome profiling of Ajay revealed 2814 and 4819 differentially expressed genes (DEGs) at the PI and GF stages, respectively, as compared to its parents. In the case of Rajalaxmi, 660 and 5264 DEGs were identified at PI and GF stages, respectively. Functionally relevant DEGs were selected for validation through qRT-PCR, which were found to be co-related with the expression patterns to RNA-seq. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis indicated significant DEGs enriched for energy metabolism pathways, such as photosynthesis, oxidative phosphorylation, and carbon fixation, at the PI stage, while carbohydrate metabolism-related pathways, such as glycolysis and starch and sucrose metabolism, were significantly involved at the GF stage. Many genes involved in energy metabolism exhibited upregulation at the PI stage, whereas the genes involved in carbohydrate biosynthesis had higher expression at the GF stage. The majority of the DEGs were successfully mapped to know yield related rice quantitative trait loci (QTLs). A set of important transcription factors (TFs) was found to be encoded by the identified DEGs. Our results indicated that a complex interplay of several genes in different pathways contributes to higher yield and vigor in rice hybrids.

Genes, pathways and transcription factors involved in seedling stage chilling stress tolerance in indica rice through RNA-Seq analysis

BACKGROUND

Rice plants show yellowing, stunting, withering, reduced tillering and utimately low productivity in susceptible varieties under low temperature stress. Comparative transcriptome analysis was performed to identify novel transcripts, gain new insights into different gene expression and pathways involved in cold tolerance in rice.

RESULTS

Comparative transcriptome analyses of 5 treatments based on chilling stress exposure revealed more down regulated genes in susceptible and higher up regulated genes in tolerant genotypes. A total of 13930 and 10599 differentially expressed genes (DEGs) were detected in cold susceptible variety (CSV) and cold tolerant variety (CTV), respectively. A continuous increase in DEGs at 6, 12, 24 and 48 h exposure of cold stress was detected in both the genotypes. Gene ontology (GO) analysis revealed 18 CSV and 28 CTV term significantly involved in molecular function, cellular component and biological process. GO classification showed a significant role of transcription regulation, oxygen, lipid binding, catalytic and hydrolase activity for tolerance response. Absence of photosynthesis related genes, storage products like starch and synthesis of other classes of molecules like fatty acids and terpenes during the stress were noticed in susceptible genotype. However, biological regulations, generation of precursor metabolites, signal transduction, photosynthesis, regulation of cellular process, energy and carbohydrate metabolism were seen in tolerant genotype during the stress. KEGG pathway annotation revealed more number of genes regulating different pathways resulting in more tolerant. During early response phase, 24 and 11 DEGs were enriched in CTV and CSV, respectively in energy metabolism pathways. Among the 1583 DEG transcription factors (TF) genes, 69 WRKY, 46 bZIP, 41 NAC, 40 ERF, 31/14 MYB/MYB-related, 22 bHLH, 17 Nin-like 7 HSF and 4C3H were involved during early response phase. Late response phase showed 30 bHLH, 65 NAC, 30 ERF, 26/20 MYB/MYB-related, 11 C3H, 12 HSF, 86 Nin-like, 41 AP2/ERF, 55 bZIP and 98 WRKY members TF genes. The recovery phase included 18 bHLH, 50 NAC, 31 ERF, 24/13 MYB/MYB-related, 4 C3H, 4 HSF, 14 Nin-like, 31 bZIP and 114 WRKY TF genes.

CONCLUSIONS

Transcriptome analysis of contrasting genotypes for cold tolerance detected the genes, pathways and transcription factors involved in the stress tolerance.

Genes, pathways and transcription factors involved in seedling stage chilling stress tolerance in indica rice through RNA-Seq analysis

BACKGROUND

Rice plants show yellowing, stunting, withering, reduced tillering and utimately low productivity in susceptible varieties under low temperature stress. Comparative transcriptome analysis was performed to identify novel transcripts, gain new insights into different gene expression and pathways involved in cold tolerance in rice.

RESULTS

Comparative transcriptome analyses of 5 treatments based on chilling stress exposure revealed more down regulated genes in susceptible and higher up regulated genes in tolerant genotypes. A total of 13930 and 10599 differentially expressed genes (DEGs) were detected in cold susceptible variety (CSV) and cold tolerant variety (CTV), respectively. A continuous increase in DEGs at 6, 12, 24 and 48 h exposure of cold stress was detected in both the genotypes. Gene ontology (GO) analysis revealed 18 CSV and 28 CTV term significantly involved in molecular function, cellular component and biological process. GO classification showed a significant role of transcription regulation, oxygen, lipid binding, catalytic and hydrolase activity for tolerance response. Absence of photosynthesis related genes, storage products like starch and synthesis of other classes of molecules like fatty acids and terpenes during the stress were noticed in susceptible genotype. However, biological regulations, generation of precursor metabolites, signal transduction, photosynthesis, regulation of cellular process, energy and carbohydrate metabolism were seen in tolerant genotype during the stress. KEGG pathway annotation revealed more number of genes regulating different pathways resulting in more tolerant. During early response phase, 24 and 11 DEGs were enriched in CTV and CSV, respectively in energy metabolism pathways. Among the 1583 DEG transcription factors (TF) genes, 69 WRKY, 46 bZIP, 41 NAC, 40 ERF, 31/14 MYB/MYB-related, 22 bHLH, 17 Nin-like 7 HSF and 4C3H were involved during early response phase. Late response phase showed 30 bHLH, 65 NAC, 30 ERF, 26/20 MYB/MYB-related, 11 C3H, 12 HSF, 86 Nin-like, 41 AP2/ERF, 55 bZIP and 98 WRKY members TF genes. The recovery phase included 18 bHLH, 50 NAC, 31 ERF, 24/13 MYB/MYB-related, 4 C3H, 4 HSF, 14 Nin-like, 31 bZIP and 114 WRKY TF genes.

CONCLUSIONS

Transcriptome analysis of contrasting genotypes for cold tolerance detected the genes, pathways and transcription factors involved in the stress tolerance.

Comparative transcriptome analysis reveals important roles of nonadditive genes in maize hybrid An'nong 591 under heat stress

BACKGROUND

Heterosis is the superior performance of F₁ hybrids relative to their parental lines for a wide range of traits. In this study, expression profiling and heterosis associated genes were analyzed by RNA sequencing (RNA-Seq) in seedlings of the maize hybrid An'nong 591 and its parental lines under control and heat stress conditions.

RESULTS

Through performing nine pairwise comparisons, the maximum number of differentially expressed genes (DEGs) was detected between the two parental lines, and the minimum number was identified between the F₁ hybrid and the paternal lines under both conditions, which suggested greater genetic contribution of the paternal line to heat stress tolerance. Gene Ontology (GO) enrichment analysis of the 4518 common DEGs indicated that GO terms associated with diverse stress responses and photosynthesis were highly overrepresented in the 76 significant terms of the biological process category. A total of 3970 and 7653 genes exhibited nonadditive expression under control and heat stress, respectively. Among these genes, 2253 (56.8%) genes overlapped, suggesting that nonadditive genes tend to be conserved in expression. In addition, 5400 nonadditive genes were found to be specific for heat stress condition, and further GO analysis indicated that terms associated with stress responses were significantly overrepresented, and 60 genes were assigned to the GO term response to heat. Pathway enrichment analysis indicated that 113 genes were involved in spliceosome metabolic pathways. Nineteen of the 33 overlapping genes assigned to the GO term response to heat showed significantly higher number of alternative splicing (AS) events under heat stress than under control conditions, suggesting that AS of these genes play an important role in response to heat stress.

CONCLUSIONS

This study reveals the transcriptomic divergence of the maize F₁ hybrid and its parental lines under control and heat stress conditions, and provides insight into the underlying molecular mechanisms of heterosis and the response to heat stress in maize.

Comparative transcriptome analysis reveals unique genetic adaptations conferring salt tolerance in a xerohalophyte

Most studies on salt tolerance in plants have been conducted using glycophytes like Arabidopsis thaliana (L.) Heynh., with limited resistance to salinity. The xerohalophyte Zygophyllum xanthoxylum (Bunge) Engl. is a salt-accumulating desert plant that efficiently transports Na+ into vacuoles to manage salt and exhibits increased growth under salinity conditions, suggesting a unique transcriptional response compared with glycophytes. We used transcriptome profiling by RNA-seq to compare gene expression in roots of Z. xanthoxylum and A. thaliana under 50 mM NaCl treatments. Gene Ontology (GO) functional annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) metabolic pathway analysis suggested that 50 mM NaCl was perceived as a stimulus for Z. xanthoxylum whereas a stress for A. thaliana. Exposure to 50 mM NaCl caused metabolic shifts towards gluconeogenesis to stimulate growth of Z. xanthoxylum, but triggered defensive systems in A. thaliana. Compared with A. thaliana, a vast array of ion transporter genes was induced in Z. xanthoxylum, revealing an active strategy to uptake Na+ and nutrients from the environment. An ascorbate-glutathione scavenging system for reactive oxygen species was also crucial in Z. xanthoxylum, based on high expression of key enzyme genes. Finally, key regulatory genes for the biosynthesis pathways of abscisic acid and gibberellin showed distinct expression patterns between the two species and auxin response genes were more active in Z. xanthoxylum compared with A. thaliana. Our results provide an important framework for understanding unique patterns of gene expression conferring salt resistance in Z. xanthoxylum.

Carbohydrate metabolism and fertility related genes high expression levels promote heterosis in autotetraploid rice harboring double neutral genes

BACKGROUND

Autotetraploid rice hybrids have great potential to increase the production, but hybrid sterility is a major hindrance in the utilization of hybrid vigor in polyploid rice, which is mainly caused by pollen abortion. Our previous study showed that double pollen fertility neutral genes, Sa-n and Sb-n, can overcome hybrid sterility in autotetraploid rice. Here, we used an autotetraploid rice line harboring double neutral genes to develop hybrids by crossing with auto- and neo-tetraploid rice, and evaluated heterosis and its underlying molecular mechanism.

RESULTS

All autotetraploid rice hybrids, which harbored double pollen fertility neutral genes, Sa-n and Sb-n, displayed high seed setting and significant positive heterosis for yield and yield-related traits. Cytological observations revealed normal chromosome behaviors and higher frequency of bivalents in the hybrid than parents during meiosis. Transcriptome analysis revealed significantly higher expressions of important saccharides metabolism and starch synthase related genes, such as OsBEIIb and OsSSIIIa, in the grains of hybrid than parents. Furthermore, many meiosis-related and specific genes, including DPW and CYP703A3, displayed up-regulation in the hybrid compared to a parent with low seed setting. Many non-additive genes were detected in the hybrid, and GO term of carbohydrate metabolic process was significantly enriched in all the transcriptome tissues except flag leaf (three days after flowering). Moreover, many differentially expressed genes (DEGs) were identified in the yield-related quantitative trait loci (QTLs) regions as possible candidate genes.

CONCLUSION

Our results revealed that increase in the number of bivalents improved the seed setting of hybrid harboring double pollen fertility neutral genes. Many important genes, including meiosis-related and meiosis-specific genes and saccharides metabolism and starch synthase related genes, exhibited heterosis specific expression patterns in polyploid rice during different development stages. The functional analysis of important genes will provide valuable information for molecular mechanisms of heterosis in polyploid rice.

Transcript analyses reveal a comprehensive role of abscisic acid in modulating fruit ripening in Chinese jujube

BACKGROUND

Chinese jujube (Ziziphus jujuba Mill.) is a non-climacteric fruit; however, the underlying mechanism of ripening and the role of abscisic acid involved in this process are not yet understood for this species.

RESULTS

In the present study, a positive correlation between dynamic changes in endogenous ABA and the onset of jujube ripening was determined. Transcript analyses suggested that the expression balance among genes encoding nine-cis-epoxycarotenoid dioxygenase (ZjNCED3), ABA-8'-hydroxylase (ZjCYP707A2), and beta-glucosidase (ZjBG4, ZjBG5, ZjBG8, and ZjBG9) has an important role in maintaining ABA accumulation, while the expression of a receptor (ZjPYL8), protein phosphatase 2C (ZjPP2C4-8), and sucrose nonfermenting 1-related protein kinase 2 (ZjSnRK2-2 and ZjSnRK2-5) is important in regulating fruit sensitivity to ABA applications. In addition, white mature 'Dongzao' fruit were harvested and treated with 50 mg L- 1 ABA or 50 mg L- 1 nordihydroguaiaretic acid (NDGA) to explore the role of ABA in jujube fruit ripening. By comparative transcriptome analyses, 1103 and 505 genes were differentially expressed in response to ABA and NDGA applications on the 1st day after treatment, respectively. These DEGs were associated with photosynthesis, secondary, lipid, cell wall, and starch and sugar metabolic processes, suggesting the involvement of ABA in modulating jujube fruit ripening. Moreover, ABA also exhibited crosstalk with other phytohormones and transcription factors, indicating a regulatory network for jujube fruit ripening.

CONCLUSIONS

Our study further elucidated ABA-associated metabolic and regulatory processes. These findings are helpful for improving strategies for jujube fruit storage and for gaining insights into understand complex non-climacteric fruit ripening processes.

Transcriptome analysis reveals that barnyard grass exudates increase the allelopathic potential of allelopathic and non-allelopathic rice (Oryza sativa) accessions

BACKGROUND

Allelopathy in rice (Oryza sativa) is a chemically induced response that is elevated by the exogenous application of chemical compounds and barnyard grass root exudates. An in-depth understanding of the response mechanisms of rice to chemical induction is necessary for the identification of target genes for increasing the allelopathic potential of rice. However, no previous studies have evaluated the transcriptomic changes associated with allelopathy in rice in response to barnyard grass exudates treatment. Thus, the aim of the present study was to reveal differentially expressed genes (DEGs) in allelopathic and non-allelopathic rice seedlings treated with barnyard grass exudates to identify target allelopathy genes.

RESULTS

The inhibitory effect of the culture solutions on the allelopathic rice accession PI312777 (PI) and the non-allelopathic rice accession Lemont (LE) significantly increased (P < 0.05) after treatment with barnyard grass root exudates. The RNA sequencing results revealed that 14,891 genes in PI(+B) vs. LE(+B), 12,505 genes in PI(+B) vs. PI(-B), and 5857 genes in LE(+B) vs. LE(-B) were differentially expressed following root exudates treatment. These DEGs were classified into three categories and 32 functional groups, i.e., 12 groups in the biological process category, 12 groups in the cellular component category, and eight groups in the molecular function category. There were 5857 and 2846 upregulated genes and 135 and 50 upregulated Gene Ontology terms (P < 0.05) in the biological process category in PI(+B) vs. PI(-B) and LE(+B) vs. LE(-B), respectively. These results indicated that the allelopathic accession PI is more sensitive than the non-allelopathic accession LE to exogenous root exudates treatment. Genes related to rice allelochemical-related biosynthesis pathways, particularly the shikimic acid and acetic acid pathways, were significantly differentially expressed in both rice accessions. These findings suggested that phenolic acids, fatty acids, and flavonoids, which constitute the downstream metabolites of the shikimic acid and acetic acid pathways, are significantly expressed in response to root exudates of barnyard grass.

CONCLUSIONS

The allelopathic potential of both rice accessions could be significantly enhanced by barnyard grass root exudates application. Furthermore, genes related to the biosynthesis pathways of reported rice allelochemicals were significantly differentially expressed in both accessions. Phenylalanine ammonia lyase was determined to be a potential target for the regulation of chemical induction.

Heterotic patterns of primary and secondary metabolites in the oilseed crop Brassica juncea

Heterosis refers to the superior performance of F1 hybrids over their respective parental inbred lines. Although the genetic and expression basis of heterosis have been previously investigated, the metabolic basis for this phenomenon is poorly understood. In a preliminary morphological study in Brassica juncea, we observed significant heterosis at the 50% flowering stage, wherein both the growth and reproduction of F1 reciprocal hybrids were greater than that of their parents. To identify the possible metabolic causes or consequences of this heterosis, we carried out targeted LC-MS analysis of 48 primary (amino acids and sugars) and secondary metabolites (phytohormones, glucosinolates, flavonoids, and phenolic esters) in five developmental tissues at 50% flowering in hybrids and inbred parents. Principal component analysis (PCA) of metabolites clearly separated inbred lines from their hybrids, particularly in the bud tissues. In general, secondary metabolites displayed more negative heterosis values in comparison to primary metabolites. The tested primary and secondary metabolites displayed both additive and non-additive modes of inheritance in F1 hybrids, wherein the number of metabolites showing an additive mode of inheritance were higher in buds and siliques (52.77-97.14%) compared to leaf tissues (47.37-80%). Partial least regression (PLS) analysis further showed that primary metabolites, in general, displayed higher association with morphological parameters in F1 hybrids. Overall, our results are consistent with a resource-cost model for heterosis in B. juncea, where metabolite allocation in hybrids appears to favor growth, at the expense of secondary metabolism.

Transcriptomic analysis of grain filling in rice inferior grains under moderate soil drying

Moderate soil drying imposed at the post-anthesis stage significantly increases starch accumulation in inferior grains of rice, but how this process is regulated at the level of gene expression remains unclear. In this study, we applied moderate drying (MD) treatments to the soil at the post-anthesis stage and followed the dynamics of the conversion process of soluble sugars to starch in inferior grains using RNA-seq analysis. An elevated level of ABA induced by MD was consistently associated with down-regulation of ABA8ox2, suggesting that lower expression of this gene may be responsible for the higher ABA content, potentially resulting in better filling in inferior grains. In addition, MD treatments up-regulated genes encoding five key enzymes involved sucrose-to-starch conversion and increased the activities of enzymes responsible for soluble-sugar reduction and starch accumulation in inferior grains. Differentially expressed transcription factors, including NAC, GATA, WRKY, and M-type MADS, were predicted to interact with other proteins in mediating filling of inferior grains as a response to MD. Transient expression analysis showed that NAC activated WAXY expression by binding to its promoter, indicating that NAC played a key role in starch synthesis of inferior grains under MD treatment. Our results provide new insights into the molecular mechanisms that regulate grain filling in inferior grains of rice under moderate soil drying.

Expression of Maize MADS Transcription Factor ZmES22 Negatively Modulates Starch Accumulation in Rice Endosperm

As major component in cereals grains, starch has been one of the most important carbohydrate consumed by a majority of world’s population. However, the molecular mechanism for regulation of biosynthesis of starch remains elusive. In the present study, ZmES22, encoding a MADS-type transcription factor, was modestly characterized from maize inbred line B73. ZmES22 exhibited high expression level in endosperm at 10 days after pollination (DAP) and peaked in endosperm at 20 DAP, indicating that ZmES22 was preferentially expressed in maize endosperm during active starch synthesis. Transient expression of ZmES22 in tobacco leaf revealed that ZmES22 protein located in nucleus. No transactivation activity could be detected for ZmES22 protein via yeast one-hybrid assay. Transformation of overexpressing plasmid 35S::ZmES22 into rice remarkedly reduced 1000-grain weight as well as the total starch content, while the soluble sugar was significantly higher in transgenic rice lines. Moreover, overexpressing ZmES22 reduced fractions of long branched starch. Scanning electron microscopy images of transverse sections of rice grains revealed that altered expression of ZmES22 also changed the morphology of starch granule from densely packed, polyhedral starch granules into loosely packed, spherical granules with larger spaces. Furthermore, RNA-seq results indicated that overexpressing ZmES22 could significantly influence mRNA expression levels of numerous key regulatory genes in starch synthesis pathway. Y1H assay illustrated that ZmES22 protein could bind to the promoter region of OsGIF1 and downregulate its mRNA expression during rice grain filling stages. These findings suggest that ZmES22 was a novel regulator during starch synthesis process in rice endosperm.

RNA-Seq Analyses Identify Additivity as the Predominant Gene Expression Pattern in F1 Chicken Embryonic Brain and Liver

The superior performance of hybrids to parents, termed heterosis, has been widely utilized in animal and plant breeding programs, but the molecular mechanism underlying heterosis remains an enigma. RNA-Seq provides a novel way to investigate heterosis at the transcriptome-wide level, because gene expression functions as an intermediate phenotype that contributes to observable traits. Here we compared embryonic gene expression between chicken hybrids and their inbred parental lines to identify inheritance patterns of gene expression. Inbred Fayoumi and Leghorn were crossed reciprocally to obtain F1 fertile eggs. RNA-Seq was carried out using 24 brain and liver samples taken from day 12 embryos, and the differentially expressed (DE) genes were identified by pairwise comparison among the hybrids, parental lines, and mid-parent expression values. Our results indicated the expression levels of the majority of the genes in the F1 cross are not significantly different from the mid-parental values, suggesting additivity as the predominant gene expression pattern in the F1. The second and third prevalent gene expression patterns are dominance and over-dominance. Additionally, we found only 7⁻20% of the DE genes exhibit allele-specific expression in the F1, suggesting that trans regulation is the main driver for differential gene expression and thus contributes to heterosis effect in the F1 crosses.

Differential transcriptome patterns associated with early seedling development in a wild and a domesticated common bean (Phaseolus vulgaris L.) accession

Genes that control "Domestication Syndrome" traits were direct targets of selection, like those controlling increased seed size in the common bean. However, selection for this trait brought about unintentional selection on genes controlling seedling growth. We hypothesized that wild and domesticated plants have different early seedling growth patterns as an indirect consequence of selection for a larger seed size during domestication, and those differences resulted from changes in gene expression patterns of the wild ancestor. Large seeds pose a challenge to reserve remobilization during early heterotrophic growth, particularly during a transition towards more fertile alluvial soils. To address our hypothesis, we characterized the patterns of gene expression of cotyledon, root, and leaf tissues of 7-day old seedlings of a wild and a landrace accession of the common bean. Differential expression analyses detected genes with contrasting patterns of expression between the two genotypes in all three tissues. Some of the differentially expressed genes with contrasting genotypic patterns are known to have domestication-related signatures of selection. Among these genes were some transcription factors associated with key roles in development. These genes may represent targets of indirect selection and ultimately explain the growth phenotypic differences between wild and domesticated seedlings.