Development and mapping of SNP assays in allotetraploid cotton

A narrow germplasm base and a complex allotetraploid genome have made the discovery of single nucleotide polymorphism (SNP) markers difficult in cotton (Gossypium hirsutum). To generate sequence for SNP discovery, we conducted a genome reduction experiment (EcoRI, BafI double digest, followed by adapter ligation, biotin–streptavidin purification, and agarose gel separation) on two accessions of G. hirsutum and two accessions of G. barbadense. From the genome reduction experiment, a total of 2.04 million genomic sequence reads were assembled into contigs with an N₅₀ of 508 bp and analyzed for SNPs. A previously generated assembly of expressed sequence tags (ESTs) provided an additional source for SNP discovery. Using highly conservative parameters (minimum coverage of 8× at each SNP and 20% minor allele frequency), a total of 11,834 and 1,679 non-genic SNPs were identified between accessions of G. hirsutum and G. barbadense in genome reduction assemblies, respectively. An additional 4,327 genic SNPs were also identified between accessions of G. hirsutum in the EST assembly. KBioscience KASPar assays were designed for a portion of the intra-specific G. hirsutum SNPs. From 704 non-genic and 348 genic markers developed, a total of 367 (267 non-genic, 100 genic) mapped in a segregating F₂ population (Acala Maxxa × TX2094) using the Fluidigm EP1 system. A G. hirsutum genetic linkage map of 1,688 cM was constructed based entirely on these new SNP markers. Of the genic-based SNPs, we were able to identify within which genome (‘A’ or ‘D’) each SNP resided using diploid species sequence data. Genetic maps generated by these newly identified markers are being used to locate quantitative, economically important regions within the cotton genome.     

Genetic Diversity and Core Collection Evaluations in Common Wheat Germplasm from the Northwestern Spring Wheat Region in China

Fluorescence microsatellite markers were employed to reveal genetic diversity of 340 wheat accessions consisting of 229 landraces and 111 modern varieties from the Northwest Spring Wheat Region in China. The 340 accessions were chosen as candidate core collections for wheat germplasm in this region. A core collection representing the genetic diversity of these accessions was identified based on a cluster dendrogram of 78 SSR loci. A total of 967 alleles were detected with a mean of 13.6 alleles (5–32) per locus. Mean PIC was 0.64, ranged from 0.05 to 0.91. All loci were distributed relatively evenly in the A, B and D wheat genomes. Mean genetic richness of A, B and D genomes for both landraces and modern varieties was B > A > D. However, mean genetic diversity indices of landraces changed to B > D > A. As a whole, genetic diversity of the landraces was considerably higher than that of the modern varieties. The big difference of genetic diversity indices in the three genomes suggested that breeding has exerted greater selection pressure in the D than the A or B genomes in this region. Changes of allelic proportions represented in the proposed core collection at different sampling scales suggested that the sampling percentage of the core collection in the Northwest Spring Wheat Region should be greater than 4% of the base collection to ensure that more than 70% of the variation is represented by the core collection. Electronic supplementary material Electronic supplementary material is available for this article at and accessible for authorised users.     

A survey of genome‐wide single nucleotide polymorphisms through genome resequencing in the Périgord black truffle (Tuber melanosporum Vittad.)

The Périgord black truffle (Tuber melanosporum Vittad.), considered a gastronomic delicacy worldwide, is an ectomycorrhizal filamentous fungus that is ecologically important in Mediterranean French, Italian and Spanish woodlands. In this study, we developed a novel resource of single nucleotide polymorphisms (SNPs) for T. melanosporum using Illumina high‐throughput resequencing. The genome from six T. melanosporum geographical accessions was sequenced to a depth of approximately 20×. These geographical accessions were selected from different populations within the northern and southern regions of the geographical species distribution. Approximately 80% of the reads for each of the six resequenced geographical accessions mapped against the reference T. melanosporum genome assembly, estimating the core genome size of this organism to be approximately 110 Mbp. A total of 442 326 SNPs corresponding to 3540 SNPs/Mbps were identified as being included in all seven genomes. The SNPs occurred more frequently in repeated sequences (85%), although 4501 SNPs were also identified in the coding regions of 2587 genes. Using the ratio of nonsynonymous mutations per nonsynonymous site (pN) to synonymous mutations per synonymous site (pS) and Tajima's D index scanning the whole genome, we were able to identify genomic regions and genes potentially subjected to positive or purifying selection. The SNPs identified represent a valuable resource for future population genetics and genomics studies.     

Quantitative structure analysis of genetic diversity among spring bread wheats (<Emphasis Type="Italic">Triticum aestivum</Emphasis> L.) from different geographical regions

Genetic diversity in spring bread wheat (T.␣aestivum L.) was studied in a total of 69 accessions. For this purpose, 52 microsatellite (SSR) markers were used and a total of 406 alleles were detected, of which 182 (44.8%) occurred at a frequency of <5% (rare alleles). The number of alleles per locus ranged from 2 to 14 with an average of 7.81. The largest number of alleles per locus occurred in the B genome (8.65) as␣compared to the A (8.43) and D (5.93) genomes, respectively. The polymorphism index content (PIC) value varied from 0.24 to 0.89 with an average of 0.68. The highest PIC for all accessions was found in the B␣genome (0.71) as compared to the A (0.68) and D␣genomes (0.63). Genetic distance-based method (standard UPGMA clustering) and a model-based method (structure analysis) were used for cluster analysis. The two methods led to analogical results. Analysis of molecular variance (AMOVA) showed that 80.6% of the total variation could be explained by the variance within the geographical groups. In comparison to the diversity detected for all accessions (H _e = 0.68), genetic diversity among European spring bread wheats was H _e = 0.65. A comparatively higher diversity was observed between wheat varieties from Southern European countries (Austria/Switzerland, Portugal/Spain) corresponding to those from other regions.     

Identification and validation of a core set of informative genic SSR and SNP markers for assaying functional diversity in barley

A ‘core set’ of 28 simple sequence repeat (SSR) and 28 single nucleotide polymorphism (SNP) markers for barley was developed after screening six diverse genotypes (DGs) representing six countries (Afghanistan, Pakistan, Algeria, Egypt, Jordan and Syria) with 50 SSR and 50 SNP markers derived from expressed sequence tags (ESTs). The markers of the core set are single locus with very high quality amplifications, high polymorphism information content (PIC) and are distributed across the barley genome. PIC values for the selected SSR and SNP markers ranged between 0.32–0.72 (average 0.58) and 0.28–0.50 (average 0.42), respectively. To make the SNP genotyping cost effective, CAPS (cleaved amplified polymorphic sequence) and indel assays were developed for 23 markers and the remaining 5 SNP markers were optimized for pyrosequencing. A high coefficient of correlations (r = 0.96, P < 0.005) between the genetic similarity matrices of SSR and SNP genotyping data of the core set on diverse genotypes (DGs) and their similar groupings according to the geographical distribution in both SSR and SNP phenograms with high bootstrap values underline the utility and reliability of the core set. A comparative allelic and sequence diversity for SSR and SNP markers between the DGs and six elite parental genotypes (PGs) of mapping populations showed comparable diverse nature of two germplasm sets. However, unique SNPs and indels were observed in both germplasm sets providing more datapoints for analysing haplotypes in a better way for the corresponding SNP marker. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Development and Evaluation of SoySNP50K,a High-Density Genotyping Array for Soybean

The objective of this research was to identify single nucleotide polymorphisms (SNPs) and to develop an Illumina Infinium BeadChip that contained over 50,000 SNPs from soybean (Glycine max L. Merr.). A total of 498,921,777 reads 35–45bp in length were obtained from DNA sequence analysis of reduced representation libraries from several soybean accessions which included six cultivated and two wild soybean (G. soja Sieb. et Zucc.) genotypes. These reads were mapped to the soybean whole genome sequence and 209,903 SNPs were identified. After applying several filters, a total of 146,161 of the 209,903 SNPs were determined to be ideal candidates for Illumina Infinium II BeadChip design. To equalize the distance between selected SNPs, increase assay success rate, and minimize the number of SNPs with low minor allele frequency, an iteration algorithm based on a selection index was developed and used to select 60,800 SNPs for Infinium BeadChip design. Of the 60,800 SNPs, 50,701 were targeted to euchromatic regions and 10,000 to heterochromatic regions of the 20 soybean chromosomes. In addition, 99 SNPs were targeted to unanchored sequence scaffolds. Of the 60,800 SNPs, a total of 52,041 passed Illumina’s manufacturing phase to produce the SoySNP50K iSelect BeadChip. Validation of the SoySNP50K chip with 96 landrace genotypes, 96 elite cultivars and 96 wild soybean accessions showed that 47,337 SNPs were polymorphic and generated successful SNP allele calls. In addition, 40,841 of the 47,337 SNPs (86%) had minor allele frequencies ≥10% among the landraces, elite cultivars and the wild soybean accessions. A total of 620 and 42 candidate regions which may be associated with domestication and recent selection were identified, respectively. The SoySNP50K iSelect SNP beadchip will be a powerful tool for characterizing soybean genetic diversity and linkage disequilibrium, and for constructing high resolution linkage maps to improve the soybean whole genome sequence assembly.     

Single nucleotide polymorphism genotyping in polyploid wheat with the Illumina GoldenGate assay

Single nucleotide polymorphisms (SNPs) are indispensable in such applications as association mapping and construction of high-density genetic maps. These applications usually require genotyping of thousands of SNPs in a large number of individuals. Although a number of SNP genotyping assays are available, most of them are designed for SNP genotyping in diploid individuals. Here, we demonstrate that the Illumina GoldenGate assay could be used for SNP genotyping of homozygous tetraploid and hexaploid wheat lines. Genotyping reactions could be carried out directly on genomic DNA without the necessity of preliminary PCR amplification. A total of 53 tetraploid and 38 hexaploid homozygous wheat lines were genotyped at 96 SNP loci. The genotyping error rate estimated after removal of low-quality data was 0 and 1% for tetraploid and hexaploid wheat, respectively. Developed SNP genotyping assays were shown to be useful for genotyping wheat cultivars. This study demonstrated that the GoldenGate assay is a very efficient tool for high-throughput genotyping of polyploid wheat, opening new possibilities for the analysis of genetic variation in wheat and dissection of genetic basis of complex traits using association mapping approach. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Identification and characterization of more than 4 million intervarietal SNPs across the group 7 chromosomes of bread wheat

Despite being a major international crop, our understanding of the wheat genome is relatively poor due to its large size and complexity. To gain a greater understanding of wheat genome diversity, we have identified single nucleotide polymorphisms between 16 Australian bread wheat varieties. Whole‐genome shotgun Illumina paired read sequence data were mapped to the draft assemblies of chromosomes 7A, 7B and 7D to identify more than 4 million intervarietal SNPs. SNP density varied between the three genomes, with much greater density observed on the A and B genomes than the D genome. This variation may be a result of substantial gene flow from the tetraploid Triticum turgidum, which possesses A and B genomes, during early co‐cultivation of tetraploid and hexaploid wheat. In addition, we examined SNP density variation along the chromosome syntenic builds and identified genes in low‐density regions which may have been selected during domestication and breeding. This study highlights the impact of evolution and breeding on the bread wheat genome and provides a substantial resource for trait association and crop improvement. All SNP data are publically available on a generic genome browser GBrowse at www.wheatgenome.info .     

Microsatellite-based molecular diversity of bread wheat germplasm and association mapping of wheat resistance to the Russian wheat aphid

Genetic diversity of a set of 71 wheat accessions, including 53 biotype 2 Russian wheat aphid (RWA2)-resistant landraces and 18 RWA2 susceptible accessions, was assessed by examining molecular variation at multiple microsatellite (SSR) loci. Fifty-one wheat SSR primer pairs were used, 81 SSR loci were determined, and 545 SSR alleles were detected. These SSR loci covered all the three genomes, 21 chromosomes, and at least 41 of the 42 chromosome arms. Diversity values averaged over SSR loci were high with mean number of SSR alleles/locus = 6.7, mean Shannon’s index (H) = 1.291, and mean Nei’s gene diversity (He) = 0.609. The three wheat genomes ranked as A > D > B and the homoeologous groups ranked as 7 > 3  > 1 > 2 > 6 > 5 > 4 based on the number of alleles per locus. Xgwm136 on chromosome arm 1AS is the most polymorphic SSR locus with the largest number of observed and effective alleles and the highest H and He. Among all 2485 pairs of wheat accessions, genetic distance (GD) ranged from 0.054 to 1.933 and averaged 0.9832. A dendrogram based on GD matrix showed that all the wheat accessions could be grouped into distinct clusters. Most of the susceptible cultivars (13/18) were clustered into groups that contains all or mostly susceptible accessions. Most of the U.S. cultivars belong to a group that is distinguishable from all the different RWA2 resistant groups. Diversity analysis was also conducted separately for subgroups containing 53 RWA2-resistant accessions and 18 RWA2-susceptible accessions. Association mapping revealed 28 SSR loci significantly associated with leaf chlorosis, and 8 with leaf rolling. New chromosome regions associated with RWA2 resistance were detected, and indicated existence of new RWA resistance genes located on chromosomes of all other homoeologous groups in addition to the groups 1 and 7 in bread wheat. This information is helpful for development of mapping populations for RWA2 resistance genes from different phylogenetic groups, and for wise utilization of the RWA-resistant germplasm in wheat breeding programs.     

Development of SNP-based CAPS and dCAPS markers in eight different genes involved in starch biosynthesis in rice

Single nucleotide polymorphisms (SNPs), which are inexhaustible, highly stable, and simply detectable sequence polymorphisms, can lead to phenotypic variations by affecting protein composition changes. Here, we report development of 25 new cleaved amplified polymorphic sequence or derived cleaved amplified polymorphic sequence markers that have discrete band sizes in relation to the SNP genotypes in eight putative gene regions. The average frequency of DNA polymorphisms was 1 per 175 bp (SNPs, 1 per 217 bp; In/dels, 1 per 906 bp). In primary statistical analysis of each marker on 55 diverse rice accessions, including different ecotypes, the mean value of the major allele frequency was 0.658 (0.509–0.927). The average polymorphism information content was 0.326 (0.126–0.375). The mean value of the inbreeding coefficient (f) was 0.950 and was positive (heterozygote deficiency) at all loci, corresponding to the inbreeding system in rice. In cluster analysis, all rice accessions clustered mainly into three groups according to the ecotypes. The association analysis showed that the SNP of Granule-bound starch synthase I and ADP-glucose pyrophosphorylase small subunit (ADPase-S) genes were highly associated with apparent amylose content variation than the others. These new SNP markers may be useful in genotyping rice germplasm, in marker-assisted selection for improving starch quality and content, and in linkage as well as association studies. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Identification of new sources of aluminum resistance in wheat

Aluminum (Al) toxicity is a major constraint for wheat production in acidic soils. An Al resistance gene on chromosome 4DL that traces to Brazilian wheat has been extensively studied, and can provide partial protection from Al damage. To identify potentially new sources of Al resistance, 590 wheat accessions, including elite wheat breeding lines from the United States and other American and European countries, landraces and commercial cultivars from East Asia, and synthetic wheat lines from CIMMYT, Mexico, were screened for Al resistance by measuring relative root elongation in culture with a nutrient solution containing Al, and by staining Al-stressed root tips with hematoxylin. Eighty-eight wheat accessions demonstrated at least moderate resistance to Al toxicity. Those selected lines were subjected to analysis of microsatellite markers linked to an Al resistance gene on 4DL and a gene marker for the Al-activated malate transporter (ALMT1) locus. Many of the selected Al-resistant accessions from East Asia did not have the Al-resistant marker alleles of ALMT1, although they showed Al resistance similar to the US Al-resistant cultivar, Atlas 66. Most of the cultivars derived from Jagger and Atlas 66 have the Al-resistant marker alleles of ALMT1. Cluster analysis separated the selected Al-resistant germplasm into two major clusters, labeled as Asian and American–European clusters. Potentially new germplasm of Al resistance different from those derived from Brazil were identified. Further investigation of Al resistance in those new germplasms may reveal alternative Al-resistance mechanisms in wheat. Electronic supplementary material The online version of this article (doi:contains supplementary material, which is available to authorized users. Responsible Editor: Thomas B. Kinraide.     

Genome-specific primer sets for starch biosynthesis genes in wheat

Common wheat (Triticum aestivum L., 2n=6x=42) is an allohexaploid composed of three closely related genomes, designated A, B, and D. Genetic analysis in wheat is complicated, as most genes are present in triplicated sets located in the same chromosomal regions of homoeologous chromosomes. The goal of this report was to use genomic information gathered from wheat–rice sequence comparison to develop genome-specific primer sets for five genes involved in starch biosynthesis. Intron locations in wheat were inferred through the alignment of wheat cDNA sequences with rice genomic sequence. Exon-anchored primers, which amplify across introns, allowed the sequencing of introns from the three genomes for each gene. Sequence variation within introns among the three wheat genomes provided the basis for genome-specific primer design. For three genes, ADP-glucose pyrophosphorylase (Agp-L), sucrose transporter (SUT), and waxy (Wx), genome-specific primer sets were developed for all three genomes. Genome-specific primers were developed for two of the three genomes for Agp-S and starch synthase I (SsI). Thus, 13 of 15 possible genome-specific primer sets were developed using this strategy. Seven genome-specific primer combinations were used to amplify alleles in hexaploid wheat lines for sequence comparison. Three single nucleotide polymorphisms (SNPs) were identified in a comparison of 5,093 bp among a minimum of ten wheat accessions. Two of these SNPs could be converted into cleaved amplified polymorphism sequence (CAPS) markers. Our results indicated that the design of genome-specific primer sets using intron-based sequence differences has a high probability of success, while the identification of polymorphism among alleles within a genome may be a challenge.     

Genotype‐specific SNP map based on whole chromosome 3B sequence information from wheat cultivars Arina and Forno

Agronomically important traits are frequently controlled by rare, genotype‐specific alleles. Such genes can only be mapped in a population derived from the donor genotype. This requires the development of a specific genetic map, which is difficult in wheat because of the low level of polymorphism among elite cultivars. The absence of sufficient polymorphism, the complexity of the hexaploid wheat genome as well as the lack of complete sequence information make the construction of genetic maps with a high density of reproducible and polymorphic markers challenging. We developed a genotype‐specific genetic map of chromosome 3B from winter wheat cultivars Arina and Forno. Chromosome 3B was isolated from the two cultivars and then sequenced to 10‐fold coverage. This resulted in a single‐nucleotide polymorphisms (SNP) database of the complete chromosome. Based on proposed synteny with the Brachypodium model genome and gene annotation, sequences close to coding regions were used for the development of 70 SNP‐based markers. They were mapped on a Arina × Forno Recombinant Inbred Lines population and found to be spread over the complete chromosome 3B. While overall synteny was well maintained, numerous exceptions and inversions of syntenic gene order were identified. Additionally, we found that the majority of recombination events occurred in distal parts of chromosome 3B, particularly in hot‐spot regions. Compared with the earlier map based on SSR and RFLP markers, the number of markers increased fourfold. The approach presented here allows fast development of genotype‐specific polymorphic markers that can be used for mapping and marker‐assisted selection.     

Characterization of genetic diversity and population structure in wheat using array based SNP markers

Genetic diversity is crucial for successful adaptation and sustained improvement in crops. India is bestowed with diverse agro-climatic conditions which makes it rich in wheat germplasm adapted to various niches. Germplasm repository consists of local landraces, trait specific genetic stocks including introgressions from wild relatives, exotic collections, released varieties, and improved germplasm. Characterization of genetic diversity is done using morpho-physiological characters as well as by analyzing variations at DNA level. However, there are not many reports on array based high throughput SNP markers having characteristics of genome wide coverage employed in Indian spring wheat germplasm. Amongst wheat SNP arrays, 35K Axiom Wheat Breeder’s Array has the highest SNP polymorphism efficiency suitable for genetic mapping and genetic diversity characterization. Therefore, genotyping was done using 35K in 483 wheat genotypes resulting in 14,650 quality filtered SNPs, that were distributed across the B (~?50%), A (~?39%), and D (~?10%) genomes. The total genetic distance coverage was 4477.85 cM with 3.27 SNP/cM and 0.49 cM/SNP as average marker density and average inter-marker distance, respectively. The PIC ranged from 0.09 to 0.38 with an average of 0.29 across genomes. Population structure and Principal Coordinate Analysis resulted in two subpopulations (SP1 and SP2). The analysis of molecular variance revealed the genetic variation of 2% among and 98% within subpopulations indicating high gene flow between SP1 and SP2. The subpopulation SP2 showed high level of genetic diversity based on genetic diversity indices viz. Shannon’s information index (I)?=?0.648, expected heterozygosity (He)?=?0.456 and unbiased expected heterozygosity (uHe)?=?0.456. To the best of our knowledge, this study is the first to include the largest set of Indian wheat genotypes studied exclusively for genetic diversity. These findings may serve as a potential source for the identification of uncharacterized QTL/gene using genome wide association studies and marker assisted selection in wheat breeding programs.

     

Genome-wide searching of single-nucleotide polymorphisms among eight distantly and closely related rice cultivars (Oryza sativa L.) and a wild accession (Oryza rufipogon Griff.).

We searched the genomes of eight rice cultivars (Oryza sativa L. ssp. japonica and ssp. indica) and a wild rice accession (Oryza rufipogon Griffith) for nucleotide polymorphisms, and identified 7805 polymorphic loci, including single-nucleotide polymorphisms (SNPs) and insertions/deletions (InDels), in predicted intergenic regions. Polymorphisms are useful as DNA markers for genetic analysis or positional cloning with segregating populations of crosses. Pairwise comparison between cultivars and a neighbor-joining tree calculated from SNPs agreed very well with relationships between rice strains predicted from pedigree data or calculated with other DNA markers such as p-SINE1 and simple sequence repeats (SSRs), suggesting that whole-genome SNP information can be used for analysis of evolutionary relationships. Using multiple SNPs to identify alleles, we drew a map to illustrate the alleles shared among the eight cultivars and the accession. The map revealed that most of the genome is mono- or di-allelic among japonica cultivars, whereas alleles well conserved among modern japonica paddy rice cultivars were often shared with indica cultivars or wild rice, suggesting that the genome structure of modern cultivars is composed of chromosomal segments from various genetic backgrounds. Use of allele-sharing analysis and association analysis were also tested and are discussed.     

Single nucleotide polymorphisms and linkage disequilibrium in sunflower

Genetic diversity in modern sunflower (Helianthus annuus L.) cultivars (elite oilseed inbred lines) has been shaped by domestication and breeding bottlenecks and wild and exotic allele introgression⁻the former narrowing and the latter broadening genetic diversity. To assess single nucleotide polymorphism (SNP) frequencies, nucleotide diversity, and linkage disequilibrium (LD) in modern cultivars, alleles were resequenced from 81 genic loci distributed throughout the sunflower genome. DNA polymorphisms were abundant; 1078 SNPs (1/45.7 bp) and 178 insertions-deletions (INDELs) (1/277.0 bp) were identified in 49.4 kbp of DNA/genotype. SNPs were twofold more frequent in noncoding (1/32.1 bp) than coding (1/62.8 bp) sequences. Nucleotide diversity was only slightly lower in inbred lines (θ = 0.0094) than wild populations (θ = 0.0128). Mean haplotype diversity was 0.74. When extraploted across the genome (~3500 Mbp), sunflower was predicted to harbor at least 76.4 million common SNPs among modern cultivar alleles. LD decayed more slowly in inbred lines than wild populations (mean LD declined to 0.32 by 5.5 kbp in the former, the maximum physical distance surveyed), a difference attributed to domestication and breeding bottlenecks. SNP frequencies and LD decay are sufficient in modern sunflower cultivars for very high-density genetic mapping and high-resolution association mapping.     

Development of a 135K SNP genotyping array for Actinidia arguta and its applications for genetic mapping and QTL analysis in kiwifruit

Kiwifruit (Actinidia spp) is a woody, perennial and deciduous vine. In this genus, there are multiple ploidy levels but the main cultivated cultivars are polyploid. Despite the availability of many genomic resources in kiwifruit, SNP genotyping is still a challenge given these different levels of polyploidy. Recent advances in SNP array technologies have offered a high-throughput genotyping platform for genome-wide DNA polymorphisms. In this study, we developed a high-density SNP genotyping array to facilitate genetic studies and breeding applications in kiwifruit. SNP discovery was performed by genome-wide DNA sequencing of 40 kiwifruit genotypes. The identified SNPs were stringently filtered for sequence quality, predicted conversion performance and distribution over the available Actinidia chinensis genome. A total of 134 729 unique SNPs were put on the array. The array was evaluated by genotyping 400 kiwifruit individuals. We performed a multidimensional scaling analysis to assess the diversity of kiwifruit germplasm, showing that the array was effective to distinguish kiwifruit accessions. Using a tetraploid F1 population, we constructed an integrated linkage map covering 3060.9 cM across 29 linkage groups and performed QTL analysis for the sex locus that has been identified on Linkage Group 3 (LG3) in Actinidia arguta. Finally, our dataset presented evidence of tetrasomic inheritance with partial preferential pairing in A. arguta. In conclusion, we developed and evaluated a 135K SNP genotyping array for kiwifruit. It has the advantage of a comprehensive design that can be an effective tool in genetic studies and breeding applications in this high-value crop.     

Whole Mitochondrial and Plastid Genome SNP Analysis of Nine Date Palm Cultivars Reveals Plastid Heteroplasmy and Close Phylogenetic Relationships among Cultivars

Date palm is a very important crop in western Asia and northern Africa, and it is the oldest domesticated fruit tree with archaeological records dating back 5000 years. The huge economic value of this crop has generated considerable interest in breeding programs to enhance production of dates. One of the major limitations of these efforts is the uncertainty regarding the number of date palm cultivars, which are currently based on fruit shape, size, color, and taste. Whole mitochondrial and plastid genome sequences were utilized to examine single nucleotide polymorphisms (SNPs) of date palms to evaluate the efficacy of this approach for molecular characterization of cultivars. Mitochondrial and plastid genomes of nine Saudi Arabian cultivars were sequenced. For each species about 60 million 100 bp paired-end reads were generated from total genomic DNA using the Illumina HiSeq 2000 platform. For each cultivar, sequences were aligned separately to the published date palm plastid and mitochondrial reference genomes, and SNPs were identified. The results identified cultivar-specific SNPs for eight of the nine cultivars. Two previous SNP analyses of mitochondrial and plastid genomes identified substantial intra-cultivar ( = intra-varietal) polymorphisms in organellar genomes but these studies did not properly take into account the fact that nearly half of the plastid genome has been integrated into the mitochondrial genome. Filtering all sequencing reads that mapped to both organellar genomes nearly eliminated mitochondrial heteroplasmy but all plastid SNPs remained heteroplasmic. This investigation provides valuable insights into how to deal with interorganellar DNA transfer in performing SNP analyses from total genomic DNA. The results confirm recent suggestions that plastid heteroplasmy is much more common than previously thought. Finally, low levels of sequence variation in plastid and mitochondrial genomes argue for using nuclear SNPs for molecular characterization of date palm cultivars.