To study the SSR loci information and develop molecular markers, a total of 435,858 unigenes in transcriptome of Polygonatum sibiricum were used to explore SSR. The distribution frequency of SSR and the basic characteristics of repeat motifs were analyzed using MISA software, and SSR primers were designed by Primer 3.0 software and then validated by PCR. Moreover, the gene function analysis of SSR Unigene was obtained by Blast. The results showed that 112,728 SSR loci were found in the transcriptome of Polygonatum sibiricum, which distributed in 435,858 unigenes with a distribution frequency of 25.86%. Mo‐nucleotide and Di‐nucleotide repeat were the main types, accounted for 83.83% of all SSRs. The repeat motifs of A/T and AC/GT were the predominant repeat types of Mo‐nucleotide and Di‐nucleotide, respectively. A total of 113,305 pairs of SSR primers with the potential to produce polymorphism were designed for maker development. One hundred and fifty‐four of the 500 randomly selected primers not only produced fragments with expected molecular size but also had high polymorphism, which could accurately separate the tested varieties. The gene function of unigenes containing SSR was mostly related to the molecular function of Polygonatum sibiricum. The SSR markers in transcriptome of Polygonatum sibiricum show rich type, strong specificity, and high potential of polymorphism, which will benefit the candidate gene mining and marker‐assisted breeding. The developed markers can also provide technical methods for molecular identification of intraspecific species of Polygonatum Mill. and maker‐assisted breeding of superior varieties of Polygonatum Mill.

Black barley is caused by phytomelanin synthesized in lemma and/or pericarp and the trait is controlled by one dominant gene Blp1. The gene is mapped on chromosome 1H by molecular markers, but it is yet to be isolated. Specific-locus amplified fragment sequencing (SLAF-seq) is an effective method for large-scale de novo single nucleotide polymorphism (SNP) discovery and genotyping. In the present study, SLAF-seq with bulked segregant analysis (BSA) was employed to obtain sufficient markers to fine mapping Blp1 gene in an F2 population derived from Hatiexi No.1 × Zhe5819. Based on SNP screening criteria, a total of 77,542 polymorphic SNPs met the requirements for association analysis. Combining two association analysis methods, the overlapped region with a size of 32.41 Mb on chromosome 1H was obtained as the candidate region of Blp1 gene. According to SLAF-seq data, markers were developed in the target region and were used for mapping the Blp1 gene. Linkage analysis showed that Blp1 co-segregated with HZSNP34 and HZSNP36, and was delimited by two markers (HZSNP35 and HZSNP39) spanning 8.1 cM in 172 homozygous yellow grain F2 plants of Hatiexi No.1 × Zhe5819. More polymorphic markers were screened in the reduced target region and were used to genotype the population. As a result, Blp1 was delimited within a 1.66 Mb on chromosome 1H by the upstream marker HZSNP63 and the downstream marker HZSNP59. Our results demonstrated the utility of SLAF-seq-BSA approach to identify the candidate region and discover polymorphic markers at the specific targeted genomic region.

Flavonoids, a class of phenolic compounds, are one of the main functional components and have a wide range of molecular structures and biological activities in Polygonatum. A few of them, including homoisoflavonoids, chalcones, isoflavones, and flavones, were identified in Polygonatum and displayed a wide range of powerful biological activities, such as anti-cancer, anti-viral, and blood sugar regulation. However, few studies have systematically been published on the flavonoid biosynthesis pathway in Polygonatum cyrtonema Hua. Therefore, in the present study, a combined transcriptome and metabolome analysis was performed on the leaf, stem, rhizome, and root tissues of P. cyrtonema to uncover the synthesis pathway of flavonoids and to identify key regulatory genes. Flavonoid-targeted metabolomics detected a total of 65 active substances from four different tissues, among which 49 substances were first study to identify in Polygonatum, and 38 substances were flavonoids. A total of 19 differentially accumulated metabolites (DAMs) (five flavonols, three flavones, two dihydrochalcones, two flavanones, one flavanol, five phenylpropanoids, and one coumarin) were finally screened by KEGG enrichment analysis. Transcriptome analysis indicated that a total of 222 unigenes encoding 28 enzymes were annotated into three flavonoid biosynthesis pathways, which were “phenylpropanoid biosynthesis”, “flavonoid biosynthesis”, and “flavone and flavonol biosynthesis”. The combined analysis of the metabolome and transcriptome revealed that 37 differentially expressed genes (DEGs) encoding 11 enzymes (C4H, PAL, 4CL, CHS, CHI, F3H, DFR, LAR, ANR, FNS, FLS) and 19 DAMs were more likely to be regulated in the flavonoid biosynthesis pathway. The expression of 11 DEGs was validated by qRT-PCR, resulting in good agreement with the RNA-Seq. Our studies provide a theoretical basis for further elucidating the flavonoid biosynthesis pathway in Polygonatum.

Barley semi-dwarf genes have been extensively explored and widely used in barley breeding programs. The semi-dwarf gene ari-e from Golden Promise is an important gene associated with some agronomic traits and salt tolerance. While ari-e has been mapped on barley chromosome 5H using traditional markers and next-generation sequencing technologies, it has not yet been finely located on this chromosome. We integrated two methods to develop molecular markers for fine-mapping the semi-dwarf gene ari-e: (1) specific-length amplified fragment sequencing (SLAF-seq) with bulked segregant analysis (BSA) to develop SNP markers, and (2) the whole-genome shotgun sequence to develop InDels. Both SNP and InDel markers were developed in the target region and used for fine-mapping the ari-e gene. Linkage analysis showed that ari-e co-segregated with marker InDel-17 and was delimited by two markers (InDel-16 and DGSNP21) spanning 6.8 cM in the doubled haploid (DH) Dash × VB9104 population. The genetic position of ari-e was further confirmed in the Hindmarsh × W1 DH population which was located between InDel-7 and InDel-17. As a result, the overlapping region of the two mapping populations flanked by InDel-16 and InDel-17 was defined as the candidate region spanning 0.58 Mb on the POPSEQ physical map. The current study demonstrated the SLAF-seq for SNP discovery and whole-genome shotgun sequencing for InDel development as an efficient approach to map complex genomic region for isolation of functional gene. The ari-e gene was fine mapped from 10 Mb to 0.58 Mb interval.