A new cacao linkage map based on codominant markers: development and integration of 201 new microsatellite markers

Interactions of Different Species of Phytophthora with Cacao Induce Genetic, Biochemical, and Morphological Plant Alterations

Diseases associated with Phytophthora cause considerable losses in cocoa production worldwide. Analyzing genes, proteins, and metabolites involved in Theobroma cacao's interaction with Phytophthora species is essential to explaining the molecular aspects of plant defense. Through a systematic literature review, this study aims to identify reports of genes, proteins, metabolites, morphological characteristics, and molecular and physiological processes of T. cacao involved in its interaction with species of Phytophthora. After the searches, 35 papers were selected for the data extraction stage, according to pre-established inclusion and exclusion criteria. In these studies, 657 genes and 32 metabolites, among other elements (molecules and molecular processes), were found to be involved in the interaction. The integration of this information resulted in the following conclusions: the expression patterns of pattern recognition receptors (PRRs) and a possible gene-to-gene interaction participate in cocoa resistance to Phytophthora spp.; the expression pattern of genes that encode pathogenesis-related (PRs) proteins is different between resistant and susceptible genotypes; phenolic compounds play an important role in preformed defenses; and proline accumulation may be involved in cell wall integrity. Only one proteomics study of T. cacao-Phytophthora spp. was found, and some genes proposed via QTL analysis were confirmed in transcriptomic studies.

Molecular characterization of national cocoa collection from the leading traditional growing areas in Ecuador

Ecuador is the leading producer and exporter of fine cocoa, with plantations over 80 years old, preserving distinctive aroma and flavor characteristics. The research objective was to screen the genetic variability of a collection of National cocoa from Ecuador's leading traditional cocoa growing areas, denominated as Centennial National Cocoa Plants (CCNC). This germplasm collection with 243 accessions was analyzed with 20 microsatellites (SSR) markers. DNA genotyping was highly informative, generating a total of 109 SSR alleles with an average of 5.5 alleles per locus. Only 0.8% of duplicate accessions were identified. The average genetic diversity obtained was 0.447, and the polymorphic content index was 0.414, which shows a high genetic diversity. The clustering, main coordinates, and population assignment analysis revealed that the samples are classified into two subpopulations (GN and GM), differentiated by their level of heterozygosity, with a fixation index value of 0.105. The results showed that microsatellite markers and statistical tools provide useful information that favors managing and conserving genetic variability in CCNC collection. Keywords: fine and aroma cocoa, Sabor Arriba, DNA genotyping, SSR markers

Sequencing and analysis of the complete mitochondrial genomes of Toona sinensis and Toona ciliata reveal evolutionary features of Toona

BACKGROUND

Toona is a critical genus in the Meliaceae, and the plants of this group are an asset for both restorative and restorative purposes, the most flexible of which are Toona sinensis and Toona ciliata. To concentrate on the advancement of mitochondrial(Mt) genome variety in T.sinensis and T.ciliata, the Mt genomes of the two species were sequenced in high throughput independently, after de novo assembly and annotation to construct a Mt genome map for comparison in genome structure. Find their repetitive sequences and analyze them in comparison with the chloroplast genome, along with Maximum-likelihood(ML) phylogenetic analysis with 16 other relatives.

RESULTS

(1) T. sinensis and T.ciliata are both circular structures with lengths of 683482 bp and 68300 bp, respectively. They share a high degree of similarity in encoding genes and have AT preferences. All of them have the largest Phe concentration and are the most frequently used codons. (2) Both of their Mt genome are highly preserved in terms of structural and functional genes, while the main variability is reflected in the length of tRNA, the number of genes, and the value of RSCU. (3) T. siniensis and T. ciliata were detected to have 94 and 87 SSRs, respectively, of which mononucleotides accounted for the absolute proportion. Besides, the vast majority of their SSRs were found to be poly-A or poly-T. (4)10 and 11 migrating fragments were identified in the comparison with the chloroplast genome, respectively. (5) In the ML evolutionary tree, T.sinensis and T.ciliata clustered individually into a small branch with 100% support, reflecting two species of Toona are very similarly related to each other.

CONCLUSIONS

This research provides a basis for the exploitation of T.sinensis and T.ciliata in terms of medicinal, edible, and timber resources to avoid confusion; at the same time, it can explore the evolutionary relationship between the Toona and related species, which does not only have an important practical value, but also provides a theoretical basis for future hybrid breeding of forest trees, molecular markers, and evolutionary aspects of plants, which has great scientific significance.

The chloroplast genome of Farsetia hamiltonii Royle, phylogenetic analysis, and comparative study with other members of Clade C of Brassicaceae

BACKGROUND

Farsetia hamiltonii Royle is a medicinally important annual plant from the Cholistan desert that belongs to the tribe Anastaticeae and clade C of the Brassicaceae family. We provide the entire chloroplast sequence of F.hamiltonii, obtained using the Illumina HiSeq2500 and paired-end sequencing. We compared F. hamiltonii to nine other clade C species, including Farsetia occidentalis, Lobularia libyca, Notoceras bicorne, Parolinia ornata, Morettia canescens, Cochlearia borzaeana, Megacarpaea polyandra, Biscutella laevigata, and Iberis amara. We conducted phylogenetic research on the 22 Brassicaceae species, which included members from 17 tribes and six clades.

RESULTS

The chloroplast genome sequence of F.hamiltonii of 154,802 bp sizes with 36.30% GC content and have a typical structure comprised of a Large Single Copy (LSC) of 83,906 bp, a Small Single Copy (SSC) of 17,988 bp, and two copies of Inverted Repeats (IRs) of 26,454 bp. The genomes of F. hamiltonii and F. occidentalis show shared amino acid frequencies and codon use, RNA editing sites, simple sequence repeats, and oligonucleotide repeats. The maximum likelihood tree revealed Farsetia as a monophyletic genus, closely linked to Morettia, with a bootstrap score of 100. The rate of transversion substitutions (Tv) was higher than the rate of transition substitutions (Ts), resulting in Ts/Tv less than one in all comparisons with F. hamiltonii, indicating that the species are closely related. The rate of synonymous substitutions (Ks) was greater than non-synonymous substitutions (Ka) in all comparisons with F. hamiltonii, with a Ka/Ks ratio smaller than one, indicating that genes underwent purifying selection. Low nucleotide diversity values range from 0.00085 to 0.08516, and IR regions comprise comparable genes on junctions with minimal change, supporting the conserved status of the selected chloroplast genomes of the clade C of the Brassicaceae family. We identified ten polymorphic regions, including rps8-rpl14, rps15-ycf1, ndhG-ndhI, psbK-psbI, ccsA-ndhD, rpl36-rps8, petA-psbJ, ndhF-rpl32, psaJ-rpl3, and ycf1 that might be exploited to construct genuine and inexpensive to solve taxonomic discrepancy and understand phylogenetic relationship amongst Brassicaceae species.

CONCLUSION

The entire chloroplast sequencing of F. hamiltonii sheds light on the divergence of genic chloroplast sequences among members of the clade C. When other Farsetia species are sequenced in the future, the full F. hamiltonii chloroplast will be used as a source for comprehensive taxonomical investigations of the genus. The comparison of F. hamiltonii and other clade C species adds new information to the phylogenetic data and evolutionary processes of the clade. The results of this study will also provide further molecular uses of clade C chloroplasts for possible plant genetic modifications and will help recognise more Brassicaceae family species.

Complete chloroplast genome of Hordeum brevisubulatum: Genome organization, synonymous codon usage, phylogenetic relationships, and comparative structure analysis

BACKGROUND

Hordeum brevisubulatum, known as fine perennial forage, is used for soil salinity improvement in northern China. Chloroplast (cp) genome is an ideal model for assessing its genome evolution and the phylogenetic relationships. We de novo sequenced and analyzed the cp genome of H. brevisubulatum, providing a fundamental reference for further studies in genetics and molecular breeding.

RESULTS

The cp genome of H. brevisubulatum was 137,155 bp in length with a typical quadripartite structure. A total of 130 functional genes were annotated and the gene of accD was lost in the process of evolution. Among all the annotated genes, 16 different genes harbored introns and the genes of ycf3 and rps12 contained two introns. Parity rule 2 (PR2) plot analysis showed that majority of genes had a bias toward T over A in the coding strand in all five Hordeum species, and a slight G over C in the other four Hordeum species except for H. bogdanil. Additionally, 52 dispersed repeat sequences and 182 simple sequence repeats were identified. Moreover, some unique SSRs of each species could be used as molecular markers for further study. Compared to the other four Hordeum species, H. brevisubulatum was most closely related to H. bogdanii and its cp genome was relatively conserved. Moreover, inverted repeat regions (IRa and IRb) were less divergent than other parts and coding regions were relatively conserved compared to non-coding regions. Main divergence was presented at the SSC/IR border.

CONCLUSIONS

This research comprehensively describes the architecture of the H. brevisubulatum cp genome and improves our understanding of its cp biology and genetic diversity, which will facilitate biological discoveries and cp genome engineering.

Initial Complete Chloroplast Genomes of Alchemilla (Rosaceae): Comparative Analysis and Phylogenetic Relationships

The genus Alchemilla L., known for its medicinal and ornamental value, is widely distributed in the Holarctic regions with a few species found in Asia and Africa. Delimitation of species within Alchemilla is difficult due to hybridization, autonomous apomixes, and polyploidy, necessitating efficient molecular-based characterization. Herein, we report the initial complete chloroplast (cp) genomes of Alchemilla. The cp genomes of two African (Afromilla) species Alchemilla pedata and Alchemilla argyrophylla were sequenced, and phylogenetic and comparative analyses were conducted in the family Rosaceae. The cp genomes mapped a typical circular quadripartite structure of lengths 152,438 and 152,427 base pairs (bp) in A. pedata and A. argyrophylla, respectively. Alchemilla cp genomes were composed of a pair of inverted repeat regions (IRa/IRb) of length 25,923 and 25,915 bp, separating the small single copy (SSC) region of 17,980 and 17,981 bp and a large single copy (LSC) region of 82,612 and 82,616 bp in A. pedata and A. argyrophylla, respectively. The cp genomes encoded 114 unique genes including 88 protein-coding genes, 37 transfer RNA (tRNA) genes, and 4 ribosomal RNA (rRNA) genes. Additionally, 88 and 95 simple sequence repeats (SSRs) and 37 and 40 tandem repeats were identified in A. pedata and A. argyrophylla, respectively. Significantly, the loss of group II intron in atpF gene in Alchemilla species was detected. Phylogenetic analysis based on 26 whole cp genome sequences and 78 protein-coding gene sequences of 27 Rosaceae species revealed a monophyletic clustering of Alchemilla nested within subfamily Rosoideae. Based on a protein-coding region, negative selective pressure (Ka/Ks < 1) was detected with an average Ka/Ks value of 0.1322 in A. argyrophylla and 0.1418 in A. pedata. The availability of complete cp genome in the genus Alchemilla will contribute to species delineation and further phylogenetic and evolutionary studies in the family Rosaceae.

The complete chloroplast genome sequence of yellow mustard (Sinapis alba L.) and its phylogenetic relationship to other Brassicaceae species

As a member of the large Brassicaceae family, yellow mustard (Sinapis alba L.) has been used as an important gene pool for the genetic improvement of cash crops in Brassicaceae. Understanding the phylogenetic relationship between Sinapis alba (S. alba) and other Brassicaceae crops can provide guidance on the introgression of its favorable alleles into related species. The chloroplast (cp) genome is an ideal model for assessing genome evolution and the phylogenetic relationships of complex angiosperm families. Herein, we de novo assembled the complete cp genome of S. alba by integrating the PacBio and Illumina sequencing platforms. A 153,760 bp quadripartite cycle without any gap was obtained, including a pair of inverted repeats (IRa and IRb) of 26,221 bp, separated by a large single copy (LSC) region of 83,506 bp and a small single copy (SSC) region of 17,821 bp. A total of 78 protein-coding genes, 30 tRNA genes, and four rRNA genes were identified in this cp genome, as were 89 simple sequence repeat (SSR) loci of 18 types. The codon usage analysis revealed a preferential use of the Leu codon with the A/U ending. The phylogenetic analysis using 82 Brassicaceae species demonstrated that S. alba had a close relationship with important Brassica and Raphanus species; moreover, it likely originated from a separate evolutionary pathway compared with the congeneric Sinapis arvensis. The synonymous (Ks) and non-synonymous (Ks) substitution rate analysis showed that genes encoding "Subunits of cytochrome b/f complex" were under the lowest purifying selection pressure, whereas those associated with "Maturase", "Subunit of acetyl-CoA", and "Subunits of NADH-dehydrogenase" underwent relatively higher purifying selection pressures. Our results provide valuable information for fully utilizing the S. alba cp genome as a potential genetic resource for the genetic improvement of Brassica and Raphanus species.

The complete chloroplast genome sequence of watercress (Nasturtium officinale R. Br.): Genome organization, adaptive evolution and phylogenetic relationships in Cardamineae

Watercress (Nasturtium officinale R. Br.), an aquatic leafy vegetable of the Brassicaceae family, is known as a nutritional powerhouse. Here, we de novo sequenced and assembled the complete chloroplast (cp) genome of watercress based on combined PacBio and Illumina data. The cp genome is 155,106 bp in length, exhibiting a typical quadripartite structure including a pair of inverted repeats (IRA and IRB) of 26,505 bp separated by a large single copy (LSC) region of 84,265 bp and a small single copy (SSC) region of 17,831 bp. The genome contained 113 unique genes, including 79 protein-coding genes, 30 tRNAs and 4 rRNAs, with 20 duplicate in the IRs. Compared with the prior cp genome of watercress deposited in GenBank, 21 single nucleotide polymorphisms (SNPs) and 27 indels were identified, mainly located in noncoding sequences. A total of 49 repeat structures and 71 simple sequence repeats (SSRs) were detected. Codon usage showed a bias for A/T-ending codons in the cp genome of watercress. Moreover, 45 RNA editing sites were predicted in 16 genes, all for C-to-U transitions. A comparative plastome study with Cardamineae species revealed a conserved gene order and high similarity of protein-coding sequences. Analysis of the Ka/Ks ratios of Cardamineae suggested positive selection exerted on the ycf2 gene in watercress, which might reflect specific adaptations of watercress to its particular living environment. Phylogenetic analyses based on complete cp genomes and common protein-coding genes from 56 species showed that the genus Nasturtium was a sister to Cardamine in the Cardamineae tribe. Our study provides valuable resources for future evolution, population genetics and molecular biology studies of watercress.

Cacao biotechnology: current status and future prospects

Theobroma cacao-The Food of the Gods, provides the raw material for the multibillion dollar chocolate industry and is also the main source of income for about 6 million smallholders around the world. Additionally, cocoa beans have a number of other nonfood uses in the pharmaceutical and cosmetic industries. Specifically, the potential health benefits of cocoa have received increasing attention as it is rich in polyphenols, particularly flavonoids. At present, the demand for cocoa and cocoa-based products in Asia is growing particularly rapidly and chocolate manufacturers are increasing investment in this region. However, in many Asian countries, cocoa production is hampered due to many reasons including technological, political and socio-economic issues. This review provides an overview of the present status of global cocoa production and recent advances in biotechnological applications for cacao improvement, with special emphasis on genetics/genomics, in vitro embryogenesis and genetic transformation. In addition, in order to obtain an insight into the latest innovations in the commercial sector, a survey was conducted on granted patents relating to T. cacao biotechnology.

A Larger Chocolate Chip-Development of a 15K Theobroma cacao L. SNP Array to Create High-Density Linkage Maps

Cacao (Theobroma cacao L.) is an important cash crop in tropical regions around the world and has a rich agronomic history in South America. As a key component in the cosmetic and confectionary industries, millions of people worldwide use products made from cacao, ranging from shampoo to chocolate. An Illumina Infinity II array was created using 13,530 SNPs identified within a small diversity panel of cacao. Of these SNPs, 12,643 derive from variation within annotated cacao genes. The genotypes of 3,072 trees were obtained, including two mapping populations from Ecuador. High-density linkage maps for these two populations were generated and compared to the cacao genome assembly. Phenotypic data from these populations were combined with the linkage maps to identify the QTLs for yield and disease resistance.

Genome of wild olive and the evolution of oil biosynthesis

Here we present the genome sequence and annotation of the wild olive tree (Olea europaea var. sylvestris), called oleaster, which is considered an ancestor of cultivated olive trees. More than 50,000 protein-coding genes were predicted, a majority of which could be anchored to 23 pseudochromosomes obtained through a newly constructed genetic map. The oleaster genome contains signatures of two Oleaceae lineage-specific paleopolyploidy events, dated at ∼28 and ∼59 Mya. These events contributed to the expansion and neofunctionalization of genes and gene families that play important roles in oil biosynthesis. The functional divergence of oil biosynthesis pathway genes, such as FAD2, SACPD, EAR, and ACPTE, following duplication, has been responsible for the differential accumulation of oleic and linoleic acids produced in olive compared with sesame, a closely related oil crop. Duplicated oleaster FAD2 genes are regulated by an siRNA derived from a transposable element-rich region, leading to suppressed levels of FAD2 gene expression. Additionally, neofunctionalization of members of the SACPD gene family has led to increased expression of SACPD2, 3, 5, and 7, consequently resulting in an increased desaturation of steric acid. Taken together, decreased FAD2 expression and increased SACPD expression likely explain the accumulation of exceptionally high levels of oleic acid in olive. The oleaster genome thus provides important insights into the evolution of oil biosynthesis and will be a valuable resource for oil crop genomics.

Report on the development of putative functional SSR and SNP markers in passion fruits

BACKGROUND

Passionflowers Passiflora edulis and Passiflora alata are diploid, outcrossing and understudied fruit bearing species. In Brazil, passion fruit cultivation began relatively recently and has earned the country an outstanding position as the world's top producer of passion fruit. The fruit's main economic value lies in the production of juice, an essential exotic ingredient in juice blends. Currently, crop improvement strategies, including those for underexploited tropical species, tend to incorporate molecular genetic approaches. In this study, we examined a set of P. edulis transcripts expressed in response to infection by Xanthomonas axonopodis, (the passion fruit's main bacterial pathogen that attacks the vines), aiming at the development of putative functional markers, i.e. SSRs (simple sequence repeats) and SNPs (single nucleotide polymorphisms).

RESULTS

A total of 210 microsatellites were found in 998 sequences, and trinucleotide repeats were found to be the most frequent (31.4%). Of the sequences selected for designing primers, 80.9% could be used to develop SSR markers, and 60.6% SNP markers for P. alata. SNPs were all biallelic and found within 15 gene fragments of P. alata. Overall, gene fragments generated 10,003 bp. SNP frequency was estimated as one SNP every 294 bp. Polymorphism rates revealed by SSR and SNP loci were 29.4 and 53.6%, respectively.

CONCLUSIONS

Passiflora edulis transcripts were useful for the development of putative functional markers for P. alata, suggesting a certain level of sequence conservation between these cultivated species. The markers developed herein could be used for genetic mapping purposes and also in diversity studies.

The Complete Chloroplast Genome of Chinese Bayberry (Morella rubra, Myricaceae): Implications for Understanding the Evolution of Fagales

Morella rubra (Myricaceae), also known as Chinese bayberry, is an economically important, subtropical, evergreen fruit tree. The phylogenetic placement of Myricaceae within Fagales and the origin of Chinese bayberry's domestication are still unresolved. In this study, we report the chloroplast (cp) genome of M. rubra and take advantage of several previously reported chloroplast genomes from related taxa to examine patterns of evolution in Fagales. The cp genomes of three M. rubra individuals were 159,478, 159,568, and 159.586 bp in length, respectively, comprising a pair of inverted repeat (IR) regions (26,014-26,069 bp) separated by a large single-copy (LSC) region (88,683-88,809 bp) and a small single-copy (SSC) region (18,676-18,767 bp). Each cp genome encodes the same 111 unique genes, consisting of 77 different protein-coding genes, 30 transfer RNA genes and four ribosomal RNA genes, with 18 duplicated in the IRs. Comparative analysis of chloroplast genomes from four representative Fagales families revealed the loss of infA and the pseudogenization of ycf15 in all analyzed species, and rpl22 has been pseudogenized in M. rubra and Castanea mollissima, but not in Juglans regia or Ostrya rehderiana. The genome size variations are detected mainly due to the length of intergenic spacers rather than gene loss, gene pseudogenization, IR expansion or contraction. The phylogenetic relationships yielded by the complete genome sequences strongly support the placement of Myricaceae as sister to Juglandaceae. Furthermore, seven cpDNA markers (trnH-psbA, psbA-trnK, rps2-rpoC2, ycf4-cemA, petD-rpoA, ndhE-ndhG, and ndhA intron) with relatively high levels of variation and variable cpSSR loci were identified within M. rubra, which will be useful in future research characterizing the population genetics of M. rubra and investigating the origin of domesticated Chinese bayberry.

Assembly and analysis of the complete Salix purpurea L. (Salicaceae) mitochondrial genome sequence

Plant mitochondrial (mt) genomes possess several complex features, including a variable size, a dynamic genome structure, and complicated patterns of gene loss and gain throughout evolutionary history. Studies of plant mt genomes can, therefore, provide unique insights into organelle evolution. We assembled the complete Salix purpurea L. mt genome by screening genomic sequence reads generated by a Roche-454 pyrosequencing platform. The pseudo-molecule obtained has a typical circular structure 598,970 bp long, with an overall GC content of 55.06%. The S. purpurea mt genome contains 52 genes: 31 protein-coding, 18 tRNAs, and three rRNAs. Eighteen tandem repeats and 404 microsatellites are distributed unevenly throughout the S. purpurea mt genome. A phylogenetic tree of 23 representative terrestrial plants strongly supports S. purpurea inclusion in the Malpighiales clade. Our analysis contributes toward understanding the organization and evolution of organelle genomes in Salicaceae species.

The Whole Genome Assembly and Comparative Genomic Research of Thellungiella parvula (Extremophile Crucifer) Mitochondrion

The complete nucleotide sequences of the mitochondrial (mt) genome of an extremophile species Thellungiella parvula (T. parvula) have been determined with the lengths of 255,773 bp. T. parvula mt genome is a circular sequence and contains 32 protein-coding genes, 19 tRNA genes, and three ribosomal RNA genes with a 11.5% coding sequence. The base composition of 27.5% A, 27.5% T, 22.7% C, and 22.3% G in descending order shows a slight bias of 55% AT. Fifty-three repeats were identified in the mitochondrial genome of T. parvula, including 24 direct repeats, 28 tandem repeats (TRs), and one palindromic repeat. Furthermore, a total of 199 perfect microsatellites have been mined with a high A/T content (83.1%) through simple sequence repeat (SSR) analysis and they were distributed unevenly within this mitochondrial genome. We also analyzed other plant mitochondrial genomes' evolution in general, providing clues for the understanding of the evolution of organelles genomes in plants. Comparing with other Brassicaceae species, T. parvula is related to Arabidopsis thaliana whose characters of low temperature resistance have been well documented. This study will provide important genetic tools for other Brassicaceae species research and improve yields of economically important plants.

Identification of candidate genes involved in Witches' broom disease resistance in a segregating mapping population of Theobroma cacao L. in Brazil

Background

Witches’ broom disease (WBD) caused by the fungus Moniliophthora perniciosa is responsible for considerable economic losses for cacao producers. One of the ways to combat WBD is to plant resistant cultivars. Resistance may be governed by a few genetic factors, mainly found in wild germplasm.

Results

We developed a dense genetic linkage map with a length of 852.8 cM that contains 3,526 SNPs and is based on the MP01 mapping population, which counts 459 trees from a cross between the resistant ‘TSH 1188’ and the tolerant ‘CCN 51’ at the Mars Center for Cocoa Science in Barro Preto, Bahia, Brazil. Seven quantitative trait loci (QTL) that are associated with WBD were identified on five different chromosomes using a multi-trait QTL analysis for outbreeders. Phasing of the haplotypes at the major QTL region on chromosome IX on a diversity panel of genotypes clearly indicates that the major resistance locus comes from a well-known source of WBD resistance, the clone ‘SCAVINA 6’. Various potential candidate genes identified within all QTL may be involved in different steps leading to disease resistance. Preliminary expression data indicate that at least three of these candidate genes may play a role during the first 12 h after infection, with clear differences between ‘CCN 51’ and ‘TSH 1188’.

Conclusions

We combined the information from a large mapping population with very distinct parents that segregate for WBD, a dense set of mapped markers, rigorous phenotyping capabilities and the availability of a sequenced genome to identify several genomic regions that are involved in WBD resistance. We also identified a novel source of resistance that most likely comes from the ‘CCN 51’ parent. Thanks to the large population size of the MP01 population, we were able to pick up QTL and markers with relatively small effects that can contribute to the creation and selection of more tolerant/resistant plant material.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2415-x) contains supplementary material, which is available to authorized users.

The complete chloroplast genome sequence of the relict woody plant Metasequoia glyptostroboides Hu et Cheng

Metasequoia glyptostroboides Hu et Cheng is the only species in the genus Metasequoia Miki ex Hu et Cheng, which belongs to the Cupressaceae family. There were around 10 species in the Metasequoia genus, which were widely spread across the Northern Hemisphere during the Cretaceous of the Mesozoic and in the Cenozoic. M. glyptostroboides is the only remaining representative of this genus. Here, we report the complete chloroplast (cp) genome sequence and the cp genomic features of M. glyptostroboides. The M. glyptostroboides cp genome is 131,887 bp in length, with a total of 117 genes comprised of 82 protein-coding genes, 31 tRNA genes and four rRNA genes. In this genome, 11 forward repeats, nine palindromic repeats, and 15 tandem repeats were detected. A total of 188 perfect microsatellites were detected through simple sequence repeat (SSR) analysis and these were distributed unevenly within the cp genome. Comparison of the cp genome structure and gene order to those of several other land plants indicated that a copy of the inverted repeat (IR) region, which was found to be IR region A (IRA), was lost in the M. glyptostroboides cp genome. The five most divergent and five most conserved genes were determined and further phylogenetic analysis was performed among plant species, especially for related species in conifers. Finally, phylogenetic analysis demonstrated that M. glyptostroboides is a sister species to Cryptomeria japonica (L. F.) D. Don and to Taiwania cryptomerioides Hayata. The complete cp genome sequence information of M. glyptostroboides will be great helpful for further investigations of this endemic relict woody plant and for in-depth understanding of the evolutionary history of the coniferous cp genomes, especially for the position of M. glyptostroboides in plant systematics and evolution.

Food Fingerprinting: Characterization of the Ecuadorean Type CCN-51 of Theobroma cacao L. Using Microsatellite Markers

The cocoa type "Colección Castro Naranjal 51" (CCN-51) is known for its resistance to specific climate conditions and its high yield, but it shows a weaker flavor profile and therefore is marketed as bulk cocoa. In a previous study, the two cocoa types Arriba and CCN-51 could easily be distinguished, but differences among the CCN-51 samples were observed. This was unexpected, as CCN-51 is reported to be a clone. To confirm whether CCN-51 is a pure clone, 10 simple sequence repeats (SSR) located on the nuclear genome were used to analyze various CCN-51 samples in comparison to the cocoa varieties Arriba and Criollo. As expected, there are differences in the SSR pattern among CCN-51, Arriba, and Criollo, but a variability within the CCN-51 sample set was detected as well. The previously described sequence variation in the chloroplast genome was confirmed by a variability in the microsatellite loci of the nuclear genome for a comprehensive cultivar collection of CCN-51 of both bean and leaf samples. In summary, beneath somaclonal variation, misidentification of plant collections and also sexual reproduction of CCN-51 can be suggested.

DNA-based differentiation of the Ecuadorian cocoa types CCN-51 and Arriba based on sequence differences in the chloroplast genome

Two cocoa types, Arriba and CCN-51, are being cultivated in Ecuador. With regard to the unique aroma, Arriba is considered a fine cocoa type, while CCN-51 is a bulk cocoa because of its weaker aroma. Because it is being assumed that Arriba is mixed with CCN-51, there is an interest in the analytical differentiation of the two types. Two methods to identify CCN-51 adulterations in Arriba cocoa were developed on the basis of differences in the chloroplast DNA. On the one hand, a different repeat of the sequence TAAAG in the inverted repeat region results in a different length of amplicons for the two cocoa types, which can be detected by agarose gel electrophoresis, capillary gel electrophoresis, and denaturing high-performance liquid chromatography. On the other hand, single nucleotide polymorphisms (SNPs) between the CCN-51 and Arriba sequences represent restriction sites, which can be used for restriction fragment length polymorphism analysis. A semi-quantitative analysis based on these SNPs is feasible. A method for an exact quantitation based on these results is not realizable. These sequence variations were confirmed for a comprehensive cultivar collection of Arriba and CCN-51, for both bean and leaf samples.

Pervasive effects of a dominant foliar endophytic fungus on host genetic and phenotypic expression in a tropical tree

It is increasingly recognized that macro-organisms (corals, insects, plants, vertebrates) consist of both host tissues and multiple microbial symbionts that play essential roles in their host's ecological and evolutionary success. Consequently, identifying benefits and costs of symbioses, as well as mechanisms underlying them are research priorities. All plants surveyed under natural conditions harbor foliar endophytic fungi (FEF) in their leaf tissues, often at high densities. Despite producing no visible effects on their hosts, experiments have nonetheless shown that FEF reduce pathogen and herbivore damage. Here, combining results from three genomic, and two physiological experiments, we demonstrate pervasive genetic and phenotypic effects of the apparently asymptomatic endophytes on their hosts. Specifically, inoculation of endophyte-free (E−) Theobroma cacao leaves with Colletotrichum tropicale (E+), the dominant FEF species in healthy T. cacao, induces consistent changes in the expression of hundreds of host genes, including many with known defensive functions. Further, E+ plants exhibited increased lignin and cellulose content, reduced maximum rates of photosynthesis (A_max), and enrichment of nitrogen-15 and carbon-13 isotopes. These phenotypic changes observed in E+ plants correspond to changes in expression of specific functional genes in related pathways. Moreover, a cacao gene (Tc00g04254) highly up-regulated by C. tropicale also confers resistance to pathogen damage in the absence of endophytes or their products in host tissues. Thus, the benefits of increased pathogen resistance in E+ plants are derived in part from up-regulation of intrinsic host defense responses, and appear to be offset by potential costs including reduced photosynthesis, altered host nitrogen metabolism, and endophyte heterotrophy of host tissues. Similar effects are likely in most plant-endophyte interactions, and should be recognized in the design and interpretation of genetic and phenotypic studies of plants.

The first genetic linkage map of Eucommia ulmoides

In accordance with pseudo-testcross strategy, the first genetic linkage map of Eucommia ulmoides Oliv. was constructed by an F1 population of 122 plants using amplified fragment length polymorphism (AFLP) markers. A total of 22 AFLP primer combinations generated 363 polymorphic markers. We selected 289 markers segregating as 1:1 and used them for constructing the parent-specific linkage maps. Among the candidate markers, 127 markers were placed on the maternal map LF and 108 markers on the paternal map Q1. The maternal map LF spanned 1116.1 cM in 14 linkage groups with a mean map distance of 8.78 cM; the paternal map Q1 spanned 929.6 cM in 12 linkage groups with an average spacing of 8.61 cM. The estimated coverage of the genome through two methods was 78.5 and 73.9% for LF, and 76.8 and 71.2% for Q1, respectively. This map is the first linkage map of E. ulmoides and provides a basis for mapping quantitative-trait loci and breeding applications.

Discovery and mapping of a new expressed sequence tag-single nucleotide polymorphism and simple sequence repeat panel for large-scale genetic studies and breeding of Theobroma cacao L

Theobroma cacao is an economically important tree of several tropical countries. Its genetic improvement is essential to provide protection against major diseases and improve chocolate quality. We discovered and mapped new expressed sequence tag-single nucleotide polymorphism (EST-SNP) and simple sequence repeat (SSR) markers and constructed a high-density genetic map. By screening 149 650 ESTs, 5246 SNPs were detected in silico, of which 1536 corresponded to genes with a putative function, while 851 had a clear polymorphic pattern across a collection of genetic resources. In addition, 409 new SSR markers were detected on the Criollo genome. Lastly, 681 new EST-SNPs and 163 new SSRs were added to the pre-existing 418 co-dominant markers to construct a large consensus genetic map. This high-density map and the set of new genetic markers identified in this study are a milestone in cocoa genomics and for marker-assisted breeding. The data are available at http://tropgenedb.cirad.fr.

The genome of Theobroma cacao

We sequenced and assembled the draft genome of Theobroma cacao, an economically important tropical-fruit tree crop that is the source of chocolate. This assembly corresponds to 76% of the estimated genome size and contains almost all previously described genes, with 82% of these genes anchored on the 10 T. cacao chromosomes. Analysis of this sequence information highlighted specific expansion of some gene families during evolution, for example, flavonoid-related genes. It also provides a major source of candidate genes for T. cacao improvement. Based on the inferred paleohistory of the T. cacao genome, we propose an evolutionary scenario whereby the ten T. cacao chromosomes were shaped from an ancestor through eleven chromosome fusions.

Mining of expressed sequence tag libraries of cacao for microsatellite markers using five computational tools

Expressed sequence tags (ESTs) provide researchers with a quick and inexpensive route for discovering new genes, data on gene expression and regulation, and also provide genic markers that help in constructing genome maps. Cacao is an important perennial crop of humid tropics. Cacao EST sequences, as available in the public domain, were downloaded and made into contigs. Microsatellites were located in these ESTs and contigs using five softwares (MISA, TRA, TROLL, SSRIT and SSR primer). MISA gave maximum coverage of SSRs in cacao ESTs and contigs, although TRA was able to detect higher order (5-mer) repeats. The frequency of SSRs was one per 26.9 kb in the known set of ESTs. One-third of the repeats in EST-contigs were found to be trimeric. A few rare repeats like 21-mer repeat were also located. A/T repeats were most abundant among the mononucleotide repeats and the AG/GA/TC/CT type was the most frequent among dimerics. Flanking primers were designed using Primer3 program and verified experimentally for PCR amplification. The results of the study are made available freely online database (http://riju.byethost31.com/cocoa/). Seven primer pairs amplified genomic DNA isolated from leaves were used to screen a representative set of 12 accessions of cacao.

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.

Development of a marker assisted selection program for cacao

ABSTRACT Production of cacao in tropical America has been severely affected by fungal pathogens causing diseases known as witches' broom (WB, caused by Moniliophthora perniciosa), frosty pod (FP, caused by M. roreri) and black pod (BP, caused by Phytophthora spp.). BP is pan-tropical and causes losses in all producing areas. WB is found in South America and parts of the Caribbean, while FP is found in Central America and parts of South America. Together, these diseases were responsible for over 700 million US dollars in losses in 2001 (4). Commercial cacao production in West Africa and South Asia are not yet affected by WB and FP, but cacao grown in these regions is susceptible to both. With the goal of providing new disease resistant cultivars the USDA-ARS and Mars, Inc. have developed a marker assisted selection (MAS) program. Quantitative trait loci have been identified for resistance to WB, FP, and BP. The potential usefulness of these markers in identifying resistant individuals has been confirmed in an experimental F(1) family in Ecuador.

Comparative analysis of expressed genes from cacao meristems infected by Moniliophthora perniciosa

BACKGROUND AND AIMS

Witches' broom disease is caused by the hemibiotrophic basidiomycete Moniliophthora perniciosa, and is one of the most important diseases of cacao in the western hemisphere. Because very little is known about the global process of such disease development, expressed sequence tags (ESTs) were used to identify genes expressed during the Theobroma cacao-Moniliophthora perniciosa interaction.

METHODS

Two cDNA libraries corresponding to the resistant (RT) and susceptible (SP) cacao-M. perniciosa interactions were constructed from total RNA, using the DB SMART Creator cDNA library kit (Clontech). Clones were randomly selected, sequenced from the 5' end and analysed using bioinformatics tools including in silico analysis of the differential gene expression.

KEY RESULTS

A total of 6884 ESTs were generated from the RT and SP cDNA libraries. These ESTs were composed of 2585 singlets and 341 contigs for a total of 2926 non-redundant sequences. The redundancy of the libraries was low and their specificity high when compared with the few other cacao libraries already published. Sequence analysis allowed the assignment of a putative functional category for 54 % of sequences, whereas approx. 22 % of sequences corresponded to unknown function and approx. 24 % of sequences did not show any significant similarity with other proteins present in the database. Despite the similar overall distribution of the sequences in functional categories between the two libraries, qualitative differences were observed. Genes involved during the defence response to pathogen infection or in programmed cell death were identified, such as pathogenesis related-proteins, trypsin inhibitor or oxalate oxidase, and some of them showed an in silico differential expression between the resistant and the susceptible interactions.

CONCLUSIONS

As far as is known this is the first EST resource from the cacao-M. perniciosa interaction and it is believed that it will provide a significant contribution to the understanding of the molecular mechanisms of the resistance and susceptibility of cacao to M. perniciosa, to develop strategies to control witches' broom, and as a source of polymorphism for molecular marker development and marker-assisted selection.

Adding value to cocoa (Theobroma cacao L.) germplasm information with domestication history and admixture mapping

A sound understanding of crop history can provide the basis for deriving novel genetic information through admixture mapping. We confirmed this, by using characterization data from an international collection of cocoa, collected 25 years ago, and from a contemporary plantation. We focus on the trees derived from three centuries of admixture between Meso-American Criollo and South American Forastero genomes. In both cacao sets of individuals, linkage disequilibrium extended over long genetic distances along chromosome regions, as expected in populations derived from recent admixture. Based on loose genome scans, genomic regions involved in useful traits were identified. Fifteen genomic regions involved in seed and fruit weight variation were highlighted. They correspond to ten previously identified QTLs and five novel ones. Admixture mapping can help to add value to genetic resources and thus, help to encourage investment in their conservation.

Microsatellite marker polymorphism and mapping in pea (Pisum sativum L.)

This paper aims at providing reliable and cost effective genotyping conditions, level of polymorphism in a range of genotypes and map position of newly developed microsatellite markers in order to promote broad application of these markers as a common set for genetic studies in pea. Optimal PCR conditions were determined for 340 microsatellite markers based on amplification in eight genotypes. Levels of polymorphism were determined for 309 of these markers. Compared to data obtained for other species, levels of polymorphism detected in a panel of eight genotypes were high with a mean number of 3.8 alleles per polymorphic locus and an average PIC value of 0.62, indicating that pea represents a rather polymorphic autogamous species. One of our main objectives was to locate a maximum number of microsatellite markers on the pea genetic map. Data obtained from three different crosses were used to build a composite genetic map of 1,430 cM (Haldane) comprising 239 microsatellite markers. These include 216 anonymous SSRs developed from enriched genomic libraries and 13 SSRs located in genes. The markers are quite evenly distributed throughout the seven linkage groups of the map, with 85% of intervals between the adjacent SSR markers being smaller than 10 cM. There was a good conservation of marker order and linkage group assignment across the three populations. In conclusion, we hope this report will promote wide application of these markers and will allow information obtained by different laboratories worldwide in diverse fields of pea genetics, such as QTL mapping studies and genetic resource surveys, to be easily aligned.

High throughput, high resolution selection of polymorphic microsatellite loci for multiplex analysis

BACKGROUND

Large-scale genetic profiling, mapping and genetic association studies require access to a series of well-characterised and polymorphic microsatellite markers with distinct and broad allele ranges. Selection of complementary microsatellite markers with non-overlapping allele ranges has historically proved to be a bottleneck in the development of multiplex microsatellite assays. The characterisation process for each microsatellite locus can be laborious and costly given the need for numerous, locus-specific fluorescent primers.

RESULTS

Here, we describe a simple and inexpensive approach to select useful microsatellite markers. The system is based on the pooling of multiple unlabelled PCR amplicons and their subsequent ligation into a standard cloning vector. A second round of amplification utilising generic labelled primers targeting the vector and unlabelled locus-specific primers targeting the microsatellite flanking region yield allelic profiles that are representative of all individuals contained within the pool. Suitability of various DNA pool sizes was then tested for this purpose. DNA template pools containing between 8 and 96 individuals were assessed for the determination of allele ranges of individual microsatellite markers across a broad population. This helped resolve the balance between using pools that are large enough to allow the detection of many alleles against the risk of including too many individuals in a pool such that rare alleles are over-diluted and so do not appear in the pooled microsatellite profile. Pools of DNA from 12 individuals allowed the reliable detection of all alleles present in the pool.

CONCLUSION

The use of generic vector-specific fluorescent primers and unlabelled locus-specific primers provides a high resolution, rapid and inexpensive approach for the selection of highly polymorphic microsatellite loci that possess non-overlapping allele ranges for use in large-scale multiplex assays.

Creation of BAC genomic resources for cocoa ( Theobroma cacao L.) for physical mapping of RGA containing BAC clones

We have constructed and validated the first cocoa ( Theobroma cacao L.) BAC library, with the aim of developing molecular resources to study the structure and evolution of the genome of this perennial crop. This library contains 36,864 clones with an average insert size of 120 kb, representing approximately ten haploid genome equivalents. It was constructed from the genotype Scavina-6 (Sca-6), a Forastero clone highly resistant to cocoa pathogens and a parent of existing mapping populations. Validation of the BAC library was carried out with a set of 13 genetically-anchored single copy and one duplicated markers. An average of nine BAC clones per probe was identified, giving an initial experimental estimation of the genome coverage represented in the library. Screening of the library with a set of resistance gene analogues (RGAs), previously mapped in cocoa and co-localizing with QTL for resistance to Phytophthora traits, confirmed at the physical level the tight clustering of RGAs in the cocoa genome and provided the first insights into the relationships between genetic and physical distances in the cocoa genome. This library represents an available BAC resource for structural genomic studies or map-based cloning of genes corresponding to important QTLs for agronomic traits such as resistance genes to major cocoa pathogens like Phytophthora spp ( palmivora and megakarya), Crinipellis perniciosa and Moniliophthora roreri.