SNP Discovery and Genetic Variation of Candidate Genes Relevant to Heat Tolerance and Agronomic Traits in Natural Populations of Sand Rice (Agriophyllum squarrosum)

Phenotypic variation in hypocotyl elongation among elite sand rice (Agriophyllum squarrosum) lines

Sand rice (Agriophyllum squarrosum), widely distributed in Central Arid Asia and prevalent in the sand dunes of northern China, presents a promising potential as a climate-resilient crop. The plasticity of hypocotyl growth is the key trait for sand rice to cope with wind erosion and sand burial, ensure seedling emergence, and determine plant architecture. In this study, we assessed the overall hypocotyl phenotype of six sand rice elite lines, which were collected from different regions of northern China, and selected by our group over past decade through common garden trials. Significant phenotypic variations were observed in thousand-seed weight (TSW), seedling emergence percentage, hypocotyl length and diameter, and seedling fresh weight among the lines. The elite line Aerxiang (AEX) exhibited excellent agronomic performance with superior and synchronous emergence, and high survival percentage, distinguishing itself as a prime candidate for further large-scale cultivation. Contrastingly, the lines from the arid regions showed markedly lower performance. Partial Least Squares Path Modeling (PLSPM) was used to assess the impact of seed provenance climate factors, including annual mean temperature (AMT) and annual mean precipitation (AMP), on trait variability among lines. The findings indicate a significant correlation between climate factors and hypocotyl length, highlighting the intricate adaptation of sand rice to local climate. The comprehensive understanding of the mechanisms behind phenotypic variations offers valuable insights for sand rice de novo domestication and innovative germplasm resources, and lays the foundation for ecological restoration in sandy areas.

Heat-induced proteomic changes in anthers of contrasting rice genotypes under variable stress regimes

Heat stress drastically affects anther tissues resulting in poor plant fertility, necessitating an urgent need to determine the key proteome regulation associated with mature anther in response to heat stress. We identified several genotype - specific protein alterations in rice anthers of Moroberekan (Japonica, heat sensitive), IR64 (Indica, moderately heat tolerant), and Nagina22 (Aus, heat tolerant) in the short-term (ST_HS; one cycle of 42°C, 4 hours before anthesis) and long-term (LT_HS; 6 cycles of 38°C, 6 hours before anthesis) heat stress. The proteins upregulated in long-term heat stress in Nagina22 were enriched in biological processes related to unfolded protein binding and carboxylic acid metabolism, including amino acid metabolism. In short-term heat stress, Nagina22 anthers were enriched in proteins associated with vitamin E biosynthesis and GTPase activator activity. In contrast, downregulated proteins were related to ribosomal proteins. The expression of different Hsp20 and DnaJ was genotype specific. Overall, the heat response in Nagina22 was associated with its capacity for adequate metabolic control and cellular homeostasis, which may be critical for its higher reproductive thermotolerance. This study improves our understanding of thermotolerance mechanisms in rice anthers during anthesis and lays a foundation for breeding thermotolerant varieties via molecular breeding.

Integrated metabolomics and transcriptomics insights on flavonoid biosynthesis of a medicinal functional forage, Agriophyllum squarrosum (L.), based on a common garden trial covering six ecotypes

Agriophyllum squarrosum (L.) Moq., well known as sandrice, is an important wild forage in sandy areas and a promising edible and medicinal resource plant with great domestication potential. Previous studies showed flavonoids are one of the most abundant medicinal ingredients in sandrice, whereby isorhamnetin and isorhamnetin-3-glycoside were the top two flavonols with multiple health benefits. However, the molecular regulatory mechanisms of flavonoids in sandrice remain largely unclear. Based on a common garden trial, in this study, an integrated transcriptomic and flavonoids-targeted metabolomic analysis was performed on the vegetative and reproductive periods of six sandrice ecotypes, whose original habitats covered a variety of environmental factor gradients. Multiple linear stepwise regression analysis unveiled that flavonoid accumulation in sandrice was positively correlated with temperature and UVB and negatively affected by precipitation and sunshine duration, respectively. Weighted co-expression network analysis (WGCNA) indicated the bHLH and MYB transcription factor (TF) families might play key roles in sandrice flavonoid biosynthesis regulation. A total of 22,778 differentially expressed genes (DEGs) were identified between ecotype DL and ecotype AEX, the two extremes in most environmental factors, whereby 85 DEGs could be related to known flavonoid biosynthesis pathway. A sandrice flavonoid biosynthesis network embracing the detected 23 flavonoids in this research was constructed. Gene families Plant flavonoid O-methyltransferase (AsPFOMT) and UDP-glucuronosyltransferase (AsUGT78D2) were identified and characterized on the transcriptional level and believed to be synthases of isorhamnetin and isorhamnetin-3-glycoside in sandrice, respectively. A trade-off between biosynthesis of rutin and isorhamnetin was found in the DL ecotype, which might be due to the metabolic flux redirection when facing environmental changes. This research provides valuable information for understanding flavonoid biosynthesis in sandrice at the molecular level and laid the foundation for precise development and utilization of this functional resource forage.

Precipitation and local environment shape the geographic variation of seed size across natural populations of sand rice (Agriophyllum squarrosum)

Sand rice (Agriophyllum squarrosum) is widely distributed on dunes in the Asian interior arid zone, and its large intraspecific trait variation makes it a very good model for investigating the ecological processes underlying its adaptation to the desert environment. In this study, seed size variation across 68 natural populations was used to establish geographic patterns and to quantify the effects of the climate, soil, and collection-year weather variables. The length of the seed major axis and thousand seed weight (TSW) both showed significant longitudinal patterns. Long-term climate variables accounted for most of the explained variances for seed major axis (57.20%) and TSW (91.54%). Specifically, annual precipitation and minimum monthly precipitation had the most significantly positive and negative effects, indicating that longitudinal clines are driven by a precipitation gradient across the species' distribution range. A substantial unique effect of soil variables (27.27%) was found for seed major axis variation, but only 3.64% of TSW variation was explained by soil variables. Two extreme groups were selected to evaluate the genetic and plastic effects on seed size in a common garden experiment. Large-seeded individuals were more competitive in semi-arid regions, and had stronger adaptive plasticity as well as better performance in early seedling establishment, and hence they have potential for use in future domestication projects.

Pre-mRNA alternative splicing as a modulator for heat stress response in plants

The molecular responses of plants to the important abiotic stress, heat stress (HS), have been extensively studied at the transcriptional level. Alternative splicing (AS) is a post-transcriptional regulatory process in which an intron-containing gene can generate more than one mRNA variant. The impact of HS on the pre-mRNA splicing process has been reported in various eukaryotes but seldom discussed in-depth, especially in plants. Here, we review AS regulation in response to HS in different plant species. We discuss potential molecular mechanisms controlling heat-inducible AS regulation in plants and hypothesize that AS regulation participates in heat-priming establishment and HS memory maintenance. We propose that the pre-mRNA splicing variation is an important regulator of plant HS responses (HSRs).

Rewilding crops for climate resilience: economic analysis and de novo domestication strategies

To match predicted population growth, annual food production should be doubled by 2050. This is not achievable by current agronomical and breeding practices, due to the impact of climate changes and associated abiotic stresses on agricultural production systems. Here, we analyze the impact of global climate trends on crop productivity and show that the overall loss in crop production from climate-driven abiotic stresses may exceed US$170 billion year-1 and represents a major threat to global food security. We also show that abiotic stress tolerance had been present in wild progenitors of modern crops but was lost during their domestication. We argue for a major shift in our paradigm of crop breeding, focusing on climate resilience, and call for a broader use of wild relatives as a major tool in this process. We argue that, while molecular tools are currently in place to harness the potential of climate-resilient genes present in wild relatives, the complex polygenic nature of tolerance traits remains a major bottleneck in this process. Future research efforts should be focused not only on finding appropriate wild relatives but also on development of efficient cell-based high-throughput phenotyping platforms allowing assessment of the in planta operation of key genes.

Organic acid metabolites involved in local adaptation to altitudinal gradient in Agriophyllum squarrosum, a desert medicinal plant

Agriophyllum squarrosum (L.) Moq., a pioneer plant endemic to the temperate deserts of Asia, could be domesticated into an ideal crop with outstanding ecological and medicinal characteristics. A previous study showed differential organic acid accumulation between two in situ altitudinal ecotypes. To verify whether this accumulation was determined by environmental or genetic factors, we conducted organic acid targeted metabolic profiling among 14 populations of A. squarrosum collected from regions with different altitudes based on a common garden experiment. Results showed that the most abundant organic acid in A. squarrosum was citric acid (96.03%, 2322.90 μg g^-1). Association analysis with in situ environmental variables showed that salicylic acid content was positively correlated with altitudinal gradient. Considering the enrichment of salicylic acid and protocatechualdehyde in high-altitude populations based on the common garden experiment, and the high expression of their biosynthesis relative genes (i.e., PAL and C4H) in the in situ high-altitude ecotype, we propose that organic acid accumulation could be involved in local adaptation to high altitudes. This study not only addresses the molecular basis of local adaptation involving the accumulation of organic acids in the desert plant A. squarrosum but also provides a method to screen wild germplasms to mitigate the impact of global climate change.

Variations in Flavonoid Metabolites Along Altitudinal Gradient in a Desert Medicinal Plant Agriophyllum squarrosum

Agriophyllum squarrosum (L.) Moq., a pioneer plant endemic to the temperate deserts of Asia, could be domesticated into an ideal crop with outstanding ecological and medicinal characteristics. A previous study showed differential flavonoid accumulation between two in situ altitudinal ecotypes. To verify whether this accumulation was determined by environmental or genetic factors, we conducted flavonoid-targeted metabolic profiling among 14 populations of A. squarrosum collected from regions with different altitudes based on a common garden experiment. Results showed that the most abundant flavonoid in A. squarrosum was isorhamnetin (48.40%, 557.45 μg/g), followed by quercetin (13.04%, 150.15 μg/g), tricin (11.17%, 128.70 μg/g), isoquercitrin (7.59%, 87.42 μg/g), isovitexin (7.20%, 82.94 μg/g), and rutin (7.00%, 80.62 μg/g). However, based on a common garden at middle-altitude environment, almost none of the flavonoids was enriched in the high-altitude populations, and even some flavonoids, such as quercetin, tricin, and rutin, were significantly enriched in low-altitude populations. This phenomenon indicated that the accumulation of flavonoids was not a result of local adaptation to high altitude. Furthermore, association analysis with in situ environmental variables showed that the contents of quercetin, tricin, and rutin were strongly positively correlated with latitude, longitude, and precipitation gradients and negatively correlated with temperature gradients. Thus, we could conclude that the accumulations of flavonoids in A. squarrosum were more likely as a result of local adaption to environmental heterogeneity combined with precipitation and temperature other than high altitude. This study not only provides an example to understand the molecular ecological basis of pharmacognosy, but also supplies methodologies for developing a new industrial crop with ecological and agricultural importance.

Evolution and functional differentiation of recently diverged phytochelatin synthase genes from Arundo donax L

Phytochelatin synthases (PCSs) play pivotal roles in the detoxification of heavy metals and metalloids in plants; however, little information on the evolution of recently duplicated PCS genes in plant species is available. Here we characterize the evolution and functional differentiation of three PCS genes from the giant reed (Arundo donax L.), a biomass/bioenergy crop with remarkable resistance to cadmium and other heavy metals. Phylogenetic reconstruction with PCS genes from fully sequenced monocotyledonous genomes indicated that the three A. donax PCSs, namely AdPCS1-3, form a monophyletic clade. The AdPCS1-3 genes were expressed at low levels in many A. donax organs and displayed different levels of cadmium-responsive expression in roots. Overexpression of AdPCS1-3 in Arabidopsis thaliana and yeast reproduced the phenotype of functional PCS genes. Mass spectrometry analyses confirmed that AdPCS1-3 are all functional enzymes, but with significant differences in the amount of the phytochelatins synthesized. Moreover, heterogeneous evolutionary rates characterized the AdPCS1-3 genes, indicative of relaxed natural selection. These results highlight the elevated functional differentiation of A. donax PCS genes from both a transcriptional and an enzymatic point of view, providing evidence of the high evolvability of PCS genes and of plant responsiveness to heavy metal stress.

Synteny-based mapping of causal point mutations relevant to sand rice (Agriophyllum squarrosum) trichomeless1 mutant by RNA-sequencing

Sand rice (Agriophyllum squarrosum), a diploid Amaranthaceae species within the order Caryophyllalles, has potential as a crop in response to concerns about ongoing climate change and future food security. Modifying the weedy traits, such as dense trichomes, is important for the domestication of sand rice. In this study, an ethyl methanesulfonate (EMS) induced Agriophyllum squarrosum trichomeless mutant (astcl1) was isolated and genetic analysis revealed that this mutant was controlled by a single recessive locus. RNA sequencing was employed to analyze sequence variations between the mutant and wild-type individuals based on allele frequencies. Synteny-based mapping against two closely related and sequenced species finally delimited the causal mutations into an approximately 18.97 mega base pair (Mb) interval on the top portion of the chromosome 6 (Bv6) of sugar beet (Beta vulgaris) and two intervals (5.56 and 2.69 Mb) on the chromosomes 14 (CqA14) and 6 (CqB06) of quinoa (Chenopodium quinoa). These two quinoa intervals were located in the orthologous chromosomal regions against sugar beet Bv6. Differential expression analysis revealed that the astcl1 mutation only affects the expression of a small subset of genes. Most significantly, 17 lipid transport and metabolism related genes, such as GDSL esterases and very-long-chain 3-oxoacyl-CoA reductase 1, and two R2R3 MYB genes, MYB39 and RAX3, were down-regulated in astcl1 mutants. These results imply that the Astcl1 protein coordinately regulates trichome initiation and cuticle biosynthesis in sand rice.

Application of Molecular Methods in the Identification of Ingredients in Chinese Herbal Medicines

There are several kinds of Chinese herbal medicines originating from diverse sources. However, the rapid taxonomic identification of large quantities of Chinese herbal medicines is difficult using traditional methods, and the process of identification itself is prone to error. Therefore, the traditional methods of Chinese herbal medicine identification must meet higher standards of accuracy. With the rapid development of bioinformatics, methods relying on bioinformatics strategies offer advantages with respect to the speed and accuracy of the identification of Chinese herbal medicine ingredients. This article reviews the applicability and limitations of biochip and DNA barcoding technology in the identification of Chinese herbal medicines. Furthermore, the future development of the two technologies of interest is discussed.

Morphological description and DNA barcoding study of sand rice (Agriophyllum squarrosum, Chenopodiaceae) collected in Kazakhstan

BACKGROUND

Sand rice (Agriophyllum squarrosum (L.) Moq.) is an annual shrub-like plant adapted to the mobile sand dunes in desert and semi-desert regions of Asia. It has a balanced nutrient composition with relatively high concentration of lipids and proteins, which results in its nutrition being similar to legumes. Sand rice's proteins contain the full range of essential amino acids. However, calories content is more similar to wheat. These features together with desert stress resistance make sand rice a potential food crop resilient to ongoing climate change. It is also an important fodder crop (on young stages of growth) for cattle in arid regions of Kazakhstan. In our work, sand rice samples were collected from two distant regions of Kazakhstan as a part of the nation-wide project to determine genetic variation of the native flora.

RESULTS

Samples were collected in western and southeastern parts of Kazakhstan separated by distances of up to 1300 km. Sequences of the nuclear ribosomal DNA ITS1-5.8S-ITS2 region and the chloroplast matK gene confirmed the identity of species defined by morphological traits. Comparison with GenBank sequences revealed polymorphic sequence positions among Kazakh populations and GenBank references, and suggested a distinction among local populations of sand rice. The phylogenetic analysis of nucleotide sequences showed a clear partition of A. squarrosum (L.) Moq. from Agriophyllum minus Fisch. & C.A. Mey, which grows in the same sand dunes environment.

CONCLUSIONS

DNA barcoding analyses of ITS and matK sequences showed a segregation of A. squarrosum from A. minus into separate clades in Maximum-Likelhood dendrograms. ITS analysis can be successfully used to characterize A. squarrosum populations growing quite distant from each other. The data obtained in this work provide the basis for further investigations on A. squarrosum population structure and may play a role in the screening of sand rice plants growing in desert and semi-desert environments of Central Asia and China.