Segregation distortion in Arabidopsis C24/Col-0 and Col-0/C24 recombinant inbred line populations is due to reduced fertility caused by epistatic interaction of two loci

A new large set of reciprocal recombinant inbred lines (RILs) was created between the Arabidopsis accessions Col-0 and C24 for quantitative trait mapping approaches, consisting of 209 Col-0 x C24 and 214 C24 x Col-0 F(7 )RI lines. Genotyping was performed using 110 evenly distributed framework single nucleotide polymorphism markers, yielding a genetic map of 425.70 cM, with an average interval of 3.87 cM. Segregation distortion (SD) was observed in several genomic regions during the construction of the genetic map. Linkage disequilibrium analysis revealed an association between a distorted region at the bottom of chromosome V and a non-distorted region on chromosome IV. A detailed analysis of the RILs for these two regions showed that an SD occurred when homozygous Col-0 alleles on chromosome IV coincided with homozygous C24 alleles at the bottom of chromosome V. Using nearly isogenic lines segregating for the distorted region we confirmed that this genotypic composition leads to reduced fertility and fitness.     

Abscisic acid antagonizes ethylene-induced hyponastic growth in Arabidopsis

Ethylene induces enhanced differential growth in petioles of Arabidopsis (Arabidopsis thaliana), resulting in an upward movement of the leaf blades (hyponastic growth). The amplitude of this effect differs between accessions, with Columbia-0 (Col-0) showing a large response, while in Landsberg erecta (Ler), hyponastic growth is minimal. Abscisic acid (ABA) was found to act as an inhibitory factor of this response in both accessions, but the relationship between ethylene and ABA differed between the two; the ability of ABA to inhibit ethylene-induced hyponasty was significantly more pronounced in Col-0. Mutations in ABI1 or ABI3 induced a strong ethylene-regulated hyponastic growth in the less responsive accession Ler, while the response was abolished in the ABA-hypersensitive era1 in Col-0. Modifications in ABA levels altered petiole angles in the absence of applied ethylene, indicating that ABA influences petiole angles also independently from ethylene. A model is proposed whereby the negative effect of ABA on hyponastic growth is overcome by ethylene in Col-0 but not in Ler. However, when ABA signaling is artificially released in Ler, this regulatory mechanism is bypassed, resulting in a strong hyponastic response in this accession.     

Using Next Generation Sequencing for Multiplexed Trait-Linked Markers in Wheat

With the advent of next generation sequencing (NGS) technologies, single nucleotide polymorphisms (SNPs) have become the major type of marker for genotyping in many crops. However, the availability of SNP markers for important traits of bread wheat ( Triticum aestivum L.) that can be effectively used in marker-assisted selection (MAS) is still limited and SNP assays for MAS are usually uniplex. A shift from uniplex to multiplex assays will allow the simultaneous analysis of multiple markers and increase MAS efficiency. We designed 33 locus-specific markers from SNP or indel-based marker sequences that linked to 20 different quantitative trait loci (QTL) or genes of agronomic importance in wheat and analyzed the amplicon sequences using an Ion Torrent Proton Sequencer and a custom allele detection pipeline to determine the genotypes of 24 selected germplasm accessions. Among the 33 markers, 27 were successfully multiplexed and 23 had 100% SNP call rates. Results from analysis of "kompetitive allele-specific PCR" (KASP) and sequence tagged site (STS) markers developed from the same loci fully verified the genotype calls of 23 markers. The NGS-based multiplexed assay developed in this study is suitable for rapid and high-throughput screening of SNPs and some indel-based markers in wheat.     

Use of VeraCode 384-plex assays for watermelon diversity analysis and integrated genetic map of watermelon with single nucleotide polymorphisms and simple sequence repeats

Watermelon (Citrullus lanatus var. lanatus) is one of the most important vegetable crops in the world. Molecular markers have become the tools of choice for resolving watermelon taxonomic relationships and evolution. Increased numbers of single nucleotide polymorphism (SNP) markers together with simple sequence repeat (SSR) markers would be useful for phylogenetic analyses of germplasm accessions and for linkage mapping for marker-assisted breeding with quantitative trait loci and single genes. We aimed to construct a genetic map based on SNPs (generated by Illumina Veracode multiplex assays for genotyping) and SSR markers and evaluate relationships inferred from SNP genotypes between 130 watermelon accessions collected throughout the world. We incorporated 282 markers (232 SNPs and 50 SSRs) into the linkage map. The genetic map consisted of 11 linkage groups spanning 924.72 cM with an average distance of 3.28 cM between markers. Because all of the SNP-containing sequences were assembled with the whole-genome sequence draft for watermelon, chromosome numbers could be readily assigned for all the linkage groups. We found that 134 SNPs were polymorphic in 130 watermelon accessions chosen for diversity studies. The current 384-plex SNP set is a powerful tool for characterizing genetic relatedness and for developing medium-resolution genetic maps.     

MSQT for choosing SNP assays from multiple DNA alignments

MOTIVATION: One challenging aspect of genotyping and association mapping projects is often the identification of markers that are informative between groups of individuals and to convert these into genotyping assays. RESULTS: The Multiple SNP Query Tool (MSQT) extracts SNP information from multiple sequence alignments, stores it in a database, provides a web interface to query the database and outputs SNP information in a format directly applicable for SNP-assay design. MSQT was applied to Arabidopsis thaliana sequence data to develop SNP genotyping assays that distinguish a recurrent parent (Col-0) from five other strains. SNPs with intermediate allele frequencies were also identified and developed into markers suitable for efficient genetic mapping among random pairs of wild strains. AVAILABILITY: The source code for MSQT is available at http://msqt.weigelworld.org, together with an online instance of MSQT containing data on 1214 sequenced fragments from 96 ecotypes (wild inbred strains) of the reference plant A. thaliana. All SNP genotyping assays are available in several formats for broad community use. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.     

SNP genotyping in melons: genetic variation, population structure, and linkage disequilibrium

Novel sequencing technologies were recently used to generate sequences from multiple melon (Cucumis melo L.) genotypes, enabling the in silico identification of large single nucleotide polymorphism (SNP) collections. In order to optimize the use of these markers, SNP validation and large-scale genotyping are necessary. In this paper, we present the first validated design for a genotyping array with 768 SNPs that are evenly distributed throughout the melon genome. This customized Illumina GoldenGate assay was used to genotype a collection of 74 accessions, representing most of the botanical groups of the species. Of the assayed loci, 91 % were successfully genotyped. The array provided a large number of polymorphic SNPs within and across accessions. This set of SNPs detected high levels of variation in accessions from this crop’s center of origin as well as from several other areas of melon diversification. Allele distribution throughout the genome revealed regions that distinguished between the two main groups of cultivated accessions (inodorus and cantalupensis). Population structure analysis showed a subdivision into five subpopulations, reflecting the history of the crop. A considerably low level of LD was detected, which decayed rapidly within a few kilobases. Our results show that the GoldenGate assay can be used successfully for high-throughput SNP genotyping in melon. Since many of the genotyped accessions are currently being used as the parents of breeding populations in various programs, this set of mapped markers could be used for future mapping and breeding efforts.     

Ecotypic variability in the metabolic response of seeds to diurnal hydration-dehydration cycles and its relationship to seed vigor

Seeds in the seed bank experience diurnal cycles of imbibition followed by complete dehydration. These conditions pose a challenge to the regulation of germination. The effect of recurring hydration-dehydration (Hy-Dh) cycles were tested on seeds from four Arabidopsis thaliana accessions [Col-0, Cvi, C24 and Ler]. Diurnal Hy-Dh cycles had a detrimental effect on the germination rate and on the final percentage of germination in Col-0, Cvi and C24 ecotypes, but not in the Ler ecotype, which showed improved vigor following the treatments. Membrane permeability measured by ion conductivity was generally increased following each Hy-Dh cycle and was correlated with changes in the redox status represented by the GSSG/GSH (oxidized/reduced glutathione) ratio. Among the ecotypes, Col-0 seeds displayed the highest membrane permeability, whilst Ler was characterized by the greatest increase in electrical conductivity following Hy-Dh cycles. Following Dh 2 and Dh 3, the respiratory activity of Ler seeds significantly increased, in contrast to the other ecotypes, indicative of a dramatic shift in metabolism. These differences were associated with accession-specific content and patterns of change of (i) cell wall-related laminaribiose and mannose; (ii) fatty acid composition, specifically of the unsaturated oleic acid and α-linoleic acid; and (iii) asparagine, ornithine and the related polyamine putrescine. Furthermore, in the Ler ecotype the content of the tricarboxylic acid (TCA) cycle intermediates fumarate, succinate and malate increased in response to dehydration, in contrast to a decrease in the other three ecotypes. These findings provide a link between seed respiration, energy metabolism, fatty acid β-oxidation, nitrogen mobilization and membrane permeability and the improved germination of Ler seeds following Hy-Dh cycles.     

Transgressive segregation reveals two Arabidopsis TIR-NB-LRR resistance genes effective against Leptosphaeria maculans, causal agent of blackleg disease

In a cross between the two resistant accessions Col-0 and Ler-0, a 15:1 segregation was found in F2, suggesting the presence of unlinked resistance loci to Leptosphaeria maculans. One hundred Col-4 x Ler-0, and 50 Ler-2 x Cvi-1 recombinant inbred lines, and seven susceptible Ler-0 x Ws-0 F2 progenies were examined to identify the two loci. Resistance in Col-4, Ws-0 and Cvi-1 (RLM1) was mapped to the marker m305 on chromosome 1. Col-4 x Ler-0 and Ler-2 x Cvi-1 mapping populations located RLM2(Ler) on the same arm of chromosome 4. A tight physical location of RLM2 was established through near-isogenic lines. This region was found to correspond to an ancient duplication event between the RLM1 and RLM2 loci. Two independent T-DNA mutants in a TIR-NB-LRR R gene (At1g64070) displayed susceptibility, and L. maculans susceptible mutant phenotypes were confirmed to be allelic for rlm1 in F1 after crosses with susceptible rlm1(Ler)rlm2(Col) plants. Complementation of rlm1(Ler)rlm2(Col) with the genomic Col-0 sequence of At1g64070 conferred resistance. In addition, two T-DNA mutants in a neighbouring homologous TIR-NB-LRR gene (At1g63880) displayed moderate susceptibility to L. maculans. Sequence analysis revealed that At1g64070 was truncated by a premature stop codon, and that At1g63880 was absent in Ler-0. RNA interference confirmed that Ler-0 resistance is dependent on genes structurally related to RLM1. Camalexin was identified as a quantitative co-dominant resistance factor of Col-0 origin, but independent of RLM1. RLM1/RLM2 resistance was, however, found to require RAR1 and partially HSP90.1.     

ERECTA, an LRR receptor-like kinase protein controlling development pleiotropically affects resistance to bacterial wilt

Bacterial wilt, one of the most devastating bacterial diseases of plants worldwide, is caused by Ralstonia solanacearum and affects many important crop species. We show that several strains isolated from solanaceous crops in Europe are pathogenic in different accessions of Arabidopsis thaliana. One of these strains, 14.25, causes wilting symptoms in A. thaliana accession Landsberg erecta (Ler) and no apparent symptoms in accession Columbia (Col-0). Disease development and bacterial multiplication in the susceptible Ler accession depend on functional hypersensitive response and pathogenicity (hrp) genes, key elements for bacterial pathogenicity. Genetic analysis using Ler x Col-0 recombinant inbred lines showed that resistance is governed by at least three loci: QRS1 (Quantitative Resistance to R. solanacearum) and QRS2 on chromosome 2, and QRS3 on chromosome 5. These loci explain about 90% of the resistance carried by the Col-0 accession. The ERECTA gene, which encodes a leucine-rich repeat receptor-like kinase (LRR-RLK) and affects development of aerial organs, is dimorphic in our population and lies close to QRS1. Susceptible Ler plants transformed with a wild-type ERECTA gene, and the LER line showed increased disease resistance to R. solanacearum as indicated by reduced wilt symptoms and impaired bacterial growth, suggesting unexpected cross-talk between resistance and developmental pathways.     

High-throughput genotyping in citrus accessions using an SNP genotyping array

We developed a 384 multiplexed SNP array, named CitSGA-1, for the genotyping of Citrus cultivars, and evaluated the performance and reliability of the genotyping. SNPs were surveyed by direct sequence comparison of the sequence tagged site (STS) fragment amplified from genomic DNA of cultivars representing the genetic diversity of citrus breeding in Japan. Among 1497 SNPs candidates, 384 SNPs for a high-throughput genotyping array were selected based on physical parameters of Illumina’s bead array criteria. The assay using CitSGA-1 was applied to a hybrid population of 88 progeny and 103 citrus accessions for breeding in Japan, which resulted in 73,726 SNP calls. A total of 351 SNPs (91 %) could call different genotypes among the DNA samples, resulting in a success rate for the assay comparable to previously reported rates for other plant species. To confirm the reliability of SNP genotype calls, parentage analysis was applied, and it indicated that the number of reliable SNPs and corresponding STSs were 276 and 213, respectively. The multiplexed SNP genotyping array reported here will be useful for the efficient construction of linkage map, for the detection of markers for marker-assisted breeding, and for the identification of cultivars.     

Construction and analysis of 2 reciprocal Arabidopsis introgression line populations

Two new large reciprocal sets of introgression lines (ILs) were created between the Arabidopsis accessions Col-0 and C24. In both sets (78 ILs with Col-0 background and 62 ILs with C24 background), the donor segments cover almost the entire genome with an average substitution size of 18.3 cM. In addition to the basic sets of ILs, further subILs were developed for 2 genomic regions allowing better mapping resolution. SubILs carrying donor segments with candidate genes for flowering time and reduced fertility were used to demonstrate the usefulness of the reciprocal ILs for quantitative trait loci detection and fine mapping. For subIL development at high resolution around the reduced fertility locus, we used modified CelI-based assays in one-well format for both marker development and genotyping. This serves as a very flexible and cost-effective approach.     

Developing Single Nucleotide Polymorphism (SNP) Markers for the Identification of Coffee Germplasm

Coffee is one of the most widely consumed beverages and represents a multibillion-dollar global industry. Accurate identification of coffee cultivars is essential for efficient management, exchange, and use of coffee genetic resources. To date, a universal platform that can allow data comparison across different laboratories and genotyping platforms has not been developed by the coffee research community. Using expressed sequence tags (EST) of Coffea arabica, C. canephora and C. racemosa from public databases, we developed 7538 single nucleotide polymorphism (SNP) markers and selected 180 for validation using 25 C. arabica and C. canephora accessions from Puerto Rico. Based on the validation result, we designated a panel of 55 SNP markers that are polymorphic across the two species. The average minor allele frequency and information index of this SNP panel are 0.281 and 0.690, respectively. This panel enabled the differentiation of all tested accessions of C. canephora, which accounts for 79.2 % of the total polymorphism in the samples. Only 21.8 % of the polymorphic SNPs were detected in the 12 C. arabica cultivars, which, nonetheless, were able to unambiguously differentiate the 12 Arabica cultivars into ten unique genotypes, including two synonymous groups. Several local Puerto Rican cultivars with partial Timor pedigree, including Limaní, Frontón, and TARS 18087, showed substantial genetic difference from the other common Arabica cultivars, such as Catuai, Borbón, and Mundo Nuevo. This coffee SNP panel provides robust and universally comparable DNA fingerprints, thus can serve as a genotyping tool to assist coffee germplasm management, propagation of planting material, and coffee cultivar authentication.     

A collection of INDEL markers for map-based cloning in seven Arabidopsis accessions

The availability of a comprehensive set of resources including an entire annotated reference genome, sequenced alternative accessions, and a multitude of marker systems makes Arabidopsis thaliana an ideal platform for genetic mapping. PCR markers based on INsertions/DELetions (INDELs) are currently the most frequently used polymorphisms. For the most commonly used mapping combination, Columbia×Landsberg erecta (Col-0×Ler-0), the Cereon polymorphism database is a valuable resource for the generation of polymorphic markers. However, because the number of markers available in public databases for accessions other than Col-0 and Ler-0 is extremely low, mapping using other accessions is far from straightforward. This issue arose while cloning mutations in the Wassilewskija (Ws-4) background. In this work, approaches are described for marker generation in Ws-4 x Col-0. Complementary strategies were employed to generate 229 INDEL markers. Firstly, existing Col-0/Ler-0 Cereon predicted polymorphisms were mined for transferability to Ws-4. Secondly, Ws-0 ecotype Illumina sequence data were analyzed to identify INDELs that could be used for the development of PCR-based markers for Col-0 and Ws-4. Finally, shotgun sequencing allowed the identification of INDELs directly between Col-0 and Ws-4. The polymorphism of the 229 markers was assessed in seven widely used Arabidopsis accessions, and PCR markers that allow a clear distinction between the diverged Ws-0 and Ws-4 accessions are detailed. The utility of the markers was demonstrated by mapping more than 35 mutations in a Col-0×Ws-4 combination, an example of which is presented here. The potential contribution of next generation sequencing technologies to more traditional map-based cloning is discussed.     

High-throughput SNP genotyping with the GoldenGate assay in maize

Single nucleotide polymorphisms (SNPs) are abundant and evenly distributed throughout the genomes of most plant species. They have become an ideal marker system for genetic research in many crops. Several high throughput platforms have been developed that allow rapid and simultaneous genotyping of up to a million SNP markers. In this study, a custom GoldenGate assay containing 1,536 SNPs was developed based on public SNP information for maize and used to genotype two recombinant inbred line (RIL) populations (Zong3 x 87-1, and B73 x By804) and a panel of 154 diverse inbred lines. Over 90% of the SNPs were successfully scored in the diversity panel and the two RIL populations, with a genotyping error rate of less than 2%. A total of 975 SNP markers detected polymorphism in at least one of the two mapping populations, with a polymorphic rate of 38.5% in Zong3 x 87-1 and 52.6% in B73 x By804. The polymorphic SNPs in B73 x By804 have been integrated with previously mapped simple sequence repeat markers to construct a high-density linkage map containing 662 markers with a total length of 1,673.7 cM and an average of 2.53 cM between two markers. The minor allelic frequency (MAF) was distributed evenly across 10 continued classes from 0.05 to 0.5, and about 16% of the SNP markers had a MAF below 10% in the diversity panel. Polymorphism rates for individual SNP markers in pair-wise comparisons of genotypes tested ranged from 0.3 to 63.8% with an average of 36.3%. Most SNPs used in this GoldenGate assay appear to be equally useful for diversity analysis, marker-trait association studies, and marker-aided breeding.     

SNP2CAPS: a SNP and INDEL analysis tool for CAPS marker development

With the influx of various SNP genotyping assays in recent years, there has been a need for an assay that is robust, yet cost effective, and could be performed using standard gel-based procedures. In this context, CAPS markers have been shown to meet these criteria. However, converting SNPs to CAPS markers can be a difficult process if done manually. In order to address this problem, we describe a computer program, SNP2CAPS, that facilitates the computational conversion of SNP markers into CAPS markers. 413 multiple aligned sequences derived from barley ESTs were analysed for the presence of polymorphisms in 235 distinct restriction sites. 282 (90%) of 314 alignments that contain sequence variation due to SNPs and InDels revealed at least one polymorphic restriction site. After reducing the number of restriction enzymes from 235 to 10, 31% of the polymorphic sites could still be detected. In order to demonstrate the usefulness of this tool for marker development, we experimentally validated some of the results predicted by SNP2CAPS.     

Indel arrays: an affordable alternative for genotyping

Natural variation and induced mutations are important resources for gene discovery and the elucidation of genetic circuits. Mapping such polymorphisms requires rapid and cost-efficient methods for genome-wide genotyping. Here we report the development of a microarray-based method that assesses 240 unique markers in a single hybridization experiment at a cost of less than US$50 in materials per line. Our genotyping array is built with 70-mer oligonucleotide elements representing insertion/deletion (indel) polymorphisms between the Arabidopsis thaliana accessions Columbia-0 (Col) and Landsberg erecta (Ler). These indel polymorphisms are recognized with great precision by comparative genomic hybridization, eliminating the need for array replicates and complex statistical analysis. Markers are present genome-wide, with an average spacing of approximately 500 kb. PCR primer information is provided for all array indels, allowing rapid single-locus inquiries. Multi-well chips allow groups of 16 lines to be genotyped in a single experiment. We demonstrate the utility of the array for accurately mapping recessive mutations, RIL populations and mixed genetic backgrounds from accessions other than Col and Ler. Given the ease of use of shotgun sequencing to generate partial genomic sequences of unsequenced species, this approach is readily transferable to non-model organisms.     

High-throughput single nucleotide polymorphism genotyping for breeding applications in rice using the BeadXpress platform

Multiplexed single nucleotide polymorphism (SNP) markers have the potential to increase the speed and cost-effectiveness of genotyping, provided that an optimal SNP density is used for each application. To test the efficiency of multiplexed SNP genotyping for diversity, mapping and breeding applications in rice (Oryza sativa L.), we designed seven GoldenGate VeraCode oligo pool assay (OPA) sets for the Illumina BeadXpress Reader. Validated markers from existing 1536 Illumina SNPs and 44?K Affymetrix SNP chips developed at Cornell University were used to select subsets of informative SNPs for different germplasm groups with even distribution across the genome. A 96-plex OPA was developed for quality control purposes and for assigning a sample into one of the five O. sativa population subgroups. Six 384-plex OPAs were designed for genetic diversity analysis, DNA fingerprinting, and to have evenly-spaced polymorphic markers for quantitative trait locus (QTL) mapping and background selection for crosses between different germplasm pools in rice: Indica/Indica, Indica/Japonica, Japonica/Japonica, Indica/O. rufipogon, and Japonica/O. rufipogon. After testing on a diverse set of rice varieties, two of the SNP sets were re-designed by replacing poor-performing SNPs. Pilot studies were successfully performed for diversity analysis, QTL mapping, marker-assisted backcrossing, and developing specialized genetic stocks, demonstrating that 384-plex SNP genotyping on the BeadXpress platform is a robust and efficient method for marker genotyping in rice.     

Development of user-friendly functional molecular markers for VvDXS gene conferring muscat flavor in grapevine

High fruit and wine quality combined with good climatic adaptation and disease resistance are essential objectives of grape breeding. While several molecular markers are available for pyramiding resistance to fungal pathogens, molecular tools for predicting fruit composition are still scarce. Muscat flavor, caused by the accumulation of monoterpenoids in the berry, is an important target trait for breeding, sought after in both table grapes and wine. Four missense mutations in the VvDXS gene in grape germplasm have been shown to be tightly linked to muscat flavor. Here we present highly reproducible and breeder-friendly functional markers for each of the targeted polymorphisms developed by using either the multiplexed minisequencing SNaPshot? method, the high-resolution melting (HRM) assay or the cleaved amplified polymorphic sequence system. A total of 242 grapevine accessions were analyzed to optimize these different genotyping methods and to provide allele-specific markers for accurate selection of muscat flavor at early stages of grape breeding programs. The HRM and the minisequencing SNaPshot multiplex assays allow for high-throughput automated screening and are suitable for large-scale breeding programs and germplasm characterization.