The 'Tommy Atkins' mango genome reveals candidate genes for fruit quality.

BACKGROUND: Mango, Mangifera indica L., an important tropical fruit crop, is grown for its sweet and aromatic fruits. Past improvement of this species has predominantly relied on chance seedlings derived from over 1000 cultivars in the Indian sub-continent with a large variation for fruit size, yield, biotic and abiotic stress resistance, and fruit quality among other traits. Historically, mango has been an orphan crop with very limited molecular information. Only recently have molecular and genomics-based analyses enabled the creation of linkage maps, transcriptomes, and diversity analysis of large collections. Additionally, the combined analysis of genomic and phenotypic information is poised to improve mango breeding efficiency. RESULTS: This study sequenced, de novo assembled, analyzed, and annotated the genome of the monoembryonic mango cultivar 'Tommy Atkins'. The draft genome sequence was generated using NRGene de-novo Magic on high molecular weight DNA of 'Tommy Atkins', supplemented by 10X Genomics long read sequencing to improve the initial assembly. A hybrid population between 'Tommy Atkins' x 'Kensington Pride' was used to generate phased haplotype chromosomes and a highly resolved phased SNP map. The final 'Tommy Atkins' genome assembly was a consensus sequence that included 20 pseudomolecules representing the 20 chromosomes of mango and included ~ 86% of the ~ 439 Mb haploid mango genome. Skim sequencing identified ~ 3.3 M SNPs using the 'Tommy Atkins' x 'Kensington Pride' mapping population. Repeat masking identified 26,616 genes with a median length of 3348 bp. A whole genome duplication analysis revealed an ancestral 65 MYA polyploidization event shared with Anacardium occidentale. Two regions, one on LG4 and one on LG7 containing 28 candidate genes, were associated with the commercially important fruit size characteristic in the mapping population. CONCLUSIONS: The availability of the complete 'Tommy Atkins' mango genome will aid global initiatives to study mango genetics.

Genetic diversity of avocado (Persea americana Mill.) germplasm using pooled sequencing.

BACKGROUND: Discovering a genome-wide set of avocado (Persea americana Mill.) single nucleotide polymorphisms and characterizing the diversity of germplasm collection is a powerful tool for breeding. However, discovery is a costly process, due to loss of loci that are proven to be non-informative when genotyping the germplasm. RESULTS: Our study on a collection of 100 accessions comprised the three race types, Guatemalan, Mexican, and West Indian. To increase the chances of discovering polymorphic loci, three pools of genomic DNA, one from each race, were sequenced and the reads were aligned to a reference transcriptome. In total, 507,917 polymorphic loci were identified in the entire collection. Of these, 345,617 were observed in all three pools, 117,692 in two pools, 44,552 in one of the pools, and only 56 (0.0001%) were homozygous in the three pools but for different alleles. The polymorphic loci were validated using 192 randomly selected SNPs by genotyping the accessions within each pool. The sensitivity of polymorphic locus prediction ranged from 0.77 to 0.94. The correlation between the allele frequency estimated from the pooled sequences and actual allele frequency from genotype calling of individual accessions was r = 0.8. A subset of 109 SNPs were then used to evaluate the genetic relationships among avocado accessions and the genetic diversity of the collection. The three races were distinctly clustered by projecting the genetic variation on a PCA plot. As expected, by estimating the kinship coefficient for all the accessions, many of the cultivars from the California breeding program were closely related to each other, especially, the Hass-like ones. The green-skin avocados, e.g., 'Bacon', 'Zutano', 'Ettinger' and 'Fuerte' were also closely related to each other. CONCLUSIONS: A framework for SNP discovery and genetically characterizing of a breeder's accessions was described. Sequencing pools of gDNA is a cost-effective approach to create a genome-wide stock of polymorphic loci for a breeding program. Reassessing the botanical and the genetic knowledge about the germplasm accessions is valuable for future breeding. Kinship analysis may be used as a first step in finding a parental candidates in a parentage analyses.

iXora: exact haplotype inferencing and trait association.

BACKGROUND: We address the task of extracting accurate haplotypes from genotype data of individuals of large F1 populations for mapping studies. While methods for inferring parental haplotype assignments on large F1 populations exist in theory, these approaches do not work in practice at high levels of accuracy. RESULTS: We have designed iXora (Identifying crossovers and recombining alleles), a robust method for extracting reliable haplotypes of a mapping population, as well as parental haplotypes, that runs in linear time. Each allele in the progeny is assigned not just to a parent, but more precisely to a haplotype inherited from the parent. iXora shows an improvement of at least 15% in accuracy over similar systems in literature. Furthermore, iXora provides an easy-to-use, comprehensive environment for association studies and hypothesis checking in populations of related individuals. CONCLUSIONS: iXora provides detailed resolution in parental inheritance, along with the capability of handling very large populations, which allows for accurate haplotype extraction and trait association. iXora is available for non-commercial use from http://researcher.ibm.com/project/3430.

Interrelationships of petiolar air canal architecture, water depth, and convective air flow in Nymphaea odorata (Nymphaeaceae).

PREMISE OF THE STUDY: Nymphaea odorata grows in water up to 2 m deep, producing fewer larger leaves in deeper water. This species has a convective flow system that moves gases from younger leaves through submerged parts to older leaves, aerating submerged parts. Petiolar air canals are the convective flow pathways. This study describes the structure of these canals, how this structure varies with water depth, and models how convective flow varies with depth. • METHODS: Nymphaea odorata plants were grown at water depths from 30 to 90 cm. Lamina area, petiolar cross-sectional area, and number and area of air canals were measured. Field-collected leaves and leaves from juvenile plants were analyzed similarly. Using these data and data from the literature, we modeled how convective flow changes with water depth. • KEY RESULTS: Petioles of N. odorata produce two central pairs of air canals; additional pairs are added peripherally, and succeeding pairs are smaller. The first three pairs account for 96% of air canal area. Air canals form 24% of petiolar cross-sectional area. Petiolar and air canal cross-sectional areas increase with water depth. Petiolar area scales with lamina area, but the slope of this relationship is lower in 90 cm water than at shallower depths. In our model, the rate of convective flow varied with depth and with the balance of influx to efflux leaves. • CONCLUSIONS: Air canals in N. odorata petioles increase in size and number in deeper water but at a decreasing amount in relation to lamina area. Convective flow also depends on the number of influx to efflux laminae.

A genetically anchored physical framework for Theobroma cacao cv. Matina 1-6.

BACKGROUND: The fermented dried seeds of Theobroma cacao (cacao tree) are the main ingredient in chocolate. World cocoa production was estimated to be 3 million tons in 2010 with an annual estimated average growth rate of 2.2%. The cacao bean production industry is currently under threat from a rise in fungal diseases including black pod, frosty pod, and witches' broom. In order to address these issues, genome-sequencing efforts have been initiated recently to facilitate identification of genetic markers and genes that could be utilized to accelerate the release of robust T. cacao cultivars. However, problems inherent with assembly and resolution of distal regions of complex eukaryotic genomes, such as gaps, chimeric joins, and unresolvable repeat-induced compressions, have been unavoidably encountered with the sequencing strategies selected. RESULTS: Here, we describe the construction of a BAC-based integrated genetic-physical map of the T. cacao cultivar Matina 1-6 which is designed to augment and enhance these sequencing efforts. Three BAC libraries, each comprised of 10× coverage, were constructed and fingerprinted. 230 genetic markers from a high-resolution genetic recombination map and 96 Arabidopsis-derived conserved ortholog set (COS) II markers were anchored using pooled overgo hybridization. A dense tile path consisting of 29,383 BACs was selected and end-sequenced. The physical map consists of 154 contigs and 4,268 singletons. Forty-nine contigs are genetically anchored and ordered to chromosomes for a total span of 307.2 Mbp. The unanchored contigs (105) span 67.4 Mbp and therefore the estimated genome size of T. cacao is 374.6 Mbp. A comparative analysis with A. thaliana, V. vinifera, and P. trichocarpa suggests that comparisons of the genome assemblies of these distantly related species could provide insights into genome structure, evolutionary history, conservation of functional sites, and improvements in physical map assembly. A comparison between the two T. cacao cultivars Matina 1-6 and Criollo indicates a high degree of collinearity in their genomes, yet rearrangements were also observed. CONCLUSIONS: The results presented in this study are a stand-alone resource for functional exploitation and enhancement of Theobroma cacao but are also expected to complement and augment ongoing genome-sequencing efforts. This resource will serve as a template for refinement of the T. cacao genome through gap-filling, targeted re-sequencing, and resolution of repetitive DNA arrays.

Sequencing of a QTL-rich region of the Theobroma cacao genome using pooled BACs and the identification of trait specific candidate genes.

BACKGROUND: BAC-based physical maps provide for sequencing across an entire genome or a selected sub-genomic region of biological interest. Such a region can be approached with next-generation whole-genome sequencing and assembly as if it were an independent small genome. Using the minimum tiling path as a guide, specific BAC clones representing the prioritized genomic interval are selected, pooled, and used to prepare a sequencing library. RESULTS: This pooled BAC approach was taken to sequence and assemble a QTL-rich region, of ~3 Mbp and represented by twenty-seven BACs, on linkage group 5 of the Theobroma cacao cv. Matina 1-6 genome. Using various mixtures of read coverages from paired-end and linear 454 libraries, multiple assemblies of varied quality were generated. Quality was assessed by comparing the assembly of 454 reads with a subset of ten BACs individually sequenced and assembled using Sanger reads. A mixture of reads optimal for assembly was identified. We found, furthermore, that a quality assembly suitable for serving as a reference genome template could be obtained even with a reduced depth of sequencing coverage. Annotation of the resulting assembly revealed several genes potentially responsible for three T. cacao traits: black pod disease resistance, bean shape index, and pod weight. CONCLUSIONS: Our results, as with other pooled BAC sequencing reports, suggest that pooling portions of a minimum tiling path derived from a BAC-based physical map is an effective method to target sub-genomic regions for sequencing. While we focused on a single QTL region, other QTL regions of importance could be similarly sequenced allowing for biological discovery to take place before a high quality whole-genome assembly is completed.

A Highly Sensitive Method to Detect Avocado Sunblotch Viroid for the Maintenance of Infection-Free Avocado Germplasm Collections.

The United States Department of Agriculture (USDA) Agricultural Research Service (ARS) Subtropical Horticulture Research Station (SHRS) in Miami, FL holds a large germplasm collection of avocado (Persea americana). The recent threat of infection by laurel wilt has encouraged the creation of a backup collection at a disease-free site. Creating the backup collection is complicated by infection of some trees in the germplasm collection with avocado sunblotch viroid (ASBVd). Infected trees are frequently asymptomatic, necessitating the use of a molecular diagnostic assay. Although a reverse-transcription based assay already exists and has been used to assay all germplasm at the station, some trees showed inconsistent results. We have developed a more sensitive and specific assay involving pre-amplification of the entire viroid cDNA followed by detection using real-time PCR and a TaqMan assay. A second screening of all germplasm identified additional ASBVd -infected trees and allowed us to confidently remove these trees from the station. This method enables avocado germplasm curators to proceed with the creation of a viroid-free backup collection.

Genome-Wide Identification of Mango (Mangifera indica L.) Polygalacturonases: Expression Analysis of Family Members and Total Enzyme Activity During Fruit Ripening.

Mango (Mangifera indica L.) is an important commercial fruit that shows a noticeable loss of firmness during ripening. Polygalacturonase (PG, E.C. 3.2.1.15) is a crucial enzyme for cell wall loosening during fruit ripening since it solubilizes pectin and its activity correlates with fruit softening. Mango PGs were mapped to a genome draft using seventeen PGs found in mango transcriptomes and 48 bonafide PGs were identified. The phylogenetic analysis suggests that they are related to Citrus sinensis, which may indicate a recent evolutive divergence and related functions with orthologs in the tree. Gene expression analysis for nine PGs showed differential expression for them during post-harvest fruit ripening, MiPG21-1, MiPG14, MiPG69-1, MiPG17, MiPG49, MiPG23-3, MiPG22-7, and MiPG16 were highly up-regulated. PG enzymatic activity also increased during maturation and these results correlate with the loss of firmness observed in mango during post-harvest ripening, between the ethylene production burst and the climacteric peak. The analysis of PGs promoter regions identified regulatory sequences associated to ripening such as MADS-box, ethylene regulation like ethylene insensitive 3 (EIN3) factors, APETALA2-like and ethylene response element factors. During mango fruit ripening the action of at least these nine PGs contribute to softening, and their expression is regulated at the transcriptional level. The prediction of the tridimensional structure of some PGs showed a conserved parallel beta-helical fold related to polysaccharide hydrolysis and a modular architecture, where exons correspond to structural elements. Further biotechnological approaches could target specific softening-related PGs to extend mango post-harvest shelf life.

Genetic Map of Mango: A Tool for Mango Breeding.

Mango (Mangifera indica) is an economically and nutritionally important tropical/subtropical tree fruit crop. Most of the current commercial cultivars are selections rather than the products of breeding programs. To improve the efficiency of mango breeding, molecular markers have been used to create a consensus genetic map that identifies all 20 linkage groups in seven mapping populations. Polyembryony is an important mango trait, used for clonal propagation of cultivars and rootstocks. In polyembryonic mango cultivars, in addition to a zygotic embryo, several apomictic embryos develop from maternal tissue surrounding the fertilized egg cell. This trait has been associated with linkage group 8 in our consensus genetic map and has been validated in two of the seven mapping populations. In addition, we have observed a significant association between trait and single nucleotide polymorphism (SNP) markers for the vegetative trait of branch habit and the fruit traits of bloom, ground skin color, blush intensity, beak shape, and pulp color.

Application of Genome Wide Association and Genomic Prediction for Improvement of Cacao Productivity and Resistance to Black and Frosty Pod Diseases.

Chocolate is a highly valued and palatable confectionery product. Chocolate is primarily made from the processed seeds of the tree species Theobroma cacao. Cacao cultivation is highly relevant for small-holder farmers throughout the tropics, yet its productivity remains limited by low yields and widespread pathogens. A panel of 148 improved cacao clones was assembled based on productivity and disease resistance, and phenotypic single-tree replicated clonal evaluation was performed for 8 years. Using high-density markers, the diversity of clones was expressed relative to 10 known ancestral cacao populations, and significant effects of ancestry were observed in productivity and disease resistance. Genome-wide association (GWA) was performed, and six markers were significantly associated with frosty pod disease resistance. In addition, genomic selection was performed, and consistent with the observed extensive linkage disequilibrium, high predictive ability was observed at low marker densities for all traits. Finally, quantitative trait locus mapping and differential expression analysis of two cultivars with contrasting disease phenotypes were performed to identify genes underlying frosty pod disease resistance, identifying a significant quantitative trait locus and 35 differentially expressed genes using two independent differential expression analyses. These results indicate that in breeding populations of heterozygous and recently admixed individuals, mapping approaches can be used for low complexity traits like pod color cacao, or in other species single gene disease resistance, however genomic selection for quantitative traits remains highly effective relative to mapping. Our results can help guide the breeding process for sustainable improved cacao productivity.

Use of capillary array electrophoresis single-strand conformational polymorphism analysis to estimate genetic diversity of candidate genes in germplasm collections.

Capillary array electrophoresis single-strand conformation polymorphism (CAE-SSCP) analysis provides a reliable high-throughput method to genotype plant germplasm collections. Primers designed for highly conserved regions of candidate genes can be used to amplify DNA from plants in the collection. These amplified DNA fragments of identical length are turned into useful markers by assaying sequence differences by CAE-SSCP analysis. Sequence differences affect the electrophoretic mobility of single-stranded DNA under non-denaturing conditions. By collecting the mobility data for both strands assayed at two temperatures, alleles can be defined by mobility alone. For a germplasm collection with an unknown number of alleles at a locus, such mobility data of homozygotes can be used to determine the number of unique alleles without the necessity of cloning and sequencing each allele.

Mango (Mangifera indica L.) cv. Kent fruit mesocarp de novo transcriptome assembly identifies gene families important for ripening.

Fruit ripening is a physiological and biochemical process genetically programmed to regulate fruit quality parameters like firmness, flavor, odor and color, as well as production of ethylene in climacteric fruit. In this study, a transcriptomic analysis of mango (Mangifera indica L.) mesocarp cv. "Kent" was done to identify key genes associated with fruit ripening. Using the Illumina sequencing platform, 67,682,269 clean reads were obtained and a transcriptome of 4.8 Gb. A total of 33,142 coding sequences were predicted and after functional annotation, 25,154 protein sequences were assigned with a product according to Swiss-Prot database and 32,560 according to non-redundant database. Differential expression analysis identified 2,306 genes with significant differences in expression between mature-green and ripe mango [1,178 up-regulated and 1,128 down-regulated (FDR ≤ 0.05)]. The expression of 10 genes evaluated by both qRT-PCR and RNA-seq data was highly correlated (R = 0.97), validating the differential expression data from RNA-seq alone. Gene Ontology enrichment analysis, showed significantly represented terms associated to fruit ripening like "cell wall," "carbohydrate catabolic process" and "starch and sucrose metabolic process" among others. Mango genes were assigned to 327 metabolic pathways according to Kyoto Encyclopedia of Genes and Genomes database, among them those involved in fruit ripening such as plant hormone signal transduction, starch and sucrose metabolism, galactose metabolism, terpenoid backbone, and carotenoid biosynthesis. This study provides a mango transcriptome that will be very helpful to identify genes for expression studies in early and late flowering mangos during fruit ripening.

Promoting utilization of Saccharum spp. genetic resources through genetic diversity analysis and core collection construction.

Sugarcane (Saccharum spp.) and other members of Saccharum spp. are attractive biofuel feedstocks. One of the two World Collections of Sugarcane and Related Grasses (WCSRG) is in Miami, FL. This WCSRG has 1002 accessions, presumably with valuable alleles for biomass, other important agronomic traits, and stress resistance. However, the WCSRG has not been fully exploited by breeders due to its lack of characterization and unmanageable population. In order to optimize the use of this genetic resource, we aim to 1) genotypically evaluate all the 1002 accessions to understand its genetic diversity and population structure and 2) form a core collection, which captures most of the genetic diversity in the WCSRG. We screened 36 microsatellite markers on 1002 genotypes and recorded 209 alleles. Genetic diversity of the WCSRG ranged from 0 to 0.5 with an average of 0.304. The population structure analysis and principal coordinate analysis revealed three clusters with all S. spontaneum in one cluster, S. officinarum and S. hybrids in the second cluster and mostly non-Saccharum spp. in the third cluster. A core collection of 300 accessions was identified which captured the maximum genetic diversity of the entire WCSRG which can be further exploited for sugarcane and energy cane breeding. Sugarcane and energy cane breeders can effectively utilize this core collection for cultivar improvement. Further, the core collection can provide resources for forming an association panel to evaluate the traits of agronomic and commercial importance.

The genome sequence of the most widely cultivated cacao type and its use to identify candidate genes regulating pod color.

BACKGROUND: Theobroma cacao L. cultivar Matina 1-6 belongs to the most cultivated cacao type. The availability of its genome sequence and methods for identifying genes responsible for important cacao traits will aid cacao researchers and breeders. RESULTS: We describe the sequencing and assembly of the genome of Theobroma cacao L. cultivar Matina 1-6. The genome of the Matina 1-6 cultivar is 445 Mbp, which is significantly larger than a sequenced Criollo cultivar, and more typical of other cultivars. The chromosome-scale assembly, version 1.1, contains 711 scaffolds covering 346.0 Mbp, with a contig N50 of 84.4 kbp, a scaffold N50 of 34.4 Mbp, and an evidence-based gene set of 29,408 loci. Version 1.1 has 10x the scaffold N50 and 4x the contig N50 as Criollo, and includes 111 Mb more anchored sequence. The version 1.1 assembly has 4.4% gap sequence, while Criollo has 10.9%. Through a combination of haplotype, association mapping and gene expression analyses, we leverage this robust reference genome to identify a promising candidate gene responsible for pod color variation. We demonstrate that green/red pod color in cacao is likely regulated by the R2R3 MYB transcription factor TcMYB113, homologs of which determine pigmentation in Rosaceae, Solanaceae, and Brassicaceae. One SNP within the target site for a highly conserved trans-acting siRNA in dicots, found within TcMYB113, seems to affect transcript levels of this gene and therefore pod color variation. CONCLUSIONS: We report a high-quality sequence and annotation of Theobroma cacao L. and demonstrate its utility in identifying candidate genes regulating traits.

Evaluation of methods for de novo genome assembly from high-throughput sequencing reads reveals dependencies that affect the quality of the results.

Recent developments in high-throughput sequencing technology have made low-cost sequencing an attractive approach for many genome analysis tasks. Increasing read lengths, improving quality and the production of increasingly larger numbers of usable sequences per instrument-run continue to make whole-genome assembly an appealing target application. In this paper we evaluate the feasibility of de novo genome assembly from short reads (≤100 nucleotides) through a detailed study involving genomic sequences of various lengths and origin, in conjunction with several of the currently popular assembly programs. Our extensive analysis demonstrates that, in addition to sequencing coverage, attributes such as the architecture of the target genome, the identity of the used assembly program, the average read length and the observed sequencing error rates are powerful variables that affect the best achievable assembly of the target sequence in terms of size and correctness.

Phylogenetic analysis of seven WRKY genes across the palm subtribe Attaleinae (Arecaceae) [corrected] identifies Syagrus as sister group of the coconut.

BACKGROUND: The Cocoseae is one of 13 tribes of Arecaceae subfam. Arecoideae, and contains a number of palms with significant economic importance, including the monotypic and pantropical Cocos nucifera L., the coconut, the origins of which have been one of the "abominable mysteries" of palm systematics for decades. Previous studies with predominantly plastid genes weakly supported American ancestry for the coconut but ambiguous sister relationships. In this paper, we use multiple single copy nuclear loci to address the phylogeny of the Cocoseae subtribe Attaleinae, and resolve the closest extant relative of the coconut. METHODOLOGY/PRINCIPAL FINDINGS: We present the results of combined analysis of DNA sequences of seven WRKY transcription factor loci across 72 samples of Arecaceae tribe Cocoseae subtribe Attaleinae, representing all genera classified within the subtribe, and three outgroup taxa with maximum parsimony, maximum likelihood, and Bayesian approaches, producing highly congruent and well-resolved trees that robustly identify the genus Syagrus as sister to Cocos and resolve novel and well-supported relationships among the other genera of the Attaleinae. We also address incongruence among the gene trees with gene tree reconciliation analysis, and assign estimated ages to the nodes of our tree. CONCLUSIONS/SIGNIFICANCE: This study represents the as yet most extensive phylogenetic analyses of Cocoseae subtribe Attaleinae. We present a well-resolved and supported phylogeny of the subtribe that robustly indicates a sister relationship between Cocos and Syagrus. This is not only of biogeographic interest, but will also open fruitful avenues of inquiry regarding evolution of functional genes useful for crop improvement. Establishment of two major clades of American Attaleinae occurred in the Oligocene (ca. 37 MYBP) in Eastern Brazil. The divergence of Cocos from Syagrus is estimated at 35 MYBP. The biogeographic and morphological congruence that we see for clades resolved in the Attaleinae suggests that WRKY loci are informative markers for investigating the phylogenetic relationships of the palm family.

SSCP markers provide a useful alternative to microsatellites in genotyping and estimating genetic diversity in populations and germplasm collections of plant specialty crops.

For well-studied plant species with whole genome sequence or extensive EST data, SNP markers are the logical choice for both genotyping and whole genome association studies. However, SNP markers may not address the needs of researchers working on specialty crops with limited available genomic information. Microsatellite markers have been frequently employed due to their robustness, but marker development can be difficult and may result in few polymorphic markers. SSCP markers, such as microsatellites, are PCR-based and scored by electrophoretic mobility but, because they are based on SNPs rather than length differences, occur more frequently and are easier to develop than microsatellites. We have examined how well correlated the estimation of genetic diversity and genetic distance are in a population or germplasm collection when measured by 13 highly polymorphic microsatellite markers or 20 SSCP markers. We observed a significant correlation in pairwise genetic distances of 82 individuals in an international cacao germplasm collection (Mantel test Rxy=0.59, p<0.0001 for 10 000 permutations). Both sets of markers could distinguish each individual in the population. These data provide strong support for the use of SSCP markers in the genotyping of plant species where development of microsatellites would be difficult or expensive.

Geographic and genetic population differentiation of the Amazonian chocolate tree (Theobroma cacao L).

Numerous collecting expeditions of Theobroma cacao L. germplasm have been undertaken in Latin-America. However, most of this germplasm has not contributed to cacao improvement because its relationship to cultivated selections was poorly understood. Germplasm labeling errors have impeded breeding and confounded the interpretation of diversity analyses. To improve the understanding of the origin, classification, and population differentiation within the species, 1241 accessions covering a large geographic sampling were genotyped with 106 microsatellite markers. After discarding mislabeled samples, 10 genetic clusters, as opposed to the two genetic groups traditionally recognized within T. cacao, were found by applying Bayesian statistics. This leads us to propose a new classification of the cacao germplasm that will enhance its management. The results also provide new insights into the diversification of Amazon species in general, with the pattern of differentiation of the populations studied supporting the palaeoarches hypothesis of species diversification. The origin of the traditional cacao cultivars is also enlightened in this study.

The potential of the WRKY gene family for phylogenetic reconstruction: an example from the Malvaceae.

The WRKY gene family of transcription factors is involved in several diverse pathways and includes components of plant-specific, ancient regulatory networks. WRKY genes contain one or two highly conserved DNA binding domains interrupted by an intron. We used partial sequences of five independent WRKY loci to assess their potential for phylogeny reconstruction. Loci were originally isolated from Theobroma cacao L. by PCR with a single pair of degenerate primers; loci-specific primers were subsequently designed. We tested those loci across the sister genera Herrania Goudot and Theobroma L., with Guazuma ulmifolia Lam. as the outgroup. Overall, the combined WRKY matrices performed as well or better than other genes in resolving the intrageneric phylogeny of Herrania and Theobroma. The ease of isolating numerous, independent WRKY loci from diverse plant species with a single pair of degenerate primers designed to the highly conserved WRKY domain, renders them extremely useful tools for generating multiple, single or low copy nuclear loci for molecular phylogenetic studies at lower taxonomic levels. This is the first demonstration of the potential for members of the WRKY gene family for phylogenetic reconstruction.

Single-strand conformation polymorphism analysis of candidate genes for reliable identification of alleles by capillary array electrophoresis.

We investigated the reliability of capillary array electrophoresis-single strand conformation polymorphism (CAE-SSCP) to determine if it can be used to identify novel alleles of candidate genes in a germplasm collection. Both strands of three different size fragments (160, 245 and 437 bp) that differed by one or more nucleotides in sequence were analyzed at four different temperatures (18 degrees C, 25 degrees C, 30 degrees C, and 35 degrees C). Mixtures of amplified fragments of either the intron interrupting the C-terminal WRKY domain of the Tc10 locus or the NBS domain of the TcRGH1 locus of Theobroma cacao were electroinjected into all 16 capillaries of an ABI 3100 Genetic Analyzer and analyzed three times at each temperature. Multiplexing of samples of different size range is possible, as intermediate and large fragments were analyzed simultaneously in these experiments. A statistical analysis of the means of the fragment mobilities demonstrated that single-stranded conformers of the fragments could be reliably identified by their mobility at all temperatures and size classes. The order of elution of fragments was not consistent over strands or temperatures for the intermediate and large fragments. If samples are only run once at a single temperature, small fragments could be identified from a single strand at a single temperature. A combination of data from both strands of a single run was needed to identify correctly all four of the intermediate fragments and no combination of data from strands or temperatures would allow the correct identification of two large fragments that differed by only a single single-nucleotide polymorphism (SNP) from a single run. Thus, to adequately assess alleles at a candidate gene locus using SSCP on a capillary array, fragments should be < or =250 bp, samples should be analyzed at two different temperatures between 18 degrees C and 30 degrees C to reduce the variability introduced by the capillaries, data should be combined from both strands and both temperatures, and undenatured double-stranded (ds)DNA molecular weight standards, such as ROX 2500, should be included as internal standards.