Genome-wide characterization of simple sequence repeats in cucumber (Cucumis sativus L.)

A high-quality cucumber genome assembly enhances computational comparative genomics

Genetic variation is expressed by the presence of polymorphisms in compared genomes of individuals that can be transferred to next generations. The aim of this work was to reveal genome dynamics by predicting polymorphisms among the genomes of three individuals of the highly inbred B10 cucumber (Cucumis sativus L.) line. In this study, bioinformatic comparative genomics was used to uncover cucumber genome dynamics (also called real-time evolution). We obtained a new genome draft assembly from long single molecule real-time (SMRT) sequencing reads and used short paired-end read data from three individuals to analyse the polymorphisms. Using this approach, we uncovered differentiation aspects in the genomes of the inbred B10 line. The newly assembled genome sequence (B10v3) has the highest contiguity and quality characteristics among the currently available cucumber genome draft sequences. Standard and newly designed approaches were used to predict single nucleotide and structural variants that were unique among the three individual genomes. Some of the variant predictions spanned protein-coding genes and their promoters, and some were in the neighbourhood of annotated interspersed repetitive elements, indicating that the highly inbred homozygous plants remained genetically dynamic. This is the first bioinformatic comparative genomics study of a single highly inbred plant line. For this project, we developed a polymorphism prediction method with optimized precision parameters, which allowed the effective detection of small nucleotide variants (SNVs). This methodology could significantly improve bioinformatic pipelines for comparative genomics and thus has great practical potential in genomic metadata handling.

A leaf shape mutant provides insight into PINOID Serine/Threonine Kinase function in cucumber (Cucumis sativus L.)

Optimizing leaf shape is a major challenge in efforts to develop an ideal plant type. Cucumber leaf shapes are diverse; however, the molecular regulatory mechanisms underlying leaf shape formation are unknown. In this study, we obtained a round leaf mutant (rl) from an ethyl methanesulfonate-induced mutagenesis population. Genetic analysis revealed that a single recessive gene, rl, is responsible for this mutation. A modified MutMap analysis combined linkage mapping identified a single nucleotide polymorphism within a candidate gene, Csa1M537400, as the mutation underlying the trait. Csa1M537400 encodes a PINOID kinase protein involved in auxin transport. Expression of Csa1M537400 was significantly lower in the rl mutant than in wild type, and it displayed higher levels of IAA (indole-3-acetic acid) in several tissues. Treatment of wild-type plants with an auxin transport inhibitor induced the formation of round leaves, similar to those in the rl mutant. Altered expression patterns of several auxin-related genes in the rl mutant suggest that rl plays a key role in auxin biosynthesis, transport, and response in cucumber. These findings provide insight into the molecular mechanism underlying the regulation of auxin signaling pathways in cucumber, and will be valuable in the development of an ideal plant type.

CsTFL1 inhibits determinate growth and terminal flower formation through interaction with CsNOT2a in cucumber

Cucumber (Cucumis sativus L.) is an important vegetable crop that carries on vegetative growth and reproductive growth simultaneously. Indeterminate growth is favourable for fresh market under protected environments, whereas determinate growth is preferred for pickling cucumber in the once-over mechanical harvest system. The genetic basis of determinacy is largely unknown in cucumber. In this study, map-based cloning of the de locus showed that the determinate growth habit is caused by a non-synonymous SNP in CsTFL1CsTFL1 is expressed in the subapical regions of the shoot apical meristem, lateral meristem and young stems. Ectopic expression of CsTFL1 rescued the terminal flower phenotype in the Arabidopsis tfl1-11 mutant and delayed flowering in wild-type Arabidopsis Knockdown of CsTFL1 resulted in determinate growth and formation of terminal flowers in cucumber. Biochemical analyses indicated that CsTFL1 interacts with a homolog of the miRNA biogenesis gene CsNOT2a; CsNOT2a interacts with FDP. Cucumber CsFT directly interacts with CsNOT2a and CsFD, and CsFD interacts with two 14-3-3 proteins. These data suggest that CsTFL1 competes with CsFT for interaction with CsNOT2a-CsFDP to inhibit determinate growth and terminal flower formation in cucumber.

QTL mapping and genome-wide association study reveal two novel loci associated with green flesh color in cucumber

BACKGROUND

Green flesh color, resulting from the accumulation of chlorophyll, is one of the most important commercial traits for the fruits. The genetic network regulating green flesh formation has been studied in tomato, melon and watermelon. However, little is known about the inheritance and molecular basis of green flesh in cucumber. This study sought to determine the main genomic regions associated with green flesh. Three F₂ and two BC₁ populations derived from the 9110Gt (cultivated cucumber, green flesh color) and PI183967 (wild cucumber, white flesh color) were used for the green flesh genetic analysis. Two F₂ populations of them were further employed to do the map construction and quantitative trait loci (QTL) study. Also, a core cucumber germplasms population was used to do the GWAS analysis.

RESULTS

We identified three indexes, flesh color (FC), flesh extract color (FEC) and flesh chlorophyll content (FCC) in three environments. Genetic analysis indicated that green flesh color in 9110Gt is controlled by a major-effect QTL. We developed two genetic maps with 192 and 174 microsatellite markers respectively. Two novel inversions in Chr1 were identified between cultivated and wild cucumbers. The major-effect QTL, qgf5.1, was identified using FC, FEC and FCC index in all different environments used. In addition, the same qgf5.1, together with qgf3.1, was identified via GWAS. Further investigation of two candidate regions using pairwise LD correlations, combined with genetic diversity of qgf5.1 in natural populations, it was found that Csa5G021320 is the candidate gene of qgf5.1. Geographical distribution revealed that green flesh color formation could be due to the high latitude, which has longer day time to produce the photosynthesis and chlorophyll synthesis during cucumber domestication and evolution.

CONCLUSIONS

We first reported the cucumber green flesh color is a quantitative trait. We detected two novel loci qgf5.1 and qgf3.1, which regulate the green flesh formation in cucumber. The QTL mapping and GWAS approaches identified several candidate genes for further validation using functional genomics or forward genetics approaches. Findings from the present study provide a new insight into the genetic control of green flesh in cucumber.

Mapping and identification of CsUp, a gene encoding an Auxilin-like protein, as a putative candidate gene for the upward-pedicel mutation (up) in cucumber

BACKGROUND

Pedicel orientation can affect the female flower orientation and seed yield in cucumber. A spontaneous mutant possessing upward growth of pedicels was identified in the wild type inbred strain 9930 and named upward-pedicel (up). The morphological and genetic analyses of up were performed in this study. In order to clone the up gene, 933 F₂ individuals and 524 BC₁ individuals derived from C-8-6 (WT) and up were used for map-based cloning.

RESULTS

up was mapped to a 35.2 kb physical interval on chromosome 1, which contains three predicted genes. Sequencing analysis revealed that a 5-bp deletion was found in the second exon of Csa1G535800, and it led to a frameshift mutation resulting in a premature stop codon. The candidate gene of CsUp (Csa1G535800) was further confirmed via genomic and cDNA sequencing in biparental and natural cucumber populations. Sequencing data showed that a 4-bp deletion was found in the sixth exon of Csa1G535800 in CGN19839, another inbred line, and there was also a mutation of an amino acid in Csa1G535800 that could contribute to the upward growth of pedicels in CGN19839. Moreover, it was found that Csa1G535800 exhibited strong expression in the pedicel of WT, suggesting its important role in development of pedicel orientation. Thus, Csa1G535800 was considered to be the candidate gene of CsUp.

CONCLUSIONS

CsUp encodes an Auxilin-like protein and controls pedicel orientation in cucumber. The identification of CsUp may help us to understand the mechanism of pedicel orientation development and allow for investigation of novel functions of Auxilin-like proteins in cucumber.

Microsatellite Borders and Micro-sequence Conservation in Juglans

Walnuts (Juglans spp.) are economically important nut and timber species with a worldwide distribution. Using the published Persian walnut genome as a reference for the assembly of short reads from six Juglans species and several interspecific hybrids, we identified simple sequence repeats in 12 Juglans nuclear and organellar genomes. The genome-wide distribution and polymorphisms of nuclear and organellar microsatellites (SSRs) for most Juglans genomes have not been previously studied. We compared the frequency of nuclear SSR motifs and their lengths across Juglans, and identified section-specific chloroplast SSR motifs. Primer pairs were designed for more than 60,000 SSR-containing sequences based on alignment against assembled scaffold sequences. Of the >60,000 loci, 39,000 were validated by e-PCR using unique primer pairs. We identified primers containing 100% sequence identity in multiple species. Across species, sequence identity in the SSR-flanking regions was generally low. Although SSRs are common and highly dispersed in the genome, their flanking sequences are conserved at about 90 to 95% identity within Juglans and within species. In a few rare cases, flanking sequences are identical across species of Juglans. This comprehensive report of nuclear and organellar SSRs in Juglans and the generation of validated SSR primers will be a useful resource for future genetic analyses, walnut breeding programs, high-level taxonomic evaluations, and genomic studies in Juglandaceae.

PolyMorphPredict: A Universal Web-Tool for Rapid Polymorphic Microsatellite Marker Discovery From Whole Genome and Transcriptome Data

Microsatellites are ubiquitously distributed, polymorphic repeat sequence valuable for association, selection, population structure and identification. They can be mined by genomic library, probe hybridization and sequencing of selected clones. Such approach has many limitations like biased hybridization and selection of larger repeats. In silico mining of polymorphic markers using data of various genotypes can be rapid and economical. Available tools lack in some or other aspects like: targeted user defined primer generation, polymorphism discovery using multiple sequence, size and number limits of input sequence, no option for primer generation and e-PCR evaluation, transferability, lack of complete automation and user-friendliness. They also lack the provision to evaluate published primers in e-PCR mode to generate additional allelic data using re-sequenced data of various genotypes for judicious utilization of previously generated data. We developed the tool (PolyMorphPredict) using Perl, R, Java and launched at Apache which is available at http://webtom.cabgrid.res.in/polypred/. It mines microsatellite loci and computes primers from genome/transcriptome data of any species. It can perform e-PCR using published primers for polymorphism discovery and across species transferability of microsatellite loci. Present tool has been evaluated using five species of different genome size having 21 genotypes. Though server is equipped with genomic data of three species for test run with gel simulation, but can be used for any species. Further, polymorphism predictability has been validated using in silico and in vitro PCR of four rice genotypes. This tool can accelerate the in silico microsatellite polymorphism discovery in re-sequencing projects of any species of plant and animal for their diversity estimation along with variety/breed identification, population structure, MAS, QTL and gene discovery, traceability, parentage testing, fungal diagnostics and genome finishing.

Novel loci fsd6.1 and Csgl3 regulate ultra-high fruit spine density in cucumber

Quantitative Trait Loci (QTL) analysis of multiple populations in multiple environments revealed that the fsd6.2 locus, which includes the candidate gene Csgl3, controls high fruit spine density in natural cucumbers. GWAS identified a novel locus fsd6.1, which regulates ultra-high fruit spine density in combination with Csgl3, and evolved during cucumber domestication. Fruit spine density, a domestication trait, largely influences the commercial value of cucumbers. However, the molecular basis of fruit spine density in cucumber remains unclear. In this study, four populations were derived from five materials, which included three with low fruit spine density, one with high fruit spine density, and one with ultra-high fruit spine density. Fruit spine densities were measured in 15 environments over a span of 6 years. The distributions were bimodal suggesting that fruit spine density is controlled by a major-effect QTL. QTL analysis determined that the same major-effect QTL, fsd6.2, is present in four populations. Fine mapping indicated that Csgl3 is the candidate gene at the fsd6.2 locus. Phylogenetic and geographical distribution analyses revealed that Csgl3 originated from China, which has the highest genetic diversity for fruit spine density. One novel minor-effect QTL, fsd6.1, was detected in the HR and HP populations derived from the cross between 65G and 02245. In addition, GWAS identified a novel locus that colocates with fsd6.1. Inspection of a candidate region of about 18 kb in size using pairwise LD correlations, combined with genetic diversity and phylogenetic analysis of fsd6.1 in natural populations, indicated that Csa6G421750 is the candidate gene responsible for ultra-high fruit spine density in cucumber. This study provides new insights into the origin of fruit spine density and the evolution of high/ultra-high fruit spine density during cucumber domestication.

Development, characterization, and cross-amplification of microsatellite markers for Psammosilene tunicoides (Caryophyllaceae)

PREMISE OF THE STUDY

Psammosilene tunicoides (Caryophyllaceae) is a narrowly distributed and endemic plant species in southwestern China. The overexploitation of natural P. tunicoides has led to the destruction of many populations. Population and genetic studies will provide crucial data for the protection and management of P. tunicoides. In this study, we develop simple sequence repeat markers of P. tunicoides to analyze population diversity.

METHODS AND RESULTS

Microsatellite loci of P. tunicoides were isolated with FIASCO. Eleven polymorphic and 10 monomorphic primers were developed. The 11 polymorphic primers were tested in three P. tunicoides populations, yielding two to nine alleles per locus. Levels of observed heterozygosity varied from 0.000 to 1.000, and levels of expected heterozygosity ranged from 0.000 to 0.615. In addition, three of these loci were successfully amplified, and showed polymorphism, in three Silene species.

CONCLUSIONS

These microsatellite markers can be valuable tools to investigate the genetic diversity and population structure of P. tunicoides.

De Novo Transcriptomic Analysis and Development of EST–SSRs for Styrax japonicus

Styrax japonicus sieb. et Zucc. is widely distributed in China with ornamental and medicinal values. However, the transcriptome of S. japonicus has not yet been reported. In this study, we carried out the first transcriptome analysis of S. japonicus and developed a set of expressed sequence tag–simple sequence repeats (EST–SSRs). We obtained 338,570,222 clean reads in total, of which the mean GC content was 41.58%. In total, 136,071 unigenes were obtained having an average length of 611 bp and 71,226 unigenes were favorably annotated in the database. In total, we identified 55,977 potential EST–SSRs from 38,611 unigenes, of which there was 1 SSR per 6.73 kb. The di-nucleotide repeats (40.40%) were the most identified SSRs. One set of 60 primer pairs was randomly selected, and the amplified products in S. japonicus were validated; 28 primer pairs successfully produced clear amplicons. A total of 21 (35%) polymorphic genic SSR markers were identified between two populations. In total, 15 alleles were detected and the average number was 6. The average of observed heterozygosity and expected heterozygosity was 0.614 and 0.552, respectively. The polymorphism information content (PIC) value fluctuated between 0.074 and 0.855, with a mean value of 0.504, which was also the middle level. This study provides useful information for diversity studies and resource assessments of S. japonicus.

Development of Polymorphic Simple Sequence Repeat Markers using High-Throughput Sequencing in Button Mushroom (Agaricus bisporus)

The white button mushroom (Agaricus bisporus) is one of the most widely cultivated species of edible mushroom. Despite its economic importance, relatively little is known about the genetic diversity of this species. Illumina paired-end sequencing produced 43,871,558 clean reads and 69,174 contigs were generated from five offspring. These contigs were subsequently assembled into 57,594 unigenes. The unigenes were annotated with reference genome in which 6,559 unigenes were associated with clusters, indicating orthologous genes. Gene ontology classification assigned many unigenes. Based on genome data of the five offspring, 44 polymorphic simple sequence repeat (SSR) markers were developed. The major allele frequency ranged from 0.42 to 0.92. The number of genotypes and the number of alleles ranged from 1 to 4, and from 2 to 4, respectively. The observed heterozygosity and the expected heterozygosity ranged from 0.00 to 1.00, and from 0.15 to 0.64, respectively. The polymorphic information content value ranged from 0.14 to 0.57. The genetic distances and UPGMA clustering discriminated offspring strains. The SSR markers developed in this study can be applied in polymorphism analyses of button mushroom and for cultivar discrimination.

De Novo assembly, characterization and development of EST-SSRs from Bletilla striata transcriptomes profiled throughout the whole growing period

Bletilla striata is an endangered orchid that has been used for millennia as a medicinal herb, in cosmetics and as a horticultural plant. To construct the first nucleotide database for this species and to develop abundant EST-SSR markers for facilitating further studies, various tissues and organs of plants in the main developmental stages were harvested for mRNA isolation and subsequent RNA sequencing. A total of 106,054,784 clean reads were generated by using Illumina paired-end sequencing technology. The reads were assembled into 127,261 unigenes by the Trinity package; the unigenes had an average length of 612 bp and an N50 of 957 bp. Of these unigenes, 67,494 (51.86%) were annotated in a series of databases. Of these annotated unigenes, 41,818 and 24,615 were assigned to gene ontology categories and clusters of orthologous groups, respectively. Additionally, 20,764 (15.96%) unigenes were mapped onto 275 pathways using the KEGG database. In addition, 25,935 high-quality EST-SSR primer pairs were developed from the 15,433 unigenes by MISA mining. To validate the accuracy of the newly designed markers, 87 of 100 randomly selected primers were effectively amplified; 63 of those yielded PCR products of the expected size, and 25 yielded products with significant amounts of polymorphism among the 4 landraces. Furthermore, the transferability test of the 25 polymorphic markers was performed in 6 individuals of two closely related genus Phalaenopsis and dendrobium. Which results showed a total of 5 markers can successfully amplified among these populations. This research provides a comprehensive nucleotide database and lays a solid foundation for functional gene mining and genomic research in B. striata. The developed EST-SSR primers could facilitate phylogenetic studies and breeding.

The Landscape of Repetitive Elements in the Refined Genome of Chilli Anthracnose Fungus Colletotrichum truncatum

The ascomycete fungus Colletotrichum truncatum is a major phytopathogen with a broad host range which causes anthracnose disease of chilli. The genome sequencing of this fungus led to the discovery of functional categories of genes that may play important roles in fungal pathogenicity. However, the presence of gaps in C. truncatum draft assembly prevented the accurate prediction of repetitive elements, which are the key players to determine the genome architecture and drive evolution and host adaptation. We re-sequenced its genome using single-molecule real-time (SMRT) sequencing technology to obtain a refined assembly with lesser and smaller gaps and ambiguities. This enabled us to study its genome architecture by characterising the repetitive sequences like transposable elements (TEs) and simple sequence repeats (SSRs), which constituted 4.9 and 0.38% of the assembled genome, respectively. The comparative analysis among different Colletotrichum species revealed the extensive repeat rich regions, dominated by Gypsy superfamily of long terminal repeats (LTRs), and the differential composition of SSRs in their genomes. Our study revealed a recent burst of LTR amplification in C. truncatum, C. higginsianum, and C. scovillei. TEs in C. truncatum were significantly associated with secretome, effectors and genes in secondary metabolism clusters. Some of the TE families in C. truncatum showed cytosine to thymine transitions indicative of repeat-induced point mutation (RIP). C. orbiculare and C. graminicola showed strong signatures of RIP across their genomes and "two-speed" genomes with extensive AT-rich and gene-sparse regions. Comparative genomic analyses of Colletotrichum species provided an insight into the species-specific SSR profiles. The SSRs in the coding and non-coding regions of the genome revealed the composition of trinucleotide repeat motifs in exons with potential to alter the translated protein structure through amino acid repeats. This is the first genome-wide study of TEs and SSRs in C. truncatum and their comparative analysis with six other Colletotrichum species, which would serve as a useful resource for future research to get insights into the potential role of TEs in genome expansion and evolution of Colletotrichum fungi and for development of SSR-based molecular markers for population genomic studies.

Pyramiding of nine transgenes in maize generates high-level resistance against necrotrophic maize pathogens

Key message Nine transgenes from different categories, viz. plant defense response genes and anti-apoptosis genes, played combined roles in maize to inhibit the necrotrophic pathogens Rhizoctonia solani and Bipolaris maydis. Maize sheath blight and southern corn leaf blight are major global threats to maize production. The management of these necrotrophic pathogens has encountered limited success due to the characteristics of their lifestyle. Here, we presented a transgenic pyramiding breeding strategy to achieve nine different resistance genes integrated in one transgenic maize line to combat different aspects of necrotrophic pathogens. These nine genes, selected from two different categories, plant defense response genes (Chi, Glu, Ace-AMP1, Tlp, Rs-AFP2, ZmPROPEP1 and Pti4), and anti-apoptosis genes (Iap and p35), were successfully transferred into maize and further implicated in resistance against the necrotrophic pathogens Rhizoctonia solani and Bipolaris maydis. Furthermore, the transgenic maize line 910, with high expression levels of the nine integrated genes, was selected from 49 lines. Under greenhouse and field trial conditions, line 910 showed significant resistance against maize sheath blight and southern corn leaf blight diseases. Higher-level resistance was obtained after the pyramiding of more resistance transgenes from different categories that function via different mechanisms. The present study provides a successful strategy for the management of necrotrophic pathogens.

Complete resistance to powdery mildew and partial resistance to downy mildew in a Cucumis hystrix introgression line of cucumber were controlled by a co-localized locus

Key message A single recessive gene for complete resistance to powdery mildew and a major-effect QTL for partial resistance to downy mildew were co-localized in a Cucumis hystrix introgression line of cucumber. Downy mildew (DM) and powdery mildew (PM) are two major foliar diseases in cucumber. DM resistance (DMR) and PM resistance (PMR) may share common components; however, the genetic relationship between them remains unclear. IL52, a Cucumis hystrix introgression line of cucumber which has been reported to possess DMR, was recently identified to exhibit PMR as well. In this study, a single recessive gene pm for PMR was mapped to an approximately 468-kb region on chromosome 5 with 155 recombinant inbred lines (RILs) and 193 F₂ plants derived from the cross between a susceptible line 'changchunmici' and IL52. Interestingly, pm was co-localized with the major-effect DMR QTL dm5.2 confirmed by combining linkage analysis and BSA-seq, which was consistent with the observed linkage of DMR and PMR in IL52. Further, phenotype-genotype correlation analysis of DMR and PMR in the RILs indicated that the co-localized locus pm/dm5.2 confers complete resistance to PM and partial resistance to DM. Seven candidate genes for DMR were identified within dm5.2 by BSA-seq analysis, of which Csa5M622800.1, Csa5M622830.1 and Csa5M623490.1 were also the same candidate genes for PMR. A single nucleotide polymorphism that is present in the 3' untranslated region (3'UTR) of Csa5M622830.1 co-segregated perfectly with PMR. The GATA transcriptional factor gene Csa5M622830.1 may be a likely candidate gene for DMR and PMR. This study has provided a clear evidence for the relationship between DMR and PMR in IL52 and sheds new light on the potential value of IL52 for cucumber DMR and PMR breeding program.

Genome-wide characterization of microsatellites and genetic diversity assessment of spinach in the Chinese germplasm collection

Spinach is a nutritional leafy green vegetable, and it also serves as a model species for studying sex chromosome evolution. Genetic marker development and genome structure analysis are important in breeding practice and theoretical evolution studies of spinach. In this study, the frequency and distribution of different microsatellites in the recently released draft spinach genome were characterized. A total of 261,002 perfect microsatellites were identified (estimated frequency: ~262.1 loci/Mbp). The most abundant microsatellites were tetranucleotide and trinucleotide, accounting for 33.2% and 27.7% of the total number of microsatellites, respectively. A total of 105 primer pairs were designed and screened, and 34 were polymorphic among the detected spinach cultivars. Combined with seven primer sets developed previously, 41 primer pairs were used to investigate genetic diversity among 43 spinach cultivars in China. The average polymorphism information content value of the 41 markers was 0.43, representing an intermediate level. The spinach cultivars had a low genetic diversity, and no detectable common factors were shared by each group in the UPGMA dendrogram. This study's findings facilitate further investigations on the organization of the microsatellites in spinach genome and provide clues for future breeding applications of spinach in China.

QTL Analysis for Downy Mildew Resistance in Cucumber Inbred Line PI 197088

Downy mildew (DM), caused by Pseudoperonospora cubensis, is one of the major foliar diseases prevailing in cucumber-growing areas. The mechanism of DM resistance in cucumber, particularly the plant introduction (PI) 197088 from India, is presently unclear. Quantitative trait locus (QTL) mapping is an efficient approach to studying DM resistance genes in cucumber. In this study, we performed QTL mapping for DM resistance in PI 197088 with 183 F₂-derived F₃ (F_2:3) families from the cross between PI 197088 (DM resistant) and Changchunmici (DM susceptible). A linkage map was constructed using 141 simple sequence repeat markers. Phenotypic data were collected from seven independent experiments. In total, five QTL were detected on chromosomes 1, 3, 4, and 5 with DM resistance contributed by PI 197088. The QTL on chromosome 4, dm4.1, was reproducibly detected in all indoor experiments, which could explain 27% of the phenotypic variance detected. Additionally, dm1.1 and dm5.2 showed moderate effects, while dm3.1 and dm5.1 were minor-effect QTL. This study revealed the unique genetic architecture of DM resistance in PI 197088, which may provide important guidance for efficient use in cucumber breeding for DM resistance.

STAYGREEN (CsSGR) is a candidate for the anthracnose (Colletotrichum orbiculare) resistance locus cla in Gy14 cucumber

Map-based cloning identified a candidate gene for resistance to the anthracnose fungal pathogen Colletotrichum orbiculare in cucumber, which reveals a novel function for the highly conserved STAYGREEN family genes for host disease resistance in plants. Colletotrichum orbiculare is a hemibiotrophic fungal pathogen that causes anthracnose disease in cucumber and other cucurbit crops. No host resistance genes against the anthracnose pathogens have been cloned in crop plants. Here, we reported fine mapping and cloning of a resistance gene to the race 1 anthracnose pathogen in cucumber inbred lines Gy14 and WI 2757. Phenotypic and QTL analysis in multiple populations revealed that a single recessive gene, cla, was underlying anthracnose resistance in both lines, but WI2757 carried an additional minor-effect QTL. Fine mapping using 150 Gy14 × 9930 recombinant inbred lines and 1043 F₂ individuals delimited the cla locus into a 32 kb region in cucumber Chromosome 5 with three predicted genes. Multiple lines of evidence suggested that the cucumber STAYGREEN (CsSGR) gene is a candidate for the anthracnose resistance locus. A single nucleotide mutation in the third exon of CsSGR resulted in the substitution of Glutamine in 9930 to Arginine in Gy14 in CsSGR protein which seems responsible for the differential anthracnose inoculation responses between Gy14 and 9930. Quantitative real-time PCR analysis indicated that CsSGR was significantly upregulated upon anthracnose pathogen inoculation in the susceptible 9930, while its expression was much lower in the resistant Gy14. Investigation of allelic diversities in natural cucumber populations revealed that the resistance allele in almost all improved cultivars or breeding lines of the U.S. origin was derived from PI 197087. This work reveals an unknown function for the highly conserved STAYGREEN (SGR) family genes for host disease resistance in plants.

Mutations in CsPID encoding a Ser/Thr protein kinase are responsible for round leaf shape in cucumber (Cucumis sativus L.)

Two round-leaf mutants, rl-1 and rl-2, were identified from EMS-induced mutagenesis. High throughput sequencing and map-based cloning suggested CsPID encoding a Ser/Thr protein kinase as the most possible candidate for rl-1. Rl-2 was allelic to Rl-1. Leaf shape is an important plant architecture trait that is affected by plant hormones, especially auxin. In Arabidopsis, PINOID (PID), a regulator for the auxin polar transporter PIN (PIN-FORMED) affects leaf shape formation, but this function of PID in crop plants has not been well studied. From an EMS mutagenesis population, we identified two round-leaf (rl) mutants, C356 and C949. Segregation analysis suggested that both mutations were controlled by single recessive genes, rl-1 and rl-2, respectively. With map-based cloning, we show that CsPID as the candidate gene of rl-1; a non-synonymous SNP in the second exon of CsPID resulted in an amino acid substitution and the round leaf phenotype. As compared in the wild type plant, CsPID had significantly lower expression in the root, leaf and female flowers in C356, which may result in the less developed roots, round leaves and abnormal female flowers, respectively in the rl-1 mutant. Among the three copies of PID genes, CsPID, CsPID2 and CSPID2L (CsPID2-like) in the cucumber genome, CsPID was the only one with significantly differential expression in adult leaves between WT and C356 suggesting CsPID plays a main role in leaf shape formation. The rl-2 mutation in C949 was also cloned, which was due to another SNP in a nearby location of rl-1 in the same CsPID gene. The two round leaf mutants and the work presented herein provide a good foundation for understanding the molecular mechanisms of CsPID in cucumber leaf development.

Fine Mapping and Candidate Gene Prediction for White Immature Fruit Skin in Cucumber (Cucumis sativus L.)

In this study, a single recessive gene (designated w₀) was identified to control the white immature fruit color. Genetic mapping with simple sequence repeats (SSR) markers located the w₀ gene in the distal region of cucumber chromosome 3 (Chr.3). Fine mapping was then conducted using the method of draft genome scaffold-assisted chromosome walking with 7304 F₂ individuals, which allowed for the assignment of the gene locus to a 100.3 kb genomic DNA region with two flanking markers, Q138 and Q193. Thirteen candidate genes were predicted in the 100.3 kb region. Quantitative real-time polymerase chain reaction (qRT-PCR) analysis showed that the expression of the Csa3G904140 gene, which encodes a two-component response regulator-like protein, was much higher in the immature fruit skin of the green parental line (Q1) than in the white parental line (H4). A coding sequence analysis suggested that a single-base insertion occurred at the ninth exon, resulting in a frameshift mutation in Csa3G904140 of H4, and the mutation was consistent with the phenotype in 17 green/white germplasms. Therefore, Csa3G904140 was taken as the likely candidate gene controlling the immature fruit color of cultivated cucumber. This study will contribute to the cloning of candidate genes and the development of white cucumber cultivars using marker-assisted breeding.

Genetic Evaluation of Natural Populations of the Endangered Conifer Thuja koraiensis Using Microsatellite Markers by Restriction-Associated DNA Sequencing

Thuja koraiensis Nakai is an endangered conifer of high economic and ecological value in Jilin Province, China. However, studies on its population structure and conservation genetics have been limited by the lack of genomic data. Here, 37,761 microsatellites (simple sequence repeat, SSR) were detected based on 875,792 de novo-assembled contigs using a restriction-associated DNA (RAD) approach. Among these SSRs, 300 were randomly selected to test for polymorphisms and 96 obtained loci were able to amplify a fragment of expected size. Twelve polymorphic SSR markers were developed to analyze the genetic diversity and population structure of three natural populations. High genetic diversity (mean NA = 5.481, HE = 0.548) and moderate population differentiation (pairwise Fst = 0.048–0.078, Nm = 2.940–4.958) were found in this species. Molecular variance analysis suggested that most of the variation (83%) existed within populations. Combining the results of STRUCTURE, principal coordinate, and neighbor-joining analysis, the 232 individuals were divided into three genetic clusters that generally correlated with their geographical distributions. Finally, appropriate conservation strategies were proposed to protect this species. This study provides genetic information for the natural resource conservation and utilization of T. koraiensis and will facilitate further studies of the evolution and phylogeography of the species.

QTL mapping of downy and powdery mildew resistances in PI 197088 cucumber with genotyping-by-sequencing in RIL population

Host resistances in PI 197088 cucumber to downy and powdery mildew pathogens are conferred by 11 (3 with major effect) and 4 (1 major effect) QTL, respectively, and three of which are co-localized. The downy mildew (DM) and powdery mildew (PM) are the two most important foliar diseases of cucurbit crops worldwide. The cucumber accession PI 197088 exhibits high-level resistances to both pathogens. Here, we reported QTL mapping results for DM and PM resistances with 148 recombinant inbred lines from a cross between PI 197088 and the susceptible line 'Coolgreen'. Phenotypic data on responses to natural DM and PM infection were collected in multi-year and multi-location replicated field trials. A high-density genetic map with 2780 single nucleotide polymorphisms (SNPs) from genotyping-by-sequencing and 55 microsatellite markers was developed, which revealed genomic regions with segregation distortion and mis-assemblies in the '9930' cucumber draft genome. QTL analysis identified 11 and 4 QTL for DM and PM resistances accounting for more than 73.5 and 63.0% total phenotypic variance, respectively. Among the 11 DM resistance QTL, dm5.1, dm5.2, and dm5.3 were major-effect contributing QTL, whereas dm1.1, dm2.1, and dm6.2 conferred susceptibility. Of the 4 QTL for PM resistance, pm5.1 was the major-effect QTL explaining 32.4% phenotypic variance and the minor-effect QTL pm6.1 contributed to disease susceptibility. Three PM QTL, pm2.1, pm5.1, and pm6.1, were co-localized with DM QTL dm2.1, dm5.2, and dm6.1, respectively, which was consistent with the observed linkage of PM and DM resistances in PI 197088. The genetic architecture of DM resistance in PI 197088 and another resistant line WI7120 (PI 330628) was compared, and the potential of using PI 197088 in cucumber breeding for downy and powdery mildew resistances is discussed.

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.

Mining and Development of Novel SSR Markers Using Next Generation Sequencing (NGS) Data in Plants

Microsatellites, or simple sequence repeats (SSRs), are one of the most informative and multi-purpose genetic markers exploited in plant functional genomics. However, the discovery of SSRs and development using traditional methods are laborious, time-consuming, and costly. Recently, the availability of high-throughput sequencing technologies has enabled researchers to identify a substantial number of microsatellites at less cost and effort than traditional approaches. Illumina is a noteworthy transcriptome sequencing technology that is currently used in SSR marker development. Although 454 pyrosequencing datasets can be used for SSR development, this type of sequencing is no longer supported. This review aims to present an overview of the next generation sequencing, with a focus on the efficient use of de novo transcriptome sequencing (RNA-Seq) and related tools for mining and development of microsatellites in plants.

Fine Mapping of CsVYL, Conferring Virescent Leaf Through the Regulation of Chloroplast Development in Cucumber

Leaf color mutants in higher plants are ideal materials for investigating the structure and function of photosynthetic system. In this study, we identified a cucumber vyl (virescent-yellow leaf) mutant in the mutant library, which exhibited reduced pigment contents and delayed chloroplast development process. F₂ and BC₁ populations were constructed from the cross between vyl mutant and cucumber inbred line 'Hazerd' to identify that the vyl trait is controlled by a simply recessive gene designated as CsVYL. The CsVYL gene was mapped to a 3.8 cM interval on chromosome 4 using these 80 F₂ individuals and BSA (bulked segregation analysis) approach. Fine genetic map was conducted with 1542 F₂ plants and narrowed down the vyl locus to an 86.3 kb genomic region, which contains a total of 11 genes. Sequence alignment between the wild type (WT) and vyl only identified one single nucleotide mutation (C→T) in the first exon of gene Csa4G637110, which encodes a DnaJ-like zinc finger protein. Gene Expression analysis confirmed the differences in transcription level of Csa4G637110 between wild type and mutant plants. Map-based cloning of the CsVYL gene could accelerate the study of chloroplast development and chlorophyll synthesis of cucumber.

Comprehensive Characterization of Simple Sequence Repeats in Eggplant (Solanum melongena L.) Genome and Construction of a Web Resource

We have characterized the simple sequence repeat (SSR) markers of the eggplant (Solanum melongena) using a recent high quality sequence of its whole genome. We found nearly 133,000 perfect SSRs, a density of 125.5 SSRs/Mbp, and also about 178,400 imperfect SSRs. Of the perfect SSRs, 15.6% were complex, with two stretches of repeats separated by an intervening block of <100 nt. Di- and trinucleotide SSRs accounted, respectively, for 43 and 37% of the total. The SSRs were classified according to their number of repeats and overall length, and were assigned to their linkage group. We found 2,449 of the perfect SSRs in 2,086 genes, with an overall density of 18.5 SSRs/Mbp across the gene space; 3,524 imperfect SSRs were present in 2,924 genes at a density of 26.7 SSRs/Mbp. Putative functions were assigned via ontology to genes containing at least one SSR. Using this data we developed an "Eggplant Microsatellite DataBase" (EgMiDB) which permits identification of SSR markers in terms of their location on the genome, type of repeat (perfect vs. imperfect), motif type, sequence, repeat number and genomic/gene context. It also suggests forward and reverse primers. We employed an in silico PCR analysis to validate these SSR markers, using as templates two CDS sets and three assembled transcriptomes obtained from diverse eggplant accessions.

Cross-species transferability of EST-SSR markers developed from the transcriptome of Melilotus and their application to population genetics research

Melilotus is one of the most important legume forages, but the lack of molecular markers has limited the development and utilization of Melilotus germplasm resources. In the present study, 151 M clean reads were generated from various genotypes of Melilotus albus using Illumina sequencing. A total of 19,263 potential EST-SSRs were identified from 104,358 unigene sequences. Moreover, 18,182 primer pairs were successfully designed, and 550 primer pairs were selected using criteria of base repeat type, fragment length and annealing temperature. In addition, 550 primer pairs were screened by using PCR amplification products and used to assess polymorphisms in 15 M. albus accessions. A total of 114 primer pairs were detected as being highly polymorphic, and the average polymorphism information content (PIC) value was 0.79. Furthermore, those 114 polymorphic primer pairs were used to evaluate the transferability to 18 species of the genus Melilotus, and 70 EST-SSR markers were found to be transferable among the 18 Melilotus species. According to the UPGMA dendrogram and STRUCTURE analysis, the 18 Melilotus species were classified into three clusters. This study offers a valuable resource for the genetic diversity and molecular assisted breeding of germplasm resources in the genus Melilotus.

Putative Microsatellite DNA Marker-Based Wheat Genomic Resource for Varietal Improvement and Management

Wheat fulfills 20% of global caloric requirement. World needs 60% more wheat for 9 billion population by 2050 but climate change with increasing temperature is projected to affect wheat productivity adversely. Trait improvement and management of wheat germplasm requires genomic resource. Simple Sequence Repeats (SSRs) being highly polymorphic and ubiquitously distributed in the genome, can be a marker of choice but there is no structured marker database with options to generate primer pairs for genotyping on desired chromosome/physical location. Previously associated markers with different wheat trait are also not available in any database. Limitations of in vitro SSR discovery can be overcome by genome-wide in silico mining of SSR. Triticum aestivum SSR database (TaSSRDb) is an integrated online database with three-tier architecture, developed using PHP and MySQL and accessible at http://webtom.cabgrid.res.in/wheatssr/. For genotyping, Primer3 standalone code computes primers on user request. Chromosome-wise SSR calling for all the three sub genomes along with choice of motif types is provided in addition to the primer generation for desired marker. We report here a database of highest number of SSRs (476,169) from complex, hexaploid wheat genome (~17 GB) along with previously reported 268 SSR markers associated with 11 traits. Highest (116.93 SSRs/Mb) and lowest (74.57 SSRs/Mb) SSR densities were found on 2D and 3A chromosome, respectively. To obtain homozygous locus, e-PCR was done. Such 30 loci were randomly selected for PCR validation in panel of 18 wheat Advance Varietal Trial (AVT) lines. TaSSRDb can be a valuable genomic resource tool for linkage mapping, gene/QTL (Quantitative trait locus) discovery, diversity analysis, traceability and variety identification. Varietal specific profiling and differentiation can supplement DUS (Distinctiveness, Uniformity, and Stability) testing, EDV (Essentially Derived Variety)/IV (Initial Variety) disputes, seed purity and hybrid wheat testing. All these are required in germplasm management as well as also in the endeavor of wheat productivity.

Development of Genome-Wide SSR Markers from Angelica gigas Nakai Using Next Generation Sequencing

Angelica gigas Nakai is an important medicinal herb, widely utilized in Asian countries especially in Korea, Japan, and China. Although it is a vital medicinal herb, the lack of sequencing data and efficient molecular markers has limited the application of a genetic approach for horticultural improvements. Simple sequence repeats (SSRs) are universally accepted molecular markers for population structure study. In this study, we found over 130,000 SSRs, ranging from di- to deca-nucleotide motifs, using the genome sequence of Manchu variety (MV) of A. gigas, derived from next generation sequencing (NGS). From the putative SSR regions identified, a total of 16,496 primer sets were successfully designed. Among them, we selected 848 SSR markers that showed polymorphism from in silico analysis and contained tri- to hexa-nucleotide motifs. We tested 36 SSR primer sets for polymorphism in 16 A. gigas accessions. The average polymorphism information content (PIC) was 0.69; the average observed heterozygosity (HO) values, and the expected heterozygosity (HE) values were 0.53 and 0.73, respectively. These newly developed SSR markers would be useful tools for molecular genetics, genotype identification, genetic mapping, molecular breeding, and studying species relationships of the Angelica genus.

Development of Highly Informative Genome-Wide Single Sequence Repeat Markers for Breeding Applications in Sesame and Construction of a Web Resource: SisatBase

The sequencing of the full nuclear genome of sesame (Sesamum indicum L.) provides the platform for functional analyses of genome components and their application in breeding programs. Although the importance of microsatellites markers or simple sequence repeats (SSR) in crop genotyping, genetics, and breeding applications is well established, only a little information exist concerning SSRs at the whole genome level in sesame. In addition, SSRs represent a suitable marker type for sesame molecular breeding in developing countries where it is mainly grown. In this study, we identified 138,194 genome-wide SSRs of which 76.5% were physically mapped onto the 13 pseudo-chromosomes. Among these SSRs, up to three primers pairs were supplied for 101,930 SSRs and used to in silico amplify the reference genome together with two newly sequenced sesame accessions. A total of 79,957 SSRs (78%) were polymorphic between the three genomes thereby suggesting their promising use in different genomics-assisted breeding applications. From these polymorphic SSRs, 23 were selected and validated to have high polymorphic potential in 48 sesame accessions from different growing areas of Africa. Furthermore, we have developed an online user-friendly database, SisatBase (http://www.sesame-bioinfo.org/SisatBase/), which provides free access to SSRs data as well as an integrated platform for functional analyses. Altogether, the reference SSR and SisatBase would serve as useful resources for genetic assessment, genomic studies, and breeding advancement in sesame, especially in developing countries.

Fine mapping identifies CsGCN5 encoding a histone acetyltransferase as putative candidate gene for tendril-less1 mutation (td-1) in cucumber

Next-generation sequencing-aided map-based cloning delimited the cucumber tendril - less1 ( td - 1 ) locus into a 190.7-kb region in chromosome 6 harboring a putative, novel-function candidate gene encoding a histone acetyltransferase ( CsGCN5 ). The tendril initiated from the lateral meristem is an important and characteristic organ for the species in the Cucurbitaceae family including cucumber (Cucumis sativus L.). While the tendril has its evolutionary significance, it also poses a nuisance in cucumber cultivation under protected environments in which tendril-less cucumber has its advantages. From an EMS mutagenesis population, we identified a tendril-less mutant B007, which was controlled by a recessive gene td-1. Through next-generation sequencing-aided map-based cloning, we show CsGCN5 (Cucumis sativus GENERAL CONTROL NONDEREPRESSIBLE 5), a cucumber gene for a histone acetyltransferase as the most possible candidate for td-1. A non-synonymous SNP in the first exon of CsGCN5 resulted in an amino-acid substitution from Asp (D) in the wild type to Asn (N) in the tendril-less mutant. The candidacy of CsGCN5 was further confirmed by multiple lines of evidence in both biparental and natural cucumber populations. Non-significant expression of CsGCN5 in multiple organs was found between the wild type and the mutant. CsGCN5 exhibited strong expression in the tendril of wild-type plants suggesting its important roles in growth and development of plant tendrils. The identification and characterization of the td-1 mutant from the present study provided a useful tool in understanding the molecular mechanisms of tendril organogenesis and investigation of novel functions of the histone acetyltransferase in cucumber.

A mutant in the CsDET2 gene leads to a systemic brassinosteriod deficiency and super compact phenotype in cucumber (Cucumis sativus L.)

A novel dwarf cucumber mutant, scp-2, displays a typical BR biosynthesis-deficient phenotype, which is due to a mutation in CsDET2 for a steroid 5-alpha-reductase. Brassinosteroids (BRs) are a group of plant hormones that play important roles in the development of plant architecture, and extreme dwarfism is a typical outcome of BR-deficiency. Most cucumber (Cucumis sativus L.) varieties have an indeterminate growth habit, and dwarfism may have its value in manipulation of plant architecture and improve production in certain production systems. In this study, we identified a spontaneous dwarf mutant, super compact-2 (scp-2), that also has dark green, wrinkle leaves. Genetic analyses indicated that scp-2 was different from two previously reported dwarf mutants: compact (cp) and super compact-1 (scp-1). Map-based cloning revealed that the mutant phenotype was due to two single nucleotide polymorphism and a single-base insertion in the CsDET2 gene that resulted in a missense mutation in a conserved amino acid and thus a truncated protein lacking the conserved catalytic domains in the predicted steroid 5α-reductase protein. Measurement of endogenous hormone levels indicated a reduced level of brassinolide (BL, a bioactive BR) in scp-2, and the mutant phenotype could be partially rescued by the application of epibrassinolide (EBR). In addition, scp-2 mutant seedlings exhibited dark-grown de-etiolation, and defects in cell elongation and vascular development. These data support that scp-2 is a BR biosynthesis-deficient mutant, and that the CsDET2 gene plays a key role in BR biosynthesis in cucumber. We also described the systemic BR responses and discussed the specific BR-related phenotypes in cucumber plants.

QTL mapping of domestication and diversifying selection related traits in round-fruited semi-wild Xishuangbanna cucumber (Cucumis sativus L. var. xishuangbannanesis)

QTL analysis revealed 11 QTL underlying flowering time and fruit size variation in the semi-wild Xishuangbanna cucumber, of which, FT6.2 and FS5.2 played the most important roles in determining photoperiod-dependent flowering time and round-fruit shape, respectively. Flowering time and fruit size are two important traits in domestication and diversifying selection in cucumber, but their genetic basis is not well understood. Here we reported QTL mapping results on flowering time and fruit size with F₂ and F_2:3 segregating populations derived from the cross between WI7200, a small fruited, early flowering primitive cultivated cucumber and WI7167, a round-fruited, later flowering semi-wild Xishuangbanna (XIS) cucumber. A linkage map with 267 microsatellite marker loci was developed with 138 F₂ plants. Phenotypic data of male and female flowering time, fruit length and diameter and three other traits (mature fruit weight and number, and seedling hypocotyl length) were collected in multiple environments. Three flowering time QTL, FT1.1, FT5.1 and FT6.2 were identified, in which FT6.2 played the most important role in conferring less photoperiod sensitive early flowering during domestication whereas FT1.1 seemed more influential in regulating flowering time within the cultivated cucumber. Eight consensus fruit size QTL distributed in 7 chromosomes were detected, each of which contributed to both longitudinal and radial growth in cucumber fruit development. Among them, FS5.2 on chromosome 5 exhibited the largest effect on the determination of round fruit shape that was characteristic of the WI7167 XIS cucumber. Possible roles of these flowering time and fruit size QTL in domestication of cucumber and crop evolution of the semi-wild XIS cucumber, as well as the genetic basis of round fruit shape in cucumber are discussed.

Characterization of the transcriptome and EST-SSR development in Boea clarkeana, a desiccation-tolerant plant endemic to China

BACKGROUND

Desiccation-tolerant (DT) plants can recover full metabolic competence upon rehydration after losing most of their cellular water (>95%) for extended periods of time. Functional genomic approaches such as transcriptome sequencing can help us understand how DT plants survive and respond to dehydration, which has great significance for plant biology and improving the drought tolerance of crops. Boea clarkeana Hemsl. (Gesneriaceae) is a DT dicotyledonous herb. Its genomic sequences characteristics remain unknown. Based on transcriptomic analyses, polymorphic EST-SSR (simple sequence repeats in expressed sequence tags) molecular primers can be designed, which will greatly facilitate further investigations of the population genetics and demographic histories of DT plants.

METHODS

In the present study, we used the platform Illumina HiSeq™2000 and de novo assembly technology to obtain leaf transcriptomes of B. clarkeana and conducted a BLASTX alignment of the sequencing data and protein databases for sequence classification and annotation. Then, based on the sequence information, the EST-SSR markers were developed, and the functional annotation of ESTs containing polymorphic SSRs were obtained through BLASTX.

RESULTS

A total of 91,449 unigenes were generated from the leaf cDNA library of B. clarkeana. Based on a sequence similarity search with a known protein database, 72,087 unigenes were annotated. Among the annotated unigenes, a total of 71,170 unigenes showed significant similarity to the known proteins of 463 popular model species in the Nr database, and 59,962 unigenes and 32,336 unigenes were assigned to Gene Ontology (GO) classifications and Cluster of Orthologous Groups (COG), respectively. In addition, 44,924 unigenes were mapped in 128 KEGG pathways. Furthermore, a total of 7,610 unigenes with 8,563 microsatellites were found. Seventy-four primer pairs were selected from 436 primer pairs designed for polymorphism validation. SSRs with higher polymorphism rates were concentrated on dinucleotides, pentanucleotides and hexanucleotides. Finally, 17 pairs with stable, highly polymorphic loci were selected for polymorphism screening. There was a total of 65 alleles, with 2-6 alleles at each locus. Primarily due to the unique biological characteristics of plants, the HE (0-0.196), HO (0.082-0.14) and PIC (0-0.155) per locus were very low. The functional annotation distribution centered on ESTs containing di- and tri-nucleotide SSRs, and the ESTs containing primers BC2, BC4 and BC12 were annotated to vegetative dehydration/desiccation pathways.

DISCUSSION

This work is the first genetic study of B. clarkeana as a new plant resource of DT genes. A substantial number of transcriptome sequences were generated in this study. These sequences are valuable resources for gene annotation and discovery as well as molecular marker development. These sequences could also provide a valuable basis for future molecular studies of B. clarkeana.

Development and Molecular Characterization of Novel Polymorphic Genomic DNA SSR Markers in Lentinula edodes

Sixteen genomic DNA simple sequence repeat (SSR) markers of Lentinula edodes were developed from 205 SSR motifs present in 46.1-Mb long L. edodes genome sequences. The number of alleles ranged from 3-14 and the major allele frequency was distributed from 0.17-0.96. The values of observed and expected heterozygosity ranged from 0.00-0.76 and 0.07-0.90, respectively. The polymorphic information content value ranged from 0.07-0.89. A dendrogram, based on 16 SSR markers clustered by the paired hierarchical clustering' method, showed that 33 shiitake cultivars could be divided into three major groups and successfully identified. These SSR markers will contribute to the efficient breeding of this species by providing diversity in shiitake varieties. Furthermore, the genomic information covered by the markers can provide a valuable resource for genetic linkage map construction, molecular mapping, and marker-assisted selection in the shiitake mushroom.

Round fruit shape in WI7239 cucumber is controlled by two interacting quantitative trait loci with one putatively encoding a tomato SUN homolog

QTL analysis revealed two interacting loci, FS1.2 and FS2.1, underlying round fruit shape in WI7239 cucumber; CsSUN , a homolog of tomato fruit shape gene SUN , was a candidate for FS1.2. Fruit size is an important quality and yield trait in cucumber, but its genetic basis remains poorly understood. Here we reported QTL mapping results on fruit size with segregating populations derived from the cross between WI7238 (long fruit) and WI7239 (round fruit) inbred cucumber lines. Phenotypic data of fruit length and diameter were collected at anthesis, immature and mature fruit stages in four environments. Ten major-effect QTL were detected for six traits; synthesis of information from these QTL supported two genes, FS1.2 and FS2.1, underlying fruit size variation in the examined populations. Under the two-gene model, deviation from expected segregation ratio in fruit length and diameter among segregating populations was observed, which could be explained mainly by the interactions between FS1.2 and FS2.1, and segregation distortion in the FS2.1 region. Genome-wide candidate gene search identified CsSUN, a homolog of the tomato fruit shape gene SUN, as the candidate for FS1.2. The round-fruited WI7239 had a 161-bp deletion in the first exon of CsSUN, and its expression in WI7239 was significantly lower than that in WI7238. A marker derived from this deletion was mapped at the peak location of FS1.2 in QTL analysis. Comparative analysis suggested the melon gene CmSUN-14, a homolog of CsSUN as a candidate of the fl2/fd2/fw2 QTL in melon. This study revealed the unique genetic architecture of round fruit shape in WI7239 cucumber. It also highlights the power of QTL analysis for traits with a simple genetic basis but their expression is complicated by other factors.

Elucidation of the molecular responses of a cucumber segment substitution line carrying Pm5.1 and its recurrent parent triggered by powdery mildew by comparative transcriptome profiling

BACKGROUND

Powdery mildew (PM) is one of the most severe fungal diseases of cucurbits, but the molecular mechanisms underlying PM resistance in cucumber remain elusive. In this study, we developed a PM resistant segment substitution line SSL508-28 that carried a segment on chromosome five representing the Pm5.1 locus from PM resistant donor Jin5-508 using marker-assisted backcrossing of an elite PM susceptible cucumber inbred line D8.

RESULTS

Whole-genome resequencing of SSL508-28, Jin5-508 and D8 was performed to identify the exact boundaries of the breakpoints for this introgression because of the low density of available single sequence repeat markers. This led to the identification of a ~6.8 Mb substituted segment predicted to contain 856 genes. RNA-seq was used to study gene expression differences in PM treated (plants harvested 48 h after inoculation) and untreated (control) SSL508-28 and D8 lines. Exactly 1,248 and 1,325 differentially expressed genes (DEGs) were identified in SSL508-28 and D8, respectively. Of those, 88 DEGs were located in the ~6.8 Mb segment interval. Based on expression data and annotation, we identified 8 potential candidate genes that may participate in PM resistance afforded by Pm5.1, including two tandemly arrayed genes encoding receptor protein kinases, two transcription factors, two genes encoding remorin proteins, one gene encoding a P-type ATPase and one gene encoding a 70 kDa heat shock protein. The transcriptome data also revealed a complex regulatory network for Pm5.1-mediated PM resistance that may involve multiple signal regulators and transducers, cell wall modifications and the salicylic acid signaling pathway.

CONCLUSION

These findings shed light on the cucumber PM defense mechanisms mediated by Pm5.1 and provided valuable information for the fine mapping of Pm5.1 and breeding of cucumber with enhanced resistance to PM.

De Novo Assembly of Transcriptome and Development of Novel EST-SSR Markers in Rhododendron rex Lévl. through Illumina Sequencing

Transcriptome sequences generated by next-generation sequencing (NGS) technologies can be utilized to rapidly detect and characterize a large number of gene-based microsatellites from different plants. Rhododendron rex Lévl. is a perennial woody species from the family Ericaceae and an endangered plant with high ornamental value endemic to Southwestern China. Nevertheless, the genetic and genomic information of R. rex remain unknown. In this study, we performed transcriptome sequencing for R. rex leaf samples, and generated large transcript sequences for functional characterization and development gene-associated SSR markers. A total of 164,242 unigenes were assembled and 115,089 (70.07%) unigenes were successfully annotated in public databases. In addition, a total of 15,314 potential EST-SSRs were identified, and the frequency of SSRs in the R. rex unigenes was 9.32%, with an average of one EST-SSR per 5.65 kb. The most abundant type was repeated di-nucleotide (54.63%), followed by mono- (26.03%) and tri-nucleotide (18.51%) repeats. Based on the SSR-containing sequence, 100 primer pairs were randomly selected and synthesized and used for assessment of the polymorphism. Thirty-six primer pairs were polymorphic and revealed polymorphism among 20 individuals from four R. rex populations. A total of 197 alleles were identified, with an average of 5.472 alleles per locus. The Polymorphism Information Content ranged from 0.154 to 0.870, with a mean of 0.482. The newly developed EST-SSR markers exhibited high transferability (58.33-83.33%) among the six subgenera. Thus, these novel EST-SSR markers developed would provide valuable sequence resources for population structure, genetic diversity analysis, and genetic resource assessments of R. rex and its related species.

Genome-Wide Analysis of Simple Sequence Repeats in Bitter Gourd (Momordica charantia)

Bitter gourd (Momordica charantia) is widely cultivated as a vegetable and medicinal herb in many Asian and African countries. After the sequencing of the cucumber (Cucumis sativus), watermelon (Citrullus lanatus), and melon (Cucumis melo) genomes, bitter gourd became the fourth cucurbit species whose whole genome was sequenced. However, a comprehensive analysis of simple sequence repeats (SSRs) in bitter gourd, including a comparison with the three aforementioned cucurbit species has not yet been published. Here, we identified a total of 188,091 and 167,160 SSR motifs in the genomes of the bitter gourd lines 'Dali-11' and 'OHB3-1,' respectively. Subsequently, the SSR content, motif lengths, and classified motif types were characterized for the bitter gourd genomes and compared among all the cucurbit genomes. Lastly, a large set of 138,727 unique in silico SSR primer pairs were designed for bitter gourd. Among these, 71 primers were selected, all of which successfully amplified SSRs from the two bitter gourd lines 'Dali-11' and 'K44'. To further examine the utilization of unique SSR primers, 21 SSR markers were used to genotype a collection of 211 bitter gourd lines from all over the world. A model-based clustering method and phylogenetic analysis indicated a clear separation among the geographic groups. The genomic SSR markers developed in this study have considerable potential value in advancing bitter gourd research.

Inheritance and quantitative trail loci mapping of adventitious root numbers in cucumber seedlings under waterlogging conditions

The hypocotyl-derived adventitious root (AR) is an important morphological acclimation to waterlogging stress; however, its genetic basis has not been adequately understood. In the present study, a mixed major gene plus polygene inheritance model was used to analyze AR numbers (ARN) 7 days after waterlogging treatment in six generations (P₁, P₂, F₁, B₁, B_2, and F₂), using cucumber waterlogging tolerant line Zaoer-N and sensitive Pepino as parents. The results showed that the genetic model D-4, mixed one negative dominance major gene and additive-dominance polygenes, is the best-fitting genetic model for waterlogging-triggered ARN phenotype. A genetic linkage map spanning 550.8 cM and consisting of 149 simple sequence repeat (SSR) markers segregating into seven linkage groups was constructed. Three QTLs (ARN3.1, ARN5.1, and ARN6.1) distributed on chromosomes 3, 5, and 6 were identified by composite interval mapping. The major-effect QTL, ARN6.1, located between SSR12898 and SSR04751, was the only locus detected in three seasons, with least likelihood (LOD) scores of 8.8, 10.4, and 9.5 and account for 17.6, 24, and 19.8% of the phenotypic variance, respectively. Using five additional single nucleotide polymorphism (SNP) makers, the ARN6.1 was narrowed down to a 0.79 Mb interval franked by SSR12898 and SNP25558853. Illumina RNA-sequencing data generated on hypocotyls of two parents 48 h after waterlogging treatment revealed 15 genes in the 0.79 Mb interval were differentially expressed, including Csa6G503880 encoding a salicylic acid methyl transferase-like protein, Csa6G504590 encoding a cytochrome P450 monooxygenase, and Csa6G505230 encoding a heavy metal-associated protein. Our findings shed light on the genetic architecture underlying adventitious rooting during waterlogging stress in cucumber, and provide a list of potential gene targets for further elucidating waterlogging tolerance in plants.

Genetic Analysis and QTL Mapping of Fruit Peduncle Length in Cucumber (Cucumis sativus L.)

Mechanized harvesting of cucumbers offers significant advantages compared to manual labor as both shortages and costs of labor increase. However the efficient use of machines depends on breeding plants with longer peduncles, but the genetic and molecular basis of fruit peduncle development in cucumber is not well understood. In this study, F2 populations were developed from a cross between two inbred lines, 1101 with a long peduncle and 1694 with a short peduncle. These were grown at two field sites, Hainan, with a tropical marine climate, in December 2014, and Beijing, with a warm temperate climate, in May 2015. Electron microscope examination of the pith cells in the peduncles of the two parental lines showed that line 1101 had significantly greater numbers of smaller cells compared to line 1694. The inheritance of cucumber fruit peduncle length (FPL) was investigated by the mixed major gene and polygene inheritance model. Genetic analysis indicated that FPL in cucumber is quantitatively inherited and controlled by one additive major gene and additive-dominant polygenes (D-2 model). A total of 1460 pairs of SSR (simple sequence repeat) primers were analyzed to identify quantitative trait loci (QTLs). Two similar genetic maps with 78 SSR markers which covered 720.6 cM in seven linkage groups were constructed based on two F2 populations. QTL analysis from the data collected at the two field sites showed that there are two minor QTLs on chromosome 1, named qfpl1.1 and qfpl1.2, and one major QTL on chromosome 6, named qfpl6.1. The marker UW021226, which was the closest one to qfpl6.1, had an accuracy rate of 79.0% when tested against plants selected from populations of the two parents. The results from this study provide insights into the inheritance and molecular mechanism of the variation of FPL in cucumber, and further research will be carried out to fine map qfpl6.1 to develop more accurate markers for MAS breeding.

QTL mapping of parthenocarpic fruit set in North American processing cucumber

Through a novel phenotyping method, four QTLs were consistently associated with increased parthenocarpic fruit set in North American processing cucumber that accounted for over 75 % of observed phenotypic variation. Parthenocarpy is a desirable trait with potential for increasing yield and quality in processing cucumber production. Although many successful parthenocarpic fresh market cucumber varieties have been developed, the genetic and molecular mechanisms behind parthenocarpic expression in cucumber remain largely unknown. Since parthenocarpy is an important yield component, it is difficult to separate the true parthenocarpic character from other yield related traits. In the present study, we developed a novel phenotypic approach for parthenocarpic fruit set focusing on early fruit development. Two hundred and five F₃ families derived from a cross between the highly parthenocarpic line 2A and low parthenocarpic line Gy8 were phenotypically evaluated in three greenhouse experiments. Seven QTLs associated with parthenocarpic fruit set were detected. Among them, one each on chromosomes 5 and 7 (parth5.1 and parth7.1) and two on chromosome 6 (parth6.1 and parth6.2) were consistently identified in all experiments, but their relative contribution to the total phenotypic variation was dependent on plant growth stages. While each of the four QTLs had almost equal contribution to the expression of the trait at commercial harvest stage, parth7.1 played an important role in early parthenocarpic fruit set. The results suggested that parthenocarpic fruit set can be accurately evaluated with as few as 20 nodes of growth. The QTLs identified in this study for parthenocarpic fruit set are a valuable resource for cucumber breeders interested in developing parthenocarpic cultivars and to researchers interested in the genetic and molecular mechanisms of parthenocarpic fruit set.

A genome-wide analysis of simple sequence repeats in Apis cerana and its development as polymorphism markers

The Asian honeybee (Apis cerana) is an important indigenous species that play an indispensable role in the ecological balance and biological diversity. Few studies have been conducted to characterize the simple sequence repeats (SSRs) derived from A. cerana, so, in this study, a genome-wide screening for SSRs were firstly performed in the genome of A. cerana by comparison with that in west honeybee (Apis mellifera). There were 20,9991 SSRs distributed throughout the genome of A. cerana (Korea strain) and di-nucleotides were the most frequent SSR type. Both total number and density of SSRs in A. cerana genome were smaller than that in A. mellifera genome. Through comparing length discrepancy of SSRs loci among several isolates based on sequence alignment, 218 potential polymorphic SSRs primers derived from A. cerana were presented. Five among these SSR markers were evaluated for amplification in twenty-eight colonies of Apis cerana cerana (Chinese honeybee), which showed highly polymorphic, with the value of Polymorphism information content (PIC) ranging from 0.47 to 0.61. All these results will contribute to further develop more effective SSRs markers derived from A. cerana, which can be used to study genetic structure and population polymorphism of Asian honeybee.

The chlorophyll-deficient golden leaf mutation in cucumber is due to a single nucleotide substitution in CsChlI for magnesium chelatase I subunit

The cucumber chlorophyll-deficient golden leaf mutation is due to a single nucleotide substitution in the CsChlI gene for magnesium chelatase I subunit which plays important roles in the chlorophyll biosynthesis pathway. The Mg-chelatase catalyzes the insertion of Mg(2+) into the protoporphyrin IX in the chlorophyll biosynthesis pathway, which is a protein complex encompassing three subunits CHLI, CHLD, and CHLH. Chlorophyll-deficient mutations in genes encoding the three subunits have played important roles in understanding the structure, function and regulation of this important enzyme. In an EMS mutagenesis population, we identified a chlorophyll-deficient mutant C528 with golden leaf color throughout its development which was viable and able to set fruits and seeds. Segregation analysis in multiple populations indicated that this leaf color mutation was recessively inherited and the green color showed complete dominance over golden color. Map-based cloning identified CsChlI as the candidate gene for this mutation which encoded the CHLI subunit of cucumber Mg-chelatase. The 1757-bp CsChlI gene had three exons and a single nucleotide change (G to A) in its third exon resulted in an amino acid substitution (G269R) and the golden leaf color in C528. This mutation occurred in the highly conserved nucleotide-binding domain of the CHLI protein in which chlorophyll-deficient mutations have been frequently identified. The mutant phenotype, CsChlI expression pattern and the mutated residue in the CHLI protein suggested the mutant allele in C528 is unique among mutations identified so far in different species. This golden leaf mutant not only has its potential in cucumber breeding, but also provides a useful tool in understanding the CHLI function and its regulation in the chlorophyll biosynthesis pathway as well as chloroplast development.

Fine Mapping of Virescent Leaf Gene v-1 in Cucumber (Cucumis sativus L.)

Leaf color mutants are common in higher plants that can be used as markers in crop breeding or as an important tool in understanding regulatory mechanisms in chlorophyll biosynthesis and chloroplast development. In virescent leaf mutants, young leaves are yellow in color, which gradually return to normal green when the seedlings grow large. In the present study, we conducted phenotypic characterization and genetic mapping of the cucumber virescent leaf mutant 9110Gt conferred by the v-1 locus. Total chlorophyll and carotenoid content in 9110Gt was reduced by 44% and 21%, respectively, as compared with its wild type parental line 9110G. Electron microscopic investigation revealed fewer chloroplasts per cell and thylakoids per chloroplast in 9110Gt than in 9110G. Fine genetic mapping allowed for the assignment of the v-1 locus to a 50.4 kb genomic DNA region in chromosome 6 with two flanking markers that were 0.14 and 0.16 cM away from v-1, respectively. Multiple lines of evidence supported CsaCNGCs as the only candidate gene for the v-1 locus, which encoded a cyclic-nucleotide-gated ion channel protein. A single nucleotide change in the promoter region of v-1 seemed to be associated with the virescent color change in 9110Gt. Real-time PCR revealed significantly lower expression of CsaCNGCs in the true leaves of 9110Gt than in 9110G. This was the first report that connected the CsaCNGCs gene to virescent leaf color change, which provided a useful tool to establish linkages among virescent leaf color change, chloroplast development, chlorophyll biosynthesis, and the functions of the CsaCNGCs gene.

A Genome-Wide Survey of the Microsatellite Content of the Globe Artichoke Genome and the Development of a Web-Based Database

The recently acquired genome sequence of globe artichoke (Cynara cardunculus var. scolymus) has been used to catalog the genome’s content of simple sequence repeat (SSR) markers. More than 177,000 perfect SSRs were revealed, equivalent to an overall density across the genome of 244.5 SSRs/Mbp, but some 224,000 imperfect SSRs were also identified. About 21% of these SSRs were complex (two stretches of repeats separated by <100 nt). Some 73% of the SSRs were composed of dinucleotide motifs. The SSRs were categorized for the numbers of repeats present, their overall length and were allocated to their linkage group. A total of 4,761 perfect and 6,583 imperfect SSRs were present in 3,781 genes (14.11% of the total), corresponding to an overall density across the gene space of 32,5 and 44,9 SSRs/Mbp for perfect and imperfect motifs, respectively. A putative function has been assigned, using the gene ontology approach, to the set of genes harboring at least one SSR. The same search parameters were applied to reveal the SSR content of 14 other plant species for which genome sequence is available. Certain species-specific SSR motifs were identified, along with a hexa-nucleotide motif shared only with the other two Compositae species (sunflower (Helianthus annuus) and horseweed (Conyza canadensis)) included in the study. Finally, a database, called “Cynara cardunculus MicroSatellite DataBase” (CyMSatDB) was developed to provide a searchable interface to the SSR data. CyMSatDB facilitates the retrieval of SSR markers, as well as suggested forward and reverse primers, on the basis of genomic location, genomic vs genic context, perfect vs imperfect repeat, motif type, motif sequence and repeat number. The SSR markers were validated via an in silico based PCR analysis adopting two available assembled transcriptomes, derived from contrasting globe artichoke accessions, as templates.