polymapR-linkage analysis and genetic map construction from F1 populations of outcrossing polyploids.

Motivation: Polyploid species carry more than two copies of each chromosome, a condition found in many of the world's most important crops. Genetic mapping in polyploids is more complex than in diploid species, resulting in a lack of available software tools. These are needed if we are to realize all the opportunities offered by modern genotyping platforms for genetic research and breeding in polyploid crops. Results: polymapR is an R package for genetic linkage analysis and integrated genetic map construction from bi-parental populations of outcrossing autopolyploids. It can currently analyse triploid, tetraploid and hexaploid marker datasets and is applicable to various crops including potato, leek, alfalfa, blueberry, chrysanthemum, sweet potato or kiwifruit. It can detect, estimate and correct for preferential chromosome pairing, and has been tested on high-density marker datasets from potato, rose and chrysanthemum, generating high-density integrated linkage maps in all of these crops. Availability and implementation: polymapR is freely available under the general public license from the Comprehensive R Archive Network (CRAN) at http://cran.r-project.org/package=polymapR. Supplementary information: Supplementary data are available at Bioinformatics online.

Genetic mapping and QTL analysis of Botrytis resistance in Gerbera hybrida.

Gerbera hybrida is an economically important cut flower. In the production and transportation of gerbera with unavoidable periods of high relative humidity, grey mould occurs and results in losses in quality and quantity of flowers. Considering the limitations of chemical use in greenhouses and the impossibility to use these chemicals in auction or after sale, breeding for resistant gerbera cultivars is considered as the best practical approach. In this study, we developed two segregating F1 populations (called S and F). Four parental linkage maps were constructed using common and parental specific SNP markers developed from expressed sequence tag sequencing. Parental genetic maps, containing 30, 29, 27 and 28 linkage groups and a consensus map covering 24 of the 25 expected chromosomes, could be constructed. After evaluation of Botrytis disease severity using three different tests, whole inflorescence, bottom (of disc florets) and ray floret, quantitative trait locus (QTL) mapping was performed using the four individual parental maps. A total of 20 QTLs (including one identical QTL for whole inflorescence and bottom tests) were identified in the parental maps of the two populations. The number of QTLs found and the explained variance of most QTLs detected reflect the complex mechanism of Botrytis disease response.

SNP-markers in Allium species to facilitate introgression breeding in onion.

BACKGROUND: Within onion, Allium cepa L., the availability of disease resistance is limited. The identification of sources of resistance in related species, such as Allium roylei and Allium fistulosum, was a first step towards the improvement of onion cultivars by breeding. SNP markers linked to resistance and polymorphic between these related species and onion cultivars are a valuable tool to efficiently introgress disease resistance genes. In this paper we describe the identification and validation of SNP markers valuable for onion breeding. RESULTS: Transcriptome sequencing resulted in 192 million RNA seq reads from the interspecific F1 hybrid between A. roylei and A. fistulosum (RF) and nine onion cultivars. After assembly, reliable SNPs were discovered in about 36 % of the contigs. For genotyping of the interspecific three-way cross population, derived from a cross between an onion cultivar and the RF (CCxRF), 1100 SNPs that are polymorphic in RF and monomorphic in the onion cultivars (RF SNPs) were selected for the development of KASP assays. A molecular linkage map based on 667 RF-SNP markers was constructed for CCxRF. In addition, KASP assays were developed for 1600 onion-SNPs (SNPs polymorphic among onion cultivars). A second linkage map was constructed for an F2 of onion x A. roylei (F2(CxR)) that consisted of 182 onion-SNPs and 119 RF-SNPs, and 76 previously mapped markers. Markers co-segregating in both the F2(CxR) and the CCxRF population were used to assign the linkage groups of RF to onion chromosomes. To validate usefulness of these SNP markers, QTL mapping was applied in the CCxRF population that segregates for resistance to Botrytis squamosa and resulted in a QTL for resistance on chromosome 6 of A. roylei. CONCLUSIONS: Our research has more than doubled the publicly available marker sequences of expressed onion genes and two onion-related species. It resulted in a detailed genetic map for the interspecific CCxRF population. This is the first paper that reports the detection of a QTL for resistance to B. squamosa in A. roylei.

Genetic mapping of resistance to Fusarium oxysporum f. sp. tulipae in tulip.

Fusarium oxysporum is a major problem in the production of tulip bulbs. Breeding for resistant cultivars through a conventional approach is a slow process due to the long life cycle of tulip. Until now, marker-assisted selection (MAS) has been hampered by the large genome size and the absence of a genetic map. This study is aimed at construction of the first genetic map for tulip and at the identification of loci associated with resistance to F. oxysporum. A cross-pollinated population of 125 individuals segregating for Fusarium resistance was obtained from Tulipa gesneriana "Kees Nelis" and T. fosteriana "Cantata." Fusarium resistance of the mapping population was evaluated through a soil infection test in two consecutive years, and a spot inoculation test in which a green fluorescent protein tagged Fusarium strain was used for inoculation. The genetic maps have been constructed for the parents separately. The genetic map of "Kees Nelis" comprised 342 markers on 27 linkage groups covering 1707 cM, while the map of "Cantata" comprised 300 markers on 21 linkage groups covering 1201 cM. Median distance between markers was 3.9 cM for "Kees Nelis" and 3.1 cM for "Cantata." Six putative quantitative trait loci (QTLs) for Fusarium resistance were identified, derived from both parents. QTL2, QTL3, and QTL6 were significant in all disease tests. For the flanking markers of the QTLs, phenotypic means of the two allelic groups, segregating from a parent for such a marker, were significantly different. These markers will be useful for the development of MAS in tulip breeding.

Using multi-locus allelic sequence data to estimate genetic divergence among four Lilium (Liliaceae) cultivars.

Next Generation Sequencing (NGS) may enable estimating relationships among genotypes using allelic variation of multiple nuclear genes simultaneously. We explored the potential and caveats of this strategy in four genetically distant Lilium cultivars to estimate their genetic divergence from transcriptome sequences using three approaches: POFAD (Phylogeny of Organisms from Allelic Data, uses allelic information of sequence data), RAxML (Randomized Accelerated Maximum Likelihood, tree building based on concatenated consensus sequences) and Consensus Network (constructing a network summarizing among gene tree conflicts). Twenty six gene contigs were chosen based on the presence of orthologous sequences in all cultivars, seven of which also had an orthologous sequence in Tulipa, used as out-group. The three approaches generated the same topology. Although the resolution offered by these approaches is high, in this case there was no extra benefit in using allelic information. We conclude that these 26 genes can be widely applied to construct a species tree for the genus Lilium.

Generation and analysis of expressed sequence tags in the extreme large genomes Lilium and Tulipa.

BACKGROUND: Bulbous flowers such as lily and tulip (Liliaceae family) are monocot perennial herbs that are economically very important ornamental plants worldwide. However, there are hardly any genetic studies performed and genomic resources are lacking. To build genomic resources and develop tools to speed up the breeding in both crops, next generation sequencing was implemented. We sequenced and assembled transcriptomes of four lily and five tulip genotypes using 454 pyro-sequencing technology. RESULTS: Successfully, we developed the first set of 81,791 contigs with an average length of 514 bp for tulip, and enriched the very limited number of 3,329 available ESTs (Expressed Sequence Tags) for lily with 52,172 contigs with an average length of 555 bp. The contigs together with singletons covered on average 37% of lily and 39% of tulip estimated transcriptome. Mining lily and tulip sequence data for SSRs (Simple Sequence Repeats) showed that di-nucleotide repeats were twice more abundant in UTRs (UnTranslated Regions) compared to coding regions, while tri-nucleotide repeats were equally spread over coding and UTR regions. Two sets of single nucleotide polymorphism (SNP) markers suitable for high throughput genotyping were developed. In the first set, no SNPs flanking the target SNP (50 bp on either side) were allowed. In the second set, one SNP in the flanking regions was allowed, which resulted in a 2 to 3 fold increase in SNP marker numbers compared with the first set. Orthologous groups between the two flower bulbs: lily and tulip (12,017 groups) and among the three monocot species: lily, tulip, and rice (6,900 groups) were determined using OrthoMCL. Orthologous groups were screened for common SNP markers and EST-SSRs to study synteny between lily and tulip, which resulted in 113 common SNP markers and 292 common EST-SSR. Lily and tulip contigs generated were annotated and described according to Gene Ontology terminology. CONCLUSIONS: Two transcriptome sets were built that are valuable resources for marker development, comparative genomic studies and candidate gene approaches. Next generation sequencing of leaf transcriptome is very effective; however, deeper sequencing and using more tissues and stages is advisable for extended comparative studies.

SNP markers retrieval for a non-model species: a practical approach.

BACKGROUND: SNP (Single Nucleotide Polymorphism) markers are rapidly becoming the markers of choice for applications in breeding because of next generation sequencing technology developments. For SNP development by NGS technologies, correct assembly of the huge amounts of sequence data generated is essential. Little is known about assembler's performance, especially when dealing with highly heterogeneous species that show a high genome complexity and what the possible consequences are of differences in assemblies on SNP retrieval. This study tested two assemblers (CAP3 and CLC) on 454 data from four lily genotypes and compared results with respect to SNP retrieval. RESULTS: CAP3 assembly resulted in higher numbers of contigs, lower numbers of reads per contig, and shorter average read lengths compared to CLC. Blast comparisons showed that CAP3 contigs were highly redundant. Contrastingly, CLC in rare cases combined paralogs in one contig. Redundant and chimeric contigs may lead to erroneous SNPs. Filtering for redundancy can be done by blasting selected SNP markers to the contigs and discarding all the SNP markers that show more than one blast hit. Results on chimeric contigs showed that only four out of 2,421 SNP markers were selected from chimeric contigs. CONCLUSION: In practice, CLC performs better in assembling highly heterogeneous genome sequences compared to CAP3, and consequently SNP retrieval is more efficient. Additionally a simple flow scheme is suggested for SNP marker retrieval that can be valid for all non-model species.