Insights into the Musa genome: syntenic relationships to rice and between Musa species

Unveiling the dynamic expression of PR-1 during Musa spp. infection by Fusarium oxysporum fsp. Cubense: a cloning and characterization study

BACKGROUND

Pathogen-related proteins (PR) are pivotal in plant defense, combating diverse biotic and abiotic stresses. While multiple gene families contribute to banana resistance against Fusarium oxysporum f sp. cubense (Foc), Pseudocercospora eumusae, and Pratylenchus coffeae, the significance of PR-1 genes in defense is paramount.

METHODS

Three PR-1 genes, up-regulated under diverse biotic stresses, were cloned from both resistant and susceptible cultivars of Foc, P. eumusae, and P. coffeae. Molecular characterization, phylogenetic analysis, and docking studies with the Foc TR4 CP gene were conducted.

RESULTS

Through transcriptomic and real-time studies, three PR-1 genes (Ma02_g15050, Ma02_g15060, and Ma04_g34800) from Musa spp. were identified. These genes exhibited significant up-regulation in resistant cultivars when exposed to Foc, P. eumusae, and P. coffeae. Cloning of these genes was successfully performed from both resistant and susceptible cultivars of Foc race 1 and TR4, P. eumusae, and P. coffeae. Distinct characteristics were observed among the PR-1 genes, with groups 1 and 2 being acidic with signal peptides, and group 3 being basic without signal peptides. All cloned PR-1 proteins belonged to the CAP superfamily (PF00188). Phylogenetic analysis revealed clustering patterns for acidic PR-1 proteins, and KEGG orthology showed associations with vital pathways, including MAPK signaling, plant hormone signal transduction, and plant-pathogen interaction. Secondary and tertiary structure analyses confirmed sequence conservation across studied species. Docking studies explored interactions between the cerato-platanin (CP) gene from Foc TR4 and Ma02_g15060 from banana, suggesting the potential hindrance of PR-1 antifungal activity through direct interaction.

CONCLUSIONS

The findings underscore the crucial role of cloned PR-1 genes in banana plant defense mechanisms against a broad spectrum of biotic stresses. These genes, especially those in groups 1 and 2, hold promise as candidates for developing stress-tolerant banana cultivars. The study provides valuable insights into the molecular aspects of banana defense strategies, emphasizing the potential applications of PR-1 genes in enhancing banana resilience.

Two haplotype-resolved genome assemblies for AAB allotriploid bananas provide insights into banana subgenome asymmetric evolution and Fusarium wilt control

Bananas (Musa spp.) are one of the world's most important fruit crops and play a vital role in food security for many developing countries. Most banana cultivars are triploids derived from inter- and intraspecific hybridizations between the wild diploid ancestor species Musa acuminate (AA) and M. balbisiana (BB). We report two haplotype-resolved genome assemblies of the representative AAB-cultivated types, Plantain and Silk, and precisely characterize ancestral contributions by examining ancestry mosaics across the genome. Widespread asymmetric evolution is observed in their subgenomes, which can be linked to frequent homologous exchange events. We reveal the genetic makeup of triploid banana cultivars and verify that subgenome B is a rich source of disease resistance genes. Only 58.5% and 59.4% of Plantain and Silk genes, respectively, are present in all three haplotypes, with >50% of genes being differentially expressed alleles in different subgenomes. We observed that the number of upregulated genes in Plantain is significantly higher than that in Silk at one-week post-inoculation with Fusarium wilt tropical race 4 (Foc TR4), which confirms that Plantain can initiate defense responses faster than Silk. Additionally, we compared genomic and transcriptomic differences among the genes related to carotenoid synthesis and starch metabolism between Plantain and Silk. Our study provides resources for better understanding the genomic architecture of cultivated bananas and has important implications for Musa genetics and breeding.

Genome-Wide Identification and Expression Analysis of DWARF53 Gene in Response to GA and SL Related to Plant Height in Banana

Dwarfing is one of the common phenotypic variations in asexually reproduced progeny of banana, and dwarfed banana is not only windproof and anti-fallout but also effective in increasing acreage yield. As a key gene in the strigolactone signalling pathway, DWARF53 (D53) plays an important role in the regulation of the height of plants. In order to gain insight into the function of the banana D53 gene, this study conducted genome-wide identification of banana D53 gene based on the M. acuminata, M. balbisiana and M. itinerans genome database. Analysis of MaD53 gene expression under high temperature, low temperature and osmotic stress based on transcriptome data and RT-qPCR was used to analyse MaD53 gene expression in different tissues as well as in different concentrations of GA and SL treatments. In this study, we identified three MaD53, three MbD53 and two MiD53 genes in banana. Phylogenetic tree analysis showed that D53 Musa are equally related to D53 Asparagales and Poales. Both high and low-temperature stresses substantially reduced the expression of the MaD53 gene, but osmotic stress treatments had less effect on the expression of the MaD53 gene. GR24 treatment did not significantly promote the height of the banana, but the expression of the MaD53 gene was significantly reduced in roots and leaves. GA treatment at 100 mg/L significantly promoted the expression of the MaD53 gene in roots, but the expression of this gene was significantly reduced in leaves. In this study, we concluded that MaD53 responds to GA and SL treatments, but "Yinniaijiao" dwarf banana may not be sensitive to GA and SL.

High-quality genome assemblies for two Australimusa bananas (Musa spp.) and insights into regulatory mechanisms of superior fiber properties

Bananas (Musa spp.) are monocotyledonous plants with high genetic diversity in the Musaceae family that are cultivated mainly in tropical and subtropical countries. The fruits are a popular food, and the plants themselves have diverse uses. Four genetic groups (genomes) are thought to have contributed to current banana cultivars: Musa acuminata (A genome), Musa balbisiana (B genome), Musa schizocarpa (S genome), and species of the Australimusa section (T genome). However, the T genome has not been effectively explored. Here, we present the high-quality TT genomes of two representative accessions, Abaca (Musa textilis), with high-quality natural fiber, and Utafun (Musa troglodytarum, Fe'i group), with abundant β-carotene. Both the Abaca and Utafun assemblies comprise 10 pseudochromosomes, and their total genome sizes are 613 Mb and 619 Mb, respectively. Comparative genome analysis revealed that the larger size of the T genome is likely attributable to rapid expansion and slow removal of transposons. Compared with those of Musa AA or BB accessions or sisal (Agava sisalana), Abaca fibers exhibit superior mechanical properties, mainly because of their thicker cell walls with a higher content of cellulose, lignin, and hemicellulose. Expression of MusaCesA cellulose synthesis genes peaks earlier in Abaca than in AA or BB accessions during plant development, potentially leading to earlier cellulose accumulation during secondary cell wall formation. The Abaca-specific expressed gene MusaMYB26, which is directly regulated by MusaMYB61, may be an important regulator that promotes precocious expression of secondary cell wall MusaCesAs. Furthermore, MusaWRKY2 and MusaNAC68, which appear to be involved in regulating expression of MusaLAC and MusaCAD, may at least partially explain the high accumulation of lignin in Abaca. This work contributes to a better understanding of banana domestication and the diverse genetic resources in the Musaceae family, thus providing resources for Musa genetic improvement.

Genome-wide characterization of 2OGD superfamily for mining of susceptibility factors responding to various biotic stresses in Musa spp

Bananas are an important staple food and cash crop, but they are vulnerable to a variety of pests and diseases that substantially reduce yield and quality. Banana diseases are challenging to control and necessitate an integrated strategy, and development of resistant cultivars is one of the effective ways of managing diseases. Lasting disease resistance is the main goal in crop improvement and resistance mediated by a single resistant (R) gene mostly lack durability. However, long-term resistance can be obtained by inactivating susceptibility factors (S), which facilitate pathogen infection and proliferation. Identification and inactivation of susceptibility factors against the major pathogens like Fusarium oxysporum f. sp. cubense (Foc), Pseudocercospora eumusae and Pratylenchus coffeae in banana will be an effective way in developing banana varieties with more durable resistance. Downy mildew resistance 6 (DMR6) and DMR-like oxygenases (DLO1) are one such susceptibility factors and they belong to 2-oxoglutarate Fe(II) dependent oxygenases (2OGD) superfamily. 2OGDs are known to catalyze a plethora of reactions and also confer resistance to different pathogens in various crops, but not much is known about the 2OGD in Musa species. Through a comprehensive genome-wide analysis, 133 and 122 potential 2OGDs were systematically identified and categorized from the A and B genomes of banana, respectively. Real time expression of dmr6 and dlo1 genes showed positive correlation with transcriptome data upon Foc race1 and TR4 infection and examination of expression pattern of Macma4_04_g22670 (Ma04_g20880) and Macma4_02_g13590 (Ma02_g12040) genes revealed their involvement in Foc race1 and TR4 infections, respectively. Further the expression profile of 2OGDs, specifically Macma4_04_g25310 (Ma04_g23390), Macma4_08_g11980 (Ma08_g12090) and Macma4_04_g38910 (Ma04_g36640) shows that they may play a significant role as a susceptibility factor, particularly against P. eumusae and P. coffeae, implying that they can be exploited as a candidate gene for editing in developing resistant cultivars against these diseases. In summary, our findings contribute to a deeper comprehension of the evolutionary and functional aspects of 2OGDs in Musa spp. Furthermore, they highlight the substantial functions of these family constituents in the progression of diseases. These insights hold significance in the context of enhancing the genetic makeup of bananas to attain extended and more durable resistance against pathogens.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01380-y.

Genome-wide identification and expression analysis of the GRAS family under low-temperature stress in bananas

Introduction

GRAS, named after GAI, RGA, and SCR, is a class of plant-specific transcription factors family that plays a crucial role in growth and development, signal transduction, and various stress responses.

Methods

To understand the biological functions of the banana GRAS gene family, a genome-wide identification and bioinformatics analysis of the banana GRAS gene family was performed based on information from the M. acuminata, M. balbisiana, and M. itinerans genomic databases.

Result

In the present study, we identified 73 MaGRAS, 59 MbGRAS, and 58 MiGRAS genes in bananas at the whole-genome scale, and 56 homologous genes were identified in the three banana genomes. Banana GRASs can be classified into 10 subfamilies, and their gene structures revealed that most banana GRAS gDNAs lack introns. The promoter sequences of GRASs had a large number of cis-acting elements related to plant growth and development, phytohormone, and adversity stress responsiveness. The expression pattern of seven key members of MaGRAS response to low-temperature stress and different tissues was also examined by quantitative reverse transcription polymerase chain reaction (qRT-PCR). The microRNAs-MaGRASs target prediction showed perfect complementarity of seven GRAS genes with the five mac-miRNAs. The expression of all seven genes was lowest in roots, and the expression of five genes was highest in leaves during low-temperature stress. The expression of MaSCL27-2, MaSCL27-3, and MaSCL6-1 was significantly lower under low-temperature stress compared to the control, except for MaSCL27-2, which was slightly higher than the 28°C control at 4 h. The expression of MaSCL27-2, MaSCL27-3, and MaSCL6-1 dropped to the lowest levels at 24 h, 12 h, and 4 h, respectively. The MaSCL27-4 and MaSCL6-2 expression was intermittently upregulated, rising to the highest expression at 24h, while the expression of MaSCL22 was less variable, remaining at the control level with small changes.

Discussion

In summary, it is tentatively hypothesized that the GRAS family has an important function in low-temperature stress in bananas. This study provides a theoretical basis for further analyzing the function of the banana GRAS gene and the resistance of bananas to cold temperatures.

Characterization of Banana SNARE Genes and Their Expression Analysis under Temperature Stress and Mutualistic and Pathogenic Fungal Colonization

SNAREs (soluble N-ethylmaleimide-sensitive-factor attachment protein receptors) are engines for almost all of the membrane fusion and exocytosis events in organism cells. In this study, we identified 84 SNARE genes from banana (Musa acuminata). Gene expression analysis revealed that the expression of MaSNAREs varied a lot in different banana organs. By analyzing their expression patterns under low temperature (4 °C), high temperature (45 °C), mutualistic fungus (Serendipita indica, Si) and fungal pathogen (Fusarium oxysporum f. sp. Cubense Tropical Race 4, FocTR4) treatments, many MaSNAREs were found to be stress responsive. For example, MaBET1d was up-regulate by both low and high temperature stresses; MaNPSN11a was up-regulated by low temperature but down-regulated by high temperature; and FocTR4 treatment up-regulated the expression of MaSYP121 but down-regulated MaVAMP72a and MaSNAP33a. Notably, the upregulation or downregulation effects of FocTR4 on the expression of some MaSNAREs could be alleviated by priorly colonized Si, suggesting that they play roles in the Si-enhanced banana wilt resistance. Foc resistance assays were performed in tobacco leaves transiently overexpressing MaSYP121, MaVAMP72a and MaSNAP33a. Results showed that transient overexpression of MaSYP121 and MaSNPA33a suppressed the penetration and spread of both Foc1 (Foc Race 1) and FocTR4 in tobacco leaves, suggesting that they play positive roles in resisting Foc infection. However, the transient overexpression of MaVAMP72a facilitated Foc infection. Our study can provide a basis for understanding the roles of MaSNAREs in the banana responses to temperature stress and mutualistic and pathogenic fungal colonization.

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.

Genome-wide analysis of pathogenesis-related protein 1 (PR-1) gene family from Musa spp. and its role in defense response during stresses

Pathogenesis related protein-1 (PR-1) is the most abundantly produced protein during defense response against many biotic and abiotic stresses. However, knowledge on PR-1 gene family and its evolutionary relationship in banana is very limited. In order to study the potential role of PR-1 genes in banana, genome wide identification, structure analysis and expressions were performed. A total of 15 and 11 PR-1 genes were identified from A and B genomes of banana and the proteins encoded by this gene family are of varying lengths and harbor conserved domains and motifs. PR-1 genes are unevenly dispersed on 11 chromosomes with segmental duplication in both A and B genome, suggesting an important contribution of duplication in expansion of PR-1 gene family in banana. qRT-PCR analysis of PR-1 gene showed positive correlation with the RNAseq data under various stresses and examination of expression pattern of selected MaPR-1 genes in banana revealed its role in biotic and abiotic stresses in general and fusarium wilt in particular. This study provides significant insight into the functions of PR-1 genes which can be further exploited as a promising candidate for developing multiple stress tolerant banana varieties.

Genome-wide identification, characterization of expansin gene family of banana and their expression pattern under various stresses

Expansin, a cell wall-modifying gene family, has been well characterized and its role in biotic and abiotic stress resistance has been proven in many monocots, but not yet studied in banana, a unique model crop. Banana is one of the staple food crops in developing countries and its production is highly influenced by various biotic and abiotic factors. Characterizing the expansin genes of the ancestor genome (M. acuminata and M. balbisiana) of present day cultivated banana will enlighten their role in growth and development, and stress responses. In the present study, 58 (MaEXPs) and 55 (MbaEXPs) putative expansin genes were identified in A and B genome, respectively, and were grouped in four subfamilies based on phylogenetic analysis. Gene structure and its duplications revealed that EXPA genes are highly conserved and are under negative selection whereas the presence of more number of introns in other subfamilies revealed that they are diversifying. Expression profiling of expansin genes showed a distinct expression pattern for biotic and abiotic stress conditions. This study revealed that among the expansin subfamilies, EXPAs contributed significantly towards stress-resistant mechanism. The differential expression of MaEXPA18 and MaEXPA26 under drought stress conditions in the contrasting cultivar suggested their role in drought-tolerant mechanism. Most of the MaEXPA genes are differentially expressed in the root lesion nematode contrasting cultivars which speculated that this expansin subfamily might be the susceptible factor. The downregulation of MaEXPLA6 in resistant cultivar during Sigatoka leaf spot infection suggested that by suppressing this gene, resistance may be enhanced in susceptible cultivar. Further, in-depth studies of these genes will lead to gain insight into their role in various stress conditions in banana.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-021-03106-x.

Genome-wide identification of FAD gene family and their contributions to the temperature stresses and mutualistic and parasitic fungi colonization responses in banana

Fatty acid desaturase (FAD) plays important roles in plant growth and development and plant defense processes. In this study, we identified 27 MaFAD genes from the banana genome. According to the amino acid sequence similarities, their encoded proteins could be classified into five subfamilies. This classification is consistently supported by their gene and protein structures, conserved motifs and subcellular localizations. Segmental duplication events were found to play predominant roles in the MaFAD gene family expansion. Thirty miRNAs targeting MaFADs were identified and many hormone- and stress-responsive cis-acting elements and transcription factor binding sites (TFBSs) were identified in their promoters, indicating that the MaFADs expression regulation was very complicated. Gene expression analysis showed that some MaFADs showed significant differential expression in response to high and low temperature. FocTR4 influenced greatly the expression of several MaFADs and greatly induced the fatty acid (FA) accumulations in roots. Although S. indica showed no significant influence on the expression of most MaFADs, it could greatly alleviate the influence of FocTR4 on several MaFADs and FA biosynthesis. Our study revealed that MaFADs contributed greatly to the responses of high and low temperature stresses and mutualistic and parasitic fungi colonization in banana.

Genome-Wide Identification and Transcript Analysis of TCP Gene Family in Banana (Musa acuminata L.)

Plant-specific TEOSINTE-BRANCHED1/CYCLOIDEA/PROLIFERATING CELL FACTOR1 (TCP) gene family has versatile functions in diverse aspects of plants. However, less research on banana TCPs was done comprehensively. Accordingly, 48 banana TCP genes were characterized on aspects of gene structure, conserved motifs, phylogenetic relationship, and expression patterns. Members of the MaTCP gene family were unevenly distributed among 11 chromosomes and purification selection was the driving force of the MaTCP gene family. Gene duplication analysis indicated that segmental duplication is the major contributor to family expansion. Promoter analysis showed that MaTCPs might be involved in banana growth, development, and abiotic stress responses. Further, the expression of 12 MaTCPs was analyzed by real-time quantitative RT-PCR, and the protein interaction analysis showed that MaPCF10 and MaPCF13 may have an important function in banana fruit development and ripening. These results lay the foundation for further study of the functions of TCP genes in banana.

Genome-wide identification and expression pattern analysis of lipoxygenase gene family in banana

The LOX genes have been identified and characterized in many plant species, but studies on the banana LOX genes are very limited. In this study, we respectively identified 18 MaLOX, 11 MbLOX, and 12 MiLOX genes from the Musa acuminata, M. balbisiana and M. itinerans genome data, investigated their gene structures and characterized the physicochemical properties of their encoded proteins. Banana LOXs showed a preference for using and ending with G/C and their encoded proteins can be classified into 9-LOX, Type I 13-LOX and Type II 13-LOX subfamilies. The expansion of the MaLOXs might result from the combined actions of genome-wide, tandem, and segmental duplications. However, tandem and segmental duplications contribute to the expansion of MbLOXs. Transcriptome data based gene expression analysis showed that MaLOX1, 4, and 7 were highly expressed in fruit and their expression levels were significantly regulated by ethylene. And 11, 12 and 7 MaLOXs were found to be low temperature-, high temperature-, and Fusarium oxysporum f. sp. Cubense tropical race 4 (FocTR4)-responsive, respectively. MaLOX8, 9 and 13 are responsive to all the three stresses, MaLOX4 and MaLOX12 are high temperature- and FocTR4-responsive; MaLOX6 and MaLOX17 are significantly induced by low temperature and FocTR4; and the expression of MaLOX7 and MaLOX16 are only affected by high temperature. Quantitative real-time PCR (qRT-PCR) analysis revealed that the expression levels of several MaLOXs are regulated by MeJA and FocTR4, indicating that they can increase the resistance of banana by regulating the JA pathway. Additionally, the weighted gene co-expression network analysis (WGCNA) of MaLOXs revealed 3 models respectively for 5 (MaLOX7-11), 3 (MaLOX6, 13, and 17), and 1 (MaLOX12) MaLOX genes. Our findings can provide valuable information for the characterization, evolution, diversity and functionality of MaLOX, MbLOX and MiLOX genes and are helpful for understanding the roles of LOXs in banana growth and development and adaptations to different stresses.

Genome-scale analyses of polyketide synthases in banana: Phylogenetics and expression profiling forecast their candidacy in specialized metabolism

Plant type III polyketide synthases (PKSs) are associated with various functions in plant growth, development and defense by providing a multitude of polyketide scaffolds for diverse specialized metabolic pathways (SMPs). To decipher banana PKSs involved in specialized metabolism, genome-wide comparative analyses were conducted with A (Musa acuminata) and B (Musa balbisiana) genomes of banana. Both genomes retained eight chalcone synthases (CHSs), seven curcumin synthases (CURSs), three diketidyl-CoA synthases (DCSs) and one anther specific CHS (ASC). Segmental (42%) and tandem (37%) duplication events majorly flourished the banana PKS family. Six of 19 PKSs of A genome (designated as MaPKSs) showed relatively a higher expression in the root, corm, sheath, leaf and embryogenic cell suspension (ECS) of banana. To determine the defense response of MaPKSs and to highlight their candidacy in various SMPs, expression profiling was conducted by qPCR in ECSs treated with 100/200 μM of jasmonic acid (JA) and salicylic acid (SA) at 24/48 h. Maximum and subordinate expression induction of MaPKSs was apparent respectively against JA and SA treatments. Notably, most MaPKSs achieved their peak expression within 24 h of JA and the total flavonoid content was reached maximum within 24 h of JA/SA elicitations. Considering the homology, phylogeny, and expression levels in each analyzed sample (n = 13), three CHSs, three DCSs along with three CURSs and one ASC were selected as most promising candidates respectively for flavonoids, phenylphenalenones and sporopollenin biosynthesis in banana. Our findings provide a first-line resource to disclose the functions of banana PKSs involved in distinct SMPs.

Badnaviruses and banana genomes: a long association sheds light on Musa phylogeny and origin

Badnaviruses are double-stranded DNA pararetroviruses of the family Caulimoviridae. Badnaviral sequences found in banana are distributed over three main clades of the genus Badnavirus and exhibit wide genetic diversity. Interestingly, the nuclear genome of many plants, including banana, is invaded by numerous badnaviral sequences although badnaviruses do not require an integration step to replicate, unlike animal retroviruses. Here, we confirm that banana streak viruses (BSVs) are restricted to clades 1 and 3. We also show that only BSVs from clade 3 encompassing East African viral species are not integrated into Musa genomes, unlike BSVs from clade 1. Finally, we demonstrate that sequences from clade 2 are definitively integrated into Musa genomes with no evidence of episomal counterparts; all are phylogenetically distant from BSVs known to date. Using different molecular approaches, we dissected the coevolution between badnaviral sequences of clade 2 and banana by comparing badnavirus integration patterns across a banana sampling representing major Musa speciation events. Our data suggest that primary viral integrations occurred millions of years ago in banana genomes under different possible scenarios. Endogenous badnaviral sequences can be used as powerful markers to better characterize the Musa phylogeny, narrowing down the likely geographical origin of the Musa ancestor.

Structural and Functional Annotation of Transposable Elements Revealed a Potential Regulation of Genes Involved in Rubber Biosynthesis by TE-Derived siRNA Interference in Hevea brasiliensis

The natural rubber biosynthetic pathway is well described in Hevea, although the final stages of rubber elongation are still poorly understood. Small Rubber Particle Proteins and Rubber Elongation Factors (SRPPs and REFs) are proteins with major function in rubber particle formation and stabilization. Their corresponding genes are clustered on a scaffold1222 of the reference genomic sequence of the Hevea brasiliensis genome. Apart from gene expression by transcriptomic analyses, to date, no deep analyses have been carried out for the genomic environment of SRPPs and REFs loci. By integrative analyses on transposable element annotation, small RNAs production and gene expression, we analysed their role in the control of the transcription of rubber biosynthetic genes. The first in-depth annotation of TEs (Transposable Elements) and their capacity to produce TE-derived siRNAs (small interfering RNAs) is presented, only possible in the Hevea brasiliensis clone PB 260 for which all data are available. We observed that 11% of genes are located near TEs and their presence may interfere in their transcription at both genetic and epigenetic level. We hypothesized that the genomic environment of rubber biosynthesis genes has been shaped by TE and TE-derived siRNAs with possible transcriptional interference on their gene expression. We discussed possible functionalization of TEs as enhancers and as donors of alternative transcription start sites in promoter sequences, possibly through the modelling of genetic and epigenetic landscapes.

The structure, functional evolution, and evolutionary trajectories of the H+-PPase gene family in plants

BACKGROUND

The H+-PPase (pyrophosphatase) gene family is an important class of proton transporters that play key roles in plant development and stress resistance. Although the physiological and biochemical functions of H+-PPases are well characterized, the structural evolution and functional differentiation of this gene family remain unclear.

RESULTS

We identified 124 H+-PPase members from 27 plant species using complete genomic data obtained from algae to angiosperms. We found that all analyzed plants carried H+-PPase genes, and members were not limited to the two main types (type I and II). Differentiation of this gene family occurred early in evolutionary history, probably prior to the emergence of algae. The type I and II H+-PPase genes were retained during the subsequent evolution of higher plants, and their copy numbers increased rapidly in some angiosperms following whole-genome duplication (WGD) events, with obvious expression pattern differentiation among the new copies. We found significant functional divergence between type I and II H+-PPase genes, with both showing evidence for positive selection pressure. We classified angiosperm type I H+-PPases into subtypes Ia and non-Ia, which probably differentiated at an early stage of angiosperm evolution. Compared with non-Ia subtype, the Ia subtype appears to confer some advantage in angiosperms, as it is highly conserved and abundantly expressed, but shows no evidence for positive selection.

CONCLUSIONS

We hypothesized that there were many types of H+-PPase genes in the plant ancestral genome, and that different plant groups retained different types of these genes. In the early stages of angiosperm evolution, the type I H+-PPase genes differentiated into various subtypes. In addition, the expression pattern varied not only among genes of different types or subtypes, but also among copies of the same subtype. Based on the expression patterns and copy numbers of H+-PPase genes in higher plants, we propose two possible evolutionary trajectories for this gene family.

Comparative and evolutionary analyses of the divergence of plant oligosaccharyltransferase STT3 isoforms

STT3 is a catalytic subunit of hetero-oligomeric oligosaccharyltransferase (OST), which is important for asparagine-linked glycosylation. In mammals and plants, OSTs with different STT3 isoforms exhibit distinct levels of enzymatic efficiency or different responses to stressors. Although two different STT3 isoforms have been identified in both plants and animals, it remains unclear whether these isoforms result from gene duplication in an ancestral eukaryote. Furthermore, the molecular mechanisms underlying the functional divergences between the two STT3 isoforms in plant have not been well elucidated. Here, we conducted phylogenetic analysis of the major evolutionary node species and suggested that gene duplications of STT3 may have occurred independently in animals and plants. Across land plants, the exon-intron structure differed between the two STT3 isoforms, but was highly conserved for each isoform. Most angiosperm STT3a genes had 23 exons with intron phase 0, while STT3b genes had 6 exons with intron phase 2. Characteristic motifs (motif 18 and 19) of STT3s were mapped to different structure domains in the plant STT3 proteins. These two motifs overlap with regions of high nonsynonymous-to-synonymous substitution rates, suggesting the regions may be related to functional difference between STT3a and STT3b. In addition, promoter elements and gene expression profiles were different between the two isoforms, indicating expression pattern divergence of the two genes. Collectively, the identified differences may result in the functional divergence of plant STT3s.

Genome-Wide Identification and Characterization of the NAC Transcription Factor Family in Musa Acuminata and Expression Analysis during Fruit Ripening

Banana (Musa acuminata, AAA group) is a representative climacteric fruit with essential nutrients and pleasant flavors. Control of its ripening determines both the fruit quality and the shelf life. NAC (NAM, ATAF, CUC2) proteins, as one of the largest superfamilies of transcription factors, play crucial roles in various functions, especially developmental processes. Thus, it is important to conduct a comprehensive identification and characterization of the NAC transcription factor family at the genomic level in M. acuminata. In this article, a total of 181 banana NAC genes were identified. Phylogenetic analysis indicated that NAC genes in M. acuminata, Arabidopsis, and rice were clustered into 18 groups (S1–S18), and MCScanX analysis disclosed that the evolution of MaNAC genes was promoted by segmental duplication events. Expression patterns of NAC genes during banana fruit ripening induced by ethylene were investigated using RNA-Seq data, and 10 MaNAC genes were identified as related to fruit ripening. A subcellular localization assay of selected MaNACs revealed that they were all localized to the nucleus. These results lay a good foundation for the investigation of NAC genes in banana toward the biological functions and evolution.

Two large reciprocal translocations characterized in the disease resistance-rich burmannica genetic group of Musa acuminata

BACKGROUND AND AIMS

Banana cultivars are derived from hybridizations involving Musa acuminata subspecies. The latter diverged following geographical isolation in distinct South-east Asian continental regions and islands. Observation of chromosome pairing irregularities in meiosis of hybrids between these subspecies suggested the presence of large chromosomal structural variations. The aim of this study was to characterize such rearrangements.

METHODS

Marker (single nucleotide polymorphism) segregation in a self-progeny of the 'Calcutta 4' accession and mate-pair sequencing were used to search for chromosomal rearrangements in comparison with the M. acuminata ssp. malaccensis genome reference sequence. Signature segment junctions of the revealed chromosome structures were identified and searched in whole-genome sequencing data from 123 wild and cultivated Musa accessions.

KEY RESULTS

Two large reciprocal translocations were characterized in the seedy banana M. acuminata ssp. burmannicoides 'Calcutta 4' accession. One consisted of an exchange of a 240 kb distal region of chromosome 2 with a 7.2 Mb distal region of chromosome 8. The other involved an exchange of a 20.8 Mb distal region of chromosome 1 with a 11.6 Mb distal region of chromosome 9. Both translocations were found only in wild accessions belonging to the burmannicoides/burmannica/siamea subspecies. Only two of the 87 cultivars analysed displayed the 2/8 translocation, while none displayed the 1/9 translocation.

CONCLUSION

Two large reciprocal translocations were identified that probably originated in the burmannica genetic group. Accurate characterization of these translocations should enhance the use of this disease resistance-rich burmannica group in 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.

Recombination and Large Structural Variations Shape Interspecific Edible Bananas Genomes

Admixture and polyploidization are major recognized eukaryotic genome evolutionary processes. Their impacts on genome dynamics vary among systems and are still partially deciphered. Many banana cultivars are triploid (sometimes diploid) interspecific hybrids between Musa acuminata (A genome) and M. balbisiana (B genome). They have no or very low fertility, are vegetatively propagated and have been classified as “AB,” “AAB,” or “ABB” based on morphological characters. We used NGS sequence data to characterize the A versus B chromosome composition of nine diploid and triploid interspecific cultivars, to compare the chromosome structures of A and B genomes and analyze A/B chromosome segregations in a polyploid context. We showed that interspecific recombination occurred frequently between A and B chromosomes. We identified two large structural variations between A and B genomes, a reciprocal translocation and an inversion that locally affected recombination and led to segregation distortion and aneuploidy in a triploid progeny. Interspecific recombination and large structural variations explained the mosaic genomes observed in edible bananas. The unprecedented resolution in deciphering their genome structure allowed us to start revisiting the origins of banana cultivars and provided new information to gain insight into the impact of interspecificity on genome evolution. It will also facilitate much more effective assessment of breeding strategies.

PISTILLATA paralogs in Tarenaya hassleriana have diverged in interaction specificity

BACKGROUND

Floral organs are specified by MADS-domain transcription factors that act in a combinatorial manner, as summarized in the (A)BCE model. However, this evolutionarily conserved model is in contrast to a remarkable amount of morphological diversity in flowers. One of the mechanisms suggested to contribute to this diversity is duplication of floral MADS-domain transcription factors. Although gene duplication is often followed by loss of one of the copies, sometimes both copies are retained. If both copies are retained they will initially be redundant, providing freedom for one of the paralogs to change function. Here, we examine the evolutionary fate and functional consequences of a transposition event at the base of the Brassicales that resulted in the duplication of the floral regulator PISTILLATA (PI), using Tarenaya hassleriana (Cleomaceae) as a model system.

RESULTS

The transposition of a genomic region containing a PI gene led to two paralogs which are located at different positions in the genome. The original PI copy is syntenic in position with most angiosperms, whereas the transposed copy is syntenic with the PI genes in Brassicaceae. The two PI paralogs of T. hassleriana have very similar expression patterns. However, they may have diverged in function, as only one of these PI proteins was able to act heterologously in the first whorl of A. thaliana flowers. We also observed differences in protein complex formation between the two paralogs, and the two paralogs exhibit subtle differences in DNA-binding specificity. Sequence analysis indicates that most of the protein sequence divergence between the two T. hassleriana paralogs emerged in a common ancestor of the Cleomaceae and the Brassicaceae.

CONCLUSIONS

We found that the PI paralogs in T. hassleriana have similar expression patterns, but may have diverged at the level of protein function. Data suggest that most protein sequence divergence occurred rapidly, prior to the origin of the Brassicaceae and Cleomaceae. It is tempting to speculate that the interaction specificities of the Brassicaceae-specific PI proteins are different compared to the PI found in other angiosperms. This could lead to PI regulating partly different genes in the Brassicaceae, and ultimately might result in change floral in morphology.

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

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

Assessment of genetic variation among wild Alpinia nigra (Zingiberaceae) population: an approach based on molecular phylogeny

Genetic structure was evaluated among wild Alpinia nigra (Gaertn.) B.L. Burtt, populations. The information of genetic relatedness was developed using random amplified polymorphic DNA (RAPD), inter-simple sequence repeat (ISSR) and barcoding loci (plastid and mitochondrial). The order (high to low) of Shannon's information index (I) and Nei's gene diversity (h) from the populations was: "IIT Guwahati" > "Amingaon" > "Saraighat". Genetic diversity decreased and genetic differentiation increased among the three populations. We observed no isolation by distance thus lower amount of gene flow was observed. Narrow range of genetic distance among the three populations and appearance of two distinct clusters strengthened the geographical isolation in dendrogram and principal component analysis. No mutation among the three populations was observed for seven plastid loci and two mitochondrial tested suggesting the taxonomic homogeneity. The phylogeny based on nine barcoding loci supported our observation that individuals of IIT Guwahati were partially isolated from the outside populations. Our study will provide a backbone for developing strategies to resist habitat fragmentation of Zingiberaceous plants.

Genome-wide identification and characterization of the CKII gene family in the cultivated banana cultivar (Musa spp. cv Tianbaojiao) and the wild banana (Musa itinerans)

Plant casein kinase II (CKII) plays an essential role in regulating plant growth and development, and responses to biotic and abiotic stresses. Here, we report the identification and characterization of the CKII family genes in Musa spp. cv. ‘Tianbaojiao’ (AAA group) and the wild banana (Musa itinerans). The 13 cDNA sequences of the CKII family members were identified both in ‘Tianbaojiao’ and wild banana, respectively. The differences between CKII α and CKII β members are corroborated through the subcellular localizations, phosphorylation sites and gene structures. The cloning of CKII β-like-2 gDNA sequences in wild banana and ‘Tianbaojiao’ and the analysis of gene structures showed MiCKIIβ-like-2b and MaCKIIβ-like-2 are likely alternatively spliced transcripts, which were derived from the alternative splicing events that involved exon deletion. The qPCR validation showed differential expression CKII family members in response to cold stress and also in all tested tissues (leaf, pseudostem and root) of wild banana. In particular, the normal transcript MiCKIIβ-like-2a was highly expressed in response to cold stress in wild banana; oppositely, the alternatively spliced transcript MiCKIIβ-like-2b was quite lowly expressed. The complex origin and long-term evolution of Musa lineage might explain the alternative splicing events of CKII β-like-2.

Expression of the SIN3 homologue from banana, MaSIN3, suppresses ABA responses globally during plant growth in Arabidopsis

The SIN3 family of co-repressors is a family of highly conserved eukaryotic repressor proteins that regulates diverse functions in yeasts and animals but remains largely uncharacterized functionally even in plants like Arabidopsis. The sole SIN3 homologue in banana, MaSIN3, was identified as a 1408 amino acids, nuclear localized protein conserved to other SIN3s in the PAH, HID and HCR domains. Interestingly, MaSIN3 over-expression in Arabidopsis mimics a state of reduced ABA responses throughout plant development affecting growth processes such as germination, root growth, stomatal closure and water loss, flowering and senescence. The reduction in ABA responses is not due to reduced ABA levels but due to suppression of expression of several transcription factors mediating ABA responses. Transcript levels of negative regulators of germination (ABI3, ABI5, PIL5, RGL2 and RGL3) are reduced post-imbibition while those responsible for GA biosynthesis are up-regulated in transgenic MaSIN3 over-expressers. ABA-associated transcription factors are also down-regulated in response to ABA treatment. The HDAC inhibitors, SAHA and sodium butyrate, in combination with ABA differentially suppress germination in control and transgenic lines suggesting the recruitment by MaSIN3 of HDACs involved in suppression of ABA responses in different processes. The studies provide an insight into the ability of MaSIN3 to specifically affect a subset of developmental processes governed largely by ABA.

Dominant Allele Phylogeny and Constitutive Subgenome Haplotype Inference in Bananas Using Mitochondrial and Nuclear Markers

Cultivated bananas (Musa spp.) have undergone domestication patterns involving crosses of wild progenitors followed by long periods of clonal propagation. Majority of cultivated bananas are polyploids with different constitutive subgenomes and knowledge on phylogenies to their progenitors at the species and subspecies levels is essential. Here, the mitochondrial (NAD1) and nuclear (CENH3) markers were used to phylogenetically position cultivated banana genotypes to diploid progenitors. The CENH3 nuclear marker was used to identify a minimum representative haplotype number in polyploids and diploid bananas based on single nucleotide polymorphisms. The mitochondrial marker NAD1 was observed to be ideal in differentiating bananas of different genomic constitutions based on size of amplicons as well as sequence. The genotypes phylogenetically segregated based on the dominant genome; AAB genotypes grouped with AA and AAA, and the ABB together with BB. Both markers differentiated banana sections, but could not differentiate subspecies within the A genomic group. On the basis of CENH3 marker, a total of 13 haplotypes (five in both diploid and triploid, three in diploids, and rest unique to triploids) were identified from the genotypes tested. The presence of haplotypes, which were common in diploids and triploids, stipulate possibility of a shared ancestry in the genotypes involved in this study. Furthermore, the presence of multiple haplotypes in some diploid bananas indicates their being heterozygous. The haplotypes identified in this study are of importance because they can be used to check the level of homozygozity in breeding lines as well as to track segregation in progenies.

Genome-wide analysis of the AP2/ERF family in Musa species reveals divergence and neofunctionalisation during evolution

AP2/ERF domain containing transcription factor super family is one of the important regulators in the plant kingdom. The involvement of AP2/ERF family members has been elucidated in various processes associated with plant growth, development as well as in response to hormones, biotic and abiotic stresses. In this study, we carried out genome-wide analysis to identify members of AP2/ERF family in Musa acuminata (A genome) and Musa balbisiana (B genome) and changes leading to neofunctionalisation of genes. Analysis identified 265 and 318 AP2/ERF encoding genes in M. acuminata and M. balbisiana respectively which were further classified into ERF, DREB, AP2, RAV and Soloist groups. Comparative analysis indicated that AP2/ERF family has undergone duplication, loss and divergence during evolution and speciation of the Musa A and B genomes. We identified nine genes which are up-regulated during fruit ripening and might be components of the regulatory machinery operating during ethylene-dependent ripening in banana. Tissue-specific expression analysis of the genes suggests that different regulatory mechanisms might be involved in peel and pulp ripening process through recruiting specific ERFs in these tissues. Analysis also suggests that MaRAV-6 and MaERF026 have structurally diverged from their M. balbisiana counterparts and have attained new functions during ripening.

Gene Evolutionary Trajectories and GC Patterns Driven by Recombination in Zea mays

Recombination occurring during meiosis is critical for creating genetic variation and plays an essential role in plant evolution. In addition to creating novel gene combinations, recombination can affect genome structure through altering GC patterns. In maize (Zea mays) and other grasses, another intriguing GC pattern exists. Maize genes show a bimodal GC content distribution that has been attributed to nucleotide bias in the third, or wobble, position of the codon. Recombination may be an underlying driving force given that recombination sites are often associated with high GC content. Here we explore the relationship between recombination and genomic GC patterns by comparing GC gene content at each of the three codon positions (GC₁, GC₂, and GC₃, collectively termed GC_x) to instances of a variable GC-rich motif that underlies double strand break (DSB) hotspots and to meiocyte-specific gene expression. Surprisingly, GC_x bimodality in maize cannot be fully explained by the codon wobble hypothesis. High GC_x genes show a strong overlap with the DSB hotspot motif, possibly providing a mechanism for the high evolutionary rates seen in these genes. On the other hand, genes that are turned on in meiosis (early prophase I) are biased against both high GC_x genes and genes with the DSB hotspot motif, possibly allowing important meiotic genes to avoid DSBs. Our data suggests a strong link between the GC-rich motif underlying DSB hotspots and high GC_x genes.

Large distribution and high sequence identity of a Copia-type retrotransposon in angiosperm families

Retrotransposons are the main component of plant genomes. Recent studies have revealed the complexity of their evolutionary dynamics. Here, we have identified Copia25 in Coffea canephora, a new plant retrotransposon belonging to the Ty1-Copia superfamily. In the Coffea genomes analyzed, Copia25 is present in relatively low copy numbers and transcribed. Similarity sequence searches and PCR analyses show that this retrotransposon with LTRs (Long Terminal Repeats) is widely distributed among the Rubiaceae family and that it is also present in other distantly related species belonging to Asterids, Rosids and monocots. A particular situation is the high sequence identity found between the Copia25 sequences of Musa, a monocot, and Ixora, a dicot species (Rubiaceae). Our results reveal the complexity of the evolutionary dynamics of the ancient element Copia25 in angiosperm, involving several processes including sequence conservation, rapid turnover, stochastic losses and horizontal transfer.

Small RNA Profiling of Two Important Cultivars of Banana and Overexpression of miRNA156 in Transgenic Banana Plants

Micro RNAs (miRNAs) are a class of non-coding, short RNAs having important roles in regulation of gene expression. Although plant miRNAs have been studied in detail in some model plants, less is known about these miRNAs in important fruit plants like banana. miRNAs have pivotal roles in plant growth and development, and in responses to diverse biotic and abiotic stress stimuli. Here, we have analyzed the small RNA expression profiles of two different economically significant banana cultivars by using high-throughput sequencing technology. We identified a total of 170 and 244 miRNAs in the two libraries respectively derived from cv. Grand Naine and cv. Rasthali leaves. In addition, several cultivar specific microRNAs along with their putative target transcripts were also detected in our studies. To validate our findings regarding the small RNA profiles, we also undertook overexpression of a common microRNA, MusamiRNA156 in transgenic banana plants. The transgenic plants overexpressing the stem-loop sequence derived from MusamiRNA156 gene were stunted in their growth together with peculiar changes in leaf anatomy. These results provide a foundation for further investigations into important physiological and metabolic pathways operational in banana in general and cultivar specific traits in particular.

Plant terpenes: defense responses, phylogenetic analysis, regulation and clinical applications

The terpenoids constitute the largest class of natural products and many interesting products are extensively applied in the industrial sector as flavors, fragrances, spices and are also used in perfumery and cosmetics. Many terpenoids have biological activities and also used for medical purposes. In higher plants, the conventional acetate-mevalonic acid pathway operates mainly in the cytosol and mitochondria and synthesizes sterols, sesquiterpenes and ubiquinones mainly. In the plastid, the non-mevalonic acid pathway takes place and synthesizes hemi-, mono-, sesqui-, and diterpenes along with carotenoids and phytol tail of chlorophyll. In this review paper, recent developments in the biosynthesis of terpenoids, indepth description of terpene synthases and their phylogenetic analysis, regulation of terpene biosynthesis as well as updates of terpenes which have entered in the clinical studies are reviewed thoroughly.

Doubling down on genomes: polyploidy and crop plants

Polyploidy, or whole genome multiplication, is ubiquitous among angiosperms. Many crop species are relatively recent allopolyploids, resulting from interspecific hybridization and polyploidy. Thus, an appreciation of the evolutionary consequences of (allo)polyploidy is central to our understanding of crop plant domestication, agricultural improvement, and the evolution of angiosperms in general. Indeed, many recent insights into plant biology have been gleaned from polyploid crops, including, but not limited to wheat, tobacco, sugarcane, apple, and cotton. A multitude of evolutionary processes affect polyploid genomes, including rapid and substantial genome reorganization, transgressive gene expression alterations, gene fractionation, gene conversion, genome downsizing, and sub- and neofunctionalization of duplicate genes. Often these genomic changes are accompanied by heterosis, robustness, and the improvement of crop yield, relative to closely related diploids. Historically, however, the genome-wide analysis of polyploid crops has lagged behind those of diploid crops and other model organisms. This lag is partly due to the difficulties in genome assembly, resulting from the genomic complexities induced by combining two or more evolutionarily diverged genomes into a single nucleus and by the significant size of polyploid genomes. In this review, we explore the role of polyploidy in angiosperm evolution, the domestication process and crop improvement. We focus on the potential of modern technologies, particularly next-generation sequencing, to inform us on the patterns and processes governing polyploid crop improvement and phenotypic change subsequent to domestication.

The common evolutionary history of badnaviruses and banana

Recent plant genome sequencing efforts have revealed myriad viral sequences suggesting a cryptic interaction between both partners. Interestingly, no integration step has ever been reported as an obligatory step in the life cycle of plant viruses. Circular dsDNA viruses belonging to the family Caulimoviridae are the most abundant among integrated plant viral sequences. In this review, we describe how this hitherto hidden interaction could inform the evolutionary history of both partners badnaviruses and banana plants.

"A draft Musa balbisiana genome sequence for molecular genetics in polyploid, inter- and intra-specific Musa hybrids"

BACKGROUND

Modern banana cultivars are primarily interspecific triploid hybrids of two species, Musa acuminata and Musa balbisiana, which respectively contribute the A- and B-genomes. The M. balbisiana genome has been associated with improved vigour and tolerance to biotic and abiotic stresses and is thus a target for Musa breeding programs. However, while a reference M. acuminata genome has recently been released (Nature 488:213-217, 2012), little sequence data is available for the corresponding B-genome.To address these problems we carried out Next Generation gDNA sequencing of the wild diploid M. balbisiana variety 'Pisang Klutuk Wulung' (PKW). Our strategy was to align PKW gDNA reads against the published A-genome and to extract the mapped consensus sequences for subsequent rounds of evaluation and gene annotation.

RESULTS

The resulting B-genome is 79% the size of the A-genome, and contains 36,638 predicted functional gene sequences which is nearly identical to the 36,542 of the A-genome. There is substantial sequence divergence from the A-genome at a frequency of 1 homozygous SNP per 23.1 bp, and a high degree of heterozygosity corresponding to one heterozygous SNP per 55.9 bp. Using expressed small RNA data, a similar number of microRNA sequences were predicted in both A- and B-genomes, but additional novel miRNAs were detected, including some that are unique to each genome. The usefulness of this B-genome sequence was evaluated by mapping RNA-seq data from a set of triploid AAA and AAB hybrids simultaneously to both genomes. Results for the plantains demonstrated the expected 2:1 distribution of reads across the A- and B-genomes, but for the AAA genomes, results show they contain regions of significant homology to the B-genome supporting proposals that there has been a history of interspecific recombination between homeologous A and B chromosomes in Musa hybrids.

CONCLUSIONS

We have generated and annotated a draft reference Musa B-genome and demonstrate that this can be used for molecular genetic mapping of gene transcripts and small RNA expression data from several allopolyploid banana cultivars. This draft therefore represents a valuable resource to support the study of metabolism in inter- and intraspecific triploid Musa hybrids and to help direct breeding programs.

Ancient cis-regulatory constraints and the evolution of genome architecture

The order of genes along metazoan chromosomes has generally been thought to be largely random, with few implications for organismal function. However, two recent studies, reporting hundreds of pairs of genes that have remained linked in diverse metazoan species over hundreds of millions of years of evolution, suggest widespread functional implications for gene order. These associations appear to largely reflect cis-regulatory constraints, with either (i) multiple genes sharing transcriptional regulatory elements, or (ii) regulatory elements for a developmental gene being found within a neighboring 'bystander' gene (known as a genomic regulatory block). We discuss implications, questions raised, and new research directions arising from these studies, as well as evidence for similar phenomena in other eukaryotic groups.

Three infectious viral species lying in wait in the banana genome

Plant pararetroviruses integrate serendipitously into their host genomes. The banana genome harbors integrated copies of banana streak virus (BSV) named endogenous BSV (eBSV) that are able to release infectious pararetrovirus. In this investigation, we characterized integrants of three BSV species-Goldfinger (eBSGFV), Imove (eBSImV), and Obino l'Ewai (eBSOLV)-in the seedy Musa balbisiana Pisang klutuk wulung (PKW) by studying their molecular structure, genomic organization, genomic landscape, and infectious capacity. All eBSVs exhibit extensive viral genome duplications and rearrangements. eBSV segregation analysis on an F1 population of PKW combined with fluorescent in situ hybridization analysis showed that eBSImV, eBSOLV, and eBSGFV are each present at a single locus. eBSOLV and eBSGFV contain two distinct alleles, whereas eBSImV has two structurally identical alleles. Genotyping of both eBSV and viral particles expressed in the progeny demonstrated that only one allele for each species is infectious. The infectious allele of eBSImV could not be identified since the two alleles are identical. Finally, we demonstrate that eBSGFV and eBSOLV are located on chromosome 1 and eBSImV is located on chromosome 2 of the reference Musa genome published recently. The structure and evolution of eBSVs suggest sequential integration into the plant genome, and haplotype divergence analysis confirms that the three loci display differential evolution. Based on our data, we propose a model for BSV integration and eBSV evolution in the Musa balbisiana genome. The mutual benefits of this unique host-pathogen association are also discussed.

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.

Cytogenetic evidence of mixed disomic and polysomic inheritance in an allotetraploid (AABB) Musa genotype

BACKGROUND AND AIMS

Edible bananas originated mainly from two wild species, Musa acuminata Colla (AA) and Musa balbisiana Colla (BB), and triploid cultivars with an AAA, AAB or ABB genome are the most widely used. In the present study, chromosome pairing affinities are investigated in a sterile AB Indian variety and in its fertile colchicine-induced allotetraploid (AABB) derivative to determine the inheritance pattern of the tetraploid genotype. The potential implications of interspecific recombination and chromosomal composition of diploid gametes for Musa improvement are presented.

METHODS

The pairing of different chromosome sets at diploid and tetraploid levels was investigated through a combination of conventional cytogenetic and genomic in-situ hybridization (GISH) analyses of meiotic chromosomes, leading to a likelihood model of the pairing behaviour. GISH analysis of mitotic chromosomes was also conducted to reveal the chromosome constitution of hybrids derived from crosses involving the allotetraploid genotype.

KEY RESULTS

Analysis of chromosome associations at both ploidy levels suggested that the newly formed allotetraploid behaves as a 'segmental allotetraploid' with three chromosome sets in a tetrasomic pattern, three sets in a likely disomic pattern and the five remaining sets in an intermediate pattern. Balanced and unbalanced diploid gametes were detected in progenies, with the chromosome constitution appearing to be more homogenous in pollen than in ovules.

CONCLUSIONS

Colchicine-induced allotetraploids in Musa provide access to the genetic background of natural AB varieties. The segmental inheritance pattern exhibited by the AABB allotetraploid genotype implies chromosome exchanges between M. acuminata and M. balbisiana species and opens new horizons for reciprocal transfer of valuable alleles.

Mining of EST-SSR markers of Musa and their transferability studies among the members of order the Zingiberales

Expressed sequence tags (ESTs) databases of 11 Musa complementary DNA libraries were retrieved from National Center of Biotechnology Information and used for mining simple sequence repeats (SSRs). Out of 21,056 unique ESTs, SSR regions were found only in 5,158 ESTs. Among these SSR containing ESTs, the occurrence of trinucleotide repeats are the most abundant followed by mono-, di-, tetra-, hexa-, and pentanucleotides. Moreover, this study showed that the rate of class II SSRs (<20 nucleotides) was higher than the class I SSRs (<20 nucleotides), and proportion of class I and II SSRs as abundant for tri-repeats. As a representative sample, primers were synthesized for 24 ESTs, carrying >12 nucleotides of SSR region, and tested among the various genomic group of Musa accessions. The result showed that 88 % of primers were functional primers, and 43 % are showing polymorphism among the Musa accessions. Transferability studies of Musa EST-SSRs among the genera of the order Zingiberales exhibited 100 and 58 % transferability in Musaceae and Zingiberaceae, respectively. The sequence comparison of SSR regions among the different Musa accessions confirmed that polymorphism is mainly due to the variation in repeat length. High percentage of cross-species, cross-genera, and cross-family transferability also suggested that these Musa EST-SSR markers will be a valuable resource for the comparative mapping by developing COS markers, in evolutionary studies and in improvement of the members of Zingiberaceae and Musaceae.

Generation and analysis of expressed sequence tags in the extreme large genomes Lilium and Tulipa

BACKGROUND

Bulbous flowers such as lily and tulip (Liliaceae family) are monocot perennial herbs that are economically very important ornamental plants worldwide. However, there are hardly any genetic studies performed and genomic resources are lacking. To build genomic resources and develop tools to speed up the breeding in both crops, next generation sequencing was implemented. We sequenced and assembled transcriptomes of four lily and five tulip genotypes using 454 pyro-sequencing technology.

RESULTS

Successfully, we developed the first set of 81,791 contigs with an average length of 514 bp for tulip, and enriched the very limited number of 3,329 available ESTs (Expressed Sequence Tags) for lily with 52,172 contigs with an average length of 555 bp. The contigs together with singletons covered on average 37% of lily and 39% of tulip estimated transcriptome. Mining lily and tulip sequence data for SSRs (Simple Sequence Repeats) showed that di-nucleotide repeats were twice more abundant in UTRs (UnTranslated Regions) compared to coding regions, while tri-nucleotide repeats were equally spread over coding and UTR regions. Two sets of single nucleotide polymorphism (SNP) markers suitable for high throughput genotyping were developed. In the first set, no SNPs flanking the target SNP (50 bp on either side) were allowed. In the second set, one SNP in the flanking regions was allowed, which resulted in a 2 to 3 fold increase in SNP marker numbers compared with the first set. Orthologous groups between the two flower bulbs: lily and tulip (12,017 groups) and among the three monocot species: lily, tulip, and rice (6,900 groups) were determined using OrthoMCL. Orthologous groups were screened for common SNP markers and EST-SSRs to study synteny between lily and tulip, which resulted in 113 common SNP markers and 292 common EST-SSR. Lily and tulip contigs generated were annotated and described according to Gene Ontology terminology.

CONCLUSIONS

Two transcriptome sets were built that are valuable resources for marker development, comparative genomic studies and candidate gene approaches. Next generation sequencing of leaf transcriptome is very effective; however, deeper sequencing and using more tissues and stages is advisable for extended comparative studies.

The banana (Musa acuminata) genome and the evolution of monocotyledonous plants

Bananas (Musa spp.), including dessert and cooking types, are giant perennial monocotyledonous herbs of the order Zingiberales, a sister group to the well-studied Poales, which include cereals. Bananas are vital for food security in many tropical and subtropical countries and the most popular fruit in industrialized countries. The Musa domestication process started some 7,000 years ago in Southeast Asia. It involved hybridizations between diverse species and subspecies, fostered by human migrations, and selection of diploid and triploid seedless, parthenocarpic hybrids thereafter widely dispersed by vegetative propagation. Half of the current production relies on somaclones derived from a single triploid genotype (Cavendish). Pests and diseases have gradually become adapted, representing an imminent danger for global banana production. Here we describe the draft sequence of the 523-megabase genome of a Musa acuminata doubled-haploid genotype, providing a crucial stepping-stone for genetic improvement of banana. We detected three rounds of whole-genome duplications in the Musa lineage, independently of those previously described in the Poales lineage and the one we detected in the Arecales lineage. This first monocotyledon high-continuity whole-genome sequence reported outside Poales represents an essential bridge for comparative genome analysis in plants. As such, it clarifies commelinid-monocotyledon phylogenetic relationships, reveals Poaceae-specific features and has led to the discovery of conserved non-coding sequences predating monocotyledon-eudicotyledon divergence.

Evaluation of different methods of protein extraction and identification of differentially expressed proteins upon ethylene-induced early-ripening in banana peels

BACKGROUND

Banana peels (Musa spp.) are a good example of a plant tissue where protein extraction is challenging due to the abundance of interfering metabolites. Sample preparation is a critical step in proteomic research and is critical for good results.

RESULTS

We sought to evaluate three methods of protein extraction: trichloroacetic acid (TCA)-acetone precipitation, phenol extraction, and TCA precipitation. We found that a modified phenol extraction protocol was the most optimal method. SDS-PAGE and two-dimensional gel electrophoresis (2-DE) demonstrated good protein separation and distinct spots of high quality protein. Approximately 300 and 550 protein spots were detected on 2-DE gels at pH values of 3-10 and 4-7, respectively. Several spots were excised from the 2-DE gels and identified by mass spectrometry.

CONCLUSIONS

The protein spots identified were found to be involved in glycolysis, the tricarboxylic acid cycle, and the biosynthesis of ethylene. Several of the identified proteins may play important roles in banana ripening.

Foundation characteristics of edible Musa triploids revealed from allelic distribution of SSR markers

BACKGROUND AND AIMS

The production of triploid banana and plantain (Musa spp.) cultivars with improved characteristics (e.g. greater disease resistance or higher yield), while still preserving the main features of current popular cultivars (e.g. taste and cooking quality), remains a major challenge for Musa breeders. In this regard, breeders require a sound knowledge of the lineage of the current sterile triploid cultivars, to select diploid parents that are able to transmit desirable traits, together with a breeding strategy ensuring final triploidization and sterility. Highly polymorphic single sequence repeats (SSRs) are valuable markers for investigating phylogenetic relationships.

METHODS

Here, the allelic distribution of each of 22 SSR loci across 561 Musa accessions is analysed.

KEY RESULTS AND CONCLUSIONS

We determine the closest diploid progenitors of the triploid 'Cavendish' and 'Gros Michel' subgroups, valuable information for breeding programmes. Nevertheless, in establishing the likely monoclonal origin of the main edible triploid banana subgroups (i.e. 'Cavendish', 'Plantain' and 'Mutika-Lujugira'), we postulated that the huge phenotypic diversity observed within these subgroups did not result from gamete recombination, but rather from epigenetic regulations. This emphasizes the need to investigate the regulatory mechanisms of genome expression on a unique model in the plant kingdom. We also propose experimental standards to compare additional and independent genotyping data for reference.

Divergent evolution of oxidosqualene cyclases in plants

Triterpenes are one of the largest classes of plant metabolites and have important functions. A diverse array of triterpenoid skeletons are synthesized via the isoprenoid pathway by enzymatic cyclization of 2,3-oxidosqualene. The genomes of the lower plants Chlamydomonas reinhardtii and moss (Physcomitrella patens) contain just one oxidosqualene cyclase (OSC) gene (for sterol biosynthesis), whereas the genomes of higher plants contain nine to 16 OSC genes. Here we carry out functional analysis of rice OSCs and rigorous phylogenetic analysis of 96 OSCs from higher plants, including Arabidopsis thaliana, Oryza sativa, Sorghum bicolor and Brachypodium distachyon. The functional analysis identified an amino acid sequence for isoarborinol synthase (OsIAS) (encoded by Os11g35710/OsOSC11) in rice. Our phylogenetic analysis suggests that expansion of OSC members in higher plants has occurred mainly through tandem duplication followed by positive selection and diversifying evolution, and consolidated the previous suggestion that dicot triterpene synthases have been derived from an ancestral lanosterol synthase instead of directly from their cycloartenol synthases. The phylogenetic trees are consistent with the reaction mechanisms of the protosteryl and dammarenyl cations which parent a wide variety of triterpene skeletal types, allowing us to predict the functions of the uncharacterized OSCs.

Individual chromosome assignment and chromosomal collinearity in Gossypium thurberi, G. trilobum and D subgenome of G. barbadense revealed by BAC-FISH

The experiment on individual chromosome assignments and chromosomal diversity was conducted using a multi-probe fluorescence in situ hybridization (FISH) system in D subgenome of tetraploid Gossypium barbadense (D(b)), G. thurberi (D(1)) and G. trilobum (D(8)), which the later two were the possible subgenome donors of tetraploid cottons. The FISH probes contained a set of bacterial artificial chromosome (BAC) clones specific to 13 individual chromosomes from D subgenome of G. hirsutum (D(h)), a D genome centromere-specific BAC clone 150D24, 45S and 5S ribosomal DNA (rDNA) clones, respectively. All tested chromosome orientations were confirmed by the centromere-specific BAC probe. In D(1) and D(8), four 45S rDNA loci were found assigning at the end of the short arm of chromosomes 03, 07, 09 and 11, while one 5S rDNA locus was successfully marked at pericentromeric region of the short arm of chromosome 09. In D(b), three 45S rDNA loci and two 5S rDNA loci were found out. Among them, two 45S rDNA loci were located at the terminal of the short arm of chromosomes D(b)07 and D(b)09, whilst one 5S rDNA locus was found situating near centromeric region of the short arm of chromosome D(b)09. The positions of the BAC clones specific to the 13 individual chromosomes from D(h) were compared between D(1), D(8) and D(b). The result showed the existence of chromosomal collinearity within D(1) and D(8), and as well between them and D(b). The results will serve as a base for understanding chromosome structure of cotton and polyploidy evolution of cotton genome and will provide bio-information for assembling the sequences of finished and the on-going cotton whole genome sequencing projects.