Genomic insights into the evolution of Echinochloa species as weed and orphan crop

Comparative metabolomics and transcriptomics provide new insights into florpyrauxifen-benzyl resistance in Echinochloa glabrescens

Echinochloa glabrescens Munro ex Hook. f. is a weed of the genus Echinocloa (Echinocloa spp.) that occurs frequently in paddy fields, causing serious harm to rice production. Florpyrauxifen-benzyl (FPB) is a foliar-applied herbicide used to control Echinocloa spp. in paddy fields. However, in recent years, with the widespread use of FPB in rice production, FPB-resistant barnyard grasses have been reported. Here, we identified an FPB-resistant E. glabrescens population with a resistance index (RI) of 10.65 and conducted a comparative analysis using untargeted metabolomics and transcriptomics to investigate the differences between an FPB-resistant E. glabrescens population and a susceptible E. glabrescens population after treatment with the recommended field dose of FPB. Our results showed that the FPB-resistant E. glabrescens had 115 differentially accumulated metabolites (DAMs; 65 up-regulated and 50 down-regulated) and 6397 differentially expressed genes (DEGs; 65 up-regulated and 50 down-regulated) compared to the susceptible E. glabrescens. The analysis of DAMs and DEGs revealed that DAMs were significantly enriched in Glutathione metabolism, Arginine and proline metabolism, and Zeatin biosynthesis pathways, while DEGs were mainly enriched in carbon fixation in photosynthetic organisms, photosynthesis, cyanoamino acid metabolism and glutathione metabolism, etc. The glutathione metabolism pathway was found to be significantly enriched for both DEGs and DAMs. Within this pathway, the metabolites (spermine) and genes (GSTU8, GSTU18, GSTF1) may play a pivotal role in the resistance mechanism of FPB-resistant E. glabrescens. Furthermore, we demonstrated the presence of GST-mediated metabolic resistance in an FPB-resistant E. glabrescens population by using NBD-Cl. Overall, our study provides new insights into the underlying mechanisms of E. glabrescens resistance to FPB through a comparative analysis of untargeted metabolomics and transcriptomics. Additionally, we identified the GST-mediated metabolic resistance in an FPB-resistant E. glabrescens population, and screened for three candidate genes (GSTU8, GSTU18, GSTF1), which has significant implications for improving the weed management efficacy of FPB in rice production and guiding judicious herbicide usage.

Current status of community resources and priorities for weed genomics research

Weeds are attractive models for basic and applied research due to their impacts on agricultural systems and capacity to swiftly adapt in response to anthropogenic selection pressures. Currently, a lack of genomic information precludes research to elucidate the genetic basis of rapid adaptation for important traits like herbicide resistance and stress tolerance and the effect of evolutionary mechanisms on wild populations. The International Weed Genomics Consortium is a collaborative group of scientists focused on developing genomic resources to impact research into sustainable, effective weed control methods and to provide insights about stress tolerance and adaptation to assist crop breeding.

Epigenetic Regulation of CYP72A385-Mediated Metabolic Resistance to Novel Auxin Herbicide Florpyrauxifen-benzyl in Echinochloa crus-galli (L.) P. Beauv

Weed's metabolic resistance to herbicides has undermined the sustainability of herbicides and global food security. Notably, we identified an Echinochloa crus-galli (L.) P. Beauv population (R) that evolved resistance to the never-used florpyrauxifen-benzyl, in which florpyrauxifen-benzyl was metabolized faster than the susceptible E. crus-galli population (S). RNA-seq identified potential metabolism-related genes, EcCYP72A385 and EcCYP85A1, whose expression in yeast exhibited the capacity to degrade florpyrauxifen-benzyl. Region-2 in the EcCYP72A385 promoter showed significant demethylation after florpyrauxifen-benzyl treatment in the R population. DNA methyltransferase inhibitors induce EcCYP72A385 overexpression in the S population and endow it with tolerance to florpyrauxifen-benzyl. Moreover, methyltransferase-like 7A (EcMETTL7A) was overexpressed in the S population and specifically bound to the EcCYP72A385 promoter. Transgenic EcCYP72A385 in Arabidopsis and Oryza sativa L. exhibited resistance to florpyrauxifen-benzyl, whereas EcMETTL7A transgenic plants were sensitive. Overall, EcCYP72A385 is the principal functional gene for conferring resistance to florpyrauxifen-benzyl and is regulated by EcMETTL7A in E. crus-galli.

Weed biology and management in the multi-omics era: Progress and perspectives

Weeds pose a significant threat to crop production, resulting in substantial yield reduction. In addition, they possess robust weedy traits that enable them to survive in extreme environments and evade human control. In recent years, the application of multi-omics biotechnologies has helped to reveal the molecular mechanisms underlying these weedy traits. In this review, we systematically describe diverse applications of multi-omics platforms for characterizing key aspects of weed biology, including the origins of weed species, weed classification, and the underlying genetic and molecular bases of important weedy traits such as crop-weed interactions, adaptability to different environments, photoperiodic flowering responses, and herbicide resistance. In addition, we discuss limitations to the application of multi-omics techniques in weed science, particularly compared with their extensive use in model plants and crops. In this regard, we provide a forward-looking perspective on the future application of multi-omics technologies to weed science research. These powerful tools hold great promise for comprehensively and efficiently unraveling the intricate molecular genetic mechanisms that underlie weedy traits. The resulting advances will facilitate the development of sustainable and highly effective weed management strategies, promoting greener practices in agriculture.

Haplotype-resolved genome of Prunus zhengheensis provides insight into its evolution and low temperature adaptation in apricot

Prunus zhengheensis, an extremely rare population of apricots, originated in warm South-East China and is an excellent material for genetic breeding. However, most apricots and two related species (P. sibirica, P. mandshurica) are found in the cold northern regions in China and the mechanism of their distribution is still unclear. In addition, the classification status of P. zhengheensis is controversial. Thus, we generated a high-quality haplotype-resolved genome for P. zhengheensis, exploring key genetic variations in its adaptation and the causes of phylogenetic incongruence. We found extensive phylogenetic discordances between the nuclear and organelle phylogenies of P. zhengheensis, which could be explained by incomplete lineage sorting. A 242.22-Mb pan-genome of the Armeniaca section was developed with 13 chromosomal genomes. Importantly, we identified a 566-bp insertion in the promoter of the HSFA1d gene in apricot and showed that the activity of the HSFA1d promoter increased under low temperatures. In addition, HSFA1d overexpression in Arabidopsis thaliana indicated that HSFA1d positively regulated plant growth under chilling. Therefore, we hypothesized that the insertion in the promoter of HSFA1d in apricot improved its low-temperature adaptation, allowing it to thrive in relatively cold locations. The findings help explain the weather adaptability of Armeniaca plants.

Syntenic analysis of ACCase loci and target-site-resistance mutations in cyhalofop-butyl resistant Echinochloa crus-galli var. crus-galli in Japan

BACKGROUND

Recently, suspected cyhalofop-butyl-resistant populations of allohexaploid weed Echinochloa crus-galli var. crus-galli were discovered in rice fields in Aichi Prefecture, Japan. Analyzing the target-site ACCase genes of cyhalofop-butyl helps understand the resistance mechanism. However, in E. crus-galli, the presence of multiple ACCase genes and the lack of detailed gene investigations have complicated the analysis of target-site genes. Therefore, in this study, we characterized the herbicide response of E. crus-galli lines and thoroughly characterized the ACCase genes, including the evaluation of gene mutations in the ACCase genes of each line.

RESULT

Four suspected resistant lines collected from Aichi Prefecture showed varying degrees of resistance to cyhalofop-butyl and other FOP-class ACCase inhibitors but were sensitive to herbicides with other modes of action. Through genomic analysis, six ACCase loci were identified in the E. crus-galli genome. We renamed each gene based on its syntenic relationship with other ACCase genes in the Poaceae species. RNA-sequencing analysis revealed that all ACCase genes, except the pseudogenized copy ACCase2A, were transcribed at a similar level in the shoots of E. crus-galli. Mutations known to confer resistance to FOP-class herbicides, that is W1999C, W2027C/S and I2041N, were found in all resistant lines in either ACCase1A, ACCase1B or ACCase2C.

CONCLUSION

In this study, we found that the E. crus-galli lines were resistant exclusively to ACCase-inhibiting herbicides, with a target-site resistance mutation in the ACCase gene. Characterization of ACCase loci in E. crus-galli provides a basis for further research on ACCase herbicide resistance in Echinochloa spp. © 2023 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

A reference genome of Commelinales provides insights into the commelinids evolution and global spread of water hyacinth (Pontederia crassipes)

Commelinales belongs to the commelinids clade, which also comprises Poales that includes the most important monocot species, such as rice, wheat, and maize. No reference genome of Commelinales is currently available. Water hyacinth (Pontederia crassipes or Eichhornia crassipes), a member of Commelinales, is one of the devastating aquatic weeds, although it is also grown as an ornamental and medical plant. Here, we present a chromosome-scale reference genome of the tetraploid water hyacinth with a total length of 1.22 Gb (over 95% of the estimated size) across 8 pseudochromosome pairs. With the representative genomes, we reconstructed a phylogeny of the commelinids, which supported Zingiberales and Commelinales being sister lineages of Arecales and shed lights on the controversial relationship of the orders. We also reconstructed ancestral karyotypes of the commelinids clade and confirmed the ancient commelinids genome having 8 chromosomes but not 5 as previously reported. Gene family analysis revealed contraction of disease-resistance genes during polyploidization of water hyacinth, likely a result of fitness requirement for its role as a weed. Genetic diversity analysis using 9 water hyacinth lines from 3 continents (South America, Asia, and Europe) revealed very closely related nuclear genomes and almost identical chloroplast genomes of the materials, as well as provided clues about the global dispersal of water hyacinth. The genomic resources of P. crassipes reported here contribute a crucial missing link of the commelinids species and offer novel insights into their phylogeny.

Target-Site and Metabolic Resistance Mechanisms to Penoxsulam in Late Watergrass (Echinochloa phyllopogon) in China

Echinochloa phyllopogon, a malignant weed in Northeast China's paddy fields, is currently presenting escalating resistance concerns. Our study centered on the HJHL-715 E. phyllopogon population, which showed heightened resistance to penoxsulam, through a whole-plant bioassay. Pretreatment with a P450 inhibitor malathion significantly increased penoxsulam sensitivity in resistant plants. In order to determine the resistance mechanism of the resistant population, we purified the resistant population from individual plants and isolated target-site resistance (TSR) and nontarget-site resistance (NTSR) materials. Pro-197-Thr and Trp-574-Leu mutations in acetolactate synthase (ALS) 1 and ALS2 of the resistant population drove reduced sensitivity of penoxsulam to the target-site ALS, the primary resistance mechanisms. To fully understand the NTSR mechanism, NTSR materials were investigated by using RNA-sequencing (RNA-seq) combined with a reference genome. High-performance liquid chromatography-tandem mass spectrometry (HPLC-MS/MS) analysis further supported the enhanced penoxsulam metabolism in NTSR materials. Gene expression data and quantitative reverse transcription polymerase chain reaction (qRT-PCR) validation confirmed 29 overexpressed genes under penoxsulam treatment, with 16 genes concurrently upregulated with quinclorac and metamifop treatment. Overall, our study confirmed coexisting TSR and NTSR mechanisms in E. phyllopogon's resistance to ALS inhibitors.

Milletdb: a multi-omics database to accelerate the research of functional genomics and molecular breeding of millets

Millets are a class of nutrient-rich coarse cereals with high resistance to abiotic stress; thus, they guarantee food security for people living in areas with extreme climatic conditions and provide stress-related genetic resources for other crops. However, no platform is available to provide a comprehensive and systematic multi-omics analysis for millets, which seriously hinders the mining of stress-related genes and the molecular breeding of millets. Here, a free, web-accessible, user-friendly millets multi-omics database platform (Milletdb, http://milletdb.novogene.com) has been developed. The Milletdb contains six millets and their one related species genomes, graph-based pan-genomics of pearl millet, and stress-related multi-omics data, which enable Milletdb to be the most complete millets multi-omics database available. We stored GWAS (genome-wide association study) results of 20 yield-related trait data obtained under three environmental conditions [field (no stress), early drought and late drought] for 2 years in the database, allowing users to identify stress-related genes that support yield improvement. Milletdb can simplify the functional genomics analysis of millets by providing users with 20 different tools (e.g., 'Gene mapping', 'Co-expression', 'KEGG/GO Enrichment' analysis, etc.). On the Milletdb platform, a gene PMA1G03779.1 was identified through 'GWAS', which has the potential to modulate yield and respond to different environmental stresses. Using the tools provided by Milletdb, we found that the stress-related PLATZs TFs (transcription factors) family expands in 87.5% of millet accessions and contributes to vegetative growth and abiotic stress responses. Milletdb can effectively serve researchers in the mining of key genes, genome editing and molecular breeding of millets.

Telomere-to-telomere genome assembly of an allotetraploid pernicious weed, Echinochloa phyllopogon

Echinochloa phyllopogon is an allotetraploid pernicious weed species found in rice fields worldwide that often exhibit resistance to multiple herbicides. An accurate genome sequence is essential to comprehensively understand the genetic basis underlying the traits of this species. Here, the telomere-to-telomere genome sequence of E. phyllopogon was presented. Eighteen chromosome sequences spanning 1.0 Gb were constructed using the PacBio highly fidelity long technology. Of the 18 chromosomes, 12 sequences were entirely assembled into telomere-to-telomere and gap-free contigs, whereas the remaining six sequences were constructed at the chromosomal level with only eight gaps. The sequences were assigned to the A and B genome with total lengths of 453 and 520 Mb, respectively. Repetitive sequences occupied 42.93% of the A genome and 48.47% of the B genome, although 32,337, and 30,889 high-confidence genes were predicted in the A and B genomes, respectively. This suggested that genome extensions and gene disruptions caused by repeated sequence accumulation often occur in the B genome before polyploidization to establish a tetraploid genome. The highly accurate and comprehensive genome sequence could be a milestone in understanding the molecular mechanisms of the pernicious traits and in developing effective weed control strategies to avoid yield loss in rice production.

Selection and validation of reliable reference genes for quantitative real-time PCR in Barnyard millet (Echinochloa spp.) under varied abiotic stress conditions

Quantitative real-time polymerase chain reaction (RT-qPCR) using a stable reference gene is widely used for gene expression research. Barnyard millet (Echinochloa spp.) is an ancient crop in Asia and Africa that is widely cultivated for food and fodder. It thrives well under drought, salinity, cold, and heat environmental conditions, besides adapting to any soil type. To date, there are no gene expression studies performed to identify the potential candidate gene responsible for stress response in barnyard millet, due to lack of reference gene. Here, 10 candidate reference genes, Actin (ACT), α-tubulin (α-TUB), β-tubulin (β-TUB), RNA pol II (RP II), elongation factor-1 alpha (EF-1α), adenine phosphoribosyltransferase (APRT), TATA-binding protein-like factor (TLF), ubiquitin-conjugating enzyme 2 (UBC2), ubiquitin-conjugating enzyme E2L5 (UBC5) and glyceraldehyde-3-phosphate dehydrogenase (GAPDH), were selected from mRNA sequences of E. crus-galli and E. colona var frumentacea. Five statistical algorithms (geNorm, NormFinder, BestKeeper, ΔCt, and RefFinder) were applied to determine the expression stabilities of these genes in barnyard millet grown under four different abiotic stress (drought, salinity, cold and heat) exposed at different time points. The UBC5 and ɑ-TUB in drought, GAPDH in salinity, GAPDH and APRT in cold, and EF-1α and RP II in heat were the most stable reference genes, whereas ß-TUB was the least stable irrespective of stress conditions applied. Further Vn/Vn + 1 analysis revealed two reference genes were sufficient to normalize gene expression across all sample sets. The suitability of identified reference genes was validated with Cu-ZnSOD (SOD1) in the plants exposed to different abiotic stress conditions. The results revealed that the relative quantification of the SOD1 gene varied according to reference genes and the number of reference genes used, thus highlighting the importance of the choice of a reference gene in such experiments. This study provides a foundational framework for standardizing RT-qPCR analyses, enabling accurate gene expression profiling in barnyard millet.

Inter-Species Investigation of Biological Traits among Eight Echinochloa Species

Due to the diversity of Echinochloa species and the limited understanding of their damage processes in rice fields, clarifying the biological properties of distinct species could help create a foundation for effective control techniques. Pot experiments and field competition trials were conducted using eight Echinochloa species to elucidate their biological differences and assess their varying levels of negative impact on rice. The survey outcomes showed that E. oryzoides had the highest 1000-grain weight (3.12 g) while E. colona had the lowest (0.90 g). The largest grain number per spikelet found in E. glabrescens (940) was 3.4 times greater than that in E. oryzoides (277). Different species responded variably to changes in temperature and photoperiod. Except for E. caudate, all Echinochloa species exhibited a shortened growth period with the delay of the sowing date. Under field competitive conditions, all Echinochloa species exhibited significantly greater net photosynthetic rates than rice, with E. crusgalli exhibiting the highest photosynthetic capacity. Moreover, in this resource-limited setting, barnyardgrass species had a decrease in tiller formation and panicle initiation but a significant increase in plant height. These findings contribute valuable insights into the biological characteristics of barnyardgrass populations and provide guidance for implementing effective control measures in rice fields.

A syntelog-based pan-genome provides insights into rice domestication and de-domestication

BACKGROUND

Asian rice is one of the world's most widely cultivated crops. Large-scale resequencing analyses have been undertaken to explore the domestication and de-domestication genomic history of Asian rice, but the evolution of rice is still under debate.

RESULTS

Here, we construct a syntelog-based rice pan-genome by integrating and merging 74 high-accuracy genomes based on long-read sequencing, encompassing all ecotypes and taxa of Oryza sativa and Oryza rufipogon. Analyses of syntelog groups illustrate subspecies divergence in gene presence-and-absence and haplotype composition and identify massive genomic regions putatively introgressed from ancient Geng/japonica to ancient Xian/indica or its wild ancestor, including almost all well-known domestication genes and a 4.5-Mbp centromere-spanning block, supporting a single domestication event in main rice subspecies. Genomic comparisons between weedy and cultivated rice highlight the contribution from wild introgression to the emergence of de-domestication syndromes in weedy rice.

CONCLUSIONS

This work highlights the significance of inter-taxa introgression in shaping diversification and divergence in rice evolution and provides an exploratory attempt by utilizing the advantages of pan-genomes in evolutionary studies.

Recurrent Selection of Echinochloa crus-galli with a Herbicide Mixture Reduces Progeny Sensitivity

Herbicide mixtures are used to increase the spectrum of weed control and to manage weeds with target-site resistance to some herbicides. However, the effect of mixtures on the evolution of herbicide resistance caused by enhanced metabolism is unknown. This study evaluated the effect of a fenoxaprop-p-ethyl and imazethapyr mixture on the evolution of herbicide resistance in Echinochloa crus-galli using recurrent selection at sublethal doses. The progeny from second generations selected with the mixture had lower control than parental plants or the unselected progeny. GR₅₀ increased 1.6- and 2.6-fold after two selection cycles with the mixture in susceptible (POP1-S) and imazethapyr-resistant (POP2-IR) biotypes, respectively. There was evidence that recurrent selection with this sublethal mixture had the potential to evolve cross-resistance to diclofop, cyhalofop, sethoxydim, and quinclorac. Mixture selection did not cause increased relative expression for a set of analyzed genes (CYP71AK2, CYP72A122, CYP72A258, CYP81A12, CYP81A14, CYP81A21, CYP81A22, and GST1). Fenoxaprop, rather than imazethapyr, is the main contributor to the decreased control in the progenies after recurrent selection with the mixture in low doses. This is the first study reporting the effect of a herbicide mixture at low doses on herbicide resistance evolution. The lack of control using the mixture may result in decreased herbicide sensitivity of the weed progenies. Using mixtures may select important detoxifying genes that have the potential to metabolize herbicides in patterns that cannot currently be predicted. The use of fully recommended herbicide rates in herbicide mixtures is recommended to reduce the risk of this type of resistance evolution.

Weed genomics: yielding insights into the genetics of weedy traits for crop improvement

Weeds cause tremendous economic and ecological damage worldwide. The number of genomes established for weed species has sharply increased during the recent decade, with some 26 weed species having been sequenced and de novo genomes assembled. These genomes range from 270 Mb (Barbarea vulgaris) to almost 4.4 Gb (Aegilops tauschii). Importantly, chromosome-level assemblies are now available for 17 of these 26 species, and genomic investigations on weed populations have been conducted in at least 12 species. The resulting genomic data have greatly facilitated studies of weed management and biology, especially origin and evolution. Available weed genomes have indeed revealed valuable weed-derived genetic materials for crop improvement. In this review, we summarize the recent progress made in weed genomics and provide a perspective for further exploitation in this emerging field.

A high-quality chromosome-level genome assembly of Pelteobagrus vachelli provides insights into its environmental adaptation and population history

Pelteobagrus vachelli is a freshwater fish with high economic value, but the lack of genome resources has severely restricted its industrial development and population conservation. Here, we constructed the first chromosome-level genome assembly of P. vachelli with a total length of approximately 662.13 Mb and a contig N50 was 14.02 Mb, and scaffolds covering 99.79% of the assembly were anchored to 26 chromosomes. Combining the comparative genome results and transcriptome data under environmental stress (high temperature, hypoxia and Edwardsiella. ictaluri infection), the MAPK signaling pathway, PI3K-Akt signaling pathway and apelin signaling pathway play an important role in environmental adaptation of P. vachelli, and these pathways were interconnected by the ErbB family and involved in cell proliferation, differentiation and apoptosis. Population evolution analysis showed that artificial interventions have affected wild populations of P. vachelli. This study provides a useful genomic information for the genetic breeding of P. vachelli, as well as references for further studies on fish biology and evolution.

Comparative Analysis of Whole Chloroplast Genomes of Three Common Species of Echinochloa (Gramineae) in Paddy Fields

Echinochloa crus-galli var. crus-galli, E. crus-galli var. zelayensis, and E. glabrescens, morphologically similar at the seedling stage, are the most pernicious barnyard grass species in paddy fields worldwide. Chloroplast (cp) genomes could be conducive to their identification. In this study, we assembled the complete cp genome sequences of Echinochloa crus-galli var. crus-galli (139,856 bp), E. crus-galli var. zelayensis (139,874 bp), and E. glabrescens (139,874 bp), which exhibited a typical circular tetramerous structure, large and small single-copy regions, and a pair of inverted repeats. In Echinochloa crus-galli var. crus-galli, there were 136 simple sequence (SSRs) and 62 long (LRs) repeats, and in the other two species, 139 SSRs and 68 LRs. Each cp genome contains 92 protein-encoding genes. In Echinochloa crus-galli var. crus-galli and E. glabrescens, 321 and 1 single-nucleotide polymorphisms were detected compared to Echinochloa crus-galli var. zelayensis. IR expansion and contraction revealed small differences between the three species. The phylogenetic tree based on cp genomes demonstrated the phylogenetic relationship between ten barnyard grass species and other common Gramineae plants, showing new genetic relationships of the genus Echinochloa. This study provides valuable information on cp genomes, useful for identifying and classifying the genus Echinochloa and studying its phylogenetic relationships and evolution.

Lateral transfers lead to the birth of momilactone biosynthetic gene clusters in grass

Momilactone A, an important plant labdane-related diterpenoid, functions as a phytoalexin against pathogens and an allelochemical against neighboring plants. The genes involved in the biosynthesis of momilactone A are found in clusters, i.e., momilactone A biosynthetic gene clusters (MABGCs), in the rice and barnyardgrass genomes. In addition, we know little about the origin and evolution of MABGCs. Here, we integrated results from comprehensive phylogeny and comparative genomic analyses of the core genes of MABGC-like clusters and MABGCs in 40 monocot plant genomes, providing convincing evidence for the birth and evolution of MABGCs in grass species. The MABGCs found in the PACMAD clade of the core grass lineage (including Panicoideae and Chloridoideae) originated from a MABGC-like cluster in Triticeae (BOP clade) via lateral gene transfer (LGT) and followed by recruitment of MAS1/2 and CYP76L1 genes. The MABGCs in Oryzoideae originated from PACMAD through another LGT event and lost CYP76L1 afterwards. The Oryza MABGC and another Oryza diterpenoid cluster c2BGC are two distinct clusters, with the latter originating from gene duplication and relocation within Oryzoideae. Further comparison of the expression patterns of the MABGC genes between rice and barnyardgrass in response to pathogen infection and allelopathy provides novel insights into the functional innovation of MABGCs in plants. Our results demonstrate LGT-mediated origination of MABGCs in grass and shed lights into the evolutionary innovation and optimization of plant biosynthetic pathways.

Horizontal transfer and evolution of the biosynthetic gene cluster for benzoxazinoids in plants

Benzoxazinoids are a class of protective and allelopathic plant secondary metabolites that have been identified in multiple grass species and are encoded by the Bx biosynthetic gene cluster (BGC) in maize. Data mining of 41 high-quality grass genomes identified complete Bx clusters (containing genes Bx1-Bx5 and Bx8) in three genera (Zea, Echinochloa, and Dichanthelium) of Panicoideae and partial clusters in Triticeae. The Bx cluster probably originated from gene duplication and chromosomal translocation of native homologs of Bx genes. An ancient Bx cluster that included additional Bx genes (e.g., Bx6) is presumed to have been present in ancestral Panicoideae. The ancient Bx cluster was putatively gained by the Triticeae ancestor via horizontal transfer (HT) from the ancestral Panicoideae and later separated into multiple segments on different chromosomes. Bx6 appears to have been under less constrained selection compared with the Bx cluster during the evolution of Panicoideae, as evidenced by the fact that it was translocated away from the Bx cluster in Zea mays, moved to other chromosomes in Echinochloa, and even lost in Dichanthelium. Further investigations indicate that purifying selection and polyploidization have shaped the evolutionary trajectory of Bx clusters in the grass family. This study provides the first candidate case of HT of a BGC between plants and sheds new light on the evolution of BGCs.