Foxtail millet (Setaria italica (L.) Beauv.) is an important gramineous grain-food and forage crop. It is grown worldwide for human and livestock consumption. Its small genome and diploid nature have led to foxtail millet fast becoming a novel model for investigating plant architecture, drought tolerance and C₄ photosynthesis of grain and bioenergy crops. Therefore, cost-effective, reliable and highly polymorphic molecular markers covering the entire genome are required for diversity, mapping and functional genomics studies in this model species. A total of 5,020 highly repetitive microsatellite motifs were isolated from the released genome of the genotype 'Yugu1’ by sequence scanning. Based on sequence comparison between S. italica and S. viridis, a set of 788 SSR primer pairs were designed. Of these primers, 733 produced reproducible amplicons and were polymorphic among 28 Setaria genotypes selected from diverse geographical locations. The number of alleles detected by these SSR markers ranged from 2 to 16, with an average polymorphism information content of 0.67. The result obtained by neighbor-joining cluster analysis of 28 Setaria genotypes, based on Nei’s genetic distance of the SSR data, showed that these SSR markers are highly polymorphic and effective. A large set of highly polymorphic SSR markers were successfully and efficiently developed based on genomic sequence comparison between different genotypes of the genus Setaria. The large number of new SSR markers and their placement on the physical map represent a valuable resource for studying diversity, constructing genetic maps, functional gene mapping, QTL exploration and molecular breeding in foxtail millet and its closely related species.

Background/Objectives: Agriophyllum squarrosum (L.) Moq. (A. squarrosum), also known as sandrice, is an important medicinal plant widely distributed in dunes across all the deserts of China. Common garden trials have shown content variations in flavonoids among the ecotypes of sandrice, which correlated with temperature heterogeneity in situ. However, there have not been any environmental control experiments to further elucidate whether the accumulation of flavonoids was triggered by cold stress; Methods: This study conducted a four-day ambient 4 °C low-temperature treatment on three ecotypes along with an in situ annual mean temperature gradient (Dulan (DL), Aerxiang (AEX), and Dengkou (DK)); Results: Target metabolomics showed that 12 out of 14 flavonoids in sandrice were driven by cold stress. Among them, several flavonoids were significantly up-regulated, such as naringenin and naringenin chalcone in all three ecotypes; isorhamnetin, quercetin, dihydroquercetin, and kaempferol in DL and AEX; and astragalin in DK. They were accompanied by 19 structural genes of flavonoid synthesis and 33 transcription factors were markedly triggered by cold stress in sandrice. The upstream genes, AsqAEX006535−CHS, AsqAEX016074−C4H, and AsqAEX004011−4CL, were highly correlated with the enrichment of naringenin, which could be fine-tuned by AsqAEX015868−bHLH62, AsqAEX001711−MYB12, and AsqAEX002220−MYB1R1; Conclusions: This study sheds light on how desert plants like sandrice adapt to cold stress by relying on a unique flavonoid biosynthesis mechanism that regulating the accumulation of naringenin. It also supports the precise development of sandrice for the medicinal industry. Specifically, quercetin and isorhamnetin should be targeted for development in DL and AEX, while astragalin should be precisely developed in DK.

The climate in arid Central Asia (ACA) has changed rapidly in recent decades, but the ecological consequences of this are far from clear. To predict the impacts of climate change on ecosystem functioning, greater attention should be given to the relationships between leaf functional traits and environmental heterogeneity. As a dominant constructive shrub widely distributed in ACA, Reaumuria soongarica provided us with an ideal model to understand how leaf functional traits of desert ecosystems responded to the heterogeneous environments of ACA. Here, to determine the influences of genetic and ecological factors, we characterized species‐wide variations in leaf traits among 30 wild populations of R. soongarica and 16 populations grown in a common garden. We found that the leaf length, width, and leaf length to width ratio (L/W) of the northern lineage were significantly larger than those of other genetic lineages, and principal component analysis based on the in situ environmental factors distinguished the northern lineage from the other lineages studied. With increasing latitude, leaf length, width, and L/W in the wild populations increased significantly. Leaf length and L/W were negatively correlated with altitude, and first increased and then decreased with increasing mean annual temperature (MAT) and mean annual precipitation (MAP). Stepwise regression analyses further indicated that leaf length variation was mainly affected by latitude. However, leaf width was uncorrelated with altitude, MAT, or MAP. The common garden trial showed that leaf width variation among the eastern populations was caused by both local adaptation and phenotypic plasticity. Our findings suggest that R. soongarica preferentially changes leaf length to adjust leaf size to cope with environmental change. We also reveal phenotypic evidence for ecological speciation of R. soongarica. These results will help us better understand and predict the consequences of climate change for desert ecosystem functioning.

As one of the hotspot regions for sympatric speciation studies, Evolution Canyon (EC) became an ideal place for its high level of microclimatic divergence interslopes. In this study, to highlight the genetic mechanisms of sympatric speciation, phenotypic variation on flowering time and transcriptomic divergence were investigated between two ecotypes of Ricotia lunaria, which inhabit the opposite temperate and tropical slopes of EC I (Lower Nahal Oren, Mount Carmel, Israel) separated by 100 m at the bottom of the slopes. Growth chamber results showed that flowering time of the ecotype from south-facing slope population # 3 (SFS 3) was significantly 3 months ahead of the north-facing slope population # 5 (NFS 5). At the same floral development stage, transcriptome analysis showed that 1,064 unigenes were differentially expressed between the two ecotypes, which enriched in the four main pathways involved in abiotic and/or biotic stresses responses, including flavonoid biosynthesis, α-linolenic acid metabolism, plant–pathogen interaction and linoleic acid metabolism. Furthermore, based on Ka/Ks analysis, nine genes were suggested to be involved in the ecological divergence between the two ecotypes, whose homologs functioned in RNA editing, ABA signaling, photoprotective response, chloroplasts protein-conducting channel, and carbohydrate metabolism in Arabidopsis thaliana. Among them, four genes, namely, SPDS1, FCLY, Tic21 and BGLU25, also showed adaptive divergence between R. lunaria and A. thaliana, suggesting that these genes could play an important role in plant speciation, at least in Brassicaceae. Based on results of both the phenotype of flowering time and comparative transcriptome, we hypothesize that, after long-time local adaptations to their interslope microclimatic environments, the molecular functions of these nine genes could have been diverged between the two ecotypes. They might differentially regulate the expression of the downstream genes and pathways that are involved in the interslope abiotic stresses, which could further diverge the flowering time between the two ecotypes, and finally induce the reproductive isolation establishment by natural selection overruling interslope gene flow, promoting sympatric speciation.