Large SNP arrays for genotyping in crop plants.

Genotyping with large numbers of molecular markers is now an indispensable tool within plant genetics and breeding. Especially through the identification of large numbers of single nucleotide polymorphism (SNP) markers using the novel high-throughput sequencing technologies, it is now possible to reliably identify many thousands of SNPs at many different loci in a given plant genome. For a number of important crop plants, SNP markers are now being used to design genotyping arrays containing thousands of markers spread over the entire genome and to analyse large numbers of samples. In this article, we discuss aspects that should be considered during the design of such large genotyping arrays and the analysis of individuals. The fact that crop plants are also often autopolyploid or allopolyploid is given due consideration. Furthermore, we outline some potential applications of large genotyping arrays including high-density genetic mapping, characterization (fingerprinting) of genetic material and breeding-related aspects such as association studies and genomic selection.

Development of a large SNP genotyping array and generation of high-density genetic maps in tomato.

The concurrent development of high-throughput genotyping platforms and next generation sequencing (NGS) has increased the number and density of genetic markers, the efficiency of constructing detailed linkage maps, and our ability to overlay recombination and physical maps of the genome. We developed an array for tomato with 8,784 Single Nucleotide Polymorphisms (SNPs) mainly discovered based on NGS-derived transcriptome sequences. Of the SNPs, 7,720 (88%) passed manufacturing quality control and could be scored in tomato germplasm. The array was used to generate high-density linkage maps for three interspecific F(2) populations: EXPEN 2000 (Solanum lycopersicum LA0925 x S. pennellii LA0716, 79 individuals), EXPEN 2012 (S. lycopersicum Moneymaker x S. pennellii LA0716, 160 individuals), and EXPIM 2012 (S. lycopersicum Moneymaker x S. pimpinellifolium LA0121, 183 individuals). The EXPEN 2000-SNP and EXPEN 2012 maps consisted of 3,503 and 3,687 markers representing 1,076 and 1,229 unique map positions (genetic bins), respectively. The EXPEN 2000-SNP map had an average marker bin interval of 1.6 cM, while the EXPEN 2012 map had an average bin interval of 0.9 cM. The EXPIM 2012 map was constructed with 4,491 markers (1,358 bins) and an average bin interval of 0.8 cM. All three linkage maps revealed an uneven distribution of markers across the genome. The dense EXPEN 2012 and EXPIM 2012 maps showed high levels of colinearity across all 12 chromosomes, and also revealed evidence of small inversions between LA0716 and LA0121. Physical positions of 7,666 SNPs were identified relative to the tomato genome sequence. The genetic and physical positions were mostly consistent. Exceptions were observed for chromosomes 3, 10 and 12. Comparing genetic positions relative to physical positions revealed that genomic regions with high recombination rates were consistent with the known distribution of euchromatin across the 12 chromosomes, while very low recombination rates were observed in the heterochromatic regions.

A large maize (Zea mays L.) SNP genotyping array: development and germplasm genotyping, and genetic mapping to compare with the B73 reference genome.

SNP genotyping arrays have been useful for many applications that require a large number of molecular markers such as high-density genetic mapping, genome-wide association studies (GWAS), and genomic selection. We report the establishment of a large maize SNP array and its use for diversity analysis and high density linkage mapping. The markers, taken from more than 800,000 SNPs, were selected to be preferentially located in genes and evenly distributed across the genome. The array was tested with a set of maize germplasm including North American and European inbred lines, parent/F1 combinations, and distantly related teosinte material. A total of 49,585 markers, including 33,417 within 17,520 different genes and 16,168 outside genes, were of good quality for genotyping, with an average failure rate of 4% and rates up to 8% in specific germplasm. To demonstrate this array's use in genetic mapping and for the independent validation of the B73 sequence assembly, two intermated maize recombinant inbred line populations - IBM (B73×Mo17) and LHRF (F2×F252) - were genotyped to establish two high density linkage maps with 20,913 and 14,524 markers respectively. 172 mapped markers were absent in the current B73 assembly and their placement can be used for future improvements of the B73 reference sequence. Colinearity of the genetic and physical maps was mostly conserved with some exceptions that suggest errors in the B73 assembly. Five major regions containing non-colinearities were identified on chromosomes 2, 3, 6, 7 and 9, and are supported by both independent genetic maps. Four additional non-colinear regions were found on the LHRF map only; they may be due to a lower density of IBM markers in those regions or to true structural rearrangements between lines. Given the array's high quality, it will be a valuable resource for maize genetics and many aspects of maize breeding.

Analysis of DNA polymorphisms in sugar beet (Beta vulgaris L.) and development of an SNP-based map of expressed genes.

A panel of 13 sugar beet lines and one genotype each of the Beta vulgaris cultivars red beet and Swiss chard, and B. vulgaris ssp. maritima were used to identify polymorphisms in alignments of genomic DNA sequences derived from 315 EST- and 43 non-coding RFLP-derived loci. In sugar beet lines, loci of expressed genes showed an average SNP frequency of 1/72 bp, 1 in 58 bp in non-coding sequences, increasing to 1/47 bp upon the addition of the remaining genotypes. Within analysed DNA fragments, alleles at different SNP positions displayed linkage disequilibrium indicative of haplotype structures. On average 2.7 haplotypes were found in sugar beet lines, and haplotype conservation in expressed genes appeared to exceed 500 bp in length. Seven different genotyping techniques including SNP detection by MALDI-TOF mass spectrometry, pyrosequencing and fluorescence scanning of labelled nucleotides were employed to perform 712 segregation analyses for 538 markers in three F(2) populations. Functions were predicted for 492 mapped sequences. Genetic maps comprised 305 loci covering 599.8 cM in population K1, 241 loci distributed over 636.6 cM in population D2, and 166 loci over 507.1 cM in population K2, respectively. Based on 156 markers common to more than one population an integrated map was constructed with 524 loci covering 664.3 cM. For 377 loci the genome positions of the most similar sequences from A. thaliana were identified, but little evidence for previously presented ancestral genome structures was found.

Minisequencing on functionalised self-assembled monolayer as a simple approach for single nucleotide polymorphism analysis of cattle.

We have developed a genetic barcode module, based on a parallel sorting facility of single nucleotide polymorphism for secure individual identification of cattle. Biotinylated allele-specific oligonucleotides were immobilized onto the predefined spots of streptavidin tethered self-assembled monolayers with long chain alkanethiols on biochips. The target DNAs for hybridization and subsequent on-chip minisequencing were produced by multiplex PCR method. After enzymatic extension, only the moiety-modified dideoxynucleotide triphosphate, when coupled to its complementary target sequence, could be detected by the corresponding antibody to the moiety in a specific and sensitive manner. The database SNPZoo was developed for storage of the sequence information consisting of cytosine/thymidine patterns This SNP chip system can further be used in the detection of any replaceable point mutations occurring in the human and animal genes.

Association of a lysine-232/alanine polymorphism in a bovine gene encoding acyl-CoA:diacylglycerol acyltransferase (DGAT1) with variation at a quantitative trait locus for milk fat content.

DGAT1 encodes diacylglycerol O-acyltransferase (EC ), a microsomal enzyme that catalyzes the final step of triglyceride synthesis. It became a functional candidate gene for lactation traits after studies indicated that mice lacking both copies of DGAT1 are completely devoid of milk secretion, most likely because of deficient triglyceride synthesis in the mammary gland. Our mapping studies placed DGAT1 close to the region of a quantitative trait locus (QTL) on bovine chromosome 14 for variation in fat content of milk. Sequencing of DGAT1 from pooled DNA revealed significant frequency shifts at several variable positions between groups of animals with high and low breeding values for milk fat content in different breeds (Holstein-Friesian, Fleckvieh, and Braunvieh). Among the variants was a nonconservative substitution of lysine by alanine (K232A), with the lysine-encoding allele being associated with higher milk fat content. Haplotype analysis indicated the lysine variant to be ancestral. Two animals that were typed heterozygous (Qq) at the QTL based on marker-assisted QTL-genotyping were heterozygous for the K232A substitution, whereas 14 animals that are most likely qq at the QTL were homozygous for the alanine-encoding allele. An independent association study in Fleckvieh animals confirmed the positive effect of the lysine variant on milk fat content. We consider the nonconservative K232A substitution to be directly responsible for the QTL variation, although our genetic studies cannot provide formal proof.