Potential of genotyping-by-sequencing for genomic selection in livestock populations

Background

Next-generation sequencing techniques, such as genotyping-by-sequencing (GBS), provide alternatives to single nucleotide polymorphism (SNP) arrays. The aim of this work was to evaluate the potential of GBS compared to SNP array genotyping for genomic selection in livestock populations.

Methods

The value of GBS was quantified by simulation analyses in which three parameters were varied: (i) genome-wide sequence read depth (x) per individual from 0.01x to 20x or using SNP array genotyping; (ii) number of genotyped markers from 3000 to 300 000; and (iii) size of training and prediction sets from 500 to 50 000 individuals. The latter was achieved by distributing the total available x of 1000x, 5000x, or 10 000x per genotyped locus among the varying number of individuals. With SNP arrays, genotypes were called from sequence data directly. With GBS, genotypes were called from sequence reads that varied between loci and individuals according to a Poisson distribution with mean equal to x. Simulated data were analyzed with ridge regression and the accuracy and bias of genomic predictions and response to selection were quantified under the different scenarios.

Results

Accuracies of genomic predictions using GBS data or SNP array data were comparable when large numbers of markers were used and x per individual was ~1x or higher. The bias of genomic predictions was very high at a very low x. When the total available x was distributed among the training individuals, the accuracy of prediction was maximized when a large number of individuals was used that had GBS data with low x for a large number of markers. Similarly, response to selection was maximized under the same conditions due to increasing both accuracy and selection intensity.

Conclusions

GBS offers great potential for developing genomic selection in livestock populations because it makes it possible to cover large fractions of the genome and to vary the sequence read depth per individual. Thus, the accuracy of predictions is improved by increasing the size of training populations and the intensity of selection is increased by genotyping a larger number of selection candidates.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-015-0102-z) contains supplementary material, which is available to authorized users.     

Prospects and limits of marker imputation in quantitative genetic studies in European elite wheat (Triticum aestivum L.)

Background

The main goal of our study was to investigate the implementation, prospects, and limits of marker imputation for quantitative genetic studies contrasting map-independent and map-dependent algorithms. We used a diversity panel consisting of 372 European elite wheat (Triticum aestivum L.) varieties, which had been genotyped with SNP arrays, and performed intensive simulation studies.

Results

Our results clearly showed that imputation accuracy was substantially higher for map-dependent compared to map-independent methods. The accuracy of marker imputation depended strongly on the linkage disequilibrium between the markers in the reference panel and the markers to be imputed. For the decay of linkage disequilibrium present in European wheat, we concluded that around 45,000 markers are needed for low cost, low-density marker profiling. This will facilitate high imputation accuracy, also for rare alleles. Genomic selection and diversity studies profited only marginally from imputing missing values. In contrast, the power of association mapping increased substantially when missing values were imputed.

Conclusions

Imputing missing values is especially of interest for an economic implementation of association mapping in breeding populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1366-y) contains supplementary material, which is available to authorized users.     

MPAgenomics: an R package for multi-patient analysis of genomic markers

Background

Last generations of Single Nucleotide Polymorphism (SNP) arrays allow to study copy-number variations in addition to genotyping measures.

Results

MPAgenomics, standing for multi-patient analysis (MPA) of genomic markers, is an R-package devoted to: (i) efficient segmentation and (ii) selection of genomic markers from multi-patient copy number and SNP data profiles. It provides wrappers from commonly used packages to streamline their repeated (sometimes difficult) manipulation, offering an easy-to-use pipeline for beginners in R.The segmentation of successive multiple profiles (finding losses and gains) is performed with an automatic choice of parameters involved in the wrapped packages. Considering multiple profiles in the same time, MPAgenomics wraps efficient penalized regression methods to select relevant markers associated with a given outcome.

Conclusions

MPAgenomics provides an easy tool to analyze data from SNP arrays in R. The R-package MPAgenomics is available on CRAN.     

A bi-filtering method for processing single nucleotide polymorphism array data improves the quality of genetic map and accuracy of quantitative trait locus mapping in doubled haploid populations of polyploid Brassica napus

Background

Single nucleotide polymorphism (SNP) markers have a wide range of applications in crop genetics and genomics. Due to their polyploidy nature, many important crops, such as wheat, cotton and rapeseed contain a large amount of repeat and homoeologous sequences in their genomes, which imposes a huge challenge in high-throughput genotyping with sequencing and/or array technologies. Allotetraploid Brassica napus (AACC, 2n = 4x = 38) comprises of two highly homoeologous sub-genomes derived from its progenitor species B. rapa (AA, 2n = 2x = 20) and B. oleracea (CC, 2n = 2x = 18), and is an ideal species to exploit methods for reducing the interference of extensive inter-homoeologue polymorphisms (mHemi-SNPs and Pseudo-simple SNPs) between closely related sub-genomes.

Results

Based on a recent B. napus 6K SNP array, we developed a bi-filtering procedure to identify unauthentic lines in a DH population, and mHemi-SNPs and Pseudo-simple SNPs in an array data matrix. The procedure utilized both monomorphic and polymorphic SNPs in the DH population and could effectively distinguish the mHemi-SNPs and Pseudo-simple SNPs that resulted from superposition of the signals from multiple SNPs. Compared with conventional procedure for array data processing, the bi-filtering method could minimize the pseudo linkage relationship caused by the mHemi-SNPs and Pseudo-simple SNPs, thus improving the quality of SNP genetic map. Furthermore, the improved genetic map could increase the accuracies of mapping of QTLs as demonstrated by the ability to eliminate non-real QTLs in the mapping population.

Conclusions

The bi-filtering analysis of the SNP array data represents a novel approach to effectively assigning the multi-loci SNP genotypes in polyploid B. napus and may find wide applications to SNP analyses in polyploid crops.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1559-4) contains supplementary material, which is available to authorized users.     

Assessment of alternative genotyping strategies to maximize imputation accuracy at minimal cost

Background

Commercial breeding programs seek to maximise the rate of genetic gain while minimizing the costs of attaining that gain. Genomic information offers great potential to increase rates of genetic gain but it is expensive to generate. Low-cost genotyping strategies combined with genotype imputation offer dramatically reduced costs. However, both the costs and accuracy of imputation of these strategies are highly sensitive to several factors. The objective of this paper was to explore the cost and imputation accuracy of several alternative genotyping strategies in pedigreed populations.

Methods

Pedigree and genotype data from a commercial pig population were used. Several alternative genotyping strategies were explored. The strategies differed in the density of genotypes used for the ancestors and the individuals to be imputed. Parents, grandparents, and other relatives that were not descendants, were genotyped at high-density, low-density, or extremely low-density, and associated costs and imputation accuracies were evaluated.

Results

Imputation accuracy and cost were influenced by the alternative genotyping strategies. Given the mating ratios and the numbers of offspring produced by males and females, an optimized low-cost genotyping strategy for a commercial pig population could involve genotyping male parents at high-density, female parents at low-density (e.g. 3000 SNP), and selection candidates at very low-density (384 SNP).

Conclusions

Among the selection candidates, 95.5 % and 93.5 % of the genotype variation contained in the high-density SNP panels were recovered using a genotyping strategy that costs respectively, $24.74 and $20.58 per candidate.     

Application of whole genome re-sequencing data in the development of diagnostic DNA markers tightly linked to a disease-resistance locus for marker-assisted selection in lupin (Lupinus angustifolius)

Background

Molecular marker-assisted breeding provides an efficient tool to develop improved crop varieties. A major challenge for the broad application of markers in marker-assisted selection is that the marker phenotypes must match plant phenotypes in a wide range of breeding germplasm. In this study, we used the legume crop species Lupinus angustifolius (lupin) to demonstrate the utility of whole genome sequencing and re-sequencing on the development of diagnostic markers for molecular plant breeding.

Results

Nine lupin cultivars released in Australia from 1973 to 2007 were subjected to whole genome re-sequencing. The re-sequencing data together with the reference genome sequence data were used in marker development, which revealed 180,596 to 795,735 SNP markers from pairwise comparisons among the cultivars. A total of 207,887 markers were anchored on the lupin genetic linkage map. Marker mining obtained an average of 387 SNP markers and 87 InDel markers for each of the 24 genome sequence assembly scaffolds bearing markers linked to 11 genes of agronomic interest. Using the R gene PhtjR conferring resistance to phomopsis stem blight disease as a test case, we discovered 17 candidate diagnostic markers by genotyping and selecting markers on a genetic linkage map. A further 243 candidate diagnostic markers were discovered by marker mining on a scaffold bearing non-diagnostic markers linked to the PhtjR gene. Nine out from the ten tested candidate diagnostic markers were confirmed as truly diagnostic on a broad range of commercial cultivars. Markers developed using these strategies meet the requirements for broad application in molecular plant breeding.

Conclusions

We demonstrated that low-cost genome sequencing and re-sequencing data were sufficient and very effective in the development of diagnostic markers for marker-assisted selection. The strategies used in this study may be applied to any trait or plant species. Whole genome sequencing and re-sequencing provides a powerful tool to overcome current limitations in molecular plant breeding, which will enable plant breeders to precisely pyramid favourable genes to develop super crop varieties to meet future food demands.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1878-5) contains supplementary material, which is available to authorized users.     

Association analysis for udder health based on SNP-panel and sequence data in Danish Holsteins

Background

The sensitivity of genome-wide association studies for the detection of quantitative trait loci (QTL) depends on the density of markers examined and the statistical models used. This study compares the performance of three marker densities to refine six previously detected QTL regions for mastitis traits: 54 k markers of a medium-density SNP (single nucleotide polymorphism) chip (MD), imputed 777 k markers of a high-density SNP chip (HD), and imputed whole-genome sequencing data (SEQ). Each dataset contained data for 4496 Danish Holstein cattle. Comparisons were performed using a linear mixed model (LM) and a Bayesian variable selection model (BVS).

Results

After quality control, 587, 7825, and 78 856 SNPs in the six targeted regions remained for MD, HD, and SEQ data, respectively. In general, the association patterns between SNPs and traits were similar for the three marker densities when tested using the same statistical model. With the LM model, 120 (MD), 967 (HD), and 7209 (SEQ) SNPs were significantly associated with mastitis, whereas with the BVS model, 43 (MD), 131 (HD), and 1052 (SEQ) significant SNPs (Bayes factor > 3.2) were observed. A total of 26 (MD), 75 (HD), and 465 (SEQ) significant SNPs were identified by both models. In addition, one, 16, and 33 QTL peaks for MD, HD, and SEQ data were detected according to the QTL intensity profile of SNP bins by post-analysis of the BVS model.

Conclusions

The power to detect significant associations increased with increasing marker density. The BVS model resulted in clearer boundaries between linked QTL than the LM model. Using SEQ data, the six targeted regions were refined to 33 candidate QTL regions for udder health. The comparison between these candidate QTL regions and known genes suggested that NPFFR2, SLC4A4, DCK, LIFR, and EDN3 may be considered as candidate genes for mastitis susceptibility.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-015-0129-1) contains supplementary material, which is available to authorized users.     

Application of in silico bulked segregant analysis for rapid development of markers linked to Bean common mosaic virus resistance in common bean

Background

Common bean was one of the first crops that benefited from the development and utilization of molecular marker-assisted selection (MAS) for major disease resistance genes. Efficiency of MAS for breeding common bean is still hampered, however, due to the dominance, linkage phase, and loose linkage of previously developed markers. Here we applied in silico bulked segregant analysis (BSA) to the BeanCAP diversity panel, composed of over 500 lines and genotyped with the BARCBEAN_3 6K SNP BeadChip, to develop codominant and tightly linked markers to the I gene controlling resistance to Bean common mosaic virus (BCMV).

Results

We physically mapped the genomic region underlying the I gene. This locus, in the distal arm of chromosome Pv02, contains seven putative NBS-LRR-type disease resistance genes. Two contrasting bulks, containing BCMV host differentials and ten BeanCAP lines with known disease reaction to BCMV, were subjected to in silico BSA for targeting the I gene and flanking sequences. Two distinct haplotypes, containing a cluster of six single nucleotide polymorphisms (SNP), were associated with resistance or susceptibility to BCMV. One-hundred and twenty-two lines, including 115 of the BeanCAP panel, were screened for BCMV resistance in the greenhouse, and all of the resistant or susceptible plants displayed distinct SNP haplotypes as those found in the two bulks. The resistant/susceptible haplotypes were validated in 98 recombinant inbred lines segregating for BCMV resistance. The closest SNP (~25-32 kb) to the distal NBS-LRR gene model for the I gene locus was targeted for conversion to codominant KASP (Kompetitive Allele Specific PCR) and CAPS (Cleaved Amplified Polymorphic Sequence) markers. Both marker systems accurately predicted the disease reaction to BCMV conferred by the I gene in all screened lines of this study.

Conclusions

We demonstrated the utility of the in silico BSA approach using genetically diverse germplasm, genotyped with a high-density SNP chip array, to discover SNP variation at a specific targeted genomic region. In common bean, many disease resistance genes are mapped and their physical genomic position can now be determined, thus the application of this approach will facilitate further development of codominant and tightly linked markers for use in MAS.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-903) contains supplementary material, which is available to authorized users.     

Genome-wide association analysis of seedling root development in maize (Zea mays L.)

Background

Plants rely on the root system for anchorage to the ground and the acquisition and absorption of nutrients critical to sustaining productivity. A genome wide association analysis enables one to analyze allelic diversity of complex traits and identify superior alleles. 384 inbred lines from the Ames panel were genotyped with 681,257 single nucleotide polymorphism markers using Genotyping-by-Sequencing technology and 22 seedling root architecture traits were phenotyped.

Results

Utilizing both a general linear model and mixed linear model, a GWAS study was conducted identifying 268 marker trait associations (p ≤ 5.3×10^-7). Analysis of significant SNP markers for multiple traits showed that several were located within gene models with some SNP markers localized within regions of previously identified root quantitative trait loci. Gene model GRMZM2G153722 located on chromosome 4 contained nine significant markers. This predicted gene is expressed in roots and shoots.

Conclusion

This study identifies putatively associated SNP markers associated with root traits at the seedling stage. Some SNPs were located within or near (<1 kb) gene models. These gene models identify possible candidate genes involved in root development at the seedling stage. These and respective linked or functional markers could be targets for breeders for marker assisted selection of seedling root traits.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1226-9) contains supplementary material, which is available to authorized users.     

Fine mapping of a QTL for ear size on porcine chromosome 5 and identification of high mobility group AT-hook 2 (HMGA2) as a positional candidate gene

Background

Ear size and shape are distinct conformation characteristics of pig breeds. Previously, we identified a significant quantitative trait locus (QTL) influencing ear surface on pig chromosome 5 in a White Duroc × Erhualian F₂ resource population. This QTL explained more than 17% of the phenotypic variance.

Methods

Four new markers on pig chromosome 5 were genotyped across this F₂ population. RT-PCR was performed to obtain expression profiles of different candidate genes in ear tissue. Standard association test, marker-assisted association test and F-drop test were applied to determine the effects of single nucleotide polymorphisms (SNP) on ear size. Three synthetic commercial lines were also used for the association test.

Results

We refined the QTL to an 8.7-cM interval and identified three positional candidate genes i.e. HMGA2, SOX5 and PTHLH that are expressed in ear tissue. Seven SNP within these three candidate genes were selected and genotyped in the F₂ population. Of the seven SNP, HMGA2 SNP ({"type":"entrez-nucleotide","attrs":{"text":"JF748727","term_id":"347311053","term_text":"JF748727"}}JF748727: g.2836 A > G) showed the strongest association with ear size in the standard association test and marker-assisted association test. With the F-drop test, F value decreased by more than 97% only when the genotypes of HMGA2 g.2836 A > G were included as a fixed effect. Furthermore, the significant association between g.2836 A > G and ear size was also demonstrated in the synthetic commercial Sutai pig line. The haplotype-based association test showed that the phenotypic variance explained by HMGA2 was similar to that explained by the QTL and at a much higher level than by SOX5. More interestingly, HMGA2 is also located within the dog orthologous chromosome region, which has been shown to be associated with ear type and size.

Conclusions

HMGA2 was the closest gene with a potential functional effect to the QTL or marker for ear size on chromosome 5. This study will contribute to identify the causative gene and mutation underlying this QTL.     

Polymorphisms of the Tissue Inhibitor of Metalloproteinase 3 Gene Are Associated with Resistance to High-Altitude Pulmonary Edema (HAPE) in a Japanese Population: A Case Control Study Using Polymorphic Microsatellite Markers

Introduction

High-altitude pulmonary edema (HAPE) is a hypoxia-induced, life-threatening, high permeability type of edema attributable to pulmonary capillary stress failure. Genome-wide association analysis is necessary to better understand how genetics influence the outcome of HAPE.

Materials and Methods

DNA samples were collected from 53 subjects susceptible to HAPE (HAPE-s) and 67 elite Alpinists resistant to HAPE (HAPE-r). The genome scan was carried out using 400 polymorphic microsatellite markers throughout the whole genome in all subjects. In addition, six single nucleotide polymorphisms (SNPs) of the gene encoding the tissue inhibitor of metalloproteinase 3 (TIMP3) were genotyped by Taqman® SNP Genotyping Assays.

Results

The results were analyzed using case-control comparisons. Whole genome scanning revealed that allele frequencies in nine markers were statistically different between HAPE-s and HAPE-r subjects. The SNP genotyping of the TIMP3 gene revealed that the derived allele C of rs130293 was associated with resistance to HAPE [odds ratio (OR) = 0.21, P = 0.0012) and recessive inheritance of the phenotype of HAPE-s (P = 0.0012). A haplotype CAC carrying allele C of rs130293 was associated with resistance to HAPE.

Discussion

This genome-wide association study revealed several novel candidate genes associated with susceptibility or resistance to HAPE in a Japanese population. Among those, the minor allele C of rs130293 (C/T) in the TIMP3 gene was linked to resistance to HAPE; while, the ancestral allele T was associated with susceptibility to HAPE.     

The SNP rs3128965 of HLA-DPB1 as a Genetic Marker of the AERD Phenotype

Background

Two common clinical syndromes of acetylsalicylic acid (aspirin) hypersensitivity, aspirin-exacerbated respiratory disease (AERD) and aspirin-exacerbated cutaneous disease (AECD), were subjected to a genome-wide association study to identify strong genetic markers for aspirin hypersensitivity in a Korean population.

Methods

A comparison of SNP genotype frequencies on an Affymetrix Genome-Wide Human SNP array of 179 AERD patients and 1989 healthy normal control subjects (NC) revealed SNPs on chromosome 6 that were associated with AERD, but not AECD. To validate the association, we enrolled a second cohort comprising AERD (n = 264), NC (n = 238) and disease-control (aspirin tolerant asthma; ATA, n = 387) groups.

Results

The minor genotype frequency (AG or AA) of a particular SNP, rs3128965, in the HLA-DPB1 region was higher in the AERD group compared to the ATA or NC group (P = 0.001, P = 0.002, in a co-dominant analysis model, respectively). Comparison of rs3128965 alleles with the clinical features of asthmatics revealed that patients harboring the A allele had increased bronchial hyperresponsiveness to inhaled aspirin and methacholine, and higher 15-HETE levels, than those without the A allele (P = 0.039, 0.037, and 0.004, respectively).

Conclusions

This implies the potential of rs3128965 as a genetic marker for diagnosis and prediction of the AERD phenotype.     

Development of EST-based SNP and InDel markers and their utilization in tetraploid cotton genetic mapping

Background

Availability of molecular markers has proven to be an efficient tool in facilitating progress in plant breeding, which is particularly important in the case of less researched crops such as cotton. Considering the obvious advantages of single nucleotide polymorphisms (SNPs) and insertion-deletion polymorphisms (InDels), expressed sequence tags (ESTs) were analyzed in silico to identify SNPs and InDels in this study, aiming to develop more molecular markers in cotton.

Results

A total of 1,349 EST-based SNP and InDel markers were developed by comparing ESTs between Gossypium hirsutum and G. barbadense, mining G. hirsutum unigenes, and analyzing 3′ untranslated region (3′UTR) sequences. The marker polymorphisms were investigated using the two parents of the mapping population based on the single-strand conformation polymorphism (SSCP) analysis. Of all the markers, 137 (10.16%) were polymorphic, and revealed 142 loci. Linkage analysis using a BC₁ population mapped 133 loci on the 26 chromosomes. Statistical analysis of base variations in SNPs showed that base transitions accounted for 55.78% of the total base variations and gene ontology indicated that cotton genes varied greatly in harboring SNPs ranging from 1.00 to 24.00 SNPs per gene. Sanger sequencing of three randomly selected SNP markers revealed discrepancy between the in silico predicted sequences and the actual sequencing results.

Conclusions

In silico analysis is a double-edged blade to develop EST-SNP/InDel markers. On the one hand, the designed markers can be well used in tetraploid cotton genetic mapping. And it plays a certain role in revealing transition preference and SNP frequency of cotton genes. On the other hand, the developmental efficiency of markers and polymorphism of designed primers are comparatively low.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1046) contains supplementary material, which is available to authorized users.     

Influence of salinity on bacterioplankton communities from the Brazilian rain forest to the coastal Atlantic Ocean

Background

Planktonic bacteria are recognized as important drivers of biogeochemical processes in all aquatic ecosystems, however, the taxa that make up these communities are poorly known. The aim of this study was to investigate bacterial communities in aquatic ecosystems at Ilha Grande, Rio de Janeiro, Brazil, a preserved insular environment of the Atlantic rain forest and how they correlate with a salinity gradient going from terrestrial aquatic habitats to the coastal Atlantic Ocean.

Methodology/Principal Findings

We analyzed chemical and microbiological parameters of water samples and constructed 16S rRNA gene libraries of free living bacteria obtained at three marine (two coastal and one offshore) and three freshwater (water spring, river, and mangrove) environments. A total of 836 sequences were analyzed by MOTHUR, yielding 269 freshwater and 219 marine operational taxonomic units (OTUs) grouped at 97% stringency. Richness and diversity indexes indicated that freshwater environments were the most diverse, especially the water spring. The main bacterial group in freshwater environments was Betaproteobacteria (43.5%), whereas Cyanobacteria (30.5%), Alphaproteobacteria (25.5%), and Gammaproteobacteria (26.3%) dominated the marine ones. Venn diagram showed no overlap between marine and freshwater OTUs at 97% stringency. LIBSHUFF statistics and PCA analysis revealed marked differences between the freshwater and marine libraries suggesting the importance of salinity as a driver of community composition in this habitat. The phylogenetic analysis of marine and freshwater libraries showed that the differences in community composition are consistent.

Conclusions/Significance

Our data supports the notion that a divergent evolutionary scenario is driving community composition in the studied habitats. This work also improves the comprehension of microbial community dynamics in tropical waters and how they are structured in relation to physicochemical parameters. Furthermore, this paper reveals for the first time the pristine bacterioplankton communities in a tropical island at the South Atlantic Ocean.     

Association of INSIG2 Polymorphism with Overweight and LDL in Children

Background

Dyslipidemia and overweight are common issues in children. Identifying genetic markers of risk could lead to targeted interventions. A polymorphism of SNP rs7566605 near insulin-induced gene 2 (INSIG2) has been identified as a strong candidate gene for obesity, through its feedback control of lipid synthesis.

Objective

To identify polymorphisms in INSIG2 which are associated with overweight (BMI ≥ 85% for age) and dyslipidemia in children. Hypothesis: The C allele of rs7566605 would be significantly associated with BMI and LDL.

Design/Methods

We genotyped 15 SNPs in/near INSIG2 in 1,058 healthy children (53% non-Hispanic white (NHW), 37% overweight) participating in a school based study. Genotype was compared with BMI and lipid markers, adjusting for age, gender, and puberty.

Results

We found a significant association between the SNP rs12464355 and LDL in NHW children, p < 0.001. The G allele is protective (lower LDL). A different SNP was associated with overweight in NHW: rs17047757. SNP rs7566605 was not associated with overweight or lipid levels.

Conclusions

We identified novel genetic associations between INSIG2 and both overweight and LDL in NHW children. Polymorphisms in INSIG2 may be important in the development of obesity through its effects on lipid regulation.     

Mapping in the era of sequencing: high density genotyping and its application for mapping TYLCV resistance in Solanum pimpinellifolium

Background

A RIL population between Solanum lycopersicum cv. Moneymaker and S. pimpinellifolium G1.1554 was genotyped with a custom made SNP array. Additionally, a subset of the lines was genotyped by sequencing (GBS).

Results

A total of 1974 polymorphic SNPs were selected to develop a linkage map of 715 unique genetic loci. We generated plots for visualizing the recombination patterns of the population relating physical and genetic positions along the genome.This linkage map was used to identify two QTLs for TYLCV resistance which contained favourable alleles derived from S. pimpinellifolium. Further GBS was used to saturate regions of interest, and the mapping resolution of the two QTLs was improved. The analysis showed highest significance on Chromosome 11 close to the region of 51.3 Mb (qTy-p11) and another on Chromosome 3 near 46.5 Mb (qTy-p3). Furthermore, we explored the population using untargeted metabolic profiling, and the most significant differences between susceptible and resistant plants were mainly associated with sucrose and flavonoid glycosides.

Conclusions

The SNP information obtained from an array allowed a first QTL screening of our RIL population. With additional SNP data of a RILs subset, obtained through GBS, we were able to perform an in silico mapping improvement to further confirm regions associated with our trait of interest. With the combination of different ~ omics platforms we provide valuable insight into the genetics of S. pimpinellifolium-derived TYLCV resistance.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1152) contains supplementary material, which is available to authorized users.