The banana (Musa acuminata) genome and the evolution of monocotyledonous plants

Telomere-to-telomere genome assembly of Eleocharis dulcis and expression profiles during corm development

Eleocharis dulcis (Burm. f.) Trin. ex Hensch., commonly known as Chinese water chestnut, is a traditional aquatic vegetable in China, and now is widely cultivated throughout the world because of its high nutritional value and unique tastes. Here, we report the assembly of a 493.24 Mb telomere-to-telomere (T2T) genome of E. dulcis accomplished by integrating ONT ultra-long reads, PacBio long reads and Hi-C data. The reference genome was anchored onto 111 gap-free chromosomes, containing 48.31% repeat elements and 33,493 predicted protein-coding genes. Whole genome duplication (WGD) and inter-genomic synteny analyses indicated that chromosome breakage and genome duplication in E. dulcis possibly occurred multiple times during genome evolution after its divergence from a common ancestor with Rhynchospora breviuscula at ca. 35.6 Mya. A comparative time-course transcriptome analysis of corm development revealed different patterns of gene expression between cultivated and wild accessions with the highest number of differentially expressed genes (DEGs, 15,870) at the middle swelling stage and some of the DEGs were significantly enriched for starch metabolic process.

Celine, a long interspersed nuclear element retrotransposon, colonizes in the centromeres of poplar chromosomes

Centromeres in most multicellular eukaryotes are composed of long arrays of repetitive DNA sequences. Interestingly, several transposable elements, including the well-known long terminal repeat centromeric retrotransposon of maize (CRM), were found to be enriched in functional centromeres marked by the centromeric histone H3 (CENH3). Here, we report a centromeric long interspersed nuclear element (LINE), Celine, in Populus species. Celine has colonized preferentially in the CENH3-associated chromatin of every poplar chromosome, with 84% of the Celine elements localized in the CENH3-binding domains. In contrast, only 51% of the CRM elements were bound to CENH3 domains in Populus trichocarpa. These results suggest different centromere targeting mechanisms employed by Celine and CRM elements. Nevertheless, the high target specificity seems to be detrimental to further amplification of the Celine elements, leading to a shorter life span and patchy distribution among plant species compared with the CRM elements. Using a phylogenetically guided approach, we were able to identify Celine-like LINE elements in tea plant (Camellia sinensis) and green ash tree (Fraxinus pennsylvanica). The centromeric localization of these Celine-like LINEs was confirmed in both species. We demonstrate that the centromere targeting property of Celine-like LINEs is of primitive origin and has been conserved among distantly related plant species.

Molecular Characteristics and Functional Identification of a Key Alpha-Amylase-Encoding Gene AMY11 in Musa acuminata

Alpha-amylase (AMY) plays a significant role in regulating the growth, development, and postharvest quality formation in plants. Nevertheless, little is known about the genome-wide features, expression patterns, subcellular localization, and functional regulation of AMY genes (MaAMYs) in the common starchy banana (Musa acuminata). Twelve MaAMY proteins from the banana genome database were clustered into two groups and contained a conserved catalytic domain. These MaAMYs formed collinear pairs with the AMYs of maize and rice. Three tandem gene pairs were found within the MaAMYs and are indicative of putative gene duplication events. Cis-acting elements of the MaAMY promoters were found to be involved in phytohormone, development, and stress responses. Furthermore, MaAMY02, 08, 09, and 11 were actively expressed during fruit development and ripening. Specifically, MaAMY11 showed the highest expression level at the middle and later stages of banana ripening. Subcellular localization showed that MaAMY02 and 11 were predominately found in the chloroplast, whereas MaAMY08 and 09 were primarily localized in the cytoplasm. Notably, transient attenuation of MaAMY11 expression resulted in an obvious increase in the starch content of banana fruit, while a significant decrease in starch content was confirmed through the transient overexpression of MaAMY11. Together, these results reveal new insights into the structure, evolution, and expression patterns of the MaAMY family, affirming the functional role of MaAMY11 in the starch degradation of banana fruit.

Phylogenomics identifies parents of naturally occurring tetraploid bananas

BACKGROUND

Triploid bananas are almost sterile. However, we succeeded in harvesting seeds from two edible triploid banana individuals (Genotype: ABB) in our conservation repository where various wild diploid bananas were also grown. The resulting rare offspring survived to seedling stages. DNA content analyses reveal that they are tetraploid. Since bananas contain maternally inherited plastids and paternally inherited mitochondria, we sequenced and assembled plastomes and mitogenomes of these seedlings to trace their hybridization history.

RESULTS

The coding sequences of both organellar genomic scaffolds were extracted, aligned, and concatenated for constructing phylogenetic trees. Our results suggest that these tetraploid seedlings be derived from hybridization between edible triploid bananas and wild diploid Musa balbisiana (BB) individuals. We propose that generating female triploid gametes via apomeiosis may allow the triploid maternal bananas to produce viable seeds.

CONCLUSIONS

Our study suggests a practical avenue towards expanding genetic recombination and increasing genetic diversity of banana breeding programs. Further cellular studies are needed to understand the fusion and developmental processes that lead to formation of hybrid embryos in banana reproduction, polyploidization, and evolution.

Striking variation in chromosome structure within Musa acuminata subspecies, diploid cultivars, and F1 diploid hybrids

The majority of cultivated bananas originated from inter- and intra(sub)specific crosses between two wild diploid species, Musa acuminata and Musa balbisiana. Hybridization and polyploidization events during the evolution of bananas led to the formation of clonally propagated cultivars characterized by a high level of genome heterozygosity and reduced fertility. The combination of low fertility in edible clones and differences in the chromosome structure among M. acuminata subspecies greatly hampers the breeding of improved banana cultivars. Using comparative oligo-painting, we investigated large chromosomal rearrangements in a set of wild M. acuminata subspecies and cultivars that originated from natural and human-made crosses. Additionally, we analyzed the chromosome structure of F1 progeny that resulted from crosses between Mchare bananas and the wild M. acuminata 'Calcutta 4' genotype. Analysis of chromosome structure within M. acuminata revealed the presence of a large number of chromosomal rearrangements showing a correlation with banana speciation. Chromosome painting of F1 hybrids was complemented by Illumina resequencing to identify the contribution of parental subgenomes to the diploid hybrid clones. The balanced presence of both parental genomes was revealed in all F1 hybrids, with the exception of one clone, which contained only Mchare-specific SNPs and thus most probably originated from an unreduced diploid gamete of Mchare.

Genome-wide identification and characterization of DTX family genes highlighting their locations, functions, and regulatory factors in banana (Musa acuminata)

The detoxification efflux carriers (DTX) are a significant group of multidrug efflux transporter family members that play diverse functions in all kingdoms of living organisms. However, genome-wide identification and characterization of DTX family transporters have not yet been performed in banana, despite its importance as an economic fruit plant. Therefore, a detailed genome-wide analysis of DTX family transporters in banana (Musa acuminata) was conducted using integrated bioinformatics and systems biology approaches. In this study, a total of 37 DTX transporters were identified in the banana genome and divided into four groups (I, II, III, and IV) based on phylogenetic analysis. The gene structures, as well as their proteins' domains and motifs, were found to be significantly conserved. Gene ontology (GO) annotation revealed that the predicted DTX genes might play a vital role in protecting cells and membrane-bound organelles through detoxification mechanisms and the removal of drug molecules from banana cells. Gene regulatory analyses identified key transcription factors (TFs), cis-acting elements, and post-transcriptional regulators (miRNAs) of DTX genes, suggesting their potential roles in banana. Furthermore, the changes in gene expression levels due to pathogenic infections and non-living factor indicate that banana DTX genes play a role in responses to both biotic and abiotic stresses. The results of this study could serve as valuable tools to improve banana quality by protecting them from a range of environmental stresses.

Evolutionary patterns of variations in chromosome counts and genome sizes show positive correlations with taxonomic diversity in tropical gingers

PREMISE

Cytogenetic traits such as an organism's chromosome number and genome size are taxonomically critical as they are instrumental in defining angiosperm diversity. Variations in these traits can be traced to evolutionary processes such as polyploidization, although geographic variations across cytogenetic traits remain underexplored. In the pantropical monocot family Zingiberaceae (~1500 species), cytogenetic traits have been well documented; however, the role of these traits in shaping taxonomic diversity and biogeographic patterns of gingers is not known.

METHODS

A time-calibrated Bayesian phylogenetic tree was constructed for 290 taxa covering three of the four subfamilies in Zingiberaceae. We tested models of chromosome number and genome size evolution within the family and whether lineage age, taxonomic diversity, and distributional range explain the variations in the cytogenetic traits. Tests were carried out at two taxonomic ranks: within Zingiberaceae and within genus Hedychium using correlations, generalized linear models and phylogenetic least square models.

RESULTS

The most frequent changes in chromosome number within Zingiberaceae were noted to be demi-polyploidization and polyploidization (~57% of the time), followed by ascending dysploidy (~27%). The subfamily Zingiberoideae showed descending dysploidy at its base, while Alpinioideae showed polyploidization at its internal nodes. Although chromosome counts and genome sizes did not corroborate with each other, suggesting that they are not equivalent; higher chromosome number variations and higher genome size variations were associated with higher taxonomic diversity and wider biogeographic distribution.

CONCLUSIONS

Within Zingiberaceae, multiple incidences of polyploidization were discovered, and cytogenetic events appear to have reduced the genome sizes and increased taxonomic diversity, distributional ranges and invasiveness.

Chromosome-level genome assemblies of Musa ornata and Musa velutina provide insights into pericarp dehiscence and anthocyanin biosynthesis in banana

Musa ornata and Musa velutina are members of the Musaceae family and are indigenous to the South and Southeast Asia. They are very popular in the horticultural market, but the lack of genomic sequencing data and genetic studies has hampered efforts to improve their ornamental value. In this study, we generated the first chromosome-level genome assemblies for both species by utilizing Oxford Nanopore long reads and Hi-C reads. The genomes of M. ornata and M. velutina were assembled into 11 pseudochromosomes with genome sizes of 427.85 Mb and 478.10 Mb, respectively. Repetitive sequences comprised 46.70% and 50.91% of the total genomes for M. ornata and M. velutina, respectively. Differentially expressed gene (DEG) and Gene Ontology (GO) enrichment analyses indicated that upregulated genes in the mature pericarps of M. velutina were mainly associated with the saccharide metabolic processes, particularly at the cell wall and extracellular region. Furthermore, we identified polygalacturonase (PG) genes that exhibited higher expression level in mature pericarps of M. velutina compared to other tissues, potentially being accountable for pericarp dehiscence. This study also identified genes associated with anthocyanin biosynthesis pathway. Taken together, the chromosomal-level genome assemblies of M. ornata and M. velutina provide valuable insights into the mechanism of pericarp dehiscence and anthocyanin biosynthesis in banana, which will significantly contribute to future genetic and molecular breeding efforts.

Genome organization and botanical diversity

The rich diversity of angiosperms, both the planet's dominant flora and the cornerstone of agriculture, is integrally intertwined with a distinctive evolutionary history. Here, we explore the interplay between angiosperm genome organization and botanical diversity, empowered by genomic approaches ranging from genetic linkage mapping to analysis of gene regulation. Commonality in the genetic hardware of plants has enabled robust comparative genomics that has provided a broad picture of angiosperm evolution and implicated both general processes and specific elements in contributing to botanical diversity. We argue that the hardware of plant genomes-both in content and in dynamics-has been shaped by selection for rather substantial differences in gene regulation between plants and animals such as maize and human, organisms of comparable genome size and gene number. Their distinctive genome content and dynamics may reflect in part the indeterminate development of plants that puts strikingly different demands on gene regulation than in animals. Repeated polyploidization of plant genomes and multiplication of individual genes together with extensive rearrangement and differential retention provide rich raw material for selection of morphological and/or physiological variations conferring fitness in specific niches, whether natural or artificial. These findings exemplify the burgeoning information available to employ in increasing knowledge of plant biology and in modifying selected plants to better meet human needs.

Ancient hybridization and repetitive element proliferation in the evolutionary history of the monocot genus Amomum (Zingiberaceae)

Genome size variation is a crucial aspect of plant evolution, influenced by a complex interplay of factors. Repetitive elements, which are fundamental components of genomic architecture, often play a role in genome expansion by selectively amplifying specific repeat motifs. This study focuses on Amomum, a genus in the ginger family (Zingiberaceae), known for its 4.4-fold variation in genome size. Using a robust methodology involving PhyloNet reconstruction, RepeatExplorer clustering, and repeat similarity-based phylogenetic network construction, we investigated the repeatome composition, analyzed repeat dynamics, and identified potential hybridization events within the genus. Our analysis confirmed the presence of four major infrageneric clades (A-D) within Amomum, with clades A-C exclusively comprising diploid species (2n = 48) and clade D encompassing both diploid and tetraploid species (2n = 48 and 96). We observed an increase in the repeat content within the genus, ranging from 84% to 89%, compared to outgroup species with 75% of the repeatome. The SIRE lineage of the Ty1-Copia repeat superfamily was prevalent in most analyzed ingroup genomes. We identified significant difference in repeatome structure between the basal Amomum clades (A, B, C) and the most diverged clade D. Our investigation revealed evidence of ancient hybridization events within Amomum, coinciding with a substantial proliferation of multiple repeat groups. This finding supports the hypothesis that ancient hybridization is a driving force in the genomic evolution of Amomum. Furthermore, we contextualize our findings within the broader context of genome size variations and repeatome dynamics observed across major monocot lineages. This study enhances our understanding of evolutionary processes within monocots by highlighting the crucial roles of repetitive elements in shaping genome size and suggesting the mechanisms that drive these changes.

Characterizing subgenome recombination and chromosomal imbalances in banana varietal lineages

BACKGROUND

Bananas and plantains (Musa spp.) are among the most important crops worldwide. The cultivated varieties are vegetatively propagated, so their genetic diversity is essentially fixed over time. Musa acuminata, M. balbisiana and M. schizocarpa have provided the named A, B and S subgenomes that predominantly constitute these varieties. Here we aimed to characterize intergenetic recombination and chromosomal imbalances between these A/B/S subgenomes, which often result in copy-number variants (CNVs) leading to changes in gene dosage and phenotype, in a diverse panel of bananas and plantains. This will allow us to characterize varietal lineages better and identify sources of genetic variation.

METHODS

We delimited population structure and clonal lineages in a diverse panel of 188 banana and plantain accessions from the most common cultivars using admixture, principal component and phylogenetic analyses. We used new scalable alignment-based methods, Relative Averaged Alignment (RAA) and Relative Coverage, to infer subgenome composition (AA, AAB, etc.) and interspecific recombination.

RESULTS

In our panel, we identified ten varietal lineages composed of somatic clones, plus three groups of tetraploid accessions. We identified chromosomal exchanges resulting in gains/losses in chromosomal segments (CNVs), particularly in AAB and ABB varieties.

CONCLUSIONS

We demonstrated alignment-based RAA and Relative Coverage can identify subgenome composition and introgressions with similar results to more complex approaches based on single nucleotide polymorphism (SNP) databases. These ab initio species-agnostic methods can be used without sequencing a panel of wild ancestors to find private SNPs, or in recently diverged pools where private SNPs are uncommon. The extensive A/B/S exchanges and the variation in the length of some introgressions between lineages further support multiple foundational events of hybridization and residual backcrossing. Imbalances between A/B/S may have resulted in CNVs and gene dosage variation. Since most edible banana genomes are fixed on time, these CNVs are stable genetic variations probably associated with phenotypic variation for future genetic studies.

Ancient Duplication and Lineage-Specific Transposition Determine Evolutionary Trajectory of ERF Subfamily across Angiosperms

AP2/ERF transcription factor family plays an important role in plant development and stress responses. Previous studies have shed light on the evolutionary trajectory of the AP2 and DREB subfamilies. However, knowledge about the evolutionary history of the ERF subfamily in angiosperms still remains limited. In this study, we performed a comprehensive analysis of the ERF subfamily from 107 representative angiosperm species by combining phylogenomic and synteny network approaches. We observed that the expansion of the ERF subfamily was driven not only by whole-genome duplication (WGD) but also by tandem duplication (TD) and transposition duplication events. We also found multiple transposition events in Poaceae, Brassicaceae, Poales, Brassicales, and Commelinids. These events may have had notable impacts on copy number variation and subsequent functional divergence of the ERF subfamily. Moreover, we observed a number of ancient tandem duplications occurred in the ERF subfamily across angiosperms, e.g., in Subgroup IX, IXb originated from ancient tandem duplication events within IXa. These findings together provide novel insights into the evolution of this important transcription factor family.

Genomic analysis of Nypa fruticans elucidates its intertidal adaptations and early palm evolution

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.

Natural Extracts for Antibacterial Applications

Bacteria-induced epidemics and infectious diseases are seriously threatening the health of people around the world. In addition, antibiotic therapy has been inducing increasingly more serious bacterial resistance, which makes it urgent to develop new treatment strategies to combat bacteria, including multidrug-resistant bacteria. Natural extracts displaying antibacterial activity and good biocompatibility have attracted much attention due to greater concerns about the safety of synthetic chemicals and emerging drug resistance. These antibacterial components can be isolated and utilized as antimicrobials, as well as transformed, combined, or wrapped with other substances by using modern assistive technologies to fight bacteria synergistically. This review summarizes recent advances in natural extracts from three kinds of sources-plants, animals, and microorganisms-for antibacterial applications. This work discusses the corresponding antibacterial mechanisms and the future development of natural extracts in antibacterial fields.

Genome-wide identification of XTH gene family in Musa acuminata and response analyses of MaXTHs and xyloglucan to low temperature

Banana (Musa spp.) production is seriously threatened by low temperature (LT) in tropical and subtropical regions. Xyloglucan endotransglycosylase/hydrolases (XTHs) are considered chief enzymes in cell wall remodelling and play a central role in stress responses. However, whether MaXTHs are involved in the low temperature stress tolerance in banana is not clear. Here, the identification and characterization of MaXTHs were carried out, followed by prediction of their cis-acting elements and protein-protein interactions. In addition, candidate MaXTHs involved in banana tolerance to LT were screened through a comparison of their responses to LT between tolerant and sensitive cultivars using RNA-Seq analysis. Moreover, immunofluorescence (IF) labelling was employed to compare changes in the temporal and spatial distribution of different types of xyloglucan components between these two cultivars upon stress. In total, 53 MaXTHs have been identified, and all were predicted to be located in the cell wall, 14 of them also in the cytoplasm. Only 11 MaXTHs have been found to interact with other proteins. Among 16 MaXTHs with LT responsiveness elements, MaXTH26/29/32/35/50 (Group I/II members) and MaXTH7/8 (Group IIIB members) might be involved in banana tolerance to LT stress. IF results suggested that the content of xyloglucan components recognized by CCRC-M87/103/104/106 antibodies might be negatively related to banana chilling tolerance. In conclusion, we have identified the MaXTH gene family and assessed cell wall re-modelling under LT stress. These results will be beneficial for banana breeding against stresses and enrich the cell wall-mediated resistance mechanism in plants to stresses.

Two haplotype-resolved genome assemblies for AAB allotriploid bananas provide insights into banana subgenome asymmetric evolution and Fusarium wilt control

Bananas (Musa spp.) are one of the world's most important fruit crops and play a vital role in food security for many developing countries. Most banana cultivars are triploids derived from inter- and intraspecific hybridizations between the wild diploid ancestor species Musa acuminate (AA) and M. balbisiana (BB). We report two haplotype-resolved genome assemblies of the representative AAB-cultivated types, Plantain and Silk, and precisely characterize ancestral contributions by examining ancestry mosaics across the genome. Widespread asymmetric evolution is observed in their subgenomes, which can be linked to frequent homologous exchange events. We reveal the genetic makeup of triploid banana cultivars and verify that subgenome B is a rich source of disease resistance genes. Only 58.5% and 59.4% of Plantain and Silk genes, respectively, are present in all three haplotypes, with >50% of genes being differentially expressed alleles in different subgenomes. We observed that the number of upregulated genes in Plantain is significantly higher than that in Silk at one-week post-inoculation with Fusarium wilt tropical race 4 (Foc TR4), which confirms that Plantain can initiate defense responses faster than Silk. Additionally, we compared genomic and transcriptomic differences among the genes related to carotenoid synthesis and starch metabolism between Plantain and Silk. Our study provides resources for better understanding the genomic architecture of cultivated bananas and has important implications for Musa genetics and breeding.

Genome-Wide Identification and Expression Analysis of DWARF53 Gene in Response to GA and SL Related to Plant Height in Banana

Dwarfing is one of the common phenotypic variations in asexually reproduced progeny of banana, and dwarfed banana is not only windproof and anti-fallout but also effective in increasing acreage yield. As a key gene in the strigolactone signalling pathway, DWARF53 (D53) plays an important role in the regulation of the height of plants. In order to gain insight into the function of the banana D53 gene, this study conducted genome-wide identification of banana D53 gene based on the M. acuminata, M. balbisiana and M. itinerans genome database. Analysis of MaD53 gene expression under high temperature, low temperature and osmotic stress based on transcriptome data and RT-qPCR was used to analyse MaD53 gene expression in different tissues as well as in different concentrations of GA and SL treatments. In this study, we identified three MaD53, three MbD53 and two MiD53 genes in banana. Phylogenetic tree analysis showed that D53 Musa are equally related to D53 Asparagales and Poales. Both high and low-temperature stresses substantially reduced the expression of the MaD53 gene, but osmotic stress treatments had less effect on the expression of the MaD53 gene. GR24 treatment did not significantly promote the height of the banana, but the expression of the MaD53 gene was significantly reduced in roots and leaves. GA treatment at 100 mg/L significantly promoted the expression of the MaD53 gene in roots, but the expression of this gene was significantly reduced in leaves. In this study, we concluded that MaD53 responds to GA and SL treatments, but "Yinniaijiao" dwarf banana may not be sensitive to GA and SL.

Dataset from transcriptome profiling of Musa resistant and susceptible cultivars in response to Fusarium oxysporum f.sp. cubense race1 and TR4 challenges using Illumina NovaSeq

In this investigation, the study focused on the RNAseq data generated in response to Fusarium oxysporum f.sp. cubense (Foc) race1 (Cavendish infecting strain VCG 0124), targeting both resistant (cv. Rose, AA) and susceptible cultivars (Namarai, AA), and Tropical Race 4 (TR4, strain VCG 01213/16), involving resistant (cv. Rose, AA) and susceptible cultivars (Matti, AA). The respective contrasting cultivars were independently challenged with Foc race1 and TR4, and the root and corm samples were collected in two replications at varying time intervals [0th (control), 2nd, 4th, 6th, and 8th days] in duplicates. The RNA samples underwent stringent quality checks, with all 80 samples meeting the primary parameters, including a satisfactory RNA integrity number (>7). Subsequent library preparation and secondary quality control steps were executed successfully for all samples, paving the way for the sequencing phase. Sequencing generated an extensive amount of data, yielding a range of 10 to 31 million paired-end raw reads per sample, resulting in a cumulative raw data size of 11-50 GB. These raw reads were aligned against the reference genome of Musa acuminata ssp. malaccensis version 2 (DH Pahang), as well as the pathogen genomes of Foc race 1 and Foc TR4, using the HISAT2 alignment tool. The focal point of this study was the investigation of differential gene expression patterns of Musa spp. upon Foc infection. In Foc race1 resistant and susceptible root samples across the designated day intervals, a significant number of genes displayed up-regulation (ranging from 1 to 228) and down-regulation (ranging from 1 to 274). In corm samples, the up-regulated genes ranged from 1 to 149, while down-regulated genes spanned from 3 to 845. For Foc TR4 resistant and susceptible root samples, the expression profiles exhibited a notable up-regulation of genes (ranging from 31 to 964), along with a down-regulation range of 316-1315. In corm samples, up-regulated genes ranged from 57 to 929, while down-regulated genes were observed in the range of 40-936. In addition to the primary analysis, a comprehensive secondary analysis was conducted, including Gene Ontology (GO), euKaryotic Orthologous Groups (KOG) classification, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, and investigations into Simple Sequence Repeats (SSRs), Single Nucleotide Polymorphisms (SNPs), and microRNA (miRNA). The complete dataset was carefully curated and housed at ICAR-NRCB, Trichy, ensuring its accuracy and accessibility for the duration of the study. Further, the raw transcriptome read datasets have been successfully submitted to the National Center for Biotechnology Information - Sequence Read Archive (NCBI-SRA) database, ensuring the accessibility and reproducibility of this valuable dataset for further research endeavors.

High-quality genome assemblies for two Australimusa bananas (Musa spp.) and insights into regulatory mechanisms of superior fiber properties

Bananas (Musa spp.) are monocotyledonous plants with high genetic diversity in the Musaceae family that are cultivated mainly in tropical and subtropical countries. The fruits are a popular food, and the plants themselves have diverse uses. Four genetic groups (genomes) are thought to have contributed to current banana cultivars: Musa acuminata (A genome), Musa balbisiana (B genome), Musa schizocarpa (S genome), and species of the Australimusa section (T genome). However, the T genome has not been effectively explored. Here, we present the high-quality TT genomes of two representative accessions, Abaca (Musa textilis), with high-quality natural fiber, and Utafun (Musa troglodytarum, Fe'i group), with abundant β-carotene. Both the Abaca and Utafun assemblies comprise 10 pseudochromosomes, and their total genome sizes are 613 Mb and 619 Mb, respectively. Comparative genome analysis revealed that the larger size of the T genome is likely attributable to rapid expansion and slow removal of transposons. Compared with those of Musa AA or BB accessions or sisal (Agava sisalana), Abaca fibers exhibit superior mechanical properties, mainly because of their thicker cell walls with a higher content of cellulose, lignin, and hemicellulose. Expression of MusaCesA cellulose synthesis genes peaks earlier in Abaca than in AA or BB accessions during plant development, potentially leading to earlier cellulose accumulation during secondary cell wall formation. The Abaca-specific expressed gene MusaMYB26, which is directly regulated by MusaMYB61, may be an important regulator that promotes precocious expression of secondary cell wall MusaCesAs. Furthermore, MusaWRKY2 and MusaNAC68, which appear to be involved in regulating expression of MusaLAC and MusaCAD, may at least partially explain the high accumulation of lignin in Abaca. This work contributes to a better understanding of banana domestication and the diverse genetic resources in the Musaceae family, thus providing resources for Musa genetic improvement.

Origin and evolution of the triploid cultivated banana genome

Most fresh bananas belong to the Cavendish and Gros Michel subgroups. Here, we report chromosome-scale genome assemblies of Cavendish (1.48 Gb) and Gros Michel (1.33 Gb), defining three subgenomes, Ban, Dh and Ze, with Musa acuminata ssp. banksii, malaccensis and zebrina as their major ancestral contributors, respectively. The insertion of repeat sequences in the Fusarium oxysporum f. sp. cubense (Foc) tropical race 4 RGA2 (resistance gene analog 2) promoter was identified in most diploid and triploid bananas. We found that the receptor-like protein (RLP) locus, including Foc race 1-resistant genes, is absent in the Gros Michel Ze subgenome. We identified two NAP (NAC-like, activated by apetala3/pistillata) transcription factor homologs specifically and highly expressed in fruit that directly bind to the promoters of many fruit ripening genes and may be key regulators of fruit ripening. Our genome data should facilitate the breeding and super-domestication of bananas.

The chromosome-level genome assembly of Fraxinus americana provides insights into the evolution of Oleaceae plants

White ash (Fraxinus americana linn.) originates from the southeastern United States. It is a tall and fast-growing tree species with strong salt-alkali resistance and cold tolerance, making it an important reforestation species and widely planted worldwide. Here, we completed the chromosome-level reference genome assembly of F. americana based on Illumina, PacBio, and Hi-C reads, with a genome size of 878.98 Mb, an N50 of 3.27 Mb, and a heterozygosity rate of 0.3 %. Based on de novo prediction, transcriptome prediction, and homology-based protein prediction, we obtained 39,538 genes. Approximately 843.21 Mb of the assembly genome was composed of 37,928 annotated protein-coding genes, with a gene function annotation rate of 95.93 %. 99.94 % of the overlap clusters (877.44 Mb) were anchored to 23 chromosomes. Synteny analysis of F. americana and other Oleaceae plants showed that F. americana underwent frequent chromosome rearrangements. The amplification of the Ale transposons effectively promoted the genome size of F. americana. Compared with other Oleaceae plants, the Glutathione S-transferase (GST) gene family in the F. americana genome has undergone significant expansion, which may help F. americana cope with adverse natural environments. Furthermore, we found that key enzyme-coding gene families related to lignin biosynthesis were expanded and highly expressed in F. americana leaves. These key genes drive lignin synthesis and benefit F. americana in fast-growing, as well as resisting biotic and abiotic stress. Overall, the F. americana genome assembly provides insights into the evolution of Oleaceae plants and provides abundant resources for breeding and germplasm conservation of white ash.

Transcriptomic and metabolomic differences between banana varieties which are resistant or susceptible to Fusarium wilt

Background

Fusarium wilt, caused by Fusarium oxysporum f. sp. cubense race 4 (Foc4), is the most lethal disease of bananas in Asia.

Methods

To better understand the defense response of banana to Fusarium wilt, the transcriptome and metabolome profiles of the roots from resistant and susceptible bananas inoculated with Foc4 were compared.

Results

After Foc4 inoculation, there were 172 and 1,856 differentially expressed genes (DEGs) in the Foc4-susceptible variety (G1) and Foc4-resistant variety (G9), respectively. In addition, a total of 800 DEGs were identified between G1 and G9, which were mainly involved in the oxidation-reduction process, cell wall organization, phenylpropanoid biosynthesis, and lipid and nitrogen metabolism, especially the DEGs of Macma4_08_g22610, Macma4_11_g19760, and Macma4_03_g06480, encoding non-classical arabinogalactan protein; GDSL-like lipase; and peroxidase. In our study, G9 showed a stronger and earlier response to Foc4 than G1. As the results of metabolomics, lipids, phenylpropanoids and polyketides, organic acids, and derivatives played an important function in response to Fusarium wilt. More importantly, Macma4_11_g19760 might be one of the key genes that gave G9 more resistance to Foc4 by a lowered expression and negative regulation of lipid metabolism. This study illustrated the difference between the transcriptomic and metabolomic profiles of resistant and susceptible bananas. These results improved the current understanding of host-pathogen interactions and will contribute to the breeding of resistant banana plants.

Genome-wide expansion and reorganization during grass evolution: from 30 Mb chromosomes in rice and Brachypodium to 550 Mb in Avena

BACKGROUND

The BOP (Bambusoideae, Oryzoideae, and Pooideae) clade of the Poaceae has a common ancestor, with similarities to the genomes of rice, Oryza sativa (2n = 24; genome size 389 Mb) and Brachypodium, Brachypodium distachyon (2n = 10; 271 Mb). We exploit chromosome-scale genome assemblies to show the nature of genomic expansion, structural variation, and chromosomal rearrangements from rice and Brachypodium, to diploids in the tribe Aveneae (e.g., Avena longiglumis, 2n = 2x = 14; 3,961 Mb assembled to 3,850 Mb in chromosomes).

RESULTS

Most of the Avena chromosome arms show relatively uniform expansion over the 10-fold to 15-fold genome-size increase. Apart from non-coding sequence diversification and accumulation around the centromeres, blocks of genes are not interspersed with blocks of repeats, even in subterminal regions. As in the tribe Triticeae, blocks of conserved synteny are seen between the analyzed species with chromosome fusion, fission, and nesting (insertion) events showing deep evolutionary conservation of chromosome structure during genomic expansion. Unexpectedly, the terminal gene-rich chromosomal segments (representing about 50 Mb) show translocations between chromosomes during speciation, with homogenization of genome-specific repetitive elements within the tribe Aveneae. Newly-formed intergenomic translocations of similar extent are found in the hexaploid A. sativa.

CONCLUSIONS

The study provides insight into evolutionary mechanisms and speciation in the BOP clade, which is valuable for measurement of biodiversity, development of a clade-wide pangenome, and exploitation of genomic diversity through breeding programs in Poaceae.

Distribution and Location of BEVs in Different Genotypes of Bananas Reveal the Coevolution of BSVs and Bananas

Members of the family Caulimoviridae contain abundant endogenous pararetroviral sequences (EPRVs) integrated into the host genome. Banana streak virus (BSV), a member of the genus Badnavirus in this family, has two distinct badnaviral integrated sequences, endogenous BSV (eBSV) and banana endogenous badnavirus sequences (BEVs). BEVs are distributed widely across the genomes of different genotypes of bananas. To clarify the distribution and location of BEVs in different genotypes of bananas and their coevolutionary relationship with bananas and BSVs, BEVs and BSVs were identified in 102 collected banana samples, and a total of 327 BEVs were obtained and categorized into 26 BEVs species with different detection rates. However, the majority of BEVs were found in Clade II, and a few were clustered in Clade I. Additionally, BEVs and BSVs shared five common conserved motifs. However, BEVs had two unique amino acids, methionine and lysine, which differed from BSVs. BEVs were distributed unequally on most of chromosomes and formed hotspots. Interestingly, a colinear relationship of BEVs was found between AA and BB, as well as AA and SS genotypes of bananas. Notably, the chromosome integration time of different BEVs varied. Based on our findings, we propose that the coevolution of bananas and BSVs is driven by BSV Driving Force (BDF), a complex interaction between BSVs, eBSVs, and BEVs. This study provides the first clarification of the relationship between BEVs and the coevolution of BSVs and bananas in China.

Narrow genetic diversity in germplasm from the Guinean and Sudano-Guinean zones in Benin indicates the need to broaden the genetic base of sweet fig banana (Musa acuminata cv Sotoumon)

Sweet fig (M. acuminata cv. Sotoumon) is an economically important dessert banana in Benin, with high nutritional, medicinal, and cultural values. Nevertheless, its productivity and yield are threatened by biotic and abiotic stresses. Relevant knowledge of the genetic diversity of this economically important crop is essential for germplasm conservation and the development of breeding programs. However, very little is known about the genetic makeup of this cultivar in Benin. To advance the understanding of genetic diversity in sweet fig banana germplasm, a Genotype-By-Sequencing (GBS) was performed on a panel of 273 accessions collected in different phytogeographical zones of Benin. GBS generated 8,457 quality SNPs, of which 1992 were used for analysis after filtering. The results revealed a low diversity in the studied germplasm (He = 0.0162). Genetic differentiation was overall very low in the collection as suggested by the negative differentiation index (Fstg = -0.003). The Analysis of Molecular Variance (AMOVA) indicated that the variation between accessions within populations accounted for 83.8% of the total variation observed (P < 0.001). The analysis of population structure and neighbor-joining tree partitioned the germplasm into three clusters out of which a predominant major one contained 98.1% of all accessions. These findings demonstrate that current sweet fig banana genotypes shared a common genetic background, which made them vulnerable to biotic and abiotic stress. Therefore, broadening the genetic base of the crop while maintaining its quality attributes and improving yield performance is of paramount importance. Moreover, the large genetic group constitutes an asset for future genomic selection studies in the crop and can guide the profiling of its conservation strategies.

Chromosome-level genomes of three key Allium crops and their trait evolution

Allium crop breeding remains severely hindered due to the lack of high-quality reference genomes. Here we report high-quality chromosome-level genome assemblies for three key Allium crops (Welsh onion, garlic and onion), which are 11.17 Gb, 15.52 Gb and 15.78 Gb in size with the highest recorded contig N50 of 507.27 Mb, 109.82 Mb and 81.66 Mb, respectively. Beyond revealing the genome evolutionary process of Allium species, our pathogen infection experiments and comparative metabolomic and genomic analyses showed that genes encoding enzymes involved in the metabolic pathway of Allium-specific flavor compounds may have evolved from an ancient uncharacterized plant defense system widely existing in many plant lineages but extensively boosted in alliums. Using in situ hybridization and spatial RNA sequencing, we obtained an overview of cell-type categorization and gene expression changes associated with spongy mesophyll cell expansion during onion bulb formation, thus indicating the functional roles of bulb formation genes.

The role of orphan crops in the transition to nutritional quality-oriented crop improvement

Micronutrient malnutrition is a persisting problem threatening global human health. Biofortification via metabolic engineering has been proposed as a cost-effective and short-term means to alleviate this burden. There has been a recent rise in the recognition of potential that underutilized, orphan crops can hold in decreasing malnutrition concerns. Here, we illustrate how orphan crops can serve as a medium to provide micronutrients to populations in need, whilst promoting and maintaining dietary diversity. We provide a roadmap, illustrating which aspects to be taken into consideration when evaluating orphan crops. Recent developments have shown successful biofortification via metabolic engineering in staple crops. This review provides guidance in the implementation of these successes to relevant orphan crop species, with a specific focus on the relevant micronutrients iron, zinc, provitamin A and folates.

Genome-wide identification and gene expression analysis of GHMP kinase gene family in banana cv. Rasthali

BACKGROUND

The GHMP kinase gene family encompasses ATP-dependent kinases, significantly involved in the biosynthesis of isoprenes, amino acids, and metabolism of carbohydrates. Banana is a staple tropical crop that is globally consumed but known for high sensitivity to salt, cold, and drought stresses. The GHMP kinases are known to play a significant role during abiotic stresses in plants. The present study emphasizes the role of GHMP kinases in various abiotic stress conditions in banana.

METHODS AND RESULTS

We identified 12 GHMP kinase (MaGHMP kinase) genes in the banana genome database and witnessed the presence of the conserved Pro-X-X-X-Gly-Leu-X-Ser-Ser-Ala domain in their protein sequences. All genes were found to be involved in ATP-binding and carried kinase activity confronting their biological roles in the isoprene (27%) and amino acid (20%) biosyntheses. The expression analysis of genes during cold, drought, and salt stress conditions in tissue culture grown banana cultivar Rasthali plants showed a significant involvement of MaGHMP kinase genes in these stress conditions. The highest expression of MaGHMP kinase3 (8.5 fold) was noted during cold stress, while MaGHMP kinase1 (25 fold and 40.01 fold) showed maximum expression during drought and salt stress conditions in leaf tissue of Rasthali.

CONCLUSION

Our findings suggested that MaGHMP kinase1 (MaHSK) and MaGHMP kinase3 (MaGlcAK) could be considered promising candidates for thwarting the abiotic stresses in banana.

Let's go bananas with paintings

We tell the story of the complex history and global distribution of the very popular banana and illustrate this with amazing depictions by artists and scientists. In addition, we discuss the major changes that occurred during domestication and the risks associated with maintaining a monoculture.

Chromosome-scale genome assembly provides insights into flower coloration mechanisms of Canna indica

Canna indica, which produce conspicuous and colorful flowers, are widely appreciated as ornamental plants. We used Pacific Biosciences sequencing (PacBio) and chromosome conformation capture (Hi-C) genome scaffolding to build a high-quality chromosome-scale genome assembly of C. indica and the genome assembly was 821Mb with a contig N50 of 48Mb assembled into nine chromosomes. The genome of C. indica was predicted to contain 31,130 genes and 30,816 genes were functionally annotated. Genome annotation identified 522 Mb (63.59 %) as repetitive sequences. Genome evolution analysis showed that whole-genome duplication occurred 53.4 million years ago. Transcriptome analysis revealed that petal coloration was linked with the expression of genes encoding enzymes involved in anthocyanin biosynthesis, carotenoid biosynthesis, and the methylerythritol phosphate (MEP) pathway. Furthermore, modules of co-expressed genes and hub genes were identified via weighted gene co-expression network analysis. These results suggested that, in Canna indica, deep red petal coloration was regulated by CHS2 and yellow petal coloration was associated with expression of ARF6 and NAC14. Considered together, the current study revealed a high-quality reference genome which may provide new insights into the molecular basis of flower coloration in Canna indica and help enhance the conservation and breeding of ornamental plants in general.

A single amino acid substitution in MdLAZY1A dominantly impairs shoot gravitropism in Malus

Plant architecture is 1 of the most important factors that determines crop yield potential and productivity. In apple (Malus domestica), genetic improvement of tree architecture has been challenging due to a long juvenile phase and growth as complex trees composed of a distinct scion and a rootstock. To better understand the genetic control of apple tree architecture, the dominant weeping growth phenotype was investigated. We report the identification of MdLAZY1A (MD13G1122400) as the genetic determinant underpinning the Weeping (W) locus that largely controls weeping growth in Malus. MdLAZY1A is 1 of the 4 paralogs in apple that are most closely related to AtLAZY1 involved in gravitropism in Arabidopsis (Arabidopsis thaliana). The weeping allele (MdLAZY1A-W) contains a single nucleotide mutation c.584T>C that leads to a leucine to proline (L195P) substitution within a predicted transmembrane domain that colocalizes with Region III, 1 of the 5 conserved regions in LAZY1-like proteins. Subcellular localization revealed that MdLAZY1A localizes to the plasma membrane and nucleus in plant cells. Overexpressing the weeping allele in apple cultivar Royal Gala (RG) with standard growth habit impaired its gravitropic response and altered the growth to weeping-like. Suppressing the standard allele (MdLAZY1A-S) by RNA interference (RNAi) in RG similarly changed the branch growth direction to downward. Overall, the L195P mutation in MdLAZY1A is genetically causal for weeping growth, underscoring not only the crucial roles of residue L195 and Region III in MdLAZY1A-mediated gravitropic response but also a potential DNA base editing target for tree architecture improvement in Malus and other crops.

Exploring Pan-Genomes: An Overview of Resources and Tools for Unraveling Structure, Function, and Evolution of Crop Genes and Genomes

The availability of multiple sequenced genomes from a single species made it possible to explore intra- and inter-specific genomic comparisons at higher resolution and build clade-specific pan-genomes of several crops. The pan-genomes of crops constructed from various cultivars, accessions, landraces, and wild ancestral species represent a compendium of genes and structural variations and allow researchers to search for the novel genes and alleles that were inadvertently lost in domesticated crops during the historical process of crop domestication or in the process of extensive plant breeding. Fortunately, many valuable genes and alleles associated with desirable traits like disease resistance, abiotic stress tolerance, plant architecture, and nutrition qualities exist in landraces, ancestral species, and crop wild relatives. The novel genes from the wild ancestors and landraces can be introduced back to high-yielding varieties of modern crops by implementing classical plant breeding, genomic selection, and transgenic/gene editing approaches. Thus, pan-genomic represents a great leap in plant research and offers new avenues for targeted breeding to mitigate the impact of global climate change. Here, we summarize the tools used for pan-genome assembly and annotations, web-portals hosting plant pan-genomes, etc. Furthermore, we highlight a few discoveries made in crops using the pan-genomic approach and future potential of this emerging field of study.

The phased telomere-to-telomere reference genome of Musa acuminata, a main contributor to banana cultivars

Musa acuminata is a main wild contributor to banana cultivars. Here, we reported a haplotype-resolved and telomere-to-telomere reference genome of M. acuminata by incorporating PacBio HiFi reads, Nanopore ultra-long reads, and Hi-C data. The genome size of the two haploid assemblies was estimated to be 469.83 Mb and 470.21 Mb, respectively. Multiple assessments confirmed the contiguity (contig N50: 16.53 Mb and 18.58 Mb; LAI: 20.18 and 19.48), completeness (BUSCOs: 98.57% and 98.57%), and correctness (QV: 45.97 and 46.12) of the genome. The repetitive sequences accounted for about half of the genome size. In total, 40,889 and 38,269 protein-coding genes were annotated in the two haploid assemblies, respectively, of which 9.56% and 3.37% were newly predicted. Genome comparison identified a large reciprocal translocation involving 3 Mb and 10 Mb from chromosomes 01 and 04 within M. acuminata. This reference genome of M. acuminata provides a valuable resource for further understanding of subgenome evolution of Musa species, and precise genetic improvement of banana.

Genome-wide characterization of 2OGD superfamily for mining of susceptibility factors responding to various biotic stresses in Musa spp

Bananas are an important staple food and cash crop, but they are vulnerable to a variety of pests and diseases that substantially reduce yield and quality. Banana diseases are challenging to control and necessitate an integrated strategy, and development of resistant cultivars is one of the effective ways of managing diseases. Lasting disease resistance is the main goal in crop improvement and resistance mediated by a single resistant (R) gene mostly lack durability. However, long-term resistance can be obtained by inactivating susceptibility factors (S), which facilitate pathogen infection and proliferation. Identification and inactivation of susceptibility factors against the major pathogens like Fusarium oxysporum f. sp. cubense (Foc), Pseudocercospora eumusae and Pratylenchus coffeae in banana will be an effective way in developing banana varieties with more durable resistance. Downy mildew resistance 6 (DMR6) and DMR-like oxygenases (DLO1) are one such susceptibility factors and they belong to 2-oxoglutarate Fe(II) dependent oxygenases (2OGD) superfamily. 2OGDs are known to catalyze a plethora of reactions and also confer resistance to different pathogens in various crops, but not much is known about the 2OGD in Musa species. Through a comprehensive genome-wide analysis, 133 and 122 potential 2OGDs were systematically identified and categorized from the A and B genomes of banana, respectively. Real time expression of dmr6 and dlo1 genes showed positive correlation with transcriptome data upon Foc race1 and TR4 infection and examination of expression pattern of Macma4_04_g22670 (Ma04_g20880) and Macma4_02_g13590 (Ma02_g12040) genes revealed their involvement in Foc race1 and TR4 infections, respectively. Further the expression profile of 2OGDs, specifically Macma4_04_g25310 (Ma04_g23390), Macma4_08_g11980 (Ma08_g12090) and Macma4_04_g38910 (Ma04_g36640) shows that they may play a significant role as a susceptibility factor, particularly against P. eumusae and P. coffeae, implying that they can be exploited as a candidate gene for editing in developing resistant cultivars against these diseases. In summary, our findings contribute to a deeper comprehension of the evolutionary and functional aspects of 2OGDs in Musa spp. Furthermore, they highlight the substantial functions of these family constituents in the progression of diseases. These insights hold significance in the context of enhancing the genetic makeup of bananas to attain extended and more durable resistance against pathogens.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01380-y.

Telomere-to-telomere haplotype-resolved reference genome reveals subgenome divergence and disease resistance in triploid Cavendish banana

Banana is one of the most important crops of the world. Cavendish-type bananas, which have a monospecific Musa acuminata origin (AAA), account for around half of the global banana production, thereby are of great significance for human societies. However, until now, the high-quality haplotype-resolved reference genome was still undecoded for banana cultivars. Here, we reported the telomere-to-telomere (T2T) and haplotype-resolved reference genome of 'Baxijiao' (Cavendish) consisting of three haploid assemblies. The sizes of the three haploid assemblies were estimated to be 477.16 Mb, 477.18 Mb, and 469.57 Mb, respectively. Although with monospecific origins, the three haploid assemblies showed great differences with low levels of sequence collinearity. Several large reciprocal translocations were identified among chromosomes 1, 4, and 7. An expansion of gene families that might affect fruit quality and aroma was detected, such as those belonging to sucrose/disaccharide/oligosaccharide catabolic processes, sucrose metabolic process, starch metabolic process, and aromatic compound biosynthetic process. Besides, an expansion of gene families related to anther and pollen development was observed, which could be associated with parthenocarpy and sterility of the Cavendish cultivar. Finally, much fewer resistance genes were identified in 'Baxijiao' than in M. acuminata, particularly in the gene clusters in chromosomes 3 and 10, providing potential targets to explore for molecular analysis of disease resistance in banana. This T2T haplotype-resolved reference genome will thus be a valuable genetic resource for biological studies, molecular breeding, and genetic improvement of banana.

Genome-wide identification and expression analysis of the GRAS family under low-temperature stress in bananas

Introduction

GRAS, named after GAI, RGA, and SCR, is a class of plant-specific transcription factors family that plays a crucial role in growth and development, signal transduction, and various stress responses.

Methods

To understand the biological functions of the banana GRAS gene family, a genome-wide identification and bioinformatics analysis of the banana GRAS gene family was performed based on information from the M. acuminata, M. balbisiana, and M. itinerans genomic databases.

Result

In the present study, we identified 73 MaGRAS, 59 MbGRAS, and 58 MiGRAS genes in bananas at the whole-genome scale, and 56 homologous genes were identified in the three banana genomes. Banana GRASs can be classified into 10 subfamilies, and their gene structures revealed that most banana GRAS gDNAs lack introns. The promoter sequences of GRASs had a large number of cis-acting elements related to plant growth and development, phytohormone, and adversity stress responsiveness. The expression pattern of seven key members of MaGRAS response to low-temperature stress and different tissues was also examined by quantitative reverse transcription polymerase chain reaction (qRT-PCR). The microRNAs-MaGRASs target prediction showed perfect complementarity of seven GRAS genes with the five mac-miRNAs. The expression of all seven genes was lowest in roots, and the expression of five genes was highest in leaves during low-temperature stress. The expression of MaSCL27-2, MaSCL27-3, and MaSCL6-1 was significantly lower under low-temperature stress compared to the control, except for MaSCL27-2, which was slightly higher than the 28°C control at 4 h. The expression of MaSCL27-2, MaSCL27-3, and MaSCL6-1 dropped to the lowest levels at 24 h, 12 h, and 4 h, respectively. The MaSCL27-4 and MaSCL6-2 expression was intermittently upregulated, rising to the highest expression at 24h, while the expression of MaSCL22 was less variable, remaining at the control level with small changes.

Discussion

In summary, it is tentatively hypothesized that the GRAS family has an important function in low-temperature stress in bananas. This study provides a theoretical basis for further analyzing the function of the banana GRAS gene and the resistance of bananas to cold temperatures.

Chromosome-level reference genome assembly provides insights into the evolution of Pennisetum alopecuroides

Pennisetum alopecuroides is an important forage grass resource, which plays a vital role in ecological environment improvement. Therefore, the acquisition of P. alopecuroides genome resources is conducive to the study of the adaptability of Pennisetum species in ecological remediation and forage breeding development. Here we assembled a P. alopecuroides cv. 'Liqiu' genome at the chromosome level with a size of approximately 845.71 Mb, contig N50 of 84.83Mb, and genome integrity of 99.13% as assessed by CEGMA. A total of 833.41-Mb sequences were mounted on nine chromosomes by Hi-C technology. In total, 60.66% of the repetitive sequences and 34,312 genes were predicted. The genomic evolution analysis showed that P. alopecuroides cv. 'Liqiu' was isolated from Setaria 7.53-13.80 million years ago and from Cenchrus 5.33-8.99 million years ago, respectively. The whole-genome event analysis showed that P. alopecuroides cv. 'Liqiu' underwent two whole-genome duplication (WGD) events in the evolution process, and the duplication events occurred at a similar time to that of Oryza sativa and Setaria viridis. The completion of the genome sequencing of P. alopecuroides cv. 'Liqiu' provides data support for mining high-quality genetic resources of P. alopecuroides and provides a theoretical basis for the origin and evolutionary characteristics of Pennisetum.

Effect of Potassium Permanganate as an Ethylene Scavenger and Physicochemical Characterization during the Shelf Life of Fresh Banana (Musa paradisiaca)

The conventional method of employing low temperatures for storage and distribution has long been the standard approach for preserving most fruits and vegetables. This practice is likewise prevalent in the retail industry, which relies on similar methods for transporting and maintaining the quality of perishable products on their shelves. The aim was to preserve bananas (Musa paradisiaca) using an ethylene scavenger, potassium permanganate, which is contained in small paper bags, to increase the storage and distribution time at low cost. The bananas were distributed in four plastic containers at a temperature of 23°C, three of the treatments contained different concentrations of potassium permanganate, and one was potassium permanganate free. The experimental period was 19 days, and the variations in weight loss, pH, titratable acidity, texture, color, and total soluble solids were analyzed. Potassium permanganate effectively reduced the changes in their physiological ripening.

Genome-Wide Identification of the MPK Gene Family and Expression Analysis under Low-Temperature Stress in the Banana

Mitogen-activated protein kinases (MAPKs and MPKs) are important in the process of resisting plant stress. In this study, 21, 12, 18, 16, and 10 MPKs were identified from Musa acuminata, Musa balbisiana, Musa itinerans, Musa schizocarpa, and Musa textilis, respectively. These MPKs were divided into Group A, B, C, and D. Phylogenetic analysis revealed that this difference in number was due to the gene shrinkage of the Group B subfamily of Musa balbisiana and Musa textilis. KEGG annotations revealed that K14512, which is involved in plant hormone signal transduction and the plant-pathogen interaction, was the most conserved pathway of the MPKs. The results of promoter cis-acting element prediction and focTR4 (Fusarium oxysporum f. sp. cubense tropical race 4) transcriptome expression analysis preliminarily confirmed that MPKs were relevant to plant hormone and biotic stress, respectively. The expression of MPKs in Group A was significantly upregulated at 4 °C, and dramatically, the MPKs in the root were affected by low temperature. miR172, miR319, miR395, miR398, and miR399 may be the miRNAs that regulate MPKs during low-temperature stress, with miR172 being the most critical. miRNA prediction and qRT-PCR results indicated that miR172 may negatively regulate MPKs. Therefore, we deduced that MPKs might coordinate with miR172 to participate in the process of the resistance to low-temperature stress in the roots of the banana. This study will provide a theoretical basis for further analysis of the mechanism of MPKs under low-temperature stress of bananas, and this study could be applied to molecular breeding of bananas in the future.

Banana disease-suppressive soil drives Bacillus assembled to defense Fusarium wilt of banana

Fusarium wilt of banana (FWB) caused by Fusarium oxysporum f. sp. cubense tropical race 4 (Foc TR4), poses a serious problem for sustainable banana production. Biological control is one of the effective measures to control this destructive disease. High-throughput sequencing of soil microorganisms could significantly improve the efficiency and accuracy of biocontrol strain screening. In this study, the soil microbial diversity of six main banana-producing areas in Yunnan was sequenced by Illumina Miseq platform. The outcome of this study showed the genus of Chujaibacter, Bacillus, and Sphingomonas were significantly enriched in microorganism community composition. Further correlation analysis with soil pathogen (Foc TR4) content showed that Bacillus was significantly negatively correlated with pathogen content. Therefore, we isolated and identified Bacillus from the disease-suppressive soils, and obtained a B. velezensis strain YN1910. In vitro and pot experiments showed that YN1910 had a significant control effect (78.43-81.76%) on banana Fusarium wilt and had a significant growth promotion effect on banana plants.

Towards the Investigation of the Adaptive Divergence in a Species of Exceptional Ecological Plasticity: Chromosome-Scale Genome Assembly of Chouardia litardierei (Hyacinthaceae)

One of the central goals of evolutionary biology is to understand the genomic basis of adaptive divergence. Different aspects of evolutionary processes should be studied through genome-wide approaches, therefore maximizing the investigated genomic space. However, in-depth genome-scale analyses often are restricted to a model or economically important species and their closely related wild congeners with available reference genomes. Here, we present the high-quality chromosome-level genome assembly of Chouardia litardierei, a plant species with exceptional ecological plasticity. By combining PacBio and Hi-C sequencing technologies, we generated a 3.7 Gbp genome with a scaffold N50 size of 210 Mbp. Over 80% of the genome comprised repetitive elements, among which the LTR retrotransposons prevailed. Approximately 86% of the 27,257 predicted genes were functionally annotated using public databases. For the comparative analysis of different ecotypes' genomes, the whole-genome sequencing of two individuals, each from a distinct ecotype, was performed. The detected above-average SNP density within coding regions suggests increased adaptive divergence-related mutation rates, therefore confirming the assumed divergence processes within the group. The constructed genome presents an invaluable resource for future research activities oriented toward the investigation of the genetics underlying the adaptive divergence that is likely unfolding among the studied species' ecotypes.

Diploid and tetraploid genomes of Acorus and the evolution of monocots

Monocots are a major taxon within flowering plants, have unique morphological traits, and show an extraordinary diversity in lifestyle. To improve our understanding of monocot origin and evolution, we generate chromosome-level reference genomes of the diploid Acorus gramineus and the tetraploid Ac. calamus, the only two accepted species from the family Acoraceae, which form a sister lineage to all other monocots. Comparing the genomes of Ac. gramineus and Ac. calamus, we suggest that Ac. gramineus is not a potential diploid progenitor of Ac. calamus, and Ac. calamus is an allotetraploid with two subgenomes A, and B, presenting asymmetric evolution and B subgenome dominance. Both the diploid genome of Ac. gramineus and the subgenomes A and B of Ac. calamus show clear evidence of whole-genome duplication (WGD), but Acoraceae does not seem to share an older WGD that is shared by most other monocots. We reconstruct an ancestral monocot karyotype and gene toolkit, and discuss scenarios that explain the complex history of the Acorus genome. Our analyses show that the ancestors of monocots exhibit mosaic genomic features, likely important for that appeared in early monocot evolution, providing fundamental insights into the origin, evolution, and diversification of monocots.

A near-complete genome assembly of Thalia dealbata Fraser (Marantaceae)

This study presents a chromosome-level, near-complete genome assembly of Thalia dealbata (Marantaceae), a typical emergent wetland plant with high ornamental and environmental value. Based on 36.99 Gb PacBio HiFi reads and 39.44 Gb Hi-C reads, we obtained a 255.05 Mb assembly, of which 251.92 Mb (98.77%) were anchored into eight pseudo-chromosomes. Five pseudo-chromosomes were completely assembled, and the other three had one to two gaps. The final assembly had a high contig N50 value (29.80 Mb) and benchmarking universal single-copy orthologs (BUSCO) recovery score (97.52%). The T. dealbata genome had 100.35 Mb repeat sequences, 24,780 protein-coding genes, and 13,679 non-coding RNAs. Phylogenetic analysis revealed that T. dealbata was closest to Zingiber officinale, whose divergence time was approximately 55.41 million years ago. In addition, 48 and 52 significantly expanded and contracted gene families were identified within the T. dealbata genome. Moreover, 309 gene families were specific to T. dealbata, and 1,017 genes were positively selected. The T. dealbata genome reported in this study provides a valuable genomic resource for further research on wetland plant adaptation and the genome evolution dynamics. This genome is also beneficial for the comparative genomics of Zingiberales species and flowering plants.

Three R2R3-MYB transcription factors from banana (Musa acuminata) activate structural anthocyanin biosynthesis genes as part of an MBW complex

OBJECTIVE

Bananas are one of the most popular fruits in the world, providing food security and employment opportunities in several developing countries. Increasing the anthocyanin content of banana fruit could improve the health-promoting properties. Anthocyanin biosynthesis is largely regulated at the transcriptional level. However, relatively little is known about the transcriptional activation of anthocyanin biosynthesis in banana.

RESULTS

We analysed the regulatory activity of three Musa acuminata MYBs that were predicted by bioinformatic analysis to transcriptionally regulate anthocyanin biosynthesis in banana. MaMYBA1, MaMYBA2 and MaMYBPA2 did not complement the anthocyanin-deficient phenotype of the Arabidopsis thaliana pap1/pap2 mutant. However, co-transfection experiments in A. thaliana protoplasts showed that MaMYBA1, MaMYBA2 and MaMYBPA2 function as components of a transcription factor complex with a bHLH and WD40 protein, the so called MBW complex, resulting in the activation of the A. thaliana ANTHOCYANIDIN SYNTHASE and DIHYDROFLAVONOL 4-REDUCTASE promoters. The activation potential of MaMYBA1, MaMYBA2 and MaMYBPA2 was increased when combined with the monocot Zea mays bHLH ZmR instead of the dicot AtEGL3. This work paves the path towards decoding the MBW complex-mediated transcriptional activation of anthocyanin biosynthesis in banana. It will also facilitate research towards increased anthocyanin content in banana and other monocot crops.

Genetic Mapping, Candidate Gene Identification and Marker Validation for Host Plant Resistance to the Race 4 of Fusarium oxysporum f. sp. cubense Using Musa acuminata ssp. malaccensis

Fusarium wilt of banana is a devastating disease that has decimated banana production worldwide. Host resistance to Fusarium oxysporum f. sp. Cubense (Foc), the causal agent of this disease, is genetically dissected in this study using two Musa acuminata ssp. Malaccensis segregating populations, segregating for Foc Tropical (TR4) and Subtropical (STR4) race 4 resistance. Marker loci and trait association using 11 SNP-based PCR markers allowed the candidate region to be delimited to a 12.9 cM genetic interval corresponding to a 959 kb region on chromosome 3 of 'DH-Pahang' reference assembly v4. Within this region, there was a cluster of pattern recognition receptors, namely leucine-rich repeat ectodomain containing receptor-like protein kinases, cysteine-rich cell-wall-associated protein kinases, and leaf rust 10 disease-resistance locus receptor-like proteins, positioned in an interspersed arrangement. Their transcript levels were rapidly upregulated in the resistant progenies but not in the susceptible F₂ progenies at the onset of infection. This suggests that one or several of these genes may control resistance at this locus. To confirm the segregation of single-gene resistance, we generated an inter-cross between the resistant parent 'Ma850' and a susceptible line 'Ma848', to show that the STR4 resistance co-segregated with marker '28820' at this locus. Finally, an informative SNP marker 29730 allowed the locus-specific resistance to be assessed in a collection of diploid and polyploid banana plants. Of the 60 lines screened, 22 lines were predicted to carry resistance at this locus, including lines known to be TR4-resistant, such as 'Pahang', 'SH-3362', 'SH-3217', 'Ma-ITC0250', and 'DH-Pahang/CIRAD 930'. Additional screening in the International Institute for Tropical Agriculture's collection suggests that the dominant allele is common among the elite 'Matooke' NARITA hybrids, as well as in other triploid or tetraploid hybrids derived from East African highland bananas. Fine mapping and candidate gene identification will allow characterization of molecular mechanisms underlying the TR4 resistance. The markers developed in this study can now aid the marker-assisted selection of TR4 resistance in breeding programs around the world.

Small cardamom genome: development and utilization of microsatellite markers from a draft genome sequence of Elettaria cardamomum Maton

Small cardamom (Elettaria cardamomum Maton), the queen of spices, is the third most expensive spice in the world after saffron and vanilla, valued highly for its aroma and taste. This perennial herbaceous plant is a native of coastal parts of Southern India and displays a significant amount of morphological diversity. Its genetic potential has not been exploited due to lack of genomic resources limiting our understanding of the genome and important metabolic pathways which give it the economic advantage in the spice industry. Here, we report upon the de novo assembled, draft whole genome sequence of cardamom variety, Njallani Green Gold. We used a hybrid assembly strategy using the reads from the Oxford Nanopore, Illumina and 10x Genomics GemCode sequencing chemistries. The assembled genome length was 1.06 Gb (gigabases) which is close to the estimated genome size of cardamom. More than 75% of the genome was captured in 8000 scaffolds with a N50 of 0.15 Mb. The genome appears to have a high repeat content and 68055 gene models were predicted. The genome is close to Musa species and displays an expansion and contraction in different gene families. The draft assembly was used for in silico mining of simple sequence repeats (SSRs). A total of 2,50,571 SSRs were identified of which 2,18,270 were perfect SSRs and 32,301 were compound SSRs. Among the perfect SSRs, trinucleotides were most abundant (1,25,329) and hexanucleotide repeats appear least (2,380). From the 2,50,571 SSRs mined, 2,27,808 primer pairs were designed based on flanking sequence information. Wet lab validation was performed for 246 SSR loci and based on their amplification profiles, 60 SSR markers were used for diversity analysis of a set of 60 diverse cardamom accessions. The average number of alleles detected per locus were 14.57 with a minimum of 4 and maximum of 30 alleles. Population structure analysis revealed the presence of high degree of admixtures which could primarily be due to cross-pollination prevalent in this species. The SSR markers identified would help in the development of gene or trait-linked markers which can be subsequently used for marker-assisted breeding for crop improvement in cardamom. The information on utilization of the SSR loci for generation of markers has been developed into a public database, 'cardamomSSRdb' that is freely available for use by the cardamom community.

Overexpression of banana GDP-L-galactose phosphorylase (GGP) modulates the biosynthesis of ascorbic acid in Arabidopsis thaliana

l-Ascorbic acid (AsA) is a potent antioxidant and essential micronutrient for the growth and development of plants and animals. AsA is predominantly synthesized by the Smirnoff-Wheeler (SW) pathway in plants where the GDP-L-galactose phosphorylase (GGP) gene encodes the rate-limiting step. In the present study, AsA was estimated in twelve banana cultivars, where Nendran carried the highest (17.2 mg/100 g) amount of AsA in ripe fruit pulp. Five GGP genes were identified from the banana genome database, and they were located at chromosome 6 (4 MaGGPs) and chromosome 10 (1 MaGGP). Based on in-silico analysis, three potential MaGGP genes were isolated from the cultivar Nendran and subsequently overexpressed in Arabidopsis thaliana. Significant enhancement in AsA (1.52 to 2.20 fold) level was noted in the leaves of all three MaGGPs overexpressing lines as compared to non-transformed control plants. Among all, MaGGP2 emerged as a potential candidate for AsA biofortification in plants. Further, the complementation assay of Arabidopsis thaliana vtc-5-1 and vtc-5-2 mutants with MaGGP genes overcome the AsA deficiency that showed improved plant growth as compared to non-transformed control plants. This study lends strong affirmation towards development of AsA biofortified plants, particularly the staples that sustain the personages in developing countries.

Metapangenomics of wild and cultivated banana microbiome reveals a plethora of host-associated protective functions

BACKGROUND

Microbiomes are critical to plants, promoting growth, elevating stress tolerance, and expanding the plant's metabolic repertoire with novel defense pathways. However, generally microbiomes within plant tissues, which intimately interact with their hosts, remain poorly characterized. These endospheres have become a focus in banana (Musa spp.)-an important plant for study of microbiome-based disease protection. Banana is important to global food security, while also being critically threatened by pandemic diseases. Domestication and clonal propagation are thought to have depleted protective microbiomes, whereas wild relatives may hold promise for new microbiome-based biological controls. The goal was to compare metapangenomes enriched from 7 Musa genotypes, including wild and cultivated varieties grown in sympatry, to assess the host associations with root and leaf endosphere functional profiles.

RESULTS

Density gradients successfully generated culture-free microbial enrichment, dominated by bacteria, with all together 24,325 species or strains distinguished, and 1.7 million metagenomic scaffolds harboring 559,108 predicted gene clusters. About 20% of sequence reads did not match any taxon databases and ~ 62% of gene clusters could not be annotated to function. Most taxa and gene clusters were unshared between Musa genotypes. Root and corm tissues had significantly richer endosphere communities that were significantly different from leaf communities. Agrobacterium and Rhizobium were the most abundant in all samples while Chitinophagia and Actinomycetia were more abundant in roots and Flavobacteria in leaves. At the bacterial strain level, there were > 2000 taxa unique to each of M. acuminata (AAA genotype) and M. balbisiana (B-genotype), with the latter 'wild' relatives having richer taxa and functions. Gene ontology functional enrichment showed core beneficial functions aligned with those of other plants but also many specialized prospective beneficial functions not reported previously. Some gene clusters with plant-protective functions showed signatures of phylosymbiosis, suggesting long-standing associations or heritable microbiomes in Musa.

CONCLUSIONS

Metapangenomics revealed key taxa and protective functions that appeared to be driven by genotype, perhaps contributing to host resistance differences. The recovery of rich novel taxa and gene clusters provides a baseline dataset for future experiments in planta or in vivo bacterization or engineering of wild host endophytes.

DNA methylation signatures of duplicate gene evolution in angiosperms

Gene duplication is a source of evolutionary novelty. DNA methylation may play a role in the evolution of duplicate genes (paralogs) through its association with gene expression. While this relationship has been examined to varying extents in a few individual species, the generalizability of these results at either a broad phylogenetic scale with species of differing duplication histories or across a population remains unknown. We applied a comparative epigenomics approach to 43 angiosperm species across the phylogeny and a population of 928 Arabidopsis (Arabidopsis thaliana) accessions, examining the association of DNA methylation with paralog evolution. Genic DNA methylation was differentially associated with duplication type, the age of duplication, sequence evolution, and gene expression. Whole genome duplicates were typically enriched for CG-only gene-body methylated or unmethylated genes, while single-gene duplications were typically enriched for non-CG methylated or unmethylated genes. Non-CG methylation, in particular, was characteristic of more recent single-gene duplicates. Core angiosperm gene families differentiated into those which preferentially retain paralogs and 'duplication-resistant' families, which convergently reverted to singletons following duplication. Duplication-resistant families that still have paralogous copies were, uncharacteristically for core angiosperm genes, enriched for non-CG methylation. Non-CG methylated paralogs had higher rates of sequence evolution, higher frequency of presence-absence variation, and more limited expression. This suggests that silencing by non-CG methylation may be important to maintaining dosage following duplication and be a precursor to fractionation. Our results indicate that genic methylation marks differing evolutionary trajectories and fates between paralogous genes and have a role in maintaining dosage following duplication.

Characterization of Banana SNARE Genes and Their Expression Analysis under Temperature Stress and Mutualistic and Pathogenic Fungal Colonization

SNAREs (soluble N-ethylmaleimide-sensitive-factor attachment protein receptors) are engines for almost all of the membrane fusion and exocytosis events in organism cells. In this study, we identified 84 SNARE genes from banana (Musa acuminata). Gene expression analysis revealed that the expression of MaSNAREs varied a lot in different banana organs. By analyzing their expression patterns under low temperature (4 °C), high temperature (45 °C), mutualistic fungus (Serendipita indica, Si) and fungal pathogen (Fusarium oxysporum f. sp. Cubense Tropical Race 4, FocTR4) treatments, many MaSNAREs were found to be stress responsive. For example, MaBET1d was up-regulate by both low and high temperature stresses; MaNPSN11a was up-regulated by low temperature but down-regulated by high temperature; and FocTR4 treatment up-regulated the expression of MaSYP121 but down-regulated MaVAMP72a and MaSNAP33a. Notably, the upregulation or downregulation effects of FocTR4 on the expression of some MaSNAREs could be alleviated by priorly colonized Si, suggesting that they play roles in the Si-enhanced banana wilt resistance. Foc resistance assays were performed in tobacco leaves transiently overexpressing MaSYP121, MaVAMP72a and MaSNAP33a. Results showed that transient overexpression of MaSYP121 and MaSNPA33a suppressed the penetration and spread of both Foc1 (Foc Race 1) and FocTR4 in tobacco leaves, suggesting that they play positive roles in resisting Foc infection. However, the transient overexpression of MaVAMP72a facilitated Foc infection. Our study can provide a basis for understanding the roles of MaSNAREs in the banana responses to temperature stress and mutualistic and pathogenic fungal colonization.