Molecular characterization of cherry varieties.
Abstract
Molecular markers are useful tools for measuring the genetic diversity among agricultural species. In plants, microsatellites are still the most used markers for germplasm characterization, conservation, and traceability purposes, while in the livestock sector, although having represented the standard for at least two decades, they are still used only for minor farm animal species. In this work, together with a review on the use of microsatellites in livestock, we also illustrate the use of these markers for the characterization of agricultural diversity and food traceability through two case studies: (i) the analysis of genetic diversity in ancient fruit tree cultivars of apple (Malus × domestica Borkh.), pear (Pyrus communis L.), sweet cherry (Prunus avium L.), and sour cherry (Prunus cerasus L.) from Northern Italy and (ii) the molecular authentication of wheat food chain. In the former case, a high genetic variability as well as the presence of different ploidy levels were detected, while in the latter microsatellite markers were shown to be useful for traceability and product authentication along the whole food chain. Overall, the presented evidence confirms the versatility of microsatellites as markers for both agrobiodiversity characterization and food traceability in cultivated plants and farm animals.
Keywords
- agrobiodiversity
- fruit tree
- livestock
- microsatellites
- traceability
1. Introduction
Molecular characterization has various purposes in plant and animal genetic resource management, such as elucidating relationships between breeds/varieties, characterizing new genotypes, monitoring shifts in population genetic structure, and exploiting associations among traits and markers [1–3]. A well-recognizable molecular profile is a key factor for the protection and conservation of any genetic resource. Researchers can properly exploit plant and animal genetic resources if the materials are well characterized. Low assay cost, affordable hardware, throughput, convenience and ease of assay development, and automation are important factors when choosing DNA-based technology.
Microsatellites, or simple sequence repeats (SSRs), are polymorphic loci that derive from the repetition of short sequence motives of one to six base pairs in length. Microsatellites are among the most useful markers mainly because they are single locus co-dominant markers [4]. In the plant field, the availability of co-dominant markers is important in the analysis of hybrids. Furthermore, with respect to some categories of multi-locus markers (e.g. RAPD), microsatellites are characterized by higher reproducibility. Microsatellites have been largely used for DNA fingerprinting in several species, both wild and domesticated, although in recent years they have been increasingly replaced by single nucleotide polymorphisms (SNPs), particularly in the livestock genetic field [5].
Microsatellites have a series of characteristics that make them ideal to analyze plant genomes: (1) co-dominance that makes possible the analysis of hybrids of plant commercial varieties; (2) the amplified fragments are usually small in size (100 and 300 base pairs) resulting in positive PCR amplifications even in highly degraded DNA; (3) because of the polyploid nature of the genome of several important crop species, a small number of selected SSRs are able to provide a high discrimination capacity, as reported in the section on plant biodiversity; (4) SSRs are automatable, reproducible between different laboratories (provided that some precautions are taken to uniform allele size scoring, such as sharing of standard samples between labs), easily multiplexed, and easy to score; (5) SSRs usually show a high level of polymorphism and several alleles can be detected for a single SSR locus. This latter aspect makes SSRs extremely useful also for organisms with limited or no information on the genomic sequence because a small number of markers can be enough to clearly discriminate between a large number of samples. Compared to SNP markers, SSRs are less numerous in the genome but present a higher number of alleles per locus (SNPs are usually bi-allelic); therefore, a small number of SSRs can result in a discrimination capacity similar to that obtained with a large number of SNPs [6].
Biological diversity—or biodiversity—is a term used to describe the variety of life on Earth. It refers to the wide variety of ecosystems and living organisms: animals, plants, their habitats, and their genes [7]. While biodiversity can be considered as the foundation of life on Earth, it is crucial for the functioning of ecosystems providing us with products and services without which we could not live. Biodiversity is also the foundation of agriculture. In presence of biodiversity, men can select the genetic material available and gradually improve varieties and breeds. Preservation of biodiversity is, therefore, recognized worldwide as a topic of great concern both in wild and agricultural species, and with respect to the latter, recently, there has been an increasing interest in preserving local plant germplasms. Local varieties as breeds, landraces, ecotypes, and ancient varieties, which have been rarely subjected to breeding, are usually characterized by high genetic variability and genotypes. These germplasm resources are well adapted to both local needs and environmental conditions with good fitness for the anthropic and natural environments in which they have evolved [2, 3].
Local germplasms, as ancient fruit tree cultivars or traditional livestock breeds, frequently face strong genetic erosion starting from the twentieth century. Genetic erosion refers to “the loss of individual genes and the loss of particular combinations of genes (i.e., gene complexes) such as those maintained in locally adapted landraces” [8]. Therefore, the term “genetic erosion” refers to both the loss of genes or alleles and the loss of varieties. Conservation of genetic materials, both using
Microsatellites have been used to evaluate crop germplasm and genetic diversity in several species, including rye [9], grape [10], sugarcane [11], rice [12], and olive [13]. Agrobiodiversity of fruit tree is of increasing concern mainly because repositories still remain a valuable source of allelic variation for many traits and can be exploited for breeding in the near future. Studying the genetic diversity of germplasm resources is not only significant for the protection of species, but also necessary for the development and utilization of germplasm resources for crop improvement and to face existing and future biotic and abiotic constraints with respect to sustainable production in the context of global environmental change [14]. Examples include apple landraces (
In the livestock sector, microsatellite markers have been widely used for more than a decade for the characterization and conservation of livestock biodiversity and for the traceability of food products. In livestock, current genotyping standards are represented by standardized SNP panels that allow the characterization of tens or hundreds of thousand markers per sample [5]; but due to the low costs and to the possibility of in-house implementation of genotyping protocols, microsatellite markers still represent a useful resource to characterize livestock breeds in several developing countries, in which the access to SNP typing or other high throughput technologies can be difficult or too expensive [17–19]. Some years ago, FAO published recommendations for standardized sets of microsatellite loci to be used for studying diversity in the major livestock species [20] in order to make possible the comparison of results across different research projects [17–19].
The good resolution power and frequent occurrence of SSR within plant and farm animal genomes make this type of marker very useful in the food sector also. Food traceability is a milestone of EU food safety policy. The European Commission has agreed to establish a ‘Reference Centre’ to combat food fraud and ensure the “authenticity and integrity” of the EU food supply chain [21]. EU enhances and supports projects related to food safety as the recently approved project Food Integrity, comprising 38 participants from 18 European countries and one from China [22]. Furthermore, the addition of products without prior declaration on the label, besides representing fraud and adulteration, can also bring health risks, in particular to allergic consumers. In recent years, food traceability has become a topical field mainly to prevent fraud, adulteration, and sophistication. A database of food ingredient fraud issues was developed by [23]. The food products more subject to fraud are, in order, olive oil, milk, honey, saffron, orange juice, coffee, apple juice and wine [23]. Most of the processed foods contain very low quality and quantity of DNA, because thermal or chemical treatments determine its degradation. Being microsatellites short repeats of 1–6 nucleotides, they are the most useful markers for DNA recovered from a treated food matrix and combined with
Apart from adulteration and fraudulent procedures, traceability is of great importance to authenticate the quality and integrity of European high value food. A biochemical and genetic approach using microsatellites was useful to discriminate the geographical origin of Italian red wines obtained from Campania region native red grape varieties [10]. Several DNA-based analytical methods have been developed and applied to identify and quantify cereal species and to fingerprint and identify varieties to verify their authenticity [24, 25] developed a microsatellite-based method to verify the presence of the four required durum wheat cultivars in “Altamura” bread, and which are cultivated in a restricted geographical area close to the town of Altamura. Altamura bread, according to its European mark of protected designation of origin (PDO), at least 80% of the total flour used for Altamura bread preparation must derive from the aforementioned traditional durum wheat cultivars used alone or in combination.
In livestock, breed discrimination is useful to detect fraud and to protect and valorize typical productions. Girgentana goat (
The aim of this paper is to highlight the utility of microsatellite markers to study both genetic diversity of domesticated plants and animals and food traceability. Some examples have been provided in the following sections.
2. Agrobiodiversity: the case study of fruit tree species in Northern Italy
Researchers [29] reported that 940 crop plants species are threatened globally and genetic erosion was described in different crop groups, such as cereals and grasses or fruits and nuts [8]. When a species, or the diversity within a species, is lost, the genes important for improving crops are also lost. Preserving local germplasms, landraces, ecotypes, and ancient varieties, means preserving not only our history and culture (such populations represented for centuries an important source of food for local people) but also an extremely useful reserve of genes usable to introduce new characteristics in modern varieties. In order to preserve the local germplasm of ancient fruit tree cultivars, a systematic recovering and characterization of the traditional material of the western part of the Emilia Romagna region was carried out. In this area the tradition of pear (
Selvaticona di Magnano | PC | 253/219 | 101/92 | 140/118 | 187/185 | 141 | |
Mora piacentina | PC | 253/219 | 101/92 | 140/118 | 187/185 | 141 | |
Picaion acc.1 | PC | 238 | 92 | 118 | 185 | 141/130 | |
Picaion acc.2 | PC | 238 | 92 | 118 | 185 | 141/130 | |
Smirne | PC | 253/219 | 92 | 118 | 185/187 | 130/123 | |
Pavesi acc. A | PC | 253/249 | 101/96 | 118 | 197/183 | 128/125 | |
Pavesi acc. C1 | PC | 253/249 | 101/96 | 118 | 197/183 | 128/125 | |
Pavesi acc. C2 | PC | 253/249 | 101/96 | 118 | 197/183 | 128/125 | |
Mori | PC | 221/219 | 92 | 118 | 187/185 | 141/132 | |
Raffaella | PC | 238 | 101 | 118 | 185 | 141 | |
Flamengo acc.A | PC | 238 | 92 | 118 | 185 | 141/130 | |
Flamengo acc.B | PC | 238 | 92 | 118 | 185 | 141/130 | |
Flamengo acc.C | PC | 238 | 92 | 118 | 185 | 141/130 | |
Duroncina della goccia | PC | 221/219 | 92 | 118 | 197/185 | 141 | |
Prima | PC | 253/249 | 96/92 | 118/114 | 185 | 145/130 | |
Mora di Vignola | PC | 221/219 | 101 | 126/114 | 197/185 | 128/123 | |
Giambella | PR | 251/219 | 96/92 | 126/118 | 187/185 | 128/123 | |
Marasca dal peduncolo lungo |
PC | 249/247/221/195 | 96 | 126/112 | 185/179/155 | 130/113 | |
Marasca Villanova | PC | 249/238/225/195 | 92 | 126/112 | 195/185/175/169 | 141/130/123 | |
Marinone I acc. A | PC | 251/225/195 | 105/92 | 140/126/114/112 | 195/185/175/167 | 141/130/123 | |
Marinone II acc. A | PC | 225/221/195 | 92 | 126/118/112 | 193/185/175/167 | 141/130/123 | |
Marinone II acc. C | PC | 251/225/195 | 105/92 | 140/126/114/112 | 195/185/175/167 | 141/130/123 | |
Amarena Piacentina | PC | 249/247/225/195 | 92 | 126/112 | 193/185/175/167 | 141/130/123 | |
Visciola | PC | 225/211/195 | 105/99 | 126/122/118/110 | 197/191/187/171 | 130/123 |
In addition, DNA analysis was carried out using SSR markers in order to obtain a preliminary fingerprint of each sampled accession and to eventually solve controversies of synonyms (different names for a single genotype) and homonyms (a single name for different genotypes). Genetic variability of the samples was evaluated using five SSR markers for each species: EMPA015, EMPA018 [30], UDP97-402 [31], UCDCH17, and UCDCH31 [32] for sweet and sour cherry; GD96, GD100 and GD162 [33] for apple; KA14, KA16 and BGT23b [34] for pear; GD142, GD147 [33] for both apple and pear (Tables 1–3). DNA extraction from young leaves and PCR amplification have been carried out as previously reported [35]. Analysis of PCR products was performed using an ABI Prism 3100 Genetic Analyzer (Applied Biosystem—Thermofisher). Expected heterozygosity and discrimination power were calculated as described in [35], while observed heterozygosity was calculated as the ratio between heterozygous genotypes over the total number of the samples (
GD96 | GD100 | GD147 | GD162 | GD142 | |||
---|---|---|---|---|---|---|---|
Ruggine acc. I | PC | 178/172 | 230/222 | 150/129 | 219/210 | 138/132 | |
Ruggine acc. II | PC | 178/172 | 230/222 | 150/129 | 219/210 | 138/132 | |
Fior d’acacia | PC | 180/172 | 224 | 146/129 | 230 | 140/138 | |
Verdone | PC | 176/172 | 234/224 | 148/135/129 | 228/210 | 144/126 | |
Rustaio | PC | 178/174 | 224 | 135/129 | 228/210 | 131 | |
Rustajò | PC | 176/174/168 | 226/224/222 | 135 | 230/228/210 | 154/144/126 | |
Restajo | PC | 174/150 | 226/219 | 142/135 | 210 | 144/138 | |
Carraia acc. I | PC | 170/168/150 | 234/230/224 | 148/146/135 | 234/230/222/210 | 140/138 | |
Carraia acc. II | PC | 174/150 | 232/230 | 148/135 | 228/222/210 | 140/138 | |
Salame | PC | 172 | 224 | 148/146 | 222/210 | 144 | |
Rosa | PR | 178/174/168 | 230/226/224 | 148/137/135 | 230/219/210 | 148/144/138 | |
Mela Rosa | PR | 187/185/164 | NA | 139 | 226/210 | 144/132/126 | |
Bella di Maggio | PR | 174 | 226 | 127 | 219/210 | 144/140 | |
Cavic | PR | 178/172 | 224 | 148/135 | 230/228 | 144 | |
Seriana | PR | 176/170/168 | 226 | 148/137/129 | 230/226/210 | 148/144/126 | |
Melo Olio | PR | 194/176 | 222 | 142/137 | 234/228 | 152/140 | |
Cucumero | PR | 172 | 224 | 148 | 222/210 | 144 | |
Ghiacciata | PR | 176/172 | 224 | 135/129 | 210 | 132/126 | |
Musona | PR | 178/172 | 234/224 | 135/129 | 230/228 | 144/142 | |
Codaro | PR | 172/166 | 224/222 | 135 | 228/210 | 152/126 |
Lauro acc. I | PC | 213/195 | 129 | 132/120 | 194/176 | 166/158 | |
Lauro acc. II | PC | 213/195 | 129 | 132/120 | 194/176 | 166/158 | |
Limone acc. I | PC | 209 | 129/115 | 118 | 184/178 | 156/152 | |
Limone acc. II | PC | 209 | 129/115 | 118 | 184/178 | 156/152 | |
Limone acc. III | PC | 209 | 129/115 | 118 | 184/178 | 156/152 | |
Rossetto | PC | 193/191 | 147 | 118 | 184 | 158/148 | |
Macagn | PC | 213/195 | 145/129 | 128/118 | 188 | 174/160 | |
Sburdacen | PC | 191 | 129/123 | 128/118 | 184 | 182/180/176/174 | |
Sburdacion acc. I | PC | NA | 129/123 | 122/118 | 222/190/184 | 178/160/156 | |
Sburdacion acc. II | PC | NA | 129/123 | 120/118 | 190/184 | 178/160/156 | |
Coda torta acc. I | PC | 505/488 | 147/129 | 124/118 | 194/176 | 174/172/146 | |
Coda torta acc. II | PC | 505/488 | 147/129 | 124/118 | 194/176 | 174/172/146 | |
Nigrò | PC | NA | 129/125 | 134/128 | 194/184/166 | 174/164/146 | |
Colar | PC | 213 | 147/129 | 126/120/118 | 194/186/184 | 174/160/146 | |
Bianchetto | PC | 543/509 | 129 | 124/120 | 184/180 | 180/158 | |
Nobile acc. I | PR | 213/195 | 129 | 132/120 | 194/176 | 166/158 | |
Nobile acc. II | PR | 213/195 | 129 | 132/120 | 194/176 | 166/158 | |
Butirra Polesine | PR | 235/231 | 145/129 | 138/118 | 190/176 | 166 | |
San Giovanni | PR | 191 | 129/125 | 125/118 | 194/184 | 168/160 | |
San Germano | PR | 209/203 | 147/131 | 118 | 186 | 160/158 | |
San Pietro | PR | 209 | 139/129 | 122/118 | 184/186 | 164/136 | |
Cipolla | PR | 209/193 | 145/129/123 | 132/118 | 186/184/176 | 166/150/136 | |
Bergamotto | PR | 203 | 131/129 | 122/118 | 184 | 160/156 | |
Nigrer | PR | 179 | 131/129 | 126/118 | 184 | 164/148 | |
Carlet | PR | 179 | 139/129 | 128/118 | 194/186 | 148/146 | |
Moscato | PR | 209 | 129 | 124/118 | 184 | 170/160 | |
Spadone | PR | 179 | 151/115 | 136/126/120 | 184/176 | 178/174/164 | |
Ingurien | PR | 169 | 129/125 | 118 | NA | 164/148 | |
Svirgolato | PR | 223/213 | 129/119 | 126/120 | 184/176 | 166/158 | |
Colar | PR | 213 | 147/129 | 120/118 | 194/186/184 | 174/160/146 | |
Pavia | PR | 209/195 | 145/131/123 | 128/118 | 186 | 158/148 | |
Ducale | PR | 209/195 | 129/125 | 118 | 184/176 | 164/136 | |
Butirra Ruggina | PR | 195 | 129/115 | 128/120/118 | NA | 174/166 |
EMPA015 | 10 | 0.881 | 0.861 | |
EMPA018 | 5 | 0.668 | 0.743 | |
UDP97/402 | 7 | 0.753 | 0.712 | |
UCDCH 17 | 13 | 0.815 | 0.712 | |
UCDCH 31 | 8 | 0.775 | 0.854 | |
Average | 8.6 | 0.778 | 0.776 | |
GD96 | 13 | 0.868 | 0.905 | |
GD100 | 8 | 0.788 | 0.867 | |
GD147 | 9 | 0.818 | 0.920 | |
GD162 | 7 | 0.779 | 0.905 | |
GD142 | 10 | 0.839 | 0.915 | |
Average | 9.4 | 0.818 | 0.902 | |
BGT23b | 16 | 0.888 | 0.915 | |
KA16 | 10 | 0.723 | 0.898 | |
GD147 | 11 | 0.778 | 0.894 | |
KA14 | 11 | 0.807 | 0.907 | |
GD 142 | 18 | 0.921 | 0.935 | |
Average | 13.2 | 0.823 | 0.909 |
Sweet cherry (
The microsatellite analysis revealed a different scenario regarding sour, sweet, and duke cherry accessions (Table 1). The number of different alleles detected is reported in Table 4, the average number of alleles is 8.6, the lowest number of alleles is 5 for EMPA018, and the highest is 13 for UCDCH17. The expected heterozygosity ranged between 0.668 (EMPA018) and 0.881 (EMPA015) (Table 4). Based on the frequencies of the different alleles, the probability to obtain a particular genotype by chance was evaluated. Despite the use of a small set of markers, we had very low probability values ranging from 10−6 to 10−9 for diploid varieties and 10−12 to 10−19 for polyploid varieties. The smallest value was obtained for the variety Visciola, this is likely a consequence of its hybrid nature (data not shown). These results confirm what had already been shown in the case of
The three accessions belonging to the sweet cherry cultivar Pavesi have the same molecular profile, indicating that they derived from a unique mother plant. The same could be noted for the accessions of the cultivars Flamengo and Picaion. Two cultivars, namely Mora piacentina and Selvaticona di Magnano, have the same SSR profile. This situation, with all the caution due to the small number of markers used, could be a typical case of synonymy and the two names could be two different local designations for plants anciently derived from the same genetic material and then vegetatively propagated. Concerning sour cherry, the cultivars Marasca and Marasca di Villanova, despite a similar name, had a different genetic profile suggesting that they belong to two different cultivars and they are a case of homonymy. A similar situation was found within the three accessions belonging to Marinone: Marinone I acc. A and Marinone II acc. C had the same profile while Marinone II acc. A was clearly different. Very likely, the first two accessions derived from the same mother plant while the last one had a different origin resulting in a case of homonymy. Comparing the profiles of the different markers in sweet and sour cherries, sweet cherries had a simple profile with the different loci having just one (homozygous) or two (heterozygous) alleles. On the contrary, sour cherries had a more complicated allelic combination and it was common to find, for each marker, the presence of single loci having three or four different alleles. This high number of alleles at the level of the single locus could be a consequence of local duplications of genomic regions or, more likely, of different ploidy levels. In this respect it is reported that sweet cherries are diploids while sour cherries are polyploids (such as tetraploids).
To have a better representation of the relationships among the different accessions analysis, principal component analysis (PCA) was carried out (Figure 4). Two clearly separated groups could be defined: the first including sweet cherry accessions and the second including sour cherry accessions. Among the sour cherry accessions, the one being closest to the sweet cherry group was the variety Visciola. The term Visciola is used to refer to a variety of duke cherry that originated by natural hybridization between a sweet and a sour cherry variety. This hybrid nature can determine the intermediate position of this sample between the sweet and sour cherry groups.
Apple and pear are among the most economically important fruit tree crops of the temperate zones. According to the FAO report on the state of world’s plant genetic resources for food and agriculture, at least 97,500 apple accessions and 1140 pear accessions are present in worldwide
The provinces of Parma and Piacenza have a long tradition of apple and pear cultivation, and a wide diversity of cultivars, well adapted to the local environmental conditions, was grown in this area since ancient times. In apple, as in cherry, the number of alleles highlighted at a single locus in the different samples, ranged from one to four supporting the presence of different ploidy levels (Table 2).
Along with cultivars having just one or two alleles at each locus, such as Ruggine, Fior d’Acacia, and Salame, there were some cultivars with three alleles per locus, such as Seriana, Rosa, and Rustajò. These results supported diploidy and triploidy as the main ploidy levels in local apple germplasm and they agree with what is generally reported in literature concerning apple varieties: most of the apples grown commercially are diploid (2n), although there are many triploid varieties (3n) [42]. The presence of four different alleles, in a single locus, was a rare event and it was found just in a single case (marker GD162, first accession of variety Carraia). While the high number of currently cultivated varieties is diploid or triploid, the presence of tetraploid forms not cultivated but useful for breeding was reported too [43]. It cannot be excluded that after centuries of vegetative propagation, some tetraploid forms could be originated and unintentionally cultivated. The number of different alleles detected by the five SSRs is reported in Table 4; the average value was 9.4, the lowest value was 7 for GD162, while the highest was 13 for GD96. The expected heterozygosity ranged between 0.779 (GD162) and 0.868 (GD96). The probabilities to obtain a particular genotype by chance were very low ranging from 10−7 to 10−10 for diploid varieties and 10−9 to 10−16 for polyploid varieties (data not shown). Also in apple there were cases of homonymy: i) the two accessions of the variety Carraia were clearly different at the genetic level and, very likely, they originated from different mother plants, ii) despite very similar denominations, the varieties Rustaio, Rustajò, and Restajo had different genetic profiles, so they can be effectively considered as different cultivated varieties.
In pear, as in the previous species, it was possible to detect the presence of loci with more than two alleles (Table 3). As for apple and cherry, this evidence suggested the presence of different ploidy levels in the local pear germplasm. Based on the results, diploid varieties were the most diffused followed by triploids. Tetraploidy was rarer, being evidenced just a single time in cultivar Sburdacen with marker GD142. The presence of varieties of pear characterized by different ploidy level, diploids, triploids, and tetraploids was already reported in the literature [44]. The number of different alleles detected by the five SSRs is reported (Table 4): the average value was 13.2, the highest among the three species, while the average expected heterozygosity and discrimination power were similar to the values of apple. The lowest allele number was 10 for KA16 while the highest was 19 for GD142. The expected heterozygosity ranged between 0.723 (KA16) and 0.921 (GD142) (Table 4). Based on the frequencies of the different alleles, we evaluated the probability to obtain any particular genotype. Once again, the probability values were very low ranging from 10−8 to 10−11 for diploid varieties and 10−10 to 10−14 for polyploid varieties (data not shown).
A clear case of synonymy was present concerning the two names Lauro and Nobile. By comparing the genetic profiles, it was possible to see that the different accessions had the same alleles showing that they derived from a common mother ancestor. In this case, the two names are linked to the different provinces, with the name Lauro diffused in the province of Piacenza and the name Nobile in the province of Parma. The three accessions belonging to the variety Limone had the same genetic profile, confirming that they derived from the same mother plant. The same was found for the two accessions of the variety Coda torta. On the contrary, the two accessions of the variety Sburdacion were slightly different, being a case of homonymy. Probably these accessions derived from a common ancestor that encountered some genetic changes (as somatic mutations). Despite the similar names, varieties Nigrò and Nigrer and varieties Butirra Polesine and Butirra ruggina had different genetic profiles and they can be considered as different cultivars. With respect to cherry and apple, in pear a higher frequency of null alleles, i.e. five cases in pear against one case in apple and none in cherry, was observed. To verify this, the amplifications were replicated at least five independent times and the amplicons were always absent. The two accessions of the variety Sburdacion with the marker BGT23b were both characterized by the absence of amplification, supporting close genetic relationships.
This study confirmed the utility of microsatellite markers for biodiversity evaluation and for all conservation actions that can follow the preliminary analysis of genetic variability. Despite the use of a small number of markers, several cases were highlighted: (1) synonymy in sweet cherry (Mora piacentina and Selvaticona di Magnano) and pear (Lauro and Nobile); (2) homonymy inside the Marinone and Marasca (sour cherry), Carraia (apple), and Sburdacion (pear); (3) accessions belonging to the same cultivated variety characterized by high genetic uniformity as a consequence of the derivation from a common ancestor; (4) high biodiversity in the old local germplasm; (5) different levels of ploidy: diploidy in sweet cherry, apple, and pear; triploidy in apple and pear; tetraploidy, rare in apple and pear, and mainly present in sour cherry.
3. Microsatellite markers in the livestock sector
For more than a decade, microsatellites have been one of the most popular types of markers used in the livestock sector for various purposes [45], e.g., the characterization and conservation of diversity [46, 47], the reconstruction of the post-domestication evolutionary history of farm animals [48, 49], parentage testing [50], mapping of quantitative trait loci (QTL) [51, 52] or other causative mutations [53], and traceability of food products [26, 54, 55]. The average number of microsatellite loci used in livestock research varied between 15 and 30 [45], even if a lower number of highly informative loci have been adopted for specific purposes. For example, the International Society for Animal Genetics has established that panels of as few as 12 microsatellite loci have enough resolution for the routine identification of individuals and parentage testing in cattle and horse [56].
A large number of national and international projects aiming at the description of farm animal species diversity have relied on the use of microsatellites. These markers have been used to estimate diversity (both within and between breeds) and genetic admixture even among closely related breeds, usually by means of clustering approaches, principal coordinate analysis, or phylogenetic inference [46]. Comprehensive microsatellite-based studies of livestock diversity have been carried out in European chicken [57], goats from Europe and the middle East [58], Eurasian sheep [59], and African cattle [48], just to mention a few.
One of the major drawbacks of microsatellite genotyping is that the use of different PCR-amplification protocols and genotyping techniques may result in different allele size scoring at the same locus in different labs or experiments, thus hampering the possibility to combine microsatellite genotypes obtained from different projects. To circumvent this, the use of the same set of markers (or at least of a common subset of markers) and genotyping of standard samples across projects has been recommended [60]. In particular, to promote the use of common marker panels, the ISAG-FAO Advisory Group on Animal Genetic Diversity has published guidelines and ranked lists of microsatellite loci to be used for studying diversity in major livestock species [20]. Using these markers in order of ranking should maximize the overlap and increase the possibility of merging data from different investigations.
Concerning allele size standardization through the inclusion of standard samples, for some species (e.g. sheep and goats) the standards adopted in the course of large-scale projects have also been shared with research initiatives in different continents to permit merging of the results. This is the case of the European project Econogene [61] whose sheep and goat standard samples have been made available to other large-scale investigations in Africa and Asia. Acknowledging the usefulness of a joint analysis of different datasets to obtain a global view of livestock diversity, as in the case of the meta-analysis performed by the EU project GlobalDiv [62, 63], a number of statistical methods have been devised that allow merging and analyzing datasets even when they have only a few breeds and/or markers in common. The method developed by [64], for example, estimates population genetics parameters (e.g., heterozygosity, allelic richness, and admixture) by means of a double regression approach and has been successfully applied to the meta-analysis of microsatellite data of cattle populations from Europe, Africa, and Asia [45]. [65], instead, have devised a method based on iterative regression to infer the contribution given by each missing allele/breed combination, which allows calculating genetic distances also on merged datasets with missing information (see [45] and figures therein).
Gaining a global view on the worldwide patterns of diversity of livestock genetic resources may allow to highlight (i) the presence of gaps, i.e., areas in which livestock characterization is incomplete or lacking, (ii) local diversity hotspots which may deserve particular attention or conservation efforts, (iii) geographical trends of clonal variation or discontinuities that can shed further light on the evolutionary history and post-domestication migration routes of farm animal species.
In livestock, current genotyping standards are represented by standardized SNP panels that allow the characterization of tens or hundreds of thousand markers per sample at the same time and at a reasonable cost. Commercial SNP chips at varying levels of marker density are already available for the major livestock species, e.g. for cattle at medium density [66] and high density [67], for sheep and goats at medium density [68, 69]. Being highly standardized, SNP panels do not suffer from allele scoring differences and thus permit an immediate comparison and merging of data produced in different labs [70]. A comparative evaluation of the effectiveness of microsatellites vs. SNP markers for individual identification and parentage assessment has recently shown that 2–3 SNPs per microsatellite were necessary to obtain a comparable exclusion power value in a highly consanguineous Angus cattle herd [71]. Therefore, in a similar context the use of, e.g. 50K SNP chip panel might be equivalent to typing of 16–25K microsatellite loci. Nevertheless, due to the low costs and to the possibility of in-house implementation of genotyping protocols, microsatellite markers still represent a useful resource, e.g. to characterize livestock breeds in several developing countries [72, 73], in which the access to SNP typing or other high throughput technologies can be difficult or just too expensive, or to set priorities for conservation at the local or regional scale [74, 75].
4. Traceability of food
Food traceability is of primary importance to avoid fraudulent procedures and to authenticate the origin of particular products. Dishonest producers may substitute, partially or totally, some food products with others less expensive to increase the profit. For this reason, certifying the origin and composition of a certain food is becoming more and more important [76–78]. Molecular analysis is one of the most recently developed methods to trace food products. Molecular traceability is useful to distinguish traditional varieties with specific high quality traits and to protect the PDO and “Protected Geographical Indication” (PGI) marks. Italian products represent 20% of protected food in Europe and the certified “made in Italy” is important for Italian product exportation. DNA is present in every food product and its analysis makes possible to recover a lot of information about the identity of the ingredients in foods and feed. It is often reported that DNA is relatively more resistant than other classes of biological molecules (e.g. proteins) to the degradation caused by food processing. Despite this, as a consequence of processes such as cooking, fermentation etc., degradation of DNA occurs anyway and, generally, the stronger the treatment the shorter the DNA fragments become. Thus, the possibility to analyze small DNA fragments is very important for traceability purposes.
An additional problem, when working with plant-derived products is that along with the DNA, a high number of different inhibitors of polymerase reactions can be recovered from a food matrix. Plants are very rich in carbohydrates and polyphenols, which tend to be co-extracted with the DNA. Their presence can prevent the activity of polymerases hindering the analysis of DNA by PCR reaction. Different commercial kits or customized protocols can be considered to tackle this problem and usually DNA extracted from most food matrixes can be analyzed using molecular tools. Molecular markers make it possible to discriminate, not only the species from which the food is originated, but also the variety (cultivar) or population of origin [79–81]. Among the different classes of markers, some are more suitable than others for traceability purposes. Recently, the two main classes of markers that have been adopted are SNPs and microsatellites. While SNPs are becoming the most used markers for animal-based product analysis and identification, microsatellites are still the election markers for genetic traceability of plant-based products.
The final goal of DNA analysis in the agro-food sector is the comparison of the molecular profile of a sample with a reference profile to evidence the presence of congruencies or discrepancies. When the SSR profile of the sample is congruent with what is expected (similar to the reference profile), the two profiles are matching and it is possible to speculate that the sample under investigation has the same origin as the reference. However, in any final conclusion that is reached in certain cases, it is also important to evaluate the probabilities that the two profiles are identical because they derive from the same genetic material and not just by chance. This requires deep knowledge of the genetic base of the species under investigation and the probability level to obtain the same marker profile, using a set of SSRs, in two independent samples just by chance. This is very important for plant species in which it is often not enough to detect the presence of a particular species in a processed product. For several plant-derived products, as for extra virgin olive oil and wine, the final price on the market is highly dependent on the cultivated variety of the species that has been used as raw material. In this situation, a possible fraud could be represented by the substitution of a declared cultivar with another one with a smaller commercial value but with similar organoleptic properties (different cultivars of olive or of grapevine).
180 | 97 | 142/96 | 129/113 | 127/150 | 167/165 | 160/129 | |
180 | 97 | 142/96 | 113 | 127 | 154/152 | 160 | |
180 | 97 | 142/96 | 117 | 127 | 154/152 | 129 | |
180 | 97 | 142/96 | 117 | 127 | 154/152 | 129 | |
180 | 97 | 142/96 | 129/113 | 127/250 | 167/165 | 160/129 | |
180 | 97 | 142/96 | 113 | 127 | 154/152 | 160 | |
180 | 97 | 142/96 | 117 | 127 | 154/152 | 129 | |
180 | 97 | 142/96 | 117 | 127/150 | 154/152 | 129 | |
180 | 97 | 142/96 | 113 | 127 | 154/152 | 160 | |
180 | 97 | 142/96 | 113 | 127 | 154/152 | 160/129 | |
180 | 97 | 142/96 | 129/113 | 127/150 | 167/165 | 160/129 | |
180 | 97 | 142/96 | 113 | 127 | 154/152 | 160 |
Correct identification and authentication of processed food is more challenging than that of fresh food mainly because of the presence of inhibitors and of DNA degradation. To face these problems, PCRs for food traceability are usually low template-DNA PCRs (LT-DNA PCRs), because increasing the amount of DNA may consequently increase the quantity of inhibitors and determine the failure of the amplification. These PCRs are usually carried out using very small amount of DNA (in the order of few dozens of picograms) and high numbers of amplification cycles (> 35) to have a visible signal. While it is reported that PCR can theoretically work even with amounts of template DNA lower than the aforementioned ones, usually LT-DNA PCRs suffer from several limitations. Concerning SSRs, LT-DNA PCRs can be characterized by marker profiles showing a higher heterozygote peak imbalance between the signals of the observed alleles in a specific sample with respect to standard PCR or by the stochastic disappearance of some allele signals (allelic drop-out, mainly a problem for the bigger size alleles). This outcome is mainly a consequence of the small amount and of the degradation of the template DNA. In these conditions, the final result of the PCR can be strongly influenced by the effect of a random selection of the template molecules during the first cycles of the amplification. Other factors that can make the interpretation of the molecular profiles difficult are the presence of: (1) stutter bands; (2) split peaks, deriving from the incomplete adenylation of the PCR products; (3) allelic drop-in, deriving often from contamination and mainly present in the multiplexing amplifications; (4) triploid profile, deriving from the unexpected amplification of three peaks (three loci) from a diploid genome.
In recent years, our laboratory dealt with the extraction and analysis of DNA from different kinds of food matrices with different purposes and different markers technologies [82–86]. In this section, as an example, the results on traceability of wheat-derived products will be provided. These SSR analyses were carried out as a work under contract for which a third party commissioned us. The samples were collected from the whole supply chain of durum wheat (
From the results obtained, it was not possible to find correspondence between the different samples within each label. As an example, seeds A did not correspond to treated seeds A, flour A, and pasta A. On the contrary, seeds A had the same profile as treated seed B and pasta B. Similarly, seeds B had the same profile as treated seeds C, flour C, and pasta C (Figure 5). Concerning the last samples, the presence of correspondence between seeds C, treated seeds A, and flour A was evidenced. Absence of correspondence was found for type A pasta whose genetic profile was more similar to the genetic profile of pasta B and for flour B whose genetic profile was unique and different from the other profiles. As previously stated, samples were received in blind without any knowledge about the origin of the different labeled samples.
Based on this, it was possible to conclude that the seeds of cultivar B (the exact name of the variety was unknown) were used to produce treated seeds C, flour C, and pasta C; seeds of cultivar A were used to produce treated seeds B and pasta B; seeds of cultivar C were used to produce treated seeds A and flour A (Figure 5). Pasta A was likely produced by mixing flour A with flour C in almost identical percentages and this was explained by the appearance of the signal corresponding to flour C allele (Figure 5). The only incongruence was about flour B. This sample had a genetic profile different from the other samples: it had the same profile of flour A with just an extra allele with SSR Xgwm577. This means that flour B was obtained from a fourth and different cultivated variety, but the possibility of contamination cannot be excluded. Concerning the sample pasta B, the amplification with marker Xgm459 was replicated four times and two times just the 113 bp allele was obtained, while the other two times both the 129 and 113 bp alleles were amplified. As reported previously, working with food-derived DNA is challenging also because of the allelic drop-out: the stochastic disappearance of one of the alleles, usually the biggest one, can be observed as a consequence of DNA degradation, which can explain the results obtained for pasta B.
The results obtained were a clear indication of the utility of SSR markers in following the whole wheat chain, despite the DNA degradation determined by processing.
5. Conclusions
The recent development of high throughput genotyping methods has prompted SNPs as desired markers for several applications in agricultural research, in particular in the livestock sector. Despite this, microsatellites, because of their characteristics, can still be considered as markers of choice for numerous studies, in particular concerning plant genomes, both for biodiversity studies and for molecular traceability of plant-derived food products. In a biodiversity study of local ancient germplasm of fruit tree species, using a small number of markers, we obtained important indications as the presence of synonymy and homonymy, high biodiversity, and different levels of ploidy. Furthermore, the high polymorphism of microsatellite loci together with the different ploidy levels detected increased the probability to link each cultivar to its corresponding genotypic profile. This is particularly interesting because it means that few properly selected SSRs can be enough to obtain robust results. In the same time, microsatellites can be very useful for molecular traceability as it was evidenced from our results of the whole production chain from durum wheat raw material to processed pasta. Indeed, despite the degradation of DNA caused by food processing, SSRs were able to find the correspondence between blind samples and genotypes highlighting some incongruences.
Acknowledgments
The study of fruit tree biodiversity was funded by Misura 214, and Azione 7 of the Emilia Romagna regional development plan.
References
- 1.
Ahrens CW, James EA. Conserving the small milkwort, Comesperma polygaloides , a vulnerable subshrub in a fragmented landscape. Conserv Genet. 2016;17:891–901. DOI 10.1007/s10592-016-0830-9 - 2.
Martins S, Simões F, MendonÇa D, Matos J, Silva AP, Carnide V. Western European wild and landraces hazelnuts evaluated by SSR markers. Plant Mol Biol Report. 2015; 33 :1712–1720. DOI 10.1007/s11105-015-0867-9 - 3.
Liang W, Dondini L, De Franceschi P, Paris R, Sansavini S, Tartarini S. Genetic diversity, population structure and construction of a core collection of apple cultivars from Italian Germplasm. Plant Mol Biol Report. 2015; 33 :458–473. DOI 10.1007/s11105-014-0754-9 - 4.
Powell W, Morgante M, Andre C, Hanafey M, Vogel J, Tingey S, Rafalski A. The comparison of RFLP, RAPD, AFLP and SSR (microsatellite) markers for germplasm analysis. Mol Breed. 1996; 2 :225–238. - 5.
Nicoloso L, Bomba L, Colli L, Negrini R, Milanesi M, Mazza R, and the Italian Goat Consortium. Genetic diversity of Italian goat breeds assessed with a medium-density SNP chip. Genet Sel Evol. 2015; 47 :62. DOI 10.1186/s12711-015-0140-6 - 6.
Scarano D, Rao R. DNA markers for food products authentication. Diversity. 2014; 6 :579–596. DOI:10.3390/d6030579 - 7.
IUCN, International Union for Conservation of Nature [Internet]. 2016. Available from: http://www.iucn.org [Accessed: 2016-07-05] - 8.
The Second Report on the state of the World’s Plant Genetic Resources for Food and Agriculture [Internet]. 2010. Available from: http://www.fao.org/agriculture/crops/thematic-sitemap/theme/seeds-pgr/sow/sow2/en/ [Accessed: 2016-07-05] - 9.
Targońska M, Bolibok-Brągoszewska H, Rakoczy-Trojanowska M. Assessment of genetic diversity in Secale cereale based on SSR markers. Plant Mol Biol Rep. 2016;34 :37–51. DOI: 10.1007/s11105-015-0896-4 - 10.
Muccillo L, Gambuti A, Frusciante L, Iorizzo M, Moio L, Raieta K, Rinaldi A, Colantuoni V, Aversano R. Biochemical features of native red wines and genetic diversity of the corresponding grape varieties from Campania region. Food Chem. 2014; 143 :506–513. DOI: 10.1016/j.foodchem.2013.07.133 - 11.
Lu X, Zhou H, Pan Y, Chen C, Zhu J, Chen P, Li Y, Cai Q, Chen R. Segregation analysis of microsatellite (SSR) markers in sugarcane polyploids. Genet Mol Res. 2015; 14 :18384–18395. DOI: 10.4238/2015.December.23.26 - 12.
Ahmad F, Hanafi MM, Hakim MA, Rafii MY, Arolu IW, Abdullah SNA. Genetic divergence and heritability of 42 coloured upland rice genotypes ( Oryza sativa ) as revealed by microsatellites marker and agro-morphological traits. PLoS One. 2015;10 . DOI: 10.1371/journal.pone.0138246 - 13.
Doveri S, Gil F, Díaz A, Reale S, Busconi M, Machado A, Martín A, Fogher C, Donini P, Lee D. Standardization of a set of microsatellite markers for use in cultivar identification studies in olive ( Olea europaea L.). Sci Hortic. 2008;116 :367–373. DOI:10.1016/j.scienta.2008.02.005 - 14.
Dwivedi SL, Ceccarelli S, Blair MW, Upadhyaya HD, Are AK, Ortiz R. Landrace germplasm for improving yield and abiotic stress adaptation. Trends Plant Sci. 2016; 21 :31–42. DOI: 10.1016/j.tplants.2015.10.012 - 15.
Lassois L, Denancé C, Ravon E, Guyader A, Guisnel R, Hibrand-Saint-Oyant L, Poncet C, Lasserre-Zuber P, Feugey L, Durel C-E. Genetic diversity, population structure, parentage analysis, and construction of core collections in the French apple germplasm based on SSR markers. Plant Mol Biol Rep. 2016:1-18. DOI: 10.1007/s11105-015-0966-7 - 16.
Gao Y, Liu F, Wang K, Wang D, Gong X, Liu L, et al. Genetic diversity of Malus cultivars and wild relatives in the Chinese National Repository of Apple Germplasm Resources. Tree Genet Genom. 2015; 11 . DOI: 10.1007/s11295-015-0913-7 - 17.
Zaman G, Shekar MC, Aziz A. Molecular characterization of Meghalaya Local pigs (Niang Megha) using microsatellite markers. Indian J Sci Technol. 2013; 6 :5302-5306. DOI: 10.17485/ijst/2013/v6i10/38777 - 18.
Colombo E, Strillacci MG, Cozzi MC, Madeddu M, Mangiagalli MG, Mosca F, Zaniboni L, Bagnato A, Cerolini S. Feasibility study on the FAO chicken microsatellite panel to assess genetic variability in the Turkey ( Meleagris gallopavo ). Ital J Anim Sci. 2014;13 :887–890. DOI: 10.4081/ijas.2014.3334 - 19.
Radhika G, Raghavan KC, Aravindakshan TV, Thirupathy V. Genetic diversity and population structure analysis of native and crossbred goat genetic groups of Kerala, India. Small Ruminant Res. 2015; 131 :50–57. DOI: 10.1016/j.smallrumres.2015.08.008 - 20.
DADÂIS Domestic Animal Diversity Information System (DADÂIS), Food and Agriculture Organization of the United Nations [Internet]. 2016. Available from: http://www.fao.org/dad-is/ [Accessed: 2016-07-05] - 21.
Traceability - EU Food Law [internet]. 2016. Available from: http://www.eurofoodlaw.com/food-safety-and-standards/traceability [Accessed: 2016-07-05] - 22.
Food Integrity [Internet] Available from: https://secure.fera.defra.gov.uk/foodintegrity [Accessed: 2016-07-05] - 23.
Moore J, Spink J, Lipp M. Development and application of a database of food ingredient fraud and economically motivated adulteration from 1980 to 2010. J Food Sci. 2012; 77 :R118–R126. DOI: 10.1111/j.1750-3841.2012.02657.x - 24.
Terzi V, Morcia C, Gorrini A, Stanca AM, Shewry PR, Faccioli P. DNA-based methods for identification and quantification of small grain cereal mixtures and fingerprinting of varieties. J Cereal Sci. 2005; 41 :213–220. DOI: 10.1016/j.jcs.2004.08.003 - 25.
Pasqualone A, Alba V, Mangini G, Blanco A, Montemurro C. Durum wheat cultivar traceability in PDO Altamura bread by analysis of DNA microsatellites. Eur Food Res Technol. 2010; 230 :723–729. DOI: 10.1007/s00217-009-1210-1 - 26.
Sardina M, Tortorici L, Mastrangelo S, Gerlando R, Tolone M, Portolano B. Application of microsatellite markers as potential tools for traceability of Girgentana goat breed dairy products. Food Res Int. 2015; 74 :115–122. DOI: 10.1016/j.foodres.2015.04.038 - 27.
Mateus JC, Russo-Almeida PA. Traceability of 9 Portuguese cattle breeds with PDO products in the market using microsatellites. Food Control. 2015; 47 :487–492. DOI: 10.1016/j.foodcont.2014.07.038 - 28.
Dalvit C, De Marchi M, Dal Zotto R, Gervaso M, Meuwissen T, Cassandro M. Breed assignment test in four Italian beef cattle breeds. Meat Sci. 2008; 80 :389–395. DOI: 10.1016/j.meatsci.2008.01.001 - 29.
Khoshbakht K, Hammer K. Threatened and rare ornamental plants. J Agr Rural Dev Trop Subtrop. 2007; 108 :19–39. - 30.
Clarke JB, Tobutt KR. Development and characterization of polymorphic microsatellites from Prunus avium ‘Napoleon’. Mol Ecol Notes. 2003;3 :578–580. DOI: 10.1046/j.1471-8286.2003.00517.x - 31.
Cipriani G, Lot G, Huang W-G, Marrazzo M T, Peterlunger E, Testolin R. AC/GT and AG/CT microsatellite repeats in peach [ Prunus persica (L) Batsch]: isolation, characterisation and cross-species amplification in Prunus. Theor Appl Genet. 1999;99 :65–72. DOI: 10.1007/s001220051209 - 32.
Turkoglu Z, Bilgener S, Ercisli S, Bakır M, Koc A, Akbulut M, Gercekcioglu R, Gunes M, Esitken A. Simple sequence repeat-based assessment of genetic relationships among Prunus rootstocks. Genet Mol Res. 2010;9 :2156–2165. DOI: 10.4238/vol9-4gmr957 - 33.
Hokanson S C, Szewc-McFadden A K, Lamboy W F, McFerson J R. Microsatellite (SSR) markers reveal genetic identities, genetic diversity and relationships in a Malus ×domestica Borkh: core subset collection. Theor Appl Genet. 1998;97 :671–683. DOI: 10.1007/s001220050943 - 34.
Yamamoto T, Kimura T, Sawamura Y, Manabe T, Kotobuki K, Hayashi T, Ban Y, Matsuta N. Simple sequence repeats for genetic analysis in pear. Euphytica. 2002; 124 :129–137. DOI: 10.1023/A:1015677505602 - 35.
Martinelli F, Busconi M, Camangi F, Fogher C, Stefani A, Sebastiani L. Ancient Pomoideae (Malus domestica Borkh. andPyrus communis L.) cultivars in “Appennino Toscano” (Tuscany, Italy): molecular (SSR) and morphological characterization. Caryologia. 2008;61 :320–331. DOI: 10.1080/00087114.2008.10589643 - 36.
Wünsch A, Hormaza JI. Molecular characterisation of sweet cherry ( Prunus avium L.) genotypes using peach [Prunus persica (L.) Batsch] SSR sequences. Heredity. 2002;89 :56–63. DOI: 10.1038/sj.hdy.6800101 - 37.
Höfer M, Peil A. Phenotypic and genotypic characterization in the collection of sour and duke cherries ( Prunus cerasus and ×P. ×gondouini ) of the Fruit Genebank in Dresden-Pillnitz, Germany. Genet Resour Crop Evol. 2015;62 :551–566. DOI: 10.1007/s10722-014-0180-8 - 38.
Pereira L, Martins-Lopes P, Batista C, Zanol GC, Clímaco P, Brazão J, Eiras-Dias JE, Guedes-Pinto H. Molecular markers for assessing must varietal origin. Food Anal Method. 2012; 5 :1252–1259. DOI: 10.1007/s12161-012-9369-7 - 39.
Velasco R, Zharkikh A, Affourtit J, Dhingra A, Cestaro A, Kalyanaraman A, Fontana P, Bhatnagar SK, Troggio M, Pruss D et al. The genome of the domesticated apple ( Malus ×domestica Borkh.). Nat Genet. 2010;42 :833–839. DOI: 10.1038/ng.654 - 40.
Cornille A, Gladieux P, Smulders M, Roldán-Ruiz I, Laurens F, Le Cam B, Nersesyan A, Clavel J, Olonova M, Feugey L, Gabrielyan I, Zhang X, Tenaillon M, Giraud T. New insight into the history of domesticated apple: secondary contribution of the European wild apple to the genome of cultivated varieties. PLoS Genet. 2012; 8 :e1002703. DOI: 10.1371/journal.pgen.1002703 - 41.
Silva GJ, Souza TM, Barbieri RL, Costa de Oliveira A. Origin, domestication, and dispersing of Pear ( Pyrus spp.). Adv Agr. 2014; 2014:541097. DOI:10.1155/2014/541097 - 42.
Höfer M, Meister A. Genome size variation in Malus species. J Bot. 2010; 2010:480873. DOI:10.1155/2010/480873. - 43.
Sedysheva GA, Gorbacheva NG. Estimation of new tetraploid apple forms as donors of diploid gametes for selection on a polyploidy level. Univ J Plant Sci. 2013; 1 :49–54. DOI: 10.13189/ujps.2013.010204 - 44.
Cao Y, Huang L, Li S, Yang Y. Genetics of ploidy and hybridized combination types for polyploid breeding in pear. Acta Hortic. 2002; 587 :207–210. DOI: 10.17660/ActaHortic.2002.587.24 - 45.
Lenstra JA, Groeneveld LF, Eding H, Kantanen J, Williams JL, Taberlet P, Nicolazzi EL, Sölkner J, Simianer H, Ciani E, Garcia JF, Bruford MW, Ajmone-Marsan P, Weigend S. Molecular tools and analytical approaches for the characterization of farm animal genetic diversity. Anim Genet. 2012; 43 :483–502. DOI: 10.1111/j.1365-2052.2011.02309.x - 46.
Groeneveld LF, Lenstra JA, Eding H, Toro MA, Scherf B, Pilling D, Negrini R, Finlay EK, Jianlin H, Groeneveld E, Weigend S; GLOBALDIV Consortium. Genetic diversity in farm animals: a review. Anim Genet. 2010; 41 S1 :6–31. DOI: 10.1111/j.1365-2052.2010.02038.x - 47.
Colli L, Perrotta G, Negrini R, Bomba L, Bigi D, Zambonelli P, Verini Supplizi A, Liotta L, Ajmone-Marsan P. Detecting population structure and recent demographic history in endangered livestock breeds: the case of the Italian autochthonous donkeys. Anim Genet. 2013; 44 :69–78. DOI: 10.1111/j.1365-2052.2012.02356.x - 48.
Hanotte O, Bradley DG, Ochieng JW, Verjee Y, Hill EW, Rege JE. African pastoralism: genetic imprints of origins and migrations. Science. 2002; 296 :336–339. DOI: 10.1126/science.1069878 - 49.
Cymbron T, Freeman AR, Isabel Malheiro M, Vigne JD, Bradley DG. Microsatellite diversity suggests different histories for Mediterranean and Northern European cattle populations. Proc Biol Sci. 2005; 272 :1837–1843. DOI: 10.1098/rspb.2005.3138 - 50.
da Silva EC, McManus CM, de Paiva Guimarães MP, Gouveia AM, Facó O, Pimentel DM, Caetano AR, Paiva SR. Validation of a microsatellite panel for parentage testing of locally adapted and commercial goats in Brazil. Genet Mol Biol. 2014; 37 :54–60. DOI: 10.1590/S1415-47572014000100010 - 51.
Uemoto Y, Sato S, Ohnishi C, Hirose K, Kameyama K, Fukawa K, Kudo O, Kobayashi E. Quantitative trait loci for leg weakness traits in a Landrace purebred population. Anim Sci J. 2010; 81 :28–33. DOI: 10.1111/j.1740-0929.2009.00713.x - 52.
Dayo GK, Gautier M, Berthier D, Poivey JP, Sidibe I, Bengaly Z, Eggen A, Boichard D, Thevenon S. Association studies in QTL regions linked to bovine trypanotolerance in a West African crossbred population. Anim Genet. 2012; 43 :123–132. DOI: 10.1111/j.1365-2052.2011.02227.x - 53.
Georges M, A B Dietz, A Mishra, D Nielsen, L S Sargeant, A Sorensen, M R Steele, X Zhao, H Leipold, J E Womack. Microsatellite mapping of the gene causing weaver disease in cattle will allow the study of an associated quantitative trait locus. Proc Natl Acad Sci U S A. 1993; 90 :1058–1062. DOI: 10.1073/pnas.90.3.1058 - 54.
Lenstra JA. Primary identification: DNA markers for animal and plant traceability. In: Smith I, Furness T (Eds.) Improving Traceability in Food Processing and Distribution. Woodhead Publ., Cambridge; 2005. p. 147–164. - 55.
Orrú L, Napolitano F, Catillo G, Moioli B. Meat molecular traceability: how to choose the best set of microsatellites? Meat Sci. 2006; 72 :312–317. DOI: 10.1016/j.meatsci.2005.07.018 - 56.
International Society of Animal Genetics. [Internet]. 2006. Available from http://www.isag.us/index.asp?autotry=true&ULnotkn=true [Accessed: 2016-07-05] - 57.
Berthouly C, Bed’Hom B, Tixier-Boichard M, Chen CF, Lee YP, Laloë D, Legros H, Verrier E, Rognon X. Using molecular markers and multivariate methods to study the genetic diversity of local European and Asian chicken breeds. Anim Genet. 2008; 39 :121–129. DOI: 10.1111/j.1365-2052.2008.01703.x - 58.
Cañón J, García D, García-Atance MA, Obexer-Ruff G, Lenstra JA, Ajmone-Marsan P, Dunner S; ECONOGENE Consortium. Geographical partitioning of goat diversity in Europe and the Middle East. Anim Genet. 2006; 37 :327–334. DOI: 10.1111/j.1365-2052.2006.01461.x - 59.
Tapio M, Ozerov M, Tapio I, Toro MA, Marzanov N, Cinkulov M, Goncharenko G, Kiselyova T, Murawski M, Kantanen J. Microsatellite-based genetic diversity and population structure of domestic sheep in northern Eurasia. BMC Genet. 2010; 11 :76. DOI: 10.1186/1471-2156-11-76 - 60.
The State of the World’s Animal Genetic Resources for Food and Agriculture. In: Rischkowsky B, Pilling D (Eds.). FAO 2007 Rome. [Internet]. 2007. http://www.fao.org/docrep/010/a1250e/a1250e00.htm [Accessed: 2016-07-05] - 61.
Sustainable conservation of animal genetic resources in margin rural areas: integrating molecular genetics socio-economic and geostatistical approaches. [Internet]. 2002. Available from http://www.econogene.eu [Accessed: 2016-07-05] - 62.
The Second Report on the State of the World’s Animal Genetic Resources for Food and Agriculture. In Scherf BD, Pilling D (Eds.). FAO Commission on Genetic Resources for Food and Agriculture Assessments. [Internet]. 2015. Available from http://www.fao.org/3/a-i4787e/index.html [Accessed: 2016-07-05] - 63.
Ajmone-Marsan P and The GLOBALDIV Consortium: a global view of livestock biodiversity and conservation – GLOBALDIV. Anim Genet. 2010; 41 :1–5. DOI: 10.1111/j.1365-2052.2010.02036.x - 64.
Freeman AR, Bradley DG, Nagda S, Gibson JP, Hanotte O. Combination of multiple microsatellite data sets to investigate genetic diversity and admixture of domestic cattle. Anim Genet. 2006; 37 :1–9. DOI: 10.1111/j.1365-2052.2005.01363.x - 65.
Taeubert H, Bradley D, Simianer H. Estimation of genetic distances from two partly overlapping microsatellite marker data sets. Proceedings of 30th International Conference on Animal Genetics, 20–25 August, 2006; Porto Seguro, Brazil. p.28 - 66.
Matukumalli LK, Lawley CT, Schnabel RD, Taylor JF, Allan MF, Heaton MP, O’Connell J, Moore SS, Smith TP, Sonstegard TS, Van Tassell CP. Development and characterization of a high density SNP genotyping assay for cattle. PLoS One. 2009; 4 :e5350. DOI: 10.1371/journal.pone.0005350 - 67.
Cañas-Álvarez JJ, González-Rodríguez A, Munilla S, Varona L, Díaz C, Baro JA, Altarriba J, Molina A, Piedrafita J. Genetic diversity and divergence among Spanish beef cattle breeds assessed by a bovine high-density SNP chip. J Anim Sci. 2015; 93 :5164–5174. DOI: 10.2527/jas.2015-9271 - 68.
Kijas JW, Lenstra JA, Hayes B, Boitard S, Porto Neto LR, San Cristobal M, Servin B, McCulloch R, Whan V, Gietzen K, Paiva S, Barendse W, Ciani E, Raadsma H, McEwan J, Dalrymple B; International Sheep Genomics Consortium Members. Genome-wide analysis of the world’s sheep breeds reveals high levels of historic mixture and strong recent selection. PLoS Biol. 2013; 10 :e1001258. DOI: 10.1371/journal.pbio.1001258 - 69.
Tosser-Klopp G, Bardou P, Bouchez O, Cabau C, Crooijmans R, Dong Y, Donnadieu-Tonon C, Eggen A, Heuven HC, Jamli S, Jiken AJ, Klopp C, Lawley CT, McEwan J, Martin P, Moreno CR, Mulsant P, Nabihoudine I, Pailhoux E, Palhière I, Rupp R, Sarry J, Sayre BL, Tircazes A, Jun Wang, Wang W, Zhang W; International Goat Genome Consortium. Design and characterization of a 52K SNP chip for goats. PLoS One. 2014; 9 :e86227. DOI: 10.1371/journal.pone.0086227 - 70.
Bruford MW, Ginja C, Hoffmann I, Joost S, Orozco-terWengel P, Alberto FJ, Amaral AJ, Barbato M, Biscarini F, Colli L, Costa M, Curik I, Duruz S, Ferenčaković M, Fischer D, Fitak R, Groeneveld LF, Hall SJG, Hanotte O, Hassan F, Helsen P, Iacolina L, Kantanen J, Leempoel K, Lenstra JA, Ajmone-Marsan P, Masembe C, Megens H-J, Miele M, Neuditschko M, Nicolazzi EL, Pompanon F, Roosen J, Sevane N, Smetko A, Štambuk A, Streeter I, Stucki S, Supakorn C, Telo Da Gama L, Tixier-Boichard M, Wegmann D, Zhan X. Prospects and challenges for the conservation of farm animal genomic resources, 2015-2025. Front Genet. 2015; 6 :314. DOI: 10.3389/fgene.2015.00314 - 71.
Fernández ME, Goszczynski DE, Lirón JP, Villegas-Castagnasso EE, Carino MH, Ripoli MV, Rogberg-Muñoz A, Posik DM, Peral-García P, Giovambattista G. Comparison of the effectiveness of microsatellites and SNP panels for genetic identification, traceability and assessment of parentage in an inbred Angus herd. Genet Mol Biol. 2013; 36 :185–191. DOI: 10.1590/S1415-47572013000200008 - 72.
Azam A, Babar ME, Firyal S, Anjum AA, Akhtar N, Asif M, Hussain T. DNA typing of Pakistani cattle breeds Tharparkar and Red Sindhi by microsatellite markers. Mol Biol Rep. 2012; 39 :845–849. DOI: 10.1007/s11033-011-0807-1 - 73.
Yadav AS, Gahlot K, Gahlot GC, Asraf M, Yadav ML. Microsatellite DNA typing for assessment of genetic variability in Marwari breed of Indian goat. Vet World. 2015; 8 :848–854. DOI: 10.14202/vetworld.2015.848-854. - 74.
Medugorac I, Veit-Kensch CE, Ramljak J, Brka M, Marković B, Stojanovic S, Bytyqi H, Kochoski L, Kume K, Grünenfelder HP, Bennewitz J, Förster M. Conservation priorities of genetic diversity in domesticated metapopulations: a study in taurine cattle breeds. Ecol Evol. 2011; 1 :408–420. DOI: 10.1002/ece3.39 - 75.
Ginja C, Gama LT, Cortes O, Delgado JV, Dunner S, García D, Landi V, Martin-Burriel I, Martinez-Martinez A, Penedo MCT, Rodellar C, Zaragoza P, Canon J and Biobovis Consortium. Analysis of conservation priorities of Iberoamerican cattle based on autosomal microsatellite markers. Genet Sel Evol. 2013; 45 :35. DOI: 10.1186/1297-9686-45-35 - 76.
Testolin R; Lain O. DNA extraction from olive oil and PCR amplification of microsatellite markers. Food Chem Toxicol. 2005; 70 :108–112. DOI: 10.1111/j.1365-2621.2005.tb09011.x - 77.
Alba V, Sabetta W, Blanco A, Pasqualone A, Montemurro C. Microsatellite marker to identify specific alleles in DNA extracted from monovarietal virgin olive oils. Eur Food Res Technol. 2009; 229 :375–382. DOI: 10.1007/s00217-009-1062-8 - 78.
Soffritti G, Busconi M, Sánchez RA, Thiercelin JM, Polissiou M, Roldán M, Fernández JA. Genetic and epigenetic approaches for the possible detection of adulteration and auto-adulteration in Saffron ( Crocus sativus L.) Spice. Molecules. 2016;21 :E343. DOI: 10.3390/molecules21030343 - 79.
Caramante M, Corrado G, Monti LM, Rao R. Simple sequence repeats are able to trace tomato cultivars in tomato food chains. Food Cont. 2011; 22 (3–4):549–554. DOI: 10.1016/j.foodcont.2010.10.002 - 80.
Corrado G, Imperato A, la Mura M, Perri E, Rao R. Genetic diversity among olive varieties of southern Italy and the traceability of olive oil using SSR markers. J Hortic Sci Biotech. 2011; 86 :461–466. DOI: 10.1080/14620316.2011.11512789 - 81.
Pasqualone A, Montemurro C, Summo C, Sabetta W, Caponio F, Blanco A. Effectiveness of microsatellite DNA markers in checking the identity of protected designation of origin extra virgin olive oil. J Agric Food Chem. 2007; 55 :3857–3862. DOI: 10.1021/jf063708r - 82.
Busconi M, Foroni C, Corradi M, Bongiorni C, Cattapan F, Fogher C. DNA extraction from olive oil and its use in the identification of the production cultivar. Food Chem. 2003; 83 :127–134. DOI: 10.1016/S0308-8146(03)00218-8 - 83.
Pafundo S, Busconi M, Agrimonti C, Fogher C, Marmiroli N. Storage-time effect on olive oil DNA assessed by amplified fragments length polymorphisms. Food Chem. 2010; 123 :787–793. DOI: 10.1016/j.foodchem.2010.05.027 - 84.
Bracci T, Busconi M, Fogher C, Sebastiani L. Molecular studies in olive ( Olea europaea L.): Overview on DNA markers applications and recent advances in genome analysis. Plant Cell Rep. 2011;30 :449–462. DOI: 10.1007/s00299-010-0991-9 - 85.
Busconi M, Reggi S, Dallolio G, Fogher C. Food microbiota diversity. In: Grillo A, Venora G (Eds.). Changing Diversity in Changing Environment. Intech, Croatia 2011; p. 17–32. DOI: 10.5772/24841 - 86.
Busconi M, Zacconi C, Scolari G. Bacterial ecology of PDO Coppa and Pancetta Piacentina at the end of ripening and after MAP storage of sliced product. Int J Food Microbiol. 2014; 172 :13–20. DOI: 10.1016/j.ijfoodmicro.2013.11.023