Chloroplast DNA evolution and phylogenetic relationships in Lycopersicon

The First Complete Chloroplast Genome of Spider Flower (Cleome houtteana) Providing a Genetic Resource for Understanding Cleomaceae Evolution

In the present study, the sequencing and analysis of the complete chloroplast genome of Cleome houtteana and its comparison with related species in the Cleomaceae family were carried out. The genome spans 157,714 base pairs (bp) and follows the typical chloroplast structure, consisting of a large single-copy (LSC) region (87,506 bp), a small single-copy (SSC) region (18,598 bp), and two inverted repeats (IRs) (25,805 bp each). We identified a total of 129 genes, including 84 protein-coding genes, 8 ribosomal RNA (rRNA) genes, and 37 transfer RNA (tRNA) genes. Our analysis of simple sequence repeats (SSRs) and repetitive elements revealed 91 SSRs, with a high number of A/T-rich mononucleotide repeats, which are common in chloroplast genomes. We also observed forward, palindromic, and tandem repeats, which are known to play roles in genome stability and evolution. When comparing C. houtteana with its relatives, we identified several highly variable regions, including ycf1, ycf2, and trnH-psbA, marking them as propitious molecular markers for the identification of species as well as phylogenetic studies. We examined the inverted repeat (IR) boundaries and found minor shifts in comparison to the other species, particularly in the ycf1 gene region, which is a known hotspot for evolutionary changes. Additionally, our analysis of selective pressures (Ka/Ks ratios) showed that most genes are under strong purifying selection, preserving their essential functions. A sliding window analysis of nucleotide diversity (Pi) identified several regions with high variability, such as trnH-psbA, ycf1, ndhI-ndhG, and trnL-ndhF, highlighting their potential for use in evolutionary and population studies. Finally, our phylogenetic analysis, using complete chloroplast genomes from species within Cleomaceae, Brassicaceae, and Capparaceae, confirmed that C. houtteana belongs within the Cleomaceae family. It showed a close evolutionary relationship with Tarenaya hassleriana and Sieruela rutidosperma, supporting previous taxonomic classifications. The findings from the current research offer invaluable insights regarding genomic structure, evolutionary adaptations, and phylogenetic relationships of C. houtteana, providing a foundation for future research on species evolution, taxonomy, and conservation within the Cleomaceae family.

GeneMiner: A tool for extracting phylogenetic markers from next-generation sequencing data

The advancement of next-generation sequencing (NGS) technologies has been revolutionary for the field of evolutionary biology. This technology has led to an abundance of available genomes and transcriptomes for researchers to mine. Specifically, researchers can mine for various types of molecular markers that are vital for phylogenetic, evolutionary and ecological studies. Numerous tools have been developed to extract these molecular markers from NGS data. However, due to an insufficient number of well-annotated reference genomes for non-model organisms, it remains challenging to obtain these markers accurately and efficiently. Here, we present GeneMiner, an improved and expanded version of our previous tool, Easy353. GeneMiner combines the reference-guided de Bruijn graph assembly with seed self-discovery and greedy extension. Additionally, it includes a verification step using a parameter-bootstrap method to reduce the pitfalls associated with using a relatively distant reference. Our results, using both experimental and simulation data, showed GeneMiner can accurately acquire phylogenetic molecular markers for plants using transcriptomic, genomic and other NGS data. GeneMiner is designed to be user-friendly, fast and memory-efficient. Further, it is compatible with Linux, Windows and macOS. All source codes are publicly available on GitHub (https://github.com/sculab/GeneMiner) and Gitee (https://gitee.com/sculab/GeneMiner) for easy accessibility and transparency.

Comparative and phylogenetic analysis of the complete chloroplast genomes of ten Pittosporum species from East Asia

Pittosporum (Pittosporaceae) is famous as the ornamental and medical values, which is distributed tropical and subtropical regions of Eastern Hemisphere. The few phylogenetic studies have included samples from the Pacific Island, but the phylogenetic relationships of Asian species has not been studied. Here, the complete chloroplast (cp) genomes of ten Pittosporum species from East Asia were first sequenced and compared with those of the published species of this genus. Our results indicated that cp genomes of these species had a typical and conserved quadripartite structure. 131 genes were identical in order and orientation and no changes of inverted repeat (IR) occurred. However, the comparative analysis of cp genomes suggested that sequence divergence mainly appeared in non-coding or intergenic regions, in which several divergence hotspots were identified. By contrast, protein-coding genes showed the lowest variance under strong purifying selection. Phylogenetic analysis based on the cp genome sequences showed that the tested Pittosporum species were clustered into two major clades, in which the Asian species formed Clade I and the remaining species from Australia and New Zealand formed Clade II with high support values, which was consistent with the results of ITS data with low support values. These results suggested that cp genome is a robust phylogenetic indicator for deep nodes in the phylogeny of Pittosporum. Meanwhile, these results will provide the valuable information to better understand the phylogeny and biogeography of Pittosporum.

First Record of Comparative Plastid Genome Analysis and Phylogenetic Relationships among Corylopsis Siebold & Zucc. (Hamamelidaceae)

Corylopsis Siebold & Zucc. (Hamamelidaceae) is widely used as a horticultural plant and comprises approximately 25 species in East Asia. Molecular research is essential to distinguish Corylopsis species, which are morphologically similar. Molecular research has been conducted using a small number of genes but not in Corylopsis. Plastid genomes of Corylopsis species (Corylopsis gotoana, Corylopsis pauciflora, and Corylopsis sinensis) were sequenced using next-generation sequencing techniques. Repeats and nucleotide diversity that could be used as DNA markers were also investigated. A phylogenetic investigation was carried out using 79 protein-coding genes to infer the evolutionary relationships within the genus Corylopsis. By including new plastomes, the overall plastid genome structure of Corylopsis was similar. Simple sequence repeats of 73-106 SSRs were identified in the protein-coding genes of the plastid genomes, and 33-40 long repeat sequences were identified in the plastomes. The Pi value of the rpl33_rps18 region, an intergenic spacer, was the highest. Phylogenetic analysis demonstrated that Corylopsis is a monophyletic group and Loropetalum is closely related to Corylopsis. C. pauciflora, C. gotoana, and C. spicata formed a clade distributed in Japan, whereas C. sinensis, C. glandulifera, and C. velutina formed a clade that was distributed in China.

Phylogenomic and syntenic data demonstrate complex evolutionary processes in early radiation of the rosids

Some of the most vexing problems of deep level relationship that remain in angiosperms involve the superrosids. The superrosid clade contains a quarter of all angiosperm species, with 18 orders in three subclades (Vitales, Saxifragales and core rosids) exhibiting remarkable morphological and ecological diversity. To help resolve deep-level relationships, we constructed a high-quality chromosome-level genome assembly for Tiarella polyphylla (Saxifragaceae) thus providing broader genomic representation of Saxifragales. Whole genome microsynteny analysis of superrosids showed that Saxifragales shared more synteny clusters with core rosids than Vitales, further supporting Saxifragales as more closely related with core rosids. To resolve the ordinal phylogeny of superrosids, we screened 122 single copy nuclear genes from genomes of 36 species, representing all 18 superrosid orders. Vitales were recovered as sister to all other superrosids (Saxifragales + core rosids). Our data suggest dramatic differences in relationships compared to earlier studies within core rosids. Fabids should be restricted to the nitrogen-fixing clade, while Picramniales, the Celastrales-Malpighiales (CM) clade, Huerteales, Oxalidales, Sapindales, Malvales and Brassicales formed an "expanded" malvid clade. The Celastrales-Oxalidales-Malpighiales (COM) clade (sensu APG IV) was not monophyletic. Crossosomatales, Geraniales, Myrtales and Zygophyllales did not belong to either of our well-supported malvids or fabids. There is strong discordance between nuclear and plastid phylogenetic hypotheses for superrosid relationships; we show that this is best explained by a combination of incomplete lineage sorting and ancient reticulation.

PhyloHerb: A high-throughput phylogenomic pipeline for processing genome skimming data

Premise

The application of high-throughput sequencing, especially to herbarium specimens, is rapidly accelerating biodiversity research. Low-coverage sequencing of total genomic DNA (genome skimming) is particularly promising and can simultaneously recover the plastid, mitochondrial, and nuclear ribosomal regions across hundreds of species. Here, we introduce PhyloHerb, a bioinformatic pipeline to efficiently assemble phylogenomic data sets derived from genome skimming.

Methods and Results

PhyloHerb uses either a built-in database or user-specified references to extract orthologous sequences from all three genomes using a BLAST search. It outputs FASTA files and offers a suite of utility functions to assist with alignment, partitioning, concatenation, and phylogeny inference. The program is freely available at https://github.com/lmcai/PhyloHerb/.

Conclusions

We demonstrate that PhyloHerb can accurately identify genes using a published data set from Clusiaceae. We also show via simulations that our approach is effective for highly fragmented assemblies from herbarium specimens and is scalable to thousands of species.

Complete pan-plastome sequences enable high resolution phylogenetic classification of sugar beet and closely related crop wild relatives

Background

As the major source of sugar in moderate climates, sugar-producing beets (Beta vulgaris subsp. vulgaris) have a high economic value. However, the low genetic diversity within cultivated beets requires introduction of new traits, for example to increase their tolerance and resistance attributes – traits that often reside in the crop wild relatives. For this, genetic information of wild beet relatives and their phylogenetic placements to each other are crucial. To answer this need, we sequenced and assembled the complete plastome sequences from a broad species spectrum across the beet genera Beta and Patellifolia, both embedded in the Betoideae (order Caryophyllales). This pan-plastome dataset was then used to determine the wild beet phylogeny in high-resolution.

Results

We sequenced the plastomes of 18 closely related accessions representing 11 species of the Betoideae subfamily and provided high-quality plastome assemblies which represent an important resource for further studies of beet wild relatives and the diverse plant order Caryophyllales. Their assembly sizes range from 149,723 bp (Beta vulgaris subsp. vulgaris) to 152,816 bp (Beta nana), with most variability in the intergenic sequences. Combining plastome-derived phylogenies with read-based treatments based on mitochondrial information, we were able to suggest a unified and highly confident phylogenetic placement of the investigated Betoideae species.

Our results show that the genus Beta can be divided into the two clearly separated sections Beta and Corollinae. Our analysis confirms the affiliation of B. nana with the other Corollinae species, and we argue against a separate placement in the Nanae section. Within the Patellifolia genus, the two diploid species Patellifolia procumbens and Patellifolia webbiana are, regarding the plastome sequences, genetically more similar to each other than to the tetraploid Patellifolia patellaris. Nevertheless, all three Patellifolia species are clearly separated.

Conclusion

In conclusion, our wild beet plastome assemblies represent a new resource to understand the molecular base of the beet germplasm. Despite large differences on the phenotypic level, our pan-plastome dataset is highly conserved. For the first time in beets, our whole plastome sequences overcome the low sequence variation in individual genes and provide the molecular backbone for highly resolved beet phylogenomics. Hence, our plastome sequencing strategy can also guide genomic approaches to unravel other closely related taxa.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08336-8.

The Generalized Robinson-Foulds Distance for Phylogenetic Trees

The Robinson-Foulds (RF) distance, one of the most widely used metrics for comparing phylogenetic trees, has the advantage of being intuitive, with a natural interpretation in terms of common splits, and it can be computed in linear time, but it has a very low resolution, and it may become trivial for phylogenetic trees with overlapping taxa, that is, phylogenetic trees that share some but not all of their leaf labels. In this article, we study the properties of the Generalized Robinson-Foulds (GRF) distance, a recently proposed metric for comparing any structures that can be described by multisets of multisets of labels, when applied to rooted phylogenetic trees with overlapping taxa, which are described by sets of clusters, that is, by sets of sets of labels. We show that the GRF distance has a very high resolution, it can also be computed in linear time, and it is not (uniformly) equivalent to the RF distance.

Whole genomes: the holy grail. A commentary on: 'Molecular phylogenomics of the tribe Shoreeae (Dipterocarpaceae) using whole plastidgenomes'

This article comments on:

Jacqueline Heckenhauer, Ovidiu Paun, Mark W. Chase, Peter S. Ashton, A.S. Kamariah and Rosabelle Samuel. 2019. Molecular phylogenomics of the tribe Shoreeae (Dipterocarpaceae) using whole plastid genomes. Annals of Botany 123(5): 857–865.

Analyses of 202 plastid genomes elucidate the phylogeny of Solanum section Petota

Our paper analyzes full plastid DNA sequence data of 202 wild and cultivated diploid potatoes, Solanum section Petota, to explore its phylogenetic utility compared to prior analyses of the same accessions using genome-wide nuclear SNPs, and plastid DNA restriction site data. The present plastid analysis discovered the same major clades as the nuclear data but with some substantial differences in topology within the clades. The considerably larger plastid and nuclear data sets add phylogenetic resolution within the prior plastid DNA restriction site data, highlight plastid/nuclear incongruence that supports hypotheses of hybridization/introgression to help explain the taxonomic difficulty in the section.

Phylogenomics of tomato chloroplasts using assembly and alignment-free method

Phylogenetics and population genetics are central disciplines in evolutionary biology. Both are based on the comparison of single DNA sequences, or a concatenation of a number of these. However, with the advent of next-generation DNA sequencing technologies, the approaches that consider large genomic data sets are of growing importance for the elucidation of evolutionary relationships among species. Among these approaches, the assembly and alignment-free methods which allow an efficient distance computation and phylogeny reconstruction are of great importance. However, it is not yet clear under what quality conditions and abundance of genomic data such methods are able to infer phylogenies accurately. In the present study we assess the method originally proposed by Fan et al. for whole genome data, in the elucidation of Tomatoes' chloroplast phylogenetics using short read sequences. We find that this assembly and alignment-free method is capable of reproducing previous results under conditions of high coverage, given that low frequency k-mers (i.e. error prone data) are effectively filtered out. Finally, we present a complete chloroplast phylogeny for the best data quality candidates of the recently published 360 tomato genomes.

A MOLECULAR REEXAMINATION OF INTROGRESSION BETWEEN HELIANTHUS ANNUUS AND H. BOLANDERI (COMPOSITAE)

Heiser (1949) hypothesized that a weedy race of Helianthus bolanderi had originated by the introgression of genes from H. annum into a serpentine race of H. bolanderi. Although Heiser's investigation of these species is frequently cited as one of the best examples of introgression in plants, definitive evidence of gene exchange is lacking (Heiser, 1973). To determine whether the weedy race of H. bolanderi actually originated via introgression, we analyzed allozyme, chloroplast-DNA (cpDNA), and nuclear-ribosomal-DNA (rDNA) variation. Evidence from enzyme electrophoresis did not support the proposed introgressive origin of weedy H. bolanderi. We detected a total of 37 low-frequency alleles distinguishing the serpentine race of H. bolanderi from H. annuus. Weedy H. bolanderi possessed only four of the 37 marker alleles. Further analysis demonstrated that serpentine H. bolanderi combined seven of the 35 alleles distinguishing H. annuus from weedy H. bolanderi, indicating that serpentine H. bolanderi shares three more marker alleles with H. annuus than does weedy H. bolanderi. These results are similar to expectations for race divergence from a single common ancestor and suggest that, if introgression occurred, the majority of marker alleles were rapidly lost following the initial hybridization event. Even more compelling evidence opposing Heiser's (1949) hypothesis, however, was from restriction-fragment analysis of cpDNA and nuclear rDNA. We detected a total of 17 cpDNA and five rDNA restriction-site mutations among the 19 populations examined. No parallel or back mutations were observed in phylogenetic trees constructed using either cpDNA or rDNA mutations, and both phylogenies were completely congruent regarding the alignment of all three taxa. In addition, the weedy race of H. bolanderi possessed a unique cpDNA, which was outside the range of variation observed among populations of either of the presumed parental species. Mean sequence divergences between the cpDNAs of weedy H. bolanderi and those of serpentine H. bolanderi and H. annuus were 0.30% and 0.35%, respectively. These estimates are comparable to sequence-divergence values observed between closely related species in other plant groups. Given the lack of parallel or convergent mutations in the cpDNA and rDNA phylogenetic trees, the complete congruence of these trees with flavonoid- and allozyme-variation patterns, and the presence of a unique and divergent chloroplast genome in the weedy race of H. bolanderi, we suggest that the weedy race of H. bolanderi was not derived recently through introgression, as hypothesized, but is relatively ancient in origin.

Cuticular hydrocarbons and sucrose esters as chemotaxonomic markers of wild and cultivated tomato species (Solanum section Lycopersicon)

The tomato (Solanum lycopersicum L.) is one of the most important vegetables worldwide. Due to the limited genetic variability, wild related species are considered as potential gene pool for breeding cultivated plants with enriched genetic basis. Taxonomic relations between tomato species at the level of single groups and taxa still remain, however, not fully resolved. Hence, in addition to already reported classification based on the morphology of the plants and molecular markers, we proposed chemotaxonomic approach to unveil some aspects of tomato taxonomy. Cuticular hydrocarbons and surface sucrose esters (SEs) were used as chemotaxonomic markers. Classification based on the cuticular hydrocarbon profile was in good agreement with other taxonomic studies as long as between-species differences were taken into account. Clear separation of the common tomato and closely related species from the majority of S. pennellii accessions was obtained. In the same time, however, S. pennellii revealed broad variation: based on the results, three highly distinct types of these plants were proposed, among them one type was very similar to cultivated tomato and its relatives. Addition of SEs profiles to the dataset did not impair the classification, but clarified the position of S. pennellii. The results suggest possible hybrid origin of some of S. pennellii and wild S. lycopersicum accessions, and the approach proposed has a potential to identify such hybrid plant lines.

Phylogenomics Reveals Three Sources of Adaptive Variation during a Rapid Radiation

Speciation events often occur in rapid bursts of diversification, but the ecological and genetic factors that promote these radiations are still much debated. Using whole transcriptomes from all 13 species in the ecologically and reproductively diverse wild tomato clade (Solanum sect. Lycopersicon), we infer the species phylogeny and patterns of genetic diversity in this group. Despite widespread phylogenetic discordance due to the sorting of ancestral variation, we date the origin of this radiation to approximately 2.5 million years ago and find evidence for at least three sources of adaptive genetic variation that fuel diversification. First, we detect introgression both historically between early-branching lineages and recently between individual populations, at specific loci whose functions indicate likely adaptive benefits. Second, we find evidence of lineage-specific de novo evolution for many genes, including loci involved in the production of red fruit color. Finally, using a "PhyloGWAS" approach, we detect environment-specific sorting of ancestral variation among populations that come from different species but share common environmental conditions. Estimated across the whole clade, small but substantial and approximately equal fractions of the euchromatic portion of the genome are inferred to contribute to each of these three sources of adaptive genetic variation. These results indicate that multiple genetic sources can promote rapid diversification and speciation in response to new ecological opportunity, in agreement with our emerging phylogenomic understanding of the complexity of both ancient and recent species radiations.

Next-generation sampling: Pairing genomics with herbarium specimens provides species-level signal in Solidago (Asteraceae)

PREMISE OF THE STUDY

The ability to conduct species delimitation and phylogeny reconstruction with genomic data sets obtained exclusively from herbarium specimens would rapidly enhance our knowledge of large, taxonomically contentious plant genera. In this study, the utility of genotyping by sequencing is assessed in the notoriously difficult genus Solidago (Asteraceae) by attempting to obtain an informative single-nucleotide polymorphism data set from a set of specimens collected between 1970 and 2010.

METHODS

Reduced representation libraries were prepared and Illumina-sequenced from 95 Solidago herbarium specimen DNAs, and resulting reads were processed with the nonreference Universal Network-Enabled Analysis Kit (UNEAK) pipeline. Multidimensional clustering was used to assess the correspondence between genetic groups and morphologically defined species.

RESULTS

Library construction and sequencing were successful in 93 of 95 samples. The UNEAK pipeline identified 8470 single-nucleotide polymorphisms, and a filtered data set was analyzed for each of three Solidago subsections. Although results varied, clustering identified genomic groups that often corresponded to currently recognized species or groups of closely related species.

DISCUSSION

These results suggest that genotyping by sequencing is broadly applicable to DNAs obtained from herbarium specimens. The data obtained and their biological signal suggest that pairing genomics with large-scale herbarium sampling is a promising strategy in species-rich plant groups.

Next-generation sampling: Pairing genomics with herbarium specimens provides species-level signal in Solidago (Asteraceae)

PREMISE OF THE STUDY

The ability to conduct species delimitation and phylogeny reconstruction with genomic data sets obtained exclusively from herbarium specimens would rapidly enhance our knowledge of large, taxonomically contentious plant genera. In this study, the utility of genotyping by sequencing is assessed in the notoriously difficult genus Solidago (Asteraceae) by attempting to obtain an informative single-nucleotide polymorphism data set from a set of specimens collected between 1970 and 2010.

METHODS

Reduced representation libraries were prepared and Illumina-sequenced from 95 Solidago herbarium specimen DNAs, and resulting reads were processed with the nonreference Universal Network-Enabled Analysis Kit (UNEAK) pipeline. Multidimensional clustering was used to assess the correspondence between genetic groups and morphologically defined species.

RESULTS

Library construction and sequencing were successful in 93 of 95 samples. The UNEAK pipeline identified 8470 single-nucleotide polymorphisms, and a filtered data set was analyzed for each of three Solidago subsections. Although results varied, clustering identified genomic groups that often corresponded to currently recognized species or groups of closely related species.

DISCUSSION

These results suggest that genotyping by sequencing is broadly applicable to DNAs obtained from herbarium specimens. The data obtained and their biological signal suggest that pairing genomics with large-scale herbarium sampling is a promising strategy in species-rich plant groups.

Genetic structure of the four wild tomato species in the Solanum peruvianum s.l. species complex

The most diverse wild tomato species Solanum peruvianum sensu lato (s.l.) has been reclassified into four separate species: Solanum peruvianum sensu stricto (s.s.), Solanum corneliomuelleri, Solanum huaylasense, and Solanum arcanum. However, reproductive barriers among the species are incomplete and this can lead to discrepancies regarding genetic identity of germplasm. We used genotyping by sequencing (GBS) of S. peruvianum s.l., Solanum neorickii, and Solanum chmielewskii to develop tens of thousands of mapped single nucleotide polymorphisms (SNPs) to analyze genetic relationships within and among species. The data set was condensed to 14,043 SNPs with no missing data across 46 sampled plants. Origins of accessions were mapped using geographical information systems (GIS). Isolation by distance, pairwise genetic distances, and number of clusters were estimated using population genetics approaches. Isolation by distance was strongly supported, especially between interspecific pairs. Eriopersicon (S. peruvianum s.s., S. corneliomuelleri, S. huaylasense) and Arcanum (S. arcanum, S. neorickii, S. chmielewskii) species groups were genetically distinct, except for S. huaylasense which showed 50% membership proportions in each group. Solanum peruvianum and S. corneliomuelleri were not significantly differentiated from each other. Many thousands of SNP markers were identified that could potentially be used to distinguish pairs of species, including S. peruvianum versus S. corneliomuelleri, if they are verified on larger numbers of samples. Diagnostic markers will be valuable for delimiting morphologically similar and interfertile species in germplasm management. Approximately 12% of the SNPs rejected a genome-wide test of selective neutrality based on differentiation among species of S. peruvianum s.l. These are candidates for more comprehensive studies of microevolutionary processes within this species complex.

Tomato breeding in the genomics era: insights from a SNP array

BACKGROUND

The major bottle neck in genetic and linkage studies in tomato has been the lack of a sufficient number of molecular markers. This has radically changed with the application of next generation sequencing and high throughput genotyping. A set of 6000 SNPs was identified and 5528 of them were used to evaluate tomato germplasm at the level of species, varieties and segregating populations.

RESULTS

From the 5528 SNPs, 1980 originated from 454-sequencing, 3495 from Illumina Solexa sequencing and 53 were additional known markers. Genotyping different tomato samples allowed the evaluation of the level of heterozygosity and introgressions among commercial varieties. Cherry tomatoes were especially different from round/beefs in chromosomes 4, 5 and 12. We were able to identify a set of 750 unique markers distinguishing S. lycopersicum 'Moneymaker' from all its distantly related wild relatives. Clustering and neighbour joining analysis among varieties and species showed expected grouping patterns, with S. pimpinellifolium as the most closely related to commercial tomatoes earlier results.

CONCLUSIONS

Our results show that a SNP search in only a few breeding lines already provides generally applicable markers in tomato and its wild relatives. It also shows that the Illumina bead array generated data are highly reproducible. Our SNPs can roughly be divided in two categories: SNPs of which both forms are present in the wild relatives and in domesticated tomatoes (originating from common ancestors) and SNPs unique for the domesticated tomato (originating from after the domestication event). The SNPs can be used for genotyping, identification of varieties, comparison of genetic and physical linkage maps and to confirm (phylogenetic) relations. In the SNPs used for the array there is hardly any overlap with the SolCAP array and it is strongly recommended to combine both SNP sets and to select a core collection of robust SNPs completely covering the entire tomato genome.

Variation revealed by SNP genotyping and morphology provides insight into the origin of the tomato

Tomato, Solanum lycopersicum, is divided into two widely distributed varieties: the cultivated S. lycopersicum var. lycopersicum, and the weedy S. lycopersicum var. cerasiforme. Solanum pimpinellifolium is the most closely related wild species of tomato.

The roles of S. pimpinellifolium and S. l. cerasiforme during the domestication of tomato are still under debate. Some authors consider S. l. cerasiforme to be the ancestor, whereas others think that S. l. cerasiforme is an admixture of S. pimpinellifolium and the cultivated S. l. lycopersicum. It is also not clear whether the domestication occurred in the Andean region or in Mesoamerica. We characterized 272 accessions (63 S. pimpinellifolium, 106 S. l. cerasiforme, 95 S. l. lycopersicum and 8 derived from hybridization processes) were morphologically and genetically using the SolCap platform (7,414 SNPs). The two species were distinguished in a PCA analysis and displayed a rich geographic structure. Solanum lycopersicum var. cerasiforme and S. l. lycopersicum were also differentiated in the PCA and Structure analyses, which supports maintaining them as different varieties. Solanum pimpinellifolium and the Andean S. l. cerasiforme were more diverse than the non-Andean S. lycopersicum. Solanum lycopersicum var. cerasiforme was morphologically and molecularly intermediate between S. pimpinellifolium and tomato. Solanum lycopersicum var. cerasiforme, with the exception of several Ecuadorian and Mexican accessions, is composed of the products of admixture processes according to the Structure analysis. The non-admixtured S. l. cerasiforme might be similar to the ancestral cultivars from which the cultivated tomato originated, and presents remarkable morphological diversity, including fruits of up to 6 cm in diameter. The data obtained would fit a model in which a pre-domestication took place in the Andean region, with the domestication being completed in Mesoamerica. Subsequently, the Spaniards took plants from Mesoamerica to Spain and from there they were exported to the rest of the world.

Complete plastome sequences of Equisetum arvense and Isoetes flaccida: implications for phylogeny and plastid genome evolution of early land plant lineages

Background

Despite considerable progress in our understanding of land plant phylogeny, several nodes in the green tree of life remain poorly resolved. Furthermore, the bulk of currently available data come from only a subset of major land plant clades. Here we examine early land plant evolution using complete plastome sequences including two previously unexamined and phylogenetically critical lineages. To better understand the evolution of land plants and their plastomes, we examined aligned nucleotide sequences, indels, gene and nucleotide composition, inversions, and gene order at the boundaries of the inverted repeats.

Results

We present the plastome sequences of Equisetum arvense, a horsetail, and of Isoetes flaccida, a heterosporous lycophyte. Phylogenetic analysis of aligned nucleotides from 49 plastome genes from 43 taxa supported monophyly for the following clades: embryophytes (land plants), lycophytes, monilophytes (leptosporangiate ferns + Angiopteris evecta + Psilotum nudum + Equisetum arvense), and seed plants. Resolution among the four monilophyte lineages remained moderate, although nucleotide analyses suggested that P. nudum and E. arvense form a clade sister to A. evecta + leptosporangiate ferns. Results from phylogenetic analyses of nucleotides were consistent with the distribution of plastome gene rearrangements and with analysis of sequence gaps resulting from insertions and deletions (indels). We found one new indel and an inversion of a block of genes that unites the monilophytes.

Conclusions

Monophyly of monilophytes has been disputed on the basis of morphological and fossil evidence. In the context of a broad sampling of land plant data we find several new pieces of evidence for monilophyte monophyly. Results from this study demonstrate resolution among the four monilophytes lineages, albeit with moderate support; we posit a clade consisting of Equisetaceae and Psilotaceae that is sister to the "true ferns," including Marattiaceae.

Chloroplast DNA evidence for the origin of the genus Heterogaura from a species of Clarkia (Onagraceae)

Restriction-site variation in chloroplast DNA was examined in the morphologically distinct and monotypic genus Heterogaura and the related speciose genus Clarkia (Onagraceae), both native to California. Of the 605 restriction sites surveyed, a total of 119 mutations were identified. Of these, 55 were shared by at least two species and were used to construct a most parsimonious phylogenetic tree. This analysis, as well as one based on a distance metric, provided evidence that Heterogaura and Clarkia dudleyana, a member of a phylogenetically advanced section, share a more recent common ancestor than either does with any other species. The two species are more closely related than nearly all paris of Clarkia tested. The origin of Heterogaura from within another genus raises important questions about the adequacy of morphological data and suggests that the relationships of other well-known monotypic plant genera should be reinvestigated.

Chloroplast DNA evolution and systematics of Phanerophlebia (Dryopteridaceae) and related fern genera

Restriction site variation in chloroplast DNA was examined in the neotropical fern genus Phanerophlebia and in selected species of the related Asiatic genus Cyrtomium and the cosmopolitan progenitor of these two, Polystichum. A total of 103 restriction site mutations was identified; these were used to construct phylogenetic networks and trees based on Wagner and Dollo parsimony and Fitch-Margoliash distance algorithms. The analyses provided evidence that Phanerophlebia did not arise from Cyrtomium. Both genera are convergent descendants from different progenitor groups in Polystichum, and Asiatic Cyrtomium is more closely related to temperate New World Polystichum than it is to neotropical Phanerophlebia. Reticulate venation, previously considered an important taxonomic character for infrageneric classification in Phanerophlebia, most likely evolved independently twice within the genus. Diploid maternal progenitors are suggested for two of four tetraploid species analyzed, and species-level distinctions for two closely related taxa of Phanerophlebia are questioned.

New World tetraploid cottons contain Old World cytoplasm

New World tetraploid cottons (Gossypium spp.) originated through hybridization of ancestral diploid species that presently have allopatric ranges in Asia-Africa (the A genome) and the New World tropics and subtropics (the D genome). Despite intensive study, the identity of the parental diploids and the antiquity of polyploidization remain unresolved. In this study, variation in the maternally inherited chloroplast genome was assessed among species representing both of the parental genomes and the tetraploids. Approximately 560 restriction sites were assayed in each accession, representing sequence information for about 3200 nucleotides. The resulting maternal phylogeny has no convergent restriction site mutations and demonstrates that the cytoplasm donor for all tetraploid species was an A genome diploid with a chloroplast genome that is similar to Gossypium arboreum and Gossypium herbaceum. No mutational differences were detected between these two species, and few mutations distinguish the chloroplast genomes of A genome diploids from those of tetraploid taxa. In contrast to expectations based on extensive taxonomic, geographic, and genetic diversity, a surprisingly low level of sequence divergence has accumulated subsequent to polyploidization. Chloroplast genomes of tetraploid species are distinguished from each other by between one and six apparent point mutations. The data suggest that tetraploid cotton originated relatively recently, perhaps within the last 1-2 million years, with subsequent rapid evolution and diversification throughout the New World tropics.

Ecological and geographic modes of species divergence in wild tomatoes

Understanding the role of geography and ecology in species divergence is central to the study of evolutionary diversification. We used climatic, geographic, and biological data from nine wild Andean tomato species to describe each species' ecological niche and to evaluate the likely ecological and geographical modes of speciation in this clade. Using data from >1000 wild accessions and publicly available data derived from geographic information systems for various environmental variables, we found most species pairs were significantly differentiated for one or more environmental variables. By comparing species' predicted niches generated by species distribution modeling (SDM), we found significant niche differentiation among three of four sister-species pairs, suggesting ecological divergence is consistently associated with recent divergence. In comparison, based on age-range correlation (ARC) analysis, there was no evidence for a predominant geographical (allopatric vs. sympatric) context for speciation in this group. Overall, our results suggest an important role for environmentally mediated differentiation, rather than simply geographical isolation, in species divergence.

A new genetic linkage map of tomato based on a Solanum lycopersicum x S. pimpinellifolium RIL population displaying locations of candidate pathogen response genes

The narrow genetic base of the cultivated tomato, Solanum lycopersicum L., necessitates introgression of new variation from related species. Wild tomato species represent a rich source of useful genes and traits. Exploitation of genetic variation within wild species can be facilitated by the use of molecular markers and genetic maps. Recently we identified an accession (LA2093) within the red-fruited wild tomato species Solanum pimpinellifolium L. with exceptionally desirable characteristics, including disease resistance, abiotic stress tolerance, and high fruit lycopene content. To facilitate genetic characterization of such traits and their exploitation in tomato crop improvement, we developed a new recombinant inbred line (RIL) population from a cross between LA2093 and an advanced tomato breeding line (NCEBR-1). Furthermore, we constructed a medium-density molecular linkage map of this population using 294 polymorphic markers, including standard RFLPs, EST sequences (used as RFLP probes), CAPS, and SSRs. The map spanned 1091 cM of the tomato genome with an average marker spacing of 3.7 cM. A majority of the EST sequences, which were mainly chosen based on the putative role of their unigenes in disease resistance, defense-related response, or fruit quality, were mapped onto the tomato chromosomes for the first time. Co-localizations of relevant EST sequences with known disease resistance genes in tomato were also examined. This map will facilitate identification, genetic exploitation, and positional cloning of important genes or quantitative trait loci in LA2093. It also will allow the elucidation of the molecular mechanism(s) underlying important traits segregating in the RIL population. The map may further facilitate characterization and exploitation of genetic variation in other S. pimpinellifolium accessions as well as in modern cultivars of tomato.

Do potatoes and tomatoes have a single evolutionary history, and what proportion of the genome supports this history?

BACKGROUND

Phylogenies reconstructed with only one or a few independently inherited loci may be unresolved or incongruent due to taxon and gene sampling, horizontal gene transfer, or differential selection and lineage sorting at individual loci. In an effort to remedy this situation, we examined the utility of conserved orthologous set (COSII) nuclear loci to elucidate the phylogenetic relationships among 29 diploid Solanum species in the sister clades that include tomato and potato, and in Datura inoxia as a far outgroup. We screened 40 COSII markers with intron content over 60% that are mapped in different chromosomes; selected a subset of 19 by the presence of single band amplification of size mostly between 600 and 1200 bp; sequenced these 19 COSII markers, and performed phylogenetic analyses with individual and concatenated datasets. The present study attempts to provide a fully resolved phylogeny among the main clades in potato and tomato that can help to identify the appropriate markers for future studies using additional species.

RESULTS

Among potatoes, when total evidence is invoked, one single predominant history is highlighted with complete resolution within and among the three main clades. It also supports the hypothesis of the North and Central American B-genome origin of the tuber-bearing members of Solanum sect. Petota and shows a clear division between A genomes in clades 3 and 4, and B genomes in clade 1+2. On the other hand, when a prior agreement approach is invoked other potato evolutionary histories are revealed but with less support. These alternative histories could be explained by past hybridization, or fast rates of speciation. In the case of tomato, the analyses with all sequence data completely resolved 19 of 21 clades, for the first time revealed the monophyly of five clades, and gave further support for the recent segregation of new species from the former Solanum peruvianum. Concordance analyses revealed and summarized the extensive discordance among COSII markers. Some potential reasons for discordance could be methodological, to include systematic errors due to using a wrong model of sequence evolution, coupled with long branches, or mixtures of branch lengths within COSII, or undetected paralogy or alignment bias. Other reasons could be biological processes such as hybridization or lineage sorting.

CONCLUSION

This study confirms and quantifies the utility of using DNA sequences from different parts of the genome in phylogenetic studies to avoid possible bias in the sampling. It shows that 11-18 loci are enough to get the dominant history in this group of Solanum, but more loci would be needed to discern the distribution of gene genealogies in more depth, and thus detect which mechanism most likely shaped the discordance.

Do plant parts compete for resources? An evolutionary viewpoint

Simultaneously growing sinks are thought to compete for plant resources. Negative correlations, for example between grain number and stem mass in cereals, indeed resemble competition; but is the notion of intra-plant competition evolutionarily justified? Here we review intra-plant competition in light of two aspects of evolutionary biology: (a) major transitions that led to the reorganization of evolutionary individuals (e.g. isolated DNA molecules and independent cells) into new units of adaptation (e.g. chromosomes and multicellular organisms) with associated constraints to intra-individual conflict; and (b) genomic conflicts within individual plants with implications for resource allocation. Against this background, we look at apparent competition among genetically identical plant parts, and conclude that plants might use competition-like mechanisms to allocate resources, but only to the extent that these proximate mechanisms enhance overall plant fitness. In dealing with apparent competition among genetically different plant structures, we emphasize developing seeds attached to the same maternal plant, and the determination of yield components in annual crops. We propose that competition-like mechanisms among genetically different plant parts have been strongly shaped by the evolution of genomic conflict between parent and offspring, between female and male parents, and among siblings. By defining the number and potential size of grain simultaneously and before fertilization, a strong maternal control of resource allocation is exerted that favours uniform offspring size and partially counteracts genomic conflict.

[Combining phylogenetic information: concept, methodology, and challenges]

The DNA sequences, morphological and other homologous characters can be used to infer the origins and histories of biological taxa. Combining all the phylogenetic information available can produce more inclusive phylogenies, improve our understanding of living organisms, and enable biologists to prompt and test hypotheses on a larger scale and with stronger statistical power. In this article, the concept of combining phylogenetic information and its comparison with traditional analysis were reviewed. The most popular approaches of supertree and supermatrix were discussed in detail, and novel ways were presented. Although the combining analysis is facing rigid challenges from data and foundation, it is currently the only approach for realization of the Tree(Net) of Life, and its development will definitely expand our knowledge of evolution on the earth and contribute to the progress of evolutionary related disciplines.

Chloroplast DNA diversity in wild and cultivated barley: implications for genetic conservation

Nine diverse lines of cultivated barley (Hordeum vulgare) and 11 lines of its wild progenitor (H. spontaneum) were assayed for variation in their chloroplast DNA by digestion with ten restriction endonucleases. The cultivated lines exhibited a single cpDNA polymorphism, whereas the wild material exhibited five. The significantly lower level of diversity among the cultivated lines was unexpected because both cultivated and wild lines had been selected for comparable levels of diversity for nuclear encoded isozyme loci. These results suggest that the level of cytoplasmic diversity was markedly restricted during the domestication of cultivated barley.

A Solanum lycopersicum x Solanum pimpinellifolium linkage map of tomato displaying genomic locations of R-genes, RGAs, and candidate resistance/defense-response ESTs

We have identified an accession (LA2093) within the tomato wild species Solanum pimpinellifolium with many desirable characteristics, including biotic and abiotic stress tolerance and good fruit quality. To utilize the full genetic potential of LA2093 in tomato breeding, we have developed a linkage map based on an F(2) population of a cross between LA2093 and a tomato breeding line, using 115 RFLP, 94 EST, and 41 RGA markers. The map spanned 1002.4 cM of the 12 tomato chromosomes with an average marker distance of 4.0 cM. The length of the map and linear order of the markers were in good agreement with the published maps of tomato. The ESTs were chosen based on their sequence similarities with known resistance or defense-response genes, signal-transduction factors, transcriptional regulators, and genes encoding pathogenesis-related proteins. Locations of several ESTs and RGAs coincided with locations of several known tomato resistance genes and quantitative resistance loci (QRLs), suggesting that candidate-gene approach may be effective in identifying and mapping new R genes. This map will be useful for marker-assisted exploitation of desirable traits in LA2093 and other S. pimpinellifolium accessions, and possibly for utilization of genetic variation within S. lycopersicum.

Poriferan mtDNA and animal phylogeny based on mitochondrial gene arrangements

Phylogenetic relationships among the metazoan phyla are the subject of an ongoing controversy. Analysis of mitochondrial gene arrangements is a powerful tool to investigate these relationships; however, its previous application outside of individual animal phyla has been hampered by the lack of informative out-group data. To address this shortcoming, we determined complete mitochondrial DNA sequences for the demosponges Geodia neptuni and Tethya actinia, two representatives of the most basal animal phylum, the Porifera. With sponges as an outgroup, we investigated phylogenetic relationships of nine bilaterian phyla using both breakpoint analysis of global mitochondrial gene arrangements and maximum parsimony analysis of mitochondrial gene adjacencies. Our results provide strong support for a group that includes protostome (but not deuterostome) coelomate, pseudocoelomate, and acoelomate animals, thus clearly rejecting the Coelomata hypothesis. Two other groups of bilaterian animals, Lophotrochozoa and Ambulacraria, are also supported by our analyses. However, due to the remarkable stability of mitochondrial gene arrangements in Deuterostomia and the Ecdysozoa, conclusions on their evolutionary history cannot be drawn.

High-resolution mapping and functional analysis of se2.1: a major stigma exsertion quantitative trait locus associated with the evolution from allogamy to autogamy in the genus Lycopersicon

The degree to which stigmas are exserted above the stamen in flowers is a key determinant of cross-pollination (and hence allogamy) in many plant species. Most species in the genus Lycopersicon are obligate or facultative outcrossers and bear flowers with highly exserted stigmas. In contrast, the cultivated tomato (Lycopersicon esculentum) bears flowers with flush or inserted stigmas promoting self-fertilization. It has been observed that a major QTL, se2.1, on chromosome 2 is responsible for a large portion of phenotypic variation for this trait and that mutation(s) at this locus were likely involved in the evolution from allogamy to autogamy in this genus. To understand the genetic and molecular basis of stigma exsertion, we have conducted a high-resolution mapping at the chromosome region harboring the se2.1 QTL. The results indicate that this is a compound locus, comprising at least five tightly linked genes, one controlling style length, three controlling stamen length, and the other affecting anther dehiscence, a taxonomic character used to distinguish Lycopersicon species from other solanaceous species. This cluster of genes may represent the vestiges of an ancient coadapted gene complex in controlling mating behavior.

Population structure and phylogeography of Solanum pimpinellifolium inferred from a nuclear gene

Phylogeographical studies are emerging as a powerful tool for understanding the population structure and evolution of wild relatives of crop species. Because of their value as genetic resources, there is great interest in exploring the distribution of variation in wild relatives of cultivated plants. In this study, we use sequence variation from the nuclear gene, fruit vacuolar invertase (Vac), to investigate the population history of Solanum pimpinellifolium. Solanum pimpinellifolium is a close relative of the cultivated tomato and has repeatedly served as a source of valuable traits for crop improvement. We sequenced the second intron of the Vac gene in 129 individuals, representing 16 populations from the northern half of Peru. Patterns of haplotype sharing among populations indicate that there is isolation by distance. However, there is no congruence between the geographical distribution of haplotypes and their genealogical relationships. Levels of outcrossing decrease towards the southernmost populations, as previously observed in an allozyme study. The geographical pattern of Vac variation supports a centre of origin in northern Peru for S. pimpinellifolium and a gradual colonization along the Pacific coast. This implies that inbreeding populations are derived from outcrossing ones and that variation present at the Vac locus predates the spread of S. pimpinellifolium. The expansion of cities and human agricultural activity in the habitat of S. pimpinellifolium currently pose a threat to the species.

Fossil evidence and phylogeny: the age of major angiosperm clades based on mesofossil and macrofossil evidence from Cretaceous deposits

The fossil record has played an important role in the history of evolutionary thought, has aided the determination of key relationships through mosaics, and has allowed an assessment of a number of ecological hypotheses. Nonetheless, expectations that it might accurately and precisely mirror the progression of taxa through time seem optimistic in light of the many factors potentially interfering with uniform preservation. In view of these limitations, attempts to use the fossil record to corroborate phylogenetic hypotheses based on extensive comparisons among extant taxa may be misplaced. Instead we suggest a method-minimum age node mapping-for combining reliable fossil evidence with hypotheses of phylogeny. We use this methodology in conjunction with a phylogeny for angiosperms to assess timing in the history of major angiosperm clades. This method places many clades both with and without fossil records in temporal perspective, reveals discrepancies among clades in propensities for preservation, and raises some interesting questions about angiosperm evolution. By providing a context for understanding the gaps in the angiosperm fossil record this technique lends credibility and support to the remainder of the angiosperm record and to its applications in understanding a variety of aspects of angiosperm history. In effect, this methodology empowers the fossil record.

Relationships, origin, and diversity of Galapagos tomatoes: implications for the conservation of natural populations

Endemic Galápagos tomatoes (Lycopersicon cheesmanii) are of great value for cultivated tomato (L. esculentum) breeding, and therefore their conservation is of significance. Although within L. cheesmanii there is heterogeneity for many traits and formal infraspecific classification is not justified, here we distinguish three forms, without taxonomic significance, of L. cheesmanii that are of interest to breeders because of their distinctive morphology and habitat preferences: L. cheesmanii 'short' (one- to two-pinnate leaves, short internodes, and coastal habitats), L. cheesmanii 'long' (one- to two-pinnate leaves, long internodes, and inland habitats), and L. cheesmanii forma minor (three- to four-pinnate leaves, short internodes, and coastal habitats). In a recent survey of tomato populations in the Galápagos Islands, we found that several populations of L. cheesmanii reported 30-50 years earlier had disappeared, mostly as a consequence of human activity. In addition, a previously unreported invasive wild red-fruited form, which we named L. esculentum 'Gal cer,' was found on the island of Santa Cruz. The total diversity (estimated with amplified fragment length polymorphisms [AFLPs]) within L. cheesmanii (H(T) = 0.051) is almost as high as that for the mainland wild species L. pimpinellifolium (H(T) = 0.072). Lycopersicon esculentum 'Gal cer,' on the other hand, has a much lower diversity (H(T) = 0.014). Comparison of AFLP fragments shared by L. esculentum 'Gal cer' with other species showed that it is closely related to weedy tomato L. esculentum var. cerasiforme and, therefore, likely of recent origin. Genetic differentiation among the three native L. cheesmanii forms is low (G(ST) = 0.235), indicating that they share a common genetic background. Nonetheless, L. cheesmanii 'short' is about twice as diverse as L. cheesmanii 'long' or L. cheesmanii f. minor. UPGMA cluster and principal components analysis distinguish four groups within Eulycopersicon: L. pimpinellifolium, cultivated L. esculentum, L. esculentum var. cerasiforme including L. esculentum 'Gal cer,' and L. cheesmanii. The geographic distance and genetic distance in the wild forms of Galápagos tomatoes were not correlated. Apart from the pressure of humans, some native L. cheesmanii populations, especially L. cheesmanii 'long,' might be displaced by invasive L. esculentum 'Gal cer' because they share a similar habitat. We did not find evidence of intercrossing of L. cheesmanii with introduced L. esculentum, but occasional hybridization that contributes to loss of genetic integrity of L. cheesmanii cannot be ruled out. Establishment of reserves of L. cheesmanii to protect this species from introduced herbivorous animals and from hybridization with L. esculentum 'Gal cer' would help to conserve L. cheesmanii. Furthermore, accessions collected by C. M. Rick and others in the 1950s-1970s and now stored in germplasm banks could be used to reinstate some extinct populations.

Comparative Sequencing in the Genus Lycopersicon: Implications for the Evolution of Fruit Size in the Domestication of Cultivated Tomatoes

Sequence variation was sampled in cultivated and related wild forms of tomato at fw2.2—a fruit weight QTL key to the evolution of domesticated tomatoes. Variation at fw2.2 was contrasted with variation at four other loci not involved in fruit weight determination. Several conclusions could be reached: (1) Fruit weight variation attributable to fw2.2 is not caused by variation in the FW2.2 protein sequence; more likely, it is due to transcriptional variation associated with one or more of eight nucleotide changes unique to the promoter of large-fruit alleles; (2) fw2.2 and loci not involved in fruit weight have not evolved at distinguishably different rates in cultivated and wild tomatoes, despite the fact that fw2.2 was likely a target of selection during domestication; (3) molecular-clock-based estimates suggest that the large-fruit allele of fw2.2, now fixed in most cultivated tomatoes, arose in tomato germplasm long before domestication; (4) extant accessions of L. esculentum var. cerasiforme, the subspecies thought to be the most likely wild ancestor of domesticated tomatoes, appear to be an admixture of wild and cultivated tomatoes rather than a transitional step from wild to domesticated tomatoes; and (5) despite the fact that cerasiforme accessions are polymorphic for large- and small-fruit alleles at fw2.2, no significant association was detected between fruit size and fw2.2 genotypes in the subspecies—as tested by association genetic studies in the relatively small sample studied—suggesting the role of other fruit weight QTL in fruit weight variation in cerasiforme.

Ascertaining maternal and paternal lineage within Musa by chloroplast and mitochondrial DNA RFLP analyses

In banana, the maternal transmission of chloroplast DNA and paternal transmission of the mitochondrial DNA provides an exceptional opportunity for studying the maternal and paternal lineage of clones. In the present study, RFLP combined with hybridization of heterologous mitochondrial and chloroplastic probes have been used to characterize 71 wild accessions and 131 diploid and 103 triploid cultivated clones. In additon to Musa acuminata and Musa balbisiana, other species from the four Musa sections were studied to investigate their contribution to the origin of cultivated bananas. These molecular analyses enable the classification of the Musa complex to be discussed. Results ascertain relationships among and between the wild accessions and the mono- and interspecific diploid and triploid bananas, particularly for the acuminata genome. Parthenocarpic varieties are shown to be linked to M. acuminata banksii and M. acuminata errans, thus suggesting that the first center of domestication was in the Philippines - New Guinea area.