Future tools for association mapping in crop plants

Association mapping currently relies on the identification of genetic markers. Several technologies have been adopted for genetic marker analysis, with single nucleotide polymorphisms (SNPs) being the most popular where a reasonable quantity of genome sequence data are available. We describe several tools we have developed for the discovery, annotation, and visualization of molecular markers for association mapping. These include autoSNPdb for SNP discovery from assembled sequence data; TAGdb for the identification of gene specific paired read Illumina GAII data; CMap3D for the comparison of mapped genetic and physical markers; and BAC and Gene Annotator for the online annotation of genes and genomic sequences.     

Construction of a high-density SNP genetic map in flue-cured tobacco based on SLAF-seq

Tobacco (Nicotiana tabacum L., 2n = 48) is an important agronomic crop and model plant. Flue-cured tobacco is the most important type and accounts for approximately 80 % of tobacco production worldwide. The low genetic diversity of flue-cured tobacco impedes the construction of a high-density genetic linkage map using simple sequence repeat (SSR) markers and warrants the exploitation of single nucleotide polymorphic (SNP) markers from genomic regions. In this article, initially using specific locus-amplified fragment sequencing, we discovered 10,891 SNPs that were subsequently used as molecular markers for genetic map construction. Combined with SSR markers, a final high-density genetic map was generated containing 4215 SNPs and 194 SSRs distributed on 24 linkage groups (LGs). The genetic map was 2662.43 cM in length, with an average distance of 0.60 cM between adjacent markers. Furthermore, by mapping the SNP markers to the ancestral genomes of Nicotiana tomentosiformis and Nicotiana sylvestris, a large number of genome rearrangements were identified as occurring after the polyploidization event. Finally, using this novel integrated map and mapping population, two major quantitative trait loci (QTLs) were identified for flue-curing and mapped to the LG6 of tobacco. This is the first report of SNP markers and a SNP-based linkage map being developed in tobacco. The high-density genetic map and QTLs related to tobacco curing will support gene/QTL fine mapping, genome sequence assembly and molecular breeding in tobacco.     

High-throughput SNP discovery and genotyping in durum wheat (<Emphasis Type="Italic">Triticum durum</Emphasis> Desf.)

We describe the application of complexity reduction of polymorphic sequences (CRoPS^®) technology for the discovery of SNP markers in tetraploid durum wheat (Triticum durum Desf.). A next-generation sequencing experiment was carried out on reduced representation libraries obtained from four durum cultivars. SNP validation and minor allele frequency (MAF) estimate were carried out on a panel of 12 cultivars, and the feasibility of genotyping these SNPs in segregating populations was tested using the Illumina Golden Gate (GG) technology. A total of 2,659 SNPs were identified on 1,206 consensus sequences. Among the 768 SNPs that were chosen irrespective of their genomic repetitiveness level and assayed on the Illumina BeadExpress genotyping system, 275 (35.8%) SNPs matched the expected genotypes observed in the SNP discovery phase. MAF data indicated that the overall SNP informativeness was high: a total of 196 (71.3%) SNPs had MAF >0.2, of which 76 (27.6%) showed MAF >0.4. Of these SNPs, 157 were mapped in one of two mapping populations (Meridiano × Claudio and Colosseo × Lloyd) and integrated into a common genetic map. Despite the relatively low genotyping efficiency of the GG assay, the validated CRoPS-derived SNPs showed valuable features for genomics and breeding applications such as a uniform distribution across the wheat genome, a prevailing single-locus codominant nature and a high polymorphism. Here, we report a new set of 275 highly robust genome-wide Triticum SNPs that are readily available for breeding purposes.     

A high density genetic map of tobacco (<Emphasis Type="Italic">Nicotiana tabacum</Emphasis> L.) obtained from large scale microsatellite marker development

Tobacco (Nicotiana tabacum L.) is a species in the large family of the Solanaceae and is important as an agronomic crop and as a model system in plant biotechnology. Despite its importance, only limited molecular marker resources are available that can be used for genome analysis, genetic mapping and breeding. We report here on the development and characterization of 5,119 new and functional microsatellite markers and on the generation of a high-resolution genetic map for the tetraploid tobacco genome. The genetic map was generated using an F2 mapping population derived from the intervarietal cross of Hicks Broadleaf × Red Russian and merges the polymorphic markers from this new set with those from a smaller set previously used to produce a lower density map. The genetic map described here contains 2,317 microsatellite markers and 2,363 loci, resulting in an average distance between mapped microsatellite markers which is less than 2 million base pairs or 1.5 cM. With this new and expanded marker resource, a sufficient number of markers are now available for multiple applications ranging from tobacco breeding to comparative genome analysis. The genetic map of tobacco is now comparable in marker density and resolution with the best characterized genomes of the Solanaceae: tomato and potato.     

Genome-Wide Discovery of DNA Polymorphisms in Mei (<Emphasis Type="Italic">Prunus mume</Emphasis> Sieb. et Zucc.), an Ornamental Woody Plant,with Contrasting Tree Architecture and their Functional Relevance for Weeping Trait

Next-generation sequencing technologies provide opportunities to ascertain the genetic basis of phenotypic differences, even in the closely related cultivars via detection of large amount of DNA polymorphisms. In this study, we performed whole-genome re-sequencing of two mei cultivars with contrasting tree architecture. 75.87 million 100 bp pair-end reads were generated, with 92 % coverage of the genome. Re-sequencing data of two former upright mei cultivars were applied for detecting DNA polymorphisms, since we were more interested in variations conferring weeping trait. Applying stringent parameters, 157,317 mutual single nucleotide polymorphisms (SNPs) and 15,064 mutual insertions-deletions (InDels) were detected and found unevenly distributed within and among the mei chromosomes, which lead to the discovery of 220 high-density, 463 low-density SNP regions together with 80 high-density InDel regions. Additionally, 322 large-effect SNPs and 433 large-effect InDels were detected, and 10.09 % of the SNPs were observed in coding regions. 5.25 % SNPs in coding regions resulted in non-synonymous changes. Ninety SNPs were chosen randomly for validation using high-resolution melt analysis. 93.3 % of the candidate SNPs contained the predicted SNPs. Pfam analysis was further conducted to better understand SNP effects on gene functions. DNA polymorphisms of two known QTL loci conferring weeping trait and their functional effect were also analyzed thoroughly. This study highlights promising functional markers for molecular breeding and a whole-genome genetic basis of weeping trait in mei.     

Pedigree‐based analysis of derivation of genome segments of an elite rice reveals key regions during its breeding

Analyses of genome variations with high‐throughput assays have improved our understanding of genetic basis of crop domestication and identified the selected genome regions, but little is known about that of modern breeding, which has limited the usefulness of massive elite cultivars in further breeding. Here we deploy pedigree‐based analysis of an elite rice, Huanghuazhan, to exploit key genome regions during its breeding. The cultivars in the pedigree were resequenced with 7.6× depth on average, and 2.1 million high‐quality single nucleotide polymorphisms (SNPs) were obtained. Tracing the derivation of genome blocks with pedigree and information on SNPs revealed the chromosomal recombination during breeding, which showed that 26.22% of Huanghuazhan genome are strictly conserved key regions. These major effect regions were further supported by a QTL mapping of 260 recombinant inbred lines derived from the cross of Huanghuazhan and a very dissimilar cultivar, Shuanggui 36, and by the genome profile of eight cultivars and 36 elite lines derived from Huanghuazhan. Hitting these regions with the cloned genes revealed they include numbers of key genes, which were then applied to demonstrate how Huanghuazhan were bred after 30 years of effort and to dissect the deficiency of artificial selection. We concluded the regions are helpful to the further breeding based on this pedigree and performing breeding by design. Our study provides genetic dissection of modern rice breeding and sheds new light on how to perform genomewide breeding by design.     

Development and evaluation of a 9K SNP array for peach by internationally coordinated SNP detection and validation in breeding germplasm

Although a large number of single nucleotide polymorphism (SNP) markers covering the entire genome are needed to enable molecular breeding efforts such as genome wide association studies, fine mapping, genomic selection and marker-assisted selection in peach [Prunus persica (L.) Batsch] and related Prunus species, only a limited number of genetic markers, including simple sequence repeats (SSRs), have been available to date. To address this need, an international consortium (The International Peach SNP Consortium; IPSC) has pursued a coordinated effort to perform genome-scale SNP discovery in peach using next generation sequencing platforms to develop and characterize a high-throughput Illumina Infinium® SNP genotyping array platform. We performed whole genome re-sequencing of 56 peach breeding accessions using the Illumina and Roche/454 sequencing technologies. Polymorphism detection algorithms identified a total of 1,022,354 SNPs. Validation with the Illumina GoldenGate® assay was performed on a subset of the predicted SNPs, verifying ∼75% of genic (exonic and intronic) SNPs, whereas only about a third of intergenic SNPs were verified. Conservative filtering was applied to arrive at a set of 8,144 SNPs that were included on the IPSC peach SNP array v1, distributed over all eight peach chromosomes with an average spacing of 26.7 kb between SNPs. Use of this platform to screen a total of 709 accessions of peach in two separate evaluation panels identified a total of 6,869 (84.3%) polymorphic SNPs.The almost 7,000 SNPs verified as polymorphic through extensive empirical evaluation represent an excellent source of markers for future studies in genetic relatedness, genetic mapping, and dissecting the genetic architecture of complex agricultural traits. The IPSC peach SNP array v1 is commercially available and we expect that it will be used worldwide for genetic studies in peach and related stone fruit and nut species.     

Genome-wide association mapping of agronomic traits in sugar beet

Recent results indicate that association mapping in populations from applied plant breeding is a powerful tool to detect QTL which are of direct relevance for breeding. The focus of this study was to unravel the genetic architecture of six agronomic traits in sugar beet. To this end, we employed an association mapping approach, based on a very large population of 924 elite sugar beet lines from applied plant breeding, fingerprinted with 677 single nucleotide polymorphism (SNP) markers covering the entire genome. We show that in this population linkage disequilibrium decays within a short genetic distance and is sufficient for the detection of QTL with a large effect size. To increase the QTL detection power and the mapping resolution a much higher number of SNPs is required. We found that for QTL detection, the mixed model including only the kinship matrix performed best, even in the presence of a considerable population structure. In genome-wide scans, main effect QTL and epistatic QTL were detected for all six traits. Our full two-dimensional epistasis scan revealed that for complex traits there appear to be epistatic master regulators, loci which are involved in a large number of epistatic interactions throughout the genome.     

Whole genome SNP discovery and analysis of genetic diversity in Turkey (Meleagris gallopavo)

ABSTRACT: BACKGROUND: The turkey (Meleagris gallopavo) is an important agricultural species and the second largest contributor to the world's poultry meat production. Genetic improvement is attributed largely to selective breeding programs that rely on highly heritable phenotypic traits, such as body size and breast muscle development. Commercial breeding with small effective population sizes and epistasis can result in loss of genetic diversity, which in turn can lead to reduced individual fitness and reduced response to selection. The presence of genomic diversity in domestic livestock species therefore, is of great importance and a prerequisite for rapid and accurate genetic improvement of selected breeds in various environments, as well as to facilitate rapid adaptation to potential changes in breeding goals. Genomic selection requires a large number of genetic markers such as e.g. single nucleotide polymorphisms (SNPs) the most abundant source of genetic variation within the genome. RESULTS: Alignment of next generation sequencing data of 32 individual turkeys from different populations was used for the discovery of 5.49 million SNPs, which subsequently were used for the analysis of genetic diversity among the different populations. All of the commercial lines branched from a single node relative to the heritage varieties and the South Mexican turkey population. Heterozygosity of all individuals from the different turkey populations ranged from 0.17-2.73 SNPs/Kb, while heterozygosity of populations ranged from 0.73-1.64 SNPs/Kb. The average frequency of heterozygous SNPs in individual turkeys was 1.07 SNPs/Kb. Five genomic regions with very low nucleotide variation were identified in domestic turkeys that showed state of fixation towards alleles different than wild alleles. CONCLUSION: The turkey genome is much less diverse with a relatively low frequency of heterozygous SNPs as compared to other livestock species like chicken and pig. The whole genome SNP discovery study in turkey resulted in the detection of 5.49 million putative SNPs compared to the reference genome. All commercial lines appear to share a common origin. Presence of different alleles/haplotypes in the SM population highlights that specific haplotypes have been selected in the modern domesticated turkey.     

Development of high-density genetic maps for barley and wheat using a novel two-enzyme genotyping-by-sequencing approach

Advancements in next-generation sequencing technology have enabled whole genome re-sequencing in many species providing unprecedented discovery and characterization of molecular polymorphisms. There are limitations, however, to next-generation sequencing approaches for species with large complex genomes such as barley and wheat. Genotyping-by-sequencing (GBS) has been developed as a tool for association studies and genomics-assisted breeding in a range of species including those with complex genomes. GBS uses restriction enzymes for targeted complexity reduction followed by multiplex sequencing to produce high-quality polymorphism data at a relatively low per sample cost. Here we present a GBS approach for species that currently lack a reference genome sequence. We developed a novel two-enzyme GBS protocol and genotyped bi-parental barley and wheat populations to develop a genetically anchored reference map of identified SNPs and tags. We were able to map over 34,000 SNPs and 240,000 tags onto the Oregon Wolfe Barley reference map, and 20,000 SNPs and 367,000 tags on the Synthetic W9784 × Opata85 (SynOpDH) wheat reference map. To further evaluate GBS in wheat, we also constructed a de novo genetic map using only SNP markers from the GBS data. The GBS approach presented here provides a powerful method of developing high-density markers in species without a sequenced genome while providing valuable tools for anchoring and ordering physical maps and whole-genome shotgun sequence. Development of the sequenced reference genome(s) will in turn increase the utility of GBS data enabling physical mapping of genes and haplotype imputation of missing data. Finally, as a result of low per-sample costs, GBS will have broad application in genomics-assisted plant breeding programs.     

Complete genome sequencing and analysis of <Emphasis Type="Italic">Capsicum annuum</Emphasis> varieties

Peppers (Capsicum spp.), economically important crops of the nightshade family, are used as vegetables, spices, and medicines. Bacterial wilt (BW), one of the most devastating bacterial plant diseases, causes heavy losses in many crops, including peppers. The aim of this study was to develop DNA-based molecular markers that could be used in future marker-assisted breeding for BW resistance. Genome-wide discovery of single-nucleotide polymorphisms (SNPs was achieved by comparing whole-genome resequencing data from two pepper varieties, Saengryeg 211 (susceptible to BW) and 82PR66 (resistant to BW) with the reference genome of Capsicum annuum cv. CM334. A total of 118,588,231 and 148,774,861 paired raw reads with an average length of 151 bp were recorded for Saengryeg 211 and 82PR66, respectively. Data mining revealed 6,804,889 and 4,293,534 SNPs for Saengryeg 211 and 82PR66, respectively. The SNPs, were assorted into 12 chromosomes. A total of 5,514,563 polymorphic SNPs were discovered between the two varieties, and 42,236 high-resolution melting marker primer sets were selected for use in future experiments. Our results reveal numerous DNA-based molecular markers that might be useful in molecular breeding and QTL mapping in pepper.     

Development of single nucleotide polymorphism markers in the large and complex rubber tree genome using next-generation sequence data

The development of single nucleotide polymorphism (SNP) markers provides the opportunity to improve many areas of plant breeding and population genetics. Unfortunately, for species such as the rubber tree (Hevea brasiliensis), the use of next-generation sequencing for genomic SNP discovery is very difficult because of the large genome size and the abundance of repeated sequences. Access to a set of validated SNP markers is a significant advantage for rubber researchers who wish to apply SNPs in scientific research. Here, we performed genomic sequencing of H. brasiliensis and generated 10,993,648 short reads, which were assembled into 10,071 contigs (N50 = 3078) by a de novo assembly strategy. A total of 2446 contigs presented no hits in the current H. brasiliensis genome assembly and may therefore be considered novel genomic sequences of rubber tree. A total of 143 putative polymorphic positions were selected, gene annotations were available for 58.7 % of the markers, and all of the sequences could be anchored to the released H. brasiliensis genome. These SNPs were validated in eight genotypes of H. brasiliensis and 15 F1 plants from a mapping population, resulting in 30 (20.9 %) positions correctly classified. The analysis revealed key candidate genes responsible for defence mechanisms and provided markers for further genetic improvement of Hevea in breeding programmes.     

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

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

Construction of a high-density genetic map of <Emphasis Type="Italic">Ziziphus jujuba</Emphasis> Mill. using genotyping by sequencing technology

The Chinese jujube (Ziziphus jujuba Mill., 2n = 2 × = 24), one of the most popular fruit trees in China, is widely cultivated and utilized in Asia. High-density genetic linkage maps are valuable resources for molecular breeding and functional genomics; however, they are still under-developed for the jujube. The genotyping by sequencing (GBS) strategy could be an efficient and cost-effective tool for single nucleotide polymorphism (SNP) discovery based on the sequenced jujube genome. Here, we report a new high-density genetic map constructed using GBS technology. An F1 population with 145 progenies and their parents (‘Dongzao’ × ‘Zhongningyuanzao’) were sequenced on the Illumina HiSeq 4000 platform. In total, 79.8 Gb of raw data containing 256,708,177 paired-end reads were generated. After data filtering and SNP genotyping, 40,372 polymorphic SNP markers were developed between the parents and 2540 (1756 non-redundant) markers were mapped onto the integrated genetic linkage map. The map spanned 1456.53 cM and was distributed among 12 linkage groups, which is consistent with the haploid chromosome number of the jujube. The average marker interval was 0.88 cM. The genetic map allowed us to anchor 224 Mb (63.7 %) of scaffolds from the sequenced ‘Junzao’ genome, containing 52 newly anchored scaffolds, which extended the genome assembly by 7 Mb. In conclusion, GBS technology was applied efficiently for SNP discovery in this study. The high-density genetic map will serve as a unique tool for molecular-assisted breeding and genomic studies, which will contribute to further research and improvement of the jujube in the near future.     

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

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

SNPMeta: SNP annotation and SNP metadata collection without a reference genome

The increase in availability of resequencing data is greatly accelerating SNP discovery and has facilitated the development of SNP genotyping assays. This, in turn, is increasing interest in annotation of individual SNPs. Currently, these data are only available through curation, or comparison to a reference genome. Many species lack a reference genome, but are still important genetic models or are significant species in agricultural production or natural ecosystems. For these species, it is possible to annotate SNPs through comparison with cDNA, or data from well‐annotated genes in public repositories. We present SNPMeta, a tool which gathers information about SNPs by comparison with sequences present in GenBank databases. SNPMeta is able to annotate SNPs from contextual sequence in SNP assay designs, and SNPs discovered through genotyping by sequencing (GBS) approaches. However, SNPs discovered through GBS occur throughout the genome, rather than only in gene space, and therefore do not annotate at high rates. SNPMeta can therefore be used to annotate SNPs in nonmodel species or species that lack a reference genome. Annotations generated by SNPMeta are highly concordant with annotations that would be obtained from a reference genome.     

新一代分子标记--SNPs及其应用

单核苷酸多态性(SNPs)是广泛存在于基因组中的一类DNA序列变异,其频率为1%或更高。它是由单个碱基的转换或颠换引起的点突变,稳定而可靠,并通常以二等位基因的形式出现。采用生物芯片和DNA微阵列技术来检测SNP,便于对基因组进行大幅度和高通量分析。因此,作为新一代分子标记,SNP在生物学诸多领域具有广阔应用前景。本文简要叙述SNPs技术的发展历史、研究动态以及相关的理论,介绍了与SNPs相关的基本术语、概念及其特点,列举了发现与检测SNPs主要技术的原理和方法,同时还根据一些具体实例介绍了SNPs在模式动、植物遗传图谱构建、品种鉴定、物种起源与亲缘关系、连锁不平衡与关联分析及其在群体遗传结构及其变化机制研究中的应用。最后展望了SNPs在群体遗传、分子育种和生物进化等研究领域中的应用前景。