Biological pathway expression complementation contributes to biomass heterosis in Arabidopsis

Insights into heterosis from histone modifications in the flag leaf of inter-subspecific hybrid rice

BACKGROUND

Inter-subspecific hybrid rice represents a significant breakthrough in agricultural genetics, offering higher yields and better resilience to various environmental stresses. While the utilization of these hybrids has shed light on the genetic processes underlying hybridization, understanding the molecular mechanisms driving heterosis remains a complex and ongoing challenge. Here, chromatin immunoprecipitation-sequencing (ChIP-seq) was used to analyze genome-wide profiles of H3K4me3 and H3K27me3 modifications in the inter-subspecific hybrid rice ZY19 and its parents, Z04A and ZHF1015, then combined them with the transcriptome and DNA methylation data to uncover the effects of histone modifications on gene expression and the contribution of epigenetic modifications to heterosis.

RESULTS

In the hybrid, there were 8,126 and 1,610 different peaks for H3K4me3 and H3K27me3 modifications when compared to its parents, respectively, with the majority of them originating from the parental lines. The different modifications between the hybrid and its parents were more frequently observed as higher levels in the hybrid than in the parents. In ZY19, there were 476 and 84 allele-specific genes with H3K4me3 and H3K27me3 modifications identified, representing 7.9% and 12% of the total analyzed genes, respectively. Only a small portion of genes that showed differences in parental H3K4me3 and H3K27me3 modifications which demonstrated allele-specific histone modifications (ASHM) in the hybrid. The H3K4me3 modification level in the hybrid was significantly lower compared to the parents. In the hybrid, DNA methylation occurs more frequently among histone modification target genes. Additionally, over 62.58% of differentially expressed genes (DEGs) were affected by epigenetic variations. Notably, there was a strong correlation observed between variations in H3K4me3 modifications and gene expression levels in the hybrid and its parents.

CONCLUSION

These findings highlight the substantial impact of histone modifications and DNA methylation on gene expression during hybridization. Epigenetic variations play a crucial role in controlling the differential expression of genes, with potential implications for heterosis.

Haplotype-resolved T2T genome assemblies and pangenome graph of pear reveal diverse patterns of allele-specific expression and the genomic basis of fruit quality traits

Hybrid crops often exhibit increased yield and greater resilience, yet the genomic mechanism(s) underlying hybrid vigor or heterosis remain unclear, hindering our ability to predict the expression of phenotypic traits in hybrid breeding. Here, we generated haplotype-resolved T2T genome assemblies of two pear hybrid varieties, 'Yuluxiang' (YLX) and 'Hongxiangsu' (HXS), which share the same maternal parent but differ in their paternal parents. We then used these assemblies to explore the genome-scale landscape of allele-specific expression (ASE) and create a pangenome graph for pear. ASE was observed for close to 6000 genes in both hybrid cultivars. A subset of ASE genes related to aspects of fruit quality such as sugars, organic acids, and cuticular wax were identified, suggesting their important contributions to heterosis. Specifically, Ma1, a gene regulating fruit acidity, is absent in the paternal haplotypes of HXS and YLX. A pangenome graph was built based on our assemblies and seven published pear genomes. Resequencing data for 139 cultivated pear genotypes (including 97 genotypes sequenced here) were subsequently aligned to the pangenome graph, revealing numerous structural variant hotspots and selective sweeps during pear diversification. As predicted, the Ma1 allele was found to be absent in varieties with low organic acid content, and this association was functionally validated by Ma1 overexpression in pear fruit and calli. Overall, these results reveal the contributions of ASE to fruit-quality heterosis and provide a robust pangenome reference for high-resolution allele discovery and association mapping.

High-quality genome assembly enables prediction of allele-specific gene expression in hybrid poplar

Poplar (Populus) is a well-established model system for tree genomics and molecular breeding, and hybrid poplar is widely used in forest plantations. However, distinguishing its diploid homologous chromosomes is difficult, complicating advanced functional studies on specific alleles. In this study, we applied a trio-binning design and PacBio high-fidelity long-read sequencing to obtain haplotype-phased telomere-to-telomere genome assemblies for the 2 parents of the well-studied F1 hybrid "84K" (Populus alba × Populus tremula var. glandulosa). Almost all chromosomes, including the telomeres and centromeres, were completely assembled for each haplotype subgenome apart from 2 small gaps on one chromosome. By incorporating information from these haplotype assemblies and extensive RNA-seq data, we analyzed gene expression patterns between the 2 subgenomes and alleles. Transcription bias at the subgenome level was not uncovered, but extensive-expression differences were detected between alleles. We developed machine-learning (ML) models to predict allele-specific expression (ASE) with high accuracy and identified underlying genome features most highly influencing ASE. One of our models with 15 predictor variables achieved 77% accuracy on the training set and 74% accuracy on the testing set. ML models identified gene body CHG methylation, sequence divergence, and transposon occupancy both upstream and downstream of alleles as important factors for ASE. Our haplotype-phased genome assemblies and ML strategy highlight an avenue for functional studies in Populus and provide additional tools for studying ASE and heterosis in hybrids.

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.

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.

Increased long-distance and homo-trans interactions related to H3K27me3 in Arabidopsis hybrids

In plants, the genome structure of hybrids changes compared with their parents, but the effects of these changes in hybrids remain elusive. Comparing reciprocal crosses between Col × C24 and C24 × Col in Arabidopsis using high-throughput chromosome conformation capture assay (Hi-C) analysis, we found that hybrid three-dimensional (3D) chromatin organization had more long-distance interactions relative to parents, and this was mainly located in promoter regions and enriched in genes with heterosis-related pathways. The interactions between euchromatin and heterochromatin were increased, and the compartment strength decreased in hybrids. In compartment domain (CD) boundaries, the distal interactions were more in hybrids than their parents. In the hybrids of CURLY LEAF (clf) mutants clfCol  × clfC24 and clfC24  × clfCol , the heterosis phenotype was damaged, and the long-distance interactions in hybrids were fewer than in their parents with lower H3K27me3. ChIP-seq data revealed higher levels of H3K27me3 in the region adjacent to the CD boundary and the same interactional homo-trans sites in the wild-type (WT) hybrids, which may have led to more long-distance interactions. In addition, the differentially expressed genes (DEGs) located in the boundaries of CDs and loop regions changed obviously in WT, and the functional enrichment for DEGs was different between WT and clf in the long-distance interactions and loop regions. Our findings may therefore propose a new epigenetic explanation of heterosis in the Arabidopsis hybrids and provide new insights into crop breeding and yield increase.

Dynamic roles of small RNAs and DNA methylation associated with heterosis in allotetraploid cotton (Gossypium hirsutum L.)

BACKGROUND

Heterosis is a complex phenomenon wherein the hybrids outperform their parents. Understanding the underlying molecular mechanism by which hybridization leads to higher yields in allopolyploid cotton is critical for effective breeding programs. Here, we integrated DNA methylation, transcriptomes, and small RNA profiles to comprehend the genetic and molecular basis of heterosis in allopolyploid cotton at three developmental stages.

RESULTS

Transcriptome analysis revealed that numerous DEGs responsive to phytohormones (auxin and salicylic acid) were drastically altered in F1 hybrid compared to the parental lines. DEGs involved in energy metabolism and plant growth were upregulated, whereas DEGs related to basal defense were downregulated. Differences in homoeologous gene expression in F1 hybrid were greatly reduced after hybridization, suggesting that higher levels of parental expression have a vital role in heterosis. Small RNAome and methylome studies showed that the degree of DNA methylation in hybrid is higher when compared to the parents. A substantial number of allele-specific expression genes were found to be strongly regulated by CG allele-specific methylation levels. The hybrid exhibited higher 24-nt-small RNA (siRNA) expression levels than the parents. The regions in the genome with increased levels of 24-nt-siRNA were chiefly related to genes and their flanking regulatory regions, demonstrating a possible effect of these molecules on gene expression. The transposable elements correlated with siRNA clusters in the F1 hybrid had higher methylation levels but lower expression levels, which suggest that these non-additively expressed siRNA clusters, reduced the activity of transposable elements through DNA methylation in the hybrid.

CONCLUSIONS

These multi-omics data provide insights into how changes in epigenetic mechanisms and gene expression patterns can lead to heterosis in allopolyploid cotton. This makes heterosis a viable tool in cotton breeding.

The mRNA-lncRNA landscape of multiple tissues uncovers key regulators and molecular pathways that underlie heterosis for feed intake and efficiency in laying chickens

BACKGROUND

Heterosis is routinely exploited to improve animal performance. However, heterosis and its underlying molecular mechanism for feed intake and efficiency have been rarely explored in chickens. Feed efficiency continues to be an important breeding goal trait since feed accounts for 60 to 70% of the total production costs in poultry. Here, we profiled the mRNA-lncRNA landscape of 96 samples of the hypothalamus, liver and duodenum mucosa from White Leghorn (WL), Beijing-You chicken (YY), and their reciprocal crosses (WY and YW) to elucidate the regulatory mechanisms of heterosis.

RESULTS

We observed negative heterosis for both feed intake and residual feed intake (RFI) in YW during the laying period from 43 to 46 weeks of age. Analysis of the global expression pattern showed that non-additivity was a major component of the inheritance of gene expression in the three tissues for YW but not for WY. The YW-specific non-additively expressed genes (YWG) and lncRNA (YWL) dominated the total number of non-additively expressed genes and lncRNA in the hypothalamus and duodenum mucosa. Enrichment analysis of YWG showed that mitochondria components and oxidation phosphorylation (OXPHOS) pathways were shared among the three tissues. The OXPHOS pathway was enriched by target genes for YWL with non-additive inheritance of expression in the liver and duodenum mucosa. Weighted gene co-expression network analysis revealed divergent co-expression modules associated with feed intake and RFI in the three tissues from WL, YW, and YY. Among the negatively related modules, the OXPHOS pathway was enriched by hub genes in the three tissues, which supports the critical role of oxidative phosphorylation. Furthermore, protein quantification of ATP5I was highly consistent with ATP5I expression in the liver, which suggests that, in crossbred YW, non-additive gene expression is down-regulated and decreases ATP production through oxidative phosphorylation, resulting in negative heterosis for feed intake and efficiency.

CONCLUSIONS

Our results demonstrate that non-additively expressed genes and lncRNA involved in oxidative phosphorylation in the hypothalamus, liver, and duodenum mucosa are key regulators of the negative heterosis for feed intake and RFI in layer chickens. These findings should facilitate the rational choice of suitable parents for producing crossbred chickens.

Grain yield improvement by genome editing of TaARF12 that decoupled peduncle and rachis development trajectories via differential regulation of gibberellin signalling in wheat

Plant breeding is constrained by trade-offs among different agronomic traits by the pleiotropic nature of many genes. Genes that contribute to two or more favourable traits with no penalty on yield are rarely reported, especially in wheat. Here, we describe the editing of a wheat auxin response factor TaARF12 by using CRISPR/Cas9 that rendered shorter plant height with larger spikes. Changes in plant architecture enhanced grain number per spike up to 14.7% with significantly higher thousand-grain weight and up to 11.1% of yield increase under field trials. Weighted Gene Co-Expression Network Analysis (WGCNA) of spatial-temporal transcriptome profiles revealed two hub genes: RhtL1, a DELLA domain-free Rht-1 paralog, which was up-regulated in peduncle, and TaNGR5, an organ size regulator that was up-regulated in rachis, in taarf12 plants. The up-regulation of RhtL1 in peduncle suggested the repression of GA signalling, whereas up-regulation of TaNGR5 in spike may promote GA response, a working model supported by differential expression patterns of GA biogenesis genes in the two tissues. Thus, TaARF12 complemented plant height reduction with larger spikes that gave higher grain yield. Manipulation of TaARF12 may represent a new strategy in trait pyramiding for yield improvement in wheat.

JGI Plant Gene Atlas: an updateable transcriptome resource to improve functional gene descriptions across the plant kingdom

Gene functional descriptions offer a crucial line of evidence for candidate genes underlying trait variation. Conversely, plant responses to environmental cues represent important resources to decipher gene function and subsequently provide molecular targets for plant improvement through gene editing. However, biological roles of large proportions of genes across the plant phylogeny are poorly annotated. Here we describe the Joint Genome Institute (JGI) Plant Gene Atlas, an updateable data resource consisting of transcript abundance assays spanning 18 diverse species. To integrate across these diverse genotypes, we analyzed expression profiles, built gene clusters that exhibited tissue/condition specific expression, and tested for transcriptional response to environmental queues. We discovered extensive phylogenetically constrained and condition-specific expression profiles for genes without any previously documented functional annotation. Such conserved expression patterns and tightly co-expressed gene clusters let us assign expression derived additional biological information to 64 495 genes with otherwise unknown functions. The ever-expanding Gene Atlas resource is available at JGI Plant Gene Atlas (https://plantgeneatlas.jgi.doe.gov) and Phytozome (https://phytozome.jgi.doe.gov/), providing bulk access to data and user-specified queries of gene sets. Combined, these web interfaces let users access differentially expressed genes, track orthologs across the Gene Atlas plants, graphically represent co-expressed genes, and visualize gene ontology and pathway enrichments.

Pervasive under-dominance in gene expression underlying emergent growth trajectories in Arabidopsis thaliana hybrids

BACKGROUND

Complex traits, such as growth and fitness, are typically controlled by a very large number of variants, which can interact in both additive and non-additive fashion. In an attempt to gauge the relative importance of both types of genetic interactions, we turn to hybrids, which provide a facile means for creating many novel allele combinations.

RESULTS

We focus on the interaction between alleles of the same locus, i.e., dominance, and perform a transcriptomic study involving 141 random crosses between different accessions of the plant model species Arabidopsis thaliana. Additivity is rare, consistently observed for only about 300 genes enriched for roles in stress response and cell death. Regulatory rare-allele burden affects the expression level of these genes but does not correlate with F1 rosette size. Non-additive, dominant gene expression in F1 hybrids is much more common, with the vast majority of genes (over 90%) being expressed below the parental average. Unlike in the additive genes, regulatory rare-allele burden in the dominant gene set is strongly correlated with F1 rosette size, even though it only mildly covaries with the expression level of these genes.

CONCLUSIONS

Our study underscores under-dominance as the predominant gene action associated with emergence of rosette growth trajectories in the A. thaliana hybrid model. Our work lays the foundation for understanding molecular mechanisms and evolutionary forces that lead to dominance complementation of rare regulatory alleles.

Genetic and molecular regulation of increased photosynthetic cell number contributes to leaf size heterosis in Arabidopsis

Heterosis is an important genetic phenomenon that has been observed and widely utilized in agriculture. However, the genetic and molecular bases of heterosis are unclear. Through transcriptome-wide association studies (TWAS) and expression quantitative trait locus (eQTL) analysis to integrate genome, transcriptome, and heterotic phenotype of a half-sibling Arabidopsis hybrid population, we report that the genetic and molecular bases of variations in leaf growth heterosis can be explained by the varied expression levels of growth-regulating genes resulting from distinct sets of heterozygous eQTLs carried by the half-sibling hybrids. In F1 versus parent, the degree of up-regulated gene expression in the cell cycle pathway in the shoot apex and the photosynthesis pathway in true leaf positively correlates with true leaf area heterosis level, and this is affected by the accumulation of superior heterozygous eQTLs. This was further corroborated by the major contribution of increased photosynthetic cell number to leaf area heterosis.

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.

Gene expression variations and allele-specific expression of two rice and their hybrid in caryopses at single-nucleus resolution

Seeds are an indispensable part of the flowering plant life cycle and a critical determinant of agricultural production. Distinct differences in the anatomy and morphology of seeds separate monocots and dicots. Although some progress has been made with respect to understanding seed development in Arabidopsis, the transcriptomic features of monocotyledon seeds at the cellular level are much less understood. Since most important cereal crops, such as rice, maize, and wheat, are monocots, it is essential to study transcriptional differentiation and heterogeneity during seed development at a finer scale. Here, we present single-nucleus RNA sequencing (snRNA-seq) results of over three thousand nuclei from caryopses of the rice cultivars Nipponbare and 9311 and their intersubspecies F₁ hybrid. A transcriptomics atlas that covers most of the cell types present during the early developmental stage of rice caryopses was successfully constructed. Additionally, novel specific marker genes were identified for each nuclear cluster in the rice caryopsis. Moreover, with a focus on rice endosperm, the differentiation trajectory of endosperm subclusters was reconstructed to reveal the developmental process. Allele-specific expression (ASE) profiling in endosperm revealed 345 genes with ASE (ASEGs). Further pairwise comparisons of the differentially expressed genes (DEGs) in each endosperm cluster among the three rice samples demonstrated transcriptional divergence. Our research reveals differentiation in rice caryopsis from the single-nucleus perspective and provides valuable resources to facilitate clarification of the molecular mechanism underlying caryopsis development in rice and other monocots.

Genetic architecture of inter-specific and -generic grass hybrids by network analysis on multi-omics data

BACKGROUND

Understanding the mechanisms underlining forage production and its biomass nutritive quality at the omics level is crucial for boosting the output of high-quality dry matter per unit of land. Despite the advent of multiple omics integration for the study of biological systems in major crops, investigations on forage species are still scarce.

RESULTS

Our results identified substantial changes in gene co-expression and metabolite-metabolite network topologies as a result of genetic perturbation by hybridizing L. perenne with another species within the genus (L. multiflorum) relative to across genera (F. pratensis). However, conserved hub genes and hub metabolomic features were detected between pedigree classes, some of which were highly heritable and displayed one or more significant edges with agronomic traits in a weighted omics-phenotype network. In spite of tagging relevant biological molecules as, for example, the light-induced rice 1 (LIR1), hub features were not necessarily better explanatory variables for omics-assisted prediction than features stochastically sampled and all available regressors.

CONCLUSIONS

The utilization of computational techniques for the reconstruction of co-expression networks facilitates the identification of key omic features that serve as central nodes and demonstrate correlation with the manifestation of observed traits. Our results also indicate a robust association between early multi-omic traits measured in a greenhouse setting and phenotypic traits evaluated under field conditions.

Transgressive and parental dominant gene expression and cytosine methylation during seed development in Brassica napus hybrids

KEY MESSAGE

Transcriptomic and epigenomic profiling of gene expression and small RNAs during seed and seedling development reveals expression and methylation dominance levels with implications on early stage heterosis in oilseed rape. The enhanced performance of hybrids through heterosis remains a key aspect in plant breeding; however, the underlying mechanisms are still not fully elucidated. To investigate the potential role of transcriptomic and epigenomic patterns in early expression of hybrid vigor, we investigated gene expression, small RNA abundance and genome-wide methylation in hybrids from two distant Brassica napus ecotypes during seed and seedling developmental stages using next-generation sequencing. A total of 31117, 344, 36229 and 7399 differentially expressed genes, microRNAs, small interfering RNAs and differentially methylated regions were identified, respectively. Approximately 70% of the differentially expressed or methylated features displayed parental dominance levels where the hybrid followed the same patterns as the parents. Via gene ontology enrichment and microRNA-target association analyses during seed development, we found copies of reproductive, developmental and meiotic genes with transgressive and paternal dominance patterns. Interestingly, maternal dominance was more prominent in hypermethylated and downregulated features during seed formation, contrasting to the general maternal gamete demethylation reported during gametogenesis in angiosperms. Associations between methylation and gene expression allowed identification of putative epialleles with diverse pivotal biological functions during seed formation. Furthermore, most differentially methylated regions, differentially expressed siRNAs and transposable elements were in regions that flanked genes without differential expression. This suggests that differential expression and methylation of epigenomic features may help maintain expression of pivotal genes in a hybrid context. Differential expression and methylation patterns during seed formation in an F₁ hybrid provide novel insights into genes and mechanisms with potential roles in early heterosis.

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.

Overdominant expression of genes plays a key role in root growth of tobacco hybrids

Heterosis has greatly improved the yield and quality of crops. However, previous studies often focused on improving the yield and quality of the shoot system, while research on the root system was neglected. We determined the root numbers of 12 F₁ hybrids, all of which showed strong heterosis, indicating that tobacco F₁ hybrids have general heterosis. To understand its molecular mechanism, we selected two hybrids with strong heterosis, GJ (G70 × Jiucaiping No.2) and KJ (K326 × Jiucaiping No.2), and their parents for transcriptome analysis. There were 84.22% and 90.25% of the differentially expressed genes were overdominantly expressed. The enrichment analysis of these overdominantly expressed genes showed that "Plant hormone signal transduction", "Phenylpropanoid biosynthesis", "MAPK signaling pathway - plant", and "Starch and sucrose metabolism" pathways were associated with root development. We focused on the analysis of the biosynthetic pathways of auxin(AUX), cytokinins(CTK), abscisic acid(ABA), ethylene(ET), and salicylic acid(SA), suggesting that overdominant expression of these hormone signaling pathway genes may enhance root development in hybrids. In addition, Nitab4.5_0011528g0020、Nitab4.5_0003282g0020、Nitab4.5_0004384g0070 may be the genes involved in root growth. Genome-wide comparative transcriptome analysis enhanced our understanding of the regulatory network of tobacco root development and provided new ideas for studying the molecular mechanisms of tobacco root development.

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.

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.

Dominant complementation of biological pathways in maize hybrid lines is associated with heterosis

Allele-specific expressed genes (ASEGs) are widespread in maize hybrid lines and play important roles of complementation of biological pathways in heterosis. Heterosis (hybrid vigor) is an important phenomenon with both theoretical and practical value. However, our understanding of the genetic and molecular mechanisms behind heterosis is still limited. Here, we analyzed a comprehensive dataset of maize (Zea mays L.), including RNA-seq data from three hybrid-parent triplets (HPTs) and acetylated protein data from one HPT. The gene expression patterns exhibited extensive variation between the hybrids and their parents, and a substantial number of allele-specific expressed genes (ASEGs) were identified in the hybrids. Notably, ASEGs from different HPTs were significantly enriched in various conserved pathways. The parental alleles of ASEGs with fewer deleterious single-nucleotide polymorphisms were more likely to be expressed in hybrid lines than other parental alleles. ASEGs were mainly enriched in the functional gene ontology terms protein biosynthesis, photosynthesis, and metabolism. In addition, the ASEGs across the three HPTs were involved in key photosynthetic pathways and might enhance the photosynthetic efficiency of the hybrids. These findings suggest that ASEGs involved in complementary biological pathways in maize hybrids contribute to heterosis, shedding new light on the molecular mechanism of heterosis.

Flexible response and rapid recovery strategies of the plateau forage Poa crymophila to cold and drought

Cold and drought stress are the two most severe abiotic stresses in alpine regions. Poa crymophila is widely grown in the Qinghai-Tibet Plateau with strong tolerance. Here, by profiling gene expression patterns and metabolomics-associated transcriptomics co-expression network, the acclimation of Poa crymophila to the two stresses was characterized. (1) The genes and metabolites with stress tolerance were induced by cold and drought, while those related with growth were inhibited, and most of them were restored faster after stresses disappeared. In particular, the genes for the photosynthesis system had strong resilience. (2) Additionally, cold and drought activated hypoxia and UV-B adaptation genes, indicating long-term life on the plateau could produce special adaptations. (3) Phenolamines, polyamines, and amino acids, especially N',N″,N'″-p-coumaroyl-cinnamoyl-caffeoyl spermidine, putrescine, and arginine, play key roles in harsh environments. Flexible response and quick recovery are strategies for adaptation to drought and cold in P. crymophila, accounting for its robust tolerance and resilience. In this study, we presented a comprehensive stress response profile of P. crymophila and provided many candidate genes or metabolites for future forage improvement.

Comparative transcriptome analysis between inbred lines and hybrids provides molecular insights into K+ content heterosis of tobacco (Nicotiana tabacum L.)

Potassium (K⁺) is essential for crop growth. Increasing the K⁺ content can often directly promote the improvement of crop yield and quality. Heterosis plays an important role in genetic improvement and leads to genetic gains. We found that the K⁺ content of tobacco showed significant heterosis, which is highly significant for cultivating tobacco varieties with high K⁺ content. However, the mechanism by which K⁺ content heterosis occurs in tobacco leaves is not clear. In this study, a comprehensive comparative transcriptome sequencing analysis of root samples from the hybrid G70 × GDH11 and its parental inbred lines G70 and GDH11 was performed to elucidate the importance of the root uptake capacity of K⁺ in the formation of heterosis. The results showed that 29.53% and 60.49% of the differentially expressed genes (DEGs) exhibited dominant and over-dominant expression patterns, respectively. These non-additive upregulated DEGs were significantly enriched in GO terms, such as metal ion transport and reaction, ion balance and homeostasis, ion channel activity, root meristem growth, and regulation of root hairs. The KEGG annotation results indicated that these genes were mainly involved in the pathways such as energy metabolism, carbohydrate formation, amino acid metabolism, and signal transduction. Further analysis showed that probable potassium transporter 17 (NtKT17) and potassium transporter 5-like (NtKT5), associated with potassium ion absorption, glutamate receptor 2.2-like and glutamate receptor 2.8-like, associated with ion channel activity, LOC107782957, protein detoxification 42-like, and probable glutamate carboxypeptidase 2, associated with root configuration, showed a significantly higher expression in the hybrids. These results indicated that the over-dominant expression pattern of DEGs played a key role in the heterosis of K⁺ content in tobacco leaves, and the overexpression of the genes related to K⁺ uptake, transport, and root development in hybrids helped to improve the K⁺ content of plants, thus showing the phenomenon of heterosis.

Epistasis Activation Contributes Substantially to Heterosis in Temperate by Tropical Maize Hybrids

Epistasis strongly affects the performance of superior maize hybrids. In this study, a multiple-hybrid population, consisting of three hybrid maize sets with varied interparental divergence, was generated by crossing 28 temperate and 23 tropical inbred lines with diverse genetic backgrounds. We obtained 1,154 tested hybrids. Among these tested hybrids, heterosis increased steadily as the heterotic genetic distance increased. Mid-parent heterosis was significantly higher in the temperate by tropical hybrids than in the temperate by temperate hybrids. Genome-wide prediction and association mapping was performed for grain weight per plant (GWPP) and days to silking (DTS) using 20K high-quality SNPs, showing that epistatic effects played a more prominent role than dominance effects in temperate by tropical maize hybrids. A total of 33 and 420 epistatic QTL were identified for GWPP and DTS, respectively, in the temperate by tropical hybrids. Protein-protein interaction network and gene-set enrichment analyses showed that epistatic genes were involved in protein interactions, which play an important role in photosynthesis, biological transcription pathways, and protein synthesis. We showed that the interaction of many minor-effect genes in the hybrids could activate the transcription activators of epistatic genes, resulting in a cascade of amplified yield heterosis. The multiple-hybrid population design enhanced our understanding of heterosis in maize, providing an insight into the acceleration of hybrid maize breeding by activating epistatic effects.

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.

Heterosis and Differential DNA Methylation in Soybean Hybrids and Their Parental Lines

Heterosis is an important biological phenomenon and is widely applied to increase agricultural productivity. However, the underlying molecular mechanisms of heterosis are still unclear. Here we constructed three combinations of reciprocal hybrids of soybean, and subsequently used MethylRAD-seq to detect CCGG and CCWGG (W = A or T) methylation in the whole genome of these hybrids and their parents at the middle development period of contemporary seed. We were able to prove that changes in DNA methylation patterns occurred in immature hybrid seeds and the parental variation was to some degree responsible for differential expression between the reciprocal hybrids. Non-additive differential methylation sites (DMSs) were also identified in large numbers in hybrids. Interestingly, most of these DMSs were hyper-methylated and were more concentrated in gene regions than the natural distribution of the methylated sites. Further analysis of the non-additive DMSs located in gene regions exhibited their participation in various biological processes, especially those related to transcriptional regulation and hormonal function. These results revealed DNA methylation reprogramming pattern in the hybrid soybean, which is associated with phenotypic variation and heterosis initiation.

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.

Inheritance of gene expression throughout fruit development in chili pepper

Gene expression is the primary molecular phenotype and can be estimated in specific organs or tissues at particular times. Here we analyzed genome-wide inheritance of gene expression in fruits of chili pepper (Capsicum annuum L.) in reciprocal crosses between a domesticated and a wild accession, estimating this parameter during fruit development. We defined a general hierarchical schema to classify gene expression inheritance which can be employed for any quantitative trait. We found that inheritance of gene expression is affected by both, the time of fruit development as well as the direction of the cross, and propose that such variations could be common in many developmental processes. We conclude that classification of inheritance patterns is important to have a better understanding of the mechanisms underlying gene expression regulation, and demonstrate that sets of genes with specific inheritance pattern at particular times of fruit development are enriched in different biological processes, molecular functions and cell components. All curated data and functions for analysis and visualization are publicly available as an R package.

Advances in Research on the Mechanism of Heterosis in Plants

Heterosis is a common biological phenomenon in nature. It substantially contributes to the biomass yield and grain yield of plants. Moreover, this phenomenon results in high economic returns in agricultural production. However, the utilization of heterosis far exceeds the level of theoretical research on this phenomenon. In this review, the recent progress in research on heterosis in plants was reviewed from the aspects of classical genetics, parental genetic distance, quantitative trait loci, transcriptomes, proteomes, epigenetics (DNA methylation, histone modification, and small RNA), and hormone regulation. A regulatory network of various heterosis-related genes under the action of different regulatory factors was summarized. This review lays a foundation for the in-depth study of the molecular and physiological aspects of this phenomenon to promote its effects on increasing the yield of agricultural production.