Temporal Regulation of the Metabolome and Proteome in Photosynthetic and Photorespiratory Pathways Contributes to Maize Heterosis

Impact of glioma metabolism-related gene ALPK1 on tumor immune heterogeneity and the regulation of the TGF-β pathway

Background

Recent years have seen persistently poor prognoses for glioma patients. Therefore, exploring the molecular subtyping of gliomas, identifying novel prognostic biomarkers, and understanding the characteristics of their immune microenvironments are crucial for improving treatment strategies and patient outcomes.

Methods

We integrated glioma datasets from multiple sources, employing Non-negative Matrix Factorization (NMF) to cluster samples and filter for differentially expressed metabolic genes. Additionally, we utilized Weighted Gene Co-expression Network Analysis (WGCNA) to identify key genes. A predictive model was developed utilizing the optimal consistency index derived from a combination of 101 machine learning techniques, and its effectiveness was confirmed through multiple datasets employing different methodologies. In-depth analyses were conducted on immune cell infiltration and tumor microenvironmental aspects. Single-cell sequencing data were employed for clustering and differential expression analysis of genes associated with glioma. Finally, the immune relevance of the model gene ALPK1 in the context of pan-cancer was explored, including its relationship with immune checkpoints.

Results

The application of NMF, coupled with differential analysis of metabolic-related genes, led to the identification of two clusters exhibiting significant differences in survival, age, and metabolic gene expression among patients. Core genes were identified through WGCNA, and a total of 101 machine learning models were constructed, with LASSO+GBM selected as the optimal model, demonstrating robust validation performance. Comprehensive analyses revealed that high-risk groups exhibited greater expression of specific genes, with ALPK1 showing significant correlations with immune regulation.

Conclusion

This research employed a multi-dataset strategy and various methods to clarify the differences in metabolic traits and immune conditions in glioma patients, while creating an innovative prognostic risk evaluation framework. These results offer fresh perspectives on the intricate biological processes that define gliomas.

Molecular concepts to explain heterosis in crops

Heterosis describes the superior performance of hybrid plants compared with their genetically distinct parents and is a pillar of global food security. Here we review the current status of the molecular dissection of heterosis. We discuss how extensive intraspecific structural genomic variation between parental genotypes leads to heterosis by genetic complementation in hybrids. Moreover, we survey how global gene expression complementation contributes to heterosis by hundreds of additionally active genes in hybrids and how overdominant single genes mediate heterosis in several species. Furthermore, we highlight the prominent role of the microbiome in improving the performance of hybrids. Taken together, the molecular understanding of heterosis will pave the way to accelerate hybrid productivity and a more sustainable agriculture.

Hybrid Prediction in Horticulture Crop Breeding: Progress and Challenges

In the context of rapidly increasing population and diversified market demands, the steady improvement of yield and quality in horticultural crops has become an urgent challenge that modern breeding efforts must tackle. Heterosis, a pivotal theoretical foundation for plant breeding, facilitates the creation of superior hybrids through crossbreeding and selection among a variety of parents. However, the vast number of potential hybrids presents a significant challenge for breeders in efficiently predicting and selecting the most promising candidates. The development and refinement of effective hybrid prediction methods have long been central to research in this field. This article systematically reviews the advancements in hybrid prediction for horticultural crops, including the roles of marker-assisted breeding and genomic prediction in phenotypic forecasting. It also underscores the limitations of some predictors, like genetic distance, which do not consistently offer reliable hybrid predictions. Looking ahead, it explores the integration of phenomics with genomic prediction technologies as a means to elevate prediction accuracy within actual breeding programs.

Synergy of gut microbiota and host genome in driving heterosis expression of chickens

Heterosis has been widely utilized in agricultural production. Despite over a century of extensive research, the underlying mechanisms of heterosis remain elusive. Most hypotheses and research have focused on the genetic basis of heterosis. However, the potential role of gut microbiota in heterosis has been largely ignored. Here, we carefully design a crossbreeding experiment with two distinct broiler breeds and conduct 16S rRNA amplicon and transcriptome sequencing to investigate the synergistic role of gut microbiota and host genes in driving heterosis. We find that the breast muscle weight of hybrids exhibits a high heterosis, 6.28% higher than the mid-parent value. A notable difference is observed in the composition and potential function of cecal microbiota between hybrids and their parents. Over 90% of differentially colonized microbiota and differentially expressed genes exhibit nonadditive patterns. Integrative analyses uncover associations between nonadditive genes and nonadditive microbiota, including a connection between the expression of cellular signaling pathways and metabolism-related genes and the abundance of Odoribacter, Oscillibacter, and Alistipes in hybrids. Moreover, higher abundances of these microbiota are related to better meat yield. In summary, these findings highlight the importance of gut microbiota in heterosis, serving as crucial factors that modulate heterosis expression in chickens.

Rhythmic lipid and gene expression responses to chilling in panicoid grasses

Chilling stress threatens plant growth and development, particularly affecting membrane fluidity and cellular integrity. Understanding plant membrane responses to chilling stress is important for unraveling the molecular mechanisms of stress tolerance. Whereas core transcriptional responses to chilling stress and stress tolerance are conserved across species, the associated changes in membrane lipids appear to be less conserved, as which lipids are affected by chilling stress varies by species. Here, we investigated changes in gene expression and membrane lipids in response to chilling stress during one 24 h cycle in chilling-tolerant foxtail millet (Setaria italica), and chilling-sensitive sorghum (Sorghum bicolor) and Urochloa (browntop signal grass, Urochloa fusca, lipids only), leveraging their evolutionary relatedness and differing levels of chilling stress tolerance. We show that most chilling-induced lipid changes are conserved across the three species, while we observed distinct, time-specific responses in chilling-tolerant foxtail millet, indicating the presence of a finely orchestrated adaptive mechanism. We detected rhythmicity in lipid responses to chilling stress in the three grasses, which were also present in Arabidopsis thaliana, suggesting the conservation of rhythmic patterns across species and highlighting the importance of accounting for time of day. When integrating lipid datasets with gene expression profiles, we identified potential candidate genes that showed corresponding transcriptional changes in response to chilling stress, providing insights into the differences in regulatory mechanisms between chilling-sensitive sorghum and chilling-tolerant foxtail millet.

Circadian and photoperiodic regulation of the vegetative to reproductive transition in plants

As sessile organisms, plants must respond constantly to ever-changing environments to complete their life cycle; this includes the transition from vegetative growth to reproductive development. This process is mediated by photoperiodic response to sensing the length of night or day through circadian regulation of light-signaling molecules, such as phytochromes, to measure the length of night to initiate flowering. Flowering time is the most important trait to optimize crop performance in adaptive regions. In this review, we focus on interplays between circadian and light signaling pathways that allow plants to optimize timing for flowering and seed production in Arabidopsis, rice, soybean, and cotton. Many crops are polyploids and domesticated under natural selection and breeding. In response to adaptation and polyploidization, circadian and flowering pathway genes are epigenetically reprogrammed. Understanding the genetic and epigenetic bases for photoperiodic flowering will help improve crop yield and resilience in response to climate change.

Unlocking allelic variation in circadian clock genes to develop environmentally robust and productive crops

MAIN CONCLUSION

Molecular mechanisms of biological rhythms provide opportunities to harness functional allelic diversity in core (and trait- or stress-responsive) oscillator networks to develop more climate-resilient and productive germplasm. The circadian clock senses light and temperature in day-night cycles to drive biological rhythms. The clock integrates endogenous signals and exogenous stimuli to coordinate diverse physiological processes. Advances in high-throughput non-invasive assays, use of forward- and inverse-genetic approaches, and powerful algorithms are allowing quantitation of variation and detection of genes associated with circadian dynamics. Circadian rhythms and phytohormone pathways in response to endogenous and exogenous cues have been well documented the model plant Arabidopsis. Novel allelic variation associated with circadian rhythms facilitates adaptation and range expansion, and may provide additional opportunity to tailor climate-resilient crops. The circadian phase and period can determine adaptation to environments, while the robustness in the circadian amplitude can enhance resilience to environmental changes. Circadian rhythms in plants are tightly controlled by multiple and interlocked transcriptional-translational feedback loops involving morning (CCA1, LHY), mid-day (PRR9, PRR7, PRR5), and evening (TOC1, ELF3, ELF4, LUX) genes that maintain the plant circadian clock ticking. Significant progress has been made to unravel the functions of circadian rhythms and clock genes that regulate traits, via interaction with phytohormones and trait-responsive genes, in diverse crops. Altered circadian rhythms and clock genes may contribute to hybrid vigor as shown in Arabidopsis, maize, and rice. Modifying circadian rhythms via transgenesis or genome-editing may provide additional opportunities to develop crops with better buffering capacity to environmental stresses. Models that involve clock gene‒phytohormone‒trait interactions can provide novel insights to orchestrate circadian rhythms and modulate clock genes to facilitate breeding of all season crops.

Comparative transcriptomic analysis reveals the molecular mechanism underlying seedling heterosis and its relationship with hybrid contemporary seeds DNA methylation in soybean

Heterosis is widely used in crop production, but phenotypic dominance and its underlying causes in soybeans, a significant grain and oil crop, remain a crucial yet unexplored issue. Here, the phenotypes and transcriptome profiles of three inbred lines and their resulting F1 seedlings were analyzed. The results suggest that F1 seedlings with superior heterosis in leaf size and biomass exhibited a more extensive recompilation in their transcriptional network and activated a greater number of genes compared to the parental lines. Furthermore, the transcriptional reprogramming observed in the four hybrid combinations was primarily non-additive, with dominant effects being more prevalent. Enrichment analysis of sets of differentially expressed genes, coupled with a weighted gene co-expression network analysis, has shown that the emergence of heterosis in seedlings can be attributed to genes related to circadian rhythms, photosynthesis, and starch synthesis. In addition, we combined DNA methylation data from previous immature seeds and observed similar recompilation patterns between DNA methylation and gene expression. We also found significant correlations between methylation levels of gene region and gene expression levels, as well as the discovery of 12 hub genes that shared or conflicted with their remodeling patterns. This suggests that DNA methylation in contemporary hybrid seeds have an impact on both the F1 seedling phenotype and gene expression to some extent. In conclusion, our study provides valuable insights into the molecular mechanisms of heterosis in soybean seedlings and its practical implications for selecting superior soybean varieties.

Toward understanding and utilizing crop heterosis in the age of biotechnology

Heterosis, a universal phenomenon in nature, mainly reflected in the superior productivity, quality, and fitness of F1 hybrids compared with their inbred parents, has been exploited in agriculture and greatly benefited human society in terms of food security. However, the flexible and efficient utilization of heterosis has remained a challenge in hybrid breeding systems because of the limitations of "three-line" and "two-line" methods. In the past two decades, rapidly developed biotechnologies have provided unprecedented conveniences for both understanding and utilizing heterosis. Notably, "third-generation" (3G) hybrid breeding technology together with high-throughput sequencing and gene editing greatly promoted the efficiency of hybrid breeding. Here, we review emerging ideas about the genetic or molecular mechanisms of heterosis and the development of 3G hybrid breeding system in the age of biotechnology. In addition, we summarized opportunities and challenges for optimal heterosis utilization in the future.

Omics-driven exploration and mining of key functional genes for the improvement of food and fiber crops

The deployment of omics technologies has obtained an incredible boost over the past few decades with the advances in next-generation sequencing (NGS) technologies, innovative bioinformatics tools, and the deluge of available biological information. The major omics technologies in the limelight are genomics, transcriptomics, proteomics, metabolomics, and phenomics. These biotechnological advances have modernized crop breeding and opened new horizons for developing crop varieties with improved traits. The genomes of several crop species are sequenced, and a huge number of genes associated with crucial economic traits have been identified. These identified genes not only provide insights into the understanding of regulatory mechanisms of crop traits but also decipher practical grounds to assist in the molecular breeding of crops. This review discusses the potential of omics technologies for the acquisition of biological information and mining of the genes associated with important agronomic traits in important food and fiber crops, such as wheat, rice, maize, potato, tomato, cassava, and cotton. Different functional genomics approaches for the validation of these important genes are also highlighted. Furthermore, a list of genes discovered by employing omics approaches is being represented as potential targets for genetic modifications by the latest genome engineering methods for the development of climate-resilient crops that would in turn provide great impetus to secure global food security.

RefMetaPlant: a reference metabolome database for plants across five major phyla

Plants are unique with tremendous chemical diversity and metabolic complexity, which is highlighted by estimates that green plants collectively produce metabolites numbering in the millions. Plant metabolites play crucial roles in all aspects of plant biology, like growth, development, stress responses, etc. However, the lack of a reference metabolome for plants, and paucity of high-quality standard compound spectral libraries and related analytical tools, have hindered the discovery and functional study of phytochemicals in plants. Here, by leveraging an advanced LC-MS platform, we generated untargeted mass spectral data from >150 plant species collected across the five major phyla. Using a self-developed computation protocol, we constructed reference metabolome for 153 plant species. A 'Reference Metabolome Database for Plants' (RefMetaPlant) was built to encompass the reference metabolome, integrated standard compound mass spectral libraries for annotation, and related query and analytical tools like 'LC-MS/MS Query', 'RefMetaBlast' and 'CompoundLibBlast' for searches and profiling of plant metabolome and metabolite identification. Analogous to a reference genome in genomic research, RefMetaPlant provides a powerful platform to support plant genome-scale metabolite analysis to promote knowledge/data sharing and collaboration in the field of metabolomics. RefMetaPlant is freely available at https://www.biosino.org/RefMetaDB/.

Phenotypic and Proteomic Insights into Differential Cadmium Accumulation in Maize Kernels

The contamination of agricultural soil with cadmium (Cd), a heavy metal, poses a significant environmental challenge, affecting crop growth, development, and human health. Previous studies have established the pivotal role of the ZmHMA3 gene, a P-type ATPase heavy metal transporter, in determining variable Cd accumulation in maize grains among 513 inbred lines. To decipher the molecular mechanism underlying mutation-induced phenotypic differences mediated by ZmHMA3, we conducted a quantitative tandem mass tag (TMT)-based proteomic analysis of immature maize kernels. This analysis aimed to identify differentially expressed proteins (DEPs) in wild-type B73 and ZmHMA3 null mutant under Cd stress. The findings demonstrated that ZmHMA3 accumulated higher levels of Cd compared to B73 when exposed to varying Cd concentrations in the soil. In comparison to soil with a low Cd concentration, B73 and ZmHMA3 exhibited 75 and 142 DEPs, respectively, with 24 common DEPs shared between them. ZmHMA3 showed a higher induction of upregulated genes related to Cd stress than B73. Amino sugar and nucleotide sugar metabolism was specifically enriched in B73, while phenylpropanoid biosynthesis, nitrogen metabolism, and glyoxylate and dicarboxylate metabolism appeared to play a more significant role in ZmHMA3. This study provides proteomics insights into unraveling the molecular mechanism underlying the differences in Cd accumulation in maize kernels.

Protein nonadditive expression and solubility contribute to heterosis in Arabidopsis hybrids and allotetraploids

Hybrid vigor or heterosis has been widely applied in agriculture and extensively studied using genetic and gene expression approaches. However, the biochemical mechanism underlying heterosis remains elusive. One theory suggests that a decrease in protein aggregation may occur in hybrids due to the presence of protein variants between parental alleles, but it has not been experimentally tested. Here, we report comparative analysis of soluble and insoluble proteomes in Arabidopsis intraspecific and interspecific hybrids or allotetraploids formed between A. thaliana and A. arenosa. Both allotetraploids and intraspecific hybrids displayed nonadditive expression (unequal to the sum of the two parents) of the proteins, most of which were involved in biotic and abiotic stress responses. In the allotetraploids, homoeolog-expression bias was not observed among all proteins examined but accounted for 17-20% of the nonadditively expressed proteins, consistent with the transcriptome results. Among expression-biased homoeologs, there were more A. thaliana-biased than A. arenosa-biased homoeologs. Analysis of the insoluble and soluble proteomes revealed more soluble proteins in the hybrids than their parents but not in the allotetraploids. Most proteins in ribosomal biosynthesis and in the thylakoid lumen, membrane, and stroma were in the soluble fractions, indicating a role of protein stability in photosynthetic activities for promoting growth. Thus, nonadditive expression of stress-responsive proteins and increased solubility of photosynthetic proteins may contribute to heterosis in Arabidopsis hybrids and allotetraploids and possibly hybrid crops.

High-quality Cymbidium mannii genome and multifaceted regulation of crassulacean acid metabolism in epiphytes

Epiphytes with crassulacean acid metabolism (CAM) photosynthesis are widespread among vascular plants, and repeated evolution of CAM photosynthesis is a key innovation for micro-ecosystem adaptation. However, we lack a complete understanding of the molecular regulation of CAM photosynthesis in epiphytes. Here, we report a high-quality chromosome-level genome assembly of a CAM epiphyte, Cymbidium mannii (Orchidaceae). The 2.88-Gb orchid genome with a contig N50 of 22.7 Mb and 27 192 annotated genes was organized into 20 pseudochromosomes, 82.8% of which consisted of repetitive elements. Recent expansions of long terminal repeat retrotransposon families have made a major contribution to the evolution of genome size in Cymbidium orchids. We reveal a holistic scenario of molecular regulation of metabolic physiology using high-resolution transcriptomics, proteomics, and metabolomics data collected across a CAM diel cycle. Patterns of rhythmically oscillating metabolites, especially CAM-related products, reveal circadian rhythmicity in metabolite accumulation in epiphytes. Genome-wide analysis of transcript and protein level regulation revealed phase shifts during the multifaceted regulation of circadian metabolism. Notably, we observed diurnal expression of several core CAM genes (especially βCA and PPC) that may be involved in temporal fixation of carbon sources. Our study provides a valuable resource for investigating post-transcription and translation scenarios in C. mannii, an Orchidaceae model for understanding the evolution of innovative traits in epiphytes.

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.

Disentangling the Heterosis in Biomass Production and Radiation Use Efficiency in Maize: A Phytomer-Based 3D Modelling Approach

Maize (Zea mays L.) benefits from heterosis in-yield formation and photosynthetic efficiency through optimizing canopy structure and improving leaf photosynthesis. However, the role of canopy structure and photosynthetic capacity in determining heterosis in biomass production and radiation use efficiency has not been separately clarified. We developed a quantitative framework based on a phytomer-based three-dimensional canopy photosynthesis model and simulated light capture and canopy photosynthetic production in scenarios with and without heterosis in either canopy structure or leaf photosynthetic capacity. The accumulated above-ground biomass of Jingnongke728 was 39% and 31% higher than its male parent, Jing2416, and female parent, JingMC01, while accumulated photosynthetically active radiation was 23% and 14% higher, correspondingly, leading to an increase of 13% and 17% in radiation use efficiency. The increasing post-silking radiation use efficiency was mainly attributed to leaf photosynthetic improvement, while the dominant contributing factor differs for male and female parents for heterosis in post-silking yield formation. This quantitative framework illustrates the potential to identify the key traits related to yield and radiation use efficiency and helps breeders to make selections for higher yield and photosynthetic efficiency.

An enhanced network of energy metabolism, lysine acetylation, and growth-promoting protein accumulation is associated with heterosis in elite hybrid rice

Heterosis refers to the superior performance of a hybrid compared with its parental lines. Although several genetic and molecular models have been proposed to explain heterosis, it remains unclear how hybrid cells integrate complementary gene expression or activity to drive heterotic growth. In this work, we show that accumulation of growth-promoting and energy metabolism proteins, enhanced energy metabolism activities, and increased protein lysine acetylation were associated with superior growth of the panicle meristem in the elite hybrid rice Shanyou 63 relative to its parental varieties. Metabolism of nuclear/cytosolic acetyl-coenzyme A was also enhanced in the hybrid, which paralleled increases in histone H3 acetylation to selectively target the expression of growth-promoting and metabolic genes. Lysine acetylation of cellular proteins, including TARGET OF RAPAMYCIN complex 1, ribosomal proteins, and energy metabolism enzymes, was also augmented and/or remodeled to modulate their activities. The data indicate that an enhanced network of energy-producing metabolic activity and growth-promoting histone acetylation/gene expression in the hybrid could contribute to its superior growth rate and may constitute a model to explain heterosis.

Lnc-AIFM2-1 promotes HBV immune escape by acting as a ceRNA for miR-330-3p to regulate CD244 expression

Chronic hepatitis B (CHB) virus infection is a major risk factor for cirrhosis and hepatocellular carcinoma (HCC). Hepatitis B virus (HBV) immune escape is regulated by the exhaustion of virus-specific CD8+ T cells, which is associated with abnormal expression of negative regulatory molecule CD244. However, the underlying mechanisms are unclear. To investigate the important roles of non-coding RNAs play in CD244 regulating HBV immune escape, we performed microarray analysis to determine the differential expression profiles of long non-coding RNAs (lncRNAs), microRNAs (miRNAs), and mRNAs in patients with CHB and patients with spontaneous clearance of HBV. Competing endogenous RNA (ceRNA) was analyzed by bioinformatics methods and confirmed by the dual-luciferase reporter assay. Furthermore, gene silencing and overexpression experiments were used to further identify the roles of lncRNA and miRNA in HBV immune escape through CD244 regulation. The results showed that the expression of CD244 on the surface of CD8+ T cells was significantly increased in CHB patients and in the co-culture system of T cells and HBV-infected HepAD38 cells, which was accompanied by the reduction of miR-330-3p and the elevation of lnc-AIFM2-1. The down-regulated miR-330-3p induced the apoptosis of T cells by lifting the inhibition of CD244, which was reversed by miR-330-3p mimic or CD244-siRNA. Lnc-AIFM2-1 promotes the accumulation of CD244, which is mediated by decreased miR-330-3p, and then reduced the clearance ability of CD8+ T cells to HBV through regulated CD244 expression. And the injury in the ability of CD8+ T cells to clear HBV can be reversed by lnc-AIFM2-1-siRNA, miR-330-3p mimic, or CD244-siRNA. Collectively, our findings indicate that lnc-AIFM2-1 on CD244 by acting as a ceRNA of miR-330-3p contributes to HBV immune escape, which may provide novel insights into the roles of interaction networks among lncRNA, miRNA, and mRNA in HBV immune escape, highlighting potential applications of lnc-AIFM2-1 and CD244 for diagnosis and treatment in CHB.

Genome-wide analysis of transcriptome and histone modifications in Brassica napus hybrid

Although utilization of heterosis has largely improved the yield of many crops worldwide, the underlying molecular mechanism of heterosis, particularly for allopolyploids, remains unclear. Here, we compared epigenome and transcriptome data of an elite hybrid and its parental lines in three assessed tissues (seedling, flower bud, and silique) to explore their contribution to heterosis in allopolyploid B. napus. Transcriptome analysis illustrated that a small proportion of non-additive genes in the hybrid compared with its parents, as well as parental expression level dominance, might have a significant effect on heterosis. We identified histone modification (H3K4me3 and H3K27me3) variation between the parents and hybrid, most of which resulted from the differences between parents. H3K4me3 variations were positively correlated with gene expression differences among the hybrid and its parents. Furthermore, H3K4me3 and H3K27me3 were rather stable in hybridization and were mainly inherited additively in the B. napus hybrid. Together, our data revealed that transcriptome reprogramming and histone modification remodeling in the hybrid could serve as valuable resources for better understanding heterosis in allopolyploid crops.

Combined transcriptome and metabolome analysis reveals the effects of light quality on maize hybrids

BACKGROUND

Heterosis, or hybrid vigor, refers to the phenotypic superiority of an F₁ hybrid relative to its parents in terms of growth rate, biomass production, grain yield, and stress tolerance. Light is an energy source and main environmental cue with marked impacts on heterosis in plants. Research into the production applications and mechanism of heterosis has been conducted for over a century and a half, but little is known about the effect of light on plant heterosis.

RESULTS

In this study, an integrated transcriptome and metabolome analysis was performed using maize (Zea mays L.) inbred parents, B73 and Mo17, and their hybrids, B73 × Mo17 (BM) and Mo17 × B73 (MB), grown in darkness or under far-red, red, or blue light. Most differentially expressed genes (73.72-92.50%) and differentially accumulated metabolites (84.74-94.32%) exhibited non-additive effects in BM and MB hybrids. Gene Ontology analysis revealed that differential genes and metabolites were involved in glutathione transfer, carbohydrate transport, terpenoid biosynthesis, and photosynthesis. The darkness, far-red, red, and blue light treatments were all associated with phenylpropanoid-flavonoid biosynthesis by Weighted Gene Co-expression Network Analysis and Kyoto Encyclopedia of Genes and Genomes enrichment analysis. Five genes and seven metabolites related to phenylpropanoid-flavonoid biosynthesis pathway were identified as potential contributors to the interactions between maize heterosis and light conditions. Consistent with the strong mid-parent heterosis observed for metabolites, significant increases in both fresh and dry weights were found in the MB and BM hybrids compared with their inbred parents. Unexpectedly, increasing light intensity resulted in higher biomass heterosis in MB, but lower biomass heterosis in BM.

CONCLUSIONS

The transcriptomic and metabolomic results provide unique insights into the effects of light quality on gene expression patterns and genotype-environment interactions, and have implications for gene mining of heterotic loci to improve maize production.

Metabolomics-centered mining of plant metabolic diversity and function: Past decade and future perspectives

Plants are natural experts in organic synthesis, being able to generate large numbers of specific metabolites with widely varying structures that help them adapt to variable survival challenges. Metabolomics is a research discipline that integrates the capabilities of several types of research including analytical chemistry, statistics, and biochemistry. Its ongoing development provides strategies for gaining a systematic understanding of quantitative changes in the levels of metabolites. Metabolomics is usually performed by targeting either a specific cell, a specific tissue, or the entire organism. Considerable advances in science and technology over the last three decades have propelled us into the era of multi-omics, in which metabolomics, despite at an earlier developmental stage than genomics, transcriptomics, and proteomics, offers the distinct advantage of studying the cellular entities that have the greatest influence on end phenotype. Here, we summarize the state of the art of metabolite detection and identification, and illustrate these techniques with four case study applications: (i) comparing metabolite composition within and between species, (ii) assessing spatio-temporal metabolic changes during plant development, (iii) mining characteristic metabolites of plants in different ecological environments and upon exposure to various stresses, and (iv) assessing the performance of metabolomics as a means of functional gene identification , metabolic pathway elucidation, and metabolomics-assisted breeding through analyzing plant populations with diverse genetic variations. In addition, we highlight the prominent contributions of joint analyses of plant metabolomics and other omics datasets, including those from genomics, transcriptomics, proteomics, epigenomics, phenomics, microbiomes, and ion-omics studies. Finally, we discuss future directions and challenges exploiting metabolomics-centered approaches in understanding plant metabolic diversity.

Vertical-scale spatial influence of radial oxygen loss on rhizosphere microbial community in constructed wetland

Complex interactions between plants and microorganisms form the basis of constructed wetlands (CWs) for pollutant removal. In the rhizosphere, radial oxygen loss (ROL) plays a key role in the activity and abundance of functional microorganisms. However, little has been done to explore how ROL would influence the niche differentiation of microbial communities at different vertical spatial scales. We demonstrate that ROL decreases with depth, promoting an oxidation-reduction rhizosphere microecosystem in CWs. The high level of ROL in the upper layer could support the oxygen supply for aerobic bacteria (Haliangium), facilitating the COD (60%) and NH₄⁺-N (50%) removal, whereas the enrichment of denitrifiers (e.g., Hydrogenophaga and Ralstonia) and methanotrophs (Methanobaterium) in the lower layer could stimulate denitrification. The function prediction results further certified that the abundance of genes catalyzing nitrifying and denitrification processes were significantly enhanced in the upper and bottom layers, respectively, which was attributed to the oxygen concentration gradient in the rhizosphere. This study contributes to further unraveling the rhizosphere effect and enables an improved understanding of the decontamination mechanisms of CWs.

The impact of epistasis in the heterosis and combining ability analyses

The current theoretical knowledge concerning the influence of epistasis on heterosis is based on a simplified multiplicative model. The objective of this study was to assess how epistasis affects the heterosis and combining ability analyses, assuming additive model, hundreds of genes, linkage disequilibrium (LD), dominance, and seven types of digenic epistasis. We developed the quantitative genetics theory for supporting the simulation of the individual genotypic values in nine populations, the selfed populations, the 36 interpopulation crosses, 180 doubled haploids (DHs), and their 16,110 crosses, assuming 400 genes on 10 chromosomes of 200 cM. Epistasis only affects population heterosis if there is LD. Only additive × additive and dominance × dominance epistasis can affect the components of the heterosis and combining ability analyses of populations. Epistasis can have a negative impact on the heterosis and combining ability analysis of populations, leading to wrong inferences regarding the identification of superior and most divergent populations. However, this depends on the type of epistasis, percentage of epistatic genes, and magnitude of their effects. Except for duplicate genes with cumulative effects and non-epistatic genic interaction, there was a decrease in the average heterosis by increasing the percentage of epistatic genes and the magnitude of their effects. The same results are generally true for the combining ability analysis of DHs. The combining ability analyses of subsets of 20 DHs showed no significant average impact of epistasis on the identification of the most divergent ones, regardless of the number of epistatic genes and magnitude of their effects. However, a negative effect on the assessment of the superior DHs can occur assuming 100% of epistatic genes, but depending on the epistasis type and the epistatic effect magnitude.

Allele-specific expression and chromatin accessibility contribute to heterosis in tea plants (Camellia sinensis)

Heterosis is extensively used to improve crop productivity, yet its allelic and chromatin regulation remains unclear. Based on our resolved genomes of the maternal TGY and paternal HD, we analyzed the contribution of allele-specific expression (ASE) and chromatin accessibility of JGY and HGY, the artificial hybrids of oolong tea with the largest cultivated area in China. The ASE genes (ASEGs) of tea hybrids with maternal-biased were mainly related to the energy and terpenoid metabolism pathways, whereas the ASEGs with paternal-biased tend to be enriched in glutathione metabolism, and these parental bias of hybrids may coordinate and lead to the acquisition of heterosis in more biological pathways. ATAC-seq results showed that hybrids have significantly higher accessible chromatin regions (ACRs) compared with their parents, which may confer broader and stronger transcriptional activity of genes in hybrids. The number of ACRs with significantly increased accessibility in hybrids was much greater than decreased, and the associated alleles were also affected by differential ACRs across different parents, suggesting enhanced positive chromatin regulation and potential genetic effects in hybrids. Core ASEGs of terpene and purine alkaloid metabolism pathways with significant positive heterosis have greater chromatin accessibility in hybrids, and were potentially regulated by several members of the MYB, DOF and TRB families. The binding motif of CsMYB85 in the promoter ACR of the rate-limiting enzyme CsDXS was verified by DAP-seq. These results suggest that higher numbers and more accessible ACRs in hybrids contribute to the regulation of ASEGs, thereby affecting the formation of heterotic metabolites.

Integrated analysis of DNA methylation, transcriptome, and global metabolites in interspecific heterotic Capsicum F1 hybrid

Hybrid breeding is one of the efficacious methods of crop improvement. Here, we report our work towards understanding the molecular basis of F1 hybrid heterosis from Capsicum chinense and C. frutescens cross. Bisulfite sequencing identified a total of 70597 CG, 108797 CHG, and 38418 CHH differentially methylated regions (DMRs) across F₁ hybrid and parents, and of these, 4891 DMRs showed higher methylation in F₁ compared to the mid-parental methylation values (MPMV). Transcriptome analysis showed higher expression of 46–55% differentially expressed genes (DE-Gs) in the F₁ hybrid. The qRT-PCR analysis of 24 DE-Gs with negative promoter methylation revealed 91.66% expression similarity with the transcriptome data. A few metabolites and 65–72% enriched genes in metabolite biosynthetic pathways showed overall increased expression in the F₁ hybrid compared to parents. These findings, taken together, provided insights into the integrated role of DNA methylation, and genes and metabolites expression in the manifestation of heterosis in Capsicum.

•

Global methylation identified significantly different proportions of mCs in hybrid

•

Of common DMRs, 33.08% showed different methylation in hybrid from the mid-parental value

•

Negatively correlated DEG pDMR-genes were enriched in metabolic pathways

•

Significant higher expression of metabolites and DE-Gs were identified in the F₁ hybrid

Evaluation of the Brewing Characteristics, Digestion Profiles, and Neuroprotective Effects of Two Typical Se-Enriched Green Teas

As a functional beverage, selenium (Se)-enriched green tea (Se-GT) has gained increasing popularity for its superior properties in promoting health. In this study, we compared the brewing characteristics, in vitro digestion profiles, and protective effects on neurotoxicity induced through the amyloid-beta (Aβ) peptide of two typical Se-GTs (Enshi Yulu (ESYL) and Ziyang Maojian (ZYMJ), representing the typical low-Se green tea and high-Se green tea, respectively). ESYL and ZYMJ showed similar chemical component leaching properties with the different brewing methods, and the optimized brewing conditions were 5 min, 90 °C, 50 mL/g, and first brewing. The antioxidant activities of the tea infusions had the strongest positive correlation with the tea polyphenols among all of the leaching substances. The tea infusions of ESYL and ZYMJ showed similar digestive behaviors, and the tea polyphenols in the tea infusions were almost totally degraded or transferred after 150 min of dynamic digestion. Studies conducted in a cell model of Alzheimer’s disease (AD) showed that the extract from the high-Se green tea was more effective for neuroprotection compared with the low-Se green tea. Overall, our results revealed the best brewing conditions and digestion behaviors of Se-GT and the great potential of Se-GT or Se-enriched green extract (Se-GTE) to be used as promising AD-preventive beverages or food ingredients.

Expression Patterns Divergence of Reciprocal F1 Hybrids Between Gossypium hirsutum and Gossypium barbadense Reveals Overdominance Mediating Interspecific Biomass Heterosis

Hybrid breeding has provided an impetus to the process and achievement of a higher yield and quality of crops. Interspecific hybridization is critical for resolving parental genetic diversity bottleneck problems. The reciprocal interspecific hybrids and their parents (Gossypium hirsutum and Gossypium barbadense) have been applied in this study to elucidate the transcription regulatory mechanism of early biomass heterosis. Phenotypically, the seed biomass, plant height over parent heterosis, leaf area over parent heterosis, and fresh and dry biomass were found to be significantly higher in hybrids than in parents. Analysis of leaf areas revealed that the one-leaf stage exhibits the most significant performance in initial vegetative growth vigor and larger leaves in hybrids, increasing the synthesis of photosynthesis compounds and enhancing photosynthesis compound synthesis. Comparative transcriptome analysis showed that transgressive down-regulation (TDR) is the main gene expression pattern in the hybrids (G. hirsutum × G. barbadense, HB), and it was found that the genes of photosystem I and Adenosine triphosphate (ATP)-binding may promote early growth vigor. Transgressive up-regulation (TUR) is the major primary gene expression pattern in the hybrids (G. barbadense × G. hirsutum, BH), and photosystem II-related genes mediated the performance of early biomass heterosis. The above results demonstrated that overdominance mediates biomass heterosis in interspecific hybrid cotton and the supervisory mechanism divergence of hybrids with different females. Photosynthesis and other metabolic process are jointly involved in controlling early biomass heterosis in interspecific hybrid cotton. The expression pattern data of transcriptome sequencing were supported using the qRT-PCR analysis. Our findings could be useful in theoretical and practical studies of early interspecific biomass heterosis, and the results provide potential resources for the theoretical and applied research on early interspecific biomass heterosis.

Comparative transcriptomic analysis reveals the molecular mechanism underlying seedling biomass heterosis in Brassica napus

BACKGROUND

Heterosis is an important biological phenomenon in which the hybrids exceed the parents in many traits. However, the molecular mechanism underlying seedling heterosis remains unclear.

RESULTS

In the present study, we analyzed the leaf transcriptomes of strong hybrids (AM, HM) and weak hybrids (CM, HW) and their parents (A, C, H, M, and W) at two periods. Phenotypically, hybrids had obvious biomass heterosis at the seedling stage, with statistically significant differences between the strong and weak hybrids. The transcriptomic analysis demonstrated that the number of differentially expressed genes (DEGs) between parents was the highest. Further analysis showed that most DEGs were biased toward parental expression. The biological processes of the two periods were significantly enriched in the plant hormone signal transduction and photosynthetic pathways. In the plant hormone signaling pathway, DEG expression was high in hybrids, with expression differences between strong and weak hybrids. In addition, DEGs related to cell size were identified. Similar changes were observed during photosynthesis. The enhanced leaf area of hybrids generated an increase in photosynthetic products, which was consistent with the phenotype of the biomass. Weighted gene co-expression network analysis of different hybrids and parents revealed that hub genes in vigorous hybrid were mainly enriched in the plant hormone signal transduction and regulation of plant hormones.

CONCLUSION

Plant hormone signaling and photosynthesis pathways, as well as differential expression of plant cell size-related genes, jointly regulate the dynamic changes between strong and weak hybrids and the generation of seedling-stage heterosis. This study may elucidate the molecular mechanism underlying early biomass heterosis and help enhance canola yield.

Inheritance of Secondary Metabolites and Gene Expression Related to Tomato Fruit Quality

Flavour and nutritional quality are important goals for tomato breeders. This study aimed to shed light upon transgressive behaviors for fruit metabolic content. We studied the metabolic contents of 44 volatile organic compounds (VOCs), 18 polyphenolics, together with transcriptome profiles in a factorial design comprising six parental lines and their 14 F1 hybrids (HF1) among which were five pairs of reciprocal HF1. After cluster analyses of the metabolome dataset and co-expression network construction of the transcriptome dataset, we characterized the mode of inheritance of each component. Both overall and per-cross mode of inheritance analyses revealed as many additive and non-additive modes of inheritance with few reciprocal effects. Up to 66% of metabolites displayed transgressions in a HF1 relative to parental values. Analysis of the modes of inheritance of metabolites revealed that: (i) transgressions were mostly of a single type whichever the cross and poorly correlated to the genetic distance between parental lines; (ii) modes of inheritance were scarcely consistent between the 14 crosses but metabolites belonging to the same cluster displayed similar modes of inheritance for a given cross. Integrating metabolome, transcriptome and modes of inheritance analyses suggested a few candidate genes that may drive important changes in fruit VOC contents.

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.

Changes in Non-Volatile and Volatile Metabolites Associated with Heterosis in Tea Plants (Camellia sinensis)

Heterosis or hybrid vigor is extensively used in plant breeding. However, the contribution of metabolites to heterosis is still elusive. Here, we systematically identified the non-volatile and volatile metabolites of two hybrids and their parents in Camellia sinensis. The metabolomics analysis showed prevalent non-additive accumulation in hybrids, among which the non-additive nucleotides, alkaloids, organic acids, and tannins contribute to the positive heterosis of hybrids, including typical inosine, guanosine, adenosine, caffeine, succinic acid, adipic acid, xylonic acid, and gallic acid. The catechins and free amino acids in hybrids showed negative heterosis compared to its maternal cultivar TGY. Furthermore, the significant accumulation of non-additive terpenes combined with the mild heterosis of other types of volatiles contributes to the aroma of tea plant hybrids. The genetics of volatiles from different parents affect the aroma of hybrids processed into oolong tea. The comprehensive heterosis of these non-additive metabolites may play an important role in the formation of desirable breeding traits for hybrids. Our results provide insights into the utilization of heterosis breeding and the regulation of heterosis metabolites in tea plants.

The multi-omics basis of potato heterosis

Heterosis is a fundamental biological phenomenon characterized by the superior performance of hybrids over their parents. Although tremendous progress has been reported in seed crops, the molecular mechanisms underlying heterosis in clonally propagated crops are largely unknown. Potato (Solanum tuberosum L.) is the most important tuber crop and an ongoing revolution is transforming potato from a clonally propagated tetraploid crop into a seed-propagated diploid hybrid potato. In our previous study, we developed the first generation of highly homozygous inbred lines of potato and hybrids with strong heterosis. Here, we integrated transcriptome, metabolome, and DNA methylation data to explore the genetic and molecular basis of potato heterosis at three developmental stages. We found that the initial establishment of heterosis in diploid potato was mainly due to dominant complementation. Flower color, male fertility, and starch and sucrose metabolism showed obvious gene dominant complementation in hybrids, and hybrids devoted more energy to primary metabolism for rapid growth. In addition, we identified ~2 700 allele-specific expression genes at each stage, which likely function in potato heterosis and might be regulated by CHH allele-specific methylation level. Our multi-omics analysis provides insight into heterosis in potato and facilitates the exploitation of heterosis in potato breeding.

Triploid Hybrid Vigor in Above-Ground Growth and Methane Fermentation Efficiency of Energy Willow

Hybrid vigor and polyploidy are genetic events widely utilized to increase the productivity of crops. Given that bioenergy usage needs to be expanded, we investigated triploid hybrid vigor in terms of the biology of biomass-related willow traits and their relevance to the control of biomethane production. To produce triploid hybrid genotypes, we crossed two female diploid Swedish cultivars (Inger, Tordis) with two male autotetraploid willow (Salix viminalis) variants (PP-E7, PP-E15). Field studies at two locations and in two successive years recorded considerable midparent heterosis (MPH%) in early shoot length that ranged between 11.14 and 68.85% and in the growth rate between 34.12 and 97.18%. The three triploid hybrids (THs) developed larger leaves than their parental cultivars, and the MPH% for their CO2 assimilation rate varied between 0.84 and 25.30%. The impact of hybrid vigor on the concentrations of plant hormones in these TH genotypes reflected essentially different hormonal statuses that depended preferentially on maternal parents. Hybrid vigor was evinced by an elevated concentration of jasmonic acid in shoot meristems of all the three THs (MPH:29.73; 67.08; 91.91%). Heterosis in auxin-type hormones, such as indole-3-acetic acid (MPH:207.49%), phenylacetic acid (MPH:223.51%), and salicylic acid (MPH:27.72%) and benzoic acid (MPH:85.75%), was detectable in the shoots of TH21/2 plants. These hormones also accumulated in their maternal Inger plants. Heterosis in cytokinin-type hormones characterized the shoots of TH3/12 and TH17/17 genotypes having Tordis as their maternal parent. Unexpectedly, we detected abscisic acid as a positive factor in the growth of TH17/17 plants with negative MPH percentages in stomatal conductance and a lower CO2 assimilation rate. During anaerobic digestion, wood raw materials from the triploid willow hybrids that provided positive MPH% in biomethane yield (6.38 and 27.87%) showed negative MPH in their acid detergent lignin contents (from -8.01 to -14.36%). Altogether, these insights into controlling factors of above-ground growth parameters of willow genotypes support the utilization of triploid hybrid vigor in willow breeding to expand the cultivation of short rotation energy trees for renewable energy production.

Plant height heterosis is quantitatively associated with expression levels of plastid ribosomal proteins

The use of hybrids is widespread in agriculture, yet the molecular basis for hybrid vigor (heterosis) remains obscure. To identify molecular components that may contribute to trait heterosis, we analyzed paired proteomic and transcriptomic data from seedling leaf and mature leaf blade tissues of maize hybrids and their inbred parents. Nuclear- and plastid-encoded subunits of complexes required for protein synthesis in the chloroplast and for the light reactions of photosynthesis were expressed above midparent and high-parent levels, respectively. Consistent with previous reports in Arabidopsis, ethylene biosynthetic enzymes were expressed below midparent levels in the hybrids, suggesting a conserved mechanism for heterosis between monocots and dicots. The ethylene biosynthesis mutant, acs2/acs6, largely phenocopied the hybrid proteome, indicating that a reduction in ethylene biosynthesis may mediate the differences between inbreds and their hybrids. To rank the relevance of expression differences to trait heterosis, we compared seedling leaf protein levels to the adult plant height of 15 hybrids. Hybrid/midparent expression ratios were most positively correlated with hybrid/midparent plant height ratios for the chloroplast ribosomal proteins. Our results show that increased expression of chloroplast ribosomal proteins in hybrid seedling leaves is mediated by reduced expression of ethylene biosynthetic enzymes and that the degree of their overexpression in seedlings can quantitatively predict adult trait heterosis.

Temperature-dependent sugar accumulation in interspecific Capsicum F1 plants showing hybrid weakness

In plants, F₁ hybrids showing hybrid weakness exhibit weaker growth than their parents. The phenotypes of hybrid weakness are often suppressed at certain temperatures. However, it is unclear whether hybrid weakness in Capsicum annuum × C. chinense is temperature-dependent or not. Our study showed that Capsicum hybrid weakness was suppressed at 30 and 35 °C and was induced at 15, 20, and 25 °C. Moreover, we investigated the time course of hybrid weakness in cell death, metabolite content, and gene expression in leaves of plants transferred to 20 °C after growing at 30 °C for 21 days. The expression of pathogen defense-related genes was upregulated at 1 day after transfer to 20 °C (DAT). Cell death was detected at 7 DAT, plant growth had almost stopped since 14 DAT, and sugars were accumulated at 42 DAT in hybrid plants. The study revealed that some sugar transporter genes, which had been upregulated since 7 DAT, were involved in sugar accumulation in Capsicum hybrid weakness. Thus, our results demonstrated that gene expression changes occur first, followed by physiological and morphological changes after induction of hybrid weakness. These responses observed in this study in Capsicum hybrid weakness are likely to be owed to plant defense responses-like reactions.

Tomato fruit quality traits and metabolite content are affected by reciprocal crosses and heterosis

Heterosis occurs when the F1s outperform their parental lines for a trait. Reciprocal hybrids are obtained by changing the cross direction of parental genotypes. Both biological phenomena could affect the external and internal attributes of fleshy fruits. This work aimed to detect reciprocal effects and heterosis in tomato (Solanum lycopersicum) fruit quality traits and metabolite content. Twelve agronomic traits and 28 metabolites identified and estimated by 1H-NMR were evaluated in five cultivars grown in two environments. Given that the genotype component was more important than the phenotype, the traits were evaluated following a full diallel mating design among those cultivars, in a greenhouse. Hybrids showed a higher phenotypic diversity than parental lines. Interestingly, the metabolites, mainly amino acids, displayed more reciprocal effects and heterosis. Agronomic traits were more influenced by general combining ability (GCA) and metabolites by specific combining ability (SCA). Furthermore, the genetic distance between parental lines was not causally related to the occurrence of reciprocal effects or heterosis. Hybrids with heterosis and a high content of metabolites linked to tomato flavour and nutritious components were obtained. Our results highlight the impact of selecting a cultivar as male or female in a cross to enhance the variability of fruit attributes through hybrids as well as the possibility to exploit heterosis for fruit composition.

The utility of metabolomics as a tool to inform maize biology

With the rise of high-throughput omics tools and the importance of maize and its products as food and bioethanol, maize metabolism has been extensively explored. Modern maize is still rich in genetic and phenotypic variation, yielding a wide range of structurally and functionally diverse metabolites. The maize metabolome is also incredibly dynamic in terms of topology and subcellular compartmentalization. In this review, we examine a broad range of studies that cover recent developments in maize metabolism. Particular attention is given to current methodologies and to the use of metabolomics as a tool to define biosynthetic pathways and address biological questions. We also touch upon the use of metabolomics to understand maize natural variation and evolution, with a special focus on research that has used metabolite-based genome-wide association studies (mGWASs).

Integrated transcript and metabolite profiling reveals coordination between biomass size and nitrogen metabolism in Arabidopsis F1 hybrids

Heterosis refers to the improved agronomic performance of F₁ hybrids relative to their parents. Although this phenomenon is widely employed to increase biomass, yield, and stress tolerance of plants, the underlying molecular mechanisms remain unclear. To dissect the metabolic fluctuations derived from genomic and/or environmental differences contributing to the improved biomass of F₁ hybrids relative to their parents, we optimized the growth condition for Arabidopsis thaliana F₁ hybrids and their parents. Modest but statistically significant increase in the biomass of F₁ hybrids was observed. Plant samples grown under the optimized condition were also utilized for integrated omics analysis to capture specific changes in the F₁ hybrids. Metabolite profiling of F₁ hybrids and parent plants was performed using gas chromatography-mass spectrometry. Among the detected 237 metabolites, 2-oxoglutarate (2-OG) and malate levels were lower and the level of aspartate was higher in the F₁ hybrids than in each parent. In addition, microarray analysis revealed that there were 44 up-regulated and 12 down-regulated genes with more than 1.5-fold changes in expression levels in the F₁ hybrid compared to each parent. Gene ontology (GO) analyses indicated that genes up-regulated in the F₁ hybrids were largely related to organic nitrogen (N) process. Quantitative PCR verified that glutamine synthetase 2 (AtGLN2) was upregulated in the F₁ hybrids, while other genes encoding enzymes in the GS-GOGAT cycle showed no significant differences between the hybrid and parent lines. These results suggested the existence of metabolic regulation that coordinates biomass and N metabolism involving AtGLN2 in F₁ hybrids.