Characterization of novel microsatellite markers in Musa acuminata subsp. burmannicoides, var. Calcutta 4

Transcriptome Profiling of the Resistance Response of Musa acuminata subsp. burmannicoides, var. Calcutta 4 to Pseudocercospora musae

Banana (Musa spp.), which is one of the world’s most popular and most traded fruits, is highly susceptible to pests and diseases. Pseudocercospora musae, responsible for Sigatoka leaf spot disease, is a principal fungal pathogen of Musa spp., resulting in serious economic damage to cultivars in the Cavendish subgroup. The aim of this study was to characterize genetic components of the early immune response to P. musae in Musa acuminata subsp. burmannicoides, var. Calcutta 4, a resistant wild diploid. Leaf RNA samples were extracted from Calcutta 4 three days after inoculation with fungal conidiospores, with paired-end sequencing conducted in inoculated and non-inoculated controls using lllumina HiSeq 4000 technology. Following mapping to the reference M. acuminata ssp. malaccensis var. Pahang genome, differentially expressed genes (DEGs) were identified and expression representation analyzed on the basis of gene ontology enrichment, Kyoto Encyclopedia of Genes and Genomes orthology and MapMan pathway analysis. Sequence data mapped to 29,757 gene transcript models in the reference Musa genome. A total of 1073 DEGs were identified in pathogen-inoculated cDNA libraries, in comparison to non-inoculated controls, with 32% overexpressed. GO enrichment analysis revealed common assignment to terms that included chitin binding, chitinase activity, pattern binding, oxidoreductase activity and transcription factor (TF) activity. Allocation to KEGG pathways revealed DEGs associated with environmental information processing, signaling, biosynthesis of secondary metabolites, and metabolism of terpenoids and polyketides. With 144 up-regulated DEGs potentially involved in biotic stress response pathways, including genes involved in cell wall reinforcement, PTI responses, TF regulation, phytohormone signaling and secondary metabolism, data demonstrated diverse early-stage defense responses to P. musae. With increased understanding of the defense responses occurring during the incompatible interaction in resistant Calcutta 4, these data are appropriate for the development of effective disease management approaches based on genetic improvement through introgression of candidate genes in superior cultivars.

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.

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.

Genome-Wide Computational Analysis of Musa Microsatellites: Classification, Cross-Taxon Transferability, Functional Annotation, Association with Transposons & miRNAs, and Genetic Marker Potential

The development of organized, informative, robust, user-friendly, and freely accessible molecular markers is imperative to the Musa marker assisted breeding program. Although several hundred SSR markers have already been developed, the number of informative, robust, and freely accessible Musa markers remains inadequate for some breeding applications. In view of this issue, we surveyed SSRs in four different data sets, developed large-scale non-redundant highly informative therapeutic SSR markers, and classified them according to their attributes, as well as analyzed their cross-taxon transferability and utility for the genetic study of Musa and its relatives. A high SSR frequency (177 per Mbp) was found in the Musa genome. AT-rich dinucleotide repeats are predominant, and trinucleotide repeats are the most abundant in transcribed regions. A significant number of Musa SSRs are associated with pre-miRNAs, and 83% of these SSRs are promising candidates for the development of therapeutic SSR markers. Overall, 74% of the SSR markers were polymorphic, and 94% were transferable to at least one Musa spp. Two hundred forty-three markers generated a total of 1047 alleles, with 2-8 alleles each and an average of 4.38 alleles per locus. The PIC values ranged from 0.31 to 0.89 and averaged 0.71. We report the largest set of non-redundant, polymorphic, new SSR markers to be developed in Musa. These additional markers could be a valuable resource for marker-assisted breeding, genetic diversity and genomic studies of Musa and related species.

Analysis of the leaf transcriptome of Musa acuminata during interaction with Mycosphaerella musicola: gene assembly, annotation and marker development

BACKGROUND

Although banana (Musa sp.) is an important edible crop, contributing towards poverty alleviation and food security, limited transcriptome datasets are available for use in accelerated molecular-based breeding in this genus. 454 GS-FLX Titanium technology was employed to determine the sequence of gene transcripts in genotypes of Musa acuminata ssp. burmannicoides Calcutta 4 and M. acuminata subgroup Cavendish cv. Grande Naine, contrasting in resistance to the fungal pathogen Mycosphaerella musicola, causal organism of Sigatoka leaf spot disease. To enrich for transcripts under biotic stress responses, full length-enriched cDNA libraries were prepared from whole plant leaf materials, both uninfected and artificially challenged with pathogen conidiospores.

RESULTS

The study generated 846,762 high quality sequence reads, with an average length of 334 bp and totalling 283 Mbp. De novo assembly generated 36,384 and 35,269 unigene sequences for M. acuminata Calcutta 4 and Cavendish Grande Naine, respectively. A total of 64.4% of the unigenes were annotated through Basic Local Alignment Search Tool (BLAST) similarity analyses against public databases.Assembled sequences were functionally mapped to Gene Ontology (GO) terms, with unigene functions covering a diverse range of molecular functions, biological processes and cellular components. Genes from a number of defense-related pathways were observed in transcripts from each cDNA library. Over 99% of contig unigenes mapped to exon regions in the reference M. acuminata DH Pahang whole genome sequence. A total of 4068 genic-SSR loci were identified in Calcutta 4 and 4095 in Cavendish Grande Naine. A subset of 95 potential defense-related gene-derived simple sequence repeat (SSR) loci were validated for specific amplification and polymorphism across M. acuminata accessions. Fourteen loci were polymorphic, with alleles per polymorphic locus ranging from 3 to 8 and polymorphism information content ranging from 0.34 to 0.82.

CONCLUSIONS

A large set of unigenes were characterized in this study for both M. acuminata Calcutta 4 and Cavendish Grande Naine, increasing the number of public domain Musa ESTs. This transcriptome is an invaluable resource for furthering our understanding of biological processes elicited during biotic stresses in Musa. Gene-based markers will facilitate molecular breeding strategies, forming the basis of genetic linkage mapping and analysis of quantitative trait loci.

Exploring wild genetic resources of Musa acuminata Colla distributed in the humid forests of southern Western Ghats of peninsular India using ISSR markers

Musa acuminata ssp. burmannica, one of the wild progenitors contributing 'A genome' to the present-day dessert bananas, has a long evolutionary history intervened by human activities. In this study, ISSR markers were used to analyze the pattern of genetic variation and differentiation in 32 individuals along with two reference samples (viz., Musa acuminata ssp. burmannicoides, var. Calcutta 4 and Musa balbisiana) of wild Musa, which corresponded to three populations across the biodiversity-rich hot spot of southern Western Ghats of India. High levels of genetic diversity were revealed both at the species and population levels, using Nei's diversity indices. The hierarchical analysis of molecular variance showed pronounced genetic differentiation, as 96% of the total variance was fixed within population and only 4% among populations. Nei's genetic differentiation coefficient (GST=0.1823) and low gene flow (Nm=1.18) further confirmed this. The positive correlation (Mantel test) between geographic distance and genetic distance (r=0.338 P<0.001) indicates geographic isolation as one of the key factors in shaping the population genetic structure. Grouping of individuals was largely in conformity with their spatial distribution, which was confirmed by UPGMA cluster analysis and PCA scatter plot clustering all 32 individuals into three major groups along a geographical gradient. The discontinuous distribution and dwindling population due to habitat fragmentation are serious threats to prevailing genetic diversity in this species. Conservation measures based on diversity pattern are suggested for long-term preservation and sustainable utilization of this precious genetic resource.

KEY MESSAGE

A diverse germplasm of Musa acuminata ssp. burmannica exists in southern Western Ghats as a possible repository of useful resistant traits, which can be effectively utilized for crop improvement.

New microsatellite markers for bananas (Musa spp)

Thirty-four microsatellite markers (SSRs) were identified in EST and BAC clones from Musa acuminata burmannicoides var. Calcutta 4 and validated in 22 Musa genotypes from the Banana Germplasm Bank of Embrapa-CNPMF, which includes wild and improved diploids. The number of alleles per locus ranged from 2 to 14. The markers were considered highly informative based on their polymorphism information content values; more than 50% were above 0.5. These SSRs will be useful for banana breeding programs, for studies of genetic diversity, germplasm characterization and selection, development of saturated genetic linkage maps, and marker assisted selection.

Development of expressed sequence tag and expressed sequence tag-simple sequence repeat marker resources for Musa acuminata

BACKGROUND AND AIMS

Banana (Musa acuminata) is a crop contributing to global food security. Many varieties lack resistance to biotic stresses, due to sterility and narrow genetic background. The objective of this study was to develop an expressed sequence tag (EST) database of transcripts expressed during compatible and incompatible banana-Mycosphaerella fijiensis (Mf) interactions. Black leaf streak disease (BLSD), caused by Mf, is a destructive disease of banana. Microsatellite markers were developed as a resource for crop improvement.

METHODOLOGY

cDNA libraries were constructed from in vitro-infected leaves from BLSD-resistant M. acuminata ssp. burmaniccoides Calcutta 4 (MAC4) and susceptible M. acuminata cv. Cavendish Grande Naine (MACV). Clones were 5'-end Sanger sequenced, ESTs assembled with TGICL and unigenes annotated using BLAST, Blast2GO and InterProScan. Mreps was used to screen for simple sequence repeats (SSRs), with markers evaluated for polymorphism using 20 diploid (AA) M. acuminata accessions contrasting in resistance to Mycosphaerella leaf spot diseases.

PRINCIPAL RESULTS

A total of 9333 high-quality ESTs were obtained for MAC4 and 3964 for MACV, which assembled into 3995 unigenes. Of these, 2592 displayed homology to genes encoding proteins with known or putative function, and 266 to genes encoding proteins with unknown function. Gene ontology (GO) classification identified 543 GO terms, 2300 unigenes were assigned to EuKaryotic orthologous group categories and 312 mapped to Kyoto Encyclopedia of Genes and Genomes pathways. A total of 624 SSR loci were identified, with trinucleotide repeat motifs the most abundant in MAC4 (54.1 %) and MACV (57.6 %). Polymorphism across M. acuminata accessions was observed with 75 markers. Alleles per polymorphic locus ranged from 2 to 8, totalling 289. The polymorphism information content ranged from 0.08 to 0.81.

CONCLUSIONS

This EST collection offers a resource for studying functional genes, including transcripts expressed in banana-Mf interactions. Markers are applicable for genetic mapping, diversity characterization and marker-assisted breeding.