The use of a candidate gene approach to study Botrytis cinerea resistance in Gerbera hybrida.

Candidate genes (CG) for Botrytis cinerea resistance described in literature were mapped on gerbera linkage maps for which several QTL for Botrytis resistance had been found previously using a rapid, low-cost platform for SNP genotyping. In total, 29 CGs were mapped in either of two mapping populations. Four CGs were mapped within the previous identified QTL intervals and three co-localized with QTL. Two of these CGs for resistance against B. cinerea, PG1 (polygalacturonase gene) and sit (sitiens, ABA-aldehyde oxidase gene) that mapped in QTL regions for the ray floret disease resistance test were studied in detail. Virus-induced gene silencing (VIGS) was used for gene function analysis to determine the CGs' role in gerbera resistance to Botrytis. Ray florets, of which the CGs were silenced, showed a significantly delayed growth of lesions upon Botrytis infection compared to controls. Combining QTL analysis, candidate gene mapping and VIGS showed to be an useful combination to identify possible causal genes and for understanding the molecular mechanisms of Botrytis resistance in gerbera. The two genes seem to act as partial S-genes and are likely among the determining genes leading to the variation observed for B. cinerea resistance in gerbera.

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.

Biotechnological advances in Lilium.

Modern powerful techniques in plant biotechnology have been developed in lilies (Lilium spp., Liliaceae) to propagate, improve and make new phenotypes. Reliable in vitro culture methods are available to multiply lilies rapidly and shorten breeding programs. Lilium is also an ideal model plant to study in vitro pollination and embryo rescue methods. Although lilies are recalcitrant to genetic manipulation, superior genotypes are developed with improved flower colour and form, disease resistance and year round forcing ability. Different DNA molecular markers have been developed for rapid indirect selection, genetic diversity evaluation, mutation detection and construction of Lilium linkage map. Some disease resistance-QTLs are already mapped on the Lilium linkage map. This review presents latest information on in vitro propagation, genetic engineering and molecular advances made in lily.

Transcriptome Analysis of Gerbera hybrida Including in silico Confirmation of Defense Genes Found.

For the ornamental crop Gerbera hybrida, breeding at the moment is done using conventional methods. As this has drawbacks in breeding speed and efficiency, especially for complex traits like disease resistance, we set out to develop genomic resources. The leaf and flower bud transcriptomes of four parents, used to generate two gerbera populations, were sequenced using Illumina paired-end sequencing. In total, 36,770 contigs with an average length of 1397 bp were generated and these have been the starting point for SNP identification and annotation. The consensus contig sequences were used to map reads of individual parents, to identify genotype specific SNPs, and to assess the presence of common SNPs between genotypes. Comparison with the non-redundant protein database (nr) showed that 29,146 contigs gave BLAST hits. Of sequences with blast results, 73.3% obtained a clear gene ontology (GO) annotation. EST contigs coding for enzymes were found in Kyoto Encyclopedia of Genes and Genomes maps (KEGG). Through, these annotated data and KEGG molecular interaction network, transcripts associated with the phenylpropanoid metabolism, other secondary metabolite biosynthesis pathways, phytohormone biosynthesis and signal transduction were analyzed in more detail. Identifying genes involved in these processes could provide genetic and genomic resources for studying the mechanism of disease resistance in gerbera.

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.

Production of polyploids and unreduced gametes in Lilium auratum × L. henryi hybrid.

Intergenomic F1 hybrids between L. auratum x L. henryi and their BC1 progeny were investigated through genomic in situ hybridization technique (GISH) to determine their potential value in lily breeding. We confirmed that F1 intergenomic hybrids possessed a set of chromosomes (x=12) from both parents and that flowers of the F1 auratum × henryi hybrid showed an intermediate morphological phenotype. Pollen size, viability and germination ability were measured through microscopic observations. F1 intergenomic hybrids produced a relevant frequency of 2n-gametes, which were successfully used to perform crosses with Oriental hybrids, resulting in the triploid Oriental Auratum Henryi (OAuH) hybrid. Twenty BC1 plants were generated by crossing between four different Oriental hybrid cultivars and F1 AuH hybrids using an in vitro embryo rescue technique, after which the genome constitution and chromosome composition were analyzed by GISH. All plants were triploid, showing 12 from female parents (diploid Oriental hybrid) and 24 from male parents (diploid F1 AuH hybrid). Overall, 16 out of 20 BC1 progeny possessed recombinant chromosomes with 1-5 crossover sites per plant. Cytological analysis of 20 BC1 plants by GISH verified that the occurrence of 2n pollen formation in all F1 AuH hybrids was derived from the FDR (first division restitution) mechanism, in which the genome composition of all BC1 plants possess 12 Oriental + 12 L. auratum + 12 L. henryi chromosomes. Allotriploids derived from the AuH hybrid were used as female for crossing with the diploid Oriental hybrid cultivar 'Sorbonne' and considerable numbers of plants (0-6.5 plants per ovary) were only obtained when female OAuH (BC1) triploids were used. Taken together, the results of this study indicate that production and analysis of F1 AuH hybrids and their progeny through sexual polyploidization can be useful for efficient creation of important horticultural traits.

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.

An assessment of chromosomal rearrangements in neopolyploids of Lilium hybrids.

Two types of newly induced polyploids (neopolyploids) of Lilium hybrids were monitored for the occurrence of chromosomal rearrangements through genomic in situ hybridization (GISH) technique. One of the populations was obtained through crossing an allotriploid Longiflorum x Oriental hybrid (LLO) with an allotetraploid Longiflorum x Trumpet hybrid (LLTT), both of which were derived from somatic chromosome doubling. The other type of allopolyploid population was derived from meiotic chromosome doubling in which numerically unreduced (2n) gametes from two different interspecific hybrids, namely, Longiflorum x Asiatic (LA) and Oriental x Asiatic (OA), were used to get backcross progeny with the Asiatic parents. GISH clearly discriminated the three constituent genomes (L, T, and O) in the complements of the progeny obtained from mitotic chromosome doubling. A total of 26 individuals were analyzed from this population and there was no evidence of chromosomal rearrangements. However, in the case of meiotically doubled allopolyploid progeny, considerable frequencies of chromosomal rearrangements were observed through GISH. The so-called chromosomal rearrangements in meiotic polyploids are the result of homoeologous recombination rather than translocations. Furthermore, evidence for the occurrence of meiotic recombination in the LA hybrids has been confirmed with GISH on meiotic chromosomes. Thus, there was evidence that neopolyploids of Lilium hybrids did not possess any noticeable chromosome rearrangements.

Construction of chromosomal recombination maps of three genomes of lilies (Lilium) based on GISH analysis.

Chromosomal recombination maps were constructed for three genomes of lily (Lilium) using GISH analyses. For this purpose, the backcross (BC) progenies of two diploid (2n = 2x = 24) interspecific hybrids of lily, viz. Longiflorum x Asiatic (LA) and Oriental x Asiatic (OA), were used. Mostly the BC progenies of LA hybrids consisted of both triploid (2n = 3x = 36) and diploid (2n = 2x = 24) with some aneuploid genotypes and those of OA hybrids consisted of triploid (2n = 3x = 36) and some aneuploid genotypes. In all cases, it was possible to identify the homoeologous recombinant chromosomes as well as accurately count the number of crossover points, which are called "recombination sites". Recombination sites were estimated in the BC progeny of 71 LA and 41 OA genotypes. In the case of BC progenies of LA hybrids, 248 recombination sites were cytologically localized on 12 different chromosomes of each genome (i.e., L and A). Similarly, 116 recombinant sites were localized on the 12 chromosomes each from the BC progenies of OA hybrids (O and A genomes). Cytological maps were constructed on the basis of the percentages of distances (micrometres) of the recombination sites from the centromeres. Since an Asiatic parent was involved in both hybrids, viz. LA and OA, two maps were constructed for the A genome that were indicated as Asiatic (L) and Asiatic (O). The other two maps were Longiflorum (A) and Oriental (A). Remarkably, the recombination sites were highly unevenly distributed among the different chromosomes of all four maps. Because the recombination sites can be unequivocally identified through GISH, they serve as reliable landmarks and pave the way for assigning molecular markers or desirable genes to chromosomes of Lilium and also monitor introgression of alien segments.

Analysis of the meiosis in the F(1) hybrids of Longiflorum x Asiatic (LA) of lilies (Lilium) using genomic in situ hybridization.

Longiflorum and Asiatic lilies of the genus Lilium of the family Liliaceae are two important groups of modern lily cultivars. One of the main trends of lily breeding is to realize introgression between these groups. With cut style pollination and embryo rescue, distant hybrids between the two groups have been obtained. However, the F(1) hybrids are highly sterile or some of them could produce a small number of 2n gametes, and their BC(1) progenies are usually triploids. Dutch lily breeders have selected many cultivars from these BC(1) progenies based on their variation. It is presumably suggested that such variation could be caused by intergenomic recombination and abnormal meiosis during gamete formation in F(1) hybrids of Longiflorum x Asiatic (LA) hybrids in Lilium. Therefore, the meiotic process of ten F(1) LA hybrids was cytologically investigated using genomic in situ hybridization and traditional cytological methods in the present research. The results showed that: at metaphase I, the homoeologous chromosome pairing among different F(1) hybrids ranged from 2.0 to 11.4 bivalents formed by homoeologous chromosomes per pollen mother cell (PMC), and very few multivalents, and even very few bivalents were formed by two chromosomes within one genome rather than homoeologous chromosomes in some PMCs; at anaphase I, all bivalents were disjoined and most univalents were divided. Both the disjoined bivalents (half-bivalents) and the divided univalents (sister chromatids) moved to the opposite poles, and then formed two groups of chromosomes; because the two resulting half-bivalents retained their axes in the cell undisturbed, many crossover types, including single crossovers, three strand double crossovers, four strand double crossovers, four strand triple crossovers, and four strand multiple crossovers between the non-sister chromatids in the tetrads of bivalents, were clearly inferred by analyzing the breakpoints on the disjoined bivalents. The present investigation not only explained the reason for sterility of the F(1) LA hybrids and the variation of their BC(1) progenies, but also provided a new method to analyze crossover types in other F(1) interspecific hybrids as well.

Ectopic expression of LLAG1, an AGAMOUS homologue from lily (Lilium longiflorum Thunb.) causes floral homeotic modifications in Arabidopsis.

The ABC model for floral development was proposed more than 10 years ago and since then many studies have been performed on model species, such as Arabidopsis thaliana, Antirrhinum majus, and many other species in order to confirm this hypothesis. This led to additional information on flower development and to more complex molecular models. AGAMOUS (AG) is the only C type gene in Arabidopsis and it is responsible for stamen and carpel development as well as floral determinacy. LLAG1, an AG homologue from lily (Lilium longiflorum Thunb.) was isolated by screening a cDNA library derived from developing floral buds. The deduced amino acid sequence revealed the MIKC structure and a high homology in the MADS-box among AG and other orthologues. Phylogenetic analysis indicated a close relationship between LLAG1 and AG orthologues from monocot species. Spatial expression data showed LLAG1 transcripts exclusively in stamens and carpels, constituting the C domain of the ABC model. Functional analysis was carried out in Arabidopsis by overexpression of LLAG1 driven by the CaMV35S promoter. Transformed plants showed homeotic changes in the two outer floral whorls with some plants presenting the second whorl completely converted into stamens. Altogether, these data strongly indicated the functional homology between LLAG1 and AG.