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

Genome-wide development of simple sequence repeats markers and genetic diversity analysis of chayote

BACKGROUND

Chayote is a high economic crop in the Cucurbitaceae family, playing an important role in food production, disease treatment and the production of degradable materials in industries. Due to the harsh environment, such as high temperature, drought and frost, some chayote resources are gradually disappearing. It is crucial to collect, characterize, and conserve chayote resources. However, the genetic diversity of chayote resources in China has not been studied so far.

RESULTS

In this study, we collected 35 individuals of chayote from 14 provinces in China. Subsequently, we found 363,156 SSR motifs from the chayote genome and designed 57 pairs of SSR primers for validation. Out of these, 48 primer pairs successfully amplified bands, with 42 of them showing polymorphism. These 42 primer pairs detected a total of 153 alleles, averaging 3.64 alleles per locus. The polymorphic information content ranged from 0.03 to 0.78, with an average value of 0.41, indicating a high level of polymorphism. Based on the analysis using STRUCTURE, PCoA, and UPGMA methods, the 35 chayote individuals were divided into two major clusters. Through further association analysis, 7 significantly associated SSR markers were identified, including four related to peel color and three related to spine.

CONCLUSIONS

These molecular markers will contribute to the analysis of genetic diversity and genetic breeding improvement of chayote in the future.

Effects of trimer repeats on Psidium guajava L. gene expression and prospection of functional microsatellite markers

Most research on trinucleotide repeats (TRs) focuses on human diseases, with few on the impact of TR expansions on plant gene expression. This work investigates TRs' effect on global gene expression in Psidium guajava L., a plant species with widespread distribution and significant relevance in the food, pharmacology, and economics sectors. We analyzed TR-containing coding sequences in 1,107 transcripts from 2,256 genes across root, shoot, young leaf, old leaf, and flower bud tissues of the Brazilian guava cultivars Cortibel RM and Paluma. Structural analysis revealed TR sequences with small repeat numbers (5-9) starting with cytosine or guanine or containing these bases. Functional annotation indicated TR-containing genes' involvement in cellular structures and processes (especially cell membranes and signal recognition), stress response, and resistance. Gene expression analysis showed significant variation, with a subset of highly expressed genes in both cultivars. Differential expression highlighted numerous down-regulated genes in Cortibel RM tissues, but not in Paluma, suggesting interplay between tissues and cultivars. Among 72 differentially expressed genes with TRs, 24 form miRNAs, 13 encode transcription factors, and 11 are associated with transposable elements. In addition, a set of 20 SSR-annotated, transcribed, and differentially expressed genes with TRs was selected as phenotypic markers for Psidium guajava and, potentially for closely related species as well.

Genome-wide identification of microsatellites for mapping, genetic diversity and cross-transferability in wheat (Triticum spp)

Wheat (Triticum aestivum L.) is a crucial global staple crop, and is consistently being improved to enhance yield, disease resistance, and quality traits. However, the development of molecular markers is a challenging task due to its hexaploid genome. Molecular marker system such as simple sequence repeat (SSR) and single nucleotide polymorphism (SNP) are helpful for breeding, but SNP has limitations due to its development cost and its conversion to breeder markers. The study proposed an in-silico approach, by utilizing the low-cost transcriptome sequencing of two parental lines, 'TAC 75' and 'WH 1105', to identify polymorphic SSRs for mapping in a recombinant inbred line (RIL) population. This study introduces a new approach to bridge wheat genetics intricacies and next-generation sequencing potential. It presents a comprehensive genome-wide SSR distribution using IWGSC CS RefSeq v2.1 genome assembly and to identify 189 polymorphic loci through in-silico strategy. Of these, 54.76% showed polymorphism between parents, surpassing the traditional low polymorphic success rate. A RIL population screening validated these markers, demonstrating the fitness of identified markers through chi-square tests. The designed SSRs were also validated for genetic diversity analysis in a subset of 37 Indian wheat genotypes and cross-transferability in the wild/relative wheat species. In diversity analysis, a subset of 38 markers revealed 95 alleles (2.5 allele/locus), indicating substantial genetic variation. Population structure analysis unveiled three distinct groups, supported by phylogenetic and PCoA analyses. Further the polymorphic SSRs were also analyzed for SSR-gene association using gene ontology analysis. By utilizing the developing seed transcriptome data within parental lines, the study has enhanced the polymorphic SSR identification precision and facilitated in the RIL population. The undertaken study pioneers the use of transcriptome sequencing and genetic mapping to overcome challenges posed by the intricate wheat genome. This approach offers a cost-effective, less labour-intensive alternative to conventional methods, providing a platform for advancing wheat breeding research.

Genome Assembly and Microsatellite Marker Development Using Illumina and PacBio Sequencing in the Carex pumila (Cyperaceae) from Korea

This study is the first to report the characterization of Carex pumila genomic information. Assembly of the genome generated a draft of C. pumila based on PacBio Sequel II and Illumina paired-end sequencing, which was assembled from 2941 contigs with an estimated genome size of 0.346 Gb. The estimate of repeats in the genome was 31.0%, and heterozygosity ranged from 0.426 to 0.441%. The integrity evaluation of the assembly revealed 1481 complete benchmarked universal single-copy orthologs (BUSCO) (91.76%), indicating the high quality of the draft assembly. A total of 23,402 protein-coding genes were successfully predicted and annotated in the protein database. UpsetR plots showed that 7481 orthogroups were shared by all species. The phylogenetic tree showed that C. pumila is a close but distant relative of Ananas comosus. C. pumila had greater contraction (3154) than expansion (392). Among the extended gene families, aquaporins have been found to be enriched. Primers for microsatellite markers determined 30 polymorphic markers out of 100. The average number of alleles amplified by these 30 polymorphic markers was 4 to 12, with an average polymorphism information content (PIC) value of 0.660. In conclusion, our study provides a useful resource for comparative genomics, phylogeny, and future population studies of C. pumila.

Genome-wide survey and characterization of microsatellites in cashew and design of a web-based microsatellite database: CMDB

The cashew is an edible tree nut crop having a wide range of food and industrial applications. Despite great economic importance, the genome-wide characterization of microsatellites [simple sequence repeats (SSRs)] in cashew is lacking. In this study, we carried out the first comprehensive genome-wide microsatellites/SSRs characterization in cashew and developed polymorphic markers and a web-based microsatellite database. A total of 54526 SSRs were discovered in the cashew genome, with a mean frequency of 153 SSRs/Mb. Among the mined genome-wide SSRs (2-6 bp size motifs), the dinucleotide repeat motifs were dominant (68.98%) followed by the trinucleotides (24.56%). The Class I type of SSRs (≥20 bp) were 45.10%, while Class II repeat motifs (≥12-<20 bp) were 54.89% of the total genomic SSRs discovered here. Further, the AT-rich SSRs occurred more frequently in the cashew genome (84%) compared to the GC-rich SSRs. The validation of the in silico-mined genome-wide SSRs by PCR screening in cashew genotypes resulted in the development of 59 polymorphic SSR markers, and the polymorphism information content (PIC) of the polymorphic SSR markers ranged from 0.19 to 0.84. Further, a web-based database, "Cashew Microsatellite Database (CMDB)," was constructed to provide access to the genome-wide SSRs mined in this study as well as transcriptome-based SSRs from our previous study to the research community through a user-friendly searchable interface. Besides, CMDB provides information on experimentally validated SSRs. CMDB permits the retrieval of SSR markers information with the customized search options. Altogether, the genome-wide SSRs characterization, the polymorphic markers and CMDB database developed in this study would serve as valuable marker resources for DNA fingerprinting, germplasm characterization, genetic studies, and molecular breeding in cashew and related Anacardium species.

Comprehensive genome-wide identification and transferability of chromosome-specific highly variable microsatellite markers from citrus species

Citrus species among the most important and widely consumed fruit in the world due to Vitamin C, essential oil glands, and flavonoids. Highly variable simple sequence repeats (SSR) markers are one of the most informative and versatile molecular markers used in perennial tree genetic research. SSR survey of Citrus sinensis and Citrus maxima were identified perfect SSRs spanning nine chromosomes. Furthermore, we categorized all SSR motifs into three major classes based on their tract lengths. We designed and validated a class I SSRs in the C. sinensis and C. maxima genome through electronic polymerase chain reaction (ePCR) and found 83.89% in C. sinensis and 78.52% in C. maxima SSRs producing a single amplicon. Then, we selected extremely variable SSRs (> 40 nt) from the ePCR-verified class I SSRs and in silico validated across seven draft genomes of citrus, which provided us a subset of 84.74% in C. sinensis and 77.53% in C. maxima highly polymorphic SSRs. Out of these, 129 primers were validated on 24 citrus genotypes through wet-lab experiment. We found 127 (98.45%) polymorphic HvSSRs on 24 genotypes. The utility of the developed HvSSRs was demonstrated by analysing genetic diversity of 181 citrus genotypes using 17 HvSSRs spanning nine citrus chromosomes and were divided into 11 main groups through 17 HvSSRs. These chromosome-specific SSRs will serve as a powerful genomic tool used for future QTL mapping, molecular breeding, investigation of population genetic diversity, comparative mapping, and evolutionary studies among citrus and other relative genera/species.

Small cardamom genome: development and utilization of microsatellite markers from a draft genome sequence of Elettaria cardamomum Maton

Small cardamom (Elettaria cardamomum Maton), the queen of spices, is the third most expensive spice in the world after saffron and vanilla, valued highly for its aroma and taste. This perennial herbaceous plant is a native of coastal parts of Southern India and displays a significant amount of morphological diversity. Its genetic potential has not been exploited due to lack of genomic resources limiting our understanding of the genome and important metabolic pathways which give it the economic advantage in the spice industry. Here, we report upon the de novo assembled, draft whole genome sequence of cardamom variety, Njallani Green Gold. We used a hybrid assembly strategy using the reads from the Oxford Nanopore, Illumina and 10x Genomics GemCode sequencing chemistries. The assembled genome length was 1.06 Gb (gigabases) which is close to the estimated genome size of cardamom. More than 75% of the genome was captured in 8000 scaffolds with a N50 of 0.15 Mb. The genome appears to have a high repeat content and 68055 gene models were predicted. The genome is close to Musa species and displays an expansion and contraction in different gene families. The draft assembly was used for in silico mining of simple sequence repeats (SSRs). A total of 2,50,571 SSRs were identified of which 2,18,270 were perfect SSRs and 32,301 were compound SSRs. Among the perfect SSRs, trinucleotides were most abundant (1,25,329) and hexanucleotide repeats appear least (2,380). From the 2,50,571 SSRs mined, 2,27,808 primer pairs were designed based on flanking sequence information. Wet lab validation was performed for 246 SSR loci and based on their amplification profiles, 60 SSR markers were used for diversity analysis of a set of 60 diverse cardamom accessions. The average number of alleles detected per locus were 14.57 with a minimum of 4 and maximum of 30 alleles. Population structure analysis revealed the presence of high degree of admixtures which could primarily be due to cross-pollination prevalent in this species. The SSR markers identified would help in the development of gene or trait-linked markers which can be subsequently used for marker-assisted breeding for crop improvement in cardamom. The information on utilization of the SSR loci for generation of markers has been developed into a public database, 'cardamomSSRdb' that is freely available for use by the cardamom community.

Genome-wide characterization and development of SSR markers for genetic diversity analysis in northwestern Himalayas Walnut (Juglans regia L.)

In the present study, we designed and validated genome-wide polymorphic SSR markers (110 SSRs) by mining the walnut genome. A total of 198,924 SSR loci were identified. Among these, successful primers were designed for 162,594 (81.73%) SSR loci. Dinucleotides were the most predominant accounting for 88.40% (175,075) of total SSRs. The SSR frequency was 377.312 SSR/Mb and it showed a decreasing trend from dinucleotide to octanucleotide motifs. We identified 20 highly polymorphic SSR markers and used them to genotype 72 walnut accessions. Over all, we obtained 118 alleles that ranged from 2 to 12 with an average value of 5.9. The higher SSR PIC values indicate their robustness in discriminating walnut genotypes. Heat map, PCA, and population structure categorized 72 walnut genotypes into 2 distinct clusters. The genetic variation within population was higher than among population as inferred by analysis of molecular variance (AMOVA). For walnut improvement, it is necessary to have a large repository of SSRs with high discriminative power. The present study reports 150,000 SSRs, which is the largest SSR repository for this important nut crop. Scientific communities may use this repository for walnut improvement such as QTL mapping, genetic studies, linkage map construction, and marker-assisted selection.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03563-6.

First Report on Development of Genome-Wide Microsatellite Markers for Stock (Matthiola incana L.)

Stock (Matthiola incana (L.) R. Br.) is a famous annual ornamental plant with important ornamental and economic value. The lack of DNA molecular markers has limited genetic analysis, genome evolution, and marker-assisted selective breeding studies of M. incana. Therefore, more DNA markers are needed to support the further elucidation of the biology and genetics of M. incana. In this study, a high-quality genome of M. incana was initially assembled and a set of effective SSR primers was developed at the whole-genome level using genome data. A total of 45,612 loci of SSRs were identified; the di-nucleotide motifs were the most abundant (77.35%). In total, 43,540 primer pairs were designed, of which 300 were randomly selected for PCR validation, and as the success rate for amplification. In addition, 22 polymorphic SSR markers were used to analyze the genetic diversity of 40 stock varieties. Clustering analysis showed that all varieties could be divided into two clusters with a genetic distance of 0.68, which were highly consistent with their flower shape (potted or cut type). Moreover, we have verified that these SSR markers are effective and transferable within the Brassicaceae family. In this study, potential SSR molecular markers were successfully developed for 40 M. incana varieties using whole genome analysis, providing an important genetic tool for theoretical and applied research on M. incana.

Genome-wide identification and expression analysis of calmodulin-like proteins in cucumber

Background

The calmodulin-like (CML) protein is a crucial Ca²⁺-binding protein that can sense and conduct the Ca²⁺ signal in response to extracellular stimuli. The CML protein families have been identified and characterized in many species. Nevertheless, scarce information on cucumber CML is retrievable.

Methods

In this study, bioinformatic analyses, including gene structure, conserved domain, phylogenetic relationship, chromosome distribution, and gene synteny, were comprehensively performed to identify and characterize CsCML gene members. Spatiotemporal expression analysis in different organs and environment conditions were assayed with real-time quantitative polymerase chain reaction (qRT-PCR).

Results

Forty-four CsCMLs family members were well characterized, and the results showed that the 44 CsCML proteins contained one to four EF-hand domains without other functional domains. Most of the CsCML proteins were intron-less and unevenly distributed on seven chromosomes; two tandemly duplicated gene pairs and three segmentally duplicated gene pairs were identified in the cucumber genome. Cis-acting element analysis showed that the hormone, stress, and plant growth and development-related elements were in the promotor regions. In addition, spatiotemporal expression analysis revealed distinctive expression patterns for CsCML genes in different tissues and environmental conditions, and a putative protein interaction network also confirmed their potential role in responding to various stimuli. These results provide a foundation for understanding CsCMLs and provide a theoretical basis for further study of the physiological functions of CsCMLs.

Fine mapping of CscpFtsY, a gene conferring the yellow leaf phenotype in cucumber (Cucumis sativus L.)

BACKGROUND

Leaf color mutants are ideal materials to study pigment metabolism and photosynthesis. Leaf color variations are mainly affected by chlorophylls (Chls) and carotenoid contents and chloroplast development in higher plants. However, the regulation of chlorophyll metabolism remains poorly understood in many plant species. The chloroplast signal-recognition particle system is responsible for the insertion of the light-harvesting chlorophyll a/b proteins (LHCPs) to thylakoid membranes, which controls the chloroplast development as well as the regulation of Chls biosynthesis post-translationally in higher plants.

RESULTS

In this study, the yellow leaf cucumber mutant, named yl, was found in an EMS-induced mutant library, which exhibited a significantly reduced chlorophyll content, abnormal chloroplast ultrastructure and decreased photosynthetic capacity. Genetic analysis demonstrated that the phenotype of yl was controlled by a recessive nuclear gene. Using BSA-seq technology combined with the map-based cloning method, we narrowed the locus to a 100 kb interval in chromosome 3. Linkage analysis and allelism test validated the candidate SNP residing in CsaV3_3G009150 encoding one homolog of chloroplast signal-recognition particle (cpSRP) receptor in Arabidopsis, cpFtsY, could be responsible for the yellow leaf phenotype of yl. The relative expression of CscpFtsY was significantly down-regulated in different organs except for the stem, of yl compared with that in the wild type (WT). Subcellular localization result showed that CscpFtsY located in the chloroplasts of mesophyll cells.

CONCLUSIONS

The yl mutant displayed Chls-deficient, impaired chloroplast ultrastructure with intermittent grana stacks and significantly decreased photosynthetic capacity. The isolation of CscpFtsY in cucumber could accelerate the progress on chloroplast development by cpSRP-dependant LHCP delivery system and regulation of Chls biosynthesis in a post-translational way.

Genome-Wide Survey and Analysis of Microsatellites in Waterlily, and Potential for Polymorphic Marker Development

Waterlily (Nymphaeaceae), a diploid dicotyledon, is an ornamental aquatic plant. In 2020, the complete draft genome for the blue-petal waterlily (Nymphaea colorata) was made available in GenBank. To date, the genome-wide mining of microsatellites or simple sequence repeats (SSRs) in waterlily is still absent. In the present study, we investigated the characteristics of genome-wide microsatellites for N. colorata and developed polymorphic SSR markers across tropical and hardy waterlilies. A total of 238,816 SSRs were identified in 14 N. colorata chromosomes with an average density of 662.60 SSRs per Mb, and the largest number of SSRs were present on chromosome 1 (n = 30,426, 705.94 SSRs per Mb). The dinucleotide was the most common type, and AT-rich repeats prevail in the N. colorata genome. The SSR occurrence frequencies decreased as the number of motif repeats increased. Among 2442 protein-coding region SSRs, trinucleotides, accounting for 63.84%, were the most abundant. Gene ontology terms for signal transduction (e.g., GO: 0045859 and GO: 0019887) and the lipoic acid metabolism (ko00785,) were overrepresented in GO and KEGG enrichment analysis, respectively. In addition, 107,152 primer pairs were identified, and 13 novel polymorphism SSR markers were employed to distinguish among nine waterlily cultivars, of which Ny-5.2 and Ny-10.1 were the most informative SSR loci. This study contributes the first detailed characterization of SSRs in N. colorata genomes and delivers 13 novel polymorphism markers, which are useful for the molecular breeding strategies, genetic diversity and population structure analysis of waterlily.

RsmA3 modulates RpoS through the RetS-Gac-Rsm signaling pathway in response to H2 O2 stress in the phytopathogen Pseudomonas syringae

Reactive oxygen species are a fatal challenge to the plant pathogenic bacterium Pseudomonas syringae. In this study, we reveal that the global regulatory protein RsmA3 from the RetS-Gac/Rsm signaling pathway modulates RpoS in the early-log growth phase in the P. syringae wild-type strain MB03, thereby regulating oxidative tolerance to H₂ O₂ and ultimately affecting pathogenicity to the host plant. Following increased H₂ O₂ by external addition or endogenous induction by menadione, the resistance of the mutant strain ΔretS to H₂ O₂ is significantly enhanced due to rapid increases in the transcription of Rsm-related non-coding small RNAs (nc sRNAs), a sigma factor RpoS, and H₂ O₂ -detoxifying enzymes. Moreover, the ΔretS mutant is significantly less pathogenic in cucumber leaves. Seven Rsm-related nc sRNAs (namely, rsmZ, rsmY, and rsmX_1-5 ) show functional redundancy in the RetS-Gac-Rsm signaling pathway. External addition of H₂ O₂ stimulates increases in the transcription of both rsmY and rsmZ. Thus, we propose a regulatory model of the RetS-Gac-Rsm signaling pathway in P. syringae MB03 for the regulation of H₂ O₂ tolerance and phytopathogenicity in the host plant.

Serial-Omics and Molecular Function Study Provide Novel Insight into Cucumber Variety Improvement

Cucumbers are rich in vitamins and minerals. The cucumber has recently become one of China’s main vegetable crops. More specifically, the adjustment of the Chinese agricultural industry’s structure and rapid economic development have resulted in increases in the planting area allocated to Chinese cucumber varieties and in the number of Chinese cucumber varieties. After complete sequencing of the “Chinese long” genome, the transcriptome, proteome, and metabolome were obtained. Cucumber has a small genome and short growing cycle, and these traits are conducive to the application of molecular breeding techniques for improving fruit quality. Here, we review the developments and applications of molecular markers and genetic maps for cucumber breeding and introduce the functions of gene families from the perspective of genomics, including fruit development and quality, hormone response, resistance to abiotic stress, epitomizing the development of other omics, and relationships among functions.

A set of sampling, preparation, and staining techniques for studying meiosis in cucumber

The development of genetic and genomic resources for biological studies in cucumber has experienced an unprecedented boom in recent years. To investigate the function of putative meiotic genes and germplasm in breeding programs, an accurate cytogenetic characterization is required. Cytological methods and reference to investigate meiosis in cucumber are limited at present. Here we provide a set of cytological techniques that have been adapted for the study of meiosis in cucumber. The meiotic stages can be identified with high precision using hierarchical criteria from developing buds, undisturbed meiocytes, and freshly stained chromosomes. A meiotic cytological atlas of all stages is presented as a reference for identifying particular stages and for comparison of meiosis between normal and mutant plants. We performed a comparative analysis of the distribution of cytoplasmic organelles between cucumber and Arabidopsis, and we described a highly nonsynchronous condensation of chromosome parts during diplotene. A simplified fluorescence in situ hybridization (FISH) protocol, using robustly spread chromosomes, were developed. In addition, we designed a single oligonucleotide probe for 5S rDNA to use in karyotyping and monitoring of homologous chromosome pairing, which will make FISH analysis of 5S rDNA easier and more economical.

Genome-Wide Survey and Development of the First Microsatellite Markers Database (AnCorDB) in Anemone coronaria L

Anemone coronaria L. (2n = 2x = 16) is a perennial, allogamous, highly heterozygous plant marketed as a cut flower or in gardens. Due to its large genome size, limited efforts have been made in order to develop species-specific molecular markers. We obtained the first draft genome of the species by Illumina sequencing an androgenetic haploid plant of the commercial line “MISTRAL^® Magenta”. The genome assembly was obtained by applying the MEGAHIT pipeline and consisted of 2 × 10⁶ scaffolds. The SciRoKo SSR (Simple Sequence Repeats)-search module identified 401.822 perfect and 188.987 imperfect microsatellites motifs. Following, we developed a user-friendly “Anemone coronaria Microsatellite DataBase” (AnCorDB), which incorporates the Primer3 script, making it possible to design couples of primers for downstream application of the identified SSR markers. Eight genotypes belonging to eight cultivars were used to validate 62 SSRs and a subset of markers was applied for fingerprinting each cultivar, as well as to assess their intra-cultivar variability. The newly developed microsatellite markers will find application in Breeding Rights disputes, developing genetic maps, marker assisted breeding (MAS) strategies, as well as phylogenetic studies.

Genome-wide SSR markers in bottle gourd: development, characterization, utilization in assessment of genetic diversity of National Genebank of India and synteny with other related cucurbits

Lagenaria siceraria (Molina) Standley is an important cultivated crop with its immense importance in pharmaceutical industry and as vegetable. Its seed, root, stem, leaves, flower, and fruit are used as an ointment for ailment of various diseases throughout Asia. Despite its worldwide importance, informative co-dominant microsatellite markers in the bottle gourd crop are very restricted, impeding genetic improvement, cultivar identification, and phylogenetic studies. Next-generation sequencing has revolutionized the approaches for discovery, assessment, and validation of molecular markers. We conducted a genome-wide analysis, for developing SSR markers by utilizing restriction site-associated DNA sequencing (RAD-Seq) data obtained from NCBI. By performing in silico mining of microsatellite repeat motifs, we developed 45,066 perfect SSR markers. Of which 207 markers were successfully validated and 120 (57.97%) polymorphic primer pairs were utilized for an in-depth genetic diversity and population structure analysis of 96 accessions from the National Genebank of India. Tetranucleotide repeats (∼34.3%) were the most prevalent followed by trinucleotide repeats (∼30.73%), further 21.03%, 9.6%, and 4.3% of di-, penta-, and hexa-nucleotide repeats in the bottle gourd genome, respectively. Synteny of SSR markers on 11 bottle gourd linkage groups was correlated with the 7 chromosomes of cucumber (93.2%), 12 chromosomes of melon (87.4%), and 11 of watermelon (90.8%). The generated SSR markers provide a valuable tool for germplasm characterization, genetic linkage map construction, studying synteny, gene discovery, and for breeding in bottle gourd and other cucurbits species. KEY MESSAGE: Development of 45,066 perfect microsatellite markers as a valuable tool for marker assisted selection (MAS) in plant breeding.

Fruit transcriptional profiling of the contrasting genotypes for shelf life reveals the key candidate genes and molecular pathways regulating post-harvest biology in cucumber

Cucumber fruits are perishable in nature and become unfit for market within 2-3 days of harvesting. A natural variant, DC-48 with exceptionally high shelf life was developed and used to dissect the genetic architecture and molecular mechanism for extended shelf life through RNA-seq for first time. A total of 1364 DEGs were identified and cell wall degradation, chlorophyll and ethylene metabolism related genes played key role. Polygalacturunase (PG), Expansin (EXP) and xyloglucan were down regulated determining fruit firmness and retention of fresh green colour was mainly attributed to the low expression level of the chlorophyll catalytic enzymes (CCEs). Gene regulatory networks revealed the hub genes and cross-talk associated with wide variety of the biological processes. Large number of SSRs (21524), SNPs (545173) and InDels (126252) identified will be instrumental in cucumber improvement. A web genomic resource, CsExSLDb developed will provide a platform for future investigation on cucumber post-harvest biology.

BigFiRSt: A Software Program Using Big Data Technique for Mining Simple Sequence Repeats From Large-Scale Sequencing Data

Background

Simple Sequence Repeats (SSRs) are short tandem repeats of nucleotide sequences. It has been shown that SSRs are associated with human diseases and are of medical relevance. Accordingly, a variety of computational methods have been proposed to mine SSRs from genomes. Conventional methods rely on a high-quality complete genome to identify SSRs. However, the sequenced genome often misses several highly repetitive regions. Moreover, many non-model species have no entire genomes. With the recent advances of next-generation sequencing (NGS) techniques, large-scale sequence reads for any species can be rapidly generated using NGS. In this context, a number of methods have been proposed to identify thousands of SSR loci within large amounts of reads for non-model species. While the most commonly used NGS platforms (e.g., Illumina platform) on the market generally provide short paired-end reads, merging overlapping paired-end reads has become a common way prior to the identification of SSR loci. This has posed a big data analysis challenge for traditional stand-alone tools to merge short read pairs and identify SSRs from large-scale data.

Results

In this study, we present a new Hadoop-based software program, termed BigFiRSt, to address this problem using cutting-edge big data technology. BigFiRSt consists of two major modules, BigFLASH and BigPERF, implemented based on two state-of-the-art stand-alone tools, FLASH and PERF, respectively. BigFLASH and BigPERF address the problem of merging short read pairs and mining SSRs in the big data manner, respectively. Comprehensive benchmarking experiments show that BigFiRSt can dramatically reduce the execution times of fast read pairs merging and SSRs mining from very large-scale DNA sequence data.

Conclusions

The excellent performance of BigFiRSt mainly resorts to the Big Data Hadoop technology to merge read pairs and mine SSRs in parallel and distributed computing on clusters. We anticipate BigFiRSt will be a valuable tool in the coming biological Big Data era.

Genome-wide identification and expression analysis of the cucumber PYL gene family

Abscisic acid (ABA) is a very important hormone in plants. It regulates growth and development of plants and plays an important role in biotic and abiotic stresses. The Pyrabactin resistance 1-like (PYR/PYL) proteins play a central role in ABA signal transduction pathways. The working system of PYL genes in cucumber, an important economical vegetable (Cucumis sativus L.), has not been fully studied yet. Through bioinformatics, a total of 14 individual PYL genes were identified in Chinese long '9930' cucumber. Fourteen PYL genes were distributed on six chromosomes of cucumber, and their encoded proteins predicted to be distributed in cytoplasm and nucleus. Based on the phylogenetic analysis, the PYL genes of cucumber, Arabidopsis, rice, apple, Brachypodium distachyon and soybeancould be classified into three groups. Genetic structures and conserved domains analysis revealed that CsPYL genes in the same group have similar exons and conserved domains. By predicting cis-elements in the promoters, we found that all CsPYL members contained hormone and stress-related elements. Additionally, the expression patterns of CsPYL genes were specific in tissues. Finally, we further examined the expression of 14 CsPYL genes under ABA, PEG, salt stress. The qRT-PCR results showed that most PYL gene expression levels were up-regulated. Furthermore, with different treatments about 3h, the relative expression of PYL8 was up-regulated and more than 20 times higher than 0h. It indicated that this gene may play an important role in abiotic stress.

Integrating physiology, genetics, and transcriptome to decipher a new thermo-sensitive and light-sensitive virescent leaf gene mutant in cucumber

Chloroplasts are the material basis of photosynthesis, and temperature and light severely affect chloroplast development and thus influence photosynthetic efficiency. This study identified a spontaneous virescent leaf mutant, SC311Y, whose cotyledons and true leaves were yellow and gradually turned green. However, temperature and light affected the process of turning green. In addition, this mutant (except at the seedling stage) had ruffled leaves with white stripes, sterile males, and poorly fertile female flowers. Genetic characteristics analysis revealed that the recessive gene controlled the virescent leaf. Two F₂ populations mapped v-3 to the interval of 33.54-35.66 Mb on chromosome 3. In this interval, BSA-Seq, RNA-Seq, and cDNA sequence analyses revealed only one nonsynonymous mutation in the Csa3G042730 gene, which encoded the RNA exosome supercomplex subunit resurrection1 (RST1). Csa3G042730 was predicted to be the candidate gene controlling the virescent leaf, and the candidate gene may regulate chloroplast development by regulating plastid division2 (PDV2). A transcriptome analysis showed that different factors caused the reduced chlorophyll and carotenoid content in the mutants. To our knowledge, this study is the first report of map-based cloning related to virescent leaf, male-sterile, and chloroplast RNA regulation in cucumber. The results could accelerate the study of the RNA exosome supercomplex for the dynamic regulation of chloroplast RNA.

Genome-Wide Analysis of Simple Sequence Repeats in Cabbage (Brassica oleracea L.)

Cabbage (Brassica oleracea L. var. capitata) accounts for a critical vegetable crop belonging to Brassicaceae family, and it has been extensively planted worldwide. Simple sequence repeats (SSRs), the markers with high polymorphism and co-dominance degrees, offer a crucial genetic research resource. The current work identified totally 64,546 perfect and 93,724 imperfect SSR motifs in the genome of the cabbage 'TO1000.' Then, we divided SSRs based on the respective overall length and repeat number into different linkage groups. Later, we characterized cabbage genomes from the perspectives of motif length, motif-type classified and SSR level, and compared them across cruciferous genomes. Furthermore, a large set of 64,546 primer pairs were successfully identified, which generated altogether 1,113 SSR primers, including 916 (82.3%) exhibiting repeated and stable amplification. In addition, there were 32 informative SSR markers screened, which might decide 32 cabbage genotypes for their genetic diversity, with level of polymorphism information of 0.14-0.88. Cultivars were efficiently identified by the new strategy designating manual diagram for identifying cultivars. Lastly, 32 cabbage accessions were clearly separately by five Bol-SSR markers. Besides, we verified whether such SSRs were available and transferable in 10 Brassicaceae relatives. Based on the above findings, those genomic SSR markers identified in the present work may facilitate cabbage research, which lay a certain foundation for further gene tagging and genetic linkage analyses, like marker-assisted selection, genetic mapping, as well as comparative genomic analysis.

Comparative genomic analysis of expansin superfamily gene members in zucchini and cucumber and their expression profiles under different abiotic stresses

Zucchini and cucumber belong to the Cucurbitaceae family, a group of economical and nutritious food plants that is consumed worldwide. Expansin superfamily proteins are generally localized in the cell wall of plants and are known to possess an effect on cell wall modification by causing the expansion of this region. Although the whole genome sequences of cucumber and zucchini plants have been resolved, the determination and characterization of expansin superfamily members in these plants using whole genomic data have not been implemented yet. In the current study, a genome-wide analysis of zucchini (Cucurbita pepo) and cucumber (Cucumis sativus) genomes was performed to determine the expansin superfamily genes. In total, 49 and 41 expansin genes were identified in zucchini and cucumber genomes, respectively. All expansin superfamily members were subjected to further bioinformatics analysis including gene and protein structure, ontology of the proteins, phylogenetic relations and conserved motifs, orthologous relations with other plants, targeting miRNAs of those genes and in silico gene expression profiles. In addition, various abiotic stress responses of zucchini and cucumber expansin genes were examined to determine their roles in stress tolerance. CsEXPB-04 and CsEXPA-11 from cucumber and CpEXPA-20 and CpEXPLA-14 from zucchini can be candidate genes for abiotic stress response and tolerance in addition to their roles in the normal developmental processes, which are supported by the gene expression analysis. This work can provide new perspectives for the roles of expansin superfamily genes and offers comprehensive knowledge for future studies investigating the modes of action of expansin proteins.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01108-w.

Molecular research progress and improvement approach of fruit quality traits in cucumber

Recent molecular studies revealed new opportunities to improve cucumber fruit quality. However, the fruit color and spine traits molecular basis remain vague despite the vast sources of genetic diversity. Cucumber is agriculturally, economically and nutritionally important vegetable crop. China produces three-fourths of the world's total cucumber production. Cucumber fruit quality depends on a number of traits such as the fruit color (peel and flesh color), spine (density, size and color), fruit shape, fruit size, defects, texture, firmness, taste, maturity stage and nutritional composition. Fruit color and spine traits determine critical quality attributes and have been the interest of researchers at the molecular level. Evaluating the molecular mechanisms of fruit quality traits is important to improve production and quality of cucumber varieties. Genes and qualitative trait locus (QTL) that are responsible for cucumber fruit color and fruit spine have been identified. The purpose of this paper is to reveal the molecular research progress of fruit color and spines as key quality traits of cucumber. The markers and genes identified so far could help for marker-assisted selection of the fruit color and spine trait in cucumber breeding and its associated nutritional improvement. Based on the previous studies, peel color and spine density as examples, we proposed a comprehensive approach for cucumber fruit quality traits improvement. Moreover, the markers and genes can be useful to facilitate cloning-mediated genetic breeding in cucumber. However, in the era of climate change, increased human population and high-quality demand of consumers, studies on molecular mechanisms of cucumber fruit quality traits are limited.

Chromosome-specific polymorphic SSR markers in tropical eucalypt species using low coverage whole genome sequences: systematic characterization and validation

Eucalyptus is one of the major plantation species with wide variety of industrial uses. Polymorphic and informative simple sequence repeats (SSRs) have broad range of applications in genetic analysis. In this study, two individuals of Eucalyptus tereticornis (ET217 and ET86), one individual each from E. camaldulensis (EC17) and E. grandis (EG9) were subjected to whole genome resequencing. Low coverage (10×) genome sequencing was used to find polymorphic SSRs between the individuals. Average number of SSR loci identified was 95,513 and the density of SSRs per Mb was from 157.39 in EG9 to 155.08 in EC17. Among all the SSRs detected, the most abundant repeat motifs were di-nucleotide (59.6%–62.5%), followed by tri- (23.7%–27.2%), tetra- (5.2%–5.6%), penta- (5.0%–5.3%), and hexa-nucleotide (2.7%–2.9%). The predominant SSR motif units were AG/CT and AAG/TTC. Computational genome analysis predicted the SSR length variations between the individuals and identified the gene functions of SSR containing sequences. Selected subset of polymorphic markers was validated in a full-sib family of eucalypts. Additionally, genome-wide characterization of single nucleotide polymorphisms, InDels and transcriptional regulators were carried out. These variations will find their utility in genome-wide association studies as well as understanding of molecular mechanisms involved in key economic traits. The genomic resources generated in this study would provide an impetus to integrate genomics in marker-trait associations and breeding of tropical eucalypts.

A mutation in CsHY2 encoding a phytochromobilin (PΦB) synthase leads to an elongated hypocotyl 1(elh1) phenotype in cucumber (Cucumis sativus L.)

The elongated hypocotyl1 (elh1) mutant in cucumber is due to a mutation in CsHY2, which is a homolog of the Arabidopsis HY2 encoding the phytochromobilin (PΦB) synthase for phytochrome biosynthesis Hypocotyl length is a critical determinant in establishing high quality seedlings for successful cucumber production, but knowledge on the molecular regulation of hypocotyl growth in cucumber is very limited. Here, we reported identification and characterization of a cucumber elongated hypocotyl 1 (elh1) mutant. We found that the longer hypocotyl in elh1 was due to longitudinal growth of hypocotyl cells. With fine mapping, the elh1 locus was delimited to a 20.9-kb region containing three annotated genes; only one polymorphism was identified in this region between two parental lines, which was a non-synonymous SNP (G28153633A) in the third exon of CsHY2 (CsGy1G030000) that encodes a phytochromobilin (PΦB) synthase. Uniqueness of the mutant allele at CsHY2 was verified in natural cucumber populations. Ectopic expression of CsHY2 in Arabidopsis hy2-1 long-hypocotyl mutant led to reduced hypocotyl length. The PΦB protein was targeted to the chloroplast. The expression levels of CsHY2 and five phytochrome genes CsPHYA1, CsPHYA2, CsPHYB, CsPHYC and CsPHYE were all significantly down-regulated while several cell elongation related genes were up-regulated in elh1 mutant compared to wild-type cucumber, which are correlated with dynamic hypocotyl elongation in the mutant. RNA-seq analysis in the WT and mutant revealed differentially expressed genes involved in porphyrin and chlorophyll metabolisms, cell elongation and plant hormone signal transduction pathways. This is the first report to characterize and clone the CsHY2 gene in cucumber. This work reveals the important of CsHY2 in regulating hypocotyl length and extends our understanding of the roles of CsHY2 in cucumber.

Identification of a putative candidate gene encoding 7-dehydrocholesterol reductase involved in brassinosteroids biosynthesis for compact plant architecture in Cucumber (Cucumis sativus L.)

By the strategy of bulked segregant analysis sequencing combined with genetic mapping, CsDWF5, which encodes 7 dehydrocholesterol reductase that involved in brassinosteroids biosynthesis, was identified as the candidate gene for cpa. Dwarf architecture is one of the most important breeding goals in crops. The biosynthesis and signal transduction of brassinosteroids (BRs) have a great impact on plant growth and development including plant architecture. Here, we identified a compact plant architecture (cpa) mutant from an EMS-induced cucumber population. cpa displayed the extremely dwarf phenotype with shortened internode and petiole, darkened and wrinkled leaf. Genetic analysis revealed that cpa was caused by a single recessive gene. By the strategy of bulked segregant analysis sequencing combined with genetic mapping, CsDWF5, encoding a 7-dehydrocholesterol reductase that involved in sterol biosynthesis, was identified as the candidate gene for cpa. One single nucleotide mutation (G→A) in splicing site causing 3-bp insertion (TAG) was found in the first base of the sixth intron of CsDWF5 in cpa, which furtherly resulted in the frameshift mutation and got a premature stop codon. The expression of CsDWF5 gene was significantly down regulated in different tissues of the cpa mutant compared with that in wild type. The phenotype of cpa could be partially recovered by exogenous BR treatment. Transcriptome analysis identified 1096 genes that exhibited differential expression between the cpa mutant and wild type. KEGG enrichment analysis indicated that differentially expressed genes were significantly enriched in BR biosynthesis and plant-pathogen interaction pathways. These results provide perspectives on the molecular mechanisms underlying the dwarfing phenotype in cucumber.

Genome Wide Characterization, Comparative and Genetic Diversity Analysis of Simple Sequence Repeats in Cucurbita Species

Simple sequence repeats (SSRs) are widely used in mapping constructions and comparative and genetic diversity analyses. Here, 103,056 SSR loci were found in Cucurbita species by in silico PCR. In general, the frequency of these SSRs decreased with the increase in the motif length, and di-nucleotide motifs were the most common type. For the same repeat types, the SSR frequency decreased sharply with the increase in the repeat number. The majority of the SSR loci were suitable for marker development (84.75% in Cucurbita moschata, 94.53% in Cucurbita maxima, and 95.09% in Cucurbita pepo). Using these markers, the cross-species transferable SSR markers between C. pepo and other Cucurbitaceae species were developed, and the complicated mosaic relationships among them were analyzed. Especially, the main syntenic relationships between C. pepo and C. moschata or C. maxima indicated that the chromosomes in the Cucurbita genomes were highly conserved during evolution. Furthermore, 66 core SSR markers were selected to measure the genetic diversity in 61 C. pepo germplasms, and they were divided into two groups by structure and unweighted pair group method with arithmetic analysis. These results will promote the utilization of SSRs in basic and applied research of Cucurbita species.

Identification of quantitative trait loci for powdery mildew resistance in highly resistant cucumber (Cucumis sativus L.) using ddRAD-seq analysis

Powdery mildew, caused by Podosphaera xanthii (syn. Sphaerotheca fuliginea ex Fr. Poll.), is one of the most economically important foliar diseases in cucumber (Cucumis sativus L.). Cucumber parental line 'Kyuri Chukanbohon Nou 5 Go', developed from weedy cucumber line CS-PMR1, is highly resistant to powdery mildew and is promising breeding material. We performed quantitative trait locus (QTL) analysis using double-digest restriction-site-associated DNA sequencing (ddRAD-Seq) in a population from a cross between 'Kyuri Chukanbohon Nou 5 Go' and the Japanese native cultivar 'Kaga-aonaga-fushinari', which is susceptible to powdery mildew. The resistance of the population and its parents was evaluated using leaf disc assays and image analysis. We detected one major QTL on Chr. 5 that was effective at both 20°C and 25°C and one minor QTL on Chr. 1 effective at 20°C. We detected two additional QTLs in subpopulation: one on Chr. 3 effective at 20°C and one on Chr. 5 effective at both 20°C and 25°C in a position different from the major QTL. The resistance alleles at all four QTLs were contributed by 'Kyuri Chukanbohon Nou 5 Go'. The results of this study can be used to develop practical DNA markers tightly linked to genes for powdery mildew resistance.

A SNP of HD-ZIP I transcription factor leads to distortion of trichome morphology in cucumber (Cucumis sativus L.)

BACKGROUND

Trichomes are excellent model systems for the analysis of cell differentiation and play essential roles in plant protection. From cucumber inbred line 'WD1', we identified an EMS-induced trichome abnormally developing mutant, nps, which exhibited smaller, denser and no pyramid-shaped head trichomes.

RESULTS

Using F₂ and BC₁ populations constructed from a cross between nps and '9930', the genetic analysis showed that the nps trait is controlled by a single recessive nuclear gene. We identified CsNps by map-based cloning with 576 individuals of the F₂ population generated from the cross of nps and inbred line '9930'. The CsNps was located at a 13.4-kb genomic region on chromosome 3, which region contains three predicted genes. Sequence analysis showed that only one single nucleotide mutation (C → T) between 9930 and nps was found in the second exon of Csa3G748220, a plant-specific class I HD-Zip gene. The result of allelism test also indicated that nps is a novel allelic mutant of Mict (Micro-trichome). Thus, nps was renamed mict-L130F. By comparing the transcriptome of mict-L130F vs WD1 and 06-2 (mict) vs 06-1 (wildtype, near-isogenic line of 06-2), several potential target genes that may be related to trichome development were identified.

CONCLUSIONS

Our results demonstrate that Mict-L130F is involved in the morphogenesis of trichomes. Map-based cloning of the Mict-L130F gene could promote the study of trichome development in cucumber.

Mapping and identification of CsSh5.1, a gene encoding a xyloglucan galactosyltransferase required for hypocotyl elongation in cucumber (Cucumis sativus L.)

CsSh5.1, which controls hypocotyl elongation under high temperature conditions in cucumber, was mapped to a 57.1 kb region on chromosome 5 containing a candidate gene encoding a xyloglucan galactosyltransferase. Hypocotyl growth is a vital process in seedling establishment. Hypocotyl elongation after germination relies more on longitudinal cell elongation than cell division. Cell elongation is largely determined by the extensibility of the cell wall. Here, we identified a spontaneous mutant in cucumber (Cucumis sativus L.), sh5.1, which exhibits a temperature-insensitive short hypocotyl phenotype. Genetic analysis showed that the phenotype of sh5.1 was controlled by a recessive nuclear gene. CsSh5.1 was mapped to a 57.1 kb interval on chromosome 5, containing eight predicted genes. Sequencing analysis revealed that the Csa5G171710 is the candidate gene of CsSh5.1, which was further confirmed via co-segregation analysis and genomic DNA sequencing in natural cucumber variations. The result indicated that hypocotyl elongation might be controlled by this gene. CsSh5.1 encodes a xyloglucan galactosyltransferase that specifically adds galactose to xyloglucan and forms galactosylated xyloglucans, which determine the strength and extensibility of the cell walls. CsSh5.1 expression in wild-type (WT) hypocotyl was significantly higher than that in sh5.1 hypocotyl under high temperature, suggesting its important role in hypocotyl cell elongation under high temperature. The identification of CsSh5.1 is helpful for elucidating the function of xyloglucan galactosyltransferase in cell wall expansion and understanding the mechanism of hypocotyl elongation in cucumber.

Genome-wide identification of CLE gene family and their potential roles in bolting and fruit bearing in cucumber (Cucumis sativus L.)

BACKGROUND

Signal peptides are essential for plant growth and development. In plants, biological processes including cell-cell communication, cellular proliferation and differentiation, cellular determination of self-incompatibility, and defensive responses, all depend heavily on peptide-signaling networks such as CLE (CLAVATA3/Embryo surrounding region-related). The CLEs are indispensable in different periods of plant growth and development, especially in maintaining the balance between proliferation and differentiation of stem cells in various meristematic tissues. The working system of CLE genes in cucumber, an important economical vegetable (Cucumis sativus L.), has not been fully studied yet. The distributional patterns of chromosome-level genome assembly in cucumber provide a fundamental basis for a genome-wide comparative analysis of CLE genes in such plants.

RESULTS

A total of 26 individual CLE genes were identified in Chinese long '9930' cucumber, the majority of which belong to unstable short alkaline and hydrophilic peptides. A comparative analysis showed a close relationship in the development of CLE genes among Arabidopsis thaliana, melon, and cucumber. Half of the exon-intron structures of all CsCLEs genes are single-exon genes, and motif 1, a typical CLE domain near the C-terminal functioning in signal pathways, is found in all cucumber CLE proteins but CsCLE9. The analysis of CREs (Cis-Regulatory Elements) in the upstream region of the 26 cucumber CLE genes indicates a possible relationship between CsCLE genes and certain functions of hormone response elements. Cucumber resulted closely related to Arabidopsis and melon, having seven and 15 orthologous CLE genes in Arabidopsis and melon, respectively. Additionally, the calculative analysis of a pair of orthologous genes in cucumber showed that as a part of the evolutionary process, CLE genes are undergoing a positive selection process which leads to functional differentiation. The specific expression of these genes was vigorous at the growth and development period and tissues. Cucumber gene CLV3 was overexpressed in Arabidopsis, more than half of the transformed plants in T1 generation showed the phenomena of obvious weakness of the development of growing point, no bolting, and a decreased ability of plant growth. Only two bolted strains showed that either the pod did not develop or the pod was short, and its development was significantly inferior to that in the wild type.

CONCLUSIONS

In this study, 26 CLE genes were identified in Chinese long '9930' cucumber genome. The CLE genes were mainly composed of alkaline hydrophilic unstable proteins. The genes of the CLE family were divided into seven classes, and shared close relationships with their homologs in Arabidopsis and melon. The specific expression of these genes was evaluated in different periods of growth and tissue development, and CLV3, which the representative gene of the family, was overexpressed in Arabidopsis, suggesting that it has a role in bolting and fruit bearing in cucumber.

Transcriptome sequencing and microsatellite marker discovery in Ailanthus altissima (Mill.) Swingle (Simaroubaceae)

Ailanthus altissima Swingle, is a tree species native to East Asia and has a great potential in decorative, bioenergy and industrial applications in many countries. To date, despite its commercial importance, the genomic and genetic resources available for this species are still insufficient. In this study, we characterized the transcriptome of A. altissima and developed thirteen EST-SSRs (expressed sequence tag-simple sequence repeats) based on Illumina paired-end RNA sequencing (RNA-seq). Besides, we developed ten polymorphic chloroplast microsatellite (cpSSR) markers using the available chloroplast genome of A. altissima. The transcriptome data produced 87,797 unigenes, of which 64,891 (73.91%) unigenes were successfully annotated in at least one protein database. For cpSSR markers the number of detected alleles (N) per marker varied from three at cpSSR12 to twelve at cpSSR8, the unbiased haploid diversity indices (uh) varied from 0.111 to 0.485, and haploid diversity indices (h) ranged from 0.101 to 0.444 with an average unbiased haploid diversity index (uh) of 0.274. Overall, a total of 65 different cpSSR alleles were identified at the ten loci among 165 individuals of A. altissima. The allele number per locus for EST-SSRs varied from 2.143 to 9.357, and the values of observed and expected heterozygosity ranged from 0.312 to 1.000 and 0.505 to 0.826, respectively. The molecular markers developed in this study will facilitate future genetic diversity, population structure, long distance-gene transfer and pollen-based gene flow analyses of A. altissima populations from its known distribution ranges in China focusing on planted and natural forest stands.

Genetic Diversity Assessment and Cultivar Identification of Cucumber (Cucumis sativus L.) Using the Fluidigm Single Nucleotide Polymorphism Assay

Genetic diversity analysis and cultivar identification were performed using a core set of single nucleotide polymorphisms (SNPs) in cucumber (Cucumis sativus L.). For the genetic diversity study, 280 cucumber accessions collected from four continents (Asia, Europe, America, and Africa) by the National Agrobiodiversity Center of the Rural Development Administration in South Korea and 20 Korean commercial F1 hybrids were genotyped using 151 Fluidigm SNP assay sets. The heterozygosity of the SNP loci per accession ranged from 4.76 to 82.76%, with an average of 32.1%. Population genetics analysis was performed using population structure analysis and hierarchical clustering (HC), which indicated that these accessions were classified mainly into four subpopulations or clusters according to their geographical origins. The subpopulations for Asian and European accessions were clearly distinguished from each other (F_ST value = 0.47), while the subpopulations for Korean F1 hybrids and Asian accessions were closely related (F_ST = 0.34). The highest differentiation was observed between American and European accessions (F_ST = 0.41). Nei's genetic distance among the 280 accessions was 0.414 on average. In addition, 95 commercial F1 hybrids of three cultivar groups (Baekdadagi-, Gasi-, and Nakhap-types) were genotyped using 82 Fluidigm SNP assay sets for cultivar identification. These 82 SNPs differentiated all cultivars, except seven. The heterozygosity of the SNP loci per cultivar ranged from 12.20 to 69.14%, with an average of 34.2%. Principal component analysis and HC demonstrated that most cultivars were clustered based on their cultivar groups. The Baekdadagi- and Gasi-types were clearly distinguished, while the Nakhap-type was closely related to the Baekdadagi-type. Our results obtained using core Fluidigm SNP assay sets provide useful information for germplasm assessment and cultivar identification, which are essential for breeding and intellectual right protection in cucumber.

Agarose-resolvable InDel markers based on whole genome re-sequencing in cucumber

Insertion and Deletion (InDel) are common features in genomes and are associated with genetic variation. The whole-genome re-sequencing data from two parents (X1 and X2) of the elite cucumber (Cucumis sativus) hybrid variety Lvmei No.1 was used for genome-wide InDel polymorphisms analysis. Obtained sequence reads were mapped to the genome reference sequence of Chinese fresh market type inbred line ‘9930’ and gaps conforming to InDel were pinpointed. Further, the level of cross-parents polymorphism among five pairs of cucumber breeding parents and their corresponding hybrid varieties were used for evaluating hybrid seeds purity test efficiency of InDel markers. A panel of 48 cucumber breeding lines was utilized for PCR amplification versatility and phylogenetic analysis of these markers. In total, 10,470 candidate InDel markers were identified for X1 and X2. Among these, 385 markers with more than 30 nucleotide difference were arbitrary chosen. These markers were selected for experimental resolvability through electrophoresis on an Agarose gel. Two hundred and eleven (211) accounting for 54.81% of markers could be validated as single and clear polymorphic pattern while 174 (45.19%) showed unclear or monomorphic genetic bands between X1 and X2. Cross-parents polymorphism evaluation recorded 68 (32.23%) of these markers, which were designated as cross-parents transferable (CPT) InDel markers. Interestingly, the marker InDel114 presented experimental transferability between cucumber and melon. A panel of 48 cucumber breeding lines including parents of Lvmei No. 1 subjected to PCR amplification versatility using CPT InDel markers successfully clustered them into fruit and common cucumber varieties based on phylogenetic analysis. It is worth noting that 16 of these markers were predominately associated to enzymatic activities in cucumber. These agarose-based InDel markers could constitute a valuable resource for hybrid seeds purity testing, germplasm classification and marker-assisted breeding in cucumber.

Inheritance of chloroplast and mitochondrial genomes in cucumber revealed by four reciprocal F1 hybrid combinations

Both genomes in chloroplasts and mitochondria of plant cell are usually inherited from maternal parent, with rare exceptions. To characterize the inheritance patterns of the organelle genomes in cucumber (Cucumis sativus var. sativus), two inbred lines and their reciprocal F₁ hybrids were analyzed using an next generation whole genome sequencing data. Their complete chloroplast genome sequences were de novo assembled, and a single SNP was identified between the parental lines. Two reciprocal F₁ hybrids have the same chloroplast genomes with their maternal parents. Meanwhile, 292 polymorphic sites were identified between mitochondrial genomes of the two parental lines, which showed the same genotypes with their paternal parents in the two reciprocal F₁ hybrids, without any recombination. The inheritance patterns of the chloroplast and mitochondria genomes were also confirmed in four additional cucumber accessions and their six reciprocal F₁ hybrids using molecular markers derived from the identified polymorphic sites. Taken together, our results indicate that the cucumber chloroplast genome is maternally inherited, as is typically observed in other plant species, whereas the large cucumber mitochondrial genome is paternally inherited. The combination of DNA markers derived from the chloroplast and mitochondrial genomes will provide a convenient system for purity test of F₁ hybrid seeds in cucumber breeding.

De novo transcriptomic analysis and identification of EST-SSR markers in Stephanandra incisa

Stephanandra incisa is a wild-type shrub with beautiful leaves and white flowers and is commonly used as a garden decoration accessory. However, the limited availability of genomic data of S. incisa has restricted its breeding process. Here, we identified EST-SSR markers using de novo transcriptome sequencing. In this study, a transcriptome database containing 35,251 unigenes, having an average length of 985 bp, was obtained from S. incisa. From these unigene sequences, we identified 5,555 EST-SSRs, with a distribution density of one SSR per 1.60 kb. Dinucleotides (52.96%) were the most detected SSRs, followed by trinucleotides (34.64%). From the EST-SSR loci, we randomly selected 100 sites for designing primer and used the DNA of 60 samples to verify the polymorphism. The average value of the effective number of alleles (Ne), Shannon's information index (I), and expective heterozygosity (He) was 1.969, 0.728, and 0.434, respectively. The polymorphism information content (PIC) value was in the range of 0.108 to 0.669, averaging 0.406, which represented a middle polymorphism level. Cluster analysis of S. incisa were also performed based on the obtained EST-SSR data in our work. As shown by structure analysis, 60 individuals could be classified into two groups. Thus, the identification of these novel EST-SSR markers provided valuable sequence information for analyzing the population structure, genetic diversity, and genetic resource assessment of S. incisa and other related species.

Comparative Genome-Wide Characterization of Microsatellites in Candida albicans and Candida dubliniensis Leading to the Development of Species-Specific Marker

BACKGROUND

Microsatellites or simple sequence repeats (SSR) are related to genomic structure, function, and certain diseases of taxonomically different organisms.

OBJECTIVE

To characterize microsatellites in two closely related Candida species by searching and comparing 1-6 bp nucleotide motifs and utilizing them to develop species-specific markers.

METHODS

Whole-genome sequence was downloaded from the public domain, microsatellites were mined and analyzed, and primers were synthesized.

RESULTS

A total of 15,821 and 7,868 microsatellites, with mono-nucleotides (8,679) and trinucleotides (3,156) as most frequent microsatellites, were mined in Candida dubliniensis and Candida albicans, respectively. Chromosome size was found positively correlated with microsatellite number in both the species, whereas it was negatively correlated with the relative abundance and density of microsatellites. A number of unique motifs were also found in both the species. Overall, microsatellite frequencies of each chromosome in C. dubliniensis were higher than in C. albicans.

CONCLUSION

The features of microsatellite distribution in the two species' genomes revealed that it is probably not conserved in the genus Candida. Data generated in this article could be used for comparative genome mapping and understanding the distribution of microsatellites and genome structure between these closely related and phenotypically misidentified species and may provide a foundation for the development of a new set of species-specific microsatellite markers. Here, we also report a novel microsatellite-based marker for C. dubliniensis-specific identification.

Comprehensive Characterization and Validation of Chromosome-Specific Highly Polymorphic SSR Markers From Pomegranate (Punica granatum L.) cv. Tunisia Genome

The simple sequence repeat (SSR) survey of 'Tunisia' genome (296.85 Mb) identified a total of 365,279 perfect SSRs spanning eight chromosomes, with a mean marker density of 1,230.6 SSRs/Mb. We found a positive trend in chromosome length and the SSR abundance as marker density enhanced with a shorter chromosome length. The highest number of SSRs (60,708) was mined from chromosome 1 (55.56 Mb), whereas the highest marker density (1,294.62 SSRs/Mb) was recorded for the shortest chromosome 8 (27.99 Mb). Furthermore, we categorized all SSR motifs into three major classes based on their tract lengths. Across the eight chromosomes, the class III had maximum number of SSR motifs (301,684, 82.59%), followed by the class II (31,056, 8.50%) and the class I (5,003, 1.37%). Examination of the distribution of SSR motif types within a chromosome suggested the abundance of hexanucleotide repeats in each chromosome followed by dinucleotides, and these results are consistent with 'Tunisia' genome features as a whole. Concerning major repeat types, AT/AG was the most frequent (14.16%), followed by AAAAAT/AAAAAG (7.89%), A/C (7.54%), AAT/AAG (5.23%), AAAT/AAAG (4.37%), and AAAAT/AAAAG (1.2%) types. We designed and validated a total of 3,839 class I SSRs in the 'Tunisia' genome through electronic polymerase chain reaction (ePCR) and found 1,165 (30.34%) SSRs producing a single amplicon. Then, we selected 906 highly variable SSRs (> 40 nt) from the ePCR-verified class I SSRs and in silico validated across multiple draft genomes of pomegranate, which provided us a subset of 265 highly polymorphic SSRs. Of these, 235 primers were validated on six pomegranate genotypes through wet-lab experiment. We found 221 (94%) polymorphic SSRs on six genotypes, and 187 of these SSRs had ≥ 0.5 PIC values. The utility of the developed SSRs was demonstrated by analyzing genetic diversity of 30 pomegranate genotypes using 16 HvSSRs spanning eight pomegranate chromosomes. In summary, we developed a comprehensive set of highly polymorphic genome-wide SSRs. These chromosome-specific SSRs will serve as a powerful genomic tool to leverage future genetic studies, germplasm management, and genomics-assisted breeding in pomegranate.

Fine Mapping and Transcriptome Analysis of Virescent Leaf Gene v-2 in Cucumber (Cucumis sativus L.)

Leaf color mutants are the ideal materials to explore the pathways of chlorophyll metabolism, chloroplast development and photosynthesis system. In this study, a new virescent leaf mutant 104Y was identified by spontaneous mutation, whose cotyledon and upper five true leaves were yellow color. The yellow true leaves gradually turned green from top to bottom with increased chlorophyll contents. Genetic analysis indicated that the virescent leaf was controlled by one single recessive gene v-2, which was accurately mapped into 36.0-39.7 Mb interval on chromosome 3 by using BSA-seq and linkage analysis. Fine mapping analysis further narrowed v-2 into 73-kb genomic region including eight genes with BC₁ and F₂ populations. Through BSA-seq and cDNA sequencing analysis, only one nonsynonymous mutation existed in the Csa3G890020 gene encoding auxin F-box protein was identified, which was predicted as the candidate gene controlling virescent leaf. Comparative transcriptome analysis and quantitative real-time PCR analysis revealed that the expression level of Csa3G890020 was not changed between EC1 and 104Y. However, RNA-seq analysis identified that the key genes involved in chlorophyll biosynthesis and auxin signaling transduction network were mainly down-regulated in 104Y compared with EC1, which indicated that the regulatory functions of Csa3G890020 could be performed at post-transcriptional level rather than transcriptional level. This is the first report to map-based clone an auxin F-box protein gene related to virescent leaf in cucumber. The results will exhibit a new insight into the chlorophyll biosynthesis regulated by auxin signaling transduction network.

The landscape of microsatellites in the enset (Ensete ventricosum) genome and web-based marker resource development

Ensete ventricosum (Musaceae, enset) is an Ethiopian food security crop. To realize the potential of enset for rural livelihoods, further knowledge of enset diversity, genetics and genomics is required to support breeding programs and conservation. This study was conducted to explore the enset genome to develop molecular markers, genomics resources, and characterize enset landraces while giving insight into the organization of the genome. We identified 233 microsatellites (simple sequence repeats, SSRs) per Mbp in the enset genome, representing 0.28% of the genome. Mono- and di-nucleotide repeats motifs were found in a higher proportion than other classes of SSR-motifs. In total, 154,586 non-redundant enset microsatellite markers (EMM) were identified and 40 selected for primer development. Marker validation by PCR and low-cost agarose gel electrophoresis revealed that 92.5% were polymorphic, showing a high PIC (Polymorphism Information Content; 0.87) and expected heterozygosity (He = 0.79-0.82). In silico analysis of genomes of closely related species showed 46.86% of the markers were transferable among enset species and 1.90% were transferable to Musa. The SSRs are robust (with basic PCR methods and agarose gel electrophoresis), informative, and applicable in measuring enset diversity, genotyping, selection and potentially breeding. Enset SSRs are available in a web-based database at https://enset-project.org/EnMom@base.html (or https://enset.aau.edu.et/index.html , downloadable from Figshare).

Genome-wide development of miRNA-based SSR markers in Cleistogenes songorica and analysis of their transferability to Gramineae/non-Gramineae species

Simple sequence repeat (SSR) markers are commonly used for many genetic applications, such as map construction, fingerprinting, and genetic diversity analyses, due to their high reproducibility, polymorphism, and abundance. Endogenous miRNAs play essential roles in plant development and gene expression under diverse biotic and abiotic stress conditions. In the present study, we predicted 110 miRNA-SSR primer pairs from 287 precursor miRNAs. Among 110 primer pairs, 85 were successfully amplified and examined for transferability to other Gramineae and non-Gramineae species. The results showed that all 82 primer pairs yielded unambiguous and strong amplification, and across the 23 studied Cleistogenes accessions, a total of 385 alleles were polymorphic. The number of alleles produced per primer varied from 3 to 11, with an average of 4.69 per locus. The expected heterozygosity (He) ranged from 0.44 to 0.88, with an average of 0.74 per locus, and the PIC (Polymorphism Information Content) values ranged from 0.34 to 0.87, with an average of 0.69 per locus. Furthermore, 1422 miRNA target genes were predicted and analyzed using the GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) databases. In conclusion, the results showed that an miRNA-based microsatellite marker system can be applicable for genetic diversity and marker-assisted breeding studies.

Development of a Large Gene-Associated SSR Marker Set and in-Depth Genetic Characterization in Scarlet Sage

Salvia splendens, scarlet or tropical sage, is a tender perennial herbaceous flowering plant popularly grown in public and private gardens all over the world. In this study, we developed a set of simple sequence repeats (SSRs) from genome-wide sequences to assess the genetic diversity and population structure among 112 cultivars. We obtained 364,379 SSRs by mining scarlet sage’s recently published whole genome sequence; 14,545 gene-associated SSR loci were identified in 2 kb gene flanking regions. Among the 768 gene-associated SSR primer sets we screened, 576 loci successfully amplified in DNA pools of 3–4 different cultivars, of which 271 remained polymorphic when tested across eight individual plants. We searched for the related gene functions attributable to these gene-associated SSRs using diverse databases, resulting in 259 Non-redundant matching sequences, 205 individual Gene Ontology (GO) terms, 236 assigned to eukaryotic orthologous groups, and 67 KEGG-annotated (Kyoto Encyclopedia of Genes and Genomes) sequences. We finally selected 41 polymorphic SSR loci to infer genetic diversity and population structure among 112 S. splendens accessions. Based on the developed gene-associated SSRs, clustering analyses consistently revealed two distinct genetic groups within the core collection of S. splendens cultivars. This work developed and characterized an exhaustive set of genome-wide gene-associated SSR markers for scarlet sage. These SSRs can provide species identification, genetic diversity and population structure information for S. splendens, and will therefore be important tools for the management and protection of S. splendens germplasm.

Genome-wide characterization of simple sequence repeats in Palmae genomes

Background

Microsatellites or simple sequence repeats (SSRs) have become the most significant DNA marker technology used in genetic research. The availability of complete draft genomes for a number of Palmae species has made it possible to perform genome-wide analysis of SSRs in these species. Palm trees are tropical and subtropical plants with agricultural and economic importance due to the nutritional value of their fruit cultivars.

Objective

This is the first comprehensive study examining and comparing microsatellites in completely-sequenced draft genomes of Palmae species.

Methods

We identified and compared perfect SSRs with 1–6 bp nucleotide motifs to characterize microsatellites in Palmae species using PERF v0.2.5. We analyzed their relative abundance, relative density, and GC content in five palm species: Phoenix dactylifera, Cocos nucifera, Calamus simplicifolius, Elaeis oleifera, and Elaeis guineensis.

Results

A total of 118241, 328189, 450753, 176608, and 70694 SSRs were identified, respectively. The six repeat types were not evenly distributed across the five genomes. Mono- and dinucleotide SSRs were the most abundant, and GC content was highest in tri- and hexanucleotide SSRs.

Conclusion

We envisage that this analysis would further substantiate more in-depth computational, biochemical, and molecular studies on the roles SSRs may play in the genome organization of the palm species. The current study contributes a detailed characterization of simple sequence repeats in palm genomes.

Electronic supplementary material

The online version of this article (10.1007/s13258-020-00924-w) contains supplementary material, which is available to authorized users.

A new SNP genotyping technology Target SNP-seq and its application in genetic analysis of cucumber varieties

To facilitate the utility of SNP-based genotyping, we developed a new method called target SNP-seq which combines the advantages of multiplex PCR amplification and high throughput sequencing. Compared with KASP, Microarrays, GBS and other SNP genotyping methods, target SNP-seq is flexible both in SNPs and samples, yields high accuracy, especially when genotyping genome wide perfect SNPs with high polymorphism and conserved flanking sequences, and is cost-effective, requiring 3 days and $7 for per DNA sample to genotype hundreds of SNP loci. The present study established a DNA fingerprint of 261 cucumber varieties by target SNP-seq with 163 perfect SNPs from 4,612,350 SNPs based on 182 cucumber resequencing datasets. Four distinct subpopulations were found in 261 Chinese cucumber varieties: the north China type, the south China type, the Europe type, and the Xishuangbanna type. The north China type and Xishuangbanna type harbored lower genetic diversity, indicating greater risk of genetic erosion in these two subpopulations. Furthermore, a core set of 24 SNPs was able to distinguish 99% of the 261 cucumber varieties. 29 core cucumber backbone varieties in China were identified. Therefore, target SNP-seq provides a new way to screen out core SNP loci from the whole genome for DNA fingerprinting of crop varieties. The high efficiency and low cost of target SNP-seq is more competitive than the current SNP genotyping methods, and it has excellent application prospects in genetic research, as well as in promoting plant breeding processes in the near future.

Development of polymorphic EST-SSR markers and their applicability in genetic diversity evaluation in Rhododendron arboreum

The genus Rhododendron, known for large impressive flowers is widely distributed throughout the world. Rhododendrons have limited genetic information, despite of comprising high species diversity, morphological overlap and weak genetic barrier. In present study, expressed sequence tag (EST) data from Rhododendron catawbiense Michx (Subgenus Hymenanthes, Section Ponticum) and Rhododendron mucronatum var. ripense (Makino) E.H. Wilson (Subgenus Tsutsusi, Section Tsutsusi) were utilized for mining and identification of the SSRs for genetic diversity analysis of R. arboreum Smith (Subgenus Tsutsusi, Section Tsutsusi). A total of 249 SSRs were developed from 1767 contigs. Di-nucleotide was found to be most abundant repeat followed by tri- and tetra-nucleotide repeats. The motif AG/CT was most common di-nucleotide motif (31.73%), whereas, AAC/GTT (8.43%), ACG/CGT (8.03%), AAG/CTT (7.23%) and AGG/CCT (6.43%) were most abundant tri-nucleotide repeat motif. Among these SSRs, 168 sequences were only fit into the criteria to design flanking primer pairs. A total of 30 randomly selected primer pairs were utilized for validation and genetic diversity study in 36 genotypes of R. arboreum collected from western Himalayan region. In aggregate, 26 SSR markers (86.66%) produced good and repeatable amplifications. Expected heterozygosity (H_E) ranged from 0.322 to 0.841 and observed heterozygosity (H_O) ranged from 0.327 to 1.000 and PIC value ranged from 0.008 to 0.786. These primers were able to distinguish the geographic differences of occurrence based on cluster analysis. These developed EST-SSRs can be useful in future population genetics analysis and micro-evolutionary studies in Rhododendron species.

An irregularly striped rind mutant reveals new insight into the function of PG1β in cucumber (Cucumis sativus L.)

Via bulked segregant analysis sequencing combined with linkage mapping, the ist gene responsible for the irregularly striped rind mutation was delimited to a 144-kb region in cucumber. Sequencing and expression analysis identified Csa1G005490 as the candidate gene. The rind appearance of cucumber is one of the most important commercial quality traits. Usually, an immature cucumber fruit has a uniform rind that varies from green to yellow to white among different cultivated varieties. In the present paper, we isolated a novel fruit appearance cucumber mutant, ist, that has an irregularly striped rind pattern. The mutant displayed green irregular stripes on a yellow-green background at the immature fruit stage. Genetic analysis revealed that a single recessive gene, ist, is responsible for this mutation. A BSA (bulked segregant analysis) sequencing approach combined with genetic mapping delimited the ist locus to an interval with a length of 144 kb, and 21 predicted genes were annotated in the region. Based on mutation site screening and expression analysis, two single-nucleotide polymorphisms within the candidate gene, Csa1G005490, were identified as constituting the mutation. Csa1G005490 encodes a polygalacturonase-1 noncatalytic subunit beta protein (PG1β) known to be involved in fruit softening. The expression of Csa1G005490 was significantly lower in the ist mutant than in the wild type. Transcriptome analysis identified 1796 differentially expressed genes (DEGs) between the ist mutant and wild type. Gene ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses indicated that these DEGs were enriched mostly in photosynthesis and chlorophyll metabolism pathways. Decreased expression patterns of several chlorophyll synthesis genes in the mutant suggest that ist plays a key role in chlorophyll biosynthesis. These results will provide new insight into the molecular mechanism underlying rind appearance polymorphisms in cucumber.

Detection and application of genome-wide variations in peach for association and genetic relationship analysis

BACKGROUND

Peach (Prunus persica L.) is a diploid species and model plant of the Rosaceae family. In the past decade, significant progress has been made in peach genetic research via DNA markers, but the number of these markers remains limited.

RESULTS

In this study, we performed a genome-wide DNA markers detection based on sequencing data of six distantly related peach accessions. A total of 650,693~1,053,547 single nucleotide polymorphisms (SNPs), 114,227~178,968 small insertion/deletions (InDels), 8386~12,298 structure variants (SVs), 2111~2581 copy number variants (CNVs) and 229,357~346,940 simple sequence repeats (SSRs) were detected and annotated. To demonstrate the application of DNA markers, 944 SNPs were filtered for association study of fruit ripening time and 15 highly polymorphic SSRs were selected to analyze the genetic relationship among 221 accessions.

CONCLUSIONS

The results showed that the use of high-throughput sequencing to develop DNA markers is fast and effective. Comprehensive identification of DNA markers, including SVs and SSRs, would be of benefit to genetic diversity evaluation, genetic mapping, and molecular breeding of peach.