Genetic diversity and population structure of Musa accessions in ex situ conservation

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.

Banana Classification Using Sanger Sequencing of the Ribosomal DNA Internal Transcribed Spacer (ITS) Region

Banana (Musa spp.) is one of the most economically important horticultural crops. There are many types of banana, with differing ploidy (usually diploid, triploid, or tetraploid) and genome types (most containing the A or/and B genome). Currently, observation and genome type detection are commonly used to identify banana germplasm resources. However, observation is tedious, while genome type detection cannot distinguish categories below genome types. It is, therefore, urgent to establish a simple and effective method for identifying banana germplasm resources. This study sequenced and analyzed the ribosomal DNA internal transcribed spacer (ITS) sequences of 62 banana germplasm resources and found that the sequencing peaks, especially the 20 bp region near the 420-bp position (referred to as the 420-bp region), exhibited relatively recognizable and repeatable polymorphism characteristics. Using the 420-bp region as a marker, we were able to quickly distinguish bananas belonging to different genome type groups or different subgroups in the same genome type group. Moreover, it appeared that Sanger sequencing of ITS could be used to identify hybrid banana offspring. In general, ITS sequencing simplifies the classification of banana germplasm resources and has potential application in several areas of Musa improvement.

Origin and evolution of the triploid cultivated banana genome

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

Pollen viabilities and gene expression profiles across Musa genomes

Banana (Musa spp.) is a major global economic fruit crop. However, cross-pollination from other Musa cultivars grown in nearby plantations results in seeded fruit that exceeds market demand. This study investigated pollen viability and germination and examined the expression profiles of pollen development-related genes across seven Musa genomes (AA, BB, AAA, BBB, AAB, ABB and ABBB). Twenty-three Musa cultivars were assessed for pollen viability using lacto-aceto-orcein and triphenyltetrazolium chloride staining methods. Results revealed that pollen viability obtained from both methods was significantly different among all the studied cultivars. Cultivars carrying BB (diploid) genomes had higher viability percentages than AA (diploid), AAA, BBB, AAB and ABB (triploid) and ABBB (tetraploid) genomes. Germination of the studied cultivars was also investigated on pollen culture medium, with results showing significant differences between the pollen of each cultivar. The best germinating cultivar was TKM (11.0 %), carrying BB genome. Expression profiles of pollen development-related genes by RT-qPCR indicated that both TPD1A and MYB80 genes were highly expressed in triploid Musa genomes but the PTC1 gene showed down-regulated expression, resulting in non-viable pollen. Pollen viability, pollen germination and pollen development-related genes differed across Musa cultivars. This knowledge will be useful for the selection of male parents for Musa cross-breeding programs. Pollen viability should also be considered when planning Musa production to avoid seeded fruit.

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.

Gene Expression, Histology and Histochemistry in the Interaction between Musa sp. and Pseudocercospora fijiensis

Bananas are the main fruits responsible for feeding more than 500 million people in tropical and subtropical countries. Black Sigatoka, caused by the fungus Pseudocercospora fijiensis, is one of the most destructive disease for the crop. This fungus is mainly controlled with the use of fungicides; however, in addition to being harmful to human health, they are associated with a high cost. The development of resistant cultivars through crosses of susceptible commercial cultivars is one of the main focuses of banana breeding programs worldwide. Thus, the objective of the present study was to investigate the interaction between Musa sp. and P. fijiensis through the relative expression of candidate genes involved in the defence response to black Sigatoka in four contrasting genotypes (resistant: Calcutta 4 and Krasan Saichon; susceptible: Grand Naine and Akondro Mainty) using quantitative real-time PCR (RT-qPCR) in addition to histological and histochemical analyses to verify the defence mechanisms activated during the interaction. Differentially expressed genes (DEGs) related to the jasmonic acid and ethylene signalling pathway, GDSL-like lipases and pathogenesis-related proteins (PR-4), were identified. The number and distance between stomata were directly related to the resistance/susceptibility of each genotype. Histochemical tests showed the production of phenolic compounds and callosis as defence mechanisms activated by the resistant genotypes during the interaction process. Scanning electron microscopy (SEM) showed pathogenic structures on the leaf surface in addition to calcium oxalate crystals. The resistant genotype Krasan Saichon stood out in the analyses and has potential for use in breeding programs for resistance to black Sigatoka in banana and plantains.

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.

Erianthus germplasm collection in Thailand: genetic structure and phylogenetic aspects of tetraploid and hexaploid accessions

BACKGROUND

The genus Erianthus, which belongs to the "Saccharum complex", includes C₄ warm-season grasses. Erianthus species are widely distributed throughout Southeast Asia, East Asia and South Asia. Erianthus arundinaceus (Retz.) Jeswiet is highly adaptable to the environment, has a high percentage of dry matter, and is highly productive. Recently, this species has attracted attention as a novel bioenergy crop and as a breeding material for sugarcane improvement. Such interest in E. arundinaceus has accelerated the collection and conservation of its genetic resources, mainly in Asian countries, and also evaluation of morphological, agricultural, and cytogenetic features in germplasm collections. In Thailand, genetic resources of E. arundinaceus have been collected over the past 20 years and their phenotypic traits have been evaluated. However, the genetic differences and relatedness of the germplasms are not fully understood.

RESULTS

A set of 41 primer pairs for nuclear simple sequence repeats (SSRs) developed from E. arundinaceus were used to assess the genetic diversity of 121 Erianthus germplasms collected in Thailand; of these primer pairs, 28 detected a total of 316 alleles. A Bayesian clustering approach with these alleles classified the accessions into four main groups, generally corresponding to the previous classification based on phenotypic analysis. The results of principal coordinate analysis and phylogenetic analysis of the 121 accessions on the basis of the SSR markers showed the same trend as Bayesian clustering, whereas sequence variations of three non-coding regions of chloroplast DNA revealed eight haplotypes among the accessions. The analysis of genetic structure and phylogenetic relationships, however, found some accessions whose classification contradicted the results of previous phenotypic classification.

CONCLUSIONS

The molecular approach used in this study characterized the genetic diversity and relatedness of Erianthus germplasms collected across Thailand. This knowledge would allow efficient maintenance and conservation of the genetic resources of this grass and would help to use Erianthus species as breeding materials for development of novel bioenergy crops and sugarcane improvement.

Maximizing genetic representation in seed collections from populations of self and cross-pollinated banana wild relatives

BACKGROUND

Conservation of plant genetic resources, including the wild relatives of crops, plays an important and well recognised role in addressing some of the key challenges faced by humanity and the planet including ending hunger and biodiversity loss. However, the genetic diversity and representativeness of ex situ collections, especially that contained in seed collections, is often unknown. This limits meaningful assessments against conservation targets, impairs targeting of future collecting and limits their use. We assessed genetic representation of seed collections compared to source populations for three wild relatives of bananas and plantains. Focal species and sampling regions were M. acuminata subsp. banksii (Papua New Guinea), M. balbisiana (Viet Nam) and M. maclayi s.l. (Bougainville, Papua New Guinea). We sequenced 445 samples using suites of 16-20 existing and newly developed taxon-specific polymorphic microsatellite markers. Samples of each species were from five populations in a region; 15 leaf samples from different individuals and 16 seed samples from one infructescence ('bunch') were analysed for each population.

RESULTS

Allelic richness of seeds compared to populations was 51, 81 and 93% (M. acuminata, M. balbisiana and M. maclayi respectively). Seed samples represented all common alleles in populations but omitted some rarer alleles. The number of collections required to achieve the 70% target of the Global Strategy for Plant Conservation was species dependent, relating to mating systems. Musa acuminata populations had low heterozygosity and diversity, indicating self-fertilization; many bunches were needed (> 15) to represent regional alleles to 70%; over 90% of the alleles from a bunch are included in only two seeds. Musa maclayi was characteristically cross-fertilizing; only three bunches were needed to represent regional alleles; within a bunch, 16 seeds represent alleles. Musa balbisiana, considered cross-fertilized, had low genetic diversity; seeds of four bunches are needed to represent regional alleles; only two seeds represent alleles in a bunch.

CONCLUSIONS

We demonstrate empirical measurement of representation of genetic material in seeds collections in ex situ conservation towards conservation targets. Species mating systems profoundly affected genetic representation in seed collections and therefore should be a primary consideration to maximize genetic representation. Results are applicable to sampling strategies for other wild species.

Comparative transcriptome analysis of unripe and ripe banana (cv. Nendran) unraveling genes involved in ripening and other related processes

Banana is one of the most important fruit crops consumed globally owing to its high nutritional value. Previously, we demonstrated that the ripe pulp of the banana cultivar (cv.) Nendran (AAB) contained a high amount of pro-vitamin A carotenoids. However, the molecular factors involved in the ripening process in Nendran fruit are unexplored. Hence, we commenced a transcriptome study by using the Illumina HiSeq 2500 at two stages i.e. unripe and ripe fruit-pulp of Nendran. Overall, 3474 up and 4727 down-regulated genes were obtained. A large number of identified transcripts were related to genes involved in ripening, cell wall degradation and aroma formation. Gene ontology analysis highlighted differentially expressed genes that play a key role in various pathways. These pathways were mainly linked to cellular, molecular and biological processes. The present transcriptome study also reveals a crucial role of up-regulated carotenoid biosynthesis pathway genes namely, lycopene beta cyclase and geranylgeranyl pyrophosphate synthase at the ripening stage. Genes related to the ripening and other processes like aroma and flavor were highly expressed in the ripe pulp. Expression of numerous transcription factor family genes was also identified. This study lays a path towards understanding the ripening, carotenoid accumulation and other related processes in banana.

Genome-Wide Novel Genic Microsatellite Marker Resource Development and Validation for Genetic Diversity and Population Structure Analysis of Banana

Trait tagging through molecular markers is an important molecular breeding tool for crop improvement. SSR markers encoded by functionally relevant parts of a genome are well suited for this task because they may be directly related to traits. However, a limited number of these markers are known for Musa spp. Here, we report 35136 novel functionally relevant SSR markers (FRSMs). Among these, 17,561, 15,373 and 16,286 FRSMs were mapped in-silico to the genomes of Musa acuminata, M. balbisiana and M. schizocarpa, respectively. A set of 273 markers was validated using eight accessions of Musa spp., from which 259 markers (95%) produced a PCR product of the expected size and 203 (74%) were polymorphic. In-silico comparative mapping of FRSMs onto Musa and related species indicated sequence-based orthology and synteny relationships among the chromosomes of Musa and other plant species. Fifteen FRSMs were used to estimate the phylogenetic relationships among 50 banana accessions, and the results revealed that all banana accessions group into two major clusters according to their genomic background. Here, we report the first large-scale development and characterization of functionally relevant Musa SSR markers. We demonstrate their utility for germplasm characterization, genetic diversity studies, and comparative mapping in Musa spp. and other monocot species. The sequences for these novel markers are freely available via a searchable web interface called Musa Marker Database.

Molecular identification and evolutionary relationships between the subspecies of Musa by DNA barcodes

BACKGROUND

The banana (Musa sp., AAA) genome is constantly increasing due to high-frequency of somaclonal variations. Due to its large diversity, a conventional numerical and morphological based taxonomic identification of banana cultivars is laborious, difficult and often leads to subject of disagreements.

RESULTS

Hence, in the present study, we used universal DNA barcode ITS2 region to identify and to find the genetic relationship between the cultivars and varieties of banana. Herein, a total of 16 banana cultivars were PCR amplified using ITS2 primer pair. In addition, 321 sequences which were retrieved from GenBank, USA, were used in this study. The sequences were then aligned using Clustal W and genetic distances were computed using MEGA V5.1. The study showed significant divergence between the intra- and inter-specific genetic distances in ITS2 region. BLAST1 and Distance methods proved that ITS2 DNA barcode region successfully identified and distinguished the cultivar and varieties of banana.

CONCLUSION

Thus, from the results of the present study, it is clear that ITS2 is not only an efficient DNA barcode to identify the banana species but also a potential candidate for enumerating the phylogenetic relationships between the subspecies and cultivars. This is the first comprehensive study to categorically distinguish the economically important banana subspecies and varieties using DNA barcodes and to understand its evolutionary relationship.

Assessment of metabolic variability and diversity present in leaf, peel and pulp tissue of diploid and triploid Musa spp

Banana (Musa spp.) plants produce many health promoting compounds in leaf, peel and pulp. For a robust metabolic analysis of these tissues, leaf at five developmental stages were compared to assess suitable sampling practices. Results confirmed that the common sampling practise of leaf 3 is applicable for metabolic comparisons. The developed work flow was applied to analyse the metabolite diversity present in 18 different Musa varieties, providing baseline levels of metabolites in leaf, peel and pulp tissue. Correlation analysis was then used to ascertain whether similar trends can be detected in the three plant tissues of the diversity panel. The genome group displayed a dominant role in the composition of the metabolome in all three tissues. This led to the conclusion that a correlation between tissues was only possible within a genome group as the different parental backgrounds caused too great a variation in the metabolomes. It also suggests the metabolome could be used to monitor the interaction/hybridisation of genomes during breeding programmes.

Musa balbisiana genome reveals subgenome evolution and functional divergence

Banana cultivars (Musa ssp.) are diploid, triploid and tetraploid hybrids derived from Musa acuminata and Musa balbisiana. We presented a high-quality draft genome assembly of M. balbisiana with 430 Mb (87%) assembled into 11 chromosomes. We identified that the recent divergence of M. acuminata (A-genome) and M. balbisiana (B-genome) occurred after lineage-specific whole-genome duplication, and that the B-genome may be more sensitive to the fractionation process compared to the A-genome. Homoeologous exchanges occurred frequently between A- and B-subgenomes in allopolyploids. Genomic variation within progenitors resulted in functional divergence of subgenomes. Global homoeologue expression dominance occurred between subgenomes of the allotriploid. Gene families related to ethylene biosynthesis and starch metabolism exhibited significant expansion at the pathway level and wide homoeologue expression dominance in the B-subgenome of the allotriploid. The independent origin of 1-aminocyclopropane-1-carboxylic acid oxidase (ACO) homoeologue gene pairs and tandem duplication-driven expansion of ACO genes in the B-subgenome contributed to rapid and major ethylene production post-harvest in allotriploid banana fruits. The findings of this study provide greater context for understanding fruit biology, and aid the development of tools for breeding optimal banana cultivars.

Metabolite profiling characterises chemotypes of Musa diploids and triploids at juvenile and pre-flowering growth stages

Bananas (Musa spp.) are consumed worldwide as dessert and cooking types. Edible banana varieties are for the most part seedless and sterile and therefore vegetatively propagated. This confers difficulties for breeding approaches against pressing biotic and abiotic threats and for the nutritional enhancement of banana pulp. A panel of banana accessions, representative of the diversity of wild and cultivated bananas, was analysed to assess the range of chemotypes available globally. The focus of this assessment was banana leaves at two growth stages (juvenile and pre-flowering), to see when during the plant growth metabolic differences can be established. The metabolic data corresponded to genomic trends reported in previous studies and demonstrated a link between metabolites/pathways and the genomes of M. acuminata and M. balbisiana. Furthermore, the vigour and resistance traits of M. balbisiana was connected to the phenolic composition and showed differences with the number of B genes in the hybrid accessions. Differences in the juvenile and pre-flowering data led to low correlation between the growth stages for prediction purposes.

Balance design for robust foliar nutrient diagnosis of "Prata" banana (Musa spp.)

The "Cavendish" and "Prata" subgroups represent respectively 47% and 24% of the world banana production. Compared to world average progressing from 10.6 to 20.6 t ha-1 between 1961 and 2016, and despite sustained domestic demand and the introduction of new cultivars, banana yield in Brazil has stagnated around 14.5 t ha-1 mainly due to nutrient and water mismanagement. "Prata" is now the dominant subgroup in N-E Brazil and is fertigated at high costs. Nutrient balances computed as isometric log-ratios (ilr) provide a comprehensive understanding of nutrient relationships in the diagnostic leaf at high yield level by combining raw concentration data. Although the most appropriate method for multivariate analysis of compositional balances may be less efficient due to non-normal data distribution and limited nutrient mobility in the plant, robustness of the nutrient balance approach could be improved using Box-Cox exponents assigned to raw foliar concentrations. Our objective was to evaluate the accuracy of nutrient balances to diagnose fertigated "Prata" orchards. The dataset comprised 609 observations on fruit yields and leaf tissue compositions collected from 2010 to 2016 in Ceará state, N-E Brazil. Raw nutrient concentration ranges were ineffective as diagnostic tool due to considerable overlapping of concentration ranges for low- and high-yielding subpopulations at cutoff yield of 40 Mg ha-1. Nutrient concentrations were combined into isometric log-ratios (ilr) and normalized by Box-Cox corrections between 0 and 1 which may also account for restricted nutrient transfer from leaf to fruit. Despite reduced ilr skewness, Box-Cox coefficients did not improve model robustness measured as the accuracy of the Cate-Nelson partition between yield and the multivariate distance across ilr values. Sensitivity was 94%, indicating that low yields are attributable primarily to nutrient imbalance. There were 148 false-positive specimens (high yield despite nutrient imbalance) likely due to suboptimal nutrition, contamination, or luxury consumption. The profitability of "Prata" orchards could be enhanced by rebalancing nutrients using ilr standards with no need for Box-Cox correction.

Comparative analysis of genetic diversity and differentiation of cauliflower (Brassica oleracea var. botrytis) accessions from two ex situ genebanks

Cauliflower (Brassica oleracea var. botrytis) is an important vegetable crop for human nutrition. We characterized 192 cauliflower accessions from the USDA and IPK genebanks with genotyping by sequencing (GBS). They originated from 26 different countries and represent about 44% of all cauliflower accessions in both genebanks. The analysis of genetic diversity revealed that accessions formed two major groups that represented the two genebanks and were not related to the country of origin. This differentiation was robust with respect to the analysis methods that included principal component analysis, ADMIXTURE and neighbor-joining trees. Genetic diversity was higher in the USDA collection and significant phenotypic differences between the two genebanks were found in three out of six traits investigated. GBS data have a high proportion of missing data, but we observed that the exclusion of single nucleotide polymorphisms (SNPs) with missing data or the imputation of missing SNP alleles produced very similar results. The results indicate that the composition and type of accessions have a strong effect on the structure of genetic diversity of ex situ collections, although regeneration procedures and local adaptation to regeneration conditions may also contribute to a divergence. F_st-based outlier tests of genetic differentiation identified only a small proportion (<1%) of SNPs that are highly differentiated between the two genebanks, which indicates that selection during seed regeneration is not a major cause of differentiation between genebanks. Seed regeneration procedures of both genebanks do not result in different levels of genetic drift and loss of genetic variation. We therefore conclude that the composition and type of accessions mainly influence the level of genetic diversity and explain the strong genetic differentiation between the two ex situ collections. In summary, GBS is a useful method for characterizing genetic diversity in cauliflower genebank material and our results suggest that it may be useful to incorporate routine genotyping into accession management and seed regeneration to monitor the diversity present in ex situ collections and to reduce the loss of genetic diversity during seed regeneration.

Genome survey of pistachio (Pistacia vera L.) by next generation sequencing: Development of novel SSR markers and genetic diversity in Pistacia species

BACKGROUND

Pistachio (Pistacia vera L.) is one of the most important nut crops in the world. There are about 11 wild species in the genus Pistacia, and they have importance as rootstock seed sources for cultivated P. vera and forest trees. Published information on the pistachio genome is limited. Therefore, a genome survey is necessary to obtain knowledge on the genome structure of pistachio by next generation sequencing. Simple sequence repeat (SSR) markers are useful tools for germplasm characterization, genetic diversity analysis, and genetic linkage mapping, and may help to elucidate genetic relationships among pistachio cultivars and species.

RESULTS

To explore the genome structure of pistachio, a genome survey was performed using the Illumina platform at approximately 40× coverage depth in the P. vera cv. Siirt. The K-mer analysis indicated that pistachio has a genome that is about 600 Mb in size and is highly heterozygous. The assembly of 26.77 Gb Illumina data produced 27,069 scaffolds at N50 = 3.4 kb with a total of 513.5 Mb. A total of 59,280 SSR motifs were detected with a frequency of 8.67 kb. A total of 206 SSRs were used to characterize 24 P. vera cultivars and 20 wild Pistacia genotypes (four genotypes from each five wild Pistacia species) belonging to P. atlantica, P. integerrima, P. chinenesis, P. terebinthus, and P. lentiscus genotypes. Overall 135 SSR loci amplified in all 44 cultivars and genotypes, 41 were polymorphic in six Pistacia species. The novel SSR loci developed from cultivated pistachio were highly transferable to wild Pistacia species.

CONCLUSIONS

The results from a genome survey of pistachio suggest that the genome size of pistachio is about 600 Mb with a high heterozygosity rate. This information will help to design whole genome sequencing strategies for pistachio. The newly developed novel polymorphic SSRs in this study may help germplasm characterization, genetic diversity, and genetic linkage mapping studies in the genus Pistacia.

DArT whole genome profiling provides insights on the evolution and taxonomy of edible Banana (Musa spp.)

BACKGROUND AND AIMS

Dessert and cooking bananas are vegetatively propagated crops of great importance for both the subsistence and the livelihood of people in developing countries. A wide diversity of diploid and triploid cultivars including AA, AB, AS, AT, AAA, AAB, ABB, AAS and AAT genomic constitutions exists. Within each of this genome groups, cultivars are classified into subgroups that are reported to correspond to varieties clonally derived from each other after a single sexual event. The number of those founding events at the basis of the diversity of bananas is a matter of debate.

METHODS

We analysed a large panel of 575 accessions, 94 wild relatives and 481 cultivated accessions belonging to the section Musa with a set of 498 DArT markers previously developed.

KEY RESULTS

DArT appeared successful and accurate to describe Musa diversity and help in the resolution of cultivated banana genome constitution and taxonomy, and highlighted discrepancies in the acknowledged classification of some accessions. This study also argues for at least two centres of domestication corresponding to South-East Asia and New Guinea, respectively. Banana domestication in New Guinea probably followed different schemes that those previously reported where hybridization underpins the emergence of edible banana. In addition, our results suggest that not all wild ancestors of bananas are known, especially in M. acuminata subspecies. We also estimate the extent of the two consecutive bottlenecks in edible bananas by evaluating the number of sexual founding events underlying our sets of edible diploids and triploids, respectively.

CONCLUSIONS

The attribution of clone identity to each sample of the sets allowed the detection of subgroups represented by several sets of clones. Although morphological characterization of some of the accessions is needed to correct potentially erroneous classifications, some of the subgroups seem polyclonal.

How endogenous plant pararetroviruses shed light on Musa evolution

BACKGROUND AND AIMS

Banana genomes harbour numerous copies of viral sequences derived from banana streak viruses (BSVs) - dsDNA viruses belonging to the family Caulimoviridae.These viral integrants (eBSVs) are mostly defective, probably as a result of 'pseudogenization' driven by host genome evolution. However, some can give rise to infection by releasing a functional viral genome following abiotic stresses. These distinct infective eBSVs correspond to the three main widespread BSV species (BSOLV, BSGFV and BSIMV), fully described within the Musa balbisiana B genomes of the seedy diploid 'Pisang Klutuk Wulung' (PKW).

METHODS

We characterize eBSV distribution among a Musa sampling including seedy BB diploids and interspecific hybrids with Musa acuminate exhibiting different levels of ploidy for the B genome (ABB, AAB, AB). We used representative samples of the two areas of sympatry between M. acuminate and M. balbisiana species representing the native area of the most widely cultivated AAB cultivars (in India and in East Asia, ranging from the Philippines to New Guinea). Seventy-seven accessions were characterized using eBSV-related PCR markers and Southern hybridization approaches. We coded both sets of results to create a common dissimilarity matrix with which to interpret eBSV distribution.

KEY RESULTS

We propose a Musa phylogeny driven by the M. balbisiana genome based on a dendrogram resulting from a joint neighbour-joining analysis of the three BSV species, showing for the first time lineages between BB and ABB/AAB hybrids. eBSVs appear to be relevant phylogenetic markers that can illustrate theM. balbisiana phylogeography story.

CONCLUSION

The theoretical implications of this study for further elucidation of the historical and geographical process of Musa domestication are numerous. Discovery of banana plants with B genome non-infective for eBSV opens the way to the introduction of new genitors in programmes of genetic banana improvement.

The triploid East African Highland Banana (EAHB) genepool is genetically uniform arising from a single ancestral clone that underwent population expansion by vegetative propagation

All East African Highland Banana varieties are genetically uniform having arisen from a single clone introduced to Africa. East African Highland bananas (EAHBs) are a subgroup of triploid (AAA genome) bananas of importance to food security in the Great Lakes region of Africa. Little is known about their genetic variation, population structure and evolutionary history. Ninety phenotypically diverse EAHB cultivars were genotyped at 100 SSR microsatellite markers to investigate population genetic diversity, the correlation of genetic variability with morphological classes, and evolutionary origins since introduction to Africa. Population-level statistics were compared to those for plantain (AAB) and dessert (AAA) cultivars representing other M. acuminata subgroups. EAHBs displayed minimal genetic variation and are largely genetically uniform, irrespective of whether they were derived from the distinct Ugandan or Kenyan germplasm collections. No association was observed between EAHB genetic diversity and currently employed morphological taxonomic systems for EAHB germplasm. Population size dynamics indicated that triploid EAHBs arose as a single hybridization event, which generated a genetic bottleneck during foundation of the EAHB genepool. As EAHB triploids are sterile, subsequent asexual vegetative propagation of EAHBs allowed a recent rapid expansion in population size. This provided a basis for emergence of genetically near-isogenic somatic mutants selected across farmers and environments in East Africa over the past 2000 years since EAHBs were first introduced to the African continent.

The striking and unexpected cytogenetic diversity of genus Tanacetum L. (Asteraceae): a cytometric and fluorescent in situ hybridisation study of Iranian taxa

BACKGROUND

Although karyologically well studied, the genus Tanacetum (Asteraceae) is poorly known from the perspective of molecular cytogenetics. The prevalence of polyploidy, including odd ploidy warranted an extensive cytogenetic study. We studied several species native to Iran, one of the most important centres of diversity of the genus. We aimed to characterise Tanacetum genomes through fluorochrome banding, fluorescent in situ hybridisation (FISH) of rRNA genes and the assessment of genome size by flow cytometry. We appraise the effect of polyploidy and evaluate the existence of intraspecific variation based on the number and distribution of GC-rich bands and rDNA loci. Finally, we infer ancestral genome size and other cytogenetic traits considering phylogenetic relationships within the genus.

RESULTS

We report first genome size (2C) estimates ranging from 3.84 to 24.87 pg representing about 11 % of those recognised for the genus. We found striking cytogenetic diversity both in the number of GC-rich bands and rDNA loci. There is variation even at the population level and some species have undergone massive heterochromatic or rDNA amplification. Certain morphometric data, such as pollen size or inflorescence architecture, bear some relationship with genome size. Reconstruction of ancestral genome size, number of CMA+ bands and number of rDNA loci show that ups and downs have occurred during the evolution of these traits, although genome size has mostly increased and the number of CMA+ bands and rDNA loci have decreased in present-day taxa compared with ancestral values.

CONCLUSIONS

Tanacetum genomes are highly unstable in the number of GC-rich bands and rDNA loci, although some patterns can be established at the diploid and tetraploid levels. In particular, aneuploid taxa and some odd ploidy species show greater cytogenetic instability than the rest of the genus. We have also confirmed a linked rDNA arrangement for all the studied Tanacetum species. The labile scenario found in Tanacetum proves that some cytogenetic features previously regarded as relatively constant, or even diagnostic, can display high variability, which is better interpreted within a phylogenetic context.

Genetic diversity in aspen and its relation to arthropod abundance

The ecological consequences of biodiversity have become a prominent public issue. Little is known on the effect of genetic diversity on ecosystem services. Here, a diversity experiment was established with European and North American aspen (Populus tremula, P. tremuloides) planted in plots representing either a single deme only or combinations of two, four and eight demes. The goals of this study were to explore the complex inter- and intraspecific genetic diversity of aspen and to then relate three measures for diversity (deme diversity, genetic diversity determined as Shannon index or as expected heterozygosity) to arthropod abundance. Microsatellite and AFLP markers were used to analyze the genetic variation patterns within and between the aspen demes and deme mixtures. Large differences were observed regarding the genetic diversity within demes. An analysis of molecular variance revealed that most of the total genetic diversity was found within demes, but the genetic differentiation among demes was also high. The complex patterns of genetic diversity and differentiation resulted in large differences of the genetic variation within plots. The average diversity increased from plots with only one deme to plots with two, four, and eight demes, respectively and separated plots with and without American aspen. To test whether intra- and interspecific diversity impacts on ecosystem services, arthropod abundance was determined. Increasing genetic diversity of aspen was related to increasing abundance of arthropods. However, the relationship was mainly driven by the presence of American aspen suggesting that species identity overrode the effect of intraspecific variation of European aspen.

Characterization of Musa sp. fruits and plantain banana ripening stages according to their physicochemical attributes

This study aimed at understanding the contribution of the fruit physicochemical parameters to Musa sp. diversity and plantain ripening stages. A discriminant analysis was first performed on a collection of 35 Musa sp. cultivars, organized in six groups based on the consumption mode (dessert or cooking banana) and the genomic constitution. A principal component analysis reinforced by a logistic regression on plantain cultivars was proposed as an analytical approach to describe the plantain ripening stages. The results of the discriminant analysis showed that edible fraction, peel pH, pulp water content, and pulp total phenolics were among the most contributing attributes for the discrimination of the cultivar groups. With mean values ranging from 65.4 to 247.3 mg of gallic acid equivalents/100 g of fresh weight, the pulp total phenolics strongly differed between interspecific and monospecific cultivars within dessert and nonplantain cooking bananas. The results of the logistic regression revealed that the best models according to fitting parameters involved more than one physicochemical attribute. Interestingly, pulp and peel total phenolic contents contributed in the building up of these models.

Genome-wide characterization of genetic diversity and population structure in Secale

BACKGROUND

Numerous rye accessions are stored in ex situ genebanks worldwide. Little is known about the extent of genetic diversity contained in any of them and its relation to contemporary varieties, since to date rye genetic diversity studies had a very limited scope, analyzing few loci and/ or few accessions. Development of high throughput genotyping methods for rye opened the possibility for genome wide characterizations of large accessions sets. In this study we used 1054 Diversity Array Technology (DArT) markers with defined chromosomal location to characterize genetic diversity and population structure in a collection of 379 rye accessions including wild species, landraces, cultivated materials, historical and contemporary rye varieties.

RESULTS

Average genetic similarity (GS) coefficients and average polymorphic information content (PIC) values varied among chromosomes. Comparison of chromosome specific average GS within and between germplasm sub-groups indicated regions of chromosomes 1R and 4R as being targeted by selection in current breeding programs. Bayesian clustering, principal coordinate analysis and Neighbor Joining clustering demonstrated that source and improvement status contributed significantly to the structure observed in the analyzed set of Secale germplasm. We revealed a relatively limited diversity in improved rye accessions, both historical and contemporary, as well as lack of correlation between clustering of improved accessions and geographic origin, suggesting common genetic background of rye accessions from diverse geographic regions and extensive germplasm exchange. Moreover, contemporary varieties were distinct from the remaining accessions.

CONCLUSIONS

Our results point to an influence of reproduction methods on the observed diversity patterns and indicate potential of ex situ collections for broadening the genetic diversity in rye breeding programs. Obtained data show that DArT markers provide a realistic picture of the genetic diversity and population structure present in the collection of 379 rye accessions and are an effective platform for rye germplasm characterization and association mapping studies.