Bananas and plantains of the genus Musa constitute the most vital fruits and staple foods. Cultivated bananas may have originated from intraspecific and interspecific hybridizations of four wild species, namely Musa acuminata (A), M. balbisiana (B), M. schizocarpa (S), and the Australimusa species (T). Here, we appraise the advances made in banana genomics, genetics, and breeding over the past few decades. The sequencing of Musa genomes has been a major breakthrough in banana research programs, presenting unprecedented possibilities for gaining deeper insights into the evolution, domestication, breeding, and genetics of indispensable agronomic traits of bananas. Also, we delve into how these genetic facets, coupled with innovative genomic-assisted tools, including genomic selection and gene editing, propel advancements in banana breeding endeavors. Ultimately, we propose the forthcoming prospects within the domain of banana genetics and breeding.

Bananas (Musa ssp.) are globally important staple crops increasingly constrained by biotic stressors, climatic instability, and the high labor demands of cultivation. The genetic improvement of dwarf phenotypes offers a strategic pathway to enhance mechanization and reduce production costs. In this study, we have carried out whole-genome resequencing of 300 Musa accessions to analyze genome-wide allelic diversity and identify loci associated with shoot architecture. Our analysis uncovered extensive genetic variation within the A subgenome, pivotal for environmental adaptability, and detected introgression from Musa itinerans (subgroup A) into cultivated varieties (subgroup F), suggesting a broadened genetic base amenable to breeding. A genome-wide association study (GWAS) pinpointed MabHLH30 as a crucial gene associated plant stature. Functional validation confirmed MabHLH30 as a critical regulator of plant stature and leaf morphology. Leveraging this finding, we developed molecular markers for MabHLH30, enabling marker-assisted selection (MAS) to accelerate the breeding of compact, high-yielding cultivars. Collectively, these results provide a genomic framework for the targeted improvement of banana architecture and represent a valuable resource for cultivar development under diverse agroecological conditions.

In recent decades, research on model organisms has significantly increased our understanding of core biological processes in plant science. However, this focus has created a substantial knowledge bottleneck due to the limited phylogenetic and ecological spectrum covered. Gymnosperms, especially conifers, represent a molecular and ecological diversity hotspot among seed plants. Despite their importance, research on these species is notably underrepresented, primarily due to a slower pace of investigation resulting from a lack of community-based resources and databases. To fill this gap, we developed the P(inus)ra(diata)-G(ene)E(xpression) (Pra-GE)-ATLAS, which consists of several tools and two main modules: transcriptomics and proteomics, presented in this work for the forestry commercial and stress-sensitive species Pinus radiata. We have summarized and centralized all the available information to provide a comprehensive view of the gene expression landscape. To illustrate how applications of the database lead to new biological insights, we have integrated multiple regulatory layers across tissues and stressors. While stress favors the retention of small introns, harmonized alternative splicing analyses reveal that genes with conifers' iconic large introns tend to be under constitutive regulation. Furthermore, the degree of convergence between stressors differed between regulatory layers, with proteomic responses remaining highly distinctive even through intergenerational memory tolerance. Overall, the Pra-GE-ATLAS aims to narrow the distance between angiosperms and gymnosperms resources, deepening our understanding of how characteristic pine features have evolved. Pra-GE-ATLAS DB is available at: http://pra-ge-atlas.valmei.es.

Pitaya is an important perennial herbaceous fruit tree. The color of fruit determines pitaya nutritive (and attractive) value, which is considered as an important objective in breeding improvement. In this study, we reported the first telomere-to-telomere (T2T) gap-free genome of “Shuangse No. 1” pitaya (Hylocereus polyrhizus; red peel). Two high-quality genomes for “Dahong” (H. polyrhizus; red peel) and “Honghuaqinglong” (H. stenopterus; stay-green) were further assembled, aiming to explore the genetic diversity of pitaya genomes. In further analysis, we noticed a high proportion of viral contamination in pitaya tissues, which hindered the efficient utilization of transcriptomic data. To address this issue, we analyzed 111 pitaya transcriptome data from different geographic regions to characterize and separate viral components. Then we developed an efficient, novel, and universal transcript purification system for pitaya transcriptomes by applying it to 27 samples from different tissues and species, thereby enhancing the utility for transcriptomic and broader biological research. Combining the purified transcriptomic data with comparative genomic analyses, we identified HuERF72, a transcription factor (TF) that potentially regulates chlorophyll degradation in pitaya. Interaction assays and plant transformation elucidated that HuERF72 acts as a repressive TF by directly binding to the promoter of HuSGR1, a key structural gene in the chlorophyll degradation pathway. This study provides high-quality genomic resources and novel methodologies for molecular investigations in pitaya. Additionally, the proposed regulatory network advances our understanding of the transcriptional regulatory mechanisms underlying chlorophyll degradation, offering valuable insights into the genetic improvement of pitaya.

Bamboo is a fast-growing and ecologically significant plant with immense economic value due to its applications in construction, textiles, and bioenergy. However, research on bamboo has been hindered by its long vegetative period, unpredictable flowering cycles, and challenges in genetic transformation. Recent developments in advanced sequencing and genetic engineering technologies have provided new insights into bamboo's evolutionary history, developmental biology, and stress resilience, paving the way for improved conservation and sustainable utilization. This review synthesizes the latest findings on bamboo's genomics, biotechnology, and the molecular mechanisms governing its growth, development, and stress response. Key genes and regulatory pathways controlling its rapid growth, internode elongation, rhizome development, culm lignification, flowering, and abiotic stress responses have been identified through multi-omics and functional studies. Complex interactions among transcription factors, epigenetic regulators, and functionally important genes shape bamboo's unique growth characteristics. Moreover, progress in genetic engineering techniques, including clustered regularly interspaced short palindromic repeats-based genome editing, has opened new avenues for targeted genetic improvements. However, technical challenges, particularly the complexity of polyploid bamboo genomes and inefficient regeneration systems, remain significant barriers to functional studies and large-scale breeding efforts. By integrating recent genomic discoveries with advancements in biotechnology, this review proposes potential strategies to overcome existing technological limitations and to accelerate the development of improved bamboo varieties. Continued efforts in multi-omics research, gene-editing applications, and sustainable cultivation practices will be essential for harnessing bamboo as a resilient and renewable resource for the future. The review presented here not only deepens our understanding of bamboo's genetic architecture but also provides a foundation for future research aimed at optimizing its ecological and industrial potential.

Deep learning-based genomic prediction (DL-based GP) has shown promising performance compared to traditional GP methods in plant breeding, particularly in handling large, complex multi-omics data sets. However, the effective development and widespread adoption of DL-based GP still face substantial challenges, including the need for large, high-quality data sets, inconsistencies in performance benchmarking, and the integration of environmental factors. Here, we summarize the key obstacles impeding the development of DL-based GP models and propose future developing directions, such as modular approaches, data augmentation, and advanced attention mechanisms.

Avena barbata, a wild oat species within the genus Avena, is a widely used model for studying plant ecological adaptation due to its strong environmental adaptability and disease resistance, serving as a valuable genetic resource for oat improvement. Here, we phased the high-quality chromosome-level genome assembly of A. barbata (6.88 Gb, contig N50 = 53.74 Mb) into A (3.57 Gb with 47,687 genes) and B (3.31 Gb with 46,029 genes) subgenomes. Comparative genomics and phylogenomic analyses clarified the evolutionary relationships and trajectories of A, B, C and D subgenomes in Avena. We inferred that the A subgenome donor of A. barbata was Avena hirtula, while the B subgenome donor was probably an extinct diploid species closely related to Avena wiestii. Genome evolution analysis revealed the dynamic transposable element (TE) content and subgenome divergence, as well as extensive structure variations across A, B, C, and D subgenomes in Avena. Population genetic analysis of 211 A. barbata accessions from distinct ecotypes identified several candidate genes related to environmental adaptability and drought resistance. Our study provides a comprehensive genetic resource for exploring the genetic basis underlying the strong environmental adaptability of A. barbata and the molecular identification of important agronomic traits for oat breeding.

Diverse heterotic groups have been developed in China over several decades, but their genomic divergences have not been systematically studied after improvement. In this study, we performed Maize6H-60K array of 5,822 maize accessions and whole-genome re-sequencing of 150 inbred lines collected in China. Using multiple population structure analysis methods, we established a genetic boundary used to categorize heterotic groups and germplasm resources. We identified three chloroplast–cytoplasmic types that evolved during adaptation to diverse climatic environments in maize through phylogenetic and haplotype analyses. Comparative analyses revealed obvious genetic differences between heterotic groups and germplasm resources at both the chloroplast and nuclear genome levels, especially in the unique heterotic groups HG1 and HG2, which exhibited distinct regionality and genetic uniqueness. The divergent differentiation of heterotic groups from germplasm resources was driven by differential selection in specific genomic regions. Genome-wide selective sweep analysis identified core selected regions and candidate selected genes associated with traits between heterotic groups, highlighting that stress response- and plant defense-related genes were selected for environmental adaptation across a broad latitudinal range in China. Meanwhile, a genome-wide association study analysis provided evidence that core selected genes served as an important candidate gene pool with a potential role in genetic improvement. Gene exchanges among heterotic groups, which avoided the predominant heterotic patterns as much as possible, occurred to achieve population improvement during modern maize breeding. This study provides insights into the population differentiation and genetic characteristics of heterotic groups, which will facilitate the utilization of germplasm resources, the creation of novel maize germplasm, and the optimization of heterotic patterns during future maize breeding in China.

Incorporating genotype-by-environment (GE) interaction effects into genomic prediction (GP) models with multi-environment climate data can improve selection accuracy to accelerate crop breeding but has received little research attention. Here, we conducted a cross-region GP study of grain moisture content (GMC) and grain yield (GY) in maize hybrids in two major Chinese growing regions using data for 19 climatic factors across 34 environments in 2020 and 2021. Predictions were conducted in 2,126 hybrids generated from 475 maize inbred lines, using 9,355 single nucleotide polymorphism markers for genotyping. Models based on genomic best linear unbiased prediction (GBLUP) incorporating GE interaction effects of 19 climatic factors associated with day length, transpiration, temperature, and radiation (GBLUP-GE_19CF) trained on whole data set outperformed the traditional GBLUP or BayesB models in predicting GMC or GY by 10-fold cross-validation, achieving prediction accuracies of 0.731 and 0.331, respectively. To refine the climate data, we examined 84 statistical features associated with these climatic factors and identified nine factors most correlated with GMC or GY. Principal component analysis of climate data yielded nine principal components responsible for 97% of the variability in the data. Incorporating these nine factors or principal components into the GBLUP-GE framework with a similarity matrix of environments (GBLUP-GE_9CF and GBLUP-GE_PCA) provided similar prediction accuracies but could reduce the computational burden. In addition, increasing the number of test set environments in the training set from 8 to 14 increased the prediction accuracy of GBLUP-GE_19CF trained with monthly average climate data for 2020–2021. Examining prediction accuracy based on concordance, the proportion of overlapping hybrids between the top 50% of predicted and observed values for GMC and GY, indicated that concordance exceeded 50% for the GBLUP-GE_19CF model, confirming the reliability of our predictions. This study can provide practical guidance for optimizing GPs for maize breeding programs in multi-environment selection.

The publication of several high-quality genomes has contributed greatly to clarifying the evolution of citrus. However, due to their complex genetic backgrounds, the origins and evolution of many citrus species remain unclear. We assembled de novo the 294-Mbp chromosome-level genome of a more than 200-year-old primitive papeda (DYC002). Comparison between the two sets of homologous chromosomes of the haplotype-resolved genome revealed 1.2% intragenomic variations, including 1.75 million SNPs, 149,471 insertions and 154,215 deletions. Using this genome as a reference, we resequenced and performed population and phylogenetic analyses of 378 representative citrus accessions. Our study confirmed that the primary origin center of core Citrus species is in South China, particularly in the Himalaya–Hengduan Mountains. Papeda species are an ancient Citrus type compared with C. ichangensis. We found that the evolution of the Citrus genus followed two radiations through two routes (to East China and Southeast Asia) along river systems. Evidence for the origin and evolution of some individual citrus species was provided. Papeda probably played an important role in the origins of Australian finger lime, citrons, Honghe papeda and pummelos; Ichang papeda originated from Yuanjiang city of Yunnan Province, China, and C. mangshanensis has a close relationship with kumquat and Ichang papeda. Moreover, the Hunan and Guangdong Provinces of China are predicted to be the origin center of mandarin, sweet orange and sour orange. Additionally, our study revealed that fruit bitterness was significantly selected against during citrus domestication. Taken together, this study provides new insight into the origin and evolution of citrus species and may serve as a valuable genomic resource for citrus breeding and improvement.

Wild species of domesticated crops provide valuable genetic resources for resistance breeding. Prunus davidiana, a wild relative of peach with high heterozygosity and diverse stress tolerance, exhibits high resistance against aphids. However, the highly heterozygous genome of P. davidiana makes determining the underlying factors influencing resistance traits challenging. Here, we present the 501.7 Mb haplotype-resolved genome assembly of P. davidiana. Genomic comparisons of the two haplotypes revealed 18,152 structural variations, 2,699 Pda_hap1-specific and 2,702 Pda_hap2-specific genes, and 1,118 allele-specific expressed genes. Genome composition indicated 4.1% of the P. davidiana genome was non-peach origin, out of which 94.5% was derived from almond. Based on the haplotype genome, the aphid resistance quantitative trait locus (QTL) was mapped at the end of Pda03. From the aphid resistance QTL, PdaWRKY4 was identified as the major dominant gene, with a 9-bp deletion in its promoter of the resistant phenotype. Specifically, PdaWRKY4 regulates aphid resistance by promoting PdaCYP716A1-mediated anti-aphid metabolite betulin biosynthesis. Moreover, we employed a genome design to develop a breeding workflow for rapidly and precisely producing aphid-resistant peaches. In conclusion, this study identifies a novel aphid resistance gene and provides insights into genome design for the development of resistant fruit cultivars.

Eleocharis vivipara, an amphibious sedge in the Cyperaceae family, has several remarkable properties, most notably its alternate use of C₃ photosynthesis underwater and C₄ photosynthesis on land. However, the absence of genomic data has hindered its utility for evolutionary and genetic research. Here, we present a high-quality genome for E. vivipara, representing the first chromosome-level genome for the Eleocharis genus, with an approximate size of 965.22 Mb mainly distributed across 10 chromosomes. Its Hi–C pattern, chromosome clustering results, and one-to-one genome synteny across two subgroups indicates a tetraploid structure with chromosome count 2n = 4x = 20. Phylogenetic analysis suggests that E. vivipara diverged from Cyperus esculentus approximately 32.96 million years ago (Mya), and underwent a whole-genome duplication (WGD) about 3.5 Mya. Numerous fusion and fission events were identified between the chromosomes of E. vivipara and its close relatives. We demonstrate that E. vivipara has holocentromeres, a chromosomal feature which can maintain the stability of such chromosomal rearrangements. Experimental transplantation and cross-section studies showed its terrestrial culms developed C₄ Kranz anatomy with increased number of chloroplasts in the bundle sheath (BS) cells. Gene expression and weighted gene co-expression network analysis (WGCNA) showed overall elevated expression of core genes associated with the C₄ pathway, and significant enrichment of genes related to modified culm anatomy and photosynthesis efficiency. We found evidence of mixed nicotinamide adenine dinucleotide - malic enzyme and phosphoenolpyruvate carboxykinase type C₄ photosynthesis in E. vivipara, and hypothesize that the evolution of C₄ photosynthesis predates the WGD event. The mixed type is dominated by subgenome A and supplemented by subgenome B. Collectively, our findings not only shed light on the evolution of E. vivipara and karyotype within the Cyperaceae family, but also provide valuable insights into the transition between C₃ and C₄ photosynthesis, offering promising avenues for crop improvement and breeding.

Pineapple is the third most crucial tropical fruit worldwide and available in five varieties. Genomes of different pineapple varieties have been released to date; however, none of them are complete, with all exhibiting substantial gaps and representing only two of the five pineapple varieties. This significantly hinders the advancement of pineapple breeding efforts. In this study, we sequenced the genomes of three varieties: a wild pineapple variety, a fiber pineapple variety, and a globally cultivated edible pineapple variety. We constructed the first gap-free reference genome (Ref) for pineapple. By consolidating multiple sources of evidence and manually revising each gene structure annotation, we identified 26,656 protein-coding genes. The BUSCO evaluation indicated a completeness of 99.2%, demonstrating the high quality of the gene structure annotations in this genome. Utilizing these resources, we identified 7,209 structural variations across the three varieties. Approximately 30.8% of pineapple genes were located within ±5 kb of structural variations, including 30 genes associated with anthocyanin synthesis. Further analysis and functional experiments demonstrated that the high expression of AcMYB528 aligns with the accumulation of anthocyanins in the leaves, both of which may be affected by a 1.9-kb insertion fragment. In addition, we developed the Ananas Genome Database, which offers data browsing, retrieval, analysis, and download functions. The construction of this database addresses the lack of pineapple genome resource databases. In summary, we acquired a seamless pineapple reference genome with high-quality gene structure annotations, providing a solid foundation for pineapple genomics and a valuable reference for pineapple breeding.

Domestication has shaped the population structure and agronomic traits of tea plants, yet the complexity of tea population structure and genetic variation that determines these traits remains unclear. We here investigated the resequencing data of 363 diverse tea accessions collected extensively from almost all tea distributions and found that the population structure of tea plants was divided into eight subgroups, which were basically consistent with their geographical distributions. The genetic diversity of tea plants in China decreased from southwest to east as latitude increased. Results also indicated that Camellia sinensis var. assamica (CSA) illustrated divergent selection signatures with Camellia sinensis var. sinensis (CSS). The domesticated genes of CSA were mainly involved in leaf development, flavonoid and alkaloid biosynthesis, while the domesticated genes in CSS mainly participated in amino acid metabolism, aroma compounds biosynthesis, and cold stress. Comparative population genomics further identified ~730 Mb novel sequences, generating 6,058 full-length protein-encoding genes, significantly expanding the gene pool of tea plants. We also discovered 217,376 large-scale structural variations and 56,583 presence and absence variations (PAVs) across diverse tea accessions, some of which were associated with tea quality and stress resistance. Functional experiments demonstrated that two PAV genes (CSS0049975 and CSS0006599) were likely to drive trait diversification in cold tolerance between CSA and CSS tea plants. The overall findings not only revealed the genetic diversity and domestication of tea plants, but also underscored the vital role of structural variations in the diversification of tea plant traits.

The Sapindaceae family, encompassing a wide range of plant forms such as herbs, vines, shrubs, and trees, is widely distributed across tropical and subtropical regions. This family includes economically important crops like litchi, longan, rambutan, and ackee. With the wide application of genomic technologies in recent years, several Sapindaceae plant genomes have been decoded, leading to an accumulation of substantial omics data in this field. This surge in data highlights the pressing need for a unified genomic data center capable of storing, sharing, and analyzing these data. Here, we introduced SapBase, that is, the Sapindaceae Genome Database. SapBase houses seven published plant genomes alongside their corresponding gene structure and functional annotations, small RNA annotations, gene expression profiles, gene pathways, and synteny block information. It offers user-friendly features for gene information mining, co-expression analysis, and inter-species comparative genomic analysis. Furthermore, we showcased SapBase's extensive capacities through a detailed bioinformatic analysis of a MYB gene in litchi. Thus, SapBase could serve as an integrative genomic resource and analysis platform for the scientific exploration of Sapinaceae species and their comparative studies with other plants.

Structural variations (SVs) are a feature of plant genomes that has been largely unexplored despite their significant impact on plant phenotypic traits and local adaptation to abiotic and biotic stress. In this study, we employed woolly grape (Vitis retordii), a species native to the tropical and subtropical regions of East Asia with both coastal and inland habitats, as a valuable model for examining the impact of SVs on local adaptation. We assembled a haplotype-resolved chromosomal reference genome for woolly grape, and conducted population genetic analyses based on whole- genome sequencing (WGS) data from coastal and inland populations. The demographic analyses revealed recent bottlenecks in all populations and asymmetric gene flow from the inland to the coastal population. In total, 1,035 genes associated with plant adaptive regulation for salt stress, radiation, and environmental adaptation were detected underlying local selection by SVs and SNPs in the coastal population, of which 37.29% and 65.26% were detected by SVs and SNPs, respectively. Candidate genes such as FSD2, RGA1, and AAP8 associated with salt tolerance were found to be highly differentiated and selected during the process of local adaptation to coastal habitats in SV regions. Our study highlights the importance of SVs in local adaptation; candidate genes related to salt stress and climatic adaptation to tropical and subtropical environments are important genomic resources for future breeding programs of grapevine and its rootstocks.

Dendrocalamus brandisii (Munro) Kurz is a sympodial bamboo species with inimitable taste and flavorful shoots. Its rapid growth and use as high-quality material make this bamboo species highly valued for both food processing and wood applications. However, genome information for D. brandisii is lacking, primarily due to its polyploidy and large genome size. Here, we assembled a high-quality genome for hexaploid D. brandisii, which comprises 70 chromosomes with a total size of 2,756 Mb, using long-read HiFi sequencing. Furthermore, we accurately separated the genome into its three constituent subgenomes. We used Oxford Nanopore Technologies long reads to construct a transcriptomic dataset covering 15 tissues for gene annotation to complement our genome assembly, revealing differential gene expression and post-transcriptional regulation. By integrating metabolome analysis, we unveiled that well-balanced lignin formation, as well as abundant flavonoid and fructose contents, contribute to the superior quality of D. brandisii shoots. Integrating genomic, transcriptomic, and metabolomic datasets provided a solid foundation for enhancing bamboo shoot quality and developing efficient gene-editing techniques. This study should facilitate research on D. brandisii and enhance its use as a food source and wood material by providing crucial genomic resources.

ChinaMu is the largest sequence-indexed Mutator (Mu) transposon insertional library in maize (Zea mays). In this study, we made significant improvements to the size and quality of the ChinaMu library. We developed a new Mu-tag isolation method Mu-Tn5-seq (MuT-seq). Compared to the previous method used by ChinaMu, MuT-seq recovered 1/3 more germinal insertions, while requiring only about 1/14 of the sequencing volume and 1/5 of the experimental time. Using MuT-seq, we identified 113,879 germinal insertions from 3,168 Mu-active F₁ families. We also assembled a high-quality genome for the Mu-active line Mu-starter, which harbors the initial active MuDR element and was used as the pollen donor for the mutation population. Using the Mu-starter genome, we recovered 33,662 (15.6%) additional germinal insertions in 3,244 (7.4%) genes in the Mu-starter line. The Mu-starter genome also improved the assignment of 117,689 (54.5%) germinal insertions. The newly upgraded ChinaMu dataset currently contains 215,889 high-quality germinal insertions. These insertions cover 32,224 pan-genes in the Mu-starter and B73Ref5 genomes, including 23,006 (80.4%) core genes shared by the two genomes. As a test model, we investigated Mu insertions in the pentatricopeptide repeat (PPR) superfamily, discovering insertions for 92% (449/487) of PPR genes in ChinaMu, demonstrating the usefulness of ChinaMu as a functional genomics resource for maize.

Nypa fruticans (Wurmb), a mangrove palm species with origins dating back to the Late Cretaceous period, is a unique species for investigating long-term adaptation strategies to intertidal environments and the early evolution of palms. Here, we present a chromosome-level genome sequence and assembly for N. fruticans. We integrated the genomes of N. fruticans and other palm family members for a comparative genomic analysis, which confirmed that the common ancestor of all palms experienced a whole-genome duplication event around 89 million years ago, shaping the distinctive characteristics observed in this clade. We also inferred a low mutation rate for the N. fruticans genome, which underwent strong purifying selection and evolved slowly, thus contributing to its stability over a long evolutionary period. Moreover, ancient duplicates were preferentially retained, with critical genes having experienced positive selection, enhancing waterlogging tolerance in N. fruticans. Furthermore, we discovered that the pseudogenization of Early Methionine-labelled 1 (EM1) and EM6 in N. fruticans underly its crypto-vivipary characteristics, reflecting its intertidal adaptation. Our study provides valuable genomic insights into the evolutionary history, genome stability, and adaptive evolution of the mangrove palm. Our results also shed light on the long-term adaptation of this species and contribute to our understanding of the evolutionary dynamics in the palm family.

Brassica napus, commonly known as rapeseed or canola, is a major oil crop contributing over 13% to the stable supply of edible vegetable oil worldwide. Identification and understanding the gene functions in the B. napus genome is crucial for genomic breeding. A group of genes controlling agronomic traits have been successfully cloned through functional genomics studies in B. napus. In this review, we present an overview of the progress made in the functional genomics of B. napus, including the availability of germplasm resources, omics databases and cloned functional genes. Based on the current progress, we also highlight the main challenges and perspectives in this field. The advances in the functional genomics of B. napus contribute to a better understanding of the genetic basis underlying the complex agronomic traits in B. napus and will expedite the breeding of high quality, high resistance and high yield in B. napus varieties.

Rice (Oryza sativa) is a significant crop worldwide with a genome shaped by various evolutionary factors. Rice centromeres are crucial for chromosome segregation, and contain some unreported genes. Due to the diverse and complex centromere region, a comprehensive understanding of rice centromere structure and function at the population level is needed. We constructed a high-quality centromere map based on the rice super pan-genome consisting of a 251-accession panel comprising both cultivated and wild species of Asian and African rice. We showed that rice centromeres have diverse satellite repeat CentO, which vary across chromosomes and subpopulations, reflecting their distinct evolutionary patterns. We also revealed that long terminal repeats (LTRs), especially young Gypsy-type LTRs, are abundant in the peripheral CentO-enriched regions and drive rice centromere expansion and evolution. Furthermore, high-quality genome assembly and complete telomere-to-telomere (T2T) reference genome enable us to obtain more centromeric genome information despite mapping and cloning of centromere genes being challenging. We investigated the association between structural variations and gene expression in the rice centromere. A centromere gene, OsMAB, which positively regulates rice tiller number, was further confirmed by expression quantitative trait loci, haplotype analysis and clustered regularly interspaced palindromic repeats (CRISPR)/CRISPR-associated protein 9 methods. By revealing the new insights into the evolutionary patterns and biological roles of rice centromeres, our finding will facilitate future research on centromere biology and crop improvement.