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

Characterization of microRNAs and Target Genes in Musa acuminata subsp. burmannicoides, var. Calcutta 4 during Interaction with Pseudocercospora musae

Endogenous microRNAs (miRNAs) are small non-coding RNAs that perform post-transcriptional regulatory roles across diverse cellular processes, including defence responses to biotic stresses. Pseudocercospora musae, the causal agent of Sigatoka leaf spot disease in banana (Musa spp.), is an important fungal pathogen of the plant. Illumina HiSeq 2500 sequencing of small RNA libraries derived from leaf material in Musa acuminata subsp. burmannicoides, var. Calcutta 4 (resistant) after inoculation with fungal conidiospores and equivalent non-inoculated controls revealed 202 conserved miRNAs from 30 miR-families together with 24 predicted novel miRNAs. Conserved members included those from families miRNA156, miRNA166, miRNA171, miRNA396, miRNA167, miRNA172, miRNA160, miRNA164, miRNA168, miRNA159, miRNA169, miRNA393, miRNA535, miRNA482, miRNA2118, and miRNA397, all known to be involved in plant immune responses. Gene ontology (GO) analysis of gene targets indicated molecular activity terms related to defence responses that included nucleotide binding, oxidoreductase activity, and protein kinase activity. Biological process terms associated with defence included response to hormone and response to oxidative stress. DNA binding and transcription factor activity also indicated the involvement of miRNA target genes in the regulation of gene expression during defence responses. sRNA-seq expression data for miRNAs and RNAseq data for target genes were validated using stem-loop quantitative real-time PCR (qRT-PCR). For the 11 conserved miRNAs selected based on family abundance and known involvement in plant defence responses, the data revealed a frequent negative correlation of expression between miRNAs and target host genes. This examination provides novel information on miRNA-mediated host defence responses, applicable in genetic engineering for the control of Sigatoka leaf spot disease.

Genome-wide identification, characterization, and evolutionary analysis of NBS genes and their association with disease resistance in Musa spp

Banana is an important food crop that is susceptible to a wide range of pests and diseases that can reduce yield and quality. The primary objective of banana breeding programs is to increase disease resistance, which requires the identification of resistance (R) genes. Despite the fact that resistant sources have been identified in bananas, the genes, particularly the nucleotide-binding site (NBS) family, which play an important role in protecting plants against pathogens, have received little attention. As a result, this study included a thorough examination of the NBS disease resistance gene family's classification, phylogenetic analysis, genome organization, evolution, cis-elements, differential expression, regulation by microRNAs, and protein-protein interaction. A total of 116 and 43 putative NBS genes from M. acuminata and M. balbisiana, respectively, were identified and characterized, and were classified into seven sub-families. Structural analysis of NBS genes revealed the presence of signal peptides, their sub-cellular localization, molecular weight and pI. Eight commonly conserved motifs were found, and NBS genes were unevenly distributed across multiple chromosomes, with the majority of NBS genes being located in chr3 and chr1 of the A and B genomes, respectively. Tandem duplication occurrences have helped bananas' NBS genes spread throughout evolution. Transcriptome analysis of NBS genes revealed significant differences in expression between resistant and susceptible cultivars of fusarium wilt, eumusae leaf spot, root lesion nematode, and drought, implying that they can be used as candidate resistant genes. Ninety miRNAs were discovered to have targets in 104 NBS genes from the A genome, providing important insights into NBS gene expression regulation. Overall, this study offers a valuable genomic resource and understanding of the function and evolution of NBS genes in relation to rapidly evolving pathogens, as well as providing breeders with selection targets for fast-tracking breeding of banana varieties with more durable resistance to pathogens.

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.

Infection by Pseudocercospora musae leads to an early reprogramming of the Musa paradisiaca defense transcriptome

Deep sequencing technologies such as RNA sequencing can help unravel mechanisms governing defense or resistance responses in plant-pathogen interactions. Several studies have been carried out to investigate the transcriptomic changes in Musa germplasm against Yellow Sigatoka disease, but the defense response of Musa paradisiaca has not been investigated so far. We carried out transcriptome sequencing of M. paradisiaca var. Kachkal infected with the pathogen Pseudocercospora musae and found that a vast set of genes were upregulated while many genes were downregulated in the resistant cultivar as a result of infection. After transcriptome assembly and differential gene expression analysis, 429 upregulated and 156 downregulated genes were filtered out (considering fold change ± 2, p < 0.01). Functional annotation of the differentially expressed genes (DEGs) enriched the upregulated genes into 49 gene ontology (GO) classes of biological processes (BP), 20 classes of molecular function (MF) and 9 classes of cellular component (CC). Similarly, the downregulated genes were classified into 35 GO classes of BP, 28 classes of MF and 6 classes of CC. The KEGG enrichment analysis revealed that the upregulated genes were most highly represented in 'metabolic' and 'biosynthesis of secondary metabolites' pathways. Additionally, 'plant hormone signal transduction', 'plant-pathogen interaction' and 'phenylpropanoid biosynthesis' pathways were also significantly enriched indicating their involvement in resistance responses against the pathogen. The RNA-seq analysis also depicts that a range of important defense-related genes are modulated as a result of infection, all of which are responsible for either mediating or activating resistance responses in the host. We studied and validated the expression profiles of ten important defense-related genes potentially involved in conferring resistance to the pathogen through qRT-PCR. Almost all the selected defense-related genes were found to be highly and significantly upregulated within 24 h post inoculation (hpi) and for some genes, the expression remained consistently high till the later time point of 72 hpi. These results, thus, indicate that the infection by P. musae leads to a rapid reprogramming of the defense transcriptome of the resistant banana cultivar. The defense-related genes identified to be modulated in response to infection are important not only for pathogen recognition and perception but also for activation and persistence of defense in the host.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03245-9.

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.

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

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

Genetic Improvement for Resistance to Black Sigatoka in Bananas: A Systematic Review

Bananas are an important staple food crop in tropical and subtropical regions in Asia, sub-Saharan Africa, and Central and South America. The plant is affected by numerous diseases, with the fungal leaf disease black Sigatoka, caused by Mycosphaerella fijiensis Morelet [anamorph: Pseudocercospora fijiensis (Morelet) Deighton], considered one of the most economically important phytosanitary problem. Although the development of resistant cultivars is recognized as most effective method for long term control of the disease, the majority of today's cultivars are susceptible. In order to gain insights into this pathosystem, this first systematic literature review on the topic is presented. Utilizing six databases (PubMed Central, Web of Science, Google Academic, Springer, CAPES and Scopus Journals) searches were performed using pre-established inclusion and exclusion criteria. From a total of 3,070 published studies examined, 24 were relevant with regard to the Musa-P. fijiensis pathosystem. Relevant papers highlighted that resistant and susceptible cultivars clearly respond differently to infection by this pathogen. M. acuminata wild diploids such as Calcutta 4 and other diploid cultivars can harbor sources of resistance genes, serving as parentals for the generation of improved diploids and subsequent gene introgression in new cultivars. From the sequenced reference genome of Musa acuminata, although the function of many genes in the genome still require validation, on the basis of transcriptome, proteome and biochemical data, numerous candidate genes and molecules have been identified for further evaluation through genetic transformation and gene editing approaches. Genes identified in the resistance response have included those associated with jasmonic acid and ethylene signaling, transcription factors, phenylpropanoid pathways, antioxidants and pathogenesis-related proteins. Papers in this study also revealed gene-derived markers in Musa applicable for downstream application in marker assisted selection. The information gathered in this review furthers understanding of the immune response in Musa to the pathogen P. fijiensis and is relevant for genetic improvement programs for bananas and plantains for control of black Sigatoka.

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.

Deciphering the Omics of Plant-Microbe Interaction: Perspectives and New Insights

Introduction

Plants do not grow in isolation, rather they are hosts to a variety of microbes in their natural environments. While, few thrive in the plants for their own benefit, others may have a direct impact on plants in a symbiotic manner. Unraveling plant-microbe interactions is a critical component in recognizing the positive and negative impacts of microbes on plants. Also, by affecting the environment around plants, microbes may indirectly influence plants. The progress in sequencing technologies in the genomics era and several omics tools has accelerated in biological science. Studying the complex nature of plant-microbe interactions can offer several strategies to increase the productivity of plants in an environmentally friendly manner by providing better insights. This review brings forward the recent works performed in building omics strategies that decipher the interactions between plant-microbiome. At the same time, it further explores other associated mutually beneficial aspects of plant-microbe interactions such as plant growth promotion, nitrogen fixation, stress suppressions in crops and bioremediation; as well as provides better insights on metabolic interactions between microbes and plants through omics approaches. It also aims to explore advances in the study of Arabidopsis as an important avenue to serve as a baseline tool to create models that help in scrutinizing various factors that contribute to the elaborate relationship between plants and microbes. Causal relationships between plants and microbes can be established through systematic gnotobiotic experimental studies to test hypotheses on biologically derived interactions.

Conclusion

This review will cover recent advances in the study of plant-microbe interactions keeping in view the advantages of these interactions in improving nutrient uptake and plant health.

Stable reference genes for RT-qPCR analysis of gene expression in the Musa acuminata-Pseudocercospora musae interaction

Leaf pathogens are limiting factors in banana (Musa spp.) production, with Pseudocercospora spp. responsible for the important Sigatoka disease complex. In order to investigate cellular processes and genes involved in host defence responses, quantitative real-time PCR (RT-qPCR) is an analytical technique for gene expression quantification. Reliable RT-qPCR data, however, requires that reference genes for normalization of mRNA levels in samples are validated under the conditions employed for expression analysis of target genes. We evaluated the stability of potential reference genes ACT1, α-TUB, UBQ1, UBQ2, GAPDH, EF1α, APT and RAN. Total RNA was extracted from leaf tissues of Musa acuminata genotypes Calcutta 4 (resistant) and Cavendish Grande Naine (susceptible), both subjected to P. musae infection. Expression stability was determined with NormFinder, BestKeeper, geNorm and RefFinder algorithms. UBQ2 and RAN were the most stable across all M. acuminata samples, whereas when considering inoculated and non-inoculated leaf samples, APT and UBQ2 were appropriate for normalization in Calcutta 4, with RAN and α-TUB most stable in Cavendish Grande Naine. This first study of reference genes for relative quantification of target gene expression in the M. acuminata-P. musae interaction will enable reliable analysis of gene expression in this pathosystem, benefiting elucidation of disease resistance mechanisms.

MusatransSSRDB (a transcriptome derived SSR database) - An advanced tool for banana improvement

Availability of transcriptome datasets for use in accelerated molecular-based breeding in Musa species is limited. Illumina Hiseq technology was employed to determine differential gene expression between the contrasting cultivars for three different stresses (Eumusae leaf spot -Mycosphaerella eumusae, root lesion nematode - Pratylenchus coffeae and moisture deficit stress) under challenged and unchallenged conditions. An average of 34.72 million of reads was assembled into ~47629 contigs, and ~5,466 simple sequence repeats (SSR) from each library were identified. GO annotation and KEGG pathway analysis were carried for all the transcripts and the SSR, SNPs were also detected. Based on this information, a MusatransSSRDB has been developed. Currently, the database consists of 32,800 SSRs with the unique information like putative function of the SSR-containing genes and their metabolic pathway and expression profiling under various stress conditions. This database provides information on in silico polymorphic SSRs (2830 SSRs) between the contrasting cultivars for each stress and within stress. Information on in silico polymorphic SSRs specific to differentially expressed genes under challenged condition for each stress can also be accessed. This database facilitates the retrieval of results by navigating the tabs for cultivars, stress and polymorphism. This database was developed using HTML, Java and PHP; datasets are stored in MySQL database and accessible in the public domain (http://bioinfnrcb.byethost7.com/nrcbbio/). This unique information facilitates the banana breeder to select the SSR primers based on specific objectives. MusatransSSRDB along with other genomics databases will facilitate the genetic dissection and breeding for complex traits in banana. Thus, this database is a step forward in economizing cost, time, manpower and other resources. Keywords.

Development, characterization, and cross-amplification of microsatellite markers for Psammosilene tunicoides (Caryophyllaceae)

PREMISE OF THE STUDY

Psammosilene tunicoides (Caryophyllaceae) is a narrowly distributed and endemic plant species in southwestern China. The overexploitation of natural P. tunicoides has led to the destruction of many populations. Population and genetic studies will provide crucial data for the protection and management of P. tunicoides. In this study, we develop simple sequence repeat markers of P. tunicoides to analyze population diversity.

METHODS AND RESULTS

Microsatellite loci of P. tunicoides were isolated with FIASCO. Eleven polymorphic and 10 monomorphic primers were developed. The 11 polymorphic primers were tested in three P. tunicoides populations, yielding two to nine alleles per locus. Levels of observed heterozygosity varied from 0.000 to 1.000, and levels of expected heterozygosity ranged from 0.000 to 0.615. In addition, three of these loci were successfully amplified, and showed polymorphism, in three Silene species.

CONCLUSIONS

These microsatellite markers can be valuable tools to investigate the genetic diversity and population structure of P. tunicoides.

Current Strategies for the Detoxification of Jatropha curcas Seed Cake: A Review

Jatropha curcas is an important oilseed plant, with considerable potential in the development of biodiesel. Although Jatropha seed cake, the byproduct of oil extraction, is a residue rich in nitrogen, phosphorus, potassium, and carbon, with high protein content suitable for application in animal feed, the presence of toxic phorbol esters limits its application in feed supplements and fertilizers. This review summarizes the current methods available for detoxification of this residue, based upon chemical, physical, biological, or combined processes. The advantages and disadvantages of each process are discussed, and future directions involving genomic and proteomic approaches for advancing our understanding of biodegradation processes involving microorganisms are highlighted.

Gene expression analysis in Musa acuminata during compatible interactions with Meloidogyne incognita

BACKGROUND AND AIMS

Endoparasitic root-knot nematodes (RKNs) ( Meloidogyne spp.) cause considerable losses in banana ( Musa spp.), with Meloidogyne incognita a predominant species in Cavendish sub-group bananas. This study investigates the root transcriptome in Musa acuminata genotypes 4297-06 (AA) and Cavendish Grande Naine (CAV; AAA) during early compatible interactions with M. incognita .

METHODS

Roots were analysed by brightfield light microscopy over a 35 d period to examine nematode penetration and morphological cell transformation. RNA samples were extracted 3, 7 and 10 days after inoculation (DAI) with nematode J2 juveniles, and cDNA libraries were sequenced using lllumina HiSeq technology. Sequences were mapped to the M. acuminata ssp. malaccensis var. Pahang genome sequence, differentially expressed genes (DEGs) identified and transcript representation determined by gene set enrichment and pathway mapping.

KEY RESULTS

Microscopic analysis revealed a life cycle of M. incognita completing in 24 d in CAV and 27 d in 4279-06. Comparable numbers of DEGs were up- and downregulated in each genotype, with potential involvement of many in early host defence responses involving reactive oxygen species and jasmonate/ethylene signalling. DEGs revealed concomitant auxin metabolism and cell wall modification processes likely to be involved in giant cell formation. Notable transcripts related to host defence included those coding for leucine-rich repeat receptor-like serine/threonine-protein kinases, peroxidases, thaumatin-like pathogenesis-related proteins, and DREB, ERF, MYB, NAC and WRKY transcription factors. Transcripts related to giant cell development included indole acetic acid-amido synthetase GH3.8 genes, involved in auxin metabolism, as well as genes encoding expansins and hydrolases, involved in cell wall modification.

CONCLUSIONS

Expression analysis in M. acuminata during compatible interactions with RKNs provides insights into genes modulated during infection and giant cell formation. Increased understanding of both defence responses to limit parasitism during compatible interactions and effector-targeted host genes in this complex interaction will facilitate the development of genetic improvement measures for RKNs.

Transcriptome sequencing of Mycosphaerella fijiensis during association with Musa acuminata reveals candidate pathogenicity genes

BACKGROUND

Mycosphaerella fijiensis, causative agent of the black Sigatoka disease of banana, is considered the most economically damaging banana disease. Despite its importance, the genetics of pathogenicity are poorly understood. Previous studies have characterized polyketide pathways with possible roles in pathogenicity. To identify additional candidate pathogenicity genes, we compared the transcriptome of this fungus during the necrotrophic phase of infection with that during saprophytic growth in medium.

RESULTS

Transcriptome analysis was conducted, and the functions of differentially expressed genes were predicted by identifying conserved domains, Gene Ontology (GO) annotation and GO enrichment analysis, Carbohydrate-Active EnZymes (CAZy) annotation, and identification of genes encoding effector-like proteins. The analysis showed that genes commonly involved in secondary metabolism have higher expression in infected leaf tissue, including genes encoding cytochrome P450s, short-chain dehydrogenases, and oxidoreductases in the 2-oxoglutarate and Fe(II)-dependent oxygenase superfamily. Other pathogenicity-related genes with higher expression in infected leaf tissue include genes encoding salicylate hydroxylase-like proteins, hydrophobic surface binding proteins, CFEM domain-containing proteins, and genes encoding secreted cysteine-rich proteins characteristic of effectors. More genes encoding amino acid transporters, oligopeptide transporters, peptidases, proteases, proteinases, sugar transporters, and proteins containing Domain of Unknown Function (DUF) 3328 had higher expression in infected leaf tissue, while more genes encoding inhibitors of peptidases and proteinases had higher expression in medium. Sixteen gene clusters with higher expression in leaf tissue were identified including clusters for the synthesis of a non-ribosomal peptide. A cluster encoding a novel fusicoccane was also identified. Two putative dispensable scaffolds were identified with a large proportion of genes with higher expression in infected leaf tissue, suggesting that they may play a role in pathogenicity. For two other scaffolds, no transcripts were detected in either condition, and PCR assays support the hypothesis that at least one of these scaffolds corresponds to a dispensable chromosome that is not required for survival or pathogenicity.

CONCLUSIONS

Our study revealed major changes in the transcriptome of Mycosphaerella fijiensis, when associating with its host compared to during saprophytic growth in medium. This analysis identified putative pathogenicity genes and also provides support for the existence of dispensable chromosomes in this fungus.

Identification of Differentially-Expressed Genes in Response to Mycosphaerella fijiensis in the Resistant Musa Accession 'Calcutta-4' Using Suppression Subtractive Hybridization

Bananas and plantains are considered an important crop around the world. Banana production is affected by several constraints, of which Black Sigatoka Disease, caused by the fungus Mycosphaerella fijiensis, is considered one of the most important diseases in banana plantations. The banana accession ‘Calcutta-4’ has a natural resistance to Black Sigatoka; however, the fruit is not valuable for commercialization. Gene identification and expression studies in ‘Calcutta-4’ might reveal possible gene candidates for resistant to the disease and elucidate mechanisms for resistance. A subtracted cDNA library was generated from leaves after 6, 9 and 12 days inoculated with M. fijiensis conidia on greenhouse banana plants of the accession ‘Calcutta-4’. Bioinformatic analysis revealed 99 good quality sequences. Blast2go analysis revealed that 31% of the sequences could not be categorized and, according to the Biological Process Category, 32 and 28 ESTs are related to general metabolic and cellular processes, respectively; while 10 ESTs response to stimulus. Seven sequences were redundant and one was similar to genes that may be involved in pathogen resistance including the putative disease resistance protein RGA1. Genes encoding zinc finger domains were identified and may play an important role in pathogen resistance by inducing the expression of downstream genes. Expression analysis of four selected genes was performed using RT-qPCR during the early stage of the disease development at 6, 9, 12 and 15 days post inoculation showing a peak of up regulation at 9 or 12 days post inoculation. Three of the four genes showed an up-regulation of expression in ‘Calcutta-4’ when compared to ‘Williams’ after inoculation with M. fijiensis, suggesting a fine regulation of specific gene candidates that may lead to a resistance response. The genes identified in early responses in a plant-pathogen interaction may be relevant for the resistance response of ‘Calcutta-4’ to Black Sigatoka. Genes with different functions may play a role in plant response to the disease. The present study suggests a fine up regulation of these genes that might be needed to perform an incompatible interaction. Further gene functional studies need to be performed to validate their use as candidate resistance genes in susceptible banana cultivars.

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.

De novo transcriptome sequence assembly from coconut leaves and seeds with a focus on factors involved in RNA-directed DNA methylation

Coconut palm (Cocos nucifera) is a symbol of the tropics and a source of numerous edible and nonedible products of economic value. Despite its nutritional and industrial significance, coconut remains under-represented in public repositories for genomic and transcriptomic data. We report de novo transcript assembly from RNA-seq data and analysis of gene expression in seed tissues (embryo and endosperm) and leaves of a dwarf coconut variety. Assembly of 10 GB sequencing data for each tissue resulted in 58,211 total unigenes in embryo, 61,152 in endosperm, and 33,446 in leaf. Within each unigene pool, 24,857 could be annotated in embryo, 29,731 could be annotated in endosperm, and 26,064 could be annotated in leaf. A KEGG analysis identified 138, 138, and 139 pathways, respectively, in transcriptomes of embryo, endosperm, and leaf tissues. Given the extraordinarily large size of coconut seeds and the importance of small RNA-mediated epigenetic regulation during seed development in model plants, we used homology searches to identify putative homologs of factors required for RNA-directed DNA methylation in coconut. The findings suggest that RNA-directed DNA methylation is important during coconut seed development, particularly in maturing endosperm. This dataset will expand the genomics resources available for coconut and provide a foundation for more detailed analyses that may assist molecular breeding strategies aimed at improving this major tropical crop.

Identification of glutathione S-transferase genes responding to pathogen infestation in Populus tomentosa

Stem blister canker, caused by Botryosphaeria dothidea, is becoming the most serious disease of poplar in China. The molecular basis of the poplar in response to stem blister canker is not well understood. To reveal the global transcriptional changes of poplar to infection by B. dothidea, Solexa paired-end sequencing of complementary DNAs (cDNAs) from control (NB) and pathogen-treated samples (WB) was performed, resulting in a total of 339,283 transcripts and 183,881 unigenes. A total of 206,586 transcripts were differentially expressed in response to pathogen stress (false discovery rate ≤0.05 and an absolute value of log2Ratio (NB/WB) ≥1). In enrichment analysis, energy metabolism and redox reaction-related macromolecules were accumulated significantly in Gene Ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways analyses, indicating components of dynamic defense against the fungus. A total of 852 transcripts (575 upregulated and 277 downregulated transcripts) potentially involved in plant-pathogen interaction were also differentially regulated, including genes encoding proteins linked to signal transduction (putative leucine-rich repeat (LRR) protein kinases and calcium-binding proteins), defense (pathogenesis-related protein 1), and cofactors (jasmonate-ZIM-domain-containing proteins and heat shock proteins). Moreover, transcripts encoding glutathione S-transferase (GST) were accumulated to high levels, revealing key genes and proteins potentially related to pathogen resistance. Poplar RNA sequence data were validated by quantitative real-time PCR (RT-qPCR), which revealed a highly reliability of the transcriptomic profiling data.

Characterization of Withania somnifera leaf transcriptome and expression analysis of pathogenesis-related genes during salicylic acid signaling

Withania somnifera (L.) Dunal is a valued medicinal plant with pharmaceutical applications. The present study was undertaken to analyze the salicylic acid induced leaf transcriptome of W. somnifera. A total of 45.6 million reads were generated and the de novo assembly yielded 73,523 transcript contig with average transcript contig length of 1620 bp. A total of 71,062 transcripts were annotated and 53,424 of them were assigned GO terms. Mapping of transcript contigs to biological pathways revealed presence of 182 pathways. Seventeen genes representing 12 pathogenesis-related (PR) families were mined from the transcriptome data and their pattern of expression post 17 and 36 hours of salicylic acid treatment was documented. The analysis revealed significant up-regulation of all families of PR genes by 36 hours post treatment except WsPR10. The relative fold expression of transcripts ranged from 1 fold to 6,532 fold. The two families of peroxidases including the lignin-forming anionic peroxidase (WsL-PRX) and suberization-associated anionic peroxidase (WsS-PRX) recorded maximum expression of 377 fold and 6532 fold respectively, while the expression of WsPR10 was down-regulated by 14 fold. Additionally, the most stable reference gene for normalization of qRT-PCR data was also identified. The effect of SA on the accumulation of major secondary metabolites of W. somnifera including withanoside V, withaferin A and withanolide A was also analyzed and an increase in content of all the three metabolites were detected. This is the first report on expression patterns of PR genes during salicylic acid signaling in W. somnifera.