KASP Assays Research Articles

Susceptibility loci for sporadic brain arteriovenous malformation; a replication study and meta-analysis

BackgroundCase–control studies have reported multiple genetic loci to be associated with sporadic brain arteriovenous malformations (AVMs) but most of these have not been replicated in independent populations. The aim of...

Fine mapping of Msv1, a major QTL for resistance to Maize Streak Virus leads to development of production markers for breeding pipelines.

Msv1 , the major QTL for MSV resistance was delimited to an interval of 0.87 cM on chromosome 1 at 87 Mb and production markers with high prediction accuracy were developed. Maize streak virus (MSV) disease is a devastating disease in the Sub-Saharan Africa (SSA), which causes significant yield loss in maize. Resistance to MSV has previously been mapped to a major QTL (Msv1) on chromosome 1 that is germplasm and environment independent and to several minor loci elsewhere in the genome. In this study, Msv1 was fine-mapped through QTL isogenic recombinant strategy using a large F 2 population of CML206 × CML312 to an interval of 0.87 cM on chromosome 1. Genome-wide association study was conducted in the DTMA (Drought Tolerant Maize for Africa)-Association mapping panel with 278 tropical/sub-tropical breeding lines from CIMMYT using the high-density genotyping-by-sequencing (GBS) markers. This study identified 19 SNPs in the region between 82 and 93 Mb on chromosome 1(B73 RefGen_V2) at a P < 1.00E-04, which coincided with the fine-mapped region of Msv1. Haplotype trend regression identified a haplotype block significantly associated with response to MSV. Three SNPs in this haplotype block at 87 Mb on chromosome 1 had an accuracy of 0.94 in predicting the disease reaction in a collection of breeding lines with known responses to MSV infection. In two biparental populations, selection for resistant Msv1 haplotype demonstrated a reduction of 1.03-1.39 units on a rating scale of 1-5, compared to the susceptible haplotype. High-throughput KASP assays have been developed for these three SNPs to enable routine marker screening in the breeding pipeline for MSV resistance.

A SNP-based consensus genetic map for synteny-based trait targeting in faba bean (Vicia faba L.).

SummaryFaba bean (Vicia faba L.) is a globally important nitrogen‐fixing legume, which is widely grown in a diverse range of environments. In this work, we mine and validate a set of 845 SNPs from the aligned transcriptomes of two contrasting inbred lines. Each V. faba SNP is assigned by BLAST analysis to a single Medicago orthologue. This set of syntenically anchored polymorphisms were then validated as individual KASP assays, classified according to their informativeness and performance on a panel of 37 inbred lines, and the best performing 757 markers used to genotype six mapping populations. The six resulting linkage maps were merged into a single consensus map on which 687 SNPs were placed on six linkage groups, each presumed to correspond to one of the six V. faba chromosomes. This sequence‐based consensus map was used to explore synteny with the most closely related crop species, lentil and the most closely related fully sequenced genome, Medicago. Large tracts of uninterrupted colinearity were found between faba bean and Medicago, making it relatively straightforward to predict gene content and order in mapped genetic interval. As a demonstration of this, we mapped a flower colour gene to a 2‐cM interval of Vf chromosome 2 which was highly colinear with Mt3. The obvious candidate gene from 78 gene models in the collinear Medicago chromosome segment was the previously characterized MtWD40‐1 gene controlling anthocyanin production in Medicago and resequencing of the Vf orthologue showed a putative causative deletion of the entire 5′ end of the gene.

Mapping resistance to the Ug99 race group of the stem rust pathogen in a spring wheat landrace.

A new gene for Ug99 resistance from wheat landrace PI 374670 was detected on the long arm of chromosome 7A. Wheat landrace PI 374670 has seedling and field resistance to stem rust caused by Puccinia graminis f. sp tritici Eriks. & E. Henn (Pgt) race TTKSK. To elucidate the inheritance of resistance, 216 BC1F2 families, 192 double haploid (DH) lines, and 185 recombinant inbred lines (RILs) were developed by crossing PI 374670 and the susceptible line LMPG-6. The parents and progeny were evaluated for seedling resistance to Pgt races TTKSK, MCCFC, and TPMKC. The DH lines were tested in field stem rust nurseries in Kenya and Ethiopia. The DH lines were genotyped with the 90K wheat iSelect SNP genotyping platform. Goodness-of-fit tests indicated that a single dominant gene in PI 374670 conditioned seedling resistance to the three Pgt races. The seedling resistance locus mapped to the long arm of chromosome 7A and this result was verified in the RIL population screened with the flanking SNP markers using KASP assays. In the same region, a major QTL for field resistance was detected in a 7.7 cM interval and explained 34-54 and 29-36% of the variation in Kenya and Ethiopia, respectively. Results from tests with specific Pgt races and the csIH81 marker showed that the resistance was not due to Sr22. Thus, a new stem rust resistance gene or allele, either closely linked or allelic to Sr15, is responsible for the seedling and field resistance of PI 374670 to Ug99.

Comprehensive analysis of the association of EGFR, CALM3 and SMARCD1 gene polymorphisms with BMD in Caucasian women.

SummaryThree genes, including EGFR (epidermal growth factor receptor), CALM3 (calmodulin 3, calcium-modulated protein 3) and SMARCD1 (SWI/SNF-related matrix-associated actin-dependent regulator of chromatin subfamily d member 1), play different roles in bone and/or fat metabolism in Caucasian women. In this population-based investigation of 870 unrelated postmenopausal Caucasian women, CALM3 polymorphisms were significantly associated with femoral neck bone mineral density (FNK BMD), hip BMD and spine BMD. Age and tobacco status also affected BMD levels and were therefore corrected for in our statistical analysis.Introduction EGFR, CALM3 and SMARCD1 play roles in bone and/or fat metabolism. However, the correlations between the polymorphisms of these three genes and body composition levels, including BMD, remain to be determined.Materials and MethodsA population-based investigation of 870 white women was conducted. Forty-four SNPs (single nucleotide polymorphisms) in EGFR, CALM3 and SMARCD1 were chosen by the software, including those of potential functional importance. The candidate SNPs were genotyped by the KASPar assay for an association analysis with body composition levels. The correlation analysis was assessed by the Pearson's product-moment correlation coefficient and Spearman rank-order correlation tests, and the family-wise error was corrected using the Wald test implemented in PLINK.ResultsThe SNP rs12461917 in the 3′-flanking region of the CALM3 gene was significantly associated with FNK BMD (P = 0.001), hip BMD (P<0.001) and spine BMD (P = 0.001); rs11083838 in the 5′-flanking region of CALM3 gene was associated with spine BMD (P = 0.009). After adjusting for multiple comparisons, rs12461917 remained significant (P-adjusted = 0.033 for FNK BMD, P-adjusted = 0.006 for hip BMD and P-adjusted = 0.018 for spine BMD).ConclusionsOur data show that polymorphisms of the CALM3 gene in Caucasian women may contribute to variations in the BMD of the hip, spine and femoral neck.

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

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

Resistance to Clover yellow vein virus in Common Bean Germplasm

ABSTRACTClover yellow vein virus (ClYVV) is a potentially devastating viral pathogen of temperate common bean (Phaseolus vulgaris L.) production. Recent research demonstrated that resistance to ClYVV is conditioned by bc‐3 and bc‐32 and that a nonsynonymous mutation in P. vulgaris eukaryotic translation initiation factor 4E (PveIF4E) is the putative functional determinant of this resistance. The objective of this research was to expand the characterization of ClYVV resistance in common bean through greenhouse‐based phenotypic evaluation of 391 accessions of the USDA P. vulgaris Core Collection, a collection of 99 snap bean cultivars and breeding lines, and an extended host differential panel of 75 genotypes with known response to other viral pathogens and to determine if previously designed Kompetitive Allele Specific Polymerase Chain Reaction (KASP) assays could be validated for future bc‐3 resistance allele mining. As a result, one accession from the Core Collection, one additional snap bean breeding line, and 19 dry bean genotypes representing most major market classes of common bean were identified to be resistant. All genotypes that remained symptomless throughout the evaluation were demonstrated to possess mutations in PveIF4E that corresponded with either bc‐3 or bc‐32. In addition, the bc‐u bc‐22 allele combination and the By‐2 allele were also demonstrated to condition resistance phenotypes to ClYVV. These results synthesize a new model for genetic resistance to ClYVV and validate the use of the KASP assays for large‐scale bc‐3 allele mining.

A major quantitative trait locus conferring adult plant partial resistance to crown rust in oat.

BackgroundCrown rust, caused by Puccinia coronata f. sp. avenae, is the most important disease of oat worldwide. Adult plant resistance (APR), based upon partial resistance, has proven to be a durable rust management strategy in other cereal rust pathosystems. The crown rust APR in the oat line MN841801 has been effective for more than 30 years. The genetic basis of this APR was studied under field conditions in three recombinant inbred line (RIL) populations: 1) AC Assiniboia/MN841801, 2) AC Medallion/MN841801, and 3) Makuru/MN841801. The populations were evaluated for crown rust resistance with the crown rust isolate CR251 (race BRBB) in multiple environments. The 6 K oat and 90 K wheat Illumina Infinium single nucleotide polymorphism (SNP) arrays were used for genotyping the AC Assiniboia/MN841801 population. KASP assays were designed for selected SNPs and genotyped on the other two populations.ResultsThis study reports a high density genetic linkage map constructed with oat and wheat SNP markers in the AC Assiniboia/MN841801 RIL population. Most wheat SNPs were monomorphic in the oat population. However the polymorphic wheat SNPs could be scored accurately and integrated well into the linkage map. A major quantitative trait locus (QTL) on oat chromosome 14D, designated QPc.crc-14D, explained up to 76% of the APR phenotypic variance. This QTL is flanked by two SNP markers, GMI_GBS_90753 and GMI_ES14_c1439_83. QPc.crc-14D was validated in the populations AC Medallion/MN841801 and Makuru/MN841801.ConclusionsWe report the first APR QTL in oat with a large and consistent effect. QPc.crc-14D was statistically significant in all environments tested in each of the three oat populations. QPc.crc-14D is a suitable candidate for use in marker-assisted breeding and also an excellent target for map-based cloning. This is also the first study to use the 90 K wheat Infinium SNP array on oat for marker development and comparative mapping. The Infinium SNP array is a useful tool for saturating oat maps with markers. Synteny with wheat suggests that QPc.crc-14D is orthologous with the stripe rust APR gene Yr16 in wheat.Electronic supplementary materialThe online version of this article (doi:10.1186/s12870-014-0250-2) contains supplementary material, which is available to authorized users.

Development of a SNP marker assay for the Lr67 gene of wheat using a genotyping by sequencing approach

The recently characterised and mapped adult plant resistance gene Lr67 contributes to the breeding of durable resistance against multiple fungal pathogens of wheat. Lack of tightly linked markers has so far prevented widespread use of Lr67 in breeding programmes. Exhaustive screening of available SSR markers and 9,000 gene-associated SNPs using bulk segregant analysis yielded several linked markers but none close enough to the gene for routine breeding applications. In contrast, genotyping by sequencing identified several better-linked SNP-based markers that were converted into single-marker KASPar assays. Marker csSNP856 co-segregated with Lr67 in two mapping populations and was significantly associated with stripe rust resistance in an association mapping population of 972 wheat lines. High frequency (43 %, 29/66) of the Lr67-associated SNP allele in landraces from the Indian subcontinent suggests the origin of Lr67 from this region. While marker csSNP856 is not diagnostic for the presence of Lr67, it represents a significant step towards assisting the transfer of this adult plant resistance gene into elite germplasm via marker-assisted breeding.

Evaluation of genetic risk loci for intracranial aneurysms in sporadic arteriovenous malformations of the brain

BackgroundIn genome-wide association studies (GWAS) five putative risk loci are associated with intracranial aneurysm. As brain arteriovenous malformations (AVM) and intracranial aneurysms are both intracranial vascular diseases and AVMs often...

SNP Discovery for mapping alien introgressions in wheat

BackgroundMonitoring alien introgressions in crop plants is difficult due to the lack of genetic and molecular mapping information on the wild crop relatives. The tertiary gene pool of wheat is a very important source of genetic variability for wheat improvement against biotic and abiotic stresses. By exploring the 5Mg short arm (5MgS) of Aegilops geniculata, we can apply chromosome genomics for the discovery of SNP markers and their use for monitoring alien introgressions in wheat (Triticum aestivum L).ResultsThe short arm of chromosome 5Mg of Ae. geniculata Roth (syn. Ae. ovata L.; 2n = 4x = 28, UgUgMgMg) was flow-sorted from a wheat line in which it is maintained as a telocentric chromosome. DNA of the sorted arm was amplified and sequenced using an Illumina Hiseq 2000 with ~45x coverage. The sequence data was used for SNP discovery against wheat homoeologous group-5 assemblies. A total of 2,178 unique, 5MgS-specific SNPs were discovered. Randomly selected samples of 59 5MgS-specific SNPs were tested (44 by KASPar assay and 15 by Sanger sequencing) and 84% were validated. Of the selected SNPs, 97% mapped to a chromosome 5Mg addition to wheat (the source of t5MgS), and 94% to 5Mg introgressed from a different accession of Ae. geniculata substituting for chromosome 5D of wheat. The validated SNPs also identified chromosome segments of 5MgS origin in a set of T5D-5Mg translocation lines; eight SNPs (25%) mapped to TA5601 [T5DL · 5DS-5MgS(0.75)] and three (8%) to TA5602 [T5DL · 5DS-5MgS (0.95)]. SNPs (gsnp_5ms83 and gsnp_5ms94), tagging chromosome T5DL · 5DS-5MgS(0.95) with the smallest introgression carrying resistance to leaf rust (Lr57) and stripe rust (Yr40), were validated in two released germplasm lines with Lr57 and Yr40 genes.ConclusionThis approach should be widely applicable for the identification of species/genome-specific SNPs. The development of a large number of SNP markers will facilitate the precise introgression and monitoring of alien segments in crop breeding programs and further enable mapping and cloning novel genes from the wild relatives of crop plants.

Genetic Diversity and Demographic History of Cajanus spp. Illustrated from Genome-Wide SNPs

Understanding genetic structure of Cajanus spp. is essential for achieving genetic improvement by quantitative trait loci (QTL) mapping or association studies and use of selected markers through genomic assisted breeding and genomic selection. After developing a comprehensive set of 1,616 single nucleotide polymorphism (SNPs) and their conversion into cost effective KASPar assays for pigeonpea (Cajanus cajan), we studied levels of genetic variability both within and between diverse set of Cajanus lines including 56 breeding lines, 21 landraces and 107 accessions from 18 wild species. These results revealed a high frequency of polymorphic SNPs and relatively high level of cross-species transferability. Indeed, 75.8% of successful SNP assays revealed polymorphism, and more than 95% of these assays could be successfully transferred to related wild species. To show regional patterns of variation, we used STRUCTURE and Analysis of Molecular Variance (AMOVA) to partition variance among hierarchical sets of landraces and wild species at either the continental scale or within India. STRUCTURE separated most of the domesticated germplasm from wild ecotypes, and separates Australian and Asian wild species as has been found previously. Among Indian regions and states within regions, we found 36% of the variation between regions, and 64% within landraces or wilds within states. The highest level of polymorphism in wild relatives and landraces was found in Madhya Pradesh and Andhra Pradesh provinces of India representing the centre of origin and domestication of pigeonpea respectively.

Single Nucleotide Polymorphism–based Genetic Diversity in the Reference Set of Peanut (Arachis spp.) by Developing and Applying Cost‐Effective Kompetitive Allele Specific Polymerase Chain Reaction Genotyping Assays

Kompetitive allele‐specific polymerase chain reaction (KASP) assays have emerged as cost‐effective marker assays especially for molecular breeding applications. Therefore, a set of 96 informative single nucleotide polymorphisms (SNPs) was used to develop KASP assays in groundnut or peanut (Arachis spp.). Developed assays were designated as groundnut KASP assay markers (GKAMs) and screened on 94 genotypes (validation set) that included parental lines of 27 mapping populations, seven synthetic autotetraploid and amphidiploid lines, and 19 wild species accessions. As a result, 90 GKAMs could be validated and 73 GKAMs showed polymorphism in the validation set. Validated GKAMs were screened on 280 diverse genotypes of the reference set for estimating diversity features and elucidating genetic relationships. Cluster analysis of marker allelic data grouped accessions according to their genome type, subspecies, and botanical variety. The subspecies Arachis hypogaea L. subsp. fastigiata Waldron and A. hypogaea subsp. hypogaea formed distinct cluster; however, some overlaps were found indicating their frequent intercrossing during the course of evolution. The wild species, having diploid genomes, were grouped into a single cluster. The average polymorphism information content value for polymorphic GKAMs was 0.32 in the validation set and 0.31 in the reference set. These validated and highly informative GKAMs may be useful for genetics and breeding applications in Arachis species.

Development of single nucleotide polymorphisms in Phaseolus vulgaris and related Phaseolus spp

In this study, new single nucleotide polymorphism (SNP) markers were developed for common bean (Phaseolus vulgaris L.) and related Phaseolus species. The applied strategy presents new and interesting aspects, such as the choice of accessions used, which was aimed at capturing a large portion of the genetic diversity present in the common bean, with particular focus on wild and domesticated materials from Mesoamerica and the identification of loci for sequencing. Indeed, the primer pairs for 34 loci were designed with the main strategy being to search for single-copy orthologous genes among the legumes (for use in other legume species and comparative analyses). The 10 remaining loci were selected as being near to domestication quantitative trait loci or detected as putatively under selection during domestication in previous studies. To provide an efficient and inexpensive genotyping platform for geneticists and breeders, we used sequence data to develop 60 new SNP markers for KASPar assay genotyping. The same sample was also genotyped with SNP markers developed for common bean in other studies for the same assay. This allowed testing for systematic bias according to the criteria chosen to select the genotypes in which the genetic diversity is surveyed during SNP discovery. Finally, we show that most of the SNP markers worked well in a set of accessions of other species belonging to the Phaseolus genus. The genetic resources developed will be very useful not only for breeding, but also for biodiversity conservation management and evolutionary studies on legumes.

A series of eIF4E alleles at the Bc-3 locus are associated with recessive resistance to Clover yellow vein virus in common bean

Clover yellow vein virus (ClYVV) is capable of causing severe damage to common bean (Phaseolus vulgaris L.) production worldwide. The snap bean market class is particularly vulnerable because infection may lead to distortion and necrosis of the fresh green pods and rejection of the harvest. Three putatively independent recessive genes (cyv, desc, bc-3) have been reported to condition resistance to ClYVV; however, their allelic relationships have not been resolved. We identified, evaluated, and characterized the phenotypic and molecular genetic variation present in 21 informative common bean genotypes for resistance to ClYVV. Allelism testing phenotypes from multiple populations provided clear evidence that the three genes were a series of recessive alleles at the Bc-3 locus that condition unique potyvirus strain- and species-specific resistance spectra. Candidate gene analysis revealed complete association between the recessive resistance alleles and unique patterns of predicted amino acid substitutions in P. vulgaris eukaryotic translation initiation factor 4E (PveIF4E). This led to the discovery and characterization of two novel PveIF4E alleles associated with resistance to ClYVV, PveIF4E (3) , and PveIF4E (4) . We developed KASPar allele-specific SNP genotyping assays and demonstrated their ability to accurately detect and differentiate all of the PveIF4E haplotypes present in the germplasm, allelism testing, and in three separate segregating populations. The results contribute to an enhanced understanding and accessibility of the important potyvirus resistance conditioned by recessive alleles at Bc-3. The KASPar assays should be useful to further enable germplasm exploration, allelic discrimination, and marker-assisted introgression of bc-3 alleles in common bean.

Genetic structure and domestication of carrot (Daucus carota subsp. sativus) (Apiaceae)

Analyses of genetic structure and phylogenetic relationships illuminate the origin and domestication of modern crops. Despite being an important worldwide vegetable, the genetic structure and domestication of carrot (Daucus carota) is poorly understood. We provide the first such study using a large data set of molecular markers and accessions that are widely dispersed around the world. • Sequencing data from the carrot transcriptome were used to develop 4000 single nucleotide polymorphisms (SNPs). Eighty-four genotypes, including a geographically well-distributed subset of wild and cultivated carrots, were genotyped using the KASPar assay. • Analysis of allelic diversity of SNP data revealed no reduction of genetic diversity in cultivated vs. wild accessions. Structure and phylogenetic analysis indicated a clear separation between wild and cultivated accessions as well as between eastern and western cultivated carrot. Among the wild carrots, those from Central Asia were genetically most similar to cultivated accessions. Furthermore, we found that wild carrots from North America were most closely related to European wild accessions. • Comparing the genetic diversity of wild and cultivated accessions suggested the absence of a genetic bottleneck during carrot domestication. In conjunction with historical documents, our results suggest an origin of domesticated carrot in Central Asia. Wild carrots from North America were likely introduced as weeds with European colonization. These results provide answers to long-debated questions of carrot evolution and domestication and inform germplasm curators and breeders on genetic substructure of carrot genetic resources.

Fine mapping and identification of candidate genes controlling the resistance to southern root-knot nematode in PI 96354

Meloidogyne incognita (Kofoid and White) Chitwood (Mi) is the most economically damaging species of the root-knot nematode to soybean and other crops in the southern USA. PI 96354 was identified to carry a high level of resistance to galling and Mi egg production. Two Quantitative Trait Locus(QTLs) were found to condition the resistance in PI 96354 including a major QTL and a minor QTL on chromosome 10 and chromosome 18, respectively. To fine map the major QTL on chromosome 10, F5:6 recombinant inbred lines from the cross between PI 96354 and susceptible genotype Bossier were genotyped with Simple Sequence Repeats(SSR) markers to identify recombinational events. Analysis of lines carrying key recombination events placed the Mi-resistant allele on chromosome 10 to a 235-kb region of the 'Williams 82' genome sequence with 30 annotated genes. Candidate gene analysis identified four genes with cell wall modification function that have several mutations in promoter, exon, 5', and 3'UTR regions. qPCR analysis showed significant difference in expression levels of these four genes in Bossier compared to PI 96354 in the presence of Mi. Thirty Mi-resistant soybean lines were found to have same SNPs in these 4 candidate genes as PI 96354 while 12 Mi-susceptible lines possess the 'Bossier' genotype. The mutant SNPs were used to develop KASP assays to detect the resistant allele on chromosome 10. The four candidate genes identified in this study can be used in further studies to investigate the role of cell wall modification genes in conferring Mi resistance in PI 96354.

4.2 Variation in Beta Defensin 1 Genotype is Associated with Preterm Birth

IntroductionHuman beta defensin 1 (HBD1) is an antimicrobial and immunomodulatory peptide present in cervical mucus. Variation in cervical antimicrobial expression is associated with preterm labour. We hypothesised that SNPs in...

Assessment of Genetic Diversity in Peruvian Amaranth (Amaranthus caudatus and A. hybridus) Germplasm using Single Nucleotide Polymorphism Markers

ABSTRACTAmaranth (Amaranthus caudatus L.) is an important seed crop in the Andes. It has excellent nutritional value and ample capacity for growth under drought, heat, and soil nutrient deficiency. The objective of this study was to characterize and quantify genetic diversity among a series of 178 mostly Peruvian amaranth individuals genotyped from among 48 accessions using 96 biallelic single nucleotide polymorphism (SNP) markers in KASPar assays on a 96.96 Fluidigm EP1 platform. The 178 seed genotypes included the predominant white‐opaque (waxy), white‐vitreous (nonwaxy), black, and brown types, all from within the Centro de Investigacion de Cultivos Andinos‐Universidad Nacional de San Antonio Abad del Cusco collection. The samples separated into two groups based on the seed color parameter (black, brown, and white‐vitreous versus white‐opaque) with STRUCTURE analysis but did not form distinct genotypic groups based on geographic origin. When analyzed as seed‐color subgroups—black vs. brown vs. white‐vitreous vs. white‐opaque—83.46% of the genetic variation was explained by genetic differences within subgroup. The highest average observed heterozygosity (Ho) within color‐based subgroups was 0.2612 in the vitreous‐seeded subgroup, and within‐group expected heterozygosity (He) was highest in the brown‐seeded subgroup (0.398). Cluster analysis and principal component analysis partitioned the amaranth accessions into two clusters with four distinct subclusters. Analyses also revealed evidence for selection of nonwaxy, vitreous‐seeded amaranth in the Andes.

SNP discovery using Paired‐End RAD‐tag sequencing on pooled genomic DNA of Sisymbrium austriacum (Brassicaceae)

Single nucleotide polymorphisms SNPs are rapidly replacing anonymous markers in population genomic studies, but their use in non model organisms is hampered by the scarcity of cost-effective approaches to uncover genome-wide variation in a comprehensive subset of individuals. The screening of one or only a few individuals induces ascertainment bias. To discover SNPs for a population genomic study of the Pyrenean rocket (Sisymbrium austriacum subsp. chrysanthum), we undertook a pooled RAD-PE (Restriction site Associated DNA Paired-End sequencing) approach. RAD tags were generated from the PstI-digested pooled genomic DNA of 12 individuals sampled across the species distribution range and paired-end sequenced using Illumina technology to produce ~24.5Mb of sequences, covering ~7% of the specie's genome. Sequences were assembled into ~76000 contigs with a mean length of 323bp (N(50) =357bp, sequencing depth=24x). In all, >15000 SNPs were called, of which 47% were annotated in putative genic regions based on homology with the Arabidopsis thaliana genome. Gene ontology (GO) slim categorization demonstrated that the identified SNPs covered extant genic variation well. The validation of 300 SNPs on a larger set of individuals using a KASPar assay underpinned the utility of pooled RAD-PE as an inexpensive genome-wide SNP discovery technique (success rate: 87%). In addition to SNPs, we discovered >600 putative SSR markers.