De novo genome assembly and analyses of 12 founder inbred lines provide insights into maize heterosis

Genome-wide chromatin accessibility reveals transcriptional regulation of heterosis in inter-subspecific hybrid rice

The utilization of rice heterosis is essential for ensuring global food security; however, its molecular mechanism remains unclear. In this study, comprehensive analyses of accessible chromatin regions (ACRs), DNA methylation, and gene expression in inter-subspecific hybrid and its parents were performed to determine the potential role of chromatin accessibility in rice heterosis. The hybrid exhibited abundant ACRs, in which the gene ACRs and proximal ACRs were directly related to transcriptional activation rather than the distal ACRs. Regarding the dynamic accessibility contribution of the parents, paternal ZHF1015 transmitted a greater number of ACRs to the hybrid. Accessible genotype-specific target genes were enriched with overrepresented transcription factors, indicating a unique regulatory network of genes in the hybrid. Compared with its parents, the differentially accessible chromatin regions with upregulated chromatin accessibility were much greater than those with downregulated chromatin accessibility, reflecting a stronger regulation in the hybrid. Furthermore, DNA methylation levels were negatively correlated with ACR intensity, and genes were strongly affected by CHH methylation in the hybrid. Chromatin accessibility positively regulated the overall expression level of each genotype. ACR-related genes with maternal Z04A-bias allele-specific expression tended to be enriched during carotenoid biosynthesis, whereas paternal ZHF1015-bias genes were more active in carbohydrate metabolism. Our findings provide a new perspective on the mechanism of heterosis based on chromatin accessibility in inter-subspecific hybrid rice.

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.

Potential Roles of the GRF Transcription Factors in Sorghum Internodes during Post-Reproductive Stages

Growth-regulating factor (GRF) is a plant-specific family of transcription factors crucial for meristem development and plant growth. Sorghum (Sorghum bicolor L. Moench) is a cereal species widely used for food, feed and fuel. While sorghum stems are important biomass components, the regulation of stem development and the carbohydrate composition of the stem tissues remain largely unknown. Here, we identified 11 SbGRF-encoding genes and found the SbGRF expansion driven by whole-genome duplication events. By comparative analyses of GRFs between rice and sorghum, we demonstrated the divergence of whole-genome duplication (WGD)-derived OsGRFs and SbGRFs. A comparison of SbGRFs' expression profiles supports that the WGD-duplicated OsGRFs and SbGRFs experienced distinct evolutionary trajectories, possibly leading to diverged functions. RNA-seq analysis of the internode tissues identified several SbGRFs involved in internode elongation, maturation and cell wall metabolism. We constructed co-expression networks with the RNA-seq data of sorghum internodes. Network analysis discovered that SbGRF1, 5 and 7 could be involved in the down-regulation of the biosynthesis of cell wall components, while SbGRF4, 6, 8 and 9 could be associated with the regulation of cell wall loosening, reassembly and/or starch biosynthesis. In summary, our genome-wide analysis of SbGRFs reveals the distinct evolutionary trajectories of WGD-derived SbGRF pairs. Importantly, expression analyses highlight previously unknown functions of several SbGRFs in internode elongation, maturation and the potential involvement in the metabolism of the cell wall and starch during post-anthesis stages.

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.

Genome size and chromosome number are critical metrics for accurate genome assembly assessment in Eukaryota

The number of genome assemblies has rapidly increased in recent history, with NCBI databases reaching over 41,000 eukaryotic genome assemblies across about 2,300 species. Increases in read length and improvements in assembly algorithms have led to increased contiguity and larger genome assemblies. While this number of assemblies is impressive, only about a third of these assemblies have corresponding genome size estimations for their respective species on publicly available databases. In this paper, genome assemblies are assessed regarding their total size compared to their respective publicly available genome size estimations. These deviations in size are assessed related to genome size, kingdom, sequencing platform, and standard assembly metrics, such as N50 and BUSCO values. A large proportion of assemblies deviate from their estimated genome size by more than 10%, with increasing deviations in size with increased genome size, suggesting nonprotein coding and structural DNA may be to blame. Modest differences in performance of sequencing platforms are noted as well. While standard metrics of genome assessment are more likely to indicate an assembly approaching the estimated genome size, much of the variation in this deviation in size is not explained with these raw metrics. A new, proportional N50 metric is proposed, in which N50 values are made relative to the average chromosome size of each species. This new metric has a stronger relationship with complete genome assemblies and, due to its proportional nature, allows for a more direct comparison across assemblies for genomes with variation in sizes and architectures.

Impacts of parental genomic divergence in non-syntenic regions on cotton heterosis

INTRODUCTION

Heterosis has revolutionized crop breeding, enhancing global agricultural production. However, the mechanisms underlying heterosis remain obscure. Xiangzamian 2# (XZM2), a super hybrid upland cotton (Gossypium hirsutum L.) characterized by high-yield heterosis, has been developed and extensively planted in China.

OBJECTIVES

We conducted a systematic analysis of CRI12 and J8891, two parents of XZM2. We aimed to reveal the precise genetic information and the role of non-syntenic divergence in shaping heterosis, laying a foundation for advancing understanding of heterosis.

METHODS

We de novo assembled high-quality genomes of CRI12 and J8891, and further uncovered abundant genetic variations and non-syntenic regions between the parents. Whole-genome comparison, association analysis, transcriptomic analysis and relative identity-by-descent (rIBD) estimation were conducted to identify structural variations (SVs) and introgressions within non-syntenic blocks and to analyze their impacts on promoting heterosis.

RESULTS

Parental genetic divergence increased in non-syntenic regions. Furthermore, these regions, accounting for only 16.71% of the total genome, contained more loci with significantly higher heterotic effects, far exceeding the syntenic background. SVs covered 97.26% of non-syntenic sequences and caused widespread gene expression differences in these regions, driving dynamic complementation of gene expression in the hybrid. A set of SVs were responsible for trait improvement and had positive effects on heterosis, contributing larger heritability than short variations. We characterized numerous parental-specific introgressions from G. barbadense. Specifically, a functional introgression segment within non-syntenic blocks introduced an elite haplotype, which significantly increased lint yield and enhanced heterosis.

CONCLUSION

Our study clarified non-syntenic regions to harbor more loci with higher heterotic effects, revealed their importance in promoting heterosis and supported the crucial role of genetic complementation in heterosis. SVs and introgressions were identified as key factors responsible for non-syntenic divergence between the parents. They had important effects on gene expression and trait improvement, positively contributing to heterosis.

Plant pangenomes for crop improvement, biodiversity and evolution

Plant genome sequences catalogue genes and the genetic elements that regulate their expression. Such inventories further research aims as diverse as mapping the molecular basis of trait diversity in domesticated plants or inquiries into the origin of evolutionary innovations in flowering plants millions of years ago. The transformative technological progress of DNA sequencing in the past two decades has enabled researchers to sequence ever more genomes with greater ease. Pangenomes - complete sequences of multiple individuals of a species or higher taxonomic unit - have now entered the geneticists' toolkit. The genomes of crop plants and their wild relatives are being studied with translational applications in breeding in mind. But pangenomes are applicable also in ecological and evolutionary studies, as they help classify and monitor biodiversity across the tree of life, deepen our understanding of how plant species diverged and show how plants adapt to changing environments or new selection pressures exerted by human beings.

Genome-Wide Association Study Identifies Quantitative Trait Loci and Candidate Genes Involved in Deep-Sowing Tolerance in Maize (Zea mays L.)

Deep sowing is an efficient strategy for maize to ensure the seedling emergence rate under adverse conditions such as drought or low temperatures. However, the genetic basis of deep-sowing tolerance-related traits in maize remains largely unknown. In this study, we performed a genome-wide association study on traits related to deep-sowing tolerance, including mesocotyl length (ML), coleoptile length (CL), plumule length (PL), shoot length (SL), and primary root length (PRL), using 255 maize inbred lines grown in three different environments. We identified 23, 6, 4, and 4 quantitative trait loci (QTLs) associated with ML, CL, PL, and SL, respectively. By analyzing candidate genes within these QTLs, we found a γ-tubulin-containing complex protein, ZmGCP2, which was significantly associated with ML, PL, and SL. Loss of function of ZmGCP2 resulted in decreased PL, possibly by affecting the cell elongation, thus affecting SL. Additionally, we identified superior haplotypes and allelic variations of ZmGCP2 with a longer PL and SL, which may be useful for breeding varieties with deep-sowing tolerance to improve maize cultivation.

A higher-yield hybrid rice is achieved by assimilating a dominant heterotic gene in inbred parental lines

The exploitation of heterosis to integrate parental advantages is one of the fastest and most efficient ways of rice breeding. The genomic architecture of heterosis suggests that the grain yield is strongly correlated with the accumulation of numerous rare superior alleles with positive dominance. However, the improvements in yield of hybrid rice have shown a slowdown or even plateaued due to the limited availability of complementary superior alleles. In this study, we achieved a considerable increase in grain yield of restorer lines by inducing an alternative splicing event in a heterosis gene OsMADS1 through CRISPR-Cas9, which accounted for approximately 34.1%-47.5% of yield advantage over their corresponding inbred rice cultivars. To achieve a higher yield in hybrid rice, we crossed the gene-edited restorer parents harbouring OsMADS1GW3p6 with the sterile lines to develop new rice hybrids. In two-line hybrid rice Guang-liang-you 676 (GLY676), the yield of modified hybrids carrying the homozygous heterosis gene OsMADS1GW3p6 significantly exceeded that of the original hybrids with heterozygous OsMADS1. Similarly, the gene-modified F1 hybrids with heterozygous OsMADS1GW3p6 increased grain yield by over 3.4% compared to the three-line hybrid rice Quan-you-si-miao (QYSM) with the homozygous genotype of OsMADS1. Our study highlighted the great potential in increasing the grain yield of hybrid rice by pyramiding a single heterosis gene via CRISPR-Cas9. Furthermore, these results demonstrated that the incomplete dominance of heterosis genes played a major role in yield-related heterosis and provided a promising strategy for breeding higher-yielding rice varieties above what is currently achievable.

Trans-crop applications of atypical R genes for multipathogen resistance

Genetic resistance to plant diseases is essential for global food security. Significant progress has been achieved for plant disease-resistance (R) genes comprising nucleotide-binding domain, leucine-rich repeat-containing receptors (NLRs), and membrane-localized receptor-like kinases or proteins (RLKs/RLPs), which we refer to as typical R genes. However, there is a knowledge gap in how non-receptor-type or atypical R genes contribute to plant immunity. Here, we summarize resources and technologies facilitating the study of atypical R genes, examine diverse atypical R proteins for broad-spectrum resistance, and outline potential approaches for trans-crop applications of atypical R genes. Studies of atypical R genes are important for a holistic understanding of plant immunity and the development of novel strategies in disease control and crop improvement.

Natural variations of heterosis-related allele-specific expression genes in promoter regions lead to allele-specific expression in maize

BACKGROUND

Heterosis has successfully enhanced maize productivity and quality. Although significant progress has been made in delineating the genetic basis of heterosis, the molecular mechanisms underlying its genetic components remain less explored. Allele-specific expression (ASE), the imbalanced expression between two parental alleles in hybrids, is increasingly being recognized as a factor contributing to heterosis. ASE is a complex process regulated by both epigenetic and genetic variations in response to developmental and environmental conditions.

RESULTS

In this study, we explored the differential characteristics of ASE by analyzing the transcriptome data of two maize hybrids and their parents under four light conditions. On the basis of allele expression patterns in different hybrids under various conditions, ASE genes were divided into three categories: bias-consistent genes involved in basal metabolic processes in a functionally complementary manner, bias-reversal genes adapting to the light environment, and bias-specific genes maintaining cell homeostasis. We observed that 758 ASE genes (ASEGs) were significantly overlapped with heterosis quantitative trait loci (QTLs), and high-frequency variations in the promoter regions of heterosis-related ASEGs were identified between parents. In addition, 10 heterosis-related ASEGs participating in yield heterosis were selected during domestication.

CONCLUSIONS

The comprehensive analysis of ASEGs offers a distinctive perspective on how light quality influences gene expression patterns and gene-environment interactions, with implications for the identification of heterosis-related ASEGs to enhance maize yield.

Enhanced pan-genomic resources at the maize genetics and genomics database

Pan-genomes, encompassing the entirety of genetic sequences found in a collection of genomes within a clade, are more useful than single reference genomes for studying species diversity. This is especially true for a species like Zea mays, which has a particularly diverse and complex genome. Presenting pan-genome data, analyses, and visualization is challenging, especially for a diverse species, but more so when pan-genomic data is linked to extensive gene model and gene data, including classical gene information, markers, insertions, expression and proteomic data, and protein structures as is the case at MaizeGDB. Here, we describe MaizeGDB's expansion to include the genic subset of the Zea pan-genome in a pan-gene data center featuring the maize genomes hosted at MaizeGDB, and the outgroup teosinte Zea genomes from the Pan-Andropoganeae project. The new data center offers a variety of browsing and visualization tools, including sequence alignment visualization, gene trees and other tools, to explore pan-genes in Zea that were calculated by the pipeline Pandagma. Combined, these data will help maize researchers study the complexity and diversity of Zea, and to use the comparative functions to validate pan-gene relationships for a selected gene model.

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.

HEMU: An integrated comparative genomics database and analysis platform for Andropogoneae grasses

This study reports an online database and analysis platform HEMU, which integrates 75 genome assemblies from 20 unique species, large amounts of multi-omics data, and six sophisticated analysis toolkits. HEMU will facilitate comparative genomics analysis within the tribe Andropogoneae.

The OsWRKY72-OsAAT30/OsGSTU26 module mediates reactive oxygen species scavenging to drive heterosis for salt tolerance in hybrid rice

Hybrid rice (Oryza sativa) generally outperforms its inbred parents in yield and stress tolerance, a phenomenon termed heterosis, but the underlying mechanism is not completely understood. Here, we combined transcriptome, proteome, physiological, and heterosis analyses to examine the salt response of super hybrid rice Chaoyou1000 (CY1000). In addition to surpassing the mean values for its two parents (mid-parent heterosis), CY1000 exhibited a higher reactive oxygen species scavenging ability than both its parents (over-parent heterosis or heterobeltiosis). Nonadditive expression and allele-specific gene expression assays showed that the glutathione S-transferase gene OsGSTU26 and the amino acid transporter gene OsAAT30 may have major roles in heterosis for salt tolerance, acting in an overdominant fashion in CY1000. Furthermore, we identified OsWRKY72 as a common transcription factor that binds and regulates OsGSTU26 and OsAAT30. The salt-sensitive phenotypes were associated with the OsWRKY72paternal genotype or the OsAAT30maternal genotype in core rice germplasm varieties. OsWRKY72paternal specifically repressed the expression of OsGSTU26 under salt stress, leading to salinity sensitivity, while OsWRKY72maternal specifically repressed OsAAT30, resulting in salinity tolerance. These results suggest that the OsWRKY72-OsAAT30/OsGSTU26 module may play an important role in heterosis for salt tolerance in an overdominant fashion in CY1000 hybrid rice, providing valuable clues to elucidate the mechanism of heterosis for salinity tolerance in hybrid rice.

Technology-enabled great leap in deciphering plant genomes

Plant genomes provide essential and vital basic resources for studying many aspects of plant biology and applications (for example, breeding). From 2000 to 2020, 1,144 genomes of 782 plant species were sequenced. In the past three years (2021-2023), 2,373 genomes of 1,031 plant species, including 793 newly sequenced species, have been assembled, representing a great leap. The 2,373 newly assembled genomes, of which 63 are telomere-to-telomere assemblies and 921 have been generated in pan-genome projects, cover the major phylogenetic clades. Substantial advances in read length, throughput, accuracy and cost-effectiveness have notably simplified the achievement of high-quality assemblies. Moreover, the development of multiple software tools using different algorithms offers the opportunity to generate more complete and complex assemblies. A database named N3: plants, genomes, technologies has been developed to accommodate the metadata associated with the 3,517 genomes that have been sequenced from 1,575 plant species since 2000. We also provide an outlook for emerging opportunities in plant genome sequencing.

Potential Regulatory Networks and Heterosis for Flavonoid and Terpenoid Contents in Pak Choi: Metabolomic and Transcriptome Analyses

Pak choi exhibits a diverse color range and serves as a rich source of flavonoids and terpenoids. However, the mechanisms underlying the heterosis and coordinated regulation of these compounds-particularly isorhamnetin-remain unclear. This study involved three hybrid combinations and the detection of 528 metabolites from all combinations, including 26 flavonoids and 88 terpenoids, through untargeted metabolomics. Analysis of differential metabolites indicated that the heterosis for the flavonoid and terpenoid contents was parent-dependent, and positive heterosis was observed for isorhamnetin in the two hybrid combinations (SZQ, 002 and HMG, ZMG). Moreover, there was a high transcription level of flavone 3'-O-methyltransferase, which is involved in isorhamnetin biosynthesis. The third group was considered the ideal hybrid combination for investigating the heterosis of flavonoid and terpenoid contents. Transcriptome analysis identified a total of 12,652 DEGs (TPM > 1) in various groups that were used for comparison, and DEGs encoding enzymes involved in various categories, including "carotenoid bio-synthesis" and "anthocyanin biosynthesis", were enriched in the hybrid combination (SZQ, 002). Moreover, the category of anthocyanin biosynthesis also was enriched in the hybrid combination (HMG, ZMG). The flavonoid pathway demonstrated more differential metabolites than the terpenoid pathway did. The WGCNA demonstrated notable positive correlations between the dark-green modules and many flavonoids and terpenoids. Moreover, there were 23 ERF genes in the co-expression network (r ≥ 0.90 and p < 0.05). Thus, ERF genes may play a significant role in regulating flavonoid and terpenoid biosynthesis. These findings enhance our understanding of the heterosis and coordinated regulation of flavonoid and terpenoid biosynthesis in pak choi, offering insights for genomics-based breeding improvements.

Striking a growth-defense balance: Stress regulators that function in maize development

Maize (Zea mays) cultivation is strongly affected by both abiotic and biotic stress, leading to reduced growth and productivity. It has recently become clear that regulators of plant stress responses, including the phytohormones abscisic acid (ABA), ethylene (ET), and jasmonic acid (JA), together with reactive oxygen species (ROS), shape plant growth and development. Beyond their well established functions in stress responses, these molecules play crucial roles in balancing growth and defense, which must be finely tuned to achieve high yields in crops while maintaining some level of defense. In this review, we provide an in-depth analysis of recent research on the developmental functions of stress regulators, focusing specifically on maize. By unraveling the contributions of these regulators to maize development, we present new avenues for enhancing maize cultivation and growth while highlighting the potential risks associated with manipulating stress regulators to enhance grain yields in the face of environmental challenges.

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.

A graph-based pan-genome of Brassica oleracea provides new insights into its domestication and morphotype diversification

The domestication of Brassica oleracea has resulted in diverse morphological types with distinct patterns of organ development. Here we report a graph-based pan-genome of B. oleracea constructed from high-quality genome assemblies of different morphotypes. The pan-genome harbors over 200 structural variant hotspot regions enriched in auxin- and flowering-related genes. Population genomic analyses revealed that early domestication of B. oleracea focused on leaf or stem development. Gene flows resulting from agricultural practices and variety improvement were detected among different morphotypes. Selective-sweep and pan-genome analyses identified an auxin-responsive small auxin up-regulated RNA gene and a CLAVATA3/ESR-RELATED family gene as crucial players in leaf-stem differentiation during the early stage of B. oleracea domestication and the BoKAN1 gene as instrumental in shaping the leafy heads of cabbage and Brussels sprouts. Our pan-genome and functional analyses further revealed that variations in the BoFLC2 gene play key roles in the divergence of vernalization and flowering characteristics among different morphotypes, and variations in the first intron of BoFLC3 are involved in fine-tuning the flowering process in cauliflower. This study provides a comprehensive understanding of the pan-genome of B. oleracea and sheds light on the domestication and differential organ development of this globally important crop species.

Genome-Wide Identification and Expression Analyses of the FAR1/FHY3 Gene Family Provide Insight into Inflorescence Development in Maize

As transcription factors derived from transposase, FAR-RED IMPAIRED RESPONSE1 (FAR1) and its homolog FHY3 play crucial roles in the regulation of light signaling and various stress responses by coordinating the expression of downstream target genes. Despite the extensive investigation of the FAR1/FHY3 family in Arabidopsis thaliana and other species, a comprehensive examination of these genes in maize has not been conducted thus far. In this study, we employed a genomic mining approach to identify 16 ZmFAR1 genes in the maize inbred line B73, which were further classified into five subgroups based on their phylogenetic relationships. The present study characterized the predicted polypeptide sequences, molecular weights, isoelectric points, chromosomal distribution, gene structure, conserved motifs, subcellular localizations, phylogenetic relationships, and cis-regulatory elements of all members belonging to the ZmFAR1 family. Furthermore, the tissue-specific expression of the 16 ZmFAR1 genes was analyzed using RNA-seq, and their expression patterns under far-red light conditions were validated in the ear and tassel through qRT-qPCR. The observed highly temporal and spatial expression patterns of these ZmFAR1 genes were likely associated with their specific functional capabilities under different light conditions. Further analysis revealed that six ZmFAR1 genes (ZmFAR1-1, ZmFAR1-10, ZmFAR1-11, ZmFAR1-12, ZmFAR1-14, and ZmFAR1-15) exhibited a response to simulated shading treatment and actively contributed to the development of maize ears. Through the integration of expression quantitative trait loci (eQTL) analyses and population genetics, we identified the presence of potential causal variations in ZmFAR1-14 and ZmFAR1-9, which play a crucial role in regulating the kernel row number and kernel volume weight, respectively. In summary, this study represents the initial identification and characterization of ZmFAR1 family members in maize, uncovering the functional variation in candidate regulatory genes associated with the improvement of significant agronomic traits during modern maize breeding.

Telomere-to-Telomere Genome Assembly of Tibetan Medicinal Mushroom Ganoderma leucocontextum and the First Copia Centromeric Retrotransposon in Macro-Fungi Genome

A complete telomere-to-telomere (T2T) genome has been a longstanding goal in the field of genomic research. By integrating high-coverage and precise long-read sequencing data using multiple assembly strategies, we present here the first T2T gap-free genome assembly of Ganoderma leucocontextum strain GL72, a Tibetan medicinal mushroom. The T2T genome, with a size of 46.69 Mb, consists 13 complete nuclear chromosomes and typical telomeric repeats (CCCTAA)n were detected at both ends of 13 chromosomes. The high mapping rate, uniform genome coverage, a complete BUSCOs of 99.7%, and base accuracy exceeding 99.999% indicate that this assembly represents the highest level of completeness and quality. Regions characterized by distinct structural attributes, including highest Hi-C interaction intensity, high repeat content, decreased gene density, low GC content, and minimal or no transcription levels across all chromosomes may represent potential centromeres. Sequence analysis revealed the first Copia centromeric retrotransposon in macro-fungi genome. Phylogenomic analysis identified that G. leucocontextum and G. tsugae diverged from the other Ganoderma species approximately 9.8-17.9 MYA. The prediction of secondary metabolic clusters confirmed the capability of this fungus to produce a substantial quantity of metabolites. This T2T gap-free genome will contribute to the genomic 'dark matter' elucidation and server as a great reference for genetics, genomics, and evolutionary studies of G. leucocontextum.

Genetic architecture of ear traits based on association mapping and co-expression networks in maize inbred lines and hybrids

Ear traits are key contributors to grain yield in maize; therefore, exploring their genetic basis facilitates the improvement of grain yield. However, the underlying molecular mechanisms of ear traits remain obscure in both inbred lines and hybrids. Here, two association panels, respectively, comprising 203 inbred lines (IP) and 246 F1 hybrids (HP) were employed to identify candidate genes for six ear traits. The IP showed higher phenotypic variation and lower phenotypic mean than the HP for all traits, except ear tip-barrenness length. By conducting a genome-wide association study (GWAS) across multiple environments, 101 and 228 significant single-nucleotide polymorphisms (SNPs) associated with six ear traits were identified in the IP and HP, respectively. Of these significant SNPs identified in the HP, most showed complete-incomplete dominance and over-dominance effects for each ear trait. Combining a gene co-expression network with GWAS results, 186 and 440 candidate genes were predicted in the IP and HP, respectively, including known ear development genes ids1 and sid1. Of these, nine candidate genes were detected in both populations and expressed in maize ear and spikelet tissues. Furthermore, two key shared genes (GRMZM2G143330 and GRMZM2G171139) in both populations were found to be significantly associated with ear traits in the maize Goodman diversity panel with high-density variations. These findings advance our knowledge of the genetic architecture of ear traits between inbred lines and hybrids and provide a valuable resource for the genetic improvement of ear traits in maize.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01426-9.

Transcriptome-wide identification and characterization of genes exhibit allele-specific imprinting in maize embryo and endosperm

BACKGROUND

Genomic imprinting refers to a subset of genes that are expressed from only one parental allele during seed development in plants. Studies on genomic imprinting have revealed that intraspecific variations in genomic imprinting expression exist in naturally genetic varieties. However, there have been few studies on the functional analysis of allele-specific imprinted genes.

RESULTS

Here, we generated three reciprocal crosses among the B73, Mo17 and CAU5 inbred lines. Based on the transcriptome-wide analysis of allele-specific expression using RNA sequencing technology, 305 allele-specific imprinting genes (ASIGs) were identified in embryos, and 655 ASIGs were identified in endosperms from three maize F1 hybrids. Of these ASIGs, most did not show consistent maternal or paternal bias between the same tissue from different hybrids or different tissues from one hybrid cross. By gene ontology (GO) analysis, five and eight categories of GO exhibited significantly higher functional enrichments for ASIGs identified in embryo and endosperm, respectively. These functional categories indicated that ASIGs are involved in intercellular nutrient transport, signaling pathways, and transcriptional regulation of kernel development. Finally, the mutation and overexpression of one ASIG (Zm305) affected the length and width of the kernel.

CONCLUSION

In this study, our data will be helpful in gaining further knowledge of genes exhibiting allele-specific imprinting patterns in seeds. The gain- and loss-of-function phenotypes of ASIGs associated with agronomically important seed traits provide compelling evidence for ASIGs as crucial targets to optimize seed traits in crop plants.

A transcription factor ZmGLK36 confers broad resistance to maize rough dwarf disease in cereal crops

Maize rough dwarf disease (MRDD), caused by maize rough dwarf virus (MRDV) or rice black-streaked dwarf virus (RBSDV), seriously threatens worldwide production of all major cereal crops, including maize, rice, wheat and barley. Here we report fine mapping and cloning of a previously reported major quantitative trait locus (QTL) (qMrdd2) for RBSDV resistance in maize. Subsequently, we show that qMrdd2 encodes a G2-like transcription factor named ZmGLK36 that promotes resistance to RBSDV by enhancing jasmonic acid (JA) biosynthesis and JA-mediated defence response. We identify a 26-bp indel located in the 5' UTR of ZmGLK36 that contributes to differential expression and resistance to RBSDV in maize inbred lines. Moreover, we show that ZmDBF2, an AP2/EREBP family transcription factor, directly binds to the 26-bp indel and represses ZmGLK36 expression. We further demonstrate that ZmGLK36 plays a conserved role in conferring resistance to RBSDV in rice and wheat using transgenic or marker-assisted breeding approaches. Our results provide insights into the molecular mechanisms of RBSDV resistance and effective strategies to breed RBSDV-resistant cereal crops.

Structure and function of rice hybrid genomes reveal genetic basis and optimal performance of heterosis

Exploitation of crop heterosis is crucial for increasing global agriculture production. However, the quantitative genomic analysis of heterosis was lacking, and there is currently no effective prediction tool to optimize cross-combinations. Here 2,839 rice hybrid cultivars and 9,839 segregation individuals were resequenced and phenotyped. Our findings demonstrated that indica-indica hybrid-improving breeding was a process that broadened genetic resources, pyramided breeding-favorable alleles through combinatorial selection and collaboratively improved both parents by eliminating the inferior alleles at negative dominant loci. Furthermore, we revealed that widespread genetic complementarity contributed to indica-japonica intersubspecific heterosis in yield traits, with dominance effect loci making a greater contribution to phenotypic variance than overdominance effect loci. On the basis of the comprehensive dataset, a genomic model applicable to diverse rice varieties was developed and optimized to predict the performance of hybrid combinations. Our data offer a valuable resource for advancing the understanding and facilitating the utilization of heterosis in rice.

Unraveling Prevalence and Effects of Deleterious Mutations in Maize Elite Lines across Decades of Modern Breeding

Future breeding is likely to involve the detection and removal of deleterious alleles, which are mutations that negatively affect crop fitness. However, little is known about the prevalence of such mutations and their effects on phenotypic traits in the context of modern crop breeding. To address this, we examined the number and frequency of deleterious mutations in 350 elite maize inbred lines developed over the past few decades in China and the United States. Our findings reveal an accumulation of weakly deleterious mutations and a decrease in strongly deleterious mutations, indicating the dominant effects of genetic drift and purifying selection for the two types of mutations, respectively. We also discovered that slightly deleterious mutations, when at lower frequencies, were more likely to be heterozygous in the developed hybrids. This is consistent with complementation as a potential explanation for heterosis. Subsequently, we found that deleterious mutations accounted for more of the variation in phenotypic traits than nondeleterious mutations with matched minor allele frequencies, especially for traits related to leaf angle and flowering time. Moreover, we detected fewer deleterious mutations in the promoter and gene body regions of differentially expressed genes across breeding eras than in nondifferentially expressed genes. Overall, our results provide a comprehensive assessment of the prevalence and impact of deleterious mutations in modern maize breeding and establish a useful baseline for future maize improvement efforts.

Gene expression and expression quantitative trait loci analyses uncover natural variations underlying the improvement of important agronomic traits during modern maize breeding

Maize (Zea mays L.) is a major staple crop worldwide, and during modern maize breeding, cultivars with increased tolerance to high-density planting and higher yield per plant have contributed significantly to the increased yield per unit land area. Systematically identifying key agronomic traits and their associated genomic changes during modern maize breeding remains a significant challenge because of the complexity of genetic regulation and the interactions of the various agronomic traits, with most of them being controlled by numerous small-effect quantitative trait loci (QTLs). Here, we performed phenotypic and gene expression analyses for a set of 137 elite inbred lines of maize from different breeding eras in China. We found four yield-related traits are significantly improved during modern maize breeding. Through gene-clustering analyses, we identified four groups of expressed genes with distinct trends of expression pattern change across the historical breeding eras. In combination with weighted gene co-expression network analysis, we identified several candidate genes regulating various plant architecture- and yield-related agronomic traits, such as ZmARF16, ZmARF34, ZmTCP40, ZmPIN7, ZmPYL10, ZmJMJ10, ZmARF1, ZmSWEET15b, ZmGLN6 and Zm00001d019150. Further, by combining expression quantitative trait loci (eQTLs) analyses, correlation coefficient analyses and population genetics, we identified a set of candidate genes that might have been under selection and contributed to the genetic improvement of various agronomic traits during modern maize breeding, including a number of known key regulators of plant architecture, flowering time and yield-related traits, such as ZmPIF3.3, ZAG1, ZFL2 and ZmBES1. Lastly, we validated the functional variations in GL15, ZmPHYB2 and ZmPYL10 that influence kernel row number, flowering time, plant height and ear height, respectively. Our results demonstrates the effectiveness of our combined approaches for uncovering key candidate regulatory genes and functional variation underlying the improvement of important agronomic traits during modern maize breeding, and provide a valuable genetic resource for the molecular breeding of maize cultivars with tolerance for high-density planting.

CIMBL55: a repository for maize drought resistance alleles

Droughts threaten crop yields worldwide. Compared to other major staple cereal crops, maize (Zea mays) is especially sensitive to drought, which can cause dramatic fluctuations in its yield potential. Natural maize populations contain many superior alleles that can enhance drought resistance through complex regulatory mechanisms. We recently de novo assembled the genome of a prominent drought-resistant maize germplasm, CIMBL55, and systematically dissected the genetic basis for its drought resistance on the genome, transcriptome, and epigenome levels. These analyses revealed 65 favorable drought resistance alleles in CIMBL55. Subsequently, we genetically verified the functions of the drought resistance genes ZmABF4, ZmNAC075, and ZmRtn16 and unraveled the function of ZmRtn16 on a molecular level.

Advances of Apetala2/Ethylene Response Factors in Regulating Development and Stress Response in Maize

Apetala2/ethylene response factor (AP2/ERF) is one of the largest families of transcription factors, regulating growth, development, and stress response in plants. Several studies have been conducted to clarify their roles in Arabidopsis and rice. However, less research has been carried out on maize. In this review, we systematically identified the AP2/ERFs in the maize genome and summarized the research progress related to AP2/ERF genes. The potential roles were predicted from rice homologs based on phylogenetic and collinear analysis. The putative regulatory interactions mediated by maize AP2/ERFs were discovered according to integrated data sources, implying that they involved complex networks in biological activities. This will facilitate the functional assignment of AP2/ERFs and their applications in breeding strategy.