Psathyrostachys juncea is a perennial forage grass which plays an important role in soil and water conservation and ecological maintenance in cold and dry areas of temperate regions. In P. juncea, a variety of biotic and abiotic stress related genes have been used in crop improvement, indicating its agronomic, economic, forage, and breeding value. To date, there have been few studies on the genetic structure of P. juncea. Here, the genetic diversity and population structure of P. juncea were analyzed by EST-SSR molecular markers to evaluate the genetic differentiation related to tillering traits in P. juncea germplasm resources. The results showed that 400 simple sequence repeat (SSR) loci were detected in 2,020 differentially expressed tillering related genes. A total of 344 scored bands were amplified using 103 primer pairs, out of which 308 (89.53%) were polymorphic. The Nei’s gene diversity of 480 individuals was between 0.092 and 0.449, and the genetic similarity coefficient was between 0.5008 and 0.9111, with an average of 0.6618. Analysis of molecular variance analysis showed that 93% of the variance was due to differences within the population, and the remaining 7% was due to differences among populations. Psathyrostachys juncea materials were clustered into five groups based on population genetic structure, principal coordinate analysis and unweighted pair-group method with arithmetic means (UPGMA) analysis. The results were similar between clustering methods, but a few individual plants were distributed differently by the three models. The clustering results, gene diversity and genetic similarity coefficients showed that the overall genetic relationship of P. juncea individuals was relatively close. A Mantel test, UPGMA and structural analysis also showed a significant correlation between genetic relationship and geographical distribution. These results provide references for future breeding programs, genetic improvement and core germplasm collection of P. juncea.

Background Tillering is a complicated process in plant and is a significant trait that affects biomass and seed yield of bunch grass Psathyrostachys juncea, a typical perennial forage species. To clarify the regulatory mechanisms of tillering in P. juncea and to explore related candidate genes could be helpful to improve the seed and forage yield of perennial gramineous forages. We selected the tiller node tissues of P. juncea for transcriptome sequencing to determine the differentially expressed genes (DEG) between dense and sparse tillering genotypes. The metabolic pathway was studied, candidate genes were screened, and reference genes stability were evaluated. Results The results showed that approximately 5466 DEGs were identified between the two genotypes with dense and sparse tillers of P. juncea, which significantly differed in tiller number. Tillering regulation pathways analysis suggested that DEGs closely related to the biosynthesis of three plant hormones, namely auxin (IAA), cytokinin (CTK), and strigolactones (SLs), while “biosynthesis of lignin” and “nitrogen metabolism” have remarkable differences between the dense and sparse tillering genotypes. Meanwhile, the reference gene Actin1, having the best stability, was screened from twelve genes with highest expression level and was used in verification of ten tillering related candidate genes. Conclusions The tillering mechanism of perennial grass P. juncea was expounded by transcriptome analysis of tiller node tissues. We demonstrated that dense-tillering genotypes may be distinguished by their low expression patterns of genes involved in SL, IAA, and high expression patterns of genes involved in CTK biosynthesis at the tillering stage, and nitrogen metabolism and lignin biosynthesis can also affect the number of tillers. Furthermore, the expression level of ten tillering related candidate genes were verified using Actin1 as reference gene. These candidate genes provide valuable breeding resources for marker assisted selection and yield traits improvement of P. juncea.