Developing Single Nucleotide Polymorphism (SNP) Markers for the Identification of Coffee Germplasm

Coffee is one of the most widely consumed beverages and represents a multibillion-dollar global industry. Accurate identification of coffee cultivars is essential for efficient management, exchange, and use of coffee genetic resources. To date, a universal platform that can allow data comparison across different laboratories and genotyping platforms has not been developed by the coffee research community. Using expressed sequence tags (EST) of Coffea arabica, C. canephora and C. racemosa from public databases, we developed 7538 single nucleotide polymorphism (SNP) markers and selected 180 for validation using 25 C. arabica and C. canephora accessions from Puerto Rico. Based on the validation result, we designated a panel of 55 SNP markers that are polymorphic across the two species. The average minor allele frequency and information index of this SNP panel are 0.281 and 0.690, respectively. This panel enabled the differentiation of all tested accessions of C. canephora, which accounts for 79.2 % of the total polymorphism in the samples. Only 21.8 % of the polymorphic SNPs were detected in the 12 C. arabica cultivars, which, nonetheless, were able to unambiguously differentiate the 12 Arabica cultivars into ten unique genotypes, including two synonymous groups. Several local Puerto Rican cultivars with partial Timor pedigree, including Limaní, Frontón, and TARS 18087, showed substantial genetic difference from the other common Arabica cultivars, such as Catuai, Borbón, and Mundo Nuevo. This coffee SNP panel provides robust and universally comparable DNA fingerprints, thus can serve as a genotyping tool to assist coffee germplasm management, propagation of planting material, and coffee cultivar authentication.     

A 34K SNP genotyping array for Populus trichocarpa: Design,application to the study of natural populations and transferability to other Populus species

Genetic mapping of quantitative traits requires genotypic data for large numbers of markers in many individuals. For such studies, the use of large single nucleotide polymorphism (SNP) genotyping arrays still offers the most cost‐effective solution. Herein we report on the design and performance of a SNP genotyping array for Populus trichocarpa (black cottonwood). This genotyping array was designed with SNPs pre‐ascertained in 34 wild accessions covering most of the species latitudinal range. We adopted a candidate gene approach to the array design that resulted in the selection of 34 131 SNPs, the majority of which are located in, or within 2 kb of, 3543 candidate genes. A subset of the SNPs on the array (539) was selected based on patterns of variation among the SNP discovery accessions. We show that more than 95% of the loci produce high quality genotypes and that the genotyping error rate for these is likely below 2%. We demonstrate that even among small numbers of samples (n = 10) from local populations over 84% of loci are polymorphic. We also tested the applicability of the array to other species in the genus and found that the number of polymorphic loci decreases rapidly with genetic distance, with the largest numbers detected in other species in section Tacamahaca. Finally, we provide evidence for the utility of the array to address evolutionary questions such as intraspecific studies of genetic differentiation, species assignment and the detection of natural hybrids.     

Single‐nucleotide polymorphism discovery and diversity in the model legume Medicago truncatula

Extensive genomic resources are available in the model legume Medicago truncatula. Here, we present the discovery and design of the first array of single‐nucleotide polymorphism (SNP) markers in M. truncatula through large‐scale Sanger resequencing of genomic fragments spanning the genome, in a diverse panel of 16 M. truncatula accessions. Both anonymous fragments and fragments targeting candidate genes for flowering phenology and symbiosis were surveyed for nucleotide variation in almost 230 kb of unique genomic regions. A set of 384 SNP markers was designed for an Illumina's GoldenGate assay, genotyped on a collection of 192 inbred lines (CC192) representing the geographical range of the species and used to survey the diversity of two natural populations. Finally, 86% of the tested SNPs were of high quality and exhibited polymorphism in the CC192 collection. Even at the population level, we detected polymorphism for more than 50% of the selected SNPs. Analysis of the allele frequency spectrum in the CC192 showed a reduced ascertainment bias, mostly limited to very rare alleles (frequency <0.01). The substantial polymorphism detected at the species and population levels, the high marker quality and the potential to survey large samples of individuals make this set of SNP markers a valuable tool to improve our understanding of the effect of demographic and selective factors that shape the natural genetic diversity within the selfing species Medicago truncatula.     

Brewing trouble: coffee invasion in relation to edges and forest structure in tropical rainforest fragments of the Western Ghats, India

While the conservation impacts of invasive plant species on tropical biodiversity is widely recognised, little is known of the potential for cultivated crops turning invasive in tropical forest regions. In the Western Ghats biodiversity hotspot, India, fragmented rainforests often adjoin coffee plantations. This study in the Anamalai hills assessed the effects of distance from edges and forest structure on the occurrence and abundance of shade-tolerant coffee (Arabica Coffea arabica and Robusta C. canephora) in four fragments (32–200 ha) using replicate line transects laid from the edges into the interiors. The coffee species cultivated in adjoining plantations was more abundant than the other coffee species inside study fragments, showing a clear decline in stem density from edge (0 m) to interior (250 m), suggesting the influence of propagule pressure of adjoining plantations, coupled with edge effects and seed dispersal by animals. Significant positive correlations of coffee density with canopy cover indicate the potential threat of coffee invasion even in closed canopy rainforests. Stem density of Coffea arabica (150–1,825 stems/ha) was higher in more disturbed fragments, whereas Coffea canephora had spread in disturbed and undisturbed sites achieving much higher densities (6.3–11,486 stems/ha). In addition, a negative relationship between C. canephora and native shrub density indicates its potential detrimental effects on native plants.     

Molecular analysis of introgressive breeding in coffee (Coffea arabica L.)

Nineteen arabica coffee introgression lines (BC₁F₄) and two accessions derived from a spontaneous interspecific cross (i.e. Timor Hybrid) between Coffea arabica (2n=4x=44) and C. canephora (2n=2x=22) were analysed for the introgression of C. canephora genetic material. The Timor Hybrid-derived genotypes were evaluated by AFLP, using 42 different primer combinations, and compared to 23 accessions of C. arabica and 8 accessions of C. canephora. A total of 1062 polymorphic fragments were scored among the 52 accessions analysed. One hundred and seventy-eight markers consisting of 109 additional bands (i.e. introgressed markers) and 69 missing bands distinguished the group composed of the Timor Hybrid-derived genotypes from the accessions of C. arabica. AFLP therefore seemed to be an extremely efficient technique for DNA marker generation in coffee as well as for the detection of introgression in C. arabica. The genetic diversity observed in the Timor Hybrid-derived genotypes appeared to be approximately double that in C. arabica. Although representing only a small proportion of the genetic diversity available in C. canephora, the Timor Hybrid obviously constitutes a considerable source of genetic diversity for arabica breeding. Analysis of genetic relationships among the Timor Hybrid-derived genotypes suggested that introgression was not restricted to chromosome substitution but also involved chromosome recombinations. Furthermore, the Timor Hybrid-derived genotypes varied considerably in the number of AFLP markers attributable to introgression. In this way, the introgressed markers identified in the analysed arabica coffee introgressed genotypes were estimated to represent from 9% to 29% of the C. canephora genome. Nevertheless, the amount of alien genetic material in the introgression arabica lines remains substantial and should justify the development of adapted breeding strategies. Received: 2 February 1999 / Accepted: 12 May 1999     

Development,validation and genetic analysis of a large soybean SNP genotyping array

Cultivated soybean (Glycine max) suffers from a narrow germplasm relative to other crop species, probably because of under‐use of wild soybean (Glycine soja) as a breeding resource. Use of a single nucleotide polymorphism (SNP) genotyping array is a promising method for dissecting cultivated and wild germplasms to identify important adaptive genes through high‐density genetic mapping and genome‐wide association studies. Here we describe a large soybean SNP array for use in diversity analyses, linkage mapping and genome‐wide association analyses. More than four million high‐quality SNPs identified from high‐depth genome re‐sequencing of 16 soybean accessions and low‐depth genome re‐sequencing of 31 soybean accessions were used to select 180 961 SNPs for creation of the Axiom^® SoyaSNP array. Validation analysis for a set of 222 diverse soybean lines showed that 170 223 markers were of good quality for genotyping. Phylogenetic and allele frequency analyses of the validation set data indicated that accessions showing an intermediate morphology between cultivated and wild soybeans collected in Korea were natural hybrids. More than 90 unanchored scaffolds in the current soybean reference sequence were assigned to chromosomes using this array. Finally, dense average spacing and preferential distribution of the SNPs in gene‐rich chromosomal regions suggest that this array may be suitable for genome‐wide association studies of soybean germplasm. Taken together, these results suggest that use of this array may be a powerful method for soybean genetic analyses relating to many aspects of soybean breeding.     

Development of an integrated 200K SNP genotyping array and application for genetic mapping,genome assembly improvement and genome wide association studies in pear (Pyrus)

Pear (Pyrus; 2n = 34), the third most important temperate fruit crop, has great nutritional and economic value. Despite the availability of many genomic resources in pear, it is challenging to genotype novel germplasm resources and breeding progeny in a timely and cost‐effective manner. Genotyping arrays can provide fast, efficient and high‐throughput genetic characterization of diverse germplasm, genetic mapping and breeding populations. We present here 200K AXIOM^® PyrSNP, a large‐scale single nucleotide polymorphism (SNP) genotyping array to facilitate genotyping of Pyrus species. A diverse panel of 113 re‐sequenced pear genotypes was used to discover SNPs to promote increased adoption of the array. A set of 188 diverse accessions and an F₁ population of 98 individuals from ‘Cuiguan’ × ‘Starkrimson’ was genotyped with the array to assess its effectiveness. A large majority of SNPs (166 335 or 83%) are of high quality. The high density and uniform distribution of the array SNPs facilitated prediction of centromeric regions on 17 pear chromosomes, and significantly improved the genome assembly from 75.5% to 81.4% based on genetic mapping. Identification of a gene associated with flowering time and candidate genes linked to size of fruit core via genome wide association studies showed the usefulness of the array in pear genetic research. The newly developed high‐density SNP array presents an important tool for rapid and high‐throughput genotyping in pear for genetic map construction, QTL identification and genomic selection.     

Haplotype‐based genotyping‐by‐sequencing in oat genome research

In a de novo genotyping‐by‐sequencing (GBS) analysis of short, 64‐base tag‐level haplotypes in 4657 accessions of cultivated oat, we discovered 164741 tag‐level (TL) genetic variants containing 241224 SNPs. From this, the marker density of an oat consensus map was increased by the addition of more than 70000 loci. The mapped TL genotypes of a 635‐line diversity panel were used to infer chromosome‐level (CL) haplotype maps. These maps revealed differences in the number and size of haplotype blocks, as well as differences in haplotype diversity between chromosomes and subsets of the diversity panel. We then explored potential benefits of SNP vs. TL vs. CL GBS variants for mapping, high‐resolution genome analysis and genomic selection in oats. A combined genome‐wide association study (GWAS) of heading date from multiple locations using both TL haplotypes and individual SNP markers identified 184 significant associations. A comparative GWAS using TL haplotypes, CL haplotype blocks and their combinations demonstrated the superiority of using TL haplotype markers. Using a principal component‐based genome‐wide scan, genomic regions containing signatures of selection were identified. These regions may contain genes that are responsible for the local adaptation of oats to Northern American conditions. Genomic selection for heading date using TL haplotypes or SNP markers gave comparable and promising prediction accuracies of up to r = 0.74. Genomic selection carried out in an independent calibration and test population for heading date gave promising prediction accuracies that ranged between r = 0.42 and 0.67. In conclusion, TL haplotype GBS‐derived markers facilitate genome analysis and genomic selection in oat.     

Development of SNP‐genotyping arrays in two shellfish species

Use of SNPs has been favoured due to their abundance in plant and animal genomes, accompanied by the falling cost and rising throughput capacity for detection and genotyping. Here, we present in vitro (obtained from targeted sequencing) and in silico discovery of SNPs, and the design of medium‐throughput genotyping arrays for two oyster species, the Pacific oyster, Crassostrea gigas, and European flat oyster, Ostrea edulis. Two sets of 384 SNP markers were designed for two Illumina GoldenGate arrays and genotyped on more than 1000 samples for each species. In each case, oyster samples were obtained from wild and selected populations and from three‐generation families segregating for traits of interest in aquaculture. The rate of successfully genotyped polymorphic SNPs was about 60% for each species. Effects of SNP origin and quality on genotyping success (Illumina functionality Score) were analysed and compared with other model and nonmodel species. Furthermore, a simulation was made based on a subset of the C. gigas SNP array with a minor allele frequency of 0.3 and typical crosses used in shellfish hatcheries. This simulation indicated that at least 150 markers were needed to perform an accurate parental assignment. Such panels might provide valuable tools to improve our understanding of the connectivity between wild (and selected) populations and could contribute to future selective breeding programmes.     

Biotechnological applications for the improvement of coffee (Coffea arabica L.)

Summary The important advances in coffee biotechnological techniques which have been made particularly during the last 10yr could benefit the coffee breeder in practice and open new perspectives for the development of new varieties. The molecular phylogeny of Coffea species has been established using DNA sequence data. The molecular markers have revealed an extremely reduced genetic diversity in Coffea arabica L. in comparison to C. canephora. However, wild accessions collected in the Ethiopian highlands appeared to constitute a valuable gene reservoir. A complete genetic linkage map of C. canephora was reported and additional ones are being constructed, particularly on C. arabica. The integration of Molecular Assisted Selection in coffee breeding promises to drastically increase the efficiency of breeding programs. Economically important genes of the caffeine biosynthetic pathway or genes encoding for seed storage proteins have been isolated. The high performance already achieved in the in vitro propagation process by somatic embryogenesis offers the possibility to mass propagate superior hybrids in different countries of both C. arabica (selected F₁ hybrids) and C. canephora (rootstock variety). Pilot productions by somatic embryogenesis currently permit preparation for commercial application. Somaclonal variation was observed. The percentage of the off-types can vary between 3 and 10% depending on the genotype. Seed cryopreservation enables a routine use for long-term conservation of coffee genetic resources. Transgenic plants have been obtained for the C. arabica and C. canephora cultivated species through Agrobacterium-mediated transformation which constitutes the technique now currently used to transfer directly genes in coffee plants.     

Genome‐wide association study for body weight in cattle populations from Siberia

Body weight is a complex trait in cattle associated with commonly used commercial breeding measurements related to growth. Although many quantitative trait loci (QTL) for body weight have been identified in cattle so far, searching for genetic determinants in different breeds or environments is promising. Therefore, we carried out a genome‐wide association study (GWAS) in two cattle populations from the Russian Federation (Siberian region) using the GGP HD150K array containing 139 376 single nucleotide polymorphism (SNP) markers. Association tests for 107 550 SNPs left after filtering revealed five statistically significant SNPs on BTA5, considering a false discovery rate of less than 0.05. The chromosomal region containing these five SNPs contains the CCND2 gene, which was previously associated with average daily weight gain and body mass index in US beef cattle populations and in humans respectively. Our study is the first GWAS for body weight in beef cattle populations from the Russian Federation. The results provided here suggest that, despite the existence of breed‐ and species‐specific QTL, the genetic architecture of body weight could be evolutionarily conserved in mammals.     

Genome‐wide detection of genetic loci and candidate genes for teat number and body conformation traits at birth in Chinese Sushan pigs

Body conformation at birth and teat number are economically important traits in the pig industry, as these traits are usually explored to evaluate the growth and reproductive potential of piglets. To detect genetic loci and candidate genes for these traits, we performed a GWAS on 269 pigs from a recently developed Chinese breed (Sushan) using 38  128 informative SNPs on the Affymetrix Porcine SNP 55K Array. In total, we detected one genome‐wide significant (P = 1.31e‐6) SNP for teat number on chromosome X and 15 chromosome‐wide significant SNPs for teat number, body weight, body length, chest circumference and cannon circumference at birth on chromosomes 1, 3, 4, 6, 7, 9, 10, 13, 14, 15, 17 and 18. The most significant SNP had an additive effect of 0.74 × total teat number, explaining 20% of phenotypic variance. Five significant SNPs resided in the previously reported quantitative trait loci for these traits and seven significant SNPs had a pleiotropic effect on multiple traits. Intriguingly, 12 of the genes nearest to the significant SNPs are functionally related to body conformation and teat number traits, including SPRED2, MKX, TMSB4X and ESR1. GO analysis revealed that candidate genes proximal to the significant SNPs were enriched in the G‐protein coupled receptor and steroid hormone‐mediated signaling pathway. Our findings shed light on the genetic basis of the measured traits and provide molecular markers especially for the genetic improvement of teat number in Sushan and related pigs.     

Identification and characterization of expressed sequence tags‐derived simple sequence repeats markers from robusta coffee variety ‘CxR’ (an interspecific hybrid of Coffea canephora × Coffea congensis)

SSR (simple sequence repeats) markers derived from ESTs (expressed sequence tags), commonly called EST‐SSRs or genic SSRs provide useful genetic markers for crop improvement. These are easy and economical to develop as by‐products of large‐scale EST resources that have become available as part of the functional genomic studies in many plant species. Here, we describe for the first time, nine genic‐SSRs of coffee that are developed from the microsatellite containing ESTs from a cDNA library of moisture‐stressed leaves of coffee variety, ‘CxR’ (a commercial interspecific hybrid between Coffea congensis and Coffea canephora). The markers show considerable allelic diversity with PIC values up to 0.70 and 0.75 for Coffea arabica and Coffea canephora, respectively, and robust cross‐species amplification in 16 other related taxa of coffee. The validation studies thus demonstrate the potential utility of the EST‐SSRs for genetic analysis of coffee germplasm.     

New resources for genetic studies in Populus nigra: genome‐wide SNP discovery and development of a 12k Infinium array

Whole genome resequencing of 51 Populus nigra (L.) individuals from across Western Europe was performed using Illumina platforms. A total number of 1 878 727 SNPs distributed along the P. nigra reference sequence were identified. The SNP calling accuracy was validated with Sanger sequencing. SNPs were selected within 14 previously identified QTL regions, 2916 expressional candidate genes related to rust resistance, wood properties, water‐use efficiency and bud phenology and 1732 genes randomly spread across the genome. Over 10 000 SNPs were selected for the construction of a 12k Infinium Bead‐Chip array dedicated to association mapping. The SNP genotyping assay was performed with 888 P. nigra individuals. The genotyping success rate was 91%. Our high success rate was due to the discovery panel design and the stringent parameters applied for SNP calling and selection. In the same set of P. nigra genotypes, linkage disequilibrium throughout the genome decayed on average within 5–7 kb to half of its maximum value. As an application test, ADMIXTURE analysis was performed with a selection of 600 SNPs spread throughout the genome and 706 individuals collected along 12 river basins. The admixture pattern was consistent with genetic diversity revealed by neutral markers and the geographical distribution of the populations. These newly developed SNP resources and genotyping array provide a valuable tool for population genetic studies and identification of QTLs through natural‐population based genetic association studies in P. nigra.     

Shift in precipitation regime promotes interspecific hybridization of introduced Coffea species

The frequency of plant species introductions has increased in a highly connected world, modifying species distribution patterns to include areas outside their natural ranges. These introductions provide the opportunity to gain new insight into the importance of flowering phenology as a component of adaptation to a new environment. Three Coffea species, C. arabica, C. canephora (Robusta), and C. liberica, native to intertropical Africa have been introduced to New Caledonia. On this archipelago, a secondary contact zone has been characterized where these species coexist, persist, and hybridize spontaneously. We investigated the impact of environmental changes undergone by each species following its introduction in New Caledonia on flowering phenology and overcoming reproductive barriers between sister species. We developed species distribution models and compared both environmental envelopes and climatic niches between native and introduced hybrid zones. Flowering phenology was monitored in a population in the hybrid zone along with temperature and precipitation sequences recorded at a nearby weather station. The extent and nature of hybridization events were characterized using chloroplast and nuclear microsatellite markers. The three Coffea species encountered weak environmental suitability compared to their native ranges when introduced to New Caledonia, especially C. arabica and C. canephora. The niche of the New Caledonia hybrid zone was significantly different from all three species' native niches based on identity tests (I Similarity and D Schoener's Similarity Indexes). This area appeared to exhibit intermediate conditions between the native conditions of the three species for temperature‐related variables and divergent conditions for precipitation‐related ones. Flowering pattern in these Coffea species was shown to have a strong genetic component that determined the time between the triggering rain and anthesis (flower opening), specific to each species. However, a precipitation regime different from those in Africa was directly involved in generating partial flowering overlap between species and thus in allowing hybridization and interspecific gene flow. Interspecific hybrids accounted for 4% of the mature individuals in the sympatric population and occurred between each pair of species with various level of introgression. Adaptation to new environmental conditions following introduction of Coffea species to New Caledonia has resulted in a secondary contact between three related species, which would not have happened in their native ranges, leading to hybridization and gene flow.     

Genome rearrangements derived from homoeologous recombination following allopolyploidy speciation in coffee

Allopolyploidization is widespread and has played a major role in flowering plant diversification. Genomic changes are common consequences of allopolyploidization, but their mechanisms of occurrence and dynamics over time are still poorly understood. Coffea arabica, a recently formed allotetraploid, was chosen as a model to investigate genetic changes in allopolyploid using an approach that exploits next‐generation sequencing technologies. Genes affected by putative homoeolog loss were inferred by comparing the numbers of single‐nucleotide polymorphisms detected using RNA‐seq in individual accessions of C. arabica, and between accessions of its two diploid progenitor species for common sequence positions. Their physical locations were investigated and clusters of genes exhibiting homoeolog loss were identified. To validate these results, genome sequencing data were generated from one accession of C. arabica and further analyzed. Genomic rearrangements involving homoeologous exchanges appear to occur in C. arabica and to be a major source of genetic diversity. At least 5% of the C. arabica genes were inferred to have undergone homoeolog loss. The detection of a large number of homoeologous exchange events (HEEs) shared by all accessions of C. arabica strongly reinforces the assumption of a single allopolyploidization event. Furthermore, HEEs were specific to one or a few accessions, suggesting that HEE accumulates gradually. Our results provide evidence for the important role of HEE in allopolyploid genome evolution.     

Molecular Weights of Glutenin- and Gliadin-polypeptides Estimated by SDS-Polyacrylamide Gel Electrophoresis

ABSTRACT

Interspecific single-nucleotide polymorphisms (SNPs) in the rbcL DNA barcode have been strictly validated and adopted as a designed SNP genotyping maker to discriminate between two major coffee species, Coffea arabica and C. canephora, and to estimate the mixing ratio of DNA from C. arabica/C. canephora in this study. The SNP genotyping is applicable to not only green (unroasted) coffee beans, but also processed coffee products (roasted coffee beans and instant coffee powder), in which genomic DNA is degraded, because the genotyping developed in this study requires only 10 copies of 63-bp-long DNA fragments of rbcL gene. The authenticity assay established in this study has several advantages: a high versatility to DNA sample conditions; simple and rapid procedures (only two steps; DNA extraction and SNP genotyping); the feasibility in coffee business for practical use to prevent false advertising and provide quality control.

Abbreviations: SNP: single-nucleotide polymorphism; SBS: single base substitution; ISR: intergenic spacer region; INDEL: insertion-deletion.     

RNA‐Seq bulked segregant analysis enables the identification of high‐resolution genetic markers for breeding in hexaploid wheat

The identification of genetic markers linked to genes of agronomic importance is a major aim of crop research and breeding programmes. Here, we identify markers for Yr15, a major disease resistance gene for wheat yellow rust, using a segregating F₂ population. After phenotyping, we implemented RNA sequencing (RNA‐Seq) of bulked pools to identify single‐nucleotide polymorphisms (SNP) associated with Yr15. Over 27 000 genes with SNPs were identified between the parents, and then classified based on the results from the sequenced bulks. We calculated the bulk frequency ratio (BFR) of SNPs between resistant and susceptible bulks, selecting those showing sixfold enrichment/depletion in the corresponding bulks (BFR > 6). Using additional filtering criteria, we reduced the number of genes with a putative SNP to 175. The 35 SNPs with the highest BFR values were converted into genome‐specific KASP assays using an automated bioinformatics pipeline (PolyMarker) which circumvents the limitations associated with the polyploid wheat genome. Twenty‐eight assays were polymorphic of which 22 (63%) mapped in the same linkage group as Yr15. Using these markers, we mapped Yr15 to a 0.77‐cM interval. The three most closely linked SNPs were tested across varieties and breeding lines representing UK elite germplasm. Two flanking markers were diagnostic in over 99% of lines tested, thus providing a reliable haplotype for marker‐assisted selection in these breeding programmes. Our results demonstrate that the proposed methodology can be applied in polyploid F₂ populations to generate high‐resolution genetic maps across target intervals.