A Genome-Wide Association Study on the Seedless Phenotype in Banana (Musa spp.) Reveals the Potential of a Selected Panel to Detect Candidate Genes in a Vegetatively Propagated Crop

Genome-wide association analyses using multilocus models on bananas (Musa spp.) reveal candidate genes related to morphology, fruit quality, and yield

Bananas (Musa spp.) are an essential fruit worldwide and rank as the fourth most significant food crop for addressing malnutrition due to their rich nutrients and starch content. The potential of their genetic diversity remains untapped due to limited molecular breeding tools. Our study examined a phenotypically diverse group of 124 accessions from the Colombian Musaceae Collection conserved in AGROSAVIA. We assessed 12 traits categorized into morphology, fruit quality, and yield, alongside sequence data. Our sequencing efforts provided valuable insights, with an average depth of about 7× per accession, resulting in 187,133 single-nucleotide polymorphisms (SNPs) against Musa acuminata (A genome) and 220,451 against Musa balbisiana (B genome). Population structure analysis grouped samples into four and five clusters based on the reference genome. By using different association models, we identified marker-trait associations (MTAs). The mixed linear model revealed four MTAs, while the Bayesian-information and linkage-disequilibrium iteratively nested keyway and fixed and random model for circulating probability unification models identified 82 and 70 MTAs, respectively. We identified 38 and 40 candidate genes in linkage proximity to significant MTAs for the A genome and B genome, respectively. Our findings provide insights into the genetic underpinnings of morphology, fruit quality, and yield. Once validated, the SNP markers and candidate genes can potentially drive advancements in genomic-guided breeding strategies to enhance banana crop improvement.

Integrative genomics reveals the polygenic basis of seedlessness in grapevine

Seedlessness is a crucial quality trait in table grape (Vitis vinifera L.) breeding. However, the development of seeds involved intricate regulations, and the polygenic basis of seed abortion remains unclear. Here, we combine comparative genomics, population genetics, quantitative genetics, and integrative genomics to unravel the evolution and polygenic basis of seedlessness in grapes. We generated the haplotype-resolved genomes for two seedless grape cultivars, "Thompson Seedless" (TS, syn. "Sultania") and "Black Monukka" (BM). Comparative genomics identified a ∼4.25 Mb hemizygous inversion on Chr10 specific in seedless cultivars, with seedless-associated genes VvTT16 and VvSUS2 located at breakpoints. Population genomic analyses of 548 grapevine accessions revealed two distinct clusters of seedless cultivars, and the identity-by-descent (IBD) results indicated that the origin of the seedlessness trait could be traced back to "Sultania." Introgression, rather than convergent selection, shaped the evolutionary history of seedlessness in grape improvement. Genome-wide association study (GWAS) analysis identified 110 quantitative trait loci (QTLs) associated with 634 candidate genes, including previously unidentified candidate genes, such as three 11S GLOBULIN SEED STORAGE PROTEIN and two CYTOCHROME P450 genes, and well-known genes like VviAGL11. Integrative genomic analyses resulted in 339 core candidate genes categorized into 13 functional categories related to seed development. Machine learning-based genomic selection achieved a remarkable prediction accuracy of 97% for seedlessness in grapevines. Our findings highlight the polygenic nature of seedlessness and provide candidate genes for molecular genetics and an effective prediction for seedlessness in grape genomic breeding.

Dissection of the Genetic Basis of Resistance to Stem Rot in Cultivated Peanuts (Arachis hypogaea L.) through Genome-Wide Association Study

Peanut (Arachis hypogaea) is an important oilseed and cash crop worldwide, contributing an important source of edible oil and protein for human nutrition. However, the incidence of stem rot disease caused by Athelia rolfsii poses a major challenge to peanut cultivation, resulting in significant yield losses. In this study, a panel of 202 peanut accessions was evaluated for their resistance to stem rot by inoculating plants in the field with A. rolfsii-infested oat grains in three environments. The mean disease index value of each environment for accessions in subsp. fasitigiate and subsp. hypogaea showed no significant difference. Accessions from southern China displayed the lowest disease index value compared to those from other ecological regions. We used whole-genome resequencing to analyze the genotypes of the accessions and to identify significant SNPs associated with stem rot resistance through genome-wide association study (GWAS). A total of 121 significant SNPs associated with stem rot resistance in peanut were identified, with phenotypic variation explained (PVE) ranging from 12.23% to 15.51%. A total of 27 candidate genes within 100 kb upstream and downstream of 23 significant SNPs were annotated, which have functions related to recognition, signal transduction, and defense response. These significant SNPs and candidate genes provide valuable information for further validation and molecular breeding to improve stem rot resistance in peanut.

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

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

A novel SNP panel developed for targeted genotyping-by-sequencing (GBS) reveals genetic diversity and population structure of Musa spp. germplasm collection

The Philippines is situated in the geographic region regarded as the center of diversity of banana and its wild relatives (Musa spp.). It holds the most extensive collection of B-genome germplasm in the world along with A-genome groups and several natural hybrids with A- and B-genome combinations. Management of this germplasm resource has relied immensely on identification using local names and morphological characters, and the extent of genetic diversity of the collection has not been achieved with molecular markers. A high-throughput and reliable genotyping method for banana and its relatives will facilitate germplasm management and support breeding initiatives toward a marker-based approach. Here, we developed a 1 K SNP genotyping panel based on filtering of high-quality genome-wide SNPs from the Musa Germplasm Information System and used it to assess the genetic diversity and population structure of 183 accessions from a Musa spp. germplasm collection containing Philippine and foreign accessions. Targeted GBS using SeqSNP™ technology generated 70,376,284 next-generation sequencing (NGS) reads with an average effective target SNP coverage of 340 × . Bioinformatics pipeline revealed 971 polymorphic SNPs containing 76.9% homozygous calls, 23.1% heterozygous calls and 4% with missing data. A final set of 952 SNPs detected 2,092 alleles. Pairwise genetic distance varied from 0.0021 to 0.3325 with most pairs of accessions distinguished with 250 to 300 loci. The SNP panel was able to detect seven (k = 7) genetically differentiated groups and its composition through Principal Component Analysis (PCA) with k-means clustering algorithm and Discriminant Analysis of Principal Components (DAPC). Accession-specific SNPs were also identified. The 1 K SNP panel effectively distinguishes between genomic groups and provides relatively good resolution of genome-wide nucleotide diversity of Musa spp. This panel is recommended for low-density genotyping for application in marker-assisted breeding and germplasm management, and could be further enhanced to increase marker density for other applications like genetic association and genomic selection in bananas.

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.

Diversity and determinants of bitterness, astringency, and fat content in cultivated Nacional and native Amazonian cocoa accessions from Ecuador

Cocoa (Theobroma cacao L.) is the only tree that can produce cocoa. Cocoa beans are highly sought after by chocolate makers to produce chocolate. Cocoa can be fine aromatic, characterized by floral and fruity notes, or it can be described as standard cocoa with a more pronounced cocoa aroma and bitterness. In this study, the genetic and biochemical determinants of sensorial notes and nonvolatile compounds related to bitterness, astringency, fat content, and protein content will be investigated in two populations: a cultivated modern Nacional population and a population of cocoa accessions collected recently in the Ecuadorian South Amazonia area of origin of the Nacional ancestral variety. For this purpose, a genome-wide association study (GWAS) was carried out on both populations, with results of biochemical compounds evaluated by near-infrared spectroscopy (NIRS) assays and with sensory evaluations. Twenty areas of associations were detected for sensorial data especially bitterness and astringency. Fifty-three areas of associations were detected linked to nonvolatile compounds. A total of 81 candidate genes could be identified in the areas of the association.

Perspective for genomic-enabled prediction against black sigatoka disease and drought stress in polyploid species

Genomic selection (GS) in plant breeding is explored as a promising tool to solve the problems related to the biotic and abiotic threats. Polyploid plants like bananas (Musa spp.) face the problem of drought and black sigatoka disease (BSD) that restrict their production. The conventional plant breeding is experiencing difficulties, particularly phenotyping costs and long generation interval. To overcome these difficulties, GS in plant breeding is explored as an alternative with a great potential for reducing costs and time in selection process. So far, GS does not have the same success in polyploid plants as with diploid plants because of the complexity of their genome. In this review, we present the main constraints to the application of GS in polyploid plants and the prospects for overcoming these constraints. Particular emphasis is placed on breeding for BSD and drought-two major threats to banana production-used in this review as a model of polyploid plant. It emerges that the difficulty in obtaining markers of good quality in polyploids is the first challenge of GS on polyploid plants, because the main tools used were developed for diploid species. In addition to that, there is a big challenge of mastering genetic interactions such as dominance and epistasis effects as well as the genotype by environment interaction, which are very common in polyploid plants. To get around these challenges, we have presented bioinformatics tools, as well as artificial intelligence approaches, including machine learning. Furthermore, a scheme for applying GS to banana for BSD and drought has been proposed. This review is of paramount impact for breeding programs that seek to reduce the selection cycle of polyploids despite the complexity of their genome.

A Flashforward Look into Solutions for Fruit and Vegetable Production

One of the most important challenges facing current and future generations is how climate change and continuous population growth adversely affect food security. To address this, the food system needs a complete transformation where more is produced in non-optimal and space-limited areas while reducing negative environmental impacts. Fruits and vegetables, essential for human health, are high-value-added crops, which are grown in both greenhouses and open field environments. Here, we review potential practices to reduce the impact of climate variation and ecosystem damages on fruit and vegetable crop yield, as well as highlight current bottlenecks for indoor and outdoor agrosystems. To obtain sustainability, high-tech greenhouses are increasingly important and biotechnological means are becoming instrumental in designing the crops of tomorrow. We discuss key traits that need to be studied to improve agrosystem sustainability and fruit yield.

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.

Cultivar-specific markers, mutations, and chimerisim of Cavendish banana somaclonal variants resistant to Fusarium oxysporum f. sp. cubense tropical race 4

Background

The selection of tissue culture–derived somaclonal variants of Giant Cavendish banana (Musa spp., Cavendish sub-group AAA) by the Taiwan Banana Research Institute (TBRI) has resulted in several cultivars resistant to Fusarium oxysporum f. sp. cubense tropical race 4 (Foc TR4), a destructive fungus threatening global banana production. However, the mutations in these somaclonal variants have not yet been determined. We performed an RNA-sequencing (RNA-seq) analysis of three TBRI Foc TR4–resistant cultivars: ‘Tai-Chiao No. 5’ (TC5), ‘Tai-Chiao No. 7’ (TC7), and ‘Formosana’ (FM), as well as their susceptible progenitor ‘Pei-Chiao’ (PC), to investigate the sequence variations among them and develop cultivar-specific markers.

Results

A group of single-nucleotide variants (SNVs) specific to one cultivar were identified from the analysis of RNA-seq data and validated using Sanger sequencing from genomic DNA. Several SNVs were further converted into cleaved amplified polymorphic sequence (CAPS) markers or derived CAPS markers that could identify the three Foc TR4–resistant cultivars among 6 local and 5 international Cavendish cultivars. Compared with PC, the three resistant cultivars showed a loss or alteration of heterozygosity in some chromosomal regions, which appears to be a consequence of single-copy chromosomal deletions. Notably, TC7 and FM shared a common deletion region on chromosome 5; however, different TC7 tissues displayed varying degrees of allele ratios in this region, suggesting the presence of chimerism in TC7.

Conclusions

This work demonstrates that reliable SNV markers of tissue culture–derived and propagated banana cultivars with a triploid genome can be developed through RNA-seq data analysis. Moreover, the analysis of sequence heterozygosity can uncover chromosomal deletions and chimerism in banana somaclonal variants. The markers obtained from this study will assist with the identification of TBRI Cavendish somaclonal variants for the quality control of tissue culture propagation, and the protection of breeders’ rights.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08692-5.

Morphoanatomy and Histochemistry of Septal Nectaries Related to Female Fertility in Banana Plants of the ‘Cavendish’ Subgroup

The objective of this study was to gain a deeper understanding of the morphoanatomical and histochemical structures that compose the nectary of pistillate flowers (female), which are involved in the female fertility of banana plants belonging to the ‘Cavendish’ subgroup. The diploid Calcutta 4 and the Grand Naine cultivar were used for the assessment. Five stages of floral development were proposed. Pistillate flower nectaries were subjected to morphological characterization, morphoanatomy, and histochemical tests (phenolic compounds, proteins, and lipids). Morphoanatomical analysis revealed a greater presence of narrow nectariferous ducts and more developed pluristratified papillae in Calcutta 4. In contrast, Grand Naine displayed cell disintegration in nectariferous ducts and pluristratified papillae, absent transmitting tissue, and greater amounts of vascular bundles at anthesis. However, Calcutta 4 displayed no changes in the nectariferous duct at any of the stages. An association was found between phenolic compounds and lipids in vacuoles adjacent to the vascular bundles, with greater amounts found in Grand Naine. The localization of phenolic compounds may suggest that these compounds play a role in nectar secretion or the oxidation of the nectary region, ultimately limiting the growth and passage of the pollen tube and preventing ovule fertilization.

Managing High-Density Genotyping Data with Gigwa

Next generation sequencing technologies enabled high-density genotyping for large numbers of samples. Nowadays SNP calling pipelines produce up to millions of such markers, but which need to be filtered in various ways according to the type of analyses. One of the main challenges still lies in the management of an increasing volume of genotyping files that are difficult to handle for many applications. Here, we provide a practical guide for efficiently managing large genomic variation data using Gigwa, a user-friendly, scalable and versatile application that may be deployed either remotely on web servers or on a local machine.

Continuous Mapping Identifies Loci Associated With Weevil Resistance [Cosmopolites sordidus (Germar)] in a Triploid Banana Population

The first step toward marker-assisted selection is linking the phenotypes to molecular markers through quantitative trait loci (QTL) analysis. While the process is straightforward in self-pollinating diploid (2x) species, QTL analysis in polyploids requires unconventional methods. In this study, we have identified markers associated with weevil Cosmopolites sordidus (Germar) resistance in bananas using 138 triploid (2n = 3x) hybrids derived from a cross between a tetraploid "Monyet" (2n = 4x) and a 2x "Kokopo" (2n = 2x) banana genotypes. The population was genotyped by Diversity Arrays Technology Sequencing (DArTSeq), resulting in 18,009 polymorphic single nucleotide polymorphisms (SNPs) between the two parents. Marker-trait association was carried out by continuous mapping where the adjusted trait means for the corm peripheral damage (PD) and total cross-section damage (TXD), both on the logit scale, were regressed on the marker allele frequencies. Forty-four SNPs that were associated with corm PD were identified on the chromosomes 5, 6, and 8, with 41 of them located on chromosome 6 and segregated in "Kokopo." Eleven SNPs associated with corm total TXD were identified on chromosome 6 and segregated in "Monyet." The additive effect of replacing one reference allele with the alternative allele was determined at each marker position. The PD QTL was confirmed using conventional QTL linkage analysis in the simplex markers segregating in "Kokopo" (AAAA × RA). We also identified 43 putative genes in the vicinity of the markers significantly associated with the two traits. The identified loci associated with resistance to weevil damage will be used in the efforts of developing molecular tools for marker-assisted breeding in bananas.

Genome-Wide Association Studies for Sex Determination and Cross-Compatibility in Water Yam (Dioscorea alata L.)

Yam (Dioscorea spp.) species are predominantly dioecious, with male and female flowers borne on separate individuals. Cross-pollination is, therefore, essential for gene flow among and within yam species to achieve breeding objectives. Understanding genetic mechanisms underlying sex determination and cross-compatibility is crucial for planning a successful hybridization program. This study used the genome-wide association study (GWAS) approach for identifying genomic regions linked to sex and cross-compatibility in water yam (Dioscorea alata L.). We identified 54 markers linked to flower sex determination, among which 53 markers were on chromosome 6 and one on chromosome 11. Our result ascertained that D. alata is characterized by the male heterogametic sex determination system (XX/XY). The cross-compatibility indices, average crossability rate (ACR) and percentage high crossability (PHC), were controlled by loci on chromosomes 1, 6 and 17. Of the significant loci, SNPs located on chromosomes 1 and 17 were the most promising for ACR and PHC, respectively, and should be validated for use in D. alata hybridization activities to predict cross-compatibility success. A total of 61 putative gene/protein families with direct or indirect influence on plant reproduction were annotated in chromosomic regions controlling the target traits. This study provides valuable insights into the genetic control of D. alata sexual reproduction. It opens an avenue for developing genomic tools for predicting hybridization success in water yam breeding programs.

Comprehensive genotyping of a Brazilian cassava (Manihot esculenta Crantz) germplasm bank: insights into diversification and domestication

KEY MESSAGE

Brazilian cassava diversity was characterized through population genetics and clustering approaches, highlighting contrasted genetic groups and spatial genetic differentiation. Cassava (Manihot esculenta Crantz) is a major staple root crop of the tropics, originating from the Amazonian region. In this study, 3354 cassava landraces and modern breeding lines from the Embrapa Cassava Germplasm Bank (CGB) were characterized. All individuals were subjected to genotyping-by-sequencing (GBS), identifying 27,045 single-nucleotide polymorphisms (SNPs). Identity-by-state and population structure analyses revealed a unique set of 1536 individuals and 10 distinct genetic groups with heterogeneous linkage disequilibrium (LD). On this basis, a density of 1300-4700 SNP markers were selected for large-effect quantitative trait loci (QTL) detection. Identified genetic groups were further characterized for population genetics parameters including minor allele frequency (MAF), observed heterozygosity [Formula: see text], effective population size estimate [Formula: see text]) and polymorphism information content (PIC). Selection footprints and introgressions of M. glaziovii were detected. Spatial population structure analysis revealed five ancestral populations related to distinct Brazilian ecoregions. Estimation of historical relationships among identified populations suggests an early population split from Amazonian to Atlantic forest and Caatinga ecoregions and active gene flows. This study provides a thorough genetic characterization of ex situ germplasm resources from cassava's center of origin, South America, with results shedding light on Brazilian cassava characteristics and its biogeographical landscape. These findings support and facilitate the use of genetic resources in modern breeding programs including implementation of association mapping and genomic selection strategies.

Comprehensive genotyping of a Brazilian cassava (Manihot esculenta Crantz) germplasm bank: insights into diversification and domestication

Brazilian cassava diversity was characterized through population genetics and clustering approaches, highlighting contrasted genetic groups and spatial genetic differentiation. Cassava (Manihot esculenta Crantz) is a major staple root crop of the tropics, originating from the Amazonian region. In this study, 3354 cassava landraces and modern breeding lines from the Embrapa Cassava Germplasm Bank (CGB) were characterized. All individuals were subjected to genotyping-by-sequencing (GBS), identifying 27,045 single-nucleotide polymorphisms (SNPs). Identity-by-state and population structure analyses revealed a unique set of 1536 individuals and 10 distinct genetic groups with heterogeneous linkage disequilibrium (LD). On this basis, a density of 1300-4700 SNP markers were selected for large-effect quantitative trait loci (QTL) detection. Identified genetic groups were further characterized for population genetics parameters including minor allele frequency (MAF), observed heterozygosity [Formula: see text], effective population size estimate [Formula: see text]) and polymorphism information content (PIC). Selection footprints and introgressions of M. glaziovii were detected. Spatial population structure analysis revealed five ancestral populations related to distinct Brazilian ecoregions. Estimation of historical relationships among identified populations suggests an early population split from Amazonian to Atlantic forest and Caatinga ecoregions and active gene flows. This study provides a thorough genetic characterization of ex situ germplasm resources from cassava's center of origin, South America, with results shedding light on Brazilian cassava characteristics and its biogeographical landscape. These findings support and facilitate the use of genetic resources in modern breeding programs including implementation of association mapping and genomic selection strategies.

Recent Large-Scale Genotyping and Phenotyping of Plant Genetic Resources of Vegetatively Propagated Crops

Several recent national and international projects have focused on large-scale genotyping of plant genetic resources in vegetatively propagated crops like fruit and berries, potatoes and woody ornamentals. The primary goal is usually to identify true-to-type plant material, detect possible synonyms, and investigate genetic diversity and relatedness among accessions. A secondary goal may be to create sustainable databases that can be utilized in research and breeding for several years ahead. Commonly applied DNA markers (like microsatellite DNA and SNPs) and next-generation sequencing each have their pros and cons for these purposes. Methods for large-scale phenotyping have lagged behind, which is unfortunate since many commercially important traits (yield, growth habit, storability, and disease resistance) are difficult to score. Nevertheless, the analysis of gene action and development of robust DNA markers depends on environmentally controlled screening of very large sets of plant material. Although more time-consuming, co-operative projects with broad-scale data collection are likely to produce more reliable results. In this review, we will describe some of the approaches taken in genotyping and/or phenotyping projects concerning a wide variety of vegetatively propagated crops.

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.

Genome-wide association studies of fruit quality traits in jujube germplasm collections using genotyping-by-sequencing

Chinese jujube (Ziziphus jujuba Mill.) is an important fruit crop and harbors many highly diverse traits of potential economic importance. Fruit size, stone size, and fruit cracking have an important influence on the commercial value of jujube. This study is the first to conduct a genome-wide association study (GWAS) on 180 accessions of jujube and focuses on locating single-nucleotide polymorphisms (SNPs) associated with nine important fruit quality traits. Genotyping was performed using genotyping-by-sequencing and 4651 high-quality SNPs were identified. A genetic diversity analysis revealed the presence of three distinct groups, and rapid linkage disequilibrium decay was observed in this jujube population. Using a mixed linear model, a total of 45 significant SNP-trait associations were detected, among which 33 SNPs had associations with fruit size-related traits, nine were associated with stone size-related traits, and three with fruit cracking-related traits. In total, 21 candidate genes involved in cell expansion, abiotic stress responses, hormone signaling, and growth development were identified from the genome sequences of jujube. These results are useful as basic data for GWAS of other jujube traits, and these significant SNP loci and candidate genes should aid marker-assisted breeding and genomic selection of improved jujube cultivars.

Gigwa v2-Extended and improved genotype investigator

Background

The study of genetic variations is the basis of many research domains in biology. From genome structure to population dynamics, many applications involve the use of genetic variants. The advent of next-generation sequencing technologies led to such a flood of data that the daily work of scientists is often more focused on data management than data analysis. This mass of genotyping data poses several computational challenges in terms of storage, search, sharing, analysis, and visualization. While existing tools try to solve these challenges, few of them offer a comprehensive and scalable solution.

Results

Gigwa v2 is an easy-to-use, species-agnostic web application for managing and exploring high-density genotyping data. It can handle multiple databases and may be installed on a local computer or deployed as an online data portal. It supports various standard import and export formats, provides advanced filtering options, and offers means to visualize density charts or push selected data into various stand-alone or online tools. It implements 2 standard RESTful application programming interfaces, GA4GH, which is health-oriented, and BrAPI, which is breeding-oriented, thus offering wide possibilities of interaction with third-party applications. The project home page provides a list of live instances allowing users to test the system on public data (or reasonably sized user-provided data).

Conclusions

This new version of Gigwa provides a more intuitive and more powerful way to explore large amounts of genotyping data by offering a scalable solution to search for genotype patterns, functional annotations, or more complex filtering. Furthermore, its user-friendliness and interoperability make it widely accessible to the life science community.

Association genetics of bunch weight and its component traits in East African highland banana (Musa spp. AAA group)

The major quantitative trait loci associated with bunch weight and its component traits in the East African highland banana-breeding population are located on chromosome 3. Bunch weight increase is one of the major objectives of banana improvement programs, but little is known about the loci controlling bunch weight and its component traits. Here we report for the first time some genomic loci associated with bunch weight and its component traits in banana as revealed through a genome-wide association study. A banana-breeding population of 307 genotypes varying in ploidy was phenotyped in three locations under different environmental conditions, and data were collected on bunch weight, number of hands and fruits; fruit length and circumference; and diameter of both fruit and pulp for three crop cycles. The population was genotyped with genotyping by sequencing and 27,178 single nucleotide polymorphisms (SNPs) were generated. The association between SNPs and the best linear unbiased predictors of traits was performed with TASSEL v5 using a mixed linear model accounting for population structure and kinship. Using Bonferroni correction, false discovery rate, and long-range linkage disequilibrium (LD), 25 genomic loci were identified with significant SNPs and most were localized on chromosome 3. Most SNPs were located in genes encoding uncharacterized and hypothetical proteins, but some mapped to transcription factors and genes involved in cell cycle regulation. Inter-chromosomal LD of SNPs was present in the population, but none of the SNPs were significantly associated with the traits. The clustering of significant SNPs on chromosome 3 supported our hypothesis that fruit filling in this population was under control of a few quantitative trait loci with major effects.

The Future of Digital Sequence Information for Plant Genetic Resources for Food and Agriculture

The recent debates on the legal status of "digital sequence information" (DSI) at the international level could have extensive consequences for the future of agriculture and food security. A large majority of recent advances in biology, medicine, or agriculture were achieved by sharing and mining of freely accessible sequencing data. It is most probably because of the tremendous success of modern genomics and advances of synthetic biology that concerns were raised about possible fair and equitable ways of sharing data. The DSI concept is relatively new, and all concerned parties agreed upon the need for a clear definition. For example, the extent to which DSI understanding is limited only to genetic sequence data has to be clarified. In this paper, I focus on a subset of DSI essential to humankind: the DSI originating from plant genetic resources for food and agriculture (PGRFA). Two international agreements shape the conservation and use of plant genetic resources: the Convention on Biodiversity and the International Treaty for Plant Genetic Resources for Food and Agriculture. In an attempt to mobilize DSI users and producers involved in research, breeding, and conservation, I describe here how the increasing amount of genomic data, information, and studies interact with the existing legal framework at the global level. Using possible scenarios, I will emphasize the complexity of the issues surrounding DSI for PGRFA and propose potential ways forward for developing an inclusive governance and fair use of these genetic resources.

Identification and characterization of cherry (Cerasus pseudocerasus G. Don) genes responding to parthenocarpy induced by GA3 through transcriptome analysis

Background

Fruit set after successful pollination is key for the production of sweet cherries, and a low fruit-setting rate is the main problem in production of this crop. As gibberellin treatment can directly induce parthenogenesis and satisfy the hormone requirement during fruit growth and development, such treatment is an important strategy for improving the fruit-setting rate of sweet cherries. Previous studies have mainly focused on physiological aspects, such as fruit quality, fruit size, and anatomical structure, whereas the molecular mechanism remains clear.

Results

In this study, we analyzed the transcriptome of ‘Meizao’ sweet cherry fruit treated with gibberellin during the anthesis and hard-core periods to identify genes associated with parthenocarpic fruit set. A total of 25,341 genes were identified at the anthesis and hard-core stages, 765 (681 upregulated, 84 downregulated) and 186 (141 upregulated, 45 downregulated) of which were significant differentially expressed genes (DEGs) at the anthesis and the hard-core stages after gibberellin 3 (GA3) treatment, respectively. Based on DEGs between the control and GA3 treatments, the GA3 response mainly involves parthenocarpic fruit set and cell division. Exogenous gibberellin stimulated sweet cherry fruit parthenocarpy and enlargement, as verified by qRT-PCR results of related genes as well as the parthenocarpic fruit set and fruit size. Based on our research and previous studies in Arabidopsis thaliana, we identified key genes associated with parthenocarpic fruit set and cell division. Interestingly, we observed patterns among sweet cherry fruit setting-related DEGs, especially those associated with hormone balance, cytoskeleton formation and cell wall modification.

Conclusions

Overall, the result provides a possible molecular mechanism regulating parthenocarpic fruit set that will be important for basic research and industrial development of sweet cherries.

Electronic supplementary material

The online version of this article (10.1186/s12863-019-0746-8) contains supplementary material, which is available to authorized users.

A genome-wide association study identified loci for yield component traits in sugarcane (Saccharum spp.)

Sugarcane (Saccharum spp.) has a complex genome with variable ploidy and frequent aneuploidy, which hampers the understanding of phenotype and genotype relations. Despite this complexity, genome-wide association studies (GWAS) may be used to identify favorable alleles for target traits in core collections and then assist breeders in better managing crosses and selecting superior genotypes in breeding populations. Therefore, in the present study, we used a diversity panel of sugarcane, called the Brazilian Panel of Sugarcane Genotypes (BPSG), with the following objectives: (i) estimate, through a mixed model, the adjusted means and genetic parameters of the five yield traits evaluated over two harvest years; (ii) detect population structure, linkage disequilibrium (LD) and genetic diversity using simple sequence repeat (SSR) markers; (iii) perform GWAS analysis to identify marker-trait associations (MTAs); and iv) annotate the sequences giving rise to SSR markers that had fragments associated with target traits to search for putative candidate genes. The phenotypic data analysis showed that the broad-sense heritability values were above 0.48 and 0.49 for the first and second harvests, respectively. The set of 100 SSR markers produced 1,483 fragments, of which 99.5% were polymorphic. These SSR fragments were useful to estimate the most likely number of subpopulations, found to be four, and the LD in BPSG, which was stronger in the first 15 cM and present to a large extension (65 cM). Genetic diversity analysis showed that, in general, the clustering of accessions within the subpopulations was in accordance with the pedigree information. GWAS performed through a multilocus mixed model revealed 23 MTAs, six, three, seven, four and three for soluble solid content, stalk height, stalk number, stalk weight and cane yield traits, respectively. These MTAs may be validated in other populations to support sugarcane breeding programs with introgression of favorable alleles and marker-assisted selection.

The Occurrence of Seedlessness in Higher Plants; Insights on Roles and Mechanisms of Parthenocarpy

Parthenocarpy in a broad sense includes those processes that allow the production of seedless fruits. Such fruits are favorable to growers, because they are set independently of successful pollination, and to processors and consumers, because they are easier to deal with and to eat. Seedless fruits however represent a biological paradox because they do not contribute to offspring production. In this work, the occurrence of parthenocarpy in Angiosperms was investigated by conducting a bibliographic survey. We distinguished monospermic (single seeded) from plurispermic (multiseeded) species and wild from cultivated taxa. Out of 96 seedless taxa, 66% belonged to plurispermic species. Of these, cultivated species were represented six times higher than wild species, suggesting a selective pressure for parthenocarpy during domestication and breeding. In monospermic taxa, wild and cultivated species were similarly represented. The occurrence of parthenocarpy in wild species suggests that seedlessness may have an adaptive role. In monospermic species, seedless fruits are proposed to reduce seed predation through deceptive mechanisms. In plurispermic fruit species, parthenocarpy may exert an adaptive advantage under suboptimal pollination regimes, when too few embryos are formed to support fruit growth. In this situation, parthenocarpy offers the opportunity to accomplish the production and dispersal of few seeds, thus representing a selective advantage. Approximately 20 sources of seedlessness have been described in tomato. Excluding the EMS induced mutation parthenocarpic fruit (pat), the parthenocarpic phenotype always emerged in biparental populations derived from wide crosses between cultivated tomato and wild relatives. Following a theory postulated for apomictic species, we argument that wide hybridization could also be the force driving parthenocarpy, following the disruption of synchrony in time and space of reproductive developmental events, from sporogenesis to fruit development. The high occurrence of polyploidy among parthenocarpic species supported this suggestion. Other commonalities between apomixis and parthenocarpy emerged from genetic and molecular studies of the two phenomena. Such insights may improve the understanding of the mechanisms underlying these two reproductive variants of great importance to modern breeding.

Using Genomic Sequence Information to Increase Conservation and Sustainable Use of Crop Diversity and Benefit-Sharing

This article describes how CGIAR centers and partners are using genomic sequence information to promote the conservation and sustainable use of crop genetic diversity, and to generate and share benefits derived from those uses. The article highlights combined institutional, and benefit-sharing-related challenges that need to be addressed to support expanded use of digital sequence information in agricultural research and development.

Three New Genome Assemblies Support a Rapid Radiation in Musa acuminata (Wild Banana)

Edible bananas result from interspecific hybridization between Musa acuminata and Musa balbisiana, as well as among subspecies in M. acuminata. Four particular M. acuminata subspecies have been proposed as the main contributors of edible bananas, all of which radiated in a short period of time in southeastern Asia. Clarifying the evolution of these lineages at a whole-genome scale is therefore an important step toward understanding the domestication and diversification of this crop. This study reports the de novo genome assembly and gene annotation of a representative genotype from three different subspecies of M. acuminata. These data are combined with the previously published genome of the fourth subspecies to investigate phylogenetic relationships. Analyses of shared and unique gene families reveal that the four subspecies are quite homogenous, with a core genome representing at least 50% of all genes and very few M. acuminata species-specific gene families. Multiple alignments indicate high sequence identity between homologous single copy-genes, supporting the close relationships of these lineages. Interestingly, phylogenomic analyses demonstrate high levels of gene tree discordance, due to both incomplete lineage sorting and introgression. This pattern suggests rapid radiation within Musa acuminata subspecies that occurred after the divergence with M. balbisiana. Introgression between M. a. ssp. malaccensis and M. a. ssp. burmannica was detected across the genome, though multiple approaches to resolve the subspecies tree converged on the same topology. To support evolutionary and functional analyses, we introduce the PanMusa database, which enables researchers to exploration of individual gene families and trees.

Plant genetic resources for food and agriculture: opportunities and challenges emerging from the science and information technology revolution

Contents Summary 1407 I. Introduction 1408 II. Technological advances and their utility for gene banks and breeding, and longer-term contributions to SDGs 1408 III. The challenges that must be overcome to realise emerging R&D opportunities 1410 IV. Renewed governance structures for PGR (and related big data) 1413 V. Access and benefit sharing and big data 1416 VI. Conclusion 1417 Acknowledgements 1417 ORCID 1417 References 1417 SUMMARY: Over the last decade, there has been an ongoing revolution in the exploration, manipulation and synthesis of biological systems, through the development of new technologies that generate, analyse and exploit big data. Users of Plant Genetic Resources (PGR) can potentially leverage these capacities to significantly increase the efficiency and effectiveness of their efforts to conserve, discover and utilise novel qualities in PGR, and help achieve the Sustainable Development Goals (SDGs). This review advances the discussion on these emerging opportunities and discusses how taking advantage of them will require data integration and synthesis across disciplinary, organisational and international boundaries, and the formation of multi-disciplinary, international partnerships. We explore some of the institutional and policy challenges that these efforts will face, particularly how these new technologies may influence the structure and role of research for sustainable development, ownership of resources, and access and benefit sharing. We discuss potential responses to political and institutional challenges, ranging from options for enhanced structure and governance of research discovery platforms to internationally brokered benefit-sharing agreements, and identify a set of broad principles that could guide the global community as it seeks or considers solutions.

The making of virgin fruit: the molecular and genetic basis of parthenocarpy

Fruit set-the commitment of an angiosperm plant to develop fruit-is a key developmental process that normally occurs following successful fertilization. Parthenocarpy arises when fruit automatically develop in the absence of fertilization. This review uses parthenocarpic fruit development as a focal device through which to recapitulate and understand the molecular effectors that mediate and regulate fruit set. The review demonstrates that studies of parthenocarpy are providing vital insight into plant development, signaling and, potentially, high-value agricultural products.

MGIS: managing banana (Musa spp.) genetic resources information and high-throughput genotyping data

Unraveling the genetic diversity held in genebanks on a large scale is underway, due to advances in Next-generation sequence (NGS) based technologies that produce high-density genetic markers for a large number of samples at low cost. Genebank users should be in a position to identify and select germplasm from the global genepool based on a combination of passport, genotypic and phenotypic data. To facilitate this, a new generation of information systems is being designed to efficiently handle data and link it with other external resources such as genome or breeding databases. The Musa Germplasm Information System (MGIS), the database for global ex situ-held banana genetic resources, has been developed to address those needs in a user-friendly way. In developing MGIS, we selected a generic database schema (Chado), the robust content management system Drupal for the user interface, and Tripal, a set of Drupal modules which links the Chado schema to Drupal. MGIS allows germplasm collection examination, accession browsing, advanced search functions, and germplasm orders. Additionally, we developed unique graphical interfaces to compare accessions and to explore them based on their taxonomic information. Accession-based data has been enriched with publications, genotyping studies and associated genotyping datasets reporting on germplasm use. Finally, an interoperability layer has been implemented to facilitate the link with complementary databases like the Banana Genome Hub and the MusaBase breeding database.

Database URL:https://www.crop-diversity.org/mgis/

Homeolog expression analysis in an allotriploid non-model crop via integration of transcriptomics and proteomics

The fate of doubled genes, from allopolyploid or autopolyploid origin, is controlled at multiple levels, resulting in the modern day cultivars. We studied the root growth of 3 different triploid banana cultivars under control and osmotic stress conditions. The root growth of the allopolyploid ABB cultivar was 42% higher under control and 61% higher under osmotic stress. By integrating transcriptomics and proteomics, we studied the gene expression of all 3 cultivars, resulting in 2,749 identified root proteins. 383 gene loci displayed genotype specific differential expression whereof 252 showed at least one Single Amino Acid Polymorphism (SAAP). In the ABB cultivar, allele expressions supposedly follow a 1/3 and 2/3 pattern for respectively the A and the B allele. Using transcriptome read alignment to assess the homeoallelic contribution we found that 63% of the allele specific genes deviated from this expectation. 32 gene loci even did not express the A allele. The identified ABB allele- specific proteins correlate well with the observed growth phenotype as they are enriched in energy related functions such as ATP metabolic processes, nicotinamide nucleotide metabolic processes, and glycolysis.

Translating the "Banana Genome" to Delineate Stress Resistance, Dwarfing, Parthenocarpy and Mechanisms of Fruit Ripening

Evolutionary frozen, genetically sterile and globally iconic fruit "Banana" remained untouched by the green revolution and, as of today, researchers face intrinsic impediments for its varietal improvement. Recently, this wonder crop entered the genomics era with decoding of structural genome of double haploid Pahang (AA genome constitution) genotype of Musa acuminata. Its complex genome decoded by hybrid sequencing strategies revealed panoply of genes and transcription factors involved in the process of sucrose conversion that imparts sweetness to its fruit. Historically, banana has faced the wrath of pandemic bacterial, fungal, and viral diseases and multitude of abiotic stresses that has ruined the livelihood of small/marginal farmers' and destroyed commercial plantations. Decoding structural genome of this climacteric fruit has given impetus to a deeper understanding of the repertoire of genes involved in disease resistance, understanding the mechanism of dwarfing to develop an ideal plant type, unraveling the process of parthenocarpy, and fruit ripening for better fruit quality. Further, injunction of comparative genomics will usher in integration of information from its decoded genome and other monocots into field applications in banana related but not limited to yield enhancement, food security, livelihood assurance, and energy sustainability. In this mini review, we discuss pre- and post-genomic discoveries and highlight accomplishments in structural genomics, genetic engineering and forward genetic accomplishments with an aim to target genes and transcription factors for translational research in banana.