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.

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.

Insights into the Key Aroma Compounds in Mango (Mangifera indica L. 'Haden') Fruits by Stable Isotope Dilution Quantitation and Aroma Simulation Experiments.

Thirty-four aroma-active compounds, previously identified with high flavor dilution factors by application of an aroma extract dilution analysis, were quantified in tree-ripened fruits of mango (Mangifera indica L. 'Haden'). From the results, the odor activity value (OAV) was calculated for each compound as the ratio of its concentration in the mangoes to its odor threshold in water. OAVs > 1 were obtained for 24 compounds, among which ethyl 2-methylbutanoate (fruity; OAV 2100), (3E,5Z)-undeca-1,3,5-triene (pineapple-like; OAV 1900), ethyl 3-methylbutanoate (fruity; OAV 1600), and ethyl butanoate (fruity; OAV 980) were the most potent, followed by (2E,6Z)-nona-2,6-dienal (cucumber-like), ethyl 2-methylpropanoate (fruity), (E)-ß-damascenone (cooked apple-like), ethyl hexanoate (fruity), 4-hydroxy-2,5-dimethyl-3(2H)-furanone (caramel-like), 3-methylbut-2-ene-1-thiol (sulfurous), γ-decalactone (peach-like), ß-myrcene (terpeny), (3Z)-hex-3-enal (green), 4-methyl-4-sulfanylpentan-2-one (tropical fruit-like), and ethyl octanoate (fruity). Aroma simulation and omission experiments revealed that these 15 compounds, when combined in a model mixture in their natural concentrations, were able to mimic the aroma of the fruits.

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.

Characterization of the major aroma-active compounds in mango (Mangifera indica L.) cultivars Haden, White Alfonso, Praya Sowoy, Royal Special, and Malindi by application of a comparative aroma extract dilution analysis.

The aroma-active compounds present in tree-ripened fruits of the five mango (Mangifera indica L.) cultivars Haden, White Alfonso, Praya Sowoy, Royal Special, and Malindi were isolated by solvent extraction followed by solvent-assisted flavor evaporation (SAFE) and analyzed by gas chromatography-olfactometery (GC-O). Application of a comparative aroma extract dilution analysis (cAEDA) afforded 54 aroma-active compounds in the flavor dilution (FD) factor range from 4 to ≥2048, 16 of which are reported for the first time in mango. The results of the identification experiments in combination with the FD factors revealed 4-hydroxy-2,5-dimethyl-3(2H)-furanone as an important aroma compound in all cultivars analyzed. Twenty-seven aroma-active compounds were present in at least one mango cultivar at an FD factor ≥128. Clear differences in the FD factors of these odorants between each of the mango cultivars suggested that they contributed to the unique sensory profiles of the individual cultivars.

Terpenoid variations within and among half-sibling avocado trees, Persea americana Mill. (Lauraceae).

Chemical analyses were conducted to determine the qualitative and quantitative differences in monoterpenes and sesquiterpenes in plant material from avocado trees, Persea americana Mill. (Lauraceae). The initial study analyzed plant material sampled from the trunk to the leaves through different branch diameters to quantify proximo-distal spatial differences within a tree. All trees were seedlings initiated from a single maternal tree. Two-way analysis of variance was conducted on 34 chemicals that comprised at least 3% of the total chemical content of at least one tree and/or location within a tree. There were significant interactions between genotype and location sampled for most chemicals. Parentage analysis using microsatellite molecular markers (SSR's) determined that the four trees had three fathers and that they represented two full-siblings and two half-sibling trees. Descriptive discriminant analysis found that both genotype and location within a tree could be separated based on chemical content, and that the chemical content from full-siblings tended to be more similar than chemical content from half-siblings. To further explore the relationship between genetic background and chemical content, samples were analyzed from leaf material from 20 trees that included two sets of full-sibling seedling trees, the maternal tree and the surviving paternal tree. Descriptive discriminant analysis found good separation between the two full-sibling groups, and that the separation was associated with chemistry of the parental trees. Six groups of chemicals were identified that explained the variation among the trees. We discuss the results in relation to the discrimination process used by wood-boring insects for site-selection on host trees, for tree selection among potential host trees, and the potential use of terpenoid chemical content in chemotaxonomy of avocado trees.

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.

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.

Isolation, characterization, and development of WRKY genes as useful genetic markers in Theobroma cacao.

There is currently an international effort in improving disease resistance and crop yield in Theobroma cacao L., an economically important crop of the tropics, using marker-assisted selection for breeding. We are developing molecular genetic markers focusing upon gene families involved with disease resistance. One such family is the WRKY proteins, which are plant-specific transcriptional factors associated with regulating defense responses to both abiotic and biotic stresses. Degenerate PCR primers were designed to the highly conserved DNA-binding domain and other conserved motifs of group I and group II, subgroups a-c, WRKY genes. Sixteen individual WRKY fragments were isolated from a mixture of T. cacao DNA using one pair of primers. Of the 16 WRKY loci investigated, seven contained single nucleotide polymorphisms within the intron as detected by sequence comparison of the PCR products. Four of these were successfully converted into molecular markers and mapped in an F2 population by capillary electrophoresis-single strand conformation polymorphism analysis. This is the first report of a pair of degenerate primers amplifying WRKY loci directly from genomic DNA and demonstrates a simple method for developing useful genetic markers from members of a large gene family.

Analysis of genetic diversity and population structure within Florida coconut (Cocos nucifera L.) germplasm using microsatellite DNA, with special emphasis on the Fiji Dwarf cultivar.

Using 15 simple sequence repeat (SSR) microsatellite DNA loci, we analyzed genetic variation within Cocos nucifera germplasm collections at two locations in south Florida, representing eight cultivars. The loci were also used in a parentage analysis of progeny of the 'Fiji Dwarf' variety at both locations. A total of 67 alleles were detected, with eight the highest number at any one locus. These loci identified 83 of the 110 individual palms. Gene diversity of the 15 loci ranged from 0.778 to 0.223, with a mean of 0.574. 'Fiji Dwarf', 'Malayan Dwarf', 'Green Niño' and 'Red Spicata' cultivars resolve as distinct clusters in a neighbor joining tree using modified Rogers distance, while the tall varieties form two aggregates. The highest gene diversity was found in the tall cultivars (H = 0.583 cumulatively), and the lowest in the 'Malayan Dwarf' (H = 0.202). After the tall coconuts, the 'Fiji Dwarf' was most genetically diverse (H = 0.436), and had the largest number of unique alleles. Genetic identity is highest among the 'Malayan Dwarf' phenotypes, and between the tall varieties. The 'Red Malayan Dwarf' is genetically distinct from the 'Green' and 'Yellow Malayan Dwarf' phenotypes, which cannot be distinguished with the SSR loci used. Off-type 'Malayan Dwarf' phenotypes (putative hybrids with talls) can be identified genotypically. Parentage analyses of 30 'Fiji Dwarf' progeny propagated from five adults surrounded by other cultivars estimate that only 20% of the progeny were out-crossed to the other varieties, while 40-46% were possible selfs. This suggests that a seed-production orchard of the variety maintained at reasonable distance from other varieties, will likely yield only 'Fiji Dwarf' genotypes. Our data are discussed in the context of hypotheses of coconut dissemination around the world.