Musa balbisiana genome reveals subgenome evolution and functional divergence

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.

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.

The Cissus quadrangularis genome reveals its adaptive features in an arid habitat

Cissus quadrangularis is a tetraploid species belonging to the Vitaceae family and is known for the Crassulacean acid metabolism (CAM) pathway in the succulent stem, while the leaves perform C3 photosynthesis. Here, we report a high-quality genome of C. quadrangularis comprising a total size of 679.2 Mb which was phased into two subgenomes. Genome annotation identified 51 857 protein-coding genes, while approximately 47.75% of the genome was composed of repetitive sequences. Gene expression ratios of two subgenomes demonstrated that the sub-A genome as the dominant subgenome played a vital role during the drought tolerance. Genome divergence analysis suggests that the tetraploidization event occurred around 8.9 million years ago. Transcriptome data revealed that pathways related to cutin, suberine, and wax metabolism were enriched in the stem during drought treatment, suggesting that these genes contributed to the drought adaption. Additionally, a subset of CAM-related genes displayed diurnal expression patterns in the succulent stems but not in leaves, indicating that stem-biased expression of existing genes contributed to the CAM evolution. Our findings provide insights into the mechanisms of drought adaptation and photosynthesis transition in plants.

Fusarium wilt pandemic: current understanding and molecular perspectives

Plant diseases pose a severe threat to the food security of the global human population. One such disease is Fusarium wilt, which affects many plant species and causes up to 100% yield losses. Fusarium pathogen has high variability in its genetic constitution; therefore, it has evolved into different physiological races to infect different plant species spread across the different geographical regions of the world. The pathogen mainly affects plant roots, leading to colonizing and blocking vascular bundle cells, specifically xylem vessels. This blocking results in chlorosis, vascular discoloration, leaf wilting, shortening of plant, and, in severe cases, premature plant death. Due to the soil-borne nature of the wilt pathogen, neither agronomic nor plant protection measures effectively reduce the incidence of the disease. Therefore, the most cost-effective management strategy for Fusarium wilt is developing varieties resistant to a particular race of the fungus wilt prevalent in a given region. This strategy requires understanding the pathogen, its disease cycle, and epidemiology with climate-changing scenarios. Hence, in the review, we will discuss the pathogenic aspect and genetics of the Fusarium wilt, including molecular interventions for developing climate-smart wilt tolerant/resistant varieties of crops. Overall, this review will add to our knowledge for advancing the breeding of resistance against the wilt pandemic.

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.

The Zinc Finger Protein MaCCCH33-Like2 Positively Regulates Banana Fruit Ripening by Modulating Genes in Starch and Cell Wall Degradation

As zinc finger protein transcription factors (TFs), the molecular mechanism of Cys-Cys-Cys-His (CCCH) TFs in regulating plant development, growth and stress response has been well studied. However, the roles of CCCH TFs in fruit ripening are still obscure. Herein, we report that MaCCCH33-like2 TF and its associated proteins modulate the fruit softening of 'Fenjiao' bananas. MaCCCH33-like2 interacts directly with the promoters of three genes: isoamylase2 (MaISA2), sugar transporter14-like (MaSUR14-like) and β-d-xylosidase23 (MaXYL23), all of which are responsible for encoding proteins involved in the degradation of starch and cell wall components. Additionally, MaCCCH33-like2 forms interactions with abscisic acid-insensitive 5 (ABI5)-like and ethylene F-box protein 1 (MaEBF1), resulting in enhanced binding and activation of promoters of genes related to starch and cell wall degradation. When MaCCCH33-like2 is transiently and ectopically overexpressed in 'Fenjiao' banana and tomato fruit, it facilitates softening and ripening processes by promoting the degradation of cell wall components and starch and the production of ethylene. Conversely, the temporary silencing of MaCCCH33-like2 using virus-induced gene silencing (VIGS) inhibits softening and ripening in the 'Fenjiao' banana by suppressing ethylene synthesis, as well as starch and cell wall degradation. Furthermore, the promoter activity of MaCCCH33-like2 is regulated by MaABI5-like. Taken together, we have uncovered a novel MaCCCH33-like2/MaEBF1/MaABI5-like module that participates in fruit softening regulation in bananas.

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.

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.

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 structure of the tetraploid sour cherry 'Schattenmorelle' (Prunus cerasus L.) genome reveals insights into its segmental allopolyploid nature

Sour cherry (Prunus cerasus L.) is an important allotetraploid cherry species that evolved in the Caspian Sea and Black Sea regions from a hybridization of the tetraploid ground cherry (Prunus fruticosa Pall.) and an unreduced pollen of the diploid sweet cherry (P. avium L.) ancestor. Details of when and where the evolution of this species occurred are unclear, as well as the effect of hybridization on the genome structure. To gain insight, the genome of the sour cherry cultivar 'Schattenmorelle' was sequenced using Illumina NovaSeqTM and Oxford Nanopore long-read technologies, resulting in a ~629-Mbp pseudomolecule reference genome. The genome could be separated into two subgenomes, with subgenome PceS_a originating from P. avium and subgenome PceS_f originating from P. fruticosa. The genome also showed size reduction compared to ancestral species and traces of homoeologous sequence exchanges throughout. Comparative analysis confirmed that the genome of sour cherry is segmental allotetraploid and evolved very recently in the past.

The Identification of the Banana Endogenous Reference Gene MaSPS1 and the Construction of Qualitative and Quantitative PCR Detection Methods

Endogenous reference genes play a crucial role in the qualitative and quantitative PCR detection of genetically modified crops. Currently, there are no systematic studies on the banana endogenous reference gene. In this study, the MaSPS1 gene was identified as a candidate gene through bioinformatics analysis. The conservation of this gene in different genotypes of banana was tested using PCR, and its specificity in various crops and fruits was also examined. Southern blot analysis showed that there is only one copy of MaSPS1 in banana. The limit of detection (LOD) test showed that the LOD of the conventional PCR method is approximately 20 copies. The real-time fluorescence quantitative PCR (qPCR) method also exhibited high specificity, with a LOD of approximately 10 copies. The standard curve of the qPCR method met the quantitative requirements, with a limit of quantification (LOQ) of 1.14 × 10-2 ng-about 20 copies. Also, the qPCR method demonstrated good repeatability and stability. Hence, the above results indicate that the detection method established in this study has strong specificity, a low detection limit, and good stability. It provides a reliable qualitative and quantitative detection system for banana.

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.

Decoupling subgenomes within hybrid lavandin provide new insights into speciation and monoterpenoid diversification of Lavandula

Polyploidization and transposon elements contribute to shape plant genome diversity and secondary metabolic variation in some edible crops. However, the specific contribution of these variations to the chemo-diversity of Lamiaceae, particularly in economic shrubs, is still poorly documented. The rich essential oils (EOs) of Lavandula plants are distinguished by monoterpenoids among the main EO-producing species, L. angustifolia (LA), L. × intermedia (LX) and L. latifolia (LL). Herein, the first allele-aware chromosome-level genome was assembled using a lavandin cultivar 'Super' and its hybrid origin was verified by two complete subgenomes (LX-LA and LX-LL). Genome-wide phylogenetics confirmed that LL, like LA, underwent two lineage-specific WGDs after the γ triplication event, and their speciation occurred after the last WGD. Chloroplast phylogenetic analysis indicated LA was the maternal source of 'Super', which produced premium EO (higher linalyl/lavandulyl acetate and lower 1,8-cineole and camphor) close to LA. Gene expression, especially the monoterpenoid biosynthetic genes, showed bias to LX-LA alleles. Asymmetric transposon insertions in two decoupling 'Super' subgenomes were responsible for speciation and monoterpenoid divergence of the progenitors. Both hybrid and parental evolutionary analysis revealed that LTR (long terminal repeat) retrotransposon associated with AAT gene loss cause no linalyl/lavandulyl acetate production in LL, and multi-BDH copies retained by tandem duplication and DNA transposon resulted in higher camphor accumulation of LL. Advances in allelic variations of monoterpenoids have the potential to revolutionize future lavandin breeding and EO production.

MaMADS1-MaNAC083 transcriptional regulatory cascade regulates ethylene biosynthesis during banana fruit ripening

The hormone ethylene is crucial in the regulation of ripening in climacteric fruit, such as bananas. The transcriptional regulation of ethylene biosynthesis throughout banana fruit ripening has received much study, but the cascaded transcriptional machinery of upstream transcriptional regulators implicated in the ethylene biosynthesis pathway is still poorly understood. Here we report that ethylene biosynthesis genes, including MaACS1, MaACO1, MaACO4, MaACO5, and MaACO8, were upregulated in ripening bananas. NAC (NAM, ATAF, CUC) transcription factor, MaNAC083, a ripening and ethylene-inhibited gene, was discovered as a potential binding protein to the MaACS1 promoter by yeast one-hybrid screening. Further in vitro and in vivo experiments indicated that MaNAC083 bound directly to promoters of the five ethylene biosynthesis genes, thereby transcriptionally repressing their expression, which was further verified by transient overexpression experiments, where ethylene production was inhibited through MaNAC083-modulated transcriptional repression of ethylene biosynthesis genes in banana fruits. Strikingly, MaMADS1, a ripening-induced MADS (MCM1, AGAMOUS, DEFICIENS, SRF4) transcription factor, was found to directly repress the expression of MaNAC083, inhibiting trans-repression of MaNAC083 to ethylene biosynthesis genes, thereby attenuating MaNAC083-repressed ethylene production in bananas. These findings collectively illustrated the mechanistic basis of a MaMADS1-MaNAC083-MaACS1/MaACOs regulatory cascade controlling ethylene biosynthesis during banana fruit ripening. These findings increase our knowledge of the transcriptional regulatory mechanisms of ethylene biosynthesis at the transcriptional level and are expected to help develop molecular approaches to control ripening and improve fruit storability.

MaBEL1 regulates banana fruit ripening by activating cell wall and starch degradation-related genes

Banana is a typical subtropical fruit, sensitive to chilling injuries and prone to softening disorder. However, the underlying regulatory mechanisms of the softening disorder caused by cold stress remain obscure. Herein, we found that BEL1-LIKE HOMEODOMAIN transcription factor 1 (MaBEL1) and its associated proteins regulate the fruit softening and ripening process. The transcript and protein levels of MaBEL1 were up-regulated with fruit ripening but severely repressed by the chilling stress. Moreover, the MaBEL1 protein interacted directly with the promoters of the cell wall and starch degradation-related genes, such as MaAMY3, MaXYL32, and MaEXP-A8. The transient overexpression of MaBEL1 alleviated fruit chilling injury and ripening disorder caused by cold stress and promoted fruit softening and ripening of "Fenjiao" banana by inducing ethylene production and starch and cell wall degradation. The accelerated ripening was also validated by the ectopic overexpression in tomatoes. Conversely, MaBEL1-silencing aggravated the chilling injury and ripening disorder and repressed fruit softening and ripening by inhibiting ethylene production and starch and cell wall degradation. MaABI5-like and MaEBF1, the two positive regulators of the fruit softening process, interacted with MaBEL1 to enhance the promoter activity of the starch and cell wall degradation-related genes. Moreover, the F-box protein MaEBF1 does not modulate the degradation of MaBEL1, which regulates the transcription of MaABI5-like protein. Overall, we report a novel MaBEL1-MaEBF1-MaABI5-like complex system that mediates the fruit softening and ripening disorder in "Fenjiao" bananas caused by cold stress.

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.

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.

Differential subgenome expression underlies biomass accumulation in allotetraploid Pennisetum giganteum

BACKGROUND

Pennisetum giganteum (AABB, 2n = 4x = 28) is a C4 plant in the genus Pennisetum with origin in Africa but currently also grown in Asia and America. It is a crucial forage and potential energy grass with significant advantages in yield, stress resistance, and environmental adaptation. However, the mechanisms underlying these advantageous traits remain largely unexplored. Here, we present a high-quality genome assembly of the allotetraploid P. giganteum aiming at providing insights into biomass accumulation.

RESULTS

Our assembly has a genome size 2.03 Gb and contig N50 of 88.47 Mb that was further divided into A and B subgenomes. Genome evolution analysis revealed the evolutionary relationships across the Panicoideae subfamily lineages and identified numerous genome rearrangements that had occurred in P. giganteum. Comparative genomic analysis showed functional differentiation between the subgenomes. Transcriptome analysis found no subgenome dominance at the overall gene expression level; however, differentially expressed homoeologous genes and homoeolog-specific expressed genes between the two subgenomes were identified, suggesting that complementary effects between the A and B subgenomes contributed to biomass accumulation of P. giganteum. Besides, C4 photosynthesis-related genes were significantly expanded in P. giganteum and their sequences and expression patterns were highly conserved between the two subgenomes, implying that both subgenomes contributed greatly and almost equally to the highly efficient C4 photosynthesis in P. giganteum. We also identified key candidate genes in the C4 photosynthesis pathway that showed sustained high expression across all developmental stages of P. giganteum.

CONCLUSIONS

Our study provides important genomic resources for elucidating the genetic basis of advantageous traits in polyploid species, and facilitates further functional genomics research and genetic improvement of P. giganteum.

Identification of annotation artifacts concerning the chalcone synthase (CHS)

OBJECTIVE

Chalcone synthase (CHS) catalyzes the initial step of the flavonoid biosynthesis. The CHS encoding gene is well studied in numerous plant species. Rapidly growing sequence databases contain hundreds of CHS entries that are the result of automatic annotation. In this study, we evaluated apparent multiplication of CHS domains in CHS gene models of four plant species.

MAIN FINDINGS

CHS genes with an apparent triplication of the CHS domain encoding part were discovered through database searches. Such genes were found in Macadamia integrifolia, Musa balbisiana, Musa troglodytarum, and Nymphaea colorata. A manual inspection of the CHS gene models in these four species with massive RNA-seq data suggests that these gene models are the result of artificial fusions in the annotation process. While there are hundreds of seemingly correct CHS records in the databases, it is not clear why these annotation artifacts appeared.

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.

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.

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.

Patterns of Chromosomal Variation, Homoeologous Exchange, and Their Relationship with Genomic Features in Early Generations of a Synthetic Rice Segmental Allotetraploid

Polyploidization is a driving force in plant evolution. Chromosomal variation often occurs at early generations following polyploid formation due to meiotic pairing irregularity that may compromise segregation fidelity and cause homoeologous exchange (HE). The trends of chromosomal variation and especially factors affecting HE remain to be fully deciphered. Here, by whole-genome resequencing, we performed nuanced analyses of patterns of chromosomal number variation and explored genomic features that affect HE in two early generations of a synthetic rice segmental allotetraploid. We found a wide occurrence of whole-chromosome aneuploidy and, to a lesser extent, also large segment gains/losses in both generations (S2 and S4) of the tetraploids. However, while the number of chromosome gains was similar between S2 and S4, that of losses in S4 was lower than in S2. HEs were abundant across all chromosomes in both generations and showed variable correlations with different genomic features at chromosomal and/or local scales. Contents of genes and transposable elements (TEs) were positively and negatively correlated with HE frequencies, respectively. By dissecting TEs into different classes, retrotransposons were found to be negatively correlated with HE frequency to a stronger extent than DNA transposons, whereas miniature terminal inverted elements (MITEs) showed a strong positive correlation. Local HE frequencies in the tetraploids and homologous recombination (HR) rates in diploids within 1 Mb sliding windows were significantly correlated with each other and showed similar overall distribution profiles. Nonetheless, non-concordant trends between HE and HR rates were found at distal regions in some chromosomes. At local scale, both shared and polymorphic retrotransposons between parents were negatively correlated with HE frequency; in contrast, both shared and polymorphic MITEs showed positive correlations with HE frequency. Our results shed new light on the patterns of chromosomal number variation and reveal genomic features influencing HE frequency in early generations following plant polyploidization.

Identification and analysis of lignin biosynthesis genes related to fruit ripening and stress response in banana (Musa acuminata L. AAA group, cv. Cavendish)

BACKGROUND

Lignin is a key component of the secondary cell wall of plants, providing mechanical support and facilitating water transport as well as having important impact effects in response to a variety of biological and abiotic stresses.

RESULTS

In this study, we identified 104 genes from ten enzyme gene families related to lignin biosynthesis in Musa acuminata genome and found the number of MaCOMT gene family was the largest, while MaC3Hs had only two members. MaPALs retained the original members, and the number of Ma4CLs in lignin biosynthesis was significantly less than that of flavonoids. Segmental duplication existed in most gene families, except for MaC3Hs, and tandem duplication was the main way to expand the number of MaCOMTs. Moreover, the expression profiles of lignin biosynthesis genes during fruit development, postharvest ripening stages and under various abiotic and biological stresses were investigated using available RNA-sequencing data to obtain fruit ripening and stress response candidate genes. Finally, a co-expression network of lignin biosynthesis genes was constructed by weighted gene co-expression network analysis to elucidate the lignin biosynthesis genes that might participate in lignin biosynthesis in banana during development and in response to stresses.

CONCLUSION

This study systematically identified the lignin biosynthesis genes in the Musa acuminata genome, providing important candidate genes for further functional analysis. The identification of the major genes involved in lignin biosynthesis in banana provides the basis for the development of strategies to improve new banana varieties tolerant to biological and abiotic stresses with high yield and high quality.

Spatiotemporal biocontrol and rhizosphere microbiome analysis of Fusarium wilt of banana

The soil-borne fungus Fusarium oxysporum f. sp. cubense tropical race 4 (Foc TR4) causes Fusarium wilt of banana (FWB), which devastates banana production worldwide. Biocontrol is considered to be the most efficient approach to reducing FWB. Here we introduce an approach that spatiotemporally applies Piriformospore indica and Streptomyces morookaensis strains according to their respective strength to increase biocontrol efficacy of FWB. P. indica successfully colonizes banana roots, promotes lateral root formation, inhibits Foc TR4 growth inside the banana plants and reduces FWB. S. morookaensis strain Sm4-1986 secretes different secondary compounds, of which xerucitrinin A (XcA) and 6-pentyl-α-pyrone (6-PP) show the strongest anti-Foc TR4 activity. XcA chelates iron, an essential nutrient in pathogen-plant interaction that determines the output of FWB. 6-PP, a volatile organic compound, inhibits Foc TR4 germination and promotes banana growth. Biocontrol trials in the field demonstrated that application of S. morookaensis lead to improvement of soil properties and increase of rhizosphere-associated microbes that are beneficial to banana growth, which significantly reduces disease incidence of FWB. Our study suggests that optimal utilization of the two biocontrol strains increases efficacy of biocontrol and that regulating iron accessibility in the rhizosphere is a promising strategy to control FWB.

Genome-Wide Characterization and Analysis of R2R3-MYB Genes Related to Fruit Ripening and Stress Response in Banana (Musa acuminata L. AAA Group, cv. 'Cavendish')

MYB is an important type of transcription factor in eukaryotes. It is widely involved in a variety of biological processes and plays a role in plant morphogenesis, growth and development, primary and secondary metabolite synthesis, and other life processes. In this study, bioinformatics methods were used to identify the R2R3-MYB transcription factor family members in the whole Musa acuminata (DH-Pahang) genome, one of the wild ancestors of banana. A total of 280 MaMYBs were obtained, and phylogenetic analysis indicated that these MaMYBs could be classified into 33 clades with MYBs from Arabidopsis thaliana. The amino acid sequences of the R2 and R3 Myb-DNA binding in all MaMYB protein sequences were quite conserved, especially Arg-12, Arg-13, Leu-23, and Leu-79. Distribution mapping results showed that 277 MaMYBs were localized on the 11 chromosomes in the Musa acuminata genome. The MaMYBs were distributed unevenly across the 11 chromosomes. More than 40.0% of the MaMYBs were located in collinear fragments, and segmental duplications likely played a key role in the expansion of the MaMYBs. Moreover, the expression profiles of MaMYBs in different fruit development and ripening stages and under various abiotic and biotic stresses were investigated using available RNA-sequencing data to obtain fruit development, ripening-specific, and stress-responsive candidate genes. Weighted gene co-expression network analysis (WGCNA) was used to analyze transcriptome data of banana from the above 11 samples. We found MaMYBs participating in important metabolic biosynthesis pathways in banana. Collectively, our results represent a comprehensive genome-wide study of the MaMYB gene family, which should be helpful in further detailed studies on MaMYBs functions related to fruit development, postharvest ripening, and the seedling response to stress in an important banana cultivar.

Genome assembly of Musa beccarii shows extensive chromosomal rearrangements and genome expansion during evolution of Musaceae genomes

BACKGROUND

Musa beccarii (Musaceae) is a banana species native to Borneo, sometimes grown as an ornamental plant. The basic chromosome number of Musa species is x = 7, 10, or 11; however, M. beccarii has a basic chromosome number of x = 9 (2n = 2x = 18), which is the same basic chromosome number of species in the sister genera Ensete and Musella. Musa beccarii is in the section Callimusa, which is sister to the section Musa. We generated a high-quality chromosome-scale genome assembly of M. beccarii to better understand the evolution and diversity of genomes within the family Musaceae.

FINDINGS

The M. beccarii genome was assembled by long-read and Hi-C sequencing, and genes were annotated using both long Iso-seq and short RNA-seq reads. The size of M. beccarii was the largest among all known Musaceae assemblies (∼570 Mbp) due to the expansion of transposable elements and increased 45S ribosomal DNA sites. By synteny analysis, we detected extensive genome-wide chromosome fusions and fissions between M. beccarii and the other Musa and Ensete species, far beyond those expected from differences in chromosome number. Within Musaceae, M. beccarii showed a reduced number of terpenoid synthase genes, which are related to chemical defense, and enrichment in lipid metabolism genes linked to the physical defense of the cell wall. Furthermore, type III polyketide synthase was the most abundant biosynthetic gene cluster (BGC) in M. beccarii. BGCs were not conserved in Musaceae genomes.

CONCLUSIONS

The genome assembly of M. beccarii is the first chromosome-scale genome assembly in the Callimusa section in Musa, which provides an important genetic resource that aids our understanding of the evolution of Musaceae genomes and enhances our knowledge of the pangenome.

Ethylene Response Factor MaERF012 Modulates Fruit Ripening by Regulating Chlorophyll Degradation and Softening in Banana

Ethylene response factors (ERFs) are one of largest plant-specific transcription factor families involved in fruit ripening. However, the regulatory mechanism by which ERFs modulate fruit yellowing and softening remains unknown in banana. We previously found that the transcription of MaERF012 was closely related to 'Fenjiao' banana fruit ripening. Herein, we found that MaERF012 was differentially expressed in the fruit pulp and peel and was closely related to fruit ripening. MaERF012 activated the promoter activity of one chlorophyll degradation gene (MaSGR1), two starch degradation genes (MaGWD1 and MaAMY3), and three cell wall degradation genes (MaPL8, MaEXP-A8, and MaXYL23-like), which were tested by EMSA, Y1H, and DLR. Transient overexpression of MaERF012 accelerates fruit ripening by promoting fruit yellowing and softening by up-regulating the transcription of chlorophyll, starch, and cell wall degradation genes. Over-expression of MaERF012 alters the transcriptome profiles of the fruit peel and pulp, and the differentially expressed genes were mainly enriched in starch and sucrose metabolism, plant hormone signal transduction, biosynthesis of secondary metabolism, and fructose and mannose metabolism. Overall, the data showed that MaERF012 acts as a transcriptional activator by regulating fruit ripening by activating the transcription of chlorophyll, starch, and cell wall degradation genes.

Elucidating the role of MaBAM9b in starch degradation

Banana is a typical starch conversion fruit. The high content of starch at harvest is quickly digested and converted to soluble sugars during the postharvest ripening process, ultimately contributing to fruit flavor. This process is regulated in a complex manner by genes and environmental factors. MaBAM9b is one of the main enzyme genes previously found by transcriptomic analysis to be highly expressed in banana fruit. However, its exact role in starch degradation remains unclear. Here, full-length MaBAM9b was isolated from banana fruit, and its subcellular localization, protein expression, and transient expression in banana fruit slices were investigated. In addition, sense and anti-sense MaBAM9b were transformed into rice (Oryza sativa L. japonica. cv. 'Nipponbare') to identify the function of MaBAM9b. MaBAM9b was 1599 bp and encoded 532 amino acids. It contained two conserved domains of PLN02803 and glycosyl hydrolase family 14 and was localized in the chloroplast. The protein expression pattern of MaBAM9b remained consistently high throughout banana fruit ripening and starch degradation. Transient overexpression or inhibition of MaBAM9b in banana fruit greatly improved or suppressed starch degradation. Genetic modification of rice indicated that overexpression of MaBAM9b greatly improved starch degradation and seed germination, while inhibition of its expression suppressed these biological processes. These results support the key role of MaBAM9b in starch degradation and provide a target gene for banana fruit quality improvement and biological breeding.

Insights to the superoxide dismutase genes and its roles in Hevea brasiliensis under abiotic stress

The superoxide dismutase (SOD) protein significantly influences the development and growth of plants and their reaction to abiotic stresses. However, little is known about the characteristics of rubber tree SOD genes and their expression changes under abiotic stresses. The present study recognized 11 SOD genes in the rubber tree genome, including 7 Cu/ZnSODs, 2 MnSODs, and 2 FeSODs. Except for HbFSD1, SODs were scattered on five chromosomes. The phylogenetic analysis of SOD proteins in rubber trees and a few other plants demonstrated that the SOD proteins contained three major subgroups. Moreover, the genes belonging to the same clade contained similar gene structures, which confirmed their classification further. The extension of the SOD gene family in the rubber tree was mainly induced by the segmental duplication events. The cis-acting components analysis showed that HbSODs were utilized in many biological procedures. The transcriptomics data indicated that the phosphorylation of the C-terminal domain of RNA polymerase II might control the cold response genes through the CBF pathway and activate the SOD system to respond to cold stress. The qRT-PCR results showed that the expression of HbCSD1 was significantly downregulated under drought and salt stresses, which might dominate the adaption capability to different stresses. Additionally, salt promoted the expression levels of HbMSD1 and HbMSD2, exhibiting their indispensable role in the salinity reaction. The study results will provide a theoretical basis for deep research on HbSODs in rubber trees.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03328-7.

Chromosome-scale genome assemblies and annotations for Poales species Carex cristatella, Carex scoparia, Juncus effusus, and Juncus inflexus

The majority of sequenced genomes in the monocots are from species belonging to Poaceae, which include many commercially important crops. Here, we expand the number of sequenced genomes from the monocots to include the genomes of 4 related cyperids: Carex cristatella and Carex scoparia from Cyperaceae and Juncus effusus and Juncus inflexus from Juncaceae. The high-quality, chromosome-scale genome sequences from these 4 cyperids were assembled by combining whole-genome shotgun sequencing of Nanopore long reads, Illumina short reads, and Hi-C sequencing data. Some members of the Cyperaceae and Juncaceae are known to possess holocentric chromosomes. We examined the repeat landscapes in our sequenced genomes to search for potential repeats associated with centromeres. Several large satellite repeat families, comprising 3.2–9.5% of our sequenced genomes, showed dispersed distribution of large satellite repeat clusters across all Carex chromosomes, with few instances of these repeats clustering in the same chromosomal regions. In contrast, most large Juncus satellite repeats were clustered in a single location on each chromosome, with sporadic instances of large satellite repeats throughout the Juncus genomes. Recognizable transposable elements account for about 20% of each of the 4 genome assemblies, with the Carex genomes containing more DNA transposons than retrotransposons while the converse is true for the Juncus genomes. These genome sequences and annotations will facilitate better comparative analysis within monocots.

The banana genome hub: a community database for genomics in the Musaceae

The Banana Genome Hub provides centralized access for genome assemblies, annotations, and the extensive related omics resources available for bananas and banana relatives. A series of tools and unique interfaces are implemented to harness the potential of genomics in bananas, leveraging the power of comparative analysis, while recognizing the differences between datasets. Besides effective genomic tools like BLAST and the JBrowse genome browser, additional interfaces enable advanced gene search and gene family analyses including multiple alignments and phylogenies. A synteny viewer enables the comparison of genome structures between chromosome-scale assemblies. Interfaces for differential expression analyses, metabolic pathways and GO enrichment were also added. A catalogue of variants spanning the banana diversity is made available for exploration, filtering, and export to a wide variety of software. Furthermore, we implemented new ways to graphically explore gene presence-absence in pangenomes as well as genome ancestry mosaics for cultivated bananas. Besides, to guide the community in future sequencing efforts, we provide recommendations for nomenclature of locus tags and a curated list of public genomic resources (assemblies, resequencing, high density genotyping) and upcoming resources-planned, ongoing or not yet public. The Banana Genome Hub aims at supporting the banana scientific community for basic, translational, and applied research and can be accessed at https://banana-genome-hub.southgreen.fr.

The genetic mechanisms underlying the convergent evolution of pollination syndromes in the Neotropical radiation of Costus L

Selection together with variation in floral traits can act to mold floral form, often driven by a plant's predominant or most effective pollinators. To investigate the evolution of traits associated with pollination, we developed a phylogenetic framework for evaluating tempo and mode of pollination shifts across the genus Costus L., known for its evolutionary toggle between traits related to bee and bird pollination. Using a target enrichment approach, we obtained 957 loci for 171 accessions to expand the phylogenetic sampling of Neotropical Costus. In addition, we performed whole genome resequencing for a subset of 20 closely related species with contrasting pollination syndromes. For each of these 20 genomes, a high-quality assembled transcriptome was used as reference for consensus calling of candidate loci hypothesized to be associated with pollination-related traits of interest. To test for the role these candidate genes may play in evolutionary shifts in pollinators, signatures of selection were estimated as dN/dS across the identified candidate loci. We obtained a well-resolved phylogeny for Neotropical Costus despite conflict among gene trees that provide evidence of incomplete lineage sorting and/or reticulation. The overall topology and the network of genome-wide single nucleotide polymorphisms (SNPs) indicate that multiple shifts in pollination strategy have occurred across Costus, while also suggesting the presence of previously undetected signatures of hybridization between distantly related taxa. Traits related to pollination syndromes are strongly correlated and have been gained and lost in concert several times throughout the evolution of the genus. The presence of bract appendages is correlated with two traits associated with defenses against herbivory. Although labellum shape is strongly correlated with overall pollination syndrome, we found no significant impact of labellum shape on diversification rates. Evidence suggests an interplay of pollination success with other selective pressures shaping the evolution of the Costus inflorescence. Although most of the loci used for phylogenetic inference appear to be under purifying selection, many candidate genes associated with functional traits show evidence of being under positive selection. Together these results indicate an interplay of phylogenetic history with adaptive evolution leading to the diversification of pollination-associated traits in Neotropical Costus.

Effects of homoeologous exchange on gene expression and alternative splicing in a newly formed allotetraploid wheat

Homoeologous exchange (HE) is a major mechanism generating post-polyploidization genetic variation with important evolutionary consequences. However, the direct impacts of HE on gene expression and transcript diversity in allopolyploids without the intertwined evolutionary processes remain to be fully understood. Here, we analyzed high-throughput RNA-seq data of young leaves from plant groups of a synthetic allotetraploid wheat (AADD), which contained variable numbers of HEs. We aimed to investigate if and to which extent HE directly impacts gene expression and alternative splicing (AS). We found that HE impacts expression of genes located within HE regions primarily via a cis-acting dosage effect, which led to significant changes in the total expression level of homoeologous gene pairs, especially for homoeologs whose original expression was biased. In parallel, HE also influences expression of a large number of genes residing in non-HE regions by trans-regulation leading to convergent expression of homoeologs. Intriguingly, when taking the original relative homoeolog expression states into account, homoeolog pairs under trans-effect are more prone to manifesting a convergent response to the HEs whereas those under cis-regulation tended to show further exacerbated subgenome-biased expression. Moreover, HE-induced quantitative, largely individual-specific, changes of AS events were detected. Similar to homoeologous expression, homoeo-AS events under trans-effect were more responsive to HE. HE therefore exerts multifaceted immediate effects on gene expression and, to a less extent, on individualized transcript diversity in nascent allopolyploidy.

Genome-wide identification, evolution analysis of LysM gene family members and their expression analysis in response to biotic and abiotic stresses in banana (Musa L.)

LysM (Lysin motif), in response to pathogenic molecular stresses, is a crucial signal recognition gene. To understand the molecular characteristics of banana LysM gene family members, we used a series of bioinformatics methods. Based on the genomic databases of Musa acuminata, Musa balbisiana and Musa itinerans, a total of 53 genes and 55 proteins were identified, with 21 genes and 23 proteins in the M.acuminata, 16 genes and 16 proteins in each of M.balbisiana and M.itinerans, respectively. According to the conserved structural domains, LysM can be divided into five classes, namely LysM&MltD, LYK, LYP, LysMn, and LysMe. The LysM gene was relatively highly conserved in the evolution of the three genomes of banana, and some differences occurred. Expression analysis revealed that MaLysM4-5 was relatively highly expressed under high-temperature stress, low-temperature stress and pathogen infection; at the same time, about one-third of the members were down-regulated under low-temperature stress and high-temperature stress, while the expression of MaLysM10-1 and MaLysM4-5 were up-regulated. After the banana wilt fungus FocTR4 infected the banana roots, MaLysM1 was down-regulated and MaLysM11-1 was up-regulated. In conclusion, our study suggests that MaLysMs may be necessary in the response to high- and low-temperature stresses, as well as the banana wilt fungus infestation. Overall, this paper found that LysM genes may be involved in biotic and abiotic stresses in banana, and provided helpful information about LysM's evolution, expression and properties, which will provide theoretical references for further studies on the functions of LysM genes and resistance breeding in the future.

The Musa troglodytarum L. genome provides insights into the mechanism of non-climacteric behaviour and enrichment of carotenoids

BACKGROUND

Karat (Musa troglodytarum L.) is an autotriploid Fe'i banana of the Australimusa section. Karat was domesticated independently in the Pacific region, and karat fruit are characterized by a pink sap, a deep yellow-orange flesh colour, and an abundance of β-carotene. Karat fruit showed non-climacteric behaviour, with an approximately 215-day bunch filling time. These features make karat a valuable genetic resource for studying the mechanisms underlying fruit development and ripening and carotenoid biosynthesis.

RESULTS

Here, we report the genome of M. troglodytarum, which has a total length of 603 Mb and contains 37,577 predicted protein-coding genes. After divergence from the most recent common ancestors, M. troglodytarum (T genome) has experienced fusion of ancestral chromosomes 8 and 9 and multiple translocations and inversions, unlike the high synteny with few rearrangements found among M. schizocarpa (S genome), M. acuminata (A genome) and M. balbisiana (B genome). Genome microsynteny analysis showed that the triplication of MtSSUIIs due to chromosome rearrangement may lead to the accumulation of carotenoids and ABA in the fruit. The expression of duplicated MtCCD4s is repressed during ripening, leading to the accumulation of α-carotene, β-carotene and phytoene. Due to a long terminal repeat (LTR)-like fragment insertion upstream of MtERF11, karat cannot produce large amounts of ethylene but can produce ABA during ripening. These lead to non-climacteric behaviour and prolonged shelf-life, which contributes to an enrichment of carotenoids and riboflavin.

CONCLUSIONS

The high-quality genome of M. troglodytarum revealed the genomic basis of non-climacteric behaviour and enrichment of carotenoids, riboflavin, flavonoids and free galactose and provides valuable resources for further research on banana domestication and breeding and the improvement of nutritional and bioactive qualities.

Chromosomal instability and phenotypic variation in a specific lineage derived from a synthetic allotetraploid wheat

Newly formed plant allopolyploids usually have meiosis defect, resulting in chromosomal instability manifested as variation in chromosome number and/or structure. However, not all nascent allopolyploids are equally unstable. The wheat group (Aegilops/Triticum) contains 13 diploid species with distinct genome types. Many of these species can be artificially hybridized to produce viable but sterile inter-specific/intergeneric F1 hybrids, which can generate fertile synthetic allotetraploid wheats after whole genome doubling. Compared with synthetic allotetraploid wheats that contain genome combinations of AADD and S*S*DD (S* refers to related S genomes of a different species), those containing an S*S*AA genome are significantly more stable. However, robustness of the relative stability of S*S*AA genomes is unknown, nor are the phenotypic and fitness consequences during occurrences of secondary chromosomal instability. Here, we report a specific lineage originated from a single individual plant of a relatively stable synthetic allotetraploid wheat with genomes S ^l S ^l AA (S ^l and A subgenomes were from Ae. longissima and T. urartu, respectively) that showed a high degree of transgenerational chromosomal instability. Both numerical chromosome variation (NCV) and structural chromosome variation (SCV) occurred widely. While substantial differences in frequencies of both NCV and SCV were detected across the different chromosomes, only NCV frequencies were significantly different between the two subgenomes. We found that NCVs and SCVs occurred primarily due to perturbed meiosis, allowing formation of multivalents and univalents as well as homoeologous exchanges. Thus, the combination of NCVs and SCVs affected multiple phenotypic traits, particularly those related to reproductive fitness.

From fruit growth to ripening in plantain: a careful balance between carbohydrate synthesis and breakdown

In this study, we aimed to investigate for the first time different fruit development stages in plantain banana in order gain insights into the order of appearance and dominance of specific enzymes and fluxes. We examined fruit development in two plantain banana cultivars during the period between 2-12 weeks after bunch emergence using high-throughput proteomics, quantification of major metabolites, and analyses of metabolic fluxes. Starch synthesis and breakdown are processes that take place simultaneously. During the first 10 weeks fruits accumulated up to 48% of their dry weight as starch, and glucose 6-phosphate and fructose were important precursors. We found a unique amyloplast transporter and hypothesize that it facilitates the import of fructose. We identified an invertase originating from the Musa balbisiana genome that would enable carbon flow back to growth and starch synthesis and maintain a high starch content even during ripening. Enzymes associated with the initiation of ripening were involved in ethylene and auxin metabolism, starch breakdown, pulp softening, and ascorbate biosynthesis. The initiation of ripening was cultivar specific, with faster initiation being particularly linked to the 1-aminocyclopropane-1-carboxylate oxidase and 4-alpha glucanotransferase disproportionating enzymes. Information of this kind is fundamental to determining the optimal time for picking the fruit in order to reduce post-harvest losses, and has potential applications for breeding to improve fruit quality.

Genomic and Bioinformatic Resources for Perennial Fruit Species

In the post-genomic era, data management and development of bioinformatic tools are critical for the adequate exploitation of genomics data. In this review, we address the actual situation for the subset of crops represented by the perennial fruit species. The agronomical singularity of these species compared to plant and crop model species provides significant challenges on the implementation of good practices generally not addressed in other species. Studies are usually performed over several years in non-controlled environments, usage of rootstock is common, and breeders heavily rely on vegetative propagation. A reference genome is now available for all the major species as well as many members of the economically important genera for breeding purposes. Development of pangenome for these species is beginning to gain momentum which will require a substantial effort in term of bioinformatic tool development. The available tools for genome annotation and functional analysis will also be presented.

Assessing Loss of Regulatory Divergence, Genome-Transcriptome Incongruence, and Preferential Expression Switching in Abaca × Banana Backcrosses

The Musa textilis var. Abuab has high fiber quality (FQ) but is susceptible to abaca bunchy top virus (AbBTV); the Musa balbisiana var. Pacol has low FQ but is resistant against AbBTV. Their backcrosses (BC₂ and BC₃) possess both desirable traits. Analysis using RNA-seq showed that the regulatory divergence of Abuab and Pacol is largely explained by cis differences with 27.4% and 22.3% if we are to assess it using BC₂ and BC₃, respectively. Cis differences between the two genotypes are significantly reduced from BC₂ to BC₃ due to changes in genomic constitution. Trans, on the other hand, is robust to changes in allelic composition. All these are attributed to the loss of heterozygosity in BC₃ relative to BC₂. Further analysis showed that both backcrosses exhibited genome-wide preferential expression of Pacol- over Abuab-specific alleles, despite the wider genetic presence of the latter in the hybrids. The ratio of the two genotype-specific expressed transcripts and the ratio of their corresponding genetic make-up are significantly disproportionate, a phenomenon that we refer to here as "genome-transcriptome incongruence". We also observed preferential expression switching in which several genes prefer the Abuab- (or Pacol-) specific allele in BC₂ but switched to the Pacol- (or Abuab-) specific allele in the BC₃ genome.

High-quality reference genome sequences of two Cannaceae species provide insights into the evolution of Cannaceae

Canna edulis Ker-Gawl and Canna indica L. are species belonging to the Cannaceae family and both have a very high economic value. Here, we aimed to assemble genomes of C. edulis and C. indica at the chromosome level to generate a reference genome for the Cannaceae family. We also comparatively analyzed the genomes of C. edulis and C. indica and examined the molecular mechanisms responsible for the remarkable differences in plant characteristics in C. edulis varieties. Our results indicated that genome-wide duplication events had recently occurred in C. edulis and C. indica. The comparative analysis of the genomes of C. edulis and C. indica revealed that C. edulis exhibited a remarkable level of replication of genes in the starch and sucrose metabolic pathways, especially during sucrose hydrolysis. This finding is consistent with the fact that the starch content of the C. edulis tuber is higher than that of C. indica. Simplified genome re-sequencing revealed the population structure of 241 C. edulis genes, and a genome-wide association study of leaf traits revealed the location of key genes related to leaf color and morphology. These findings extend our understanding of Cannaceade at the molecular level, and provide an effective theoretical basis for further study and utilization of Cannaceae plants.

Genome-wide analysis of pathogenesis-related protein 1 (PR-1) gene family from Musa spp. and its role in defense response during stresses

Pathogenesis related protein-1 (PR-1) is the most abundantly produced protein during defense response against many biotic and abiotic stresses. However, knowledge on PR-1 gene family and its evolutionary relationship in banana is very limited. In order to study the potential role of PR-1 genes in banana, genome wide identification, structure analysis and expressions were performed. A total of 15 and 11 PR-1 genes were identified from A and B genomes of banana and the proteins encoded by this gene family are of varying lengths and harbor conserved domains and motifs. PR-1 genes are unevenly dispersed on 11 chromosomes with segmental duplication in both A and B genome, suggesting an important contribution of duplication in expansion of PR-1 gene family in banana. qRT-PCR analysis of PR-1 gene showed positive correlation with the RNAseq data under various stresses and examination of expression pattern of selected MaPR-1 genes in banana revealed its role in biotic and abiotic stresses in general and fusarium wilt in particular. This study provides significant insight into the functions of PR-1 genes which can be further exploited as a promising candidate for developing multiple stress tolerant banana varieties.

A chromosome-level reference genome of Ensete glaucum gives insight into diversity and chromosomal and repetitive sequence evolution in the Musaceae

BACKGROUND

Ensete glaucum (2n = 2x = 18) is a giant herbaceous monocotyledonous plant in the small Musaceae family along with banana (Musa). A high-quality reference genome sequence assembly of E. glaucum is a resource for functional and evolutionary studies of Ensete, Musaceae, and the Zingiberales.

FINDINGS

Using Oxford Nanopore Technologies, chromosome conformation capture (Hi-C), Illumina and RNA survey sequence, supported by molecular cytogenetics, we report a high-quality 481.5 Mb genome assembly with 9 pseudo-chromosomes and 36,836 genes. A total of 55% of the genome is composed of repetitive sequences with predominantly LTR-retroelements (37%) and DNA transposons (7%). The single 5S ribosomal DNA locus had an exceptionally long monomer length of 1,056 bp, more than twice that of the monomers at multiple loci in Musa. A tandemly repeated satellite (1.1% of the genome, with no similar sequence in Musa) was present around all centromeres, together with a few copies of a long interspersed nuclear element (LINE) retroelement. The assembly enabled us to characterize in detail the chromosomal rearrangements occurring between E. glaucum and the x = 11 species of Musa. One E. glaucum chromosome has the same gene content as Musa acuminata, while others show multiple, complex, but clearly defined evolutionary rearrangements in the change between x= 9 and 11.

CONCLUSIONS

The advance towards a Musaceae pangenome including E. glaucum, tolerant of extreme environments, makes a complete set of gene alleles, copy number variation, and a reference for structural variation available for crop breeding and understanding environmental responses. The chromosome-scale genome assembly shows the nature of chromosomal fusion and translocation events during speciation, and features of rapid repetitive DNA change in terms of copy number, sequence, and genomic location, critical to understanding its role in diversity and evolution.

Musa Germplasm A and B Genomic Composition Differentially Affects Their Susceptibility to Banana Bunchy Top Virus and Its Aphid Vector, Pentalonia nigronervosa

Banana bunchy top disease (BBTD), caused by the banana bunchy top virus (BBTV, genus Babuvirus), is the most destructive viral disease of banana and plantain (Musa spp.). The virus is transmitted persistently by the banana aphid, Pentalonia nigronervosa Coquerel (Hemiptera: Aphididae). While research efforts have focused on screening Musa genotypes for BBTD resistance, comparatively little work has been carried out to identify resistance to banana aphids. This study assessed 44 Musa germplasm of different A and B genome composition for the performance of banana aphids under semicontrolled environmental screenhouse conditions and in a field trial established in a BBTD endemic location. In the screenhouse, the AA diploid Calcutta 4 had the lowest apterous aphid density per plant (9.7 ± 4.6) compared with AAB triploid Waema, which had the highest aphid densities (395.6 ± 20.8). In the field, the highest apterous aphid density per plant (29.2 ± 6.7) occurred on the AAB triploid Batard and the lowest (0.4 ± 0.2) on the AA diploid Pisang Tongat. The AA diploid Tapo was highly susceptible to BBTD (100% infection) compared with the genotypes Balonkawe (ABB), PITA 21 (AAB), Calcutta 4 (AA), and Balbisiana Los Banos (BB), which remained uninfected. The Musa genotypes with apparent resistance to BBTD and least susceptibility to aphid population growth provide options for considering aphid and BBTD resistance in banana and plantain breeding programs.

Genetic Characteristics and Metabolic Interactions between Pseudocercospora fijiensis and Banana: Progress toward Controlling Black Sigatoka

The international importance of banana and severity of black Sigatoka disease have led to extensive investigations into the genetic characteristics and metabolic interactions between the Dothideomycete Pseudocercospora fijiensis and its banana host. P. fijiensis was shown to have a greatly expanded genome compared to other Dothideomycetes, due to the proliferation of retrotransposons. Genome analysis suggests the presence of dispensable chromosomes that may aid in fungal adaptation as well as pathogenicity. Genomic research has led to the characterization of genes and metabolic pathways involved in pathogenicity, including: secondary metabolism genes such as PKS10-2, genes for mitogen-activated protein kinases such as Fus3 and Slt2, and genes for cell wall proteins such as glucosyl phosphatidylinositol (GPI) and glycophospholipid surface (Gas) proteins. Studies conducted on resistance mechanisms in banana have documented the role of jasmonic acid and ethylene pathways. With the development of banana transformation protocols, strategies for engineering resistance include transgenes expressing antimicrobial peptides or hydrolytic enzymes as well as host-induced gene silencing (HIGS) targeting pathogenicity genes. Pseudocercospora fijiensis has been identified as having high evolutionary potential, given its large genome size, ability to reproduce both sexually and asexually, and long-distance spore dispersal. Thus, multiple control measures are needed for the sustainable control of black Sigatoka disease.

RNA-Seq Reveals Differentially Expressed Genes Associated with High Fiber Quality in Abaca (Musa textilis Nee)

Despite the importance of and current demand for abaca (Musa textilis Nee) fiber, there has been limited study that capitalizes on RNA-seq to identify candidate genes associated with high fiber quality and bunchy top virus (AbBTV) resistance. Three varieties (Abuab, Inosa, and Tangongon), one wild banana variety (Musa balbisiana Colla) Pacol, and two developed backcrosses (Abuab × Pacol BC₂ and BC₃) were grown at the Institute of Plant Breeding (IPB), Laguna, Philippines. The pseudostems of 3-month-old suckers of each genotype were sampled for RNA-seq. Datasets were analyzed for differential expression (DE) implementing various model frameworks, including pairwise, genotypic and non-DE models. Results indicate that Abuab and BC₃ induce the highest proportion (70%) of abaca-specific genes. Gene ontology (GO) enrichment analysis showed several genes associated with cellulose synthase activity, callose synthase, ß-glucosidase activity, glucan biosynthetic process, etc. KEGG pathway analysis showed several genes encoding for enzymes involved in the lignin biosynthetic pathway. Analysis using genotypic DE (GDE) between abaca bunchy top virus (AbBTV)-resistant and -susceptible groups revealed genes such as pathogenesis-related protein and NBS-LRR. As the genotypes were not infected with the pathogen, these genes are yet to be confirmed for their roles in disease resistance and are an interesting subject for further investigation.

The pangenome of banana highlights differences between genera and genomes

Banana (Musaceae family) has a complex genetic history and includes a genus Musa with a variety of cultivated clones with edible fruits, Ensete species that are grown for their edible corm, and monospecific Musella whose generic status has been questioned. The most commonly exported banana cultivars belong to Cavendish, a subgroup of Musa triploid cultivars, which is under threat by fungal pathogens, though there are also related species M. balbisiana Colla (B genome), M. textilis Née (T genome), and M. schizocarpa N. W. Simmonds (S genome), along with hybrids of these genomes, which potentially host genes of agronomic interest. Here we present the first cross-genus pangenome of banana, which contains representatives of the Musa and Ensete genera. Clusters based on gene presence-absence variation (PAV) clearly separate Musa and Ensete, while Musa is split further based on species. These results present the first pangenome study across genus boundaries and identifies genes that differentiate between Musaceae species, information that may support breeding programs in these crops.