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

Magali Lescot, Pietro Piffanelli, Ana Ciampi, Manuel Ruiz, Guillaume Blanc, Jim Leebens-Mack, Felipe da Silva, Candice Mr Santos, Angélique d'Hont, Olivier Garsmeur, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

Magali Lescot, Pietro Piffanelli, Ana Ciampi, Manuel Ruiz, Guillaume Blanc, et al. Insights into the Musa genome: Syntenic relationships to rice and between Musa species. BMC Genomics, 2008, 9 (1), 10.1186/1471-2164-9-58. doi:10.1186/1471-2164-9-58

HAL Id: hal-02978134

https://hal.science/hal-02978134

Submitted on 28 Apr 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>
Research article

Open Access

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

Magali Lescot*†¹,⁹, Pietro Piffanelli†¹,¹⁰, Ana Y Ciampi³, Manuel Ruiz¹, Guillaume Blanc⁵, Jim Leebens-Mack⁶, Felipe R da Silva³, Candice MR Santos³, Angélique D'Hont¹, Olivier Garsmeur¹, Alberto D Vilarinhos¹,⁷, Hiroyuki Kanamori⁸, Takashi Matsumoto⁸, Catherine M Ronning², Foo Cheung², Brian J Haas², Ryan Althoff², Tammy Arbogast², Erin Hine², Georgios J Pappas Jr⁴, Takuji Sasaki⁸, Manoel T Souza Jr³, Robert NG Miller⁴, Jean-Christophe Glaszmann¹ and Christopher D Town²

Address: ¹French Agricultural Research Center for International Development (CIRAD), UMR 1096, Avenue Agropolis, TA40/03, FR-34398, Montpellier, Cedex 5, France, ²The J. Craig Venter Institute, 9704 Medical Center Drive, Rockville MD 20850, USA, ³Embrapa Genetic Resources and Biotechnology (CENARGEN), P B I., Final Av. W/5 Norte, Asa Norte 70770-900, Caixa-Postal: 02372, Brasília, Brazil, ⁴Genomic Sciences and Biotechnology program, Catholic University of Brasília (UCB), SGAN 916, Modulo B, Asa Norte 70790-160, Brasília DF, Brazil, ⁵Structural and Genomic Information Laboratory (IGS), C.N.R.S. UPR 2589, Institute of Structural Biology and Microbiology (IBSM), Parc Scientifique de Luminy, 163 avenue de Luminy, FR-13288 Marseille Cedex 9, France, ⁶Department of Plant Biology, University of Georgia, Athens, GA 30602, USA, ⁷Brazilian Agricultural Research Corporation (EMBRAPA), National Research Center of Cassava and Fruit Crops (CNPMF), P.O Box 007, Zip Code 44380.000, Cruz das Almas BA, Brazil, ⁸Rice Genome Research Program (RGP), National Institute of Agrobiological Sciences (NIAS)/Institute of the Society for Techno-innovation of Agriculture, Forestry and Fisheries, Tsukuba, Ibaraki 305-8602, Japan, ⁹Structural and Genomic Information Laboratory (IGS), C.N.R.S. UPR 2589, Institute of Structural Biology and Microbiology (IBSM), Parc Scientifique de Luminy, 163 avenue de Luminy, FR-13288 Marseille Cedex 9, France and ¹⁰Parco Tecnologico Padano, Via Einstein, Lodi 26900, Italy

Email: Magali Lescot* - magali.lescot@igs.cnrs-mrs.fr; Pietro Piffanelli - pietro.piffanelli@tecnoparco.org; Ana Y Ciampi - aciampi@cenargen.embrapa.br; Manuel Ruiz - manuel.ruiz@cirad.fr; Guillaume Blanc - blanc@igs.cnrs-mrs.fr; Jim Leebens-Mack - jleebensmack@plantbio.uga.edu; Felipe R da Silva - felipes@cenargen.embrapa.br; Candice MR Santos - candice@cenargen.embrapa.br; Angélique D'Hont - dhont@cirad.fr; Olivier Garsmeur - garsmeur@cirad.fr; Alberto D Vilarinhos - vila@cnpmf.embrapa.br; Hiroyuki Kanamori - kan@staff.or.jp; Takashi Matsumoto - mat@nias.affrc.go.jp; Catherine M Ronning - cronning@jcvi.org; Foo Cheung - FCheung@jcvi.org; Brian J Haas - bhaas@broad.mit.edu; Ryan Althoff - ryan_w_althoff@hotmail.com; Tammy Arbogast - TArbogast@jcvi.org; Erin Hine - ehine@jcvi.org; Georgios J Pappas - gpappas@cenargen.embrapa.br; Takuji Sasaki - tsasaki@nias.affrc.go.jp; Manoel T Souza - msouza@cenargen.embrapa.br; Robert NG Miller - rmiller@pos.ucb.br; Jean-Christophe Glaszmann - glaszmann@cirad.fr; Christopher D Town - cdtown@jcvi.org

Published: 30 January 2008

Received: 11 June 2007

BMC Genomics 2008, 9:58 doi:10.1186/1471-2164-9-58

Accepted: 30 January 2008

This article is available from: http://www.biomedcentral.com/1471-2164/9/58

© 2008 Lescot et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## Abstract

**Background:** *Musa* species (Zingiberaceae, Zingiberales) including bananas and plantains are collectively the fourth most important crop in developing countries. Knowledge concerning *Musa* genome structure and the origin of distinct cultivars has greatly increased over the last few years. Until now, however, no large-scale analyses of *Musa* genomic sequence have been conducted. This study compares genomic sequence in two *Musa* species with orthologous regions in the rice genome.

<footnote>* Corresponding author †Equal contributors</footnote>
**Results:** We produced 1.4 Mb of *Musa* sequence from 13 BAC clones, annotated and analyzed them along with 4 previously sequenced BACs. The 443 predicted genes revealed that Zingiberales genes share GC content and distribution characteristics with eudicot and Poaceae genomes. Comparison with rice revealed microsynteny regions that have persisted since the divergence of the Commelinid orders Poales and Zingiberales at least 117 Mya. The previously hypothesized large-scale duplication event in the common ancestor of major cereal lineages within the Poaceae was verified. The divergence time distributions for *Musa*-Zingiber (Zingiberaceae, Zingiberales) orthologs and paralogs provide strong evidence for a large-scale duplication event in the *Musa* lineage after its divergence from the Zingiberaceae approximately 61 Mya. Comparisons of genomic regions from *M. acuminata* and *M. balbisiana* revealed highly conserved genome structure, and indicated that these genomes diverged circa 4.6 Mya.

**Conclusion:** These results point to the utility of comparative analyses between distantly-related monocot species such as rice and *Musa* for improving our understanding of monocot genome evolution. Sequencing the genome of *M. acuminata* would provide a strong foundation for comparative genomics in the monocots. In addition a genome sequence would aid genomic and genetic analyses of cultivated *Musa* polyploid genotypes in research aimed at localizing and cloning genes controlling important agronomic traits for breeding purposes.

## Background

Taken together, *Musa* species (bananas and plantains) comprise the fourth most important crop in developing countries [1]. The fruit is a staple food in sub-Saharan Africa, South and Central America and much of Asia, while the leaves are used for sheltering and wrapping food and the male bud can be eaten as a vegetable. *Musa* is a member of the monocot order Zingiberales, a Commelinid lineage that diverged from the line leading to rice (Poales) in the mid-cretaceous period over 100 million years ago (Figure 1) [2,3]. The *Musa* species *Musa acuminata* (AA genome) and *Musa balbisiana* (BB genome), both with 2n = 22 chromosomes, represent the two main progenitors of cultivated banana varieties. Table bananas are sterile, parthenocarpic and diploids AA or triploid with the AAA genome constitution, and represent only a fraction of world production, although they are an important cash crop. Cooking bananas and plantain cultivars, mostly consumed in the countries of production, generally have an AAB or ABB genome constitutions [4]; these are boiled, fried, dried, or sometimes ground into flour.

Knowledge concerning the genetic diversity, the origin of cultivars [5-12] and *Musa* genome structure [13-15] has greatly increased over the last few years. The haploid genome of *Musa* species was estimated as varying between 560 to 600 Mb in size [16,17], just four times larger than that of the model plant Arabidopsis (125 Mb) [18] and 30% larger than that of rice (390 Mb) [19]. Genetic maps have been developed [20-23] and recently, BAC resources were generated for both *M. acuminata* [24,25] and *M. balbisiana* [26]. A cytogenetic map based on BAC-FISH is being anchored to genetic maps in order to better characterize structural variation among *M. acuminata* genomes [22]. These resources will pave the way for studies of *Musa*

genome structure and evolution through comparisons with other monocot and eudicot genomes.

The utility of genomic comparisons of monocot and eudicot plants (e.g. [27-30]) is growing with the availability of the complete genome sequences of rice [19], Arabidopsis [18] and poplar [31], and active genome sequencing projects for a growing number of other angiosperms [32]. Most genome-scale comparative investigations within the monocots have focused on analyses of closely-related species of monocots belonging to the family of Poaceae [27,33-36]. Numerous papers have described extensive microsynteny between rice, barley, wheat, maize, *Sorghum* and sugarcane [27,35,37-42], although the degree of conservation varies between different chromosomal locations. Fewer attempts have been made to investigate the synteny between distantly-related plants. In addition, whereas extensive genomic resources have been developed for rice and other cereal species in the grass family (Poaceae), there is relatively little data on gene content or genome structure for non-grass monocots (Figure 1). Recently, the first two BAC clones genomic sequences [43], and a BAC end sequencing study of the *M. acuminata* genome [44] have been published. Here we present data on the genomic structure and organization of 1.8 Mb of *Musa* genomic nuclear DNA (including the two BAC sequenced previously [43]), show for the first time the existence of microsynteny between *Musa*, rice and Arabidopsis, characterize the extent of microsynteny between the two *Musa* species representing the progenitors of most cultivated genotypes, analyze monocot EST sequences and discuss the evolutionary implication of these results. The BAC clones sequenced in this study were identified by hybridization with gene sequences previously selected to correspond to one or a few loci in *Musa*, rice and Arabi-
<figure><img src="image_3.png" /><figcaption>**Figure 1**</figcaption></figure>

**Current understanding of relationships among monocot orders [118]**. Families are shown in bold caps and genera with EST sequences in dbEST [119]. The number of sequences in GenBank (as of 10/08/07) are shown in parentheses for each order and the shaded box highlights Commelinid orders. The nodes with < 75% bootstrap support are grey.
dopsis, thus possibly contain orthologous sequences with
these distantly related plant species.

## Results

### Selection of Musa BAC clones using broad-spectrum Sorghum cDNA and Musa RFLP probes

As part of a program aiming at selecting conserved probes from monocotyledons and *Arabidopsis thaliana* towards comparative genetic mapping studies, genomic and cDNA probes from various species were tested by Southern hybridization on DNA of various monocotyledons including *Musa* and rice. Among the probes found conserved between rice, *M. acuminata* cv. Madang, *M. balbisiana* cv. Pisang Klutuk Wulung (PKW) and Arabidopsis that revealed a single or low copy locus hybridization pattern, nine were selected. These nine probes that correspond to *Sorghum bicolor* cDNAs (SbRPG) were used to screen a *M. acuminata* cv. Calcutta-4 bacterial artificial chromosome (BAC) library (Table 1 and Additional file 1). Of these nine SbRPG genes, four encoded nuclear genes targeted to the chloroplast and/or implicated in photosynthetic-related functions, supporting the notion that this class of genes is under strong pressure for functional conservation. All *Musa* BACs identified were subjected to HindIII fingerprinting. This enabled us to separate the *Musa* BACs into groups likely to be derived from different regions of the *Musa* genome. Overall, a good correlation was observed between the number of loci identified in rice, *Sorghum* and *Musa* by Southern blot, BAC fingerprint (for *Musa*) and analysis of whole genome sequence (for rice) for these nine SbRPG genes (Table 1), all of which were found to be in single or low-copy in both *Musa* and rice.

One BAC clone was selected for sequencing for probe SbRPG132. For probes SbRPG373, SbRPG661 and SbRPG851, which were found to be present in one or two copies in rice, two *Musa* BACs with distinct HindIII fingerprints that might be derived from homeologous regions were selected for sequencing with the aim of studying the evolution of lineage-specific duplications in both *Musa* and rice (Table 1). Two BAC clones from *M. acuminata* cv. Calcutta-4 (*Musa* A) and two BACs from *M. balbisiana* cv. PKW (*Musa* B) isolated using the genetically-mapped RFLP single-copy probes CIR560 and CIR257 [23] were also fully sequenced with the objective of studying the extent of synteny between *Musa* A and B species as well as against the rice genome. These RFLP probes were selected because they corresponded to genomic clones encoding genes of known function, CIR257 for a GA-20 oxidase and CIR560 for a beta 1-3 glucanase, previously shown to be associated to traits of agronomic importance in controlling plant height [45,46] and stress response [47-50], respectively.

### Analysis of 1.8 Mb of *Musa* genomic sequences reveals particular features for the *Musa* genes

#### *Musa* genome statistics

A total of 13 BACs (Table 1) were sequenced, generating over 1.4 Mb of unique *Musa* sequence. In order to provide a uniform set, data four additional BAC sequences (see Additional file 2 and [43]) were included in our annotation pipeline. These analyses revealed 443 predicted genes (on a total of 1.8 Mb of *Musa* genomic sequence from 17 BAC inserts), after elimination of all putative protein coding genes smaller than 100 amino acid residues. Approximately half of the gene models had matches in GenBank. Their classification based on similarities to genes found in the public sequence databases is presented in Additional file 3 and an annotation overview of the *Musa* genes in Additional file 4. Gene models were also compared against the *Musa* EST database donated to the Global *Musa* Genomics Consortium by Syngenta and maintained at the MIPS (Munich Information Center for Protein Sequences, Munich Germany), revealing that at least 10% of the predicted genes had a perfect match with EST sequences, thus probably being expressed in *Musa* tissues. Analysis of gene size, exon-intron structure and base composition for these 443 predicted genes is summarized in Table 2. The annotation revealed that, with the exception of MA4_78I12, the BACs analyzed were gene-rich (an average density of one gene per 4.1 kb). Our annotation of MuH9 revealed a total of 23 gene models for an average gene density of one gene per 3.6 kb compared with one gene per 6.9 kb based upon the earlier annotation [43]. In the case of MuG9, our pipeline predicted a total of 14 gene models in the first 52 kb of this BAC followed by a region of ~21 kb containing only transposons (5). Thus the gene density in the non-transposon containing region is one gene per 3.7 kb, very similar to MuH9. Previous annotation of this BAC [43] predicted 7 genes in the same 52 kb region for a density of one gene per 7 kb, with the remainder being transposon-related. The difference between the two annotations is due mainly to the larger number of hypothetical genes identified by the TIGR pipeline as well as some gene splits (e.g. MuH9-5 is split into three genes). Like the last ~21 kb of MuG9, BAC MA4_78I12 found to be mainly composed of class II transposable elements and also contains 7 interspersed predicted genes of which only the homolog of the SbRPG661 probe had a match in public databases. BAC-FISH experiments showed that BAC MA4_78I12 hybridized to all *M. acuminata* chromosomes except their extremities (Figure 2A). This pattern of hybridization is similar to what we observed by genomic in situ hybridization (GISH) using total genomic DNA as a probe and suggested that the extremities of the chromosomes are poor in repeated sequences [14]. Two gene-rich BACs (MA4_54N07 and MA4_82I11) were also analyzed by BAC-FISH and each hybridized at the extremity of one chromosome (see Figure 2B for MA4_54N07).
<table><caption>Table 1: List of probes used to identify the *Musa* BAC clones sequenced as part of the present study. Estimated copy numbers of these sequences in rice, *Sorghum* and *Musa* are indicated for SbRPG (*Sorghum bicolor*) sequences. MA4 are BAC clones from *M. acuminata* cv. Calcutta-4 and MBP are BAC clones from *M. balbisiana* cv. Pisang Klutuk Wulung.</caption><thead><tr><th>Probe name and AC number*</th><th>Putative function</th><th>Estimated copy number in rice by Blast analysis (Rice genes locus identifier)</th><th>Estimated copy number in Sorghum by Southern blot analysis</th><th>Estimated copy number in Musa by Southern blot analysis</th><th>Number of identified Musa BAC clones</th><th>Number of Musa BAC fingerprint groups</th><th>Musa BAC clones sequenced (size) and AC number*</th></tr></thead><tbody><tr><td>SbRPG132<br>DQ185891</td><td>chlorophyll A-B binding protein type I</td><td>6<br>Os01g41710.1<br>Os09g17740.1<br>Os01g52240.1<br>Os03g39610.1<br>Os07g37550.1<br>Os11g13890.2</td><td>4</td><td>more than 4</td><td>23</td><td>6</td><td>MA4_25J11<br>(105019 bp)<br>AC186746</td></tr><tr><td>SbRPG373<br>DQ185892</td><td>hypothetical protein</td><td>1<br>Os07g02340.1</td><td>1</td><td>1</td><td>29</td><td>2</td><td>MA4_64C22<br>(80932 bp)<br>AC186752<br>MA4_8L21 (115790 bp)<br>AC186748</td></tr><tr><td>SbRPG661<br>DQ185893</td><td>thioredoxin</td><td>2<br>Os10g34520.1<br>Os07g10250.1</td><td>2</td><td>2-3</td><td>20</td><td>2</td><td>MA4_54B05<br>(54106 bp)<br>AC186753<br>MA4_78I12 (150982 bp)<br>AC186750</td></tr><tr><td>SbRPG748<br>DQ185894</td><td>porphobilinogen deaminase</td><td>1<br>Os07g10250.1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>MA4_42M13<br>(29567 bp)<br>AC186749</td></tr><tr><td>SbRPG851<br>DQ185895</td><td>phosphoglycerate kinase</td><td>2<br>Os05g41640.1<br>Os01g58610.1</td><td>2</td><td>2-3</td><td>12</td><td>2</td><td>MA4_11I1I0<br>(102441 bp)<br>AC186756<br>MA4_106017<br>(143796 bp) AC186747</td></tr><tr><td>SbRPG854<br>DQ185896</td><td>mitochondrial rieske protein</td><td>2<br>Os04g32660.1<br>Os02g32120.1</td><td>2</td><td>2</td><td>5</td><td>1</td><td>MA4_11I1B14<br>(146821 bp)<br>AC186954</td></tr><tr><td>CIR257<br>DQ334868</td><td>GA-20 oxidase</td><td>-</td><td>-</td><td>-</td><td>6</td><td>1</td><td>MA4_82I11<br>(102232 bp)<br>AC186955</td></tr><tr><td>CIR560<br>DQ334869</td><td>beta 1-3 glucanase</td><td>-</td><td>-</td><td>-</td><td>21</td><td>1</td><td>MBP_81C12 (142973 bp)<br>AC186754<br>MA4_54N07<br>(96443 bp)<br>AC186751</td></tr></tbody></table>

<footnote>*AC number : accession number.</footnote>
<table><caption>**Table 2:** Features of *Musa* genes in comparison with those of Arabidopsis and rice.</caption><thead><tr><th></th><th>Musa<sup>1</sup></th><th>Arabidopsis<sup>2</sup></th><th>Rice<sup>3</sup></th></tr></thead><tbody><tr><td>GC content: overall (%)</td><td>39.0</td><td>36.0</td><td>43.5</td></tr><tr><td>Exons (%)</td><td>48.4</td><td>44.2</td><td>53.1</td></tr><tr><td>Introns (%)</td><td>38.4</td><td>32.3</td><td>38.7</td></tr><tr><td>Intergenic (%)</td><td>38.5</td><td>31.2</td><td>41.4</td></tr><tr><td>Exon length (bp)</td><td>252</td><td>276</td><td>312</td></tr><tr><td>Intron length (bp)</td><td>366</td><td>169</td><td>364</td></tr><tr><td>Number of exons/gene</td><td>4.8</td><td>5.4</td><td>4.2</td></tr><tr><td>Gene length (bp)</td><td>2,504</td><td>2,232</td><td>2,519</td></tr><tr><td>Protein length (aa)</td><td>411</td><td>417</td><td>437</td></tr><tr><td>Gene density (kb/gene)</td><td>4.1</td><td>4.5</td><td>6.2</td></tr></tbody></table>

### Base Composition and GC Distribution along the *Musa* genes

The GC content of *Musa* coding sequences was compared with those of other monocots (rice, onion, asparagus) and dicots (Arabidopsis) using two data sets -unigene clusters and -singleton ESTs found in the TIGR plant transcript assembly database (TC/ESTs; [51]) and the 443 annotated genes (CDS) from the 17 *Musa* spp. sequenced BAC clones (Figure 3). The GC distributions of TC/EST (Figure 3A) and CDS regions of the *Musa* BACs (Figure 3B) were found to be asymmetrical and bimodal as compared to Arabidopsis and onion which are clearly symmetrical and unimodal (this report, [52,53]). The *Musa* GC content distribution resembles that of rice and other Poales with higher average GC content than eudicots (see also Table 2) and a long tail towards high GC values. We next examined GC content along the direction of transcription from the ATG start codon for each predicted *Musa* CDS using a sliding window of 129 bases (Figure 4). By manual inspection of the data, we were able to identify two categories of GC profiles from the *Musa* CDS: the first set shows a marked "rice-like" gradient of GC composition from 5' to 3' end and a higher GC content than Arabidopsis all along the CDS (Figure 4A), and the second set is "Arabidopsis-like" lacking a significant GC gradient from 5' to 3' (Figure 4B).

### Analysis of *Musa* repetitive elements

Several approaches were used to characterize the genomic sequence with respect to repeats. Database searches of the predicted genes against a non-redundant protein database (see Methods section) revealed a total of 78 transposable element (TE)-related sequences. Excluding TE-rich BAC MA4_78112, there are on average ~2.6 retrotransposons (TE of class I) per 100 kb. Only one TE of class II encoded protein was detected. BAC sequences were also screened for previously characterized *Musa* RADKA repeats [15]; an average of 1.8 RADKA-related repeats (GenBank Acces-

<figure><img src="image_4.png" /><figcaption>**Figure 2**
Chromosome preparations of *M. acuminata* cv. Calcutta-4 (2n = 22) stained with DAPI after FISH of BAC. (A) MA4_78112 (detected with Texas red). (B) MA4_54N07 (detected FITC). Scale bar = 10 microns.</figcaption></figure>

sions AF399938-AF399941, AF399943-AF399946 and AF399948) per 100 kb were identified. In an attempt to identify as yet uncharacterized repeats, the BAC sequences were also analyzed by RepeatScout [54]. After removing repeats having similarity to Arabidopsis or rice proteins, *Musa* CDS, RADKA sequences and transposable elements, six repeats with at least three copies were identified (data not shown). Five of these sequences have no significant hits to genes in GenBank while the sixth matches Gen-Bank accession X99496 with a strong similarity to a part of the *Musa ycf2* chloroplast gene. Analysis of individual BACs with PrintRepeats [55] shows that each BAC contains only a small number of regions that are repeated within the BAC, an observation that is supported by the relative ease with which the BAC sequences could be closed and finished.

### Microsynteny analysis between *Musa* and either rice or Arabidopsis

The 443 *Musa* predicted proteins were aligned against the rice and Arabidopsis proteomes. The results showed that 268 and 224 *Musa* proteins have hits with an E-value threshold of 1e-10 against the rice and Arabidopsis proteomes, respectively. The relative positions of the homologous genes identified in the rice and Arabidopsis genomes were compared to the order of the corresponding *Musa* genes with i-ADHoRe software [56]. Using this stringent approach, we were able to identify nine *Musa* BAC sequences showing microsynteny among the 17 *Musa* BACs analyzed: eight cases with rice and one case with Arabidopsis (Additional file 5).

The i-ADHoRe analyses identified syntenic blocks of three to ten genes (Additional file 5). We then refined the analyses by conducting reciprocal BLASTP searches between

<footnote n="1">¹Based on 1.8 Mb of genomic sequences</footnote>

<footnote n="2">²[18]</footnote>

<footnote n="3">³[120]</footnote>
A

<figure><img src="image_5.png" /></figure>

B

<figure><img src="image_6.png" /><figcaption>**Figure 3**
Distribution of GC content in *Musa* and its comparison with other plant species. (A) All TCs/ESTs from the named species. (B) All annotated CDS from 17 *Musa* BACs (this data set) and the complete genomes of Arabidopsis and rice.</figcaption></figure>

the genes in the orthologous regions. This analysis extended the number of genes included in these syntenic blocks. The most interesting cases of synteny conservation were found between BAC MBP_91N22 and rice chromosome 1 (Figure 5A), BAC MA4_25J11 and rice chromosomes 1 and 5 (Figure 5B), BAC MA4_8L21 and rice chromosome 3 (Figure 5C), BAC MuH9 and rice chromosome 4 (Additional file 6), and BAC MA4_42M13 and rice chromosome 2 (Additional file 7). Between five and eleven genes were found in common between the syntenic *Musa* and rice orthologous regions. In most cases the common genes were found in the same order and orientation in rice and *Musa*. However, several additional genes were typically found between the shared orthologs. Interestingly, the number of genes without orthologs within otherwise syntenic regions is much higher in rice as compared to *Musa*. This could be explained by differences between the rice and *Musa* lineages in the rate of translocation, duplication and gene death. Note also that *Musa* BAC

MBP_91N22 displays conservation of synteny with two very distant segments on rice chromosome 1.

In the case of BAC MA4_25J11, two rice orthologous regions were found with i-ADHoRe analyses and reciprocal BLASTP searches (rice chromosomes 1 and 5), revealing a duplication of this region in the rice genome. It is interesting to note that the two rice orthologous regions on chromosomes 1 and 5 have lost different sets of genes compared to Musa, as has been observed previously in other duplicated regions in angiosperms [37,57,58]. Phylogenetic analyses on the 10 Musa genes from BAC MA4_25J11 and co-orthologs [59] found on rice chromosomes 1 and 5 revealed that these regions were the product of the genome-wide duplication that has been hypothesized to have occurred early in the history of the Poaceae [56,60-62]. Duplicate maize, sugarcane, Sorghum, wheat, and barley genes occur in two separate clades in trees for loci 4, 6 and 7 (Figure 6 and Additional
<figure><img src="image_7.png" /><figcaption>**Figure 4**
Mean GC content from 5' to 3' across 129 bp sliding windows. (A) for 77 *Musa* genes with a "rice-like" gradient. (B) for 180 *Musa* genes with an "Aradidopsis-like" pattern</figcaption></figure>

file 8) indicating that the duplication occurred before the divergence of the major grain lineages including rice, maize and wheat.

The only significant case of microcolinearity found between *Musa* and Arabidopsis involved three consecutive genes (Additional file 9). Interestingly, this *Musa*-Arabidopsis syntenic block was not found to be conserved in rice.

**Syntenic relationships between two regions of M. acuminata and M. balbisiana**

We also investigated conservation of synteny between two regions of the genomes of *M. acuminata* and *M. balbisiana* species. Hybrids between these two species represent the majority of cultivated *Musa* genotypes worldwide. To carry out this pilot study, we selected BACs from orthologous regions of two single-copy, genetically-mapped RFLP

probes (CIR560 and CIR257) encoding genes of agro-nomic interest. In both cases a high level of sequence conservation was found (see Figure 7 and Additional file 10) over the entire length of the sequenced regions in common between the two pairs of BACs analyzed (82.9% of nucleotide sequence identity for the MA4_82111-MBP_81C12 pair and 87.6% for the MA4_54N7-MBP_91N22 pair). The overall levels of sequence identity in genic regions were similar between the two pairs of orthologous BACs: 96.0% for the MA4_82111-MBP_81C12 pair and 96.3% for the MA4_54N7-MBP_91N22 pair (based on the aligned orthologous gene pairs defined in Table 3; see the following paragraph for further details). A high degree of synteny was found between the orthologous sequences in both gene content and gene orientation. However, we observed some incongruence between the gene predictions of the orthologous BACs whose protein products have no match in public databases (i.e. hypothetical protein genes). In contrast, the predicted structures of genes that are homologous to known sequences were largely congruent between the orthologous BACs. Given the high levels of sequence conservation between the two *Musa* species, such variation of gene structure and exon/intron boundaries is unlikely for most functional genes. Hence, this analysis supports that further validation of gene models through additional EST sequencing or targeted RT-PCR is required.

**Divergence between M. acuminata and M. balbisiana**

In order to evaluate the degree of divergence between the two *Musa* genomes, we obtained maximum likelihood estimates for Ks values comparing pairs of orthologous genes identified in the *M. acuminata* and *M. balbisiana* BACs. We restricted our analysis to those genes (detailed in Table 3) that were intact and matched known gene sequences. For example the gene model for the 14th locus in the *M. acuminata* genome (Figure 7; L14) is similar to a pectinesterase related protein, but the gene model was excluded from the analysis because the predicted coding sequence contained several in-frame stop codons indicating that this sequence is a pseudogene. The estimated Ks values ranged from 0.0231 (Additional file 10; L19) to 0.0960 (L17), far below saturation levels (i.e. Ks << 1), with an average of 0.0410 (Table 3). Applying an average synonymous substitution rate of 4.5 per 10⁹ years for nuclear genes in the Zingiberales (see below), this suggests that *M. acuminata* and *M. balbisiana* diverged approximately 4.6 Mya ago.

**Evidence for a large-scale duplication event in the *Musa* ancestor**

We also estimated the Ks values between 1,446 pairs of duplicated *Musa* genes identified among 15,661 EST-derived unigenes found to be part of the paralog sets [63]. The distribution of Ks values was estimated in order to
A

<figure><img src="image_8.png" /></figure>

B

<figure><img src="image_9.png" /></figure>

C

<figure><img src="image_10.png" /><figcaption>**Figure 5**</figcaption></figure>

**Musa-rice syntenic regions.** Predicted genes and their orientation are shown as boxed areas. Genes annotated such as hypothetical genes are represented in white. The probes used to identify the *Musa* BAC clones are indicated in brackets. Conserved genes between *Musa* and rice regions are connected by shaded areas. (A) Syntenic relationship between *Musa* MBP_91N22 BAC clone and rice chromosome 1. (B) Syntenic relationship between *Musa* MA4_25J11 BAC clone and rice chromosome 1 and 5. The numbers above the genes correspond to the locus numbers used for phylogenetic analyses. (C) Syntenic relationship between *Musa* MA4_8L21 BAC clone and rice chromosome 3.
Locus 4

<figure><img src="image_11.png" /></figure>

Locus 6

<figure><img src="image_12.png" /></figure>

Locus 7

<figure><img src="image_13.png" /><figcaption>**Figure 6**</figcaption></figure>

**Phylogenetic analyses on three of the ten M. acuminata genes from MA4_25J11 BAC clone.** These three *Musa* genes have homologous genes in rice chromosomes 1 and 5 and the locus numbers are taken from Figure 5B. Stars indicate duplication events in the most recent common ancestor of major grain lineages (i.e. rice, wheat and maize). MA4_25J11 BAC clone was isolated by the SbRPG132 probe.
<table><caption>Table 3: Level of synonymous substitution (Ks) between homologous sequences in M. acuminata and M. balbisiana.</caption><thead><tr><th>Locus</th><th>Length</th><th>K<sub>s</sub></th><th>Predicted function</th></tr></thead><tbody><tr><td>L4</td><td>1065</td><td>0.0496</td><td>GDSL-motif lipase/hydrolase family protein</td></tr><tr><td>L5</td><td>1401</td><td>0.0311</td><td>protein kinase family protein</td></tr><tr><td>L7</td><td>1341</td><td>0.0318</td><td>hypothetical protein</td></tr><tr><td>L9</td><td>1647</td><td>0.0252</td><td>protein kinase family protein</td></tr><tr><td>L11</td><td>825</td><td>0.0323</td><td>protein kinase-related</td></tr><tr><td>L13</td><td>2901</td><td>0.0420</td><td>leucine-rich repeat-containing protein kinase family protein</td></tr><tr><td>L15</td><td>1140</td><td>0.0435</td><td>gibberellin 20-oxidase family protein</td></tr><tr><td>L16</td><td>1899</td><td>0.0342</td><td>glucose-inhibited division A family protein</td></tr><tr><td>L18*</td><td>1371</td><td>0.0960</td><td>leucine rich repeat family protein</td></tr><tr><td>L20</td><td>4119</td><td>0.0231</td><td>transcriptional repressor protein-related</td></tr><tr><td>L23</td><td>1941</td><td>0.0369</td><td>protein kinase family protein</td></tr><tr><td>L24</td><td>1737</td><td>0.0413</td><td>exostosin family protein</td></tr><tr><td>L26</td><td>2145</td><td>0.0461</td><td>kinesin light chain-related</td></tr><tr><td>K<sub>s</sub> mean</td><td>23532</td><td><strong>0.0410</strong></td><td></td></tr><tr><td>Concatenation K<sub>s</sub></td><td>23532</td><td><strong>0.0349</strong></td><td></td></tr></tbody></table>

assess spikes in the accumulation of duplicated genes [64]. If we assume that gene duplications and gene deletions are random and have relatively steady rates during the course of evolution, such a distribution is expected to show an L shape [64-68]. The distribution of Ks values for duplicated *Musa* genes exhibits a large peak centered around Ks = 0.55 (Figure 8) indicating an increase in the number of gene duplications that occurred in the *Musa* ancestor *circa* 61 Mya (assuming a synonymous substitution rate of 4.5 per 10⁹ years; see below). This ancient burst of duplications is likely the result of one or more large-scale duplication events. Alternatively, the observed duplications could be associated with a burst of transposition activity as has been hypothesized for some duplicate gene pairs in Arabidopsis [69]. However, analyses of Ks plots for duplicated rice genes were unable to detect the 60 Mya duplication event in the Poaceae that is evident in analyses of gene trees and duplicated blocks in the rice genome (e.g. [70,61]; this study). This may be due in part to the slower substitution rate we estimate for the Zingiberales relative to the Poaceae (see below).

We also analyzed the 18,612 ginger (*Zingiber officinale*; Zingiberaceae, Zingiberales) EST-derived unigenes available on the TIGR Plant Transcript Assemblies web site [71] (sequences generated by David Gang, University of Arizona) and found no evidence of large-scale duplication in the Ks distribution for paralogous pairs (Figure 8). Moreover, the modal Ks for reciprocal best matches between the *Musa* and *Zingiber* unigene sets is 0.78 (Figure 8), larger than the mode for *Musa* paralogous pairs. The age of the most recent common ancestor for the Musacaceae and Zingiberaceae is estimated at 87 Mya [3,72,73]. This implies an average synonymous substitution rate of 4.5 per 10⁹ years (0.78 synonymous substitutions per site/

(2*87,000,000 years)), intermediate between rates estimated for the Poaceae (6.1–6.5 per 10⁹ years) and palms in the order Arecales (2.61 per 10⁹ years; [74]). We must emphasize that all of these rate estimates are approximate, based on rough estimates of minimum divergence times. However, regardless of ambiguities in substitution rate calibrations, our results indicate that the predicted large-scale duplication that occurred in the *Musa* lineage (Ks = 0.55) post-dates the divergence of lineages leading to *Zingiber* and *Musa* (Ks = 0.78), but occurred well before the separation of *Musa* A and *Musa* B (Ks = 0.0410).

Ks values were also computed on 1,034 pairs of homologous genes identified between the *Musa* ESTs and the rice genome sequences. As expected, the distribution of Ks values between rice-*Musa* homologs form a single peak centered around Ks = 1.7 (Figure 8). Using this Ks value to estimate the age of the Poales-Zingiberales split is less straightforward than described above for the *Musa-Zingiber* split, because synonymous substitution rates clearly vary between these Commelinid monocot lineages.

BAC fingerprint analyses revealed that whereas SbRPG854 hybridized to a single locus in the *Musa* genome, SbRPG probes SbRPG132 hybridized to 6 regions, SbRPG663 hybridized to 5 loci, and two loci were identified for SbRPG373, SbRPG661 and SbRPG851 (Table 1 and Additional file 1). BACs representing both distinct loci hybridizing to probes SbRPG661, SbRPG373 and SbRPG851 were sequenced with the aim of dating the time of duplication relative to the divergence of the *Musa* and rice lineages. Pair-wise estimations of Ks, the number of synonymous substitutions per synonymous site, were 0.93 (± 0.25), 1.39 (± 0.19) and 1.43 (± 0.60) for *Musa* homologs of the coding regions of SbRPG661 (thiore-

<footnote>* Alignments between the *M. acuminata* gene and the orthologous sequence on the *M. balbisiana* sequence were identified by BLAST</footnote>
<figure><img src="image_14.png" /><figcaption>**Figure 7**

**Comparison between *M. acuminata* MA4_82I11** BAC clone and *M. balbisiana* MBP_81C12 BAC clone.** Predicted genes and their orientation are shown as boxed areas. Genes annotated such as hypothetical genes are represented in white. The probe used to identify the *Musa* BAC clones is indicated in brackets. Conserved genes between the two *Musa* regions are connected by shaded areas. (A) Dot plot analysis of the two pairs of homeologous BACs from *M. acuminata* and *M. balbisiana*. (B) Diagram of the syntenic regions between the two BAC clones.</figcaption></figure>
<figure><img src="image_15.png" /><figcaption>**Figure 8**
Frequency of synonymous substitution (Ks) in different pairwise species comparisons. These results reveal the existence of whole genome duplications within *Musa* and revealed an extensive event pre-dating the ginger-*Musa* or rice-*Musa* divergences.</figcaption></figure>

doxin), SbRPG851 (phosphoglycerate kinase) and SbRPG373 (hypothetical protein), respectively. Phylogenetic analyses suggested that the SbRPG851 *Musa* homologs duplicated prior to the divergence of the Poales and the Zingiberales, (probably independent from the large-scale duplication described above), and the SbRPG661 and SbRPG373 *Musa* homologs are sister to each other in the gene tree, suggesting the duplications arose after the divergence of the Poales and the Zingiberales (data not shown).

We also analyzed the degree of conservation between genomic regions surrounding SbRPG661, SbRPG851 and SbRPG373 duplicated genes in *Musa* and rice and found no synteny in regions anchored by these homologs. This absence of synteny could be explained by duplication events and subsequent gene losses or by the translocation of the focal genes.

**Discussion**

*Analysis of *Musa* genes reveals some particular features*

Sequencing and annotation of ~1.8 Mb of *Musa* genomic sequence indicated that most of the BACs analyzed were gene rich with a low content of transposable element. Our analyses of 443 *Musa* genes predicted revealed that *Musa* genes generally have a "rice-like" bimodal GC distribution with a very asymmetrical and long tail towards high GC content as in previous studies [43,44]. However, a second class of "Arabidopsis-like" genes was found with an overall low GC content and no significant gradient along the coding sequence. In contrast to a previous comparison of grass and non-grass monocots [52,53], our analyses suggest that Zingiberales genes share some characteristics with the genomes of both eudicots and members of the Poaceae. This result suggests that the *Musa* genome is more similar to cereal genomes relative to onion, asparagus and the basal-most monocot lineage, *Acorus*.
**Syntenic relationships between distantly-related monocots**

Whereas widespread conservation of synteny has been
well established for members of the grass family
(Poaceae), gene order has not been generally conserved
between rice and Arabidopsis (e.g. reviewed by [75]). Few
studies have compared genome structure between the
members of the Poaceae and other monocot families, but
recent comparisons between onion, garden asparagus and
rice have failed to find evidence of conservation of macro-
or micro-synteny [76,77]. However the genomic tag
approach developed by [78] has allowed detecting anchor
points between grasses and monocots. In this study we
were able to identify microsyntenic regions in the *Musa*
and rice genomes that have persisted over some 117 mil-
lion years of evolution since these two lineages diverged
[2]. However, in all syntenic regions detected, the shared
genes were separated by intervening genes reflecting the
occurrence of numerous insertions and deletion of genes
in both rice and *Musa*. Insertions and deletions have been
observed between rice and Arabidopsis regions showing
micro-colinearity [58] and to a much lower extent
between colinear regions among Poacea genomes [37,79].
Further sequencing of the *Musa* and other monocot
genomes will provide more insight on the extent of line-
age-specific gene gain and loss in otherwise syntenic
regions.

A first insight into syntenic relationships between Musa A and B

We focused our pilot study on two genomic regions con-
taining genes of agronomic importance for *Musa* and rice
to gain insight into the extent of conservation between the
two cultivated species, *M. acuminata* (A genome) and *M.
balbisiana* (B genome). Our data revealed an extremely
high level of colinearity between the two *Musa* genomes
in both regions. However several insertions and deletions
occurred during the period of divergence (~4.6 Mya) of
the two *Musa* species. The high level of microsynteny
between the two genomes is likely to accelerate gene iso-
lation in *M. balbisiana* once the construction of the whole
genome physical map of *M. acuminata* has been com-
pleted by the Global *Musa* Genomics Consortium.

Unveiling the paleopolyploid nature of Musa species

There is accumulating data supporting that polyploidy is
one of the most important evolutionary mechanisms
influencing the structure and content of angiosperm
genomes [80]. Our work indicates ancient polyploidiza-
tion in the lineage leading to *Musa* approximately 60 Mya.
Similar lineage-specific events were described in the
Poaceae [81,82], Brassicaceae [56,83,84], *Populus* [31],
Solanaceae, Leguminoceae [64], Papaveraceae, *Acorus*, the
Magnoliids and the Nymphaceae [65]. Polyploidy has
clearly been an important source of genetic variation
across the angiosperms as retained duplicate genes typi-

cally show divergent patterns of gene expression [85,86].
In *Musa*, as in other plant species, novel phenotypes can
emerge from this genomic amalgam, including some with
high visibility to natural selection, such as organ size and
disease resistance.

Of particular interest is the "composite" nature of the
duplicated rice regions relative to the syntenic *Musa* BAC
MA4_25J11; different sets of genes were lost in rice chro-
mosome 1 and 5, respectively as compared to *Musa*. This
type of evolution is likely to reflect a dynamic of duplica-
tion [62] and independent evolution in both monocot
lineages including recurrent cycles of genome duplication
followed by diploidization. This phenomenon was also
identified by [58] in their analysis of differential gene loss
following duplication events in rice and Arabidopsis. Fur-
thermore, our phylogenetic analyses of gene sets includ-
ing the genes on *Musa* BAC MA4_25J11, rice orthologs
and related genes found in the Arabidopsis genome and
TIGR gene indices corroborate previous results suggesting
that a genome-wide duplication in the common ancestor
of all major cereal lineages is responsible for the large
duplicated segments observed in the rice genome
[61,62,87]. This finding illustrates how comparative anal-
yses of distantly-related monocot species can complement
studies on cereal genomes.

Is rice a good model to study the structure and evolution of
Musa genomes?

The use of rice as a reference species to accelerate map-
based cloning projects by extrapolating marker position
data and increasing marker density in targeted regions has
a proven efficiency among cereal crops (e.g. barley, wheat,
*Sorghum*), with a perceivable trend towards decreased effi-
ciency when phylogenetic distance increases. Our analy-
ses of the amount of microsynteny between rice and *Musa*
suggest that there are cases in which predictions based
upon microsynteny are useful but also this may not be
general. In addition although our data showed that *Musa*
genome is more similar to grain genomes relative to
onion, asparagus and the basal monocot, *Acorus*, the dif-
ferences observed confirmed that cereal genomes are not
representative of all monocots [52,53,76,77]. This work
also highlight that comparative analyses between dis-
tantly-related species such as rice and *Musa* are very
important to improve our understanding of monocot
genomes and more generally of angiosperms genome evo-
lution.

Conclusion

In conclusion, this study represents the first effort to
investigate the existence and extent of microsynteny
between rice and *Musa*, two-distantly related monocot
species. Our analyses revealed a higher degree of synteny
than has been reported for other comparisons between
the rice and species outside of the grass family. In addition, we identified evidence for an extensive microsynteny between the two *Musa* species representing the progenitors of most cultivated genotypes. In addition, we identified evidences for an ancient genome-scale duplication event in the lineage leading to *Musa* and highlighted the complexity of analyzing the structure and evolution of plant genomes following independent cycles of genome duplication and diploidization.

## Methods

### Selection of *Musa* BAC clones

Nine probes known from previous data to be conserved between rice, *Musa acuminata* cv. Madang, *Musa balbisiana* cv. PKW and Arabidopsis and revealing single or very low copy number locus were selected. These nine probes (SbRPG) correspond to *Sorghum* cDNA developed by Rus-tica Prograin Génétique and CIRAD [88]. These cDNAs and two Musa genomic probes CIR257 and CIR560 [20] were used to screen high density filters of the *M. acuminata* Calcutta-4 BAC library [25] according to standard protocols [89](see Table 1). The probes CIR257 and CIR560 were also used to screen *M. balbisiana* cv. PKW BAC library [26]. BAC DNA of positive clones was isolated using a Qiagen Robot 9600 and Qiagen 96-well BAC DNA isola-tion kit and digested with the restriction enzyme HindIII. The HindIII fingerprints were then hybridized with the corresponding probe to determine the number of loci.

### BAC-FISH analysis

Chromosome preparations were made as described in D'Hont et al [14] from root tips of *M. acuminata* cv. Calcutta-4 cultivated in glasshouse. Fluorescent *in situ* hybridizations (FISH) were performed as described in D'Hont et al [14], with 30 ng of BAC DNA labeled with digoxigenin or biotin as probes and 50 ng/µl of sheared salmon sperm DNA. The chromosomes were counterstained with DAPI (4'.6-diamidino-2-phenylindole).

### BAC sequencing

Selected BAC clones were sequenced by similar shotgun approaches at The Institute for Genomic Research (TIGR), Empresa Brasileira de Pesquisa Agropecuária (EMBRAPACENARGEN), Universidade Catolica de Brasilia (UCB) and National Institute of Agrobiological Sciences (NIAS). At TIGR, purified BAC DNA was sheared by nebulization, size-selected (2–3 kb) and ligated into a pUC-derived vector, pHOS1, using BstXI linkers. BAC DNA sent for sequencing to EMBRAPA and UCB was fragmented at Genoscope Centre (Evry, Paris, France) using a hydros hearer, size-selected (5 kb) and ligated into pcDNA2.1 vector using BstXI linkers. Clones were sequenced from both ends using ABI Big Dye terminator chemistry on ABI 3730 sequencing machines at TIGR and using a DYEnamicTM ET Terminator Sequencing Kit (Amersham

Pharmacia Biotech) on Applied Biosystems 377 sequencers at EMBRAPA and UCB. Sequences were assembled using TIGR assembler and additional directed sequencing reactions performed as necessary to complete the sequence to high quality. BAC shotgun sequencing from NIAS were performed using shotgun (2 kb and 5–7 kb) clones of 10x coverage and Big Dye Terminator Kit (ABI) on ABI 3700 sequencers, assembled with Phred/phrap software [90,91], and contig gaps were filled if necessary.

### Sequence annotation

Annotation of the BAC assemblies was carried out using the TIGR annotation pipeline, a collection of software known as Eukaryotic Genome Control (EGC) that serves as the central data management system. Each BAC sequence was processed through a series of algorithms for predicting genes (Genscan+, Genemark.hmm, Glimmer) [92-94], splice sites [95,96] and tRNAs [97]. The AAT package [98] was used for homology search against nucleotide and protein databases, that include plant-specific cDNA and EST sequences, TIGR plant gene indices [99], a non-redundant amino acid database filtered from public sources, and SwissProt [100]. Protein models generated by the searches and predictions are further searched against Markov model (HMM) databases, including PFAM [101], and automatically assigned a putative name based on domain hits or homology to previously identified proteins. Gene structures and names were manually inspected and refined as necessary. Annotated gene models were scanned for *Musa* transposable element nucleotide sequences downloaded from GenBank and then compared to a curated database of transposable element-encoded proteins [102]. The top match from each hit was used to classify the transposable element.

### Comparison of BACs with one another

In order to determine whether the BACs selected by hybridization actually arose from duplicated regions of the *M. acuminata* (A) genome or homeologous regions of the *M. balbisiana* (B) genome, or to identify duplicated regions in the *M. acuminata* (A) genome (pairs of BACs hybridizing with the same probes), each BAC was compared against all other BACs using MUMmer [103]; Dotter [104]; [105] and an all-by-all BLASTP search [106]. The sequence identity of the overlapping sequences between BACs: MA4_82111 and MBP_81C12 or MA4_54N7 and MBP_91N22, was computed with Stretcher from the EMBOSS package [107].

### Synteny search

The 443 *Musa* predicted proteins were aligned against the rice and Arabidopsis proteomes using the BLASTP program (e-value < 1e-10) [108]. The i-ADHoRe software [56] which looks for regions where the gene order is similar between two genomic sequences was used with the
following parameters: a gap size and a cluster gap of 40-40, a q value of 0.9, three anchor points and a probability cutoff of 0.001. For four BACs (MA4_25J11, MA4_8L21, MBP_91N22 and MuH9), we tried to extend the regions of synteny between *Musa* and rice found by i-ADHoRe by conducting reciprocal BLASTP searches between the genes corresponding on the homologous regions.

### Phylogenetic analyses

*Musa* genes were used in BLASTX searches to query a database of rice and Arabidopsis gene family clusters [109]. Translated blast searches (tBLASTX) against the TIGR plant gene indices [110] were also performed and inferred protein sequences with e-values < 1e-30 were compiled with homologous *Musa*, rice and Arabidopsis sequences. Amino acid alignments of the compiled sequences were constructed using MUSCLE [111] and manually adjusted. Parsimony analyses were performed on the amino acid alignments using PAUP* v4.0b10 [112].

### Construction of unigenes

*Musa* EST sequences were provided by the Global *Musa* Genomics Consortium [113]. These sequences were first assembled into unigenes using the TGICL package [114] to eliminate sequence redundancy. Because unigenes are derived from EST sequences and so have no annotated open reading frames and may contain frameshift sequencing errors, the following approach was taken. Each unigene was aligned against the rice proteome (downloaded from GenBank) using BLASTX. The best match was considered significant if the alignment length was >100 amino acids and the Expect value (E) was <1e-15. The open reading frame was then extracted from the unigene sequence using the Genewise program (which can infer frameshift sites; [115] with the corresponding best match protein as a guide.

### Estimation of the level of synonymous substitution between two sequences

For each pair of coding sequence, the two translation products were aligned using the MUSCLE program [111] and the resulting alignment was used as a guide to align the nucleotide sequences. After removing gaps and N-containing codons, the level of synonymous substitution ($K_s$) was estimated using the maximum likelihood method implemented in CODEML [116] under the F3x4 model [117].

### Distribution of the age of duplication of *Musa* genes

All-against-all nucleotide sequence similarity searches were done among the open reading frame extracted from the unigene sequences using BLASTN [106]. Sequences aligned over >300 bp and showing at least 40% identity were defined as pairs of paralogs. Then we estimated $K_s$ for each pair of paralogs. We systematically discarded one

sequence from a pair of paralogs showing no synonymous substitutions ($K_s$ = 0) as well as all $K_s$ values involving this sequence to avoid the inclusion of redundant entries of the same gene in the analysis (see [64] for further details). A gene family of *n* members results from $n-1$ gene duplication events. However, the number of possible pairwise comparisons within a gene family ($n \times (n-1)/2$) can be substantially larger than the number of gene duplications, which results in multiple estimates of the ages of some duplications. To eliminate the redundant $K_s$ values, pairs of duplicated sequences were grouped into gene families using a single linkage clustering method. Then we used the hierarchical clustering method described in [64] to reconstruct the approximated phylogeny of each gene family: (1) Initially, all sequences in the family were treated as a separate clusters. (2) Then, the $K_s$ values for all possible pairs of clusters were compared. (3) The pair of clusters having the smallest $K_s$ value was replaced by a single new cluster containing all their sequences. (4) The median $K_s$ value was chosen to represent the duplication event that gave rise to the two merged clusters. (5) Steps 2 to 4 were repeated until all sequences were contained in a single cluster. When two clusters A and B contained more than one sequence, their associated $K_s$ value in step 2 corresponded to the median $K_s$ obtained for all possible pairs of any sequence from A and any sequence from B.

### Abbreviations

aa – amino acid

AC number – accession number

CDS – coding DNA sequence

EST – expressed sequence tag

FISH – fluorescent *in situ* hybridization

GISH – genomic *in situ* hybridization

$K_s$ – synonymous substitution rate

Mya – million years ago

PKW – Pisang Klutuk Wulung

RFLP – restriction fragment length polymorphism

TE – transposable element

TC – tentative consensi.

### Authors' contributions

PP, AD, GJP, TS, MTSJ, RNGM and JCG conceived of the study and participated in its design. PP, AYC, CMRS, OG,
ADV, HK, TM, RA, TA, EH, GJP, RNGM and CDT performed the experiments. ML, MR, GB, JLM, FRDS, CMR, FC, BJH and CDT analysed and interpreted the data. RNGM contributed reagents. ML, GB, JLM, AD and CDT wrote the manuscript. All authors read and approved the final manuscript.

## Additional material

### Additional file 1

**Supplementary Table 1. Additional list of probes used to identify the Musa BAC clones.** Estimated copy numbers of these sequences in rice, Sorghum and Musa are indicated for SbRPG (Sorghum bicolor) sequences.

Click here for file
[http://www.biomedcentral.com/content/supplementary/1471-2164-9-58-S1.doc]

### Additional file 2

**Supplementary Table 2. Additional BAC clones analyzed to define Musa gene features and syntenic relationships with rice.**

Click here for file
[http://www.biomedcentral.com/content/supplementary/1471-2164-9-58-S2.doc]

### Additional file 3

**Supplementary Table 3. Statistics of the 17 Musa BAC clones analyzed in the present study.**

Click here for file
[http://www.biomedcentral.com/content/supplementary/1471-2164-9-58-S3.doc]

### Additional file 4

**Supplementary Table 4. Annotation overview of the Musa genes.** The column «Pseudogene» indicates by [1] if the gene is a pseudogene and [0] if not. Closest sequence homolog is the first similar protein sequence found by BLASTP after the sequence itself.

Click here for file
[http://www.biomedcentral.com/content/supplementary/1471-2164-9-58-S4.doc]

### Additional file 5

**Supplementary Table 5. List of genes involved in synteny relationship between Musa and rice based on i-ADHoRE results.** Multiplicon is a BAC genomic sequence on which the baseclusters are isolated and represents a cluster of 3 genes minima.

Click here for file
[http://www.biomedcentral.com/content/supplementary/1471-2164-9-58-S5.doc]

### Additional file 6

**Supplementary Figure 1.** Musa-rice syntenic region between MuH9 BAC clone and rice chromosome 4. Homologous genes between Musa and rice are indicated by shaded areas. Genes annotated such as hypothetical genes are white.

Click here for file
[http://www.biomedcentral.com/content/supplementary/1471-2164-9-58-S6.ppt]

### Additional file 7

**Supplementary Figure 2.** The Musa-rice syntenic region around the highly conserved porphobilinogen deaminase gene. Shaded areas connect homologous genes conserved between chromosome 2 of rice and Musa MA4_42M13 BAC clone (isolated by SbRPG748 probe). Genes annotated such as hypothetical genes are white. In Musa, a recent local duplication of the porphobilinogen deaminase gene occurred (genes MA4_42M13.6 and MA4_42M13.8).

Click here for file
[http://www.biomedcentral.com/content/supplementary/1471-2164-9-58-S7.ppt]

### Additional file 8

**Supplementary Figure 3.** Phylogenetic analyses on the seven of the ten M. acuminata genes from MA4_25J11 BAC clone. These seven Musa genes have homologous genes in rice chromosomes 1 and 5 and the locus numbers are available on Figure 5B. MA4_25J11 BAC clone was isolated by SbRPG132 probe.

Click here for file
[http://www.biomedcentral.com/content/supplementary/1471-2164-9-58-S8.ppt]

### Additional file 9

**Supplementary Figure 4.** Musa-Arabidopsis syntenic region between Musa MA4_54B05 BAC clone and Arabidopsis chromosome 5. Homologous genes between Musa and Arabidopsis are indicated by shaded areas. Genes annotated such as hypothetical genes are white. MA4_54B05 BAC clone was isolated by SbRPG661 probe.

Click here for file
[http://www.biomedcentral.com/content/supplementary/1471-2164-9-58-S9.ppt]

### Additional file 10

**Supplementary Figure 5.** Collinearity between M. acuminata (MA4_54N07) and M. balbisiana (MBP_91N22) around the CIR560 marker. The shaded areas connecting the two genomic regions represent conserved genes. Predicted genes and their orientation in each Musa BAC clone are shown as boxed areas. The genes for which the name is in bold hybridize with the marker. Genes annotated such as hypothetical genes are white. (A) Dot plot analysis of the two pairs of homeologous BACs from M. acuminata and M. balbisiana. (B) Diagram of the syntenic regions between the two BAC clones.

Click here for file
[http://www.biomedcentral.com/content/supplementary/1471-2164-9-58-S10.ppt]

## Acknowledgements

We thank the Montpellier Languedoc-Roussillon Genopole® for hosting the BAC library production and screening. We thank the Genoscope Centre in Evry, Paris, France for assisting AYC to carry out the subcloning of the five BAC clones sequenced at EMBRAPA and UCB.

Access to the Syngenta Musa EST database, donated by Syngenta to the International Network for the Improvement of Banana and Plantain (INIBAP) for use within the framework of the Global Musa Genomics Consortium is acknowledged.

This work was supported by CIRAD, INIBAP, NIAS, EMBRAPA, UCB, the National Council for Scientific and Technological Development (CNPq) in Brazil, TIGR and Generation Challenge program.
## References

1. Arias P, Dankers C, Liu P, Pilkauskas P: The world banana economy 1985–2002. FAO 2003 [http://www.fao.org/documents/show_cdr.asp?url_file=/docrep/007/y5102e/y5102e00.htm].

2. Janssen T, Bremer K: The age of major monocot groups inferred from 800 + rbcL sequences. *Botanical Journal of the Linnean Society* 2004, 146:385-398.

3. Sanderson MJ, Thorne JL, Wikström N, Bremer K: Molecular evidence on plant divergence times. *American Journal of Botany* 2004, 91(1656–1665):.

4. Simmonds N, Shepherd K: The taxonomy and origins of the cultivated bananas. *Bot J Linn Soc* 1955, 55:302-312.

5. Bartos J, Alkhimova O, Dolezelova M, De Langhe E, Dolezel J: Nuclear genome size and genomic distribution of ribosomal DNA in Musa and Ensete (Musaceae): taxonomic implications. *Cytogenet Genome Res* 2005, 109(1–3):50-57.

6. Carreel F, Fauré S, Gonzalez de Leon D, Lagoda PJL, Perrier X, Bakry F, Tezenas du Montcel H, Lanaud C, Horry JP: Evaluation of the genetic diversity in diploid bananas (Musa sp.). *Genetics, Selection, Evolution* 1994, 26:125s-136s.

7. Carreel F, Gonzalez de Leon D, Lagoda P, Lanaud C, Jenny C, Horry JP, Tezenas du Montcel H: Ascertaining maternal and paternal lineage within Musa by chloroplast and mitochondrial DNA RFLP analyses. *Genome* 2002, 45(4):679-692.

8. Grapin A, Noyer JL, Dambier D, Carreel F, Lanaud C, Baurens F-C, Lagoda PJL: Diploid Musa acuminata genetic diversity with Sequence Tagged Microsatellite Sites. *Electrophoresis* 1998, 19:1374-1380.

9. Noyer JL, Causse S, Tomekpke K, Bouet A, Baurens FC: A new image of plantain diversity assessed by SSR, AFLP and MSAP markers. *Genetica* 2005, 124(1):61-69.

10. Raboin LM, Carreel F, Noyer J-L, Baurens F-C, Horry J-P, Bakry F, Tezenas Du Montcel H, Ganry J, Lanaud C, Lagoda PJL: Diploid Ancestors of Triploid Export Banana Cultivars: Molecular Identification of 2n Restitution Gamete Donors and n Gamete Donors. *Molecular breeding* 2005, 16(4):333-341.

11. Ude G, Pillay M, Nwakanma D, Tenkouano A: Genetic Diversity in Musa acuminata Colla and Musa balbisiana Colla and some of their natural hybrids using AFLP Markers. *Theor Appl Genet* 2002, 104(8):1246-1252.

12. Ge XJ, Liu MH, Wang K, Schaal BA, Chiang TY: Population structure of wild bananas, Musa balbisiana, in China determined by SSR fingerprinting and cpDNA PCR-RFLP. *Molecular ecol-ogy* 2005, 14:933-944.

13. Baurens FC, Noyer JL, Lanaud C, Lagoda PJ: Use of competitive PCR to assay copy number of repetitive elements in banana. *Mol Gen Genet* 1996, 253(1–2):57-64.

14. D'Hont A, Paget-Goy A, Escoute J, Carreel F: The interspecific genome structure of cultivated banana, Musa spp. revealed by genomic DNA in situ hybridization. *Theor Appl Genet* 2000, 100:177-183.

15. Valarik M, Simkova H, Hribova E, Safar J, Dolezelova M, Dolezel J: Isolation, characterization and chromosome localization of repetitive DNA sequences in bananas (Musa spp.). *Chromosome Res* 2002, 10(2):89-100.

16. Kamate K, Brown S, Durand P, Bureau JM, De Nay D, Trinh TH: Nuclear DNA content and base composition in 28 taxa of Musa. *Genome* 2001, 44(4):622-627.

17. Lysak M, Dolezelova M, Horry J, Swennen R, Dolezel J: Flow cytometric analysis of nuclear DNA content in Musa. *Theor Appl Genet* 1999, 98:1344-1350.

18. Arabidopsis Genome Initiative: Analysis of the genome sequence of the flowering plant Arabidopsis thaliana. *Nature* 2000, 408(6814):796-815.

19. International Rice Genome Sequencing Project: The map-based sequence of the rice genome. *Nature* 2005, 436(7052):793-800.

20. Fauré S, Noyer J, Horry J, Bakry F, Lanaud C, Gonzalez D, Leon D: A molecular marker-based linkage map of diploid bananas (Musa acuminata). *Theor Appl Genet* 1993, 87:517-526.

21. Noyer J, Dambier D, Lanaud C, Lagoda P: The saturated map of diploid banana (Musa acuminata). *Abstract Plant & Animal Genome V Conference* 1997.

22. Vilarinhos A, Carreel F, Rodier M, Hippolyte I, Benabdelmouna A, Triaire D, Bakry F, Courtois B, D'Hont A: Characterization Of Translocations In Banana By FISH Of BAC Clones Anchored

**To A Genetic Map.** In *Plant & Animal Genomes XIV Conference San Diego, CA: 2006*. January 14–18, 2006

23. Tropgenedb [http://tropgenedb.cirad.fr/en/banana.html]

24. Ortiz-Vazquez E, Kaemmer D, Zhang HB, Muth J, Rodriguez-Mendíola M, Arias-Castro C, James A: Construction and characterization of a plant transformation-competent BIBAC library of the black Sigatoka-resistant banana Musa acuminata cv. Tuu Gia (AA). *Theor Appl Genet* 2005, 110(4):706-713.

25. Vilarinhos AD, Piffanelli P, Lagoda P, Thibivilliers S, Sabau X, Carreel F, D'Hont A: Construction and characterization of a bacterial artificial chromosome library of banana (Musa acuminata Colla). *Theor Appl Genet* 2003, 106(6):1102-1106.

26. Safar J, Noa-Carrazana JC, Vrana J, Bartos J, Alkhimova O, Sabau X, Simkova H, Lheureux F, Caruana ML, Dolezel J, et al.: Creation of a BAC resource to study the structure and evolution of the banana (Musa balbisiana) genome. *Genome* 2004, 47(6):1182-1191.

27. Devos KM: Updating the 'crop circle'. *Curr Opin Plant Biol* 2005, 8(2):155-162.

28. Mudge J, Cannon SB, Kalo P, Oldroyd GE, Roe BA, Town CD, Young ND: Highly syntenic regions in the genomes of soybean, Medicago truncatula, and Arabidopsis thaliana. *BMC Plant Biol* 2005, 5(1):15.

29. Town CD, Cheung F, Maiti R, Crabtree J, Haas BJ, Wortman JR, Hine EE, Althoff R, Arbogast TS, Tallon LJ, et al.: Comparative Genomics of Brassica oleracea and Arabidopsis thaliana Reveal Gene Loss, Fragmentation, and Dispersal after Polyploidy. *Plant Cell* 2006, 18(6):1348-1359.

30. Zhu H, Choi HK, Cook DR, Shoemaker RC: Bridging model and crop legumes through comparative genomics. *Plant Physiol* 2005, 137(4):1189-1196.

31. Tuskan GA, Difazio S, Jansson S, Bohlmann J, Grigoriev I, Hellsten U, Putnam N, Ralph S, Rombauts S, Salamov A, et al.: The genome of black cottonwood, Populus trichocarpa (Torr. & Gray). *Science* 2006, 313(5793):1596-1604.

32. Jackson S, Rounsley S, Purugganan M: Comparative sequencing of plant genomes: choices to make. *Plant Cell* 2006, 18(5):1100-1104.

33. Bowers JE, Arias MA, Asher R, Avise JA, Ball RT, Brewer GA, Buss RW, Chen AH, Edwards TM, Estill JC, et al.: Comparative physical mapping links conservation of microsynteny to chromosome structure and recombination in grasses. *Proc Natl Acad Sci USA* 2005, 102(37):13206-13211.

34. Buell CR, Yuan Q, Ouyang S, Liu J, Zhu W, Wang A, Maiti R, Haas B, Wortman J, Pertea M, et al.: Sequence, annotation, and analysis of synteny between rice chromosome 3 and diverged grass species. *Genome Res* 2005, 15(9):1284-1291. Epub 2005 Aug 1218

35. La Rota M, Sorrells ME: Comparative DNA sequence analysis of mapped wheat ESTs reveals the complexity of genome relationships between rice and wheat. *Funct Integr Genomics* 2004, 4(1):34-46.

36. Singh NK, Raghuvanshi S, Srivastava SK, Gaur A, Pal AK, Dalal V, Singh A, Ghazi IA, Bhargav A, Yadav M, et al.: Sequence analysis of the long arm of rice chromosome 11 for rice-wheat synteny. *Funct Integr Genomics* 2004, 4(2):102-117. Epub 2004 Apr 2014

37. Ilic K, SanMiguel PJ, Bennetzen JL: A complex history of rearrangement in an orthologous region of the maize, sorghum, and rice genomes. *Proc Natl Acad Sci USA* 2003, 100(21):12265-12270. Epub 12003 Oct 12266

38. Gu Y, Coleman-Derr D, Kong X, Anderson O: Rapid genome evolution revealed by comparative sequence analysis of orthologous regions from four triticeae genomes. *Plant Physiol* 2004, 135(1):459-470.

39. Jaillon O, Aury J, Brunet F, Petit J, Stange-Thomann N, Mauceli E, Bouneau L, Fischer C, Ozouf-Costaz C, Bernot A, et al.: Genome duplication in the teleost fish Tetraodon nigroviridis reveals the early vertebrate proto-karyotype. *Nature* 2004, 431(7011):946-957.

40. Jannoo N, Grivet L, Chantret N, Garsmeur O, Glaszmann JC, Arruda P, D'Hont A: Orthologous comparison in a gene-rich region among grasses reveals stability in the sugarcane polyploid genome. *Plant Journal* 2007 in press.

41. Salse J, Piegu B, Cooke R, Delseny M: Synteny between Arabidopsis thaliana and rice at the genome level: a tool to identify conservation in the ongoing rice genome sequencing project. *Nucleic Acids Res* 2002, 30(11):2316-2328.
42. Salse J, Piegu B, Cooke R, Delseny M: New in silico insight into the synteny between rice (Oryza sativa L.) and maize (Zea mays L.) highlights reshuffling and identifies new duplications in the rice genome. *Plant J* 2004, **38**(3):396-409. Erratum in: *Plant J*. 2004 Jun;2038(2005):2873

43. Aert R, Sagi L, Volckaert G: Gene content and density in banana (Musa acuminata) as revealed by genomic sequencing of BAC clones. *Theor Appl Genet* 2004, **109**(1):129-139.

44. Cheung F, Town CD: A BAC end view of the Musa acuminata genome. *BMC Plant Biol* 2007, **7**(29):29.

45. Oikawa T, Koshioka M, Kojima K, Yoshida H, Kawata M: A role of OsGA20ox1, encoding an isoform of gibberellin 20-oxidase, for regulation of plant stature in rice. *Plant Mol Biol* 2004, **55**(5):687-700.

46. Spielmeyer W, Ellis MH, Chandler PM: Semidwarf (sd-1), "green revolution" rice, contains a defective gibberellin 20-oxidase gene. *Proc Natl Acad Sci USA* 2002, **99**(13):9043-9048.

47. Nishizawa Y, Saruta M, Nakazono K, Nishio Z, Soma M, Yoshida T, Nakajima E, Hibi T: Characterization of transgenic rice plants over-expressing the stress-inducible beta-glucanase gene Gns1. *Plant Mol Biol* 2003, **51**(1):143-152.

48. Thomas BR, Romero GO, Nevins DJ, Rodriguez RL: New perspectives on the endo-beta-glucanases of glycosyl hydrolase Family 17. *Int J Biol Macromol* 2000, **27**(2):139-144.

49. Romero GO, Simmons C, Yaneshita M, Doan M, Thomas BR, Rodriguez RL: Characterization of rice endo-beta-glucanase genes (Gns2-Gns14) defines a new subgroup within the gene family. *Gene* 1998, **223**(1-2):311-320.

50. Simmons CR, Litts JC, Huang N, Rodriguez RL: Structure of a rice beta-glucanase gene regulated by ethylene, cytokinin, wounding, salicylic acid and fungal elicitors. *Plant Mol Biol* 1992, **18**(1):33-45.

51. Childs KL, Hamilton JP, Zhu W, Ly E, Cheung F, Wu H, Rabinowicz PD, Town CD, Buell CR, Chan AP: The TIGR Plant Transcript Assemblies database. *Nucleic Acids Res* 2007:D846-851.

52. Kuhl JC, Cheung F, Yuan Q, Martin W, Zewdie Y, McCallum J, Cataacha N, Rutherford P, Sink KC, Jenderek M, et al.: A unique set of 11,008 onion expressed sequence tags reveals expressed sequence and genomic differences between the monocot orders Asparagales and Poales. *Plant Cell* 2004, **16**(1):114-125. Epub 2003 Dec 11

53. Kuhl JC, Havey MJ, Martin WJ, Cheung F, Yuan Q, Landherr L, Hu Y, Leebens-Mack J, Town CD, Sink KC: Comparative genomic analyses in Asparagus. *Genome* 2005, **48**:1052-1060.

54. Price AL, Jones NC, Pevzner PA: De novo identification of repeat families in large genomes. *Bioinformatics* 2005, **21**(Suppl 1):i351-358.

55. Parsons JD: Miropeats: graphical DNA sequence comparisons. *Comput Appl Biosci* 1995, **11**(6):615-619.

56. Simmillion C, Vandepoele K, Saeys Y, Van de Peer Y: Building genomic profiles for uncovering segmental homology in the twilight zone. *Genome Res* 2004, **14**(6):1095-1106.

57. Sampedro J, Lee Y, Carey RE, dePamphilis C, Cosgrove DJ: Use of genomic history to improve phylogeny and understanding of births and deaths in a gene family. *Plant J* 2005, **44**(3):409-419.

58. Vandepoele K, Simmillion C, Van de Peer Y: Detecting the undetectable: uncovering duplicated segments in Arabidopsis by comparison with rice. *Trends Genet* 2002, **18**(12):606-608.

59. Sonnhammer EL, Koonin EV: Orthology, paralogy and proposed classification for paralog subtypes. *Trends Genet* 2002, **18**(12):619-620.

60. Blanc G, Barakat A, Guyot R, Cooke R, Delseny M: Extensive duplication and reshuffling in the Arabidopsis genome. *Plant Cell* 2000, **12**(7):1093-1101.

61. Paterson AH, Bowers JE, Chapman BA: Ancient polyploidization predating divergence of the cereals, and its consequences for comparative genomics. *Proc Natl Acad Sci USA* 2004, **101**(26):9903-9908.

62. Yu J, Wang J, Lin W, Li S, Li H, Zhou J, Ni P, Dong W, Hu S, Zeng C, et al.: The Genomes of Oryza sativa: a history of duplications. *PLoS Biol* 2005, **3**(2):e38.

63. A Global Programme for Musa Improvement [http://www.promusa.org/].

64. Blanc G, Wolfe KH: Widespread paleopolyploidy in model plant species inferred from age distributions of duplicate genes. *Plant Cell* 2004, **16**(7):1667-1678. Epub 2004 Jun 1618

65. Cui L, Wall PK, Leebens-Mack JH, Lindsay BG, Soltis DE, Doyle JJ, Soltis PS, Carlson JE, Arumuganathan K, Barakat A, et al.: Widespread genome duplications throughout the history of flowering plants. *Genome Res* 2006, **16**(6):738-749.

66. Lynch M, Conery JS: The evolutionary fate and consequences of duplicate genes. *Science* 2000, **290**(5494):1151-1155.

67. Maere S, De Bodt S, Raes J, Casneuf T, Van Montagu M, Kuiper M, Van de Peer Y: Modeling gene and genome duplications in eukaryotes. *Proc Natl Acad Sci USA* 2005, **102**(15):5454-5459.

68. Schlueter SD, Wilkerson MD, Huala E, Rhee SY, Brendel V: Community-based gene structure annotation. *Trends Plant Sci* 2005, **10**(1):9-14.

69. Hughes AL, Friedman R, Ekollu V, Rose JR: Non-random association of transposable elements with duplicated genomic blocks in Arabidopsis thaliana. *Mol Phylogenet Evol* 2003, **29**(3):410-416.

70. Vandepoele K, Simmillion C, Van de Peer Y: Evidence that rice and other cereals are ancient aneuploids. *Plant Cell* 2003, **15**(9):2192-2202.

71. The TIGR Plant Transcript Assemblies database [http://plantta.tiger.org/].

72. Kress WJ, Prince LM, Hahn WJ, Zimmer EA: Unraveling the evolutionary radiation of the families of the Zingiberales using morphological and molecular evidence. *Syst Biol* 2001, **50**(6):926-944.

73. Bremer K: Early Cretaceous lineages of monocot flowering plants. *Proc Natl Acad Sci USA* 2000, **97**(9):4707-4711.

74. Gaut BS, Morton BR, McCaig BC, Clegg MT: Substitution rate comparisons between grasses and palms: synonymous rate differences at the nuclear gene Adh parallel rate differences at the plastid gene rbcL. *Proc Natl Acad Sci USA* 1996, **93**(19):10274-10279.

75. Bennetzen JL, Ma J, Devos KM: Mechanisms of recent genome size variation in flowering plants. *Ann Bot* (Lond) 2005, **95**(1):127-132.

76. Martin WJ, McCallum J, Shigyo M, Jakse J, Kuhl JC, Yamane N, Pither-Joyce M, Gokce AF, Sink KC, Town CD, et al.: Genetic mapping of expressed sequences in onion and in silico comparisons with rice show scant colinearity. *Mol Genet Genomics* 2005:1-8.

77. Jakse J, Telgmann A, Jung C, Khar A, Melgar S, Cheung F, Town CD, Havey MJ: Comparative sequence and genetic analyses of asparagus BACs reveal no microsynteny with onion or rice. *Theor Appl Genet* 2006, **114**(1):31-39.

78. Lohithaswa HC, Feltus FA, Singh HP, Bacon CD, Bailey CD, Paterson AH: Leveraging the rice genome sequence for monocot comparative and translational genomics. *Theor Appl Genet* 2007, **115**(2):237-243.

79. Dubcovsky J, Ramakrishna W, SanMiguel PJ, Busso CS, Yan L, Shiloff BA, Bennetzen JL: Comparative sequence analysis of colinear barley and rice bacterial artificial chromosomes. *Plant Physiol* 2001, **125**(3):1342-1353.

80. Adams KL, Wendel JF: Polyploidy and genome evolution in plants. *Curr Opin Plant Biol* 2005, **8**(2):135-141.

81. Paterson AH, Bowers JE, Peterson DG, Estill JC, Chapman BA: Structure and evolution of cereal genomes. *Curr Opin Genet Dev* 2003, **13**(6):644-650.

82. Wang X, Shi X, Hao B, Ge S, Luo J: Duplication and DNA segmental loss in the rice genome: implications for diploidization. *New Phytol* 2005, **165**(3):937-946.

83. Vision TJ, Brown DG, Tanksley SD: The origins of genomic duplications in Arabidopsis. *Science* 2000, **290**(5499):2114-2117.

84. Bowers JE, Abbey C, Anderson S, Chang C, Draye X, Hoppe AH, Jesup R, Lemke C, Lennington J, Li Z, et al.: A high-density genetic recombination map of sequence-tagged sites for sorghum, as a framework for comparative structural and evolutionary genomics of tropical grains and grasses. *Genetics* 2003, **165**(1):367-386.

85. Blanc G, Wolfe KH: Functional divergence of duplicated genes formed by polyploidy during Arabidopsis evolution. *Plant Cell* 2004, **16**(7):1679-1691.

86. Duarte JM, Cui L, Wall PK, Zhang Q, Zhang X, Leebens-Mack J, Ma H, Altman N, dePamphilis CW: Expression pattern shifts following duplication indicative of subfunctionalization and neofunctionalization in regulatory genes of Arabidopsis. *Mol Biol Evol* 2006, **23**(2):469-478.
87. Rong J, Abbey C, Bowers JE, Brubaker CL, Chang C, Chee PW, Delmonte TA, Ding X, Garza JJ, Marler BS, et al.: A 3347-locus genetic recombination map of sequence-tagged sites reveals features of genome organization, transmission and evolution of cotton (Gossypium). *Genetics* 2004, **166**(1):389-417.

88. Boivin K, Deu M, Rami J-F, Trouche G, Hamon P: Towards a saturated sorghum map using RFLP and AFLP markers. *Theor Appl Genet* 1999, **98**:320.

89. Luo M, Wang YH, Frisch D, Joobeur T, Wing RA, Dean RA: Melon bacterial artificial chromosome (BAC) library construction using improved methods and identification of clones linked to the locus conferring resistance to melon Fusarium wilt (Fom-2). *Genome* 2001, **44**:154-116.

90. Ewing B, Green P: Base-calling of automated sequencer traces using phred. II. Error probabilities. *Genome Res* 1998, **8**(3):186-194.

91. Ewing B, Hillier L, Wendl MC, Green P: Base-calling of automated sequencer traces using phred. I. Accuracy assessment. *Genome Res* 1998, **8**(3):175-185.

92. Burge C, Karlin S: Prediction of complete gene structures in human genomic DNA. *J Mol Biol* 1997, **268**(1):78-94.

93. Lukashin AV, Borodovsky M: GeneMark.hmm: new solutions for gene finding. *Nucleic Acids Res* 1998, **26**(4):1107-1115.

94. Salzberg SL, Pertea M, Delcher AL, Gardner MJ, Tettelin H: Interpolated Markov models for eukaryotic gene finding. *Genomics* 1999, **59**(1):24-31.

95. Pertea M, Lin X, Salzberg SL: GeneSplicer: a new computational method for splice site prediction. *Nucleic Acids Res* 2001, **29**(5):1185-1190.

96. Hebsgaard SM, Korning PG, Tolstrup N, Engelbrecht J, Rouze P, Brunak S: Splice site prediction in Arabidopsis thaliana pre-mRNA by combining local and global sequence information. *Nucleic Acids Res* 1996, **24**(17):3439-3452.

97. Lowe TM, Eddy SR: tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. *Nucleic Acids Res* 1997, **25**(5):955-964.

98. Huang X, Adams MD, Zhou H, Kerlavage AR: A tool for analyzing and annotating genomic sequences. *Genomics* 1997, **46**(1):37-45.

99. Quackenbush J, Liang F, Holt I, Pertea G, Upton J: The TIGR gene indices: reconstruction and representation of expressed gene sequences. *Nucleic Acids Res* 2000, **28**(1):141-145.

100. Bairoch A, Apweiler R: The SWISS-PROT protein sequence database and its supplement TrEMBL in 2000. *Nucleic Acids Res* 2000, **28**(1):45-48.

101. Bateman A, Birney E, Cerruti L, Durbin R, Etwiller L, Eddy SR, Griffiths-Jones S, Howe KL, Marshall M, Sonnhammer EL: The Pfam protein families database. *Nucleic Acids Res* 2002, **30**(1):276-280.

102. Transposable elements database on TIGR FTP site [ftp://ftp.tigr.org/pub/data/TransposableElements/transposon_db.pep]

103. Delcher AL, Phillippy A, Carlton J, Salzberg SL: Fast algorithms for large-scale genome alignment and comparison. *Nucleic Acids Res* 2002, **30**(11):2478-2483.

104. Sonnhammer EL, Durbin R: A dot-matrix program with dynamic threshold control suited for genomic DNA and protein sequence analysis. *Gene* 1995, **167**(1-2):GCI-10.

105. Dotter web site [http://bioinfo.hku.hk/doc/dotter.html]

106. Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. *Nucleic Acids Res* 1997, **25**(17):3389-3402.

107. Rice P, Longden I, Bleasby A: EMBOSS: the European Molecular Biology Open Software Suite. *Trends Genet* 2000, **16**(6):276-277.

108. Altschul S, Gish W, Miller W, Myers E, Lipman D: Basic local alignment search tool. *J Mol Biol* 1990, **215**(3):403-410.

109. The Floral Genome Project - PlantTribes [http://fgp.dev.huck.psu.edu/tribe.pl]

110. The TIGR plant gene indices web site [http://www.tigr.org/tdb/tgi/plant.shtml]

111. Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. *Nucleic Acids Res* 2004, **32**(5):1792-1797. Print 2004

112. Swofford DL: PAUP*. Phylogenetic Analysis Using Parsimony (*and Other Methods*). Version 4. *Sinauer Associates, Sunderland, Massachusetts* 2003.

113. The Global Musa Genomics Consortium [http://www.musage.nomics.org]

114. Pertea G, Huang X, Liang F, Antonescu V, Sultana R, Karamycheva S, Lee Y, White J, Cheung F, Parvizi B, et al.: TIGR Gene Indices clustering tools (TGICL): a software system for fast clustering of large EST datasets. *Bioinformatics* 2003, **19**(5):651-652.

115. Birney E, Thompson JD, Gibson TJ: PairWise and SearchWise: finding the optimal alignment in a simultaneous comparison of a protein profile against all DNA translation frames. *Nucleic Acids Res* 1996, **24**(14):2730-2739.

116. Yang Z: Phylogenetic analysis by maximum likelihood (PAML), version 2. University College London, England 1999.

117. Goldman NAZY: A codon-based model of nucleotide substitution for protein-coding DNA sequences. *Mol Biol Evol* 1994, **11**(5):725-736.

118. Chase MW: Monocot relationships: an overview. *American Journal of Botany* 2004, **91**:1645-1655.

119. The Expressed Sequence Tags Database [http://www.ncbi.nlm.nih.gov/dbEST/]

120. Yuan Q, Ouyang S, Wang A, Zhu W, Maiti R, Lin H, Hamilton J, Haas B, Sultana R, Cheung F, et al.: The institute for genomic research Osal rice genome annotation database. *Plant Physiol* 2005, **138**(1):18-26.

Publish with **BioMed Central** and every scientist can read your work free of charge

"BioMed Central will be the most significant development for disseminating the results of biomedical research in our lifetime."

Sir Paul Nurse, Cancer Research UK

Your research papers will be:

* available free of charge to the entire biomedical community

* peer reviewed and published immediately upon acceptance

* cited in PubMed and archived on PubMed Central

* yours — you keep the copyright

Submit your manuscript here:
http://www.biomedcentral.com/info/publishing_adv.asp

<figure><img src="image_16.png" /></figure>

BioMedcentral