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

Development of EST-SSR markers based on transcriptome sequencing for germplasm evaluation of 65 lilies (Lilium)

BACKGROUND

Lilium genus consists of approximately 100 species and numerous varieties, many of which are interspecific hybrids, which result in a complicated genetic background. The germplasm identification, genetic relationship analysis, and breeding of Lilium rely on exploiting genetic information among different accessions. Hence, an attempt was made to develop new EST-SSR markers and study the molecular divergence among 65 genotypes of Lilium.

METHODS AND RESULTS

A total of 5509 EST-SSRs were identified from the high-throughput sequencing database of L. 'Elodie'. After primer screening, six primer pairs with the most abundant polymorphic bands were selected from 100 primer pairs. Combined with the other 10 reported SSR primers, a total of 16 pairs detected genetic information with an average PIC value of 0.7583. The number of alleles per locus varied from four to 33, the expected heterozygosity varied from 0.3289 to 0.9231, and the observed heterozygosity varied from 0.2857 to 0.5122. Based on the phylogenic results, 22 Asiatic hybrids (A), seven Longiflorum × Asiatic hybrids (LA), as well as two native species were grouped. Eighteen Oriental hybrids (O) and nine Oriental × Trumpet (OT) hybrids, four native species, one LO, and one LL (L. pardalinum × L. longiflorum) variety were grouped.

CONCLUSIONS

Two major clusters were reported and a large number of genotypes were grouped based on UPGMA and STRUCTURE analysis methods. The PIC value as well as other parameters revealed that the EST-SSR markers selected were informative. In addition, the clustering pattern displayed better agreement with the cultivar's pedigree. The newly identified SSRs in this study provides molecular markers for germplasm characterization and genetic diversity for Lilium.

Fire Blight Susceptibility in Lilium spp. Correlates to Sensitivity to Botrytis elliptica Secreted Cell Death Inducing Compounds

Fire blight represents a widespread disease in Lilium spp. and is caused by the necrotrophic Ascomycete Botrytis elliptica. There are >100 Lilium species that fall into distinct phylogenetic groups and these have been used to generate the contemporary commercial genotypes. It is known among lily breeders and growers that different groups of lilies differ in susceptibility to fire blight, but the genetic basis and mechanisms of susceptibility to fire blight are unresolved. The aim of this study was to quantify differences in fire blight susceptibility between plant genotypes and differences in virulence between fungal isolates. To this end we inoculated, in four biological replicates over 2 years, a set of 12 B. elliptica isolates on a panel of 18 lily genotypes representing seven Lilium hybrid groups. A wide spectrum of variation in symptom severity was observed in different isolate-genotype combinations. There was a good correlation between the lesion diameters on leaves and flowers of the Lilium genotypes, although the flowers generally showed faster expanding lesions. It was earlier postulated that B. elliptica pathogenicity on lily is conferred by secreted proteins that induce programmed cell death in lily cells. We selected two aggressive isolates and one mild isolate and collected culture filtrate (CF) samples to compare the cell death inducing activity of their secreted compounds in lily. After leaf infiltration of the CFs, variation was observed in cell death responses between the diverse lilies. The severity of cell death responses upon infiltration of the fungal CF observed among the diverse Lilium hybrid groups correlated well to their fire blight susceptibility. These results support the hypothesis that susceptibility to fire blight in lily is mediated by their sensitivity to B. elliptica effector proteins in a quantitative manner. Cell death-inducing proteins may provide an attractive tool to predict fire blight susceptibility in lily breeding programs.

Transcriptome wide SSR discovery cross-taxa transferability and development of marker database for studying genetic diversity population structure of Lilium species

Lily belongs to family liliaceae, which mainly propagates vegetatively. Therefore, sufficient number of polymorphic, informative, and functional molecular markers are essential for studying a wide range of genetic parameters in Lilium species. We attempted to develop, characterize and design SSR (simple sequence repeat) markers using online genetic resources for analyzing genetic diversity and population structure of Lilium species. We found di-nucleotide repeat motif were more frequent (4684) within 0.14 gb (giga bases) transcriptome than other repeats, of which was two times higher than tetra-repeat motifs. Frequency of di-(AG/CT), tri-(AGG/CTT), tetra-(AAAT), penta-(AGAGG), and hexa-(AGAGGG) repeats was 34.9%, 7.0%, 0.4%, 0.3%, and 0.2%, respectively. A total of 3607 non-redundant SSR primer pairs was designed based on the sequences of CDS, 5'-UTR and 3'-UTR region covering 34%, 14%, 23%, respectively. Among them, a sub set of primers (245 SSR) was validated using polymerase chain reaction (PCR) amplification, of which 167 primers gave expected PCR amplicon and 101 primers showed polymorphism. Each locus contained 2 to 12 alleles on average 0.82 PIC (polymorphic information content) value. A total of 87 lily accessions was subjected to genetic diversity analysis using polymorphic SSRs and found to separate into seven groups with 0.73 to 0.79 heterozygosity. Our data on large scale SSR based genetic diversity and population structure analysis may help to accelerate the breeding programs of lily through utilizing different genomes, understanding genetics and characterizing germplasm with efficient manner.

Transcriptome Sequencing Analysis Provides Insights Into the Response to Fusarium oxysporum in Lilium pumilum

Lily basal rot, caused by Fusarium oxysporum f. sp. lilii, is one of the most serious diseases of lily. Although the lily germplasm which is resistant to F. oxysporum has been used in disease-resistant breeding, few studies on its molecular mechanism of disease resistance have been reported. To comprehensively study the mechanism of resistance to F. oxysporum, transcriptome sequencings of root tissues from Lilium pumilum inoculated with F. oxysporum or sterile water for 6, 12, or 24 h were performed. A total of 50 GB of data were obtained from the transcriptome sequencings of the 6 L. pumilum samples, and 217 098 Unigenes were obtained after the de novo assembly, of which 38.36% Unigenes were annotated. The sequencing results showed that the numbers of differentially expressed genes at 6, 12, and 24 h after inoculation compared with the control were 111, 254, and 2500, respectively. The functional enrichment analysis of the differentially expressed genes showed that several pathways were involved in responses of L. pumilum, mainly including starch and sucrose metabolism, glycolysis/gluconeogenesis, phenylpropanoid biosynthesis, plant hormone signal transduction, flavonoid biosynthesis, vitamin B6 (VB6) biosynthesis, acid biosynthesis, proteasome, and ribosome. Transcription factor analysis revealed that the WRKY and ERF families played important roles in responses of L. pumilum to F. oxysporum. The results of this study elucidate the molecular responses to F. oxysporum in lily and lay a theoretical foundation for improving lily breeding and strategies for lily basal rot resistance.

Exploration and Exploitation of Novel SSR Markers for Candidate Transcription Factor Genes in Lilium Species

Lilies (Lilium sp.) are commercially important horticultural crops widely cultivated for their flowers and bulbs. Here, we conducted large-scale data mining of the lily transcriptome to develop transcription factor (TF)-associated microsatellite markers (TFSSRs). Among 216,768 unigenes extracted from our sequence data, 6966 unigenes harbored simple sequence repeats (SSRs). Seventy-one SSRs were associated with TF genes, and these were used to design primers and validate their potential as markers. These 71 SSRs were accomplished with 31 transcription factor families; including bHLH, MYB, C2H2, ERF, C3H, NAC, bZIP, and so on. Fourteen highly polymorphic SSRs were selected based on Polymorphic Information Content (PIC) values and used to study genetic diversity and population structure in lily accessions. Higher genetic diversity was observed in Longiflorum compared to Oriental and Asiatic populations. Lily accessions were divided into three sub-populations based in our structure analysis, and an un-rooted neighbor-joining tree effectively separated the accessions according to Asiatic, Oriental, and Longiflorum subgroups. Finally, we showed that 46 of the SSR-associated genes were differentially expressed in response to Botrytis elliptica infection. Thus, our newly developed TFSSR markers represent a powerful tool for large-scale genotyping, high-density and comparative mapping, marker-aided backcrossing, and molecular diversity analysis of Lilium sp.

Spatiotemporal Transcriptome Analysis Provides Insights into Bicolor Tepal Development in Lilium "Tiny Padhye"

The bicolor Asiatic hybrid lily cultivar "Tiny Padhye" is an attractive variety because of its unique color pattern. During its bicolor tepal development, the upper tepals undergo a rapid color change from green to white, while the tepal bases change from green to purple. However, the molecular mechanisms underlying these changes remain largely uncharacterized. To systematically investigate the dynamics of the lily bicolor tepal transcriptome during development, we generated 15 RNA-seq libraries from the upper tepals (S2-U) and basal tepals (S1-D, S2-D, S3-D, and S4-D) of Lilium "Tiny Padhye." Utilizing the Illumina platform, a total of 295,787 unigenes were obtained from 713.12 million high-quality paired-end reads. A total of 16,182 unigenes were identified as differentially expressed genes during tepal development. Using Kyoto Encyclopedia of Genes and Genomes pathway analysis, candidate genes involved in the anthocyanin biosynthetic pathway (61 unigenes), and chlorophyll metabolic pathway (106 unigenes) were identified. Further analyses showed that most anthocyanin biosynthesis genes were transcribed coordinately in the tepal bases, but not in the upper tepals, suggesting that the bicolor trait of "Tiny Padhye" tepals is caused by the transcriptional regulation of anthocyanin biosynthetic genes. Meanwhile, the high expression level of chlorophyll degradation genes and low expression level of chlorophyll biosynthetic genes resulted in the absence of chlorophylls from "Tiny Padhye" tepals after flowering. Transcription factors putatively involved in the anthocyanin biosynthetic pathway and chlorophyll metabolism in lilies were identified using a weighted gene co-expression network analysis and their possible roles in lily bicolor tepal development were discussed. In conclusion, these extensive transcriptome data provide a platform for elucidating the molecular mechanisms of bicolor tepals in lilies and provide a basis for similar research in other closely related species.

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.