Unravelling the complex story of intergenomic recombination in ABB allotriploid bananas

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.

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.

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.

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.

Interploidy Introgression Shaped Adaptation during the Origin and Domestication History of Brassica napus

Polyploidy is recurrent across the tree of life and known as an evolutionary driving force in plant diversification and crop domestication. How polyploid plants adapt to various habitats has been a fundamental question that remained largely unanswered. Brassica napus is a major crop cultivated worldwide, resulting from allopolyploidy between unknown accessions of diploid B. rapa and B. oleracea. Here, we used whole-genome resequencing data of accessions representing the majority of morphotypes and ecotypes from the species B. rapa, B. oleracea, and B. napus to investigate the role of polyploidy during domestication. To do so, we first reconstructed the phylogenetic history of B. napus, which supported the hypothesis that the emergence of B. napus derived from the hybridization of European turnip of B. rapa and wild B. oleracea. These analyses also showed that morphotypes of swede and Siberian kale (used as vegetable and fodder) were domesticated before rapeseed (oil crop). We next observed that frequent interploidy introgressions from sympatric diploids were prominent throughout the domestication history of B. napus. Introgressed genomic regions were shown to increase the overall genetic diversity and tend to be localized in regions of high recombination. We detected numerous candidate adaptive introgressed regions and found evidence that some of the genes in these regions contributed to phenotypic diversification and adaptation of different morphotypes. Overall, our results shed light on the origin and domestication of B. napus and demonstrate interploidy introgression as an important mechanism that fuels rapid diversification in polyploid species.

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.

Shared pedigree relationships and transmission of unreduced gametes in cultivated banana

BACKGROUND AND AIMS

Cultivated bananas resulted from inter(sub)specific hybridizations involving Musa species and subspecies (M. acuminata subspecies, M. schizocarpa, M. balbisiana) and the subsequent selection, centuries ago, of hybrids with parthenocarpic, seedless fruits. Cultivars have low fertility and are vegetatively propagated, forming groups of somaclones. Relatively few of them, mainly triploids, are grown on a large scale and characterization of their parental relationships may be useful for breeding strategies. Here we investigate parental relationships and gamete-type contributions among diploid and polyploid banana cultivars.

METHODS

We used SNP genotyping data from whole-genome sequencing of 178 banana individuals, including 111 cultivars, 55 wild bananas and 12 synthetic F1 hybrids. We analysed the proportion of SNP sites in accordance with direct parentage with a global statistic and along chromosomes for selected individuals.

KEY RESULTS

We characterized parentage relationships for 7 diploid cultivars, 11 triploid cultivars and 1 tetraploid cultivar. Results showed that both diploid and triploid cultivars could have contributed gametes to other banana cultivars. Diploids may have contributed 1x or 2x gametes and triploids 1x to 3x gametes. The Mchare diploid cultivar group, nowadays only found in East Africa, was found as parent of two diploid and eight triploid cultivars. In five of its identified triploid offspring, corresponding to main export or locally popular dessert bananas, Mchare contributed a 2x gamete with full genome restitution without recombination. Analyses of remaining haplotypes in these Mchare offspring suggested ancestral pedigree relationships between different interspecific banana cultivars.

CONCLUSIONS

The current cultivated banana resulted from different pathways of formation, with implication of recombined or un-recombined unreduced gametes produced by diploid or triploid cultivars. Identification of dessert banana's parents and the types of gametes they contributed should support the design of breeding strategies.

Unravelling the diversity in water usage among wild banana species in response to vapour pressure deficit

The rise in global temperature is not only affecting plant functioning directly, but is also increasing air vapour pressure deficit (VPD). The yield of banana is heavily affected by water deficit but so far breeding programs have never addressed the issue of water deficit caused by high VPD. A reduction in transpiration at high VPD has been suggested as a key drought tolerance breeding trait to avoid excessive water loss, hydraulic failure and to increase water use efficiency. In this study, stomatal and transpiration responses under increasing VPD at the leaf and whole-plant level of 8 wild banana (sub)species were evaluated, displaying significant differences in stomatal reactivity. Three different phenotypic groups were identified under increasing VPD. While (sub)species of group III maintained high transpiration rates under increasing VPD, M. acuminata ssp. errans (group I), M. acuminata ssp. zebrina (group II) and M. balbisiana (group II) showed the highest transpiration rate limitations to increasing VPD. In contrast to group I, group II only showed strong reductions at high VPD levels, limiting the cost of reduced photosynthesis and strongly increasing their water use efficiency. M. acuminata ssp. zebrina and M. balbisiana thus show the most favourable responses. This study provides a basis for the identification of potential parent material in gene banks for breeding future-proof bananas that cope better with lack of water.

Patterns, mechanisms, and consequences of homoeologous exchange in allopolyploid angiosperms: a genomic and epigenomic perspective

Allopolyploids result from hybridization between different evolutionary lineages coupled with genome doubling. Homoeologous chromosomes (chromosomes with common shared ancestry) may undergo recombination immediately after allopolyploid formation and continue over successive generations. The outcome of this meiotic pairing behavior is dynamic and complex. Homoeologous exchanges (HEs) may lead to the formation of unbalanced gametes, reduced fertility, and selective disadvantage. By contrast, HEs could act as sources of novel evolutionary substrates, shifting the relative dosage of parental gene copies, generating novel phenotypic diversity, and helping the establishment of neo-allopolyploids. However, HE patterns vary among lineages, across generations, and even within individual genomes and chromosomes. The causes and consequences of this variation are not fully understood, though interest in this evolutionary phenomenon has increased in the last decade. Recent technological advances show promise in uncovering the mechanistic basis of HEs. Here, we describe recent observations of the common patterns among allopolyploid angiosperm lineages, underlying genomic and epigenomic features, and consequences of HEs. We identify critical research gaps and discuss future directions with far-reaching implications in understanding allopolyploid evolution and applying them to the development of important phenotypic traits of polyploid crops.

Analysis of genetic diversity and agronomic variation in banana sub-populations for genomic selection under drought stress in southern Benin

In the perspective of investigating genomic selection (GS) among Musa genotypes in West and Central Africa, banana accessions were phenotyped under natural drought stress in Benin and genotyped using genotyping by sequencing. Sixty-one (61) accessions grouped into three major genomic groups AAA, AAB and ABB and those without genomic affiliation information were used. Variation within the population was determined by phenotypic variables while population structure and clustering analysis were carried out to understand the genetic diversity at the molecular level. Among the genomic groups evaluated, the group AAB showed the best performance for fruit weight at maturity, (3.41 ± 1.99 kg) and for plant height (198.46 ± 12.66 cm). At the accession level, HD 117 S1 and NIA 27 showed the best plant height (263.16 ± 20.98 cm) and the best fruit weight at maturity (9.43 ± 0.0 kg) respectively. Phenotypic data did not reveal clear genetic diversity among accessions; however, the genetic diversity was conspicuous at the molecular level using 5000 markers. The affiliations of local accessions in genomic groups were determined for the first time based on the phenotypic and molecular data obtained in this study. The knowledge generated allows the possibility to apply GS in banana.

Interspecific introgression patterns reveal the origins of worldwide cultivated bananas in New Guinea

Hybridizations between Musa species and subspecies, enabled by their transport via human migration, were proposed to have played an important role in banana domestication. We exploited sequencing data of 226 Musaceae accessions, including wild and cultivated accessions, to characterize the inter(sub)specific hybridization pattern that gave rise to cultivated bananas. We identified 11 genetic pools that contributed to cultivars, including two contributors of unknown origin. Informative alleles for each of these genetic pools were pinpointed and used to obtain genome ancestry mosaics of accessions. Diploid and triploid cultivars had genome mosaics involving three up to possibly seven contributors. The simplest mosaics were found for some diploid cultivars from New Guinea, combining three contributors, i.e., banksii and zebrina representing Musa acuminata subspecies and, more unexpectedly, the New Guinean species Musa schizocarpa. Breakpoints of M. schizocarpa introgressions were found to be conserved between New Guinea cultivars and the other analyzed diploid and triploid cultivars. This suggests that plants bearing these M. schizocarpa introgressions were transported from New Guinea and gave rise to currently cultivated bananas. Many cultivars showed contrasted mosaics with predominant ancestry from their geographical origin across Southeast Asia to New Guinea. This revealed that further diversification occurred in different Southeast Asian regions through hybridization with other Musa (sub)species, including two unknown ancestors that we propose to be M. acuminata ssp. halabanensis and a yet to be characterized M. acuminata subspecies. These results highlighted a dynamic crop formation process that was initiated in New Guinea, with subsequent diversification throughout Southeast Asia.

Hybridization, missing wild ancestors and the domestication of cultivated diploid bananas

Hybridization and introgressions are important evolutionary forces in plants. They contribute to the domestication of many species, including understudied clonal crops. Here, we examine their role in the domestication of a clonal crop of outmost importance, banana (Musa ssp.). We used genome-wide SNPs generated for 154 diploid banana cultivars and 68 samples of the wild M. acuminata to estimate and geo-localize the contribution of the different subspecies of M. acuminata to cultivated banana. We further investigated the wild to domesticate transition in New Guinea, an important domestication center. We found high levels of admixture in many cultivars and confirmed the existence of unknown wild ancestors with unequal contributions to cultivated diploid. In New Guinea, cultivated accessions exhibited higher diversity than their direct wild ancestor, the latter recovering from a bottleneck. Introgressions, balancing selection and positive selection were identified as important mechanisms for banana domestication. Our results shed new lights on the radiation of M. acuminata subspecies and on how they shaped banana domestication. They point candidate regions of origin for two unknown ancestors and suggest another contributor in New Guinea. This work feed research on the evolution of clonal crops and has direct implications for conservation, collection, and breeding.

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.

GeMo: a web-based platform for the visualization and curation of genome ancestry mosaics

In silico chromosome painting is a technique by which contributions of distinct genetic groups are represented along chromosomes of hybrid individuals. This type of analysis is used to study the mechanisms by which these individuals were formed. Such techniques are well adapted to identify genetic groups contributing to these individuals as well as hybridization events. It can also be used to follow chromosomal recombinations that occurred naturally or were generated by selective breeding. Here, we present GeMo, a novel interactive web-based and user-oriented interface to visualize in a linear-based fashion results of in silico chromosome painting. To facilitate data input generation, a script to execute analytical commands is provided and an interactive data curation mode is supported to ensure consistency of the automated procedure. GeMo contains preloaded datasets from published studies on crop domestication but can be applied to other purposes, such as breeding programs Although only applied so far on plants, GeMo can handle data from animals as well.

Database URL: https://gemo.southgreen.fr/

Application of Restriction Site-Associated DNA Sequencing (RAD-Seq) for Copy Number Variation and Triploidy Detection in Human

At present, low-pass whole-genome sequencing (WGS) is frequently used in clinical research and in the screening of copy number variations (CNVs). However, there are still some challenges in the detection of triploids. Restriction site-associated DNA sequencing (RAD-Seq) technology is a reduced-representation genome sequencing technology developed based on next-generation sequencing. Here, we verified whether RAD-Seq could be employed to detect CNVs and triploids. In this study, genomic DNA of 11 samples was extracted employing a routine method and used to build libraries. Five cell lines of known karyotypes and 6 triploid abortion tissue samples were included for RAD-Seq testing. The triploid samples were confirmed by STR analysis and also tested by low-pass WGS. The accuracy and efficiency of detecting CNVs and triploids by RAD-Seq were then assessed, compared with low-pass WGS. In our results, RAD-Seq detected 11 out of 11 (100%) chromosomal abnormalities, including 4 deletions and 1 aneuploidy in the purchased cell lines and all triploid samples. By contrast, these triploids were missed by low-pass WGS. Furthermore, RAD-Seq showed a higher resolution and more accurate allele frequency in the detection of triploids than low-pass WGS. Our study shows that, compared with low-pass WGS, RAD-Seq has relatively higher accuracy in CNV detection at a similar cost and is capable of identifying triploids. Therefore, the application of this technique in medical genetics has a significant potential value.

Reassessing Banana Phylogeny and Organelle Inheritance Modes Using Genome Skimming Data

Bananas (Musa spp.) are some of the most important fruit crops in the world, contributing up to US$10 billion in export values annually. In this study, we use high-throughput sequencing to obtain genomic resources of high-copy DNA molecules in bananas. We sampled 13 wild species and eight cultivars that represent the three genera (Ensete, Musa, and Musella) of the banana family (Musaceae). Their plastomic, 45S rDNA, and mitochondrial scaffolds were recovered from genome skimming data. Two major clades (Clades I & II) within Musa are strongly supported by the three genomic compartment data. We document, for the first time, that the plastomes of Musaceae have expanded inverted repeats (IR) after they diverged from their two close relatives, Heliconiaceae (the lobster-claws) and Strelitziaceae (the traveler's bananas). The presence/absence of rps19 within IR regions reinforces the two intra-generic clades within Musa. Our comparisons of the bananas' plastomic and mitochondrial DNA sequence trees aid in identifying hybrid bananas' parentage. As the mitochondrial genes of Musa have elevated substitution rates, paternal inheritance likely plays an influential role on the Musa mitogenome evolution. We propose genome skimming as a useful method for reliable genealogy tracing and phylogenetics in bananas.

Understanding polyploid banana origins. A commentary on: 'Unravelling the complex story of intergenomic recombination in ABB allotriploid bananas'

This article comments on: Alberto Cenci, Julie Sardos, Yann Hueber, Guillaume Martin, Catherine Breton, Nicolas Roux, Rony Swennen, Sebastien Christian Carpentier and Mathieu Rouard, Unravelling the complex story of intergenomic recombination in ABB allotriploid bananas, Annals of Botany, Volume 127, Issue 1, 01 January 2021, Pages 7–20, https://doi.org/10.1093/aob/mcaa032

Chromosome reciprocal translocations have accompanied subspecies evolution in bananas

Chromosome rearrangements and the way that they impact genetic differentiation and speciation have long raised questions from evolutionary biologists. They are also a major concern for breeders because of their bearing on chromosome recombination. Banana is a major crop that derives from inter(sub)specific hybridizations between various once geographically isolated Musa species and subspecies. We sequenced 155 accessions, including banana cultivars and representatives of Musa diversity, and genotyped-by-sequencing 1059 individuals from 11 progenies. We precisely characterized six large reciprocal translocations and showed that they emerged in different (sub)species of Musa acuminata, the main contributor to currently cultivated bananas. Most diploid and triploid cultivars analyzed were structurally heterozygous for 1 to 4 M. acuminata translocations, highlighting their complex origin. We showed that all translocations induced a recombination reduction of variable intensity and extent depending on the translocations, involving only the breakpoint regions, a chromosome arm, or an entire chromosome. The translocated chromosomes were found preferentially transmitted in many cases. We explore and discuss the possible mechanisms involved in this preferential transmission and its impact on translocation colonization.